Genomic Diversity of Bacteriophages Infecting the Genus Acinetobacter

Comprehensive Approaches to Combatting Acinetobacter baumannii Biofilms: From Biofilm Structure to Phage-Based Therapies

Acinetobacter baumannii-a multidrug-resistant (MDR) pathogen that causes, for example, skin and soft tissue wounds; urinary tract infections; pneumonia; bacteremia; and endocarditis, particularly due to its ability to form robust biofilms-poses a significant challenge in clinical settings. This structure protects the bacteria from immune responses and antibiotic treatments, making infections difficult to eradicate. Given the rise in antibiotic resistance, alternative therapeutic approaches are urgently needed. Bacteriophage-based strategies have emerged as a promising solution for combating A. baumannii biofilms. Phages, which are viruses that specifically infect bacteria, offer a targeted and effective means of disrupting biofilm and lysing bacterial cells. This review explores the current advancements in bacteriophage therapy, focusing on its potential for treating A. baumannii biofilm-related infections. We described the mechanisms by which phages interact with biofilms, the challenges in phage therapy implementation, and the strategies being developed to enhance its efficacy (phage cocktails, engineered phages, combination therapies with antibiotics). Understanding the role of bacteriophages in both biofilm disruption and in inhibition of its forming could pave the way for innovative treatments in combating MDR A. baumannii infections as well as the prevention of their development.

Host population structure and species resolution reveal prophage transmission dynamics

Much knowledge about bacteriophages has been obtained via genomics and metagenomics over the last decades. However, most studies dealing with prophage diversity have rarely conducted phage species delimitation (aspect 1) and have hardly integrated the population structure of the host (aspect 2). Yet, these two aspects are essential in assessing phage diversity. Here, we implemented an operational definition of phage species (clustering at 95% identity, 90% coverage) and integrated the host's population structure to understand prophage diversity better. Gathering the most extensive data set of Acinetobacter baumannii phages (4,152 prophages + 122 virulent phages, distributed in 46 countries in the world), we show that 91% (875 out of 963) of the prophage species have four or fewer prophages per species, and just five prophage species have more than 100 prophages. Most prophage species have a narrow host range and are geographically restricted; yet, very few have a broad host range being well spread in distant lineages of A. baumannii. These few broad host range prophage species are not only cosmopolitan but also the most abundant species. We also noted that polylysogens had very divergent prophages, belonging to different prophage species, and prophages can easily be gained and lost within the bacterial lineages. Finally, even with this extensive data set, the prophage diversity has not been fully grasped. Our study highlights how integrating the host population structure and a solid operational definition of phage species allows us to better appreciate phage diversity and its transmission dynamics.

IMPORTANCE

Much knowledge about bacteriophages has been obtained via genomics and metagenomics over the last decades. However, most studies dealing with prophage diversity have rarely conducted phage species delimitation (aspect 1) and have hardly integrated the population structure of the host (aspect 2). Yet, these two aspects are essential in assessing phage diversity. Here, we implemented an operational definition of phage species (clustering at 95% identity, 90% coverage) and integrated the host's population structure to understand prophage diversity better. Gathering the most extensive data set of Acinetobacter baumannii phages, we show that most prophage species have four or fewer prophages per species, and just five prophage species have more than 100 prophages. Most prophage species have a narrow host range and are geographically restricted; yet, very few have a broad host range being well spread in distant lineages of A. baumannii. These few broad host range prophage species are cosmopolitan and the most abundant species. Prophages in the same bacterial genome are very divergent, and prophages can easily be gained and lost within the bacterial lineages. Finally, even with this extensive data set, the prophage diversity has not been fully grasped. This study shows how integrating the host population structure and clustering at the species level allows us to better appreciate phage diversity and its transmission dynamics.

Phage therapy combats pan drug-resistant Acinetobacter baumannii infection safely and efficiently

