Evaluation of the genomic diversity of viruses infecting bacteria, archaea and eukaryotes using a common bioinformatic platform: steps towards a unified taxonomy

How to introduce a new bacteriophage on the block: a short guide to phage classification

Bacteriophage (phage) studies established the field of molecular biology and continue to propel life science research forward due to their diversity, abundance, and potential applications. In this Gem article, we orient newcomers to four common ways phages are currently classified: infection cycle, morphology, taxonomy, and supergroup. By using these classifications, researchers can determine where any novel phage fits into the scheme of the known "phage-verse".

Exploring the historical roots, advantages and efficacy of phage therapy in plant diseases management

In the drug-resistance era, phage therapy has received considerable attention from worldwide researchers. Phage therapy has been given much attention in public health but is rarely applied to control plant diseases. Herein, we discuss phage therapy as a biocontrol approach against several plant diseases. The emergence of antibiotic resistance in agriculturally important pathogenic bacteria and the toxic nature of different synthetic compounds used to control microbes has driven researchers to rethink the century-old strategy of phage therapy''. Compared to other treatment strategies, phage therapy offers remarkable advantages such as high specificity, less chances of drug resistance, non-harmful nature, and benefit to soil microbial flora. The optimizations and protective formulations of phages are significant accomplishments; however, steps towards a better understanding of the physiologic characteristics of phages need to be preceded to commercialize their use. The future of phage therapy in the context of plant disease management is promising and could play a significant role in sustainable agriculture. Ongoing research will likely affirm the safety of phage therapy, ensuring that it does not harm non-target organisms, including beneficial soil microbes. Phage therapy could become vital in addressing global food security challenges, particularly in regions heavily impacted by plant bacterial diseases. Efforts to create formulations that enhance the stability and shelf-life of phages will be crucial, especially for their use in varied environmental conditions.

Novel metagenomics analysis of stony coral tissue loss disease

Stony coral tissue loss disease (SCTLD) has devastated coral reefs off the coast of Florida and continues to spread throughout the Caribbean. Although a number of bacterial taxa have consistently been associated with SCTLD, no pathogen has been definitively implicated in the etiology of SCTLD. Previous studies have predominantly focused on the prokaryotic community through 16S rRNA sequencing of healthy and affected tissues. Here, we provide a different analytical approach by applying a bioinformatics pipeline to publicly available metagenomic sequencing samples of SCTLD lesions and healthy tissues from 4 stony coral species. To compensate for the lack of coral reference genomes, we used data from apparently healthy coral samples to approximate a host genome and healthy microbiome reference. These reads were then used as a reference to which we matched and removed reads from diseased lesion tissue samples, and the remaining reads associated only with disease lesions were taxonomically classified at the DNA and protein levels. For DNA classifications, we used a pathogen identification protocol originally designed to identify pathogens in human tissue samples, and for protein classifications, we used a fast protein sequence aligner. To assess the utility of our pipeline, a species-level analysis of a candidate genus, Vibrio, was used to demonstrate the pipeline's effectiveness. Our approach revealed both complementary and unique coral microbiome members compared with a prior metagenome analysis of the same dataset.

Characterization and genomic analysis of a lytic Stenotrophomonas maltophilia short-tailed phage A1432 revealed a new genus of the family Mesyanzhinovviridae

Stenotrophomonas maltophilia (S. maltophilia) is an emerging opportunistic pathogen that exhibits resistant to a majority of commonly used antibiotics. Phages have the potential to serve as an alternative treatment for S. maltophilia infections. In this study, a lytic phage, A1432, infecting S. maltophilia YCR3A-1, was isolated and characterized from a karst cave. Transmission electron microscopy revealed that phage A1432 possesses an icosahedral head and a shorter tail. Phage A1432 demonstrated a narrow host range, with an optimal multiplicity of infection of 0.1. The one-step growth curve indicated a latent time of 10 min, a lysis period of 90 min, a burst size of 43.2 plaque-forming units per cell. In vitro bacteriolytic activity test showed that phage A1432 was capable to inhibit the growth of S. maltophilia YCR3A-1 in an MOI-dependent manner after 2 h of co-culture. BLASTn analysis showed that phage A1432 genome shares the highest similarity (81.46%) with Xanthomonas phage Xoo-sp2 in the NCBI database, while the query coverage was only 37%. The phage contains double-stranded DNA with a genome length of 61,660 bp and a GC content of 61.92%. It is predicted to have 79 open reading frames and one tRNA, with no virulence or antibiotic resistance genes. Phylogenetic analysis using terminase large subunit and DNA polymerase indicated that phage A1432 clustered with members of the Bradleyvirinae subfamily but diverged into a distinct branch. Further phylogenetic comparison analysis using Average Nucleotide Identity, proteomic phylogenetic analysis, genomic network analysis confirmed that phage A1432 belongs to a novel genus within the Bradleyvirinae subfamily, Mesyanzhinovviridae family. Additionally, phylogenetic analysis of the so far isolated S. maltophilia phages revealed significant genetic diversity among these phages. The results of this research will contribute valuable information for further studies on their morphological and genetic diversity, will aid in elucidating the evolutionary mechanisms that give rise to them.

Novel metagenomics analysis of stony coral tissue loss disease

Stony coral tissue loss disease (SCTLD) has devastated coral reefs off the coast of Florida and continues to spread throughout the Caribbean. Although a number of bacterial taxa have consistently been associated with SCTLD, no pathogen has been definitively implicated in the etiology of SCTLD. Previous studies have predominantly focused on the prokaryotic community through 16S rRNA sequencing of healthy and affected tissues. Here, we provide a different analytical approach by applying a bioinformatics pipeline to publicly available metagenomic sequencing samples of SCTLD lesions and healthy tissues from four stony coral species. To compensate for the lack of coral reference genomes, we used data from apparently healthy coral samples to approximate a host genome and healthy microbiome reference. These reads were then used as a reference to which we matched and removed reads from diseased lesion tissue samples, and the remaining reads associated only with disease lesions were taxonomically classified at the DNA and protein levels. For DNA classifications, we used a pathogen identification protocol originally designed to identify pathogens in human tissue samples, and for protein classifications, we used a fast protein sequence aligner. To assess the utility of our pipeline, a species-level analysis of a candidate genus, Vibrio, was used to demonstrate the pipeline's effectiveness. Our approach revealed both complementary and unique coral microbiome members compared to a prior metagenome analysis of the same dataset.

The virtue of training: extending phage host spectra against vancomycin-resistant Enterococcus faecium strains using the Appelmans method

