From Trees to Clouds: PhageClouds for Fast Comparison of ∼640,000 Phage Genomic Sequences and Host-Centric Visualization Using Genomic Network Graphs

Weberviruses are gut-associated phages that infect Klebsiella spp

Weberviruses are bacteriophages (phages) that can infect and lyse clinically relevant, multidrug-resistant (MDR) strains of Klebsiella. They are an attractive therapeutic option to tackle Klebsiella infections due to their high burst sizes, long shelf life, and associated depolymerases. In this study, we isolated and characterized seven new lytic phages and compared their genomes with those of their closest relatives. Gene-sharing network, ViPTree proteome, and terL gene-sequence-based analyses incorporating all publicly available webervirus genomes [n = 258 from isolates, n = 65 from metagenome-assembled genome (MAG) datasets] confirmed the seven phages as members of the genus Webervirus and identified a novel genus (Defiantjazzvirus) within the family Drexlerviridae. Using our curated database of 265 isolated phage genomes and 65 MAGs (n = 330 total), we found that weberviruses are distributed globally and primarily associated with samples originating from the gut: sewage (154/330, 47%), wastewater (83/330, 25%), and human faeces (66/330, 20%). We identified three distinct clusters of potential depolymerases encoded within the 330 genomes. Due to their global distribution, frequency of isolation and lytic activity against the MDR clinical Klebsiella strains used in this study, we conclude that weberviruses and their depolymerases show promise for development as therapeutic agents against Klebsiella spp.

Characterization and anti-biofilm potentiality of an isolated novel Aeromonas hydrophila-infecting bacteriophage AHPMCC11, belonging to the genus Ahphunavirus

Aeromonas hydrophila is a major aquatic habitat pathogen responsible for huge economic losses in the aquaculture and food industries. In this study, a lytic bacteriophage AHPMCC11 was isolated by using A. hydrophila MTCC 1739. AHPMCC11 showed a short latent period of 10 min and the burst size was 215 PFU/cell. AHPMCC11 had potent bacteriolytic activity within 2 h in liquid culture inhibition assay and exhibited biofilm scavenging activity against A. hydrophila MTCC 1739. AHPMCC11 was found stable at a wide range of pH levels (3-12), temperature ranges (4-37 °C), and salinity conditions (0-40 ppt). The AHPMCC11 genome was determined to be 42,439 bp in length with 58.9 % G + C content, 51 CDS, and no tRNA. Comparative genome study suggested that AHPMCC11 may represent a novel species within the Autographiviridae family, belonging to the Ahphunavirus genus. In conclusion, AHPMCC11 might be used as a biocontrol agent in aquaculture and the food industry.

Characterization and Comparative Genomic Analysis of vB_BceM_CEP1: A Novel Temperate Bacteriophage Infecting Burkholderia cepacia Complex

The increasing prevalence of multidrug-resistant bacteria imminently threatens public health and jeopardizes nearly all aspects of modern medicine. The Burkholderia cepacia complex (Bcc) comprises Burkholderia cepacia and the related species of Gram-negative bacteria. Members of the Bcc group are opportunistic pathogens responsible for various chronic illnesses, including cystic fibrosis and chronic granulomatous disease. Phage therapy is emerging as a potential solution to combat the antimicrobial resistance crisis. In this study, a temperate phage vB_BceM_CEP1 was isolated from sewage and fully characterized. Transmission electron microscopy indicated that vB_BceM_CEP1 belongs to the family Peduoviridae. The isolated phage demonstrated enhanced environmental stability and antibiofilm potential. One-step growth analysis revealed a latent period of 30 min and an average burst size of 139 plaque-forming units per cell. The genome of vB_BceM_CEP1 consists of 32,486 bp with a GC content of 62.05%. A total of 40 open reading frames were annotated in the phage genome, and none of the predicted genes was annotated as tRNA. Notably, genes associated with antibiotic resistance, host virulence factors, and toxins were absent from the vB_BceM_CEP1 genome. Based on its unique phenotype and phylogeny, the isolated phage vB_BceM_CEP1 is classified as a new temperate phage with lytic activity. The findings of this study enhance our understanding of the diversity of Bcc phages.