Phage therapy offers a promising approach to combat the growing threat of antimicrobial resistance. Yet, key questions remain regarding dosage, administration routes, combination therapy, and the causes of therapeutic failure. In this study, we focused on a novel lytic phage, ФAb4B, which specifically targeted the Acinetobacter baumannii strains with KL160 capsular polysaccharide, including the pan-drug resistant A. baumannii YQ4. ФAb4B exhibited the ability to effectively inhibit biofilm formation and eradicate mature biofilms independently of dosage. Additionally, it demonstrated a wide spectrum of antibiotic-phage synergy and did not show any cytotoxic or haemolytic effects. Continuous phage injections, both intraperitoneally and intravenously over 7 d, showed no acute toxicity in vivo. Importantly, phage therapy significantly improved neutrophil counts, outperforming ciprofloxacin. However, excessive phage injections suppressed neutrophil levels. The combinatorial treatment of phage-ciprofloxacin rescued 91% of the mice, a superior outcome compared to phage alone (67%). The efficacy of the combinatorial treatment was independent of phage dosage. Notably, prophylactic administration of the combinatorial regimen provided no protection, but even when combined with a delayed therapeutic regimen, it saved all the mice. Bacterial resistance to the phage was not a contributing factor to treatment failure. Our preclinical study systematically describes the lytic phage's effectiveness in both in vitro and in vivo settings, filling in crucial details about phage treatment against bacteriemia caused by A. baumannii, which will provide a robust foundation for the future of phage therapy.

In vivo efficacy of phage cocktails against carbapenem resistance Acinetobacter baumannii in the rat pneumonia model

Acinetobacter baumannii, an opportunistic pathogen, poses a significant threat in intensive care units, leading to severe nosocomial infections. The rise of multi-drug-resistant strains, particularly carbapenem-resistant A. baumannii, has created formidable challenges for effective treatment. Given the prolonged development cycle and high costs associated with antibiotics, phages have garnered clinical attention as an alternative for combating infections caused by drug-resistant bacteria. However, the utilization of phage therapy encounters notable challenges, including the narrow host spectrum, where each phage targets a limited subset of bacteria, increasing the risk of phage resistance development. Additionally, uncertainties in immune system dynamics during treatment hinder tailoring symptomatic interventions based on patient-specific states. In this study, we isolated two A. baumannii phages from wastewater and conducted a comprehensive assessment of their potential applications. This evaluation included sequencing analysis, genome classification, pH and temperature stability assessments, and in vitro bacterial inhibition assays. Further investigations involved analyzing histological and cytokine alterations in rats undergoing phage cocktail treatment for pneumonia. The therapeutic efficacy of the phages was validated, and transcriptomic studies of rat lung tissue during phage treatment revealed crucial changes in the immune system. The findings from our study underscore the potential of phages for future development as a treatment strategy and offer compelling evidence regarding immune system dynamics throughout the treatment process.IMPORTANCEDue to the growing problem of multi-drug-resistant bacteria, the use of phages is being considered as an alternative to antibiotics, and the genetic safety and application stability of phages determine the potential of phage application. The absence of drug resistance genes and virulence genes in the phage genome can ensure the safety of phage application, and the fact that phage can remain active in a wide range of temperatures and pH is also necessary for application. In addition, the effect evaluation of preclinical studies is especially important for clinical application. By simulating the immune response situation during the treatment process through mammalian models, the changes in animal immunity can be observed, and the effect of phage therapy can be further evaluated. Our study provides compelling evidence that phages hold promise for further development as therapeutic agents for Acinetobacter baumannii infections.

Lytic Capsule-Specific Acinetobacter Bacteriophages Encoding Polysaccharide-Degrading Enzymes

The genus Acinetobacter comprises both environmental and clinically relevant species associated with hospital-acquired infections. Among them, Acinetobacter baumannii is a critical priority bacterial pathogen, for which the research and development of new strategies for antimicrobial treatment are urgently needed. Acinetobacter spp. produce a variety of structurally diverse capsular polysaccharides (CPSs), which surround the bacterial cells with a thick protective layer. These surface structures are primary receptors for capsule-specific bacteriophages, that is, phages carrying tailspikes with CPS-depolymerizing/modifying activities. Phage tailspike proteins (TSPs) exhibit hydrolase, lyase, or esterase activities toward the corresponding CPSs of a certain structure. In this study, the data on all lytic capsule-specific phages infecting Acinetobacter spp. with genomes deposited in the NCBI GenBank database by January 2024 were summarized. Among the 149 identified TSPs encoded in the genomes of 143 phages, the capsular specificity (K specificity) of 46 proteins has been experimentally determined or predicted previously. The specificity of 63 TSPs toward CPSs, produced by various Acinetobacter K types, was predicted in this study using a bioinformatic analysis. A comprehensive phylogenetic analysis confirmed the prediction and revealed the possibility of the genetic exchange of gene regions corresponding to the CPS-recognizing/degrading parts of different TSPs between morphologically and taxonomically distant groups of capsule-specific Acinetobacter phages.

Resistance to bacteriophage incurs a cost to virulence in drug-resistant Acinetobacter baumannii