Phage therapy has (re)emerged as a serious possibility for combating multidrug-resistant bacterial infections, including those caused by vancomycin-resistant Enterococcus faecium strains. These opportunistic pathogens belong to a specific clonal complex 17, against which relatively few phages have been screened. We isolated a collection of 21 virulent phages growing on these vancomycin-resistant isolates. Each of these phages harbored a typical narrow plaquing host range, lysing at most 5 strains and covering together 10 strains of our panel of 14 clinical isolates. To enlarge the host spectrum of our phages, the Appelmans protocol was used. We mixed four out of our most complementary phages in a cocktail that we iteratively grew on eight naive strains from our panel, of which six were initially refractory to at least three of the combined phages. Fifteen successive passages permitted to significantly improve the lytic activity of the cocktail, from which phages with extended host ranges within the E. faecium species could be isolated. A single evolved phage able to kill up to 10 of the 14 initial E. faecium strains was obtained, and it barely infected nearby species. All evolved phages had acquired point mutations or a recombination event in the tail fiber genetic region, suggesting these genes might have driven phage evolution by contributing to their extended host spectra.

Efficient Recovery of Complete Gut Viral Genomes by Combined Short- and Long-Read Sequencing

Current metagenome assembled human gut phage catalogs contained mostly fragmented genomes. Here, comprehensive gut virome detection procedure is developed involving virus-like particle (VLP) enrichment from ≈500 g feces and combined sequencing of short- and long-read. Applied to 135 samples, a Chinese Gut Virome Catalog (CHGV) is assembled consisting of 21,499 non-redundant viral operational taxonomic units (vOTUs) that are significantly longer than those obtained by short-read sequencing and contained ≈35% (7675) complete genomes, which is ≈nine times more than those in the Gut Virome Database (GVD, ≈4%, 1,443). Interestingly, the majority (≈60%, 13,356) of the CHGV vOTUs are obtained by either long-read or hybrid assemblies, with little overlap with those assembled from only the short-read data. With this dataset, vast diversity of the gut virome is elucidated, including the identification of 32% (6,962) novel vOTUs compare to public gut virome databases, dozens of phages that are more prevalent than the crAssphages and/or Gubaphages, and several viral clades that are more diverse than the two. Finally, the functional capacities are also characterized of the CHGV encoded proteins and constructed a viral-host interaction network to facilitate future research and applications.

Accumulation of defense systems in phage-resistant strains of Pseudomonas aeruginosa

Prokaryotes encode multiple distinct anti-phage defense systems in their genomes. However, the impact of carrying a multitude of defense systems on phage resistance remains unclear, especially in a clinical context. Using a collection of antibiotic-resistant clinical strains of Pseudomonas aeruginosa and a broad panel of phages, we demonstrate that defense systems contribute substantially to defining phage host range and that overall phage resistance scales with the number of defense systems in the bacterial genome. We show that many individual defense systems target specific phage genera and that defense systems with complementary phage specificities co-occur in P. aeruginosa genomes likely to provide benefits in phage-diverse environments. Overall, we show that phage-resistant phenotypes of P. aeruginosa with at least 19 phage defense systems exist in the populations of clinical, antibiotic-resistant P. aeruginosa strains.

Bacteriophage Taxonomy: A Continually Evolving Discipline

While taxonomy is an often underappreciated branch of science, it serves very important roles. Bacteriophage taxonomy has evolved from a discipline based mainly on morphology, characterized by the work of David Bradley and Hans-Wolfgang Ackermann, to the sequence-based approach that is taken today. The Bacterial Viruses Subcommittee of the International Committee on Taxonomy of Viruses (ICTV) takes a holistic approach to classifying prokaryote viruses by measuring overall DNA and protein similarity and phylogeny before making decisions about the taxonomic position of a new virus. The huge number of complete genomes being deposited with the National Center for Biotechnology Information (NCBI) and other public databases has resulted in a reassessment of the taxonomy of many viruses, and the future will see the introduction of new viral families and higher orders.

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.

Identification of Structural and Morphogenesis Genes of Sulfitobacter Phage ΦGT1 and Placement within the Evolutionary History of the Podoviruses

ΦGT1 is a lytic podovirus of an alphaproteobacterial Sulfitobacter species, with few closely matching sequences among characterized phages, thus defying a useful description by simple sequence clustering methods. The history of the ΦGT1 core structure module was reconstructed using timetrees, including numerous related prospective prophages, to flesh out the evolutionary lineages spanning from the origin of the ejectosomal podovirus >3.2 Gya to the present genes of ΦGT1 and its closest relatives. A peculiarity of the ΦGT1 structural proteome is that it contains two paralogous tubular tail A (tubeA) proteins. The origin of the dual tubeA arrangement was traced to a recombination between two more ancient podoviral lineages occurring ~0.7 Gya in the alphaproteobacterial order Rhizobiales. Descendants of the ancestral dual A recombinant were tracked forward forming both temperate and lytic phage clusters and exhibiting both vertical transmission with patchy persistence and horizontal transfer with respect to host taxonomy. The two ancestral lineages were traced backward, making junctions with a major metagenomic podoviral family, the LUZ24-like gammaproteobacterial phages, and Myxococcal phage Mx8, and finally joining near the origin of podoviruses with P22. With these most conservative among phage genes, deviations from uncomplicated vertical and nonrecombinant descent are numerous but countable. The use of timetrees allowed conceptualization of the phage's evolution in the context of a sequence of ancestors spanning the time of life on Earth.

VirClust-A Tool for Hierarchical Clustering, Core Protein Detection and Annotation of (Prokaryotic) Viruses

Recent years have seen major changes in the classification criteria and taxonomy of viruses. The current classification scheme, also called "megataxonomy of viruses", recognizes six different viral realms, defined based on the presence of viral hallmark genes (VHGs). Within the realms, viruses are classified into hierarchical taxons, ideally defined by the phylogeny of their shared genes. To enable the detection of shared genes, viruses have first to be clustered, and there is currently a need for tools to assist with virus clustering and classification. Here, VirClust is presented. It is a novel, reference-free tool capable of performing: (i) protein clustering, based on BLASTp and Hidden Markov Models (HMMs) similarities; (ii) hierarchical clustering of viruses based on intergenomic distances calculated from their shared protein content; (iii) identification of core proteins and (iv) annotation of viral proteins. VirClust has flexible parameters both for protein clustering and for splitting the viral genome tree into smaller genome clusters, corresponding to different taxonomic levels. Benchmarking on a phage dataset showed that the genome trees produced by VirClust match the current ICTV classification at family, sub-family and genus levels. VirClust is freely available, as a web-service and stand-alone tool.

Isolation, screening and characterization of phage

Bacterial resistance threatens public health due to a lack of novel antibacterial classes since the 21st century. Bacteriophages, the most ubiquitous microorganism on Earth and natural predators of bacteria, have the potential to save the world from the post-antibiotic era. Therefore, phage isolation and characterization are in high demand to find suitable phages for therapeutic and bacterial control applications. The chapter presents brief guidance supported by recommendations on the isolation of phages, and initial screening of phage antimicrobial efficacy, in addition to, conducting comprehensive characterization addressing morphological, biological, genomic, and taxonomic features.

Genome Characterization and Infectivity Potential of Vibriophage-ϕLV6 with Lytic Activity against Luminescent Vibrios of Penaeus vannamei Shrimp Aquaculture

