Assessing the diversity and specificity of two freshwater viral communities through metagenomics

Virome Analysis of Paddy Floodwater in Two Cropping Seasons

Paddy fields include two interconnected ecosystems-soil and floodwater. Microbes and viruses are an integral component of these ecosystems, yet the viral communities have not been extensively studied. We present an analysis of the viromes of paddy floodwater collected during the two cropping seasons in India, the kharif and rabi seasons respectively. The overall taxonomic and functional characteristics appeared to be similar in both seasons, suggesting stability of the viral community. Taxonomically, the families of tailed bacteriophages dominated. The predominance of functional roles related to lytic phages further confirmed this. We reconstructed two complete and several partial viral genomes from the assembled data. The genomes did not align with any known sequences, thus representing novel viruses of the floodwater ecosystem.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12088-024-01292-9.

Sediment microbial diversity, functional potentials, and antibiotic resistance pattern: a case study of Cochin Estuary core sediment

Marine sediments are an important part of the marine environment and the world's greatest organic carbon source. Sediment microorganisms are important regulators of major geochemical and eco-environmental processes in marine environments, especially nutrient dynamics and biogeochemical cycles. Despite their importance, core marine microorganisms are virtually unknown due to a lack of consensus on how to identify them. Most core microbiotas have been characterized thus far based on species abundance and occurrence. The combined effects of habitat and depth on benthic bacterial communities and ecological functions were studied using "Next-Generation sequencing (NGS) and Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) predictive functional profiling" at the surface (0.2 cm) and bottom depth (250 cm) in a sediment core sample from Cochin Estuary, Kerala, India. The results showed that bacterial diversity and richness were significantly higher in the surface sediment sample with the most abundant phyla being Proteobacteria, Acidobacteria, Chloroflexi, and Bacteroidetes. The major metabolic functions were metabolism, followed by environmental information processing and genetic information processing. Antibiotic resistance genes between the surface and bottom samples help to understand the resistance pattern among multidrug resistance is the most prominent one. Among viruses, Siphoviridae is the dominant family, followed by Myoviridae. In the case of Archea, Crenarchaeota is dominant, whereas among eukaryotes phyla Streptophyta and Chordata were dominant in the surface and the bottom samples respectively.

From shells to sequences: A proof-of-concept study for on-site analysis of hemolymphatic circulating cell-free DNA from sentinel mussels using Nanopore technology

Blue mussels are often abundant and widely distributed in polar marine coastal ecosystems. Because of their wide distribution, ecological importance, and relatively stationary lifestyle, bivalves have long been considered suitable indicators of ecosystem health and changes. Monitoring the population dynamics of blue mussels can provide information on the overall biodiversity, species interactions, and ecosystem functioning. In the present work, we combined the concept of liquid biopsy (LB), an emerging concept in medicine based on the sequencing of free circulating DNA, with the Oxford Nanopore Technologies (ONT) platform using a portable laboratory in a remote area. Our results demonstrate that this platform is ideally suited for sequencing hemolymphatic circulating cell-free DNA (ccfDNA) fragments found in blue mussels. The percentage of non-self ccfDNA accounted for >50 % of ccfDNA at certain sampling Sites, allowing the quick, on-site acquisition of a global view of the biodiversity of a coastal marine ecosystem. These ccfDNA fragments originated from viruses, bacteria, plants, arthropods, algae, and multiple Chordata. Aside from non-self ccfDNA, we found DNA fragments from all 14 blue mussel chromosomes, as well as those originating from the mitochondrial genomes. However, the distribution of nuclear and mitochondrial DNA was significantly different between Sites. Similarly, analyses between various sampling Sites showed that the biodiversity varied significantly within microhabitats. Our work shows that the ONT platform is well-suited for LB in sentinel blue mussels in remote and challenging conditions, enabling faster fieldwork for conservation strategies and resource management in diverse settings.

Benchmarking informatics approaches for virus discovery: caution is needed when combining in silico identification methods

Understanding the ecological impacts of viruses on natural and engineered ecosystems relies on the accurate identification of viral sequences from community sequencing data. To maximize viral recovery from metagenomes, researchers frequently combine viral identification tools. However, the effectiveness of this strategy is unknown. Here, we benchmarked combinations of six widely used informatics tools for viral identification and analysis (VirSorter, VirSorter2, VIBRANT, DeepVirFinder, CheckV, and Kaiju), called "rulesets." Rulesets were tested against mock metagenomes composed of taxonomically diverse sequence types and diverse aquatic metagenomes to assess the effects of the degree of viral enrichment and habitat on tool performance. We found that six rulesets achieved equivalent accuracy [Matthews Correlation Coefficient (MCC) = 0.77, Padj ≥ 0.05]. Each contained VirSorter2, and five used our "tuning removal" rule designed to remove non-viral contamination. While DeepVirFinder, VIBRANT, and VirSorter were each found once in these high-accuracy rulesets, they were not found in combination with each other: combining tools does not lead to optimal performance. Our validation suggests that the MCC plateau at 0.77 is partly caused by inaccurate labeling within reference sequence databases. In aquatic metagenomes, our highest MCC ruleset identified more viral sequences in virus-enriched (44%-46%) than in cellular metagenomes (7%-19%). While improved algorithms may lead to more accurate viral identification tools, this should be done in tandem with careful curation of sequence databases. We recommend using the VirSorter2 ruleset and our empirically derived tuning removal rule. Our analysis provides insight into methods for in silico viral identification and will enable more robust viral identification from metagenomic data sets.

IMPORTANCE

The identification of viruses from environmental metagenomes using informatics tools has offered critical insights in microbial ecology. However, it remains difficult for researchers to know which tools optimize viral recovery for their specific study. In an attempt to recover more viruses, studies are increasingly combining the outputs from multiple tools without validating this approach. After benchmarking combinations of six viral identification tools against mock metagenomes and environmental samples, we found that these tools should only be combined cautiously. Two to four tool combinations maximized viral recovery and minimized non-viral contamination compared with either the single-tool or the five- to six-tool ones. By providing a rigorous overview of the behavior of in silico viral identification strategies and a pipeline to replicate our process, our findings guide the use of existing viral identification tools and offer a blueprint for feature engineering of new tools that will lead to higher-confidence viral discovery in microbiome studies.

Discovering and exploring the hidden diversity of human gut viruses using highly enriched virome samples

Viruses are an abundant and crucial component of the human microbiome, but accurately discovering them via metagenomics is still challenging. Currently, the available viral reference genomes poorly represent the diversity in microbiome samples, and expanding such a set of viral references is difficult. As a result, many viruses are still undetectable through metagenomics even when considering the power of de novo metagenomic assembly and binning, as viruses lack universal markers. Here, we describe a novel approach to catalog new viral members of the human gut microbiome and show how the resulting resource improves metagenomic analyses. We retrieved >3,000 viral-like particles (VLP) enriched metagenomic samples (viromes), evaluated the efficiency of the enrichment in each sample to leverage the viromes of highest purity, and applied multiple analysis steps involving assembly and comparison with hundreds of thousands of metagenome-assembled genomes to discover new viral genomes. We reported over 162,000 viral sequences passing quality control from thousands of gut metagenomes and viromes. The great majority of the retrieved viral sequences (~94.4%) were of unknown origin, most had a CRISPR spacer matching host bacteria, and four of them could be detected in >50% of a set of 18,756 gut metagenomes we surveyed. We included the obtained collection of sequences in a new MetaPhlAn 4.1 release, which can quantify reads within a metagenome matching the known and newly uncovered viral diversity. Additionally, we released the viral database for further virome and metagenomic studies of the human microbiome.

Metagenomic analysis of antibiotic-resistance genes and viruses released from glaciers into downstream habitats

Glaciers serve as effective reservoirs of antibiotic resistance genes (ARGs) and viruses for millions of years. Climate change and anthropogenic activity have accelerated the melting of glaciers, but the patterns of release of ARGs and viruses from melting glaciers into downstream habitats remain unknown. We analyzed 171 metagenomic samples from glaciers and their downstream habitats and found that the abundance and diversity of ARGs were higher in glaciers (polar and plateau glaciers) than downstream habitats (Arctic Ocean, Qinghai Lake, and Yangtze River Basin), with the diversity of viruses having the opposite pattern. Proteobacteria and Actinobacteria were the main potential hosts of ARGs and viruses, and the richness of ARGs carried by the hosts was positively correlated with viral abundance, suggesting that the transmission of viruses in the hosts could disseminate ARGs. Source tracking indicated that >18 % of the ARGs and >25 % of the viruses detected in downstream habitats originated from glaciers, demonstrating that glaciers could be one of the potential sources of ARGs and viruses in downstream habitats. Increased solar radiation and emission of carbon dioxide mainly influenced the release of the ARGs and viruses from glaciers into downstream habitats. This study provides a systematic insight demonstrating the release of ARGs and viruses from the melting glaciers, potentially increasing ecological pressure.

