Isolation of the Novel Phage PHB09 and Its Potential Use against the Plant Pathogen Pseudomonas syringae pv. actinidiae

Bacteriophages are viruses that specifically infect target bacteria. Recently, bacteriophages have been considered potential biological control agents for bacterial pathogens due to their host specificity. Pseudomonas syringae pv. actinidiae (Psa) is a reemerging pathogen that causes bacterial canker of kiwifruit (Actinidia sp.). The economic impact of this pest and the development of resistance to antibiotics and copper sprays in Psa and other pathovars have led to investigation of alternative management strategies. Phage therapy may be a useful alternative to conventional treatments for controlling Psa infections. Although the efficacy of bacteriophage φ6 was evaluated for the control of Psa, the characteristics of other DNA bacteriophages infecting Psa remain unclear. In this study, the PHB09 lytic bacteriophage specific to Psa was isolated from kiwifruit orchard soil. Extensive host range testing using Psa isolated from kiwifruit orchards and other Pseudomonas strains showed PHB09 has a narrow host range. It remained stable over a wide range of temperatures (4–50 °C) and pH values (pH 3–11) and maintained stability for 50 min under ultraviolet irradiation. Complete genome sequence analysis indicated PHB09 might belong to a new myovirus genus in Caudoviricetes. Its genome contains a total of 94,844 bp and 186 predicted genes associated with phage structure, packaging, host lysis, DNA manipulation, transcription, and additional functions. The isolation and identification of PHB09 enrich the research on Pseudomonas phages and provide a promising biocontrol agent against kiwifruit bacterial canker.


Introduction
Pseudomonas syringae is a species complex of bacterial plant pathogens infecting more than 180 plant species [1], including phytopathogens that have significant effects on the agriculture sector by causing severe economic losses worldwide. Many globally important crops fall within the host range of P. syringae, making this group one of the most economically destructive pathogens [2]. Some kiwifruit cultivars, especially red-fleshed kiwifruit "Hongyang" and "Enza-Red" (Actinidia chinensis var. chinensis), are highly susceptible to P. syringae pv. actinidiae (Psa) [3,4], with infection causing bacterial canker and severe economic losses in many countries [5,6]. In recent studies, Psa can be divided into at least six different biovars [7]. Biovar 1 is associated with outbreaks of the disease in Japan and Italy [8]. Biovar 2 has only been reported in South Korea [9]. Biovar 3 derived from a Psa population from China, is responsible for substantial economic losses worldwide [10]. Biovar 4 is an avirulent strain [11], and Biovars 5 and 6 are identified in Japan [12,13]. Copper pesticide combined with antibiotics and strict orchard management have been used to control the disease and spread of Psa [14]. However, these strategies are not completely effective, and their widespread use can induce resistance in the pathogen and further environmental damage [15]. Therefore, a phage-based alternative strategy is required.

Bacterial Strains and Culture Conditions
The bacterial strain PsaBJ530 (CGMCC1.19157) used in this study was isolated previously from the branches of canker-infected, red-fleshed kiwifruit "Hongyang" cultivar in Cangxi County, Sichuan Province, and its pathogenicity has been confirmed in vivo (data not shown). The biovar of PsaBJ530 was determined by PCR using primer CAGGAATTCAT-GACTTCTCA and TAGTCTCGAAGATTCAATGG. Preservation of the strain was achieved using 15% glycerol stocks stored at − 80 • C. PsaBJ530 was subjected to culturing at 25 • C with shaking in tryptic soy broth (TSB) prior to phage infection.

