Next Article in Journal
Forced Swimming-Induced Depressive-like Behavior and Anxiety Are Reduced by Chlorpheniramine via Suppression of Oxidative and Inflammatory Mediators and Activating the Nrf2-BDNF Signaling Pathway
Next Article in Special Issue
Recent Advances in Molecular Mechanisms of Nucleoside Antivirals
Previous Article in Journal
Antidiabetic Effect of Fermented Mesembryanthemum crystallinum L. in db/db Mice Involves Regulation of PI3K-Akt Pathway
Previous Article in Special Issue
Plant Antimicrobial Peptides: Insights into Structure-Function Relationships for Practical Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Phenotypic Characterization and Comparative Genomic Analyses of Mycobacteriophage WIVsmall as A New Member Assigned to F1 Subcluster

Institute of Biomedical Informatics, School of Basic Medical Sciences, Henan University, Kaifeng 475004, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Curr. Issues Mol. Biol. 2023, 45(8), 6432-6448; https://doi.org/10.3390/cimb45080406
Submission received: 28 June 2023 / Revised: 28 July 2023 / Accepted: 31 July 2023 / Published: 3 August 2023
(This article belongs to the Special Issue Advanced Research in Antimicrobial and Antiviral Drugs)

Abstract

:
In this study, we conducted the morphological observation, biological and genomic characterization, evolutionary analysis, comparative genomics description, and proteome identification of a recently isolated mycobacteriophage, WIVsmall. Morphologically, WIVsmall is classified as a member of the Siphoviridae family, characterized by a flexible tail, measuring approximately 212 nm in length. The double-stranded phage genome DNA of WIVsmall spans 53,359 base pairs, and exhibits a G + C content of 61.01%. The genome of WIVsmall comprises 103 protein-coding genes, while no tRNA genes were detected. The genome annotation unveiled the presence of functional gene clusters responsible for mycobacteriophage assembly and maturation, replication, cell lysis, and functional protein synthesis. Based on the analysis of the phylogenetic tree, the genome of WIVsmall was classified as belonging to subgroup F1. A comparative genomics analysis indicated that the WIVsmall genome exhibited the highest similarity to the phage SG4, with a percentage of 64%. The single-step growth curve analysis of WIVsmall revealed a latent period of 120 min, and an outbreak period of 200 min.

1. Introduction

As a crucial milestone in modern medicine, antibiotics have played a pivotal role in controlling pathogenic infections, and saving innumerable lives, over the past century. However, as bacteria have evolved multiple mechanisms to impede the efficacy of antimicrobial agents, it is imperative that we discover innovative approaches to combating drug-resistant bacterial infections [1]. In fact, there are now numerous alternative therapies available to combat drug resistance, including monoclonal antibodies, and microbiota therapy [2]. However, bacteriophages have emerged as a highly promising alternative to antibiotics (Figure 1A). Compared with antibiotics, bacteriophages offer several significant advantages against bacterial infections, including host specificity, the absence of cross immunogenicity with antibiotics [3], abundance in nature [4], and self-replication [5], improving the speed and quality of treatment. Thus far, bacteriophages have proven effective in treating various clinical cases of antibiotic-resistant bacterial infections. For example, a 63-year-old patient diagnosed with a urinary tract infection caused by K. pneumoniae was successfully treated through a phage cocktail therapy, comprising six phages in combination with anti-infective drugs, resulting in clinical cure [6]. Additionally, bacteriophage therapy has demonstrated efficacy in treating infections caused by S. aureus, E. coli, A. bambini, and M. abscessus [7,8]. Although the clinical application of mycobacteriophages for treating M. abscessus infections has been successful, their widespread implementation is challenging, due to their limited efficacy against other clinically isolated M. abscessus strains. The clinical application of phage therapy may be perceived as a form of personalized therapy, which may limit its accessibility to those with the financial means.
The World Health Organization (WHO) released a global priority list of pathogens in 2017 [9]. The pathogens were categorized as critical-, high-, and medium-priority bacteria, based on the urgency of developing new antibiotics to combat them. Mycobacterium tuberculosis has been identified as one of the high-priority bacteria. M. tuberculosis is the causative agent of pulmonary tuberculosis (TB), a disease that poses a significant threat to human health. It is estimated that approximately one quarter of the world’s population is infected with M. tuberculosis, with over 1.4 million deaths attributed to tuberculosis each year [10]. The emergence of drug-resistant strains of M. tuberculosis poses a significant obstacle to the effective control of tuberculosis [11]. Based on the 2022 World Health Organization (WHO) report on tuberculosis (TB), the proportion of multi-drug resistant cases was estimated as being 20% in 2021, and almost 43% of MDR-TB cases were treated between 2018 and 2021 [12]. Mycobacteriophages, either in isolation or as a cocktail, or as part of a synergistic therapy with antibiotics, have potential in the treatment of mycobacterial infections [13]. It is reported that 20 patients with antibiotic-refractory mycobacterial infections received treatment with a mycobacteriophage cocktail [14]. Favorable clinical or microbiological responses were observed in 11 of these patients. Additionally, Jessica S. Little and her team reported a case of a 56-year-old patient with a refractory disseminated cutaneous Mycobacterium chelonae infection, who received mycobacteriophage therapy, and experienced a significant improvement in their skin lesions [15]. Mycobacteriophages can also be utilized as a tool for engineering shuttle plasmids to detect drug-resistant M. tuberculosis strains, including DS6A and TM4-based shuttle plasmids [13,16].

2. Materials and Methods

2.1. Phage Isolation and Preparation

M. smegmatis mc2155 served as the host bacterium for the isolation of the mycobacteriophages. The isolation of WIVsmall was carried out in accordance with a previously described methodology [17]. For the phage isolation, 10 g of soil sample was mixed with 10–15 mL of buffer (50 mM Tris-HCl, pH 7.5, 150 mM NaCl, 10 mM MgSO4, and 2 mM CaCl2), and 1 mL of M. smegmatis mc2155, in a 50 mL Erlenmeyer flask. The mixture was incubated at 37 °C for 24 h. The supernatant was filtered using a 0.22-μm membrane filter, to eliminate residual bacteria and phytoplankton. The conventional double-layer agar technique and spot assay were employed to isolate individual phages. The process of single plaque-picking was repeated thrice. The phage lysates were subsequently stored at 4 °C for further experimentation

2.2. DNA Extraction

Genomic DNA was extracted from the purified phage using SDS-proteinase K protocols [18]. The structural protein of the phage WIVsmall was digested, using a combined reagent consisting of 50 mg of protease K (Invitrogen, Shanghai, China), 20 mM EDTA, and 1% SDS, for 4 h at 56 °C. The phage DNA was extracted using the conventional phenol–chloroform method, with an equal volume of phenol/chloroform/isoamyl alcohol (25:24:1, v/v). The aforementioned step was repeated, and an equal volume of chloroform was employed to further refine the supernatant. The combined reagent, consisting of one volume of 2-propanol, and 0.4 volumes of 3M sodium acetate (pH 4.6), was added to the supernatant. The resulting mixture was then centrifuged at 7000× g for 10 min at 4 °C, to precipitate the phage DNA.