In vitro and in vivo assessment of the competence of a novel lytic phage vB_EcoS_UTEC10 targeting multidrug resistant Escherichia coli with a robust biofilm eradication activity

Escherichia coli (E. coli) is a leading cause of human infections worldwide and is considered a major cause of nosocomial infections, sepsis, meningitis and diarrhea. Lately, there has been an alarming increase in the incidence of antimicrobial resistance among clinical E. coli isolates. In the current study, a novel bacteriophage (phage) vB_EcoS_UTEC10 was isolated and characterized. The isolated phage showed high stability over wide temperature and pH ranges beside its promising bacteriolytic activity against multidrug resistant (MDR) E. coli isolates. In addition, vB_EcoS_UTEC10 showed a marked antibiofilm capability against mature E. coli biofilms. Genomic investigation revealed that vB_EcoS_UTEC10 has a double stranded DNA genome that consists of 44,772 bp comprising a total of 73 open reading frames (ORFs), out of which 35 ORFs were annotated as structural or functional proteins, and none were related to antimicrobial resistance or lysogeny. In vivo investigations revealed a promising bacteriolytic activity of vB_EcoS_UTEC10 against MDR E. coli which was further supported by a significant reduction in bacterial load in specimens collected from the phage-treated mice. Histopathology examination demonstrated minimal signs of inflammation and necrosis in the tissues of phage-treated mice compared to the degenerative tissue damage observed in untreated mice. In summary, the present findings suggest that vB_EcoS_UTEC10 has a remarkable ability to eradicate MDR E. coli infections and biofilms. These findings could be further invested for the development of targeted phage therapies that offer a viable alternative to traditional antibiotics against resistant E. coli.

Isolation and characterization of Yersinia phage fMtkYen3-01

Yersinia enterocolitica causes yersiniosis, the third most common gastrointestinal infection in humans throughout Europe. The emergence of multidrug resistance and the lack of effective new antibiotics have drawn attention to phage therapy as a treatment option. Here, we report the complete genome sequence of phage fMtkYen3-01, which infects Y. enterocolitica serotype O:3 strains. This phage has a genome 40,415 bp in length with 45.1% GC content and 49 predicted genes. fMtkYen3-01 infected 9.5% of the 42 Y. enterocolitica strains tested and showed stability at 25-40 °C, as well as pH 5.0-10.0. These results suggest the therapeutic potential of this phage.

Phenotypic characterization and genomic analysis of a Salmonella phage L223 for biocontrol of Salmonella spp. in poultry

The escalating incidence of foodborne salmonellosis poses a significant global threat to food safety and public health. As antibiotic resistance in Salmonella continues to rise, there is growing interest in bacteriophages as potential alternatives. In this study, we isolated, characterized, and evaluated the biocontrol efficacy of lytic phage L223 in chicken meat. Phage L223 demonstrated robust stability across a broad range of temperatures (20-70 °C) and pH levels (2-11) and exhibited a restricted host range targeting Salmonella spp., notably Salmonella Typhimurium and Salmonella Enteritidis. Characterization of L223 revealed a short latent period of 30 min and a substantial burst size of 515 PFU/cell. Genomic analysis classified L223 within the Caudoviricetes class, Guernseyvirinae subfamily and Jerseyvirus genus, with a dsDNA genome size of 44,321 bp and 47.9% GC content, featuring 72 coding sequences devoid of antimicrobial resistance, virulence factors, toxins, and tRNA genes. Application of L223 significantly (p < 0.005) reduced Salmonella Typhimurium ATCC 14,028 counts by 1.24, 2.17, and 1.55 log CFU/piece after 2, 4, and 6 h of incubation, respectively, in experimentally contaminated chicken breast samples. These findings highlight the potential of Salmonella phage L223 as a promising biocontrol agent for mitigating Salmonella contamination in food products, emphasizing its relevance for enhancing food safety protocols.

