PHROG: families of prokaryotic virus proteins clustered using remote homology

phage therapy candidates from Sphae: An automated toolkit for predicting sequencing data

Motivation

Phage therapy is a viable alternative for treating bacterial infections amidst the escalating threat of antimicrobial resistance. However, the therapeutic success of phage therapy depends on selecting safe and effective phage candidates. While experimental methods focus on isolating phages and determining their lifecycle and host range, comprehensive genomic screening is critical to identify markers that indicate potential risks, such as toxins, antimicrobial resistance, or temperate lifecycle traits. These analyses are often labor-intensive and time-consuming, limiting the rapid deployment of phage in clinical settings.

Results

We developed Sphae, an automated bioinformatics pipeline designed to streamline therapeutic potential of a phage in under ten minutes. Using Snakemake workflow manager, Sphae integrates tools for quality control, assembly, genome assessment, and annotation tailored specifically for phage biology. Sphae automates the detection of key genomic markers, including virulence factors, antimicrobial resistance genes, and lysogeny indicators like integrase, recombinase, and transposase, which could preclude therapeutic use. Benchmarked on 65 phage sequences, 28 phage samples showed therapeutic potential, 8 failed during assembly due to low sequencing depth, 22 samples included prophage or virulent markers, and the remaining 23 samples included multiple phage genomes per sample. This workflow outputs a comprehensive report, enabling rapid assessment of phage safety and suitability for phage therapy under these criteria. Sphae is scalable, portable, facilitating efficient deployment across most high-performance computing (HPC) and cloud platforms, expediting the genomic evaluation process.

Availability

Sphae is source code and freely available at https://github.com/linsalrob/sphae, with installation supported on Conda, PyPi, Docker containers.

Isolation and Characterization of Infection of Four New Bacteriophages Infecting a Vibrio parahaemolyticus Strain

Bacteria from genus Vibrio continue to be one of the most common threats to aquaculture sustainability. Vibrio spp. have been associated with infectious outbreaks in fish, shrimp, bivalves and even algae farms worldwide. Moreover, several Vibrio spp. are also pathogens that impact human health and are a threat to public health when transferred to consumers through contaminated seafood products. The use of bacteriophages is an evolving technology that could be applied in the treatment of Vibrio spp. either to protect aquaculture farms or to decontaminate seafood, namely bivalves during their depuration. In the present study, bacteriophages vB_VpS_LMAVpS1 (S1) vB_VpS_LMAVpVPP (VPP), vB_VpS_LMAVpSH (SH) and vB_VpS_LMAVpH (H) infecting V. parahaemolyticus were isolated and characterized. All phages presented fast adsorption rates and were able to control V. parahaemolyticus at all multiplicity of infections (MOIs) tested (MOI of 1, 10 and 100), with reductions of more than 4 log CFU/mL being recorded, but only in the presence of divalent cation calcium. The rate of emergence of phage-resistant mutants was very low (1.8 × 10-6 to 3.1 × 10-6). Bacterial phage resistance was not permanent and led to a loss of bacterial fitness. All four phages presented with lysins encoded in their genomes. The results presented provide valuable insights for future studies in the application of these bacteriophages in different scenarios to control, decontaminate or treat bacterial infections or contaminations of V. parahaemolyticus.

Complete genome sequences of three Pseudomonas aeruginosa jumbo bacteriophages discovered in Kenya

The genomes of three Pseudomonas aeruginosa Phikzvirus bacteriophages isolated in Kenya are described. The genomes of phages vB_PaePAO1-KEN19, vB_Pae3705-KEN49, and vB_Pae10145-KEN51, respectively, had lengths of 278,921, 280,231, and 280,173 bp, with 36.93%, 36.84%, and 36.86% GC content, containing 419, 417, and 417 coding sequences (including seven tRNAs in each genome).

Pipolins are bimodular platforms that maintain a reservoir of defense systems exchangeable with various bacterial genetic mobile elements

Defense genes gather in diverse types of genomic islands in bacteria and provide immunity against viruses and other genetic mobile elements. Here, we disclose pipolins, previously found in diverse bacterial phyla and encoding a primer-independent PolB, as a new category of widespread defense islands. The analysis of the occurrence and structure of pipolins revealed that they are commonly integrative elements flanked by direct repeats in Gammaproteobacteria genomes, mainly Escherichia, Vibrio or Aeromonas, often taking up known mobile elements integration hotspots. Remarkably, integrase dynamics correlates with alternative integration spots and enables diverse lifestyles, from integrative to mobilizable and plasmid pipolins, such as in members of the genera Limosilactobacillus, Pseudosulfitobacter or Staphylococcus. Pipolins harbor a minimal core and a large cargo module enriched for defense factors. In addition, analysis of the weighted gene repertoire relatedness revealed that many of these defense factors are actively exchanged with other mobile elements. These findings indicate pipolins and, potentially other defense islands, act as orthogonal reservoirs of defense genes, potentially transferable to immune autonomous MGEs, suggesting complementary exchange mechanisms for defense genes in bacterial populations.

Isolation and characterization of the new Streptomyces phages Kamino, Geonosis, Abafar, and Scarif infecting a broad range of host species

Streptomyces, a multifaceted genus of soil-dwelling bacteria belonging to the phylum Actinomycetota, features intricate phage-host interactions shaped by its complex life cycle and the synthesis of a diverse array of specialized metabolites. Here, we describe the isolation and characterization of four novel Streptomyces phages infecting a variety of different host species. While phage Kamino, isolated on Streptomyces kasugaensis, is predicted to be temperate and encodes a serine integrase in its genome, phages Geonosis (isolated on Streptomyces griseus) and Abafar and Scarif, isolated on Streptomyces albidoflavus, are virulent phages. Phages Kamino and Geonosis were shown to amplify well in liquid culture leading to a pronounced culture collapse already at low titers. Determination of the host range by testing >40 different Streptomyces species identified phages Kamino, Abafar, and Scarif as broad host-range phages. Overall, the phages described in this study expand the publicly available portfolio of phages infecting Streptomyces and will be instrumental in advancing the mechanistic understanding of the intricate antiviral strategies employed by these multicellular bacteria.IMPORTANCEThe actinobacterial genus Streptomyces is characterized by multicellular, filamentous growth and the synthesis of a diverse range of bioactive molecules. These characteristics also play a role in shaping their interactions with the most abundant predator in the environment, bacteriophages-viruses infecting bacteria. In this study, we characterize four new phages infecting Streptomyces. Out of those, three phages feature a broad host range infecting up to 15 different species. The isolated phages were characterized with respect to plaque and virion morphology, host range, and amplification in liquid culture. In summary, the phages reported in this study contribute to the broader collection of publicly available phages infecting Streptomyces, playing a crucial role in advancing our mechanistic understanding of phage-host interactions of these multicellular bacteria.

A long-context language model for deciphering and generating bacteriophage genomes

Inspired by the success of large language models (LLMs), we develop a long-context generative model for genomes. Our multiscale transformer model, megaDNA, is pre-trained on unannotated bacteriophage genomes with nucleotide-level tokenization. We demonstrate the foundational capabilities of our model including the prediction of essential genes, genetic variant effects, regulatory element activity and taxonomy of unannotated sequences. Furthermore, it generates de novo sequences up to 96 K base pairs, which contain potential regulatory elements and annotated proteins with phage-related functions.

V- and VL-Scores Uncover Viral Signatures and Origins of Protein Families

Viruses are key drivers of microbial diversity, nutrient cycling, and co-evolution in ecosystems, yet their study is hindered due to challenges in culturing. Traditional gene-centric methods, which focus on a few hallmark genes like for capsids, miss much of the viral genome, leaving key viral proteins and functions undiscovered. Here, we introduce two powerful annotation-free metrics, V-score and VL-score, designed to quantify the "virus-likeness" of protein families and genomes and create an open-access searchable database, 'V-Score-Search'. By applying V- and VL-scores to public databases (KEGG, Pfam, and eggNOG), we link 38-77% of protein families with viruses, a 9-16x increase over current estimates. These metrics outperform existing approaches, enabling precise detection of viral genomes, prophages, and host-derived auxiliary viral genes (AVGs) from fragmented sequences, and significantly improving genome binning. Remarkably, we identify up to 17x more AVGs, dominated by non-metabolic proteins of unknown function. This innovation unlocks new insights into virus signatures and host interactions, with wide-ranging implications from genomics to biotechnology.

Complete genomes of Asgard archaea reveal diverse integrated and mobile genetic elements