Introduction . Acinetobacter baumannii is a critical priority pathogen for novel antimicrobials (World Health Organization) because of the rise in nosocomial infections and its ability to evolve resistance to last resort antibiotics. A. baumannii is thus a priority target for phage therapeutics. Two strains of a novel, virulent bacteriophage (LemonAid and Tonic) able to infect carbapenem-resistant A. baumannii (strain NCTC 13420), were isolated from environmental water samples collected through a citizen science programme.Gap statement. Phage-host coevolution can lead to emergence of host resistance, with a concomitant reduction in the virulence of host bacteria; a potential benefit to phage therapy applications.Methodology. In vitro and in vivo assays, genomics and microscopy techniques were used to characterize the phages; determine mechanisms and impact of phage resistance on host virulence, and the efficacy of the phages against A. baumannii.Results. A. baumannii developed resistance to both viruses, LemonAid and Tonic. Resistance came at a cost to virulence, with the resistant variants causing significantly reduced mortality in a Galleria mellonella larval in vivo model. A replicated 8 bp insertion increased in frequency (~40 % higher frequency than in the wild-type) within phage-resistant A. baumannii mutants, putatively resulting in early truncation of a protein of unknown function. Evidence from comparative genomics and an adsorption assay suggests this protein acts as a novel phage receptor site in A. baumannii. We find no evidence linking resistance to changes in capsule structure, a known virulence factor. LemonAid efficiently suppressed growth of A. baumanni in vitro across a wide range of titres. However, in vivo, while survival of A. baumannii infected larvae significantly increased with both remedial and prophylactic treatment with LemonAid (107 p.f.u. ml-1), the effect was weak and not sufficient to save larvae from morbidity and mortality.Conclusion. While LemonAid and Tonic did not prove effective as a treatment in a Galleria larvae model, there is potential to harness their ability to attenuate virulence in drug-resistant A. baumannii.

Whole-Genome Sequencing and Phenotypic Analysis of Streptococcus equi subsp. zooepidemicus Sequence Type 147 Isolated from China

Streptococcus equi subsp. zooepidemicus (S. zooepidemicus) is one of the important zoonotic and opportunistic pathogens. In recent years, there has been growing evidence that supports the potential role of S. zooepidemicus in severe diseases in horses and other animals, including humans. Furthermore, the clinical isolation and drug resistance rates of S. zooepidemicus have been increasing yearly, leading to interest in its in-depth genomic analysis. In order to deepen the understanding of the S. zooepidemicus characteristics and genomic features, we investigated the genomic islands, mobile genetic elements, virulence and resistance genes, and phenotype of S. zooepidemicus strain ZHZ 211 (ST147), isolated from an equine farm in China. We obtained a 2.18 Mb, high-quality chromosome and found eight genomic islands. According to a comparative genomic investigation with other reference strains, ZHZ 211 has more virulence factors, like an iron uptake system, adherence, exoenzymes, and antiphagocytosis. More interestingly, ZHZ 211 has acquired a mobile genetic element (MGE), prophage Ph01, which was found to be in the chromosome of this strain and included two hyaluronidase (hyl) genes, important virulence factors of the strain. Moreover, two transposons and two virulence (virD4) genes were found to be located in the same genome island of ZHZ 211. In vitro phenotypic results showed that ZHZ 211 grows faster and is resistant to clarithromycin, enrofloxacin, and sulfonamides. The higher biofilm-forming capabilities of ZHZ 211 may provide a competitive advantage for survival in its niche. The results expand our understanding of the genomic, pathogenicity, and resistance characterization of Streptococcus zooepidemicus and facilitate further exploration of its molecular pathogenic mechanism.

New Obolenskvirus Phages Brutus and Scipio: Biology, Evolution, and Phage-Host Interaction

Two novel virulent phages of the genus Obolenskvirus infecting Acinetobacter baumannii, a significant nosocomial pathogen, have been isolated and studied. Phages Brutus and Scipio were able to infect A. baumannii strains belonging to the K116 and K82 capsular types, respectively. The biological properties and genomic organization of the phages were characterized. Comparative genomic, phylogenetic, and pangenomic analyses were performed to investigate the relationship of Brutus and Scipio to other bacterial viruses and to trace the possible origin and evolutionary history of these phages and other representatives of the genus Obolenskvirus. The investigation of enzymatic activity of the tailspike depolymerase encoded in the genome of phage Scipio, the first reported virus infecting A. baumannii of the K82 capsular type, was performed. The study of new representatives of the genus Obolenskvirus and mechanisms of action of depolymerases encoded in their genomes expands knowledge about the diversity of viruses within this taxonomic group and strategies of Obolenskvirus-host bacteria interaction.

