Evolutionarily conserved orthologous families in phages are relatively rare in their prokaryotic hosts

A novel genus of Pectobacterium bacteriophages display broad host range by targeting several species of Danish soft rot isolates

The bacterial diseases black leg and soft rot in potatoes cause heavy losses of potatoes worldwide. Bacteria within the genus Pectobacteriaceae are the causative agents of black leg and soft rot. The use of antibiotics in agriculture is heavily regulated and no other effective treatment currently exists, but bacteriophages (phages) have shown promise as potential biocontrol agents. In this study we isolated soft rot bacteria from potato tubers and plant tissue displaying soft rot or black leg symptoms collected in Danish fields. We then used the isolated bacterial strains as hosts for phage isolation. Using organic waste, we isolated phages targeting different species within Pectobacterium. Here we focus on seven of these phages representing a new genus primarily targeting P. brasiliense; phage Ymer, Amona, Sabo, Abuela, Koroua, Taid and Pappous. TEM image of phage Ymer showed siphovirus morphotype, and the proposed Ymer genus belongs to the class Caudoviricetes, with double-stranded DNA genomes varying from 39 kb to 43 kb. In silico host range prediction using a CRISPR-Cas spacer database suggested both P. brasiliense, P. polaris and P. versatile as natural hosts for phages within the proposed Ymer genus. A following host range experiment, using 47 bacterial isolates from Danish tubers and plants symptomatic with soft rot or black leg disease verified the in silico host range prediction, as the genus as a group were able to infect all three Pectobacterium species. Phages did, however, primarily target P. brasiliense isolates and displayed differences in host range even within the species level. Two of the phages were able to infect two or more Pectobacterium species. Despite no nucleotide similarity with any phages in the NCBI database, the proposed Ymer genus did share some similarity at the protein level, as well as gene synteny, with currently known phages. None of the phages encoded integrases or other genes typically associated with lysogeny. Similarly, no virulence factors nor antimicrobial resistance genes were found, and combined with their ability to infect several soft rot-causing Pectobacterium species from Danish fields, demonstrates their potential as biocontrol agents against soft rot and black leg diseases in potatoes.

Use of an Integrated Approach Involving AlphaFold Predictions for the Evolutionary Taxonomy of Duplodnaviria Viruses

The evaluation of the evolutionary relationships is exceptionally important for the taxonomy of viruses, which is a rapidly expanding area of research. The classification of viral groups belonging to the realm Duplodnaviria, which include tailed bacteriophages, head-tailed archaeal viruses and herpesviruses, has undergone many changes in recent years and continues to improve. One of the challenging tasks of Duplodnaviria taxonomy is the classification of high-ranked taxa, including families and orders. At the moment, only 17 of 50 families have been assigned to orders. The evaluation of the evolutionary relationships between viruses is complicated by the high level of divergence of viral proteins. However, the development of structure prediction algorithms, including the award-winning AlphaFold, encourages the use of the results of structural predictions to clarify the evolutionary history of viral proteins. In this study, the evolutionary relationships of two conserved viral proteins, the major capsid protein and terminase, representing different viruses, including all classified Duplodnaviria families, have been analysed using AlphaFold modelling. This analysis has been undertaken using structural comparisons and different phylogenetic methods. The results of the analyses mainly indicated the high quality of AlphaFold modelling and the possibility of using the AlphaFold predictions, together with other methods, for the reconstruction of the evolutionary relationships between distant viral groups. Based on the results of this integrated approach, assumptions have been made about refining the taxonomic classification of bacterial and archaeal Duplodnaviria groups, and problems relating to the taxonomic classification of Duplodnaviria have been discussed.

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.

Potential Solutions Using Bacteriophages against Antimicrobial Resistant Bacteria

Bacteriophages are viruses that specifically infect a bacterial host. They play a great role in the modern biotechnology and antibiotic-resistant microbe era. Since the discovery of phages, their application as a control agent has faced challenges that made antibiotics a better fit for combating pathogenic bacteria. Recently, with the novel sequencing technologies providing new insight into the nature of bacteriophages, their application has a second chance to be used. However, novel challenges need to be addressed to provide proper strategies for their practical application. This review focuses on addressing these challenges by initially introducing the nature of bacteriophages and describing the phage-host-dependent strategies for phage application. We also describe the effect of the long-term application of phages in natural environments and other bacterial communities. Overall, this review gathered crucial information for the future application of phages. We predict the use of phages will not be the only control strategy against pathogenic bacteria. Therefore, more studies must be done for low-risk control methods against antimicrobial-resistant bacteria.

Informative Regions In Viral Genomes

Viruses, far from being just parasites affecting hosts' fitness, are major players in any microbial ecosystem. In spite of their broad abundance, viruses, in particular bacteriophages, remain largely unknown since only about 20% of sequences obtained from viral community DNA surveys could be annotated by comparison with public databases. In order to shed some light into this genetic dark matter we expanded the search of orthologous groups as potential markers to viral taxonomy from bacteriophages and included eukaryotic viruses, establishing a set of 31,150 ViPhOGs (Eukaryotic Viruses and Phages Orthologous Groups). To do this, we examine the non-redundant viral diversity stored in public databases, predict proteins in genomes lacking such information, and used all annotated and predicted proteins to identify potential protein domains. The clustering of domains and unannotated regions into orthologous groups was done using cogSoft. Finally, we employed a random forest implementation to classify genomes into their taxonomy and found that the presence or absence of ViPhOGs is significantly associated with their taxonomy. Furthermore, we established a set of 1457 ViPhOGs that given their importance for the classification could be considered as markers or signatures for the different taxonomic groups defined by the ICTV at the order, family, and genus levels.

The LUCA and its complex virome

The last universal cellular ancestor (LUCA) is the most recent population of organisms from which all cellular life on Earth descends. The reconstruction of the genome and phenotype of the LUCA is a major challenge in evolutionary biology. Given that all life forms are associated with viruses and/or other mobile genetic elements, there is no doubt that the LUCA was a host to viruses. Here, by projecting back in time using the extant distribution of viruses across the two primary domains of life, bacteria and archaea, and tracing the evolutionary histories of some key virus genes, we attempt a reconstruction of the LUCA virome. Even a conservative version of this reconstruction suggests a remarkably complex virome that already included the main groups of extant viruses of bacteria and archaea. We further present evidence of extensive virus evolution antedating the LUCA. The presence of a highly complex virome implies the substantial genomic and pan-genomic complexity of the LUCA itself.

Long-read metagenomic exploration of extrachromosomal mobile genetic elements in the human gut

BACKGROUND

Elucidating the ecological and biological identity of extrachromosomal mobile genetic elements (eMGEs), such as plasmids and bacteriophages, in the human gut remains challenging due to their high complexity and diversity.

RESULTS

Here, we show efficient identification of eMGEs as complete circular or linear contigs from PacBio long-read metagenomic data. De novo assembly of PacBio long reads from 12 faecal samples generated 82 eMGE contigs (2.5~666.7-kb), which were classified as 71 plasmids and 11 bacteriophages, including 58 novel plasmids and six bacteriophages, and complete genomes of five diverse crAssphages with terminal direct repeats. In a dataset of 413 gut metagenomes from five countries, many of the identified plasmids were highly abundant and prevalent. The ratio of gut plasmids by our plasmid data is more than twice that in the public database. Plasmids outnumbered bacterial chromosomes three to one on average in this metagenomic dataset. Host prediction suggested that Bacteroidetes-associated plasmids predominated, regardless of microbial abundance. The analysis found several plasmid-enriched functions, such as inorganic ion transport, while antibiotic resistance genes were harboured mostly in low-abundance Proteobacteria-associated plasmids.

CONCLUSIONS

Overall, long-read metagenomics provided an efficient approach for unravelling the complete structure of human gut eMGEs, particularly plasmids.

Defining a Core Genome for the Herpesvirales and Exploring their Evolutionary Relationship with the Caudovirales

The order Herpesvirales encompasses a wide variety of important and broadly distributed human pathogens. During the last decades, similarities in the viral cycle and the structure of some of their proteins with those of the order Caudovirales, the tailed bacterial viruses, have brought speculation regarding the existence of an evolutionary relationship between these clades. To evaluate such hypothesis, we used over 600 Herpesvirales and 2000 Caudovirales complete genomes to search for the presence or absence of clusters of orthologous protein domains and constructed a dendrogram based on their compositional similarities. The results obtained strongly suggest an evolutionary relationship between the two orders. Furthermore, they allowed to propose a core genome for the Herpesvirales, composed of 4 proteins, including the ATPase subunit of the DNA-packaging terminase, the only protein with previously verified conservation. Accordingly, a phylogenetic tree constructed with sequences derived from the clusters associated to these proteins grouped the Herpesvirales strains accordingly to the established families and subfamilies. Overall, this work provides results supporting the hypothesis that the two orders are evolutionarily related and contributes to the understanding of the history of the Herpesvirales.

Expanding the Diversity of Myoviridae Phages Infecting Lactobacillus plantarum-A Novel Lineage of Lactobacillus Phages Comprising Five New Members

Lactobacillus plantarum is a bacterium with probiotic properties and promising applications in the food industry and agriculture. So far, bacteriophages of this bacterium have been moderately addressed. We examined the diversity of five new L. plantarum phages via whole genome shotgun sequencing and in silico protein predictions. Moreover, we looked into their phylogeny and their potential genomic similarities to other complete phage genome records through extensive nucleotide and protein comparisons. These analyses revealed a high degree of similarity among the five phages, which extended to the vast majority of predicted virion-associated proteins. Based on these, we selected one of the phages as a representative and performed transmission electron microscopy and structural protein sequencing tests. Overall, the results suggested that the five phages belong to the family Myoviridae, they have a long genome of 137,973-141,344 bp, a G/C content of 36.3-36.6% that is quite distinct from their host's, and surprisingly, 7 to 15 tRNAs. Only an average 41/174 of their predicted genes were assigned a function. The comparative analyses unraveled considerable genetic diversity for the five L. plantarum phages in this study. Hence, the new genus "Semelevirus" was proposed, comprising exclusively of the five phages. This novel lineage of Lactobacillus phages provides further insight into the genetic heterogeneity of phages infecting Lactobacillus sp. The five new Lactobacillus phages have potential value for the development of more robust starters through, for example, the selection of mutants insensitive to phage infections. The five phages could also form part of phage cocktails, which producers would apply in different stages of L. plantarum fermentations in order to create a range of organoleptic outputs.

Phage Therapy in the Twenty-First Century: Facing the Decline of the Antibiotic Era; Is It Finally Time for the Age of the Phage?

Burgeoning problems of antimicrobial resistance dictate that new solutions be developed to combat old foes. Use of lytic bacteriophages (phages) for the treatment of drug-resistant bacterial infections is one approach that has gained significant traction in recent years. Fueled by reports of experimental phage therapy cases with very positive patient outcomes, several early-stage clinical trials of therapeutic phage products have been launched in the United States. Eventual licensure enabling widespread access to phages is the goal; however, new paths to regulatory approval and mass-market distribution, distinct from those of small-molecule antibiotics, must be forged first. This review highlights unique aspects related to the clinical use of phages, including advantages to be reaped as well as challenges to be overcome.

Evolution of the Large Nucleocytoplasmic DNA Viruses of Eukaryotes and Convergent Origins of Viral Gigantism

The Nucleocytoplasmic Large DNA Viruses (NCLDV) of eukaryotes (proposed order "Megavirales") comprise an expansive group of eukaryotic viruses that consists of the families Poxviridae, Asfarviridae, Iridoviridae, Ascoviridae, Phycodnaviridae, Marseilleviridae, Pithoviridae, and Mimiviridae, as well as Pandoraviruses, Molliviruses, and Faustoviruses that so far remain unaccounted by the official virus taxonomy. All these viruses have double-stranded DNA genomes that range in size from about 100 kilobases (kb) to more than 2.5 megabases. The viruses with genomes larger than 500kb are informally considered "giant," and the largest giant viruses surpass numerous bacteria and archaea in both particle and genome size. The discovery of giant viruses has been highly unexpected and has changed the perception of viral size and complexity, and even, arguably, the entire concept of a virus. Given that giant viruses encode multiple proteins that are universal among cellular life forms and are components of the translation system, the quintessential cellular molecular machinery, attempts have been made to incorporate these viruses in the evolutionary tree of cellular life. Moreover, evolutionary scenarios of the origin of giant viruses from a fourth, supposedly extinct domain of cellular life have been proposed. However, despite all the differences in the genome size and gene repertoire, the NCLDV can be confidently defined as monophyletic group, on the strength of the presence of about 40 genes that can be traced back to their last common ancestor. Using several most strongly conserved genes from this ancestral set, a well-resolved phylogenetic tree of the NCLDV was built and employed as the scaffold to reconstruct the history of gene gain and loss throughout the course of the evolution of this group of viruses. This reconstruction reveals extremely dynamic evolution that involved extensive gene gain and loss in many groups of viruses and indicates that giant viruses emerged independently in several clades of the NCLDV. Thus, these giants of the virus world evolved repeatedly from smaller and simpler viruses, rather than from a fourth domain of cellular life, and captured numerous genes, including those for translation system components, from eukaryotes, along with some bacterial genes. Even deeper evolutionary reconstructions reveal apparent links between the NCLDV and smaller viruses of eukaryotes, such as adenoviruses, and ultimately, derive all these viruses from tailless bacteriophages.

A comprehensive and quantitative exploration of thousands of viral genomes

The complete assembly of viral genomes from metagenomic datasets (short genomic sequences gathered from environmental samples) has proven to be challenging, so there are significant blind spots when we view viral genomes through the lens of metagenomics. One approach to overcoming this problem is to leverage the thousands of complete viral genomes that are publicly available. Here we describe our efforts to assemble a comprehensive resource that provides a quantitative snapshot of viral genomic trends – such as gene density, noncoding percentage, and abundances of functional gene categories – across thousands of viral genomes. We have also developed a coarse-grained method for visualizing viral genome organization for hundreds of genomes at once, and have explored the extent of the overlap between bacterial and bacteriophage gene pools. Existing viral classification systems were developed prior to the sequencing era, so we present our analysis in a way that allows us to assess the utility of the different classification systems for capturing genomic trends.

Sequence analysis of malacoherpesvirus proteins: Pan-herpesvirus capsid module and replication enzymes with an ancient connection to "Megavirales"

The order Herpesvirales includes animal viruses with large double-strand DNA genomes replicating in the nucleus. The main capsid protein in the best-studied family Herpesviridae contains a domain with HK97-like fold related to bacteriophage head proteins, and several virion maturation factors are also homologous between phages and herpesviruses. The origin of herpesvirus DNA replication proteins is less well understood. While analyzing the genomes of herpesviruses in the family Malacohepresviridae, we identified nearly 30 families of proteins conserved in other herpesviruses, including several phage-related domains in morphogenetic proteins. Herpesvirus DNA replication factors have complex evolutionary history: some are related to cellular proteins, but others are closer to homologs from large nucleocytoplasmic DNA viruses. Phylogenetic analyses suggest that the core replication machinery of herpesviruses may have been recruited from the same pool as in the case of other large DNA viruses of eukaryotes.

To Be or Not To Be T4: Evidence of a Complex Evolutionary Pathway of Head Structure and Assembly in Giant Salmonella Virus SPN3US

Giant Salmonella phage SPN3US has a 240-kb dsDNA genome and a large complex virion composed of many proteins for which the functions of most are undefined. We recently determined that SPN3US shares a core set of genes with related giant phages and sequenced and characterized 18 amber mutants to facilitate its use as a genetic model system. Notably, SPN3US and related giant phages contain a bolus of ejection proteins within their heads, including a multi-subunit virion RNA polymerase (vRNAP), that enter the host cell with the DNA during infection. In this study, we characterized the SPN3US virion using mass spectrometry to gain insight into its head composition and the features that its head shares with those of related giant phages and with T4 phage. SPN3US has only homologs to the T4 proteins critical for prohead shell formation, the portal and major capsid proteins, as well as to the major enzymes essential for head maturation, the prohead protease and large terminase subunit. Eight of ~50 SPN3US head proteins were found to undergo proteolytic processing at a cleavage motif by the prohead protease gp245. Gp245 undergoes auto-cleavage of its C-terminus, suggesting this is a conserved activation and/or maturation feature of related phage proteases. Analyses of essential head gene mutants showed that the five subunits of the vRNAP must be assembled for any subunit to be incorporated into the prohead, although the assembled vRNAP must then undergo subsequent major conformational rearrangements in the DNA packed capsid to allow ejection through the ~30 Å diameter tail tube for transcription from the injected DNA. In addition, ejection protein candidate gp243 was found to play a critical role in head assembly. Our analyses of the vRNAP and gp243 mutants highlighted an unexpected dichotomy in giant phage head maturation: while all analyzed giant phages have a homologous protease that processes major capsid and portal proteins, processing of ejection proteins is not always a stable/defining feature. Our identification in SPN3US, and related phages, of a diverged paralog to the prohead protease further hints toward a complicated evolutionary pathway for giant phage head structure and assembly.

Genome analysis of three novel lytic Vibrio coralliilyticus phages isolated from seawater, Okinawa, Japan

