Ecogenomics and potential biogeochemical impacts of globally abundant ocean viruses

RaFAH: Host prediction for viruses of Bacteria and Archaea based on protein content

Culture-independent approaches have recently shed light on the genomic diversity of viruses of prokaryotes. One fundamental question when trying to understand their ecological roles is: which host do they infect? To tackle this issue we developed a machine-learning approach named Random Forest Assignment of Hosts (RaFAH), that uses scores to 43,644 protein clusters to assign hosts to complete or fragmented genomes of viruses of Archaea and Bacteria. RaFAH displayed performance comparable with that of other methods for virus-host prediction in three different benchmarks encompassing viruses from RefSeq, single amplified genomes, and metagenomes. RaFAH was applied to assembled metagenomic datasets of uncultured viruses from eight different biomes of medical, biotechnological, and environmental relevance. Our analyses led to the identification of 537 sequences of archaeal viruses representing unknown lineages, whose genomes encode novel auxiliary metabolic genes, shedding light on how these viruses interfere with the host molecular machinery. RaFAH is available at https://sourceforge.net/projects/rafah/.

•

RaFAH was developed to predict the hosts of viruses of Bacteria and Archaea

•

RaFAH displayed comparable or superior performance to other host-prediction tools

•

RaFAH performed well across viromes from eight different ecosystems

•

RaFAH identified hundreds of genomic sequences as derived from viruses of Archaea

Viruses that infect Bacteria and Archaea are ubiquitous and extremely abundant. Recent advances have led to the discovery of many thousands of complete and partial genomes of these biological entities. Understanding the biology of these viruses and how they influence their ecosystems depends on knowing which hosts they infect. We developed a tool that uses data from complete or fragmented genomes to predict the hosts of viruses using a machine-learning approach. Our tool, RaFAH, displayed performance comparable with or superior to that of other host-prediction tools. In addition, it identified hundreds of sequences as derived from the genomes of viruses of Archaea, which are one of the least characterized fractions of the global virosphere.

Fishing for phages in metagenomes: what do we catch, what do we miss?

Metagenomics and metatranscriptomics have become the principal approaches for discovery of novel bacteriophages and preliminary characterization of their ecology and biology. Metagenomic sequencing dramatically expanded the known diversity of tailed and non-tailed phages with double-stranded DNA genomes and those with single-stranded DNA genomes, whereas metatranscriptomics led to the discovery of thousands of new single-stranded RNA phages. Apart from expanding phage diversity, metagenomics studies discover major novel groups of phages with unique features of genome organization, expression strategy and virus-host interaction, such as the putative order 'crAssvirales', which includes the most abundant human-associated viruses. The continued success of metagenomics hinges on the combination of the most powerful computational methods for phage genome assembly and analysis including harnessing CRISPR spacers for the discovery of novel phages and host assignment. Together, these approaches could make a comprehensive characterization of the earth phageome a realistic goal.

Expanding standards in viromics: in silico evaluation of dsDNA viral genome identification, classification, and auxiliary metabolic gene curation

BACKGROUND

Viruses influence global patterns of microbial diversity and nutrient cycles. Though viral metagenomics (viromics), specifically targeting dsDNA viruses, has been critical for revealing viral roles across diverse ecosystems, its analyses differ in many ways from those used for microbes. To date, viromics benchmarking has covered read pre-processing, assembly, relative abundance, read mapping thresholds and diversity estimation, but other steps would benefit from benchmarking and standardization. Here we use in silico-generated datasets and an extensive literature survey to evaluate and highlight how dataset composition (i.e., viromes vs bulk metagenomes) and assembly fragmentation impact (i) viral contig identification tool, (ii) virus taxonomic classification, and (iii) identification and curation of auxiliary metabolic genes (AMGs).

RESULTS

The in silico benchmarking of five commonly used virus identification tools show that gene-content-based tools consistently performed well for long (≥3 kbp) contigs, while k-mer- and blast-based tools were uniquely able to detect viruses from short (≤3 kbp) contigs. Notably, however, the performance increase of k-mer- and blast-based tools for short contigs was obtained at the cost of increased false positives (sometimes up to ∼5% for virome and ∼75% bulk samples), particularly when eukaryotic or mobile genetic element sequences were included in the test datasets. For viral classification, variously sized genome fragments were assessed using gene-sharing network analytics to quantify drop-offs in taxonomic assignments, which revealed correct assignations ranging from ∼95% (whole genomes) down to ∼80% (3 kbp sized genome fragments). A similar trend was also observed for other viral classification tools such as VPF-class, ViPTree and VIRIDIC, suggesting that caution is warranted when classifying short genome fragments and not full genomes. Finally, we highlight how fragmented assemblies can lead to erroneous identification of AMGs and outline a best-practices workflow to curate candidate AMGs in viral genomes assembled from metagenomes.

CONCLUSION

Together, these benchmarking experiments and annotation guidelines should aid researchers seeking to best detect, classify, and characterize the myriad viruses 'hidden' in diverse sequence datasets.

Global overview and major challenges of host prediction methods for uncultivated phages

Bacterial communities play critical roles across all of Earth's biomes, affecting human health and global ecosystem functioning. They do so under strong constraints exerted by viruses, that is, bacteriophages or 'phages'. Phages can reshape bacterial communities' structure, influence long-term evolution of bacterial populations, and alter host cell metabolism during infection. Metagenomics approaches, that is, shotgun sequencing of environmental DNA or RNA, recently enabled large-scale exploration of phage genomic diversity, yielding several millions of phage genomes now to be further analyzed and characterized. One major challenge however is the lack of direct host information for these phages. Several methods and tools have been proposed to bioinformatically predict the potential host(s) of uncultivated phages based only on genome sequence information. Here we review these different approaches and highlight their distinct strengths and limitations. We also outline complementary experimental assays which are being proposed to validate and refine these bioinformatic predictions.

Ecology of inorganic sulfur auxiliary metabolism in widespread bacteriophages

Microbial sulfur metabolism contributes to biogeochemical cycling on global scales. Sulfur metabolizing microbes are infected by phages that can encode auxiliary metabolic genes (AMGs) to alter sulfur metabolism within host cells but remain poorly characterized. Here we identified 191 phages derived from twelve environments that encoded 227 AMGs for oxidation of sulfur and thiosulfate (dsrA, dsrC/tusE, soxC, soxD and soxYZ). Evidence for retention of AMGs during niche-differentiation of diverse phage populations provided evidence that auxiliary metabolism imparts measurable fitness benefits to phages with ramifications for ecosystem biogeochemistry. Gene abundance and expression profiles of AMGs suggested significant contributions by phages to sulfur and thiosulfate oxidation in freshwater lakes and oceans, and a sensitive response to changing sulfur concentrations in hydrothermal environments. Overall, our study provides fundamental insights on the distribution, diversity, and ecology of phage auxiliary metabolism associated with sulfur and reinforces the necessity of incorporating viral contributions into biogeochemical configurations.

Exploring Viral Diversity in a Gypsum Karst Lake Ecosystem Using Targeted Single-Cell Genomics

Little is known about the diversity and distribution of viruses infecting green sulfur bacteria (GSB) thriving in euxinic (sulfuric and anoxic) habitats, including gypsum karst lake ecosystems. In this study, we used targeted cell sorting combined with single-cell sequencing to gain insights into the gene content and genomic potential of viruses infecting sulfur-oxidizing bacteria Chlorobium clathratiforme, obtained from water samples collected during summer stratification in gypsum karst Lake Kirkilai (Lithuania). In total, 82 viral contigs were bioinformatically identified in 62 single amplified genomes (SAGs) of C. clathratiforme. The majority of viral gene and protein sequences showed little to no similarity with phage sequences in public databases, uncovering the vast diversity of previously undescribed GSB viruses. We observed a high level of lysogenization in the C. clathratiforme population, as 87% SAGs contained intact prophages. Among the thirty identified auxiliary metabolic genes (AMGs), two, thiosulfate sulfurtransferase (TST) and thioredoxin-dependent phosphoadenosine phosphosulfate (PAPS) reductase (cysH), were found to be involved in the oxidation of inorganic sulfur compounds, suggesting that viruses can influence the metabolism and cycling of this essential element. Finally, the analysis of CRISPR spacers retrieved from the consensus C. clathratiforme genome imply persistent and active virus–host interactions for several putative phages prevalent among C. clathratiforme SAGs. Overall, this study provides a glimpse into the diversity of phages associated with naturally occurring and highly abundant sulfur-oxidizing bacteria.

Temporal Changes of Virus-Like Particle Abundance and Metagenomic Comparison of Viral Communities in Cropland and Prairie Soils