Asgard archaea are of great interest as the progenitors of Eukaryotes, but little is known about the mobile genetic elements (MGEs) that may shape their ongoing evolution. Here, we describe MGEs that replicate in Atabeyarchaeia, a wetland Asgard archaea lineage represented by two complete genomes. We used soil depth-resolved population metagenomic data sets to track 18 MGEs for which genome structures were defined and precise chromosome integration sites could be identified for confident host linkage. Additionally, we identified a complete 20.67 kbp circular plasmid and two family-level groups of viruses linked to Atabeyarchaeia, via CRISPR spacer targeting. Closely related 40 kbp viruses possess a hypervariable genomic region encoding combinations of specific genes for small cysteine-rich proteins structurally similar to restriction-homing endonucleases. One 10.9 kbp integrative conjugative element (ICE) integrates genomically into the Atabeyarchaeum deiterrae-1 chromosome and has a 2.5 kbp circularizable element integrated within it. The 10.9 kbp ICE encodes an expressed Type IIG restriction-modification system with a sequence specificity matching an active methylation motif identified by Pacific Biosciences (PacBio) high-accuracy long-read (HiFi) metagenomic sequencing. Restriction-modification of Atabeyarchaeia differs from that of another coexisting Asgard archaea, Freyarchaeia, which has few identified MGEs but possesses diverse defense mechanisms, including DISARM and Hachiman, not found in Atabeyarchaeia. Overall, defense systems and methylation mechanisms of Asgard archaea likely modulate their interactions with MGEs, and integration/excision and copy number variation of MGEs in turn enable host genetic versatility.

MVP: a modular viromics pipeline to identify, filter, cluster, annotate, and bin viruses from metagenomes

While numerous computational frameworks and workflows are available for recovering prokaryote and eukaryote genomes from metagenome data, only a limited number of pipelines are designed specifically for viromics analysis. With many viromics tools developed in the last few years alone, it can be challenging for scientists with limited bioinformatics experience to easily recover, evaluate quality, annotate genes, dereplicate, assign taxonomy, and calculate relative abundance and coverage of viral genomes using state-of-the-art methods and standards. Here, we describe Modular Viromics Pipeline (MVP) v.1.0, a user-friendly pipeline written in Python and providing a simple framework to perform standard viromics analyses. MVP combines multiple tools to enable viral genome identification, characterization of genome quality, filtering, clustering, taxonomic and functional annotation, genome binning, and comprehensive summaries of results that can be used for downstream ecological analyses. Overall, MVP provides a standardized and reproducible pipeline for both extensive and robust characterization of viruses from large-scale sequencing data including metagenomes, metatranscriptomes, viromes, and isolate genomes. As a typical use case, we show how the entire MVP pipeline can be applied to a set of 20 metagenomes from wetland sediments using only 10 modules executed via command lines, leading to the identification of 11,656 viral contigs and 8,145 viral operational taxonomic units (vOTUs) displaying a clear beta-diversity pattern. Further, acting as a dynamic wrapper, MVP is designed to continuously incorporate updates and integrate new tools, ensuring its ongoing relevance in the rapidly evolving field of viromics. MVP is available at https://gitlab.com/ccoclet/mvp and as versioned packages in PyPi and Conda.IMPORTANCEThe significance of our work lies in the development of Modular Viromics Pipeline (MVP), an integrated and user-friendly pipeline tailored exclusively for viromics analyses. MVP stands out due to its modular design, which ensures easy installation, execution, and integration of new tools and databases. By combining state-of-the-art tools such as geNomad and CheckV, MVP provides high-quality viral genome recovery and taxonomy and host assignment, and functional annotation, addressing the limitations of existing pipelines. MVP's ability to handle diverse sample types, including environmental, human microbiome, and plant-associated samples, makes it a versatile tool for the broader microbiome research community. By standardizing the analysis process and providing easily interpretable results, MVP enables researchers to perform comprehensive studies of viral communities, significantly advancing our understanding of viral ecology and its impact on various ecosystems.

Improving viral annotation with artificial intelligence

Viruses of bacteria, "phages," are fundamental, poorly understood components of microbial community structure and function. Additionally, their dependence on hosts for replication positions phages as unique sensors of ecosystem features and environmental pressures. High-throughput sequencing approaches have begun to give us access to the diversity and range of phage populations in complex microbial community samples, and metagenomics is currently the primary tool with which we study phage populations. The study of phages by metagenomic sequencing, however, is fundamentally limited by viral diversity, which results in the vast majority of viral genomes and metagenome-annotated genomes lacking annotation. To harness bacteriophages for applications in human and environmental health and disease, we need new methods to organize and annotate viral sequence diversity. We recently demonstrated that methods that leverage self-supervised representation learning can supplement statistical sequence representations for remote viral protein homology detection in the ocean virome and propose that consideration of the functional content of viral sequences allows for the identification of similarity in otherwise sequence-diverse viruses and viral-like elements for biological discovery. In this review, we describe the potential and pitfalls of large language models for viral annotation. We describe the need for new approaches to annotate viral sequences in metagenomes, the fundamentals of what protein language models are and how one can use them for sequence annotation, the strengths and weaknesses of these models, and future directions toward developing better models for viral annotation more broadly.

Viromes vs. mixed community metagenomes: choice of method dictates interpretation of viral community ecology

BACKGROUND

Viruses, the majority of which are uncultivated, are among the most abundant biological entities on Earth. From altering microbial physiology to driving community dynamics, viruses are fundamental members of microbiomes. While the number of studies leveraging viral metagenomics (viromics) for studying uncultivated viruses is growing, standards for viromics research are lacking. Viromics can utilize computational discovery of viruses from total metagenomes of all community members (hereafter metagenomes) or use physical separation of virus-specific fractions (hereafter viromes). However, differences in the recovery and interpretation of viruses from metagenomes and viromes obtained from the same samples remain understudied.

RESULTS

Here, we compare viral communities from paired viromes and metagenomes obtained from 60 diverse samples across human gut, soil, freshwater, and marine ecosystems. Overall, viral communities obtained from viromes had greater species richness and total viral genome abundances than those obtained from metagenomes, although there were some exceptions. Despite this, metagenomes still contained many viral genomes not detected in viromes. We also found notable differences in the predicted lytic state of viruses detected in viromes vs metagenomes at the time of sequencing. Other forms of variation observed include genome presence/absence, genome quality, and encoded protein content between viromes and metagenomes, but the magnitude of these differences varied by environment.

CONCLUSIONS

Overall, our results show that the choice of method can lead to differing interpretations of viral community ecology. We suggest that the choice of whether to target a metagenome or virome to study viral communities should be dependent on the environmental context and ecological questions being asked. However, our overall recommendation to researchers investigating viral ecology and evolution is to pair both approaches to maximize their respective benefits. Video Abstract.

Differentiation of the Anammox core microbiome: Unraveling the evolutionary impetus of scalable gene flow

Anaerobic ammonium oxidation bacteria (AAOB), distinguished by their unique autotrophic nitrogen metabolism, hold pivotal positions in the global nitrogen cycle and environmental biotechnologies. However, the ecophysiology and evolution of AAOB remain poorly understood, attributed to the absence of monocultures. Hence, a comprehensive elucidation of the AAOB-dominated core microbiome, anammox core, is imperative to further completing the theory of engineered nitrogen removal and ecological roles of anammox. Performing taxonomic and phylogenetic analyses on collected genome repertoires, we show here that Candidatus Brocadia and Candidatus Kuenenia possesses a more compact core than Candidatus Jettenia, which partly explains why the latter has a less common ecological presence. Evidence of gene flow is particularly striking in functions related to biosynthesis and oxygen detoxification, underscoring the evolutionary forces driving lineage and core differentiation. Furthermore, CRISPR spacer traceback of the AAOB metagenome-assembled genomes (MAGs) reveals a series of genetic traces for the concealed phages. By reconceptualizing the functional divergence of AAOB with the historical role of phages, we ultimately propose a coevolutionary framework to understand the evolutionary trajectory of anammox microecology. The discoveries provided in this study offer new insights into understanding the evolution of AAOB and the ecology of anammox.

Targeted phage hunting to specific Klebsiella pneumoniae clinical isolates is an efficient antibiotic resistance and infection control strategy