Phage Isolation and Purification
Soil samples were collected from "Hongyang" cultivar kiwifruit orchards in Sichuan Province and used to screen phages specific to the isolated bacterial strain PsaBJ530. For phage isolation, 2 g soil was placed in a 50-mL Falcon tube containing 20 mL TSB and 4 sterile tungsten beads (Qiagen, Manchester, UK). The tube was vortexed for 5 min and centrifuged at 4500 rpm for 5 min. A 5-mL aliquot of the resulting supernatant was filtered through a 0.22-µm filter (Millipore, Sigma-Aldrich, Gillingham, UK) to remove any bacteria. The optical density (OD) of overnight bacterial cultures were measured at 600 nm and adjusted to an optical density of 0.8. Then, 100 µL of the bacterial host and 2 mL of a serially diluted environmental sample were added to a soft top agar overlay and incubated overnight at 25 • C. Three successive single-plaque isolations were performed to achieve pure phage isolates. The phage stocks were stored at 4 • C until use. DNA integrity and size distribution were assessed on a 1% (w/v) agarose gel and visualized with ethidium bromide [30].

Transmission Electron Microscopy
Transmission electron microscopy was used for determination of phage morphology. A filtered high-titer (~10 11 PFU/mL) phage stock solution was placed on the surface of carbon-coated copper grids (Agar Scientific, Essex, UK). The phage was negatively stained using 2% (w/v) uranyl acetate, which was removed after 2 min. The grids with adsorbed phages were air dried, and the morphology of the phage was observed by transmission electron microscopy (Hitachi H-7650, Tokyo, Japan) at an acceleration voltage of 80 kV using a charge-coupled device camera (AMT400, Woburn, MA, USA).

Determination of the Optimal Multiplicity of Infection and One-Step Growth Curve
The double agar overlay plaque assay was used with some modifications to determine phage titer. Briefly, 100 µL serially diluted phage culture was mixed with 100 µL PsaBJ530 cells (early log phase, OD 600 0.8). This mixture was added to 10 mL top agar (TSB with 0.7% agar), poured onto a plate (TSB with 1.5% agar), and incubated at 25 • C for 12-16 h to form phage plaques. Three parallel assessments were performed.
To determine the optimal multiplicity of infection (MOI), exponential-phase Psa cells were diluted to 10 9 CFU/mL. Phage solution (1 mL) was added to 1 mL of diluted bacterial cells at ratios of 1000, 100, 10, 1, 0.1, 0.01, 0.001, and 0.0001, respectively, and the mixture was shaken at 25 • C for 4 h at 200 rpm. Subsequently, the culture was centrifuged at 12,000× g for 10 min to remove bacterial cells, and the supernatant was filtered through a 0.22-µm pore-size membrane filter. The titer of the phage filtrate was determined using the double agar overlay plaque assay, as described previously. The dilution that generated the highest phage titer was considered the optimal MOI. Three parallel assessments were performed.
In the one-step growth curve assay, 15 mL bacterial cells (10 9 CFU/ mL) were mixed with 15 mL phage solution at an MOI of 0.001. The mixture was incubated at 25 • C for 5 min and centrifuged at 12,000× g for 30 s to remove unabsorbed free phages. The precipitate was suspended in TSB after two washes and then mixed with 30 mL TSB, followed by shaking at 25 • C and 180 rpm. This time was defined as t = 0 and subsequent time points as t = 15, 30, 45, 60, 75, 90, 105, 120, 135, and 150 min; at each time point, we collected 500 µL samples. The phage titer was determined using the double-layer agar method. The experiment was repeated three times, and three parallel tests were performed to measure phage titers at each time point.

Phage Host Range
The host range of phage PHB09 was investigated via infection of 22 bacterial strains using the spot assay method. Six biovar 2/3 Psa strains, including BJ530, BJ9, and BST isolated from kiwifruit orchards in Sichuan Province and three Psa strains from the Korean Agriculture Culture Collection, as well as other Pseudomonas sp. strains and three bacterial strains of other genera from China General Microbiological Culture Collection were tested. Bacteria in the log phase (100 µL) were mixed with 10 mL soft TSB agar (0.7% agar) and poured on top of a bottom layer containing 1.5% agar (15 mL). Then, the phage suspension was spotted onto the surface of double-layer agar plates containing lawns of the target bacterial strains. The plates were incubated at 25 • C for 12-18 h, and plaque formation was observed. Bacterial sensitivity to a phage was established based on the lysis-cleared zone around the spot. The results were classified into two categories: clear lysis zone and no lysis zone.