Shrimp aquaculture, especially during the hatchery phase, is prone to economic losses due to infections caused by luminescent vibrios. In the wake of antimicrobial resistance (AMR) in bacteria and the food safety requirements of farmed shrimp, aqua culturists are seeking alternatives to antibiotics for shrimp health management, and bacteriophages are fast emerging as natural and bacteria-specific antimicrobial agents. This study analyzed the whole genome of vibriophage-ϕLV6 that showed lytic activity against six luminescent vibrios isolated from the larval tanks of P. vannamei shrimp hatcheries. The Vibriophage-ϕLV6 genome was 79,862 bp long with 48% G+C content and 107 ORFs that coded for 31 predicted protein functions, 75 hypothetical proteins, and a tRNA. Pertinently, the vibriophage-ϕLV6 genome harbored neither AMR determinants nor virulence genes, indicating its suitability for phage therapy. There is a paucity of whole genome-based information on vibriophages that lyse luminescent vibrios, and this study adds pertinent data to the database of V. harveyi infecting phage genomes and, to our knowledge, is the first vibriophage genome report from India. Transmission electron microscopy (TEM) of vibriophage-ϕLV6 revealed an icosahedral head (~73 nm) and a long, flexible tail (~191 nm) suggesting siphovirus morphology. The vibriophage-ϕLV6 phage at a multiplicity of infection (MOI) of 80 inhibited the growth of luminescent V. harveyi at 0.25%, 0.5%, 1%, 1.5%, 2%, 2.5%, and 3% salt gradients. In vivo experiments conducted with post-larvae of shrimp showed that vibriophage-ϕLV6 reduced luminescent vibrio counts and post-larval mortalities in the phage-treated tank compared to the bacteria-challenged tank, suggesting the potentiality of vibriophage-ϕLV6 as a promising candidate in treating luminescent vibriosis in shrimp aquaculture. The vibriophage-ϕLV6 survived for 30 days in salt (NaCl) concentrations ranging from 5 ppt to 50 ppt and was stable at 4 °C for 12 months.

Putative Mitoviruses without In-Frame UGA(W) Codons: Evolutionary Implications

Mitoviruses are small vertically transmitted RNA viruses found in fungi, plants and animals. Taxonomically, a total of 105 species and 4 genera have been formally recognized by ICTV, and recently, 18 new putative species have been included in a new proposed genus. Transcriptomic and metatranscriptomic studies are a major source of countless new virus-like sequences that are continually being added to open databases and these may be good sources for identifying new putative mitoviruses. The search for mitovirus-like sequences in the NCBI databases resulted in the discovery of more than one hundred new putative mitoviruses, with important implications for taxonomy and also for the evolutionary scenario. Here, we propose the inclusion of four new putative members to the genus Kvaramitovirus, and the existence of a new large basally divergent lineage composed of 144 members that lack internal UGA codons (subfamily "Arkeomitovirinae"), a feature not shared by the vast majority of mitoviruses. Finally, a taxonomic categorization proposal and a detailed description of the evolutionary history of mitoviruses were carried out. This in silico study supports the hypothesis of the existence of a basally divergent lineage that could have had an impact on the early evolutionary history of mitoviruses.

Abolishment of morphology-based taxa and change to binomial species names: 2022 taxonomy update of the ICTV bacterial viruses subcommittee

This article summarises the activities of the Bacterial Viruses Subcommittee of the International Committee on Taxonomy of Viruses for the period of March 2021-March 2022. We provide an overview of the new taxa proposed in 2021, approved by the Executive Committee, and ratified by vote in 2022. Significant changes to the taxonomy of bacterial viruses were introduced: the paraphyletic morphological families Podoviridae, Siphoviridae, and Myoviridae as well as the order Caudovirales were abolished, and a binomial system of nomenclature for species was established. In addition, one order, 22 families, 30 subfamilies, 321 genera, and 862 species were newly created, promoted, or moved.

The Emerging Role of the Gut Virome in Health and Inflammatory Bowel Disease: Challenges, Covariates and a Viral Imbalance

Virome research is a rapidly growing area in the microbiome field that is increasingly associated with human diseases, such as inflammatory bowel disease (IBD). Although substantial progress has been made, major methodological challenges limit our understanding of the virota. In this review, we describe challenges that must be considered to accurately report the virome composition and the current knowledge on the virome in health and IBD. First, the description of the virome shows strong methodological biases related to wetlab (e.g., VLP enrichment) and bioinformatics approaches (viral identification and classification). Second, IBD patients show consistent viral imbalances characterized by a high relative abundance of phages belonging to the Caudovirales and a low relative abundance of phages belonging to the Microviridae. Simultaneously, a sporadic contraction of CrAss-like phages and a potential expansion of the lysogenic potential of the intestinal virome are observed. Finally, despite numerous studies that have conducted diversity analysis, it is difficult to draw firm conclusions due to methodological biases. Overall, we present the many methodological and environmental factors that influence the virome, its current consensus in health and IBD, and a contributing hypothesis called the "positive inflammatory feedback loop" that may play a role in the pathophysiology of IBD.

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.

Isolation and characterization of novel Fusobacterium nucleatum bacteriophages

Fusobacterium nucleatum is a strictly anaerobic, Gram-negative bacterial species that is a member of the commensal flora in the oral cavity and gut. Recent studies suggested that the increase of abundance is associated with the development of various diseases, among which colorectal cancer is of the biggest concerns. Phage therapy is regarded as a potential approach to control the number of F. nucleatum, which may contribute to the prevention and treatment of related diseases. In this study, we isolated five isolates of bacteriophage targeting F. nucleatum. The morphological, biological, genomic and functional characteristics of five bacteriophages were investigated. Transmission electron microscopy revealed that JD-Fnp1 ~ JD-Fnp5 are all myoviruses. The size of the JD-Fnp1 ~ JD-Fnp5 genomes was 180,066 bp (JD-Fnp1), 41,329 bp (JD-Fnp2), 38,962 bp (JD-Fnp3), 180,231 bp (JD-Fnp4), and 41,353 bp (JD-Fnp5) respectively. The biological features including pH and heat stability, host range, growth characteristics of JD-Fnp1 ~ JD-Fnp5 displayed different patterns. Among them, JD-Fnp4 is considered to have the greatest clinical application value. The identification and characterization of JD-Fnp1 ~ JD-Fnp5 provides a basis for subsequent therapeutic strategy exploration of F. nucleatum-related diseases.

Ayka, a Novel Curtobacterium Bacteriophage, Provides Protection against Soybean Bacterial Wilt and Tan Spot

Diseases caused by the Gram-positive bacterium Curtobacteriumflaccumfaciens pv. flaccumfaciens (Cff) inflict substantial economic losses in soybean cultivation. Use of specific bacterial viruses (bacteriophages) for treatment of seeds and plants to prevent the development of bacterial infections is a promising approach for bioprotection in agriculture. Phage control has been successfully tested for a number of staple crops. However, this approach has never been applied to treat bacterial diseases of legumes caused by Cff, and no specific bacteriophages have been known to date. This paper presents detailed characteristics of the first lytic bacteriophage infecting this pathogen. Phage Ayka, related to φ29-like (Salasmaviridae) viruses, but representing a new subfamily, was shown to control the development of bacterial wilt and tan spot in vitro and in greenhouse plants.

