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

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.

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.

Prophages and plasmids can display opposite trends in the types of accessory genes they carry

Mobile genetic elements (MGEs), such as phages and plasmids, often possess accessory genes encoding bacterial functions, facilitating bacterial evolution. Are there rules governing the arsenal of accessory genes MGEs carry? If such rules exist, they might be reflected in the types of accessory genes different MGEs carry. To test this hypothesis, we compare prophages and plasmids with respect to the frequencies at which they carry antibiotic resistance genes (ARGs) and virulence factor genes (VFGs) in the genomes of 21 pathogenic bacterial species using public databases. Our results indicate that prophages tend to carry VFGs more frequently than ARGs in three species, whereas plasmids tend to carry ARGs more frequently than VFGs in nine species, relative to genomic backgrounds. In Escherichia coli, where this prophage-plasmid disparity is detected, prophage-borne VFGs encode a much narrower range of functions than do plasmid-borne VFGs, typically involved in damaging host cells or modulating host immunity. In the species where the above disparity is not detected, ARGs and VFGs are barely found in prophages and plasmids. These results indicate that MGEs can differentiate in the types of accessory genes they carry depending on their infection strategies, suggesting a rule governing horizontal gene transfer mediated by MGEs.

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.

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.

Going through phages: a computational approach to revealing the role of prophage in Staphylococcus aureus

Prophages have important roles in virulence, antibiotic resistance, and genome evolution in Staphylococcus aureus . Rapid growth in the number of sequenced S. aureus genomes allows for an investigation of prophage sequences at an unprecedented scale. We developed a novel computational pipeline for phage discovery and annotation. We combined PhiSpy, a phage discovery tool, with VGAS and PROKKA, genome annotation tools to detect and analyse prophage sequences in nearly 10 011  S . aureus genomes, discovering thousands of putative prophage sequences with genes encoding virulence factors and antibiotic resistance. To our knowledge, this is the first large-scale application of PhiSpy on a large-scale set of genomes (10 011  S . aureus ). Determining the presence of virulence and resistance encoding genes in prophage has implications for the potential transfer of these genes/functions to other bacteria via transduction and thus can provide insight into the evolution and spread of these genes/functions between bacterial strains. While the phage we have identified may be known, these phages were not necessarily known or characterized in S. aureus and the clustering and comparison we did for phage based on their gene content is novel. Moreover, the reporting of these genes with the S. aureus genomes is novel.

Burkholderia contaminans Bacteriophage CSP3 Requires O-Antigen Polysaccharides for Infection

The Burkholderia cepacia complex is a group of opportunistic pathogens that cause both severe acute and chronic respiratory infections. Due to their large genomes containing multiple intrinsic and acquired antimicrobial resistance mechanisms, treatment is often difficult and prolonged. One alternative to traditional antibiotics for treatment of bacterial infections is bacteriophages. Therefore, the characterization of bacteriophages infective for the Burkholderia cepacia complex is critical to determine their suitability for any future use. Here, we describe the isolation and characterization of novel phage, CSP3, infective against a clinical isolate of Burkholderia contaminans. CSP3 is a new member of the Lessievirus genus that targets various Burkholderia cepacia complex organisms. Single nucleotide polymorphism (SNP) analysis of CSP3-resistant B. contaminans showed that mutations to the O-antigen ligase gene, waaL, consequently inhibited CSP3 infection. This mutant phenotype is predicted to result in the loss of cell surface O-antigen, contrary to a related phage that requires the inner core of the lipopolysaccharide for infection. Additionally, liquid infection assays showed that CSP3 provides suppression of B. contaminans growth for up to 14 h. Despite the inclusion of genes that are typical of the phage lysogenic life cycle, we saw no evidence of CSP3's ability to lysogenize. Continuation of phage isolation and characterization is crucial in developing large and diverse phage banks for global usage in cases of antibiotic-resistant bacterial infections. IMPORTANCE Amid the global antibiotic resistance crisis, novel antimicrobials are needed to treat problematic bacterial infections, including those from the Burkholderia cepacia complex. One such alternative is the use of bacteriophages; however, a lot is still unknown about their biology. Bacteriophage characterization studies are of high importance for building phage banks, as future work in developing treatments such as phage cocktails should require well-characterized phages. Here, we report the isolation and characterization of a novel Burkholderia contaminans phage that requires the O-antigen for infection, a distinct phenotype seen among other related phages. Our findings presented in this article expand on the ever-evolving phage biology field, uncovering unique phage-host relationships and mechanisms of infection.

Characterization of prophages in bacterial genomes from the honey bee (Apis mellifera) gut microbiome