Stability of the Phage under Various Thermal, Ultraviolet, and pH Conditions
To evaluate the thermal stability of bacteriophages, phage preparations (10 10 PFU/mL) in TSB broth were incubated in a 4 • C, 25 • C, 37 • C, or 50 • C water bath. To evaluate the ultraviolet (UV) stability of bacteriophages, phage preparations (10 10 PFU/mL) in TSB broth were illuminated with a UV lamp (365 nm, 18 µW/cm 2 ) for 0, 5, 15, 30, 45, or 60 min. Three samples were serially diluted. Subsequently, their titers were determined via the double-layer agar method, and the samples were placed in an incubator at 25 • C for 12-18 h. Three parallel experiments were conducted. To estimate pH stability, TSB was adjusted to various pH values (3.0, 4.0, 5.0, 6.0 7.0, 8.0, 9.0, 10.0, and 11.0). Next, 100 µL phage (10 10 PFU/mL) was added to 900 µL TSB broth at each pH and incubated at 25 • C for 1 h.

In Vitro and In Vivo Phage Efficacy in Psa Control
To examine the lytic activity of PHB09 in vitro against a target bacterium, a bacteriophage aliquot was added to bacterial solution in the exponential phase (PsaBJ530 strain, 10 9 CFU/mL) to achieve an MOI of 1. The solutions were mixed and incubated at 25 • C for 48 h with shaking (200 rpm). As the negative control, bacterial culture without phage was inoculated with the same volume of TSB medium. Bacterial growth was monitored using a TECAN microplate reader (TECAN, Männedorf, Switzerland). The OD 600 was measured for 48 h. The lytic activity of PHB09 was determined based on the phage titer at the same time points. Three parallel tests were performed, and each experiment was conducted in triplicate.
In vivo experimentation to evaluate phage efficacy was performed with leaf discs (≈8 cm diameter). The discs were obtained from healthy kiwifruit plants. Leaf surface was disinfected with sodium hypochlorite 1% and washed twice with sterile distilled water. Subsequently, the leaf discs were placed in humidity chambers and inoculated individually with three drops of 100 µL of PsaBJ530 (10 9 CFU/mL). In Psa + phage group, PHB09 was added two hours after Psa infection with a MOI of 1. Leaf discs were deposited in a plate containing 20 mL of sterile distilled water supplemented with cycloheximide (100 µg/mL) to avoid fungus growth [25]. Leaf samples from each group were collected at 0, 12, 24, 48 and 72 h after infection. For each sample, 0.5 cm 2 leaf tissue was homogenized in 1 mL sterile TSB medium to determine Psa concentration (CFU/mL) and phage titer (PFU/mL). All experiments were performed in biological and technical triplicate.

DNA Extraction
The phage lysates were centrifuged at 8000× g for 5 min to remove cellular debris. The supernatant, containing the majority of viral particles, was filtered through a 0.22-µm syringe filter to remove cellular debris. The treated lysate was concentrated by centrifuging in a 100-kDa molecular weight cutoff ultrafiltration centrifuge tube (Amicon Ultra-15 centrifugal filter units; Millipore, MA, USA) at 5000× g to a final volume of 1 mL. The Nucleic Acid Extraction Kit II (Geneaid Biotech Ltd., Taiwan, China) was used to extract phage DNA from a high-titer plate lysate (minimum of 10 9 PFU/mL). The extracted DNA was stored in a 1.5-mL ultracentrifuge tube at −20 • C until needed. DNA quality was evaluated by agarose gel electrophoresis.