Cyanophages from a less virulent clade dominate over their sister clade in global oceans

Environmental virus communities are highly diverse. However, the infection physiology underlying the evolution of diverse phage lineages and their ecological consequences are largely unknown. T7-like cyanophages are abundant in nature and infect the marine unicellular cyanobacteria, Synechococcus and Prochlorococcus, important primary producers in the oceans. Viruses belonging to this genus are divided into two distinct phylogenetic clades: clade A and clade B. These viruses have narrow host-ranges with clade A phages primarily infecting Synechococcus genotypes, while clade B phages are more diverse and can infect either Synechococcus or Prochlorococcus genotypes. Here we investigated infection properties (life history traits) and environmental abundances of these two clades of T7-like cyanophages. We show that clade A cyanophages have more rapid infection dynamics, larger burst sizes and greater virulence than clade B cyanophages. However, clade B cyanophages were at least 10-fold more abundant in all seasons, and infected more cyanobacteria, than clade A cyanophages in the Red Sea. Models predicted that steady-state cyanophage abundances, infection frequency, and virus-induced mortality, peak at intermediate virulence values. Our findings indicate that differences in infection properties are reflected in virus phylogeny at the clade level. They further indicate that infection properties, together with differences in subclade diversity and host repertoire, have important ecological consequences with the less aggressive, more diverse virus clade having greater ecological impacts.

Complete Genome Sequence of a Novel Hafnia alvei Bacteriophage

A bacteriophage that is able to infect Hafnia alvei and can also infect two other hosts, Klebsiella pneumoniae and Salmonella enterica, was isolated from wastewater. The complete genome sequence was determined by long-read PacBio sequencing. Based on sequence similarity, the bacteriophage is assigned to the viral genus Kolesnikvirus within the family Myoviridae.

Highly diverse flavobacterial phages isolated from North Sea spring blooms

It is generally recognized that phages are a mortality factor for their bacterial hosts. This could be particularly true in spring phytoplankton blooms, which are known to be closely followed by a highly specialized bacterial community. We hypothesized that phages modulate these dense heterotrophic bacteria successions following phytoplankton blooms. In this study, we focused on Flavobacteriia, because they are main responders during these blooms and have an important role in the degradation of polysaccharides. A cultivation-based approach was used, obtaining 44 lytic flavobacterial phages (flavophages), representing twelve new species from two viral realms. Taxonomic analysis allowed us to delineate ten new phage genera and ten new families, from which nine and four, respectively, had no previously cultivated representatives. Genomic analysis predicted various life styles and genomic replication strategies. A likely eukaryote-associated host habitat was reflected in the gene content of some of the flavophages. Detection in cellular metagenomes and by direct-plating showed that part of these phages were actively replicating in the environment during the 2018 spring bloom. Furthermore, CRISPR/Cas spacers and re-isolation during two consecutive years suggested that, at least part of the new flavophages are stable components of the microbial community in the North Sea. Together, our results indicate that these diverse flavophages have the potential to modulate their respective host populations.

Bioinformatics of virus taxonomy: foundations and tools for developing sequence-based hierarchical classification

The genome sequence is the only characteristic readily obtainable for all known viruses, underlying the growing role of comparative genomics in organizing knowledge about viruses in a systematic evolution-aware way, known as virus taxonomy. Overseen by the International Committee on Taxonomy of Viruses (ICTV), development of virus taxonomy involves taxa demarcation at 15 ranks of a hierarchical classification, often in host-specific manner. Outside the ICTV remit, researchers assess fitting numerous unclassified viruses into the established taxa. They employ different metrics of virus clustering, basing on conserved domain(s), separation of viruses in rooted phylogenetic trees and pair-wise distance space. Computational approaches differ further in respect to methodology, number of ranks considered, sensitivity to uneven virus sampling, and visualization of results. Advancing and using computational tools will be critical for improving taxa demarcation across the virosphere and resolving rank origins in research that may also inform experimental virology.

Multiple novel filamentous phages detected in the cloacal swab samples of birds using viral metagenomics approach

Members of the family Inoviridae (inoviruses) are characterized by their unique filamentous morphology and infection cycle. The viral genome of inovirus is able to integrate into the host genome and continuously releases virions without lysing the host, establishing chronic infection. A large number of inoviruses have been obtained from microbial genomes and metagenomes recently, but putative novel inoviruses remaining to be identified. Here, using viral metagenomics, we identified four novel inoviruses from cloacal swab samples of wild and breeding birds. The circular genome of those four inoviruses are 6732 to 7709 nt in length with 51.4% to 56.5% GC content and encodes 9 to 13 open reading frames, respectively. The zonula occludens toxin gene implicated in the virulence of pathogenic host bacteria were identified in all four inoviruses and shared the highest amino acid sequences identity (< 37.3%) to other reference strains belonging to different genera of the family Inoviridae and among themselves. Phylogenetic analysis indicated that all the four inoviruses were genetically far away from other strains belonging to the family Inoviridae and formed an independent clade. According to the genetic distance-based criteria, all the four inoviruses identified in the present study respectively belong to four novel putative genera in the family Inoviridae.

Diversity, taxonomy, and evolution of archaeal viruses of the class Caudoviricetes

The archaeal tailed viruses (arTV), evolutionarily related to tailed double-stranded DNA (dsDNA) bacteriophages of the class Caudoviricetes, represent the most common isolates infecting halophilic archaea. Only a handful of these viruses have been genomically characterized, limiting our appreciation of their ecological impacts and evolution. Here, we present 37 new genomes of haloarchaeal tailed virus isolates, more than doubling the current number of sequenced arTVs. Analysis of all 63 available complete genomes of arTVs, which we propose to classify into 14 new families and 3 orders, suggests ancient divergence of archaeal and bacterial tailed viruses and points to an extensive sharing of genes involved in DNA metabolism and counterdefense mechanisms, illuminating common strategies of virus-host interactions with tailed bacteriophages. Coupling of the comparative genomics with the host range analysis on a broad panel of haloarchaeal species uncovered 4 distinct groups of viral tail fiber adhesins controlling the host range expansion. The survey of metagenomes using viral hallmark genes suggests that the global architecture of the arTV community is shaped through recurrent transfers between different biomes, including hypersaline, marine, and anoxic environments.

FlexTaxD: flexible modification of taxonomy databases for improved sequence classification

SUMMARY

The Flexible Taxonomy Database (FlexTaxD) framework provides a method for modification and merging official and custom taxonomic databases to create improved databases. Using such databases will increase accuracy and precision of existing methods to classify sequence reads.

AVAILABILITY AND IMPLEMENTATION

Source code is freely available at https://github.com/FOI-Bioinformatics/flextaxd and installable through Bioconda.

