INfrastructure for a PHAge REference Database: Identification of Large-Scale Biases in the Current Collection of Cultured Phage Genomes

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.

Deeper Exploration of Gut Microbiome: Profile of Resistome, Virome and Viral Auxiliary Metabolic Genes of Three Ethnic Indian Groups

The current study explored the resistomes and viromes of three Indian ethnic populations: Jaisalmer, Khargone, and Ladakh. These three groups had different dietary habits and antibiotic consumption rates. A resistome analysis indicated that compared to the Jaisalmer (n = 10) group, the burden of antibiotic resistance genes in the gut microbiome was higher in the Khargone (n = 12) and Ladakh (n = 9) groups. However, correlational analysis factoring in food habits, healthcare, and economic status was not statistically significant due to the limited number of samples. A considerable number of antibiotic resistance genes (ARGs) were present in well-known gut commensals such as Bifidobacteriaceae, Acidomonococcaceae, etc., as retrieved directly by mapping to the Resfinder database using the Groot tool. Further, the raw reads were assembled using MEGAHIT, and putative bacteriophages were retrieved using the VIBRANT tool. Many of the classified bacteriophages of the virome revealed that bacteria belonging to the families Bifidobacteriaceae and Enterocococcaceae were their hosts. The prophages identified in these groups primarily contained auxiliary metabolic genes (AMGs) for primary amino acid metabolism. However, there were significantly fewer AMGs in the Ladakh group than in the Jaisalmer group (p < 0.05). None of the classified bacteriophages or prophages contained ARGs. This indicates that phages do not normally carry antibiotic resistance genes.

An enterococcal phage protein inhibits type IV restriction enzymes involved in antiphage defense

The prevalence of multidrug resistant (MDR) bacterial infections continues to rise as the development of antibiotics needed to combat these infections remains stagnant. MDR enterococci are a major contributor to this crisis. A potential therapeutic approach for combating MDR enterococci is bacteriophage (phage) therapy, which uses lytic viruses to infect and kill pathogenic bacteria. While phages that lyse some strains of MDR enterococci have been identified, other strains display high levels of resistance and the mechanisms underlying this resistance are poorly defined. Here, we use a CRISPR interference (CRISPRi) screen to identify a genetic locus found on a mobilizable plasmid from Enterococcus faecalis involved in phage resistance. This locus encodes a putative serine recombinase followed by a Type IV restriction enzyme (TIV-RE) that we show restricts the replication of phage phi47 in vancomycin-resistant E. faecalis. We further find that phi47 evolves to overcome restriction by acquiring a missense mutation in a TIV-RE inhibitor protein. We show that this inhibitor, termed type IV restriction inhibiting factor A (tifA), binds and inactivates diverse TIV-REs. Overall, our findings advance our understanding of phage defense in drug-resistant E. faecalis and provide mechanistic insight into how phages evolve to overcome antiphage defense systems.

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.

An enterococcal phage protein broadly inhibits type IV restriction enzymes involved in antiphage defense

The prevalence of multidrug resistant (MDR) bacterial infections continues to rise as the development of antibiotics needed to combat these infections remains stagnant. MDR enterococci are a major contributor to this crisis. A potential therapeutic approach for combating MDR enterococci is bacteriophage (phage) therapy, which uses lytic viruses to infect and kill pathogenic bacteria. While phages that lyse some strains of MDR enterococci have been identified, other strains display high levels of resistance and the mechanisms underlying this resistance are poorly defined. Here, we use a CRISPR interference (CRISPRi) screen to identify a genetic locus found on a mobilizable plasmid from Enterococcus faecalis involved in phage resistance. This locus encodes a putative serine recombinase followed by a Type IV restriction enzyme (TIV-RE) that we show restricts the replication of phage phi47 in E. faecalis. We further find that phi47 evolves to overcome restriction by acquiring a missense mutation in a TIV-RE inhibitor protein. We show that this inhibitor, termed type IV restriction inhibiting factor A (tifA), binds and inactivates diverse TIV-REs. Overall, our findings advance our understanding of phage defense in drug-resistant E. faecalis and provide mechanistic insight into how phages evolve to overcome antiphage defense systems.

Characterization and Abundance of Plasmid-Dependent Alphatectivirus Bacteriophages

Antimicrobial resistance (AMR) is a major public health threat, exacerbated by the ability of bacteria to rapidly disseminate antimicrobial resistance genes (ARG). Since conjugative plasmids of the incompatibility group P (IncP) are ubiquitous mobile genetic elements that often carry ARG and are broad-host-range, they are important targets to prevent the dissemination of AMR. Plasmid-dependent phages infect plasmid-carrying bacteria by recognizing components of the conjugative secretion system as receptors. We sought to isolate plasmid-dependent phages from wastewater using an avirulent strain of Salmonella enterica carrying the conjugative IncP plasmid pKJK5. Irrespective of the site, we only obtained bacteriophages belonging to the genus Alphatectivirus. Eleven isolates were sequenced, their genomes analyzed, and their host range established using S. enterica, Escherichia coli, and Pseudomonas putida carrying diverse conjugative plasmids. We confirmed that Alphatectivirus are abundant in domestic and hospital wastewater using culture-dependent and culture-independent approaches. However, these results are not consistent with their low or undetectable occurrence in metagenomes. Therefore, overall, our results emphasize the importance of performing phage isolation to uncover diversity, especially considering the potential of plasmid-dependent phages to reduce the spread of ARG carried by conjugative plasmids, and to help combat the AMR crisis.

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.

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.

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.