The gut of the European honey bee (Apis mellifera) possesses a relatively simple bacterial community, but little is known about its community of prophages (temperate bacteriophages integrated into the bacterial genome). Although prophages may eventually begin replicating and kill their bacterial hosts, they can also sometimes be beneficial for their hosts by conferring protection from other phage infections or encoding genes in metabolic pathways and for toxins. In this study, we explored prophages in 17 species of core bacteria in the honey bee gut and two honey bee pathogens. Out of the 181 genomes examined, 431 putative prophage regions were predicted. Among core gut bacteria, the number of prophages per genome ranged from zero to seven and prophage composition (the compositional percentage of each bacterial genome attributable to prophages) ranged from 0 to 7%. Snodgrassella alvi and Gilliamella apicola had the highest median prophages per genome (3.0 ± 1.46; 3.0 ± 1.59), as well as the highest prophage composition (2.58% ± 1.4; 3.0% ± 1.59). The pathogen Paenibacillus larvae had a higher median number of prophages (8.0 ± 5.33) and prophage composition (6.40% ± 3.08) than the pathogen Melissococcus plutonius or any of the core bacteria. Prophage populations were highly specific to their bacterial host species, suggesting most prophages were acquired recently relative to the divergence of these bacterial groups. Furthermore, functional annotation of the predicted genes encoded within the prophage regions indicates that some prophages in the honey bee gut encode additional benefits to their bacterial hosts, such as genes in carbohydrate metabolism. Collectively, this survey suggests that prophages within the honey bee gut may contribute to the maintenance and stability of the honey bee gut microbiome and potentially modulate specific members of the bacterial community, particularly S. alvi and G. apicola.

Isolation and in vitro characterization of novel S. epidermidis phages for therapeutic applications

S. epidermidis is an important opportunistic pathogen causing chronic prosthetic joint infections associated with biofilm growth. Increased tolerance to antibiotic therapy often requires prolonged treatment or revision surgery. Phage therapy is currently used as compassionate use therapy and continues to be evaluated for its viability as adjunctive therapy to antibiotic treatment or as an alternative treatment for infections caused by S. epidermidis to prevent relapses. In the present study, we report the isolation and in vitro characterization of three novel lytic S. epidermidis phages. Their genome content analysis indicated the absence of antibiotic resistance genes and virulence factors. Detailed investigation of the phage preparation indicated the absence of any prophage-related contamination and demonstrated the importance of selecting appropriate hosts for phage development from the outset. The isolated phages infect a high proportion of clinically relevant S. epidermidis strains and several other coagulase-negative species growing both in planktonic culture and as a biofilm. Clinical strains differing in their biofilm phenotype and antibiotic resistance profile were selected to further identify possible mechanisms behind increased tolerance to isolated phages.

Highly diverse and unknown viruses may enhance Antarctic endoliths' adaptability

BACKGROUND

Rock-dwelling microorganisms are key players in ecosystem functioning of Antarctic ice free-areas. Yet, little is known about their diversity and ecology, and further still, viruses in these communities have been largely unexplored despite important roles related to host metabolism and nutrient cycling. To begin to address this, we present a large-scale viral catalog from Antarctic rock microbial communities.

RESULTS

We performed metagenomic analyses on rocks from across Antarctica representing a broad range of environmental and spatial conditions, and which resulted in a predicted viral catalog comprising > 75,000 viral operational taxonomic units (vOTUS). We found largely undescribed, highly diverse and spatially structured virus communities which had predicted auxiliary metabolic genes (AMGs) with functions indicating that they may be potentially influencing bacterial adaptation and biogeochemistry.

CONCLUSION

This catalog lays the foundation for expanding knowledge of virosphere diversity, function, spatial ecology, and dynamics in extreme environments. This work serves as a step towards exploring adaptability of microbial communities in the face of a changing climate. Video Abstract.

iPHoP: An integrated machine learning framework to maximize host prediction for metagenome-derived viruses of archaea and bacteria

The extraordinary diversity of viruses infecting bacteria and archaea is now primarily studied through metagenomics. While metagenomes enable high-throughput exploration of the viral sequence space, metagenome-derived sequences lack key information compared to isolated viruses, in particular host association. Different computational approaches are available to predict the host(s) of uncultivated viruses based on their genome sequences, but thus far individual approaches are limited either in precision or in recall, i.e., for a number of viruses they yield erroneous predictions or no prediction at all. Here, we describe iPHoP, a two-step framework that integrates multiple methods to reliably predict host taxonomy at the genus rank for a broad range of viruses infecting bacteria and archaea, while retaining a low false discovery rate. Based on a large dataset of metagenome-derived virus genomes from the IMG/VR database, we illustrate how iPHoP can provide extensive host prediction and guide further characterization of uncultivated viruses.

Dynamic Development of Viral and Bacterial Diversity during Grass Silage Preservation

