The marine viromes of four oceanic regions

Viral communities locked in high elevation permafrost up to 100 m in depth on the Tibetan Plateau

Permafrost serves as a natural cold reservoir for viral communities. However, little is known about the viromes in deep permafrost soil, as most studies of permafrost were restricted to shallow areas. Here, permafrost soil samples of up to 100 m in depth were collected from two sites in the Tuotuo River permafrost area on the Tibetan Plateau. We investigated the viral composition in these permafrost soil samples and analyzed the relationship of viral composition and diversity along with depths. Our study revealed that greater permafrost thickness corresponds to higher diversity within the viral community. Bacteriophages were found to be the dominant viral communities, with "kill the winner" dynamics observed within the Siphoviridae and Myoviridae. The abundance and diversity of viral communities may follow a potential pattern along soil layers and depths, influenced by pH, trace elements, and permafrost thickness. Notably, strong correlations were discovered between the content of inorganic elements, including B, Mg, Cr, Bi, Ti, Na, Ni, and Cu, and the viral composition. Moreover, we discovered highly conserved sequences of giant viruses at depth of 10, 20, and 50 m in permafrost, which play a crucial role in evolutionary processes. These findings provide valuable insights into the viral community patterns from shallow to 100-m-depth in high-elevation permafrost, offering crucial data support for the formulation of strategies for permafrost thaw caused by climate change and anthropogenic activities.

Long-read powered viral metagenomics in the oligotrophic Sargasso Sea

Dominant microorganisms of the Sargasso Sea are key drivers of the global carbon cycle. However, associated viruses that shape microbial community structure and function are not well characterised. Here, we combined short and long read sequencing to survey Sargasso Sea phage communities in virus- and cellular fractions at viral maximum (80 m) and mesopelagic (200 m) depths. We identified 2,301 Sargasso Sea phage populations from 186 genera. Over half of the phage populations identified here lacked representation in global ocean viral metagenomes, whilst 177 of the 186 identified genera lacked representation in genomic databases of phage isolates. Viral fraction and cell-associated viral communities were decoupled, indicating viral turnover occurred across periods longer than the sampling period of three days. Inclusion of long-read data was critical for capturing the breadth of viral diversity. Phage isolates that infect the dominant bacterial taxa Prochlorococcus and Pelagibacter, usually regarded as cosmopolitan and abundant, were poorly represented.

Distribution patterns and functional diversity of DNA viruses determined by ecological niches in huge river ecosystems

While the vast number of DNA and RNA viruses participate in biogeochemical cycles in natural systems, little is known about virome in river ecosystems. Here, we analyzed the DNA viral composition and its metabolic potential in the Yangtze River, including freshwater (FW) and freshwater sediments (FWS). A total of 1237 river-derived virus contigs (RVCs) were obtained following de novo assembly from 62 metagenomics. We found that the viral diversity is significantly positively correlated longitudinally. Moreover, FW exhibited a greater viral variety and significantly different composition than FWS. The viral co-occurrence network suggested that positive correlations predominate between RVCs. Lastly, 1657 viral functions were predicted by gene ontology. Notably, 96 of 150 RVCs with higher weights identified by random-forest classier were more abundant in FW, which most engage organic cyclic compound metabolic processes and hydrolase activity. Together, this study highlights the previously unrecognized viruses and the importance of their distributions and functions in major river systems.

Environmental viromes reveal the global distribution signatures of deep-sea DNA viruses

INTRODUCTION

Viruses are abundant and ecologically significant in marine ecosystems. However, the virome of deep-sea sediments is not extensively investigated.

OBJECTIVES

To explore the distribution pattern of deep-sea viruses on a global scale, the viromes of DNA viruses isolated from 138 sediments of 5 deep-sea ecosystems were characterized.

METHODS

The viral particles were purified from each sediment sample. Then the viral DNAs were extracted and subjected to viral metagenomic analysis.

RESULTS

Here, we constructed a global deep-sea environmental virome dataset by analyzing the viral DNA of 138 sediment samples. A total of 347,737 viral operational taxonomic units (vOTUs) were identified, of which 84.94% were hitherto unknown, indicating that deep sea was a reservoir of novel DNA viruses. Furthermore, circular viral genome analysis revealed 98,581 complete genomes. The classified vOTUs included eukaryotic (44.55%) and prokaryotic (25.75%) viruses, and were taxonomically assigned to 63 viral families. The composition and abundance of the deep-sea sediment viromes were dependent on the deep-sea ecosystem as opposed to geographical region. Further analysis revealed that the viral community differentiation in different deep-sea ecosystems was driven by the virus-mediated energy metabolism.

CONCLUSION

Our findings showed that deep-sea ecosystems are a reservoir of novel DNA viruses and the viral community is shaped by the environmental characteristics of deep-sea ecosystems, thus presenting critical information for determining the ecological significance of viruses in global deep-sea ecosystems.

Aeromicrobiology: A global review of the cycling and relationships of bioaerosols with the atmosphere

Airborne microorganisms and biological matter (bioaerosols) play a key role in global biogeochemical cycling, human and crop health trends, and climate patterns. Their presence in the atmosphere is controlled by three main stages: emission, transport, and deposition. Aerial survival rates of bioaerosols are increased through adaptations such as ultra-violet radiation and desiccation resistance or association with particulate matter. Current research into modern concerns such as climate change, global gene transfer, and pathogenicity often neglects to consider atmospheric involvement. This comprehensive review outlines the transpiring of bioaerosols across taxa in the atmosphere, with significant focus on their interactions with environmental elements including abiotic factors (e.g., atmospheric composition, water cycle, and pollution) and events (e.g., dust storms, hurricanes, and wildfires). The aim of this review is to increase understanding and shed light on needed research regarding the interplay between global atmospheric phenomena and the aeromicrobiome. The abundantly documented bacteria and fungi are discussed in context of their cycling and human health impacts. Gaps in knowledge regarding airborne viral community, the challenges and importance of studying their composition, concentrations and survival in the air are addressed, along with understudied plant pathogenic oomycetes, and archaea cycling. Key methodologies in sampling, collection, and processing are described to provide an up-to-date picture of ameliorations in the field. We propose optimization to microbiological methods, commonly used in soil and water analysis, that adjust them to the context of aerobiology, along with other directions towards novel and necessary advancements. This review offers new perspectives into aeromicrobiology and calls for advancements in global-scale bioremediation, insights into ecology, climate change impacts, and pathogenicity transmittance.