Expression of a recombinant endolysin from bacteriophage CAP 10-3 with lytic activity against Cutibacterium acnes

The bacteriophage CAP 10-3 forming plaques against Cutibacterium acnes which causes skin acne was previously isolated from human skin acne lesion. Incomplete whole genome sequence (WGS) of the bacteriophage CAP 10-3 was obtained and it had 29,643 bp long nucleotide with 53.86% GC content. The sequence was similar to C. acnes phage PAP 1-1 with a nucleotide sequence identity of 89.63% and the bacteriophage belonged to Pahexavirus. Bioinformatic analysis of the WGS predicted 147 ORFs and functions of 40 CDSs were identified. The predicted endolysin gene of bacteriophage CAP 10-3 was 858 bp long which was deduced as 285 amino acids (~ 31 kDa). The protein had the highest similarity with amino acid sequence of the endolysin from Propionibacterium phage PHL071N05 with 97.20% identity. The CAP 10-3 endolysin gene was amplified by PCR with primer pairs based on the gene sequence, cloned into an expression vector pET-15b and transformed into Escherichia coli BL21(DE3) strain. The predicted protein band (~ 33 kDa) for the recombinant endolysin was detected in an SDS-PAGE gel and western blot assay. The concentrated supernatant of cell lysate from E. coli BL21(DE3) (pET-15b_CAP10-3 end) and a partially purified recombinant CAP 10-3 endolysin showed antibacterial activity against C. acnes KCTC 3314 in a dose-dependent manner. In conclusion, the recombinant CAP 10-3 endolysin was successfully produced in E. coli strain and it can be considered as a therapeutic agent candidate for treatment of human skin acne.

Characterization of a novel and active temperate phage vB_AbaM_ABMM1 with antibacterial activity against Acinetobacter baumannii infection

Acinetobacter baumannii is an opportunistic pathogen that significantly causes hospital-acquired infections. Due to its multidrug resistance, treating infections caused by this pathogen is challenging. Recently, phages have gained attention as a potential alternative to antibiotics in treating bacterial infections. While lytic phages are preferred in therapy, the use of temperate phages for this purpose has received less attention. This study characterized a novel temperate phage vB_AbaM_ABMM1 (ABMM1) with antibacterial activity toward A. baumannii. ABMM1 adsorbs quickly, has short latent periods, and is relatively stable at various temperatures and neutral pH. ABMM1 has an icosahedral head and a contractile tail. It has a 75,731 kb circular permuted dsDNA genome containing 86 gene products with 37.3% G + C content and a mosaic arrangement typical of temperate phages. Genomic analysis confirmed that ABMM1 does not have antibiotic-resistance genes or virulence-related factors. The packaging strategy was predicted in silico, suggesting that ABMM1 represents a headful phage. Only truncated ABMM1 prophage was detected and has similarity in the genome of several A. baumannii strains. Despite its ability to integrate into the host chromosome, the high MOI of ABMM1 (MOI 10) effectively killed the host bacterial cells and reduced the fatality rate of bacterial infection in the zebrafish model. These findings indicate that ABMM1 can be an alternative treatment for A. baumannii infection.

Characterization and Diversity of Klebsiella pneumoniae Prophages

Klebsiella pneumoniae is a common human commensal and opportunistic pathogen. In recent years, the clinical isolation and resistance rates of K. pneumoniae have shown a yearly increase, leading to a special interest in mobile genetic elements. Prophages are a representative class of mobile genetic elements that can carry host-friendly genes, transfer horizontally between strains, and coevolve with the host's genome. In this study, we identified 15,946 prophages from the genomes of 1437 fully assembled K. pneumoniae deposited in the NCBI database, with 9755 prophages on chromosomes and 6191 prophages on plasmids. We found prophages to be notably diverse and widely disseminated in the K. pneumoniae genomes. The K. pneumoniae prophages encoded multiple putative virulence factors and antibiotic resistance genes. The comparison of strain types with prophage types suggests that the two may be related. The differences in GC content between the same type of prophages and the genomic region in which they were located indicates the alien properties of the prophages. The overall distribution of GC content suggests that prophages integrated on chromosomes and plasmids may have different evolutionary characteristics. These results suggest a high prevalence of prophages in the K. pneumoniae genome and highlight the effect of prophages on strain characterization.