Genome Sequencing and Bioinformatics Analysis
Extracted PHB09 genomic DNA was sequenced on an Illumina sequencer (Illumina, San Diego, CA, USA). Raw reads were trimmed using Trimmomatic version 0.36 (parameters: version 0.36, illuminaclip: TruSeq3-PE.fa:2:30:10 leading:3 trailing:3 slidingwindow:4:15 minlen:40) [31] to obtain clean reads, which constituted more than 90% of the raw reads. Bowtie2 version 2.3.4 [32] was used to remove sequences of the host bacterial genome; high-quality clean reads were then assembled using IDBA-UD version 1.1.3 (parameters: kmer min 21, max 91, and Step 10) [33]. The final assemblies were filtered to obtain 2923 contigs. The filtered contigs were compared with clean reads using samtools v1.8. The GC content and average sequencing depth of the contigs were calculated. Contigs with avgDepth ≥ 100 and length ≥ 5 kb were selected, followed by the alignment to viral genomes with covering ≥50%, identity ≥80% by NCBI BLASTn (https://blast.ncbi.nlm.nih.gov/Blast.cgi, accessed on 20 October 2021). The circular viral genome was obtained according to the overlap of the two terminals of the sequence. The genes from assembled genomic sequences were predicted using Gene-MarkS [34] (http://topaz.gatech.edu/GeneMark/genemarks.cgi, accessed 9 April 2021) and RAST [35] (http://rast.nmpdr.org/, accessed 9 April 2021). The program tRNAscan-SE was used to predict tRNA sequences [36]. The putative protein function associated with each open reading frame (ORF) was manually verified by searching the NCBI nonredundant protein sequence and conserved domain databases using the BLASTp, with the e-value to <1.0 × 10 −5 . The whole genome was compared with other nucleotide sequences using NCBI BLASTn (https://blast.ncbi.nlm.nih.gov/Blast.cgi, accessed on 20 October 2021). Subsequently, the average nucleotide identity was determined using the BLASTn alignment tool in the pyani package, and an interactive heatmap was constructed using heatmaply [37,38]. A genomic map was generated using CGview (http://cgview.ca/, accessed on 20 October 2021). The complete annotated genome sequence of phage PHB09, displayed in Table S1, was deposited in the NCBI nucleotide database under accession number OK040171. Easyfig v.2.0 (https://mjsull.github.io/Easyfig, accessed on 21 October 2021) was used to assess genomic organization and identify core genes in a given group of bacteriophages [39]. ViPTree (https://www.genome.jp/viptree, accessed on 6 November 2021) was used to generate viral proteomic trees for classification of viruses based on genome-wide similarities [40].

Phylogenetic Analysis
To elucidate the genome of PHB09, the large subunit of terminase and major capsid protein were compared phylogenetically with those of other bacteriophages (Table 1) using Mega-X software (version 10.1.6). ClustalX was used to align the inferred amino acid sequences using the default parameters. Based on the multiple-sequence alignment, the Jones-Taylor-Thornton model was employed for construction of a maximum likelihood tree with 100 bootstrap replicates.

Statistical Analysis
Data analysis was performed using GraphPad Prism 8.0 (GraphPad Software, San Diego, CA, USA). Multiple t-tests were performed to determine the differences between groups at a significance level of p < 0.05.

Isolation and Morphology of Phage PHB09
A new Pseudomonas phage, vB_PsyM-PHB09 (designated phage PHB09 hereafter), was isolated from a kiwifruit orchard in Sichuan, China. The phage was purified by successive single-plaque isolations using the double-layer agar technique. Phage PHB09 formed clear plaques approximately 2 mm in diameter on a bacterial lawn of the host strain PsaBJ530 (Figure 1a). Enzyme digestion showed that PHB09 is a DNA phage, as it could be digested with DNase but not RNase (Figure 1b). Electron microscopy showed that PHB01 has an icosahedral head (height, 55.2 nm ± 1.0 nm; width, 54.5 nm ± 1.5 nm) and contractile tail (length, 145.0 nm ± 1.0 nm; width, 14.0 ± 0.5 nm) (Figure 1c).
Based on these morphological characteristics and according to the latest International Committee on Taxonomy of Viruses classification, PHB09 exhibited typical head and tail morphologies associated with the class Caudoviricetes (tailed ds-DNA phages). Viruses in all of these families have icosahedral or oblate heads but differ in the length and contractile abilities of their tails [41]. Electron microscopy showed that PHB09 was a myovirus [42]. Numerous studies have examined phages that infect Pseudomonas species, and more than 97% of those observed using electron microscopy belong to the Caudovirales [43]. Rombouts et al. tested a phage cocktail of five Myoviridae phages against 41 different strains of P. syringae pv. porri on leek leaves [44]. Frampton et al. characterized 24 Psa bacteriophages isolated from soil, water, and leaf samples collected in infected kiwifruit orchards, among which 22 belonged to Myoviridae, 1 to Podoviridae, and 1 to Siphoviridae [14]. Di Lallo et al. reported two new bacteriophages, phiPSA1 (Siphoviridae) and phiPSA2 (Podoviridae), which were isolated from kiwifruit leaves infected with Psa [45]. Based on these morphological characteristics and according to the latest International Committee on Taxonomy of Viruses classification, PHB09 exhibited typical head and tail morphologies associated with the class Caudoviricetes (tailed ds-DNA phages). Viruses in all of these families have icosahedral or oblate heads but differ in the length and contractile abilities of their tails [41]. Electron microscopy showed that PHB09 was a myovirus [42]. Numerous studies have examined phages that infect Pseudomonas species, and more than 97% of those observed using electron microscopy belong to the Caudovirales [43]. Rombouts et al. tested a phage cocktail of five Myoviridae phages against 41 different strains of P. syringae pv. porri on leek leaves [44]. Frampton et al. characterized 24 Psa bacteriophages isolated from soil, water, and leaf samples collected in infected kiwifruit orchards, among which 22 belonged to Myoviridae, 1 to Podoviridae, and 1 to Siphoviridae [14]. Di Lallo et al. reported two new bacteriophages, phiPSA1 (Siphoviridae) and phiPSA2 (Podoviridae), which were isolated from kiwifruit leaves infected with Psa [45].

Biological Characterization of Phage PHB09
Phage PHB09 generated a maximum titer of 1.91 × 10 11 PFU/mL at an optimal MOI of 0.001 (Figure 2a). To determine the latent period and burst size of phage PHB09, a onestep growth curve was plotted at 25 °C. The latent period was approximately 60 min and the rise period approximately 40 min (Figure 2b). Furthermore, the burst size of phage PHB09 was estimated at 182 PFU/infected cell, indicating that phage PHB09 replicates efficiently in Psa but requires a long latency period.

Biological Characterization of Phage PHB09
Phage PHB09 generated a maximum titer of 1.91 × 10 11 PFU/mL at an optimal MOI of 0.001 (Figure 2a). To determine the latent period and burst size of phage PHB09, a one-step growth curve was plotted at 25 • C. The latent period was approximately 60 min and the rise period approximately 40 min (Figure 2b). Furthermore, the burst size of phage PHB09 was estimated at 182 PFU/infected cell, indicating that phage PHB09 replicates efficiently in Psa but requires a long latency period.

Thermal, Ultraviolet, and pH Stability Tests
The stability of lytic activity against target bacteria is important for phage therapy. Thermal stability tests showed that PHB09 activity was high following water incubation

Thermal, Ultraviolet, and pH Stability Tests
The stability of lytic activity against target bacteria is important for phage therapy. Thermal stability tests showed that PHB09 activity was high following water incubation from 4 • C to 37 • C for 12 h, after which it began to decrease. However, the maximum decrease was only 1.47 log PFU/mL after 24 h when the samples were held at a temperature of 37 • C. When the temperature was increased to 50 • C, the phage titer decreased by 2.08 log PFU/mL after 24 h of incubation. Compared with the samples at 4 • C and 25 • C, phage titers decreased significantly at 37 • C and 50 • C after 6 h treatment (Figure 3a). UV stability tests showed that PHB09 retained high levels of activity from 0 to 60 min treatment. With increasing UV exposure time, the phage titer gradually decreased. When phage PHB09 was exposed to UV for 60 min, a titer reduction of 2.13 log PFU/mL was observed (Figure 3b). To evaluate pH stability, the phage was treated with acidic, neutral, and basic TSB medium for 1 h. PHB09 remained stable over a pH range of 3-11. However, it showed significant reductions at pH 3 and 11, with titer decreases of 0.97 and 0.88 log PFU/mL, respectively (Figure 3c). Environmental factors exert significant effects on phage stability, impeding the efficacy of phage treatment. For bacteriophages to be used as biocontrol agents, their lytic activity must be stable under various environmental conditions such as temperature, UV light, and soil pH [16]. Temperature plays a fundamental role in the attachment, penetration, and amplification of phage particles within their host cells [46]. At low temperatures, genetic material from only a few phages enters bacterial host cells, and therefore fewer  Environmental factors exert significant effects on phage stability, impeding the efficacy of phage treatment. For bacteriophages to be used as biocontrol agents, their lytic activity must be stable under various environmental conditions such as temperature, UV light, and soil pH [16]. Temperature plays a fundamental role in the attachment, penetration, and amplification of phage particles within their host cells [46]. At low temperatures, genetic material from only a few phages enters bacterial host cells, and therefore fewer phage particles are involved in the multiplication phase. On the other hand, high temperatures promote an extended phage latency period. In our study, Phage PHB09 was stable over a wide range of temperatures from 4 • C to 37 • C, which is a significant advantage for future applications. Thus, in summer, when temperatures occasionally rise to 37 • C, the phage maintains its activity. As the most critical period for plants is from autumn through winter to early spring, and the infective ability is reduced at temperatures above 25 • C [47], temperature is not a problem for implementation of this phage in therapy. Solar radiation or UV irradiation can directly affect free viruses via degradation of proteins, alteration of the viral structure, and reduction of infectivity [48]. The abundance of phage PHB09 particles decreased when exposed to UV radiation. However, sensitivity to UV wavelengths in solar radiation can be overcome by applying the phage at high titers and late in the day or at night, a period during which radiation is limited [49]. pH is another important factor affecting phage attachment, infectivity, intracellular replication, and amplification [50]. Unfavorable pH values can interfere with the lysozyme enzyme or with other phage capsid proteins, thereby preventing phage attachment to receptor sites on the host cell. Phage PHB09 exhibited high pH tolerance, as indicated by maintenance of a stable titer across a pH range of 3-11. In addition, kiwifruit grow optimally at pH 5.5-6.0, and the temperature at which kiwifruit canker grows in Sichuan Province is below 30 • C. Based on these environmental conditions, PHB09 bacteriophages are likely to be stable in the natural growth environment of kiwifruit.

Host Range of Phage PHB09
The host range of phage PHB09 was investigated by infecting 22 bacterial strains of different genera and species using the spot assay method. The results, shown in Table 2, indicated that phage PHB09 could lyse not only biovar 3 Psa strains BJ530, BJ9, and BST isolated in the laboratory but also other biovar 2 Psa strains from the Korean Agriculture Culture Collection that are considered highly virulent plant pathogens and causative agents of kiwifruit losses worldwide. Other tested Pseudomonas sp. strains were not infected.  Lytic phages can kill the host bacterium and play an important role in maintaining biodiversity [14]. With advantages such as target specificity and rapid self-replication, phages are being explored as an alternative strategy for the treatment of Pseudomonas infections [51]. In the host range experiment, PHB09 showed a narrow host range, with a limited ability to infect Psa and an inability to infect other Pseudomonas sp. strains. This result indicates that PHB09 is a promising biocontrol agent that can specifically inactivate Psa strains.

In Vitro and In Vivo Efficacy of Phage PHB09
To further assess the application potential of phage PHB09, the phage was used to inhibit the growth of Psa in vitro. To examine the duration of PHB09 lytic activity against the target bacterium in vitro, the cell density (OD 600 ) was measured over 48 h (Figure 4a). The density of PsaBJ530 without phage increased by 0.49 after 48 h. The density of bacteria treated with the phage decreased gradually, to 0.29 at 12 h, and then slowly increased until 48 h. After 48 h of phage treatment, the bacterial concentration was significantly lower than that in non-treated cultures.
To test the efficacy of phage PHB09 against bacterial canker in vivo, kiwifruit leaves were treated by PHB09 two hours after Psa infection. The results showed that PHB09 could reduce the Psa load over kiwifruit leaves 24-72 h post-infection (Figure 4b). Phage PHB09 was detectable throughout the entire experiment, showing a significant titer increase at 12 h (Figure 4b). After 72 h, leaf infected with Psa had typical symptoms (browning and cankering), while leaf infected with Psa and treated with PHB09 showed no symptoms (Figure 4c). These results indicate that the treatment of phage PHB09 can efficiently control bacterial canker in kiwifruit.
Based on these results, PHB09 can maintain stable lytic activity against PsaBJ530 both in vitro and in vivo, highlighting the potential of PHB09 in the biological control of Psa infection. The lytic effect pattern is very similar to those in previous reports of Psa phages (KHUϕ34, KHUϕ38, and PPPL-1) [11,28]. This similarity indicates that PHB09 can be considered a promising candidate for phage therapy, and that its characterization in this study may provide a starting point for further exploration of its potential in the biological control of bacterial canker of kiwifruit. Due to the systemic nature of bacterial canker of kiwifruit [52], phage therapy is best used as a preventative strategy. As phage particles can move in moist environments such as plant tissues in addition to infecting bacteria that are on the plant's surface, they can be used to assess and control bacteria within a plant's tissues during an infection. However, the natural environment includes water, wind, rain, and sunlight, which influence the success of phage treatment, and no studies to date have reported the effects of these factors. In the future, protective formulations and carrier bacteria [53] must be identified to improve viral survival. Further studies are needed for this purpose.

In Vitro and In Vivo Efficacy of Phage PHB09
To further assess the application potential of phage PHB09, the phage was used to inhibit the growth of Psa in vitro. To examine the duration of PHB09 lytic activity against the target bacterium in vitro, the cell density (OD600) was measured over 48 h (Figure 4a). The density of PsaBJ530 without phage increased by 0.49 after 48 h. The density of bacteria treated with the phage decreased gradually, to 0.29 at 12 h, and then slowly increased until 48 h. After 48 h of phage treatment, the bacterial concentration was significantly lower than that in non-treated cultures.

Genome Characterization and Comparative Genomic Analysis
The 94,844-bp genome of PHB09 is circular, with a GC content of 57.61% (accession number: OK040171). In total, 186 genes were predicted, and no tRNA genes were identified. The majority of the predicted genes were detected on the forward strand, accounting for 74.2% of the total (138 of 186). Thirty-four genes had significant similarity to genes with known functions, while over 80% of all predicted genes were hypothetical proteins that could not be annotated to any homologs (BLASTp; e-value cutoff 10 −5 ). Based on the results of BLASTp and Conserved Domain Database analyses, 34 genes were annotated to encode functional proteins (Table 3), and their arrangement at the whole-genome level was mapped ( Figure 5). Although the complete genome of phage PHB09 has been analyzed, many of its proteins have unknown functions. Genome sequencing suggested that no previously described antibiotic resistance or virulence factor genes are present in the genome of phage PHB09; this characteristic increases the safety of this phage for agricultural application (Table S1). Proteomic tree of PHB09 was generated based on genomic similarity scores (denoted by SG) derived from tBLASTx scores (Figure 6a). The viral proteomic tree revealed that phage PHB09 has a high sequence identity to Pseudomonas phage VCM, belonging to genus Otagovirus (query coverage 37%; identity 78.75%; accession number LN887844.1), and to Pseudomonas phage vB_PsyM_KIL5, belonging to genus Flaumdravirus (query coverage 51%; identity 88.29%; accession number KU130130.1) (Figure 6b). The whole genome of PHB09 was colinear with those two Psa phage genomes, with low similarity in terms of genomic architecture (Figure 7). The heatmap results were similar to the NCBI search results based on the complete genome of phage PHB09, indicating that phage PHB09 has low genomic similarity to known phages ( Figure 8). Proteomic tree of PHB09 was generated based on genomic similarity scores (denoted by SG) derived from tBLASTx scores (Figure 6a). The viral proteomic tree revealed that phage PHB09 has a high sequence identity to Pseudomonas phage VCM, belonging to genus Otagovirus (query coverage 37%; identity 78.75%; accession number LN887844.1), and to Pseudomonas phage vB_PsyM_KIL5, belonging to genus Flaumdravirus (query coverage 51%; identity 88.29%; accession number KU130130.1) (Figure 6b). The whole genome of PHB09 was colinear with those two Psa phage genomes, with low similarity in terms of genomic architecture (Figure 7). The heatmap results were similar to the NCBI search results based on the complete genome of phage PHB09, indicating that phage PHB09 has low genomic similarity to known phages ( Figure 8).
The Bacterial and Archaeal Viruses Subcommittee [54] defines a genus as a cohesive group of viruses sharing a high degree of nucleotide identity (>50%) and being distinct from viruses in other genera. Thus, phage PHB09 might be classified as belonging to a novel genus.     The Bacterial and Archaeal Viruses Subcommittee [54] defines a genus as a cohesive group of viruses sharing a high degree of nucleotide identity (> 50%) and being distinct October 2021). The average nucleotide identity (ANI) was valued. Values range from 0 (0%) ANI to 1 (100% ANI): gray represents 0% ANI; clusters of highly similar phages are highlighted in pink and red.