Novel phages of Pseudomonas syringae unveil numerous potential auxiliary metabolic genes

Relatively few phages that infect plant pathogens have been isolated and investigated. The Pseudomonas syringae species complex is present in various environments, including plants. It can cause major crop diseases, such as bacterial canker on apricot trees. This study presents a collection of 25 unique phage genomes that infect P. syringae. These phages were isolated from apricot orchards with bacterial canker symptoms after enrichment with 21 strains of P. syringae. This collection comprises mostly virulent phages, with only three being temperate. They belong to 14 genera, 11 of which are newly discovered, and 18 new species, revealing great genetic diversity within this collection. Novel DNA packaging systems have been identified bioinformatically in one of the new phage species, but experimental confirmation is required to define the precise mechanism. Additionally, many phage genomes contain numerous potential auxiliary metabolic genes with diversified putative functions. At least three phages encode genes involved in bacterial tellurite resistance, a toxic metalloid. This suggests that viruses could play a role in bacterial stress tolerance. This research emphasizes the significance of continuing the search for new phages in the agricultural ecosystem to unravel novel ecological diversity and new gene functions. This work contributes to the foundation for future fundamental and applied research on phages infecting phytopathogenic bacteria.

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.

Phage predation, disease severity, and pathogen genetic diversity in cholera patients

Despite an increasingly detailed picture of the molecular mechanisms of bacteriophage (phage)-bacterial interactions, we lack an understanding of how these interactions evolve and impact disease within patients. In this work, we report a year-long, nationwide study of diarrheal disease patients in Bangladesh. Among cholera patients, we quantified Vibrio cholerae (prey) and its virulent phages (predators) using metagenomics and quantitative polymerase chain reaction while accounting for antibiotic exposure using quantitative mass spectrometry. Virulent phage (ICP1) and antibiotics suppressed V. cholerae to varying degrees and were inversely associated with severe dehydration depending on resistance mechanisms. In the absence of antiphage defenses, predation was "effective," with a high predator:prey ratio that correlated with increased genetic diversity among the prey. In the presence of antiphage defenses, predation was "ineffective," with a lower predator:prey ratio that correlated with increased genetic diversity among the predators. Phage-bacteria coevolution within patients should therefore be considered in the deployment of phage-based therapies and diagnostics.

Discovery and description of novel phage genomes from urban microbiomes sampled by the MetaSUB consortium

Bacteriophages are recognized as the most abundant members of microbiomes and have therefore a profound impact on microbial communities through the interactions with their bacterial hosts. The International Metagenomics and Metadesign of Subways and Urban Biomes Consortium (MetaSUB) has sampled mass-transit systems in 60 cities over 3 years using metagenomics, throwing light into these hitherto largely unexplored urban environments. MetaSUB focused primarily on the bacterial community. In this work, we explored MetaSUB metagenomic data in order to recover and analyze bacteriophage genomes. We recovered and analyzed 1714 phage genomes with size at least 40 kbp, from the class Caudoviricetes, the vast majority of which (80%) are novel. The recovered genomes were predicted to belong to temperate (69%) and lytic (31%) phages. Thirty-three of these genomes have more than 200 kbp, and one of them reaches 572 kbp, placing it among the largest phage genomes ever found. In general, the phages tended to be site-specific or nearly so, but 194 genomes could be identified in every city from which phage genomes were retrieved. We predicted hosts for 48% of the phages and observed general agreement between phage abundance and the respective bacterial host abundance, which include the most common nosocomial multidrug-resistant pathogens. A small fraction of the phage genomes are carriers of antibiotic resistance genes, and such genomes tended to be particularly abundant in the sites where they were found. We also detected CRISPR-Cas systems in five phage genomes. This study expands the previously reported MetaSUB results and is a contribution to the knowledge about phage diversity, global distribution, and phage genome content.

Nanopore and Illumina sequencing reveal different viral populations from human gut samples

The advent of viral metagenomics, or viromics, has improved our knowledge and understanding of global viral diversity. High-throughput sequencing technologies enable explorations of the ecological roles, contributions to host metabolism, and the influence of viruses in various environments, including the human intestinal microbiome. However, bacterial metagenomic studies frequently have the advantage. The adoption of advanced technologies like long-read sequencing has the potential to be transformative in refining viromics and metagenomics. Here, we examined the effectiveness of long-read and hybrid sequencing by comparing Illumina short-read and Oxford Nanopore Technology (ONT) long-read sequencing technologies and different assembly strategies on recovering viral genomes from human faecal samples. Our findings showed that if a single sequencing technology is to be chosen for virome analysis, Illumina is preferable due to its superior ability to recover fully resolved viral genomes and minimise erroneous genomes. While ONT assemblies were effective in recovering viral diversity, the challenges related to input requirements and the necessity for amplification made it less ideal as a standalone solution. However, using a combined, hybrid approach enabled a more authentic representation of viral diversity to be obtained within samples.