Bacteriophage classification for assembled contigs using graph convolutional network

Motivation

Bacteriophages (aka phages), which mainly infect bacteria, play key roles in the biology of microbes. As the most abundant biological entities on the planet, the number of discovered phages is only the tip of the iceberg. Recently, many new phages have been revealed using high-throughput sequencing, particularly metagenomic sequencing. Compared to the fast accumulation of phage-like sequences, there is a serious lag in taxonomic classification of phages. High diversity, abundance and limited known phages pose great challenges for taxonomic analysis. In particular, alignment-based tools have difficulty in classifying fast accumulating contigs assembled from metagenomic data.

Results

In this work, we present a novel semi-supervised learning model, named PhaGCN, to conduct taxonomic classification for phage contigs. In this learning model, we construct a knowledge graph by combining the DNA sequence features learned by convolutional neural network and protein sequence similarity gained from gene-sharing network. Then we apply graph convolutional network to utilize both the labeled and unlabeled samples in training to enhance the learning ability. We tested PhaGCN on both simulated and real sequencing data. The results clearly show that our method competes favorably against available phage classification tools.

Availability and implementation

The source code of PhaGCN is available via: https://github.com/KennthShang/PhaGCN.

Engineered Bacteriophage Therapeutics: Rationale, Challenges and Future

The current problems with increasing bacterial resistance to antibacterial therapies, resulting in a growing frequency of incurable bacterial infections, necessitates the acceleration of studies on antibacterials of a new generation that could offer an alternative to antibiotics or support their action. Bacteriophages (phages) can kill antibiotic-sensitive as well as antibiotic-resistant bacteria, and thus are a major subject of such studies. Their efficacy in curing bacterial infections has been demonstrated in in vivo experiments and in the clinic. Unlike antibiotics, phages have a narrow range of specificity, which makes them safe for commensal microbiota. However, targeting even only the most clinically relevant strains of pathogenic bacteria requires large collections of well characterized phages, whose specificity would cover all such strains. The environment is a rich source of diverse phages, but due to their complex relationships with bacteria and safety concerns, only some naturally occurring phages can be considered for therapeutic applications. Still, their number and diversity make a detailed characterization of all potentially promising phages virtually impossible. Moreover, no single phage combines all the features required of an ideal therapeutic agent. Additionally, the rapid acquisition of phage resistance by bacteria may make phages already approved for therapy ineffective and turn the search for environmental phages of better efficacy and new specificity into an endless race. An alternative strategy for acquiring phages with desired properties in a short time with minimal cost regarding their acquisition, characterization, and approval for therapy could be based on targeted genome modifications of phage isolates with known properties. The first example demonstrating the potential of this strategy in curing bacterial diseases resistant to traditional therapy is the recent successful treatment of a progressing disseminated Mycobacterium abscessus infection in a teenage patient with the use of an engineered phage. In this review, we briefly present current methods of phage genetic engineering, highlighting their advantages and disadvantages, and provide examples of genetically engineered phages with a modified host range, improved safety or antibacterial activity, and proven therapeutic efficacy. We also summarize novel uses of engineered phages not only for killing pathogenic bacteria, but also for in situ modification of human microbiota to attenuate symptoms of certain bacterial diseases and metabolic, immune, or mental disorders.

Thousands of previously unknown phages discovered in whole-community human gut metagenomes

BACKGROUND

Double-stranded DNA bacteriophages (dsDNA phages) play pivotal roles in structuring human gut microbiomes; yet, the gut virome is far from being fully characterized, and additional groups of phages, including highly abundant ones, continue to be discovered by metagenome mining. A multilevel framework for taxonomic classification of viruses was recently adopted, facilitating the classification of phages into evolutionary informative taxonomic units based on hallmark genes. Together with advanced approaches for sequence assembly and powerful methods of sequence analysis, this revised framework offers the opportunity to discover and classify unknown phage taxa in the human gut.

RESULTS

A search of human gut metagenomes for circular contigs encoding phage hallmark genes resulted in the identification of 3738 apparently complete phage genomes that represent 451 putative genera. Several of these phage genera are only distantly related to previously identified phages and are likely to found new families. Two of the candidate families, "Flandersviridae" and "Quimbyviridae", include some of the most common and abundant members of the human gut virome that infect Bacteroides, Parabacteroides, and Prevotella. The third proposed family, "Gratiaviridae," consists of less abundant phages that are distantly related to the families Autographiviridae, Drexlerviridae, and Chaseviridae. Analysis of CRISPR spacers indicates that phages of all three putative families infect bacteria of the phylum Bacteroidetes. Comparative genomic analysis of the three candidate phage families revealed features without precedent in phage genomes. Some "Quimbyviridae" phages possess Diversity-Generating Retroelements (DGRs) that generate hypervariable target genes nested within defense-related genes, whereas the previously known targets of phage-encoded DGRs are structural genes. Several "Flandersviridae" phages encode enzymes of the isoprenoid pathway, a lipid biosynthesis pathway that so far has not been known to be manipulated by phages. The "Gratiaviridae" phages encode a HipA-family protein kinase and glycosyltransferase, suggesting these phages modify the host cell wall, preventing superinfection by other phages. Hundreds of phages in these three and other families are shown to encode catalases and iron-sequestering enzymes that can be predicted to enhance cellular tolerance to reactive oxygen species.

CONCLUSIONS

Analysis of phage genomes identified in whole-community human gut metagenomes resulted in the delineation of at least three new candidate families of Caudovirales and revealed diverse putative mechanisms underlying phage-host interactions in the human gut. Addition of these phylogenetically classified, diverse, and distinct phages to public databases will facilitate taxonomic decomposition and functional characterization of human gut viromes. Video abstract.

A Roadmap for Genome-Based Phage Taxonomy

Bacteriophage (phage) taxonomy has been in flux since its inception over four decades ago. Genome sequencing has put pressure on the classification system and recent years have seen significant changes to phage taxonomy. Here, we reflect on the state of phage taxonomy and provide a roadmap for the future, including the abolition of the order Caudovirales and the families Myoviridae, Podoviridae, and Siphoviridae. Furthermore, we specify guidelines for the demarcation of species, genus, subfamily and family-level ranks of tailed phage taxonomy.

The human oral phageome

Oral bacteriophages (or phages), especially periodontal ones, constitute a growing area of interest, but research on oral phages is still in its infancy. Phages are bacterial viruses that may persist as intracellular parasitic deoxyribonucleic acid (DNA) or use bacterial metabolism to replicate and cause bacterial lysis. The microbiomes of saliva, oral mucosa, and dental plaque contain active phage virions, bacterial lysogens (ie, carrying dormant prophages), and bacterial strains containing short fragments of phage DNA. In excess of 2000 oral phages have been confirmed or predicted to infect species of the phyla Actinobacteria (>300 phages), Bacteroidetes (>300 phages), Firmicutes (>1000 phages), Fusobacteria (>200 phages), and Proteobacteria (>700 phages) and three additional phyla (few phages only). This article assesses the current knowledge of the diversity of the oral phage population and the mechanisms by which phages may impact the ecology of oral biofilms. The potential use of phage-based therapy to control major periodontal pathogens is also discussed.