Three novel Vibrio phages were isolated from seawater in Okinawa. The Vibrio phage RYC infected Vibrio coralliilyticus SWA 07, while Vibrio phages CKB-S1 and CKB-S2 infected the coral pathogen V. coralliilyticus P1 (LMG 23696). The Vibrio phages CKB-S1 and CKB-S2 displayed head-tail structures whereas the Vibrio phage RYC showed a tailless non-enveloped capsid. All these Vibrio phages contained linear and double-stranded DNA. The whole genome sequencing revealed that Vibrio phage RYC has a larger genome size compared to Vibrio phages CKB-S1 and CKB-S2, and six tRNAs genes were found only in Vibrio phage RYC. Genome-wide comparison showed that Vibrio phage CKB-S1 was closely related, but was not identical, to Vibrio parahaemolyticus phages VP16T and VP16C. Meanwhile, the Vibrio phages RYC and CKB-S2 did not show high genome-wide similarity to any phages. These results suggest that the Vibrio phages CKB-S1, CKB-S2 and RYC are novel phages, which need further exploration, especially for their potential applications in phage therapy.

Jumbo Bacteriophages: An Overview

Tailed bacteriophages with genomes larger than 200 kbp are classified as Jumbo phages, and are rarely isolated by conventional methods. These phages are designated “jumbo” owing to their most notable features of a large phage virion and large genome size. However, in addition to these, jumbo phages also exhibit several novel characteristics that have not been observed for phages with smaller genomes, which differentiate jumbo phages in terms of genome organization, virion structure, progeny propagation, and evolution. In this review, we summarize available reports on jumbo phages and discuss the differences between jumbo phages and small-genome phages. We also discuss data suggesting that jumbo phages might have evolved from phages with smaller genomes by acquiring additional functional genes, and that these additional genes reduce the dependence of the jumbo phages on the host bacteria.

Scaffolding bacterial genomes and probing host-virus interactions in gut microbiome by proximity ligation (chromosome capture) assay

The biochemical activities of microbial communities, or microbiomes, are essential parts of environmental and animal ecosystems. The dynamics, balance, and effects of these communities are strongly influenced by phages present in the population. Being able to characterize bacterium-phage relationships is therefore essential to investigate these ecosystems to the full extent of their complexity. However, this task is currently limited by (i) the ability to characterize complete bacterial and viral genomes from a complex mix of species and (ii) the difficulty to assign phage sequences to their bacterial hosts. We show that both limitations can be circumvented using meta3C, an experimental and computational approach that exploits the physical contacts between DNA molecules to infer their proximity. In a single experiment, dozens of bacterial and phage genomes present in a complex mouse gut microbiota were assembled and scaffolded de novo. The phage genomes were then assigned to their putative bacterial hosts according to the physical contacts between the different DNA molecules, opening new perspectives for a comprehensive picture of the genomic structure of the gut flora. Therefore, this work holds far-reaching implications for human health studies aiming to bridge the virome to the microbiome.

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

Viruses are the most abundant and diverse biological entities on earth, and while most of this diversity remains completely unexplored, advances in genome sequencing have provided unprecedented glimpses into the virosphere. The Prokaryotic Virus Orthologous Groups (pVOGs, formerly called Phage Orthologous Groups, POGs) resource has aided in this task over the past decade by using automated methods to keep pace with the rapid increase in genomic data. The uses of pVOGs include functional annotation of viral proteins, identification of genes and viruses in uncharacterized DNA samples, phylogenetic analysis, large-scale comparative genomics projects, and more. The pVOGs database represents a comprehensive set of orthologous gene families shared across multiple complete genomes of viruses that infect bacterial or archaeal hosts (viruses of eukaryotes will be added at a future date). The pVOGs are constructed within the Clusters of Orthologous Groups (COGs) framework that is widely used for orthology identification in prokaryotes. Since the previous release of the POGs, the size has tripled to nearly 3000 genomes and 300 000 proteins, and the number of conserved orthologous groups doubled to 9518. User-friendly webpages are available, including multiple sequence alignments and HMM profiles for each VOG. These changes provide major improvements to the pVOGs database, at a time of rapid advances in virus genomics. The pVOGs database is hosted jointly at the University of Iowa at http://dmk-brain.ecn.uiowa.edu/pVOGs and the NCBI at ftp://ftp.ncbi.nlm.nih.gov/pub/kristensen/pVOGs/home.html.

The Double-Stranded DNA Virosphere as a Modular Hierarchical Network of Gene Sharing

Virus genomes are prone to extensive gene loss, gain, and exchange and share no universal genes. Therefore, in a broad-scale study of virus evolution, gene and genome network analyses can complement traditional phylogenetics. We performed an exhaustive comparative analysis of the genomes of double-stranded DNA (dsDNA) viruses by using the bipartite network approach and found a robust hierarchical modularity in the dsDNA virosphere. Bipartite networks consist of two classes of nodes, with nodes in one class, in this case genomes, being connected via nodes of the second class, in this case genes. Such a network can be partitioned into modules that combine nodes from both classes. The bipartite network of dsDNA viruses includes 19 modules that form 5 major and 3 minor supermodules. Of these modules, 11 include tailed bacteriophages, reflecting the diversity of this largest group of viruses. The module analysis quantitatively validates and refines previously proposed nontrivial evolutionary relationships. An expansive supermodule combines the large and giant viruses of the putative order "Megavirales" with diverse moderate-sized viruses and related mobile elements. All viruses in this supermodule share a distinct morphogenetic tool kit with a double jelly roll major capsid protein. Herpesviruses and tailed bacteriophages comprise another supermodule, held together by a distinct set of morphogenetic proteins centered on the HK97-like major capsid protein. Together, these two supermodules cover the great majority of currently known dsDNA viruses. We formally identify a set of 14 viral hallmark genes that comprise the hubs of the network and account for most of the intermodule connections.

IMPORTANCE

Viruses and related mobile genetic elements are the dominant biological entities on earth, but their evolution is not sufficiently understood and their classification is not adequately developed. The key reason is the characteristic high rate of virus evolution that involves not only sequence change but also extensive gene loss, gain, and exchange. Therefore, in the study of virus evolution on a large scale, traditional phylogenetic approaches have limited applicability and have to be complemented by gene and genome network analyses. We applied state-of-the art methods of such analysis to reveal robust hierarchical modularity in the genomes of double-stranded DNA viruses. Some of the identified modules combine highly diverse viruses infecting bacteria, archaea, and eukaryotes, in support of previous hypotheses on direct evolutionary relationships between viruses from the three domains of cellular life. We formally identify a set of 14 viral hallmark genes that hold together the genomic network.

Genomes of Abundant and Widespread Viruses from the Deep Ocean

The deep sea is a massive, largely oligotrophic ecosystem, stretched over nearly 65% of the planet’s surface. Deep-sea planktonic communities are almost completely dependent upon organic carbon sinking from the productive surface, forming a vital component of global biogeochemical cycles. However, despite their importance, viruses from the deep ocean remain largely unknown. Here, we describe the first complete genomes of deep-sea viruses assembled from metagenomic fosmid libraries. “Candidatus Pelagibacter” (SAR11) phage HTVC010P and Puniceispirillum phage HMO-2011 are considered the most abundant cultured marine viruses known to date. Remarkably, some of the viruses described here recruited as many reads from deep waters as these viruses do in the photic zone, and, considering the gigantic scale of the bathypelagic habitat, these genomes provide information about what could be some of the most abundant viruses in the world at large. Their role in the viral shunt in the global ocean could be very significant. Despite the challenges encountered in inferring the identity of their hosts, we identified one virus predicted to infect members of the globally distributed SAR11 cluster. We also identified a number of putative proviruses from diverse taxa, including deltaproteobacteria, bacteroidetes, SAR11, and gammaproteobacteria. Moreover, our findings also indicate that lysogeny is the preferred mode of existence for deep-sea viruses inhabiting an energy-limited environment, in sharp contrast to the predominantly lytic lifestyle of their photic-zone counterparts. Some of the viruses show a widespread distribution, supporting the tenet “everything is everywhere” for the deep-ocean virome.

The deep sea is among the largest known habitats and a critical cog in biogeochemical cycling but remains underexplored in its microbiology. Even more than is the case for its prokaryotic community, our knowledge of its viral component has remained limited by the paucity of information provided by studies dependent upon short sequence fragments. In this work, we attempt to fill this existing gap by using a combination of classical fosmid libraries with next-generation sequencing and assembly to recover long viral genomic fragments. We have sequenced ca. 6,000 fosmids from two metagenomics libraries made from prokaryotic biomass from the deep Mediterranean Sea and recovered twenty-eight complete viral genomes, all of them novel and quite distinct from all previously described viral genomes. They are preferentially found in deeper waters and are widely distributed all over the oceans. To our knowledge, this is the first report on complete and cosmopolitan viral genomes from the bathypelagic habitat.

A census of α-helical membrane proteins in double-stranded DNA viruses infecting bacteria and archaea

Background

Viruses are the most abundant and genetically diverse biological entities on earth, yet the repertoire of viral proteins remains poorly explored. As the number of sequenced virus genomes grows into the thousands, and the number of viral proteins into the hundreds of thousands, we report a systematic computational analysis of the point of first-contact between viruses and their hosts, namely viral transmembrane (TM) proteins.

Results

The complement of α-helical TM proteins in double-stranded DNA viruses infecting bacteria and archaea reveals large-scale trends that differ from those of their hosts. Viruses typically encode a substantially lower fraction of TM proteins than archaea or bacteria, with the notable exception of viruses with virions containing a lipid component such as a lipid envelope, internal lipid core, or inner membrane vesicle. Compared to bacteriophages, archaeal viruses are substantially enriched in membrane proteins. However, this feature is not always stable throughout the evolution of a viral lineage; for example, TM proteins are not part of the common heritage shared between Lipothrixviridae and Rudiviridae. In contrast to bacteria and archaea, viruses almost completely lack proteins with complicated membrane topologies composed of more than 4 TM segments, with the few detected exceptions being obvious cases of relatively recent horizontal transfer from the host.

Conclusions

The dramatic differences between the membrane proteomes of cells and viruses stem from the fact that viruses do not depend on essential membranes for energy transformation, ion homeostasis, nutrient transport and signaling.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0817-4) contains supplementary material, which is available to authorized users.

Viruses and cells intertwined since the dawn of evolution

Many attempts have been made to define nature of viruses and to uncover their origin. Our aim within this work was to show that there are different perceptions of viruses and many concepts to explain their emergence: the virus-first concept (also called co-evolution), the escape and the reduction theories. Moreover, a relatively new concept of polyphyletic virus origin called “three RNA cells, three DNA viruses” proposed by Forterre is described herein. In this paper, not only is each thesis supported by a body of evidence but also counter-argued in the light of various findings to give more insightful considerations to the readers. As the origin of viruses and that of living cells are most probably interdependent, we decided to reveal ideas concerning nature of cellular last universal common ancestor (LUCA). Furthermore, we discuss monophyletic ancestry of cellular domains and their relationships at the molecular level of membrane lipids and replication strategies of these three types of cells. In this review, we also present the emergence of DNA viruses requiring an evolutionary transition from RNA to DNA and recently discovered giant DNA viruses possibly involved in eukaryogenesis. In the course of evolution viruses emerged many times. They have always played a key role through horizontal gene transfer in evolutionary events and in formation of the tree of life or netlike routes of evolution providing a great deal of genetic diversity. In our opinion, future findings are crucial to better understand past relations between viruses and cells and the origin of both.

Analysis of virus genomes from glacial environments reveals novel virus groups with unusual host interactions

Microbial communities in glacial ecosystems are diverse, active, and subjected to strong viral pressures and infection rates. In this study we analyse putative virus genomes assembled from three dsDNA viromes from cryoconite hole ecosystems of Svalbard and the Greenland Ice Sheet to assess the potential hosts and functional role viruses play in these habitats. We assembled 208 million reads from the virus-size fraction and developed a procedure to select genuine virus scaffolds from cellular contamination. Our curated virus library contained 546 scaffolds up to 230 Kb in length, 54 of which were circular virus consensus genomes. Analysis of virus marker genes revealed a wide range of viruses had been assembled, including bacteriophages, cyanophages, nucleocytoplasmic large DNA viruses and a virophage, with putative hosts identified as Cyanobacteria, Alphaproteobacteria, Gammaproteobacteria, Actinobacteria, Firmicutes, eukaryotic algae and amoebae. Whole genome comparisons revealed the majority of circular genome scaffolds (CGS) formed 12 novel groups, two of which contained multiple phage members with plasmid-like properties, including a group of phage-plasmids possessing plasmid-like partition genes and toxin-antitoxin addiction modules to ensure their replication and a satellite phage-plasmid group. Surprisingly we also assembled a phage that not only encoded plasmid partition genes, but a clustered regularly interspaced short palindromic repeat (CRISPR)/Cas adaptive bacterial immune system. One of the spacers was an exact match for another phage in our virome, indicating that in a novel use of the system, the lysogen was potentially capable of conferring immunity on its bacterial host against other phage. Together these results suggest that highly novel and diverse groups of viruses are present in glacial environments, some of which utilize very unusual life strategies and genes to control their replication and maintain a long-term relationship with their hosts.

Effects of actin-like proteins encoded by two Bacillus pumilus phages on unstable lysogeny, revealed by genomic analysis

We characterized two newly isolated myoviruses, Bp8p-C and Bp8p-T, infecting the ginger rhizome rot disease pathogen Bacillus pumilus GR8. The plaque of Bp8p-T exhibited a clear center with a turbid rim, suggesting that Bp8p-T could transform into latent phage. Lysogeny assays showed that both the two phages could form latent states, while Bp8p-T could form latent phage at a higher frequency and stability than Bp8p-C. The genomes of Bp8p-C and Bp8p-T were 151,417 and 151,419 bp, respectively; both encoded 212 putative proteins, and only differed by three nucleotides. Moreover, owing to this difference, Bp8p-C encoded a truncated, putative actin-like plasmid segregation protein Gp27-C. Functional analysis of protein Gp27 showed that Gp27-T encoded by Bp8p-T exhibited higher ATPase activity and assembly ability than Gp27-C. The results indicate that the difference in Gp27 affected the phage lysogenic ability. Structural proteome analysis of Bp8p-C virion resulted in the identification of 14 structural proteins, among which a pectin lyase-like protein, a putative poly-gamma-glutamate hydrolase, and three proteins with unknown function, were firstly identified as components of the phage virion. Both phages exhibited specific lytic ability to the host strain GR8. Bp8p-C showed better control effect on the pathogen in ginger rhizome slices than Bp8p-T, suggesting that Bp8p-C has a potential application in bio-control of ginger rhizome rot disease.

Extensive gene remodeling in the viral world: new evidence for nongradual evolution in the mobilome network

Complex nongradual evolutionary processes such as gene remodeling are difficult to model, to visualize, and to investigate systematically. Despite these challenges, the creation of composite (or mosaic) genes by combination of genetic segments from unrelated gene families was established as an important adaptive phenomena in eukaryotic genomes. In contrast, almost no general studies have been conducted to quantify composite genes in viruses. Although viral genome mosaicism has been well-described, the extent of gene mosaicism and its rules of emergence remain largely unexplored. Applying methods from graph theory to inclusive similarity networks, and using data from more than 3,000 complete viral genomes, we provide the first demonstration that composite genes in viruses are 1) functionally biased, 2) involved in key aspects of the arm race between cells and viruses, and 3) can be classified into two distinct types of composite genes in all viral classes. Beyond the quantification of the widespread recombination of genes among different viruses of the same class, we also report a striking sharing of genetic information between viruses of different classes and with different nucleic acid types. This latter discovery provides novel evidence for the existence of a large and complex mobilome network, which appears partly bound by the sharing of genetic information and by the formation of composite genes between mobile entities with different genetic material. Considering that there are around 10E31 viruses on the planet, gene remodeling appears as a hugely significant way of generating and moving novel sequences between different kinds of organisms on Earth.

Identification and assembly of genomes and genetic elements in complex metagenomic samples without using reference genomes

Most current approaches for analyzing metagenomic data rely on comparisons to reference genomes, but the microbial diversity of many environments extends far beyond what is covered by reference databases. De novo segregation of complex metagenomic data into specific biological entities, such as particular bacterial strains or viruses, remains a largely unsolved problem. Here we present a method, based on binning co-abundant genes across a series of metagenomic samples, that enables comprehensive discovery of new microbial organisms, viruses and co-inherited genetic entities and aids assembly of microbial genomes without the need for reference sequences. We demonstrate the method on data from 396 human gut microbiome samples and identify 7,381 co-abundance gene groups (CAGs), including 741 metagenomic species (MGS). We use these to assemble 238 high-quality microbial genomes and identify affiliations between MGS and hundreds of viruses or genetic entities. Our method provides the means for comprehensive profiling of the diversity within complex metagenomic samples.

Bacteriophages and Biofilms

Biofilms are an extremely common adaptation, allowing bacteria to colonize hostile environments. They present unique problems for antibiotics and biocides, both due to the nature of the extracellular matrix and to the presence within the biofilm of metabolically inactive persister cells. Such chemicals can be highly effective against planktonic bacterial cells, while being essentially ineffective against biofilms. By contrast, bacteriophages seem to have a greater ability to target this common form of bacterial growth. The high numbers of bacteria present within biofilms actually facilitate the action of bacteriophages by allowing rapid and efficient infection of the host and consequent amplification of the bacteriophage. Bacteriophages also have a number of properties that make biofilms susceptible to their action. They are known to produce (or to be able to induce) enzymes that degrade the extracellular matrix. They are also able to infect persister cells, remaining dormant within them, but re-activating when they become metabolically active. Some cultured biofilms also seem better able to support the replication of bacteriophages than comparable planktonic systems. It is perhaps unsurprising that bacteriophages, as the natural predators of bacteria, have the ability to target this common form of bacterial life.

Comparison of different assembly and annotation tools on analysis of simulated viral metagenomic communities in the gut

BACKGROUND

The main limitations in the analysis of viral metagenomes are perhaps the high genetic variability and the lack of information in extant databases. To address these issues, several bioinformatic tools have been specifically designed or adapted for metagenomics by improving read assembly and creating more sensitive methods for homology detection. This study compares the performance of different available assemblers and taxonomic annotation software using simulated viral-metagenomic data.

RESULTS

We simulated two 454 viral metagenomes using genomes from NCBI's RefSeq database based on the list of actual viruses found in previously published metagenomes. Three different assembly strategies, spanning six assemblers, were tested for performance: overlap-layout-consensus algorithms Newbler, Celera and Minimo; de Bruijn graphs algorithms Velvet and MetaVelvet; and read probabilistic model Genovo. The performance of the assemblies was measured by the length of resulting contigs (using N50), the percentage of reads assembled and the overall accuracy when comparing against corresponding reference genomes. Additionally, the number of chimeras per contig and the lowest common ancestor were estimated in order to assess the effect of assembling on taxonomic and functional annotation. The functional classification of the reads was evaluated by counting the reads that correctly matched the functional data previously reported for the original genomes and calculating the number of over-represented functional categories in chimeric contigs. The sensitivity and specificity of tBLASTx, PhymmBL and the k-mer frequencies were measured by accurate predictions when comparing simulated reads against the NCBI Virus genomes RefSeq database.

CONCLUSIONS

Assembling improves functional annotation by increasing accurate assignations and decreasing ambiguous hits between viruses and bacteria. However, the success is limited by the chimeric contigs occurring at all taxonomic levels. The assembler and its parameters should be selected based on the focus of each study. Minimo's non-chimeric contigs and Genovo's long contigs excelled in taxonomy assignation and functional annotation, respectively.tBLASTx stood out as the best approach for taxonomic annotation for virus identification. PhymmBL proved useful in datasets in which no related sequences are present as it uses genomic features that may help identify distant taxa. The k-frequencies underperformed in all viral datasets.

Expanding the marine virosphere using metagenomics

Viruses infecting prokaryotic cells (phages) are the most abundant entities of the biosphere and contain a largely uncharted wealth of genomic diversity. They play a critical role in the biology of their hosts and in ecosystem functioning at large. The classical approaches studying phages require isolation from a pure culture of the host. Direct sequencing approaches have been hampered by the small amounts of phage DNA present in most natural habitats and the difficulty in applying meta-omic approaches, such as annotation of small reads and assembly. Serendipitously, it has been discovered that cellular metagenomes of highly productive ocean waters (the deep chlorophyll maximum) contain significant amounts of viral DNA derived from cells undergoing the lytic cycle. We have taken advantage of this phenomenon to retrieve metagenomic fosmids containing viral DNA from a Mediterranean deep chlorophyll maximum sample. This method allowed description of complete genomes of 208 new marine phages. The diversity of these genomes was remarkable, contributing 21 genomic groups of tailed bacteriophages of which 10 are completely new. Sequence based methods have allowed host assignment to many of them. These predicted hosts represent a wide variety of important marine prokaryotic microbes like members of SAR11 and SAR116 clades, Cyanobacteria and also the newly described low GC Actinobacteria. A metavirome constructed from the same habitat showed that many of the new phage genomes were abundantly represented. Furthermore, other available metaviromes also indicated that some of the new phages are globally distributed in low to medium latitude ocean waters. The availability of many genomes from the same sample allows a direct approach to viral population genomics confirming the remarkable mosaicism of phage genomes.

Narrow-host-range bacteriophages that infect Rhizobium etli associate with distinct genomic types

In this work, we isolated and characterized 14 bacteriophages that infect Rhizobium etli. They were obtained from rhizosphere soil of bean plants from agricultural lands in Mexico using an enrichment method. The host range of these phages was narrow but variable within a collection of 48 R. etli strains. We obtained the complete genome sequence of nine phages. Four phages were resistant to several restriction enzymes and in vivo cloning, probably due to nucleotide modifications. The genome size of the sequenced phages varied from 43 kb to 115 kb, with a median size of ≈ 45 to 50 kb. A large proportion of open reading frames of these phage genomes (65 to 70%) consisted of hypothetical and orphan genes. The remainder encoded proteins needed for phage morphogenesis and DNA synthesis and processing, among other functions, and a minor percentage represented genes of bacterial origin. We classified these phages into four genomic types on the basis of their genomic similarity, gene content, and host range. Since there are no reports of similar sequences, we propose that these bacteriophages correspond to novel species.

Genome sequence analysis indicates that the model eukaryote Nematostella vectensis harbors bacterial consorts

