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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Giant Virus Global Proteomics Innovation: Comparative Evaluation of In-Gel and In-Solution Digestion Methods

With their unusually large genome and particle sizes, giant viruses (GVs) defy the conventional definition of viruses. Although most GVs isolated infect unicellular protozoans, such as amoeba, studies in the last decade have established their much wider prevalence infecting most eukaryotic supergroups and some giant viral families with the potential to be human pathogens. Their complexity, almost autonomous life cycle, and enigmatic evolution necessitate the study of GVs. The accurate assessment of GV proteome is a veritable challenge. We have compared the coverage of global protein identification using different methods for GVs isolated in Mumbai, Mimivirus Bombay (MVB), Powai Lake Megavirus (PLMV), and Kurlavirus (KV), along with two previously studied GVs, Acanthamoeba polyphaga Mimivirus (APMV) and Marseillevirus (MV). Our study shows that the simultaneous use of in-gel and in-solution digestion methods can significantly increase the coverage of protein identification in the global proteome analysis of purified GV particles. Combining the two methods of analyses, we identified an additional 72 proteins in APMV and 114 in MV compared with what have been previously reported. Similarly, proteomes of MVB, PLMV, and KV were analyzed, and a total of 242 proteins in MVB, 287 proteins in PLMV, and 174 proteins in KV were identified. Our results suggest that a combined methodology of in-gel and in-solution methods is more efficient and opens up new avenues for innovation in global proteome analysis of GVs. Future planetary health research on GVs can benefit from consideration of a broader range of proteomics methodologies as illustrated by the present study.

Dynamic coexistence driven by physiological transitions in microbial communities

Microbial ecosystems are commonly modeled by fixed interactions between species in steady exponential growth states. However, microbes often modify their environments so strongly that they are forced out of the exponential state into stressed or non-growing states. Such dynamics are typical of ecological succession in nature and serial-dilution cycles in the laboratory. Here, we introduce a phenomenological model, the Community State model, to gain insight into the dynamic coexistence of microbes due to changes in their physiological states. Our model bypasses specific interactions (e.g., nutrient starvation, stress, aggregation) that lead to different combinations of physiological states, referred to collectively as "community states", and modeled by specifying the growth preference of each species along a global ecological coordinate, taken here to be the total community biomass density. We identify three key features of such dynamical communities that contrast starkly with steady-state communities: increased tolerance of community diversity to fast growth rates of species dominating different community states, enhanced community stability through staggered dominance of different species in different community states, and increased requirement on growth dominance for the inclusion of late-growing species. These features, derived explicitly for simplified models, are proposed here to be principles aiding the understanding of complex dynamical communities. Our model shifts the focus of ecosystem dynamics from bottom-up studies based on idealized inter-species interaction to top-down studies based on accessible macroscopic observables such as growth rates and total biomass density, enabling quantitative examination of community-wide characteristics.

Comparative Analyses of Bacteriophage Genomes

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

Viral Prevalence and Genomic Xenology in the Coevolution of HzNV-2 (Nudiviridae) with Host Helicoverpa zea (Lepidoptera: Noctuidae)

Insect viruses have been described from numerous lineages, yet patterns of genetic exchange and viral prevalence, which are essential to understanding host-virus coevolution, are rarely studied. In Helicoverpa zea, the virus HzNV-2 can cause deformity of male and female genitalia, resulting in sterility. Using ddPCR, we found that male H. zea with malformed genitalia (agonadal) contained high levels of HzNV-2 DNA, confirming previous work. HzNV-2 was found to be prevalent throughout the United States, at more than twice the rate of the baculovirus HaSNPV, and that it contained several host-acquired DNA sequences. HzNV-2 possesses four recently endogenized lepidopteran genes and several more distantly related genes, including one gene with a bacteria-like sequence found in both host and virus. Among the recently acquired genes is cytosolic serine hydroxymethyltransferase (cSHMT). In nearly all tested H. zea, cSHMT contained a 200 bp transposable element (TE) that was not found in cSHMT of the sister species H. armigera. No other virus has been found with host cSHMT, and the study of this shared copy, including possible interactions, may yield new insights into the function of this gene with possible applications to insect biological control, and gene editing.