A microbial knowledge graph-based deep learning model for predicting candidate microbes for target hosts

Predicting interactions between microbes and hosts plays critical roles in microbiome population genetics and microbial ecology and evolution. How to systematically characterize the sophisticated mechanisms and signal interplay between microbes and hosts is a significant challenge for global health risks. Identifying microbe-host interactions (MHIs) can not only provide helpful insights into their fundamental regulatory mechanisms, but also facilitate the development of targeted therapies for microbial infections. In recent years, computational methods have become an appealing alternative due to the high risk and cost of wet-lab experiments. Therefore, in this study, we utilized rich microbial metagenomic information to construct a novel heterogeneous microbial network (HMN)-based model named KGVHI to predict candidate microbes for target hosts. Specifically, KGVHI first built a HMN by integrating human proteins, viruses and pathogenic bacteria with their biological attributes. Then KGVHI adopted a knowledge graph embedding strategy to capture the global topological structure information of the whole network. A natural language processing algorithm is used to extract the local biological attribute information from the nodes in HMN. Finally, we combined the local and global information and fed it into a blended deep neural network (DNN) for training and prediction. Compared to state-of-the-art methods, the comprehensive experimental results show that our model can obtain excellent results on the corresponding three MHI datasets. Furthermore, we also conducted two pathogenic bacteria case studies to further indicate that KGVHI has excellent predictive capabilities for potential MHI pairs.

Preclinical characterization and in silico safety assessment of three virulent bacteriophages targeting carbapenem-resistant uropathogenic Escherichia coli

Phage therapy has recently been revitalized in the West with many successful applications against multi-drug-resistant bacterial infections. However, the lack of geographically diverse bacteriophage (phage) genomes has constrained our understanding of phage diversity and its genetics underpinning host specificity, lytic capability, and phage-bacteria co-evolution. This study aims to locally isolate virulent phages against uropathogenic Escherichia coli (E. coli) and study its phenotypic and genomic features. Three obligately virulent Escherichia phages (øEc_Makalu_001, øEc_Makalu_002, and øEc_Makalu_003) that could infect uropathogenic E. coli were isolated and characterized. All three phages belonged to Krischvirus genus. One-step growth curve showed that the latent period of the phages ranged from 15 to 20 min, the outbreak period ~ 50 min, and the burst size ranged between 74 and 127 PFU/bacterium. Moreover, the phages could tolerate a pH range of 6 to 9 and a temperature range of 25-37 °C for up to 180 min without significant loss of phage viability. All phages showed a broad host spectrum and could lyse up to 30% of the 35 tested E. coli isolates. Genomes of all phages were approximately ~ 163 kb with a gene density of 1.73 gene/kbp and an average gene length of ~ 951 bp. The coding density in all phages was approximately 95%. Putative lysin, holin, endolysin, and spanin genes were found in the genomes of all three phages. All phages were strictly virulent with functional lysis modules and lacked any known virulence or toxin genes and antimicrobial resistance genes. Pre-clinical experimental and genomic analysis suggest these phages may be suitable candidates for therapeutic applications.

Collateral sensitivity increases the efficacy of a rationally designed bacteriophage combination to control Salmonella enterica

The ability of virulent bacteriophages to lyse bacteria influences bacterial evolution, fitness, and population structure. Knowledge of both host susceptibility and resistance factors is crucial for the successful application of bacteriophages as biological control agents in clinical therapy, food processing, and agriculture. In this study, we isolated 12 bacteriophages termed SPLA phage which infect the foodborne pathogen Salmonella enterica. To determine phage host range, a diverse collection of Enterobacteriaceae and Salmonella enterica was used and genes involved in infection by six SPLA phages were identified using Salmonella Typhimurium strain ST4/74. Candidate host receptors included lipopolysaccharide (LPS), cellulose, and BtuB. Lipopolysaccharide was identified as a susceptibility factor for phage SPLA1a and mutations in LPS biosynthesis genes spontaneously emerged during culture with S. Typhimurium. Conversely, LPS was a resistance factor for phage SPLA5b which suggested that emergence of LPS mutations in culture with SPLA1a represented collateral sensitivity to SPLA5b. We show that bacteria-phage co-culture with SPLA1a and SPLA5b was more successful in limiting the emergence of phage resistance compared to single phage co-culture. Identification of host susceptibility and resistance genes and understanding infection dynamics are critical steps in the rationale design of phage cocktails against specific bacterial pathogens.IMPORTANCEAs antibiotic resistance continues to emerge in bacterial pathogens, bacterial viruses (phage) represent a potential alternative or adjunct to antibiotics. One challenge for their implementation is the predisposition of bacteria to rapidly acquire resistance to phages. We describe a functional genomics approach to identify mechanisms of susceptibility and resistance for newly isolated phages that infect and lyse Salmonella enterica and use this information to identify phage combinations that exploit collateral sensitivity, thus increasing efficacy. Collateral sensitivity is a phenomenon where resistance to one class of antibiotics increases sensitivity to a second class of antibiotics. We report a functional genomics approach to rationally design a phage combination with a collateral sensitivity dynamic which resulted in increased efficacy. Considering such evolutionary trade-offs has the potential to manipulate the outcome of phage therapy in favor of resolving infection without selecting for escape mutants and is applicable to other virus-host interactions.

