Characterization of a bacteriophage with broad host range against strains of Pseudomonas aeruginosa isolated from domestic animals

Background Pseudomonas aeruginosa is an opportunistic pathogen and one of the leading causes of nosocomial infections. Moreover, the species can cause severe infections in cystic fibrosis patients, in burnt victims and cause disease in domestic animals. The control of these infections is often difficult due to its vast repertoire of mechanisms for antibiotic resistance. Phage therapy investigation with P. aeruginosa bacteriophages has aimed mainly the control of human diseases. In the present work, we have isolated and characterized a new bacteriophage, named Pseudomonas phage BrSP1, and investigated its host range against 36 P. aeruginosa strains isolated from diseased animals and against P. aeruginosa ATCC strain 27853. Results We have isolated a Pseudomonas aeruginosa phage from sewage. We named this virus Pseudomonas phage BrSP1. Our electron microscopy analysis showed that phage BrSP1 had a long tail structure found in members of the order Caudovirales. “In vitro” biological assays demonstrated that phage BrSP1 was capable of maintaining the P. aeruginosa population at low levels for up to 12 h post-infection. However, bacterial growth resumed afterward and reached levels similar to non-treated samples at 24 h post-infection. Host range analysis showed that 51.4% of the bacterial strains investigated were susceptible to phage BrSP1 and efficiency of plating (EOP) investigation indicated that EOP values in the strains tested varied from 0.02 to 1.72. Analysis of the phage genome revealed that it was a double-stranded DNA virus with 66,189 bp, highly similar to the genomes of members of the genus Pbunavirus, a group of viruses also known as PB1-like viruses. Conclusion The results of our “in vitro” bioassays and of our host range analysis suggested that Pseudomonas phage BrSP1 could be included in a phage cocktail to treat veterinary infections. Our EOP investigation confirmed that EOP values differ considerably among different bacterial strains. Comparisons of complete genome sequences indicated that phage BrSP1 is a novel species of the genus Pbunavirus. The complete genome of phage BrSP1 provides additional data that may help the broader understanding of pbunaviruses genome evolution. Electronic supplementary material The online version of this article (10.1186/s12866-019-1481-z) contains supplementary material, which is available to authorized users.


Background
Pseudomonas aeruginosa is a Gram-negative free-living and parasitic bacterium widespread in the environment. The species is an opportunistic pathogen of great concern to human health due to its broad spectrum of antibiotic resistance and diverse virulence factors. P. aeruginosa is the primary cause of respiratory infection in cystic fibrosis patients and one of the most important causes of severe infections in burn victims [1][2][3]. Moreover, it is a frequent etiological agent of hospital-acquired diseases, such as pneumonia, urinary tract and surgical site infections [4,5].
The versatile metabolic capacity of P. aeruginosa allows it to multiply in a variety of environments and hosts. A study published in 1952 showed a high occurrence in sewage water samples (90.4%) and a much lower incidence in human feces (11%) and soil samples (3%) [6]. Subsequent investigations confirmed the presence of P. aeruginosa in soil [7,8], in human feces [9] and water [10,11]. One of the consequences of the presence of P. aeruginosa in soil and water is that plants are often in contact with these bacteria. In fact, several plant species can be a source of infection with P. aeruginosa [12]. P. aeruginosa is also able to infect distantly related species of animals, such as insects [13], nematodes [14], fishes [15] and birds [16]. Moreover, a variety of domestic animals are susceptible to its infection. For instance, P. aeruginosa may cause endometritis in mares [17,18], mastitis in dairy animals [19][20][21], ear infection, pneumonia, septicemia and enteritis in chinchilla [22] and hemorrhagic pneumonia in confined mink [23][24][25]. In addition, the species is associated with ear, urinary and skin infection in dogs [26][27][28][29][30][31][32]. These infections are often difficult to heal due to multi-antibiotic resistance mechanisms, including the formation of biofilm [33,34].
The difficulties related to the antibiotic treatment of P. aeruginosa and the often damaging nature of its infection have prompted investigations of viruses (bacteriophages) as antimicrobial agents against this species [35,36]. Of prime interest is the possibility of using bacteriophages to treat P. aeruginosa infections in cystic fibrosis patients and burn victims [37,38]. On the other hand, there are much fewer investigations aiming for veterinary applications. One of these studies used antibiotic-resistant strains isolated from dogs to screen phages with potential therapeutic value [39]. Another investigation showed evidence of successful application of bacteriophages in the treatment of canine otitis [40].
Here we present our work with a new P. aeruginosa bacteriophage, named Pseudomonas phage BrSP1. We studied the infectivity and host range of phage BrSP1 against P. aeruginosa strains isolated from domestic animals (including dogs, cattle, swine, horse, and chinchilla) and P. aeruginosa ATCC strain 27853. Our experiments demonstrated that this phage was able to infect 51.4% of the bacterial strains analyzed. The EOP values for phage BrSP1 varied greatly, reaching approximately one hundred fold difference between the strains with the highest and lowest EOPs. In the conditions used in our "in vitro" bioassays, infection with phage BrSP1 maintained low levels of the bacterial population for at least 12 h post-infection. Sequence analysis indicated that phage BrSP1 is a putative new virus species with 66,189 bp, belonging to the genus Pbunavirus of the family Myoviridae. Members of this genus are also referred to as PB1-like viruses [41,42]. Our results indicated that BrSP1 could be part of a phage cocktail to control P. aeruginosa infection in domestic animals. Moreover, the complete genome of phage BrSP1 provided additional data that may help the broader understanding of pbunavirus genome.

