Next Article in Journal
Perspectives of Healthcare Professionals Regarding Factors Associated with Antimicrobial Resistance (AMR) and Their Consequences: A Cross Sectional Study in Eastern Province of Saudi Arabia
Next Article in Special Issue
Successful Intratracheal Treatment of Phage and Antibiotic Combination Therapy of a Multi-Drug Resistant Pseudomonas aeruginosa Murine Model
Previous Article in Journal
Blood Cultures and Appropriate Antimicrobial Administration after Achieving Sustained Return of Spontaneous Circulation in Adults with Nontraumatic Out-of-Hospital Cardiac Arrest
Previous Article in Special Issue
Basics for Improved Use of Phages for Therapy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pseudomonas aeruginosa PAO 1 In Vitro Time–Kill Kinetics Using Single Phages and Phage Formulations—Modulating Death, Adaptation, and Resistance

1
INL-International Iberian Nanotechnology Laboratory, Av. Mestre José Veiga, 4715-330 Braga, Portugal
2
Centre of Biological Engineering, LIBRO—Laboratório de Investigação em Biofilmes Rosário Oliveira, University of Minho, 4710-057 Braga, Portugal
3
Clinical Pathology Department, Hospital de Braga, 4710-243 Braga, Portugal
*
Author to whom correspondence should be addressed.
Submission received: 17 June 2021 / Revised: 2 July 2021 / Accepted: 8 July 2021 / Published: 19 July 2021
(This article belongs to the Special Issue Phage Therapy to Control Pathogenic Bacteria)

Abstract

:
Pseudomonas aeruginosa is responsible for nosocomial and chronic infections in healthcare settings. The major challenge in treating P. aeruginosa-related diseases is its remarkable capacity for antibiotic resistance development. Bacteriophage (phage) therapy is regarded as a possible alternative that has, for years, attracted attention for fighting multidrug-resistant infections. In this work, we characterized five phages showing different lytic spectrums towards clinical isolates. Two of these phages were isolated from the Russian Microgen Sextaphage formulation and belong to the Phikmvviruses, while three Pbunaviruses were isolated from sewage. Different phage formulations for the treatment of P. aeruginosa PAO1 resulted in diversified time–kill outcomes. The best result was obtained with a formulation with all phages, prompting a lower frequency of resistant variants and considerable alterations in cell motility, resulting in a loss of 73.7% in swimming motility and a 79% change in swarming motility. These alterations diminished the virulence of the phage-resisting phenotypes but promoted their growth since most became insensitive to a single or even all phages. However, not all combinations drove to enhanced cell killings due to the competition and loss of receptors. This study highlights that more caution is needed when developing cocktail formulations to maximize phage therapy efficacy. Selecting phages for formulations should consider the emergence of phage-resistant bacteria and whether the formulations are intended for short-term or extended antibacterial application.

1. Introduction

Pseudomonas aeruginosa is an opportunistic pathogen and a leading cause of severe nosocomial infections [1,2]. P. aeruginosa infections in healthcare settings include pneumonia, surgical, wound, and urinary tract infections and are associated with risk factors such as mechanical ventilation, immunosuppression, catheterization [2,3,4,5]. Effective treatment of P. aeruginosa infections is challenging [6], and multidrug-resistant and extensively drug-resistant P. aeruginosa have increased in prevalence and augmented morbidity and mortality [2,7,8]. P. aeruginosa is listed since 2017 as one of the critical priority pathogens listed by the World Health Organization to encourage research and development into new antibacterials [9].
Bacteriophage (phage) therapy is in use under the umbrella of article §37 (Unproven Interventions in Clinical Practice) of the Helsinki Declaration of Ethical Principles for Medical Research Involving Human Subjects [10,11]. In addition, the compassionate treatment of patients using phages is set out by different regulatory agencies, such as the Food and Drug Administration and the European Medicines Agency [12]. Commercial phage formulations are available in Russian and Georgian pharmacies and target a particular species or a panoply of species. P. aeruginosa phages are present in the Pyo- and Intesti-bacteriophage (Eliava BioPreparations, <105 plaque forming units (PFU)/mL), and in the complex Pyobacteriophage, Sextaphage, and Intesti-bacteriophage preparations (Microgen, at unspecified concentrations). Although each formulation target is acknowledged, the number of phages for each is not detailed. Metagenomics of Russian and Georgian pyophage cocktails evidenced the existence of 41 full-length phage genomes. In both formulations, P. aeruginosa phages shared homology with the Luz24virus, but diverged in the other contributing phages [13]. Three PYO phage cocktails, produced from 1997 to 2014, included 30 (PYO 1997) and 29 (PYO 2000) draft genomes [14]. Only 11 were common in all preparations, evidencing that the cocktail formulations are not stagnant through times but instead tweaked to fit the predominant strains isolated.
In this work, we fully characterized five P. aeruginosa PAO1 phages and evaluated the antimicrobial efficacy of different formulations, assessing their impact on killing, the emergence of insensitive mutants, and evaluated the motility and virulence changes.

2. Results

2.1. Selection of Phages Based on Host Range

In this work, we tested 30 phages against P. aeruginosa PAO1 and clinical isolates. Isolates were collected from hospitalized patients with different infections. The patients were aged from 3 (male, ear infection) to 90 (female, urinary tract infection) years old. The majority were isolated from blood (42.1%), followed by urine (26.3%) and sputum (21.1%). The antimicrobial susceptibility of the clinical isolates varied from sensitive to all antibiotics (I500546) up to resistant to six antibiotics (I93488) (Table 1). All isolates were sensitive to aztreonam, cefepime, colistin sulfate, and ticarcilin with clavulanic acid. In addition, maximum resistance was observed for seven different isolates to gentamicin and ciprofloxacin.
Phages’ lytic spectra varied from 15 to 55%. Lysis from without was perceived following tests in two urine isolates (U572569, I97824) and one blood isolate (I60026) (Table 2). In these clinical isolates, lysis occurred due to the interaction of multiple phages with the outer membrane of the isolates and not due to a lytic infection cycle and consequent release of progeny phages. Four clinical isolates (I499131, I29074, I41151, and U14706) were utterly insensitive to all phages (Table 2). These four did not share similarities in antibiotic susceptibilities changing in resistance towards one antibiotic (aminoglycosides) up to four antibiotic classes (isolate I41151) (Table 1).
Phages were selected were based on differences in the spectrum of activity. Phages SPCB and SPCG, included in the Sextaphage formulation, lysed 50 and 55% of the tested isolates. The spectra showed SPCB’s ability to infect C364224, I60026, and I60584, and SPCG’s killing isolates I97824 and C80117. SMS12, SMS21, and SMS29 isolated from raw sewage, lysed between 30 and 45%. SMS12 infected U570696 that both Sextaphages did not lyse. SMS21 also killed this isolate and further lysed H73832; however, it could not infect I202628.

2.2. Virion Particle and Plaque Morphologies

The characteristics of two phages previously isolated [15] and three additional phages were studied in terms of virion and plaque morphologies (Figure 1). Phages SPCB and SPCG had short tails resembling members of the Autographiviridae family. On the other hand, phages SMS12, SMS21, and SMS29 had long contractile tails resembling phages of the Myoviridae family.
Plaque morphologies varied from tiny (1 mm) to big (5 mm), and the plaque + halo (p + h) diameters ranged between 3 and 35 mm (Figure 2). SPCB and SPCG had considerably larger plaques (p) which did not alter with time, and haloes that increased, since the start (Figure 2A), at average speeds of 0.41 ± 0.24 and 0.43 ± 0.18 cm/day until 168 h, respectively (Figure 2C). In addition, the lysis zones of SPCB and SPCG had significant numbers of colonies already after 24 h. The haloes of SMS12, SMS21, and SMS29 were smaller (p + h, Figure 2B) and started to increase after 96 h at average speeds of 0.51 ± 0.33, 0.51 ± 0.48, and 0.69 ± 0.35 cm/day (Figure 2D), respectively.

2.3. One-Step Growth Characteristics