Diverse Prophage Elements of Salmonella enterica Serovars Show Potential Roles in Bacterial Pathogenicity

Nontyphoidal salmonellosis is an important foodborne and zoonotic infection that causes significant global public health concern. Diverse serovars are multidrug-resistant and encode several virulence indicators; however, little is known on the role prophages play in driving these traits. Here, we extracted prophages from seventy-five Salmonella genomes which represent the fifteen important serovars in the United Kingdom. We analyzed the intact prophages for the presence of virulence genes and established their genomic relationships. We identified 615 prophages from the Salmonella strains, from which 195 prophages are intact, 332 are incomplete, while 88 are questionable. The average prophage carriage was found to be 'extreme' in S. Heidelberg, S. Inverness, and S. Newport (10.2-11.6 prophages/strain), 'high' in S. Infantis, S. Stanley, S. Typhimurium, and S. Virchow (8.2-9.0 prophages/strain), 'moderate' in S. Agona, S. Braenderup, S. Bovismorbificans, S. Choleraesuis, S. Dublin, and S. Java (6.0-7.8 prophages/strain), and 'low' in S. Javiana and S. Enteritidis (5.8 prophages/strain). Cumulatively, 61 virulence genes (1500 gene copies) were detected from representative intact prophages and linked to Salmonella delivery/secretion system (42.62%), adherence (32.7%), magnesium uptake (3.88%), regulation (5%), stress/survival (1.6%), toxins (10%), and antivirulence (1.6%). Diverse clusters were formed among the intact prophages and with bacteriophages of other enterobacteria, suggesting different lineages and associations. Our work provides a strong body of data to support the contributions diverse prophages make to the pathogenicity of Salmonella, including thirteen previously unexplored serovars.

The long and short of it: benchmarking viromics using Illumina, Nanopore and PacBio sequencing technologies

Viral metagenomics has fuelled a rapid change in our understanding of global viral diversity and ecology. Long-read sequencing and hybrid assembly approaches that combine long- and short-read technologies are now being widely implemented in bacterial genomics and metagenomics. However, the use of long-read sequencing to investigate viral communities is still in its infancy. While Nanopore and PacBio technologies have been applied to viral metagenomics, it is not known to what extent different technologies will impact the reconstruction of the viral community. Thus, we constructed a mock bacteriophage community of previously sequenced phage genomes and sequenced them using Illumina, Nanopore and PacBio sequencing technologies and tested a number of different assembly approaches. When using a single sequencing technology, Illumina assemblies were the best at recovering phage genomes. Nanopore- and PacBio-only assemblies performed poorly in comparison to Illumina in both genome recovery and error rates, which both varied with the assembler used. The best Nanopore assembly had errors that manifested as SNPs and INDELs at frequencies 41 and 157 % higher than found in Illumina only assemblies, respectively. While the best PacBio assemblies had SNPs at frequencies 12 and 78 % higher than found in Illumina-only assemblies, respectively. Despite high-read coverage, long-read-only assemblies recovered a maximum of one complete genome from any assembly, unless reads were down-sampled prior to assembly. Overall the best approach was assembly by a combination of Illumina and Nanopore reads, which reduced error rates to levels comparable with short-read-only assemblies. When using a single technology, Illumina only was the best approach. The differences in genome recovery and error rates between technology and assembler had downstream impacts on gene prediction, viral prediction, and subsequent estimates of diversity within a sample. These findings will provide a starting point for others in the choice of reads and assembly algorithms for the analysis of viromes.

Viral metagenomics reveals diverse virus-host interactions throughout the soil depth profile

IMPORTANCE

Soil viruses can moderate the roles that their host microbes play in global carbon cycling. However, given that most studies investigate the surface layer (i.e., top 20 cm) of soil, the extent to which this occurs in subsurface soil (i.e., below 20 cm) is unknown. Here, we leveraged public sequencing data to investigate the interactions between viruses and their hosts at soil depth intervals, down to 115 cm. While most viruses were detected throughout the soil depth profile, their adaptation to host microbes varied. Nonetheless, we uncovered evidence for the potential of soil viruses to encourage their hosts to recycle plant-derived carbon in both surface and subsurface soils. This work reasons that our understanding of soil viral functions requires us to continue to dig deeper and compare viruses existing throughout soil ecosystems.