Phage isolation, morphology and storage stability
One lytic phage against P. aeruginosa strain Lfar01 was isolated from sewage and a homogeneous stock prepared after three rounds of purification. We named this virus Pseudomonas phage BrSP1. Electron microscopy of phage BrSP1 showed the long tail structure typical of members of the order Caudovirales (Fig. 1). The titer of a Phage BrSP1 stock had a 12.4 fold reduction upon storage, under refrigerated temperature (approximately 5°C), for 104 days (Additional file 1: Table S1).

Host range investigation and EOP analysis
We evaluated host range using, at first, spot test with undiluted phage stock and 37 P. aeruginosa strains, which included 36 strains isolated from domestic animals and P. aeruginosa ATCC 27853 (Table 1 and Additional file 2: Table S2). This preliminary analysis resulted in clear or turbid zones in 24 of the tested strains (Table 1). Subsequently, we inoculated these 24 strains with doses of diluted phage. These later experiments indicated that five of the samples were, in fact, resistant to phage BrSP1 (Table 1). The clear zones produced initially in these five strains were likely due to the disruption of cell membranes [43] or to bacteriocins produced during phage stock production. The formation of individual lysis plaques, due to the dilution of the phage stock, confirmed the susceptibility of 19 strains (51.4%), including P. aeruginosa ATCC 27853 (Table 1). These results, therefore, confirmed the importance of using low titer viral stocks when investigating host range, as has been previously noted [44].
We calculated EOP values as the ratio of lysis plaques produced in each susceptible strain (for a fixed dose of phages) divided by the number of plaques produced in P. aeruginosa strain ATCC 27853 (Fig. 2). In these assays, we inoculated most of the samples with ten microliters of dilution 10 4 of our phage stock. As this dose did not produce plaques on strains CASM03, BOIJ03 and CASM11, we used a dose that was ten times higher (dilution 10 3 ) in their analysis.

Molecular characterization of bacterial strains
The 16S rRNA gene of the strains that were not susceptible to phage BrSP1 infection were PCR amplified and sequenced. Blast analysis of the sequences confirmed the resistant strains were indeed P. aeruginosa (data not shown). The Genbank accession numbers of these sequences are available in Additional file 3: Table S3.
"In vitro" lytic activity of phage BrSP1 and origin of the bacteriophage resistant mutants We investigated the lytic activity of phage BrSP1 in assays with P. aeruginosa strains Lfar01, ATCC 27853 and BOIJ02, a strain with high EOP. In these assays, there was a remarkable reduction in bacterial cell population and bacterial growth remained checked until approximately 12 h p.i. (Fig. 3). At 24 h p.i., however, the cell density of infected and non-infected samples was almost the same. The tenfold difference in the viral concentration used in some of the assays did not produce any noticeable variation (Fig. 3).
In order to understand the nature of the emergence of resistance observed in our assays, we investigated the frequency of occurrence of bacteriophage insensitive mutants (BIMs) in three P. aeruginosa strains. Our results indicated that BIMs arose independently of phage infection at frequencies that ranged from 7.07 × 10 − 9 (minimum frequency for strain BOIJ02) to 2.10 × 10 − 7 (maximum frequency for strain Lfar01; Table 2).

BrSP1 genome size, GC content and accession number
The sequencing of the Pseudomonas phage BrSP1 genome showed that it was 66,189 bp in length with a G + C content of 55.7%. The Genbank accession number of the annotated genome of phage BrSP1 is MF623055.