One-step growth experiments with the five phages were performed (Figure 3).
Phages SPCB and SPCG reached 141.5 ± 70.5 and 169.0 ± 52.9 PFU/infected bacteria having similar latent periods (15 min). The other three phages resulted in burst sizes of 138.5 ± 40.3 (SMS12), 209.0 ± 41.0 (SMS21), and 199.0 ± 25.5 (SMS29) PFU/infected bacteria, respectively. Despite this resemblance in burst size, SMS29 presented a 10 min shorter latent period than SMS12 and SMS21 (Figure 3).

2.4. Phage Genomes and Comparative Analysis

SPCB and SPCG resembled Phikmvviruses from the Krylovirinae sub-family of the Autographiviridae family. Pairwise identity at the nucleotide level was 97.1%, with minor differences in the regions between ORFs 4–6 of SPCG and ORF 25 of SPCB (Figure 4A, Table S1). In addition, both SPCB and SPCG showed homology to phages PT5 (EU056923) and vB_Pae_TbilisiM32 (NC_017865) but at different percentages (Supplementary Tables S2 and S3 and Supplementary Figures S1 and S2).
Pairwise identity, at the nucleotide level, between SMS12, SMS21, and SMS29, was above 98%, with only a few unique regions (Figure 4B, Table S1). The most considerable differences started following the DNA primase (ORF_0074 in phages SMS21 and SMS2, ORF_0075 in SMS12) until ORF_0082. All three phages resembled Pbunaviruses, more specifically JG024 (NC_017674, and GU815091), phage 14-1 (NC_011703), respectively (Supplementary Tables S4–S6 and Supplementary Figures S3–S5).
No transmembrane domains and no tRNAs were present in the genomes, indicating their sole dependence on the host tRNA molecules. Promoter numbers varied between 1 and 3. Rho-independent terminators found no terminators for SPCB, one for SPCG, eight for SMS12, seven for SMS21, and six for SMS29, respectively.

2.5. Time–Kill of Single Phage and Phage Cocktail Formulations

Time–kill experiments were performed with single phages or phage formulations (Figure 5A,B).
Time–kill experiments showed that Pbunaviruses (SMS12, SMS21, and SMS29) produced the best antibacterial effect (reducing approximately 4 log10 CFU/mL after 3 h). The Phikmvviruses SPCB and SPCG led to 3 log10 CFU/mL reductions (Figure 5A). After 3 h post-infection, a growth of P. aeruginosa cells was observed, which was faster for Phikmvviruses than Pbunaviruses. All phages, except SPCB, exhibited after 24 h post-infection significantly lower (p ≤ 0.05) viable cell counts compared to the controls.
The number of phages used in the multi-phage combination experiments (Figure 5B) varied from 2 to 5 phages. The five-phage cocktail formulation (5PCF) led to a faster lysis at 1 h (4.8 log10 CFU/mL reduction) that further increased after 3 h (6.4 log10 CFU/mL) (Figure 5B). Removing just phage SPCB from this 5PCF reduced the antibacterial efficacy of the formulation. In fact, this four-phage cocktail formulation (4PCF) (SPCG + SMS12 + SMS21 + SMS29), and the two-phage formulations SMS12 + SMS21, SMS12 + SMS29, and SMS21 + SMS29 showed similar CFU/mL reductions (p > 0.05), oscillating between 4.0 and 4.5 log10 at 3 h. The combination of the two Sextaphage Phikmvviruses caused the lowest reduction (2.7 log10 at 3 h). A statistically significant CFU/mL difference (p ≤ 0.05), after 24 h, was observed between the control and the following cocktails: the three-phage cocktail formulation (3PCF) (SMS12 + SMS21 + SMS29, 1.7 log10), SMS12 + SMS29 (1.6 log10), and SMS12 + SMS21 (1.5 log10). As in single-phage experiments, a rapid increase in cells after a few hours of treatment also took place in all combination assays.
Due to the increase in P. aeruginosa after being challenged following a single and a multi-phage approach, surviving cells were recovered post-infection. Their susceptibility to all phages and possible changes in motility were evaluated to understand this growth phenomenon.

2.6. Assessment of the Survivor’s Susceptibility and Motility

Survivors were isolated following 24 h post-infection and assessed for their susceptibility to all phages used in this study (Table 3). Surviving cells were challenged with all phages and grouped into 18 specific resistance patterns, where R stands for resistance and S for susceptible to a given phage (Table 3). Even though we registered different susceptibility patterns following phage treatment, some recovered survivors continued to show susceptibility to the phage used but had acquired resistance towards other phages (Table 3). These susceptible survivors are highlighted in green, and an example of this is what was observed with phage SMS12. For example, survivors of SMS12 treatment showing patterns 7, 8, 10, 12, and 15, continued to be susceptible to this phage but had become resistant to one (patterns 12 and 15), two (patterns 8 and 10), or three (pattern 7) phages, respectively. It is also worth mentioning that P. aeruginosa PAO1 is susceptible to all phages, and these replicate to produce progeny. However, after 24 h, colonies obtained from the control samples (not challenged with any phage) presented resistance towards phages SPCB and SPCG (patterns 15—SRSSS (30%), and 8—RRSSS (30%)).
A higher number of patterns showed that survivors had acquired resistance towards SPCB (79.6%) and SPCG (64.7%) (see the last two rows in the susceptibility profile columns). The 5PCF resulted in the highest number of patterns (nine patterns), followed by treatment with SPCB and SMS12 alone (six). Conversely, the least varying susceptibility profiles were achieved with individual phage treatments using SPCG, SMS21, and SMS29, and the two-phage SMS12 + SMS29 formulation. Although 5PCF presented the most patterns, the removal of a phage SPCB (4PCF) and the two Phikmvviruses phages (3PCF) increased the percentage of bacteria resistant to four (R4) and five (R5) phages. R4 phage survivors after 5PCF were 30% (sum of patterns 2, 3, 5, 9) while after 4PCF and 3PCF, these increased to 57.1% (4PCF, sum of patterns 2, 3, and 5), and 50% (3PCF, patterns 2 plus 5), respectively. Additionally, R5 phage survivors (pattern 1) increased from 10% in 5PCF to 14.3% (4PCF) and 25% (3PCF) (Table 3). The removal of all Pbunaviruses from the 5PCF increased R5 phage survivors from 10% (5PCF) to 33.3% (SPCB + SPCG), respectively.
The results between any two-phage formulation comprising two of the following phages: SMS12, SMS21, and SMS29, vary considerably compared to the 3PCF consisting of the same phages. For instance, in terms of R5, the two-phage formulations had either lower (7.5%, SMS12 + SMS21), equal (25%, SMS21 + SMS29), or higher percentages (80%, SMS21 + SMS29).
In addition to the susceptibilities of surviving cells towards the different phages, survivors were also evaluated for possible changes in motility (Table 4). Several motility differences were registered, particularly considering the predominant swimming and swarming motilities (blue and grey highlighted values). Although the control cells (non-phage challenged) showed changes in susceptibility to phages, 100% of the survivors maintained excellent swimming characteristics and dendritic swarming motility. Additionally, survivors from single SMS21 and SMS29 phage experiments remained mostly good swimmers (87.5 and 79%, respectively), while those following treatment with phages SPCB, SPCG, and SMS12 changed their swimming predominance to nonswimmers. Compared to the non-phage-treated survivors, changes in swarming were only perceived for single phage treatments with SMS12 and SMS29.
The use of phage cocktails also resulted in motility shifts compared to non-phage exposed P. aeruginosa PAO1 (Table 4). Only survivors from SPCB + SPCG remained good swimmers (50%), although in a fairly similar amount to nonswimmers (45.5%). Therefore, the other cocktail formulations will be compared, regarding swimming characteristics, only in terms of the percentage changes in nonswimming survivors. The 5PCF treatment resulted in 73.7% of nonswimmers and a predominant smooth edge swarming (79.0%). The removal of phage SPCB from the 5PCF decreased nonswimmers to 38.4% (4PCF). However, this removal increased nonswarmers from 15.8% (5PCF) to 57.1% (4PCF) (Table 4). The removal of both SPCB and SPCG from the 5PCF caused a decrease in nonswimmers (42.9%, 3PCF) and increased nonswarmers (87.5%, 3PCF).
Comparing the 3PCF (SMS12 + SMS21 + SMS29) survivors and any of the two-phage formulations prepared with these phages showed an increase in nonswimmers from 42.9% (3PCF) to 71.4% (SMS12 + SMS21), 55.6% (SMS12 + SMS29), and 60.0% (SMS21 + SMS29), respectively. Additionally, the predominant swarming motility changed considerably from dendritic (87.5%, 3PCF) to nonswarmers (SMS12 + SMS21, 71.4%), with a smooth edge (SMS12 + SMS29, 57.1%). Only SMS21 + SMS29 kept the same swarming motility as the 3PCF, although this was reduced from 87.5 to 50.0%.
A few examples of swimming and swarming profiles are present in Figure 6 and Figure 7. In terms of swimming, some survivors had no motility (Figure 6A) while others presented reduced (Figure 6B,C), moderate (Figure 6D–F), and excellent swimming competencies (Figure 6G–I). A few survivors also produced flares that swam beyond the uniform swimming zone (Figure 6J,K), and a minor fraction presented, after 48 h, a red/brown pigmented area indicative of pyorubrin production (Figure 6I).
Differences in swarming behaviors were perceived, particularly in the swarm zone coverage, tendril formation, and alterations related to the swarm zone edges (Figure 7). Some survivors lost their swarming motility (Figure 7A). Furthermore, most cells presented smooth edges with wandering colonies (Figure 7D–F) or multiple fronts radiating outwards that swarmed fast due to the developed tendrils (Figure 7I,J). In addition, some survivors showed deep creases connected to the central swarm zone (Figure 7C,H), while others showed suppressors emerging from the central colony as motility flares (Figure 7B,G).