Ensilaging is one of the most common feed preservation processes using lactic acid bacteria to stabilize feed and save feed quality. The silage bacterial community is well known but the role of the virome and its relationship with the bacterial community is scarce. In the present study, metagenomics and amplicon sequencing were used to describe the composition of the bacterial and viral community during a 40-day grass silage preservation. During the first two days, we observed a rapid decrease in the pH and a shift in the bacterial and viral composition. The diversity of the dominant virus operational taxonomic units (vOTUs) decreased throughout the preservation. The changes in the bacterial community resembled the predicted putative host of the recovered vOTUs during each sampling time. Only 10% of the total recovered vOTUs clustered with a reference genome. Different antiviral defense mechanisms were found across the recovered metagenome-assembled genomes (MAGs); however, only a history of bacteriophage infection with Lentilactobacillus and Levilactobacillus was observed. In addition, vOTUs harbored potential auxiliary metabolic genes related to carbohydrate metabolism, organic nitrogen, stress tolerance, and transport. Our data suggest that vOTUs are enriched during grass silage preservation, and they could have a role in the establishment of the bacterial community.

Expanding known viral diversity in the healthy infant gut

The gut microbiome is shaped through infancy and impacts the maturation of the immune system, thus protecting against chronic disease later in life. Phages, or viruses that infect bacteria, modulate bacterial growth by lysis and lysogeny, with the latter being especially prominent in the infant gut. Viral metagenomes (viromes) are difficult to analyse because they span uncharted viral diversity, lacking marker genes and standardized detection methods. Here we systematically resolved the viral diversity in faecal viromes from 647 1-year-olds belonging to Copenhagen Prospective Studies on Asthma in Childhood 2010, an unselected Danish cohort of healthy mother-child pairs. By assembly and curation we uncovered 10,000 viral species from 248 virus family-level clades (VFCs). Most (232 VFCs) were previously unknown, belonging to the Caudoviricetes viral class. Hosts were determined for 79% of phage using clustered regularly interspaced short palindromic repeat spacers within bacterial metagenomes from the same children. Typical Bacteroides-infecting crAssphages were outnumbered by undescribed phage families infecting Clostridiales and Bifidobacterium. Phage lifestyles were conserved at the viral family level, with 33 virulent and 118 temperate phage families. Virulent phages were more abundant, while temperate ones were more prevalent and diverse. Together, the viral families found in this study expand existing phage taxonomy and provide a resource aiding future infant gut virome research.

Lytic bacteriophage vB_KmiS-Kmi2C disrupts biofilms formed by members of the Klebsiella oxytoca complex, and represents a novel virus family and genus

AIMS

This study aimed to characterize the lytic phage vB_KmiS-Kmi2C, isolated from sewage water on a GES-positive strain of Klebsiella michiganensis.

METHODS AND RESULTS

Comparative phylogenetic and network-based analyses were used to characterize the genome of phage vB_KmiS-Kmi2C (circular genome of 42 234 bp predicted to encode 55 genes), demonstrating it shared little similarity with other known phages. The phage was lytic on clinical strains of K. oxytoca (n = 2) and K. michiganensis (n = 4), and was found to both prevent biofilm formation and disrupt established biofilms produced by these strains.

CONCLUSIONS

We have identified a phage capable of killing clinically relevant members of the K. oxytoca complex (KoC). The phage represents a novel virus family (proposed name Dilsviridae) and genus (proposed name Dilsvirus).

Genetic determinants of host tropism in Klebsiella phages

Bacteriophages play key roles in bacterial ecology and evolution and are potential antimicrobials. However, the determinants of phage-host specificity remain elusive. Here, we isolate 46 phages to challenge 138 representative clinical isolates of Klebsiella pneumoniae, a widespread opportunistic pathogen. Spot tests show a narrow host range for most phages, with <2% of 6,319 phage-host combinations tested yielding detectable interactions. Bacterial capsule diversity is the main factor restricting phage host range. Consequently, phage-encoded depolymerases are key determinants of host tropism, and depolymerase sequence types are associated with the ability to infect specific capsular types across phage families. However, all phages with a broader host range found do not encode canonical depolymerases, suggesting alternative modes of entry. These findings expand our knowledge of the complex interactions between bacteria and their viruses and point out the feasibility of predicting the first steps of phage infection using bacterial and phage genome sequences.

What Lies Beneath? Taking the Plunge into the Murky Waters of Phage Biology

The sequence revolution revealed that bacteria-infecting viruses, known as phages, are Earth's most abundant biological entities. Phages have far-reaching impacts on the form and function of microbial communities and play a fundamental role in ecological processes. However, even well into the sequencing revolution, we have only just begun to explore the murky waters around the phage biology iceberg. Many viral reads cannot be assigned to a culturable isolate, and reference databases are biased toward more easily collectible samples, which likely distorts our conclusions. This minireview points out alternatives to mapping reads to reference databases and highlights innovative bioinformatic and experimental approaches that can help us overcome some of the challenges in phage research and better decipher the impact of phages on microbial communities. Moving beyond the identification of novel phages, we highlight phage metabolomics as an important influencer of bacterial host cell physiology and hope to inspire the reader to consider the effects of phages on host metabolism and ecosystems at large. We encourage researchers to report unassigned/unknown sequencing reads and contigs and to continue developing alternative methods to investigate phages within sequence data.

Expansion of Kuravirus-like Phage Sequences within the Past Decade, including Escherichia Phage YF01 from Japan, Prompt the Creation of Three New Genera

Bacteriophages, viruses that infect bacteria, are currently receiving significant attention amid an ever-growing global antibiotic resistance crisis. In tandem, a surge in the availability and affordability of next-generation and third-generation sequencing technologies has driven the deposition of a wealth of phage sequence data. Here, we have isolated a novel Escherichia phage, YF01, from a municipal wastewater treatment plant in Yokohama, Japan. We demonstrate that the YF01 phage shares a high similarity to a collection of thirty-five Escherichia and Shigella phages found in public databases, six of which have been previously classified into the Kuravirus genus by the International Committee on Taxonomy of Viruses (ICTV). Using modern phylogenetic approaches, we demonstrate that an expansion and reshaping of the current six-membered Kuravirus genus is required to accommodate all thirty-six member phages. Ultimately, we propose the creation of three additional genera, Vellorevirus, Jinjuvirus, and Yesanvirus, which will allow a more organized approach to the addition of future Kuravirus-like phages.

Identification of integrons and gene cassette-associated recombination sites in bacteriophage genomes

Bacteriophages are versatile mobile genetic elements that play key roles in driving the evolution of their bacterial hosts through horizontal gene transfer. Phages co-evolve with their bacterial hosts and have plastic genomes with extensive mosaicism. In this study, we present bioinformatic and experimental evidence that temperate and virulent (lytic) phages carry integrons, including integron-integrase genes, attC/attI recombination sites and gene cassettes. Integrons are normally found in Bacteria, where they capture, express and re-arrange mobile gene cassettes via integron-integrase activity. We demonstrate experimentally that a panel of attC sites carried in virulent phage can be recognized by the bacterial class 1 integron-integrase (IntI1) and then integrated into the paradigmatic attI1 recombination site using an attC x attI recombination assay. With an increasing number of phage genomes projected to become available, more phage-associated integrons and their components will likely be identified in the future. The discovery of integron components in bacteriophages establishes a new route for lateral transfer of these elements and their cargo genes between bacterial host cells.

Three Phages One Host: Isolation and Characterization of Pantoea agglomerans Phages from a Grasshopper Specimen

The bacterial genus Pantoea comprises species found in a variety of different environmental sources. Pantoea spp. are often recovered from plant material and are capable of both benefitting the plants and acting like phytopathogens. Some species of Pantoea (including P. agglomerans) are considered opportunistic human pathogens capable of causing various infections in immunocompromised subjects. In this study, a strain of P. agglomerans (identified by 16S rRNA gene sequencing) was isolated from a dead specimen of an unidentified Latvian grasshopper species. The retrieved strain of P. agglomerans was then used as a host for the potential retrieval of phages from the same source material. After rounds of plaque purification and propagation, three high-titer lysates corresponding to putatively distinct phages were acquired. Transmission electron microscopy revealed that one of the phages was a myophage with an unusual morphology, while the two others were typical podophages. Whole-genome sequencing (WGS) was performed for each of these isolated phages. Genome de novo assembly and subsequent functional annotation confirmed that three different strictly lytic phages were isolated. Elaborate genomic characterization of the acquired phages was performed to elucidate their place within the so-far-uncovered phage diversity.

Significant Differences in Planktonic Virus Communities Between "Cellular Fraction" (0.22 ~ 3.0 µm) and "Viral Fraction" (< 0.22 μm) in the Ocean

Compared to free-living viruses (< 0.22 m) in the ocean, planktonic viruses in the "cellular fraction" (0.22 ~ 3.0 μm) are now far less well understood, and the differences between them remain largely unexplored. Here, we revealed that even in the same seawater samples, the "cellular fraction" comprised significantly distinct virus communities from the free virioplankton, with only 13.87% overlap in viral contigs at the species level. Compared to the viral genomes deposited in NCBI RefSeq database, 99% of the assembled viral genomes in the "cellular fraction" represented novel genera. Notably, the assembled (near-) complete viral genomes within the "cellular fraction" were significantly larger than that in the "viral fraction," and the "cellular fraction" contained three times more species of giant viruses or jumbo phages with genomes > 200 kb than the "viral fraction." The longest complete genomes of jumbo phage (~ 252 kb) and giant virus (~ 716 kb) were both detected only in the "cellular fraction." Moreover, a relatively higher proportion of proviruses were predicted within the "cellular fraction" than "viral fraction." Besides the substantial divergence in viral community structure, the different fractions also contained their unique viral auxiliary metabolic genes; e.g., those potentially participating in inorganic carbon fixation in deep sea were detected only in the "cellular-fraction" viromes. In addition, there was a considerable divergence in the community structure of both "cellular fraction" and "viral fraction" viromes between the surface and deep-sea habitats, suggesting that they might have similar environmental adaptation properties. The findings deepen our understanding of the complexity of viral community structure and function in the ocean.

Comprehensive Genome Analysis of Neisseria meningitidis from South America Reveals a Distinctive Pathogenicity-Related Prophage Repertoire

Neisseria meningitidis, a bacterium that colonizes in the human nasopharynx, occasionally causes invasive meningococcal disease leading to meningitis or septicemia. Different serogroups and lineages (clonal complexes) are related to the occurrence and epidemiology of N. meningitidis. Despite vaccines for most serogroups, N. meningitidis lineages causing unusual clinical manifestations and a higher fatality rate compared to other lineages have been reported in South America. The present study focused on exploring the diversity of N. meningitidis prophages from South America and their relationship with the epidemiological variables of these strains. We found a high diversity of prophages among the different clonal complexes. By comparing them with previously described N. meningitidis phages and prophages, we revealed groups of prophages sharing similar compositions, which could be useful for prophage comparison in N. meningitidis. Furthermore, we observed a high correlation between the prophage content and epidemiological features, e.g., pathogenicity or clonal complex. Additionally, a distinctive filamentous prophage named here as IMSAR-11 (Invasive Meningococci from South America Related to cc11) was identified. Interestingly, two versions of IMSAR-11, circular and chromosomally integrated, were found. Overall, this study reinforces the importance of the genomic characterization of circulating N. meningitidis lineages to generate new targets for lineage monitoring, diagnosis, or appropriateness of vaccine development. Further studies are necessary to understand the role of these prophages in the persistence, dispersal, and virulence of N. meningitidis in the world.

Genomic Profiling of Non-O157 Shiga Toxigenic Escherichia coli-Infecting Bacteriophages from South Africa

Background

Non-O157 Shiga toxigenic Escherichia coli (STEC) are one of the most important food and waterborne pathogens worldwide. Although bacteriophages (phages) have been used for the biocontrol of these pathogens, a comprehensive understanding of the genetic characteristics and lifestyle of potentially effective candidate phages is lacking.

Materials and Methods

In this study, 10 non-O157-infecting phages previously isolated from feedlot cattle and dairy farms in the North-West province of South Africa were sequenced, and their genomes were analyzed.

Results

Comparative genomics and proteomics revealed that the phages were closely related to other E. coli-infecting Tunaviruses, Seuratviruses, Carltongylesviruses, Tequatroviruses, and Mosigviruses from the National Center for Biotechnology Information GenBank database. Phages lacked integrases associated with a lysogenic cycle and genes associated with antibiotic resistance and Shiga toxins.

Conclusions

Comparative genomic analysis identified a diversity of unique non-O157-infecting phages, which could be used to mitigate the abundance of various non-O157 STEC serogroups without safety concerns.

Tyroviruses are a new group of temperate phages that infect Bacillus species in soil environments worldwide

BACKGROUND

Bacteriophages are widely considered to be highly abundant and genetically diverse, with their role in the evolution and virulence of many pathogens becoming increasingly clear. Less attention has been paid on phages preying on Bacillus, despite the potential for some of its members, such as Bacillus anthracis, to cause serious human disease.

RESULTS

We have isolated five phages infecting the causative agent of anthrax, Bacillus anthracis. Using modern phylogenetic approaches we place these five new Bacillus phages, as well as 21 similar phage genomes retrieved from publicly available databases and metagenomic datasets into the Tyrovirus group, a newly proposed group named so due to the conservation of three distinct tyrosine recombinases. Genomic analysis of these large phages (~ 160-170 kb) reveals their DNA packaging mechanism and genomic features contributing to virion morphogenesis, host cell lysis and phage DNA replication processes. Analysis of the three tyrosine recombinases suggest Tyroviruses undergo a prophage lifecycle that may involve both host integration and plasmid stages. Further we show that Tyroviruses rely on divergent invasion mechanisms, with a subset requiring host S-layer for infection.

CONCLUSIONS

Ultimately, we expand upon our understanding on the classification, phylogeny, and genomic organisation of a new and substantial phage group that prey on critically relevant Bacillus species. In an era characterised by a rapidly evolving landscape of phage genomics the deposition of future Tyroviruses will allow the further unravelling of the global spread and evolutionary history of these Bacillus phages.

Tracking the phage trends: A comprehensive review of applications in therapy and food production

In the present scenario, the challenge of emerging antimicrobial resistance is affecting human health globally. The increasing incidences of multidrug-resistant infections have become harder to treat, causing high morbidity, and mortality, and are posing extensive financial loss. Limited discovery of new antibiotic molecules has further complicated the situation and has forced researchers to think and explore alternatives to antibiotics. This has led to the resurgence of the bacteriophages as an effective alternative as they have a proven history in the Eastern world where lytic bacteriophages have been used since their first implementation over a century ago. To help researchers and clinicians towards strengthening bacteriophages as a more effective, safe, and economical therapeutic alternative, the present review provides an elaborate narrative about the important aspects of bacteriophages. It abridges the prerequisite essential requirements of phage therapy, the role of phage biobank, and the details of immune responses reported while using bacteriophages in the clinical trials/compassionate grounds by examining the up-to-date case reports and their effects on the human gut microbiome. This review also discusses the potential of bacteriophages as a biocontrol agent against food-borne diseases in the food industry and aquaculture, in addition to clinical therapy. It finishes with a discussion of the major challenges, as well as phage therapy and phage-mediated biocontrols future prospects.

Expanding the genomic encyclopedia of Actinobacteria with 824 isolate reference genomes

The phylum Actinobacteria includes important human pathogens like Mycobacterium tuberculosis and Corynebacterium diphtheriae and renowned producers of secondary metabolites of commercial interest, yet only a small part of its diversity is represented by sequenced genomes. Here, we present 824 actinobacterial isolate genomes in the context of a phylum-wide analysis of 6,700 genomes including public isolates and metagenome-assembled genomes (MAGs). We estimate that only 30%-50% of projected actinobacterial phylogenetic diversity possesses genomic representation via isolates and MAGs. A comparison of gene functions reveals novel determinants of host-microbe interaction as well as environment-specific adaptations such as potential antimicrobial peptides. We identify plasmids and prophages across isolates and uncover extensive prophage diversity structured mainly by host taxonomy. Analysis of >80,000 biosynthetic gene clusters reveals that horizontal gene transfer and gene loss shape secondary metabolite repertoire across taxa. Our observations illustrate the essential role of and need for high-quality isolate genome sequences.

Genome Sequences of Gordonia Phages GrootJr and NovumRegina

Two Gordonia bacteriophages, GrootJr and NovumRegina, were discovered, sequenced, and annotated. These phages were isolated from distinct soil and water samples, respectively, on Gordonia terrae 3612. These phages belong to the CR2 subcluster and are similar in genome size and gene content.

Characteristics and Comparative Genomic Analysis of a Novel Virus, VarioGold, the First Bacteriophage of Variovorax

Variovorax represents a widespread and ecologically significant genus of soil bacteria. Despite the ecological importance of these bacteria, our knowledge about the viruses infecting Variovorax spp. is quite poor. This study describes the isolation and characterization of the mitomycin-induced phage, named VarioGold. To the best of our knowledge, VarioGold represents the first characterized virus for this genus. Comparative genomic analyses suggested that VarioGold is distinct from currently known bacteriophages at both the nucleotide and protein levels; thus, it could be considered a new virus genus. In addition, another 37 prophages were distinguished in silico within the complete genomic sequences of Variovorax spp. that are available in public databases. The similarity networking analysis highlighted their general high diversity, which, despite clustering with previously described phages, shows their unique genetic load. Therefore, the novelty of Variovorax phages warrants the great enrichment of databases, which could, in turn, improve bioinformatic strategies for finding (pro)phages.

MetaPhage: an Automated Pipeline for Analyzing, Annotating, and Classifying Bacteriophages in Metagenomics Sequencing Data

Phages are the most abundant biological entities on the planet, and they play an important role in controlling density, diversity, and network interactions among bacterial communities through predation and gene transfer. To date, a variety of bacteriophage identification tools have been developed that differ in the phage mining strategies used, input files requested, and results produced. However, new users attempting bacteriophage analysis can struggle to select the best methods and interpret the variety of results produced. Here, we present MetaPhage, a comprehensive reads-to-report pipeline that streamlines the use of multiple phage miners and generates an exhaustive report. The report both summarizes and visualizes the key findings and enables further exploration of key results via interactive filterable tables. The pipeline is implemented in Nextflow, a widely adopted workflow manager that enables an optimized parallelization of tasks in different locations, from local server to the cloud; this ensures reproducible results from containerized packages. MetaPhage is designed to enable scalability and reproducibility; also, it can be easily expanded to include new miners and methods as they are developed in this continuously growing field. MetaPhage is freely available under a GPL-3.0 license at https://github.com/MattiaPandolfoVR/MetaPhage. IMPORTANCE Bacteriophages (viruses that infect bacteria) are the most abundant biological entities on earth and are increasingly studied as members of the resident microbiota community in many environments, from oceans to soils and the human gut. Their identification is of great importance to better understand complex bacterial dynamics and microbial ecosystem function. A variety of metagenome bacteriophage identification tools have been developed that differ in the phage mining strategies used, input files requested, and results produced. To facilitate the management and the execution of such a complex workflow, we developed MetaPhage (MP), a comprehensive reads-to-report pipeline that streamlines the use of multiple phage miners and generates an exhaustive report. The pipeline is implemented in Nextflow, a widely adopted workflow manager that enables an optimized parallelization of tasks. MetaPhage is designed to enable scalability and reproducibility and offers an installation-free, dependency-free, and conflict-free workflow execution.

Crop management shapes the diversity and activity of DNA and RNA viruses in the rhizosphere

BACKGROUND

The rhizosphere is a hotspot for microbial activity and contributes to ecosystem services including plant health and biogeochemical cycling. The activity of microbial viruses, and their influence on plant-microbe interactions in the rhizosphere, remains undetermined. Given the impact of viruses on the ecology and evolution of their host communities, determining how soil viruses influence microbiome dynamics is crucial to build a holistic understanding of rhizosphere functions.

RESULTS

Here, we aimed to investigate the influence of crop management on the composition and activity of bulk soil, rhizosphere soil, and root viral communities. We combined viromics, metagenomics, and metatranscriptomics on soil samples collected from a 3-year crop rotation field trial of oilseed rape (Brassica napus L.). By recovering 1059 dsDNA viral populations and 16,541 ssRNA bacteriophage populations, we expanded the number of underexplored Leviviricetes genomes by > 5 times. Through detection of viral activity in metatranscriptomes, we uncovered evidence of "Kill-the-Winner" dynamics, implicating soil bacteriophages in driving bacterial community succession. Moreover, we found the activity of viruses increased with proximity to crop roots, and identified that soil viruses may influence plant-microbe interactions through the reprogramming of bacterial host metabolism. We have provided the first evidence of crop rotation-driven impacts on soil microbial communities extending to viruses. To this aim, we present the novel principal of "viral priming," which describes how the consecutive growth of the same crop species primes viral activity in the rhizosphere through local adaptation.

CONCLUSIONS

Overall, we reveal unprecedented spatial and temporal diversity in viral community composition and activity across root, rhizosphere soil, and bulk soil compartments. Our work demonstrates that the roles of soil viruses need greater consideration to exploit the rhizosphere microbiome for food security, food safety, and environmental sustainability. Video Abstract.

Antimicrobial Treatment of Pseudomonas aeruginosa Severe Sepsis

Pseudomonas aeruginosa is a pathogen often encountered in a healthcare setting. It has consistently ranked among the most frequent pathogens seen in nosocomial infections, particularly bloodstream and respiratory tract infections. Aside from having intrinsic resistance to many antibiotics, it rapidly acquires resistance to novel agents. Given the high mortality of pseudomonal infections generally, and pseudomonal sepsis particularly, and with the rise of resistant strains, treatment can be very challenging for the clinician. In this paper, we will review the latest evidence for the optimal treatment of P. aeruginosa sepsis caused by susceptible as well as multidrug-resistant strains including the difficult to treat pathogens. We will also discuss the mode of drug infusion, indications for combination therapy, along with the proper dosing and duration of treatment for various conditions with a brief discussion of the use of non-antimicrobial agents.

Phylogenic analysis of new viral cluster of large phages with unusual DNA genomes containing uracil in place of thymine in gene-sharing network, using phages S6 and PBS1 and relevant uncultured phages derived from sewage metagenomics

Bacteriophages (phages) are the most diverse and abundant life-form on Earth. Jumbophages are phages with double-stranded DNA genomes longer than 200 kbp. Among these, some jumbophages with uracil in place of thymine as a nucleic acid base, which we have tentatively termed "dU jumbophages" in this study, have been reported. Because the dU jumbophages are considered to be a living fossil from the RNA world, the evolutionary traits of dU jumbophages are of interest. In this study, we examined the phylogeny of dU jumbophages. First, tBLASTx analysis of newly sequenced dU jumbophages such as Bacillus phage PBS1 and previously isolated Staphylococcus phage S6 showed similarity to the other dU jumbophages. Second, we detected the two partial genome sequences of uncultured phages possibly relevant to dU jumbophages, scaffold_002 and scaffold_007, from wastewater metagenomics. Third, according to the gene-sharing network analysis, the dU jumbophages, including phages PBS1 and S6, and uncultured phage scaffold_002 formed a cluster, which suggested a new viral subfamily/family. Finally, analyses of the phylogenetic relationship with other phages showed that the dU jumbophage cluster, which had two clades of phages infecting Gram-negative and Gram-positive bacteria, diverged from the single ancestral phage. These findings together with previous reports may imply that dU jumbophages evolved from the same origin before divergence of Gram-negative and Gram-positive bacteria.

Membrane lipid renovation in Pseudomonas aeruginosa - implications for phage therapy?

Pseudomonas aeruginosa is an important Gram-negative pathogen with intrinsic resistance to many clinically used antibiotics. It is particularly troublesome in nosocomial infections, immunocompromised patients, and individuals with cystic fibrosis. Antimicrobial resistance (AMR) is a huge threat to global health, with a predicted 10 million people dying from resistant infections by 2050. A promising therapy for combatting AMR infections is phage therapy. However, more research is required to investigate mechanisms that may influence the efficacy of phage therapy. An important overlooked aspect is the impact of membrane lipid remodelling on phage binding ability. P. aeruginosa undergoes changes in membrane lipids when it encounters phosphorus stress, an environmental perturbation that is likely to occur during infection. Lipid changes include the substitution of glycerophospholipids with surrogate glycolipids and the over-production of ornithine-containing aminolipids. Given that membrane lipids are known to influence the structure and function of membrane proteins, we propose that changes in the composition of membrane lipids during infection may alter phage binding and subsequent phage infection dynamics. Consideration of such effects needs to be urgently prioritised in order to develop the most effective phage therapy strategies for P. aeruginosa infections.

Validation of the application of gel beads-based single-cell genome sequencing platform to soil and seawater

Single-cell genomics is applied to environmental samples as a method to solve the problems of current metagenomics. However, in the fluorescence-activated cell sorting-based cell isolation and subsequent whole genome amplification, the sorting efficiency and the sequence quality are greatly affected by the type of target environment, limiting its adaptability. Here, we developed an improved single-cell genomics platform, named SAG-gel, which utilizes gel beads for single-cell isolation, lysis, and whole genome amplification. To validate the versatility of SAG-gel, single-cell genome sequencing was performed with model bacteria and microbial samples collected from eight environmental sites, including soil and seawater. Gel beads enabled multiple lysis treatments. The genome coverage with model bacteria was improved by 9.1-25%. A total of 734 single amplified genomes were collected from the diverse environmental samples, and almost full-length 16S rRNA genes were recovered from 57.8% of them. We also revealed two marine Rhodobacter strains harboring nearly identical 16S rRNA genes but having different genome contents. In addition, searching for viral sequences elucidated the virus-host linkage over the sampling sites, revealing the geographic distribution and diverse host range of viruses.

The transcriptional regulator CtrA controls gene expression in Alphaproteobacteria phages: Evidence for a lytic deferment pathway

Pilitropic and flagellotropic phages adsorb to bacterial pili and flagella. These phages have long been used to investigate multiple aspects of bacterial physiology, such as the cell cycle control in the Caulobacterales. Targeting cellular appendages for adsorption effectively constrains the population of infectable hosts, suggesting that phages may have developed strategies to maximize their infective yield. Brevundimonas phage vB_BsubS-Delta is a recently characterized pilitropic phage infecting the Alphaproteobacterium Brevundimonas subvibrioides. Like other Caulobacterales, B. subvibrioides divides asymmetrically and its cell cycle is governed by multiple transcriptional regulators, including the master regulator CtrA. Genomic characterization of phage vB_BsubS-Delta identified the presence of a large intergenic region with an unusually high density of putative CtrA-binding sites. A systematic analysis of the positional distribution of predicted CtrA-binding sites in complete phage genomes reveals that the highly skewed distribution of CtrA-binding sites observed in vB_BsubS-Delta is an unequivocal genomic signature that extends to other pilli- and flagellotropic phages infecting the Alphaproteobacteria. Moreover, putative CtrA-binding sites in these phage genomes localize preferentially to promoter regions and have higher scores than those detected in other phage genomes. Phylogenetic and comparative genomics analyses show that this genomic signature has evolved independently in several phage lineages, suggesting that it provides an adaptive advantage to pili/flagellotropic phages infecting the Alphaproteobacteria. Experimental results demonstrate that CtrA binds to predicted CtrA-binding sites in promoter regions and that it regulates transcription of phage genes in unrelated Alphaproteobacteria-infecting phages. We propose that this focused distribution of CtrA-binding sites reflects a fundamental new aspect of phage infection, which we term lytic deferment. Under this novel paradigm, pili- and flagellotropic phages exploit the CtrA transduction pathway to monitor the host cell cycle state and synchronize lysis with the presence of infectable cells.