Gene content analysis and genome features
We identified 94 open reading frames (ORFs; Additional file 4: Table S4) in the genome of phage BrSP1. These ORFs were assigned based mainly on previously annotated genomes of closely related phages. We also checked ORFs with greater than 100 nucleotides that initiated with CTG, TTG, and GTG and identified eight putative genes using this criterion (Additional file 4: Table S4). Blast analysis showed that two of the hypothetical genes having alternative start codons (ORFs 48 and 56) had 100% similarities to putative proteins of closely related pbunaviruses (data not shown). The annotation of the genome of phage BrSP1 showed the typical block organization, with groups of genes in the same orientation (Additional file 4: Table S4). Other members of the genus Pbunavirus display similar genomic organization [37,41]. We analyzed codon usage based on the relative synonymous codon usage (RSCU), as described by Sharp and Li [45]. RSCU is the ratio of the actual number of times a codon occurs by the number of times that codon was expected to be used supposing an equal distribution of all codons for that particular amino acid [45]. RSCU values higher than one indicate a codon is used more frequently than expected (positive bias), whereas values smaller than one indicate a codon that is used less frequently than expected (negative bias). Our analysis showed that phage BrSP1 codon preferences were, in general, similar to those of the host (Additional file 5: Table S5). Overt differences were observed in few cases as, for instance, with the TAA stop codon and the ACT codon for threonine. Both of these codons display a strong negative bias in P. aeruginosa [46], but showed positive bias (TAA) or no bias (ACT) in phage BrSP1 (Additional file 5: Table S5).
We performed a whole genome comparison among all genomes selected for phylogenetic analysis using the MAUVE method (Fig. 4). In this analysis, the genomes were oriented to start at the terminase gene to overcome the uncertainties in the genome start of several isolates. The genomic architecture of all viruses showed a conserved collinear region flanked by less conserved areas towards the extremities (Fig. 4). All viruses analyzed here showed limited genomic diversity, but some long indels (> 1 kb) were observed, such as the indel resulting in the loss or gain of PB1-ORF070 (depending on the genome), which codes for a hypothetical protein related to chaperone DnaJ (Fig. 4, highlighted in black). Overall, the genome of Pseudomonas phage BrSP1 was similar in gene content and length to all members of the genus Pbunavirus, suggesting that horizontal gene transfer (HGT) events from other phages and host genomes are not frequent in the members of this genus. This observation is compatible with the low gene flux (LGCF) evolutionary mode proposed by Mavrich and Hatfull [47] for the evolution of lytic phages.
Classification and phylogeny of phage BrSP1 based on the analysis of the whole genome and of the terminase gene The International Committee on Taxonomy of Viruses (ICTV) has chosen 95% whole genome DNA sequence identity as the criterion for demarcation of phage species based on the BLASTN algorithm [48]. Whole genome analysis using the Basic Local Alignment Search Tool (BLAST) [49] showed that the phage BrSP1 BLAST best hit was Pseudomonas phage NH-4 (Genbank accession number: JN254800) with the highest percent identity of 97% (coverage of 96%, data not show). In order to classify phage BrSP1 accurately, pairwise genome alignments were performed between all members of Pbunavirus Results were evaluated as "Not clear" when no inhibition zone was observed; "Clear" when complete clearing was observed; "Turbid" when the inhibition zone was mostly turbid 3 Spot test done with diluted phage stock (up to 10 6 dilution). The formation of individual lysis plaques is indicated 4 Strain Lfar01 was used as a control genus using PASC [50]. As shown in Table 3, the pairwise analysis yielded lower identity values than those obtained with BLAST, as the percentage identity in the later was based on partial coverage. Pairwise genomic alignment revealed that BrSP1 phage whole genome identity varied from 86% (compared to Pseudomonas phage PaMx13, Table 3) to 94% (compared to Pseudomonas phage JG024, Table 3). As the genome of BrSP1 differed from all the other species by more than the 5% limit proposed by the ICTV phage BrSP1 is a putative novel species of the genus. The terminase gene, which encodes an enzyme involved in the packaging of phage DNA into capsids, is a phylogenetic marker used in the investigation of several phage groups [51]. We inferred the phylogenetic relationship of phage BrSP1 to other pbunaviruses using both the terminase gene ( Fig. 5a) and whole genome alignments (Fig. 5b). As shown in Fig.  5a, the phage BrSP1 terminase gene clustered together with a group of three ICTV recognized species (Pseudomonas virus LBL3, Pseudomonas virus SN, and Pseudomonas virus DL68), with a pairwise identity of 98, 99, and 99, respectively (data not shown). The phylogenetic tree inferred using whole-genome alignment ( Fig. 5b) was partially incongruent with the terminase-based phylogeny (Fig. 5a)