During the last several decades, viruses have been increasingly recognized for their abundance, ubiquity, and important roles in different ecosystems. Despite known contributions to aquatic systems, few studies examine viral abundance and community structure over time in terrestrial ecosystems. The effects of land conversion and land management on soil microbes have been previously investigated, but their effects on virus population are not well studied. This study examined annual dynamics of viral abundance in soils from a native tallgrass prairie and two croplands, conventional till winter wheat and no-till canola, in Oklahoma. Virus-like particle (VLP) abundance varied across sites, and showed clear seasonal shifts. VLP abundance significantly correlated with environmental variables that were generally reflective of land use, including air temperature, soil nitrogen, and plant canopy coverage. Structural equation modeling supported the effects of land use on soil communities by emphasizing interactions between management, environmental factors, and viral and bacterial abundance. Between the viral metagenomes from the prairie and tilled wheat field, 1,231 unique viral operational taxonomic units (vOTUs) were identified, and only five were shared that were rare in the contrasting field. Only 13% of the vOTUs had similarity to previously identified viruses in the RefSeq database, with only 7% having known taxonomic classification. Together, our findings indicated land use and tillage practices influence virus abundance and community structure. Analyses of viromes over time and space are vital to viral ecology in providing insight on viral communities and key information on interactions between viruses, their microbial hosts, and the environment.

IMPORTANCE Conversion of land alters the physiochemical and biological environments by not only changing the aboveground community, but also modifying the soil environment for viruses and microbes. Soil microbial communities are critical to nutrient cycling, carbon mineralization, and soil quality; and viruses are known for influencing microbial abundance, community structure, and evolution. Therefore, viruses are considered an important part of soil functions in terrestrial ecosystems. In aquatic environments, virus abundance generally exceeds bacterial counts by an order of magnitude, and they are thought to be one of the greatest genetic reservoirs on the planet. However, data are extremely limited on viruses in soils, and even less is known about their responses to the disturbances associated with land use and management. The study provides important insights into the temporal dynamics of viral abundance and the structure of viral communities in response to the common practice of turning native habitats into arable soils.

Host population diversity as a driver of viral infection cycle in wild populations of green sulfur bacteria with long standing virus-host interactions

Temperate phages are viruses of bacteria that can establish two types of infection: a lysogenic infection in which the virus replicates with the host cell without producing virions, and a lytic infection where the host cell is eventually destroyed, and new virions are released. While both lytic and lysogenic infections are routinely observed in the environment, the ecological and evolutionary processes regulating these viral dynamics are still not well understood, especially for uncultivated virus-host pairs. Here, we characterized the long-term dynamics of uncultivated viruses infecting green sulfur bacteria (GSB) in a model freshwater lake (Trout Bog Lake, TBL). As no GSB virus has been formally described yet, we first used two complementary approaches to identify new GSB viruses from TBL; one in vitro based on flow cytometry cell sorting, the other in silico based on CRISPR spacer sequences. We then took advantage of existing TBL metagenomes covering the 2005-2018 period to examine the interactions between GSB and their viruses across years and seasons. From our data, GSB populations in TBL were constantly associated with at least 2-8 viruses each, including both lytic and temperate phages. The dominant GSB population in particular was consistently associated with two prophages with a nearly 100% infection rate for >10 years. We illustrate with a theoretical model that such an interaction can be stable given a low, but persistent, level of prophage induction in low-diversity host populations. Overall, our data suggest that lytic and lysogenic viruses can readily co-infect the same host population, and that host strain-level diversity might be an important factor controlling virus-host dynamics including lytic/lysogeny switch.

The bacterial sulfur cycle in expanding dysoxic and euxinic marine waters

Dysoxic marine waters (DMW, < 1 μM oxygen) are currently expanding in volume in the oceans, which has biogeochemical, ecological and societal consequences on a global scale. In these environments, distinct bacteria drive an active sulfur cycle, which has only recently been recognized for open-ocean DMW. This review summarizes the current knowledge on these sulfur-cycling bacteria. Critical bottlenecks and questions for future research are specifically addressed. Sulfate-reducing bacteria (SRB) are core members of DMW. However, their roles are not entirely clear, and they remain largely uncultured. We found support for their remarkable diversity and taxonomic novelty by mining metagenome-assembled genomes from the Black Sea as model ecosystem. We highlight recent insights into the metabolism of key sulfur-oxidizing SUP05 and Sulfurimonas bacteria, and discuss the probable involvement of uncultivated SAR324 and BS-GSO2 bacteria in sulfur oxidation. Uncultivated Marinimicrobia bacteria with a presumed organoheterotrophic metabolism are abundant in DMW. Like SRB, they may use specific molybdoenzymes to conserve energy from the oxidation, reduction or disproportionation of sulfur cycle intermediates such as S0 and thiosulfate, produced from the oxidation of sulfide. We expect that tailored sampling methods and a renewed focus on cultivation will yield deeper insight into sulfur-cycling bacteria in DMW.

Efficient dilution-to-extinction isolation of novel virus-host model systems for fastidious heterotrophic bacteria

Microbes and their associated viruses are key drivers of biogeochemical processes in marine and soil biomes. While viruses of phototrophic cyanobacteria are well-represented in model systems, challenges of isolating marine microbial heterotrophs and their viruses have hampered experimental approaches to quantify the importance of viruses in nutrient recycling. A resurgence in cultivation efforts has improved the availability of fastidious bacteria for hypothesis testing, but this has not been matched by similar efforts to cultivate their associated bacteriophages. Here, we describe a high-throughput method for isolating important virus-host systems for fastidious heterotrophic bacteria that couples advances in culturing of hosts with sequential enrichment and isolation of associated phages. Applied to six monthly samples from the Western English Channel, we first isolated one new member of the globally dominant bacterial SAR11 clade and three new members of the methylotrophic bacterial clade OM43. We used these as bait to isolate 117 new phages, including the first known siphophage-infecting SAR11, and the first isolated phage for OM43. Genomic analyses of 13 novel viruses revealed representatives of three new viral genera, and infection assays showed that the viruses infecting SAR11 have ecotype-specific host ranges. Similar to the abundant human-associated phage ɸCrAss001, infection dynamics within the majority of isolates suggested either prevalent lysogeny or chronic infection, despite a lack of associated genes, or host phenotypic bistability with lysis putatively maintained within a susceptible subpopulation. Broader representation of important virus-host systems in culture collections and genomic databases will improve both our understanding of virus-host interactions, and accuracy of computational approaches to evaluate ecological patterns from metagenomic data.

A Novel Broad Host Range Phage Infecting Alteromonas

Bacteriophages substantially contribute to bacterial mortality in the ocean and play critical roles in global biogeochemical processes. Alteromonas is a ubiquitous bacterial genus in global tropical and temperate waters, which can cross-protect marine cyanobacteria and thus has important ecological benefits. However, little is known about the biological and ecological features of Alteromonas phages (alterophages). Here, we describe a novel alterophage vB_AmeP-R8W (R8W), which belongs to the Autographiviridae family and infects the deep-clade Alteromonas mediterranea. R8W has an equidistant and icosahedral head (65 ± 1 nm in diameter) and a short tail (12 ± 2 nm in length). The genome size of R8W is 48,825 bp, with a G + C content of 40.55%. R8W possesses three putative auxiliary metabolic genes encoding proteins involved in nucleotide metabolism and DNA binding: thymidylate synthase, nucleoside triphosphate pyrophosphohydrolase, and PhoB. R8W has a rapid lytic cycle with a burst size of 88 plaque-forming units/cell. Notably, R8W has a wide host range, such that it can infect 35 Alteromonas strains; it exhibits a strong specificity for strains isolated from deep waters. R8W has two specific receptor binding proteins and a compatible holin-endolysin system, which contribute to its wide host range. The isolation of R8W will contribute to the understanding of alterophage evolution, as well as the phage-host interactions and ecological importance of alterophages.

Ecology and molecular targets of hypermutation in the global microbiome

Changes in the sequence of an organism's genome, i.e., mutations, are the raw material of evolution. The frequency and location of mutations can be constrained by specific molecular mechanisms, such as diversity-generating retroelements (DGRs). DGRs have been characterized from cultivated bacteria and bacteriophages, and perform error-prone reverse transcription leading to mutations being introduced in specific target genes. DGR loci were also identified in several metagenomes, but the ecological roles and evolutionary drivers of these DGRs remain poorly understood. Here, we analyze a dataset of >30,000 DGRs from public metagenomes, establish six major lineages of DGRs including three primarily encoded by phages and seemingly used to diversify host attachment proteins, and demonstrate that DGRs are broadly active and responsible for >10% of all amino acid changes in some organisms. Overall, these results highlight the constraints under which DGRs evolve, and elucidate several distinct roles these elements play in natural communities.

Prophage Genomics and Ecology in the Family Rhodobacteraceae