Comparative genomics and proteomics analysis of phages infecting multi-drug resistant Escherichia coli O177 isolated from cattle faeces

The increasing prevalence of antimicrobial-resistant (AMR) pathogens has become a major global health concern. To address this challenge, innovative strategies such as bacteriophage therapy must be optimised. Genomic characterisation is a crucial step in identifying suitable phage candidates for combating AMR pathogens. The aim of this study was to characterise seven phages that infect the Escherichia coli O177 strain using a whole genome sequencing. The analysis of genome sequences revealed that these phages had linear dsDNA, with genome sizes spanning from 136, 483 to 166,791 bp and GC content varying from 35.39 to 43.63%. Taxonomically, the phages were classified under three different subfamilies (Stephanstirmvirinae, Tevenvirinae, and Vequintavirinae) and three genera (Phapecoctavirus, Tequatrovirus, and Vequintavirus) within the class Caudoviricetes. In silico PhageAI analysis predicted that all the phages were virulent, with confidence levels between 96.07 and 97.26%. The phage genomes contained between 66 and 82 ORFs, which encode hypothetical and putative functional proteins. In addition, the phage genomes contained core genes associated with molecular processes such as DNA replication, transcription modulation, nucleotide metabolism, phage structure (capsid and tail), and lysis. None of the genomes carried genes associated with undesirable traits such as integrase, antimicrobial resistance, virulence, and toxins. The study revealed high genome and proteome homology among E. coli O177 phages and other known Escherichia phages. The results suggest that the seven phages are new members of the genera Phapecoctavirus, Tequatrovirus, and Vequintavirus under the subfamilies Stephanstirmvirinae, Tevenvirinae, and Vequintavirinae, respectively.

Diverse Durham collection phages demonstrate complex BREX defense responses

Bacteriophages (phages) outnumber bacteria ten-to-one and cause infections at a rate of 1025 per second. The ability of phages to reduce bacterial populations makes them attractive alternative antibacterials for use in combating the rise in antimicrobial resistance. This effort may be hindered due to bacterial defenses such as Bacteriophage Exclusion (BREX) that have arisen from the constant evolutionary battle between bacteria and phages. For phages to be widely accepted as therapeutics in Western medicine, more must be understood about bacteria-phage interactions and the outcomes of bacterial phage defense. Here, we present the annotated genomes of 12 novel bacteriophage species isolated from water sources in Durham, UK, during undergraduate practical classes. The collection includes diverse species from across known phylogenetic groups. Comparative analyses of two novel phages from the collection suggest they may be founding members of a new genus. Using this Durham phage collection, we determined that particular BREX defense systems were likely to confer a varied degree of resistance against an invading phage. We concluded that the number of BREX target motifs encoded in the phage genome was not proportional to the degree of susceptibility. IMPORTANCE Bacteriophages have long been the source of tools for biotechnology that are in everyday use in molecular biology research laboratories worldwide. Phages make attractive new targets for the development of novel antimicrobials. While the number of phage genome depositions has increased in recent years, the expected bacteriophage diversity remains underrepresented. Here we demonstrate how undergraduates can contribute to the identification of novel phages and that a single City in England can provide ample phage diversity and the opportunity to find novel technologies. Moreover, we demonstrate that the interactions and intricacies of the interplay between bacterial phage defense systems such as Bacteriophage Exclusion (BREX) and phages are more complex than originally thought. Further work will be required in the field before the dynamic interactions between phages and bacterial defense systems are fully understood and integrated with novel phage therapies.

The Lytic Activity of Bacteriophage ZCSE9 against Salmonella enterica and Its Synergistic Effects with Kanamycin