Friunavirus Phage-Encoded Depolymerases Specific to Different Capsular Types of Acinetobacter baumannii

Acinetobacter baumannii is a critical priority nosocomial pathogen that produces a variety of capsular polysaccharides (CPSs), the primary receptors for specific depolymerase-carrying phages. In this study, the tailspike depolymerases (TSDs) encoded in genomes of six novel Friunaviruses, APK09, APK14, APK16, APK86, APK127v, APK128, and one previously described Friunavirus phage, APK37.1, were characterized. For all TSDs, the mechanism of specific cleavage of corresponding A. baumannii capsular polysaccharides (CPSs) was established. The structures of oligosaccharide fragments derived from K9, K14, K16, K37/K3-v1, K86, K127, and K128 CPSs degradation by the recombinant depolymerases have been determined. The crystal structures of three of the studied TSDs were obtained. A significant reduction in mortality of Galleria mellonella larvae infected with A. baumannii of K9 capsular type was shown in the example of recombinant TSD APK09_gp48. The data obtained will provide a better understanding of the interaction of phage-bacterial host systems and will contribute to the formation of principles of rational usage of lytic phages and phage-derived enzymes as antibacterial agents.

Genomic Analysis of Two Novel Bacteriophages Infecting Acinetobacter beijerinckii and halotolerans Species

Bacteriophages are the most diverse genetic entities on Earth. In this study, two novel bacteriophages, nACB1 (Podoviridae morphotype) and nACB2 (Myoviridae morphotype), which infect Acinetobacter beijerinckii and Acinetobacter halotolerans, respectively, were isolated from sewage samples. The genome sequences of nACB1 and nACB2 revealed that their genome sizes were 80,310 bp and 136,560 bp, respectively. Comparative analysis showed that both genomes are novel members of the Schitoviridae and the Ackermannviridae families, sharing ≤ 40% overall nucleotide identities with any other phages. Interestingly, among other genetic features, nACB1 encoded a very large RNA polymerase, while nACB2 displayed three putative depolymerases (two capsular depolymerases and one capsular esterase) encoded in tandem. This is the first report of phages infecting A. halotolerans and beijerinckii human pathogenic species. The findings regarding these two phages will allow us to further explore phage-Acinetobacter interactions and the genetic evolution for this group of phages.

Loss of a Branch Sugar in the Acinetobacter baumannii K3-Type Capsular Polysaccharide Due To Frameshifts in the gtr6 Glycosyltransferase Gene Leads To Susceptibility To Phage APK37.1

The type of capsular polysaccharide (CPS) on the cell surface of Acinetobacter baumannii can determine the specificity of lytic bacteriophage under consideration for therapeutic use. Here, we report the isolation of a phage on an extensively antibiotic resistant ST2 A. baumannii isolate AB5001 that carries the KL3 CPS biosynthesis gene cluster predicting a K3-type CPS. As the phage did not infect isolates carrying KL3 or KL22 and known to produce K3 CPS, the structure of the CPS isolated from A. baumannii AB5001 was determined. AB5001 produced a variant CPS form, K3-v1, that lacks the β-d-GlсpNAc side chain attached to the d-Galp residue in the K3 structure. Inspection of the KL3 sequence in the genomes of AB5001 and other phage-susceptible isolates with a KL3 locus revealed single-base deletions in gtr6, causing loss of the Gtr6 glycosyltransferase that adds the missing d-GlсpNAc side chain to the K3 CPS. Hence, the presence of this sugar profoundly restricts the ability of the phage to digest the CPS. The 41-kb linear double-stranded DNA (dsDNA) phage genome was identical to the genome of a phage isolated on a K37-producing isolate and thus was named APK37.1. APK37.1 also infected isolates carrying KL116. Consistent with this, K3-v1 resembles the K37 and K116 structures. APK37.1 is a Friunavirus belonging to the Autographiviridae family. The phage-encoded tail spike depolymerase DpoAPK37.1 was not closely related to Dpo encoded by other sequenced Friunaviruses, including APK37 and APK116. IMPORTANCE Lytic bacteriophage have potential for the treatment of otherwise untreatable extensively antibiotic-resistant bacteria. For Acinetobacter baumannii, most phage exhibit specificity for the type of capsular polysaccharide (CPS) produced on the cell surface. However, resistance can arise via mutations in CPS genes that abolish this phage receptor. Here, we show that single-base deletions in a CPS gene result in alteration of the final structure rather than deletion of the capsule layer and hence affect the ability of a newly reported podophage to infect strains producing the K3 CPS.