2.3. Sequence Data and Phylogenetic Analysis

The complete genome was sequenced, using 454 technology, with the GS Junior 454 system platform (Roche Diagnostics, Indianapolis, IN, USA), yielding 145-fold coverage of the phage genome. The 15 available nucleotide sequences of the mycobacteriophage genomes in the F cluster were retrieved from the NCBI nucleotide database. Based on the phage coding sequence, a phylogenetic tree was drawn, using the neighbor linkage (NJ) method. The bootstrap value was 1000 replicates on MEGA X [19].

2.4. Optimal Multiplicity of Infection (MOI)

M. smegmatis mc2155 was cultured to the exponential phase, and subsequently washed thrice with 7H9 medium, to eliminate Tween-80. The bacterial pellets were resuspended in 7H9 medium, and the cell concentration was adjusted to 106 colony-forming units per milliliter (CFU/mL). The phages were mixed with M. smegmatis mc2155 at varying dilutions, including the ratios of 10:1, 1:1, 1:10, 1:100, and 1:1000. The aforementioned mixtures were incubated at 37 °C, with agitation for 12 h, at 177 revolutions per minute. Subsequently, the 1 mL mixture was centrifuged at 10,000× g for 20 min, to remove the precipitated bacteria. The resulting supernatant was then filtered through a 0.22 µm pore-size membrane filter. The phage filtrate was detected through the employment of the gradient dilution method, and a double-layer agar (DLA) assay, as previously described [20]. The dilution that generated the highest phage titer was considered as the optimal MOI.

2.5. Single-Step Growth Experiment

To determine the single-step growth curve, 2 mL of M. smegmatis mc2155, with a titer of 106 CFU/mL, was mixed with 2 mL of the WIVsmall phage solution, at the optimal MOI, and incubated at 37 °C for 25 min. The mixture was centrifuged at 12,500× g for 2 min, to remove unabsorbed phage particles. The precipitated M. smegmatis mc2155 was suspended in mycobacteriophage buffer, and incubated at 37 °C and 177 rpm. The supernatant was collected to detect the phage titer, using the DLA method, every 30 min. The experiment was repeated three times, and a single-step growth curve was plotted.

2.6. Stability of the Phage under Various Conditions

The stability of the phage under several conditions was detected, as previously described, with some modifications [21]. To determine the phage stability at different temperatures, phage suspensions with a titer of 109 PFU/mL were incubated at 4 °C, 10 °C, 20 °C, 30 °C, 40 °C, 50 °C, 60 °C, 70 °C, 80 °C, and 90 °C for 2 h. To evaluate the ultraviolet (UV) stability of the phage WIVsmall, phage preparations with a titer of 109 PFU/mL were illuminated using a UV lamp (365 nm, 18µW/cm2) for 20, 40, 60, 80, 100, and 120 min. Subsequently, the phage titers were determined, using the DLA method. Three parallel experiments were conducted. To assess the pH stability of the phage WIVsmall, 200 µL of phage lysates was diluted in 1.8 mL of MP buffers, at various pH values ranging from 3–12, and incubated at 37 °C for 2 h.

2.7. Genomic Bioinformatics Analysis