Xanthomonas Phage PBR31: Classifying the Unclassifiable

The ability of bacteriophages to destroy bacteria has made them the subject of extensive research. Interest in bacteriophages has recently increased due to the spread of drug-resistant bacteria, although genomic research has not kept pace with the growth of genomic data. Genomic analysis and, especially, the taxonomic description of bacteriophages are often difficult due to the peculiarities of the evolution of bacteriophages, which often includes the horizontal transfer of genes and genomic modules. The latter is particularly pronounced for temperate bacteriophages, which are capable of integration into the bacterial chromosome. Xanthomonas phage PBR31 is a temperate bacteriophage, which has been neither described nor classified previously, that infects the plant pathogen Xanthomonas campestris pv. campestris. Genomic analysis, including phylogenetic studies, indicated the separation of phage PBR31 from known classified bacteriophages, as well as its distant relationship with other temperate bacteriophages, including the Lederbervirus group. Bioinformatic analysis of proteins revealed distinctive features of PBR31, including the presence of a protein similar to the small subunit of D-family DNA polymerase and advanced lysis machinery. Taxonomic analysis showed the possibility of assigning phage PBR31 to a new taxon, although the complete taxonomic description of Xanthomonas phage PBR31 and other related bacteriophages is complicated by the complex evolutionary history of the formation of its genome. The general biological features of the PBR31 phage were analysed for the first time. Due to its presumably temperate lifestyle, there is doubt as to whether the PBR31 phage is appropriate for phage control purposes. Bioinformatics analysis, however, revealed the presence of cell wall-degrading enzymes that can be utilised for the treatment of bacterial infections.

Phage-inspired strategies to combat antibacterial resistance

Antimicrobial resistance (AMR) in clinically priority pathogensis now a major threat to public health worldwide. Phages are bacterial parasites that efficiently infect or kill specific strains and represent the most abundant biological entities on earth, showing great attraction as potential antibacterial therapeutics in combating AMR. This review provides a summary of phage-inspired strategies to combat AMR. We firstly cover the phage diversity, and then explain the biological principles of phage therapy that support the use of phages in the post-antimicrobial era. Furthermore, we state the versatility methods of phage therapy both from direct access as well as collateral access. Among the direct access approaches, we discuss the use of phage cocktail therapy, phage-encoded endolysins and the bioengineering for function improvement of used phages or endolysins. On the other hand, we introduce the collateral access, including the phages antimicrobial immunity combined therapy and phage-based novel antibacterial mimic molecules. Nowadays, more and more talented and enthusiastic scientist, doctors, pharmacists, media, authorities, and industry are promoting the progress of phage therapy, and proposed more phages-inspired strategy to make them more tractable to combat AMR and benefit more people, more animal and diverse environment in "one health" framework.

Exploring the gut DNA virome in fecal immunochemical test stool samples reveals associations with lifestyle in a large population-based study

Stool samples for fecal immunochemical tests (FIT) are collected in large numbers worldwide as part of colorectal cancer screening programs. Employing FIT samples from 1034 CRCbiome participants, recruited from a Norwegian colorectal cancer screening study, we identify, annotate and characterize more than 18000 DNA viruses, using shotgun metagenome sequencing. Only six percent of them are assigned to a known taxonomic family, with Microviridae being the most prevalent viral family. Linking individual profiles to comprehensive lifestyle and demographic data shows 17/25 of the variables to be associated with the gut virome. Physical activity, smoking, and dietary fiber consumption exhibit strong and consistent associations with both diversity and relative abundance of individual viruses, as well as with enrichment for auxiliary metabolic genes. We demonstrate the suitability of FIT samples for virome analysis, opening an opportunity for large-scale studies of this enigmatic part of the gut microbiome. The diverse viral populations and their connections to the individual lifestyle uncovered herein paves the way for further exploration of the role of the gut virome in health and disease.

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

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

The Citizen Phage Library: Rapid Isolation of Phages for the Treatment of Antibiotic Resistant Infections in the UK

