PHROG: families of prokaryotic virus proteins clustered using remote homology

Isolation and genomic characterization of Staphylococcus aureus bacteriophages from Chennai, India

We isolated and characterized two lytic bacteriophages against Staphylococcus aureus named TANUVAS_MVC-VPHSA1 and TANUVAS_MVC-VPHSA2, with the aim of investigating their genomic and structural features. The bacteriophages belong to the Caudoviricetes, and their genomes have sizes of 50,505 and 50,516 base pairs with a GC content of 41.4%.

MetaCerberus: distributed highly parallelized HMM-based processing for robust functional annotation across the tree of life

MOTIVATION

MetaCerberus is a massively parallel, fast, low memory, scalable annotation tool for inference gene function across genomes to metacommunities. MetaCerberus provides an elusive HMM/HMMER-based tool at a rapid scale with low memory. It offers scalable gene elucidation to major public databases, including KEGG (KO), COGs, CAZy, FOAM, and specific databases for viruses, including VOGs and PHROGs, from single genomes to metacommunities.

RESULTS

MetaCerberus is 1.3× as fast on a single node than eggNOG-mapper v2 on 5× less memory using an exclusively HMM/HMMER mode. In a direct comparison, MetaCerberus provides better annotation of viruses, phages, and archaeal viruses than DRAM, Prokka, or InterProScan. MetaCerberus annotates more KOs across domains when compared to DRAM, with a 186× smaller database, and with 63× less memory. MetaCerberus is fully integrated for automatic analysis of statistics and pathways using differential statistic tools (i.e. DESeq2 and edgeR), pathway enrichment (GAGE R), and pathview R. MetaCerberus provides a novel tool for unlocking the biosphere across the tree of life at scale.

AVAILABILITY AND IMPLEMENTATION

MetaCerberus is written in Python and distributed under a BSD-3 license. The source code of MetaCerberus is freely available at https://github.com/raw-lab/metacerberus compatible with Python 3 and works on both Mac OS X and Linux. MetaCerberus can also be easily installed using bioconda: mamba create -n metacerberus -c bioconda -c conda-forge metacerberus.

Phage-plasmids promote recombination and emergence of phages and plasmids

Phages and plasmids are regarded as distinct types of mobile genetic elements that drive bacterial evolution by horizontal gene transfer. However, the distinction between both types is blurred by the existence of elements known as prophage-plasmids or phage-plasmids, which transfer horizontally between cells as viruses and vertically within cellular lineages as plasmids. Here, we study gene flow between the three types of elements. We show that the gene repertoire of phage-plasmids overlaps with those of phages and plasmids. By tracking recent recombination events, we find that phage-plasmids exchange genes more frequently with plasmids than with phages, and that direct gene exchange between plasmids and phages is less frequent in comparison. The results suggest that phage-plasmids can mediate gene flow between plasmids and phages, including exchange of mobile element core functions, defense systems, and antibiotic resistance. Moreover, a combination of gene transfer and gene inactivation may result in the conversion of elements. For example, gene loss turns P1-like phage-plasmids into integrative prophages or into plasmids (that are no longer phages). Remarkably, some of the latter have acquired conjugation-related functions to became mobilisable by conjugation. Thus, our work indicates that phage-plasmids can play a key role in the transfer of genes across mobile elements within their hosts, and can act as intermediates in the conversion of one type of element into another.

Potential of an Isolated Bacteriophage to Inactivate Klebsiella pneumoniae: Preliminary Studies to Control Urinary Tract Infections

Urinary tract infections (UTIs) caused by resistant Klebsiella pneumoniae can lead to severe clinical complications and even death. An alternative treatment option for infected patients is using bacteriophages. In the present study, we isolated phage VB_KPM_KP1LMA (KP1LMA) from sewage water using a K. pneumoniae strain as a host. Whole-genome analysis indicated that the genome was a double-stranded linear 176,096-bp long DNA molecule with 41.8% GC content and did not contain virulence or antibiotic resistance genes. The inactivation potential of phage KP1LMA was assessed in broth at an MOI of 1 and 10, and a maximum inactivation of 4.9 and 5.4 log CFU/mL, respectively, was observed after 9 h. The efficacy at an MOI of 10 was also assessed in urine to evaluate the phage's performance in an acidic environment. A maximum inactivation of 3.8 log CFU/mL was observed after 9 h. The results suggest that phage KP1LMA could potentially control a UTI caused by this strain of K. pneumoniae, indicating that the same procedure can be used to control UTIs caused by other strains if new specific phages are isolated. Although phage KP1LMA has a narrow host range, in the future, efforts can be made to expand its spectrum of activity and also to combine this phage with others, potentially enabling its use against other K. pneumoniae strains involved in UTIs.

Dynamic nature of viral and bacterial communities in human faeces

Bacteriophages are a major component of the gut microbiome and are believed to play a role in establishment and stabilization of microbial communities by influencing taxonomic and functional diversity. We show that the activity of lytic and temperate phages can also significantly affect bacterial community structure in a model of extended colonic retention. Intact fresh human feces were incubated anaerobically at 37°C without homogenization and subjected to metagenomic sequencing. We observed subject-specific blooms and collapses of selected bacteriophage and bacterial populations within some individuals. Most notable were striking collapses of Prevotella populations accompanied by increases in specific bacteriophages. In a number of cases, we even observed a shift from one bacterial "enterotype" to another within 48 h. These results confirm that intact feces represents a highly dynamic ecological system and suggests that colonic retention time could have a profound effect on microbiome composition, including a significant impact by bacteriophages.

Prophage enhances the ability of deep-sea bacterium Shewanella psychrophila WP2 to utilize D-amino acid

Prophages are prevalent in the marine bacterial genomes and reshape the physiology and metabolism of their hosts. However, whether and how prophages influence the microbial degradation of D-amino acids (D-AAs), which is one of the widely distributed recalcitrant dissolved organic matters (RDOMs) in the ocean, remain to be explored. In this study, we addressed this issue in a representative marine bacterium, Shewanella psychrophila WP2 (WP2), and its integrated prophage SP1. Notably, compared to the WP2 wild-type strain, the SP1 deletion mutant of WP2 (WP2ΔSP1) exhibited a significantly lower D-glutamate (D-Glu) consumption rate and longer lag phase when D-Glu was used as the sole nitrogen source. The subsequent transcriptome analysis identified 1,523 differentially expressed genes involved in diverse cellular processes, especially that multiple genes related to inorganic nitrogen metabolism were highly upregulated. In addition, the dynamic profiles of ammonium, nitrate, and nitrite were distinct between the culture media of WP2 and WP2ΔSP1. Finally, we provide evidence that SP1 conferred a competitive advantage to WP2 when D-Glu was used as the sole nitrogen source and SP1-like phages may be widely distributed in the global ocean. Taken together, these findings offer novel insight into the influences of prophages on host metabolism and RDOM cycling in marine environments.IMPORTANCEThis work represents the first exploration of the impact of prophages on the D-amino acid (D-AA) metabolism of deep-sea bacteria. By using S. psychrophila WP2 and its integrated prophage SP1 as a representative system, we found that SP1 can significantly increase the catabolism rate of WP2 to D-glutamate and produce higher concentrations of ammonium, resulting in faster growth and competitive advantages. Our findings not only deepen our understanding of the interaction between deep-sea prophages and hosts but also provide new insights into the ecological role of prophages in refractory dissolved organic matter and the nitrogen cycle in deep oceans.