Discussion
The work here presented comprises the isolation and the characterization of the Pseudomonas aeruginosa BrSP1 bacteriophage. Our study involved the sequencing of the phage genome and the investigation of biological traits, such as host-range, "in vitro" biocontrol capacity, EOP estimation and the appearance of BIMs. These latter studies were done as a preliminary investigation to evaluate the potential of using BrSP1 in phage therapy against P. aeruginosa strains that infect domestic animals.
The sequencing of the BrSP1 genome revealed that its G + C content was 55.7%, which is very similar to those of other pbunaviruses [41] and lower than the GC content of P. aeruginosa, which is 66.6% [52]. Several other lytic phages that infect P. aeruginosa also display lower levels of GC content compared to the host [53][54][55]. Likewise, lytic phages that infect other members of the genus Pseudomonas, such as P. fluorescens [56][57][58] P. putida [59] and P. stutzeri [60], also have lower GC content levels compared to their host. Kwan and coworkers [61] hypothesized that differences in GC content between strains of P. aeruginosa and their phages could be due either to lateral gene transfer, to the recent invasion of these viruses to P. aeruginosa or some other characteristic feature of the viruses, such as variations on DNA polymerase and other replication related genes. Recent studies have shown that there is a relationship between the occurrence of mutator genes and genome GC variation [62,63], suggesting that enzymes involved in The GC content of bacteria is considered a result of a complex interaction of endogenous and exogenous factors [64]. Analysis of genome databases has shown that the length of bacterial genomes is directly related to GC content [64]. Moreover, the presence of certain enzymes involved in replication and repair processes correlates with GC content [63,65]. In their recent analysis, Almpanis and coworkers [64] concluded that phages, in general, have lower GC content than their hosts. The evolutionary advantage for lower GC content in these obligatory intracellular parasites is not well understood [64]. The higher energetic costs involved in the synthesis of GTP and CTP may be an influence [66]. Also, the fact that ATP pools are more abundant than that of other nucleotides may also affect this feature [66].
Our host range investigation, done mostly with P. aeruginosa strains isolated from domestic animals, showed that phage BrSP1 infected 51.4% of the strains tested (Table 1 and Additional file 2: Table S2). The susceptible strains included bacteria isolated from cows, dogs, and horses from different regions of Southern Brazil (Table 1 and Additional file 2: Table S2). The ability to infect a broad diversity of P. aeruginosa strains seems to be a trait of the pbunaviruses. Garbe and coworkers [42], for instance, found that phage JG024 infected 50% of 100 environmental strains collected from rivers and 84% of 19 clinical strains. Likewise, an investigation with phage  KPP12 showed that this virus was able to infect 52.6% of 38 clinical strains [67]. Our EOP analysis revealed up to nearly one hundred fold EOP differences among the susceptible strains (Fig.  2). Changes on surface molecules that serve as bacteriophage receptors can significantly affect EOP and reduce phage adsorption [68][69][70]. It is, therefore, likely that the substantial EOP differences we found are a result of differences on the receptors of the strains. The majority of the susceptible strains were of medium or high efficiency according to the parameters proposed by Viazis and coworkers [71]. That is, their EOP values were above 0,2 [71]. Only three of the susceptible strains (CASM03, BOIJ03, and CASM11) displayed EOP values below 0,2 and were, thus, considered "low efficiency" (Fig. 2). The phage stock solution used in the EOP assay for these three strains had to be increased tenfold in order to produce lysis plaques. This need suggests that the use of phage BrSP1 as a control agent for the low-efficiency strains would require higher doses in order to be effective.
Our "in vitro" biological assays showed that bacterial growth remained checked for the first 12 h after phage infection, but that after 24 h there were little differences between infected and non-infected samples (Fig. 3). To Fig. 4 Whole genome synteny analysis of pbunaviruses genomes. Schematic genome alignment diagram obtained using the Mauve software package and plotted using genoPlotR. The analysis was performed using 20 complete genomes available in Genbank, as described in Table 3, and with the phage BrSP1 genome. To facilitate synteny analysis, the genomes were oriented to start at the terminase gene (highlighted in black) to overcome the uncertainties in the genome start of several isolates. Gray shading indicates homologous regions among genomes, colored bars represent nucleotide syntenic blocks, and the white spaces between blocks indicate missing regions. All viruses were collinear. The PB1-ORF070 is annotated and highlighted in black find out the origin of this rapid development of resistance, we plated P. aeruginosa on a double layer agar plate which had phage BrSP1 on the overlay. The growth of colonies under these conditions showed that resistant variants arose independently of the presence of phage BrSP1 (Table 2). This result indicated that the growth observed at 24 h post-infection was attributed, at least partially, to the multiplication of resistant variants that occurred independently of the presence of the phage. Thus, these resistant bacteria were of the type Luria and Delbruck [72] termed "original variants." Part of the resistant colonies that grew on the phage overlay produced a red pigment, similar to what Le and colleagues [69] observed in their study of phage-resistant P. aeruginosa. Their results also indicated that resistance arose spontaneously [69]. Other recent studies, both with P. aeruginosa and with other species, have investigated the spontaneous emergence of BIMs [68,73,74].
The development of resistance is a critical concern for phage therapy [75]. Thus, it is essential to consider the rapid recovery observed in our "in vitro" assays if phage BrSP1 is to be used to control P. aeruginosa. However, it is important to point out that the rapid development of resistance "in vitro" is not necessarily a barrier to the application of bacteriophages. Two recently published studies with P. aeruginosa support this assertion [76,77]. Oeshslin and colleagues [76], for instance, showed that a combination of phages and antibiotics cleared heart valve infections in a rat model. Importantly, they investigated the growth of phage resistant bacteria in fibrin clots and heart tissues of infected rats. The rapid growth of resistant bacteria in the "in vitro" fibrin-clot model contrasted with the absence of BIMs in the tissues of the experimental animals [76]. Similarly, effective infection treatment in a cancer patient was observed using antibiotics and a phage cocktail which included phages that promoted rapid growth of resistant P. aeruginosa in broth cultures [77]. In both of the above cases, the "in vitro" assays displayed patterns similar to that observed in our investigation regarding the growth of phage-resistant variants (Fig. 3).
In fact, "in vitro" control patterns similar to that of BrSP1 has been observed in other studies done with Gram-negative bacteria [78,79] (Jumpei Uchiyama, Ken Fukuda and Shigenobu Matsuzaki, personal communication, January, 2019) and the rapid regrowth of the bacterial populations seems to be a general phenomenon. The study with phages that infect Vibrio anguillarum, a pathogen of fish larvae, also showed that the fast recovery of the bacterial population "in vitro" is not necessarily an impediment to the successful application of a phage cocktail as a control measure [79].  Evidence from studies with other pbunaviruses indicates that resistance development is not a barrier for their application in phage therapy. Fukuda and coworkers [67], for instance, successfully treated keratitis in mice with a single application of Pseudomonas phage KPP12, a pbunavirus closely related to phage BrSP1 ( Table 3). The treatment reduced bacterial load to almost undetectable levels and the cornea of treated mice had a remarkable improvement [67]. In another study, the use of a phage cocktail containing a pbunavirus and a podovirus practically eliminated P. aeruginosa from the lungs of infected mice [37].
The use of phage cocktails and the concomitant use of phages and antibiotics are strategies that can be used to reduce the growth of BIMs [76][77][78]. As changes in bacteriophage receptor is a common cause of resistance, the phages in the cocktail should ideally use different receptor molecules [68]. Crucially, mutations that render the bacteria resistant to bacteriophages often reduce the fitness of the new variants [68,69]. This fitness cost may be exploited to promote the successful application of phage therapy [75].