3. Discussion

The clinical use of phages has witnessed significant advances after several compassionate reports and the successful healing of antibiotic-resistant infections. However, despite the enthusiasm for phage therapy, phage treatments cause the emergence of phage-insensitive phenotypes, which can compromise the therapeutic outcome. A recent review showed that phage-insensitive variants occurred in 80% of studies targeting the intestinal milieu, 50% of studies using sepsis models, and 75% of humans [15].
This study focuses on the impact of different phage formulations on the killing, survival, and resistance of P. aeruginosa. P. aeruginosa is a member of the challenging ESKAPE pathogens group, which shows an excellent ability to “escape” killing by antibiotics [16,17,18]. The virulence of P. aeruginosa has been attributed to several cell-associated factors such as LPS, flagellum, as well as pilus and non-pilus adhesins, and to exoenzymes or different secretory virulence factors [19].
Five P. aeruginosa PAO1 phages were characterized and used alone and in combination to control this strain. The Phikmvviruses (SPCB and SPCG) showed high homology to many Phikmvviruses, with burst sizes (116 (phages phiKMV and PT2) to 225 (phage LKA1) PFU/infected cell), and genome lengths in agreement with those reported (41,593 bp (phage LKA1)–43,548 bp (LUZ19)) [20]. The three other phages were Myoviruses, resembling other Pbunaviruses. Their genome lengths agreed with other Pbunaviruses (64,423 bp for phage LBL3 up to 68,871 bp for phage vB_Pae_PS44). However, SMS21 and SMS29 had slightly superior burst sizes to those reported (64 (phage KT28) [21] to 180 (phage JG024) PFU/infected bacteria [20,22]), with 209 and 194 PFU/infected bacteria, while SMS12 presented a burst in that range (138 PFU/infected bacteria). In addition, SMS12 and SMS21 had latent periods of 35 min, which are typically observed with other Pbunaviruses (35 min (KTN6 and KT28) [23] up to 50 min (phage JG024) [22]).
Single phage application significantly reduced the living cell population, and the best antibacterial efficacy was perceived 3 h post-treatment with the Pbunaviruses (Figure 5A). Despite Phikmvviruses’ (SPCB and SPCG) shorter latent period, this factor did not significantly outweigh the antibacterial performance of other phages. In fact, SPCB and SPCG caused the lowest viable cell reductions and the most rapid increase in survivors. On the other hand, Myoviruses use in the experiments maintained the viable cell counts low until 7 h of treatment.
In terms of combinatory phage experiments, only six phage combinations were tested (Figure 5B). The selection of the combinations tested relied on the limited antibacterial action of the dual combination of SPCB + SPCG. Different phage formulations gave rise to startling results at the population level. Cell death was much faster after one hour and may hypothetically be due: (i) to lysis-from-without phenomena; or (ii) cooperation between Autographiviridae and Myoviridae, which use different host cell receptors for adsorption (discussed below). Lysis was not further enhanced between 1 and 3 h of treatment (Figure 5B), as had previously occurred in single-phage treatments (Figure 5A). The loss of further activity may be due to competition between the access of the phages to the host receptors. Although the multiplicity of infection remains the same, taking into account each phage and number of bacteria present in the culture, the combination of phages in a cocktail can double the number of phages targeting a specific receptor when the two Phikmvviruses or two of the three Pbunaviruses are combined or even triple the phages available for adsorption to a receptor when SMS12, SMS21, and SMS29 are combined. Competition for receptors in the cell wall delays adsorption and concomitantly delays phage lysis [24,25], as recently demonstrated using fluorescently labeled phages [24]. In that work, the authors showed observed that the lysis only derived due to one phage relative to the mixed lysis fluorescence, suggesting a direct or indirect suppression of one of the phages at some point. Phage dominance was due to the blocking of DNA replication through resource sequestration, and the impotence of phage was due to an incomplete ejection of the DNA into the cell.
The selective pressure from phages on their host population can result in potential alterations in phage receptors that may hinder the phage adsorption step [26,27] and change the virulence of the emerged variants [28]. Many have reported lower frequencies of phage-resistant mutants using cocktail formulations rather than monotherapy against Klebsiella pneumoniae [29] and Escherichia coli [30,31]. The results presented herein are in agreement with these articles. The 4PCF and 3PCF gave higher percentages of phage-resistant bacteria compared to 5PCF. Nevertheless, this was not universal, since this hypothesis did not apply to all two-phage formulations that were prepared using the same phages used to produce the 3PCF.
Survivors from control experiments became insensitive to two phages, but this can be due to the phenotypic variations of P. aeruginosa PAO1 itself that are known to occur at a relatively high frequency [32].
Phages can use different host receptors for adsorption, such as regions of the LPS (O-antigen (phage P22 [33]), outer core (phage SSU5 [34])); outer membrane proteins (e.g., FepA and TonB (phage H8 [35]), OmpC (phage T4 [36])); type IV pili (phages DLP1 and DLP2 [37]), flagella (phage 7-7-1 [38]); among others. Phages related to phiKMV, such as SPCB and SPCG, are pili-dependent [20,39], and Pbunaviruses, such as phages SMS12, SMS21, and SMS29 are LPS-dependent [39].
P. aeruginosa possesses two surface structures, a single polar flagellum (flagella) and a polar type IV pili (TFP), that facilitate its motility [40], and both can serve as phage receptors. Flagella are mostly virulence factors for the establishment, persistence, and inflammatory profile and a common cause of acute and chronic of P. aeruginosa infections [23]. Flagella are not permanent cellular structures; instead, the cell’s probability of having a flagellum differs across different growth phases [41]. Bacteria without flagella generally cause less inflammation and mortality [42,43]. The swimming motility of the survivors varied according to the phage(s) used, and higher losses in motility were achieved in treatments with the 5PCF and the SMS12 + SMS21 cocktails, presenting 73.7 and 71.4% of nonswimmers, respectively. In theory, these nonswimmers may be defective in flagella, and, as a positive consequence, these cells will have a decline in virulence and be less prone to form biofilms on surfaces and tissues.
The swarming of P. aeruginosa PAO1 is typically characterized by a dendritic colonial appearance. However, when changes in flagellar quantity and placement or both occur, swarming motility can be compromised. Individual Myoviruses, the 3PCF, and SMS12 + SMS21 revealed the highest losses in swarming motility. Diminished or absent swarming motility can result from a loss of the signal recognition particle-like protein FlhF, resulting in the assembly of flagella at nonpolar locations on the cell resulting in defective swimming and swarming motilities in P. aeruginosa [44]. Additionally, mutations in LPS can influence the exposition of flagella and pili on the bacterial surface, such as that observed with PAO1 ΔwaaL mutants which encoded a functional O-antigen ligase, which showed drastic alterations in swimming and twitching motilities [45,46]. This may justify the loss of infectivity of phages in some of these experiments performed with the Myoviruses.