Salmonella, the causative agent of several diseases in humans and animals, including salmonellosis, septicemia, typhoid fever, and fowl typhoid, poses a serious threat to global public health and food safety. Globally, reports of therapeutic failures are increasing because of the increase in bacterial antibiotic resistance. Thus, this work highlights the combined phage–antibiotic therapy as a promising approach to combating bacterial resistance. In this manner, the phage ZCSE9 was isolated, and the morphology, host infectivity, killing curve, combination with kanamycin, and genome analysis of this phage were all examined. Morphologically, phage ZCSE9 is a siphovirus with a relatively broad host range. In addition, the phage can tolerate high temperatures until 80 °C with one log reduction and a basic environment (pH 11) without a significant decline. Furthermore, the phage prevents bacterial growth in the planktonic state, according to the results of the time-killing curve. Moreover, using the phage at MOI 0.1 with kanamycin against five different Salmonella serotypes reduces the required antibiotics to inhibit the growth of the bacteria. Comparative genomics and phylogenetic analysis suggested that phage ZCSE9, along with its close relatives Salmonella phages vB_SenS_AG11 and wksl3, belongs to the genus Jerseyvirus. In conclusion, phage ZCSE9 and kanamycin form a robust heterologous antibacterial combination that enhances the effectiveness of a phage-only approach for combating Salmonella.

Four principles to establish a universal virus taxonomy

A universal taxonomy of viruses is essential for a comprehensive view of the virus world and for communicating the complicated evolutionary relationships among viruses. However, there are major differences in the conceptualisation and approaches to virus classification and nomenclature among virologists, clinicians, agronomists, and other interested parties. Here, we provide recommendations to guide the construction of a coherent and comprehensive virus taxonomy, based on expert scientific consensus. Firstly, assignments of viruses should be congruent with the best attainable reconstruction of their evolutionary histories, i.e., taxa should be monophyletic. This fundamental principle for classification of viruses is currently included in the International Committee on Taxonomy of Viruses (ICTV) code only for the rank of species. Secondly, phenotypic and ecological properties of viruses may inform, but not override, evolutionary relatedness in the placement of ranks. Thirdly, alternative classifications that consider phenotypic attributes, such as being vector-borne (e.g., "arboviruses"), infecting a certain type of host (e.g., "mycoviruses," "bacteriophages") or displaying specific pathogenicity (e.g., "human immunodeficiency viruses"), may serve important clinical and regulatory purposes but often create polyphyletic categories that do not reflect evolutionary relationships. Nevertheless, such classifications ought to be maintained if they serve the needs of specific communities or play a practical clinical or regulatory role. However, they should not be considered or called taxonomies. Finally, while an evolution-based framework enables viruses discovered by metagenomics to be incorporated into the ICTV taxonomy, there are essential requirements for quality control of the sequence data used for these assignments. Combined, these four principles will enable future development and expansion of virus taxonomy as the true evolutionary diversity of viruses becomes apparent.

Characterization and comprehensive genome analysis of novel bacteriophage, vB_Kpn_ZCKp20p, with lytic and anti-biofilm potential against clinical multidrug-resistant Klebsiella pneumoniae

Introduction

The rise of infections by antibiotic-resistant bacterial pathogens is alarming. Among these, Klebsiella pneumoniae is a leading cause of death by hospital-acquired infections, and its multidrug-resistant strains are flagged as a global threat to human health, which necessitates finding novel antibiotics or alternative therapies. One promising therapeutic alternative is the use of virulent bacteriophages, which specifically target bacteria and coevolve with them to overcome potential resistance. Here, we aimed to discover specific bacteriophages with therapeutic potential against multiresistant K. pneumoniae clinical isolates.

Methods and Results