Klebsiella pneumoniae is one of the most threatening multi-drug-resistant pathogens today, with phage therapy being a promising alternative for personalized treatments. However, the intrinsic capsule diversity in Klebsiella spp. poses a substantial barrier to the phage host range, complicating the development of broad-spectrum phage-based treatments. Here, we have isolated and genomically characterized phages capable of infecting each of the acquired 77 reference serotypes of Klebsiella spp., including capsular types widespread among high-risk K. pneumoniae clones causing nosocomial infections. We demonstrated the possibility of isolating phages for all capsular types in the collection, revealing high capsular specificity among taxonomically related phages, in contrast to a few phages that exhibited broad-spectrum infection capabilities. To decipher the determinants of the specificity of these phages, we focused on their receptor-binding proteins, with particular attention to depolymerases. We also explored the possibility of designing a broad-spectrum phage cocktail based on phages isolated in reference capsular-type strains and determining the ability to lyse relevant clinical isolates. A combination of 12 phages capable of infecting 55% of the reference Klebsiella spp. serotypes was tested on a panel of carbapenem-resistant K. pneumoniae clinical isolates. Thirty-one percent of isolates were susceptible to the phage cocktail. However, our results suggest that in a highly variable encapsulated bacterial host, phage hunting must be directed to the specific Klebsiella isolates. This work is a step forward in the understanding of the complexity of phage-host interactions and highlights the importance of implementing precise and phage-specific strategies to treat K. pneumoniae infections worldwide.IMPORTANCEThe emergence of resistant bacteria is a serious global health problem. In the absence of effective treatments, phages are a personalized and effective therapeutic alternative. However, little is still known about phage-host interactions, which are key to implementing effective strategies. Here, we focus on the study of Klebsiella pneumoniae, a highly pathogenic encapsulated bacterium. The complexity and variability of the capsule, where in most cases phage receptors are found, make it difficult for phage-based treatments. Here, we isolated a large collection of Klebsiella phages against all the reference strains and in a cohort of clinical isolates. Our results suggest that clinical isolates represent a challenge, especially high-risk clones. Thus, we propose targeted phage hunting as an effective strategy to implement phage-derived therapies. Our results are a step forward for new phage-based strategies to control K. pneumoniae infections, highlighting the importance of understanding phage-host interactions to design personalized treatments against Klebsiella spp.

C. Rocha E. Intragenomic conflicts with plasmids and chromosomal mobile genetic elements drive the evolution of natural transformation within species

Natural transformation is the only mechanism of genetic exchange controlled by the recipient bacteria. We quantified its rates in 786 clinical strains of the human pathogens Legionella pneumophila (Lp) and 496 clinical and environmental strains of Acinetobacter baumannii (Ab). The analysis of transformation rates in the light of phylogeny revealed they evolve by a mixture of frequent small changes and a few large quick jumps across 6 orders of magnitude. In standard conditions close to half of the strains of Lp and a more than a third in Ab are below the detection limit and thus presumably non-transformable. Ab environmental strains tend to have higher transformation rates than the clinical ones. Transitions to non-transformability were frequent and usually recent, suggesting that they are deleterious and subsequently purged by natural selection. Accordingly, we find that transformation decreases genetic linkage in both species, which might accelerate adaptation. Intragenomic conflicts with chromosomal mobile genetic elements (MGEs) and plasmids could explain these transitions and a GWAS confirmed systematic negative associations between transformation and MGEs: plasmids and other conjugative elements in Lp, prophages in Ab, and transposable elements in both. In accordance with the hypothesis of modulation of transformation rates by genetic conflicts, transformable strains have fewer MGEs in both species and some MGEs inactivate genes implicated in the transformation with heterologous DNA (in Ab). Innate defense systems against MGEs are associated with lower transformation rates, especially restriction-modification systems. In contrast, CRISPR-Cas systems are associated with higher transformation rates suggesting that adaptive defense systems may facilitate cell protection from MGEs while preserving genetic exchanges by natural transformation. Ab and Lp have different lifestyles, gene repertoires, and population structure. Nevertheless, they exhibit similar trends in terms of variation of transformation rates and its determinants, suggesting that genetic conflicts could drive the evolution of natural transformation in many bacteria.

Protein Set Transformer: A protein-based genome language model to power high diversity viromics

Exponential increases in microbial and viral genomic data demand transformational advances in scalable, generalizable frameworks for their interpretation. Standard homology-based functional analyses are hindered by the rapid divergence of microbial and especially viral genomes and proteins that significantly decreases the volume of usable data. Here, we present Protein Set Transformer (PST), a protein-based genome language model that models genomes as sets of proteins without considering sparsely available functional labels. Trained on >100k viruses, PST outperformed other homology- and language model-based approaches for relating viral genomes based on shared protein content. Further, PST demonstrated protein structural and functional awareness by clustering capsid-fold-containing proteins with known capsid proteins and uniquely clustering late gene proteins within related viruses. Our data establish PST as a valuable method for diverse viral genomics, ecology, and evolutionary applications. We posit that the PST framework can be a foundation model for microbial genomics when trained on suitable data.

PhageScanner: a reconfigurable machine learning framework for bacteriophage genomic and metagenomic feature annotation

Bacteriophages are the most prolific organisms on Earth, yet many of their genomes and assemblies from metagenomic sources lack protein sequences with identified functions. While most bacteriophage proteins are structural proteins, categorized as Phage Virion Proteins (PVPs), a considerable number remain unclassified. Complicating matters further, traditional lab-based methods for PVP identification can be tedious. To expedite the process of identifying PVPs, machine-learning models are increasingly being employed. Existing tools have developed models for predicting PVPs from protein sequences as input. However, none of these efforts have built software allowing for both genomic and metagenomic data as input. In addition, there is currently no framework available for easily curating data and creating new types of machine learning models. In response, we introduce PhageScanner, an open-source platform that streamlines data collection for genomic and metagenomic datasets, model training and testing, and includes a prediction pipeline for annotating genomic and metagenomic data. PhageScanner also features a graphical user interface (GUI) for visualizing annotations on genomic and metagenomic data. We further introduce a BLAST-based classifier that outperforms ML-based models and an efficient Long Short-Term Memory (LSTM) classifier. We then showcase the capabilities of PhageScanner by predicting PVPs in six previously uncharacterized bacteriophage genomes. In addition, we create a new model that predicts phage-encoded toxins within bacteriophage genomes, thus displaying the utility of the framework.

Isolation and characterization of Septuagintavirus; a novel clade of Escherichia coli phages within the subfamily Vequintavirinae

Escherichia coli is a commensal inhabitant of the mammalian gut microbiota, frequently associated with various gastrointestinal diseases. There is increasing interest in comprehending the variety of bacteriophages (phages) that target this bacterium, as such insights could pave the way for their potential use in therapeutic applications. Here, we report the isolation and characterization of four newly identified E. coli infecting tailed phages (W70, A7-1, A5-4, and A73) that were found to constitute a novel genus, Septuagintavirus, within the subfamily Vequintavirinae. Genomes of these phages ranged from 137 kbp to 145 kbp, with a GC content of 41 mol%. They possess a maximum nucleotide similarity of 30% with phages of the closest phylogenetic genus, Certrevirus, while displaying limited homology to other genera of the Vequintavirinae family. Host range analysis showed that these phages have limited activity against a panel of E. coli strains, infecting 6 out of 16 tested isolates, regardless of their phylotype. Electrospray ionization-tandem mass spectrometry (ESI-MS/MS) was performed on the virion of phage W70, allowing the identification of 28 structural proteins, 19 of which were shared with phages of other genera of Vequintavirinae family. The greatest diversity was identified with proteins forming tail fiber structures, likely indicating the adaptation of virions of each phage genus of this subfamily for the recognition of their target receptor on host cells. The findings of this study provide greater insights into the phages of the subfamily Vequintavirinae, contributing to the pool of knowledge currently known about these phages.

IMPORTANCE

Escherichia coli is a well-known bacterium that inhabits diverse ecological niches, including the mammalian gut microbiota. Certain strains are associated with gastrointestinal diseases, and there is a growing interest in using bacteriophages, viruses that infect bacteria, to combat bacterial infections. Here, we describe the isolation and characterization of four novel E. coli bacteriophages that constitute a new genus, Septuagintavirus, within the subfamily Vequintavirinae. We conducted mass spectrometry on virions of a representative phage of this novel clade and compared it to other phages within the subfamily. Our analysis shows that virion structure is highly conserved among all phages, except for proteins related to tail fiber structures implicated in the host range. These findings provide greater insights into the phages of the subfamily Vequintavirinae, contributing to the existing pool of knowledge about these phages.

Expanding the Phage Galaxy: Isolation and Characterization of Five Novel Streptomyces Siphoviruses Ankus, Byblos, DekoNeimoidia, Mandalore, and Naboo

Background

Key features of the actinobacterial genus Streptomyces are multicellular, filamentous growth, and production of a broad portfolio of bioactive molecules. These characteristics appear to play an important role in phage-host interactions and are modulated by phages during infection. To accelerate research of such interactions and the investigation of novel immune systems in multicellular bacteria, phage isolation, sequencing, and characterization are needed. This is a prerequisite for establishing systematic collections that appropriately cover phage diversity for comparative analyses.