Antimicrobial resistance poses one of the greatest threats to global health and there is an urgent need for new therapeutic options. Phages are viruses that infect and kill bacteria and phage therapy could provide a valuable tool for the treatment of multidrug-resistant infections. In this study, water samples collected by citizen scientists as part of the Citizen Phage Library (CPL) project, and wastewater samples from the Environment Agency yielded phages with activity against clinical strains Klebsiella pneumoniae BPRG1484 and Enterobacter cloacae BPRG1482. A total of 169 and 163 phages were found for K. pneumoniae and E. cloacae, respectively, within four days of receiving the strains. A third strain (Escherichia coli BPRG1486) demonstrated cross-reactivity with 42 E. coli phages already held in the CPL collection. Seed lots were prepared for four K. pneumoniae phages and a cocktail combining these phages was found to reduce melanisation in a Galleria mellonella infection model. The resources and protocols utilised by the Citizen Phage Library enabled the rapid isolation and characterisation of phages targeted against multiple strains. In the future, within a clearly defined regulatory framework, phage therapy could be made available on a named-patient basis within the UK.

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.

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.

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.

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.

Phylogenetic diversity and functional potential of large and cell-associated viruses in the Bay of Bengal

IMPORTANCE

The BoB, the world's largest bay, is of significant economic importance to surrounding countries, particularly Bangladesh, which heavily relies on its coastal resources. Concurrently, the BoB holds substantial ecological relevance due to the region's high vulnerability to climate change-induced impacts. Yet, our understanding of the BoB's microbiome in relation to marine food web and biogeochemical cycling remains limited. Particularly, there are little or no data on the viral diversity and host association in the BoB. We examined the viral community in two distinct BoB coastal regions to reveal a multitude of viral species interacting with a wide range of microbial hosts, some of which play key roles in coastal biogeochemical cycling or potential pathogens. Furthermore, we demonstrate that the BoB coast harbors a diverse community of large and giant viruses, underscoring the importance of investigating understudied environments to discover novel viral lineages with complex metabolic capacities.

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.

Genome sequence and characterization of a novel Pseudomonas putida phage, MiCath

Pseudomonads are ubiquitous bacteria with importance in medicine, soil, agriculture, and biomanufacturing. We report a novel Pseudomonas putida phage, MiCath, which is the first known phage infecting P. putida S12, a strain increasingly used as a synthetic biology chassis. MiCath was isolated from garden soil under a tomato plant using P. putida S12 as a host and was also found to infect four other P. putida strains. MiCath has a ~ 61 kbp double-stranded DNA genome which encodes 97 predicted open reading frames (ORFs); functions could only be predicted for 48 ORFs using comparative genomics. Functions include structural phage proteins, other common phage proteins (e.g., terminase), a queuosine gene cassette, a cas4 exonuclease, and an endosialidase. Restriction digestion analysis suggests the queuosine gene cassette encodes a pathway capable of modification of guanine residues. When compared to other phage genomes, MiCath shares at most 74% nucleotide identity over 2% of the genome with any sequenced phage. Overall, MiCath is a novel phage with no close relatives, encoding many unique gene products.

Investigating the Human Intestinal DNA Virome and Predicting Disease-Associated Virus-Host Interactions in Severe Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS)

Understanding how the human virome, and which of its constituents, contributes to health or disease states is reliant on obtaining comprehensive virome profiles. By combining DNA viromes from isolated virus-like particles (VLPs) and whole metagenomes from the same faecal sample of a small cohort of healthy individuals and patients with severe myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS), we have obtained a more inclusive profile of the human intestinal DNA virome. Key features are the identification of a core virome comprising tailed phages of the class Caudoviricetes, and a greater diversity of DNA viruses including extracellular phages and integrated prophages. Using an in silico approach, we predicted interactions between members of the Anaerotruncus genus and unique viruses present in ME/CFS microbiomes. This study therefore provides a framework and rationale for studies of larger cohorts of patients to further investigate disease-associated interactions between the intestinal virome and the bacteriome.

Genomic and proteomic characterization of vB_SauM-UFV_DC4, a novel Staphylococcus jumbo phage

Staphylococcus aureus is one of the most relevant mastitis pathogens in dairy cattle, and the acquisition of antimicrobial resistance genes presents a significant health issue in both veterinary and human fields. Among the different strategies to tackle S. aureus infection in livestock, bacteriophages have been thoroughly investigated in the last decades; however, few specimens of the so-called jumbo phages capable of infecting S. aureus have been described. Herein, we report the biological, genomic, and structural proteomic features of the jumbo phage vB_SauM-UFV_DC4 (DC4). DC4 exhibited a remarkable killing activity against S. aureus isolated from the veterinary environment and stability at alkaline conditions (pH 4 to 12). The complete genome of DC4 is 263,185 bp (GC content: 25%), encodes 263 predicted CDSs (80% without an assigned function), 1 tRNA (Phe-tRNA), multisubunit RNA polymerase, and an RNA-dependent DNA polymerase. Moreover, comparative analysis revealed that DC4 can be considered a new viral species belonging to a new genus DC4 and showed a similar set of lytic proteins and depolymerase activity with closely related jumbo phages. The characterization of a new S. aureus jumbo phage increases our understanding of the diversity of this group and provides insights into the biotechnological potential of these viruses. KEY POINTS: • vB_SauM-UFV_DC4 is a new viral species belonging to a new genus within the class Caudoviricetes. • vB_SauM-UFV_DC4 carries a set of RNA polymerase subunits and an RNA-directed DNA polymerase. • vB_SauM-UFV_DC4 and closely related jumbo phages showed a similar set of lytic proteins.