Conclusion
In this work, we have isolated and characterized a P. aeruginosa bacteriophage, which we named Pseudomonas phage BrSP1. The virus was able to infect 51.4% of our collection of P. aeruginosa strains. The susceptible strains were clinical isolates from domestic animals and the reference strain ATCC 27853. "In vitro" bioassays using phage BrSP1 resulted in a rapid reduction of bacterial load. These results indicated that phage BrSP1 is a broad host range virus that could be included in a phage cocktail to treat P. aeruginosa infections in domestic animals. The efficiency of plating (EOP) values for phage BrSP1 varied considerably among the tested strains. We suggest the use of publicly available strains as the reference in EOP analysis, as it allows comparison of results obtained in different laboratories. Based on comparisons of complete genome sequences, we propose that Pseudomonas phage BrSP1 is a novel species of the genus Pbunavirus. In our phylogenetic analysis, we found topological incongruences between the commonly used terminase gene and the tree constructed using the complete genome. The complete genome of BrSP1 phage provides additional data that may help the broader understanding of pbunaviruses genome evolution.

P. aeruginosa origin and maintenance
Except for P. aeruginosa Lfar01 and P. aeruginosa ATCC 27853, the strains used in this study were isolated from samples sent from various sources for microbiological analysis at Dr. Agueda Vargas laboratory (Department of Preventive Veterinary Medicine, Federal University of Santa Maria, Santa Maria, Brazil; Additional file 2: Table S2). The bacteria isolated from animals are registered in the SISGEN platform under the registration code A28EE0E. SISGEN (Sistema Nacional de Gestão do Patrimônio Genético e do Conhecimento Tradicional Associado) is a national registration system of Brazilian biological material. Samples were cultivated on blood agar containing 5% sheep blood and on MacConkey agar (Difco). Cultures were incubated aerobically for 24-48 h at 36°C. After incubation, isolates were identified according to conventional microbiological methods [80]. The origin of strain Lfar01 is unknown. P. aeruginosa was propagated in Tryptic Soy Broth (TSB) and Agar (TSA) (Fluka Analytical). TSA supplemented with 15% glycerol was used for long-term storage at − 80°C.