Large language models improve annotation of prokaryotic viral proteins

Viral genomes are poorly annotated in metagenomic samples, representing an obstacle to understanding viral diversity and function. Current annotation approaches rely on alignment-based sequence homology methods, which are limited by the paucity of characterized viral proteins and divergence among viral sequences. Here we show that protein language models can capture prokaryotic viral protein function, enabling new portions of viral sequence space to be assigned biologically meaningful labels. When applied to global ocean virome data, our classifier expanded the annotated fraction of viral protein families by 29%. Among previously unannotated sequences, we highlight the identification of an integrase defining a mobile element in marine picocyanobacteria and a capsid protein that anchors globally widespread viral elements. Furthermore, improved high-level functional annotation provides a means to characterize similarities in genomic organization among diverse viral sequences. Protein language models thus enhance remote homology detection of viral proteins, serving as a useful complement to existing approaches.

The gut virome is associated with stress-induced changes in behaviour and immune responses in mice

The microbiota-gut-brain axis has been shown to play an important role in the stress response, but previous work has focused primarily on the role of the bacteriome. The gut virome constitutes a major portion of the microbiome, with bacteriophages having the potential to remodel bacteriome structure and activity. Here we use a mouse model of chronic social stress, and employ 16S rRNA and whole metagenomic sequencing on faecal pellets to determine how the virome is modulated by and contributes to the effects of stress. We found that chronic stress led to behavioural, immune and bacteriome alterations in mice that were associated with changes in the bacteriophage class Caudoviricetes and unassigned viral taxa. To determine whether these changes were causally related to stress-associated behavioural or physiological outcomes, we conducted a faecal virome transplant from mice before stress and autochthonously transferred it to mice undergoing chronic social stress. The transfer of the faecal virome protected against stress-associated behaviour sequelae and restored stress-induced changes in select circulating immune cell populations, cytokine release, bacteriome alterations and gene expression in the amygdala. These data provide evidence that the virome plays a role in the modulation of the microbiota-gut-brain axis during stress, indicating that these viral populations should be considered when designing future microbiome-directed therapies.

Ultraconserved bacteriophage genome sequence identified in 1300-year-old human palaeofaeces

Bacteriophages are widely recognised as rapidly evolving biological entities. However, knowledge about ancient bacteriophages is limited. Here, we analyse DNA sequence datasets previously generated from ancient palaeofaeces and human gut-content samples, and identify an ancient phage genome nearly identical to present-day Mushuvirus mushu, a virus that infects gut commensal bacteria. The DNA damage patterns of the genome are consistent with its ancient origin and, despite 1300 years of evolution, the ancient Mushuvirus genome shares 97.7% nucleotide identity with its modern counterpart, indicating a long-term relationship between the prophage and its host. In addition, we reconstruct and authenticate 297 other phage genomes from the last 5300 years, including those belonging to unknown families. Our findings demonstrate the feasibility of reconstructing ancient phage genome sequences, thus expanding the known virosphere and offering insights into phage-bacteria interactions spanning several millennia.

Molecular Characterization and Genome Mechanical Features of Two Newly Isolated Polyvalent Bacteriophages Infecting Pseudomonas syringae pv. garcae

Coffee plants have been targeted by a devastating bacterial disease, a condition known as bacterial blight, caused by the phytopathogen Pseudomonas syringae pv. garcae (Psg). Conventional treatments of coffee plantations affected by the disease involve frequent spraying with copper- and kasugamycin-derived compounds, but they are both highly toxic to the environment and stimulate the appearance of bacterial resistance. Herein, we report the molecular characterization and mechanical features of the genome of two newly isolated (putative polyvalent) lytic phages for Psg. The isolated phages belong to class Caudoviricetes and present a myovirus-like morphotype belonging to the genuses Tequatrovirus (PsgM02F) and Phapecoctavirus (PsgM04F) of the subfamilies Straboviridae (PsgM02F) and Stephanstirmvirinae (PsgM04F), according to recent bacterial viruses' taxonomy, based on their complete genome sequences. The 165,282 bp (PsgM02F) and 151,205 bp (PsgM04F) genomes do not feature any lysogenic-related (integrase) genes and, hence, can safely be assumed to follow a lytic lifestyle. While phage PsgM02F produced a morphogenesis yield of 124 virions per host cell, phage PsgM04F produced only 12 virions per host cell, indicating that they replicate well in Psg with a 50 min latency period. Genome mechanical analyses established a relationship between genome bendability and virion morphogenesis yield within infected host cells.

Complete genome sequence of Escherichia phage iGC_PHA_EC001

Antimicrobial resistance (AMR) in bacteria poses a global health emergency due to limited treatment options. Here, we report a lytic bacteriophage belonging to Stephanstirmvirinae family against an AMR Escherichia coli (ST2089). Escherichia phage iGC_PHA_EC001 is of genus Phapecoctavirus and 148,445 bp in length, encoding 269 predicted protein-coding sequences and 10 tRNAs. The phage encodes two lytic proteins containing phage_lysozyme (PF00959.22) and cell wall hydrolase_2 (PF07486.15) as catalytic domains, respectively.

Genome sequence of the Klebsiella quasipneumoniae bacteriophage EKq1 with activity against Klebsiella pneumoniae

We describe the genome of a lytic phage EKq1 isolated on Klebsiella quasipneumoniae, with activity against Klebsiella pneumoniae. EKq1 is an unclassified representative of the class Caudoviricetes, similar to Klebsiella phages VLCpiS8c, phiKp_7-2, and vB_KleS-HSE3. The 48,244-bp genome has a GC content of 56.43% and 63 predicted protein-coding genes.

Sheaths are diverse and abundant cell surface layers in archaea