Narrow host range phages infect essential bacteria for water purification reactions in groundwater-fed rapid sand filters

Biofiltration is used worldwide to provide safe potable water due to its low energy demand and excellent treatment performance. For instance, in Denmark, over 95% of drinking water is supplied through groundwater-fed rapid sand filters (RSF). Bacteriophages, viruses that infect bacteria, have been shown to shape the taxonomic and functional composition of microbial communities across a range of natural and engineering systems. However, phages in the biofiltration systems are rarely studied, despite the central role microbes play in water purification. To probe this, metagenomic data from surface water, groundwater and mixed source water biofiltration units (n = 26 from China, Europe and USA) for drinking water production were analysed to characterize prokaryotic viruses and to identify their potential microbial hosts. The source water type and geographical location are found to exert influence on the composition of the phageome in biofilters. Although the viral abundance (71,676 ± 17,841 RPKM) in biofilters is only 14.4% and 17.0% lower than those of the nutrient-rich wastewater treatment plants and fresh surface waters, the richness (1,441 ± 1,046) and diversity (Inverse Simpson: 91 ± 61) in biofiltration units are significantly less by a factor of 2-5 and 3-4, respectively. In depth analysis of data from 24 groundwater-fed RSFs in Denmark revealed a core phageome shared by most RSFs, which was consistently linked to dominant microbial hosts involved in key biological reactions for water purification. Finally, the high number of specific links detected between phages and bacterial species and the large proportion of lytic phages (77%) led to the conjecture that phages regulate bacterial populations through predation, preventing the proliferation of dominant species and contributing to the established functional redundancy among the dominant microbial groups. In conclusion, bacteriophages are likely to play a significant role in water treatment within biofilters, particularly through interactions with key bacterial species.

Marine viruses disperse bidirectionally along the natural water cycle

Marine viruses in seawater have frequently been studied, yet their dispersal from neuston ecosystems at the air-sea interface towards the atmosphere remains a knowledge gap. Here, we show that 6.2% of the studied virus population were shared between air-sea interface ecosystems and rainwater. Virus enrichment in the 1-mm thin surface microlayer and sea foams happened selectively, and variant analysis proved virus transfer to aerosols collected at ~2 m height above sea level and rain. Viruses detected in rain and these aerosols showed a significantly higher percent G/C base content compared to marine viruses. CRISPR spacer matches of marine prokaryotes to foreign viruses from rainwater prove regular virus-host encounters at the air-sea interface. Our findings on aerosolization, adaptations, and dispersal support transmission of viruses along the natural water cycle.

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.

Bacteriophage communities are a reservoir of unexplored microbial diversity in neonatal health and disease

Acquisition and development of the gut microbiome are vital for immune education in neonates, especially those born preterm. As such, microbial communities have been extensively studied in the context of postnatal health and disease. Bacterial communities have been the focus of research in this area due to the relative ease of targeted bacterial sequencing and the availability of databases to align and validate sequencing data. Recent increases in high-throughput metagenomic sequencing accessibility have facilitated research to investigate bacteriophages within the context of neonatal gut microbial communities. Focusing on unexplored viral diversity, has identified novel bacteriophage species and previously uncharacterised viral diversity. In doing so, studies have highlighted links between bacteriophages and bacterial community structure in the context of health and disease. However, much remains unknown about the complex relationships between bacteriophages, the bacteria they infect and their human host. With a particular focus on preterm infants, this review highlights opportunities to explore the influence of bacteriophages on developing microbial communities and the tripartite relationships between bacteriophages, bacteria and the neonatal human host. We suggest a focus on expanding collections of isolated bacteriophages that will further our understanding of the growing numbers of bacteriophages identified in metagenomes.

Novel pelagiphage isolate Polarivirus skadi is a polar specialist that dominates SAR11-associated bacteriophage communities at high latitudes

The SAR11 clade are the most abundant members of surface marine bacterioplankton and a critical component of global biogeochemical cycles. Similarly, pelagiphages that infect SAR11 are ubiquitous and highly abundant in the oceans. Pelagiphages are predicted to shape SAR11 community structures and increase carbon turnover throughout the oceans. Yet, ecological drivers of host and niche specificity of pelagiphage populations are poorly understood. Here we report the global distribution of a novel pelagiphage called "Polarivirus skadi", which is the sole representative of a novel genus. P. skadi was isolated from the Western English Channel using a cold-water ecotype of SAR11 as bait. P. skadi is closely related to the globally dominant pelagiphage HTVC010P. Along with other HTVC010P-type viruses, P. skadi belongs to a distinct viral family within the order Caudovirales, for which we propose the name Ubiqueviridae. Metagenomic read recruitment identified P. skadi as one of the most abundant pelagiphages on Earth. P. skadi is a polar specialist, replacing HTVC010P at high latitudes. Experimental evaluation of P. skadi host range against cold- and warm-water SAR11 ecotypes supported cold-water specialism. Relative abundance of P. skadi in marine metagenomes correlated negatively with temperature, and positively with nutrients, available oxygen, and chlorophyll concentrations. In contrast, relative abundance of HTVC010P correlated negatively with oxygen and positively with salinity, with no significant correlation to temperature. The majority of other pelagiphages were scarce in most marine provinces, with a few representatives constrained to discrete ecological niches. Our results suggest that pelagiphage populations persist within a global viral seed bank, with environmental parameters and host availability selecting for a few ecotypes that dominate ocean viromes.