16S rRNA gene sequencing
Sequencing analysis of the 16S rRNA genes were done for the strain Lfar01 and for those isolates that were not susceptible to BrSP1 infection. DNA extraction was done with the kit Brazol (LGC Biotecnologia), according to the manufacturer's instructions. Primers 27f (5′-CAGGCC TAA CAC ATG CAA GTC-3′) and 1387r (5′-GGG CGG WGT GTA CAA GGC-3′) [81] were used for the amplification and the sequencing reaction of the purified amplicon. The reactions were carried out in a T100™ Thermal Cycler (Bio-Rad®) following the conditions described by Marchesi, et al. [81]. The PCR product was purified using E.Z.N.A.® Cycle-Pure Kit 200 (Omega Bio-Tek®) and sequenced with the BigDye Terminator Kit® (v. 3.1 RR-100) according to the manufacture's instruction (PE Applied Biosystems). Electrophoresis was carried out on ABI 3130 Genetic Analyzer (Applied biosystems) at Genomic Engenharia Molecular Ltda (São Paulo, Brazil). The sequences obtained were assembled in the CodonCode Aligner, version 3 (CodonCode Corporation) and compared with those of the Genbank using the Blast program [49].

Phage isolation
Phage BrSP1 was isolated from samples of untreated sewage collected at SABESP Pinheiros (Water treatment utility located in Sao Paulo, Brazil). Six samples of 50 mL were transferred to conical centrifuge tubes and subjected to centrifugation (centrifuge Hitachi®, model CR22G) at 4000 x g for 5 min. A total of 45 mL of supernatant was transferred to borosilicate flasks containing 45 mL of doubled concentrated TSB media. Approximately 5 mL of P. aeruginosa strain Lfar01 (Additional file 2: Table S2) at exponential growth phase was added to the flask, and incubated at 36°C in a benchtop incubator shaker (model Classic C24, New Brunswick Scien-tific®) at 120 rpm for 24 h. After this period, a 10 mL sample was removed from the flask, transferred to a 15 mL conical tube and centrifuged for 10 min at 4000 x g (centrifuge Hitachi®, model CR22G). Approximately 5 mL of the supernatant were removed and passed through a sterile polyvinylidene fluoride syringe filter (JET BIOFIL®) with pore size of 0.22 μm.
Phage presence was evaluated using the small drop assay. The double agar method [82] was used to prepare the bacterial lawn. TSA media was plated in the bottom layer and the top layer was prepared adding 0.25% of agarose (Agargen®) to 5 mL of TSB. The agarose was dissolved by heating the broth in a microwave oven and then kept at 50°C in a water bath. Approximately 1 mL of P. aeruginosa at exponential growth phase was added to the TSB/agarose mixture immediately before pouring it over the bottom layer. Aliquots of 20 μL of each of the six samples, processed as described above, were spotted on the top layer. The agar plates were incubated for 24 h at 36°C. After this period, the plates were checked for the presence of clear lysis zones, indicating the presence of bacteriophages.
Phage purification, titration of phage stocks, storage stability and electron microscopy The plaque purification method described by Miller [83] was used to obtain a homogeneous phage stock. The viral sample that produced plaques was serially diluted using SM buffer (100 mM NaCl; 8 mM MgSO4•7H2O; 50 mM Tris-Cl pH 7.5; 0.002% gelatin) up to 10 million fold. A total of 20 μL of the serially diluted phage was added to 200 μL of P. aeruginosa cells at stationary phase of growth. The samples were homogenized, kept at room temperature for 20 min., mixed to 6 mL of molten TSB supplement with salts and 0.25% agarose, poured over a TSA plate, and incubated for a period of 48 h. at 36°C.
After the incubation period, a slab of the double layer with a single plaque was removed from the plate with a blade and transferred to a microtube containing 1 mL of SM buffer. The microtube was subjected to agitation in a benchtop shaker (Quimis®) for 1 h at 500 rpm. Phage released in the buffer was then purified with 200 μL chloroform. After 10 min of agitation, the sample was centrifuged for 10 min at 5000 x g (Centrifuge 5418, Eppendorf®). The supernatant was transferred to a new microtube. The whole process, starting from the serial dilution, was repeated twice. After the third round of purification, 10 μL of the chloroform-purified sample was used to infect 10 mL of P. aeruginosa cells at early exponential growth phase. The infected cells were incubated at 36°C in a benchtop incubator shaker (model Classic C24, New Brunswick Scientific®) at 120 rpm for 24 h. The culture was then centrifuged and the supernatant passed through a sterile Polyvinylidene fluoride syringe filter (JET BIO-FIL®) with pore size of 0.22 μm. This was the phage stock that was used in the subsequent analysis. The titer of this stock was evaluated. For that, three series of serial dilutions were done with a sample of the stock. The serially diluted phages were used to infect cells as described above. Plaques formed at specific dilutions were counted and used to estimate the plaque forming units per mL.
We investigated the storage stability of phage BrSP1 under refrigeration (approximately 5°C). We estimated the titer of a BrSP1 stock in triplicates before and after storage for 104 days in the refrigerator (Additional file 6).
We also examined Phage BrSP1 by electron microscopy. Two hundred microliters of a high titer chloroform purified phage stock were placed on a Parafilm® surface. Copper grids (200 mesh) covered with carboncoated Formvar film were floated on them for 10 min, washed with droplets of distilled water and floated on droplets of 1% aqueous uranyl acetate for 10 min. After that, grids were removed and excess of liquid eliminated with a filter paper. These negatively stained preparations were examined in a JEOL JEM 1011 transmission electron microscope (installed at the Laboratório de Microscopia Eletrônica, Departamento de Fitopatologia e Entomologia, Escola Superior de Agricultura Luiz de Queiroz, Universidade de São Paulo, campus Piracicaba) at 60 KV and images registered digitally.