Roseobacters are globally abundant bacteria with critical roles in carbon and sulfur biogeochemical cycling. Here, we identified 173 new putative prophages in 79 genomes of Rhodobacteraceae. These prophages represented 1.3 ± 0.15% of the bacterial genomes and had no to low homology with reference and metagenome-assembled viral genomes from aquatic and terrestrial ecosystems. Among the newly identified putative prophages, 35% encoded auxiliary metabolic genes (AMGs), mostly involved in secondary metabolism, amino acid metabolism, and cofactor and vitamin production. The analysis of integration sites and gene homology showed that 22 of the putative prophages were actually gene transfer agents (GTAs) similar to a GTA of Rhodobacter capsulatus. Twenty-three percent of the predicted prophages were observed in the TARA Oceans viromes generated from free viral particles, suggesting that they represent active prophages capable of induction. The distribution of these prophages was significantly associated with latitude and temperature. The prophages most abundant at high latitudes encoded acpP, an auxiliary metabolic gene involved in lipid synthesis and membrane fluidity at low temperatures. Our results show that prophages and gene transfer agents are significant sources of genomic diversity in roseobacter, with potential roles in the ecology of this globally distributed bacterial group.

Deep ocean metagenomes provide insight into the metabolic architecture of bathypelagic microbial communities

The deep sea, the largest ocean's compartment, drives planetary-scale biogeochemical cycling. Yet, the functional exploration of its microbial communities lags far behind other environments. Here we analyze 58 metagenomes from tropical and subtropical deep oceans to generate the Malaspina Gene Database. Free-living or particle-attached lifestyles drive functional differences in bathypelagic prokaryotic communities, regardless of their biogeography. Ammonia and CO oxidation pathways are enriched in the free-living microbial communities and dissimilatory nitrate reduction to ammonium and H2 oxidation pathways in the particle-attached, while the Calvin Benson-Bassham cycle is the most prevalent inorganic carbon fixation pathway in both size fractions. Reconstruction of the Malaspina Deep Metagenome-Assembled Genomes reveals unique non-cyanobacterial diazotrophic bacteria and chemolithoautotrophic prokaryotes. The widespread potential to grow both autotrophically and heterotrophically suggests that mixotrophy is an ecologically relevant trait in the deep ocean. These results expand our understanding of the functional microbial structure and metabolic capabilities of the largest Earth aquatic ecosystem.

Bacteriophage-Host Association in the Phytoplasma Insect Vector Euscelidius variegatus

Insect vectors transmit viruses and bacteria that can cause severe diseases in plants and economic losses due to a decrease in crop production. Insect vectors, like all other organisms, are colonized by a community of various microorganisms, which can influence their physiology, ecology, evolution, and also their competence as vectors. The important ecological meaning of bacteriophages in various ecosystems and their role in microbial communities has emerged in the past decade. However, only a few phages have been described so far in insect microbiomes. The leafhopper Euscelidius variegatus is a laboratory vector of the phytoplasma causing Flavescence dorée, a severe grapevine disease that threatens viticulture in Europe. Here, the presence of a temperate bacteriophage in E. variegatus (named Euscelidius variegatus phage 1, EVP-1) was revealed through both insect transcriptome analyses and electron microscopic observations. The bacterial host was isolated in axenic culture and identified as the bacterial endosymbiont of E. variegatus (BEV), recently assigned to the genus Candidatus Symbiopectobacterium. BEV harbors multiple prophages that become active in culture, suggesting that different environments can trigger different mechanisms, finely regulating the interactions among phages. Understanding the complex relationships within insect vector microbiomes may help in revealing possible microbe influences on pathogen transmission, and it is a crucial step toward innovative sustainable strategies for disease management in agriculture.

Spinal Cord Injury Changes the Structure and Functional Potential of Gut Bacterial and Viral Communities

Emerging data indicate that gut dysbiosis contributes to many human diseases, including several comorbidities that develop after traumatic spinal cord injury (SCI). To date, all analyses of SCI-induced gut dysbiosis have used 16S rRNA amplicon sequencing. This technique has several limitations, including being susceptible to taxonomic "blind spots," primer bias, and an inability to profile microbiota functions or identify viruses. Here, SCI-induced gut dysbiosis was assessed by applying genome- and gene-resolved metagenomic analysis of murine stool samples collected 21 days after an experimental SCI at the 4th thoracic spine (T4) or 10th thoracic spine (T10) spinal level. These distinct injuries partially (T10) or completely (T4) abolish sympathetic tone in the gut. Among bacteria, 105 medium- to high-quality metagenome-assembled genomes (MAGs) were recovered, with most (n = 96) representing new bacterial species. Read mapping revealed that after SCI, the relative abundance of beneficial commensals (Lactobacillus johnsonii and CAG-1031 spp.) decreased, while potentially pathogenic bacteria (Weissella cibaria, Lactococcus lactis _A, Bacteroides thetaiotaomicron) increased. Functionally, microbial genes encoding proteins for tryptophan, vitamin B6, and folate biosynthesis, essential pathways for central nervous system function, were reduced after SCI. Among viruses, 1,028 mostly novel viral populations were recovered, expanding known murine gut viral species sequence space ∼3-fold compared to that of public databases. Phages of beneficial commensal hosts (CAG-1031, Lactobacillus, and Turicibacter) decreased, while phages of pathogenic hosts (Weissella, Lactococcus, and class Clostridia) increased after SCI. Although the microbiomes and viromes were changed in all SCI mice, some of these changes varied as a function of spinal injury level, implicating loss of sympathetic tone as a mechanism underlying gut dysbiosis.IMPORTANCE To our knowledge, this is the first article to apply metagenomics to characterize changes in gut microbial population dynamics caused by a clinically relevant model of central nervous system (CNS) trauma. It also utilizes the most current approaches in genome-resolved metagenomics and viromics to maximize the biological inferences that can be made from these data. Overall, this article highlights the importance of autonomic nervous system regulation of a distal organ (gut) and its microbiome inhabitants after traumatic spinal cord injury (SCI). By providing information on taxonomy, function, and viruses, metagenomic data may better predict how SCI-induced gut dysbiosis influences systemic and neurological outcomes after SCI.

Analysis of Different Size Fractions Provides a More Complete Perspective of Viral Diversity in a Freshwater Embayment

Inspired by recent discoveries of the prevalence of large viruses in the environment, we reassessed the longstanding approach of filtering water through small-pore-size filters to separate viruses from cells before metagenomic analysis. We collected samples from three sites in Hamilton Harbour, an embayment of Lake Ontario, and studied 6 data sets derived from <0.45-μm- and >0.45-μm-size fractions to compare the diversity of viruses in these fractions. At the level of virus order/family, we observed highly diverse and distinct virus communities in the >0.45-μm-size fractions, whereas the <0.45-μm-size fractions were composed primarily of Caudovirales The relative abundances of Caudovirales for which hosts could be inferred varied widely between size fractions, with higher relative abundances of cyanophages in the >0.45-μm-size fractions, potentially indicating replication within cells during ongoing infections. Many viruses of eukaryotes, such as Mimiviridae, Phycodnaviridae, Iridoviridae, and Poxviridae, were detected exclusively in the often-disregarded >0.45-μm-size fractions. In addition to observing unique virus communities associated with each size fraction from every site we examined, we detected viruses common to both fractions, suggesting that these are candidates for further exploration because they could be the product of ongoing or recent lytic events. Most importantly, our observations indicate that analysis of either fraction alone provides only a partial perspective of double-stranded DNA (dsDNA) viruses in the environment, highlighting the need for more comprehensive approaches for analyzing virus communities inferred from metagenomic sequencing.IMPORTANCE Most studies of aquatic virus communities analyze DNA sequences derived from the smaller-size "free-virus" fraction. Our study demonstrates that analysis of virus communities using only the smaller-size fraction can lead to erroneously low diversity estimates for many of the larger viruses such as Mimiviridae, Phycodnaviridae, Iridoviridae, and Poxviridae, whereas analyzing only the larger->0.45-μm-size fraction can lead to underestimates of Caudovirales diversity and relative abundance. Similarly, our data show that examining only the smaller-size fraction can lead to underestimations of virophage and cyanophage relative abundances that could, in turn, cause researchers to assume their limited ecological importance. Given the considerable differences we observed in this study, we recommend cautious interpretations of environmental virus community assemblages and dynamics when based on metagenomic data derived from different size fractions.

Comparison of ultrafiltration and iron chloride flocculation in the preparation of aquatic viromes from contrasting sample types