Phylogenetic Analysis
Two phylogenetic trees were constructed using the amino acid sequences of the predicted large subunit of terminase and the major capsid protein, which are often the most conserved sequences in phage genomes. Phage PHB09 is highly homologous to Pseudomonas phage VCM based on the terminase large subunit, and this sequence was distinct from those of other groups (Figure 9a). Phylogenetic analysis using the sequences of the major capsid protein placed PHB09 on a separate branch, displaying a distant phylogenetic relationship to the Otagovirus genus (Figure 9b). According to phylogenetic analysis of the complete genome sequence and amino acid sequences of conserved proteins (major capsid protein and terminase large subunit), phage PHB09 belongs to a new phage lineage, indicating that PHB09 is a novel genus in the class Caudoviricete [55]. from those of other groups (Figure 9a). Phylogenetic analysis using the sequences of the major capsid protein placed PHB09 on a separate branch, displaying a distant phylogenetic relationship to the Otagovirus genus (Figure 9b). According to phylogenetic analysis of the complete genome sequence and amino acid sequences of conserved proteins (major capsid protein and terminase large subunit), phage PHB09 belongs to a new phage lineage, indicating that PHB09 is a novel genus in the class Caudoviricete [55].

Conclusions
In conclusion, this study describes the isolation, characterization, and application of vB_PsyM-PHB09, a novel DNA phage against Psa. It exhibits a short latent period, large burst size, and stability across a broad range of pH and temperature conditions. Both morphological and genetic analyses indicated that PHB09 is a member of class Caudoviricete.

Conclusions
In conclusion, this study describes the isolation, characterization, and application of vB_PsyM-PHB09, a novel DNA phage against Psa. It exhibits a short latent period, large burst size, and stability across a broad range of pH and temperature conditions. Both morphological and genetic analyses indicated that PHB09 is a member of class Caudoviricete. However, PHB09 does not belong to any genera previously identified in myovirus and should be assigned to a new genus. This study of PHB09 enriches the research on Pseudomonas phages. Moreover, PHB09 exhibited a strong ability to kill Psa strains both in vitro and in vivo, suggesting the potential for biocontrol of Psa in kiwifruit agriculture.