Microviruses: A World Beyond phiX174

Two decades of metagenomic analyses have revealed that in many environments, small (∼5 kb), single-stranded DNA phages of the family Microviridae dominate the virome. Although the emblematic microvirus phiX174 is ubiquitous in the laboratory, most other microviruses, particularly those of the gokushovirus and amoyvirus lineages, have proven to be much more elusive. This puzzling lack of representative isolates has hindered insights into microviral biology. Furthermore, the idiosyncratic size and nature of their genomes have resulted in considerable misjudgments of their actual abundance in nature. Fortunately, recent successes in microvirus isolation and improved metagenomic methodologies can now provide us with more accurate appraisals of their abundance, their hosts, and their interactions. The emerging picture is that phiX174 and its relatives are rather rare and atypical microviruses, and that a tremendous diversity of other microviruses is ready for exploration.

Diverse Durham collection phages demonstrate complex BREX defense responses

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

PTBGRP: predicting phage-bacteria interactions with graph representation learning on microbial heterogeneous information network

Identifying the potential bacteriophages (phage) candidate to treat bacterial infections plays an essential role in the research of human pathogens. Computational approaches are recognized as a valid way to predict bacteria and target phages. However, most of the current methods only utilize lower-order biological information without considering the higher-order connectivity patterns, which helps to improve the predictive accuracy. Therefore, we developed a novel microbial heterogeneous interaction network (MHIN)-based model called PTBGRP to predict new phages for bacterial hosts. Specifically, PTBGRP first constructs an MHIN by integrating phage-bacteria interaction (PBI) and six bacteria-bacteria interaction networks with their biological attributes. Then, different representation learning methods are deployed to extract higher-level biological features and lower-level topological features from MHIN. Finally, PTBGRP employs a deep neural network as the classifier to predict unknown PBI pairs based on the fused biological information. Experiment results demonstrated that PTBGRP achieves the best performance on the corresponding ESKAPE pathogens and PBI dataset when compared with state-of-art methods. In addition, case studies of Klebsiella pneumoniae and Staphylococcus aureus further indicate that the consideration of rich heterogeneous information enables PTBGRP to accurately predict PBI from a more comprehensive perspective. The webserver of the PTBGRP predictor is freely available at http://120.77.11.78/PTBGRP/.

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.

Ecogenomics and cultivation reveal distinctive viral-bacterial communities in the surface microlayer of a Baltic Sea slick

Visible surface films, termed slicks, can extensively cover freshwater and marine ecosystems, with coastal regions being particularly susceptible to their presence. The sea-surface microlayer (SML), the upper 1-mm at the air-water interface in slicks (herein slick SML) harbors a distinctive bacterial community, but generally little is known about SML viruses. Using flow cytometry, metagenomics, and cultivation, we characterized viruses and bacteria in a brackish slick SML in comparison to non-slick SML as well as seawater below slick and non-slick areas (subsurface water = SSW). Size-fractionated filtration of all samples distinguished viral attachment to hosts and particles. The slick SML contained higher abundances of virus-like particles, prokaryotic cells, and dissolved organic carbon compared to non-slick SML and SSW. The community of 428 viral operational taxonomic units (vOTUs), 426 predicted as lytic, distinctly differed across all size fractions in the slick SML compared to non-slick SML and SSW. Specific metabolic profiles of bacterial metagenome-assembled genomes and isolates in the slick SML included a prevalence of genes encoding motility and carbohydrate-active enzymes (CAZymes). Several vOTUs were enriched in slick SML, and many virus variants were associated with particles. Nine vOTUs were only found in slick SML, six of them being targeted by slick SML-specific clustered-regularly interspaced short palindromic repeats (CRISPR) spacers likely originating from Gammaproteobacteria. Moreover, isolation of three previously unknown lytic phages for Alishewanella sp. and Pseudoalteromonas tunicata, abundant and actively replicating slick SML bacteria, suggests that viral activity in slicks contributes to biogeochemical cycling in coastal ecosystems.

A novel Queuovirinae lineage of Pseudomonas aeruginosa phages encode dPreQ0 DNA modifications with a single GA motif that provide restriction and CRISPR Cas9 protection in vitro