Material & Methods

As part of a public outreach program within the priority program SPP 2330, involving local schools, we describe the isolation and characterization of five novel Streptomyces siphoviruses infecting S. griseus, S. venezuelae, and S. olivaceus.

Results

All isolates are virulent members of two existing genera and, additionally, establish a new genus in the Stanwilliamsviridae family. In addition to an extensive set of tRNAs and proteins involved in phage replication, about 80% of phage genes encode hypothetical proteins, underlining the yet underexplored phage diversity and genomic dark matter still found in bacteriophages infecting actinobacteria.

Conclusions

Taken together, phages Ankus, Byblos, DekoNeimoidia, Mandalore, and Naboo expand the phage diversity and contribute to ongoing research in the field of Streptomyces phage-host interactions.

Discovery and characterisation of new phage targeting uropathogenic Escherichia coli

Antimicrobial resistance is an escalating threat with few new therapeutic options in the pipeline. Urinary tract infections (UTIs) are one of the most prevalent bacterial infections globally and are prone to becoming recurrent and antibiotic resistant. We discovered and characterized six novel Autographiviridae and Guernseyvirinae bacterial viruses (phage) against uropathogenic Escherichia coli (UPEC), a leading cause of UTIs. The phage genomes were between 39,471 bp - 45,233 bp, with 45.0%-51.0% GC%, and 57-84 predicted coding sequences per genome. We show that tail fiber domain structure, predicted host capsule type, and host antiphage repertoire correlate with phage host range. In vitro characterisation of phage cocktails showed synergistic improvement against a mixed UPEC strain population and when sequentially dosed. Together, these phage are a new set extending available treatments for UTI from UPEC, and phage vM_EcoM_SHAK9454 represents a promising candidate for further improvement through engineering.

A novel genus of Pectobacterium bacteriophages display broad host range by targeting several species of Danish soft rot isolates

The bacterial diseases black leg and soft rot in potatoes cause heavy losses of potatoes worldwide. Bacteria within the genus Pectobacteriaceae are the causative agents of black leg and soft rot. The use of antibiotics in agriculture is heavily regulated and no other effective treatment currently exists, but bacteriophages (phages) have shown promise as potential biocontrol agents. In this study we isolated soft rot bacteria from potato tubers and plant tissue displaying soft rot or black leg symptoms collected in Danish fields. We then used the isolated bacterial strains as hosts for phage isolation. Using organic waste, we isolated phages targeting different species within Pectobacterium. Here we focus on seven of these phages representing a new genus primarily targeting P. brasiliense; phage Ymer, Amona, Sabo, Abuela, Koroua, Taid and Pappous. TEM image of phage Ymer showed siphovirus morphotype, and the proposed Ymer genus belongs to the class Caudoviricetes, with double-stranded DNA genomes varying from 39 kb to 43 kb. In silico host range prediction using a CRISPR-Cas spacer database suggested both P. brasiliense, P. polaris and P. versatile as natural hosts for phages within the proposed Ymer genus. A following host range experiment, using 47 bacterial isolates from Danish tubers and plants symptomatic with soft rot or black leg disease verified the in silico host range prediction, as the genus as a group were able to infect all three Pectobacterium species. Phages did, however, primarily target P. brasiliense isolates and displayed differences in host range even within the species level. Two of the phages were able to infect two or more Pectobacterium species. Despite no nucleotide similarity with any phages in the NCBI database, the proposed Ymer genus did share some similarity at the protein level, as well as gene synteny, with currently known phages. None of the phages encoded integrases or other genes typically associated with lysogeny. Similarly, no virulence factors nor antimicrobial resistance genes were found, and combined with their ability to infect several soft rot-causing Pectobacterium species from Danish fields, demonstrates their potential as biocontrol agents against soft rot and black leg diseases in potatoes.

Birth of protein folds and functions in the virome

The rapid evolution of viruses generates proteins that are essential for infectivity and replication but with unknown functions, due to extreme sequence divergence1. Here, using a database of 67,715 newly predicted protein structures from 4,463 eukaryotic viral species, we found that 62% of viral proteins are structurally distinct and lack homologues in the AlphaFold database2,3. Among the remaining 38% of viral proteins, many have non-viral structural analogues that revealed surprising similarities between human pathogens and their eukaryotic hosts. Structural comparisons suggested putative functions for up to 25% of unannotated viral proteins, including those with roles in the evasion of innate immunity. In particular, RNA ligase T-like phosphodiesterases were found to resemble phage-encoded proteins that hydrolyse the host immune-activating cyclic dinucleotides 3',3'- and 2',3'-cyclic GMP-AMP (cGAMP). Experimental analysis showed that RNA ligase T homologues encoded by avian poxviruses similarly hydrolyse cGAMP, showing that RNA ligase T-mediated targeting of cGAMP is an evolutionarily conserved mechanism of immune evasion that is present in both bacteriophage and eukaryotic viruses. Together, the viral protein structural database and analyses presented here afford new opportunities to identify mechanisms of virus-host interactions that are common across the virome.

The Potential of Phage Treatment to Inactivate Planktonic and Biofilm-Forming Pseudomonas aeruginosa

Pseudomonas aeruginosa is a common cause of hospital-acquired infections and exhibits a strong resistance to antibiotics. An alternative treatment option for bacterial infections is the use of bacteriophages (or phages). In this study, two distinct phages, VB_PaD_phPA-G (phPA-G) and VB_PaN_phPA-Intesti (phPA-Intesti), were used as single suspensions or in a phage cocktail to inactivate the planktonic cells and biofilms of P. aeruginosa. Preliminary experiments in culture medium showed that phage phPA-Intesti (reductions of 4.5-4.9 log CFU/mL) outperformed phPA-G (reductions of 0.6-2.6 log CFU/mL) and the phage cocktail (reduction of 4.2 log CFU/mL). Phage phPA-Intesti caused a maximum reduction of 5.5 log CFU/cm2 in the P. aeruginosa biofilm in urine after 4 h of incubation. The combination of phage phPA-Intesti and ciprofloxacin did not improve the efficacy of bacterial inactivation nor reduce the development of resistant mutants. However, the development of resistant bacteria was lower in the combined treatment with the phage and the antibiotic compared to treatment with the antibiotic alone. This phage lacks known toxins, virulence, antibiotic resistance, and integrase genes. Overall, the results suggest that the use of phage phPA-Intesti could be a potential approach to control urinary tract infections (UTIs), namely those caused by biofilm-producing and multidrug-resistant strains of P. aeruginosa.

Isolation and characterization of two novel bacteriophages against carbapenem-resistant Klebsiella pneumoniae

The increase of antibiotic-resistant bacteria has become a global health emergency and the need to explore alternative therapeutic options arises. Phage therapy uses bacteriophages to target specific bacterial strains. Phages are highly specific and can target resistant bacteria. Currently, research in this regard is focused on ensuring reliability and safety to bring this tool into clinical practice. The first step is to conduct comprehensive preclinical research. In this work, we present two novel bacteriophages vB_Kpn_F13 and vB_Kpn_F14 isolated against clinical carbapenem-resistant Klebsiella pneumoniae strains obtained from hospital sewage. Multiple studies in vitro were conducted, such as sequencing, electron microscopy, stability, host range infectivity, planktonic effect and biofilm inhibition in order to discover their ability to be used against carbapenem-resistant K. pneumoniae pathogens causing difficult-to-treat infections.

Complete genome sequence of Escherichia coli O157:H7 phage Φ241

We report the genome sequence of phage Φ241 infecting Escherichia coli O157:H7. Phage Φ241 was isolated from an industrial cucumber fermentation at high acidity (pH 3.7) and high salinity (5% NaCl). The phage genome consists of a 157,291 bp circular double-stranded DNA with 203 coding regions and 44.96% GC content.

An enterococcal phage-derived enzyme suppresses graft-versus-host disease