The role of bacteriophages in shaping bacterial composition and diversity in the human gut

The microbiota of the gut has continued to co-evolve alongside their human hosts conferring considerable health benefits including the production of nutrients, drug metabolism, modulation of the immune system, and playing an antagonistic role against pathogen invasion of the gastrointestinal tract (GIT). The gut is said to provide a habitat for diverse groups of microorganisms where they all co-habit and interact with one another and with the immune system of humans. Phages are bacterial parasites that require the host metabolic system to replicate via the lytic or lysogenic cycle. The phage and bacterial populations are regarded as the most dominant in the gut ecosystem. As such, among the various microbial interactions, the phage-bacteria interactions, although complex, have been demonstrated to co-evolve over time using different mechanisms such as predation, lysogenic conversion, and phage induction, alongside counterdefense by the bacterial population. With the help of models and dynamics of phage-bacteria interactions, the complexity behind their survival in the gut ecosystem was demystified, and their roles in maintaining gut homeostasis and promoting the overall health of humans were elucidated. Although the treatment of various gastrointestinal infections has been demonstrated to be successful against multidrug-resistant causative agents, concerns about this technique are still very much alive among researchers owing to the potential for phages to evolve. Since a dearth of knowledge exists regarding the use of phages for therapeutic purposes, more studies involving experimental models and clinical trials are needed to widen the understanding of bacteria-phage interactions and their association with immunological responses in the gut of humans.

VIRify: An integrated detection, annotation and taxonomic classification pipeline using virus-specific protein profile hidden Markov models

The study of viral communities has revealed the enormous diversity and impact these biological entities have on various ecosystems. These observations have sparked widespread interest in developing computational strategies that support the comprehensive characterisation of viral communities based on sequencing data. Here we introduce VIRify, a new computational pipeline designed to provide a user-friendly and accurate functional and taxonomic characterisation of viral communities. VIRify identifies viral contigs and prophages from metagenomic assemblies and annotates them using a collection of viral profile hidden Markov models (HMMs). These include our manually-curated profile HMMs, which serve as specific taxonomic markers for a wide range of prokaryotic and eukaryotic viral taxa and are thus used to reliably classify viral contigs. We tested VIRify on assemblies from two microbial mock communities, a large metagenomics study, and a collection of publicly available viral genomic sequences from the human gut. The results showed that VIRify could identify sequences from both prokaryotic and eukaryotic viruses, and provided taxonomic classifications from the genus to the family rank with an average accuracy of 86.6%. In addition, VIRify allowed the detection and taxonomic classification of a range of prokaryotic and eukaryotic viruses present in 243 marine metagenomic assemblies. Finally, the use of VIRify led to a large expansion in the number of taxonomically classified human gut viral sequences and the improvement of outdated and shallow taxonomic classifications. Overall, we demonstrate that VIRify is a novel and powerful resource that offers an enhanced capability to detect a broad range of viral contigs and taxonomically classify them.

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.

From Trees to Clouds: PhageClouds for Fast Comparison of ∼640,000 Phage Genomic Sequences and Host-Centric Visualization Using Genomic Network Graphs

Background: Fast and computationally efficient strategies are required to explore genomic relationships within an increasingly large and diverse phage sequence space. Here, we present PhageClouds, a novel approach using a graph database of phage genomic sequences and their intergenomic distances to explore the phage genomic sequence space. Methods: A total of 640,000 phage genomic sequences were retrieved from a variety of databases and public virome assemblies. Intergenomic distances were calculated with dashing, an alignment-free method suitable for handling massive data sets. These data were used to build a Neo4j^® graph database. Results: PhageClouds supported the search of related phages among all complete phage genomes from GenBank for a single query phage in just 10 s. Moreover, PhageClouds expanded the number of closely related phage sequences detected for both finished and draft phage genomes, in comparison with searches exclusively targeting phage entries from GenBank. Conclusions: PhageClouds is a novel resource that will facilitate the analysis of phage genomic sequences and the characterization of assembled phage genomes.