Viral metagenomes (viromes) are a valuable untargeted tool for studying viral diversity and the central roles viruses play in host disease, ecology, and evolution. Establishing effective methods to concentrate and purify viral genomes prior to sequencing is essential for high quality viromes. Using virus spike-and-recovery experiments, we stepwise compared two common approaches for virus concentration, ultrafiltration and iron chloride flocculation, across diverse matrices: wastewater influent, wastewater secondary effluent, river water, and seawater. Viral DNA was purified by removing cellular DNA via chloroform cell lysis, filtration, and enzymatic degradation of extra-viral DNA. We found that viral genomes were concentrated 1-2 orders of magnitude more with ultrafiltration than iron chloride flocculation for all matrices and resulted in higher quality DNA suitable for amplification-free and long-read sequencing. Given its widespread use and utility as an inexpensive field method for virome sampling, we nonetheless sought to optimize iron flocculation. We found viruses were best concentrated in seawater with five-fold higher iron concentrations than the standard used, inhibition of DNase activity reduced purification effectiveness, and five-fold more iron was needed to flocculate viruses from freshwater than seawater—critical knowledge for those seeking to apply this broadly used method to freshwater virome samples. Overall, our results demonstrated that ultrafiltration and purification performed better than iron chloride flocculation and purification in the tested matrices. Given that the method performance depended on the solids content and salinity of the samples, we suggest spike-and-recovery experiments be applied when concentrating and purifying sample types that diverge from those tested here.

Interaction dynamics and virus-host range for estuarine actinophages captured by epicPCR

Viruses impact microbial diversity, gene flow and function through virus-host interactions. Although metagenomics surveys are rapidly cataloguing viral diversity, methods are needed to capture specific virus-host interactions in situ. Here, we leveraged metagenomics and repurposed emulsion paired isolation-concatenation PCR (epicPCR) to investigate viral diversity and virus-host interactions in situ over time in an estuarine environment. The method fuses a phage marker, the ribonucleotide reductase gene, with the host 16S rRNA gene of infected bacterial cells within emulsion droplets providing single-cell resolution for dozens of samples. EpicPCR captured in situ virus-host interactions for viral clades with no closely related database representatives. Abundant freshwater Actinobacteria lineages, in particular Rhodoluna sp., were the most common hosts for these poorly characterized viruses, with interactions correlated with environmental factors. Multiple methods used to identify virus-host interactions, including epicPCR, identified different and largely non-overlapping interactions within the vast virus-host interaction space. Tracking virus-host interaction dynamics also revealed that multi-host viruses had significantly longer periods with observed virus-host interactions, whereas single-host viruses were observed interacting with hosts at lower minimum abundances, suggesting more efficient interactions. Capturing in situ interactions with epicPCR revealed environmental and ecological factors shaping virus-host interactions, highlighting epicPCR as a valuable technique in viral ecology.

Abundance and Diversity of Phages, Microbial Taxa, and Antibiotic Resistance Genes in the Sediments of the River Ganges Through Metagenomic Approach

In this study, we have analyzed the metagenomic DNA from the pooled sediment sample of the river Ganges to explore the abundance and diversity of phages, microbial community, and antibiotic resistance genes (ARGs). Utilizing data from Illumina platform, 4,174 (∼0.0013%) reads were classified for the 285 different DNA viruses largely dominated by the group of 260 distinctive phages (3,602 reads, ∼86.3%). Among all, Microcystis (782 hits), Haemophilus (403), Synechococcus (386), Pseudomonas (279), Enterococcus (232), Bacillus (196), Rhodococcus (166), Caulobacter (163), Salmonella (146), Enterobacteria (143), Mycobacterium and (128) phages show the highest abundance and account for ∼90% of the total identified phages. In addition, we have also identified corresponding host pertaining to these phages. Mainly, Proteobacteria (∼69.3%) dominates the microbial population structure. Primarily, orders such as Caulobacterales (∼28%), Burkholderiales (∼13.9%), Actinomycetales (∼13.7%), and Pseudomonadales (∼7.5%) signify the core section. Furthermore, 21,869 (∼0.00695%) reads were classified in 20 ARG types (classes) and 240 ARGs (subtypes), among which 4 ARG types, namely multidrug resistance (12,041 reads, ∼55%), bacitracin (3,202 reads, ∼15%), macrolide-lincosamide-streptogramin (1,744 reads, ∼7.98%), and fosmidomycin (990 reads, ∼4.53%), have the highest abundance. Simultaneously, six resistance mechanisms were also recognized with the dominance of antibiotic efflux (72.8%, 15,919 reads). The results unveil the distribution of (pro)-phages; microbial community; and various ARGs in the Ganges river sediments.

Nanopore Sequencing Is a Credible Alternative to Recover Complete Genomes of Geminiviruses

Next-generation sequencing (NGS), through the implementation of metagenomic protocols, has led to the discovery of thousands of new viruses in the last decade. Nevertheless, these protocols are still laborious and costly to implement, and the technique has not yet become routine for everyday virus characterization. Within the context of CRESS DNA virus studies, we implemented two alternative long-read NGS protocols, one that is agnostic to the sequence (without a priori knowledge of the viral genome) and the other that use specific primers to target a virus (with a priori). Agnostic and specific long read NGS-based assembled genomes of two capulavirus strains were compared to those obtained using the gold standard technique of Sanger sequencing. Both protocols allowed the detection and accurate full genome characterization of both strains. Globally, the assembled genomes were very similar (99.5-99.7% identity) to the Sanger sequences consensus, but differences in the homopolymeric tracks of these sequences indicated a specific lack of accuracy of the long reads NGS approach that has yet to be improved. Nevertheless, the use of the bench-top sequencer has proven to be a credible alternative in the context of CRESS DNA virus study and could offer a new range of applications not previously accessible.

Extended Evaluation of Viral Diversity in Lake Baikal through Metagenomics

Lake Baikal is a unique oligotrophic freshwater lake with unusually cold conditions and amazing biological diversity. Studies of the lake's viral communities have begun recently, and their full diversity is not elucidated yet. Here, we performed DNA viral metagenomic analysis on integral samples from four different deep-water and shallow stations of the southern and central basins of the lake. There was a strict distinction of viral communities in areas with different environmental conditions. Comparative analysis with other freshwater lakes revealed the highest similarity of Baikal viromes with those of the Asian lakes Soyang and Biwa. Analysis of new data, together with previously published data allowed us to get a deeper insight into the diversity and functional potential of Baikal viruses; however, the true diversity of Baikal viruses in the lake ecosystem remains still unknown. The new metaviromic data will be useful for future studies of viral composition, distribution, and the dynamics associated with global climatic and anthropogenic impacts on this ecosystem.

Potential virus-mediated nitrogen cycling in oxygen-depleted oceanic waters

Viruses play an important role in the ecology and biogeochemistry of marine ecosystems. Beyond mortality and gene transfer, viruses can reprogram microbial metabolism during infection by expressing auxiliary metabolic genes (AMGs) involved in photosynthesis, central carbon metabolism, and nutrient cycling. While previous studies have focused on AMG diversity in the sunlit and dark ocean, less is known about the role of viruses in shaping metabolic networks along redox gradients associated with marine oxygen minimum zones (OMZs). Here, we analyzed relatively quantitative viral metagenomic datasets that profiled the oxygen gradient across Eastern Tropical South Pacific (ETSP) OMZ waters, assessing whether OMZ viruses might impact nitrogen (N) cycling via AMGs. Identified viral genomes encoded six N-cycle AMGs associated with denitrification, nitrification, assimilatory nitrate reduction, and nitrite transport. The majority of these AMGs (80%) were identified in T4-like Myoviridae phages, predicted to infect Cyanobacteria and Proteobacteria, or in unclassified archaeal viruses predicted to infect Thaumarchaeota. Four AMGs were exclusive to anoxic waters and had distributions that paralleled homologous microbial genes. Together, these findings suggest viruses modulate N-cycling processes within the ETSP OMZ and may contribute to nitrogen loss throughout the global oceans thus providing a baseline for their inclusion in the ecosystem and geochemical models.

The potential of viruses to influence the magnitude of greenhouse gas emissions in an inland wetland

Wetlands are among the earth's most efficient ecosystems for carbon sequestration, but can also emit potent greenhouse gases (GHGs) depending on how they are managed. Global management strategies have sought to maximize carbon drawdown by wetlands by manipulating wetland hydrology to inhibit bacterially-mediated emissions. However, it has recently been hypothesized within wetlands that viruses have the potential to dictate the magnitude and direction of GHG emissions by attacking prokaryotes involved in the carbon cycle. Here we tested this hypothesis in a whole-ecosystem manipulation by hydrologically-restoring a degraded wetland ('rewetting') and investigated the changes in GHG emissions, prokaryotes, viruses, and virus-host interactions. We found that hydrological restoration significantly increased prokaryotic diversity, methanogenic Methanomicrobia, as well as putative iron/sulfate-cyclers (Geobacteraceae), nitrogen-cyclers (Nitrosomonadaceae), and fermentative bacteria (Koribacteraceae). These results provide insights into successional microbial community shifts during rehabilitation. Additionally, in response to watering, viral-induced prokaryotic mortality declined by 77%, resulting in limited carbon released by viral shunt that was significantly correlated with the 2.8-fold reduction in wetland carbon emissions. Our findings highlight, for the first time, the potential implications of viral infections in soil prokaryotes on wetland greenhouse gas dynamics and confirm the importance of wetland rehabilitation as a tool to offset carbon emissions.