Changes in the gut microbiome have pivotal roles in the pathogenesis of acute graft-versus-host disease (aGVHD) after allogenic haematopoietic cell transplantation (allo-HCT)1-6. However, effective methods for safely resolving gut dysbiosis have not yet been established. An expansion of the pathogen Enterococcus faecalis in the intestine, associated with dysbiosis, has been shown to be a risk factor for aGVHD7-10. Here we analyse the intestinal microbiome of patients with allo-HCT, and find that E. faecalis escapes elimination and proliferates in the intestine by forming biofilms, rather than by acquiring drug-resistance genes. We isolated cytolysin-positive highly pathogenic E. faecalis from faecal samples and identified an anti-E. faecalis enzyme derived from E. faecalis-specific bacteriophages by analysing bacterial whole-genome sequencing data. The antibacterial enzyme had lytic activity against the biofilm of E. faecalis in vitro and in vivo. Furthermore, in aGVHD-induced gnotobiotic mice that were colonized with E. faecalis or with patient faecal samples characterized by the domination of Enterococcus, levels of intestinal cytolysin-positive E. faecalis were decreased and survival was significantly increased in the group that was treated with the E. faecalis-specific enzyme, compared with controls. Thus, administration of a phage-derived antibacterial enzyme that is specific to biofilm-forming pathogenic E. faecalis-which is difficult to eliminate with existing antibiotics-might provide an approach to protect against aGVHD.

Identification of mobile genetic elements with geNomad

Identifying and characterizing mobile genetic elements in sequencing data is essential for understanding their diversity, ecology, biotechnological applications and impact on public health. Here we introduce geNomad, a classification and annotation framework that combines information from gene content and a deep neural network to identify sequences of plasmids and viruses. geNomad uses a dataset of more than 200,000 marker protein profiles to provide functional gene annotation and taxonomic assignment of viral genomes. Using a conditional random field model, geNomad also detects proviruses integrated into host genomes with high precision. In benchmarks, geNomad achieved high classification performance for diverse plasmids and viruses (Matthews correlation coefficient of 77.8% and 95.3%, respectively), substantially outperforming other tools. Leveraging geNomad's speed and scalability, we processed over 2.7 trillion base pairs of sequencing data, leading to the discovery of millions of viruses and plasmids that are available through the IMG/VR and IMG/PR databases. geNomad is available at https://portal.nersc.gov/genomad .

DepoScope: Accurate phage depolymerase annotation and domain delineation using large language models

Bacteriophages (phages) are viruses that infect bacteria. Many of them produce specific enzymes called depolymerases to break down external polysaccharide structures. Accurate annotation and domain identification of these depolymerases are challenging due to their inherent sequence diversity. Hence, we present DepoScope, a machine learning tool that combines a fine-tuned ESM-2 model with a convolutional neural network to identify depolymerase sequences and their enzymatic domains precisely. To accomplish this, we curated a dataset from the INPHARED phage genome database, created a polysaccharide-degrading domain database, and applied sequential filters to construct a high-quality dataset, which is subsequently used to train DepoScope. Our work is the first approach that combines sequence-level predictions with amino-acid-level predictions for accurate depolymerase detection and functional domain identification. In that way, we believe that DepoScope can greatly enhance our understanding of phage-host interactions at the level of depolymerases.

Optimizing bacteriophage treatment of resistant Pseudomonas

Phage therapy is increasing in relevance as an alternative treatment to combat antibiotic resistant bacteria. Phage cocktails are the state-of-the-art method of administering phages in clinical settings, preferred over monophage treatment because of their ability to eliminate multiple bacterial strains and reduce resistance formation. In our study, we compare monophage applications and phage cocktails to our chosen method of phage sequential treatments. To do so, we isolated four novel bacteriophages capable of infecting Pseudomonas alcaligenes T3, a close relative of P. aeruginosa, and characterized them using sequencing and transmission electron microscopy. While investigating monophage treatments, we observed that different phage concentrations had a strong impact on the timing and amount of resistance formation. When using phage cocktails, we observed that P. alcaligenes were capable of forming resistance in the same timespan it took them to become resistant to single phages. We isolated mutants resistant to each single phage as well as mutants exposed to phage cocktails, resulting in bacteria resistant to all four phages at once. Sequencing these mutants showed that different treatments yielded unique single nucleotide polymorphism mutation patterns. In order to combat resistance formation, we added phages one by one in intervals of 24 h, thus managing to delay resistance development and keeping bacterial growth significantly lower compared to phage cocktails.IMPORTANCEWHO declared antimicrobial resistance a top threat to global health; while antibiotics have stood at the forefront in the fight against bacterial infection, the increasing number of multidrug-resistant bacteria highlights a need to branch out in order to address the threat of antimicrobial resistance. Bacteriophages, viruses solely infecting bacteria, could present a solution due to their abundance, versatility, and adaptability. For this study, we isolated new phages infecting a fast-mutating Pseudomonas alcaligenes strain capable of forming resistance within 30 h. By using a sequential treatment approach of adding one phage after another, we were able to curb bacterial growth significantly more compared to state-of-the-art phage cocktails.

Protein Set Transformer: A protein-based genome language model to power high diversity viromics

Exponential increases in microbial and viral genomic data demand transformational advances in scalable, generalizable frameworks for their interpretation. Standard homology-based functional analyses are hindered by the rapid divergence of microbial and especially viral genomes and proteins that significantly decreases the volume of usable data. Here, we present Protein Set Transformer (PST), a protein-based genome language model that models genomes as sets of proteins without considering sparsely available functional labels. Trained on >100k viruses, PST outperformed other homology- and language model-based approaches for relating viral genomes based on shared protein content. Further, PST demonstrated protein structural and functional awareness by clustering capsid-fold-containing proteins with known capsid proteins and uniquely clustering late gene proteins within related viruses. Our data establish PST as a valuable method for diverse viral genomics, ecology, and evolutionary applications. We posit that the PST framework can be a foundation model for microbial genomics when trained on suitable data.

VOGDB-Database of Virus Orthologous Groups

Computational models of homologous protein groups are essential in sequence bioinformatics. Due to the diversity and rapid evolution of viruses, the grouping of protein sequences from virus genomes is particularly challenging. The low sequence similarities of homologous genes in viruses require specific approaches for sequence- and structure-based clustering. Furthermore, the annotation of virus genomes in public databases is not as consistent and up to date as for many cellular genomes. To tackle these problems, we have developed VOGDB, which is a database of virus orthologous groups. VOGDB is a multi-layer database that progressively groups viral genes into groups connected by increasingly remote similarity. The first layer is based on pair-wise sequence similarities, the second layer is based on the sequence profile alignments, and the third layer uses predicted protein structures to find the most remote similarity. VOGDB groups allow for more sensitive homology searches of novel genes and increase the chance of predicting annotations or inferring phylogeny. VOGD B uses all virus genomes from RefSeq and partially reannotates them. VOGDB is updated with every RefSeq release. The unique feature of VOGDB is the inclusion of both prokaryotic and eukaryotic viruses in the same clustering process, which makes it possible to explore old evolutionary relationships of the two groups. VOGDB is freely available at vogdb.org under the CC BY 4.0 license.

Distinct prophage gene profiles of Staphylococcus aureus strains from atopic dermatitis patients and healthy individuals

Staphylococcus aureus strains exhibit varying associations with atopic dermatitis (AD), but the genetic determinants underpinning the pathogenicity are yet to be fully characterized. To reveal the genetic differences between S. aureus strains from AD patients and healthy individuals (HE), we developed and employed a random forest classifier to identify potential marker genes responsible for their phenotypic variations. The classifier was able to effectively distinguish strains from AD and HE. We also uncovered strong links between certain marker genes and phage functionalities, with phage holin emerging as the most pivotal differentiating factor. Further examination of S. aureus gene content highlighted the genetic diversity and functional implications of prophages in driving differentiation between strains from AD and HE. The HE group exhibited greater gene content diversity, largely influenced by their prophages. While strains from both AD and HE universally housed prophages, those in the HE group were distinctively higher at the strain level. Moreover, although prophages in the HE group exhibited variously higher enrichment of differential functions, the AD group displayed a notable enrichment of virulence factors within their prophages, underscoring the important contribution of prophages to the pathogenesis of AD-associated strains. Overall, prophages significantly shape the genetic and functional profiles of S. aureus strains, shedding light on their pathogenic potential and elucidating the mechanisms behind the phenotypic variations in AD and HE environments.

IMPORTANCE

Through a nuanced exploration of Staphylococcus aureus strains obtained from atopic dermatitis (AD) patients and healthy controls (HE), our research unveils pivotal genetic determinants influencing their pathogenic associations. Utilizing a random forest classifier, we illuminate distinct marker genes, with phage holin emerging as a critical differential factor, revealing the profound impact of prophages on genetic and pathogenic profiles. HE strains exhibited a diverse gene content, notably shaped by unique, heightened prophages. Conversely, AD strains emphasized a pronounced enrichment of virulence factors within prophages, signifying their key role in AD pathogenesis. This work crucially highlights prophages as central architects of the genetic and functional attributes of S. aureus strains, providing vital insights into pathogenic mechanisms and phenotypic variations, thereby paving the way for targeted AD therapeutic approaches and management strategies by demystifying specific genetic and pathogenic mechanisms.