Breaking the Ice: A Review of Phages in Polar Ecosystems

Bacteriophages, or phages, are viruses that infect and replicate within bacterial hosts, playing a significant role in regulating microbial populations and ecosystem dynamics. However, phages from extreme environments such as polar regions remain relatively understudied due to challenges such as restricted ecosystem access and low biomass. Understanding the diversity, structure, and functions of polar phages is crucial for advancing our knowledge of the microbial ecology and biogeochemistry of these environments. In this review, we will explore the current state of knowledge on phages from the Arctic and Antarctic, focusing on insights gained from -omic studies, phage isolation, and virus-like particle abundance data. Metagenomic studies of polar environments have revealed a high diversity of phages with unique genetic characteristics, providing insights into their evolutionary and ecological roles. Phage isolation studies have identified novel phage-host interactions and contributed to the discovery of new phage species. Virus-like particle abundance and lysis rate data, on the other hand, have highlighted the importance of phages in regulating bacterial populations and nutrient cycling in polar environments. Overall, this review aims to provide a comprehensive overview of the current state of knowledge about polar phages, and by synthesizing these different sources of information, we can better understand the diversity, dynamics, and functions of polar phages in the context of ongoing climate change, which will help to predict how polar ecosystems and residing phages may respond to future environmental perturbations.

Phylogenomic analysis expands the known repertoire of single-stranded DNA viruses in benthic zones of the South Indian Ocean

Single-stranded (ss) DNA viruses are ubiquitous and constitute some of the most diverse entities on Earth. Most studies have focused on ssDNA viruses from terrestrial environments resulting in a significant deficit in benthic ecosystems including aphotic zones of the South Indian Ocean (SIO). Here, we assess the diversity and phylogeny of ssDNA in deep waters of the SIO using a combination of established viral taxonomy tools and a Hidden Markov Model based approach. Replication initiator protein-associated (Rep) phylogenetic reconstruction and sequence similarity networks were used to show that the SIO hosts divergent and as yet unknown circular Rep-encoding ssDNA viruses. Several sequences appear to represent entirely novel families, expanding the repertoire of known ssDNA viruses. Results suggest that a small proportion of these viruses may be circular genetic elements, which may strongly influence the diversity of both eukaryotes and prokaryotes in the SIO. Taken together, our data show that the SIO harbours a diverse assortment of previously unknown ssDNA viruses. Due to their potential to infect a variety of hosts, these viruses may be crucial for marine nutrient recycling through their influence of the biological carbon pump.

Choice of Ultrafilter Affects Recovery Rate of Bacteriophages

Studies into the viral fraction of complex microbial communities, like in the mammalian gut, have recently garnered much interest. Yet there is still no standardized protocol for extracting viruses from such samples, and the protocols that exist employ procedures that skew the viral community of the sample one way or another. The first step of the extraction pipeline often consists of the basic filtering of macromolecules and bacteria, yet even this affects the viruses in a strain-specific manner. In this study, we investigate a protocol for viral extraction based on ultrafiltration and how the choice of ultrafilter might influence the extracted viral community. Clinical samples (feces, vaginal swabs, and tracheal suction samples) were spiked with a mock community of known phages (T4, c2, Φ6, Φ29, Φx174, and Φ2972), filtered, and quantified using spot and plaque assays to estimate the loss in recovery. The enveloped Φ6 phage is especially severely affected by the choice of filter, but also tailed phages such as T4 and c2 have a reduced infectivity after ultrafiltration. We conclude that the pore size of ultrafilters may affect the recovery of phages in a strain- and sample-dependent manner, suggesting the need for greater thought when selecting filters for virus extraction.

Microviruses: A World Beyond phiX174

Two decades of metagenomic analyses have revealed that in many environments, small (∼5 kb), single-stranded DNA phages of the family Microviridae dominate the virome. Although the emblematic microvirus phiX174 is ubiquitous in the laboratory, most other microviruses, particularly those of the gokushovirus and amoyvirus lineages, have proven to be much more elusive. This puzzling lack of representative isolates has hindered insights into microviral biology. Furthermore, the idiosyncratic size and nature of their genomes have resulted in considerable misjudgments of their actual abundance in nature. Fortunately, recent successes in microvirus isolation and improved metagenomic methodologies can now provide us with more accurate appraisals of their abundance, their hosts, and their interactions. The emerging picture is that phiX174 and its relatives are rather rare and atypical microviruses, and that a tremendous diversity of other microviruses is ready for exploration.

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.

Metavirome Profiling and Dynamics of the DNA Viral Community in Seawater in Chuuk State, Federated States of Micronesia

Despite their abundance and ecological importance, little is known about the diversity of marine viruses, in part because most cannot be cultured in the laboratory. Here, we used high-throughput viral metagenomics of uncultivated viruses to investigate the dynamics of DNA viruses in tropical seawater sampled from Chuuk State, Federated States of Micronesia, in March, June, and December 2014. Among the identified viruses, 71-79% were bacteriophages belonging to the families Myoviridae, Siphoviridae, and Podoviridae (Caudoviriales), listed in order of abundance at all sampling times. Although the measured environmental factors (temperature, salinity, and pH) remained unchanged in the seawater over time, viral dynamics changed. The proportion of cyanophages (34.7%) was highest in June, whereas the proportion of mimiviruses, phycodnaviruses, and other nucleo-cytoplasmic large DNA viruses (NCLDVs) was higher in March and December. Although host species were not analysed, the dramatic viral community change observed in June was likely due to changes in the abundance of cyanophage-infected cyanobacteria, whereas that in NCLDVs was likely due to the abundance of potential eukaryote-infected hosts. These results serve as a basis for comparative analyses of other marine viral communities, and guide policy-making when considering marine life care in Chuuk State.

DNA virome of ticks in the Northeast and Hubei provinces of China reveals diverse single-stranded circular DNA viruses

BACKGROUND

Ticks are medically important vectors capable of transmitting a variety of pathogens to and between host species. Although the spectrum of tick-borne RNA viruses has been frequently investigated, the diversity of tick-borne DNA viruses remains largely unknown.

METHODS

A total of 1571 ticks were collected from forests and infested animals, and the diversity of the viruses they harbored was profiled using a DNA-specific virome method. The viromic data were phylogenetically analyzed and validated by PCR assays.

RESULTS