VPF-Class: Taxonomic assignment and host prediction of uncultivated viruses based on viral protein families

Motivation

Two key steps in the analysis of uncultured viruses recovered from metagenomes are the taxonomic classification of the viral sequences and the identification of putative host(s). Both steps rely mainly on the assignment of viral proteins to orthologs in cultivated viruses. Viral Protein Families (VPFs) can be used for the robust identification of new viral sequences in large metagenomics datasets. Despite the importance of VPF information for viral discovery, VPFs have not yet been explored for determining viral taxonomy and host targets.

Results

In this work, we classified the set of VPFs from the IMG/VR database and developed VPF-Class. VPF-Class is a tool that automates the taxonomic classification and host prediction of viral contigs based on the assignment of their proteins to a set of classified VPFs. Applying VPF-Class on 731K uncultivated virus contigs from the IMG/VR database, we were able to classify 363K contigs at the genus level and predict the host of over 461K contigs. In the RefSeq database, VPF-class reported an accuracy of nearly 100% to classify dsDNA, ssDNA and retroviruses, at the genus level, considering a membership ratio and a confidence score of 0.2. The accuracy in host prediction was 86.4%, also at the genus level, considering a membership ratio of 0.3 and a confidence score of 0.5. And, in the prophages dataset, the accuracy in host prediction was 86% considering a membership ratio of 0.6 and a confidence score of 0.8. Moreover, from the Global Ocean Virome dataset, over 817K viral contigs out of 1 million were classified.

Availability and implementation

The implementation of VPF-Class can be downloaded from https://github.com/biocom-uib/vpf-tools.

Supplementary information

Supplementary data are available at Bioinformatics online.

Binomial nomenclature for virus species: a long view

On several occasions over the past century it has been proposed that Latinized (Linnaean) binomial names (LBs) should be used for the formal names of virus species, and the opinions expressed in the early debates are still valid. The use of LBs would be sensible for the current Taxonomy if confined to the names of the specific and generic taxa of viruses of which some basic biological properties are known (e.g. ecology, hosts and virions); there is no advantage in filling the literature with formal names for partly described viruses or virus-like gene sequences. The ICTV should support the time-honoured convention that LBs are only used with biological (phylogenetic) classifications. Recent changes have left the ICTV Taxonomy and its Code uncoordinated, and its aims and audience uncertain.

Structure and assembly of archaeal viruses

Viruses of archaea represent one of the most enigmatic parts of the virosphere. Most of the characterized archaeal viruses infect extremophilic hosts and display remarkable diversity of virion morphotypes, many of which have never been observed among bacteriophages or viruses of eukaryotes. However, recent environmental studies have shown that archaeal viruses are widespread also in moderate ecosystems, where they play an important ecological role by influencing the turnover of microbial communities, with a global impact on the carbon and nitrogen cycles. In this review, we summarize recent advances in understanding the molecular details of virion organization and assembly of archaeal viruses. We start by briefly introducing the 20 officially recognized families of archaeal viruses and then outline the similarities and differences of archaeal virus assembly with the morphogenesis pathways used by bacterial and eukaryotic viruses, and discuss the evolutionary implications of these observations. Generally, the assembly of the icosahedral archaeal viruses closely follows the mechanisms employed by evolutionarily related bacterial and eukaryotic viruses with the HK97 fold and double jelly-roll major capsid proteins, emphasizing the overall conservation of these pathways over billions of years of evolution. By contrast, archaea-specific viruses employ unique virion assembly mechanisms. We also highlight some of the molecular adaptations underlying the stability of archaeal viruses in extreme environments. Despite considerable progress during the past few years, the archaeal virosphere continues to represent one of the least studied parts of the global virome, with many molecular features awaiting to be discovered and characterized.

The Application of Porphyrins and Their Analogues for Inactivation of Viruses

The problem of treating viral infections is extremely relevant due to both the emergence of new viral diseases and to the low effectiveness of existing approaches to the treatment of known viral infections. This review focuses on the application of porphyrin, chlorin, and phthalocyanine series for combating viral infections by chemical and photochemical inactivation methods. The purpose of this review paper is to summarize the main approaches developed to date in the chemical and photodynamic inactivation of human and animal viruses using porphyrins and their analogues and to analyze and discuss the information on viral targets and antiviral activity of porphyrins, chlorins, of their conjugates with organic/inorganic compounds obtained in the last 10–15 years in order to identify the most promising areas.

Shining Light on Human Gut Bacteriophages

The human gut is a complex environment that contains a multitude of microorganisms that are collectively termed the microbiome. Multiple factors have a role to play in driving the composition of human gut bacterial communities either toward homeostasis or the instability that is associated with many disease states. One of the most important forces are likely to be bacteriophages, bacteria-infecting viruses that constitute by far the largest portion of the human gut virome. Despite this, bacteriophages (phages) are the one of the least studied residents of the gut. This is largely due to the challenges associated with studying these difficult to culture entities. Modern high throughput sequencing technologies have played an important role in improving our understanding of the human gut phageome but much of the generated sequencing data remains uncharacterised. Overcoming this requires database-independent bioinformatic pipelines and even those phages that are successfully characterized only provide limited insight into their associated biological properties, and thus most viral sequences have been characterized as “viral dark matter.” Fundamental to understanding the role of phages in shaping the human gut microbiome, and in turn perhaps influencing human health, is how they interact with their bacterial hosts. An essential aspect is the isolation of novel phage-bacteria host pairs by direct isolation through various screening methods, which can transform in silico phages into a biological reality. However, this is also beset with multiple challenges including culturing difficulties and the use of traditional methods, such as plaquing, which may bias which phage-host pairs that can be successfully isolated. Phage-bacteria interactions may be influenced by many aspects of complex human gut biology which can be difficult to reproduce under laboratory conditions. Here we discuss some of the main findings associated with the human gut phageome to date including composition, our understanding of phage-host interactions, particularly the observed persistence of virulent phages and their hosts, as well as factors that may influence these highly intricate relationships. We also discuss current methodologies and bottlenecks hindering progression in this field and identify potential steps that may be useful in overcoming these hurdles.

Exploring the Remarkable Diversity of Culturable Escherichia coli Phages in the Danish Wastewater Environment

Phages drive bacterial diversity, profoundly influencing microbial communities, from microbiomes to the drivers of global biogeochemical cycling. Aiming to broaden our understanding of Escherichiacoli (MG1655, K-12) phages, we screened 188 Danish wastewater samples and isolated 136 phages. Ninety-two of these have genomic sequences with less than 95% similarity to known phages, while most map to existing genera several represent novel lineages. The isolated phages are highly diverse, estimated to represent roughly one-third of the true diversity of culturable virulent dsDNA Escherichia phages in Danish wastewater, yet almost half (40%) are not represented in metagenomic databases, emphasising the importance of isolating phages to uncover diversity. Seven viral families, Myoviridae, Siphoviridae, Podoviridae,Drexlerviridae,Chaseviridae,Autographviridae, and Microviridae, are represented in the dataset. Their genomes vary drastically in length from 5.3 kb to 170.8 kb, with a guanine and cytosine (GC) content ranging from 35.3% to 60.0%. Hence, even for a model host bacterium, substantial diversity remains to be uncovered. These results expand and underline the range of coliphage diversity and demonstrate how far we are from fully disclosing phage diversity and ecology.