Combinations of Bacteriophage Are Efficacious against Multidrug-Resistant Pseudomonas aeruginosa and Enhance Sensitivity to Carbapenem Antibiotics

The Gram-negative ESKAPE bacterium Pseudomonas aeruginosa has become a pathogen of serious concern due its extensive multi-drug resistance (MDR) profile, widespread incidences of hospital-acquired infections throughout the United States, and high occurrence in wound infections suffered by warfighters serving abroad. Bacteriophage (phage) therapy has received renewed attention as an alternative therapeutic option against recalcitrant bacterial infections, both as multi-phage cocktails and in combination with antibiotics as synergistic pairings. Environmental screening and phage enrichment has yielded three lytic viruses capable of infecting the MDR P. aeruginosa strain PAO1. Co-administration of each phage with the carbapenem antibiotics ertapenem, imipenem, and meropenem generated enhanced overall killing of bacteria beyond either phage or drug treatments alone. A combination cocktail of all three phages was completely inhibitory to growth, even without antibiotics. The same 3× phage cocktail also disrupted PAO1 biofilms, reducing biomass by over 75% compared to untreated biofilms. Further, the phage cocktail demonstrated broad efficacy as well, capable of infecting 33 out of 100 diverse clinical isolate strains of P. aeruginosa. Together, these results indicate a promising approach for designing layered medical countermeasures to potentiate antibiotic activity and possibly overcome resistance against recalcitrant, MDR bacteria such as P. aeruginosa. Combination therapy, either by synergistic phage-antibiotic pairings, or by phage cocktails, presents a means of controlling mutations that can allow for bacteria to gain a competitive edge.

Complete genome sequence of Enterococcus faecalis phage EF_RCK

A lytic bacteriophage EF_RCK infecting Enterococcus faecalis was isolated from a water sample collected in a raw cockle storage container at Taman Ria market, Sungai Petani, Malaysia. The phage has a 57,848-bp double-stranded DNA genome harboring 107 protein-encoding genes and shares 90.9% nucleotide similarity with Enterococcus phage EFKL (Saphexavirus genus).

Profiling the interplay and coevolution of Microcystis aeruginosa and cyanosiphophage Mic1

The cyanosiphophage Mic1 specifically infects the bloom-forming Microcystis aeruginosa FACHB 1339 from Lake Chaohu, China. Previous genomic analysis showed that its 92,627 bp double-stranded DNA genome consists of 98 putative open reading frames, 63% of which are of unknown function. Here, we investigated the transcriptome dynamics of Mic1 and its host using RNA sequencing. In the early, middle, and late phases of the 10 h lytic cycle, the Mic1 genes are sequentially expressed and could be further temporally grouped into two distinct clusters in each phase. Notably, six early genes, including gp49 that encodes a TnpB-like transposase, immediately reach the highest transcriptional level in half an hour, representing a pioneer cluster that rapidly regulates and redirects host metabolism toward the phage. An in-depth analysis of the host transcriptomic profile in response to Mic1 infection revealed significant upregulation of a polyketide synthase pathway and a type III-B CRISPR system, accompanied by moderate downregulation of the photosynthesis and key metabolism pathways. The constant increase of phage transcripts and relatively low replacement rate over the host transcripts indicated that Mic1 utilizes a unique strategy to gradually take over a small portion of host metabolism pathways after infection. In addition, genomic analysis of a less-infective Mic1 and a Mic1-resistant host strain further confirmed their dynamic interplay and coevolution via the frequent horizontal gene transfer. These findings provide insights into the mutual benefit and symbiosis of the highly polymorphic cyanobacteria M. aeruginosa and cyanophages.

IMPORTANCE

The highly polymorphic Microcystis aeruginosa is one of the predominant bloom-forming cyanobacteria in eutrophic freshwater bodies and is infected by diverse and abundant cyanophages. The presence of a large number of defense systems in M. aeruginosa genome suggests a dynamic interplay and coevolution with the cyanophages. In this study, we investigated the temporal gene expression pattern of Mic1 after infection and the corresponding transcriptional responses of its host. Moreover, the identification of a less-infective Mic1 and a Mic1-resistant host strain provided the evolved genes in the phage-host coevolution during the multiple-generation cultivation in the laboratory. Our findings enrich the knowledge on the interplay and coevolution of M. aeruginosa and its cyanophages and lay the foundation for the future application of cyanophage as a potential eco-friendly and bio-safe agent in controlling the succession of harmful cyanobacterial blooms.

Seed coating with phages for sustainable plant biocontrol of plant pathogens and influence of the seed coat mucilage

Pathogens resistant to classical control strategies pose a significant threat to crop yield, with seeds being a major transmission route. Bacteriophages, viruses targeting bacteria, offer an environmentally sustainable biocontrol solution. In this study, we isolated and characterized two novel phages, Athelas and Alfirin, which infect Pseudomonas syringae and Agrobacterium fabrum, respectively, and included the recently published Pfeifenkraut phage infecting Xanthomonas translucens. Using a simple immersion method, phages coated onto seeds successfully lysed bacteria post air-drying. The seed coat mucilage (SCM), a polysaccharide-polymer matrix exuded by seeds, plays a critical role in phage binding. Seeds with removed mucilage formed five to 10 times less lysis zones compared to those with mucilage. The podovirus Athelas showed the highest mucilage dependency. Phages from the Autographiviridae family also depended on mucilage for seed adhesion. Comparative analysis of Arabidopsis SCM mutants suggested the diffusible cellulose as a key component for phage binding. Long-term activity tests demonstrated high phage stability on seed surfaces and significantly increasing seedling survival rates in the presence of pathogens. Using non-virulent host strains enhanced phage presence on seeds but also has potential limitations. These findings highlight phage-based interventions as promising, sustainable strategies for combating pathogen resistance and improving crop yield.

VIBES: a workflow for annotating and visualizing viral sequences integrated into bacterial genomes

Bacteriophages are viruses that infect bacteria. Many bacteriophages integrate their genomes into the bacterial chromosome and become prophages. Prophages may substantially burden or benefit host bacteria fitness, acting in some cases as parasites and in others as mutualists. Some prophages have been demonstrated to increase host virulence. The increasing ease of bacterial genome sequencing provides an opportunity to deeply explore prophage prevalence and insertion sites. Here we present VIBES (Viral Integrations in Bacterial genomES), a workflow intended to automate prophage annotation in complete bacterial genome sequences. VIBES provides additional context to prophage annotations by annotating bacterial genes and viral proteins in user-provided bacterial and viral genomes. The VIBES pipeline is implemented as a Nextflow-driven workflow, providing a simple, unified interface for execution on local, cluster and cloud computing environments. For each step of the pipeline, a container including all necessary software dependencies is provided. VIBES produces results in simple tab-separated format and generates intuitive and interactive visualizations for data exploration. Despite VIBES's primary emphasis on prophage annotation, its generic alignment-based design allows it to be deployed as a general-purpose sequence similarity search manager. We demonstrate the utility of the VIBES prophage annotation workflow by searching for 178 Pf phage genomes across 1072 Pseudomonas spp. genomes.

Isolation and characterization of a phage collection against Pseudomonas putida

The environmental bacterium, Pseudomonas putida, possesses a broad spectrum of metabolic pathways. This makes it highly promising for use in biotechnological production as a cell factory, as well as in bioremediation strategies to degrade various aromatic pollutants. For P. putida to flourish in its environment, it must withstand the continuous threats posed by bacteriophages. Interestingly, until now, only a handful of phages have been isolated for the commonly used laboratory strain, P. putida KT2440, and no phage defence mechanisms have been characterized. In this study, we present a new Collection of Environmental P. putida Phages from Estonia, or CEPEST. This collection comprises 67 double-stranded DNA phages, which belong to 22 phage species and 9 phage genera. Our findings reveal that most phages in the CEPEST collection are more infectious at lower temperatures, have a narrow host range, and require an intact lipopolysaccharide for P. putida infection. Furthermore, we show that cryptic prophages present in the P. putida chromosome provide strong protection against the infection of many phages. However, the chromosomal toxin-antitoxin systems do not play a role in the phage defence of P. putida. This research provides valuable insights into the interactions between P. putida and bacteriophages, which could have significant implications for biotechnological and environmental applications.

Tradeoffs Between Evolved Phage Resistance and Antibiotic Susceptibility in a Highly Drug-Resistant Cystic Fibrosis-Derived Pseudomonas aeruginosa Strain

Background