4. Materials and Methods

4.1. Bacteria, Phages, Growth Conditions

This work used P. aeruginosa PAO1 (DSM 22644) and different clinical isolates (Table 1). All strains were grown at 37 °C in liquid tryptic soy broth (TSB) or solid TSB medium (TSB + 1.2% (w/v) of agar). The Pseudomonas phages used were isolated from the Sextaphage Cocktail (SPC) (Microgen, ImBio Nizhny Novgorod, Russia, series N789 date of issue 03/2016) [47] and also from raw sewage, as previously described [48].

4.2. Phage Host Range Determination

Phages were tested against a panel of isolates using the standard spot test, with phages being serially diluted in SM buffer (5.8 g/L NaCl, 2 g/L MgSO4.7H2O, and 50 mL/L of 50 mM Tris/HCl (pH 7.5)) to investigate lysis from within and from without phenomena [49].

4.3. Phage Propagation and Titration

Phage amplification was performed using the plate lysis and elution method [50], and phage titrations were performed according to a previously described method [51].

4.4. Phage Plaque Morphology and Replication Characteristics

Ten different plaques were analyzed in terms of plaque and halo widths using a high-performance imaging apparatus (Chemi XT4, GBOX-CHEMI-XT4-E, AlphaMetrix Biotech, Rödermark, Germany), coupled with a 4.2 MP imaging 16-bit CCD camera. Phage growth parameters were determined, as previously described [48]. In brief, for single-step experiments, 10 mL of a mid-exponential-phase culture was harvested by centrifugation (7000× g, 5 min, 4 °C) and the pellet resuspended in 5 mL fresh TSB (OD600 of 1.0). Phages (5 mL) were added at a MOI of 0.001, homogenized and allowed to adsorb for 5 min at room temperature. The samples were centrifuged (7000× g, 5 min, 4 °C) and the pellet resuspended in 10 mL of fresh TSB. Samples were taken during a period of 70 min. The samples were serially 10-fold diluted and plated immediately. Plaque forming units (PFU) were determined following 16 h incubation at 37 °C.

4.5. Transmission Electron Microscopy (TEM) Analysis of Phages

Phage particles were centrifuged for one hour (25,000× g, 4 °C) and the pellet was washed twice in tap water and centrifuged once more (1 h, 25,000× g, 4 °C). Phages were deposited on copper grids (400 mesh), stained with 2% (w/v) uranyl acetate (pH 4.0), and imaged at 200 kV using a JEM-2100-HT electron microscope (JEOL, Tokyo, Japan).

4.6. Phage DNA Extraction, Genome Sequencing, and Annotation