VirION2: a short- and long-read sequencing and informatics workflow to study the genomic diversity of viruses in nature

Microbes play fundamental roles in shaping natural ecosystem properties and functions, but do so under constraints imposed by their viral predators. However, studying viruses in nature can be challenging due to low biomass and the lack of universal gene markers. Though metagenomic short-read sequencing has greatly improved our virus ecology toolkit—and revealed many critical ecosystem roles for viruses—microdiverse populations and fine-scale genomic traits are missed. Some of these microdiverse populations are abundant and the missed regions may be of interest for identifying selection pressures that underpin evolutionary constraints associated with hosts and environments. Though long-read sequencing promises complete virus genomes on single reads, it currently suffers from high DNA requirements and sequencing errors that limit accurate gene prediction. Here we introduce VirION2, an integrated short- and long-read metagenomic wet-lab and informatics pipeline that updates our previous method (VirION) to further enhance the utility of long-read viral metagenomics. Using a viral mock community, we first optimized laboratory protocols (polymerase choice, DNA shearing size, PCR cycling) to enable 76% longer reads (now median length of 6,965 bp) from 100-fold less input DNA (now 1 nanogram). Using a virome from a natural seawater sample, we compared viromes generated with VirION2 against other library preparation options (unamplified, original VirION, and short-read), and optimized downstream informatics for improved long-read error correction and assembly. VirION2 assemblies combined with short-read based data (‘enhanced’ viromes), provided significant improvements over VirION libraries in the recovery of longer and more complete viral genomes, and our optimized error-correction strategy using long- and short-read data achieved 99.97% accuracy. In the seawater virome, VirION2 assemblies captured 5,161 viral populations (including all of the virus populations observed in the other assemblies), 30% of which were uniquely assembled through inclusion of long-reads, and 22% of the top 10% most abundant virus populations derived from assembly of long-reads. Viral populations unique to VirION2 assemblies had significantly higher microdiversity means, which may explain why short-read virome approaches failed to capture them. These findings suggest the VirION2 sample prep and workflow can help researchers better investigate the virosphere, even from challenging low-biomass samples. Our new protocols are available to the research community on protocols.io as a ‘living document’ to facilitate dissemination of updates to keep pace with the rapid evolution of long-read sequencing technology.

A persistent giant algal virus, with a unique morphology, encodes an unprecedented number of genes involved in energy metabolism

Viruses have long been viewed as entities possessing extremely limited metabolic capacities. Over the last decade, however, this view has been challenged, as metabolic genes have been identified in viruses possessing large genomes and virions-the synthesis of which is energetically demanding. Here, we unveil peculiar phenotypic and genomic features of Prymnesium kappa virus RF01 (PkV RF01), a giant virus of the Mimiviridae family. We found that this virus encodes an unprecedented number of proteins involved in energy metabolism, such as all four succinate dehydrogenase (SDH) subunits (A-D) as well as key enzymes in the β-oxidation pathway. The SDHA gene was transcribed upon infection, indicating that the viral SDH is actively used by the virus- potentially to modulate its host's energy metabolism. We detected orthologous SDHA and SDHB genes in numerous genome fragments from uncultivated marine Mimiviridae viruses, which suggests that the viral SDH is widespread in oceans. PkV RF01 was less virulent compared with other cultured prymnesioviruses, a phenomenon possibly linked to the metabolic capacity of this virus and suggestive of relatively long co-evolution with its hosts. It also has a unique morphology, compared to other characterized viruses in the Mimiviridae family. Finally, we found that PkV RF01 is the only alga-infecting Mimiviridae virus encoding two aminoacyl-tRNA synthetases and enzymes corresponding to an entire base-excision repair pathway, as seen in heterotroph-infecting Mimiviridae These Mimiviridae encoded-enzymes were found to be monophyletic and branching at the root of the eukaryotic tree of life. This placement suggests that the last common ancestor of Mimiviridae was endowed with a large, complex genome prior to the divergence of known extant eukaryotes.IMPORTANCE Viruses on Earth are tremendously diverse in terms of morphology, functionality, and genomic composition. Over the last decade, the conceptual gap separating viruses and cellular life has tightened because of the detection of metabolic genes in viral genomes that express complex virus phenotypes upon infection. Here, we describe Prymnesium kappa virus RF01, a large alga-infecting virus with a unique morphology, an atypical infection profile, and an unprecedented number of genes involved in energy metabolism (such as the tricarboxylic (TCA) cycle and the β-oxidation pathway). Moreover, we show that the gene corresponding to one of these enzymes (the succinate dehydrogenase subunit A) is transcribed during infection and is widespread among marine viruses. This discovery provides evidence that a virus has the potential to actively regulate energy metabolism with its own gene.

Genomic Characterization and Distribution Pattern of a Novel Marine OM43 Phage

Bacteriophages have a significant impact on the structure and function of marine microbial communities. Phages of some major bacterial lineages have recently been shown to dominate the marine viral communities. However, phages that infect many important bacterial clades still remained unexplored. Members of the marine OM43 clade are methylotrophs that play important roles in C1 metabolism. OM43 phages (phages that infect the OM43 bacteria) represent an understudied viral group with only one known isolate. In this study, we describe the genomic characterization and biogeography of an OM43 phage that infects the strain HTCC2181, designated MEP301. MEP301 has a genome size of 34,774 bp. We found that MEP301 is genetically distinct from other known phage isolates and only displays significant sequence similarity with some metagenomic viral genomes (MVGs). A total of 12 MEP301-type MVGs were identified from metagenomic datasets. Comparative genomic and phylogenetic analyses revealed that MEP301-type phages can be separated into two subgroups (subgroup I and subgroup II). We also performed a metagenomic recruitment analysis to determine the relative abundance of reads mapped to these MEP301-type phages, which suggested that subgroup I MEP301-type phages are present predominantly in the cold upper waters with lower salinity. Notably, subgroup II phages have an inverse different distribution pattern, implying that they may infect hosts from a distinct OM43 subcluster. Our study has expanded the knowledge about the genomic diversity of marine OM43 phages and identified a new phage group that is widespread in the ocean.

Metagenomic Data Assembly - The Way of Decoding Unknown Microorganisms

Metagenomics is a segment of conventional microbial genomics dedicated to the sequencing and analysis of combined genomic DNA of entire environmental samples. The most critical step of the metagenomic data analysis is the reconstruction of individual genes and genomes of the microorganisms in the communities using metagenomic assemblers - computational programs that put together small fragments of sequenced DNA generated by sequencing instruments. Here, we describe the challenges of metagenomic assembly, a wide spectrum of applications in which metagenomic assemblies were used to better understand the ecology and evolution of microbial ecosystems, and present one of the most efficient microbial assemblers, SPAdes that was upgraded to become applicable for metagenomics.

Hybrid assembly of an agricultural slurry virome reveals a diverse and stable community with the potential to alter the metabolism and virulence of veterinary pathogens

BACKGROUND

Viruses are the most abundant biological entities on Earth, known to be crucial components of microbial ecosystems. However, there is little information on the viral community within agricultural waste. There are currently ~ 2.7 million dairy cattle in the UK producing 7-8% of their own bodyweight in manure daily, and 28 million tonnes annually. To avoid pollution of UK freshwaters, manure must be stored and spread in accordance with guidelines set by DEFRA. Manures are used as fertiliser, and widely spread over crop fields, yet little is known about their microbial composition. We analysed the virome of agricultural slurry over a 5-month period using short and long-read sequencing.

RESULTS

Hybrid sequencing uncovered more high-quality viral genomes than long or short-reads alone; yielding 7682 vOTUs, 174 of which were complete viral genomes. The slurry virome was highly diverse and dominated by lytic bacteriophage, the majority of which represent novel genera (~ 98%). Despite constant influx and efflux of slurry, the composition and diversity of the slurry virome was extremely stable over time, with 55% of vOTUs detected in all samples over a 5-month period. Functional annotation revealed a diverse and abundant range of auxiliary metabolic genes and novel features present in the community, including the agriculturally relevant virulence factor VapE, which was widely distributed across different phage genera that were predicted to infect several hosts. Furthermore, we identified an abundance of phage-encoded diversity-generating retroelements, which were previously thought to be rare on lytic viral genomes. Additionally, we identified a group of crAssphages, including lineages that were previously thought only to be found in the human gut.

CONCLUSIONS