The evolution of coevolution in the study of species interactions

The study of reciprocal adaptation in interacting species has been an active and inspiring area of evolutionary research for nearly 60 years. Perhaps owing to its great natural history and potential consequences spanning population divergence to species diversification, coevolution continues to capture the imagination of biologists. Here we trace developments following Ehrlich and Raven's classic paper, with a particular focus on the modern influence of two studies by Dr. May Berenbaum in the 1980s. This series of classic work presented a compelling example exhibiting the macroevolutionary patterns predicted by Ehrlich and Raven and also formalized a microevolutionary approach to measuring selection, functional traits, and understanding reciprocal adaptation between plants and their herbivores. Following this breakthrough was a wave of research focusing on diversifying macroevolutionary patterns, mechanistic chemical ecology, and natural selection on populations within and across community types. Accordingly, we breakdown coevolutionary theory into specific hypotheses at different scales: reciprocal adaptation between populations within a community, differential coevolution among communities, lineage divergence, and phylogenetic patterns. We highlight progress as well as persistent gaps, especially the link between reciprocal adaptation and diversification.

Application of machine learning in bacteriophage research

Phages are one of the key components in the structure, dynamics, and interactions of microbial communities in different bins. It has a clear impact on human health and the food industry. Bacteriophage characterization using in vitro approaches are time/cost consuming and laborious tasks. On the other hand, with the advent of new high-throughput sequencing technology, the development of a powerful computational framework to characterize the newly identified bacteriophages is inevitable for future research. Machine learning includes powerful techniques that enable the analysis of complex datasets for knowledge discovery and pattern recognition. In this study, we have conducted a comprehensive review of machine learning methods application using different types of features were applied in various aspects of bacteriophage research including, automated curation, identification, classification, host species recognition, virion protein identification, and life cycle prediction. Moreover, potential limitations and advantages of the developed frameworks were discussed.

Bacteriophage-Mediated Control of Phytopathogenic Xanthomonads: A Promising Green Solution for the Future

Xanthomonads, members of the family Xanthomonadaceae, are economically important plant pathogenic bacteria responsible for infections of over 400 plant species. Bacteriophage-based biopesticides can provide an environmentally friendly, effective solution to control these bacteria. Bacteriophage-based biocontrol has important advantages over chemical pesticides, and treatment with these biopesticides is a minor intervention into the microflora. However, bacteriophages’ agricultural application has limitations rooted in these viruses’ biological properties as active substances. These disadvantageous features, together with the complicated registration process of bacteriophage-based biopesticides, means that there are few products available on the market. This review summarizes our knowledge of the Xanthomonas-host plant and bacteriophage-host bacterium interaction’s possible influence on bacteriophage-based biocontrol strategies and provides examples of greenhouse and field trials and products readily available in the EU and the USA. It also details the most important advantages and limitations of the agricultural application of bacteriophages. This paper also investigates the legal background and industrial property right issues of bacteriophage-based biopesticides. When appropriately applied, bacteriophages can provide a promising tool against xanthomonads, a possibility that is untapped. Information presented in this review aims to explore the potential of bacteriophage-based biopesticides in the control of xanthomonads in the future.

The impact of quorum sensing on the modulation of phage-host interactions

Bacteriophages are the most diverse and abundant biological entities on the Earth and require host bacteria to replicate. Because of this obligate relationship, in addition to the challenging conditions of surrounding environments, phages must integrate information about extrinsic and intrinsic factors when infecting their host. This integration helps to determine whether the infection becomes lytic or lysogenic, which likely influences phage spreading and long-term survival. Although a variety of environmental and physiological clues are known to modulate lysis-lysogeny decisions, the social interplay among phages and host populations has been overlooked until recently. A growing body of evidence indicates that cell-cell communication in bacteria and, more recently, peptide-based communication among phage-phage populations, affect phage-host interactions by controlling phage lysis-lysogeny decisions and phage counter-defensive strategies in bacteria. Here, we explore and discuss the role of signal molecules as well as quorum sensing and quenching factors that mediate phage-host interactions. Our aim is to provide an overview of population-dependent mechanisms that influence phage replication, and how social communication may affect the dynamics and evolution of microbial communities, including their implications in phage therapy.

Distribution of Prophages in the Oenococcus oeni Species

Oenococcus oeni is the most exploited lactic acid bacterium in the wine industry and drives the malolactic fermentation of wines. Although prophage-like sequences have been identified in the species, many are not characterized, and a global view of their integration and distribution amongst strains is currently lacking. In this work, we analyzed the complete genomes of 231 strains for the occurrence of prophages, and analyzed their size and positions of insertion. Our data show the limited variation in the number of prophages in O. oeni genomes, and that six sites of insertion within the bacterial genome are being used for site-specific recombination. Prophage diversity patterns varied significantly for different host lineages, and environmental niches. Overall, the findings highlight the pervasive presence of prophages in the O. oeni species, their role as a major source of within-species bacterial diversity and drivers of horizontal gene transfer. Our data also have implications for enhanced understanding of the prophage recombination events which occurred during evolution of O. oeni, as well as the potential of prophages in influencing the fitness of these bacteria in their distinct niches.

A Place for Viruses on the Tree of Life

Viruses are ubiquitous. They infect almost every species and are probably the most abundant biological entities on the planet, yet they are excluded from the Tree of Life (ToL). However, there can be no doubt that viruses play a significant role in evolution, the force that facilitates all life on Earth. Conceptually, viruses are regarded by many as non-living entities that hijack living cells in order to propagate. A strict separation between living and non-living entities places viruses far from the ToL, but this may be theoretically unsound. Advances in sequencing technology and comparative genomics have expanded our understanding of the evolutionary relationships between viruses and cellular organisms. Genomic and metagenomic data have revealed that co-evolution between viral and cellular genomes involves frequent horizontal gene transfer and the occasional co-option of novel functions over evolutionary time. From the giant, ameba-infecting marine viruses to the tiny Porcine circovirus harboring only two genes, viruses and their cellular hosts are ecologically and evolutionarily intertwined. When deciding how, if, and where viruses should be placed on the ToL, we should remember that the Tree functions best as a model of biological evolution on Earth, and it is important that models themselves evolve with our increasing understanding of biological systems.

Bacteriophages that infect Gram-negative bacteria as source of signal-arrest-release motif lysins

Treatment of infections caused by multidrug-resistant (MDR) Gram-negative bacteria is challenging, a potential solution for which is the use of bacteriophage-derived lytic enzymes. However, the exogenous action of bacteriophage lysins against Gram-negative bacteria is hindered due to the presence of an impermeable outer membrane in these bacteria. Nevertheless, recent research has demonstrated that some lysins are capable of permeating the outer membrane of Gram-negative bacteria with the help of signal peptides. In the present study, we investigated the genomes of 309 bacteriophages that infect Gram-negative pathogens of clinical interest in order to determine the evolutionary markers of signal peptide-containing lysins. Complete genomes displayed 265 putative lysins, of which 17 (6.41%) contained signal-arrest-release motifs and 41 (15.47%) contained cleavable signal peptides. There was no apparent relationship between host specificity and lysin diversity. Nevertheless, the evolution of lysin genes might not be independent of the rest of the bacteriophage genome once pan-genome clustering and lysin diversity appear to be correlated. In addition, signal peptide- and signal-arrest-release-containing lysins were monophyletically distributed in the protein cladogram, suggesting that the natural selection of holin-independent lysins is divergent. Our study screened 58 (21.89%) out of 265 potential candidates for in vitro experimentation against MDR bacteria.

Silent rain: does the atmosphere-mediated connectivity between microbiomes influence bacterial evolutionary rates?

Air carries a vast number of bacteria and viruses over great distances all the time. This leads to continuous introduction of foreign genetic material to local, established microbial communities. In this perspective, I ask whether this silent rain may have a slowing effect on the overall evolutionary rates in the microbial biosphere. Arguably, the greater the genetic divergence between gene 'donors' and 'recipients', the greater the chance that the gene product has a deleterious epistatic interaction with other gene products in its genetic environment. This is due to the long-term absence of check for mutual compatibility. As such, if an organism is extensively different from other bacteria, genetic innovations are less probable to fit to the genome. Here, genetic innovation would be anything that elevates the fitness of the gene vehicle (e.g. bacterium) over its contemporaries. Adopted innovations increase the fitness of the compatible genome over incompatible ones, thus possibly tempering the pace at which mutations accumulate in existing genomes over generations. I further discuss the transfer of bacteriophages through atmosphere and potential effects that this may have on local dynamics and perhaps phage survival.

The network structure and eco-evolutionary dynamics of CRISPR-induced immune diversification

As a heritable sequence-specific adaptive immune system, CRISPR-Cas is a powerful force shaping strain diversity in host-virus systems. While the diversity of CRISPR alleles has been explored, the associated structure and dynamics of host-virus interactions have not. We explore the role of CRISPR in mediating the interplay between host-virus interaction structure and eco-evolutionary dynamics in a computational model and compare the results with three empirical datasets from natural systems. We show that the structure of the networks describing who infects whom and the degree to which strains are immune, are respectively modular (containing groups of hosts and viruses that interact strongly) and weighted-nested (specialist hosts are more susceptible to subsets of viruses that in turn also infect the more generalist hosts with many spacers matching many viruses). The dynamic interplay between these networks influences transitions between dynamical regimes of virus diversification and host control. The three empirical systems exhibit weighted-nested immunity networks, a pattern our theory shows is indicative of hosts able to suppress virus diversification. Previously missing from studies of microbial host-pathogen systems, the immunity network plays a key role in the coevolutionary dynamics.

Two Lineages of Pseudomonas aeruginosa Filamentous Phages: Structural Uniformity over Integration Preferences

Pseudomonas aeruginosa filamentous (Pf) bacteriophages are important factors contributing to the pathogenicity of this opportunistic bacterium, including biofilm formation and suppression of bacterial phagocytosis by macrophages. In addition, the capacity of Pf phages to form liquid crystal structures and their high negative charge density makes them potent sequesters of cationic antibacterial agents, such as aminoglycoside antibiotics or host antimicrobial peptides. Therefore, Pf phages have been proposed as a potential biomarker for risk of antibiotic resistance development. The majority of studies describing biological functions of Pf viruses have been performed with only three of them: Pf1, Pf4, and Pf5. However, our analysis revealed that Pf phages exist as two evolutionary lineages (I and II), characterized by substantially different structural/morphogenesis properties, despite sharing the same integration sites in the host chromosomes. All aforementioned model Pf phages are members of the lineage I. Hence, it is reasonable to speculate that their interactions with P. aeruginosa and impact on its pathogenicity may be not completely extrapolated to the lineage II members. Furthermore, in order to organize the present numerical nomenclature of Pf phages, we propose a more informative approach based on the insertion sites, that is, Pf-tRNA-Gly, -Met, -Sec, -tmRNA, and -DR (direct repeats), which are fully compatible with one of five types of tyrosine integrases/recombinases XerC/D carried by these viruses. Finally, we discuss possible evolutionary mechanisms behind this division and consequences from the perspective of virus-virus, virus-bacterium, and virus-human interactions.

Identification of a novel archaea virus, detected in hydrocarbon polluted Hungarian and Canadian samples

Metagenomics is a helpful tool for the analysis of unculturable organisms and viruses. Viruses that target bacteria and archaea play important roles in the microbial diversity of various ecosystems. Here we show that Methanosarcina virus MV (MetMV), the second Methanosarcina sp. virus with a completely determined genome, is characteristic of hydrocarbon pollution in environmental (soil and water) samples. It was highly abundant in Hungarian hydrocarbon polluted samples and its genome was also present in the NCBI SRA database containing reads from hydrocarbon polluted samples collected in Canada, indicating the stability of its niche and the marker feature of this virus. MetMV, as the only currently identified marker virus for pollution in environmental samples, could contribute to the understanding of the complicated network of prokaryotes and their viruses driving the decomposition of environmental pollutants.

Impact of phages on soil bacterial communities and nitrogen availability under different assembly scenarios

BACKGROUND

Bacteriophages, the viruses infecting bacteria, are biological entities that can control their host populations. The ecological relevance of phages for microbial systems has been widely explored in aquatic environments, but the current understanding of the role of phages in terrestrial ecosystems remains limited. Here, our objective was to quantify the extent to which phages drive the assembly and functioning of soil bacterial communities. We performed a reciprocal transplant experiment using natural and sterilized soil incubated with different combinations of two soil microbial communities, challenged against native and non-native phage suspensions as well as against a cocktail of phage isolates. We tested three different community assembly scenarios by adding phages: (a) during soil colonization, (b) after colonization, and (c) in natural soil communities. One month after inoculation with phage suspensions, bacterial communities were assessed by 16S rRNA amplicon gene sequencing.

RESULTS

By comparing the treatments inoculated with active versus autoclaved phages, our results show that changes in phage pressure have the potential to impact soil bacterial community composition and diversity. We also found a positive effect of active phages on the soil ammonium concentration in a few treatments, which indicates that increased phage pressure may also be important for soil functions.

CONCLUSIONS

Overall, the present work contributes to expand the current knowledge about soil phages and provide some empirical evidence supporting their relevance for soil bacterial community assembly and functioning. Video Abstract.

Bacteriophages of the lower urinary tract

The discovery of bacteria in the female urinary bladder has fundamentally changed current dogma regarding the urinary tract and related urinary disorders. Previous research characterized many of the bacterial components of the female urinary tract, but the viral fraction of this community is largely unknown. Viruses within the human microbiota far outnumber bacterial cells, with the most abundant viruses being those that infect bacteria (bacteriophages). Similar to observations within the microbiota of the gut and oral cavity, preliminary surveys of the urinary tract and bladder microbiota indicate a rich diversity of uncharacterized bacteriophage (phage) species. Phages are vital members of the microbiota, having critical roles in shaping bacterial metabolism and community structure. Although phages have been discovered in the urinary tract, such as phages that infect Escherichia coli, sampling them is challenging owing to low biomass, possible contamination when using non-invasive methods and the invasiveness of methods that reduce the potential for contamination. Phages could influence bladder health, but an understanding of the association between phage communities, bacterial populations and bladder health is in its infancy. However, evidence suggests that phages can defend the host against pathogenic bacteria and, therefore, modulation of the microbiome using phages has therapeutic potential for lower urinary tract symptoms. Furthermore, as natural predators of bacteria, phages have garnered renewed interest for their use as antimicrobial agents, for instance, in the treatment of urinary tract infections.

Origin of a Core Bacterial Gene via Co-option and Detoxification of a Phage Lysin

Temperate phages constitute a potentially beneficial genetic reservoir for bacterial innovation despite being selfish entities encoding an infection cycle inherently at odds with bacterial fitness. These phages integrate their genomes into the bacterial host during infection, donating new but deleterious genetic material: the phage genome encodes toxic genes, such as lysins, that kill the bacterium during the phage infection cycle. Remarkably, some bacteria have exploited the destructive properties of phage genes for their own benefit by co-opting them as toxins for functions related to bacterial warfare, virulence, and secretion. However, do toxic phage genes ever become raw material for functional innovation? Here, we report on a toxic phage gene whose product has lost its toxicity and has become a domain of a core cellular factor, SpmX, throughout the bacterial order Caulobacterales. Using a combination of phylogenetics, bioinformatics, structural biology, cell biology, and biochemistry, we have investigated the origin and function of SpmX and determined that its occurrence is the result of the detoxification of a phage peptidoglycan hydrolase gene. We show that the retained, attenuated activity of the phage-derived domain plays an important role in proper cell morphology and developmental regulation in representatives of this large bacterial clade. To our knowledge, this is the first observation of a phage gene domestication event in which a toxic phage gene has been co-opted for core cellular function at the root of a large bacterial clade.

Bacterial 'Grounded' Prophages: Hotspots for Genetic Renovation and Innovation