Genomic nucleotide-based distance analysis for delimiting old world monkey derived herpes simplex virus species

BACKGROUND

Herpes simplex viruses form a genus within the alphaherpesvirus subfamily, with three identified viral species isolated from Old World monkeys (OWM); Macacine alphaherpesvirus 1 (McHV-1; herpes B), Cercopithecine alphaherpesvirus 2 (SA8), and Papiine alphaherpesvirus 2 (PaHV-2; herpes papio). Herpes B is endemic to macaques, while PaHV-2 and SA8 appear endemic to baboons. All three viruses are genetically and antigenically similar, with SA8 and PaHV-2 thought to be avirulent in humans, while herpes B is a biosafety level 4 pathogen. Recently, next-generation sequencing (NGS) has resulted in an increased number of published OWM herpes simplex genomes, allowing an encompassing phylogenetic analysis.

RESULTS

In this study, phylogenetic networks, in conjunction with a genome-based genetic distance cutoff method were used to examine 27 OWM monkey herpes simplex isolates. Genome-based genetic distances were calculated, resulting in distances between lion and pig-tailed simplex viruses themselves, and versus herpes B core strains that were higher than those between PaHV-2 and SA8 (approximately 14 and 10% respectively). The species distance cutoff was determined to be 8.94%, with the method recovering separate species status for PaHV-2 and SA8 and showed that lion and pig-tailed simplex viruses (vs core herpes B strains) were well over the distance species cutoff.

CONCLUSIONS

We propose designating lion and pig-tailed simplex viruses as separate, individual viral species, and that this may be the first identification of viral cryptic species.

Domain-centric database to uncover structure of minimally characterized viral genomes

Protein domain-based approaches to analyzing sequence data are valuable tools for examining and exploring genomic architecture across genomes of different organisms. Here, we present a complete dataset of domains from the publicly available sequence data of 9,051 reference viral genomes. The data provided contain information such as sequence position and neighboring domains from 30,947 pHMM-identified domains from each reference viral genome. Domains were identified from viral whole-genome sequence using automated profile Hidden Markov Models (pHMM). This study also describes the framework for constructing "domain neighborhoods", as well as the dataset representing it. These data can be used to examine shared and differing domain architectures across viral genomes, to elucidate potential functional properties of genes, and potentially to classify viruses.

Global Organization and Proposed Megataxonomy of the Virus World

Viruses and mobile genetic elements are molecular parasites or symbionts that coevolve with nearly all forms of cellular life. The route of virus replication and protein expression is determined by the viral genome type. Comparison of these routes led to the classification of viruses into seven "Baltimore classes" (BCs) that define the major features of virus reproduction. However, recent phylogenomic studies identified multiple evolutionary connections among viruses within each of the BCs as well as between different classes. Due to the modular organization of virus genomes, these relationships defy simple representation as lines of descent but rather form complex networks. Phylogenetic analyses of virus hallmark genes combined with analyses of gene-sharing networks show that replication modules of five BCs (three classes of RNA viruses and two classes of reverse-transcribing viruses) evolved from a common ancestor that encoded an RNA-directed RNA polymerase or a reverse transcriptase. Bona fide viruses evolved from this ancestor on multiple, independent occasions via the recruitment of distinct cellular proteins as capsid subunits and other structural components of virions. The single-stranded DNA (ssDNA) viruses are a polyphyletic class, with different groups evolving by recombination between rolling-circle-replicating plasmids, which contributed the replication protein, and positive-sense RNA viruses, which contributed the capsid protein. The double-stranded DNA (dsDNA) viruses are distributed among several large monophyletic groups and arose via the combination of distinct structural modules with equally diverse replication modules. Phylogenomic analyses reveal the finer structure of evolutionary connections among RNA viruses and reverse-transcribing viruses, ssDNA viruses, and large subsets of dsDNA viruses. Taken together, these analyses allow us to outline the global organization of the virus world. Here, we describe the key aspects of this organization and propose a comprehensive hierarchical taxonomy of viruses.

Enteric Virome and Carcinogenesis in the Gut

Colorectal cancer (CRC) is a leading cause of cancer-related deaths in both the USA and the world. Recent research has demonstrated the involvement of the gut microbiota in CRC development and progression. Microbial biomarkers of disease have focused primarily on the bacterial component of the microbiome; however, the viral portion of the microbiome, consisting of both bacteriophages and eukaryotic viruses, together known as the virome, has been lesser studied. Here we review the recent advancements in high-throughput sequencing (HTS) technologies and bioinformatics, which have enabled scientists to better understand how viruses might influence the development of colorectal cancer. We discuss the contemporary findings revealing modulations in the virome and their correlation with CRC development and progression. While a variety of challenges still face viral HTS detection in clinical specimens, we consider herein numerous next steps for future basic and clinical research. Clinicians need to move away from a single infectious agent model for disease etiology by grasping new, more encompassing etiological paradigms, in which communities of various microbial components interact with each other and the host. The reporting and indexing of patient health information, socioeconomic data, and other relevant metadata will enable identification of predictive variables and covariates of viral presence and CRC development. Altogether, the virome has a more profound role in carcinogenesis and cancer progression than once thought, and viruses, specific for either human cells or bacteria, are clinically relevant in understanding CRC pathology, patient prognosis, and treatment development.

Iroki: automatic customization and visualization of phylogenetic trees

Phylogenetic trees are an important analytical tool for evaluating community diversity and evolutionary history. In the case of microorganisms, the decreasing cost of sequencing has enabled researchers to generate ever-larger sequence datasets, which in turn have begun to fill gaps in the evolutionary history of microbial groups. However, phylogenetic analyses of these types of datasets create complex trees that can be challenging to interpret. Scientific inferences made by visual inspection of phylogenetic trees can be simplified and enhanced by customizing various parts of the tree. Yet, manual customization is time-consuming and error prone, and programs designed to assist in batch tree customization often require programming experience or complicated file formats for annotation. Iroki, a user-friendly web interface for tree visualization, addresses these issues by providing automatic customization of large trees based on metadata contained in tab-separated text files. Iroki’s utility for exploring biological and ecological trends in sequencing data was demonstrated through a variety of microbial ecology applications in which trees with hundreds to thousands of leaf nodes were customized according to extensive collections of metadata. The Iroki web application and documentation are available at https://www.iroki.net or through the VIROME portal http://virome.dbi.udel.edu. Iroki’s source code is released under the MIT license and is available at https://github.com/mooreryan/iroki.