Although diverse and abundant prokaryotic viruses were identified in the collected ticks, only eukaryotic DNA viruses with single-stranded circular genomes covering the anelloviruses and circular replication-associated (Rep) protein-encoding single-stranded (CRESS) DNA viruses were recovered from ticks. Anelloviruses were detected only in two tick pools, but CRESS DNA viruses were prevalent across these ticks except in one pool of Dermacentor spp. ticks. Phylogenetic analyses revealed that these tick-borne CRESS DNA viruses were related to viruses recovered from animal feces, tissues and even environmental samples, suggesting that their presence may be largely explained by environmental factors rather than by tick species and host blood meals.

CONCLUSIONS

Based on the results, tick-borne eukaryotic DNA viruses appear to be much less common than eukaryotic RNA viruses. Investigations involving a wider collection area and more diverse tick species are required to further support this speculation.

Covariance of Marine Nucleocytoplasmic Large DNA Viruses with Eukaryotic Plankton Communities in the Sub-Arctic Kongsfjorden Ecosystem: A Metagenomic Analysis of Marine Microbial Ecosystems

Nucleocytoplasmic large DNA viruses (NCLDVs) infect various marine eukaryotes. However, little is known about NCLDV diversity and their relationships with eukaryotic hosts in marine environments, the elucidation of which will advance the current understanding of marine ecosystems. This study characterizes the interplay between NCLDVs and the eukaryotic plankton community (EPC) in the sub-Arctic area using metagenomics and metabarcoding to investigate NCLDVs and EPC, respectively, in the Kongsfjorden ecosystem of Svalbard (Norway) in April and June 2018. Gyrodinium helveticum (Dinophyceae) is the most prevalent eukaryotic taxon in the EPC in April, during which time Mimiviridae (31.8%), Poxviridae (25.1%), Phycodnaviridae (14.7%) and Pandoraviridae (13.1%) predominate. However, in June, the predominant taxon is Aureococcus anophagefferens (Pelagophyceae), and the NCLDVs, Poxviridae (32.9%), Mimiviridae (29.1%), and Phycodnaviridae (18.5%) appear in higher proportions with an increase in Pelagophyceae, Bacillariophyceae, and Chlorophyta groups. Thus, differences in NCLDVs may be caused by changes in EPC composition in response to environmental changes, such as increases in water temperature and light intensity. Taken together, these findings are particularly relevant considering the anticipated impact of NCLDV-induced EPC control mechanisms on polar regions and, therefore, improve the understanding of the Sub-Arctic Kongsfjorden ecosystem.

Widespread yet persistent low abundance of TIM5-like cyanophages in the oceans

Ocean ecosystems are inhabited by a diverse set of viruses that impact microbial mortality and evolution. However, the distribution and abundances of specific viral lineages, particularly those from the large bank of rare viruses, remains largely unknown. Here, we assessed the diversity and abundance of the TIM5-like cyanophages. The sequencing of three new TIM5-like cyanophage genomes and environmental amplicons of a signature gene from the Red Sea revealed highly conserved gene content and sequence similarity. We adapted the polony method, a solid-phase polymerase chain reaction assay, to quantify TIM5-like cyanophages during three 2000 km expeditions in the Pacific Ocean and four annual cycles in the Red Sea. TIM5-like cyanophages were widespread, detected at all latitudes and seasons surveyed throughout the photic zone. Yet they were generally rare, ranging between <100 and 4000 viruses·ml^-1 . Occasional peaks in abundance of 10- to 100-fold were observed, reaching 71,000 viruses·ml^-1 . These peaks were ephemeral and seasonally variable in the Red Sea. Infection levels, quantified during one such peak, were very low. These characteristics of low diversity and abundance, as well as variable outbreaks, distinguishes the TIM5-like lineage from other major cyanophage lineages and illuminates that rare virus lineages can be persistent and widespread in the oceans.

Ecogenomics reveals viral communities across the Challenger Deep oceanic trench

Despite the environmental challenges and nutrient scarcity, the geographically isolated Challenger Deep in Mariana trench, is considered a dynamic hotspot of microbial activity. Hadal viruses are the least explored microorganisms in Challenger Deep, while their taxonomic and functional diversity and ecological impact on deep-sea biogeochemistry are poorly described. Here, we collect 13 sediment cores from slope and bottom-axis sites across the Challenger Deep (down to ~11 kilometers depth), and identify 1,628 previously undescribed viral operational taxonomic units at species level. Community-wide analyses reveals 1,299 viral genera and distinct viral diversity across the trench, which is significantly higher at the bottom-axis vs. slope sites of the trench. 77% of these viral genera have not been previously identified in soils, deep-sea sediments and other oceanic settings. Key prokaryotes involved in hadal carbon and nitrogen cycling are predicted to be potential hosts infected by these viruses. The detected putative auxiliary metabolic genes suggest that viruses at Challenger Deep could modulate the carbohydrate and sulfur metabolisms of their potential hosts, and stabilize host's cell membranes under extreme hydrostatic pressures. Our results shed light on hadal viral metabolic capabilities, contribute to understanding deep sea ecology and on functional adaptions of hadal viruses for future research.

Viruses in Subsurface Environments

Over the past 20 years, our knowledge of virus diversity and abundance in subsurface environments has expanded dramatically through application of quantitative metagenomic approaches. In most subsurface environments, viral diversity and abundance rival viral diversity and abundance observed in surface environments. Most of these viruses are uncharacterized in terms of their hosts and replication cycles. Analysis of accessory metabolic genes encoded by subsurface viruses indicates that they evolved to replicate within the unique features of their environments. The key question remains: What role do these viruses play in the ecology and evolution of the environments in which they replicate? Undoubtedly, as more virologists examine the role of viruses in subsurface environments, new insights will emerge.

A metagenomic insight into the Yangtze finless porpoise virome

The Yangtze finless porpoise (Neophocaena phocaenoides asiaeorientalis) inhabiting the Yantze River, China is critically endangered because of the influences of infectious disease, human activity, and water contamination. Viral diseases are one of the crucial factors that threatening the health of Yangtze finless porpoise. However, there are few studies which elaborate the viral diversity of Yangtze finless. Therefore, this study was performed to investigate the viral diversity of Yangtze finless by metagenomics. Results indicated that a total of 12,686,252 high-quality valid sequences were acquired and 2,172 virus reads were recognized. Additionally, we also obtained a total of 10,600 contigs. Phages was the most abundant virus in the samples and the ratio of DNA and RNA viruses were 69.75 and 30.25%, respectively. Arenaviridae, Ackermannviridae and Siphoviridae were the three most predominant families in all the samples. Moreover, the majority of viral genus were Mammarenavirus, Limestonevirus and Lambdavirus. The results of gene prediction indicated that these viruses play vital roles in biological process, cellular component, molecular function, and disease. To the best of our knowledge, this is the first report on the viral diversity of Yangtze finless porpoise, which filled the gaps in its viral information. Meanwhile, this study can also provide a theoretical basis for the establishment of the prevention and protection system for virus disease of Yangtze finless porpoise.