Revealing viral diversity in the Napahai plateau wetland based on metagenomics

We focused on exploring the diversity of viruses in the Napahai plateau wetland, a unique ecosystem located in Yunnan, China. While viruses in marine environments have been extensively studied for their influence on microbial metabolism and biogeochemical cycles, little is known about their composition and function in plateau wetlands. Metagenomic analysis was employed to investigate the viral diversity and biogeochemical impacts in the Napahai wetland. It revealed that the Caudoviricetes and Malgrandaviricetes class level was the most abundant viral category based on phylogenetic analysis. Additionally, a gene-sharing network highlighted the presence of numerous unexplored viruses and demonstrated their unique characteristics and significant variation within the viral community of the Napahai wetland. Furthermore, the study identified the auxiliary metabolic genes (AMGs). AMGs provide phages with additional functions, such as protection against host degradation and involvement in metabolic pathways, such as the pentose phosphate pathway and DNA biosynthesis. The viruses in the Napahai wetland were found to influence carbon, nitrogen, sulfur, and amino acid metabolism, indirectly contributing to biogeochemical cycling through these AMGs. Overall, the research sheds light on the diverse and unique viral communities in the Napahai plateau wetland and emphasizes the significant roles of viruses in microbial ecology. The findings contribute to a deeper understanding of the characteristics and ecological functions of viral communities in plateau wetland ecosystems.

Circulation of adenovirus and other viruses in urban drainage channels: an environmental surveillance in Belém, Amazon region, Brazil

Urban channels in amazon cities are very polluted, with garbage and sewage disposal in these aquatic environments, favoring the high dissemination of waterborne viruses such as human adenovirus (HAdV). The aim of this study was to perform the detection and molecular characterization of adenovirus in urban channels and in a wastewater treatment plant located in a metropolitan city in the Amazon. Additionally, metagenomic analyses were performed to assess viral diversity. Samples were concentrated by organic flocculation, analyzed by quantitative real time PCR (qPCR) and sequenced (Sanger e next generation sequencing). Cell culture was performed to verify the viability of HAdV particles. A total of 104 samples were collected, being the HAdV positivity of 76% (79/104). Among the positive samples, 29.1% (23/79) were characterized as HAdV-F40 (87%, 20/23), HAdV-F41 (8.7%, 2/23), and HAdV-B (4.3%, 1/23). Average precipitation rates ranged from 163 to 614 mm, while the pH ranged from 6.9 to 7.6. Eight positive samples were inoculated into A549 cells and in 4 of these, was observed changes in the structure of the cell monolayer, alteration in the structure of the cell monolayer was observed, but without amplification when analyzed by PCR. The metagenomic data demonstrated the presence of 14 viral families, being the most abundant: Myoviridae (41% of available reads), Siphoviridae (24.5%), Podoviridae (14.1%), and Autographiviridae (6.9%) with more than 85% of the total number of identified reads. This study reinforcing that continuous surveillance may contribute to monitoring viral diversity in aquatic environments and provide early warning of potential outbreaks.

Extracellular vesicles are the main contributor to the non-viral protected extracellular sequence space

Environmental virus metagenomes, commonly referred to as "viromes", are typically generated by physically separating virus-like particles (VLPs) from the microbial fraction based on their size and mass. However, most methods used to purify VLPs, enrich extracellular vesicles (EVs) and gene transfer agents (GTAs) simultaneously. Consequently, the sequence space traditionally referred to as a "virome" contains host-associated sequences, transported via EVs or GTAs. We therefore propose to call the genetic material isolated from size-fractionated (0.22 µm) and DNase-treated samples protected environmental DNA (peDNA). This sequence space contains viral genomes, DNA transduced by viruses and DNA transported in EVs and GTAs. Since there is no genetic signature for peDNA transported in EVs, GTAs and virus particles, we rely on the successful removal of contaminating remaining cellular and free DNA when analyzing peDNA. Using marine samples collected from the North Sea, we generated a thoroughly purified peDNA dataset and developed a bioinformatic pipeline to determine the potential origin of the purified DNA. This pipeline was applied to our dataset as well as existing global marine "viromes". Through this pipeline, we identified known GTA and EV producers, as well as organisms with actively transducing proviruses as the source of the peDNA, thus confirming the reliability of our approach. Additionally, we identified novel and widespread EV producers, and found quantitative evidence suggesting that EV-mediated gene transfer plays a significant role in driving horizontal gene transfer (HGT) in the world's oceans.

Lower viral evolutionary pressure under stable versus fluctuating conditions in subzero Arctic brines

BACKGROUND

Climate change threatens Earth's ice-based ecosystems which currently offer archives and eco-evolutionary experiments in the extreme. Arctic cryopeg brine (marine-derived, within permafrost) and sea ice brine, similar in subzero temperature and high salinity but different in temporal stability, are inhabited by microbes adapted to these extreme conditions. However, little is known about their viruses (community composition, diversity, interaction with hosts, or evolution) or how they might respond to geologically stable cryopeg versus fluctuating sea ice conditions.

RESULTS

We used long- and short-read viromics and metatranscriptomics to study viruses in Arctic cryopeg brine, sea ice brine, and underlying seawater, recovering 11,088 vOTUs (~species-level taxonomic unit), a 4.4-fold increase of known viruses in these brines. More specifically, the long-read-powered viromes doubled the number of longer (≥25 kb) vOTUs generated and recovered more hypervariable regions by >5-fold compared to short-read viromes. Distribution assessment, by comparing to known viruses in public databases, supported that cryopeg brine viruses were of marine origin yet distinct from either sea ice brine or seawater viruses, while 94% of sea ice brine viruses were also present in seawater. A virus-encoded, ecologically important exopolysaccharide biosynthesis gene was identified, and many viruses (~half of metatranscriptome-inferred "active" vOTUs) were predicted as actively infecting the dominant microbial genera Marinobacter and Polaribacter in cryopeg and sea ice brines, respectively. Evolutionarily, microdiversity (intra-species genetic variations) analyses suggested that viruses within the stable cryopeg brine were under significantly lower evolutionary pressures than those in the fluctuating sea ice environment, while many sea ice brine virus-tail genes were under positive selection, indicating virus-host co-evolutionary arms races.

CONCLUSIONS

Our results confirmed the benefits of long-read-powered viromics in understanding the environmental virosphere through significantly improved genomic recovery, expanding viral discovery and the potential for biological inference. Evidence of viruses actively infecting the dominant microbes in subzero brines and modulating host metabolism underscored the potential impact of viruses on these remote and underexplored extreme ecosystems. Microdiversity results shed light on different strategies viruses use to evolve and adapt when extreme conditions are stable versus fluctuating. Together, these findings verify the value of long-read-powered viromics and provide foundational data on viral evolution and virus-microbe interactions in Earth's destabilized and rapidly disappearing cryosphere. Video Abstract.

Metagenomic comparisons reveal a highly diverse and unique viral community in a seasonally fluctuating hypersaline microbial mat

In spring 2016, a shallow hypersaline pond (50×25 m) was found in the Cuatro Cienegas Basin (CCB). This pond, known as Archaean Domes (AD) because of its elastic microbial mats that form dome-shaped structures due to the production of reducing gases reminiscent of the Archaean eon, such as methane and hydrogen sulfide, harbour a highly diverse microbial community, rich in halophilic and methanogenic archaea. AD is a seasonally fluctuating hypersaline site, with salinity ranging from low hypersaline (5.3%) during the wet season to high hypersaline (saturation) during the dry season. To characterize the viral community and to test whether it resembles those of other hypersaline sites (whose diversity is conditioned by salinity), or if it is similar to other CCB sites (with which it shares a common geological history), we generated 12 metagenomes from different seasons and depths over a 4 year period and compared them to 35 metagenomes from varied environments. Haloarchaeaviruses were detected, but were never dominant (average of 15.37 % of the total viral species), and the viral community structure and diversity were not affected by environmental fluctuations. In fact, unlike other viral communities at hypersaline sites, AD remained more diverse than other environments regardless of season. β-Diversity analyses show that AD is closely related to other CCB sites, although it has a unique viral community that forms a cluster of its own. The similarity of two surface samples to the 30 and 50 cm depth samples, as well as the observed increase in diversity at greater depths, supports the hypothesis that the diversity of CCB has evolved as a result of a long time environmental stability of a deep aquifer that functions as a 'seed bank' of great microbial diversity that is transported to the surface by sporadic groundwater upwelling events.

