The Bacteriophage Pf-10—A Component of the Biopesticide “Multiphage” Used to Control Agricultural Crop Diseases Caused by Pseudomonas syringae

Phytopathogenic pseudomonads are widespread in the world and cause a wide range of plant diseases. In this work, we describe the Pseudomonas phage Pf-10, which is a part of the biopesticide “Multiphage” used for bacterial diseases of agricultural crops caused by Pseudomonas syringae. The Pf-10 chromosome is a dsDNA molecule with two direct terminal repeats (DTRs). The phage genomic DNA is 39,424 bp long with a GC-content of 56.5%. The Pf-10 phage uses a packaging mechanism based on T7-like short DTRs, and the length of each terminal repeat is 257 bp. Electron microscopic analysis has shown that phage Pf-10 has the podovirus morphotype. Phage Pf-10 is highly stable at pH values from 5 to 10 and temperatures from 4 to 60 °C and has a lytic activity against Pseudomonas strains. Phage Pf-10 is characterized by fast adsorption rate (80% of virions attach to the host cells in 10 min), but has a relatively small number of progeny (37 ± 8.5 phage particles per infected cell). According to the phylogenetic analysis, phage Pf-10 can be classified as a new phage species belonging to the genus Pifdecavirus, subfamily Studiervirinae, family Autographiviridae, order Caudovirales.


Introduction
Pseudomonas is a genus of Gram-negative bacteria including, inter alia, the known plant pathogen species. The phytopathogenic pseudomonads cause a wide range of plant diseases: necrotic lesions of fruit [1,2], stems [1,3], and leaves [1,2,4], hyperplasias (galls) [5], tissue macerations (rots) [6,7], cankers [1,8,9], blights [10], and vascular infections (wilts) [11]. Plant diseases caused by Pseudomonas are widespread throughout the world and infect a large number of higher plants, including agriculturally important crops, resulting in the loss of agribusiness profits and other consequential damages that ultimately have a negative impact on the world food supplies and economy. Moreover, the phytopathogeninfested plants can be a source of toxic molecules [12,13].
Pseudomonas syringae is one of the most widespread and best studied Pseudomonas plant pathogens. Depending on the range of hosts and the type of symptoms they cause, these pathogens are assigned to different intraspecific taxa called pathovars [14]. Within Pseudomonas syringae, there are at least 60 pathovars [15]. In particular, Pseudomonas syringae pv. tomato is the causal agent of tomato speck disease [16]; P. syringae pv. syringae cause blossom blast of apple [17] and pear [2], apical necrosis of mango [18], canker of cherry [19]

Phage Isolation and Propagation
Phage Pf-10 was isolated from a tissue sample of green bean (Phaseolus vulgaris) infected with P. syringae pv. syringae in 2014 [40]. Briefly, the plant material was collected in the Vitebsk Oblast, Belarus. The sensitive strain P. fluorescens BIM B-582 was used as the host strain for initial propagation. P. fluorescens BIM B-582 culture was grown in GRM broth [8.0 g/L pancreatic hydrolysate of fish-meal, 8.0 g/L enzymatic peptone, 4.0 g/L NaCl] (Obolensk, Russia) at 28 • C with shaking (250 rpm) overnight. The culture was then diluted 10-fold in a fresh medium and incubated at 28 • C with shaking (250 rpm) to OD 590 of 0.8. The plant sample was pounded in a mortar with 30 mL of GRM broth, followed by centrifugation at 5000× g for 30 min. Then, 25 mL of the supernatant was mixed with 25 mL of a BIM B-582 culture and incubated for 18 h at 28 • C with shaking (250 rpm). The cell debris was removed by centrifugation at 6000× g for 30 min, and the supernatant was analyzed for the presence of phages. Phage assay was carried out by using the double agar layer method. After overnight incubation at 28 • C, the plates were examined for plaques. A separate plaque was selected and used for further phage propagation, and the resultant phage preparation was used to obtain separate plaques. The propagation-plaque selection cycle was repeated five times in order to exclude the presence of other phages. Phage Pf-10 was deposited into the Belarusian collection of non-pathogenic microorganisms (BIM) Viruses 2022, 14,42 3 of 20 of the Institute of Microbiology, National Academy of Sciences, Minsk, Belarus under the collection number BIM BV-61-D [39].
The sensitive strain P. syringae BIM B-268 [41] was used for phage propagation to obtain high-titer suspension. Briefly, 500 µL of the overnight bacterial culture was transferred into 50 mL of LB with 10 mM CaCl 2 and 10 mM MgCl 2 and incubated at 28 • C with shaking at 120 rpm to the OD 590 of 0.2 (at approximately 1 × 10 8 colony-forming units (CFU)/mL). Then, the 25 µL phage sample (2 × 10 8 plaque-forming units (PFU)/mL) was mixed with 50 mL of the BIM B-268 culture (OD 590 of 0.2) with a multiplicity of infection (MOI) of 0.001 and incubated at 28 • C with shaking at 120 rpm for 6 h until the complete lysis of the bacterial culture. After the lysis, PEG 8000 (polyethylene glycol 8000) precipitation was performed as described previously [42]. The PEG-precipitated phage sample was dissolved in SM buffer [50 mM Tris-HCl, pH 7.5; 100 mM NaCl; 1 mM MgSO 4 ; 0.01% gelatin]. The resultant phage preparation was filtered through a 0.22 µm sterile filter and stored at 4 • C. Phage titer was determined by the double agar overlay plaque assay method. The Pf-10 phage titer in the final high-titer preparation was 5 × 10 11 PFU/mL.