Comprehensive Evaluation of RNA and DNA Viromic Methods Based on Species Richness and Abundance Analyses Using Marmot Rectal Samples

Viral metagenomics is the most powerful tool to profile viromic composition for a given sample. Different viromic methods, including amplification-free ones, have been developed, but choosing them for different purposes requires comprehensive benchmarks. Here, we assessed the performance of four routinely used methods, i.e., multiple displacement amplification (MDA), direct metagenomic sequencing (MTG), sequence-independent single-primer amplification (SIA), and metatranscriptomic sequencing (MTT), using marmot rectal samples as the templates spiked with five known viruses of different genome types. The obtained clean data were differently contaminated by host and bacterial genomes, resulting in MDA having the most, with ~72.1%, but MTT had only ~7.5% data, useful for follow-up viromic analysis. MDA showed a broader spectrum with higher efficiency to profile the DNA virome, and MTT captured almost all RNA viruses with extraordinary sensitivity; hence, they are advisable in richness-based viromic studies. MTG was weak in capturing single-stranded DNA viruses, and SIA could detect both RNA and DNA viruses but with high randomness. Due to biases to certain types of viruses, the four methods caused different alterations to species abundance compared to the initial virus composition. SIA and MDA introduced greater stochastic errors to relative abundances of species, genus, and family taxa, whereas the two amplification-free methods were more tolerant toward such errors and thus are recommendable in abundance-based analyses. In addition, genus taxon is a compromising analytic level that ensures technically supported and biologically and/or ecologically meaningful viromic conclusions. IMPORTANCE Viral metagenomics can be roughly divided into species richness-based studies and species abundance-based analyses. Viromic methods with different principles have been developed, but rational selection of these techniques according to different purposes requires comprehensive understanding of their properties. By assessing the four most widely used methods using template samples, we found that multiple displacement amplification (MDA) and metatranscriptomic sequencing (MTT) are advisable for species richness-based viromic studies, as they show excellent efficiency to detect DNA and RNA viruses. Meanwhile, metagenomic sequencing (MTG) and MTT are more compatible with stochastic errors of methods introduced into relative abundance of viromic taxa and hence are rational choices in species abundance-based analyses. This study also highlights that MTG needs to tackle host genome contamination and ameliorate the capacity to detect single-stranded DNA viruses in the future, and the MTT method requires an improvement in bacterial rRNA depletion prior to library preparation.

Virioplankton assemblages from challenger deep, the deepest place in the oceans

Hadal ocean biosphere, that is, the deepest part of the world's oceans, harbors a unique microbial community, suggesting a potential uncovered co-occurring virioplankton assemblage. Herein, we reveal the unique virioplankton assemblages of the Challenger Deep, comprising 95,813 non-redundant viral contigs from the surface to the hadal zone. Almost all of the dominant viral contigs in the hadal zone were unclassified, potentially related to Alteromonadales and Oceanospirillales. 2,586 viral auxiliary metabolic genes from 132 different KEGG orthologous groups were mainly related to the carbon, nitrogen, sulfur, and arsenic metabolism. Lysogenic viral production and integrase genes were augmented in the hadal zone, suggesting the prevalence of viral lysogenic life strategy. Abundant rve genes in the hadal zone, which function as transposase in the caudoviruses, further suggest the prevalence of viral-mediated horizontal gene transfer. This study provides fundamental insights into the virioplankton assemblages of the hadal zone, reinforcing the necessity of incorporating virioplankton into the hadal biogeochemical cycles.

Comparative genomic analysis of five freshwater cyanophages and reference-guided metagenomic data mining

BACKGROUND

As important producers using photosynthesis on Earth, cyanobacteria contribute to the oxygenation of atmosphere and the primary production of biosphere. However, due to the eutrophication of urban waterbodies and global warming, uncontrollable growth of cyanobacteria usually leads to the seasonal outbreak of cyanobacterial blooms. Cyanophages, a group of viruses that specifically infect and lyse cyanobacteria, are considered as potential environment-friendly agents to control the harmful blooms. Compared to the marine counterparts, only a few freshwater cyanophages have been isolated and genome sequenced to date, largely limiting their characterizations and applications.

RESULTS

Here, we isolated five freshwater cyanophages varying in tail morphology, termed Pam1~Pam5, all of which infect the cyanobacterium Pseudanabaena mucicola Chao 1806 that was isolated from the bloom-suffering Lake Chaohu in Anhui, China. The whole-genome sequencing showed that cyanophages Pam1~Pam5 all contain a dsDNA genome, varying in size from 36 to 142 Kb. Phylogenetic analyses suggested that Pam1~Pam5 possess different DNA packaging mechanisms and are evolutionarily distinct from each other. Notably, Pam1 and Pam5 have lysogeny-associated gene clusters, whereas Pam2 possesses 9 punctuated DNA segments identical to the CRISPR spacers in the host genome. Metagenomic data-based calculation of the relative abundance of Pam1~Pam5 at the Nanfei estuary towards the Lake Chaohu revealed that the short-tailed Pam1 and Pam5 account for the majority of the five cyanophages. Moreover, comparative analyses of the reference genomes of Pam1~Pam5 and previously reported cyanophages enabled us to identify three circular and seven linear contigs of virtual freshwater cyanophages from the metagenomic data of the Lake Chaohu.

CONCLUSIONS

We propose a high-throughput strategy to systematically identify cyanophages based on the currently available metagenomic data and the very limited reference genomes of experimentally isolated cyanophages. This strategy could be applied to mine the complete or partial genomes of unculturable bacteriophages and viruses. Transformation of the synthesized whole genomes of these virtual phages/viruses to proper hosts will enable the rescue of bona fide viral particles and eventually enrich the library of microorganisms that exist on Earth. Video abstract.

New Microviridae isolated from Sulfitobacter reveals two cosmopolitan subfamilies of single-stranded DNA phages infecting marine and terrestrial Alphaproteobacteria

The Microviridae family represents one of the major clades of single-stranded DNA (ssDNA) phages. Their cultivated members are lytic and infect Proteobacteria, Bacteroidetes, and Chlamydiae. Prophages have been predicted in the genomes from Bacteroidales, Hyphomicrobiales, and Enterobacteriaceae and cluster within the 'Alpavirinae', 'Amoyvirinae', and Gokushovirinae. We have isolated 'Ascunsovirus oldenburgi' ICBM5, a novel phage distantly related to known Microviridae. It infects Sulfitobacter dubius SH24-1b and uses both a lytic and a carrier-state life strategy. Using ICBM5 proteins as a query, we uncovered in publicly available resources sixty-five new Microviridae prophages and episomes in bacterial genomes and retrieved forty-seven environmental viral genomes (EVGs) from various viromes. Genome clustering based on protein content and phylogenetic analysis showed that ICBM5, together with Rhizobium phages, new prophages, episomes, and EVGs cluster within two new phylogenetic clades, here tentatively assigned the rank of subfamily and named 'Tainavirinae' and 'Occultatumvirinae'. They both infect Rhodobacterales. Occultatumviruses also infect Hyphomicrobiales, including nitrogen-fixing endosymbionts from cosmopolitan legumes. A biogeographical assessment showed that tainaviruses and occultatumviruses are spread worldwide, in terrestrial and marine environments. The new phage isolated here sheds light onto new and diverse branches of the Microviridae tree, suggesting that much of the ssDNA phage diversity remains in the dark.

Applying the concept of liquid biopsy to monitor the microbial biodiversity of marine coastal ecosystems

Liquid biopsy (LB) is a concept that is rapidly gaining ground in the biomedical field. Its concept is largely based on the detection of circulating cell-free DNA (ccfDNA) fragments that are mostly released as small fragments following cell death in various tissues. A small percentage of these fragments are from foreign (nonself) tissues or organisms. In the present work, we applied this concept to mussels, a sentinel species known for its high filtration capacity of seawater. We exploited the capacity of mussels to be used as natural filters to capture environmental DNA fragments of different origins to provide information on the biodiversity of marine coastal ecosystems. Our results showed that hemolymph of mussels contains DNA fragments that varied considerably in size, ranging from 1 to 5 kb. Shotgun sequencing revealed that a significant amount of DNA fragments had a nonself microbial origin. Among these, we found DNA fragments derived from bacteria, archaea, and viruses, including viruses known to infect a variety of hosts that commonly populate coastal marine ecosystems. Taken together, our study shows that the concept of LB applied to mussels provides a rich and yet unexplored source of knowledge regarding the microbial biodiversity of a marine coastal ecosystem.

Viral Eco-Genomic Tools: Development and Implementation for Aquatic Biomonitoring

Enteric viruses (EVs) occurrence within aquatic environments varies and leads to significant risk on public health of humans, animals, and diversity of aquatic taxa. Early and efficacious recognition of cultivable and fastidious EVs in aquatic systems are important to ensure the sanitary level of aquatic water and implement required treatment strategies. Herein, we provided a comprehensive overview of the conventional and up-to-date eco-genomic tools for aquatic biomonitoring of EVs, aiming to develop better water pollution monitoring tools. In combination with bioinformatics techniques, genetic tools including cloning sequencing analysis, DNA microarray, next-generation sequencing (NGS), and metagenomic sequencing technologies are implemented to make informed decisions about the global burden of waterborne EVs-associated diseases. The data presented in this review are helpful to recommend that: (1) Each viral pollution detection method has its own merits and demerits; therefore, it would be advantageous for viral pollution evaluation to be integrated as a complementary platform. (2) The total viral genome pool extracted from aquatic environmental samples is a real reflection of pollution status of the aquatic eco-systems; therefore, it is recommended to conduct regular sampling through the year to establish an updated monitoring system for EVs, and quantify viral peak concentrations, viral typing, and genotyping. (3) Despite that conventional detection methods are cheaper, it is highly recommended to implement molecular-based technologies to complement aquatic ecosystems biomonitoring due to numerous advantages including high-throughput capability. (4) Continuous implementation of the eco-genetic detection tools for monitoring the EVs in aquatic ecosystems is recommended.

The Case for Studying New Viruses of New Hosts

Virology has largely focused on viruses that are pathogenic to humans or to the other species that we care most about. There is no doubt that this has been a worthwhile investment. But many transformative advances have been made through the in-depth study of relatively obscure viruses that do not appear on lists of prioritized pathogens. In this review, I highlight the benefits that can accrue from the study of viruses and hosts off the beaten track. I take stock of viral sequence diversity across host taxa as an estimate of the bias that exists in our understanding of host-virus interactions. I describe the gains that have been made through the metagenomic discovery of thousands of new viruses in previously unsampled hosts as well as the limitations of metagenomic surveys. I conclude by suggesting that the study of viruses that naturally infect existing and emerging model organisms represents an opportunity to push virology forward in useful and hard to predict ways.Expected final online publication date for the Annual Review of Virology, Volume 9 is September 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

A Novel and Ubiquitous Marine Methylophage Provides Insights into Viral-Host Coevolution and Possible Host-Range Expansion in Streamlined Marine Heterotrophic Bacteria