Analysis of the genome sequence of the starlet sea anemone, Nematostella vectensis, reveals many genes whose products are phylogenetically closer to proteins encoded by bacteria or bacteriophages than to any metazoan homologs. One explanation for such sequence affinities could be that these genes have been horizontally transferred from bacteria to the Nematostella lineage. We show, however, that bacterium-like and phage-like genes sequenced by the N. vectensis genome project tend to cluster on separate scaffolds, which typically do not include eukaryotic genes and differ from the latter in their GC contents. Moreover, most of the bacterium-like genes in N. vectensis either lack introns or the introns annotated in such genes are false predictions that, when translated, often restore the missing portions of their predicted protein products. In a freshwater cnidarian, Hydra, for which a proteobacterial endosymbiont is known, these gene features have been used to delineate the DNA of that endosymbiont sampled by the genome sequencing project. We predict that a large fraction of bacterium-like genes identified in the N. vectensis genome similarly are drawn from the contemporary bacterial consorts of the starlet sea anemone. These uncharacterized bacteria associated with N. vectensis are a proteobacterium and a representative of the phylum Bacteroidetes, each represented in the database by an apparently random sample of informational and operational genes. A substantial portion of a putative bacteriophage genome was also detected, which would be especially unlikely to have been transferred to a eukaryote.

The global virome: not as big as we thought?

Viruses likely infect all organisms, serving to unknown extent as genetic vectors in complex networks of organisms. Environmental virologists have revealed that these abundant nanoscale entities are global players with critical roles in every ecosystem investigated. Curiously, novel genes dominate viral genomes and metagenomes, which has led to the suggestion that viruses represent the largest reservoir of unexplored genetic material on Earth with literature estimates, extrapolating from 14 mycobacteriophage genomes, suggesting that two billion phage-encoded ORFs remain to be discovered. Here we examine (meta)genomic data available in the decade since this provocative assertion, and use 'protein clusters' to evaluate whether sampling technologies have advanced to the point that we may be able to sample 'all' of viral diversity in nature.

Orthologous gene clusters and taxon signature genes for viruses of prokaryotes

Viruses are the most abundant biological entities on earth and encompass a vast amount of genetic diversity. The recent rapid increase in the number of sequenced viral genomes has created unprecedented opportunities for gaining new insight into the structure and evolution of the virosphere. Here, we present an update of the phage orthologous groups (POGs), a collection of 4,542 clusters of orthologous genes from bacteriophages that now also includes viruses infecting archaea and encompasses more than 1,000 distinct virus genomes. Analysis of this expanded data set shows that the number of POGs keeps growing without saturation and that a substantial majority of the POGs remain specific to viruses, lacking homologues in prokaryotic cells, outside known proviruses. Thus, the great majority of virus genes apparently remains to be discovered. A complementary observation is that numerous viral genomes remain poorly, if at all, covered by POGs. The genome coverage by POGs is expected to increase as more genomes are sequenced. Taxon-specific, single-copy signature genes that are not observed in prokaryotic genomes outside detected proviruses were identified for two-thirds of the 57 taxa (those with genomes available from at least 3 distinct viruses), with half of these present in all members of the respective taxon. These signatures can be used to specifically identify the presence and quantify the abundance of viruses from particular taxa in metagenomic samples and thus gain new insights into the ecology and evolution of viruses in relation to their hosts.

Antibiotic resistance: from Darwin to Lederberg to Keynes

The emergence and spread of antibiotic-resistant bacteria reflects both, a gradual, completely Darwinian evolution, which mostly yields slight decreases in antibiotic susceptibility, along with phenotypes that are not precisely characterized as "resistance"; and sudden changes, from full susceptibility to full resistance, which are driven by a vast array of horizontal gene transfer mechanisms. Antibiotics select for more than just antibiotic resistance (i.e., increased virulence and enhanced gene exchange abilities); and many non-antibiotic agents or conditions select for or maintain antibiotic resistance traits as a result of a complex network of underlying and often overlapping mechanisms. Thus, the development of new antibiotics and thoughtful, integrated anti-infective strategies is needed to address the immediate and long-term threat of antibiotic resistance. Since the biology of resistance is complex, these new drugs and strategies will not come from free-market forces, or from "incentives" for pharmaceutical companies.

Contribution of phage-derived genomic islands to the virulence of facultative bacterial pathogens

Facultative pathogens have extremely dynamic pan-genomes, to a large extent derived from bacteriophages and other mobile elements. We developed a simple approach to identify phage-derived genomic islands and apply it to show that pathogens from diverse bacterial genera are significantly enriched in clustered phage-derived genes compared with related benign strains. These findings show that genome expansion by integration of prophages containing virulence factors is a major route of evolution of facultative bacterial pathogens.

Computational methods for Gene Orthology inference

Accurate inference of orthologous genes is a pre-requisite for most comparative genomics studies, and is also important for functional annotation of new genomes. Identification of orthologous gene sets typically involves phylogenetic tree analysis, heuristic algorithms based on sequence conservation, synteny analysis, or some combination of these approaches. The most direct tree-based methods typically rely on the comparison of an individual gene tree with a species tree. Once the two trees are accurately constructed, orthologs are straightforwardly identified by the definition of orthology as those homologs that are related by speciation, rather than gene duplication, at their most recent point of origin. Although ideal for the purpose of orthology identification in principle, phylogenetic trees are computationally expensive to construct for large numbers of genes and genomes, and they often contain errors, especially at large evolutionary distances. Moreover, in many organisms, in particular prokaryotes and viruses, evolution does not appear to have followed a simple 'tree-like' mode, which makes conventional tree reconciliation inapplicable. Other, heuristic methods identify probable orthologs as the closest homologous pairs or groups of genes in a set of organisms. These approaches are faster and easier to automate than tree-based methods, with efficient implementations provided by graph-theoretical algorithms enabling comparisons of thousands of genomes. Comparisons of these two approaches show that, despite conceptual differences, they produce similar sets of orthologs, especially at short evolutionary distances. Synteny also can aid in identification of orthologs. Often, tree-based, sequence similarity- and synteny-based approaches can be combined into flexible hybrid methods.