Phage DNA was extracted using the quick-DNA Viral kit (Zymo Research, D3015) according to the manufacturer’s specifications, and sequencing performed on an Illumina HiSeq platform after DNA library preparations using the Illumina Nextera XT library preparation kit. The reads obtained were trimmed to remove adapters, contaminations, or low-quality sequences. Contigs were assembled with PATRIC 3.5.43, running BayesHammer on reads, assembling with Velvet, IDBA and SPAdes, and sorting assemblies by ARAST quality score (https://www.patricbrc.org/, 15 September 2019) and manual inspection. Genomes were auto-annotated using MyRAST [52] and manually checked (Geneious Prime 2019, Biomatters, Newark, NJ, USA). Frameshifts were checked with BLASTX [53], and functions of translated open reading frames (ORFs) were searched using BLASTP [54] (E-value ≤ 10–5) and HHPRED server [55] (September–October 2019). Transmembrane domains were checked using TMHMM [56] and Phobius [57], tRNAs searched using tRNAscan-SE [58], and putative promoters using Promoter 2.0 [59]. ARNold was used to predict Rho-independent terminators [60], and the energy was calculated using Mfold [61]. Whole-genome comparisons using OrthoVenn [62] and Progressive MAUVE were performed [63]. The complete phage genome sequences were submitted to GenBank under the accession numbers: vB_PaeP_SPCB (MN615698); vB_PaeP_SPCG (MN615699); vB_PaeM_SMS12 (MN615700); vB_PaeM_SMS21 (MN615701); and vB_PaeM_SMS29 (MN615702).

4.7. Time–Kill Experiments with Different Formulations

The M26-A document of the Clinical & Laboratory Standards Institute was adopted to carry the time–kill experiments [64]. In brief, P. aeruginosa (1 mL, 5 × 108 CFU/mL) was diluted in 9 mL of TSB, and 100 µL of phage (5 × 109 PFU/mL) or 100 µL of SM buffer (control) were added. The mixtures were incubated at 37 °C (120 rpm), and samples were taken at 0, 1, 3, 5, 7, and 24 h post-infection. Serial 10-fold dilutions of P. aeruginosa cells were performed in saline containing 10 mM ferrous ammonium sulphate. Three independent experiments conducted in triplicate were performed. After 24 h of each independent time–kill experiment, 25 surviving colonies were isolated, and tested for their susceptibility not only to the phage they had been challenged with but also all other phages. The susceptibility assay was performed according to a previously described protocol [65].

4.8. Characterization of the Motility Properties of Survivor Cells

The swimming and swarming motilities of survivors were analyzed as previously described [66,67], and Petri dishes observed using Chemi XT4 coupled with a 4.2 MP imaging 16-bit CCD camera.

4.9. Statistical Analysis

Statistical analysis was performed using two-way ANOVA followed by Tukey’s multiple comparison statistical tests, using GraphPad Prism 6 (GraphPad Software, La Jolla, CA, USA). At least three independent experiments were performed, and the results are presented as mean ± standard deviation (SD). Differences were considered as statistically different if p ≤ 0.05 (95% confidence interval).

5. Conclusions

The results of this study show that universal assumptions regarding the decrease of phage-resistant variants using cocktails are not valid for any given phage combination. Many studies found in literature isolate phages, briefly characterize them, and combine a few in cocktails failing to properly identify the host receptors which they target. Therefore, caution is necessary when combining phages from the same genus due to potentially adverse outcomes. Ideally, the phage formulations should combine different genera, more than just the two used in this study, to better understand the phage-phage and phage-bacteria interactions and produce better antibacterial solutions than the commonly available ones.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/antibiotics10070877/s1, Table S1: Clusters and number of singletons formed using OrthoVenn2 by the different phages used in this work, Table S2: Percentage of identity after MAFFT alignment (65% similarity (5.0/-4.0)) and tree building using PHYML (Hasegawa-Kishino-Yano substitution model, no bootstrapping/likelihood) of phage SPCB with other homologous phages, Table S3: Percentage of identity after MAFFT alignment (65% similarity (5.0/-4.0)) and tree building using PHYML (Hasegawa-Kishino-Yano substitution model, no bootstrapping/likelihood) of phage SPCG with other homologous phages, Table S4: Percentage of identity after MAFFT alignment (65% similarity (5.0/-4.0)) and tree building using PHYML (Hasegawa-Kishino-Yano substitution model, no bootstrapping/likelihood) of phage SMS12 with other homologous phages, Table S5: Percentage of identity after MAFFT alignment (65% similarity (5.0/-4.0)) and tree building using PHYML (Hasegawa-Kishino-Yano substitution model, no bootstrapping/likelihood) of phage SMS21 with other homologous phages, Table S6: Percentage of identity after MAFFT alignment (65% similarity (5.0/-4.0)) and tree building using PHYML (Hasegawa-Kishino-Yano substitution model, no bootstrapping/likelihood) of phage SMS29 with other homologous phages, Figure S1: Progressive MAUVE whole-genome alignment (using MUSCLE 3.6) of phage SPCB (upper genome) with five Phikmvviruses based on MAFFT alignment, Figure S2: Figure S2. Progressive MAUVE whole-genome alignment (using MUSCLE 3.6) of phage SPCG (upper genome) with five Phikmvviruses based on MAFFT alignment, Figure S3: Progressive MAUVE whole-genome alignment (using MUSCLE 3.6) of phage SMS12 (upper genome) with various Pbunaviruses, Figure S4. Progressive MAUVE whole-genome alignment (using MUSCLE 3.6) of phage SMS21 (upper genome) with various Pbunaviruses, Figure S5. Progressive MAUVE whole-genome alignment (using MUSCLE 3.6) of phage SMS29 (upper genome) with various Pbunaviruses.

Author Contributions

Conceptualization, S.S.; methodology, S.S., A.M.P. and A.F.; writing—original draft preparation, A.M.P. and S.S.; writing—review and editing, L.M.P., M.B.-L., S.S.; funding acquisition, M.B.-L. and S.S. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Portuguese Foundation for Science and Technology (FCT) under the scope of the strategic funding of UID/BIO/04469/2020 unit. S.S. acknowledges funding by FCT through the individual scientific employment program contract (2020.03171.CEECIND).

Informed Consent Statement

Patient consent was waived since the bacteria used were only to confirm the spectrum of activity of the isolated phages.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Breidenstein, E.B.M.; de la Fuente-Núñez, C.; Hancock, R.E.W. Pseudomonas aeruginosa: All roads lead to resistance. Trends Microbiol. 2011, 19, 419–426. [Google Scholar] [CrossRef] [PubMed]
  2. Lister, P.D.; Wolter, D.J.; Hanson, N.D. Antibacterial-resistant Pseudomonas aeruginosa: Clinical impact and complex regulation of chromosomally encoded resistance mechanisms. Clin. Microbiol. Rev. 2009, 22, 582–610. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Bassetti, M.; Vena, A.; Croxatto, A.; Righi, E.; Guery, B. How to manage Pseudomonas aeruginosa infections. Drugs Context 2018, 7, 212527. [Google Scholar] [CrossRef]
  4. Serra, R.; Grande, R.; Butrico, L.; Rossi, A.; Settimio, U.F.; Caroleo, B.; Amato, B.; Gallelli, L.; De Franciscis, S. Chronic wound infections: The role of Pseudomonas aeruginosa and Staphylococcus aureus. Expert Rev. Anti. Infect. Ther. 2015, 13, 605–613. [Google Scholar] [CrossRef] [PubMed]
  5. Mulcahy, L.R.; Isabella, V.M.; Lewis, K. Pseudomonas aeruginosa Biofilms in Disease. Microb. Ecol. 2014, 108, 329–335. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Pang, Z.; Raudonis, R.; Glick, B.R.; Lin, T.J.; Cheng, Z. Antibiotic resistance in Pseudomonas aeruginosa: Mechanisms and alternative therapeutic strategies. Biotechnol. Adv. 2019, 37, 177–192. [Google Scholar] [CrossRef] [PubMed]
  7. Tümmler, B.; Kiewitz, C. Cystic fibrosis: An inherited susceptibility to bacterial respiratory infections. Mol. Med. Today 1999, 5, 351–358. [Google Scholar] [CrossRef]
  8. Oliver, A.; Mulet, X.; López-Causapé, C.; Juan, C. The increasing threat of Pseudomonas aeruginosa high-risk clones. Drug Resist. Updates 2015, 21–22, 41–59. [Google Scholar] [CrossRef] [PubMed]
  9. WHO. Global Priority List of Antibiotic-Resistant Bacteria to Guide Research, Discovery, and Development of New Antibiotics; WHO: Geneva, Switzerland, 2017. [Google Scholar]
  10. Dedrick, R.M.; Guerrero-Bustamante, C.A.; Garlena, R.A.; Russell, D.A.; Ford, K.; Harris, K.; Gilmour, K.C.; Soothill, J.; Jacobs-Sera, D.; Schooley, R.T.; et al. Engineered bacteriophages for treatment of a patient with a disseminated drug-resistant Mycobacterium abscessus. Nat. Med. 2019, 25, 730–733. [Google Scholar] [CrossRef]
  11. Wright, A.; Hawkins, C.H.; Änggård, E.E.; Harper, D.R. 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. [Google Scholar] [CrossRef]
  12. McCallin, S.; Sacher, J.C.; Zheng, J.; Chan, B.K. Current state of compassionate phage therapy. Viruses 2019, 11, 343. [Google Scholar] [CrossRef] [Green Version]
  13. McCallin, S.; Sarker, S.A.; Sultana, S.; Oechslin, F.; Brüssow, H. Metagenome analysis of Russian and Georgian Pyophage cocktails and a placebo-controlled safety trial of single phage versus phage cocktail in healthy Staphylococcus aureus carriers. Environ. Microbiol. 2018, 20, 3278–3293. [Google Scholar] [CrossRef] [PubMed]
  14. Villarroel, J.; Larsen, M.V.; Kilstrup, M.; Nielsen, M. Metagenomic analysis of therapeutic PYO phage cocktails from 1997 to 2014. Viruses 2017, 9, 328. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Kohler, T.; Curty, L.K.; Barja, F.; Van Delden, C.; Pechere, J.C. Swarming of Pseudomonas aeruginosa is dependent on cell-to-cell signaling and requires flagella and pili. J. Bacteriol. 2000, 182, 5990–5996. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Oechslin, F. Resistance development to bacteriophages occurring during bacteriophage therapy. Viruses 2018, 10, 351. [Google Scholar] [CrossRef] [Green Version]
  17. Peterson, L.R. Bad Bugs, No Drugs: No ESCAPE Revisited. Clin. Infect. Dis. 2009, 49, 992–993. [Google Scholar] [CrossRef] [PubMed]
  18. Payne, D.J.; Gwynn, M.N.; Holmes, D.J.; Pompliano, D.L. Drugs for bad bugs: Confronting the challenges of antibacterial discovery. Nat. Rev. Drug Discov. 2007, 6, 29–40. [Google Scholar] [CrossRef] [PubMed]
  19. Opar, A. Bad bugs need more drugs. Nat. Rev. Drug Discov. 2007, 6, 943–944. [Google Scholar] [CrossRef]
  20. Nadal Jimenez, P.; Koch, G.; Thompson, J.A.; Xavier, K.B.; Cool, R.H.; Quax, W.J. The multiple signaling systems regulating virulence in Pseudomonas aeruginosa. Microbiol. Mol. Biol. Rev. 2012, 76, 46–65. [Google Scholar] [CrossRef] [Green Version]
  21. Ceyssens, P.J.; Glonti, T.; Kropinski, N.M.; Lavigne, R.; Chanishvili, N.; Kulakov, L.; Lashkhi, N.; Tediashvili, M.; Merabishvili, M. Phenotypic and genotypic variations within a single bacteriophage species. Virol. J. 2011, 8, 134. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Danis-Wlodarczyk, K.; Olszak, T.; Arabski, M.; Wasik, S.; Majkowska-Skrobek, G.; Augustyniak, D.; Gula, G.; Briers, Y.; Jang, H.B.; Vandenheuvel, D.; et al. Characterization of the newly isolated lytic bacteriophages KTN6 and KT28 and their efficacy against Pseudomonas aeruginosa biofilm. PLoS ONE 2015, 10, e0127603. [Google Scholar] [CrossRef]
  23. Murray, T.S.; Ledizet, M.; Kazmierczak, B.I. Swarming motility, secretion of type 3 effectors and biofilm formation phenotypes exhibited within a large cohort of Pseudomonas aeruginosa clinical isolates. J. Med. Microbiol. 2010, 59, 511–520. [Google Scholar] [CrossRef] [PubMed]
  24. Garbe, J.; Wesche, A.; Bunk, B.; Kazmierczak, M.; Selezska, K.; Rohde, C.; Sikorski, J.; Rohde, M.; Jahn, D.; Schobert, M. Characterization of JG024, a Pseudomonas aeruginosa PB1-like broad host range phage under simulated infection conditions. BMC Microbiol. 2010, 10, 301. [Google Scholar] [CrossRef] [Green Version]
  25. Trinh, J.T.; Székely, T.; Shao, Q.; Balázsi, G.; Zeng, L. Cell fate decisions emerge as phages cooperate or compete inside their host. Nat. Commun. 2017, 8, 14341. [Google Scholar] [CrossRef]
  26. Refardt, D. Within-host competition determines reproductive success of temperate bacteriophages. ISME J. 2011, 5, 1451–1460. [Google Scholar] [CrossRef] [PubMed]
  27. Wayne, R.; Neilands, J.B. Evidence for common binding sites for ferrichrome compounds and bacteriophage phi 80 in the cell envelope of Escherichia coli. J. Bacteriol. 1975, 121, 497–503. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Hyman, P.; Abedon, S.T. Bacteriophage host range and bacterial resistance. Adv. Appl. Microbiol. 2010, 70, 217–248. [Google Scholar] [CrossRef] [PubMed]
  29. Labrie, S.J.; Samson, J.E.; Moineau, S. Bacteriophage resistance mechanisms. Nat. Rev. Microbiol. 2010, 8, 317–327. [Google Scholar] [CrossRef] [PubMed]
  30. Friman, V.P.; Hiltunen, T.; Jalasvuori, M.; Lindstedt, C.; Laanto, E.; Örmälä, A.M.; Laakso, J.; Mappes, J.; Bamford, J.K.H. High temperature and bacteriophages can indirectly select for bacterial pathogenicity in environmental reservoirs. PLoS ONE 2011, 6, e17651. [Google Scholar] [CrossRef]
  31. Gu, J.; Liu, X.; Li, Y.; Han, W.; Lei, L.; Yang, Y.; Zhao, H.; Gao, Y.; Song, J.; Lu, R.; et al. A method for generation phage cocktail with great therapeutic potential. PLoS ONE 2012, 7, e31698. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. O’Flynn, G.; Ross, R.P.; Fitzgerald, G.F.; Coffey, A. Evaluation of a cocktail of three bacteriophages for biocontrol of Escherichia coli O157:H7. Appl. Environ. Microbiol. 2004, 70, 3417–3424. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Tanji, Y.; Shimada, T.; Yoichi, M.; Miyanaga, K.; Hori, K.; Unno, H. Toward rational control of Escherichia coli O157:H7 by a phage cocktail. Appl. Microbiol. Biotechnol. 2004, 64, 270–274. [Google Scholar] [CrossRef] [PubMed]
  34. Chandler, C.E.; Horspool, A.M.; Hill, P.J.; Wozniak, D.J.; Schertzer, J.W.; Rasko, D.A.; Ernsta, R.K. Genomic and phenotypic diversity among ten laboratory isolates of Pseudomonas aeruginosa PAO1. J. Bacteriol. 2019, 201, e00595-18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Andres, D.; Hanke, C.; Baxa, U.; Seul, A.; Barbirz, S.; Seckler, R. Tailspike interactions with lipopolysaccharide effect DNA ejection from phage P22 particles in vitro. J. Biol. Chem. 2010, 286, 4165–4172. [Google Scholar] [CrossRef] [Green Version]
  36. Kim, M.; Kim, S.; Park, B.; Ryu, S. Core lipopolysaccharide-specific phage SSU5 as an auxiliary component of a phage cocktail for Salmonella biocontrol. Appl. Environ. Microbiol. 2014, 80, 1026–1034. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Rabsch, W.; Ma, L.; Wiley, G.; Najar, F.Z.; Kaserer, W.; Schuerch, D.W.; Klebba, J.E.; Roe, B.A.; Laverde Gomez, J.A.; Schallmey, M.; et al. FepA- and TonB-dependent bacteriophage H8: Receptor binding and genomic sequence. J. Bacteriol. 2007, 189, 5658–5674. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Washizaki, A.; Yonesaki, T.; Otsuka, Y. Characterization of the interactions between Escherichia coli receptors, LPS and OmpC, and bacteriophage T4 long tail fibers. Microbiologyopen 2016, 5, 1003–1015. [Google Scholar] [CrossRef] [Green Version]
  39. McCutcheon, J.G.; Peters, D.L.; Dennis, J.J. Identification and characterization of type IV pili as the cellular receptor of broad host range Stenotrophomonas maltophilia bacteriophages DLP1 and DLP2. Viruses 2018, 10, 338. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Gonzalez, F.; Helm, R.F.; Broadway, K.M.; Scharf, B.E. More than rotating flagella: Lipopolysaccharide as a secondary receptor for flagellotropic phage 7-7-1. J. Bacteriol. 2018, 200, e00363-18. [Google Scholar] [CrossRef] [Green Version]
  41. Lorenz, A.; Preuße, M.; Bruchmann, S.; Pawar, V.; Grahl, N.; Pils, M.C.; Nolan, L.M.; Filloux, A.; Weiss, S.; Häussler, S. Importance of flagella in acute and chronic Pseudomonas aeruginosa infections. Environ. Microbiol. 2019, 21, 883–897. [Google Scholar] [CrossRef] [PubMed]
  42. Zhuang, X.-Y.; Guo, S.; Li, Z.; Zhao, Z.; Kojima, S.; Homma, M.; Wang, P.; Lo, C.-J.; Bai, F. Dynamic production and loss of flagellar filaments during the bacterial life cycle. bioRxiv 2019, 767319. [Google Scholar] [CrossRef] [Green Version]
  43. Feldman, M.; Bryan, R.; Rajan, S.; Scheffler, L.; Brunnert, S.; Tang, H.; Prince, A. Role of flagella in pathogenesis of Pseudomonas aeruginosa pulmonary infection. Infect. Immun. 1998, 66, 43–51. [Google Scholar] [CrossRef] [Green Version]
  44. Patankar, Y.R.; Lovewell, R.R.; Poynter, M.E.; Jyot, J.; Kazmierczak, B.I.; Berwin, B. Flagellar motility is a key determinant of the magnitude of the inflammasome response to Pseudomonas aeruginosa. Infect. Immun. 2013, 81, 2043–2052. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Murray, T.S.; Kazmierczak, B.I. FlhF Is required for swimming and swarming in Pseudomonas aeruginosa. J. Bacteriol. 2006, 188, 6995–7004. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Abeyrathne, P.D.; Daniels, C.; Poon, K.K.H.; Matewish, M.J.; Lam, J.S. Functional characterization of WaaL, a ligase associated with linking O-antigen polysaccharide to the core of Pseudomonas aeruginosa lipopolysaccharide. J. Bacteriol. 2005, 187, 3002–3012. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Oliveira, A.; Sousa, J.C.; Silva, A.C.; Melo, L.D.R.; Sillankorva, S. Chestnut honey and bacteriophage application to control Pseudomonas aeruginosa and Escherichia coli biofilms: Evaluation in an ex vivo wound model. Front. Microbiol. 2018, 9, 1725. [Google Scholar] [CrossRef] [Green Version]
  48. Sillankorva, S.; Neubauer, P.; Azeredo, J. Isolation and characterization of a T7-like lytic phage for Pseudomonas fluorescens. BMC Biotechnol. 2008, 8, 80. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Milho, C.; Andrade, M.; Vilas Boas, D.; Alves, D.; Sillankorva, S. Antimicrobial assessment of phage therapy using a porcine model of biofilm infection. Int. J. Pharm. 2019, 557, 112–123. [Google Scholar] [CrossRef] [Green Version]
  50. Sambrook, J.; Russel, D. Molecular Cloning: A Laboratory Manual, 4th ed.; Cold Spring Harbor Laboratory Press: Cold Spring Harbor, NY, USA, 2001. [Google Scholar]
  51. Adams, M. Bacteriophages; Interscience Publishers: New York, NY, USA, 1959; ISBN 9781588296825. [Google Scholar]
  52. Aziz, R.K.; Bartels, D.; Best, A.; DeJongh, M.; Disz, T.; Edwards, R.A.; Formsma, K.; Gerdes, S.; Glass, E.M.; Kubal, M.; et al. The RAST server: Rapid annotations using subsystems technology. BMC Genom. 2008, 9, 75. [Google Scholar] [CrossRef] [Green Version]
  53. Altschul, S.F.; Madden, T.L.; Schäffer, A.A.; Zhang, J.; Zhang, Z.; Miller, W.; Lipman, D.J. Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res. 1997, 25, 3389–3402. [Google Scholar] [CrossRef] [Green Version]
  54. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  55. Söding, J.; Biegert, A.; Lupas, A.N. The HHpred interactive server for protein homology detection and structure prediction. Nucleic Acids Res. 2005, 33, W244–W248. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Käll, L.; Sonnhammer, E.L.L. Reliability of transmembrane predictions in whole-genome data. FEBS Lett. 2002, 532, 415–418. [Google Scholar] [CrossRef] [Green Version]
  57. Käll, L.; Krogh, A.; Sonnhammer, E.L.L. A combined transmembrane topology and signal peptide prediction method. J. Mol. Biol. 2004, 338, 1027–1036. [Google Scholar] [CrossRef]
  58. Schattner, P.; Brooks, A.N.; Lowe, T.M. The tRNAscan-SE, snoscan and snoGPS web servers for the detection of tRNAs and snoRNAs. Nucleic Acids Res. 2005, 33, W686–W689. [Google Scholar] [CrossRef] [PubMed]
  59. Knudsen, S. Promoter2.0: For the recognition of PolII promoter sequences. Bioinformatics 1999, 15, 356–361. [Google Scholar] [CrossRef] [PubMed]
  60. Naville, M.; Ghuillot-Gaudeffroy, A.; Marchais, A.; Gautheret, D. ARNold: A web tool for the prediction of Rho-independent transcription terminators. RNA Biol. 2011, 8, 11–13. [Google Scholar] [CrossRef] [Green Version]
  61. Zuker, M. Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31, 3406–3415. [Google Scholar] [CrossRef]
  62. Wang, Y.; Coleman-Derr, D.; Chen, G.; Gu, Y.Q. OrthoVenn: A web server for genome wide comparison and annotation of orthologous clusters across multiple species. Nucleic Acids Res. 2015, 43, W78–W84. [Google Scholar] [CrossRef] [PubMed]
  63. Darling, A.E.; Mau, B.; Perna, N.T. ProgressiveMauve: Multiple genome alignment with gene gain, loss and rearrangement. PLoS ONE 2010, 5, e11147. [Google Scholar] [CrossRef] [Green Version]
  64. Barry, A.L.; Craig, W.A.; Nadler, H.; Reller, L.B.; Sanders, C.C.; Swenson, J.M. M26-A Methods for Determining Bactericidal Activity of Antimicrobial Agents, Approved Guideline, 1st ed.; Clinical and Laboratory Standards Institute: Annapolis Junction, MD, USA, 1999. [Google Scholar]
  65. Moons, P.; Faster, D.; Aertsen, A. Lysogenic conversion and phage resistance development in phage exposed Escherichia coli biofilms. Viruses 2013, 5, 150–161. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Vater, S.M.; Weiße, S.; Maleschlijski, S.; Lotz, C.; Koschitzki, F.; Schwartz, T.; Obst, U.; Rosenhahn, A. Swimming behavior of Pseudomonas aeruginosa studied by holographic 3D tracking. PLoS ONE 2014, 9, e87765. [Google Scholar] [CrossRef] [PubMed]
  67. Ha, D.G.; Kuchma, S.L.; O’Toole, G.A. Plate-based assay for swarming motility in Pseudomonas aeruginosa. In Pseudomonas Methods and Protocols. Methods in Molecular Biology (Methods and Protocols); Humana Press: New York, NY, USA, 2014; Volume 1149. [Google Scholar] [CrossRef]
Figure 1. Virion particle morphology of P. aeruginosa PAO1 phages observed by TEM.
Figure 1. Virion particle morphology of P. aeruginosa PAO1 phages observed by TEM.
Antibiotics 10 00877 g001
Figure 2. Plaque characteristics of the five phages used in this work after 24, 72, 120, and 168 h: (A) morphological characteristics (scale bar 1 cm); (B) microscopy of phages SMS12, SMS21, and SMS29 (5× objective) (scale bar 500 µm); (C) plaque diameters of phages SPCB and SPCG); (D) plaque diameters of phages SMS12, SMS21, and SMS2. Data in (C) and (D) are shown either as mean of plaque diameter (p) ± SD or mean of plaque and halo diameter (p + h) ± SD.
Figure 2. Plaque characteristics of the five phages used in this work after 24, 72, 120, and 168 h: (A) morphological characteristics (scale bar 1 cm); (B) microscopy of phages SMS12, SMS21, and SMS29 (5× objective) (scale bar 500 µm); (C) plaque diameters of phages SPCB and SPCG); (D) plaque diameters of phages SMS12, SMS21, and SMS2. Data in (C) and (D) are shown either as mean of plaque diameter (p) ± SD or mean of plaque and halo diameter (p + h) ± SD.
Antibiotics 10 00877 g002
Figure 3. One-step growth characteristics of the five phages used. Data are shown as mean ± SD.
Figure 3. One-step growth characteristics of the five phages used. Data are shown as mean ± SD.
Antibiotics 10 00877 g003
Figure 4. Comparison of phages from the two families: (A) Progressive MAUVE genome alignment between phages SPCB (above) and SPCG (below); (B) Progressive MAUVE genome alignment between phages SMS12, SMS21, and SMS29. White gap areas show regions that were not aligned and contain sequence elements specific to a particular genome. The height of the similarity profile presents an average level of conservation in a particular region of the genome sequence.
Figure 4. Comparison of phages from the two families: (A) Progressive MAUVE genome alignment between phages SPCB (above) and SPCG (below); (B) Progressive MAUVE genome alignment between phages SMS12, SMS21, and SMS29. White gap areas show regions that were not aligned and contain sequence elements specific to a particular genome. The height of the similarity profile presents an average level of conservation in a particular region of the genome sequence.
Antibiotics 10 00877 g004
Figure 5. Time–kill experiment of each phage alone or in cocktail against P. aeruginosa PAO1. MOI of 10, 120 rpm, 37 °C: (A) phages alone; and (B) phages in different cocktail formulations. 5PCF (phages SPCB + SPCG + SMS12 + SMS21 + SMS29); 4PCF (SPCG + SMS12 + SMS21 + SMS29); 3PCF (SMS12 + SMS21 + SMS29). Data are shown as mean ± SD.
Figure 5. Time–kill experiment of each phage alone or in cocktail against P. aeruginosa PAO1. MOI of 10, 120 rpm, 37 °C: (A) phages alone; and (B) phages in different cocktail formulations. 5PCF (phages SPCB + SPCG + SMS12 + SMS21 + SMS29); 4PCF (SPCG + SMS12 + SMS21 + SMS29); 3PCF (SMS12 + SMS21 + SMS29). Data are shown as mean ± SD.
Antibiotics 10 00877 g005
Figure 6. Different swimming motilities observed in surviving cells: (AI) images acquired using Chemi XT4; and (JL) digital photographs.
Figure 6. Different swimming motilities observed in surviving cells: (AI) images acquired using Chemi XT4; and (JL) digital photographs.
Antibiotics 10 00877 g006
Figure 7. Different swarming motilities observed in surviving cells. (AJ) images acquired using Chemi XT4.
Figure 7. Different swarming motilities observed in surviving cells. (AJ) images acquired using Chemi XT4.
Antibiotics 10 00877 g007
Table 1. Antibiogram of P. aeruginosa clinical isolates used for assessing the lytic spectrum of phages.
Table 1. Antibiogram of P. aeruginosa clinical isolates used for assessing the lytic spectrum of phages.
IDSourceYearSex
(M/F) *
Age
(Years)
Antibiotic
AKTOBCNATMFEPCAZCIPCSIMIMRPPRLTZPTCTTC
I500546Cutaneous2013M828≤1≤18≤240.25≤1228≤284
I499131Sputum2013F59>16≤1>48≤240.25≤1244≤284
C364224Sputum2013M524≤1≤18≤242≤1228≤28≤2
U570696Urine2013M892≤1≤18≤240.25≤14>168884
U572569Urine2013M82>16>4>48≤244≤12>164>328≤2
I66897Blood2018M632≤1≤116≤240.25≤1228≤2324
I29074Blood2017M682≤1≤116≤2161≤12283288
I41151Blood2017M80168816≤222≤1442≤2324
H73832Blood2018M612≤1≤116≤240.25≤1244≤2>164
U88885Blood2017F744≤1≤116≤220.25≤1424≤284
I60026Blood2018F832≤1≤116≤240.25≤1224≤2324
I41152Blood2018F524≤1≤116≤240.25≤1284≤284
I60584Blood2018F732≤1≤116≤240.25≤1428≤2324
I97824Urine2019M882≤188≤2160.25≤11628≤288
I93488Sputum2009M7432888≤2161≤11688≤288
I92986Urine2019F7332≤1168≤222≤1228≤288
C80117Ear2009F316≤188≤240.25≤1224≤284
U14706Urine2009M908≤1≤18≤241≤1448≤284
I202628Sputum2019M682≤1≤18≤220.25≤1228≤284
MIC BreakpointR > 16R > 4R > 4R > 16R > 8R > 8R > 0.5R > 2R > 8R > 8R > 16R > 16R > 16R > 16
* M: male; F: female. Identification of P. aeruginosa clinical isolates using VITEK2, MicroScan WalkAway as well as MALDI-TOF (for strains isolated after 2017). Antibiograms performed using the Kirby–Bauer method using Müeller–Hinton agar according to CLSI rules (for strains isolated until January 2014) and according to EUCAST afterward using VITEK2 and MicroScan WalkAway. Metallo-beta-lactamases search was done using Etest (imipenem/imipenem + EDTA) and the Kirby–Bauer method using Müeller–Hinton agar. Antibiotics and concentrations tested: AK: amikacin (30 µg); TOB: tobramycin (10 µg); CN: gentamicin (30 µg); ATM: aztreonam (30 µg); FEP: cefepime (30 µg); CAZ: ceftazidime (10 µg); CIP: ciprofloxacin (5 µg); CS: colistin sulfate (25 µg); IMI: imipenem (10 µg); MRP: meropenem (10 µg); PRL: piperacillin (30 µg); TZP: piperacillin + tazobactam (36 µg); TC: ticarcillin (75 µg); TTC: ticarcillin + clavulanic acid (85 µg). Antibiotics 10 00877 i002 Isolates showing resistance to a particular antibiotic according to the MIC Breakpoint (EUCAST, http://www.eucast.org/clinical_breakpoints/, accessed 20 January 2020).
Table 2. Lytic spectra against P. aeruginosa PAO1 and different clinical isolates.
Table 2. Lytic spectra against P. aeruginosa PAO1 and different clinical isolates.
Isolation SourcePhageBacterial Isolates
PAO1I5005461499131C364224U570696U572569I66897I29074I41151H73832U88885I60026I51359I60584I97824I93488I92986C80117U14706I202628Efficacy (%) *
SPC (Microgen)SPCA 55
SPCB 55
SPCC 45
SPCE 45
SPCF 55
SPCG 50
Sewage (2013)SMS9 45
SMS10 30
SMS11 20
Sewage (2019)SMS12 35
SMS13 30
SMS14 45
SMS15 30
SMS16 55
SMS17 30
SMS18 20
SMS19 35
SMS20 40
SMS21 30
SMS22 45
SMS23 30
SMS24 35
SMS25 30
SMS26 40
SMS27 15
SMS28 45
SMS29 45
SMS30 45
Antibiotics 10 00877 i001 lysis from within; Antibiotics 10 00877 i002 lysis from without; Antibiotics 10 00877 i003 no lysis detected; Antibiotics 10 00877 i004 phages selected; * percentage of isolates that are susceptible to a given phage.
Table 3. Susceptibility of P. aeruginosa surviving cells recovered after 24 h from the different time–kill experiments with different phage combinations.
Table 3. Susceptibility of P. aeruginosa surviving cells recovered after 24 h from the different time–kill experiments with different phage combinations.
Susceptibility ProfileSusceptibility of Surviving Cells (%)
PatternSPCBSPCGSMS12SMS21SMS29Control5PCF4PCF3PCFSPCB +
SPCG
SMS12 + SMS21SMS12 +
SMS29
SMS21 + SMS29SPCBSPCGSMS12SMS21SMS29
1RRRRR 10.014.325.033.37.580.025.019.1 6.3
2RRRRS 5.028.525.0 12.5
3RRRSR 5.014.3 11.1
4RRRSS 10.0 22.215.0 6.2
5RRSRR 15.014.325.0 20.012.5
6RRSRS 12.5 12.54.8
7RRSSR 11.15.5 12.5 25.0
8RRSSS30.020.0 4.8 18.8
9RSRRR 5.0
10RSSRS 15.0 10.0 9.5 18.8
11RSSSR 14.3
12RSSSS 15.014.3 28.0 42.9100.012.5
13SRRRR 4.8
14SRSRR 12.5 12.5
15SRSSS30.0 22.2 12.514.3
16SSSRS 35.0 12.5100.0100.0
17SSSSS40.0
R *79.664.735.352.947.0100.0100.085.7100.0100.066.5100.087.580.90.012.5100.00.0
S *29.435.364.747.153.00.00.014.30.00.033.50.012.519.1100.071.50.0100.0
Antibiotics 10 00877 i003 surviving P. aeruginosa cells that became resistant; Antibiotics 10 00877 i001 surviving cells that remained susceptible to the phage used in the treatment. When cocktail formulations were used, the value is highlighted in green only if the survivors remained susceptible to all phages present in the specific formulation. * R and S refer to the percentage of survivors showing resistance (R) or susceptibility (S) to a particular phage. 5PCF (phages SPCB + SPCG + SMS12 + SMS21 + SMS29); 4PCF (SPCG + SMS12 + SMS21 + SMS29); 3PCF (SMS12 + SMS21 + SMS29).
Table 4. Swimming and swarming motilities of P. aeruginosa cells surviving 24 h of treatment with different phage formulations.
Table 4. Swimming and swarming motilities of P. aeruginosa cells surviving 24 h of treatment with different phage formulations.
Swimming (%)Swarming (%)
Cells Surviving Specific TreatmentNoReduced to Moderate *Good NoDendriticSmooth EdgeSuppressor
Control 100.0 100.0
SPCB48.611.440.0 50.037.512.5
SPCG50.09.140.98.458.333.3
SMS1244.414.840.766.7 11.122.2
SMS21 12.587.522.266.711.1
SMS2920.010.070.062.525.012.5
5PCF73.7 26.315.85.279.0
4PCF38.430.830.857.142.9
3PCF42.921.435.787.512.5
SPCB + SPCG45.54.550.0 63.627.39.1
SMS12 + SMS2171.414.314.371.428.6
SMS12 + SMS2955.633.311.128.614.357.1
SMS21 + SMS2960.026.713.325.050.025.0
* diameter varying between 0 and 10 mm; diameter above 10 mm; 5PCF (phages SPCB + SPCG + SMS12 + SMS21 + SMS29); 4PCF (SPCG + SMS12 + SMS21 + SMS29); 3PCF (SMS12 + SMS21 + SMS29); Antibiotics 10 00877 i005 predominant swimming pattern; Antibiotics 10 00877 i006 predominant swarming pattern.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pinto, A.M.; Faustino, A.; Pastrana, L.M.; Bañobre-López, M.; Sillankorva, S. Pseudomonas aeruginosa PAO 1 In Vitro Time–Kill Kinetics Using Single Phages and Phage Formulations—Modulating Death, Adaptation, and Resistance. Antibiotics 2021, 10, 877. https://0-doi-org.brum.beds.ac.uk/10.3390/antibiotics10070877

AMA Style

Pinto AM, Faustino A, Pastrana LM, Bañobre-López M, Sillankorva S. Pseudomonas aeruginosa PAO 1 In Vitro Time–Kill Kinetics Using Single Phages and Phage Formulations—Modulating Death, Adaptation, and Resistance. Antibiotics. 2021; 10(7):877. https://0-doi-org.brum.beds.ac.uk/10.3390/antibiotics10070877

Chicago/Turabian Style

Pinto, Ana Mafalda, Alberta Faustino, Lorenzo M. Pastrana, Manuel Bañobre-López, and Sanna Sillankorva. 2021. "Pseudomonas aeruginosa PAO 1 In Vitro Time–Kill Kinetics Using Single Phages and Phage Formulations—Modulating Death, Adaptation, and Resistance" Antibiotics 10, no. 7: 877. https://0-doi-org.brum.beds.ac.uk/10.3390/antibiotics10070877

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop