Prokaryotic Virus Orthologous Groups (pVOGs): a resource for comparative genomics and protein family annotation

PCR Assay for Rapid Taxonomic Differentiation of Virulent Staphylococcus aureus and Klebsiella pneumoniae Bacteriophages

Phage therapy is now seen as a promising way to overcome the current global crisis in the spread of multidrug-resistant bacteria. However, phages are highly strain-specific, and in most cases one will have to isolate a new phage or search for a phage suitable for a therapeutic application in existing libraries. At an early stage of the isolation process, rapid screening techniques are needed to identify and type potential virulent phages. Here, we propose a simple PCR approach to differentiate between two families of virulent Staphylococcus phages (Herelleviridae and Rountreeviridae) and eleven genera of virulent Klebsiella phages (Przondovirus, Taipeivirus, Drulisvirus, Webervirus, Jiaodavirus, Sugarlandvirus, Slopekvirus, Jedunavirus, Marfavirus, Mydovirus and Yonseivirus). This assay includes a thorough search of a dataset comprising S. aureus (n = 269) and K. pneumoniae (n = 480) phage genomes available in the NCBI RefSeq/GenBank database for specific genes that are highly conserved at the taxonomic group level. The selected primers showed high sensitivity and specificity for both isolated DNA and crude phage lysates, which permits circumventing DNA purification protocols. Our approach can be extended and applied to any group of phages, given the large number of available genomes in the databases.

High Level of Interaction between Phages and Bacteria in an Artisanal Raw Milk Cheese Microbial Community

Microbial starter cultures are used in the production of many cheeses around the world, such as Parmigiano-Reggiano, in Italy, Époisses, in France, and Canastra, in Brazil, providing many of the unique features of these cheeses. Bacteriophages (phages) are ubiquitous and well known to modulate the structure of bacterial communities, and recent data indicate that cheeses contain a high abundance of naturally occurring phages. Here, we analyze the viral and bacterial metagenomes of Canastra cheese: a traditional artisanal Brazilian cheese produced using an endogenous starter culture and raw milk. Over 1,200 viral operational taxonomic units were recovered using both isolated viral-like particles and complete metagenomic DNA. Common viral families identified included Siphoviridae and Myoviridae, with 40% of putative phage genomes unidentified at the family level of classification. We observed very high phage diversity, which varied greatly across different cheese producers, with 28% of phage genomes detected in only one producer. Several metagenome-assembled genomes were recovered for lactic acid-producing bacteria, as well as nonstarter bacterial species, and we identified several phage-bacterium interactions, at the strain level of resolution, varying across distinct cheese producers. We postulate that at least one bacterial strain detected could be endogenous and unique to the Canastra cheese-producing region in Brazil and that its growth seems to be modulated by autochthonous phages present in this artisanal production system. This phage-host relationship is likely to influence the fermentation dynamics and ultimately the sensorial profile of these cheeses, with implications for other similar cheese production systems around the world. IMPORTANCE Our work demonstrated a dynamic yet stable microbial ecosystem during cheese production using an endogenous starter culture. This was observed across several distinct producers and was marked by genomic evidence of continued phage-bacterium interactions, such as the presence of bacterial defense mechanisms. Furthermore, we provide evidence of unique microbial signatures for each individual cheese producer studied in the region, a fact that may have profound consequences on product traceability. This was the first effort to describe and understand the bacteriophage composition and ecological dynamics within the Brazilian Canastra cheese production system. The study of this prototypical backslopping production system provides a solid background for further mechanistic studies of the production of many cheeses around the world.

Awaking the dormant virome in the rhizosphere

The rhizosphere is a vital soil compartment providing key plant-beneficial functions. However, little is known about the mechanisms driving viral diversity in the rhizosphere. Viruses can establish lytic or lysogenic interactions with their bacterial hosts. In the latter, they assume a dormant state integrated in the host genome and can be awakened by different perturbations that impact host cell physiology, triggering a viral bloom, which is potentially a fundamental mechanism driving soil viral diversity, as 22%-68% of soil bacteria are predicted to harbour dormant viruses. Here we assessed the viral bloom response in rhizospheric viromes by exposing them to three contrasting soil perturbation agents: earthworms, herbicide and antibiotic pollutant. The viromes were next screened for rhizosphere-relevant genes and also used as inoculant on microcosms incubations to test their impacts on pristine microbiomes. Our results show that while post-perturbation viromes diverged from control conditions, viral communities exposed to both herbicide and antibiotic pollutant were more similar to each other than those influenced by earthworms. The latter also favoured an increase in viral populations harbouring genes involved in plant-beneficial functions. Post-perturbation viromes inoculated on soil microcosms changed the diversity of pristine microbiomes, suggesting that viromes are important components of the soil ecological memory driving eco-evolutionary processes that determine future microbiome trajectories according to past events. Our findings demonstrate that viromes are active players in the rhizosphere and need to be considered in efforts to understand and control the microbial processes towards sustainable crop production.

Rational Design of Profile HMMs for Sensitive and Specific Sequence Detection with Case Studies Applied to Viruses, Bacteriophages, and Casposons

Profile hidden Markov models (HMMs) are a powerful way of modeling biological sequence diversity and constitute a very sensitive approach to detecting divergent sequences. Here, we report the development of protocols for the rational design of profile HMMs. These methods were implemented on TABAJARA, a program that can be used to either detect all biological sequences of a group or discriminate specific groups of sequences. By calculating position-specific information scores along a multiple sequence alignment, TABAJARA automatically identifies the most informative sequence motifs and uses them to construct profile HMMs. As a proof-of-principle, we applied TABAJARA to generate profile HMMs for the detection and classification of two viral groups presenting different evolutionary rates: bacteriophages of the Microviridae family and viruses of the Flavivirus genus. We obtained conserved models for the generic detection of any Microviridae or Flavivirus sequence, and profile HMMs that can specifically discriminate Microviridae subfamilies or Flavivirus species. In another application, we constructed Cas1 endonuclease-derived profile HMMs that can discriminate CRISPRs and casposons, two evolutionarily related transposable elements. We believe that the protocols described here, and implemented on TABAJARA, constitute a generic toolbox for generating profile HMMs for the highly sensitive and specific detection of sequence classes.

When Plaquing Is Not Possible: Computational Methods for Detecting Induced Phages

High-throughput sequencing of microbial communities has uncovered a large, diverse population of phages. Frequently, phages found are integrated into their bacterial host genome. Distinguishing between phages in their integrated (lysogenic) and unintegrated (lytic) stage can provide insight into how phages shape bacterial communities. Here we present the Prophage Induction Estimator (PIE) to identify induced phages in genomic and metagenomic sequences. PIE takes raw sequencing reads and phage sequence predictions, performs read quality control, read assembly, and calculation of phage and non-phage sequence abundance and completeness. The distribution of abundances for non-phage sequences is used to predict induced phages with statistical confidence. In silico tests were conducted to benchmark this tool finding that PIE can detect induction events as well as phages with a relatively small burst size (10×). We then examined isolate genome sequencing data as well as a mock community and urinary metagenome data sets and found instances of induced phages in all three data sets. The flexibility of this software enables users to easily include phage predictions from their preferred tool of choice or phage sequences of interest. Thus, genomic and metagenomic sequencing now not only provides a means for discovering and identifying phage sequences but also the detection of induced prophages.

Bacteriophages playing nice: Lysogenic bacteriophage replication stable in the human gut microbiota

Bacteriophages, viruses specific to bacteria, coexist with their bacterial hosts with limited diversity fluctuations in the guts of healthy individuals where they replicate mostly via lysogenic replication. This favors 'piggy-back-the-winner' over 'kill-the-winner' dynamics which are driven by lytic bacteriophage replication. Revisiting the deep-viral sequencing data of a healthy individual sampled over 2.4 years, we explore how these dynamics occur. Prophages found in assembled bacterial metagenomes were also found extra-cellularly, as induced phage particles (iPPs), likely derived from prophage activation. These iPPs were diverse and continually present in low abundance, relative to the highly abundant but less diverse lytic phage population. The continuous detection of low levels of iPPs suggests that spontaneous induction regularly occurs in this healthy individual, possibly allowing prophages to maintain their ability to replicate and avoiding degradation and loss from the gut microbiota.

Temperate bacteriophages infecting the mucin-degrading bacterium Ruminococcus gnavus from the human gut

Ruminococcus gnavus is a prevalent gut microbe reported to occur in higher abundance among individuals with inflammatory bowel disease (IBD). This study reports the isolation and characterization of six bacteriophages (phages) isolated from human fecal material and environmental samples that infect this species. Isolated phages have a siphovirus morphology, with genomes ranging between 36.5 and 37.8 kbp. Genome analysis indicates that the phages have a temperate lifestyle, which was confirmed by their ability to form lysogens on their host bacterial species. In contrast to the finding that phages lyse their host in liquid medium, results from a mouse trial indicate these phages can co-exist with the host bacterium in the gut without causing a significant reduction of R. gnavus. The bacterial counts in the feces of phage-treated mice did not significantly differ in the presence of phage. Furthermore, analysis of publicly available gut virome sequence data indicates a high abundance of these phages among individuals suffering from IBD. This work provides the first insight into how phages interact with R. gnavus in the human gut microbiome.

Multiomic spatial analysis reveals a distinct mucosa-associated virome

The human gut virome has been increasingly explored in recent years. However, nearly all virome-sequencing efforts rely solely on fecal samples and few studies leverage multiomic approaches to investigate phage-host relationships. Here, we combine metagenomics, metaviromics, and metatranscriptomics to study virome-bacteriome interactions at the colonic mucosal-luminal interface in a cohort of three individuals with inflammatory bowel disease; non-IBD controls were not included in this study. We show that the mucosal viral population is distinct from the stool virome and houses abundant crAss-like phages that are undetectable by fecal sampling. Through viral protein prediction and metatranscriptomic analysis, we explore viral gene transcription, prophage activation, and the relationship between the presence of integrase and temperate phages in IBD subjects. We also show the impact of deep sequencing on virus recovery and offer guidelines for selecting optimal sequencing depths in future metaviromic studies. Systems biology approaches such as those presented in this report will enhance our understanding of the human virome and its interactions with our microbiome and our health.

Metagenomic analysis reveals unexplored diversity of archaeal virome in the human gut

The human gut microbiome has been extensively explored, while the archaeal viruses remain largely unknown. Here, we present a comprehensive analysis of the archaeal viruses from the human gut metagenomes and the existing virus collections using the CRISPR spacer and viral signature-based approach. This results in 1279 viral species, of which, 95.2% infect Methanobrevibacteria_A, 56.5% shared high identity (>95%) with the archaeal proviruses, 37.2% have a host range across archaeal species, and 55.7% are highly prevalent in the human population (>1%). A methanogenic archaeal virus-specific gene for pseudomurein endoisopeptidase (PeiW) frequently occurs in the viral sequences (n = 150). Analysis of 33 Caudoviricetes viruses with a complete genome often discovers the genes (integrase, n = 29; mazE, n = 10) regulating the viral lysogenic-lytic cycle, implying the dominance of temperate viruses in the archaeal virome. Together, our work uncovers the unexplored diversity of archaeal viruses, revealing the novel facet of the human gut microbiome.

Genomic Insights into the Mobilome and Resistome of Sentinel Microorganisms Originating from Farms of Two Different Swine Production Systems

Antimicrobial resistance (AMR) is a threat to public health due to long-term antimicrobial use (AMU), which promotes the bacterial acquisition of antimicrobial resistance determinants (ARDs). Within food-producing animals, organic and extensive Iberian swine production is based on sustainable and eco-friendly management systems, providing an excellent opportunity to evaluate how sustained differences in AMU impact the development and spread of AMR. Here, through a whole-genome sequencing approach, we provide an in-depth characterization of the resistome and mobilome and their interaction in 466 sentinel bacteria, namely, Escherichia coli, Enterococcus spp., Campylobacter coli, and Staphylococcus spp., recovered from 37 intensive and organic-extensive pig farms. Both ARDs and mobile genetic elements (MGEs) were primarily taxon-associated, with higher similarities among bacteria which were closely phylogenetically related. E. coli exhibited the most diverse resistome and mobilome, with 85.4% mobilizable ARDs, 50.3% of which were plasmid-associated. Staphylococcus spp. exhibited a broad repertoire of ARDs and MGEs, with 52.3% of its resistome being mobilizable. Although Enterococcus spp. carried the highest number of ARDs per isolate and its plasmidome was similar in size to that of E. coli, 43.7% of its resistome was mobilizable. A narrow spectrum of ARDs constituted the C. coli resistome, with point mutations as its main AMR driver. A constrained AMU, as observed in organic-extensive herds, determined a reduction in the quantitative composition of the resistome and the complexity of the resistome-mobilome interaction. These results demonstrate taxon-associated AMR-MGE interactions and evidence that responsible AMU can contribute to reducing AMR pressure in the food chain. IMPORTANCE This study provides the first integral genomic characterization of the resistome and mobilome of sentinel microorganisms for antimicrobial resistance (AMR) surveillance from two different swine production systems. Relevant differences were observed among taxa in the resistomes and mobilomes they harbored, revealing their distinctive risk in AMR dissemination and spread. Thus, Escherichia coli and, to a lesser extent, Staphylococcus spp. constituted the main reservoirs of mobilizable antimicrobial resistance genes, which were predominantly plasmid-associated; in contrast to Campylobacter coli, whose resistome was mainly determined by point mutations. The reduced complexity of mobilome-resistome interaction in Enterococcus spp. suggested its limited role in AMR dissemination from swine farms. The significant differences in antimicrobial use among the studied farms allowed us to assess the suitability of whole-genome sequencing as a rapid and efficient technique for the assessment of mid- to long-term on-farm interventions for the reduction of antimicrobial use and the evaluation of AMR status.

vHULK, a New Tool for Bacteriophage Host Prediction Based on Annotated Genomic Features and Neural Networks

Background

The experimental determination of a bacteriophage host is a laborious procedure. Thus, there is a pressing need for reliable computational predictions of bacteriophage hosts.

Materials and Methods

We developed the program vHULK for phage host prediction based on 9504 phage genome features, which consider alignment significance scores between predicted proteins and a curated database of viral protein families. The features were fed to a neural network, and two models were trained to predict 77 host genera and 118 host species.

Results

In controlled random test sets with 90% redundancy reduction in terms of protein similarity, vHULK obtained on average 83% precision and 79% recall at the genus level, and 71% precision and 67% recall at the species level. The performance of vHULK was compared against three other tools on a test data set with 2153 phage genomes. On this data set, vHULK achieved better performance at both the genus and the species levels than the other tools.

Conclusions

Our results suggest that vHULK represents an advance on the state of art in phage host prediction.

Prokaryotic and viral genomes recovered from 787 Japanese gut metagenomes revealed microbial features linked to diets, populations, and diseases

We reconstructed 19,084 prokaryotic and 31,395 viral genomes from 787 Japanese gut metagenomes as Japanese metagenome-assembled genomes (JMAG) and Japanese Virus Database (JVD), which are large microbial genome datasets for a single population. Population-specific enrichment of the Bacillus subtilis and β-porphyranase among the JMAG could derive from the Japanese traditional food natto (fermented soybeans) and nori (laver), respectively. Dairy-related Enterococcus_B lactis and Streptococcus thermophilus were nominally associated with the East Asian-specific missense variant rs671:G>A in ALDH2, which was associated with dairy consumption. Of the species-level viral genome clusters in the JVD, 62.9% were novel. The β crAss-like phage composition was low among the Japanese but relatively high among African and Oceanian peoples. Evaluations of the association between crAss-like phages and diseases showed significant disease-specific associations. Our large catalog of virus-host pairs identified the positive correlation between the abundance of the viruses and their hosts.

Comparing In Vitro Faecal Fermentation Methods as Surrogates for Phage Therapy Application

The human microbiome and its importance in health and disease have been the subject of numerous research articles. Most microbes reside in the digestive tract, with up to 10¹² cells per gram of faecal material found in the colon. In terms of gene number, it has been estimated that the gut microbiome harbours >100 times more genes than the human genome. Several human intestinal diseases are strongly associated with disruptions in gut microbiome composition. Less studied components of the gut microbiome are the bacterial viruses called bacteriophages that may be present in numbers equal to or greater than the prokaryotes. Their potential to lyse their bacterial hosts, or to act as agents of horizontal gene transfer makes them important research targets. In this study in vitro faecal fermentation systems were developed and compared for their ability to act as surrogates for the human colon. Changes in bacterial and viral composition occurred after introducing a high-titre single phage preparation both with and without a known bacterial host during the 24 h-long fermentation. We also show that during this timeframe 50 mL plastic tubes can provide data similar to that generated in a sophisticated faecal fermenter system. This knowledge can guide us to a better understanding of the short-term impact of bacteriophage transplants on the bacteriomes and viromes of human recipients.

What the Phage: a scalable workflow for the identification and analysis of phage sequences

Phages are among the most abundant and diverse biological entities on earth. Phage prediction from sequence data is a crucial first step to understanding their impact on the environment. A variety of bacteriophage prediction tools have been developed over the years. They differ in algorithmic approach, results, and ease of use. We, therefore, developed "What the Phage" (WtP), an easy-to-use and parallel multitool approach for phage prediction combined with an annotation and classification downstream strategy, thus supporting the user's decision-making process by summarizing the results of the different prediction tools in charts and tables. WtP is reproducible and scales to thousands of datasets through a workflow manager (Nextflow). WtP is freely available under a GPL-3.0 license (https://github.com/replikation/What_the_Phage).

Multi-kingdom gut microbiota analyses define COVID-19 severity and post-acute COVID-19 syndrome

Our knowledge of the role of the gut microbiome in acute coronavirus disease 2019 (COVID-19) and post-acute COVID-19 is rapidly increasing, whereas little is known regarding the contribution of multi-kingdom microbiota and host-microbial interactions to COVID-19 severity and consequences. Herein, we perform an integrated analysis using 296 fecal metagenomes, 79 fecal metabolomics, viral load in 1378 respiratory tract samples, and clinical features of 133 COVID-19 patients prospectively followed for up to 6 months. Metagenomic-based clustering identifies two robust ecological clusters (hereafter referred to as Clusters 1 and 2), of which Cluster 1 is significantly associated with severe COVID-19 and the development of post-acute COVID-19 syndrome. Significant differences between clusters could be explained by both multi-kingdom ecological drivers (bacteria, fungi, and viruses) and host factors with a good predictive value and an area under the curve (AUC) of 0.98. A model combining host and microbial factors could predict the duration of respiratory viral shedding with 82.1% accuracy (error ± 3 days). These results highlight the potential utility of host phenotype and multi-kingdom microbiota profiling as a prognostic tool for patients with COVID-19.

Viral community-wide auxiliary metabolic genes differ by lifestyles, habitats, and hosts

BACKGROUND

Viral-encoded auxiliary metabolic genes (AMGs) are important toolkits for modulating their hosts' metabolisms and the microbial-driven biogeochemical cycles. Although the functions of AMGs have been extensively reported in numerous environments, we still know little about the drivers that shape the viral community-wide AMG compositions in natural ecosystems. Exploring the drivers of viral community-wide AMG compositions is critical for a deeper understanding of the complex interplays among viruses, hosts, and the environments.

RESULTS

Here, we investigated the impact of viral lifestyles (i.e., lytic and lysogenic), habitats (i.e., water, particle, and sediment), and prokaryotic hosts on viral AMG profiles by utilizing metagenomic and metatranscriptomic techniques. We found that viral lifestyles were the most important drivers, followed by habitats and host identities. Specifically, irrespective of what habitats viruses came from, lytic viruses exhibited greater AMG diversity and tended to encode AMGs for chaperone biosynthesis, signaling proteins, and lipid metabolism, which could boost progeny reproduction, whereas temperate viruses were apt to encode AMGs for host survivability. Moreover, the lytic and temperate viral communities tended to mediate the microbial-driven biogeochemical cycles, especially nitrogen metabolism, in different manners via AMGs. When focusing on each lifestyle, we further found clear dissimilarity in AMG compositions between water and sediment, as well the divergent AMGs encoded by viruses infecting different host orders.

CONCLUSIONS

Overall, our study provides a first systematic characterization of the drivers of viral community-wide AMG compositions and further expands our knowledge of the distinct interactions of lytic and temperate viruses with their prokaryotic hosts from an AMG perspective, which is critical for understanding virus-host-environment interactions in natural conditions. Video Abstract.

PhaMMseqs: a new pipeline for constructing phage gene phamilies using MMseqs2

The diversity and mosaic architecture of phage genomes present challenges for whole-genome phylogenies and comparative genomics. There are no universally conserved core genes, ∼70% of phage genes are of unknown function, and phage genomes are replete with small (<500 bp) open reading frames. Assembling sequence-related genes into "phamilies" ("phams") based on amino acid sequence similarity simplifies comparative phage genomics and facilitates representations of phage genome mosaicism. With the rapid and substantial increase in the numbers of sequenced phage genomes, computationally efficient pham assembly is needed, together with strategies for including newly sequenced phage genomes. Here, we describe the Python package PhaMMseqs, which uses MMseqs2 for pham assembly, and we evaluate the key parameters for optimal pham assembly of sequence- and functionally related proteins. PhaMMseqs runs efficiently with only modest hardware requirements and integrates with the pdm_utils package for simple genome entry and export of datasets for evolutionary analyses and phage genome map construction.

Dynamic impact of virome on colitis and colorectal cancer: Immunity, inflammation, prevention and treatment

The gut microbiome includes a series of microorganism genomes, such as bacteriome, virome, mycobiome, etc. The gut microbiota is critically involved in intestine immunity and diseases, including inflammatory bowel disease (IBD) and colorectal cancer (CRC); however, the underlying mechanism remains incompletely understood. Clarifying the relationship between microbiota and inflammation may profoundly improve our understanding of etiology, disease progression, patient management, and the development of prevention and treatment. In this review, we discuss the latest studies of the influence of enteric viruses (i.e., commensal viruses, pathogenic viruses, and bacteriophages) in the initiation, progression, and complication of colitis and colorectal cancer, and their potential for novel preventative approaches and therapeutic application. We explore the interplay between gut viruses and host immune systems for its effects on the severity of inflammatory diseases and cancer, including both direct and indirect interactions between enteric viruses with other microbes and microbial products. Furthermore, the underlying mechanisms of the virome's roles in gut inflammatory response have been explained to infer potential therapeutic targets with examples in specific clinical trials. Given that very limited literature has thus far discussed these various topics with the gut virome, we believe these extensive analyses may provide insight into the understanding of the molecular pathogenesis of IBD and CRC, which could help add the design of improved therapies for these important human diseases.

Phage-Plasmids Spread Antibiotic Resistance Genes through Infection and Lysogenic Conversion

Antibiotic resistance is rapidly spreading via the horizontal transfer of resistance genes in mobile genetic elements. While plasmids are key drivers of this process, few integrative phages encode antibiotic resistance genes. Here, we find that phage-plasmids, elements that are both phages and plasmids, often carry antibiotic resistance genes. We found 60 phage-plasmids with 184 antibiotic resistance genes, providing resistance for broad-spectrum-cephalosporins, carbapenems, aminoglycosides, fluoroquinolones, and colistin. These genes are in a few hot spots, seem to have been cotranslocated with transposable elements, and are often in class I integrons, which had not been previously found in phages. We tried to induce six phage-plasmids with resistance genes (including four with resistance integrons) and succeeded in five cases. Other phage-plasmids and integrative prophages were coinduced in these experiments. As a proof of concept, we focused on a P1-like element encoding an extended spectrum β-lactamase, blaCTX-M-55. After induction, we confirmed that it is capable of infecting and converting four other E. coli strains. Its reinduction led to the further conversion of a sensitive strain, confirming that it is a fully functional phage. This study shows that phage-plasmids carry a large diversity of clinically relevant antibiotic resistance genes that they can transfer across bacteria. As plasmids, these elements seem plastic and capable of acquiring genes from other plasmids. As phages, they may provide novel paths of transfer for resistance genes because they can infect bacteria that are distant in time and space from the original host. As a matter of alarm, they may also mediate transfer to other types of phages. IMPORTANCE The dissemination of antimicrobial resistance is a major threat to global health. Here, we show that a group of temperate bacterial viruses (phages), termed phage-plasmids, commonly encode different and multiple types of resistance genes of high clinical importance, often in integrons. This is unexpected, as phages typically do not carry resistance genes and, hence, do not confer upon their hosts resistance via infection and genome integration. Our experiments with phage-plasmids isolated from clinical settings confirmed that they infect sensitive strains and render them antibiotic resistant. The spread of antibiotic resistance genes by phage-plasmids is worrisome because it dispenses cell-to-cell contact, which is necessary for canonical plasmid transfer (conjugation). Furthermore, their integrons become genetic platforms for the acquisition of novel resistance genes.

Diversity and Ecology of Caudoviricetes Phages with Genome Terminal Repeats in Fecal Metagenomes from Four Dutch Cohorts

The human gut harbors numerous viruses infecting the human host, microbes, and other inhabitants of the gastrointestinal tract. Most of these viruses remain undiscovered, and their influence on human health is unknown. Here, we characterize viral genomes in gut metagenomic data from 1950 individuals from four population and patient cohorts. We focus on a subset of viruses that is highly abundant in the gut, remains largely uncharacterized, and allows confident complete genome identification—phages that belong to the class Caudoviricetes and possess genome terminal repeats. We detect 1899 species-level units belonging to this subset, 19% of which do not have complete representative genomes in major public gut virome databases. These units display diverse genomic features, are predicted to infect a wide range of microbial hosts, and on average account for <1% of metagenomic reads. Analysis of longitudinal data from 338 individuals shows that the composition of this fraction of the virome remained relatively stable over a period of 4 years. We also demonstrate that 54 species-level units are highly prevalent (detected in >5% of individuals in a cohort). Finally, we find 34 associations between highly prevalent phages and human phenotypes, 24 of which can be explained by the relative abundance of potential hosts.

A functional selection reveals previously undetected anti-phage defence systems in the E. coli pangenome

The ancient, ongoing coevolutionary battle between bacteria and their viruses, bacteriophages, has given rise to sophisticated immune systems including restriction-modification and CRISPR-Cas. Many additional anti-phage systems have been identified using computational approaches based on genomic co-location within defence islands, but these screens may not be exhaustive. Here we developed an experimental selection scheme agnostic to genomic context to identify defence systems in 71 diverse E. coli strains. Our results unveil 21 conserved defence systems, none of which were previously detected as enriched in defence islands. Additionally, our work indicates that intact prophages and mobile genetic elements are primary reservoirs and distributors of defence systems in E. coli, with defence systems typically carried in specific locations or hotspots. These hotspots encode dozens of additional uncharacterized defence system candidates. Our findings reveal an extended landscape of antiviral immunity in E. coli and provide an approach for mapping defence systems in other species.

Development and Validation of a Bioinformatic Workflow for the Rapid Detection of Viruses in Biosecurity

The field of biosecurity has greatly benefited from the widespread adoption of high-throughput sequencing technologies, for its ability to deeply query plant and animal samples for pathogens for which no tests exist. However, the bioinformatics analysis tools designed for rapid analysis of these sequencing datasets are not developed with this application in mind, limiting the ability of diagnosticians to standardise their workflows using published tool kits. We sought to assess previously published bioinformatic tools for their ability to identify plant- and animal-infecting viruses while distinguishing from the host genetic material. We discovered that many of the current generation of virus-detection pipelines are not adequate for this task, being outperformed by more generic classification tools. We created synthetic MinION and HiSeq libraries simulating plant and animal infections of economically important viruses and assessed a series of tools for their suitability for rapid and accurate detection of infection, and further tested the top performing tools against the VIROMOCK Challenge dataset to ensure that our findings were reproducible when compared with international standards. Our work demonstrated that several methods provide sensitive and specific detection of agriculturally important viruses in a timely manner and provides a key piece of ground truthing for method development in this space.

Insights into the global freshwater virome

Viruses are by far the most abundant life forms on this planet. Yet, the full viral diversity remains mostly unknown, especially in environments like freshwater. Therefore, we aimed to study freshwater viruses in a global context. To this end, we downloaded 380 publicly available viral metagenomes (>1 TB). More than 60% of these metagenomes were discarded based on their levels of cellular contamination assessed by ribosomal DNA content. For the remaining metagenomes, assembled contigs were decontaminated using two consecutive steps, eventually yielding 273,365 viral contigs longer than 1,000 bp. Long enough contigs (≥ 10 kb) were clustered to identify novel genomes/genome fragments. We could recover 549 complete circular and high-quality draft genomes, out of which 10 were recognized as being novel. Functional annotation of these genomes showed that most of the annotated coding sequences are DNA metabolic genes or phage structural genes. On the other hand, taxonomic analysis of viral contigs showed that most of the assigned contigs belonged to the order Caudovirales, particularly the families of Siphoviridae, Myoviridae, and Podoviridae. The recovered viral contigs contained several auxiliary metabolic genes belonging to several metabolic pathways, especially carbohydrate and amino acid metabolism in addition to photosynthesis as well as hydrocarbon degradation and antibiotic resistance. Overall, we present here a set of prudently chosen viral contigs, which should not only help better understanding of freshwater viruses but also be a valuable resource for future virome studies.

mobileOG-db: a Manually Curated Database of Protein Families Mediating the Life Cycle of Bacterial Mobile Genetic Elements

Bacterial mobile genetic elements (MGEs) encode functional modules that perform both core and accessory functions for the element, the latter of which are often only transiently associated with the element. The presence of these accessory genes, which are often close homologs to primarily immobile genes, incur high rates of false positives and, therefore, limits the usability of these databases for MGE annotation. To overcome this limitation, we analyzed 10,776,849 protein sequences derived from eight MGE databases to compile a comprehensive set of 6,140 manually curated protein families that are linked to the “life cycle” (integration/excision, replication/recombination/repair, transfer, stability/transfer/defense, and phage-specific processes) of plasmids, phages, integrative, transposable, and conjugative elements. We overlay experimental information where available to create a tiered annotation scheme of high-quality annotations and annotations inferred exclusively through bioinformatic evidence. We additionally provide an MGE-class label for each entry (e.g., plasmid or integrative element), and assign to each entry a major and minor category. The resulting database, mobileOG-db (for mobile orthologous groups), comprises over 700,000 deduplicated sequences encompassing five major mobileOG categories and more than 50 minor categories, providing a structured language and interpretable basis for an array of MGE-centered analyses. mobileOG-db can be accessed at mobileogdb.flsi.cloud.vt.edu/, where users can select, refine, and analyze custom subsets of the dynamic mobilome.

IMPORTANCE The analysis of bacterial mobile genetic elements (MGEs) in genomic data is a critical step toward profiling the root causes of antibiotic resistance, phenotypic or metabolic diversity, and the evolution of bacterial genera. Existing methods for MGE annotation pose high barriers of biological and computational expertise to properly harness. To bridge this gap, we systematically analyzed 10,776,849 proteins derived from eight databases of MGEs to identify 6,140 MGE protein families that can serve as candidate hallmarks, i.e., proteins that can be used as “signatures” of MGEs to aid annotation. The resulting resource, mobileOG-db, provides a multilevel classification scheme that encompasses plasmid, phage, integrative, and transposable element protein families categorized into five major mobileOG categories and more than 50 minor categories. mobileOG-db thus provides a rich resource for simple and intuitive element annotation that can be integrated seamlessly into existing MGE detection pipelines and colocalization analyses.

Analysis of the anaerobic digestion metagenome under environmental stresses stimulating prophage induction

BACKGROUND

The viral community has the potential to influence the structure of the microbiome and thus the yield of the anaerobic digestion process. However, the virome composition in anaerobic digestion is still under-investigated. A viral induction experiment was conducted on separate batches undergoing a series of DNA-damaging stresses, in order to coerce temperate viruses to enter the lytic cycle.

RESULTS

The sequencing of the metagenome revealed a viral community almost entirely composed of tailed bacteriophages of the order Caudovirales. Following a binning procedure 1,092 viral and 120 prokaryotic genomes were reconstructed, 64 of which included an integrated prophage in their sequence. Clustering of coverage profiles revealed the presence of species, both viral and microbial, sharing similar reactions to shocks. A group of viral genomes, which increase under organic overload and decrease under basic pH, uniquely encode the yopX gene, which is involved in the induction of temperate prophages. Moreover, the in-silico functional analysis revealed an enrichment of sialidases in viral genomes. These genes are associated with tail proteins and, as such, are hypothesised to be involved in the interaction with the host. Archaea registered the most pronounced changes in relation to shocks and featured behaviours not shared with other species. Subsequently, data from 123 different samples of the global anaerobic digestion database was used to determine coverage profiles of host and viral genomes on a broader scale.

CONCLUSIONS

Viruses are key components in anaerobic digestion environments, shaping the microbial guilds which drive the methanogenesis process. In turn, environmental conditions are pivotal in shaping the viral community and the rate of induction of temperate viruses. This study provides an initial insight into the complexity of the anaerobic digestion virome and its relation with the microbial community and the diverse environmental parameters. Video Abstract.

New Microviridae isolated from Sulfitobacter reveals two cosmopolitan subfamilies of single-stranded DNA phages infecting marine and terrestrial Alphaproteobacteria

The Microviridae family represents one of the major clades of single-stranded DNA (ssDNA) phages. Their cultivated members are lytic and infect Proteobacteria, Bacteroidetes, and Chlamydiae. Prophages have been predicted in the genomes from Bacteroidales, Hyphomicrobiales, and Enterobacteriaceae and cluster within the 'Alpavirinae', 'Amoyvirinae', and Gokushovirinae. We have isolated 'Ascunsovirus oldenburgi' ICBM5, a novel phage distantly related to known Microviridae. It infects Sulfitobacter dubius SH24-1b and uses both a lytic and a carrier-state life strategy. Using ICBM5 proteins as a query, we uncovered in publicly available resources sixty-five new Microviridae prophages and episomes in bacterial genomes and retrieved forty-seven environmental viral genomes (EVGs) from various viromes. Genome clustering based on protein content and phylogenetic analysis showed that ICBM5, together with Rhizobium phages, new prophages, episomes, and EVGs cluster within two new phylogenetic clades, here tentatively assigned the rank of subfamily and named 'Tainavirinae' and 'Occultatumvirinae'. They both infect Rhodobacterales. Occultatumviruses also infect Hyphomicrobiales, including nitrogen-fixing endosymbionts from cosmopolitan legumes. A biogeographical assessment showed that tainaviruses and occultatumviruses are spread worldwide, in terrestrial and marine environments. The new phage isolated here sheds light onto new and diverse branches of the Microviridae tree, suggesting that much of the ssDNA phage diversity remains in the dark.

Impact of a phage cocktail targeting Escherichia coli and Enterococcus faecalis as members of a gut bacterial consortium in vitro and in vivo

Escherichia coli and Enterococcus faecalis have been implicated as important players in human gut health that have been associated with the onset of inflammatory bowel disease (IBD). Bacteriophage (phage) therapy has been used for decades to target pathogens as an alternative to antibiotics, but the ability of phage to shape complex bacterial consortia in the lower gastrointestinal tract is not clearly understood. We administered a cocktail of six phages (either viable or heat-inactivated) targeting pro-inflammatory Escherichia coli LF82 and Enterococcus faecalis OG1RF as members of a defined community in both a continuous fermenter and a murine colitis model. The two target strains were members of a six species simplified human microbiome consortium (SIHUMI-6). In a 72-h continuous fermentation, the phage cocktail caused a 1.1 and 1.5 log (log10 genome copies/mL) reduction in E. faecalis and E. coli numbers, respectively. This interaction was accompanied by changes in the numbers of other SIHUMI-6 members, with an increase of Lactiplantibacillus plantarum (1.7 log) and Faecalibacterium prausnitzii (1.8 log). However, in germ-free mice colonized by the same bacterial consortium, the same phage cocktail administered twice a week over nine weeks did not cause a significant reduction of the target strains. Mice treated with active or inactive phage had similar levels of pro-inflammatory cytokines (IFN-y/IL12p40) in unstimulated colorectal colonic strip cultures. However, histology scores of the murine lower GIT (cecum and distal colon) were lower in the viable phage-treated mice, suggesting that the phage cocktail did influence the functionality of the SIHUMI-6 consortium. For this study, we conclude that the observed potential of phages to reduce host populations in in vitro models did not translate to a similar outcome in an in vivo setting, with this effect likely brought about by the reduction of phage numbers during transit of the mouse GIT.

Isolation and Molecular Characterization of a Novel Lytic Bacteriophage That Inactivates MDR Klebsiella pneumoniae Strains

The worldwide increase in serious infections caused by multidrug-resistant (MDR) K. pneumoniae emphasizes the urgent need of new therapeutic strategies for the control of this pathogen. There is growing interest in the use of bacteriophages (or phages) to treat K. pneumoniae infections, and newly isolated phages are needed. Here, we report the isolation and physical/biological/molecular characterization of a novel lytic phage and its efficacy in the control of MDR K. pneumoniae. The phage vB_KpnS_Uniso31, referred to hereafter as phage Kpn31, was isolated from hospital wastewater using K. pneumoniae CCCD-K001 as the host. Phage Kpn31 presents a siphovirus-like morphotype and was classified as Demerecviridae; Sugarlandvirus based on its complete genome sequence. The 113,444 bp Kpn31 genome does not encode known toxins or antimicrobial resistance genes, nor does it encode depolymerases related sequences. Phage Kpn31 showed an eclipse time of 15 min and a burst size of 9.12 PFU/host cell, allowing us to conclude it replicates well in K. pneumoniae CCCD-K001 with a latency period of 30 min. Phage Kpn31 was shown to be effective against at least six MDR K. pneumoniae clinical isolates in in vitro antibacterial activity assays. Based on its features, phage Kpn31 has potential for controlling infections caused by MDR K. pneumoniae.

Cell-free production of personalized therapeutic phages targeting multidrug-resistant bacteria

Bacteriophages are potent therapeutics against biohazardous bacteria, which rapidly develop multidrug resistance. However, routine administration of phage therapy is hampered by a lack of rapid production, safe bioengineering, and detailed characterization of phages. Thus, we demonstrate a comprehensive cell-free platform for personalized production, transient engineering, and proteomic characterization of a broad spectrum of phages. Using mass spectrometry, we validated hypothetical and non-structural proteins and could also monitor the protein expression during phage assembly. Notably, a few microliters of a one-pot reaction produced effective doses of phages against enteroaggregative Escherichia coli (EAEC), Yersinia pestis, and Klebsiella pneumoniae. By co-expressing suitable host factors, we could extend the range of cell-free production to phages targeting gram-positive bacteria. We further introduce a non-genomic phage engineering method, which adds functionalities for only one replication cycle. In summary, we expect this cell-free methodology to foster reverse and forward phage engineering and customized production of clinical-grade bacteriophages.

Pangenome Evolution in Environmentally Transmitted Symbionts of Deep-Sea Mussels Is Governed by Vertical Inheritance

Microbial pangenomes vary across species; their size and structure are determined by genetic diversity within the population and by gene loss and horizontal gene transfer (HGT). Many bacteria are associated with eukaryotic hosts where the host colonization dynamics may impact bacterial genome evolution. Host-associated lifestyle has been recognized as a barrier to HGT in parentally transmitted bacteria. However, pangenome evolution of environmentally acquired symbionts remains understudied, often due to limitations in symbiont cultivation. Using high-resolution metagenomics, here we study pangenome evolution of two co-occurring endosymbionts inhabiting Bathymodiolus brooksi mussels from a single cold seep. The symbionts, sulfur-oxidizing (SOX) and methane-oxidizing (MOX) gamma-proteobacteria, are environmentally acquired at an early developmental stage and individual mussels may harbor multiple strains of each symbiont species. We found differences in the accessory gene content of both symbionts across individual mussels, which are reflected by differences in symbiont strain composition. Compared with core genes, accessory genes are enriched in genome plasticity functions. We found no evidence for recent HGT between both symbionts. A comparison between the symbiont pangenomes revealed that the MOX population is less diverged and contains fewer accessory genes, supporting that the MOX association with B. brooksi is more recent in comparison to that of SOX. Our results show that the pangenomes of both symbionts evolved mainly by vertical inheritance. We conclude that genome evolution of environmentally transmitted symbionts that associate with individual hosts over their lifetime is affected by a narrow symbiosis where the frequency of HGT is constrained.

Phage-host coevolution in natural populations

Coevolution between bacteriophages (phages) and their bacterial hosts occurs through changes in resistance and counter-resistance mechanisms. To assess phage-host evolution in wild populations, we isolated 195 Vibrio crassostreae strains and 243 vibriophages during a 5-month time series from an oyster farm and combined these isolates with existing V. crassostreae and phage isolates. Cross-infection studies of 81,926 host-phage pairs delineated a modular network where phages are best at infecting co-occurring hosts, indicating local adaptation. Successful propagation of phage is restricted by the ability to adsorb to closely related bacteria and further constrained by strain-specific defence systems. These defences are highly diverse and predominantly located on mobile genetic elements, and multiple defences are active within a single genome. We further show that epigenetic and genomic modifications enable phage to adapt to bacterial defences and alter host range. Our findings reveal that the evolution of bacterial defences and phage counter-defences is underpinned by frequent genetic exchanges with, and between, mobile genetic elements.

Genomes of six viruses that infect Asgard archaea from deep-sea sediments

Asgard archaea are globally distributed prokaryotic microorganisms related to eukaryotes; however, viruses that infect these organisms have not been described. Here, using metagenome sequences recovered from deep-sea hydrothermal sediments, we characterize six relatively large (up to 117 kb) double-stranded DNA (dsDNA) viral genomes that infected two Asgard archaeal phyla, Lokiarchaeota and Helarchaeota. These viruses encode Caudovirales-like structural proteins, as well as proteins distinct from those described in known archaeal viruses. Their genomes contain around 1-5% of genes associated with eukaryotic nucleocytoplasmic large DNA viruses (NCLDVs) and appear to be capable of semi-autonomous genome replication, repair, epigenetic modifications and transcriptional regulation. Moreover, Helarchaeota viruses may hijack host ubiquitin systems similar to eukaryotic viruses. Genomic analysis of these Asgard viruses reveals that they contain features of both prokaryotic and eukaryotic viruses, and provides insights into their potential infection and host interaction mechanisms.

UG/Abi: a highly diverse family of prokaryotic reverse transcriptases associated with defense functions

Reverse transcriptases (RTs) are enzymes capable of synthesizing DNA using RNA as a template. Within the last few years, a burst of research has led to the discovery of novel prokaryotic RTs with diverse antiviral properties, such as DRTs (Defense-associated RTs), which belong to the so-called group of unknown RTs (UG) and are closely related to the Abortive Infection system (Abi) RTs. In this work, we performed a systematic analysis of UG and Abi RTs, increasing the number of UG/Abi members up to 42 highly diverse groups, most of which are predicted to be functionally associated with other gene(s) or domain(s). Based on this information, we classified these systems into three major classes. In addition, we reveal that most of these groups are associated with defense functions and/or mobile genetic elements, and demonstrate the antiphage role of four novel groups. Besides, we highlight the presence of one of these systems in novel families of human gut viruses infecting members of the Bacteroidetes and Firmicutes phyla. This work lays the foundation for a comprehensive and unified understanding of these highly diverse RTs with enormous biotechnological potential.

Computational Tools for the Analysis of Uncultivated Phage Genomes

Over a century of bacteriophage research has uncovered a plethora of fundamental aspects of their biology, ecology, and evolution. Furthermore, the introduction of community-level studies through metagenomics has revealed unprecedented insights on the impact that phages have on a range of ecological and physiological processes. It was not until the introduction of viral metagenomics that we began to grasp the astonishing breadth of genetic diversity encompassed by phage genomes. Novel phage genomes have been reported from a diverse range of biomes at an increasing rate, which has prompted the development of computational tools that support the multilevel characterization of these novel phages based solely on their genome sequences. The impact of these technologies has been so large that, together with MAGs (Metagenomic Assembled Genomes), we now have UViGs (Uncultivated Viral Genomes), which are now officially recognized by the International Committee for the Taxonomy of Viruses (ICTV), and new taxonomic groups can now be created based exclusively on genomic sequence information. Even though the available tools have immensely contributed to our knowledge of phage diversity and ecology, the ongoing surge in software programs makes it challenging to keep up with them and the purpose each one is designed for. Therefore, in this review, we describe a comprehensive set of currently available computational tools designed for the characterization of phage genome sequences, focusing on five specific analyses: (i) assembly and identification of phage and prophage sequences, (ii) phage genome annotation, (iii) phage taxonomic classification, (iv) phage-host interaction analysis, and (v) phage microdiversity.

Metagenomic assembled plasmids of the human microbiome vary across disease cohorts

We compiled a human metagenome assembled plasmid (MAP) database and interrogated differences across multiple studies that were originally designed to investigate the composition of the human microbiome across various lifestyles, life stages and events. This was performed as plasmids enable bacteria to rapidly expand their functional capacity through mobilisation, yet their contribution to human health and disease is poorly understood. We observed that inter-sample β-diversity differences of plasmid content (plasmidome) could distinguish cohorts across a multitude of conditions. We also show that reduced intra-sample plasmidome α-diversity is consistent amongst patients with inflammatory bowel disease (IBD) and Clostridioides difficile infections. We also show that faecal microbiota transplants can restore plasmidome diversity. Overall plasmidome diversity, specific plasmids, and plasmid-encoded functions can all potentially act as biomarkers of IBD or its severity. The human plasmidome is an overlooked facet of the microbiome and should be integrated into investigations regarding the role of the microbiome in promoting health or disease. Including MAP databases in analyses will enable a greater understanding of the roles of plasmid-encoded functions within the gut microbiome and will inform future human metagenome analyses.

Complete Genome Sequence of a Phage Infecting Sphingomonadaceae

We isolated a phage infecting a member of the Sphingomonadaceae family from a freshwater lake. The phage has a DNA genome of 41,771 bp, with a GC content of 61.7%. The genome harbors 50 predicted protein-coding genes and an auxiliary metabolic gene, which encodes a protein belonging to the radical S-adenosylmethionine superfamily.

RNN-VirSeeker: A Deep Learning Method for Identification of Short Viral Sequences From Metagenomes

Viruses are the most abundant biological entities on earth, and play vital roles in many aspects of microbial communities. As major human pathogens, viruses have caused huge mortality and morbidity to human society in history. Metagenomic sequencing methods could capture all microorganisms from microbiota, with sequences of viruses mixed with these of other species. Therefore, it is necessary to identify viral sequences from metagenomes. However, existing methods perform poorly on identifying short viral sequences. To solve this problem, a deep learning based method, RNN-VirSeeker, is proposed in this paper. RNN-VirSeeker was trained by sequences of 500bp sampled from known Virus and Host RefSeq genomes. Experimental results on the testing set have shown that RNN-VirSeeker exhibited AUROC of 0.9175, recall of 0.8640 and precision of 0.9211 for sequences of 500bp, and outperformed three widely used methods, VirSorter, VirFinder, and DeepVirFinder, on identifying short viral sequences. RNN-VirSeeker was also used to identify viral sequences from a CAMI dataset and a human gut metagenome. Compared with DeepVirFinder, RNN-VirSeeker identified more viral sequences from these metagenomes and achieved greater values of AUPRC and AUROC. RNN-VirSeeker is freely available at https://github.com/crazyinter/RNN-VirSeeker.

More Positive or More Negative? Metagenomic Analysis Reveals Roles of Virome in Human Disease-Related Gut Microbiome

Viruses are increasingly viewed as vital components of the human gut microbiota, while their roles in health and diseases remain incompletely understood. Here, we first sequenced and analyzed the 37 metagenomic and 18 host metabolomic samples related to irritable bowel syndrome (IBS) and found that some shifted viruses between IBS and controls covaried with shifted bacteria and metabolites. Especially, phages that infect beneficial lactic acid bacteria depleted in IBS covaried with their hosts. We also retrieved public whole-genome metagenomic datasets of another four diseases (type 2 diabetes, Crohn’s disease, colorectal cancer, and liver cirrhosis), totaling 438 samples including IBS, and performed uniform analysis of the gut viruses in diseases. By constructing disease-specific co-occurrence networks, we found viruses actively interacting with bacteria, negatively correlated with possible dysbiosis-related and inflammation-mediating bacteria, increasing the connectivity between bacteria modules, and contributing to the robustness of the networks. Functional enrichment analysis showed that phages interact with bacteria through predation or expressing genes involved in the transporter and secretion system, metabolic enzymes, etc. We further built a viral database to facilitate systematic functional classification and explored the functions of viral genes on interacting with bacteria. Our analyses provided a systematic view of the gut virome in the disease-related microbial community and suggested possible positive roles of viruses concerning gut health.

Phage-encoded ribosomal protein S21 expression is linked to late-stage phage replication

The ribosomal protein S21 (bS21) gene has been detected in diverse viruses with a large range of genome sizes, yet its in situ expression and potential significance have not been investigated. Here, we report five closely related clades of bacteriophages (phages) represented by 47 genomes (8 curated to completion and up to 331 kbp in length) that encode a bS21 gene. The bS21 gene is on the reverse strand within a conserved region that encodes the large terminase, major capsid protein, prohead protease, portal vertex proteins, and some hypothetical proteins. Based on CRISPR spacer targeting, the predominance of bacterial taxonomic affiliations of phage genes with those from Bacteroidetes, and the high sequence similarity of the phage bS21 genes and those from Bacteroidetes classes of Flavobacteriia, Cytophagia and Saprospiria, these phages are predicted to infect diverse Bacteroidetes species that inhabit a range of depths in freshwater lakes. Thus, bS21 phages have the potential to impact microbial community composition and carbon turnover in lake ecosystems. The transcriptionally active bS21-encoding phages were likely in the late stage of replication when collected, as core structural genes and bS21 were highly expressed. Thus, our analyses suggest that the phage bS21, which is involved in translation initiation, substitutes into the Bacteroidetes ribosomes and selects preferentially for phage transcripts during the late-stage replication when large-scale phage protein production is required for assembly of phage particles.

Isolation and Characterization of Novel Lytic Phages Infecting Multidrug-Resistant Escherichia coli

Urinary tract infections (UTIs) are the second most frequent bacterial infections worldwide, with Escherichia coli being the main causative agent. The increase of antibiotic-resistance determinants among isolates from clinical samples, including UTIs, makes the development of novel therapeutic strategies a necessity. In this context, the use of bacteriophages as a therapeutic alternative has been proposed, due to their ability to efficiently kill bacteria. In this work, we isolated and characterized three novel bacteriophages, microbes laboratory phage 1 (MLP1), MLP2, and MLP3, belonging to the Chaseviridae, Myoviridae, and Podoviridae families, respectively. These phages efficiently infect and kill laboratory reference strains and multidrug-resistant clinical E. coli isolates from patients with diagnosed UTIs. Interestingly, these phages are also able to infect intestinal pathogenic Escherichia coli strains, such as enteroaggregative E. coli and diffusely adherent E. coli. Our data show that the MLP phages recognize different regions of the lipopolysaccharide (LPS) molecule, an important virulence factor in bacteria that is also highly variable among different E. coli strains. Altogether, our results suggest that these phages may represent an interesting alternative for the treatment of antibiotic-resistant E. coli. IMPORTANCE Urinary tract infections affect approximately 150 million people annually. The current antibiotic resistance crisis demands the development of novel therapeutic alternatives. Our results show that three novel phages, MLP1, MLP2, and MLP3 are able to infect both laboratory and multidrug-resistant clinical isolates of Escherichia coli. Since these phages (i) efficiently kill antibiotic-resistant clinical isolates of uropathogenic Escherichia coli (UPEC), (ii) recognize different portions of the LPS molecule, and (iii) are able to efficiently infect intestinal pathogenic Escherichia coli hosts, we believe that these novel phages are good candidates to be used as a therapeutic alternative to treat antibiotic-resistant E. coli strains generating urinary tract and/or intestinal infections.

Transkingdom Analysis of the Female Reproductive Tract Reveals Bacteriophages form Communities

The female reproductive tract (FRT) microbiome plays a vital role in maintaining vaginal health. Viruses are key regulators of other microbial ecosystems, but little is known about how the FRT viruses (virome), particularly bacteriophages that comprise the phageome, impact FRT health and dysbiosis. We hypothesize that bacterial vaginosis (BV) is associated with altered FRT phageome diversity, transkingdom interplay, and bacteriophage discriminate taxa. Here, we conducted a retrospective, longitudinal analysis of vaginal swabs collected from 54 BV-positive and 46 BV-negative South African women. Bacteriome analysis revealed samples clustered into five distinct bacterial community groups (CGs), and further, bacterial alpha diversity was significantly associated with BV. Virome analysis on a subset of baseline samples showed FRT bacteriophages clustering into novel viral state types (VSTs), a viral community clustering system based on virome composition and abundance. Distinct BV bacteriophage signatures included increased alpha diversity along with discriminant Bacillus, Burkholderia, and Escherichia bacteriophages. Bacteriophage-bacteria transkingdom associations were also identified between Bacillus and Burkholderia viruses and BV-associated bacteria, providing key insights for future studies elucidating the transkingdom interactions driving BV-associated microbiome perturbations. In this cohort, bacteriophage-bacterial associations suggest complex interactions, which may play a role in the establishment and maintenance of BV.

In Silico Characterisation of Putative Prophages in Lactobacillaceae Used in Probiotics for Vaginal Health

While live biotherapeutics offer a promising approach to optimizing vaginal microbiota, the presence of functional prophages within introduced Lactobacillaceae strains could impact their safety and efficacy. We evaluated the presence of prophages in 895 publicly available Lactobacillaceae genomes using Phaster, Phigaro, Phispy, Prophet and Virsorter. Prophages were identified according to stringent (detected by ≥4 methods) or lenient criteria (detected by ≥2 methods), both with >80% reciprocal sequence overlap. The stringent approach identified 448 prophages within 359 genomes, with 40.1% genomes harbouring at least one prophage, while the lenient approach identified 1671 prophages within 83.7% of the genomes. To confirm our in silico estimates in vitro, we tested for inducible prophages in 57 vaginally-derived and commercial Lactobacillaceae isolates and found inducible prophages in 61.4% of the isolates. We characterised the in silico predicted prophages based on weighted gene repertoire relatedness and found that most belonged to the Siphoviridae or Myoviridae families. ResFam and eggNOG identified four potential antimicrobial resistance genes within the predicted prophages. Our results suggest that while Lactobacillaceae prophages seldomly carry clinically concerning genes and thus unlikely a pose a direct risk to human vaginal microbiomes, their high prevalence warrants the characterisation of Lactobacillaceae prophages in live biotherapeutics.

Genomic Diversity of Bacteriophages Infecting the Genus Acinetobacter

