Characterization and Full Genome Sequence of Novel KPP-5 Lytic Phage against Klebsiella pneumoniae Responsible for Recalcitrant Infection

Klebsiella pneumoniae is a hazardous opportunistic pathogen that is involved in many serious human diseases and is considered to be an important foodborne pathogen found in many food types. Multidrug resistance (MDR) K. pneumoniae strains have recently spread and increased, making bacteriophage therapy an effective alternative to multiple drug-resistant pathogens. As a consequence, this research was conducted to describe the genome and basic biological characteristics of a novel phage capable of lysing MDR K. pneumoniae isolated from food samples in Egypt. The host range revealed that KPP-5 phage had potent lytic activity and was able to infect all selected MDR K. pneumoniae strains from different sources. Electron microscopy images showed that KPP-5 lytic phage was a podovirus morphology. The one-step growth curve exhibited that KPP-5 phage had a relatively short latent period of 25 min, and the burst size was about 236 PFU/infected cells. In addition, KPP-5 phage showed high stability at different temperatures and pH levels. KPP-5 phage has a linear dsDNA genome with a length of 38,245 bp with a GC content of 50.8% and 40 predicted open reading frames (ORFs). Comparative genomics and phylogenetic analyses showed that KPP-5 is most closely associated with the Teetrevirus genus in the Autographviridae family. No tRNA genes have been identified in the KPP-5 phage genome. In addition, phage-borne virulence genes or drug resistance genes were not present, suggesting that KPP-5 could be used safely as a phage biocontrol agent.


Introduction
Klebsiella pneumoniae, a Gram-negative bacterium, a facultative anaerobic, is the primary cause of serious infections in humans around the world [1]. It is well-known that many nosocomial infections are caused by it, leading to increased morbidity and mortality rates [2,3]. In addition, K. pneumoniae has been considered to be an important food-borne pathogen found in many food types [4][5][6]. Unfortunately, K. pneumoniae exhibits a high resistance to many antimicrobials due to the presence of resistance genes, and many other resistances are encoded in plasmids that make treatment difficult [7].
In many habitats, including the human body, bacteria can form biofilms to face the adverse impacts of environmental challenges [8]. Biofilms are known as surface-attached aggregates embedded in an extracellular polymeric substance matrix, consisting of polysaccharides, proteins, enzymes, DNA, lipids, and water [8,9]. The ability of K. pneumoniae to the replating of captured plaques was frequently performed, and all dilutions for the phage titer determination were performed. The experiment was conducted in triplicate.

Host Range Determination
Host spectrum of isolated bacteriophages was determined by the spot-test assay according to Huang et al. [28] and Jun et al. [29] on 31 bacterial strains from different sources. The bacteriolytic activity was assayed against 19 strains of Klebsiella pneumoniae (clinical and food isolates), six of Escherichia coli (clinical isolates), three of Salmonella Typhimurium (food isolates), and 3 of Pseudomonas aeruginosa (clinical and food isolates). Three milliliters of soft agar (0.7% agar) medium was inoculated with 100 µL of test strain and poured onto Trypticase Soy Agar (TSA, Difco TM , BD, Franklin Lakes, NJ, USA) plates. The plates were left to solidify, and 10 µL of phage lysate was dropped on the surface of the previous plates and incubated at 37 • C for 24 h. After the incubation period, areas of clear spots were observed on the bacterial lawn. The host range of the most bacteriolytic phage (broadest host range) was also confirmed by the efficiency of plating method (EOP) in accordance with Mirzaei and Nilsson [30].

Transmission Electron Microscopy (TEM)
Stock KPP-5 phage was centrifuged at 16,000× g for 60 min and 4 • C after washing in a CM buffer. Supernatant was discarded, and the pellet was gently suspended in 20 µL of CM buffer; then, 5 mL was taken and placed on the carbon grids (200 mesh) coated with formvar and left for 2 min. The phage was stained by 2% uranyl acetate (negative stain), and the excess of the dye was removed by filter paper [31]. Sample examination in accordance with Accolas and Spillmann [32] was done at the Regional Center of Mycology and Biotechnology, Al-Azhar University, Cairo, Egypt using electron microscopy (Model Beckman 1010, Operated at 80 KV).

Determination of KPP-5 Phage Bacterial Culture Clearance Tendency
The potential effect of KPP-5 phage on the reduction of bacterial growth in the broth medium was evaluated using the Khawaja et al. [33] method. A host culture was incubated for 24 h at 37 • C; then, the culture (1 × 10 8 CFU/mL) was added to 2 sterile clean flasks (100 mL) after incubation. KPP-5 phage (1 × 10 8 PFU/mL) was added to one flask (test flask), while the second flask remained the same (positive control flask). In addition to the presence of a third flask containing only TSB media (negative control flask), it is indicated during the observations of the experiment. Flasks were incubated overnight in shaking (220 rpm) at 37 • C, while optical density (OD 620 ) was recorded during the incubation period at 2, 4, 6, 8, 10, 12, and 24 h. The experiment was conducted in triplicate.