Softberry (http://linux1.softberry.com/berry.phtmltopic=virus0&group=programs&subgroup=gfindv, accessed on 20 May 2023) and GeneMarkTM (http://exon.gatech.edu/GeneMark/, accessed on 20 May 2023) were used to predict the ORFs of the WIVsmall genome [22]. The starting codons were selected as ATG, TTG, and GTG. Sequences with bases less than 150 bp, and ORFs without SD sequences were manually deleted. The active domains and isoelectric points of ORFs with protein functions were predicted using Pfam (http://pfam.sanger.ac.uk/, accessed on 20 May 2023) [23] and Isoelectric Point Finder (http://greengene.uml.edu/programs/Find MW.html, accessed on 20 May 2023). tRNAscan-SE (http://lowe lab.ucsc.edu/tRNAscan-SE/, accessed on 20 May 2023) [24] and Aragorn (http://mbio-serv2.Mbioekol.lu.se/ARAGORN/, accessed on 20 May 2023) [25] were used to search for the tRNA-encoding genes in the genome. BLASTn (http://blast.ncbi.nlm.nih.gov/Blast.cgi, accessed on 20 May 2023) and FASTA (http://www.ebi.ac.uk/Tools/fasta33/index.html, accessed on 20 May 2023) were used to perform the DNA sequence alignment. Multiple genome alignments were performed using the software easyfig. CGview (http://wishart.biology.ualberta.ca/cgview/, accessed on 20 May 2023) [26] and GenomeVx (http://wolfe.ucd.ie/GenomeVx/, accessed on 20 May 2023) were used in the visual analysis of the genomic maps. MEGA-X was used to perform neighbor-joining evolutionary tree analysis via 1000 bootstrap replications [19]. DNAPlotter was used to detect the classification of the bacteriophage clusters [27].

2.8. Identification of Phage Structural Proteins

Polyethylene glycol 8000 powder was added to the bacteriophage solution to achieve a final concentration of 20%, followed by overnight precipitation at 4 °C. The phage precipitate was collected via centrifugation at 12,000 rpm for 30 min, and then re-suspended in a dithiothreitol (DTT) solution. The lysate was placed on ice, and subjected to pulse ultrasound treatment. The ultrasonically treated samples were electrophoresed ona 12% SDS-PAGE, and subsequently stained with Coomassie Brilliant Blue R250. The protein bands were excised, digested with trypsin, and subsequently analyzed via mass spectrometry.

3. Results

3.1. Phage Isolation and Morphology Analysis of Plaque

The mycobacteriophage WIVsmall was isolated from soil samples collected in Henan Province, China, using M. smegmatis mc2155 as the host organism. The plaques of WIVsmall manifested as pinpricks accompanied by a slight turbidity on an M. smegmatis mc2155 lawn (Figure 2A), indicating a low frequency of lysogeny [28]. This result is consistent with the plaque formation characteristics of temperate mycobacteriophages [29]. The TEM analysis of purified phage particles revealed that the phage WIVsmall, a member of the Siphoviridae family, comprised an isometric head with a diameter of approximately 72 nm, and a long, flexible tail with a diameter of approximately 212 nm (Figure 2B,C).

3.2. Optimal Multiplicity of Infection (MOI) and One-Step Growth Curve

A familiarity with the growth characteristics of bacteriophages is a prerequisite to investigating their life cycle. We conducted an analysis of the optimal multiplicity of infection, and the one-step growth curve of the phage WIVsmall. To ascertain the optimal multiplicity of infection (MOI), the phage WIVsmall was incubated with its host, M. smegmatis mc2155, under varying MOI conditions. The MOI of 0.1 could yield the highest titer of the progeny phage (Figure 3B), and could be considered the optimal MOI. The single-step growth curve analysis exhibited that WIVsmall had a latent period of approximately 120 min, and a burst size of approximately 12.8 PFU per infected cell after 200 min of incubation at 37 °C (Figure 3A). Like other temperate phages, the extended latency period and limited burst size of WIVsmall can likely be attributed to the low activity of its DNA polymerase [30].

3.3. General Genome Analysis

The genomic sequence has been deposited in the NCBI database, under the accession number GenBank: KC736071.1. The genome size of WIVsmall was 53,359 bases, with a G + C content of 61.01%, which is comparable to that of its host mycobacteria, and other sequenced mycobacteriophages. The genome harbored 103 open reading frames (ORFs), spanning from 99 bp to 3639 bp, and no putative genes encoding tRNA or tmRNA were identified. WIVsmall demonstrated a condensed genome organization (Figure 4), and the 3′-terminal genome’s sticky end sequence was CGGACGGCGC. Based on sequence comparisons within the Actinobacteriophage Database (phaseDB.org), all the ORFs exhibited sequence homology with other entries in the database. Among them,48 ORFs have been functionally characterized, while the remaining 55 ORFs exhibit homology to proteins that have not yet been characterized. There are 82 genes exhibiting mosaic structures in the genome of WIVsmall, suggesting that these genes may have undergone horizontal gene transfer, which is comparable to other F1 cluster mycobacteriophages. Compared to the majority of the phages in the F1 cluster, the genome of WIVsmall exhibits the conspicuous absence of eight coding genes, three of which encode proteins with well-defined functions, namely carboxypeptidase, DNA methylase, and endonuclease. Among the three functional proteins, endonuclease was found to have multiple homologs with an identical function in the genome of the phage WIVsmall, while no homologs with an identical function were identified for the remaining two proteins. However, other F1 cluster phages, including Cerasum and BigPhil, with genome lengths comparable to WIVsmall, were also found to lack coding genes for carboxypeptidase and methylase. Table 1 presents the anticipated dimensions, location, transcriptional alignment, and the nearest phage protein analogue.

3.4. Putative Functions of the Predicted ORFs

The bioinformatics analysis revealed that the WIVsmall genome presented a functional mosaic structure. The WIVsmall genome can be partitioned into four distinct functional modules: phage assembly and maturation, replication, cell lysis, and functional proteins. The genes associated with phage assembly and maturation are located on the left arm of the genome, comprising five ORFs. ORF80 and ORF79 exhibit a significant sequence similarity with the large and small subunits of the terminal enzymes, respectively. These enzymes are accountable for the assembly of the phage DNA into the capsid. Three site-specific DNA endonucleases, encoded by ORF67, ORF48, and ORF23, are implicated in the horizontal transfer of phage genes.
The lysis module of WIVsmall comprises three contiguous open reading frames (ORF49, ORF50, and ORF51). The protein encoded by ORF50 exhibits a significant similarity to LysinA, and demonstrates catalytic activity toward the peptidoglycan layer of mycobacterial cell walls. ORF51 encodes the LysinB protein, which functions as a mycolylarabinogalactan hydrolase enzyme, capable of cleaving the linkage of the PG–AG polymer and, thereby, facilitating the detachment of the unique mycolic acid layer of the mycobacterial cell wall [31]. The deletion of the lysinB gene has been shown to decrease both the plaque and burst sizes [32]. ORF49 encodes the drilling protein holin, which facilitates the formation of pores in the cytoplasmic membrane, thereby enabling the release of LysinA and LysinB from the cytoplasm to the target cell wall [33].
The DNA replication module comprises ORFs encoding proteins involved in DNA replication and transcription. ORF7 encodes a DNA polymerase, whereas ORF38 and ORF25 exhibit significant sequence homology with transcriptional regulators belonging to the Xre family. ORF28 and ORF26 encode putative transcriptional regulatory proteins of the WhiB family, which recognize promoter regions, and regulate gene expression in the phage WIVsmall [34]. In general, primer enzymes and helicases are essential to the initiation of DNA replication. However, the absence of primases and helicases in the genome of the phage WIVsmall suggests its dependence on the host for genome replication, repair, and transcription. A putative integrase (ORF41) was also predicted, which can determine whether phages undergo lysogenic or lytic cycles.
The morphology module encompasses genes that encode the structural proteins of WIVsmall. Three adjacent ORFs, namely ORF77, ORF78, and ORF79, were identified as encoding the major capsid protein, head-scaffolding protein, and head-maturation protease, respectively. These proteins are believed to play crucial roles in stabilizing the condensed form of DNA, and facilitating head development. ORF63 and ORF44 encode large and small tail proteins, respectively. The ORFs involved in the tail structure formation comprise minor tail proteins, major tail proteins, and tail-assembly chaperones.ORF67 encodes a tape-measure protein that precisely measures the length of the bacteriophage tail.ORF78 encodes a portal protein capable of facilitating the transfer of the phage DNA into host cells through the formation of a channel.

3.5. Phylogenetic Relationships

The phage genome was partitioned into clusters and subclusters, using dot mapping. As per Hatfull’s work [35,36], two genomes can be assigned to the same cluster if their sequences in the dot plot exhibit a similarity higher than 50%. Moreover, a cluster can be partitioned into subclusters if the relationships within the cluster are heterogeneous. All the genomes within the same subcluster typically exhibit a greater sequence similarity. The dot-plot results for WIVsmall and other F1 cluster phages indicate that WIVsmall is likely a member of the F1 subcluster (Figure 5A). To investigate the evolutionary position of WIVsmall within the mycobacteriophage family, we conducted the phylogenetic analysis of the tape-measure protein (ORF67), which is the longest gene in the mycobacteriophage genomes. The tape-measure-protein-encoding gene is highly conserved, and serves as a typical phylogenetic marker for mycobacteriophages [37]. The amino acid sequences of the tape-measure protein from WIVsmall, and 14 closely related mycobacteriophages (F1cluster) were aligned. A phylogenetic tree was then constructed, using MEGA X software. WIVsmall was found to be most closely related to SG4, which belongs to the F1 subcluster of mycobacteriophages (Figure 5B).

3.6. Comparative Genomics Analysis

Based on BLASTn analysis, 15 phages belonging to the mycobacteriophage F cluster exhibiteda significant similarity to the WIVsmall genome. Among them, the mycobacteriophage SG4 exhibited the highest similarity, 64%, with the genome (Table 2). Meanwhile, a comparative genome analysis was conducted among the phages WIVsmall, SG4, Bobi, and Boomer, using CGview. As depicted in Figure 6, the genes implicated in phage structure and assembly exhibited a high degree of similarity among the three genomes. By conducting a BLASTp comparison against the Actinobacteriophage Database (phaseDB.org), we identified three proteins (Gp58, Gp59, and Gp61) in the WIVsmall genome that exhibited no similarities with the mycobacteriophage belonging to cluster F1 (Table 3). Among them, Gp58 exhibited homology with proteins found in non-mycobacteriophages, including those of the Tsukamurella phage, Gordonia phage, and Rhodococcus phage. Notably, all of the host bacteria for these phages, as well as mycobacteria, belong to the Actinomycetales taxonomic order. Gp59 exhibiteda similarity to proteins found in mycobacteriophages, including IdentityCrisis, Shweta, Ruthiejr, Willsammy, and Taquito, none of which are classified as members of the F1 cluster. Gp61 displayed homology with the protein found in the phage Moosehead, which was obtained utilizing Gordonia terra 3612 as the host. However, Gp61 also exhibited a significant homology with proteins originating from mycobacteriophages, including Wilder, MkaliMitinis3, LilDestination, and Lewan. Notably, none of these mycobacteriophages are affiliated with the F1 cluster. After performing BLASTpanalysis against the Actinobacteriophage Database, we found that the remaining 100 proteins exhibited a homology with other mycobacteriophages belonging to the F1 cluster. However, the coverage and identity values of two proteins, Gp34.1 and Gp60, were below 40%, casting doubt on the functional analysis results obtained through sequence alignment.

3.7. Biological Characteristics of the WIVsmall Phage

When cultured at temperatures of 50 °C or lower, the bacteriophage WIVsmall remained infectious to M. smegmatis mc2155, demonstrating an exceptional thermal stability. However, the phage titer gradually decreased in the water bath above 50 °C. Consistently, the titer of the phage WIVsmall remained undetectable at temperatures exceeding 70 °C (Figure 7A). The bacteriophage WIVsmall exhibited stability within the pH range of 4 to 11, over a duration of 1 h. However, under the acidic conditions of pH 2 to 3, the survival rate of the phages was negligible (Figure 7C). Furthermore, chloroform did not significantly impact WIVsmall’s infectivity toward M. smegmatis mc2155, indicating that WIVsmall is a lipid-free bacteriophage (Figure 7D).

3.8. Mass-Spectrometric Identification of Phage Proteins

The LC-MS/MS analysis identified 10 out of the 103 predicted ORF expression products. Among them, eight were phage structural proteins, including three types of phage minor tail protein, a putative structural protein, a tail-length tape measure protein, a major tail protein, a major capsid subunit, and a putative portal protein (Table 4). The absence of a significantly similar sequence in the database led to the classification of two additional proteins, ORF61 and ORF66, as uncharacterized proteins. However, the LC-MS/MS analysis failed to detect other predictions as hypothetical and functional proteins. This could be attributed to the low expression levels of the aforementioned proteins, which are beyond the detection limit of LC-MS/MS analysis. Additionally, these structural proteins exhibit a higher degree of similarity to mycobacteriophages possessing comparable genome sizes, such as SG4, Ramsey, Squirty, and Job42, than to other phages.

4. Discussion

In recent years, the irrational use of antibiotics has resulted in the continuous emergence of drug-resistant bacteria, and even “superbugs”, ushering in a “post-antibiotic era”, where in effective antibiotics are no longer available to combat infections in humans. M. tuberculosis is widely recognized as one of the most formidable drug-resistant bacterial pathogens in clinical settings. Multidrug-resistant tuberculosis (MDR-TB) and extensively drug-resistant tuberculosis (XDR-TB) constitute 20% of all cases of primary tuberculosis [38]. Several studies have tackled this challenge by exploring the immense potential of mycobacteriophages in TB detection and treatment [39]. In this study, we have identified and sequenced a novel phage, WIVsmall, that infects M. smegmatis mc2155. Dot-plot and phylogenetic analyses indicate that WIVsmall belongs to the F1 subcluster, as a novel member. According to the Actinobacteriophage Database (PhageDB.org), the genome length of the phages belonging to the F1 cluster ranges from 52,141 bp to 61,164 bp, while the number of genes encoded by the phages in this cluster ranges from 88 to 113. The total length of the WIVsmall genome is 53,359 bp, comprising 103 coding genes, which falls within the range of the phage F1 cluster. Comparative genomic analysis revealed that WIVsmall shares a maximum sequence similarity of only 64% with other mycobacteriophages. WIVsmall exhibited remarkable stability across diverse conditions, encompassing pH and temperature, which suggests its potential suitability in clinical settings. These findings provide new insights into the diversity and evolution of mycobacteriophages, and highlight the potential of WIVsmall as a model system for studying phage–host interactions. Recently, the focus of engineered phage development has primarily been on phages targeting S. aureus [40], E. coli [41], and P. aeruginosa [42], with comparatively fewer studies on engineered mycobacteriophages. Hence, it is imperative to delve deeper into bioengineering research concerning mycobacteriophages in this context.

Author Contributions

Conceptualization, T.T. and Q.L. methodology, T.T.; writing—original draft preparation, X.G.; writing—review and editing, J.Z.; supervision, F.Z. and Y.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by The Key R and D and Promotion Projects of Henan Province, grant number 222102310621; China Postdoctoral Science Foundation, grant number 2021m690095; Natural Science Foundation of Henan Province, grant number 202300410052.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The GenBank accession number for the phage WIVsmall is KC736071.1.

Acknowledgments

The authors wish to thank Yanzhang Li and Sixuan Lv for valuable help on the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Darby, E.M.; Trampari, E. Molecular mechanisms of antibiotic resistance revisited. Nat. Rev. Microbiol. 2023, 21, 280–295. [Google Scholar] [CrossRef] [PubMed]
  2. Wang, C.H.; Hsieh, Y.H.; Powers, Z.M.; Kao, C.Y. Defeating Antibiotic-Resistant Bacteria: Exploring Alternative Therapies for a Post-Antibiotic Era. Int. J. Mol. Sci. 2020, 21, 1061. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Carrigy, N.B.; Larsen, S.E.; Reese, V.; Pecor, T.; Harrison, M.; Kuehl, P.J.; Hatfull, G.F.; Sauvageau, D.; Baldwin, S.L.; Finlay, W.H.; et al. Prophylaxis of Mycobacterium tuberculosis H37Rv Infection in a Preclinical Mouse Model via Inhalation of Nebulized Bacteriophage D29. Antimicrob. Agents Chemother. 2019, 63, e00871. [Google Scholar] [CrossRef]
  4. Li, X.; He, Y.; Wang, Z.; Wei, J.; Hu, T.; Si, J.; Tao, G.; Zhang, L.; Xie, L.; Abdalla, A.E.; et al. A combination therapy of Phages and Antibiotics: Two is better than one. Int. J. Biol. Sci. 2021, 17, 3573–3582. [Google Scholar] [CrossRef]
  5. Jariah, R.O.A.; Hakim, M.S. Interaction of phages, bacteria, and the human immune system: Evolutionary changes in phage therapy. Rev. Med. Virol. 2019, 29, e2055. [Google Scholar] [CrossRef] [PubMed]
  6. Bao, J.; Wu, N.; Zeng, Y.; Chen, L.; Li, L.; Yang, L.; Zhang, Y.; Guo, M.; Li, L.; Li, J.; et al. Non-active antibiotic and bacteriophage synergism to successfully treat recurrent urinary tract infection caused by extensively drug-resistant Klebsiella pneumoniae. Emerg. Microbes Infect. 2020, 9, 771–774. [Google Scholar] [CrossRef] [Green Version]
  7. Nick, J.A.; Dedrick, R.M.; Gray, A.L.; Vladar, E.K.; Smith, B.E.; Freeman, K.G.; Malcolm, K.C.; Epperson, L.E.; Hasan, N.A.; Hendrix, J.; et al. Host and pathogen response to bacteriophage engineered against Mycobacterium abscessus lung infection. Cell 2022, 185, 1860–1874.e12. [Google Scholar] [CrossRef]
  8. Dedrick, R.M.; Freeman, K.G.; Nguyen, J.A.; Bahadirli-Talbott, A.; Cardin, M.E.; Cristinziano, M.; Smith, B.E.; Jeong, S.; Ignatius, E.H.; Lin, C.T.; et al. Nebulized Bacteriophage in a Patient With Refractory Mycobacterium abscessus Lung Disease. Open Forum Infect. Dis. 2022, 9, ofac194. [Google Scholar] [CrossRef]
  9. Tarín-Pelló, A.; Suay-García, B.; Pérez-Gracia, M.T. Antibiotic resistant bacteria: Current situation and treatment options to accelerate the development of a new antimicrobial arsenal. Expert Rev. Anti-Infect. Ther. 2022, 20, 1095–1108. [Google Scholar] [CrossRef]
  10. Chakaya, J.; Petersen, E.; Nantanda, R.; Mungai, B.N.; Migliori, G.B.; Amanullah, F.; Lungu, P.; Ntoumi, F.; Kumarasamy, N.; Maeurer, M.; et al. The WHO Global Tuberculosis 2021 Report—Not so good news and turning the tide back to End TB. Int. J. Infect. Dis. 2022, 124 (Suppl. 1), S26–S29. [Google Scholar] [CrossRef]
  11. Traoré, A.N.; Rikhotso, M.C. Effectiveness of the Novel Anti-TB Bedaquiline against Drug-Resistant TB in Africa: A Systematic Review of the Literature. Pathogens 2022, 11, 636. [Google Scholar] [CrossRef] [PubMed]
  12. Walker, T.M.; Miotto, P.; Köser, C.U.; Fowler, P.W.; Knaggs, J.; Iqbal, Z.; Hunt, M.; Chindelevitch, L.; Farhat, M.; Cirillo, D.M.; et al. The 2021 WHO catalogue of Mycobacterium tuberculosis complex mutations associated with drug resistance: A genotypic analysis. Lancet Microbe 2022, 3, e265–e273. [Google Scholar] [CrossRef] [PubMed]
  13. Hosseiniporgham, S.; Sechi, L.A. A Review on Mycobacteriophages: From Classification to Applications. Pathogens 2022, 11, 777. [Google Scholar] [CrossRef]
  14. Dedrick, R.M.; Smith, B.E.; Cristinziano, M.; Freeman, K.G.; Jacobs-Sera, D.; Belessis, Y.; Whitney Brown, A.; Cohen, K.A.; Davidson, R.M.; van Duin, D.; et al. Phage Therapy of Mycobacterium Infections: Compassionate Use of Phages in 20 Patients With Drug-Resistant Mycobacterial Disease. Clin. Infect. Dis. 2023, 76, 103–112. [Google Scholar] [CrossRef] [PubMed]
  15. Little, J.S.; Dedrick, R.M.; Freeman, K.G.; Cristinziano, M.; Smith, B.E.; Benson, C.A.; Jhaveri, T.A.; Baden, L.R.; Solomon, D.A.; Hatfull, G.F. Bacteriophage treatment of disseminated cutaneous Mycobacterium chelonae infection. Nat. Commun. 2022, 13, 2313. [Google Scholar] [CrossRef] [PubMed]
  16. Rondón, L.; Urdániz, E.; Latini, C.; Payaslian, F.; Matteo, M.; Sosa, E.J.; Do Porto, D.F.; Turjanski, A.G.; Nemirovsky, S.; Hatfull, G.F.; et al. Fluoromycobacteriophages Can Detect Viable Mycobacterium tuberculosis and Determine Phenotypic Rifampicin Resistance in 3-5 Days From Sputum Collection. Front. Microbiol. 2018, 9, 1471. [Google Scholar] [CrossRef]
  17. Li, X.; Hu, T.; Wei, J.; He, Y.; Abdalla, A.E.; Wang, G.; Li, Y.; Teng, T. Characterization of a Novel Bacteriophage Henu2 and Evaluation of the Synergistic Antibacterial Activity of Phage-Antibiotics. Antibiotics 2021, 10, 174. [Google Scholar] [CrossRef]
  18. Teng, T.; Li, Q.; Liu, Z.; Li, X.; Liu, Z.; Liu, H.; Liu, F.; Xie, L.; Wang, H.; Zhang, L.; et al. Characterization and genome analysis of novel Klebsiella phage Henu1 with lytic activity against clinical strains of Klebsiella pneumoniae. Arch. Virol. 2019, 164, 2389–2393. [Google Scholar] [CrossRef]
  19. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
  20. Teng, T.; Zhang, G.; Fan, X.; Zhang, Z.; Zhang, L.; Wu, D.; Chen, S.; Li, Y.; Jin, J. Complete genome sequence analysis of PS2, a novel T4-like bacteriophage that infects Serratia marcescens clinical isolates. Arch. Virol. 2018, 163, 1997–2000. [Google Scholar] [CrossRef]
  21. Wdowiak, M.; Paczesny, J.; Raza, S. Enhancing the Stability of Bacteriophages Using Physical, Chemical, and Nano-Based Approaches: A Review. Pharmaceutics 2022, 14, 1936. [Google Scholar] [CrossRef] [PubMed]
  22. Besemer, J.; Borodovsky, M. GeneMark: Web software for gene finding in prokaryotes, eukaryotes and viruses. Nucleic Acids Res. 2005, 33, W451–W454. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Mistry, J.; Chuguransky, S.; Williams, L.; Qureshi, M.; Salazar, G.A.; Sonnhammer, E.L.L.; Tosatto, S.C.E.; Paladin, L.; Raj, S.; Richardson, L.J.; et al. Pfam: The protein families database in 2021. Nucleic Acids Res. 2021, 49, D412–D419. [Google Scholar] [CrossRef]
  24. Chan, P.P.; Lin, B.Y.; Mak, A.J.; Lowe, T.M. tRNAscan-SE 2.0: Improved detection and functional classification of transfer RNA genes. Nucleic Acids Res. 2021, 49, 9077–9096. [Google Scholar] [CrossRef]
  25. Laslett, D.; Canback, B. ARAGORN, a program to detect tRNA genes and tmRNA genes in nucleotide sequences. Nucleic Acids Res. 2004, 32, 11–16. [Google Scholar] [CrossRef] [PubMed]
  26. Stothard, P.; Grant, J.R.; Van Domselaar, G. Visualizing and comparing circular genomes using the CGView family of tools. Brief. Bioinform. 2019, 20, 1576–1582. [Google Scholar] [CrossRef] [Green Version]
  27. Carver, T.; Thomson, N.; Bleasby, A.; Berriman, M.; Parkhill, J. DNAPlotter: Circular and linear interactive genome visualization. Bioinformatics 2009, 25, 119–120. [Google Scholar] [CrossRef] [Green Version]
  28. Gudlavalleti, B.S.; Phung, T.; Barton, C.L.; Becker, A.; Graul, B.L.; Griffin, J.T.; Hays, C.J.; Horn, B.; Liang, D.R.; Rutledge, L.M.; et al. Whole genome sequencing identifies an allele responsible for clear vs. turbid plaque morphology in a Mycobacteriophage. BMC Microbiol. 2020, 20, 148. [Google Scholar] [CrossRef]
  29. Kumar, V.; Loganathan, P.; Sivaramakrishnan, G.; Kriakov, J.; Dusthakeer, A.; Subramanyam, B.; Chan, J.; Jacobs, W.R., Jr.; Paranji Rama, N. Characterization of temperate phage Che12 and construction of a new tool for diagnosis of tuberculosis. Tuberculosis 2008, 88, 616–623. [Google Scholar] [CrossRef] [Green Version]
  30. Lo, C.Y.; Gao, Y. DNA Helicase-Polymerase Coupling in Bacteriophage DNA Replication. Viruses 2021, 13, 1739. [Google Scholar] [CrossRef]
  31. Abouhmad, A.; Korany, A.H.; Hatti-Kaul, R.; Grey, C.; Dishisha, T. Exploring the Enzymatic and Antibacterial Activities of Novel Mycobacteriophage Lysin B Enzymes. Int. J. Mol. Sci. 2020, 21, 3176. [Google Scholar] [CrossRef] [PubMed]
  32. Catalão, M.J.; Pimentel, M. Mycobacteriophage Lysis Enzymes: Targeting the Mycobacterial Cell Envelope. Viruses 2018, 10, 428. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Bavda, V.R.; Yadav, A.; Jain, V. Decoding the molecular properties of mycobacteriophage D29 Holin provides insights into Holin engineering. J. Virol. 2021, 95, e02173-20. [Google Scholar] [CrossRef]
  34. Ghai, R.; Mehrshad, M.; Mizuno, C.M.; Rodriguez-Valera, F. Metagenomic recovery of phage genomes of uncultured freshwater actinobacteria. ISME J. 2017, 11, 304–308. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Hatfull, G.F. Mycobacteriophages: Genes and genomes. Annu. Rev. Microbiol. 2010, 64, 331–356. [Google Scholar] [CrossRef] [PubMed]
  36. Hatfull, G.F.; Jacobs-Sera, D.; Lawrence, J.G.; Pope, W.H.; Russell, D.A.; Ko, C.C.; Weber, R.J.; Patel, M.C.; Germane, K.L.; Edgar, R.H.; et al. Comparative genomic analysis of 60 Mycobacteriophage genomes: Genome clustering, gene acquisition, and gene size. J. Mol. Biol. 2010, 397, 119–143. [Google Scholar] [CrossRef] [Green Version]
  37. Gissendanner, C.R.; Wiedemeier, A.M.; Wiedemeier, P.D.; Minton, R.L.; Bhuiyan, S.; Harmson, J.S.; Findley, A.M. A web-based restriction endonuclease tool for mycobacteriophage cluster prediction. J. Basic Microbiol. 2014, 54, 1140–1145. [Google Scholar] [CrossRef] [PubMed]
  38. Paton, N.I.; Cousins, C.; Suresh, C.; Burhan, E.; Chew, K.L.; Dalay, V.B.; Lu, Q.; Kusmiati, T.; Balanag, V.M.; Lee, S.L.; et al. Treatment Strategy for Rifampin-Susceptible Tuberculosis. N. Engl. J. Med. 2023, 388, 873–887. [Google Scholar] [CrossRef]
  39. Hatfull, G.F. Mycobacteriophages: From Petri dish to patient. PLoS Pathog. 2022, 18, e1010602. [Google Scholar] [CrossRef]
  40. Peng, H.; Borg, R.E.; Dow, L.P.; Pruitt, B.L.; Chen, I.A. Controlled phage therapy by photothermal ablation of specific bacterial species using gold nanorods targeted by chimeric phages. Proc. Natl. Acad. Sci. USA 2020, 117, 1951–1961. [Google Scholar] [CrossRef] [Green Version]
  41. Yehl, K.; Lemire, S.; Yang, A.C.; Ando, H.; Mimee, M.; Torres, M.T.; de la Fuente-Nunez, C.; Lu, T.K. Engineering Phage Host-Range and Suppressing Bacterial Resistance through Phage Tail Fiber Mutagenesis. Cell 2019, 179, 459–469.e9. [Google Scholar] [CrossRef] [PubMed]
  42. Kot, W.; Hammer, K.; Neve, H.; Vogensen, F.K. Identification of the receptor-binding protein in lytic Leuconostoc pseudomesenteroides bacteriophages. Appl. Environ. Microbiol. 2013, 79, 3311–3314. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. (A) Antibiotic replacement therapy strategies; (B) WHO list of critical pathogens for the development of new antibiotics.
Figure 1. (A) Antibiotic replacement therapy strategies; (B) WHO list of critical pathogens for the development of new antibiotics.
Cimb 45 00406 g001
Figure 2. (A) Plaque morphology of WIVsmall; (B,C) TEM showed that the phage WIVsmall had an equal-length head attached to a retractable tail.
Figure 2. (A) Plaque morphology of WIVsmall; (B,C) TEM showed that the phage WIVsmall had an equal-length head attached to a retractable tail.
Cimb 45 00406 g002
Figure 3. Single-step growth experiment of the mycobacteriophage WIVsmall (A); and analysis of the optimal MOI (B).
Figure 3. Single-step growth experiment of the mycobacteriophage WIVsmall (A); and analysis of the optimal MOI (B).
Cimb 45 00406 g003
Figure 4. Map of the phage WIVsmall genome.
Figure 4. Map of the phage WIVsmall genome.
Cimb 45 00406 g004
Figure 5. (A) Dot plot of parts of the F cluster mycobacteriophage genomes, displayed using Gepard. All genome sequences are connected into one sequence, so that related genomes are adjacent to each other. (B) Phylogenetic tree of the mycobacteriophage WIVsmall isolates. The neighbor-joining tree is based on the amino acid sequence alignment of the available sequences in GenBank and the mycobacteriophage database.
Figure 5. (A) Dot plot of parts of the F cluster mycobacteriophage genomes, displayed using Gepard. All genome sequences are connected into one sequence, so that related genomes are adjacent to each other. (B) Phylogenetic tree of the mycobacteriophage WIVsmall isolates. The neighbor-joining tree is based on the amino acid sequence alignment of the available sequences in GenBank and the mycobacteriophage database.
Cimb 45 00406 g005
Figure 6. Genetic and physical maps were prepared using CGview. Blast 1–3 revealed the Genomic sequence similarity of WIVsmall with SG4, Bobi, and Boomer, respectively. Three annular trajectories were described (from inside to outside): GC tilt ([G − C]/[G + C]), with inward peaks indicating a larger proportion of G; GC content (the inner peak is lower than the average GC content); and ORFs and transcription direction.
Figure 6. Genetic and physical maps were prepared using CGview. Blast 1–3 revealed the Genomic sequence similarity of WIVsmall with SG4, Bobi, and Boomer, respectively. Three annular trajectories were described (from inside to outside): GC tilt ([G − C]/[G + C]), with inward peaks indicating a larger proportion of G; GC content (the inner peak is lower than the average GC content); and ORFs and transcription direction.
Cimb 45 00406 g006
Figure 7. (A) Thermal stability test; (B) UV stability test; (C) pH stability test; (D) and chloroform sensitivity.
Figure 7. (A) Thermal stability test; (B) UV stability test; (C) pH stability test; (D) and chloroform sensitivity.
Cimb 45 00406 g007
Table 1. Predicted molecular functions of the gene products of the phage WIVsmall.
Table 1. Predicted molecular functions of the gene products of the phage WIVsmall.
ORF
Number
Start and Stop
Position
StrandLength (bp)MW
(kDa)
Molecular FunctionE-ValueCoverage
(%)
Identity (%)
183–69761522.83Glucosyltransferase2E-14210095.59
2694–9182258.08Hypothetical protein2E-4310095.59
3989–12012137.78Hypothetical protein2E-4210097.14
41204–2823162061.18Glycosyltransferase0.00 10093.70
52820–30292107.46Hypothetical protein7E-43100100
63059–346040214.36Hypothetical protein4E-8910099.25
6.13457–398452820.44Hypothetical protein3E-697186.40
6.23702–38391384.86Hypothetical protein1E-24100100
73866–40301656.07GIY-YIG nuclease8E-259095.92
7.13936–4181+2468.75Hypothetical protein1E-277795
84217–481960322.39HtrLYibB protein2E-14310098.50
94816–549668125.77Hypothetical protein2E-162100100
105499–579529711.10Hypothetical protein3E-5384100
115832–6020+1896.66Hypothetical protein2E-38100100
11.16167–646930310.65Hypothetical protein5E-216780.36
126531–67161867.12Hypothetical protein9E-3610098.36
136719–707535713.68Mobile element MPME2E-799699.12
147158–73852288.38Hypothetical protein9E-439795.89
157444–793248917.83Hypothetical protein2E-11310099.38
15.17929–80661384.71Hypothetical protein8E-2210097.78
168063–847040814.38Hypothetical protein1E-379568.09
16.18393–8683+29110.59Hypothetical protein3E-408591
16.28467–86611957.07Hypothetical protein1E-3810098.44
178661–88852258.37Hypothetical protein3E-4010093.24
17.19041–93042649.31Hypothetical protein1E-5210095.40
189301–10,26396334.65Replication initiation
protein
0.010098.44
18.110,474–10,7162439.25Hypothetical protein4E-4210086.25
1910,632–10,7781475.40Hypothetical protein5E-27100100
2010,778–10,8821053.79Hypothetical protein2E-1410097.06
20.110,879–11,24136313.72HNH endonuclease6E-8310098.33
2111,238–11,81657920.60Single strand annealing protein3E-13510099.48
2211,813–12,29548318.13HNH endonuclease5E-11310097.50
2312,292–12,66637513.94DNA methyltransferase2E-6710089.68
2412,663–12,8932318.26Hypothetical protein3E-4310097.37
2513,048–13,41336613.39DNA binding protein3E-8310099.17
2613,410–13,6432348.57WhiB transcriptional factor3E-449498.63
2713,692–14,12943816.57Hypothetical protein2E-98100100
2814,174–14,66849518.60WhiB transcriptional factor6E-11710098.78
2914,668–15,01534812.35Hypothetical protein3E-629979.82
3015,012–15,29027910.67Hypothetical protein7E-5610093.48
3115,369–16,16980127.87Hypothetical protein2E-10910076.69
3215,379–16,041+66322.30Hypothetical protein1E-139100100
3316,199–16,76556720.32Hypothetical protein5E-2810056.12
33.116,873–17,0311596.10Hypothetical protein2E-27100100
33.217,058–17,48943216.18HNH endonuclease9E-9110078.57
3417,489–17,7462589.82Excisionase1E-5110096.51
34.117,727–18,215995.87Hypothetical protein2E-1110096.88
34.217,743–17,84132411.63Hypothetical protein1E-426698.59
3517,838–18,16132411.63Hypothetical protein1E-426698.59
3618,076–18,89181630.12Hypothetical protein2E-269532
3718,907–19,50660021.83Hypothetical protein2E-859867.30
3819,729–20,03430611.50Transcriptional repressor2E-6710099.01
3920,218–20,721+50418.73Transcriptional repressor3E-11210097.60
4020,984–21,430+44715.55Pin protein2E-8310091.33
4121,522–22,640111941.52Integrase0.010099.73
4221,960–22,634+67522.65Hypothetical protein2E-148100100
4323,014–23,250+2378.49Hypothetical protein5E-229462.67
4423,250–23,417+1688.49Hypothetical protein6E-3110098.18
4523,399–23,75836013.37Hypothetical protein2E-779498.23
4623,928–24,161+2348.39Hypothetical protein1E-4610098.70
4724,234–25,05582231.03DNA polymerase exonuclease subunit0.010099.63
47.125,042–25,3022619.92Hypothetical protein4E-56100100
4825,224–25,93170825.33Hypothetical protein8E-169100100
4925,299–25,67337514.21Minor tail protein2E-8010099.19
49.125,670–25,9032347.78Holin7E-4310098.70
5025,920–26,921100236.63Endolysin0.010098.80
5126,921–28,198127847.42Endolysin0.010095.29
5228,195–28,4282348.29Hypothetical protein2E-4110098.70
5328,499–28,6391414.94Hypothetical protein1E-2210095.65
5428,636–28,7851505.38Hypothetical protein7E-219585.11
5528,792–29,08829710.31Hypothetical protein4E-4110084.69
55.129,232–29,444+2138.16Hypothetical protein4E-3610092.86
5629,452–30,29184026.26Minor tail protein2E-9810090.32
5730,288–30,81552817.50Minor tail protein2E-508461.49
5830,830–31,28245315.63Hypothetical protein7E-4610062.00
5931,299–31,745+44715.54Hypothetical protein7E-083256.00
6031,468–33,264179759.21Hypothetical proteinMinor tail protein7E-5038
6132,291–33,682139245.56Minor tail protein7E-076725
61.133,395–33,6102168.00Hypothetical protein6E-1610052.44
6233,660–35,447178863.70Minor tail protein6E-946349.00
6335,444–36,28984628.88Minor tail protein0.09998.21
6436,332–38,065+173464.58Minor tail protein0.09898.77
6538,124–39,833171063.44Putative structural protein0.010098.95
6639,819–42,185+236788.37Hypothetical protein0.0100100
6739,834–43,472+3639122.34Taillength tape measure protein0.010098.51
6842,040–42,63359420.41Hypothetical protein7E-136100100
6942,931–43,46453418.59Hypothetical protein1E-113100100
7043,472–44,01154020.12HNH endonuclease1E-12910098.88
7144,035–44,53850419.21Tail assembly chaperone4E-11410098.20
7244,420–44,97155220.50Tail assembly chaperone2E-12910099.45
7345,090–45,89981029.78Major tail protein0.010099.26
7446,008–46,41240514.70Head–tail adaptor6E-8910099.25
7546,402–46,6262257.83Head–tail connector
protein
8E-4610098.65
75.146,735–47,06433011.74Head–tail adaptor3E-7010098.17
7647,061–47,43237213.21Head–tail adaptor Ad15E-8110099.19
7747,445–48,35390931.71Major capsid subunit5E-1629974.67
7848,416–48,95253719.55Head scaffolding protein8E-12210099.44
7949,031–49,58255220.14Head maturation
protease
1E-13110098.91
8049,740–51,161142251.43Portal protein0.010099.15
8151,202–52,689148853.90Terminase0.010099.60
8252,661–52,9152559.36Terminase small subunit4E-5110098.81
8352,988–53,32033312.13HNH endonuclease3E-7310099.09
83.153,146–53,3582137.59Hypothetical protein9E-369598.51
“+” respresent the sense strand; “−” respresent the nonsense strand.
Table 2. Summary of similar genomic sequences with the phage WIVsmall.
Table 2. Summary of similar genomic sequences with the phage WIVsmall.
Phage NameQuery
Cover
IdentityAccession
Number
Genome Size (bp)
SG4 64%95.74%NC_026593.159,419
Ramsey63%96%NC_011289.158,578
Coco1262%95.34%NC_051644.157,693
Job4262%97.47%NC_021538.159,626
BuzzLyseyear56%88.73%NC_023699.159,016
ShiLan48%93.36%NC_041988.159,794
Pacc4047%93.53%NC_011287.158,554
Squirty37%95%NC_026588.160,285
Jabbawokkie22%95.03%NC_022069.155,213
Yoshi18%80.90%NC_042030.158,714
Table 3. Summary of unique proteins with the phage WIVsmall.
Table 3. Summary of unique proteins with the phage WIVsmall.
ProteinStrandLocationAA
Length
Best Match
Phage_GeneHostPhage
Cluster
Predicted
Function
Cov. %/
Ident. %
E
Value
Gp58-30,830–
31,282
453TPA4_27Tsukamurella paurometabola CON55SingletonHypothetical protein100/747E-62
GMA1_27Gordonia malaquae G239SingletonHypothetical protein100/627E-46
REQ3_58Rhodococcus equi Requ28SingletonHypothetical protein100/582E-45
PatrickStar_38Gordonia terrae 3612CXHypothetical protein96/583E-43
Keelan_56Gordonia terrae 3612DPHypothetical protein100/552E-39
Gp59-31,299–
31,745
447IdentityCrisis_23Mycobacterium smegmatis mc2155SingletonHypothetical protein32/567E-08
Shweta_22Mycobacterium smegmatis mc2155NHypothetical protein34/498E-07
Ruthiejr_30Mycobacterium smegmatis mc2155K4Hypothetical protein34/498E-07
Willsammy_23Mycobacterium smegmatis mc2155P1Hypothetical protein34/467E-06
Taquito_28Mycobacterium smegmatis mc2155K4Hypothetical protein34/467E-06
Gp61-32,291–
33,682
1392Moosehead_28Gordonia terrae
3612
CZ6Hypothetical protein 67/257E-07
Wilder_24Mycobacterium smegmatis mc2155L2Minor tail protein40/231E-06
MkaliMitinis3_24Mycobacterium smegmatis mc2155L2Minor tail protein50/231E-06
LilDestine_24Mycobacterium smegmatis mc2155L2Minor tail protein50/231E-06
Lewan_24Mycobacterium smegmatis mc2155SingletonMinor tail protein50/231E-06
Table 4. Mass spectrometry analysis of the mycobacteriophage WIVsmall.
Table 4. Mass spectrometry analysis of the mycobacteriophage WIVsmall.
No.Detected
Proteins
Predicated FunctionMolecular Mass
1Gp60Phage minor tail protein59.30 kDa
2Gp62Phage minor tail protein63.75 kDa
3Gp64Phage minor tail protein64.63 kDa
4Gp65Putative structural protein63.52 kDa
5Gp67Tail-length tape measure protein122.46 kDa
6Gp73Major tail protein29.82 kDa
7Gp77Major capsid subunit31.73 kDa
8Gp80Putative portal protein51.45 kDa
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Guo, X.; Zhang, J.; Wang, Y.; Zhou, F.; Li, Q.; Teng, T. Phenotypic Characterization and Comparative Genomic Analyses of Mycobacteriophage WIVsmall as A New Member Assigned to F1 Subcluster. Curr. Issues Mol. Biol. 2023, 45, 6432-6448. https://doi.org/10.3390/cimb45080406

AMA Style

Guo X, Zhang J, Wang Y, Zhou F, Li Q, Teng T. Phenotypic Characterization and Comparative Genomic Analyses of Mycobacteriophage WIVsmall as A New Member Assigned to F1 Subcluster. Current Issues in Molecular Biology. 2023; 45(8):6432-6448. https://doi.org/10.3390/cimb45080406

Chicago/Turabian Style

Guo, Xinge, Jing Zhang, Yuhan Wang, Fang Zhou, Qiming Li, and Tieshan Teng. 2023. "Phenotypic Characterization and Comparative Genomic Analyses of Mycobacteriophage WIVsmall as A New Member Assigned to F1 Subcluster" Current Issues in Molecular Biology 45, no. 8: 6432-6448. https://doi.org/10.3390/cimb45080406

APA Style

Guo, X., Zhang, J., Wang, Y., Zhou, F., Li, Q., & Teng, T. (2023). Phenotypic Characterization and Comparative Genomic Analyses of Mycobacteriophage WIVsmall as A New Member Assigned to F1 Subcluster. Current Issues in Molecular Biology, 45(8), 6432-6448. https://doi.org/10.3390/cimb45080406

Article Metrics

Back to TopTop