The cattle slurry virome is complex, diverse and dominated by novel genera, many of which are not recovered using long or short-reads alone. Phages were found to encode a wide range of AMGs that are not constrained to particular groups or predicted hosts, including virulence determinants and putative ARGs. The application of agricultural slurry to land may therefore be a driver of bacterial virulence and antimicrobial resistance in the environment. Video abstract.

Uncultivated Viral Populations Dominate Estuarine Viromes on the Spatiotemporal Scale

Viruses are ubiquitous and abundant in the oceans, and viral metagenomes (viromes) have been investigated extensively via several large-scale ocean sequencing projects. However, there have not been any systematic viromic studies in estuaries. Here, we investigated the viromes of the Delaware Bay and Chesapeake Bay, two Mid-Atlantic estuaries. Deep sequencing generated a total of 48,190 assembled viral sequences (>5 kb) and 26,487 viral populations (9,204 virus clusters and 17,845 singletons), including 319 circular viral contigs between 7.5 kb and 161.8 kb. Unknown viruses represented the vast majority of the dominant populations, while the composition of known viruses, such as pelagiphage and cyanophage, appeared to be relatively consistent across a wide range of salinity gradients and in different seasons. A difference between estuarine and ocean viromes was reflected by the proportions of Myoviridae, Podoviridae, Siphoviridae, Phycodnaviridae, and a few well-studied virus representatives. The difference in viral community between the Delaware Bay and Chesapeake Bay is significantly more pronounced than the difference caused by temperature or salinity, indicating strong local profiles caused by the unique ecology of each estuary. Interestingly, a viral contig similar to phages infecting Acinetobacter baumannii (“Iraqibacter”) was found to be highly abundant in the Delaware Bay but not in the Chesapeake Bay, the source of which is yet to be identified. Highly abundant viruses in both estuaries have close hits to viral sequences derived from the marine single-cell genomes or long-read single-molecule sequencing, suggesting that important viruses are still waiting to be discovered in the estuarine environment.

IMPORTANCE This is the first systematic study about spatial and temporal variation of virioplankton communities in estuaries using deep metagenomics sequencing. It is among the highest-quality viromic data sets to date, showing remarkably consistent sequencing depth and quality across samples. Our results indicate that there exists a large pool of abundant and diverse viruses in estuaries that have not yet been cultivated, their genomes only available thanks to single-cell genomics or single-molecule sequencing, demonstrating the importance of these methods for viral discovery. The spatiotemporal pattern of these abundant uncultivated viruses is more variable than that of cultured viruses. Despite strong environmental gradients, season and location had surprisingly little impact on the viral community within an estuary, but we saw a significant distinction between the two estuaries and also between estuarine and open ocean viromes.

The communities and functional profiles of virioplankton along a salinity gradient in a subtropical estuary

Viruses are the major drivers shaping microorganismal communities, and impact marine biogeochemical cycling. They are affected by various environmental parameters, such as salinity. Although the spatiotemporal distribution and dynamics of virioplankton have been extensively studied in saline environments, few detailed studies of community structure and function of viruses along salinity gradients have been conducted. Here, we used the 16S and 18S rRNA gene amplicon and metagenomic sequencing from a subtropical estuary (Pearl River Estuary, PRE; located in Shenzhen, Guangdong Province, China) to explore how viral community composition and function vary along a salinity gradient. Results showed that the detected viruses were mainly bacteriophages. The double-stranded DNA viruses were the most abundant (especially Siphoviridae, Myoviridae, Mimiviridae, Phycodnaviridae, and Podoviridae), followed by a small number of single-stranded DNA (Circoviridae) and RNA (Retroviridae) viruses. Viral biodiversity significantly declined and community structure varied greatly along the salinity gradient. The salinity, ammonium and dissolved oxygen were dominated factors influencing the community composition of viruses. Association network analysis showed that viruses had a negative effect on multiple host taxa (prokaryotic and eukaryotic species). Metagenomic data revealed that the main viral functional potential was involved in organic matter metabolism by carbohydrate-active enzymes (CAZymes). Deeper comparative functional analyses showed that viruses in the low-salinity environment had more carbohydrate-binding module and glycosidase hydrolases activities than those under high-salinity conditions. However, an opposite pattern was observed for carbohydrate esterases. These results suggest that virus-encoded CAZyme genes may alter the bacterial metabolism in estuaries. Overall, our results demonstrate that there is a spatial heterogeneity in the composition and function of virioplankton along a salinity gradient. This study enhances our understanding of viral distribution and their contribution to regulating carbon degradation throughout environments with varying salinities in subtropical estuaries.

The Fennoscandian Shield deep terrestrial virosphere suggests slow motion 'boom and burst' cycles

The deep biosphere contains members from all three domains of life along with viruses. Here we investigate the deep terrestrial virosphere by sequencing community nucleic acids from three groundwaters of contrasting chemistries, origins, and ages. These viromes constitute a highly unique community compared to other environmental viromes and sequenced viral isolates. Viral host prediction suggests that many of the viruses are associated with Firmicutes and Patescibacteria, a superphylum lacking previously described active viruses. RNA transcript-based activity implies viral predation in the shallower marine water-fed groundwater, while the deeper and more oligotrophic waters appear to be in ‘metabolic standby’. Viral encoded antibiotic production and resistance systems suggest competition and antagonistic interactions. The data demonstrate a viral community with a wide range of predicted hosts that mediates nutrient recycling to support a higher microbial turnover than previously anticipated. This suggests the presence of ‘kill-the-winner’ oscillations creating slow motion ‘boom and burst’ cycles.

Karin Holmfeldt et al. sequence metagenomes and metatranscriptomes of viruses in deep groundwaters down to 448 m below the surface. The results reveal ecological dynamics of viruses including slow motion ‘boom and burst’ cycles and a ‘kill the winner’ model potentially driven by viral predation.

MetaHiC phage-bacteria infection network reveals active cycling phages of the healthy human gut

Bacteriophages play important roles in regulating the intestinal human microbiota composition, dynamics, and homeostasis, and characterizing their bacterial hosts is needed to understand their impact. We applied a metagenomic Hi-C approach on 10 healthy human gut samples to unveil a large infection network encompassing more than 6000 interactions bridging a metagenomic assembled genomes (MAGs) and a phage sequence, allowing to study in situ phage-host ratio. Whereas three-quarters of these sequences likely correspond to dormant prophages, 5% exhibit a much higher coverage than their associated MAG, representing potentially actively replicating phages. We detected 17 sequences of members of the crAss-like phage family, whose hosts diversity remained until recently relatively elusive. For each of them, a unique bacterial host was identified, all belonging to different genus of Bacteroidetes. Therefore, metaHiC deciphers infection network of microbial population with a high specificity paving the way to dynamic analysis of mobile genetic elements in complex ecosystems.

Metaviromics coupled with phage-host identification to open the viral 'black box'

Viruses are found in almost all biomes on Earth, with bacteriophages (phages) accounting for the majority of viral particles in most ecosystems. Phages have been isolated from natural environments using the plaque assay and liquid medium-based dilution culturing. However, phage cultivation is restricted by the current limitations in the number of culturable bacterial strains. Unlike prokaryotes, which possess universally conserved 16S rRNA genes, phages lack universal marker genes for viral taxonomy, thus restricting cultureindependent analyses of viral diversity. To circumvent these limitations, shotgun viral metagenome sequencing (i.e., metaviromics) has been developed to enable the extensive sequencing of a variety of viral particles present in the environment and is now widely used. Using metaviromics, numerous studies on viral communities have been conducted in oceans, lakes, rivers, and soils, resulting in many novel phage sequences. Furthermore, auxiliary metabolic genes such as ammonic monooxygenase C and β-lactamase have been discovered in viral contigs assembled from viral metagenomes. Current attempts to identify putative bacterial hosts of viral metagenome sequences based on sequence homology have been limited due to viral sequence variations. Therefore, culture-independent approaches have been developed to predict bacterial hosts using single-cell genomics and fluorescentlabeling. This review focuses on recent viral metagenome studies conducted in natural environments, especially in aquatic ecosystems, and their contributions to phage ecology. Here, we concluded that although metaviromics is a key tool for the study of viral ecology, this approach must be supplemented with phage-host identification, which in turn requires the cultivation of phage-bacteria systems.

Housekeeping in the Hydrosphere: Microbial Cooking, Cleaning, and Control under Stress