Novel Virus Identification through Metagenomics: A Systematic Review

Metagenomic Next Generation Sequencing (mNGS) allows the evaluation of complex microbial communities, avoiding isolation and cultivation of each microbial species, and does not require prior knowledge of the microbial sequences present in the sample. Applications of mNGS include virome characterization, new virus discovery and full-length viral genome reconstruction, either from virus preparations enriched in culture or directly from clinical and environmental specimens. Here, we systematically reviewed studies that describe novel virus identification through mNGS from samples of different origin (plant, animal and environment). Without imposing time limits to the search, 379 publications were identified that met the search parameters. Sample types, geographical origin, enrichment and nucleic acid extraction methods, sequencing platforms, bioinformatic analytical steps and identified viral families were described. The review highlights mNGS as a feasible method for novel virus discovery from samples of different origins, describes which kind of heterogeneous experimental and analytical protocols are currently used and provides useful information such as the different commercial kits used for the purification of nucleic acids and bioinformatics analytical pipelines.

[A better understanding of Earth's viruses thanks to metagenomes]

Despite their large number, viruses present in the environment remain largely unknown. Metagenomic approaches, targeting viruses specifically or not, have allowed us a better understanding of the composition of natural viral communities, with Caudoviricetes, Microviridae, Cressdnaviricota or Phycodnaviridae being the most frequently found viral groups. Metagenomes are gradually revealing the extent of the diversity of these groups and their structure, highlighting the large number of species, genera and even viral families, most of which being seen for the first time. Within these groups, the gene content, infected hosts and inhabited ecosystems are often consistent with the evolutionary history traced with marker genes. Thus, the diversity of viruses and their genes is more a reflection of their ancient origin and long coevolution with their hosts than of their ability to mutate rapidly.

Insights into the global freshwater virome

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

Nutrient levels and prokaryotes affect viral communities in plateau lakes

Viruses are the most abundant organisms in aquatic environments. Recent advances of viral metagenomic have greatly expanded our understanding of aquatic viral communities. However, little is known about the difference of viral communities and driving factors in freshwater lake. This study seeks to understand the spatio-temporal variation, differences, and driving factors of viral communities in two plateau lakes (Dianchi and Fuxian Lakes) with significant nutritional differences. The viral communities exhibited apparent seasonal variation in Dianchi Lake, while seasonal influences on the viral communities were greater than location-based influences. Two-thirds of all detected viral taxa were shared in two lakes, but there was variation in the composition of viral communities. Correlations between prokaryotic communities, environmental factors and viral communities were analyzed. The nutrients, chlorophyll a were primarily environmental parameters affecting viral communities, and the prokaryotic community was significantly correlated with the viral community. In addition, several viruses infecting humans were identified in two lakes, with the most abundant being Herpesviridae and Poxviridae. Overall, these findings provide information on the dynamics, composition, and differences of viral and prokaryotic communities in plateau lakes with different nutrient levels. These results suggest that nutritional levels and prokaryotic communities could play an important role in shaping viral communities in freshwater lakes.

Diverse viromes in polar regions: A retrospective study of metagenomic data from Antarctic animal feces and Arctic frozen soil in 2012-2014

Antarctica and the Arctic are the coldest places, containing a high diversity of microorganisms, including viruses, which are important components of polar ecosystems. However, owing to the difficulties in obtaining access to animal and environmental samples, the current knowledge of viromes in polar regions is still limited. To better understand polar viromes, this study performed a retrospective analysis using metagenomic sequencing data of animal feces from Antarctica and frozen soil from the Arctic collected during 2012-2014. The results reveal diverse communities of DNA and RNA viruses from at least 23 families from Antarctic animal feces and 16 families from Arctic soils. Although the viral communities from Antarctica and the Arctic show a large diversity, they have genetic similarities with known viruses from different ecosystems and organisms with similar viral proteins. Phylogenetic analysis of Microviridae, Parvoviridae, and Larvidaviridae was further performed, and complete genomic sequences of two novel circular replication-associated protein (rep)-encoding single-stranded (CRESS) DNA viruses closely related to Circoviridae were identified. These results reveal the high diversity, complexity, and novelty of viral communities from polar regions, and suggested the genetic similarity and functional correlations of viromes between the Antarctica and Arctic. Variations in viral families in Arctic soils, Arctic freshwater, and Antarctic soils are discussed. These findings improve our understanding of polar viromes and suggest the importance of performing follow-up in-depth investigations of animal and environmental samples from Antarctica and the Arctic, which would reveal the substantial role of these viruses in the global viral community.

A snapshot of the global drinking water virome: Diversity and metabolic potential vary with residual disinfectant use

Viruses are important drivers of microbial community ecology and evolution, influencing microbial mortality, metabolism, and horizontal gene transfer. However, the effects of viruses remain largely unknown in many environments, including in drinking water systems. Drinking water metagenomic studies have offered a whole community perspective of bacterial impacts on water quality, but have not yet considered the influences of viruses. In this study, we address this gap by mining viral DNA sequences from publicly available drinking water metagenomes from distribution systems in six countries around the world. These datasets provide a snapshot of the taxonomic diversity and metabolic potential of the global drinking water virome; and provide an opportunity to investigate the effects of geography, climate, and drinking water treatment practices on viral diversity. Both environmental conditions and differences in sample processing were found to influence the viral composition. Using free chlorine as the residual disinfectant was associated with clear differences in viral taxonomic diversity and metabolic potential, with significantly fewer viral populations and less even viral community structures than observed in distribution systems without residual disinfectant. Additionally, drinking water viruses carry antibiotic resistance genes (ARGs), as well as genes to survive oxidative stress and nitrogen limitation. Through this study, we have demonstrated that viral communities are diverse across drinking water systems and vary with the use of residual disinfectant. Our findings offer directions for future research to develop a more robust understanding of how virus-bacteria interactions in drinking water distribution systems affect water quality.

Viral Metagenomics Reveals Widely Diverse Viral Community of Freshwater Amazonian Lake

Despite the importance of understanding the ecology of freshwater viruses, there are not many studies on the subject compared to marine viruses. The microbiological interactions in these environments are still poorly known, especially between bacteriophages and their host bacteria and between cyanophages and cyanobacteria. Lake Bologna, Belém, capital of the Brazilian State of Pará, is a water source that supplies the city and its metropolitan region. However, it remains unexplored regarding the contents of its virome and viral diversity composition. Therefore, this work aims to explore the taxonomic diversity of DNA viruses in this lake, especially bacteriophages and cyanophages, since they can act as transducers of resistance genes and reporters of water quality for human consumption. We used metagenomic sequencing data generated by previous studies. We analyzed it at the taxonomic level using the tools Kraken2, Bracken, and Pavian; later, the data was assembled using Genome Detective, which performs the assembly of viruses. The results observed here suggest the existence of a widely diverse viral community and established microbial phage-regulated dynamics in Lake Bolonha. This work is the first ever to describe the virome of Lake Bolonha using a metagenomic approach based on high-throughput sequencing, as it contributes to the understanding of water-related public health concerns regarding the spreading of antibiotic resistance genes and population control of native bacteria and cyanobacteria.

Record-Breaking Rain Event Altered Estuarine Viral Assemblages

Viruses are the dominant biological entity in the ocean, play a vital role in biogeochemical cycles, and provide their hosts with novel metabolic capabilities through auxiliary metabolic genes (AMGs). Hurricane Harvey was a category 4 hurricane that made landfall on the Texas coast in 2017 and lashed the Houston area with 1.4–1.7 × 1010 m3 of rainfall. In this paper, we aim to characterize how the changes in abiotic conditions brought by Hurricane Harvey altered the viral assemblages of Galveston Bay at the taxonomic level and determine how viral ecosystem functions were altered. Metagenomes of the viruses and their hosts were sequenced from a transect in Galveston Bay over the five weeks following the storm. Our results show that the viral assemblages of Galveston Bay dramatically changed following Hurricane Harvey’s landfall. Of the abiotic parameters measured, salinity had the strongest effect on shaping the viral assemblages. In the five weeks following Hurricane Harvey, there was a steady increase of metabolic genes and putative viral infections. Our study provides the first in-depth look at how marine viral assemblages respond and recover from extreme rainfall events, which models predict will become more frequent and intense with climate change.

Particle foraging strategies promote microbial diversity in marine environments

Microbial foraging in patchy environments, where resources are fragmented into particles or pockets embedded in a large matrix, plays a key role in natural environments. In the oceans and freshwater systems, particle-associated bacteria can interact with particle surfaces in different ways: some colonize only during short transients, while others form long-lived, stable colonies. We do not yet understand the ecological mechanisms by which both short- and long-term colonizers can coexist. Here, we address this problem with a mathematical model that explains how marine populations with different detachment rates from particles can stably coexist. In our model, populations grow only while on particles, but also face the increased risk of mortality by predation and sinking. Key to coexistence is the idea that detachment from particles modulates both net growth and mortality, but in opposite directions, creating a trade-off between them. While slow-detaching populations show the highest growth return (i.e., produce more net offspring), they are more susceptible to suffer higher rates of mortality than fast-detaching populations. Surprisingly, fluctuating environments, manifesting as blooms of particles (favoring growth) and predators (favoring mortality) significantly expand the likelihood that populations with different detachment rates can coexist. Our study shows how the spatial ecology of microbes in the ocean can lead to a predictable diversification of foraging strategies and the coexistence of multiple taxa on a single growth-limiting resource.

Microbial Diversity in the Indian Ocean Sediments: An Insight into the Distribution and Associated Factors

Indian Ocean is the third largest oceanic division of the world and shelter to a huge microbial diversity. These microbes play an important role in the metabolism of carbon, sulfur, nitrogen, and phosphorus in the ocean water. They are also major contributors of carbon fixing and sequestration, as much as terrestrial plants to achieve CO₂ emissions reduction. The prokaryotic community in the East Indian Ocean primarily comprises of heterotrophic bacteria like Alphaproteobacteria and Gammaproteobacteria, followed by Firmicutes and Actinobacteria. The Arabian Sea and the Bay of Bengal are typically characterized by presence of vast areas of oxygen minimum zones (OMZs) and have been witnessing a shift in the microbial diversity due to the changing conditions in the ocean water. Several canonical correspondence analyses reveal temperature, salinity, and phosphate levels as crucial environmental factors in propelling the distribution of diazotrophs. The viral consortia are dominated by the Caudovirales, an order of tailed bacteriophages. Due to the rapid change in the environmental factors such as topography, temperature, and sunlight contributing toward climate change, their role in sustaining the chemical composition of the ocean can be drastically affected especially with the evidence of several bacterial and fungal communities responding to latitudinal and temperature change. Therefore, we aim to critically review the status of microbial diversity in Indian Ocean to predict their response toward climate change as they are the sentinels of change in marine life and to understand the dynamics of microbial communities in the various locations of Indian Ocean.

Genomic Characterization of a Novel Freshwater Cyanophage Reveals a New Lineage of Cyanopodovirus

Cyanobacteria are one of the dominant autotrophs in tropical freshwater communities, yet phages infecting them remain poorly characterized. Here we present the characterization of cyanophage S-SRP02, isolated from a tropical freshwater lake in Singapore, which infects Synechococcus sp. Strain SR-C1 isolated from the same lake. S-SRP02 represents a new evolutionary lineage of cyanophage. Out of 47 open reading frames (ORFs), only 20 ORFs share homology with genes encoding proteins of known function. There is lack of auxiliary metabolic genes which was commonly found as core genes in marine cyanopodoviruses. S-SRP02 also harbors unique structural genes highly divergent from other cultured phages. Phylogenetic analysis and viral proteomic tree further demonstrate the divergence of S-SRP02 from other sequenced phage isolates. Nonetheless, S-SRP02 shares synteny with phage genes of uncultured phages obtained from the Mediterranean Sea deep chlorophyll maximum fosmids, indicating the ecological importance of S-SRP02 and its related viruses. This is further supported by metagenomic mapping of environmental viral metagenomic reads onto the S-SRP02 genome.

Metagenomic Analysis of Viral Community in the Yangtze River Expands Known Eukaryotic and Prokaryotic Virus Diversity in Freshwater

Viruses in aquatic ecosystems are characterized by extraordinary abundance and diversity. Thus far, there have been limited studies focused on viral communities in river water systems. Here, we investigated the virome of the Yangtze River Delta using viral metagenomic analysis. The compositions of viral communities from six sampling sites were analyzed and compared. By using library construction and next generation sequencing, contigs and singlet reads similar to viral sequences were classified into 17 viral families, including nine dsDNA viral families, four ssDNA viral families and four RNA viral families. Statistical analysis using Friedman test suggested that there was no significant difference among the six sampling sites (P ​> ​0.05). The viromes in this study were all dominated by the order Caudovirales, and a group of Freshwater phage uvFW species were particularly prevalent among all the samples. The virome from Nanjing presented a unique pattern of viral community composition with a relatively high abundance of family Parvoviridae. Phylogenetic analyses based on virus hallmark genes showed that the Caudovirales order and CRESS-DNA viruses presented high genetic diversity, while viruses in the Microviridae and Parvoviridae families and the Riboviria realm were relatively conservative. Our study provides the first insight into viral community composition in large river ecosystem, revealing the diversity and stability of river water virome, contributing to the proper utilization of freshwater resource.

•

First insight into viral community composition in large river ecosystem.

•

Virus hallmark genes present both diverse and conservative characteristics.

•

The composition of viral communities is similar on the whole.

•

Slight regional variation of virome is existed in individual areas.

Unique T4-like phages in high-altitude lakes above 4500 m on the Tibetan Plateau

Viruses are the most abundant biological entities in the biosphere; however, little is known about viral ecology in high altitude lakes. Here, we characterized viruses from 13 lakes, nine of which located ≥4500 m above sea level, on the Tibetan Plateau, the highest plateau on Earth. The abundance of virus-like particle (VLP) in Tibetan lakes ranged from 4.8 ± 0.2 × 10⁵ VLPs mL^-1 to 6.0 ± 0.2 × 10⁷ VLPs mL^-1 and the virus-to-bacterium ratio was in the lower range of values reported for other lakes. The viral population size was positively correlated with turbidity and negatively correlated with particulate organic carbon concentration. Highly diverse VLP morphologies, including large (~300 nm) morphotypes, were observed. Phylogenetic analysis of T4-like bacteriophages based on major capsid gene (g23) identified a novel viral group, which were detected in abundance in hyposaline and mesosaline Tibetan lakes. Adaptation to lake evolution, water source (glacier-fed or non-glacier-fed) and environmental conditions (e.g., salinity, phosphorus concentration and productivity) are likely responsible for the variation in T4-like myovirus community composition in contrasting Tibetan lakes. This first investigation of viruses in high-altitude alpine lakes above 4500 m could contribute to our understanding of viral ecology in global alpine lakes.

Bacteriophages in Biological Wastewater Treatment Systems: Occurrence, Characterization, and Function

Phage bacteria interactions can affect structure, dynamics, and function of microbial communities. In the context of biological wastewater treatment (BWT), the presence of phages can alter the efficiency of the treatment process and influence the quality of the treated effluent. The active role of phages in BWT has been demonstrated, but many questions remain unanswered regarding the diversity of phages in these engineered environments, the dynamics of infection, the determination of bacterial hosts, and the impact of their activity in full-scale processes. A deeper understanding of the phage ecology in BWT can lead the improvement of process monitoring and control, promote higher influent quality, and potentiate the use of phages as biocontrol agents. In this review, we highlight suitable methods for studying phages in wastewater adapted from other research fields, provide a critical overview on the current state of knowledge on the effect of phages on structure and function of BWT bacterial communities, and highlight gaps, opportunities, and priority questions to be addressed in future research.

Novel Freshwater Cyanophages Provide New Insights into Evolutionary Relationships between Freshwater and Marine Cyanophages

Cyanobacteria and cyanophages are present widely in both freshwater and marine environments. However, freshwater cyanophages remain unknown largely due to the small numbers of cyanophage isolates despite their ecological and environmental significance. In this study, we present the characterization of two novel lytic freshwater cyanophages isolated from a tropical inland lake in Singapore, namely, cyanopodovirus S-SRP01 and cyanomyovirus S-SRM01, infecting two different strains of Synechococcus spp. Functional annotation of S-SRP01 and S-SRM01 genomes revealed a high degree of homology with marine cyanophages. Phylogenetic trees of concatenated genes and whole-genome alignment provided further evidence that S-SRP01 is close evolutionarily to marine cyanopodoviruses, while S-SRM01 is evolutionarily close to marine cyanomyoviruses. Few genetic similarities between freshwater and marine cyanophages have been identified in previous studies. The isolation of S-SRP01 and S-SRM01 expand current knowledge on freshwater cyanophages infecting Synechococcus spp. Their high degree of gene sharing provides new insights into the evolutionary relationships between freshwater and marine cyanophages. This relatedness is further supported by the discovery of similar phenomenon from other freshwater viral metagenomes. IMPORTANCE This study expands the current knowledge on freshwater cyanophage isolates and cyanophage genetic diversity, indicating that freshwater and marine cyanophages infecting Synechococcus spp. may share close genetic similarity and evolutionary relationships.

Adaptation of the polony technique to quantify Gokushovirinae, a diverse group of single-stranded DNA phage

Advances in metagenomics have revealed the ubiquity of single-stranded DNA (ssDNA) phage belonging to the subfamily Gokushovirinae in the oceans; however, the abundance and ecological roles of this group are unknown. Here, we quantify gokushoviruses through adaptation of the polony method, in which viral template DNA is immobilized in a gel, amplified by PCR, and subsequently detected by hybridization. Primers and probes for this assay were designed based on PCR amplicon diversity of gokushovirus major capsid protein gene sequences from a depth profile in the Gulf of Aqaba, Red Sea sampled in September 2015. At ≥95% identity, these 87 gokushovirus sequences formed 14 discrete clusters with the largest clades showing distinct depth distributions. The application of the polony method enabled the first quantification of gokushoviruses in any environment. The gokushoviruses were most abundant in the upper 40 m of the stratified water column, with a subsurface peak in abundance of 1.26 × 10⁵ viruses ml^-1 . These findings suggest that discrete gokushovirus genotypes infect bacterial hosts that differentially partition in the water column. Since the designed primers and probe are conserved across marine ecosystems, this polony method can be applied broadly for the quantification of gokushoviruses throughout the global oceans.

Legionella pneumophila CRISPR-Cas Suggests Recurrent Encounters with One or More Phages in the Family Microviridae

Legionella pneumophila is a ubiquitous freshwater pathogen and the causative agent of Legionnaires' disease. L. pneumophila growth within protists provides a refuge from desiccation, disinfection, and other remediation strategies. One outstanding question has been whether this protection extends to phages. L. pneumophila isolates are remarkably devoid of prophages and to date no Legionella phages have been identified. Nevertheless, many L. pneumophila isolates maintain active CRISPR-Cas defenses. So far, the only known target of these systems is an episomal element that we previously named Legionella mobile element 1 (LME-1). The continued expansion of publicly available genomic data promises to further our understanding of the role of these systems. We now describe over 150 CRISPR-Cas systems across 600 isolates to establish the clearest picture yet of L. pneumophila's adaptive defenses. By searching for targets of 1,500 unique CRISPR-Cas spacers, LME-1 remains the only identified CRISPR-Cas targeted integrative element. We identified 3 additional LME-1 variants-all targeted by previously and newly identified CRISPR-Cas spacers-but no other similar elements. Notably, we also identified several spacers with significant sequence similarity to microviruses, specifically those within the subfamily Gokushovirinae. These spacers are found across several different CRISPR-Cas arrays isolated from geographically diverse isolates, indicating recurrent encounters with these phages. Our analysis of the extended Legionella CRISPR-Cas spacer catalog leads to two main conclusions: current data argue against CRISPR-Cas targeted integrative elements beyond LME-1, and the heretofore unknown L. pneumophila phages are most likely lytic gokushoviruses. IMPORTANCE Legionnaires' disease is an often-fatal pneumonia caused by Legionella pneumophila, which normally grows inside amoebae and other freshwater protists. L. pneumophila trades diminished access to nutrients for the protection and isolation provided by the host. One outstanding question is whether L. pneumophila is susceptible to phages, given the protection provided by its intracellular lifestyle. In this work, we use Legionella CRISPR spacer sequences as a record of phage infection to predict that the "missing" L. pneumophila phages belong to the microvirus subfamily Gokushovirinae. Gokushoviruses are known to infect another intracellular pathogen, Chlamydia. How do gokushoviruses access L. pneumophila (and Chlamydia) inside their "cozy niches"? Does exposure to phages happen during a transient extracellular period (during cell-to-cell spread) or is it indicative of a more complicated environmental lifestyle? One thing is clear, 100 years after their discovery, phages continue to hold important secrets about the bacteria upon which they prey.

Glacier ice archives nearly 15,000-year-old microbes and phages

BACKGROUND

Glacier ice archives information, including microbiology, that helps reveal paleoclimate histories and predict future climate change. Though glacier-ice microbes are studied using culture or amplicon approaches, more challenging metagenomic approaches, which provide access to functional, genome-resolved information and viruses, are under-utilized, partly due to low biomass and potential contamination.

RESULTS

We expand existing clean sampling procedures using controlled artificial ice-core experiments and adapted previously established low-biomass metagenomic approaches to study glacier-ice viruses. Controlled sampling experiments drastically reduced mock contaminants including bacteria, viruses, and free DNA to background levels. Amplicon sequencing from eight depths of two Tibetan Plateau ice cores revealed common glacier-ice lineages including Janthinobacterium, Polaromonas, Herminiimonas, Flavobacterium, Sphingomonas, and Methylobacterium as the dominant genera, while microbial communities were significantly different between two ice cores, associating with different climate conditions during deposition. Separately, ~355- and ~14,400-year-old ice were subject to viral enrichment and low-input quantitative sequencing, yielding genomic sequences for 33 vOTUs. These were virtually all unique to this study, representing 28 novel genera and not a single species shared with 225 environmentally diverse viromes. Further, 42.4% of the vOTUs were identifiable temperate, which is significantly higher than that in gut, soil, and marine viromes, and indicates that temperate phages are possibly favored in glacier-ice environments before being frozen. In silico host predictions linked 18 vOTUs to co-occurring abundant bacteria (Methylobacterium, Sphingomonas, and Janthinobacterium), indicating that these phages infected ice-abundant bacterial groups before being archived. Functional genome annotation revealed four virus-encoded auxiliary metabolic genes, particularly two motility genes suggest viruses potentially facilitate nutrient acquisition for their hosts. Finally, given their possible importance to methane cycling in ice, we focused on Methylobacterium viruses by contextualizing our ice-observed viruses against 123 viromes and prophages extracted from 131 Methylobacterium genomes, revealing that the archived viruses might originate from soil or plants.

CONCLUSIONS

Together, these efforts further microbial and viral sampling procedures for glacier ice and provide a first window into viral communities and functions in ancient glacier environments. Such methods and datasets can potentially enable researchers to contextualize new discoveries and begin to incorporate glacier-ice microbes and their viruses relative to past and present climate change in geographically diverse regions globally. Video Abstract.

Ecogenomics of Groundwater Phages Suggests Niche Differentiation Linked to Specific Environmental Tolerance

Viruses are ubiquitous microbiome components, shaping ecosystems via strain-specific predation, horizontal gene transfer and redistribution of nutrients through host lysis. Viral impacts are important in groundwater ecosystems, where microbes drive many nutrient fluxes and metabolic processes; however, little is known about the diversity of viruses in these environments. We analyzed four groundwater plasmidomes (the entire plasmid content of an environment) and identified 200 viral sequences, which clustered into 41 genus-level viral clusters (approximately equivalent to viral genera) including 9 known and 32 putative new genera. We used publicly available bacterial whole-genome sequences (WGS) and WGS from 261 bacterial isolates from this groundwater environment to identify potential viral hosts. We linked 76 of the 200 viral sequences to a range of bacterial phyla, the majority associated with Proteobacteria, followed by Firmicutes, Bacteroidetes, and Actinobacteria. The publicly available WGS enabled mapping bacterial hosts to several viral sequences. The WGS of groundwater isolates increased the depth of host prediction by allowing host identification at the strain level. The latter included 4 viruses that were almost entirely (>99% query coverage, >99% identity) identified as integrated in the genomes of Pseudomonas, Acidovorax, and Castellaniella strains, resulting in high-confidence host assignments. Lastly, 21 of these viruses carried putative auxiliary metabolite genes for metal and antibiotic resistance, which might drive their infection cycles and/or provide selective advantage to infected hosts. Exploring the groundwater virome provides a necessary foundation for integration of viruses into ecosystem models where they are key players in microbial adaption to environmental stress. IMPORTANCE To our knowledge, this is the first study to identify the bacteriophage distribution in a groundwater ecosystem shedding light on their prevalence and distribution across metal-contaminated and background sites. Our study is uniquely based on selective sequencing of solely the extrachromosomal elements of a microbiome followed by analysis for viral signatures, thus establishing a more focused approach for phage identifications. Using this method, we detected several novel phage genera along with those previously established. Our approach of using the whole-genome sequences of hundreds of bacterial isolates from the same site enabled us to make host assignments with high confidence, several at strain levels. Certain phage genes suggest that they provide an environment-specific selective advantage to their bacterial hosts. Our study lays the foundation for future research on directed phage isolations using specific bacterial host strains to further characterize groundwater phages, their life cycles, and their effects on groundwater microbiome and biogeochemistry.

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.

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.

Ten Years of Collaborative Progress in the Quest for Orthologs

Accurate determination of the evolutionary relationships between genes is a foundational challenge in biology. Homology-evolutionary relatedness-is in many cases readily determined based on sequence similarity analysis. By contrast, whether or not two genes directly descended from a common ancestor by a speciation event (orthologs) or duplication event (paralogs) is more challenging, yet provides critical information on the history of a gene. Since 2009, this task has been the focus of the Quest for Orthologs (QFO) Consortium. The sixth QFO meeting took place in Okazaki, Japan in conjunction with the 67th National Institute for Basic Biology conference. Here, we report recent advances, applications, and oncoming challenges that were discussed during the conference. Steady progress has been made toward standardization and scalability of new and existing tools. A feature of the conference was the presentation of a panel of accessible tools for phylogenetic profiling and several developments to bring orthology beyond the gene unit-from domains to networks. This meeting brought into light several challenges to come: leveraging orthology computations to get the most of the incoming avalanche of genomic data, integrating orthology from domain to biological network levels, building better gene models, and adapting orthology approaches to the broad evolutionary and genomic diversity recognized in different forms of life and viruses.

Occurrence and diversity of viruses associated with cyanobacterial communities in a Brazilian freshwater reservoir

As part of the phytoplankton of marine and freshwater environments around the world, cyanobacteria interact with viruses (cyanophages) that affect their abundance and diversity. Investigations focusing on cyanophages co-occurring with freshwater cyanobacteria are scarce, particularly in Brazil. The aim of this study was to assess the diversity of cyanophages associated with a Microcystis-dominated cyanobacterial bloom in a tropical reservoir. Samples were processed as viral fractions of water and cellular fractions, and temporal fluctuations in the abundance of Ma-LMM01-type cyanophages and their Microcystis hosts were determined by qPCR. We applied shotgun metagenomics to obtain a wider characterization of the cyanophage community. During the study period, Microcystis gene copies were quantified in all cellular fractions, and the copy number of the Ma-LMM01 phage gene tended to increase with host abundance. Metagenomic analysis demonstrated that Caudovirales was the major viral order associated with the cyanophage families Myoviridae (34-88%), Podoviridae (3-42%), and Siphoviridae (6-23%). The metagenomic analysis results confirmed the presence of Microcystis cyanophages in both viral and cellular fractions and demonstrated a high relative abundance of picocyanobacteria-related viruses and Prochlorococcus (36-52%) and Synechococcus (37-50%) phages. For other main cyanobacterial genera, no related cyanophages were identified, which was probably due to the scarce representation of cyanophage sequences in databanks. Thus, the studied reservoir hosted a diverse cyanophage community with a remarkable contribution of phages related to picoplanktonic cyanobacteria. These results provide insights that motivate future sequencing efforts to assess cyanophage diversity and recover complete genomes.

Genome Sequences of Microviruses Identified in Gila Monster Feces

The complete genome sequences of 33 microviruses were determined from fecal samples collected from 14 Arizona-dwelling Gila monsters using high-throughput sequencing. These microviruses with genomes 4,383 to 6,782 nucleotides (nt) long were broadly distributed across the 14 samples.

The complete genome sequences of 33 microviruses were determined from fecal samples collected from 14 Arizona-dwelling Gila monsters using high-throughput sequencing. These microviruses with genomes 4,383 to 6,782 nucleotides (nt) long were broadly distributed across the 14 samples.

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.

Sediment Metagenomes as Time Capsules of Lake Microbiomes

The reconstruction of ecological time series from lake sediment archives can retrace the environmental impact of human activities. Molecular genetic approaches in paleolimnology have provided unprecedented access to DNA time series, which record evidence of the microbial ecologies that underlaid historical lake ecosystems. Such studies often rely on single-gene surveys, and consequently, the full diversity of preserved microorganisms remains unexplored. In this study, we probed the diversity archived in contemporary and preindustrial sediments by comparative shotgun metagenomic analysis of surface water and sediment samples from three eastern Canadian lakes. In a strategy that was aimed at disentangling historical DNA from the indigenous sediment background, microbial preservation signals were captured by mapping sequence similarities between sediment metagenome reads and reference surface water metagenome assemblies. We detected preserved Cyanobacteria, diverse bacterioplankton, microeukaryotes, and viruses in sediment metagenomes. Among the preserved microorganisms were important groups never before reported in paleolimnological reconstructions, including bacteriophages (Caudovirales) and ubiquitous freshwater Betaproteobacteria (Polynucleobacter and Limnohabitans). In contrast, ultramicroscopic Actinobacteria ("Candidatus Nanopelagicales") and Alphaproteobacteria (Pelagibacterales) were apparently not well preserved in sediment metagenomes even though they were numerically dominant in surface water metagenomes. Overall, our study explored a novel application of whole-metagenome shotgun sequencing for discovering the DNA remains of a broad diversity of microorganisms preserved in lake sediments. The recovery of diverse microbial time series supports the taxonomic expansion of microbiome reconstructions and the development of novel microbial paleoindicators.IMPORTANCE Lakes are critical freshwater resources under mounting pressure from climate change and other anthropogenic stressors. The reconstruction of ecological time series from sediment archives with paleolimnological techniques has been shown to be an effective means of understanding how humans are modifying lake ecosystems over extended timescales. In this study, we combined shotgun DNA sequencing with a novel comparative analysis of surface water and sediment metagenomes to expose the diversity of microorganisms preserved in lake sediments. The detection of DNA from a broad diversity of preserved microbes serves to more fully reconstruct historical microbiomes and describe preimpact lake conditions.

Seasonal dynamics in taxonomy and function within bacterial and viral metagenomic assemblages recovered from a freshwater agricultural pond

BACKGROUND

Ponds are important freshwater habitats that support both human and environmental activities. However, relative to their larger counterparts (e.g. rivers, lakes), ponds are understudied, especially with regard to their microbial communities. Our study aimed to fill this knowledge gap by using culture-independent, high-throughput sequencing to assess the dynamics, taxonomy, and functionality of bacterial and viral communities in a freshwater agricultural pond.

RESULTS

Water samples (n = 14) were collected from a Mid-Atlantic agricultural pond between June 2017 and May 2018 and filtered sequentially through 1 and 0.2 μm filter membranes. Total DNA was then extracted from each filter, pooled, and subjected to 16S rRNA gene and shotgun sequencing on the Illumina HiSeq 2500 platform. Additionally, on eight occasions water filtrates were processed for viral metagenomes (viromes) using chemical concentration and then shotgun sequenced. A ubiquitous freshwater phylum, Proteobacteria was abundant at all sampling dates throughout the year. However, environmental characteristics appeared to drive the structure of the community. For instance, the abundance of Cyanobacteria (e.g. Nostoc) increased with rising water temperatures, while a storm event appeared to trigger an increase in overall bacterial diversity, as well as the relative abundance of Bacteroidetes. This event was also associated with an increase in the number of antibiotic resistance genes. The viral fractions were dominated by dsDNA of the order Caudovirales, namely Siphoviridae and Myovirdae.

CONCLUSIONS

Overall, this study provides one of the largest datasets on pond water microbial ecology to date, revealing seasonal trends in the microbial taxonomic composition and functional potential.

Metagenomic Insights into the Sewage RNA Virosphere of a Large City

Sewage-associated viruses can cause several human and animal diseases, such as gastroenteritis, hepatitis, and respiratory infections. Therefore, their detection in wastewater can reflect current infections within the source population. To date, no viral study has been performed using the sewage of any large South American city. In this study, we used viral metagenomics to obtain a single sample snapshot of the RNA virosphere in the wastewater from Santiago de Chile, the seventh largest city in the Americas. Despite the overrepresentation of dsRNA viruses, our results show that Santiago’s sewage RNA virosphere was composed mostly of unknown sequences (88%), while known viral sequences were dominated by viruses that infect bacteria (60%), invertebrates (37%) and humans (2.4%). Interestingly, we discovered three novel genogroups within the Picobirnaviridae family that can fill major gaps in this taxa’s evolutionary history. We also demonstrated the dominance of emerging Rotavirus genotypes, such as G8 and G6, that have displaced other classical genotypes, which is consistent with recent clinical reports. This study supports the usefulness of sewage viral metagenomics for public health surveillance. Moreover, it demonstrates the need to monitor the viral component during the wastewater treatment and recycling process, where this virome can constitute a reservoir of human pathogens.

Virus and Potential Host Microbes from Viral-Enriched Metagenomic Characterization in the High-Altitude Wetland, Salar de Huasco, Chile

Salar de Huasco is a wetland in the Andes mountains, located 3800 m above sea level at the Chilean Altiplano. Here we present a study aimed at characterizing the viral fraction and the microbial communities through metagenomic analysis. Two ponds (H0 and H3) were examined in November 2015. Water samples were processed using tangential flow filtration to obtain metagenomes from which the DNA fraction of the sample was amplified and sequenced (HiSeq system, Illumina). The ponds were characterized by freshwater and the viral-like particles to picoplankton ratio was 12.1 and 2.3 for H0 and H3, respectively. A great number of unassigned viral sequences were found in H0 (55.8%) and H3 (32.8%), followed by the family Fuselloviridae 20.8% (H0) and other less relatively abundant groups such as Microviridae (H0, 11.7% and H3, 3.3%) and Inoviridae (H3, 2.7%). The dominant viral sequences in both metagenomes belong to the order Caudovirales, with Siphoviridae being the most important family, especially in H3 (32.7%). The most important bacteria phyla were Proteobacteria, Bacteroidetes and Firmicutes in both sites, followed by Cyanobacteria (H0). Genes encoding lysogenic and lytic enzymes (i.e., recombinases and integrases) were found in H0 and H3, indicating a potential for active viral replication at the time of sampling; this was supported by the presence of viral metabolic auxiliary genes at both sites (e.g., cysteine hydrolase). In total, our study indicates a great novelty of viral groups, differences in taxonomic diversity and replication pathways between sites, which contribute to a better understanding of how viruses balance the cycling of energy and matter in this extreme environment.

Estimate of the diversity of viral and bacterial assemblage in the coastal water of Lake Baikal

In this study, we analysed the diversity and composition of double-stranded DNA viral and bacterial communities within the sample of surface coastal water of Southern Baikal through metagenomics and deep sequencing of the 16S ribosomal RNA gene, respectively. The 16S rRNA gene analysis has revealed 14 phyla and dominance of the 'Actinobacteria' (43.6%), 'Proteobacteria' (25.2%) and 'Bacteroidetes' (11.5%). The bacterial composition was similar to that obtained previously in Lake Baikal littoral zone. Out of 1 030 169 processed virome reads, 37.4% of sequences (385 421) were identified as viral; 15.1% were identified as nonviral and related to the domains Eukarya, Bacteria and Archaea; and 47.5% had no matches in the databases. The identified virotypes belonged to different families and were predicted to infect a wide range of organisms, from bacteria to mammals. Six families (Myoviridae, Poxviridae, Mimiviridae, Siphoviridae, Phycodnaviridae and Podoviridae) were dominant accounting for more than 90% of the identified sequences (48.3%, 17.4%, 8.3%, 6.8%, 5.8% and 4.1%, respectively). In contrast to other freshwater systems, high percentage of the Poxviridae and Mimiviridae was recorded in the water sample of Lake Baikal.

Freshwater viral metagenome reveals novel and functional phage-borne antibiotic resistance genes

BACKGROUND

Antibiotic resistance developed by bacteria is a significant threat to global health. Antibiotic resistance genes (ARGs) spread across different bacterial populations through multiple dissemination routes, including horizontal gene transfer mediated by bacteriophages. ARGs carried by bacteriophages are considered especially threatening due to their prolonged persistence in the environment, fast replication rates, and ability to infect diverse bacterial hosts. Several studies employing qPCR and viral metagenomics have shown that viral fraction and viral sequence reads in clinical and environmental samples carry many ARGs. However, only a few ARGs have been found in viral contigs assembled from metagenome reads, with most of these genes lacking effective antibiotic resistance phenotypes. Owing to the wide application of viral metagenomics, nevertheless, different classes of ARGs are being continuously found in viral metagenomes acquired from diverse environments. As such, the presence and functionality of ARGs encoded by bacteriophages remain up for debate.

RESULTS

We evaluated ARGs excavated from viral contigs recovered from urban surface water viral metagenome data. In virome reads and contigs, diverse ARGs, including polymyxin resistance genes, multidrug efflux proteins, and β-lactamases, were identified. In particular, when a lenient threshold of e value of ≤ 1 × e-5 and query coverage of ≥ 60% were employed in the Resfams database, the novel β-lactamases blaHRV-1 and blaHRVM-1 were found. These genes had unique sequences, forming distinct clades of class A and subclass B3 β-lactamases, respectively. Minimum inhibitory concentration analyses for E. coli strains harboring blaHRV-1 and blaHRVM-1 and catalytic kinetics of purified HRV-1 and HRVM-1 showed reduced susceptibility to penicillin, narrow- and extended-spectrum cephalosporins, and carbapenems. These genes were also found in bacterial metagenomes, indicating that they were harbored by actively infecting phages.

CONCLUSION

Our results showed that viruses in the environment carry as-yet-unreported functional ARGs, albeit in small quantities. We thereby suggest that environmental bacteriophages could be reservoirs of widely variable, unknown ARGs that could be disseminated via virus-host interactions. Video abstract.

Extreme Viral Partitioning in a Marine-Derived High Arctic Lake

The Arctic is warming at an accelerating pace, and the rise in temperature has increasing impacts on the Arctic biome. Lakes are integrators of their surroundings and thus excellent sentinels of environmental change. Despite their importance in the regulation of key microbial processes, viruses remain largely uncharacterized in Arctic lacustrine environments. We sampled a highly stratified meromictic lake near the northern limit of the Canadian High Arctic, a region in rapid transition due to climate change. We found that the different layers of the lake harbored viral communities that were strikingly dissimilar and highly divergent from known viruses. Viruses were more abundant in the deepest part of the lake containing ancient Arctic Ocean seawater that was trapped during glacial retreat and were genomically unlike any viruses previously described. This research demonstrates the complexity and novelty of viral communities in an environment that is vulnerable to ongoing perturbation.

High-latitude, perennially stratified (meromictic) lakes are likely to be especially vulnerable to climate warming because of the importance of ice in maintaining their water column structure and associated distribution of microbial communities. This study aimed to characterize viral abundance, diversity, and distribution in a meromictic lake of marine origin on the far northern coast of Ellesmere Island, in the Canadian High Arctic. We collected triplicate samples for double-stranded DNA (dsDNA) viromics from five depths that encompassed the major features of the lake, as determined by limnological profiling of the water column. Viral abundance and virus-to-prokaryote ratios were highest at greater depths, while bacterial and cyanobacterial counts were greatest in the surface waters. The viral communities from each zone of the lake defined by salinity, temperature, and dissolved oxygen concentrations were markedly distinct, suggesting that there was little exchange of viral types among lake strata. Ten viral assembled genomes were obtained from our libraries, and these also segregated with depth. This well-defined structure of viral communities was consistent with that of potential hosts. Viruses from the monimolimnion, a deep layer of ancient Arctic Ocean seawater, were more diverse and relatively abundant, with few similarities to available viral sequences. The Lake A viral communities also differed from published records from the Arctic Ocean and meromictic Ace Lake in Antarctica. This first characterization of viral diversity from this sentinel environment underscores the microbial richness and complexity of an ecosystem type that is increasingly exposed to major perturbations in the fast-changing Arctic.

IMPORTANCE The Arctic is warming at an accelerating pace, and the rise in temperature has increasing impacts on the Arctic biome. Lakes are integrators of their surroundings and thus excellent sentinels of environmental change. Despite their importance in the regulation of key microbial processes, viruses remain largely uncharacterized in Arctic lacustrine environments. We sampled a highly stratified meromictic lake near the northern limit of the Canadian High Arctic, a region in rapid transition due to climate change. We found that the different layers of the lake harbored viral communities that were strikingly dissimilar and highly divergent from known viruses. Viruses were more abundant in the deepest part of the lake containing ancient Arctic Ocean seawater that was trapped during glacial retreat and were genomically unlike any viruses previously described. This research demonstrates the complexity and novelty of viral communities in an environment that is vulnerable to ongoing perturbation.

Network Analysis of the Papaya Orchard Virome from Two Agroecological Regions of Chiapas, Mexico

Virus-virus interactions in plants can modify host symptoms. As a result, disease management strategies may be unsuccessful if they are based solely on visual assessment and diagnostic assays for known individual viruses. Papaya ringspot virus is an important limiting factor for papaya production and likely has interactions with other viruses that are not yet known. Using high-throughput sequencing, we recovered known and novel RNA and DNA viruses from papaya orchards in Chiapas, Mexico, and categorized them by host and, in the case of papaya, symptom type: asymptomatic papaya, papaya with ringspot virus symptoms, papaya with nonringspot symptoms, weeds, and insects. Using network analysis, we demonstrated virus associations within and among host types and described the ecological community patterns. Recovery of viruses from weeds and asymptomatic papaya suggests the need for additional management attention. These analyses contribute to the understanding of the community structure of viruses in the agroecological landscape.

The study of complex ecological interactions, such as those among host, pathogen, and vector communities, can help to explain host ranges and the emergence of novel pathogens. We evaluated the viromes of papaya orchards, including weed and insect viromes, to identify common viruses in intensive production of papaya in the Pacific Coastal Plain and the Central Depression of Chiapas, Mexico. Samples of papaya cultivar Maradol, susceptible to papaya ringspot virus (PRSV), were categorized by symptoms by local farmers (papaya ringspot symptoms, non-PRSV symptoms, or asymptomatic). These analyses revealed the presence of 61 viruses, where only 4 species were shared among both regions, 16 showed homology to known viruses, and 36 were homologous with genera including Potyvirus, Comovirus, and Tombusvirus (RNA viruses) and Begomovirus and Mastrevirus (DNA viruses). We analyzed the network of associations between viruses and host-location combinations, revealing ecological properties of the network, such as an asymmetric nested pattern, and compared the observed network to null models of network association. Understanding the network structure informs management strategies, for example, revealing the potential role of PRSV in asymptomatic papaya and that weeds may be an important pathogen reservoir. We identify three key management implications: (i) each region may need a customized management strategy; (ii) visual assessment of papaya may be insufficient for PRSV, requiring diagnostic assays; and (iii) weed control within orchards may reduce the risk of virus spread to papaya. Network analysis advances understanding of host-pathogen interactions in the agroecological landscape.

IMPORTANCE Virus-virus interactions in plants can modify host symptoms. As a result, disease management strategies may be unsuccessful if they are based solely on visual assessment and diagnostic assays for known individual viruses. Papaya ringspot virus is an important limiting factor for papaya production and likely has interactions with other viruses that are not yet known. Using high-throughput sequencing, we recovered known and novel RNA and DNA viruses from papaya orchards in Chiapas, Mexico, and categorized them by host and, in the case of papaya, symptom type: asymptomatic papaya, papaya with ringspot virus symptoms, papaya with nonringspot symptoms, weeds, and insects. Using network analysis, we demonstrated virus associations within and among host types and described the ecological community patterns. Recovery of viruses from weeds and asymptomatic papaya suggests the need for additional management attention. These analyses contribute to the understanding of the community structure of viruses in the agroecological landscape.

Metagenomic Analysis of Virioplankton from the Pelagic Zone of Lake Baikal

This study describes two viral communities from the world’s oldest lake, Lake Baikal. For the analysis, we chose under-ice and late spring periods of the year as the most productive for Lake Baikal. These periods show the maximum seasonal biomass of phytoplankton and bacterioplankton, which are targets for viruses, including bacteriophages. At that time, the main group of viruses were tailed bacteriophages of the order Caudovirales that belong to the families Myoviridae, Siphoviridae and Podoviridae. Annotation of functional genes revealed that during the under-ice period, the “Phages, Prophages, Transposable Elements and Plasmids” (27.4%) category represented the bulk of the virome. In the late spring period, it comprised 9.6% of the virome. We assembled contigs by two methods: Separately assembled in each virome or cross-assembled. A comparative analysis of the Baikal viromes with other aquatic environments indicated a distribution pattern by soil, marine and freshwater groups. Viromes of lakes Baikal, Michigan, Erie and Ontario form the joint World’s Largest Lakes clade.

Phage-centric ecological interactions in aquatic ecosystems revealed through ultra-deep metagenomics

The persistent inertia in the ability to culture environmentally abundant microbes from aquatic ecosystems represents an obstacle in disentangling the complex web of ecological interactions spun by a diverse assortment of participants (pro- and eukaryotes and their viruses). In aquatic microbial communities, the numerically most abundant actors, the viruses, remain the most elusive, and especially in freshwaters their identities and ecology remain unknown. Here, using ultra-deep metagenomic sequencing from pelagic freshwater habitats, we recovered complete genomes of > 2000 phages, including small "miniphages" and large "megaphages" infecting iconic freshwater prokaryotic lineages. For instance, abundant freshwater Actinobacteria support infection by a very broad size range of phages (13-200 Kb). We describe many phages encoding genes that likely afford protection to their host from reactive oxygen species (ROS) in the aquatic environment and in the oxidative burst in protist phagolysosomes (phage-mediated ROS defense). Spatiotemporal abundance analyses of phage genomes revealed evanescence as the primary dynamic in upper water layers, where they displayed short-lived existences. In contrast, persistence was characteristic for the deeper layers where many identical phage genomes were recovered repeatedly. Phage and host abundances corresponded closely, with distinct populations displaying preferential distributions in different seasons and depths, closely mimicking overall stratification and mixis.

Phage-centric ecological interactions in aquatic ecosystems revealed through ultra-deep metagenomics

The persistent inertia in the ability to culture environmentally abundant microbes from aquatic ecosystems represents an obstacle in disentangling the complex web of ecological interactions spun by a diverse assortment of participants (pro- and eukaryotes and their viruses). In aquatic microbial communities, the numerically most abundant actors, the viruses, remain the most elusive, and especially in freshwaters their identities and ecology remain unknown. Here, using ultra-deep metagenomic sequencing from pelagic freshwater habitats, we recovered complete genomes of > 2000 phages, including small "miniphages" and large "megaphages" infecting iconic freshwater prokaryotic lineages. For instance, abundant freshwater Actinobacteria support infection by a very broad size range of phages (13-200 Kb). We describe many phages encoding genes that likely afford protection to their host from reactive oxygen species (ROS) in the aquatic environment and in the oxidative burst in protist phagolysosomes (phage-mediated ROS defense). Spatiotemporal abundance analyses of phage genomes revealed evanescence as the primary dynamic in upper water layers, where they displayed short-lived existences. In contrast, persistence was characteristic for the deeper layers where many identical phage genomes were recovered repeatedly. Phage and host abundances corresponded closely, with distinct populations displaying preferential distributions in different seasons and depths, closely mimicking overall stratification and mixis.