Prokaryotic cells employ multiple protective layers crucial for defense, structural integrity, and cellular interactions in the environment. Archaea often feature an S-layer, with some species possessing additional and remarkably resistant sheaths. The archaeal sheath has been studied in Methanothrix and Methanospirillum, revealing a complex structure consisting of amyloid proteins organized into rings. Here, we conducted a comprehensive survey of sheath-forming proteins (SH proteins) across archaeal genomes. Structural modeling reveals a rich diversity of SH proteins, indicating the presence of a sheath in members of the TACK superphylum (Thermoprotei), as well as in the methanotrophic ANME-1. SH proteins are present in up to 40 copies per genome and display diverse domain arrangements suggesting multifunctional roles within the sheath, and potential involvement in cell-cell interaction with syntrophic partners. We uncover a complex evolutionary dynamic, indicating active exchange of SH proteins in archaeal communities. We find that viruses infecting sheathed archaea encode a diversity of SH-like proteins and we use them as markers to identify 580 vOTUs potentially associated with sheathed archaea. Structural modeling suggests that viral SH proteins can form complexes with the host SH proteins. We propose a previously unreported egress strategy where the expression of viral SH-like proteins may disrupt the integrity of the host sheath and facilitate viral exit during lysis. Together, our results significantly expand knowledge of the diversity and evolution of the archaeal sheath, which has been largely understudied but might have an important role in shaping microbial communities.

Simple Porifera holobiont reveals complex interactions between the host, an archaeon, a bacterium, and a phage

The basal metazoan phylum Porifera (sponges) is increasingly used as a model to investigate ecological and evolutionary features of microbe-animal symbioses. However, sponges often host complex microbiomes, which has hampered our understanding of their interactions with their microbial symbionts. Here, we describe the discovery and characterization of the simplest sponge holobiont reported to date, consisting of the deep-sea glass sponge Aphrocallistes beatrix and two newly-described microbial symbionts: an autotrophic ammonia-oxidizing archaeon and a bacterial heterotroph. Omics analyses and metabolic modeling revealed the dependency of the ammonia-oxidizing archaea on sponge-derived ammonia to drive primary production, which in turn supports the bacterium's growth by providing the dicarboxylate fumarate. Furthermore, virus-mediated archaeal lysis appears crucial to overcome the bacterium's vitamin B12 auxotrophy. These findings reveal that the exchanges of vitamin B12 and dicarboxylate may be evolutionarily conserved features of symbiosis as they can also be found in interactions between free-living marine bacteria, and between microbes and plants or diatoms.

Coordinated transcriptional response to environmental stress by a Synechococcus virus

Viruses are a major control on populations of microbes. Often, their virulence is examined in controlled laboratory conditions. Yet, in nature, environmental conditions lead to changes in host physiology and fitness that may impart both costs and benefits on viral success. Phosphorus (P) is a major abiotic control on the marine cyanobacterium Synechococcus. Some viruses infecting Synechococcus have acquired, from their host, a gene encoding a P substrate binding protein (PstS), thought to improve virus replication under phosphate starvation. Yet, pstS is uncommon among cyanobacterial viruses. Thus, we asked how infections with viruses lacking PstS are affected by P scarcity. We show that the production of infectious virus particles of such viruses is reduced in low P conditions. However, this reduction in progeny is not caused by impaired phage genome replication, thought to be a major sink for cellular phosphate. Instead, transcriptomic analysis showed that under low P conditions, a PstS-lacking cyanophage increased the expression of a specific gene set that included mazG, hli2, and gp43 encoding a pyrophosphatase, a high-light inducible protein and DNA polymerase, respectively. Moreover, several of the upregulated genes were controlled by the host's phoBR two-component system. We hypothesize that recycling and polymerization of nucleotides liberates free phosphate and thus allows viral morphogenesis, albeit at lower rates than when phosphate is replete or when phages encode pstS. Altogether, our data show how phage genomes, lacking obvious P-stress-related genes, have evolved to exploit their host's environmental sensing mechanisms to coordinate their own gene expression in response to resource limitation.

dbAPIS: a database of anti-prokaryotic immune system genes