One-Step Growth Curve
The KPP-5 phage one-step growth curve against the host was carried out in accordance with Wang et al. [34]. KPP-5 phage (1 × 10 8 PFU/mL) was added to a host fresh culture (1 × 10 6 CFU/mL) at a multiplicity of infection (MOI) of 100. The mixture was placed in the incubator at 37 • C for 3 min, then centrifuged at 12,000× g for one minute to remove the free phages. The pellet was resuspended (time zero) in 100 mL of TSB medium and incubated at 37 • C. For 3 h, 5-min interval samples were taken, centrifuged, and the KPP-5 phage titer was calculated using the agar overlay method [27]. Burst size was detected by determining the ratio of the mean of virions released after bacterial infection to the mean of virions used at the beginning of the host infection. The experiment was conducted in triplicate.

Phage DNA Isolation and Genome Sequencing of KPP-5 Phage
Highly concentrated purified lysate was treated with RNase and DNase I Sigma Aldrich to digest any RNA or exogenous DNA. Then phage lysate was concentrated, and DNA was extracted using the Gene JET genomic DNA purification kit (Cat No. K0721, Thermo Scientific) according to the manufacturer's instructions. Phage was whole-genome fragmented using restriction enzymes from the Ion Xpress plus library preparation kit (Thermo Fisher, USA); then, fragments were Nick repaired and adaptors ligated, and prepared libraries were barcoded with Ion Xpress Barcode Adaptors. After barcoding, a 200-bp library size was selected and purified using Agencourt beads (Beckman Coulter, Brea, CA, USA), and libraries were quantified using the Ion Library TaqMan