Transmission Electron Microscopy
Three mL of high-titer phage preparation were centrifuged using the preformed CsCl gradient (1.25 g/mL, 1.4 g/mL, 1.5 g/mL and 1.7 g/mL, 2.5 mL each) at 10 • C, 25,000 rpm, for 2.5 h in a Beckman Coulter ultracentrifuge, L7-55, using a SW 41 Ti rotor. Further, ten µL of the concentrated phage suspension (10 9 PFU/mL) were applied to carbon-coated copper grids (400 mesh) and negatively stained with 1% uranyl acetate. The grids were analyzed using a JEM 1200EX (JEOL, Tokyo Japan) transmission electron microscope at 80 kV accelerating voltage. Images were taken on Kodak film SO-163 (Kodak, Cat. \# 74144, Hatfield, PA, USA). Phage particle dimensions were measured using ImageJ version 1.53e in relation to the scale bar generated by the microscope.

Host Range
Phage host range analysis was performed by spot test using 45 strains of the genus Pseudomonas (Supplementary Information, Table S1). Five µL of the phage suspension (10 9 PFU/mL) were dripped onto agar plates with the different strains, followed by incubation at 28 • C for 20 h.

Thermal and pH Stability Tests
Phage stability tests were conducted as described previously [43], with some modifications. Briefly, the phage stability at various temperatures (4,30,40,50,60,70,80 and 90 • C) was investigated by incubating the phage (at approximately 5 × 10 9 PFU/mL) at the respective temperatures for 1 h. Phage stability at various pH values (ranging from 2.2 to 10) was evaluated using four buffers: Glycine-HCl buffer (pH values 2.2 and 3), sodium acetate buffer (pH 4 and 5), phosphate buffer (pH 6, 7 and 8) and Glycine-NaOH buffer (pH 9 and 10). The phage sample was mixed with each buffer to a final phage concentration of 5 × 10 9 PFU/mL and incubated at 28 • C for 1 h. After incubation, the phage titer was determined by double-layer plate titration. Three independent trials of the experiment were performed. The results were processed using GraphPad Prism 8.4.3 as the mean of three observations ± standard deviation.
2.6. The Effects of Ca 2+ and Mg 2+ on the Killing Activity of Pf-10 In order to assess the killing activity of Pf-10 and the calcium and magnesium effects on this process, the BIM B-268 culture was infected with the phage at MOI values of 10, 1 and 0.1 in the presence or absence of Ca 2+ and/or Mg 2+ ions in solution. Briefly, 50 µL of phage preparations (at approximately 8 × 10 9 , 8 × 10 8 , 8 × 10 7 PFU/mL) were mixed separately with 450 µL of the log-phase BIM B-268 culture (at approximately 8 × 10 7 CFU/mL) in a 48-well microplate. A non-infected BIM B-268 culture was used as a control to assess the killing activity. In order to estimate the effects of Ca 2+ and/or Mg 2+ ions, CaCl 2 and/or MgCl 2 were added to the phage-host mixtures (to the final concentration of 10 mM each), SM buffer [50 mM Tris-HCl, pH 7.5; 100 mM NaCl; 1 mM MgSO 4 ; 0.01% gelatin] was added as a control. Then the microplate was incubated at 28 • C for 3 h in a FilterMax F5 microplate reader (Molecular Devices, San Jose, CA, USA) with the following shaking settings: shake mode-orbital, shake intensity-medium, with OD 595 being measured every 10 min. Three independent trials of the experiment were performed. The results were reported as the mean of three observations ± standard deviation. The resultant growth curves were visualized using GraphPad Prism 8.4.3.
The Pf-10 lytic activity in the presence or absence of Ca 2+ and/or Mg 2+ ions was evaluated using the PhageScore method developed by Konopacki [44]. Briefly, a fitting of the data on the growth kinetic of P. syringae BIM B-268 upon Pf-10 infection at different MOI was performed by non-linear least-squares method using the LMFIT library (https: //lmfit.github.io/lmfit-py/intro.html, accessed 29 September 2021). To fit the growth curves of the control (non-infected BIM B-268 culture) and test (the BIM B-268 culture infected by the Pf-10 phage) samples, formulas 1 and 2 described by Konopacki [44] were used, respectively. To calculate the PhageScore values (P S ) of the Pf-10 phage at different MOI (P Si , PhageScore for given MOI i , i-the number of different MOI values for the same bacteriophage), we used the area ratios obtained by the mathematical fitting of functional dependencies, formula 8 in Konopacki's article [44]. The areas under the curve were calculated using Simpson's rule. The value of the Phage Total Score factor (P TS ) was calculated using formula 12 from the article [44].

Adsorption Assay
In order to determine the attachment time of the phage to a bacterial cell, an adsorption assay was performed according to the protocol developed by Kropinski [45]. Briefly, Eppendorf tubes containing 0.95 mL of LB broth with three drops of chloroform were placed on ice to chill for 10 min. Nine mL of the log-phase BIM B-268 culture (OD 590 of 0.2, approximately 1 × 10 8 CFU/mL) were transferred into a 100-mL laboratory flask and incubated in a shaking water bath at 28 • C, 60 rpm for 5 min. Nine mL of LB broth were used as a control. Further, 1 mL of the phage suspension, preheated for 5 min at 28 • C, with a concentration of 1 × 10 7 PFU/mL was added into the test flask to provide MOI values of 0.01. The phage suspension was added into the control flask in the similar way. Then 50-µL aliquots were taken from both test and control flasks every minute, transferred into the prepared Eppendorf tubes and vortexed vigorously. The obtained mixtures were serially diluted in SM buffer and quantified using double agar overlay plaque assay. Five independent trials of the experiment were performed. The results were presented as percentages of the initial phage number and processed using GraphPad Prism 8.4.3 with error bars representing standard deviation for five trials. The adsorption rate was calculated using the equation described by Kropinski [45].

One-Step Growth Curve
In order to determine the average burst size of the phage, a one-step growth experiment was carried out as described by Hyman and Abedon [46]. One mL of phage suspension (1 × 10 7 PFU/mL), preheated for 5 min at 28 • C, was added to 9 mL of the bacterial culture (1 × 10 8 CFU/mL) to provide MOI values of 0.01. The resultant mixture was incubated in a shaking water bath at 28 • C, 60 rpm for 6 min for phage adsorption. After that, a 1-mL aliquot was transferred into a 1.5 mL Eppendorf tube and centrifuged at 3500× g for 10 min at 4 • C to precipitate the cells. The pellet was resuspended in 1 mL of fresh LB broth followed by addition to 9 mL of LB into a 100-mL laboratory flask and incubation in a shaking water bath at 28 • C, 60 rpm. The 50-µL aliquots of the mixture were collected at 10-min intervals for 1 h and diluted in SM buffer followed by phage quantification using double agar overlay plaque assay. Three independent trials of the experiment were performed. The PFU/mL values were calculated and plotted against time. The results were processed using GraphPad Prism 8.4.3 with error bars representing standard deviation for three trials. The latent period was determined as the interval between the adsorption of the phages to the bacterial cells and the release of phage progeny. The burst size of the phage Pf-10 was determined as the ratio of the average number of free phage particles after the release phase (plateau average [PFU/mL]) to their number during the latency phase (latent average [PFU/mL]).

Genome Sequencing, Assembly and Sequence Analysis
Phage preparation was treated with DNase I and RNase A. Phage DNA was extracted using DNeasy Blood & Tissue Kit (Qiagen, Hilden, Germany) according to the manufacturer's instructions. The library for genome sequencing was constructed using a Nextera XT DNA Library Prep Kit (Illumina Inc., San Diego, CA, USA) and sequenced on the Illumina MiSeq platform using a MiSeq reagent kit v.3 (2 × 300 bp). The genome assembly was performed by SPAdes v. 3.1 [47], and resulted in a single contig with the average coverage of 788×. Open reading frames (ORFs) were identified using RASTtk v.2.0 [48] followed by functional annotation using BLASTp v2.13.0+ (NCBI) [49] and HHpred (accessed 12 May 2021) [50]. ARAGORN v1.2.41 [51] was used to identify putative tRNAs and transfer-messenger RNAs. BRIG software v. 0.95 [52] was used for genome visualization. Phage protein were compared using the PSI-BLAST (Position-Specific Iterated BLAST) algorithm [49].

The Genome Packaging Strategy
Standard restriction analysis was used to confirm the correctness of genome assembly and identification of the bacteriophage genome termini of Pf-10 as described previously [43]. In order to determine the ends of the Pf-10 genome more accurately, rapid amplification of the ends of the genome (RAGE) was performed, as described previously [43], with minor changes to the protocol. Phage DNA was used for a typical DNA tailing reaction with the terminal transferase (New England Biolabs, Ipswich MA, USA, Cat. # M0315L) according to the protocol provided by the enzyme manufacturer. Further, two PCRs were carried out sequentially using the TaqSE DNA polymerase (SibEnzyme, Novosibirsk, Russia, Cat. # E314) and the pairs of oligonucleotides designed for the right and left ends of the Pf-10 genome (Supplementary Information, Table S3). The 3'-end tailing fragments were used as the DNA template for the first PCR, and the product of the first PCR was used as the template for the second reaction. The final PCR products were extracted from electrophoresis gel and used for Sanger sequencing with primers Pf-10_R_for2 5 -CCTTTAGGGTGCAGCACATC-3 and Pf-10_L_rev2 5 -TGATGGTCCTCTATGGGCCT-3 for the right and left ends of the genome, respectively. Due to the limitations of the RAGE method in locating the genome ends in the sequences containing "A" nucleotides, the right end of the Pf-10 genome was re-sequenced using the genomic DNA of phage Pf-10 as a template. The additional Sanger sequencing reaction for the right genome end was performed using a LiCor4300 DNA analyzer (LI-COR, Lincoln, NE, USA) with the primer Pf-10-t 5 -TGAGAATCATGTGCTATCTG-3 .
The genome packaging strategy of Pf-10 was predicted by phylogenetic analysis of the large terminase subunit as described by Merrill [53]. 57 large terminase proteins were aligned using MAFFT v7.308 [54] with default settings. The phylogenetic tree was computed using the neighbor-joining method by MEGA X [55] with a bootstrapping set to 100 and was visualized using FigTree v1.4.4 [56] with midpoint rooting.

Comparative Genomics
In order to find the related phages, ViPTree server version 1.9 [57] was used to generate a proteomic tree based on the genome-wide sequence similarities computed by tBLASTx. In addition to the reference virus genomes from the ViPTree server database, 10 phage genomes found by the BLASTn search using the whole Pf-10 genome sequence as the query were added to create the proteomic tree (Supplementary Information, Table S5). A linear comparison diagram showing the genomic identity between Pf-10 and the most closely related phages was also created with the ViPTree server version 1.9 [57]. The number of shared proteins was computed using the GET_HOMOLOGUES software v3.3.3 [58] with the COGtriangles algorithm [59] (−t 0 −C 75 −e).

Accession Number
The phage genome sequence was deposited into the GenBank database under the accession number KP025626 and NCBI Reference Sequence database under the accession number NC_027292.1 (BioProject accession number PRJNA263998, BioSample accession number SAMN21529907). The raw sequence dataset is available at the Sequence Read Archive under the accession number SRR15990561.

Statistical Analysis
All data comparisons performed in this study were verified using statistical analyses in order to ascertain the significance of the results. The statistical analyses were performed using GraphPad Prism 8.4.3. One-way ANOVA with repeated measures was used to assess the significance of differences in phage concentrations between control and tested samples during phage stability experiments (Section 3.2. pH and Thermal Stability Tests). Two-way ANOVA with repeated measures was used to analyze the statistical differences between the growth curves of the non-infected and infected P. syringae BIM B-268 cultures upon Pf-10 infections at different MOI (Section 3.3. Killing Assay and Effect of Ca 2+ and Mg 2+ on the Killing Activity of Pf-10). Two-way ANOVA with repeated measures was also used to analyze the statistical differences between the growth curves of the P. syringae BIM B-268 upon Pf-10 infections at different MOI in the presence and absence of Ca 2+ and Mg 2+ ions (Section 3.3. Killing Assay and Effect of Ca 2+ and Mg 2+ on the Killing Activity of Pf-10). A value of p ≤ 0.05 was considered statistically significant.

Phage Isolation, Host Range and Morphology
Phage Pf-10 was isolated in 2014 from a sample of green bean (Phaseolus vulgaris) tissue infected with P. syringae pv. syringae [40]. The Pf-10 phage was kindly provided by the Belarusian collection of non-pathogenic microorganisms (BIM) of the Institute of Microbiology, National Academy of Sciences, Minsk, Belarus. On the lawn of the sensitive strain P. syringae BIM B-268, the phage Pf-10 formed clear plaques approximately 4-8 mm in diameter with a characteristic small zone of incomplete lysis ( Figure 1). Host range analysis revealed that Pf-10 was able to lyse 13 (28.9%) out of 45 strains (Supplementary Information, Table S1). Seven of the thirteen strains susceptible to the Pf-10 phage (Supplementary Information, Table S1) are considered non-pathogenic, including the P. fluorescens BIM B-582 strain used for the industrial production of the "Multiphage" phages [39].  Table S1). Moreover, the latter three pathogenic strains infect cucumber fruits (Cucumis sativus L.). It should be noted that the "Multiphage" containing the Pf-10 phage is registered as a biopesticide against Pseudomonas-associated cucumber infections [39].
Morphological analysis via TEM showed that phage Pf-10 had an icosahedral capsid of approximately 55.5 ± 2.8 nm in diameter and a very short, noncontractile tail of approximately 11.2 ± 0.8 nm in length ( Figure 2). Electron microscopic analyses demonstrated that phage Pf-10 has the typical morphological features of the podovirus morphotype ( Figure 2).

pH and Thermal Stability Tests
The loss of phage activity during storage, transfer, or phage inactivation due to in vivo environment (e.g., pH) may result in poor outcomes in case of using phage-based preparations in phage therapy or bacteriophage-based biocontrol. This issue requires attention in the area of phage formulation to ensure phage viability and stability during storage and use.

pH and Thermal Stability Tests
The loss of phage activity during storage, transfer, or phage inactivation due to in vivo environment (e.g., pH) may result in poor outcomes in case of using phage-based preparations in phage therapy or bacteriophage-based biocontrol. This issue requires attention in the area of phage formulation to ensure phage viability and stability during storage and use.

pH and Thermal Stability Tests
The loss of phage activity during storage, transfer, or phage inactivation due to in vivo environment (e.g., pH) may result in poor outcomes in case of using phage-based preparations in phage therapy or bacteriophage-based biocontrol. This issue requires attention in the area of phage formulation to ensure phage viability and stability during storage and use.
The thermal stability tests showed that Pf-10 was highly stable at temperatures of 30-60 °C (ANOVA, p > 0.05), as the phage titer was similar to those of the control sample incubated at 4 °C (Figure 3b). Phage Pf-10 completely lost its activity after 1-h incubation at 70 °C and higher temperatures (ANOVA, p < 0.05), as no plaques were detected after the incubation under these temperature conditions (Figure 3b).

Killing Assay and the Effect of Ca 2+ and Mg 2+ on the Killing Activity of Pf-10
One of the important phage characteristics that can be useful for choosing phages for phage-based preparations is the killing activity, also called the lytic activity of a phage. Analysis of the growth kinetics of a bacterial host upon phage infection is the simplest method for assessing the effectiveness of the killing activity of a phage. The rate of the OD curve slope is directly associated with the lytic activity and the concentration of the phage used. Various phages are characterized by a relatively high or low efficiency of their lytic activity, depending on the time and the magnitude of the decrease in OD of the bacterial culture upon infection.
The killing activity of phage Pf-10 was assessed at MOI values of 0.1, 1 and 10 by measuring OD595 of infected P. syringae BIM B-268 cultures ( Figure 4). The non-infected cultures were used as a control ( Figure 4). Simultaneously, the effects of Ca 2+ and Mg 2+ on the killing activity of Pf-10 was assessed by adding CaCl2 and/or MgCl2 to the P. syringae BIM B-268 (Figure 4b-f) and comparing with the culture without the addition of the ions (Figure 4a,e,f). Figure 4 shows that the inhibitory effect of phage Pf-10 on bacterial growth increases significantly in a MOI-dependent manner compared to the control (Figure 4a-d; for each graph ANOVA, p < 0.05). The statistical analyses showed significant differences between bacterial growth curves under Pf-10 infection in the presence and absence of Mg 2+ and/or Ca 2+ ions (for each MOI ANOVA, p < 0.05). As shown in Figure 4, the lysis of bacterial cells is more efficient without the addition of the ions at each MOI. The PhageScore method based on bacterial biomass concentration changes, which are directly related to the lytic potential of the tested phage, was used to assess the lytic activity of phage Pf-10 under various salt conditions. The PhageScore method is a simple method recently developed by Konopacki and co-authors for comparative evaluation of lytic activity of phages, which The thermal stability tests showed that Pf-10 was highly stable at temperatures of 30-60 • C (ANOVA, p > 0.05), as the phage titer was similar to those of the control sample incubated at 4 • C (Figure 3b). Phage Pf-10 completely lost its activity after 1-h incubation at 70 • C and higher temperatures (ANOVA, p < 0.05), as no plaques were detected after the incubation under these temperature conditions (Figure 3b).

Killing Assay and the Effect of Ca 2+ and Mg 2+ on the Killing Activity of Pf-10
One of the important phage characteristics that can be useful for choosing phages for phage-based preparations is the killing activity, also called the lytic activity of a phage. Analysis of the growth kinetics of a bacterial host upon phage infection is the simplest method for assessing the effectiveness of the killing activity of a phage. The rate of the OD curve slope is directly associated with the lytic activity and the concentration of the phage used. Various phages are characterized by a relatively high or low efficiency of their lytic activity, depending on the time and the magnitude of the decrease in OD of the bacterial culture upon infection.
The killing activity of phage Pf-10 was assessed at MOI values of 0.1, 1 and 10 by measuring OD 595 of infected P. syringae BIM B-268 cultures ( Figure 4). The non-infected cultures were used as a control ( Figure 4). Simultaneously, the effects of Ca 2+ and Mg 2+ on the killing activity of Pf-10 was assessed by adding CaCl 2 and/or MgCl 2 to the P. syringae BIM B-268 (Figure 4b-f) and comparing with the culture without the addition of the ions (Figure 4a,e,f). Figure 4 shows that the inhibitory effect of phage Pf-10 on bacterial growth increases significantly in a MOI-dependent manner compared to the control (Figure 4a-d; for each graph ANOVA, p < 0.05). The statistical analyses showed significant differences between bacterial growth curves under Pf-10 infection in the presence and absence of Mg 2+ and/or Ca 2+ ions (for each MOI ANOVA, p < 0.05). As shown in Figure 4, the lysis of bacterial cells is more efficient without the addition of the ions at each MOI. The PhageScore method based on bacterial biomass concentration changes, which are directly related to the lytic potential of the tested phage, was used to assess the lytic activity of phage Pf-10 under various salt conditions. The PhageScore method is a simple method recently developed by Konopacki and co-authors for comparative evaluation of lytic activity of phages, which is measured in PhageScore factor values (P S ) [44]. The proposed method makes it possible to accurately compare and analyze the activity of bacteriophages, therefore it can be applied for comparative analysis and selection of a proper phage for specific phage-based preparations. We believe that this method will find application for assessing the phage lytic activity under various medium conditions. Estimation of phage killing activity in an environment, similar to that where the phage-based preparation will be used, is an important step in the selection of candidate phages. The PhageScore values (P S ) of Pf-10 upon infection at MOI of 0.1 and 1 are clearly reduced in the presence of MgCl 2 and CaCl 2 compared to the infection without additional Ca 2+ and Mg 2+ ions (Figure 4e). At MOI of 10, the Pf-10 P S values decrease slightly under various conditions (Figure 4e). In addition, the Pf-10 Total PhageScore values have also been calculated (Figure 4f). As described by Konopacki [44], the Total PhageScore factor (P ST ) is a parameter that characterizes the total lytic activity of a specific bacteriophage across the MOI range, being sensitive to the efficiency of the lowest dilution of the phage. Figure 4f shows that the Pf-10 P ST decreases in the presence of: Mg 2+ by 1.5 times, Ca 2+ by 1.6 times and both Mg 2+ and Ca 2+ by 2.1 times. Thus, the presence of Ca 2+ and Mg 2+ ions can reduce the killing activity of the phage Pf-10.

Adsorption Assay and One-Step Growth Curve
The data on the adsorption constant, adsorption rate, latent period, and burst size are the most important parameters that can affect the success of phage-based preparations. Understanding the interaction of phages and bacteria (including the parameters of phage infection, as well as the conditions under which the infection occurs) can help us optimize phage production, develop effective protocols for the use of phage-based preparations,

Adsorption Assay and One-Step Growth Curve
The data on the adsorption constant, adsorption rate, latent period, and burst size are the most important parameters that can affect the success of phage-based preparations. Understanding the interaction of phages and bacteria (including the parameters of phage infection, as well as the conditions under which the infection occurs) can help us optimize phage production, develop effective protocols for the use of phage-based preparations, determine the optimal dose concentrations and prognosticate the results of phage therapy [60][61][62].
The one-step growth curve assay showed the main parameters of the phage infection: the Pf-10 latent period (virion maturation) is 10-20 min, the duration of the virion release phase is about 30 min. The average burst size of the phage per infected cell is 37 ± 8.5 PFU. The latency period of the Pf-10 phage was nearly similar to Pseudomonas phage gh-1 (21 min) [65], and shorter than that of most Pseudomonas phages of the Autographiviridae family, namely Pseudomonas phage Andromeda (25-30 min) [63], Pseudomonas phage PAXYB1 (30 min) [64] and Pseudomonas phage RLP (25 min) [66], and especially of Pseudomonas phage BF7, which shows an extremely long latent period (140 min) [67]. Interestingly, the burst size of Pf-10 was smaller than that of other reported Pseudomonas phages of the Autographiviridae family such as Pseudomonas phage gh-1 (103 PFU) [65], Pseudomonas phage Andromeda (70 PFU) [63], Pseudomonas phage PAXYB1 (141 PFU) [64], Pseudomonas phage RLP (154 PFU) [66] and Pseudomonas phage BF7 (237 PFU) [67]. However, it should also be noted that all time and kinetic parameters of phage adsorption and one-step growth curve may differ depending on the bacterial host. It is not impossible that a small number of Pf-10 progeny particles are associated with the intracellular metabolism of the P. syringae BIM B-268 strain.

General Genome Organization
Genome sequencing and annotation have revealed that the Pf-10 genome is a linear dsDNA with the length of 39,167 bp (not including the second copy of the DTR) and the GC-content of 56.5%. The genome contains 49 genes (Supplementary Information, Table  S2). Thirty two out of 49 predicted CDSs (65.3%) were functionally assigned using BLASTp (NCBI) [49] and HHpred [50] ( Supplementary Information, Table S2). The schematic map of the phage genome is presented in Figure 6. The encoded proteins of phage The one-step growth curve assay showed the main parameters of the phage infection: the Pf-10 latent period (virion maturation) is 10-20 min, the duration of the virion release phase is about 30 min. The average burst size of the phage per infected cell is 37 ± 8.5 PFU. The latency period of the Pf-10 phage was nearly similar to Pseudomonas phage gh-1 (21 min) [65], and shorter than that of most Pseudomonas phages of the Autographiviridae family, namely Pseudomonas phage Andromeda (25-30 min) [63], Pseudomonas phage PAXYB1 (30 min) [64] and Pseudomonas phage RLP (25 min) [66], and especially of Pseudomonas phage BF7, which shows an extremely long latent period (140 min) [67]. Interestingly, the burst size of Pf-10 was smaller than that of other reported Pseudomonas phages of the Autographiviridae family such as Pseudomonas phage gh-1 (103 PFU) [65], Pseudomonas phage Andromeda (70 PFU) [63], Pseudomonas phage PAXYB1 (141 PFU) [64], Pseudomonas phage RLP (154 PFU) [66] and Pseudomonas phage BF7 (237 PFU) [67]. However, it should also be noted that all time and kinetic parameters of phage adsorption and one-step growth curve may differ depending on the bacterial host. It is not impossible that a small number of Pf-10 progeny particles are associated with the intracellular metabolism of the P. syringae BIM B-268 strain.

General Genome Organization
Genome sequencing and annotation have revealed that the Pf-10 genome is a linear dsDNA with the length of 39,167 bp (not including the second copy of the DTR) and the GC-content of 56.5%. The genome contains 49 genes ( Supplementary Information, Table S2). Thirty two out of 49 predicted CDSs (65.3%) were functionally assigned using BLASTp (NCBI) [49] and HHpred [50] ( Supplementary Information, Table S2). The schematic map of the phage genome is presented in Figure 6. The encoded proteins of phage Pf-10 are divided into several functional modules, the main of them being: DNA packaging and structure/morphogenesis module, lytic module, DNA replication and recombination module. It should be noted that the arrangement and structure of the Pf-10 modules resemble those in the Enterobacteria phage T7 ( Supplementary Information, Figure S3). Pf-10 are divided into several functional modules, the main of them being: DNA packaging and structure/morphogenesis module, lytic module, DNA replication and recombination module. It should be noted that the arrangement and structure of the Pf-10 modules resemble those in the Enterobacteria phage T7 ( Supplementary Information, Figure S3).

DNA Packaging Genes and the Genome Packaging Strategy
The DNA packaging module includes two proteins: a small terminase subunit (ORF45; protein_id YP_009145640.1) and a large terminase subunit (ORF47; protein_id YP_009145642.1). The critical step of the assembly of phages belonging to Caudovirales is known to be genome packaging performed by the large terminase subunit. In order to predict the type of genome termini and the packaging mechanism of phage Pf-10, the amino acid sequence of the large terminase protein was compared with the terminase protein from the published phages with well-studied DNA packaging mechanisms [68] (Supplementary Information, Table S4). As shown in Figure 7, the terminases of the phages: Pseudomonas phage Pf-10, Pseudomonas phage PFP1, Pseudomonas phage phi15, Enterobacteria phage T3, Yersinia phage phiYeO3-12 and Enterobacteria phage T7 form a separate clade, where the packaging protein of the phage Pf-10 is the closest to that of the phage PFP1 (Figure 7). Li M. et al. (2018) have shown that the phage PFP1 DNA has short direct terminal repeats of 270 bp in length [69]. Restriction analysis was performed to identify the type of genome termini of phage Pf-10, as well as to confirm the completeness and

DNA Packaging Genes and the Genome Packaging Strategy
The DNA packaging module includes two proteins: a small terminase subunit (ORF45; protein_id YP_009145640.1) and a large terminase subunit (ORF47; protein_id YP_009145642.1). The critical step of the assembly of phages belonging to Caudovirales is known to be genome packaging performed by the large terminase subunit. In order to predict the type of genome termini and the packaging mechanism of phage Pf-10, the amino acid sequence of the large terminase protein was compared with the terminase protein from the published phages with well-studied DNA packaging mechanisms [68] ( Supplementary Information, Table S4). As shown in Figure 7, the terminases of the phages: Pseudomonas phage Pf-10, Pseudomonas phage PFP1, Pseudomonas phage phi15, Enterobacteria phage T3, Yersinia phage phiYeO3-12 and Enterobacteria phage T7 form a separate clade, where the packaging protein of the phage Pf-10 is the closest to that of the phage PFP1 (Figure 7). Li, M. et al. (2018) have shown that the phage PFP1 DNA has short direct terminal repeats of 270 bp in length [69]. Restriction analysis was performed to identify the type of genome termini of phage Pf-10, as well as to confirm the completeness and correctness of the genome assembly. In silico restriction analysis of the Pf-10 genome was performed using the NEBcutter V2.0 tool [70] (Figure 8b). As shown in Figure 8, the restriction patterns for in vitro restriction analysis (Figure 8a) and in silico restriction analysis (Figure 8b) have some differences. It is clearly seen that in the case of restriction with HindIII, AfeI and NdeI, the DNA fragments containing the left end of the phage DNA ( Figure 8, indicated with yellow arrow) are longer in vitro compared to their size in the in silico restriction analysis. In order to determine the ends of the Pf-10 genome more accurately, RAGE method was used, followed by Sanger sequencing. The analysis of the sequencing results revealed the presence of two short DTRs at the ends of the Pf-10 DNA (Figure 8c). Due to limitations of the RAGE method, it was difficult to determine the final nucleotide at the right end of the Pf-10 DNA using the method, since the sequence in the DTR-containing region had two consecutive adenine nucleotides (Figure 8c and Supplementary Information, Figure S5). Sanger sequencing of the right end of the Pf-10 DNA clarified the location of the right terminus (shown in Figure 8c by the second arrow). According to the results of phylogenetic and restriction analyzes, as well as sequencing of Pf-10 DNA ends, the phage uses a packaging mechanism based on T7-like short DTRs (Figures 7 and 8), and the length of each terminal repeat in the phage chromosome is 257 bp (Figure 8d). Thus, the total length of the Pf-10 phage chromosome, including both DTRs, is 39,424 bp (Figure 8d).

Comparative Genomics
In order to assess the phylogenetic relationship between the phage Pf-10 and the known phages, the viral proteomic tree was constructed using the ViPTree server version 1.9. The viral proteomic tree showed that the closest relative of the phage Pf-10 was the Pseudomonas phage BIM BV-46, which, as with Pf-10, is part of the biopesticide "Multiphage" (Figure 9). According to the BLASTn (NCBI), the whole genome identity of Pf-10 and BIM BV-46 is 93.5% ( Supplementary Information, Table S5), indicating that they are the different phage species in accordance with the official ICTV classification (the species demarcation criterion is the genome nucleotide identity of 95%) [75]. As shown in Figure 9, phage Pf-10 forms a common giant branch with the Pseudomonas phages. However, Pf-10 has more than 70% of the whole genome identity with only eight phages from this huge branch, namely, phage BIM BV-46, phage Phi-S1, phage phiIBB-PF7A, phage PFP1, phage 22PfluR64PP, phage 67PfluR64PP, phage 71PfluR64PP, and phage UNO-SLW1 ( Figure 9; Supplementary Information, Table S5), which indicates that these phages belong to the same phage genus in accordance with the genus demarcation criterion of the official ICTV classification [75]. The phage DNA comparison diagram visualized using the ViPTree server version 1.9 shows the tBLASTx pairwise similarities between the most closely related phage genomic DNA ( Figure 10).

Conclusions
In summary, we have described a new bacteriophage Pf-10 isolated from a tissue sample of green bean (Phaseolus vulgaris) infected with the P. syringae pv. syringae. The Pseudomonas phage Pf-10, along with five other bacteriophages of the Pseudomonas genus, is a part of the biopesticide "Multiphage" [39], which is used in Belarus to control the bacterial diseases of agricultural crops caused by Pseudomonas syringae. The phage Pf-10 is the first phage from the "Multiphage" that has been sequenced and subjected to a thorough analysis, including TEM analysis, stability test, genome analyses (the genome organization and the genome packaging strategy), phylogenetic analysis and determination of phage infection parameters such as the rate of adsorption, the life cycle and size of phage progeny on the P. syringae BIM B-268. The results of the study contribute to the accumulation of data on phage biology that can help in formulation of the most efficient phage combinations in phage-based preparations to combat plant diseases caused by Pseudomonas. Thus, phage Pf-10 can be classified as a new phage species belonging to the genus Pifdecavirus, the subfamily Studiervirinae, the family Autographiviridae, and the order Caudovirales.

Conclusions
In summary, we have described a new bacteriophage Pf-10 isolated from a tissue sample of green bean (Phaseolus vulgaris) infected with the P. syringae pv. syringae. The Pseudomonas phage Pf-10, along with five other bacteriophages of the Pseudomonas genus, is a part of the biopesticide "Multiphage" [39], which is used in Belarus to control the bacterial diseases of agricultural crops caused by Pseudomonas syringae. The phage Pf-10 is the first phage from the "Multiphage" that has been sequenced and subjected to a thorough analysis, including TEM analysis, stability test, genome analyses (the genome organization and the genome packaging strategy), phylogenetic analysis and determination of phage infection parameters such as the rate of adsorption, the life cycle and size of phage progeny on the P. syringae BIM B-268. The results of the study contribute to the accumulation of data on phage biology that can help in formulation of the most efficient phage combinations in phage-based preparations to combat plant diseases caused by Pseudomonas.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/v14010042/s1. Table S1: Pseudomonas strains used in this study and infectivity of phage Pf-10; Figure S1: Transmission electron microscopy of the Pseudomonas phage Pf-10; Figure S2: Transmission electron microscopy of the Pseudomonas phage Pf-10; Table S2: Annotation of Pseudomonas phage Pf-10; Table S3: Primers used to locate the ends of the Pf-10 genome; Table S4: Large terminase subunit of Pf-10 and large terminase subunits of phages with well-known packaging mechanisms used for the phylogenetic tree; Figure S3: The pairwise whole-genome tBLASTx comparison of Pseudomonas phage Pf-10 and Enterobacteria phage T7 visualized with ViPTree server version 1.9; Table S5: Phage genomes used for phylogenetic inference; Figure S4: Restriction analysis of the phage genomic DNA; Figure S5: Schematic representation of the DTR-containing region of the Pf-10 genome after whole genome sequencing and genome assembly.