Deazaguanine DNA modifications are widespread in phages, particularly in those with pathogenic hosts. Pseudomonas phage iggy substitutes ∼16.5% of its genomic 2'-deoxyguanosine (G) with dPreQ0, and the iggy deazaguanine transglycosylase (DpdA) is unique in having a strict GA target motif, not observed previously. The iggy PreQ0 modification is shown to provide protection against both restriction endonucleases and Cas9 (when present in PAM), thus expanding our understanding of the deazaguanine modification system, its potential, and diversity. Phage iggy represents a new genus of Pseudomonas phages within the Queuovirinae subfamily; which have very little in common with other published phage genomes in terms of nucleotide similarity (<10%) and common proteins (<2%). Interestingly, shared similarity is concentrated in dpdA and preQ0 biosynthesis genes. TEM imaging confirmed a siphovirus morphology with a prolate icosahedral head and a non-contractile flexible tail with one long central tail spike. The observed protective effect of the deazaguanine modification on the iggy DNA may contribute to its broad within-species host range. Phage iggy was isolated on Pseudomonas aeruginosa PAO1, but also infects PDO300, PAK, PA14, as well as 10 of 27 tested environmental isolates and 13 of 20 tested clinical isolates of P. aeruginosa from patients with cystic fibrosis.

A predicted CRISPR-mediated symbiosis between uncultivated archaea

CRISPR-Cas systems defend prokaryotic cells from invasive DNA of viruses, plasmids and other mobile genetic elements. Here, we show using metagenomics, metatranscriptomics and single-cell genomics that CRISPR systems of widespread, uncultivated archaea can also target chromosomal DNA of archaeal episymbionts of the DPANN superphylum. Using meta-omics datasets from Crystal Geyser and Horonobe Underground Research Laboratory, we find that CRISPR spacers of the hosts Candidatus Altiarchaeum crystalense and Ca. A. horonobense, respectively, match putative essential genes in their episymbionts' genomes of the genus Ca. Huberiarchaeum and that some of these spacers are expressed in situ. Metabolic interaction modelling also reveals complementation between host-episymbiont systems, on the basis of which we propose that episymbionts are either parasitic or mutualistic depending on the genotype of the host. By expanding our analysis to 7,012 archaeal genomes, we suggest that CRISPR-Cas targeting of genomes associated with symbiotic archaea evolved independently in various archaeal lineages.

Functional biology and biotechnology of thermophilic viruses

Viruses have developed sophisticated biochemical and genetic mechanisms to manipulate and exploit their hosts. Enzymes derived from viruses have been essential research tools since the first days of molecular biology. However, most viral enzymes that have been commercialized are derived from a small number of cultivated viruses, which is remarkable considering the extraordinary diversity and abundance of viruses revealed by metagenomic analysis. Given the explosion of new enzymatic reagents derived from thermophilic prokaryotes over the past 40 years, those obtained from thermophilic viruses should be equally potent tools. This review discusses the still-limited state of the art regarding the functional biology and biotechnology of thermophilic viruses with a focus on DNA polymerases, ligases, endolysins, and coat proteins. Functional analysis of DNA polymerases and primase-polymerases from phages infecting Thermus, Aquificaceae, and Nitratiruptor has revealed new clades of enzymes with strong proofreading and reverse transcriptase capabilities. Thermophilic RNA ligase 1 homologs have been characterized from Rhodothermus and Thermus phages, with both commercialized for circularization of single-stranded templates. Endolysins from phages infecting Thermus, Meiothermus, and Geobacillus have shown high stability and unusually broad lytic activity against Gram-negative and Gram-positive bacteria, making them targets for commercialization as antimicrobials. Coat proteins from thermophilic viruses infecting Sulfolobales and Thermus strains have been characterized, with diverse potential applications as molecular shuttles. To gauge the scale of untapped resources for these proteins, we also document over 20,000 genes encoded by uncultivated viral genomes from high-temperature environments that encode DNA polymerase, ligase, endolysin, or coat protein domains.

Phage diversity in cell-free DNA identifies bacterial pathogens in human sepsis cases

Bacteriophages, viruses that infect bacteria, have great specificity for their bacterial hosts at the strain and species level. However, the relationship between the phageome and associated bacterial population dynamics is unclear. Here we generated a computational pipeline to identify sequences associated with bacteriophages and their bacterial hosts in cell-free DNA from plasma samples. Analysis of two independent cohorts, including a Stanford Cohort of 61 septic patients and 10 controls and the SeqStudy cohort of 224 septic patients and 167 controls, reveals a circulating phageome in the plasma of all sampled individuals. Moreover, infection is associated with overrepresentation of pathogen-specific phages, allowing for identification of bacterial pathogens. We find that information on phage diversity enables identification of the bacteria that produced these phages, including pathovariant strains of Escherichia coli. Phage sequences can likewise be used to distinguish between closely related bacterial species such as Staphylococcus aureus, a frequent pathogen, and coagulase-negative Staphylococcus, a frequent contaminant. Phage cell-free DNA may have utility in studying bacterial infections.