The methylotrophic OM43 clade are Gammaproteobacteria that comprise some of the smallest free-living cells known and have highly streamlined genomes. OM43 represents an important microbial link between marine primary production and remineralization of carbon back to the atmosphere. Bacteriophages shape microbial communities and are major drivers of mortality and global marine biogeochemistry. Recent cultivation efforts have brought the first viruses infecting members of the OM43 clade into culture. Here, we characterize a novel myophage infecting OM43 called Melnitz. Melnitz was isolated independently from water samples from a subtropical ocean gyre (Sargasso Sea) and temperate coastal (Western English Channel) systems. Metagenomic recruitment from global ocean viromes confirmed that Melnitz is globally ubiquitous, congruent with patterns of host abundance. Bacteria with streamlined genomes such as OM43 and the globally dominant SAR11 clade use riboswitches as an efficient method to regulate metabolism. Melnitz encodes a two-piece tmRNA (ssrA), controlled by a glutamine riboswitch, providing evidence that riboswitch use also occurs for regulation during phage infection of streamlined heterotrophs. Virally encoded tRNAs and ssrA found in Melnitz were phylogenetically more closely related to those found within the alphaproteobacterial SAR11 clade and their associated myophages than those within their gammaproteobacterial hosts. This suggests the possibility of an ancestral host transition event between SAR11 and OM43. Melnitz and a related myophage that infects SAR11 were unable to infect hosts of the SAR11 and OM43, respectively, suggesting host transition rather than a broadening of host range. IMPORTANCE Isolation and cultivation of viruses are the foundations on which the mechanistic understanding of virus-host interactions and parameterization of bioinformatic tools for viral ecology are based. This study isolated and characterized the first myophage known to infect the OM43 clade, expanding our knowledge of this understudied group of microbes. The nearly identical genomes of four strains of Melnitz isolated from different marine provinces and the global abundance estimations from metagenomic data suggest that this viral population is globally ubiquitous. Genome analysis revealed several unusual features in Melnitz and related genomes recovered from viromes, such as a curli operon and virally encoded tmRNA controlled by a glutamine riboswitch, neither of which are found in the host. Further phylogenetic analysis of shared genes indicates that this group of viruses infecting the gammaproteobacterial OM43 shares a recent common ancestor with viruses infecting the abundant alphaproteobacterial SAR11 clade. Host ranges are affected by compatible cell surface receptors, successful circumvention of superinfection exclusion systems, and the presence of required accessory proteins, which typically limits phages to singular narrow groups of closely related bacterial hosts. This study provides intriguing evidence that for streamlined heterotrophic bacteria, virus-host transitioning may not be necessarily restricted to phylogenetically related hosts but is a function of shared physical and biochemical properties of the cell.

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.

Unravelling the ecological complexity of soil viromes: Challenges and opportunities

Viruses are extremely abundant and ubiquitous in soil, and significantly contribute to various terrestrial ecosystem processes such as biogeochemical nutrient cycling, microbiome regulation and community assembly, and host evolutionary dynamics. Despite their numerous dominance and functional importance, understanding soil viral ecology is a formidable challenge, because of the technological challenges to characterize the abundance, diversity and community compositions of viruses, and their interactions with other organisms in the complex soil environment. Viruses may engage in a myriad of biological interactions within soil food webs across a broad range of spatiotemporal scales and are exposed to various biotic and abiotic disturbances. Current studies on the soil viromes, however, often describe the complexity of their tremendous diversity, but lack of exploring their potential ecological roles. In this article, we summarized the major methods to decipher the ecology of soil viruses, discussed biotic and abiotic factors and global change factors that shape the diversity and composition of soil viromes, and the ecological roles of soil viruses. We also proposed a new framework to understand the ecological complexity of viruses from micro to macro ecosystem scales and to predict and unravel their activities in terrestrial ecosystems.

Oral Toxicity Study for Salmonella Killing Lytic Bacteriophage NINP13076 in BALB/c Mice and Its Effect on Probiotic Microbiota

Viruses that infect bacteria are emerging as attractive biocontrol agents and biopreservatives for foods. Since these bacteriophages kill the target pathogens by lysis and are also consumed along with food, it is essential to evaluate their collateral toxicity on the probiotic gut microbiota. In this study, we examined the acute oral toxicity of a Salmonella phage isolated from sewage in mice. Acute oral administration of the Salmonella phage for five consecutive days did not show any significant pathological changes in the vital organs like lung, kidneys, heart, liver, and intestine. In addition, growth of typical probiotic microbiota remained unaffected even after incubation up to 24 h with the Salmonella phage. The results of this study clearly showed that oral administration of the lytic Salmonella phage did not have any significant adverse effects on the animals, may not harm the probiotic gut microbiota, and are likely to be safe for use in food preservation.

Transcriptional Organization of the Salmonella Typhimurium Phage P22 pid ORFan Locus

Many phage genes lack sequence similarity to any other open reading frame (ORF) in current databases. These enigmatic ORFan genes can have a tremendous impact on phage propagation and host interactions but often remain experimentally unexplored. We previously revealed a novel interaction between phage P22 and its Salmonella Typhimurium host, instigated by the ORFan gene pid (for phage P22 encoded instigator of dgo expression) and resulting in derepression of the host dgoRKAT operon. The pid gene is highly expressed in phage carrier cells that harbor a polarly located P22 episome that segregates asymmetrically among daughter cells. Here, we discovered that the pid locus is fitted with a weak promoter, has an exceptionally long 5′ untranslated region that is instructive for a secondary pid mRNA species, and has a 3′ Rho-independent termination loop that is responsible for stability of the pid transcript.

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.

Viruses infecting a warm water picoeukaryote shed light on spatial co-occurrence dynamics of marine viruses and their hosts

The marine picoeukaryote Bathycoccus prasinos has been considered a cosmopolitan alga, although recent studies indicate two ecotypes exist, Clade BI (B. prasinos) and Clade BII. Viruses that infect Bathycoccus Clade BI are known (BpVs), but not that infect BII. We isolated three dsDNA prasinoviruses from the Sargasso Sea against Clade BII isolate RCC716. The BII-Vs do not infect BI, and two (BII-V2 and BII-V3) have larger genomes (~210 kb) than BI-Viruses and BII-V1. BII-Vs share ~90% of their proteins, and between 65% to 83% of their proteins with sequenced BpVs. Phylogenomic reconstructions and PolB analyses establish close-relatedness of BII-V2 and BII-V3, yet BII-V2 has 10-fold higher infectivity and induces greater mortality on host isolate RCC716. BII-V1 is more distant, has a shorter latent period, and infects both available BII isolates, RCC716 and RCC715, while BII-V2 and BII-V3 do not exhibit productive infection of the latter in our experiments. Global metagenome analyses show Clade BI and BII algal relative abundances correlate positively with their respective viruses. The distributions delineate BI/BpVs as occupying lower temperature mesotrophic and coastal systems, whereas BII/BII-Vs occupy warmer temperature, higher salinity ecosystems. Accordingly, with molecular diagnostic support, we name Clade BII Bathycoccus calidus sp. nov. and propose that molecular diversity within this new species likely connects to the differentiated host-virus dynamics observed in our time course experiments. Overall, the tightly linked biogeography of Bathycoccus host and virus clades observed herein supports species-level host specificity, with strain-level variations in infection parameters.

Viral Characteristics of the Warm Atlantic and Cold Arctic Water Masses in the Nordic Seas