Biological Properties of 12 Newly Isolated Acinetobacter baumannii-Specific Bacteriophages

Infections with the opportunistic Gram-negative bacterium Acinetobacter baumannii pose a serious threat today, which is aggravated by the growing problem of multi-drug resistance among bacteria, caused by the overuse of antibiotics. Treatment of infections caused by antibiotic-resistant A. baumannii strains with the use of phage therapy is not only a promising alternative, but sometimes the only option. Therefore, phages specific for clinical multi-drug resistant A. baumannii were searched for in environmental, municipal, and hospital wastewater samples collected from different locations in Poland. The conducted research allowed us to determine the biological properties and morphology of the tested phages. As a result of our research, 12 phages specific for A. baumannii, 11 of which turned out to be temperate and only one lytic, were isolated. Their lytic spectra ranged from 11 to 75%. The plaques formed by most phages were small and transparent, while one of them formed relatively large plaques with a clearly marked 'halo' effect. Based on Transmission Electron Microscopy (TEM), most of our phages have been classified as siphoviruses (only one phage was classified as a podovirus). All phages have icosahedral capsid symmetry, and 11 of them have a long tail. Optimal multiplicity of infections (MOIs) and the adsorption rate were also determined. MOI values varied depending on the phage-from 0.001 to 10. Based on similarities to known bacteriophages, our A. baumannii-specific phages have been proposed to belong to the Beijerinckvirinae and Junivirinae subfamilies. This study provides an additional tool in the fight against this important pathogen and may boost the interest in phage therapy as an alternative and supplement to the current antibiotics.

Comparative genomics of Acinetobacter baumannii and therapeutic bacteriophages from a patient undergoing phage therapy

In 2016, a 68-year-old patient with a disseminated multidrug-resistant Acinetobacter baumannii infection was successfully treated using lytic bacteriophages. Here we report the genomes of the nine phages used for treatment and three strains of A. baumannii isolated prior to and during treatment. The phages used in the initial treatment are related, T4-like myophages. Analysis of 19 A. baumannii isolates collected before and during phage treatment shows that resistance to the T4-like phages appeared two days following the start of treatment. We generate complete genomic sequences for three A. baumannii strains (TP1, TP2 and TP3) collected before and during treatment, supporting a clonal relationship. Furthermore, we use strain TP1 to select for increased resistance to five of the phages in vitro, and identify mutations that are also found in phage-insensitive isolates TP2 and TP3 (which evolved in vivo during phage treatment). These results support that in vitro investigations can produce results that are relevant to the in vivo environment.

Whole-Genome Analysis Reveals That Bacteriophages Promote Environmental Adaptation of Staphylococcus aureus via Gene Exchange, Acquisition, and Loss

The study of bacteriophages is experiencing a resurgence owing to their antibacterial efficacy, lack of side effects, and low production cost. Nonetheless, the interactions between Staphylococcus aureus bacteriophages and their hosts remain unexplored. In this study, whole-genome sequences of 188 S. aureus bacteriophages—20 Podoviridae, 56 Herelleviridae, and 112 Siphoviridae—were obtained from the National Center for Biotechnology Information (NCBI, USA) genome database. A phylogenetic tree was constructed to estimate their genetic relatedness using single-nucleotide polymorphism analysis. Comparative analysis was performed to investigate the structural diversity and ortholog groups in the subdividing clusters. Mosaic structures and gene content were compared in relation to phylogeny. Phylogenetic analysis revealed that the bacteriophages could be distinguished into three lineages (I–III), including nine subdividing clusters and seven singletons. The subdividing clusters shared similar mosaic structures and core ortholog clusters, including the genes involved in bacteriophage morphogenesis and DNA packaging. Notably, several functional modules of bacteriophages 187 and 2368A shared more than 95% nucleotide sequence identity with prophages in the S. aureus strain RJ1267 and the Staphylococcus pseudintermedius strain SP_11306_4, whereas other modules exhibited little nucleotide sequence similarity. Moreover, the cluster phages shared similar types of holins, lysins, and DNA packaging genes and harbored diverse genes associated with DNA replication and virulence. The data suggested that the genetic diversity of S. aureus bacteriophages was likely due to gene replacement, acquisition, and loss among staphylococcal phages, which may have crossed species barriers. Moreover, frequent module exchanges likely occurred exclusively among the subdividing cluster phages. We hypothesize that during evolution, the S. aureus phages enhanced their DNA replication in host cells and the adaptive environment of their host.