Bacterial genomes are highly plastic allowing the generation of variants through mutations and acquisition of genetic information. The fittest variants are then selected by the econiche thereby allowing the bacterial adaptation and colonization of the habitat. Larger genomes, however, may impose metabolic burden and hence bacterial genomes are optimized by the loss of frivolous genetic information. The activity of temperate bacteriophages has acute consequences on the bacterial population as well as the bacterial genome through lytic and lysogenic cycles. Lysogeny is a selective advantage as the prophage provides immunity to the lysogen against secondary phage attack. Since the non-lysogens are eliminated by the lytic phages, lysogens multiply and colonize the habitat. Nevertheless, all lysogens have an imminent risk of lytic cycle activation and cell lysis. However, a mutation in the attachment sites or in the genes that encode the specific recombinase responsible for prophage excision could result in 'grounding' of the prophage. Since the lysogens with grounded prophage are immune to respective phage infection as well as dodge the induction of lytic cycle, we hypothesize that the selection of these mutant lysogens is favored relative to their normal lysogenic counterparts. These grounded prophages offer several advantages to the bacterial genome evolution through propensity for genetic variations including inversions, deletions, and insertions via horizontal gene transfer. We propose that the grounded prophages expedite bacterial genome evolution by acting as 'genetic buffer zones' thereby increasing the frequency as well as the diversity of variations on which natural selection favors the beneficial variants. The grounded prophages are also hotspots for horizontal gene transfer wherein several ecologically significant genes such as those involved in stress tolerance, antimicrobial resistance, and novel metabolic pathways, are integrated. Moreover, the high frequency of genetic changes within prophages also allows proportionate probability for the de novo genesis of genetic information. Through sequence analyses of well-characterized E. coli prophages we exemplify various roles of grounded prophages in E. coli ecology and evolution. Therefore, the temperate prophages are one of the most significant drivers of bacterial genome evolution and sites of biogenesis of genetic information.

MARVEL, a Tool for Prediction of Bacteriophage Sequences in Metagenomic Bins

Here we present MARVEL, a tool for prediction of double-stranded DNA bacteriophage sequences in metagenomic bins. MARVEL uses a random forest machine learning approach. We trained the program on a dataset with 1,247 phage and 1,029 bacterial genomes, and tested it on a dataset with 335 bacterial and 177 phage genomes. We show that three simple genomic features extracted from contig sequences were sufficient to achieve a good performance in separating bacterial from phage sequences: gene density, strand shifts, and fraction of significant hits to a viral protein database. We compared the performance of MARVEL to that of VirSorter and VirFinder, two popular programs for predicting viral sequences. Our results show that all three programs have comparable specificity, but MARVEL achieves much better performance on the recall (sensitivity) measure. This means that MARVEL should be able to identify many more phage sequences in metagenomic bins than heretofore has been possible. In a simple test with real data, containing mostly bacterial sequences, MARVEL classified 58 out of 209 bins as phage genomes; other evidence suggests that 57 of these 58 bins are novel phage sequences. MARVEL is freely available at https://github.com/LaboratorioBioinformatica/MARVEL.

MARVEL, a Tool for Prediction of Bacteriophage Sequences in Metagenomic Bins

Here we present MARVEL, a tool for prediction of double-stranded DNA bacteriophage sequences in metagenomic bins. MARVEL uses a random forest machine learning approach. We trained the program on a dataset with 1,247 phage and 1,029 bacterial genomes, and tested it on a dataset with 335 bacterial and 177 phage genomes. We show that three simple genomic features extracted from contig sequences were sufficient to achieve a good performance in separating bacterial from phage sequences: gene density, strand shifts, and fraction of significant hits to a viral protein database. We compared the performance of MARVEL to that of VirSorter and VirFinder, two popular programs for predicting viral sequences. Our results show that all three programs have comparable specificity, but MARVEL achieves much better performance on the recall (sensitivity) measure. This means that MARVEL should be able to identify many more phage sequences in metagenomic bins than heretofore has been possible. In a simple test with real data, containing mostly bacterial sequences, MARVEL classified 58 out of 209 bins as phage genomes; other evidence suggests that 57 of these 58 bins are novel phage sequences. MARVEL is freely available at https://github.com/LaboratorioBioinformatica/MARVEL.