Multi-drug resistant pathogens pose significant challenges towards the effective resolution of bacterial infections. A promising alternative strategy is phage therapy in which limited applications has afforded lifesaving resolution from drug resistant pathogens. However, adoption of this strategy is hampered by narrow bacteriophage host ranges, and as with antibiotics, bacteria can acquire resistance to phage.

Methods

To address these issues, we isolated 25 broad-host range phages against multiple cystic fibrosis (CF)-derived P. aeruginosa clinical strains thus promoting their application against conspecific pathogens. To investigate evolved resistance to phage in relation to antibiotic resistance, one CF-derived P. aeruginosa strain was exposed to a lytic phage over a short time scale.

Results

Trade-offs were observed in which evolved phage resistant P. aeruginosa strains showed decreased resistance to antibiotics. These traits that likely reflect single nucleotide polymorphisms.

Conclusion

Results suggest phage and antibiotics may be a combined approach to treat bacterial infections.

Genomic and taxonomic evaluation of 38 Treponema prophage sequences

BACKGROUND

Despite Spirochetales being a ubiquitous and medically important order of bacteria infecting both humans and animals, there is extremely limited information regarding their bacteriophages. Of the genus Treponema, there is just a single reported characterised prophage.

RESULTS

We applied a bioinformatic approach on 24 previously published Treponema genomes to identify and characterise putative treponemal prophages. Thirteen of the genomes did not contain any detectable prophage regions. The remaining eleven contained 38 prophage sequences, with between one and eight putative prophages in each bacterial genome. The prophage regions ranged from 12.4 to 75.1 kb, with between 27 and 171 protein coding sequences. Phylogenetic analysis revealed that 24 of the prophages formed three distinct sequence clusters, identifying putative myoviral and siphoviral morphology. ViPTree analysis demonstrated that the identified sequences were novel when compared to known double stranded DNA bacteriophage genomes.

CONCLUSIONS

In this study, we have started to address the knowledge gap on treponeme bacteriophages by characterising 38 prophage sequences in 24 treponeme genomes. Using bioinformatic approaches, we have been able to identify and compare the prophage-like elements with respect to other bacteriophages, their gene content, and their potential to be a functional and inducible bacteriophage, which in turn can help focus our attention on specific prophages to investigate further.

Exploring the roles of ribosomal peptides in prokaryote-phage interactions through deep learning-enabled metagenome mining

BACKGROUND

Microbial secondary metabolites play a crucial role in the intricate interactions within the natural environment. Among these metabolites, ribosomally synthesized and post-translationally modified peptides (RiPPs) are becoming a promising source of therapeutic agents due to their structural diversity and functional versatility. However, their biosynthetic capacity and ecological functions remain largely underexplored.

RESULTS

Here, we aim to explore the biosynthetic profile of RiPPs and their potential roles in the interactions between microbes and viruses in the ocean, which encompasses a vast diversity of unique biomes that are rich in interactions and remains chemically underexplored. We first developed TrRiPP to identify RiPPs from ocean metagenomes, a deep learning method that detects RiPP precursors in a hallmark gene-independent manner to overcome the limitations of classic methods in processing highly fragmented metagenomic data. Applying this method to metagenomes from the global ocean microbiome, we uncover a diverse array of previously uncharacterized putative RiPP families with great novelty and diversity. Through correlation analysis based on metatranscriptomic data, we observed a high prevalence of antiphage defense-related and phage-related protein families that were co-expressed with RiPP families. Based on this putative association between RiPPs and phage infection, we constructed an Ocean Virus Database (OVD) and established a RiPP-involving host-phage interaction network through host prediction and co-expression analysis, revealing complex connectivities linking RiPP-encoding prokaryotes, RiPP families, viral protein families, and phages. These findings highlight the potential of RiPP families involved in prokaryote-phage interactions and coevolution, providing insights into their ecological functions in the ocean microbiome.

CONCLUSIONS

This study provides a systematic investigation of the biosynthetic potential of RiPPs from the ocean microbiome at a global scale, shedding light on the essential insights into the ecological functions of RiPPs in prokaryote-phage interactions through the integration of deep learning approaches, metatranscriptomic data, and host-phage connectivity. This study serves as a valuable example of exploring the ecological functions of bacterial secondary metabolites, particularly their associations with unexplored microbial interactions. Video Abstract.

Utilizing profile hidden Markov model databases for discovering viruses from metagenomic data: a comprehensive review

Profile hidden Markov models (pHMMs) are able to achieve high sensitivity in remote homology search, making them popular choices for detecting novel or highly diverged viruses in metagenomic data. However, many existing pHMM databases have different design focuses, making it difficult for users to decide the proper one to use. In this review, we provide a thorough evaluation and comparison for multiple commonly used profile HMM databases for viral sequence discovery in metagenomic data. We characterized the databases by comparing their sizes, their taxonomic coverage, and the properties of their models using quantitative metrics. Subsequently, we assessed their performance in virus identification across multiple application scenarios, utilizing both simulated and real metagenomic data. We aim to offer researchers a thorough and critical assessment of the strengths and limitations of different databases. Furthermore, based on the experimental results obtained from the simulated and real metagenomic data, we provided practical suggestions for users to optimize their use of pHMM databases, thus enhancing the quality and reliability of their findings in the field of viral metagenomics.

An autonomous plasmid as an inovirus phage satellite

Bacterial viruses (phages) are potent agents of lateral gene transfer and thus are important drivers of evolution. A group of mobile genetic elements, referred to as phage satellites, exploits phages to disseminate their own genetic material. Here, we isolated a novel member of the family Inoviridae, Shewanella phage Dolos, along with an autonomously replicating plasmid, pDolos. Dolos causes a chronic infection in its host Shewanella oneidensis by phage production with only minor effects on the host cell proliferation. When present, plasmid pDolos hijacks Dolos functions to be predominantly packaged into phage virions and released into the environment and, thus, acts as a phage satellite. pDolos can disseminate further genetic material encoding, e.g., resistances or fluorophores to host cells sensitive to Dolos infection. Given the rather simple requirements of a plasmid for takeover of an inovirus and the wide distribution of phages of this group, we speculate that similar phage-satellite systems are common among bacteria.IMPORTANCEPhage satellites are mobile genetic elements, which hijack phages to be transferred to other host cells. The vast majority of these phage satellites integrate within the host's chromosome, and they all carry remaining phage genes. Here, we identified a novel phage satellite, pDolos, which uses an inovirus for dissemination. pDolos (i) remains as an autonomously replicating plasmid within its host, (ii) does not carry recognizable phage genes, and (iii) is smaller than any other phage satellites identified so far. Thus, pDolos is the first member of a new class of phage satellites, which resemble natural versions of phagemids.

Time-series metaproteogenomics of a high-CO2 aquifer reveals active viruses with fluctuating abundances and broad host ranges

Ecosystems subject to mantle degassing are of particular interest for understanding global biogeochemistry, as their microbiomes are shaped by prolonged exposure to high CO2 and have recently been suggested to be highly active. While the genetic diversity of bacteria and archaea in these deep biosphere systems have been studied extensively, little is known about how viruses impact these microbial communities. Here, we show that the viral community in a high-CO2 cold-water geyser (Wallender Born, Germany) undergoes substantial fluctuations over a period of 12 days, although the corresponding prokaryotic community remains stable, indicating a newly observed "infect to keep in check" strategy that maintains prokaryotic community structure. We characterized the viral community using metagenomics and metaproteomics, revealing 8 654 viral operational taxonomic units (vOTUs). CRISPR spacer-to-protospacer matching linked 278 vOTUs to 32 hosts, with many vOTUs sharing hosts from different families. High levels of viral structural proteins present in the metaproteome (several structurally annotated based on AlphaFold models) indicate active virion production at the time of sampling. Viral genomes expressed many proteins involved in DNA metabolism and manipulation, and encoded for auxiliary metabolic genes, which likely bolster phosphate and sulfur metabolism of their hosts. The active viral community encodes genes to facilitate acquisition and transformation of host nutrients, and appears to consist of many nutrient-demanding members, based on abundant virion proteins. These findings indicate viruses are inextricably linked to the biogeochemical cycling in this high-CO2 environment and substantially contribute to prokaryotic community stability in the deep biosphere hotspots.

Frequent nonhomologous replacement of replicative helicase loaders by viruses in Vibrionaceae