Who’s cooking, who’s cleaning, and who’s got the remote control within the waters blanketing Earth? Anatomically tiny, numerically dominant microbes are the crucial “homemakers” of the watery household. Phytoplankton’s culinary abilities enable them to create food by absorbing sunlight to fix carbon and release oxygen, making microbial autotrophs top-chefs in the aquatic kitchen. However, they are not the only bioengineers that balance this complex household. Ubiquitous heterotrophic microbes including prokaryotic bacteria and archaea (both “bacteria” henceforth), eukaryotic protists, and viruses, recycle organic matter and make inorganic nutrients available to primary producers. Grazing protists compete with viruses for bacterial biomass, whereas mixotrophic protists produce new organic matter as well as consume microbial biomass. When viruses press remote-control buttons, by modifying host genomes or lysing them, the outcome can reverberate throughout the microbial community and beyond. Despite recognition of the vital role of microbes in biosphere housekeeping, impacts of anthropogenic stressors and climate change on their biodiversity, evolution, and ecological function remain poorly understood. How trillions of the smallest organisms in Earth’s largest ecosystem respond will be hugely consequential. By making the study of ecology personal, the “housekeeping” perspective can provide better insights into changing ecosystem structure and function at all scales.

Mitochondrial Genome Evolution in Pelagophyte Algae

The Pelagophyceae are marine stramenopile algae that include Aureoumbra lagunensis and Aureococcus anophagefferens, two microbial species notorious for causing harmful algal blooms. Despite their ecological significance, relatively few genomic studies of pelagophytes have been carried out. To improve understanding of the biology and evolution of pelagophyte algae, we sequenced complete mitochondrial genomes for A. lagunensis (CCMP1510), Pelagomonas calceolata (CCMP1756), and five strains of Aureoc. anophagefferens (CCMP1707, CCMP1708, CCMP1850, CCMP1984, and CCMP3368) using Nanopore long-read sequencing. All pelagophyte mitochondrial genomes assembled into single, circular mapping contigs between 39,376 bp (P. calceolata) and 55,968 bp (A. lagunensis) in size. Mitochondrial genomes for the five Aureoc. anophagefferens strains varied slightly in length (42,401–42,621 bp) and were 99.4–100.0% identical. Gene content and order were highly conserved between the Aureoc. anophagefferens and P. calceolata genomes, with the only major difference being a unique region in Aureoc. anophagefferens containingDNA adenine and cytosine methyltransferase (dam/dcm) genes that appear to be the product of lateral gene transfer from a prokaryotic or viral donor. Although the A. lagunensis mitochondrial genome shares seven distinct syntenic blocks with the other pelagophyte genomes, it has a tandem repeat expansion comprising ∼40% of its length, and lacks identifiable rps19 and glycine tRNA genes. Laterally acquired self-splicing introns were also found in the 23S rRNA (rnl) gene of P. calceolata and the coxI gene of the five Aureoc. anophagefferens genomes. Overall, these data provide baseline knowledge about the genetic diversity of bloom-forming pelagophytes relative to nonbloom-forming species.

VirSorter2: a multi-classifier, expert-guided approach to detect diverse DNA and RNA viruses

BACKGROUND

Viruses are a significant player in many biosphere and human ecosystems, but most signals remain "hidden" in metagenomic/metatranscriptomic sequence datasets due to the lack of universal gene markers, database representatives, and insufficiently advanced identification tools.

RESULTS

Here, we introduce VirSorter2, a DNA and RNA virus identification tool that leverages genome-informed database advances across a collection of customized automatic classifiers to improve the accuracy and range of virus sequence detection. When benchmarked against genomes from both isolated and uncultivated viruses, VirSorter2 uniquely performed consistently with high accuracy (F1-score > 0.8) across viral diversity, while all other tools under-detected viruses outside of the group most represented in reference databases (i.e., those in the order Caudovirales). Among the tools evaluated, VirSorter2 was also uniquely able to minimize errors associated with atypical cellular sequences including eukaryotic genomes and plasmids. Finally, as the virosphere exploration unravels novel viral sequences, VirSorter2's modular design makes it inherently able to expand to new types of viruses via the design of new classifiers to maintain maximal sensitivity and specificity.

CONCLUSION