Bioinformatic Analysis of KPP-5 Phage Genome
Sequence reads were de novo assembled using trial Geneious prime software (version 11.0.6) using Geneious assembler with the default settings. The resulting scaffolds were filtered based on the length of the scaffolds. The whole assembly processes were repeated using unicycler workflow on galaxy (Galaxy Version 0.4.6.0) to ensure the reproducibility of the results with another assembler (SPAdes assembler 3.13.0). Prediction of genes was achieved by predicting open reading frames (ORFs) using ATG as starting codons with a minimum nucleotide length of 300 bp using the ORF finder in the National Center for Biotechnology Information (NCBI) database. The ORFs were translated and analyzed by SmartBlast and BLASTP to describe the potential functions and were further confirmed using the Pfam database [37]. SnapGene software was used to create a genome map of KPP-5 phage [38]. All ORFs were checked for genes encoding putative tRNAs using tRNA scan-SE version 2.0 (https://lowelab.ucsc.edu/tRNAscan-SE/) [39] and checked for genes encoding the potential protein-containing transmembrane domain using the TMHMM server v. 2.0 (https://www.cbs.dtu.dk) [40] with the default settings. ORFs were also analyzed using the databases of Virulence Searcher (http://www.mgc.ac.cn/ VFs/) [41], ResFinder (https://cge.cbs.dtu.dk/services/ResFinder/) [42], and Antibiotic Resistance Genes (https://card.mcmaster.ca/) [43] to examine the virulence factor and drug resistance gene of the KPP-5 phage. The full-genome sequence of the KPP-5 phage was compared to the phage genome sequences in GenBank using BLASTN (somewhat similar sequences) in the NCBI database, and the most closely related phages were identified. The average nucleotide identity (ANI) values were calculated using the ANI calculator (http://enve-omics.ce.gatech.edu/ani/). The Circoletto program (http://tools.bat.infspire. org/circoletto/) [44] was used to visualize the phage comparative genome. CoreGenes 3.5 (http://binf.gmu.edu:8080/CoreGenes3.5/custdata.html) was used for the analysis of the phage-core genes where genes with a score greater than 75 were assumed to be the core genes [45]. The full-genome phylogenetic tree was created using the genome-BLAST distance phylogeny method in the Virus Classification and Tree Building Online Resource (VICTOR) [46]. In order to determine the taxonomic position of KPP-5 and to investigate its evolutionary history, a phylogenetic tree was constructed on the basis of conserved proteins frequently used in the analysis of the diversity, such as DNA polymerase, major capsid protein, and terminase large subunit. The clustalW option and the neighbor-joining method with 1000 bootstrap replicates implemented in Molecular Evolutionary Genetic Analysis (MEGA) software version 7.0 were used for multiple sequences alignment and phylogenetic tree construction, respectively.

Broad Host Range Bacteriophage Isolation and Morphology
Six bacteriophages against multidrug resistance (MDR) Klebsiella pneumoniae were isolated from sewage water samples of the Wastewater Treatment Plant (Kafr El-Sheikh Governorate, Egypt). Only one bacteriophage isolate named KPP-5 showed clear plaques small in size (1-1.5 mm in diameter) that looked like dots ( Figure 1A). The KPP-5 bacteriophage was purified, propagated, and titrated for further analysis. KPP-5 phage showed potent lytic activity, being able to infect all K. pneumoniae strains of different sources, with a high efficiency of plating (EOP) activity on 14 out of 19 strains (Table 1). On the other hand, KPP-5 could not infect any of the other Escherichia coli, Salmonella Typhimurium, and Pseudomonas aeruginosa strains tested (Table 1). Klebsiella pneumoniae is a dangerous opportunistic pathogen that is involved in many serious human diseases [47,48]. In addition, K. pneumoniae has been considered to be an important food-borne pathogen found in many food types [4][5][6]. Recently, MDR K. pneumoniae strains have recently spread and increased so that they have emerged as a cause of recalcitrant infections worldwide [49][50][51]. Bacteriophage therapy is therefore an effective alternative to multiple drug-resistant pathogens [52]. KPP-5 phage showed clear plaques, and this can be attributed to the production of soluble, polysaccharide-degrading enzymes [53]. Interestingly, the KPP-5 phage showed a broad host range where it could infect all the K. pneumoniae strains tested, although there was no lytic effect against the unrelated genes. On the contrary, the KP1513 and KP-34 K. pneumoniae phages showed a narrow host range [54,55]. The effect of the KPP-5 phage depends on the presence of specific b receptors located in the bacterial cell and other factors that determine the ability of the bacteriophages to multiply within their hosts [56]. An explanation of why the KPP-5 phage was able to infect all K. pneumoniae strains selectively in terms of the host range is possibly due to infecting the capsular type of K. pneumoniae [57]. This reinforces the suggestion that capsule depolymerase enzymes are contained in the KPP-5 phage, which enables the identification and digestion of particular capsular types. Previous studies have found that most K. pneumoniae phages encode depolymerases that can digest the Klebsiella capsule [58]. Polysaccharide depolymerase degrades the macromolecular carbohydrate that industrializes the capsule enclosing the bacterial cell wall [54], followed by peptidoglycan hydrolases, which break down the layer of peptidoglycan to penetrate the cell wall and enter the cytoplasm to enable the phage to transfer its genetic material [54].

One-Step Growth Curve and Thermal and pH Stability of KPP-5 Phage
The one-step growth curve exhibited that the KPP-5 phage had a relatively short latent period of 25 min, and the burst size was about 236 PFU/infected cells (Figure 2A). The short latent period and the high burst size are most likely due to the high processivity of phage DNA polymerase and, therefore, the high replication efficiency [60,61]. According to our findings, the latency period of the KPP-5 phage (25 min) was longer than that of the Klebsiella phage ϕBO1E (Podoviridae; 10 min) [53]. Fortunately, the burst size of KPP-5 (236 phage particles/infected cells) was higher than that of other Klebsiella phages from the Siphoviridae [62] and Podoviridae [53] families. In comparison to a previous study by Tabassum et al. [63], the latent period for the TSK1 phage was 30 min, with a moderate burst size of 113 PFU/infected cells, indicating that the KPP-5 phage is more effective. Our results indicate that the replication properties of the KPP-5 phage are similar to that of phage B5055 [64]. The bacterial growth reduction tendency of the KPP-5 phage was measured over 24 h of incubation, and the log growth reduction was from the start to six hours of incubation. Then, the bacterial growth remained lower than the undetectable range over the remainder of the experiment time ( Figure 2B). This result clearly indicated that there was no development of bacteriophage-resistant mutants due to lack of resistance in bacteria against the phage. Our explanation for the lack of resistance to bacteria may be that the bacteria have not behaved violently, like the formation of proteins that mask the recognition of phage receptor sites, the digestion of the phage genome, or inhibit the action of phage enzymes. It may have other interpretations that require further investigation [65].
On the other hand, the KPP-5 phage showed high titer stability at different temperatures and pH levels, where the phage titer showed high stability after incubation at 4, 25, 37, and 50 °C for all periods of time while phage titer reduction (2.38 log10 PFU/mL) was observed at 60 °C after incubation for 24 h. However, more than half of the phage titer (3.76 log10 PFU/mL) was reduced for one hour of incubation at 75 °C, and no stable titer was observed at 8, 12, and 24 h of incubation ( Figure 2C). Since temperature plays a key role in the bacteriophage survival capacity of the attachment and the length of the latent period [66], a previous study in line with our results showed that the Caudovirales bacteriophages were thermostable [67]. Similarly, the acidity and alkalinity of the environment are other important factors affecting the phage stability. The KPP-5 phage showed the highest titer stability at pH 7, while a phage titer reduction was detected at pH 5 and pH 9. Likewise, at pH 11, the phage titer showed a severe reduction of approximately 4.45   The TEM image showed that the KPP-5 phage had an icosahedral head (36.25 nm in diameter) with a short tail (12.59 nm in length), so that the KPP-5 phage had a podovirus morphology ( Figure 1B). In line with our findings, a recent study by Tan et al. [59] showed that two phages (phage 117 and phage 31) with short tails and icosahedral heads were isolated against all MDR K. pneumoniae.

One-Step Growth Curve and Thermal and pH Stability of KPP-5 Phage
The one-step growth curve exhibited that the KPP-5 phage had a relatively short latent period of 25 min, and the burst size was about 236 PFU/infected cells (Figure 2A). The short latent period and the high burst size are most likely due to the high processivity of phage DNA polymerase and, therefore, the high replication efficiency [60,61]. According to our findings, the latency period of the KPP-5 phage (25 min) was longer than that of the Klebsiella phage φBO1E (Podoviridae; 10 min) [53]. Fortunately, the burst size of KPP-5 (236 phage particles/infected cells) was higher than that of other Klebsiella phages from the Siphoviridae [62] and Podoviridae [53] families. In comparison to a previous study by Tabassum et al. [63], the latent period for the TSK1 phage was 30 min, with a moderate burst size of 113 PFU/infected cells, indicating that the KPP-5 phage is more effective.
Our results indicate that the replication properties of the KPP-5 phage are similar to that of phage B5055 [64]. The bacterial growth reduction tendency of the KPP-5 phage was measured over 24 h of incubation, and the log growth reduction was from the start to six hours of incubation. Then, the bacterial growth remained lower than the undetectable range over the remainder of the experiment time ( Figure 2B). This result clearly indicated that there was no development of bacteriophage-resistant mutants due to lack of resistance in bacteria against the phage. Our explanation for the lack of resistance to bacteria may be that the bacteria have not behaved violently, like the formation of proteins that mask the recognition of phage receptor sites, the digestion of the phage genome, or inhibit the action of phage enzymes. It may have other interpretations that require further investigation [65].
On the other hand, the KPP-5 phage showed high titer stability at different temperatures and pH levels, where the phage titer showed high stability after incubation at 4, 25, 37, and 50 • C for all periods of time while phage titer reduction (2.38 log 10 PFU/mL) was observed at 60 • C after incubation for 24 h. However, more than half of the phage titer (3.76 log 10 PFU/mL) was reduced for one hour of incubation at 75 • C, and no stable titer was observed at 8, 12, and 24 h of incubation ( Figure 2C). Since temperature plays a key role in the bacteriophage survival capacity of the attachment and the length of the latent period [66], a previous study in line with our results showed that the Caudovirales bacteriophages were thermostable [67]. Similarly, the acidity and alkalinity of the environment are other important factors affecting the phage stability. The KPP-5 phage showed the highest titer stability at pH 7, while a phage titer reduction was detected at pH 5 and pH 9. Likewise, at pH 11, the phage titer showed a severe reduction of approximately 4.45 and 5.68 (log 10 PFU/mL) after 8 and 12 h of incubation, while the complete loss of a phage titer was observed at pH 4 and pH 11 after incubation periods of 4 and 24 h, respectively ( Figure 2D). Our study was consistent with Jamalludeen et al. [36], who reported that most phages are able to survive well over a wide pH range (5-9) under physiological conditions that maintain a normal virion structure and stability. Extreme pH conditions affect bacteriophage activity by irreversible precipitation, coagulation, or shaking of the phages [68]. and 5.68 (log10 PFU/mL) after 8 and 12 h of incubation, while the complete loss of a phage titer was observed at pH 4 and pH 11 after incubation periods of 4 and 24 h, respectively ( Figure 2D). Our study was consistent with Jamalludeen et al. [36], who reported that most phages are able to survive well over a wide pH range (5-9) under physiological conditions that maintain a normal virion structure and stability. Extreme pH conditions affect bacteriophage activity by irreversible precipitation, coagulation, or shaking of the phages [68].

Genomic Analysis and Annotation of the KPP-5 Phage
In general, a bioinformatics analysis can predict both the biological properties and the safety of phage to be used for medical purposes. As a result, the complete KPP-5 phage genome was sequenced and analyzed using various bioinformatics tools. The KPP-5 phage has a linear dsDNA genome of 38,245 bp in length with a GC content of 50.8%. The genome contains 40 predicted open reading frames (ORFs), with the primary starting codon for all ORFs being ATG, all of which have been shown to be located on the KPP-5 phage genome′s positive strand ( Figure 3). Of the 40 ORFs predicted, only 31 are identical to the actual functional genes that have been shown to encode proteins, while the remaining nine ORFs were predicted to encode hypothetical proteins with no assigned function. No tRNA genes have been identified in the KPP-5 phage genome. The characteristics of the KPP-5 phage gene products, including ORF positions, predicted products (aa), and homologs with predicted protein motifs, are shown in Table S2. The total annotated ORFs can be further categorized into five modules (Table S3): (i) replication, regulation, transcription, and translation (14 ORFs); (ii) host cell lysis (3 ORFs); (iii) phage structure (12 ORFs); (iv) DNA packaging (2 ORFs); and (v) unknown functions (9 ORFs). Interestingly, phage-borne virulence genes or drug resistance genes were not present, suggesting that KPP-5 could be safely used as a phage biocontrol agent where phage genomes used for therapy should be absolutely safe [69].

Genomic Analysis and Annotation of the KPP-5 Phage
In general, a bioinformatics analysis can predict both the biological properties and the safety of phage to be used for medical purposes. As a result, the complete KPP-5 phage genome was sequenced and analyzed using various bioinformatics tools. The KPP-5 phage has a linear dsDNA genome of 38,245 bp in length with a GC content of 50.8%. The genome contains 40 predicted open reading frames (ORFs), with the primary starting codon for all ORFs being ATG, all of which have been shown to be located on the KPP-5 phage genome s positive strand ( Figure 3). Of the 40 ORFs predicted, only 31 are identical to the actual functional genes that have been shown to encode proteins, while the remaining nine ORFs were predicted to encode hypothetical proteins with no assigned function. No tRNA genes have been identified in the KPP-5 phage genome. The characteristics of the KPP-5 phage gene products, including ORF positions, predicted products (aa), and homologs with predicted protein motifs, are shown in Table S2. The total annotated ORFs can be further categorized into five modules (Table S3): (i) replication, regulation, transcription, and translation (14 ORFs); (ii) host cell lysis (3 ORFs); (iii) phage structure (12 ORFs); (iv) DNA packaging (2 ORFs); and (v) unknown functions (9 ORFs). Interestingly, phage-borne virulence genes or drug resistance genes were not present, suggesting that KPP-5 could be safely used as a phage biocontrol agent where phage genomes used for therapy should be absolutely safe [69]. Biomedicines 2021, 9, x FOR PEER REVIEW 9 of 21

Comparative Genomic Analysis
For a comparative genomic analysis, the genome of the KPP-5 phage was searched against the phage genome sequences in the GenBank using BLASTN in NCBI, where 33 phages belonging to the Autographiviridae family and four genera (Teetrevirus (n = 18), Teseptimavirus (n = 5), Berlinvirus (n = 5), and Przondovirus (n = 5)) were used for analysis (Table S4). The comparative genomic analysis identified KPP-5 phage as belonging to the Teetrevirus genus with a high nucleotide sequence similar to the phages in this genus (identity, 85.74%-96.41% and query coverage, 79-99% with E value 0) (Table S4). Of these 18 phages, Klebsiella phage vB_KpnP_Emp27 (GenBank accession no. MN013074), Klebsiella phage Patroon (GenBank accession no. MK608335), and Klebsiella phage NL_ZS_3 (GenBank accession no. MT813142) are closely related to KPP-5 with 99%, 95%, and 94% query coverage, respectively, and 95.19%, 94.85%, and 96.41% complete genome nucleotide sequence identity, respectively ( Figure 4A-C). On the other hand, Klebsiella phage 31 (GenBank accession no. MN149904) was similar to KPP-5 with 93% query coverage and 85.74% complete genome nucleotide sequence identity ( Figure 4D). The average nucleotide identity (ANI) values were calculated using the reciprocal best hits (two-way ANI) between two genomic datasets, as calculated by Goris et al. [70] to further evaluate the relationship between KPP-5 and phage vB_KpnP_Emp27, phage Patroon, phage NL_ZS_3, and phage 31. Since the ANI values were 94.85%, 95.20%, 95.89%, and 95.52%, respectively, this shows that KPP-5 phage is a member of the new species. A phylogenetic analysis of the whole KPP-5 genome using the VICTOR web service revealed five clusters of four genera (Teetrevirus, Teseptimavirus, Berlinvirus, and Przondovirus) within the Autographviridae family and outgroup within the Siphoviridae family (Enterobacteria phage DE3). KPP-5 phage is most closely associated with the Teetrevirus genus in the Autographviridae family, where it was found to be the closest relative to Klebsiella phage vB_KpnP_Emp27 (GenBank accession no. MN013074) ( Figure 5). It is therefore suggested

Comparative Genomic Analysis
For a comparative genomic analysis, the genome of the KPP-5 phage was searched against the phage genome sequences in the GenBank using BLASTN in NCBI, where 33 phages belonging to the Autographiviridae family and four genera (Teetrevirus (n = 18), Teseptimavirus (n = 5), Berlinvirus (n = 5), and Przondovirus (n = 5)) were used for analysis (Table S4). The comparative genomic analysis identified KPP-5 phage as belonging to the Teetrevirus genus with a high nucleotide sequence similar to the phages in this genus (identity, 85.74%-96.41% and query coverage, 79-99% with E value 0) (Table S4). Of these 18 phages, Klebsiella phage vB_KpnP_Emp27 (GenBank accession no. MN013074), Klebsiella phage Patroon (GenBank accession no. MK608335), and Klebsiella phage NL_ZS_3 (GenBank accession no. MT813142) are closely related to KPP-5 with 99%, 95%, and 94% query coverage, respectively, and 95.19%, 94.85%, and 96.41% complete genome nucleotide sequence identity, respectively ( Figure 4A-C). On the other hand, Klebsiella phage 31 (GenBank accession no. MN149904) was similar to KPP-5 with 93% query coverage and 85.74% complete genome nucleotide sequence identity ( Figure 4D). The average nucleotide identity (ANI) values were calculated using the reciprocal best hits (two-way ANI) between two genomic datasets, as calculated by Goris et al. [70] to further evaluate the relationship between KPP-5 and phage vB_KpnP_Emp27, phage Patroon, phage NL_ZS_3, and phage 31. Since the ANI values were 94.85%, 95.20%, 95.89%, and 95.52%, respectively, this shows that KPP-5 phage is a member of the new species. A phylogenetic analysis of the whole KPP-5 genome using the VICTOR web service revealed five clusters of four genera (Teetrevirus, Teseptimavirus, Berlinvirus, and Przondovirus) within the Autographviridae family and outgroup within the Siphoviridae family (Enterobacteria phage DE3). KPP-5 phage is most closely associated with the Teetrevirus genus in the Autographviridae family, where it was found to be the closest relative to Klebsiella phage vB_KpnP_Emp27 (GenBank accession no. MN013074) ( Figure 5). It is therefore suggested that the KPP-5 phage be classified into the Autographviridae family and the Teetrevirus genus. lation; 3 predicted to encode proteins for host cell lysis; 12 predicted to encode proteins for the phage structure; 2 predicted to encode proteins for DNA packaging; and 6 predicted to encode proteins for unknown functions (hypothetical proteins) (Tables S3 and  S5). Three distinct genes in the KPP-5 phage genome were shown in the comparison between the KPP-5 and vB-KpnP-Emp27 phage genomes, including gp2 (GenBank No. QSJ04722), gp9 (GenBank No. QSJ04729), and gp35 (GenBank No. QSJ04755), which were predicted to encode hypothetical proteins (Tables S3 and S5).  The core genes of KPP-5, vB_KpnP_Emp27, Patroon, NL_ZS_3, and 31 phages were examined, and 34 genes were conserved in their genomes (Table S5), while 36 genes were found to be conserved in the genomes of four high-identity phages, including phage KPP-5, vB_KpnP_Emp27, Patroon, and NL_ZS_3 (Table S5). Interestingly, 37 genes were found to be conserved in the genomes of the two KPP-5 and vB_KpnP_Emp27 phages that exhibited high identity and query coverage ( Table S). The gene products of these 37 genes were 14 predicted to encode proteins for replication, regulation, transcription, and translation; 3 predicted to encode proteins for host cell lysis; 12 predicted to encode proteins for the phage structure; 2 predicted to encode proteins for DNA packaging; and 6 predicted to encode proteins for unknown functions (hypothetical proteins) (Tables S3 and S5). Three distinct genes in the KPP-5 phage genome were shown in the comparison between the KPP-5 and vB-KpnP-Emp27 phage genomes, including gp2 (GenBank No. QSJ04722), gp9 (GenBank No. QSJ04729), and gp35 (GenBank No. QSJ04755), which were predicted to encode hypothetical proteins (Tables S3 and S5).
DNA polymerase is a conserved protein frequently used in the analysis of the diversity and global spread of podophages [61]. Therefore, to observe their relationships in greater detail, we compared KPP-5 phage DNA polymerase with the same phages used in the entire genome phylogeny (Table S4) Table S5 shows the similarity of 14 genes of Klebsiella phage KPP-5 with genes   (Table S3).
DNA polymerase is a conserved protein frequently used in the analysis of the diversity and global spread of podophages [61]. Therefore, to observe their relationships in greater detail, we compared KPP-5 phage DNA polymerase with the same phages used in the entire genome phylogeny (Table S4). A DNA polymerase phylogenetic analysis showed that KPP-5 is most closely related to the Teetrevirus genus in the Autographviridae family ( Figure 6)

Host Cell Lysis Related Genes
Three genes were identified as involving enzymatic cleavage, one encoding endolysin (gp7, GenBank No. QSJ04727), the second encoding holin (gp30, GenBank No. QSJ04750), and the third encoding endopeptidase (gp32, GenBank No. QSJ04752), typically a three-component lysis system for host cell phage release out of the 40 genes predicted in the KPP-5 phage genome. KPP-5 phage endolysin showed a high similarity to     [71]. Four classes are classified as endolysins: muramidase, endopeptidases, amidase, and lytic transglycosylases [71,72]. The endolysin of phage KPP-5 was found to be 16.88 KDa (Table S2), where most Gram-negative endolysininfecting bacteria typically represent globular proteins of 15-20-kDa single domains [71,73]. A T7 phage endolysin (N-acetylmuramoyl-L-alanine amidase) is a 17-kDa protein that, by hydrolyzing the amide bond between the L-alanine and N-acetylmuramoyl residues of the peptidoglycan layer, lyses a variety of Gram-negative [71]. Holins regulate the access of phage-encoded endolysins to peptidoglycan by accumulating and developing lesions in the cytoplasmic membrane, thus triggering host cell lysis at a specific time point [71,72,74]. It was noted that KPP-5 phage holin was classified as class II and contained a single hydrophobic transmembrane domain (TMD) region via the prediction of the TMHMM server (Figure 7). Depending on their topology, holins are divided into three categories (classes I-III), in which all holins have at least one TMD [75].  (Tables S2 and S5). Endolysins are distinct peptidoglycan hydrolases and, in various infectious models, have tremendous potential as effective enzybiotics [71]. Four classes are classified as endolysins: muramidase, endopeptidases, amidase, and lytic transglycosylases [71,72]. The endolysin of phage KPP-5 was found to be 16.88 KDa (Table  S2), where most Gram-negative endolysin-infecting bacteria typically represent globular proteins of 15-20-kDa single domains [71,73]. A T7 phage endolysin (N-acetylmuramoyl-L-alanine amidase) is a 17-kDa protein that, by hydrolyzing the amide bond between the L-alanine and N-acetylmuramoyl residues of the peptidoglycan layer, lyses a variety of Gram-negative [71]. Holins regulate the access of phage-encoded endolysins to peptidoglycan by accumulating and developing lesions in the cytoplasmic membrane, thus triggering host cell lysis at a specific time point [71,72,74]. It was noted that KPP-5 phage holin was classified as class II and contained a single hydrophobic transmembrane domain (TMD) region via the prediction of the TMHMM server (Figure 7). Depending on their topology, holins are divided into three categories (classes I-III), in which all holins have at least one TMD [75].

Phage Structure Related Genes
Six protein-encoding genes associated with the structural tail, including tail assembly protein (gp17, GenBank No. QSJ04739), head-to-tail joining protein (gp18, GenBank No. QSJ04737), major tail protein (gp21, GenBank No. QSJ04741), tail tubular protein A (gp22, GenBank No. QSJ04742), tail tubular protein B (gp23, GenBank No. QSJ04743), and tail fiber protein (gp28, GenBank No. QSJ04748), were annotated among the 12 actual phage structure genes of phage KPP-5 (Table S3). Most likely, they are involved in tail assembly or phage penetration during infection through the host cell s external membrane [76]. KPP-5 phage tail fiber protein showed high similarity to Klebsiella phage vB_KpnP_Emp27 (94.38% identity, and 98% query coverage, GenBank No. QEG11854) compared to this in Klebsiella phage NL_ZS_3 (70.49% identity, and 51% query coverage, GenBank No. QNN97333), Klebsiella phage Patroon (51.89% identity, and 89% query coverage, GenBank No. QBQ72909), and Klebsiella phage 31 (51.89% identity, and 89% query coverage, GenBank No. QGH73753) (Table S5). Interestingly, there was a match (E-value 9.24e-140) to the superfamily of phage T7-like tail fiber protein (PHA00430) and to the superfamily of phage T7 tail fiber protein with E-value 2.36e-63 (pfam03906) [77,78]. Further analysis of KPP-5 phage tail fiber protein using the Phyre2 server [79] revealed that it exhibited structural similarity to the phage T7 tail fiber protein (confidence, 99.1% and identity, 34%) [80]. The tail fiber protein can accurately recognize and bind to the host surface receptors and can mutate throughout evolution, contributing to a shift in the phage host s range [81]. Interestingly, KPP-5 phage tail tubular protein A (TTPA) showed a high similarity to Yersinia phage phiYeO3-12 (100% identity and 100% query coverage, GenBank No. NP_052110). In view of the hydrolytic activity of the Yersinia phage phiYeO3-12 TTPA towards Red starch, it has been shown that this enzyme could be classified within the alpha-1, 4-glucosidase family [82]. Further analysis of KPP-5 phage TTPA using the Phyre2 server [79] revealed that it was structurally similar to the T7 gatekeeper protein (gp11) structure (confidence, 100% and identity, 80%) [83] and to klebsiella pneumoniae phage kp32 TTPA, which has lytic activity towards capsular exopolysaccharide (EPS) (confidence, 100% and identity, 61%) [84], since the use of phages is promising for the destruction of bacterial biofilm via depolymerases [84][85][86]. Depolymerases may be found in the tail of the phage or may exist as an extracellular enzyme [84]. Klebsiella pneumoniae phage kp32 TTPA exhibited lytic activity towards capsular exopolysaccharide of the multi-resistant clinical strain of Klebsiella pneumoniae, PCM2713 where the enzymatic activity of TTPA may reflect the presence of a peptidoglycan hydrolase domain in the α-helical region [84]. Furthermore, protein folding and aggregation tests have shown that Yersinia phage phiYeO3-12 TTPAgp11 is a single-domain protein that may be aggregated by maltose or N-acetylglucosamine [82].

DNA Packaging-Related Genes
Two genes, including DNA packaging protein A (terminase small subunit, gp31, Gen-Bank No. QSJ04751) and DNA packaging protein B (terminase large subunit, gp33, Gen-Bank No. QSJ04753), were annotated as phage genome packaging-associated genes in the KPP-5 phage genome (Table S3). During the phage replication process, the dsDNA needs to be packaged, and this requires the ATP-guided terminase proteins responsible for slicing dsDNA into the final genome size sequences, followed by placing them into an empty capsid [91,92]. The phage KPP-5 terminase small subunit is highly similar to the terminase small subunit of Klebsiella phage Patroon and Klebsiella phage 31 (100% identity and 100% query coverage) (Table S5), while the phage KPP-5 terminase large subunit is highly similar to the terminase large subunit of Klebsiella phage vB_KpnP_Emp27 and Klebsiella phage NL_ZS_3 (100% identity, 100% query coverage) ( Table S5). The terminase large subunit is often used to classify phages, because it is generally well-conserved in tailed

DNA Packaging-Related Genes
Two genes, including DNA packaging protein A (terminase small subunit, gp31, GenBank No. QSJ04751) and DNA packaging protein B (terminase large subunit, gp33, GenBank No. QSJ04753), were annotated as phage genome packaging-associated genes in the KPP-5 phage genome (Table S3). During the phage replication process, the dsDNA needs to be packaged, and this requires the ATP-guided terminase proteins responsible for slicing dsDNA into the final genome size sequences, followed by placing them into an empty capsid [91,92]. The phage KPP-5 terminase small subunit is highly similar to the terminase small subunit of Klebsiella phage Patroon and Klebsiella phage 31 (100% identity and 100% query coverage) (Table S5), while the phage KPP-5 terminase large subunit is highly similar to the terminase large subunit of Klebsiella phage vB_KpnP_Emp27 and Klebsiella phage NL_ZS_3 (100% identity, 100% query coverage) ( Table S5). The terminase large subunit is often used to classify phages, because it is generally well-conserved in tailed phages where phages with a similar gene would assume that their terminase large subunit would be clustered together [93]. Therefore, on the basis of the alignment of the terminase large subunit amino acid sequences, the phylogenetic relationship between Klebsiella phage KPP-5 and the same phages used in the entire genome phylogeny (Table S4), the terminase large subunit phylogenetic analysis showed that KPP-5 is most closely related to the Teetrevirus genus in the Autographviridae family ( Figure 9). edicines 2021, 9, x FOR PEER REVIEW 16 of 21 phages where phages with a similar gene would assume that their terminase large subunit would be clustered together [93]. Therefore, on the basis of the alignment of the terminase large subunit amino acid sequences, the phylogenetic relationship between Klebsiella phage KPP-5 and the same phages used in the entire genome phylogeny (Table S4), the terminase large subunit phylogenetic analysis showed that KPP-5 is most closely related to the Teetrevirus genus in the Autographviridae family ( Figure 9).

Related Unknown Function Genes
It was predicted that nine genes (gp2, gp3, gp4, gp9, gp10, gp16, gp35, gp36, and gp39 with GenBank Nos. QSJ04722, QSJ04723, QSJ04724, QSJ04729, QSJ04730, QSJ04736, QSJ04755, QSJ04756, and QSJ04759, respectively) would encode hypothetical unassigned function proteins among the 40 Klebsiella phage KPP-5 genes (Table S3). Of the nine genes predicted, only seven were homologous to other phage-encoded uncharacterized proteins, and two (gp2 and gp35) were detected only in the KPP-5 phage genome (Table S2). Furthermore, six out of seven genes were homologous to the hypothetical proteins of Klebsiella phages, where five genes (gp3, gp4, gp10, gp16, and gp39) were homologous to Klebsiella phage vB_KpnP_Emp27, Klebsiella phage Patroon, Klebsiella phage NL_ZS_3, and Klebsiella phage 31, while gp36 was homologous to Klebsiella phage vB_KpnP_Emp27, Klebsiella phage Patroon, and Klebsiella phage NL_ZS_3 (Table S5). Interestingly, gp9 was highly similar to the hypothetical protein of Citrobacter phage SH1 (97.25% identity and 100% query coverage, GenBank No. YP_009286649) contained double TMD regions via the prediction of the TMHMM server ( Figure 10). In phage genomes, 70% of the predicted genes are annotated to encode hypothetical proteins, most of which range from 30 to 200 amino acids with unknown functions [18,[94][95][96]. The high potential for the use of these proteins as toxic proteins has already been demonstrated, where they can be toxic through unknown mechanisms and have the potential to develop into novel antimicrobial peptides [18].  (Table S5). Interestingly, gp9 was highly similar to the hypothetical protein of Citrobacter phage SH1 (97.25% identity and 100% query coverage, GenBank No. YP_009286649) contained double TMD regions via the prediction of the TMHMM server ( Figure 10). In phage genomes, 70% of the predicted genes are annotated to encode hypothetical proteins, most of which range from 30 to 200 amino acids with unknown functions [18,[94][95][96]. The high potential for the use of these proteins as toxic proteins has already been demonstrated, where they can be toxic through unknown mechanisms and have the potential to develop into novel antimicrobial peptides [18].

Conclusions
In conclusion, we identified a new phage named KPP-5 capable of lysing MDR K. pneumoniae isolated from food samples in Egypt. Based on our findings, the KPP-5 phage is a promising candidate for the control of Klebsiella strains and could be safely used as a biocontrol agent for phage-borne virulence genes or for drug resistance genes in their genomes.
Supplementary Materials: The following are available online at www.mdpi.com/xxx/s1: Table S1: Antibiotics resistance profiles of the K. pneumoniae host strain. Table S2: Genome feature of phage KPP-5, gene products, and their predicted functions. Table S3: KPP-5 phage gene products are involved in five different modules. Table S4: Comparison of phages with BLASTN scores against the KPP-5 phage genome.

Conclusions
In conclusion, we identified a new phage named KPP-5 capable of lysing MDR K. pneumoniae isolated from food samples in Egypt. Based on our findings, the KPP-5 phage is a promising candidate for the control of Klebsiella strains and could be safely used as a biocontrol agent for phage-borne virulence genes or for drug resistance genes in their genomes.