Host range investigation
In total, 36 P. aeruginosa strains isolated from domestic animals and P. aeruginosa strain ATCC 27853 were investigated regarding their susceptibility to bacteriophage BrSP1. Spot test on a double agar overlay [82] was used in this analysis. Plates contained each of the 37 bacterial isolates on the top layer were prepared adding 1 mL of a culture at early exponential phase to 5 mL of molten TSB media and 0.25% agarose (Agargen). A first round of analysis was done with the undiluted phage stock. Aliquots of 20 μL of phage stock containing approximately 4.5 × 10 5 PFU were spotted with a precision pipette on the surface of the agar overlay. The plates were incubated for a period of up to 48 h at 36°C. After this period, the plates were checked for the appearance of lytic inhibition zones, which were classified as clear, turbid or not present.
Bacterial isolates that displayed clear or turbid zones of inhibition on this test were used in a second round of analysis following the same procedure described above, except that serial diluted phage stock was used (up to 10 6 dilutions). Again, the plates were incubated for a period of up to 48 h at 36°C. This time, the plates were checked for the appearance of individual inhibition plaques.

Efficiency of plating (EOP) analysis
For the efficiency of plating (EOP) evaluation, the double agar plate assays were done according to the procedure described by Mirzaei and Nilsson [84]. The phage stock used in these experiments was titrated in triplicates using P. aeruginosa strain ATCC 27853. A volume of 10 μL of diluted phage stock, containing approximately 24 PFU, was added to 200 μL of susceptible P. aeruginosa culture at the stationary phase of growth. The samples were homogenized, kept at room temperature for 20 min., mixed to 6 mL of molten TSB and 0.25% agarose, poured over a TSA plate, and incubated for a period of up to 48 h at 36°C. The number of lysis plaques formed was recorded after this period. A few samples did not display lysis plaques with the standard phage inoculum used in the assay. For those samples, the experiment was redone using a phage stock that was 10 times larger. EOP was estimated as the ratio of the number of lysis plaques produced in the P. aeruginosa strain ATCC 27853 to the number of plaques produced in each of the other strains analyzed (Additional file 7). These assays were done in triplicates.

Biological assays
In order to evaluate the capacity of phage BrSP1 to control P. aeruginosa, we carried out "in vitro" experiments with strains Lfar01, ATCC 27853 and BOIJ02. TSB medium supplemented with CaCl2 (40 mg/L) and MgSO4 (40 mg/L) was used in these analyses. The phage concentration we used were either 10 6 or 10 7 plaque forming units per mL of bacterial culture media (PFU/mL). Test tubes with 10 mL of P. aeruginosa cells at early exponential phase were inoculated with phage BrSP1 and cultured at 36°C in a benchtop incubator shaker (model Classic C24, New Brunswick Scientific®) at 120 rpm for 24 h. We removed samples from the culture at the time of viral inoculation (t = 0) and at 2, 4, 6, 8, 10, 12, and 24 h post-infection (h p.i.) (t = 2, t = 4, t = 6, t = 8, t = 10, t = 12, and t = 24, respectively). Control tubes were inoculated with SM buffer. We evaluated cell concentrations at each of these times measuring the absorbance at 600 nm (OD 600 ) (NanoDrop model 2000c, Thermo Scien-tific®) (Additional file 8). Each set of assays was done in triplicate.