The number of sequenced Acinetobacter phage genomes in the International Nucleotide Sequence Database Collaboration has increased significantly in recent years, from 37 in 2017 to a total of 139 as of January 2021 with genome sizes ranging from 31 to 378 kb. Here, we explored the genetic diversity of the Acinetobacter phages using comparative genomics approaches that included assessment of nucleotide similarity, shared gene content, single gene phylogeny, and the network-based classification tool vConTACT2. Phages infecting Acinetobacter sp. are genetically diverse and can be grouped into 8 clusters (subfamilies) and 46 sub-clusters (genera), of which 8 represent genomic singletons (additional genera). We propose the creation of five new subfamilies and suggest a reorganisation of the genus Obolenskvirus. These results provide an updated view of the viruses infecting Acinetobacter species, providing insights into their diversity.

Resolving the structure of phage-bacteria interactions in the context of natural diversity

Microbial communities are shaped by viral predators. Yet, resolving which viruses (phages) and bacteria are interacting is a major challenge in the context of natural levels of microbial diversity. Thus, fundamental features of how phage-bacteria interactions are structured and evolve in the wild remain poorly resolved. Here we use large-scale isolation of environmental marine Vibrio bacteria and their phages to obtain estimates of strain-level phage predator loads, and use all-by-all host range assays to discover how phage and host genomic diversity shape interactions. We show that lytic interactions in environmental interaction networks (as observed in agar overlay) are sparse-with phage predator loads being low for most bacterial strains, and phages being host-strain-specific. Paradoxically, we also find that although overlap in killing is generally rare between tailed phages, recombination is common. Together, these results suggest that recombination during cryptic co-infections is an important mode of phage evolution in microbial communities. In the development of phages for bioengineering and therapeutics it is important to consider that nucleic acids of introduced phages may spread into local phage populations through recombination, and that the likelihood of transfer is not predictable based on lytic host range.

Hidden diversity of double-stranded DNA phages in symbiotic Rhizobium species

In this study, we addressed the extent of diversification of phages associated with nitrogen-fixing symbiotic Rhizobium species. Despite the ecological and economic importance of the Rhizobium genus, little is known about the diversity of the associated phages. A thorough assessment of viral diversity requires investigating both lytic phages and prophages harboured in diverse Rhizobium genomes. Protein-sharing networks identified 56 viral clusters (VCs) among a set of 425 isolated phages and predicted prophages. The VCs formed by phages had more proteins in common and a higher degree of synteny, and they group together in clades in the associated phylogenetic tree. By contrast, the VCs of prophages showed significant genetic variation and gene loss, with selective pressure on the remaining genes. Some VCs were found in various Rhizobium species and geographical locations, suggesting that they have wide host ranges. Our results indicate that the VCs represent distinct taxonomic units, probably representing taxa equivalent to genera or even species. The finding of previously undescribed phage taxa indicates the need for further exploration of the diversity of phages associated with Rhizobium species. This article is part of the theme issue 'The secret lives of microbial mobile genetic elements'.

To catch a hijacker: abundance, evolution and genetic diversity of P4-like bacteriophage satellites

Bacteriophages (phages) are bacterial parasites that can themselves be parasitized by phage satellites. The molecular mechanisms used by satellites to hijack phages are sometimes understood in great detail, but the origins, abundance, distribution and composition of these elements are poorly known. Here, we show that P4-like elements are present in more than 30% of the genomes of Enterobacterales, and in almost half of those of Escherichia coli, sometimes in multiple distinct copies. We identified over 1000 P4-like elements with very conserved genetic organization of the core genome and a few hotspots with highly variable genes. These elements are never found in plasmids and have very little homology to known phages, suggesting an independent evolutionary origin. Instead, they are scattered across chromosomes, possibly because their integrases are often exchanged with other elements. The rooted phylogenies of hijacking functions are correlated and suggest longstanding coevolution. They also reveal broad host ranges in P4-like elements, as almost identical elements can be found in distinct bacterial genera. Our results show that P4-like phage satellites constitute a very distinct, widespread and ancient family of mobile genetic elements. They pave the way for studying the molecular evolution of antagonistic interactions between phages and their satellites. This article is part of the theme issue 'The secret lives of microbial mobile genetic elements'.

Selective Isolation of Eggerthella lenta from Human Faeces and Characterisation of the Species Prophage Diversity

Eggerthella lenta is an anaerobic, high GC, Gram-positive bacillus commonly found in the human digestive tract that belongs to the class Coriobacteriia of the phylum Actinobacteria. This species has been of increasing interest as an important player in the metabolism of xenobiotics and dietary compounds. However, little is known regarding its susceptibility to bacteriophage predation and how this may influence its fitness. Here, we report the isolation of seven novel E. lenta strains using cefotaxime and ceftriaxone as selective agents. We conducted comparative and pangenome analyses of these strains and those publicly available to investigate the diversity of prophages associated with this species. Prophage gene products represent a minimum of 5.8% of the E. lenta pangenome, comprising at least ten distantly related prophage clades that display limited homology to currently known bacteriophages. All clades possess genes implicated in virion structure, lysis, lysogeny and, to a limited extent, DNA replication. Some prophages utilise tyrosine recombinases and diversity generating retroelements to generate phase variation among targeted genes. The prophages have differing levels of sensitivity to the CRISPR/cas systems of their hosts, with spacers from 44 E. lenta isolates found to target only five out of the ten identified prophage clades. Furthermore, using a PCR-based approach targeting the prophage attP site, we were able to determine that several of these elements can excise from the host chromosome, thus supporting the notion that these are active prophages. The findings of this study provide further insights into the diversity of prophages infecting species of the phylum Actinobacteria.

Characterization and in vitro testing of newly isolated lytic bacteriophages for the biocontrol of Pseudomonas aeruginosa

Aim: Two lytic phages were isolated using P. aeruginosa DSM19880 as host and fully characterized. Materials & methods: Phages were characterized physicochemically, biologically and genomically. Results & conclusion: Host range analysis revealed that the phages also infect some multidrug-resistant (MDR) P. aeruginosa clinical isolates. Increasing MOI from 1 to 1000 significantly increased phage efficiency and retarded bacteria regrowth, but phage ph0034 (reduction of 7.5 log CFU/ml) was more effective than phage ph0031 (reduction of 5.1 log CFU/ml) after 24 h. Both phages belong to Myoviridae family. Genome sequencing of phages ph0031 and ph0034 showed that they do not carry toxin, virulence, antibiotic resistance and integrase genes. The results obtained are highly relevant in the actual context of bacterial resistance to antibiotics.

OUP accepted manuscript

Motivation

Phage–host associations play important roles in microbial communities. But in natural communities, as opposed to culture-based lab studies where phages are discovered and characterized metagenomically, their hosts are generally not known. Several programs have been developed for predicting which phage infects which host based on various sequence similarity measures or machine learning approaches. These are often based on whole viral and host genomes, but in metagenomics-based studies, we rarely have whole genomes but rather must rely on contigs that are sometimes as short as hundreds of bp long. Therefore, we need programs that predict hosts of phage contigs on the basis of these short contigs. Although most existing programs can be applied to metagenomic datasets for these predictions, their accuracies are generally low. Here, we develop ContigNet, a convolutional neural network-based model capable of predicting phage–host matches based on relatively short contigs, and compare it to previously published VirHostMatcher (VHM) and WIsH.

Results

On the validation set, ContigNet achieves 72–85% area under the receiver operating characteristic curve (AUROC) scores, compared to the maximum of 68% by VHM or WIsH for contigs of lengths between 200 bps to 50 kbps. We also apply the model to the Metagenomic Gut Virus (MGV) catalogue, a dataset containing a wide range of draft genomes from metagenomic samples and achieve 60–70% AUROC scores compared to that of VHM and WIsH of 52%. Surprisingly, ContigNet can also be used to predict plasmid-host contig associations with high accuracy, indicating a similar genetic exchange between mobile genetic elements and their hosts.

Availability and implementation

The source code of ContigNet and related datasets can be downloaded from https://github.com/tianqitang1/ContigNet.