Analysis of Spounaviruses as a Case Study for the Overdue Reclassification of Tailed Phages

Tailed bacteriophages are the most abundant and diverse viruses in the world, with genome sizes ranging from 10 kbp to over 500 kbp. Yet, due to historical reasons, all this diversity is confined to a single virus order-Caudovirales, composed of just four families: Myoviridae, Siphoviridae, Podoviridae, and the newly created Ackermannviridae family. In recent years, this morphology-based classification scheme has started to crumble under the constant flood of phage sequences, revealing that tailed phages are even more genetically diverse than once thought. This prompted us, the Bacterial and Archaeal Viruses Subcommittee of the International Committee on Taxonomy of Viruses (ICTV), to consider overall reorganization of phage taxonomy. In this study, we used a wide range of complementary methods-including comparative genomics, core genome analysis, and marker gene phylogenetics-to show that the group of Bacillus phage SPO1-related viruses previously classified into the Spounavirinae subfamily, is clearly distinct from other members of the family Myoviridae and its diversity deserves the rank of an autonomous family. Thus, we removed this group from the Myoviridae family and created the family Herelleviridae-a new taxon of the same rank. In the process of the taxon evaluation, we explored the feasibility of different demarcation criteria and critically evaluated the usefulness of our methods for phage classification. The convergence of results, drawing a consistent and comprehensive picture of a new family with associated subfamilies, regardless of method, demonstrates that the tools applied here are particularly useful in phage taxonomy. We are convinced that creation of this novel family is a crucial milestone toward much-needed reclassification in the Caudovirales order.

Charting the diversity of uncultured viruses of Archaea and Bacteria

BACKGROUND

Viruses of Archaea and Bacteria are among the most abundant and diverse biological entities on Earth. Unraveling their biodiversity has been challenging due to methodological limitations. Recent advances in culture-independent techniques, such as metagenomics, shed light on the unknown viral diversity, revealing thousands of new viral nucleotide sequences at an unprecedented scale. However, these novel sequences have not been properly classified and the evolutionary associations between them were not resolved.

RESULTS

Here, we performed phylogenomic analysis of nearly 200,000 viral nucleotide sequences to establish GL-UVAB: Genomic Lineages of Uncultured Viruses of Archaea and Bacteria. The pan-genome content of the identified lineages shed light on some of their infection strategies, potential to modulate host physiology, and mechanisms to escape host resistance systems. Furthermore, using GL-UVAB as a reference database for annotating metagenomes revealed elusive habitat distribution patterns of viral lineages and environmental drivers of community composition.

CONCLUSIONS

These findings provide insights about the genomic diversity and ecology of viruses of prokaryotes. The source code used in these analyses is freely available at https://sourceforge.net/projects/gluvab/.

Charting the diversity of uncultured viruses of Archaea and Bacteria

Background

Viruses of Archaea and Bacteria are among the most abundant and diverse biological entities on Earth. Unraveling their biodiversity has been challenging due to methodological limitations. Recent advances in culture-independent techniques, such as metagenomics, shed light on the unknown viral diversity, revealing thousands of new viral nucleotide sequences at an unprecedented scale. However, these novel sequences have not been properly classified and the evolutionary associations between them were not resolved.

Results

Here, we performed phylogenomic analysis of nearly 200,000 viral nucleotide sequences to establish GL-UVAB: Genomic Lineages of Uncultured Viruses of Archaea and Bacteria. The pan-genome content of the identified lineages shed light on some of their infection strategies, potential to modulate host physiology, and mechanisms to escape host resistance systems. Furthermore, using GL-UVAB as a reference database for annotating metagenomes revealed elusive habitat distribution patterns of viral lineages and environmental drivers of community composition.

Conclusions

These findings provide insights about the genomic diversity and ecology of viruses of prokaryotes. The source code used in these analyses is freely available at https://sourceforge.net/projects/gluvab/.

Re-assessing the diversity of negative strand RNA viruses in insects

The spectrum of viruses in insects is important for subjects as diverse as public health, veterinary medicine, food production, and biodiversity conservation. The traditional interest in vector-borne diseases of humans and livestock has drawn the attention of virus studies to hematophagous insect species. However, these represent only a tiny fraction of the broad diversity of Hexapoda, the most speciose group of animals. Here, we systematically probed the diversity of negative strand RNA viruses in the largest and most representative collection of insect transcriptomes from samples representing all 34 extant orders of Hexapoda and 3 orders of Entognatha, as well as outgroups, altogether representing 1243 species. Based on profile hidden Markov models we detected 488 viral RNA-directed RNA polymerase (RdRp) sequences with similarity to negative strand RNA viruses. These were identified in members of 324 arthropod species. Selection for length, quality, and uniqueness left 234 sequences for analyses, showing similarity to genomes of viruses classified in Bunyavirales (n = 86), Articulavirales (n = 54), and several orders within Haploviricotina (n = 94). Coding-complete genomes or nearly-complete subgenomic assemblies were obtained in 61 cases. Based on phylogenetic topology and the availability of coding-complete genomes we estimate that at least 20 novel viral genera in seven families need to be defined, only two of them monospecific. Seven additional viral clades emerge when adding sequences from the present study to formerly monospecific lineages, potentially requiring up to seven additional genera. One long sequence may indicate a novel family. For segmented viruses, cophylogenies between genome segments were generally improved by the inclusion of viruses from the present study, suggesting that in silico misassembly of segmented genomes is rare or absent. Contrary to previous assessments, significant virus-host codivergence was identified in major phylogenetic lineages based on two different approaches of codivergence analysis in a hypotheses testing framework. In spite of these additions to the known spectrum of viruses in insects, we caution that basing taxonomic decisions on genome information alone is challenging due to technical uncertainties, such as the inability to prove integrity of complete genome assemblies of segmented viruses.

Gut Bacteriophage: Current Understanding and Challenges

The gut microbiome is widely accepted to have a significant impact on human health yet, despite years of research on this complex ecosystem, the contributions of different forces driving microbial population structure remain to be fully elucidated. The viral component of the human gut microbiome is dominated by bacteriophage, which are known to play crucial roles in shaping microbial composition, driving bacterial diversity, and facilitating horizontal gene transfer. Bacteriophage are also one of the most poorly understood components of the human gut microbiome, with the vast majority of viral sequences sharing little to no homology to reference databases. If we are to understand the dynamics of bacteriophage populations, their interaction with the human microbiome and ultimately their influence on human health, we will depend heavily on sequence based approaches and in silico tools. This is complicated by the fact that, as with any research field in its infancy, methods of analyses vary and this can impede our ability to compare the outputs of different studies. Here, we discuss the major findings to date regarding the human virome and reflect on our current understanding of how gut bacteriophage shape the microbiome. We consider whether or not the virome field is built on unstable foundations and if so, how can we provide a solid basis for future experimentation. The virome is a challenging yet crucial piece of the human microbiome puzzle. In order to develop our understanding, we will discuss the need to underpin future studies with robust research methods and suggest some solutions to existing challenges.