Frequency of occurrence of phage resistant mutants
We investigated the frequency of phage-resistant mutants in P. aeruginosa strains Lfar01, ATCC 27853 and BOIJ02 following the method employed by Filippov and colleagues [68]. We inoculated the strains on TSA plates using the steak plate technique to isolate individual colonies. We inoculated bacterial cells from 5 colonies of each strain to a test tube containing 5 mL of TSA. After overnight growth at 36°C, samples of the bacterial cultures were serially diluted. These samples were inoculated on TSA plates with or without a double layer containing approximately 1 × 10 9 PFU of phage BrSP1. One hundred microliters of undiluted bacterial culture and of dilutions up to 10 − 2 were spread plated on the surfaces of TSA plates containing phages on the top layer. The same volumes of dilutions 10 − 7 to 10 − 9 were spread plated on TSA plates without phages. Plates were incubated for a period of up to 48 h at 36°C. Individual colonies were counted and the concentration of phage resistant bacteria and total bacteria was estimated in terms of colony forming units per mL (CFU/mL). The frequency of spontaneous resistant mutations was estimated as the ratio between the numbers of resistant to total bacterial counts per mL of the overnight culture (Additional file 9).

Phage DNA extraction
Phage DNA extraction was done following the procedure described by Su et al. [85]. A total of 10 mL of P. aeruginosa cells (strain Lfar01) at exponential growth phase were inoculated with 50 μL of diluted phage stock. The inoculated cells were cultured at 36°C in a benchtop incubator shaker (model Classic C24, New Brunswick Scientific®) at 120 rpm for 24 h. The culture was then transferred to a 15 mL conical tube and centrifuged (centrifuge Hitachi®, model CR22G) at 4000 x g for a period of 5 min. The supernatant was transferred to a new tube and treated with 10 units of amplification grade DNase I (Sigma-Aldrich) for a period of 30 min. at 36°C. Following this treatment, 200 μL of a 2 M ZnCl 2 (filter sterilized) solution, which promotes phage particles precipitation, was added. The mixture was allowed to stand for 5 min at 36°C and was centrifuged (centrifuge Hitachi®, model CR22G) at 4000 x g for a period of 5 min. The phage pellet was dissolved in TENS buffer (50 mM Tris-HCl, pH 8, 100 mM EDTA, 100 mM NaCl and 0.3% SDS) and proteinase K was added to a final concentration of 100 μg/mL. The solution was vigorously mixed and incubated at 65°C overnight. The DNA was then extracted using a phenol/chloroform/isoamyl alcohol mixture and ethanol precipitated. The DNA was resuspended in 100 μL of ultrapure water.

Sequencing, assembly and analysis of the phage BrSP1 genome
The sequencing library was prepared using TruSeq Nano DNA Kit and sequenced with HiSeq 2500 System (2 × 100 paired-end) at Macrogen Inc. (Seoul, Republic of Korea). The resulting reads were filtered and trimmed to remove short and low-quality regions/reads and then assembled using CLC Genomics Workbench (CLC Bio, Aarhus, Denmark) with default parameters. The read mapping and genome annotation was performed using Geneious 7.1.8 [86]. ORFs were annotated using BLASTx search against the NCBI non-redundant protein database. Maximum likelihood (ML) tree reconstruction was performed with FastTree using the default options [87] on alignments of terminase genes and complete Funding This work was partially supported by CNPq (Conselho Nacional de Desenvolvimento Científico e Tecnológico) grant number 407908/2013-7 to BMR and FAP-DF (Fundação de Apoio à Pesquisa do Distrito Federal), grant number 193.001.532/2016 to BMR. The funding agencies had no role in the design of the study, in the collection, analysis, and interpretation of data and in writing the manuscript.

Availability of data and materials
All data generated or analysed during this study are included in this published article and its supplementary information files.
Authors' contributions BMR and JLCW conceived the project and designed the research experiments. APCV and VLE isolated and prepared the P. aeruginosa stocks. AMG and JLCW isolated and purified the bacteriophage. AMG carried out the "in vitro" bioassays. ACCM carried out the host range, the Efficiency of Plating experiments, the frequency of mutations experiments and the phage stability investigation. ACCM and JLCW sequenced, assembled and submitted the 16S rRNA genes from the P. aeruginosa isolates. DMPAA and FLM were responsible for the phage genome sequencing, assembly and annotation. DMPAA and FLM carried out the phylogenetic tree and syntheny analysis. EWK prepared the transmission electron micrograph of the bacteriophage. BMR, DMPAA, FLM, and JLCW wrote and revised the manuscript. BMR and JLCW provided research funds. All authors read and approved the final manuscript.
Ethics approval and consent to participate Not applicable.

Consent for publication
Not applicable.

Competing interests
The authors declare that they have no competing interests.

Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.