Several microbial genomes lack textbook-defined essential genes. If an essential gene is absent from a genome, then an evolutionarily independent gene of unknown function complements its function. Here, we identified frequent nonhomologous replacement of an essential component of DNA replication initiation, a replicative helicase loader gene, in Vibrionaceae. Our analysis of Vibrionaceae genomes revealed two genes with unknown function, named vdhL1 and vdhL2, that were substantially enriched in genomes without the known helicase-loader genes. These genes showed no sequence similarities to genes with known function but encoded proteins structurally similar with a viral helicase loader. Analyses of genomic syntenies and coevolution with helicase genes suggested that vdhL1/2 encodes a helicase loader. The in vitro assay showed that Vibrio harveyi VdhL1 and Vibrio ezurae VdhL2 promote the helicase activity of DnaB. Furthermore, molecular phylogenetics suggested that vdhL1/2 were derived from phages and replaced an intrinsic helicase loader gene of Vibrionaceae over 20 times. This high replacement frequency implies the host's advantage in acquiring a viral helicase loader gene.

The virtue of training: extending phage host spectra against vancomycin-resistant Enterococcus faecium strains using the Appelmans method

Phage therapy has (re)emerged as a serious possibility for combating multidrug-resistant bacterial infections, including those caused by vancomycin-resistant Enterococcus faecium strains. These opportunistic pathogens belong to a specific clonal complex 17, against which relatively few phages have been screened. We isolated a collection of 21 virulent phages growing on these vancomycin-resistant isolates. Each of these phages harbored a typical narrow plaquing host range, lysing at most 5 strains and covering together 10 strains of our panel of 14 clinical isolates. To enlarge the host spectrum of our phages, the Appelmans protocol was used. We mixed four out of our most complementary phages in a cocktail that we iteratively grew on eight naive strains from our panel, of which six were initially refractory to at least three of the combined phages. Fifteen successive passages permitted to significantly improve the lytic activity of the cocktail, from which phages with extended host ranges within the E. faecium species could be isolated. A single evolved phage able to kill up to 10 of the 14 initial E. faecium strains was obtained, and it barely infected nearby species. All evolved phages had acquired point mutations or a recombination event in the tail fiber genetic region, suggesting these genes might have driven phage evolution by contributing to their extended host spectra.

Decoding huge phage diversity: a taxonomic classification of Lak megaphages

High-throughput sequencing for uncultivated viruses has accelerated the understanding of global viral diversity and uncovered viral genomes substantially larger than any that have so far been cultured. Notably, the Lak phages are an enigmatic group of viruses that present some of the largest known phage genomes identified in human and animal microbiomes, and are dissimilar to any cultivated viruses. Despite the wealth of viral diversity that exists within sequencing datasets, uncultivated viruses have rarely been used for taxonomic classification. We investigated the evolutionary relationships of 23 Lak phages and propose a taxonomy for their classification. Predicted protein analysis revealed the Lak phages formed a deeply branching monophyletic clade within the class Caudoviricetes which contained no other phage genomes. One of the interesting features of this clade is that all current members are characterised by an alternative genetic code. We propose the Lak phages belong to a new order, the 'Grandevirales'. Protein and nucleotide-based analyses support the creation of two families, three sub-families, and four genera within the order 'Grandevirales'. We anticipate that the proposed taxonomy of Lak megaphages will simplify the future classification of related viral genomes as they are uncovered. Continued efforts to classify divergent viruses are crucial to aid common analyses of viral genomes and metagenomes.

Bacteriophages to control Vibrio alginolyticus in live feeds prior to their administration in larviculture

AIMS

This study aimed to evaluate the efficiency of two phages [VB_VaC_TDDLMA (phage TDD) and VB_VaC_SRILMA (phage SRI)] alone and in a cocktail to control Vibrio alginolyticus in brine shrimp before their administration in larviculture.

METHODS AND RESULTS

Phages were isolated from seawater samples and characterized by host spectrum, growth parameters, adsorption rate, genomic analysis, and inactivation efficiency. Both phages belong to the Caudoviricetes class and lack known virulence or antibiotic-resistance genes. They exhibit specificity, infecting only their host, V. alginolyticus CECT 521. Preliminary experiments in a culture medium showed that phage TDD (reduction of 5.8 log CFU ml-1 after 10 h) outperformed phage SRI (reduction of 4.6 log CFU ml-1 after 6 h) and the cocktail TDD/SRI (reduction of 5.2 log CFU ml-1 after 8 h). In artificial marine water experiments with Artemia franciscana, both single phage suspensions and the phage cocktail, effectively inactivated V. alginolyticus in culture water (reduction of 4.3, 2.1, and 1.9 log CFU ml-1 for phages TDD, SRI, and the phage cocktail, respectively, after 12 h) and in A. franciscana (reduction of 51.6%, 87.3%, and 85.3% for phages TDD, SRI, and the phage cocktail, respectively, after 24 h). The two phages and the phage cocktail did not affect A. franciscana natural microbiota or other Vibrio species in the brine shrimp.

CONCLUSIONS

The results suggest that phages can safely and effectively control V. alginolyticus in A. franciscana prior to its administration in larviculture.

Systematic genome-wide discovery of host factors governing bacteriophage infectivity

Bacterial host factors regulate the infection cycle of bacteriophages. Except for some well-studied host factors (e.g., receptors or restriction-modification systems), the contribution of the rest of the host genome on phage infection remains poorly understood. We developed PHAGEPACK, a pooled assay that systematically and comprehensively measures each host-gene impact on phage fitness. PHAGEPACK combines CRISPR interference with phage packaging to link host perturbation to phage fitness during active infection. Using PHAGEPACK, we constructed a genome-wide map of genes impacting T7 phage fitness in permissive E. coli, revealing pathways previously unknown to affect phage packaging. When applied to the non-permissive E. coli O121, PHAGEPACK identified pathways leading to host resistance; their removal increased phage susceptibility up to a billion-fold. Bioinformatic analysis indicates phage genomes carry homologs or truncations of key host factors, potentially for fitness advantage. In summary, PHAGEPACK offers valuable insights into phage-host interactions, phage evolution, and bacterial resistance.

Benchmarking bioinformatic virus identification tools using real-world metagenomic data across biomes

BACKGROUND

As most viruses remain uncultivated, metagenomics is currently the main method for virus discovery. Detecting viruses in metagenomic data is not trivial. In the past few years, many bioinformatic virus identification tools have been developed for this task, making it challenging to choose the right tools, parameters, and cutoffs. As all these tools measure different biological signals, and use different algorithms and training and reference databases, it is imperative to conduct an independent benchmarking to give users objective guidance.

RESULTS

We compare the performance of nine state-of-the-art virus identification tools in thirteen modes on eight paired viral and microbial datasets from three distinct biomes, including a new complex dataset from Antarctic coastal waters. The tools have highly variable true positive rates (0-97%) and false positive rates (0-30%). PPR-Meta best distinguishes viral from microbial contigs, followed by DeepVirFinder, VirSorter2, and VIBRANT. Different tools identify different subsets of the benchmarking data and all tools, except for Sourmash, find unique viral contigs. Performance of tools improved with adjusted parameter cutoffs, indicating that adjustment of parameter cutoffs before usage should be considered.

CONCLUSIONS

Together, our independent benchmarking facilitates selecting choices of bioinformatic virus identification tools and gives suggestions for parameter adjustments to viromics researchers.

Diverse and abundant phages exploit conjugative plasmids

Phages exert profound evolutionary pressure on bacteria by interacting with receptors on the cell surface to initiate infection. While the majority of phages use chromosomally encoded cell surface structures as receptors, plasmid-dependent phages exploit plasmid-encoded conjugation proteins, making their host range dependent on horizontal transfer of the plasmid. Despite their unique biology and biotechnological significance, only a small number of plasmid-dependent phages have been characterized. Here we systematically search for new plasmid-dependent phages targeting IncP and IncF plasmids using a targeted discovery platform, and find that they are common and abundant in wastewater, and largely unexplored in terms of their genetic diversity. Plasmid-dependent phages are enriched in non-canonical types of phages, and all but one of the 65 phages we isolated were non-tailed, and members of the lipid-containing tectiviruses, ssDNA filamentous phages or ssRNA phages. We show that plasmid-dependent tectiviruses exhibit profound differences in their host range which is associated with variation in the phage holin protein. Despite their relatively high abundance in wastewater, plasmid-dependent tectiviruses are missed by metaviromic analyses, underscoring the continued importance of culture-based phage discovery. Finally, we identify a tailed phage dependent on the IncF plasmid, and find related structural genes in phages that use the orthogonal type 4 pilus as a receptor, highlighting the evolutionarily promiscuous use of these distinct contractile structures by multiple groups of phages. Taken together, these results indicate plasmid-dependent phages play an under-appreciated evolutionary role in constraining horizontal gene transfer via conjugative plasmids.