Out of six bacteriophages that we isolated from urban and medical sewage, phage vB_Kpn_ZCKp20p had the broadest host range and was thus characterized in detail. Transmission electron microscopy suggests vB_Kpn_ZCKp20p to be a tailed phage of the siphoviral morphotype. In vitro evaluation indicated a high lytic efficiency (30 min latent period and burst size of ∼100 PFU/cell), and extended stability at temperatures up to 70°C and a wide range of (2-12) pH. Additionally, phage vB_Kpn_ZCKp20p possesses antibiofilm activity that was evaluated by the crystal violet assay and was not cytotoxic to human skin fibroblasts. The whole genome was sequenced and annotated, uncovering one tRNA gene and 33 genes encoding proteins with assigned functions out of 85 predicted genes. Furthermore, comparative genomics and phylogenetic analysis suggest that vB_Kpn_ZCKp20p most likely represents a new species, but belongs to the same genus as Klebsiella phages ZCKP8 and 6691. Comprehensive genomic and bioinformatics analyses substantiate the safety of the phage and its strictly lytic lifestyle.

Conclusion

Phage vB_Kpn_ZCKp20p is a novel phage with potential to be used against biofilm-forming K. pneumoniae and could be a promising source for antibacterial and antibiofilm products, which will be individually studied experimentally in future studies.

Diversity, Dynamics and Therapeutic Application of Clostridioides difficile Bacteriophages

Clostridioides difficile causes antibiotic-induced diarrhoea and pseudomembranous colitis in humans and animals. Current conventional treatment relies solely on antibiotics, but C. difficile infection (CDI) cases remain persistently high with concomitant increased recurrence often due to the emergence of antibiotic-resistant strains. Antibiotics used in treatment also induce gut microbial imbalance; therefore, novel therapeutics with improved target specificity are being investigated. Bacteriophages (phages) kill bacteria with precision, hence are alternative therapeutics for the targeted eradication of the pathogen. Here, we review current progress in C. difficile phage research. We discuss tested strategies of isolating C. difficile phages directly, and via enrichment methods from various sample types and through antibiotic induction to mediate prophage release. We also summarise phenotypic phage data that reveal their morphological, genetic diversity, and various ways they impact their host physiology and pathogenicity during infection and lysogeny. Furthermore, we describe the therapeutic development of phages through efficacy testing in different in vitro, ex vivo and in vivo infection models. We also discuss genetic modification of phages to prevent horizontal gene transfer and improve lysis efficacy and formulation to enhance stability and delivery of the phages. The goal of this review is to provide a more in-depth understanding of C. difficile phages and theoretical and practical knowledge on pre-clinical, therapeutic evaluation of the safety and effectiveness of phage therapy for CDI.

Extended genomic analyses of the broad-host-range phages vB_KmiM-2Di and vB_KmiM-4Dii reveal slopekviruses have highly conserved genomes

High levels of antimicrobial resistance among members of the Klebsiella oxytoca complex (KoC) have led to renewed interest in the use of bacteriophage (phage) therapy to tackle infections caused by these bacteria. In this study we characterized two lytic phages, vB_KmiM-2Di and vB_KmiM-4Dii, that were isolated from sewage water against two GES-5-positive Klebsiella michiganensis strains (PS_Koxy2 and PS_Koxy4, respectively). ViPTree analysis showed both phages belonged to the genus Slopekvirus. rpoB gene-based sequence analysis of 108 presumptive K. oxytoca isolates (n=59 clinical, n=49 veterinary) found K. michiganensis to be more prevalent (46 % clinical and 43 % veterinary, respectively) than K. oxytoca (40 % clinical and 6 % veterinary, respectively). Host range analysis against these 108 isolates found both vB_KmiM-2Di and vB_KmiM-4Dii showed broad lytic activity against KoC species. Several hypothetical homing endonuclease genes were encoded within the genomes of both phages, which may contribute to their broad host range. Differences in the tail fibre protein may explain the non-identical host range of the two phages. Pangenome analysis of 24 slopekviruses found that genomes within this genus are highly conserved, with more than 50 % of all predicted coding sequences representing core genes at ≥95 % identity and ≥70 % coverage. Given their broad host ranges, our results suggest vB_KmiM-2Di and vB_KmiM-4Dii represent attractive potential therapeutics. In addition, current recommendations for phage-based pangenome analyses may require revision.