With multi-classifier and modular design, VirSorter2 demonstrates higher overall accuracy across major viral groups and will advance our knowledge of virus evolution, diversity, and virus-microbe interaction in various ecosystems. Source code of VirSorter2 is freely available ( https://bitbucket.org/MAVERICLab/virsorter2 ), and VirSorter2 is also available both on bioconda and as an iVirus app on CyVerse ( https://de.cyverse.org/de ). Video abstract.

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.

mMGE: a database for human metagenomic extrachromosomal mobile genetic elements

Extrachromosomal mobile genetic elements (eMGEs), including phages and plasmids, that can move across different microbes, play important roles in genome evolution and shaping the structure of microbial communities. However, we still know very little about eMGEs, especially their abundances, distributions and putative functions in microbiomes. Thus, a comprehensive description of eMGEs is of great utility. Here we present mMGE, a comprehensive catalog of 517 251 non-redundant eMGEs, including 92 492 plasmids and 424 759 phages, derived from diverse body sites of 66 425 human metagenomic samples. About half the eMGEs could be further grouped into 70 074 clusters using relaxed criteria (referred as to eMGE clusters below). We provide extensive annotations of the identified eMGEs including sequence characteristics, taxonomy affiliation, gene contents and their prokaryotic hosts. We also calculate the prevalence, both within and across samples for each eMGE and eMGE cluster, enabling users to see putative associations of eMGEs with human phenotypes or their distribution preferences. All eMGE records can be browsed or queried in multiple ways, such as eMGE clusters, metagenomic samples and associated hosts. The mMGE is equipped with a user-friendly interface and a BLAST server, facilitating easy access/queries to all its contents easily. mMGE is freely available for academic use at: https://mgedb.comp-sysbio.org.

The International Virus Bioinformatics Meeting 2020

The International Virus Bioinformatics Meeting 2020 was originally planned to take place in Bern, Switzerland, in March 2020. However, the COVID-19 pandemic put a spoke in the wheel of almost all conferences to be held in 2020. After moving the conference to 8-9 October 2020, we got hit by the second wave and finally decided at short notice to go fully online. On the other hand, the pandemic has made us even more aware of the importance of accelerating research in viral bioinformatics. Advances in bioinformatics have led to improved approaches to investigate viral infections and outbreaks. The International Virus Bioinformatics Meeting 2020 has attracted approximately 120 experts in virology and bioinformatics from all over the world to join the two-day virtual meeting. Despite concerns being raised that virtual meetings lack possibilities for face-to-face discussion, the participants from this small community created a highly interactive scientific environment, engaging in lively and inspiring discussions and suggesting new research directions and questions. The meeting featured five invited and twelve contributed talks, on the four main topics: (1) proteome and RNAome of RNA viruses, (2) viral metagenomics and ecology, (3) virus evolution and classification and (4) viral infections and immunology. Further, the meeting featured 20 oral poster presentations, all of which focused on specific areas of virus bioinformatics. This report summarizes the main research findings and highlights presented at the meeting.

Genome-resolved viral ecology in a marine oxygen minimum zone

Oxygen minimum zones (OMZs) are critical to marine nitrogen cycling and global climate change. While OMZ microbial communities are relatively well-studied, little is known about their viruses. Here, we assess the viral community ecology of 22 deeply sequenced viral metagenomes along a gradient of oxygenated to anoxic waters (<0.02 μmol/l O₂ ) in the Eastern Tropical South Pacific (ETSP) OMZ. We identified 46 127 viral populations (≥5 kb), which augments the known viruses from ETSP by 10-fold. Viral communities clustered into six groups that correspond to oceanographic features. Oxygen concentration was the predominant environmental feature driving viral community structure. Alpha and beta diversity of viral communities in the anoxic zone were lower than in surface waters, which parallels the low microbial diversity seen in other studies. ETSP viruses were largely endemic, with the majority of shared viruses (87%) also present in other OMZ samples. We detected 543 putative viral-encoded auxiliary metabolic genes (AMGs), of which some have a distribution that reflects physico-chemical characteristics across depth. Together these findings provide an ecological baseline for viral community structure, drivers and population variability in OMZs that will help future studies assess the role of viruses in these climate-critical environments.

Widespread endogenization of giant viruses shapes genomes of green algae

Endogenous viral elements (EVEs)-viruses that have integrated their genomes into those of their hosts-are prevalent in eukaryotes and have an important role in genome evolution^1,2. The vast majority of EVEs that have been identified to date are small genomic regions comprising a few genes², but recent evidence suggests that some large double-stranded DNA viruses may also endogenize into the genome of the host¹. Nucleocytoplasmic large DNA viruses (NCLDVs) have recently become of great interest owing to their large genomes and complex evolutionary origins^3-6, but it is not yet known whether they are a prominent component of eukaryotic EVEs. Here we report the widespread endogenization of NCLDVs in diverse green algae; these giant EVEs reached sizes greater than 1 million base pairs and contained as many as around 10% of the total open reading frames in some genomes, substantially increasing the scale of known viral genes in eukaryotic genomes. These endogenized elements often shared genes with host genomic loci and contained numerous spliceosomal introns and large duplications, suggesting tight assimilation into host genomes. NCLDVs contain large and mosaic genomes with genes derived from multiple sources, and their endogenization represents an underappreciated conduit of new genetic material into eukaryotic lineages that can substantially impact genome composition.

The Gut Virome Database Reveals Age-Dependent Patterns of Virome Diversity in the Human Gut

The gut microbiome profoundly affects human health and disease, and their infecting viruses are likely as important, but often missed because of reference database limitations. Here, we (1) built a human Gut Virome Database (GVD) from 2,697 viral particle or microbial metagenomes from 1,986 individuals representing 16 countries, (2) assess its effectiveness, and (3) report a meta-analysis that reveals age-dependent patterns across healthy Westerners. The GVD contains 33,242 unique viral populations (approximately species-level taxa) and improves average viral detection rates over viral RefSeq and IMG/VR nearly 182-fold and 2.6-fold, respectively. GVD meta-analyses show highly personalized viromes, reveal that inter-study variability from technical artifacts is larger than any "disease" effect at the population level, and document how viral diversity changes from human infancy into senescence. Together, this compact foundational resource, these standardization guidelines, and these meta-analysis findings provide a systematic toolkit to help maximize our understanding of viral roles in health and disease.

The Microbial Composition in Circumneutral Thermal Springs from Chignahuapan, Puebla, Mexico Reveals the Presence of Particular Sulfur-Oxidizing Bacterial and Viral Communities

Terrestrial thermal springs are widely distributed globally, and these springs harbor a broad diversity of organisms of biotechnological interest. In Mexico, few studies exploring this kind of environment have been described. In this work, we explore the microbial community in Chignahuapan hot springs, which provides clues to understand these ecosystems' diversity. We assessed the diversity of the microorganism communities in a hot spring environment with a metagenomic shotgun approach. Besides identifying similarities and differences with other ecosystems, we achieved a systematic comparison against 11 metagenomic samples from diverse localities. The Chignahuapan hot springs show a particular prevalence of sulfur-oxidizing bacteria from the genera Rhodococcus, Thermomonas, Thiomonas, Acinetobacter, Sulfurovum, and Bacillus, highlighting those that are different from other recovered bacterial populations in circumneutral hot springs environments around the world. The co-occurrence analysis of the bacteria and viruses in these environments revealed that within the Rhodococcus, Thiomonas, Thermonas, and Bacillus genera, the Chignahuapan samples have specific species of bacteria with a particular abundance, such as Rhodococcus erytropholis. The viruses in the circumneutral hot springs present bacteriophages within the order Caudovirales (Siphoviridae, Myoviridae, and Podoviridae), but the family of Herelleviridae was the most abundant in Chignahuapan samples. Furthermore, viral auxiliary metabolic genes were identified, many of which contribute mainly to the metabolism of cofactors and vitamins as well as carbohydrate metabolism. Nevertheless, the viruses and bacteria present in the circumneutral environments contribute to the sulfur cycle. This work represents an exhaustive characterization of a community structure in samples collected from hot springs in Mexico and opens opportunities to identify organisms of biotechnological interest.

DRAM for distilling microbial metabolism to automate the curation of microbiome function

Microbial and viral communities transform the chemistry of Earth's ecosystems, yet the specific reactions catalyzed by these biological engines are hard to decode due to the absence of a scalable, metabolically resolved, annotation software. Here, we present DRAM (Distilled and Refined Annotation of Metabolism), a framework to translate the deluge of microbiome-based genomic information into a catalog of microbial traits. To demonstrate the applicability of DRAM across metabolically diverse genomes, we evaluated DRAM performance on a defined, in silico soil community and previously published human gut metagenomes. We show that DRAM accurately assigned microbial contributions to geochemical cycles and automated the partitioning of gut microbial carbohydrate metabolism at substrate levels. DRAM-v, the viral mode of DRAM, established rules to identify virally-encoded auxiliary metabolic genes (AMGs), resulting in the metabolic categorization of thousands of putative AMGs from soils and guts. Together DRAM and DRAM-v provide critical metabolic profiling capabilities that decipher mechanisms underpinning microbiome function.

Metatranscriptomic Identification of Diverse and Divergent RNA Viruses in Green and Chlorarachniophyte Algae Cultures

Our knowledge of the diversity and evolution of the virosphere will likely increase dramatically with the study of microbial eukaryotes, including the microalgae within which few RNA viruses have been documented. By combining total RNA sequencing with sequence and structural-based homology detection, we identified 18 novel RNA viruses in cultured samples from two major groups of microbial algae: the chlorophytes and the chlorarachniophytes. Most of the RNA viruses identified in the green algae class Ulvophyceae were related to the Tombusviridae and Amalgaviridae viral families commonly associated with land plants. This suggests that the evolutionary history of these viruses extends to divergence events between algae and land plants. Seven Ostreobium sp-associated viruses exhibited sequence similarity to the mitoviruses most commonly found in fungi, compatible with horizontal virus transfer between algae and fungi. We also document, for the first time, RNA viruses associated with chlorarachniophytes, including the first negative-sense (bunya-like) RNA virus in microalgae, as well as a distant homolog of the plant virus Virgaviridae, potentially signifying viral inheritance from the secondary chloroplast endosymbiosis that marked the origin of the chlorarachniophytes. More broadly, these data suggest that the scarcity of RNA viruses in algae results from limited investigation rather than their absence.

Culturing novel and abundant pelagiphages in the ocean

Viruses play a key role in biogeochemical cycling and host mortality, metabolism, physiology and evolution in the ocean. Viruses that infect the globally abundant SAR11 bacteria (pelagiphages) were reported to be an important component of the marine viral communities. Our current knowledge of pelagiphages is based on a few studies and therefore is limited. In this study, 10 new pelagiphages were isolated and genomically characterized. These pelagiphages represent the first cultivated representatives of four viral lineages only found in metagenomic sequencing datasets previously. Many abundant environmental viral sequences, i.e., single-virus vSAG 37-F6 and several Global Ocean Viromes (GOV) viral populations, are now further confirmed with these pelagiphages. Viromic read mapping reveals that these new pelagiphages are globally distributed in the ocean and can be detected throughout the water column. Remarkably, isolation of these pelagiphages contributed up to 12% of all viromic reads annotated in the analysed viromes. Altogether, this study has greatly broadened our understanding of pelagiphages regarding their morphology, genetic diversity, infection strategies, and distribution pattern. The availability of these newly isolated pelagiphages and their genome sequences will allow us to further explore their infectivities and ecological strategies.

Metatranscriptomic Identification of Diverse and Divergent RNA Viruses in Green and Chlorarachniophyte Algae Cultures

Our knowledge of the diversity and evolution of the virosphere will likely increase dramatically with the study of microbial eukaryotes, including the microalgae within which few RNA viruses have been documented. By combining total RNA sequencing with sequence and structural-based homology detection, we identified 18 novel RNA viruses in cultured samples from two major groups of microbial algae: the chlorophytes and the chlorarachniophytes. Most of the RNA viruses identified in the green algae class Ulvophyceae were related to the Tombusviridae and Amalgaviridae viral families commonly associated with land plants. This suggests that the evolutionary history of these viruses extends to divergence events between algae and land plants. Seven Ostreobium sp-associated viruses exhibited sequence similarity to the mitoviruses most commonly found in fungi, compatible with horizontal virus transfer between algae and fungi. We also document, for the first time, RNA viruses associated with chlorarachniophytes, including the first negative-sense (bunya-like) RNA virus in microalgae, as well as a distant homolog of the plant virus Virgaviridae, potentially signifying viral inheritance from the secondary chloroplast endosymbiosis that marked the origin of the chlorarachniophytes. More broadly, these data suggest that the scarcity of RNA viruses in algae results from limited investigation rather than their absence.

Viral metagenomes of Lake Soyang, the largest freshwater lake in South Korea

A high number of viral metagenomes have revealed countless genomes of putative bacteriophages that have not yet been identified due to limitations in bacteriophage cultures. However, most virome studies have been focused on marine or gut environments, thereby leaving the viral community structure of freshwater lakes unclear. Because the lakes located around the globe have independent ecosystems with unique characteristics, viral community structures are also distinctive but comparable. Here, we present data on viral metagenomes that were seasonally collected at a depth of 1 m from Lake Soyang, the largest freshwater reservoir in South Korea. Through shotgun metagenome sequencing using the Illumina MiSeq platform, 3.08 to 5.54-Gbps of reads per virome were obtained. To predict the viral genome sequences within Lake Soyang, contigs were constructed and 648 to 1,004 putative viral contigs were obtained per sample. We expect that both viral metagenome reads and viral contigs would contribute in comparing and understanding of viral communities among different freshwater lakes depending on seasonal changes.