Anti-prokaryotic immune system (APIS) proteins, typically encoded by phages, prophages, and plasmids, inhibit prokaryotic immune systems (e.g. restriction modification, toxin-antitoxin, CRISPR-Cas). A growing number of APIS genes have been characterized and dispersed in the literature. Here we developed dbAPIS (https://bcb.unl.edu/dbAPIS), as the first literature curated data repository for experimentally verified APIS genes and their associated protein families. The key features of dbAPIS include: (i) experimentally verified APIS genes with their protein sequences, functional annotation, PDB or AlphaFold predicted structures, genomic context, sequence and structural homologs from different microbiome/virome databases; (ii) classification of APIS proteins into sequence-based families and construction of hidden Markov models (HMMs); (iii) user-friendly web interface for data browsing by the inhibited immune system types or by the hosts, and functions for searching and batch downloading of pre-computed data; (iv) Inclusion of all types of APIS proteins (except for anti-CRISPRs) that inhibit a variety of prokaryotic defense systems (e.g. RM, TA, CBASS, Thoeris, Gabija). The current release of dbAPIS contains 41 verified APIS proteins and ∼4400 sequence homologs of 92 families and 38 clans. dbAPIS will facilitate the discovery of novel anti-defense genes and genomic islands in phages, by providing a user-friendly data repository and a web resource for an easy homology search against known APIS proteins.

PhageScope: a well-annotated bacteriophage database with automatic analyses and visualizations

Bacteriophages are viruses that infect bacteria or archaea. Understanding the diverse and intricate genomic architectures of phages is essential to study microbial ecosystems and develop phage therapy strategies. However, the existing phage databases are short of meticulous annotations. To this end, we propose PhageScope (https://phagescope.deepomics.org), an online phage database with comprehensive annotations. PhageScope harbors a collection of 873 718 phage sequences from various sources. Applying fifteen state-of-the-art tools to perform systematic annotations and analyses, PhageScope provides annotations on genome completeness, host range, lifestyle information, taxonomy classification, nine types of structural and functional genetic elements, and three types of comparative genomic studies for curated phages. Additionally, PhageScope incorporates automatic analyses and visualizations for curated and customized phages, serving as an efficient platform for phage study.

Large-scale genomic survey with deep learning-based method reveals strain-level phage specificity determinants

BACKGROUND

Phage therapy, reemerging as a promising approach to counter antimicrobial-resistant infections, relies on a comprehensive understanding of the specificity of individual phages. Yet the significant diversity within phage populations presents a considerable challenge. Currently, there is a notable lack of tools designed for large-scale characterization of phage receptor-binding proteins, which are crucial in determining the phage host range.

RESULTS

In this study, we present SpikeHunter, a deep learning method based on the ESM-2 protein language model. With SpikeHunter, we identified 231,965 diverse phage-encoded tailspike proteins, a crucial determinant of phage specificity that targets bacterial polysaccharide receptors, across 787,566 bacterial genomes from 5 virulent, antibiotic-resistant pathogens. Notably, 86.60% (143,200) of these proteins exhibited strong associations with specific bacterial polysaccharides. We discovered that phages with identical tailspike proteins can infect different bacterial species with similar polysaccharide receptors, underscoring the pivotal role of tailspike proteins in determining host range. The specificity is mainly attributed to the protein's C-terminal domain, which strictly correlates with host specificity during domain swapping in tailspike proteins. Importantly, our dataset-driven predictions of phage-host specificity closely match the phage-host pairs observed in real-world phage therapy cases we studied.

CONCLUSIONS

Our research provides a rich resource, including both the method and a database derived from a large-scale genomics survey. This substantially enhances understanding of phage specificity determinants at the strain level and offers a valuable framework for guiding phage selection in therapeutic applications.

Driving through stop signs: predicting stop codon reassignment improves functional annotation of bacteriophages

The majority of bacteriophage diversity remains uncharacterized, and new intriguing mechanisms of their biology are being continually described. Members of some phage lineages, such as the Crassvirales, repurpose stop codons to encode an amino acid by using alternate genetic codes. Here, we investigated the prevalence of stop codon reassignment in phage genomes and its subsequent impacts on functional annotation. We predicted 76 genomes within INPHARED and 712 vOTUs from the Unified Human Gut Virome Catalogue (UHGV) that repurpose a stop codon to encode an amino acid. We re-annotated these sequences with modified versions of Pharokka and Prokka, called Pharokka-gv and Prokka-gv, to automatically predict stop codon reassignment prior to annotation. Both tools significantly improved the quality of annotations, with Pharokka-gv performing best. For sequences predicted to repurpose TAG to glutamine (translation table 15), Pharokka-gv increased the median gene length (median of per genome median) from 287 to 481 bp for UHGV sequences (67.8% increase) and from 318 to 550 bp for INPHARED sequences (72.9% increase). The re-annotation increased median coding capacity from 66.8% to 90.0% and from 69.0% to 89.8% for UHGV and INPHARED sequences predicted to use translation table 15. Furthermore, the proportion of genes that could be assigned functional annotation increased, including an increase in the number of major capsid proteins that could be identified. We propose that automatic prediction of stop codon reassignment before annotation is beneficial to downstream viral genomic and metagenomic analyses.

Bacteriophage Taxonomy: A Continually Evolving Discipline

While taxonomy is an often underappreciated branch of science, it serves very important roles. Bacteriophage taxonomy has evolved from a discipline based mainly on morphology, characterized by the work of David Bradley and Hans-Wolfgang Ackermann, to the sequence-based approach that is taken today. The Bacterial Viruses Subcommittee of the International Committee on Taxonomy of Viruses (ICTV) takes a holistic approach to classifying prokaryote viruses by measuring overall DNA and protein similarity and phylogeny before making decisions about the taxonomic position of a new virus. The huge number of complete genomes being deposited with the National Center for Biotechnology Information (NCBI) and other public databases has resulted in a reassessment of the taxonomy of many viruses, and the future will see the introduction of new viral families and higher orders.

Comparative Analyses of Bacteriophage Genomes

Bacterial viruses (bacteriophages or phages) are the most abundant and diverse biological entities on Earth. There is a renewed worldwide interest in phage-centered research motivated by their enormous potential as antimicrobials to cope with multidrug-resistant pathogens. An ever-growing number of complete phage genomes are becoming available, derived either from newly isolated phages (cultivated phages) or recovered from metagenomic sequencing data (uncultivated phages). Robust comparative analysis is crucial for a comprehensive understanding of genotypic variations of phages and their related evolutionary processes, and to investigate the interaction mechanisms between phages and their hosts. In this chapter, we present a protocol for phage comparative genomics employing tools selected out of the many currently available, focusing on complete genomes of phages classified in the class Caudoviricetes. This protocol provides accurate identification of similarities, differences, and patterns among new and previously known complete phage genomes as well as phage clustering and taxonomic classification.

Bacteria and bacteriophage consortia are associated with protective intestinal metabolites in patients receiving stem cell transplantation

The microbiome is a predictor of clinical outcome in patients receiving allogeneic hematopoietic stem cell transplantation (allo-SCT). Microbiota-derived metabolites can modulate these outcomes. How bacteria, fungi and viruses contribute to the production of intestinal metabolites is still unclear. We combined amplicon sequencing, viral metagenomics and targeted metabolomics from stool samples of patients receiving allo-SCT (n = 78) and uncovered a microbiome signature of Lachnospiraceae and Oscillospiraceae and their associated bacteriophages, correlating with the production of immunomodulatory metabolites (IMMs). Moreover, we established the IMM risk index (IMM-RI), which was associated with improved survival and reduced relapse. A high abundance of short-chain fatty acid-biosynthesis pathways, specifically butyric acid via butyryl-coenzyme A (CoA):acetate CoA-transferase (BCoAT, which catalyzes EC 2.8.3.8) was detected in IMM-RI low-risk patients, and virome genome assembly identified two bacteriophages encoding BCoAT as an auxiliary metabolic gene. In conclusion, our study identifies a microbiome signature associated with protective IMMs and provides a rationale for considering metabolite-producing consortia and metabolite formulations as microbiome-based therapies.

Anaerobiosis, a neglected factor in phage-bacteria interactions

IMPORTANCE

Many parameters affect phage-bacteria interaction. Some of these parameters depend on the environment in which the bacteria are present. Anaerobiosis effect on phage infection in facultative anaerobic bacteria has not yet been studied. The absence of oxygen triggers metabolic changes in facultative bacteria and this affects phage infection and viral life cycle. Understanding how an anaerobic environment can alter the behavior of phages during infection is relevant for the phage therapy success.

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 64 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 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.

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.

Predicting stop codon reassignment improves functional annotation of bacteriophages

The majority of bacteriophage diversity remains uncharacterised, and new intriguing mechanisms of their biology are being continually described. Members of some phage lineages, such as the Crassvirales, repurpose stop codons to encode an amino acid by using alternate genetic codes. Here, we investigated the prevalence of stop codon reassignment in phage genomes and subsequent impacts on functional annotation. We predicted 76 genomes within INPHARED and 712 vOTUs from the Unified Human Gut Virome catalogue (UHGV) that repurpose a stop codon to encode an amino acid. We re-annotated these sequences with modified versions of Pharokka and Prokka, called Pharokka-gv and Prokka-gv, to automatically predict stop codon reassignment prior to annotation. Both tools significantly improved the quality of annotations, with Pharokka-gv performing best. For sequences predicted to repurpose TAG to glutamine (translation table 15), Pharokka-gv increased the median gene length (median of per genome medians) from 287 to 481 bp for UHGV sequences (67.8% increase) and from 318 to 550 bp for INPHARED sequences (72.9% increase). The re-annotation increased mean coding density from 66.8% to 90.0%, and from 69.0% to 89.8% for UHGV and INPHARED sequences. Furthermore, the proportion of genes that could be assigned functional annotation increased, including an increase in the number of major capsid proteins that could be identified. We propose that automatic prediction of stop codon reassignment before annotation is beneficial to downstream viral genomic and metagenomic analyses.

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 were more abundant and species rich 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.

Genome Sequence and Characterization of Coliphage Môr Ffagbaw

Background

The quality of coastal waters around the United Kingdom is an area of increasing concern following sewer overflow, where wastewater is discharged into the environment. Coliphages, viruses that infect coliform bacteria, are associated with water quality in aquatic systems, yet remain largely uncharacterized at the genomic level.

Materials and Methods

Phage môr ffagbaw was isolated from seawater against Escherichia coli by enrichment and plaque assays. Whole genome sequencing, transmission electron microscopy, and host range analysis against the E. coli reference (ECOR) collection were used to characterize the phage.

Results

The virion had a siphovirus morphology and genomic analysis placed it within the family Drexlerviridae, subfamily Tempevirinae, and forms a new species within the genus Hanrivervirus. Spot assays revealed that phage môr ffagbaw could form plaques on 6 out of 72 ECOR strains (8%).

Conclusions

Môr ffagbaw represents a new species of phage within the genus Hanrivervirus, with a narrow host range.

Identification of Huge Phages from Wastewater Metagenomes

Huge phages have genomes larger than 200 kilobases, which are particularly interesting for their genetic inventory and evolution. We screened 165 wastewater metagenomes for the presence of viral sequences. After identifying over 600 potential huge phage genomes, we reduced the dataset using manual curation by excluding viral contigs that did not contain viral protein-coding genes or consisted of concatemers of several small phage genomes. This dataset showed seven fully annotated huge phage genomes. The phages grouped into distinct phylogenetic clades, likely forming new genera and families. A phylogenomic analysis between our huge phages and phages with smaller genomes, i.e., less than 200 kb, supported the hypothesis that huge phages have undergone convergent evolution. The genomes contained typical phage protein-coding genes, sequential gene cassettes for metabolic pathways, and complete inventories of tRNA genes covering all standard and rare amino acids. Our study showed a pipeline for huge phage analyses that may lead to new enzymes for therapeutic or biotechnological applications.

Ongoing shuffling of protein fragments diversifies core viral functions linked to interactions with bacterial hosts

Biological modularity enhances evolutionary adaptability. This principle is vividly exemplified by bacterial viruses (phages), which display extensive genomic modularity. Phage genomes are composed of independent functional modules that evolve separately and recombine in various configurations. While genomic modularity in phages has been extensively studied, less attention has been paid to protein modularity-proteins consisting of distinct building blocks that can evolve and recombine, enhancing functional and genetic diversity. Here, we use a set of 133,574 representative phage proteins and highly sensitive homology detection to capture instances of domain mosaicism, defined as fragment sharing between two otherwise unrelated proteins, and to understand its relationship with functional diversity in phage genomes. We discover that unrelated proteins from diverse functional classes frequently share homologous domains. This phenomenon is particularly pronounced within receptor-binding proteins, endolysins, and DNA polymerases. We also identify multiple instances of recent diversification via domain shuffling in receptor-binding proteins, neck passage structures, endolysins and some members of the core replication machinery, often transcending distant taxonomic and ecological boundaries. Our findings suggest that ongoing diversification via domain shuffling is reflective of a co-evolutionary arms race, driven by the need to overcome various bacterial resistance mechanisms against phages.

A compendium of viruses from methanogenic archaea reveals their diversity and adaptations to the gut environment

Methanogenic archaea are major producers of methane, a potent greenhouse gas and biofuel, and are widespread in diverse environments, including the animal gut. The ecophysiology of methanogens is likely impacted by viruses, which remain, however, largely uncharacterized. Here we carried out a global investigation of viruses associated with all current diversity of methanogens by assembling an extensive CRISPR database consisting of 156,000 spacers. We report 282 high-quality (pro)viral and 205 virus-like/plasmid sequences assigned to hosts belonging to ten main orders of methanogenic archaea. Viruses of methanogens can be classified into 87 families, underscoring a still largely undiscovered genetic diversity. Viruses infecting gut-associated archaea provide evidence of convergence in adaptation with viruses infecting gut-associated bacteria. These viruses contain a large repertoire of lysin proteins that cleave archaeal pseudomurein and are enriched in glycan-binding domains (Ig-like/Flg_new) and diversity-generating retroelements. The characterization of this vast repertoire of viruses paves the way towards a better understanding of their role in regulating methanogen communities globally, as well as the development of much-needed genetic tools.

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, and 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, 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' 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 1,072 Pseudomonas spp. genomes. VIBES software is available at https://github.com/TravisWheelerLab/VIBES.

Unveiling the hidden role of aquatic viruses in hydrocarbon pollution bioremediation

Hydrocarbon pollution poses substantial environmental risks to water and soil. Bioremediation, which utilizes microorganisms to manage pollutants, offers a cost-effective solution. However, the role of viruses, particularly bacteriophages (phages), in bioremediation remains unexplored. This study examines the diversity and activity of hydrocarbon-degradation genes encoded by environmental viruses, focusing on phages, within public databases. We identified 57 high-quality phage-encoded auxiliary metabolic genes (AMGs) related to hydrocarbon degradation, which we refer to as virus-encoded hydrocarbon degradation genes (vHYDEGs). These genes are encoded by taxonomically diverse aquatic phages and highlight the under-characterized global virosphere. Six protein families involved in the initial alkane hydroxylation steps were identified. Phylogenetic analyses revealed the diverse evolutionary trajectories of vHYDEGs across habitats, revealing previously unknown biodegraders linked evolutionarily with vHYDEGs. Our findings suggest phage AMGs may contribute to alkane and aromatic hydrocarbon degradation, participating in the initial, rate-limiting hydroxylation steps, thereby aiding hydrocarbon pollution bioremediation and promoting their propagation. To support future research, we developed vHyDeg, a database containing identified vHYDEGs with comprehensive annotations, facilitating the screening of hydrocarbon degradation AMGs and encouraging their bioremediation applications.

Phables: from fragmented assemblies to high-quality bacteriophage genomes

MOTIVATION

Microbial communities have a profound impact on both human health and various environments. Viruses infecting bacteria, known as bacteriophages or phages, play a key role in modulating bacterial communities within environments. High-quality phage genome sequences are essential for advancing our understanding of phage biology, enabling comparative genomics studies and developing phage-based diagnostic tools. Most available viral identification tools consider individual sequences to determine whether they are of viral origin. As a result of challenges in viral assembly, fragmentation of genomes can occur, and existing tools may recover incomplete genome fragments. Therefore, the identification and characterization of novel phage genomes remain a challenge, leading to the need of improved approaches for phage genome recovery.

RESULTS

We introduce Phables, a new computational method to resolve phage genomes from fragmented viral metagenome assemblies. Phables identifies phage-like components in the assembly graph, models each component as a flow network, and uses graph algorithms and flow decomposition techniques to identify genomic paths. Experimental results of viral metagenomic samples obtained from different environments show that Phables recovers on average over 49% more high-quality phage genomes compared to existing viral identification tools. Furthermore, Phables can resolve variant phage genomes with over 99% average nucleotide identity, a distinction that existing tools are unable to make.

AVAILABILITY AND IMPLEMENTATION

Phables is available on GitHub at https://github.com/Vini2/phables.

Isolation and genome sequencing of five lytic bacteriophages from hospital wastewater in the Philippines

Here, we report the genome sequences of five bacteriophages isolated from hospital wastewater, including two new species and two candidates for therapeutic application. No virulence, temperate marker and antibiotic resistance genes were found in the genomes of Escherichia phage vB_VIPECOOM03 and Klebsiella phage vB_VIPKPNUMC01, making them suitable candidate for therapy.

Complete genome sequence of the broad host range Acinetobacter baumannii phage EAb13

We describe the genome of a lytic phage EAb13 isolated from sewage, with broad activity against multidrug-resistant Acinetobacter baumannii. EAb13 is an unclassified siphovirus. Its genome consists of 82,411 bp, with 40.15% GC content, 126 protein-coding sequences, 1 tRNA, and 2,177 bp-long direct terminal repeats.

Phables: from fragmented assemblies to high-quality bacteriophage genomes

Microbial communities influence both human health and different environments. Viruses infecting bacteria, known as bacteriophages or phages, play a key role in modulating bacterial communities within environments. High-quality phage genome sequences are essential for advancing our understanding of phage biology, enabling comparative genomics studies, and developing phage-based diagnostic tools. Most available viral identification tools consider individual sequences to determine whether they are of viral origin. As a result of the challenges in viral assembly, fragmentation of genomes can occur, leading to the need for new approaches in viral identification. Therefore, the identification and characterisation of novel phages remain a challenge. We introduce Phables, a new computational method to resolve phage genomes from fragmented viral metagenome assemblies. Phables identifies phage-like components in the assembly graph, models each component as a flow network, and uses graph algorithms and flow decomposition techniques to identify genomic paths. Experimental results of viral metagenomic samples obtained from different environments show that Phables recovers on average over 49% more high-quality phage genomes compared to existing viral identification tools. Furthermore, Phables can resolve variant phage genomes with over 99% average nucleotide identity, a distinction that existing tools are unable to make. Phables is available on GitHub at https://github.com/Vini2/phables.

Novel Escherichia coli-Infecting Bacteriophages Isolated from Uganda That Target Human Clinical Isolates

Background

The antimicrobial resistance catastrophe is a growing global health threat and predicted to be worse in developing countries. Phages for Global Health (PGH) is training scientists in these regions to isolate relevant therapeutic phages for pathogenic bacteria within their locality, and thus contributing to making phage technology universally available.

Materials and Methods

During the inaugural PGH workshop in East Africa, samples from Ugandan municipal sewage facilities were collected and two novel Escherichia coli lytic phages were isolated and characterized.

Results

The phages, UP19 (capsid diameter ∼100 nm, contractile tail ∼120/20 nm) and UP30 (capsid diameter ∼70 nm, noncontractile tail of ∼170/20 nm), lysed ∼82% and ∼36% of the 11 clinical isolates examined, respectively. The genomes of UP19 (171.402 kb, 282 CDS) and UP30 (49.834 kb, 75 CDS) closely match the genera Dhakavirus and Tunavirus, respectively.

Conclusion

The phages isolated have therapeutic potential for further development against E. coli infections.

Identifying and tracking mobile elements in evolving compost communities yields insights into the nanobiome

Microbial evolution is driven by rapid changes in gene content mediated by horizontal gene transfer (HGT). While mobile genetic elements (MGEs) are important drivers of gene flux, the nanobiome-the zoo of Darwinian replicators that depend on microbial hosts-remains poorly characterised. New approaches are necessary to increase our understanding beyond MGEs shaping individual populations, towards their impacts on complex microbial communities. A bioinformatic pipeline (xenoseq) was developed to cross-compare metagenomic samples from microbial consortia evolving in parallel, aimed at identifying MGE dissemination, which was applied to compost communities which underwent periodic mixing of MGEs. We show that xenoseq can distinguish movement of MGEs from demographic changes in community composition that otherwise confounds identification, and furthermore demonstrate the discovery of various unexpected entities. Of particular interest was a nanobacterium of the candidate phylum radiation (CPR) which is closely related to a species identified in groundwater ecosystems (Candidatus Saccharibacterium), and appears to have a parasitic lifestyle. We also highlight another prolific mobile element, a 313 kb plasmid hosted by a Cellvibrio lineage. The host was predicted to be capable of nitrogen fixation, and acquisition of the plasmid coincides with increased ammonia production. Taken together, our data show that new experimental strategies combined with bioinformatic analyses of metagenomic data stand to provide insight into the nanobiome as a driver of microbial community evolution.

VAPEX: an interactive web server for the deep exploration of natural virus and phage genomes

MOTIVATION

Studying the genetic makeup of viruses and phages through genome analysis is crucial for comprehending their function in causing diseases, progressing medicine, tracing their evolutionary history, monitoring the environment, and creating innovative biotechnologies. However, accessing the necessary data can be challenging due to a lack of dedicated comparative genomic tools and viral and phage databases, which are often outdated. Moreover, many wet bench experimentalists may not have the computational proficiency required to manipulate large amounts of genomic data.

RESULTS

We have developed VAPEX (Virus And Phage EXplorer), a web server which is supported by a database and features a user-friendly web interface. This tool enables users to easily perform various genomic analysis queries on all natural viruses and phages that have been fully sequenced and are listed in the NCBI compendium. VAPEX therefore excels in producing visual depictions of fully resolved synteny maps, which is one of its key strengths. VAPEX has the ability to exhibit a vast array of orthologous gene classes simultaneously through the use of symbolic representation. Additionally, VAPEX can fully analyze user-submitted viral and phage genomes, including those that have not yet been annotated.

AVAILABILITY AND IMPLEMENTATION

VAPEX can be accessed from all current web browsers such as Chrome, Firefox, Edge, Safari, and Opera. VAPEX is freely accessible at https://archaea.i2bc.paris-saclay.fr/vapex/.

Phages are unrecognized players in the ecology of the oral pathogen Porphyromonas gingivalis

BACKGROUND

Porphyromonas gingivalis (hereafter "Pg") is an oral pathogen that has been hypothesized to act as a keystone driver of inflammation and periodontal disease. Although Pg is most readily recovered from individuals with actively progressing periodontal disease, healthy individuals and those with stable non-progressing disease are also colonized by Pg. Insights into the factors shaping the striking strain-level variation in Pg, and its variable associations with disease, are needed to achieve a more mechanistic understanding of periodontal disease and its progression. One of the key forces often shaping strain-level diversity in microbial communities is infection of bacteria by their viral (phage) predators and symbionts. Surprisingly, although Pg has been the subject of study for over 40 years, essentially nothing is known of its phages, and the prevailing paradigm is that phages are not important in the ecology of Pg.

RESULTS

Here we systematically addressed the question of whether Pg are infected by phages-and we found that they are. We found that prophages are common in Pg, they are genomically diverse, and they encode genes that have the potential to alter Pg physiology and interactions. We found that phages represent unrecognized targets of the prevalent CRISPR-Cas defense systems in Pg, and that Pg strains encode numerous additional mechanistically diverse candidate anti-phage defense systems. We also found that phages and candidate anti-phage defense system elements together are major contributors to strain-level diversity and the species pangenome of this oral pathogen. Finally, we demonstrate that prophages harbored by a model Pg strain are active in culture, producing extracellular viral particles in broth cultures.

CONCLUSION

This work definitively establishes that phages are a major unrecognized force shaping the ecology and intra-species strain-level diversity of the well-studied oral pathogen Pg. The foundational phage sequence datasets and model systems that we establish here add to the rich context of all that is already known about Pg, and point to numerous avenues of future inquiry that promise to shed new light on fundamental features of phage impacts on human health and disease broadly. Video Abstract.

Protein embeddings improve phage-host interaction prediction

With the growing interest in using phages to combat antimicrobial resistance, computational methods for predicting phage-host interactions have been explored to help shortlist candidate phages. Most existing models consider entire proteomes and rely on manual feature engineering, which poses difficulty in selecting the most informative sequence properties to serve as input to the model. In this paper, we framed phage-host interaction prediction as a multiclass classification problem that takes as input the embeddings of a phage's receptor-binding proteins, which are known to be the key machinery for host recognition, and predicts the host genus. We explored different protein language models to automatically encode these protein sequences into dense embeddings without the need for additional alignment or structural information. We show that the use of embeddings of receptor-binding proteins presents improvements over handcrafted genomic and protein sequence features. The highest performance was obtained using the transformer-based protein language model ProtT5, resulting in a 3% to 4% increase in weighted F1 and recall scores across different prediction confidence thresholds, compared to using selected handcrafted sequence features.

Newly Established Genetic System for Functional Analysis of MetSV

The linear chromosome of the Methanosarcina spherical virus with 10,567 bp exhibits 22 ORFs with mostly unknown functions. Annotation using common tools and databases predicted functions for a few genes like the type B DNA polymerase (MetSVORF07) or the small (MetSVORF15) and major (MetSVORF16) capsid proteins. For verification of assigned functions of additional ORFs, biochemical or genetic approaches were found to be essential. Consequently, we established a genetic system for MetSV by cloning its genome into the E. coli plasmid pCR-XL-2. Comparisons of candidate plasmids with the MetSV reference based on Nanopore sequencing revealed several mutations of yet unknown provenance with an impact on protein-coding sequences. Linear MetSV inserts were generated by BamHI restriction, purified and transformed in Methanosarcina mazei by an optimized liposome-mediated transformation protocol. Analysis of resulting MetSV virions by TEM imaging and infection experiments demonstrated no significant differences between plasmid-born viruses and native MetSV particles regarding their morphology or lytic behavior. The functionality of the genetic system was tested by the generation of a ΔMetSVORF09 mutant that was still infectious. Our genetic system of MetSV, the first functional system for a virus of methanoarchaea, now allows us to obtain deeper insights into MetSV protein functions and virus-host interactions.

Unexplored diversity and ecological functions of transposable phages

Phages are prevalent in diverse environments and play major ecological roles attributed to their tremendous diversity and abundance. Among these viruses, transposable phages (TBPs) are exceptional in terms of their unique lifestyle, especially their replicative transposition. Although several TBPs have been isolated and the life cycle of the representative phage Mu has been extensively studied, the diversity distribution and ecological functions of TBPs on the global scale remain unknown. Here, by mining TBPs from enormous microbial genomes and viromes, we established a TBP genome dataset (TBPGD), that expands the number of accessible TBP genomes 384-fold. TBPs are prevalent in diverse biomes and show great genetic diversity. Based on taxonomic evaluations, we propose the categorization of TBPs into four viral groups, including 11 candidate subfamilies. TBPs infect multiple bacterial phyla, and seem to infect a wider range of hosts than non-TBPs. Diverse auxiliary metabolic genes (AMGs) are identified in the TBP genomes, and genes related to glycoside hydrolases and pyrimidine deoxyribonucleotide biosynthesis are highly enriched. Finally, the influences of TBPs on their hosts are experimentally examined by using the marine bacterium Shewanella psychrophila WP2 and its infecting transposable phage SP2. Collectively, our findings greatly expand the genetic diversity of TBPs, and comprehensively reveal their potential influences in various ecosystems.

New Bacteriophages with Podoviridal Morphotypes Active against Yersinia pestis: Characterization and Application Potential

Phages of highly pathogenic bacteria represent an area of growing interest for bacterial detection and identification and subspecies typing, as well as for phage therapy and environmental decontamination. Eight new phages-YpEc56, YpEc56D, YpEc57, YpEe58, YpEc1, YpEc2, YpEc11, and YpYeO9-expressing lytic activity towards Yersinia pestis revealed a virion morphology consistent with the Podoviridae morphotype. These phages lyse all 68 strains from 2 different sets of Y. pestis isolates, thus limiting their potential application for subtyping of Y. pestis strains but making them rather promising in terms of infection control. Two phages-YpYeO9 and YpEc11-were selected for detailed studies based on their source of isolation and lytic cross activity towards other Enterobacteriaceae. The full genome sequencing demonstrated the virulent nature of new phages. Phage YpYeO9 was identified as a member of the Teseptimavirus genus and YpEc11 was identified as a member of the Helsettvirus genus, thereby representing new species. A bacterial challenge assay in liquid microcosm with a YpYeO9/YpEc11 phage mixture showed elimination of Y. pestis EV76 during 4 h at a P/B ratio of 1000:1. These results, in combination with high lysis stability results of phages in liquid culture, the low frequency of formation of phage resistant mutants, and their viability under different physical-chemical factors indicate their potential for their practical use as an antibacterial mean.

Genomic analysis of Anderson typing phages of Salmonella Typhimrium: towards understanding the basis of bacteria-phage interaction

The Anderson phage typing scheme has been successfully used worldwide for epidemiological surveillance of Salmonella enterica serovar Typhimurium. Although the scheme is being replaced by whole genome sequence subtyping methods, it can provide a valuable model system for study of phage-host interaction. The phage typing scheme distinguishes more than 300 definitive types of Salmonella Typhimurium based on their patterns of lysis to a unique collection of 30 specific Salmonella phages. In this study, we sequenced the genomes of 28 Anderson typing phages of Salmonella Typhimurium to begin to characterize the genetic determinants that are responsible for the differences in these phage type profiles. Genomic analysis of typing phages reveals that Anderson phages can be classified into three different groups, the P22-like, ES18-like and SETP3-like clusters. Most Anderson phages are short tailed P22-like viruses (genus Lederbergvirus); but phages STMP8 and STMP18 are very closely related to the lambdoid long tailed phage ES18, and phages STMP12 and STMP13 are related to the long noncontractile tailed, virulent phage SETP3. Most of these typing phages have complex genome relationships, but interestingly, two phage pairs STMP5 and STMP16 as well as STMP12 and STMP13 differ by a single nucleotide. The former affects a P22-like protein involved in DNA passage through the periplasm during its injection, and the latter affects a gene whose function is unknown. Using the Anderson phage typing scheme would provide insights into phage biology and the development of phage therapy for the treatment of antibiotic resistant bacterial infections.

The Outer Membrane Cytochrome OmcA Is Essential for Infection of Shewanella oneidensis by a Zebrafish-Associated Bacteriophage

The microbiota-the mixture of microorganisms in the intestinal tract of animals-plays an important role in host biology. Bacteriophages are a prominent, though often overlooked, component of the microbiota. The mechanisms that phage use to infect susceptible cells associated with animal hosts, and the broader role they could play in determining the substituents of the microbiota, are poorly understood. In this study, we isolated a zebrafish-associated bacteriophage, which we named Shewanella phage FishSpeaker. This phage infects Shewanella oneidensis strain MR-1, which cannot colonize zebrafish, but it is unable to infect Shewanella xiamenensis strain FH-1, a strain isolated from the zebrafish gut. Our data suggest that FishSpeaker uses the outer membrane decaheme cytochrome OmcA, which is an accessory component of the extracellular electron transfer (EET) pathway in S. oneidensis, as well as the flagellum to recognize and infect susceptible cells. In a zebrafish colony that lacks detectable FishSpeaker, we found that most Shewanella spp. are sensitive to infection and that some strains are resistant to infection. Our results suggest that phage could act as a selectivity filter for zebrafish-associated Shewanella and show that the EET machinery can be targeted by phage in the environment. IMPORTANCE Phage exert selective pressure on bacteria that influences and shapes the composition of microbial populations. However, there is a lack of native, experimentally tractable systems for studying how phage influence microbial population dynamics in complex communities. Here, we show that a zebrafish-associated phage requires both the outer membrane-associated extracellular electron transfer protein OmcA and the flagellum to infect Shewanella oneidensis strain MR-1. Our results suggest that the newly discovered phage-FishSpeaker-could exert selective pressure that restricts which Shewanella spp. colonize zebrafish. Moreover, the requirement of OmcA for infection by FishSpeaker suggests that the phage preferentially infects cells that are oxygen limited, a condition required for OmcA expression and an ecological feature of the zebrafish gut.

Genome Sequence of Staphylococcus aureus Phage ESa2

We describe the genome of a lytic phage, ESa2, isolated from environmental water and specific for Staphylococcus aureus. ESa2 belongs to the family Herelleviridae and genus Kayvirus. Its genome consists of 141,828 bp, with 30.25% GC content, 253 predicted protein-coding sequences, 3 tRNAs, and 10,130-bp-long terminal repeats.

Deciphering Phage-Host Specificity Based on the Association of Phage Depolymerases and Bacterial Surface Glycan with Deep Learning

Phage tailspike proteins are depolymerases that target diverse bacterial surface glycans with high specificity, determining the host-specificity of numerous phages. To address the challenge of identifying tailspike proteins due to their sequence diversity, we developed SpikeHunter, an approach based on the ESM-2 protein language model. Using SpikeHunter, we successfully identified 231,965 tailspike proteins from a dataset comprising 8,434,494 prophages found within 165,365 genomes of five common pathogens. Among these proteins, 143,035 tailspike proteins displayed strong associations with serotypes. Moreover, we observed highly similar tailspike proteins in species that share closely related serotypes. We found extensive domain swapping in all five species, with the C-terminal domain being significantly associated with host serotype highlighting its role in host range determination. Our study presents a comprehensive cross-species analysis of tailspike protein to serotype associations, providing insights applicable to phage therapy and biotechnology.

Isolation and Characterization of a Novel Phage Collection against Avian-Pathogenic Escherichia coli

The increase in antibiotic-resistant avian-pathogenic Escherichia coli (APEC), the causative agent of colibacillosis in poultry, warrants urgent research and the development of alternative therapies. This study describes the isolation and characterization of 19 genetically diverse, lytic coliphages, 8 of which were tested in combination for their efficacy in controlling in ovo APEC infections. Genome homology analysis revealed that the phages belong to nine different genera, one of them being a novel genus (Nouzillyvirus). One phage, REC, was derived from a recombination event between two Phapecoctavirus phages (ESCO5 and ESCO37) isolated in this study. Twenty-six of the 30 APEC strains tested were lysed by at least one phage. Phages exhibited varying infectious capacities, with narrow to broad host ranges. The broad host range of some phages could be partially explained by the presence of receptor-binding protein carrying a polysaccharidase domain. To demonstrate their therapeutic potential, a phage cocktail consisting of eight phages belonging to eight different genera was tested against BEN4358, an APEC O2 strain. In vitro, this phage cocktail fully inhibited the growth of BEN4358. In a chicken lethality embryo assay, the phage cocktail enabled 90% of phage-treated embryos to survive infection with BEN4358, compared with 0% of nontreated embryos, indicating that these novel phages are good candidates to successfully treat colibacillosis in poultry. IMPORTANCE Colibacillosis, the most common bacterial disease affecting poultry, is mainly treated by antibiotics. Due to the increased prevalence of multidrug-resistant avian-pathogenic Escherichia coli, there is an urgent need to assess the efficacy of alternatives to antibiotherapy, such as phage therapy. Here, we have isolated and characterized 19 coliphages that belong to nine phage genera. We showed that a combination of 8 of these phages was efficacious in vitro to control the growth of a clinical isolate of E. coli. Used in ovo, this phage combination allowed embryos to survive APEC infection. Thus, this phage combination represents a promising treatment for avian colibacillosis.