Nordic Seas are the subarctic seas connecting the Arctic Ocean and North Atlantic Ocean with complex water masses, experiencing an abrupt climate change. Though knowledge of the marine virosphere has expanded rapidly, the diversity of viruses and their relationships with host cells and water masses in the Nordic Seas remain to be fully revealed. Here, we establish the Nordic Sea DNA virome (NSV) data set of 55,315 viral contigs including 1,478 unique viral populations from seven stations influenced by both the warm Atlantic and cold Arctic water masses. Caudovirales dominated in the seven NSVs, especially in the warm Atlantic waters. The major giant nucleocytoplasmic large DNA viruses (NCLDVs) contributed a significant proportion of the classified viral contigs in the NSVs (32.2%), especially in the cold Arctic waters (44.9%). The distribution patterns of Caudovirales and NCLDVs were a reflection of the community structure of their hosts in the corresponding water masses and currents. Latitude, pH, and flow speed were found to be key factors influencing the microbial communities and coinfluencing the variation of viral communities. Network analysis illustrated the tight coupling between the variation of viral communities and microbial communities in the Nordic Seas. This study suggests a probable linkage between viromes, host cells, and surface water masses from both the cool Arctic and warm Atlantic Oceans. IMPORTANCE This is a systematic study of Nordic Sea viromes using metagenomic analysis. The viral diversity, community structure, and their relationships with host cells and the complex water masses from both the cool Arctic and the warm Atlantic oceans were illustrated. The NCLDVs and Caudovirales are proposed as the viral characteristics of the cold Arctic and warm Atlantic waters, respectively. This study provides an important background for the viromes in the subarctic seas connecting the Arctic Ocean and North Atlantic Ocean and sheds light on their responses to abrupt climate change in the future.

Viruses and Their Interactions With Bacteria and Archaea of Hypersaline Great Salt Lake

Viruses play vital biogeochemical and ecological roles by (a) expressing auxiliary metabolic genes during infection, (b) enhancing the lateral transfer of host genes, and (c) inducing host mortality. Even in harsh and extreme environments, viruses are major players in carbon and nutrient recycling from organic matter. However, there is much that we do not yet understand about viruses and the processes mediated by them in the extreme environments such as hypersaline habitats. The Great Salt Lake (GSL) in Utah, United States is a hypersaline ecosystem where the biogeochemical role of viruses is poorly understood. This study elucidates the diversity of viruses and describes virus–host interactions in GSL sediments along a salinity gradient. The GSL sediment virosphere consisted of Haloviruses (32.07 ± 19.33%) and members of families Siphoviridae (39.12 ± 19.8%), Myoviridae (13.7 ± 6.6%), and Podoviridae (5.43 ± 0.64%). Our results demonstrate that salinity alongside the concentration of organic carbon and inorganic nutrients (nitrogen and phosphorus) governs the viral, bacteria, and archaeal diversity in this habitat. Computational host predictions for the GSL viruses revealed a wide host range with a dominance of viruses that infect Proteobacteria, Actinobacteria, and Firmicutes. Identification of auxiliary metabolic genes for photosynthesis (psbA), carbon fixation (rbcL, cbbL), formaldehyde assimilation (SHMT), and nitric oxide reduction (NorQ) shed light on the roles played by GSL viruses in biogeochemical cycles of global relevance.

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.

CoCoNet: an efficient deep learning tool for viral metagenome binning

MOTIVATION

Metagenomic approaches hold the potential to characterize microbial communities and unravel the intricate link between the microbiome and biological processes. Assembly is one of the most critical steps in metagenomics experiments. It consists of transforming overlapping DNA sequencing reads into sufficiently accurate representations of the community's genomes. This process is computationally difficult and commonly results in genomes fragmented across many contigs. Computational binning methods are used to mitigate fragmentation by partitioning contigs based on their sequence composition, abundance or chromosome organization into bins representing the community's genomes. Existing binning methods have been principally tuned for bacterial genomes and do not perform favorably on viral metagenomes.

RESULTS

We propose Composition and Coverage Network (CoCoNet), a new binning method for viral metagenomes that leverages the flexibility and the effectiveness of deep learning to model the co-occurrence of contigs belonging to the same viral genome and provide a rigorous framework for binning viral contigs. Our results show that CoCoNet substantially outperforms existing binning methods on viral datasets.

AVAILABILITY AND IMPLEMENTATION

CoCoNet was implemented in Python and is available for download on PyPi (https://pypi.org/). The source code is hosted on GitHub at https://github.com/Puumanamana/CoCoNet and the documentation is available at https://coconet.readthedocs.io/en/latest/index.html. CoCoNet does not require extensive resources to run. For example, binning 100k contigs took about 4 h on 10 Intel CPU Cores (2.4 GHz), with a memory peak at 27 GB (see Supplementary Fig. S9). To process a large dataset, CoCoNet may need to be run on a high RAM capacity server. Such servers are typically available in high-performance or cloud computing settings.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Exploration of the Diversity of Clustered Regularly Interspaced Short Palindromic Repeats-Cas Systems in Clostridium novyi sensu lato

Classified as the genospecies Clostridium novyi sensu lato and distributed into four lineages (I-IV), Clostridium botulinum (group III), Clostridium novyi, and Clostridium haemolyticum are clostridial pathogens that cause animal diseases. Clostridium novyi sensu lato contains a large mobilome consisting of plasmids and circular bacteriophages. Here, we explored clustered regularly interspaced short palindromic repeats (CRISPR) arrays and their associated proteins (Cas) to shed light on the link between evolution of CRISPR-Cas systems and the plasmid and phage composition in a study of 58 Clostridium novyi sensu lato genomes. In 55 of these genomes, types I-B (complete or partial), I-D, II-C, III-B, III-D, or V-U CRISPR-Cas systems were detected in chromosomes as well as in mobile genetic elements (MGEs). Type I-B predominated (67.2%) and was the only CRISPR type detected in the Ia, III, and IV genomic lineages. Putative type V-U CRISPR Cas14a genes were detected in two different cases: next to partial type-IB CRISPR loci on the phage encoding the botulinum neurotoxin (BoNT) in lineage Ia and in 12 lineage II genomes, as part of a putative integrative element related to a phage-inducible chromosomal island (PICI). In the putative PICI, Cas14a was associated with CRISPR arrays and restriction modification (RM) systems as part of an accessory locus. This is the first time a PICI containing such locus has been detected in C. botulinum. Mobilome composition and dynamics were also investigated based on the contents of the CRISPR arrays and the study of spacers. A large proportion of identified protospacers (20.2%) originated from Clostridium novyi sensu lato (p1_Cst, p4_BKT015925, p6_Cst, CWou-2020a, p1_BKT015925, and p2_BKT015925), confirming active exchanges within this genospecies and the key importance of specific MGEs in Clostridium novyi sensu lato.

Hostile Takeover: How Viruses Reprogram Prokaryotic Metabolism

To reproduce, prokaryotic viruses must hijack the cellular machinery of their hosts and redirect it toward the production of viral particles. While takeover of the host replication and protein synthesis apparatus has long been considered an essential feature of infection, recent studies indicate that extensive reprogramming of host primary metabolism is a widespread phenomenon among prokaryotic viruses that is required to fulfill the biosynthetic needs of virion production. In this review we provide an overview of the most significant recent findings regarding virus-induced reprogramming of prokaryotic metabolism and suggest how quantitative systems biology approaches may be used to provide a holistic understanding of metabolic remodeling during lytic viral infection. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Illuminating the Virosphere Through Global Metagenomics

Viruses are the most abundant biological entity on Earth, infect cellular organisms from all domains of life, and are central players in the global biosphere. Over the last century, the discovery and characterization of viruses have progressed steadily alongside much of modern biology. In terms of outright numbers of novel viruses discovered, however, the last few years have been by far the most transformative for the field. Advances in methods for identifying viral sequences in genomic and metagenomic datasets, coupled to the exponential growth of environmental sequencing, have greatly expanded the catalog of known viruses and fueled the tremendous growth of viral sequence databases. Development and implementation of new standards, along with careful study of the newly discovered viruses, have transformed and will continue to transform our understanding of microbial evolution, ecology, and biogeochemical cycles, leading to new biotechnological innovations across many diverse fields, including environmental, agricultural, and biomedical sciences.

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.

Viral footprints across Gulfs of Kathiawar Peninsula and Arabian Sea: Unraveled from pelagic sediment metagenomic data

Marine biosphere is one of the largest, diverse and dynamic system hosting numerous of microorganisms. Viruses being the most abundant under explored lifeforms in ocean, represent a reservoir of great genetic diversity. We report the metagenomic insights on the viral communities in the deep sediments of the two Gulfs of Gujarat i.e. Gulf of Khambhat and Gulf of Kutch, with one sample from Arabian Sea, treated as open sea control. The viral reads were filtered from the whole dataset, assembled and studied for viral diversity, which was visualized by Pavian. The sequences were checked for the viral abundance, diversity and functionality. The resulting viral taxonomic classification contained 6 orders, 8 families and 47 genera. The results revealed that the phages infecting Cyanobacterium, Bacillus and Vibrio dominated the sediments. Further, it was observed that majority of viral sequences belonged to double-stranded DNA phages. The present study attempts to provide a primary insight of the viral signals and potential genetic content in the Gulfs of Kathiawar.

Going to extremes - a metagenomic journey into the dark matter of life

The Virus-X-Viral Metagenomics for Innovation Value-project was a scientific expedition to explore and exploit uncharted territory of genetic diversity in extreme natural environments such as geothermal hot springs and deep-sea ocean ecosystems. Specifically, the project was set to analyse and exploit viral metagenomes with the ultimate goal of developing new gene products with high innovation value for applications in biotechnology, pharmaceutical, medical, and the life science sectors. Viral gene pool analysis is also essential to obtain fundamental insight into ecosystem dynamics and to investigate how viruses influence the evolution of microbes and multicellular organisms. The Virus-X Consortium, established in 2016, included experts from eight European countries. The unique approach based on high throughput bioinformatics technologies combined with structural and functional studies resulted in the development of a biodiscovery pipeline of significant capacity and scale. The activities within the Virus-X consortium cover the entire range from bioprospecting and methods development in bioinformatics to protein production and characterisation, with the final goal of translating our results into new products for the bioeconomy. The significant impact the consortium made in all of these areas was possible due to the successful cooperation between expert teams that worked together to solve a complex scientific problem using state-of-the-art technologies as well as developing novel tools to explore the virosphere, widely considered as the last great frontier of life.

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.

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.

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.

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.

A library preparation optimized for metagenomics of RNA viruses

Our understanding of the viral communities associated to animals has not yet reached the level attained on the bacteriome. This situation is due to, among others, technical challenges in adapting metagenomics using high-throughput sequencing to the study of RNA viromes in animals. Although important developments have been achieved in most steps of viral metagenomics, there is yet a key step that has received little attention: the library preparation. This situation differs from bacteriome studies in which developments in library preparation have largely contributed to the democratisation of metagenomics. Here, we present a library preparation optimized for metagenomics of RNA viruses from insect vectors of viral diseases. The library design allows a simple PCR-based preparation, such as those routinely used in bacterial metabarcoding, that is adapted to shotgun sequencing as required in viral metagenomics. We first optimized our library preparation using mock viral communities and then validated a full metagenomic approach incorporating our preparation in two pilot studies with field-caught insect vectors; one including a comparison with a published metagenomic protocol. Our approach provided a fold increase in virus-like sequences compared to other studies, and nearly-full genomes from new virus species. Moreover, our results suggested conserved trends in virome composition within a population of a mosquito species. Finally, the sensitivity of our approach was compared to a commercial diagnostic PCR for the detection of an arbovirus in field-caught insect vectors. Our approach could facilitate studies on viral communities from animals and the democratization of metagenomics in community ecology of viruses.

The Problem with ‘Microbiome’

The term “microbiome” is currently applied predominantly to assemblages of organisms with 16S rRNA genes. In this context, “microbiome” is a misnomer that has been conferred a wide-ranging primacy over terms for community members lacking such genes, e.g., mycobiome, eukaryome, and virome, yet these are also important subsets of microbial communities. Widespread convenient and affordable 16S rRNA sequencing pipelines have accelerated continued use of such a “microbiome”, but at what intellectual and practical costs? Here we show that the use of “microbiome” in ribosomal gene-based studies has been egregiously misapplied, and discuss potential impacts. We argue that the current focus of “microbiome” research, predominantly on only ‘bacteria’, presents a dangerous narrowing of scope which encourages dismissal and even ignorance of other organisms’ contributions to microbial diversity, sensu stricto, and as etiologic agents; we put this in context by discussing cases in both marine microbial diversity and the role of pathogens in global amphibian decline. Fortunately, the solution is simple. We must use descriptive nouns that strictly reflect the outcomes attainable by the methods used. “Microbiome”, as a descriptive noun, should only be used when diversity in the three recognized domains is explored.

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 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.