Metagenomic analysis of coastal RNA virus communities

Novel RNA viruses from the Atlantic Ocean: Ecogenomics, biogeography, and total virioplankton mass contribution from surface to the deep ocean

Marine viruses play a major role in the energy and nutrient cycle and affect the evolution of their hosts. Despite their importance, there is still little knowledge about RNA viruses. Here, we have explored the Atlantic Ocean, from surface to deep (4.296 m), and used viromics and quantitative methods to unveil the genomics, biogeography, and the mass contribution of RNA viruses to the total viroplankton. A total of 2481 putative RNA viral contigs (>500 bp) and 107 larger bona fide RNA viral genomes (>2.5 kb) were identified; 88 of them representing novel viruses belonging mostly to two clades: Yangshan assemblage (sister clade to the class Alsuviricetes) and Nodaviridae. These viruses were highly endemic and locally abundant, with little or no presence in other oceans since only ≈15% of them were found in at least one of the Tara sampling metatranscriptomes. Quantitative data indicated that the abundance of RNA viruses in the surface and deep chlorophyll maximum zone was within ≈106  VLP/mL representing a potential contribution of 5.2%-24.4% to the total viroplankton community (DNA and RNA viruses), with DNA viruses being the predominant members (≈107  VLP/mL). However, for the deep sample, the observed trend was the opposite, although as further discussed, several biases should be considered. Together these results contribute to our understanding of the diversity, abundance, and distribution of RNA viruses in the oceans and provide a basis for further investigation into their ecological roles and biogeography.

Metagenomic surveillance uncovers diverse and novel viral taxa in febrile patients from Nigeria

Effective infectious disease surveillance in high-risk regions is critical for clinical care and pandemic preemption; however, few clinical diagnostics are available for the wide range of potential human pathogens. Here, we conduct unbiased metagenomic sequencing of 593 samples from febrile Nigerian patients collected in three settings: i) population-level surveillance of individuals presenting with symptoms consistent with Lassa Fever (LF); ii) real-time investigations of outbreaks with suspected infectious etiologies; and iii) undiagnosed clinically challenging cases. We identify 13 distinct viruses, including the second and third documented cases of human blood-associated dicistrovirus, and a highly divergent, unclassified dicistrovirus that we name human blood-associated dicistrovirus 2. We show that pegivirus C is a common co-infection in individuals with LF and is associated with lower Lassa viral loads and favorable outcomes. We help uncover the causes of three outbreaks as yellow fever virus, monkeypox virus, and a noninfectious cause, the latter ultimately determined to be pesticide poisoning. We demonstrate that a local, Nigerian-driven metagenomics response to complex public health scenarios generates accurate, real-time differential diagnoses, yielding insights that inform policy.

In Situ Molecular Ecological Analyses Illuminate Distinct Factors Regulating Formation and Demise of a Harmful Dinoflagellate Bloom

The development and demise of a harmful algal bloom (HAB) are generally regulated by multiple processes; identifying specific critical drivers for a specific bloom is important yet challenging. Here, we conducted a whole-assemblage molecular ecological study on a dinoflagellate bloom to address the hypothesis that energy and nutrient acquisition, defense against grazing and microbial attacks, and sexual reproduction are critical to the rise and demise of the bloom. Microscopic and molecular analyses identified the bloom-causing species as Karenia longicanalis and showed that the ciliate Strombidinopsis sp. was dominant in a nonbloom plankton community, whereas the diatom Chaetoceros sp. dominated the after-bloom community, along with remarkable shifts in the community structure for both eukaryotes and prokaryotes. Metatranscriptomic analysis indicated that heightened energy and nutrient acquisition in K. longicanalis significantly contributed to bloom development. In contrast, active grazing by the ciliate Strombidinopsis sp. and attacks by algicidal bacteria (Rhodobacteracea, Cryomorphaceae, and Rhodobacteracea) and viruses prevented (at nonbloom stage) or collapsed the bloom (in after-bloom stage). Additionally, nutrition competition by the Chaetoceros diatoms plausibly contributed to bloom demise. The findings suggest the importance of energy and nutrients in promoting this K. longicanalis bloom and the failure of antimicrobial defense and competition of diatoms as the major bloom suppressor and terminator. This study provides novel insights into bloom-regulating mechanisms and the first transcriptomic data set of K. longicanalis, which will be a valuable resource and essential foundation for further elucidation of bloom regulators of this and related species of Kareniaceae in the future. IMPORTANCE HABs have increasingly occurred and impacted human health, aquatic ecosystems, and coastal economies. Despite great efforts, the factors that drive the development and termination of a bloom are poorly understood, largely due to inadequate in situ data about the physiology and metabolism of the causal species and the community. Using an integrative molecular ecological approach, we determined that heightened energy and nutrient acquisition promoted the bloom, while resource allocation in defense and failure to defend against grazing and microbial attacks likely prevented or terminated the bloom. Our findings reveal the differential roles of multiple abiotic and biotic environmental factors in driving the formation or demise of a toxic dinoflagellate bloom, suggesting the importance of a balanced biodiverse ecosystem in preventing a dinoflagellate bloom. The study also demonstrates the power of whole-assemblage metatranscriptomics coupled to DNA barcoding in illuminating plankton ecological processes and the underlying species and functional diversities.

Expanding diversity and ecological roles of RNA viruses

While the diversity of global environmental RNA viruses has remained largely unexplored, recent advances have reported on the discovery of over 10⁶ RNA viral contigs from both terrestrial and marine ecosystems that will help us to better understand the diversity, evolution, ecological roles, and transmission of RNA viruses.

Virome analysis in lakes of the South Shetland Islands, Antarctica - 2020

Polar freshwater ecosystems are characterized by a distinct microbiota. However, little is known about viral diversity and abundance, especially regarding the ecology of RNA viruses. We used shotgun metagenomic analysis on samples from Antarctic ecosystems, and report here the characterization of the virome fraction, from different lakes located in the South Shetland Islands (Penguin, Ardley, Deception and King George Island) in the Peninsula Antarctica, in the summer season 2020. DNA viruses (99.4 %) prevailed over RNA viruses (0.6 %) in the lake samples. Six viral orders were identified in the metagenomic libraries: Caudovirales (dsDNA), which was prevalent in most lakes; Picornavirales (ssRNA+); Sobelivirales (ssRNA+); Tolivirales (ssRNA+); Petitvirales (ssDNA) and Baphyvirales (ssDNA), including eight viral families (Herelleviridae, Siphoviridae, Myoviridae, Microviridae, Marnaviridae, Bacilladnaviridae, Barnaviridae and Tombusviridae) and several other, mainly non-classified ssRNA⁽⁺⁾ viruses in the lakes of Ardley Island. Bacteriophages (dsDNA) (Herelleviridae family) infecting the phylum Firmicutes and Siphoviridae were predominant in most lakes evaluated. Functional analysis demonstrated a prevalence of unknown proteins (68 %) in the virome. Our prospective study provides virome analysis data from different lakes in the South Shetland Islands, Antarctica, opening exploratory lines for future research related to the biodiversity and viral ecology in this extreme 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.

Natural Populations from the Phytophthora palustris Complex Show a High Diversity and Abundance of ssRNA and dsRNA Viruses

We explored the virome of the “Phytophthora palustris complex”, a group of aquatic specialists geographically limited to Southeast and East Asia, the native origin of many destructive invasive forest Phytophthora spp. Based on high-throughput sequencing (RNAseq) of 112 isolates of “P. palustris” collected from rivers, mangroves, and ponds, and natural forests in subtropical and tropical areas in Indonesia, Taiwan, and Japan, 52 putative viruses were identified, which, to varying degrees, were phylogenetically related to the families Botybirnaviridae, Narnaviridae, Tombusviridae, and Totiviridae, and the order Bunyavirales. The prevalence of all viruses in their hosts was investigated and confirmed by RT-PCR. The rich virus composition, high abundance, and distribution discovered in our study indicate that viruses are naturally infecting taxa from the “P. palustris complex” in their natural niche, and that they are predominant members of the host cellular environment. Certain Indonesian localities are the viruses’ hotspots and particular “P. palustris” isolates show complex multiviral infections. This study defines the first bi-segmented bunya-like virus together with the first tombus-like and botybirna-like viruses in the genus Phytophthora and provides insights into the spread and evolution of RNA viruses in the natural populations of an oomycete species.

Diverse RNA Viruses Associated with Diatom, Eustigmatophyte, Dinoflagellate, and Rhodophyte Microalgae Cultures

Unicellular microalgae are of immense ecological importance with growing commercial potential in industries such as renewable energy, food, and pharmacology. Viral infections can have a profound impact on the growth and evolution of their hosts. However, very little is known of the diversity within, and the effect of, unicellular microalgal RNA viruses. In addition, identifying RNA viruses in these organisms that could have originated more than a billion years ago constitutes a robust data set to dissect molecular events and address fundamental questions in virus evolution. We assessed the diversity of RNA viruses in eight microalgal cultures, including representatives from the diatom, eustigmatophyte, dinoflagellate, red algae, and euglenid groups. Using metatranscriptomic sequencing combined with bioinformatic approaches optimized to detect highly divergent RNA viruses, we identified 10 RNA virus sequences, with nine constituting new viral species. Most of the newly identified RNA viruses belonged to the double-stranded Totiviridae, Endornaviridae, and Partitiviridae, greatly expanding the reported host range for these families. Two new species belonging to the single-stranded RNA viral clade Marnaviridae, commonly associated with microalgal hosts, were also identified. This study highlights that a substantial diversity of RNA viruses likely exists undetected within the unicellular microalgae. It also highlights the necessity for RNA viral characterization and for investigation of the effects of viral infections on microalgal physiology, biology, and growth, considering their environmental and industrial roles. IMPORTANCE Our knowledge of the diversity of RNA viruses infecting microbial algae-the microalgae-is minimal. However, describing the RNA viruses infecting these organisms is of primary importance at both the ecological and economic scales because of the fundamental roles these organisms play in aquatic environments and their growing value across a range of industrial fields. Using metatranscriptomic sequencing, we aimed to reveal the RNA viruses present in cultures of eight microalgae species belonging to the diatom, dinoflagellate, eustigmatophyte, rhodophyte, and euglena major clades of algae. Accordingly, we identified 10 new divergent RNA virus species belonging to RNA virus families as diverse as the double-stranded Totiviridae, Endornaviridae, and Partitiviridae and the single-stranded Marnaviridae. By expanding the known diversity of RNA viruses infecting unicellular eukaryotes, this study contributes to a better understanding of the early evolution of the virosphere and will inform the use of microalgae in industrial applications.

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.

RdRp-Scan: A Bioinformatic Resource to Identify and Annotate Divergent RNA Viruses in Metagenomic Sequence Data

Despite a rapid expansion in the number of documented viruses following the advent of metagenomic sequencing, the identification and annotation of highly divergent RNA viruses remains challenging, particularly from poorly characterized hosts and environmental samples. Protein structures are more conserved than primary sequence data, such that structure-based comparisons provide an opportunity to reveal the viral “dusk matter”: viral sequences with low, but detectable, levels of sequence identity to known viruses with available protein structures. Here, we present a new open computational and resource – RdRp-scan – that contains a standardized bioinformatic toolkit to identify and annotate divergent RNA viruses in metagenomic sequence data based on the detection of RNA dependent RNA polymerase (RdRp) sequences. By combining RdRp-specific Hidden Markov models (HMM) and structural comparisons we show that RdRp-scan can efficiently detect RdRp sequences with identity levels as low as 10% to those from known viruses and not identifiable using standard sequence-to-sequence comparisons. In addition, to facilitate the annotation and placement of newly detected and divergent virus-like sequences into the diversity of RNA viruses, RdRp-scan provides new custom and curated databases of viral RdRp sequences and core motifs, as well as pre-built RdRp multiple sequence alignments. In parallel, our analysis of the sequence diversity detected by RdRp-scan revealed that while most of the taxonomically unassigned RdRps fell into pre-established clusters, with some falling into potentially new orders of RNA viruses related to the Wolframvirales and Tolivirales. Finally, a survey of the conserved A, B and C RdRp motifs within the RdRp-scan sequence database revealed additional variations of both sequence and position that might provide new insights into the structure, function and evolution of viral polymerases.

ViBE: a hierarchical BERT model to identify eukaryotic viruses using metagenome sequencing data

Viruses are ubiquitous in humans and various environments and continually mutate themselves. Identifying viruses in an environment without cultivation is challenging; however, promoting the screening of novel viruses and expanding the knowledge of viral space is essential. Homology-based methods that identify viruses using known viral genomes rely on sequence alignments, making it difficult to capture remote homologs of the known viruses. To accurately capture viral signals from metagenomic samples, models are needed to understand the patterns encoded in the viral genomes. In this study, we developed a hierarchical BERT model named ViBE to detect eukaryotic viruses from metagenome sequencing data and classify them at the order level. We pre-trained ViBE using read-like sequences generated from the virus reference genomes and derived three fine-tuned models that classify paired-end reads to orders for eukaryotic deoxyribonucleic acid viruses and eukaryotic ribonucleic acid viruses. ViBE achieved higher recall than state-of-the-art alignment-based methods while maintaining comparable precision. ViBE outperformed state-of-the-art alignment-free methods for all test cases. The performance of ViBE was also verified using real sequencing datasets, including the vaginal virome.

Diversity and ecological footprint of Global Ocean RNA viruses

DNA viruses are increasingly recognized as influencing marine microbes and microbe-mediated biogeochemical cycling. However, little is known about global marine RNA virus diversity, ecology, and ecosystem roles. In this study, we uncover patterns and predictors of marine RNA virus community- and "species"-level diversity and contextualize their ecological impacts from pole to pole. Our analyses revealed four ecological zones, latitudinal and depth diversity patterns, and environmental correlates for RNA viruses. Our findings only partially parallel those of cosampled plankton and show unexpectedly high polar ecological interactions. The influence of RNA viruses on ecosystems appears to be large, as predicted hosts are ecologically important. Moreover, the occurrence of auxiliary metabolic genes indicates that RNA viruses cause reprogramming of diverse host metabolisms, including photosynthesis and carbon cycling, and that RNA virus abundances predict ocean carbon export.

Cryptic and abundant marine viruses at the evolutionary origins of Earth's RNA virome

Whereas DNA viruses are known to be abundant, diverse, and commonly key ecosystem players, RNA viruses are insufficiently studied outside disease settings. In this study, we analyzed ≈28 terabases of Global Ocean RNA sequences to expand Earth's RNA virus catalogs and their taxonomy, investigate their evolutionary origins, and assess their marine biogeography from pole to pole. Using new approaches to optimize discovery and classification, we identified RNA viruses that necessitate substantive revisions of taxonomy (doubling phyla and adding >50% new classes) and evolutionary understanding. "Species"-rank abundance determination revealed that viruses of the new phyla "Taraviricota," a missing link in early RNA virus evolution, and "Arctiviricota" are widespread and dominant in the oceans. These efforts provide foundational knowledge critical to integrating RNA viruses into ecological and epidemiological models.

RNA Viruses in Aquatic Ecosystems through the Lens of Ecological Genomics and Transcriptomics

Massive amounts of data from nucleic acid sequencing have changed our perspective about diversity and dynamics of marine viral communities. Here, we summarize recent metatranscriptomic and metaviromic studies targeting predominantly RNA viral communities. The analysis of RNA viromes reaffirms the abundance of lytic (+) ssRNA viruses of the order Picornavirales, but also reveals other (+) ssRNA viruses, including RNA bacteriophages, as important constituents of extracellular RNA viral communities. Sequencing of dsRNA suggests unknown diversity of dsRNA viruses. Environmental metatranscriptomes capture the dynamics of ssDNA, dsDNA, ssRNA, and dsRNA viruses simultaneously, unravelling the full complexity of viral dynamics in the marine environment. RNA viruses are prevalent in large size fractions of environmental metatranscriptomes, actively infect marine unicellular eukaryotes larger than 3 µm, and can outnumber bacteriophages during phytoplankton blooms. DNA and RNA viruses change abundance on hourly timescales, implying viral control on a daily temporal basis. Metatranscriptomes of cultured protists host a diverse community of ssRNA and dsRNA viruses, often with multipartite genomes and possibly persistent intracellular lifestyles. We posit that RNA viral communities might be more diverse and complex than formerly anticipated and that the influence they exert on community composition and global carbon flows in aquatic ecosystems may be underestimated.

Evidence for a novel partitivirus isolated from the entomopathogenic nematode Steinernema ceratophorum

Nematodes are abundant, but little is known about their viruses. In this study, we report a novel partitivirus isolated from the entomopathogenic nematode species Steinernema ceratophorum, named "Steinernema ceratophorum partitivirus 1" (ScPV-1). The complete genome of ScPV-1 comprises two dsRNA segments, dsRNA1 (2352 bp) and dsRNA2 (2196 bp). Each dsRNA contains a single open reading frame (ORF), encoding a putative RNA-dependent RNA polymerase (RdRp) and a coat protein (CP), respectively. The sequences of the RdRp and CP showed the highest similarity (47% and 33% identity, respectively) to Plasmopara viticola associated partitivirus 7 (PvAP-7). A multiple sequence alignment and phylogenetic analysis of the RdRp of ScPV-1 and other selected viruses indicated that ScPV-1 is a new member of the genus Betapartitivirus in the family Partitiviridae.

Eukaryotic Microbial RNA Viruses-Acute or Persistent? Insights into Their Function in the Aquatic Ecosystem

Isolated RNA viruses mainly parasitize eukaryotes. RNA viruses either expand horizontally by infecting hosts (acute type) or coexist with the host and are vertically inherited (persistent type). The significance of persistent-type RNA viruses in environmental viromes (the main hosts are expected to be microbes) was only recently reported because they had previously been overlooked in virology. In this review, we summarize the host-virus relationships of eukaryotic microbial RNA viruses. Picornavirales and Reoviridae are recognized as representative acute-type virus families, and most of the microbial viruses in Narnaviridae, Totiviridae, and Partitiviridae are categorized as representative persistent-type viruses. Acute-type viruses have only been found in aquatic environments, while persistent-type viruses are present in various environments, including aquatic environments. Moreover, persistent-type viruses are potentially widely spread in the RNA viral sequence space. This emerging evidence provides novel insights into RNA viral diversity, host-virus relationships, and their history of co-evolution.

RNA Virosphere in a Marine Zooplankton Community in the Subtropical Western North Pacific

Zooplankton and viruses play a key role in marine ecosystems; however, their interactions have not been examined in detail. In the present study, the diversity of viruses associated with zooplankton collected using a plankton net (mesh size: 100‍ ‍μm) in the subtropical western North Pacific was investigated by fragmented and primer ligated dsRNA sequencing. We obtained 21 and 168 operational taxonomic units (OTUs) of ssRNA and dsRNA viruses, respectively, containing RNA-dependent RNA polymerase (RdRp). These OTUs presented average amino acid similarities of 43.5 and 44.0% to the RdRp genes of known viruses in ssRNA viruses and dsRNA viruses, respectively. Dominant OTUs mainly belonged to narna-like and picorna-like ssRNA viruses and chryso-like, partiti-like, picobirna-like, reo-like, and toti-like dsRNA viruses. Phylogenetic ana-lyses of the RdRp gene revealed that OTUs were phylogenetically diverse and clustered into distinct clades from known viral groups. The community structure of the same zooplankton sample was investigated using small subunit (SSU) rRNA sequences assembled from the metatranscriptome of single-stranded RNA. More than 90% of the sequence reads were derived from metazoan zooplankton; copepods comprised approximately 70% of the sequence reads. Although this ana-lysis provided no direct evidence of the host species of RNA viruses, these dominant zooplankton are expected to be associated with the RNA viruses detected in the present study. The present results indicate that zooplankton function as a reservoir of diverse RNA viruses and suggest that investigations of zooplankton viruses will provide a more detailed understanding of the role of viruses in marine ecosystems.

Putative Invertebrate, Plant, and Wastewater Derived ssRNA Viruses in Plankton of the Anthropogenically Impacted Anacostia River, District of Columbia, USA

The Anacostia River is a highly impacted watershed in the Northeastern United States which experiences combined sewage outfall in downstream waters. We examined the composition of RNA viruses at three sites in the river using viral metagenomics. Viromes had well represented Picornaviruses, Tombusviruses, Wolframviruses, Nodaviruses, with fewer Tobamoviruses, Sobemoviruses, and Densoviruses (ssDNA). Phylogenetic ana-lyses of detected viruses provide evidence for putatively autochthonous and allochthonous invertebrate, plant, and vertebrate host origin. The number of viral genomes matching Ribovaria increased downstream, and assemblages were most disparate between distant sites, suggesting impacts of the combined sewage overflows at these sites. Additionally, we recovered a densovirus genome fragment which was highly similar to the Clinch ambidensovirus 1, which has been attributed to mass mortality of freshwater mussels in Northeastern America. Taken together, these data suggest that RNA viromes of the Anacostia River reflect autochthonous production of virus particles by benthic metazoan and plants, and inputs from terrestrial habitats including sewage.

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.

Virus-derived sequences from the transcriptomes of two snail vectors of schistosomiasis, Biomphalaria pfeifferi and Bulinus globosus from Kenya

Schistosomiasis, which infects more than 230 million people, is vectored by freshwater snails. We identified viral sequences in the transcriptomes of Biomphalaria pfeifferi (BP) and Bulinus globosus (BuG), two of the world's most important schistosomiasis vectors in Africa. Sequences from 26 snails generated using Illumina Hi-Seq or 454 sequencing were assembled using Trinity and CAP3 and putative virus sequences were identified using a bioinformatics pipeline. Phylogenetic analyses were performed using viral RNA-dependent RNA polymerase and coat protein sequences to establish relatedness between virus sequences identified and those of known viruses. Viral sequences were identified from the entire snail holobiont, including symbionts, ingested material and organisms passively associated with the snails. Sequences derived from more than 17 different viruses were found including five near full-length genomes, most of which were small RNA viruses with positive sense RNA genomes (i.e., picorna-like viruses) and some of which are likely derived from adherent or ingested diatoms. Based on phylogenetic analysis, five of these viruses (including BPV2 and BuGV2) along with four Biomphalaria glabrata viruses reported previously, cluster with known invertebrate viruses and are putative viruses of snails. The presence of RNA sequences derived from four of these novel viruses in samples was confirmed. Identification of the genome sequences of candidate snail viruses provides a first step toward characterization of additional gastropod viruses, including from species of biomedical significance.

A High Rate Algal Pond Hosting a Dynamic Community of RNA Viruses

Despite a surge of RNA virome sequencing in recent years, there are still many RNA viruses to uncover-as indicated by the relevance of viral dark matter to RNA virome studies (i.e., putative viruses that do not match to taxonomically identified viruses). This study explores a unique site, a high-rate algal pond (HRAP), for culturing industrially microalgae, to elucidate new RNA viruses. The importance of viral-host interactions in aquatic systems are well documented, and the ever-expanding microalgae industry is no exception. As the industry becomes a more important source of sustainable plastic manufacturing, a producer of cosmetic pigments and alternative protein sources, and a means of CO2 remediation in the face of climate change, studying microalgal viruses becomes a vital practice for proactive management of microalgae cultures at the industrial level. This study provides evidence of RNA microalgal viruses persisting in a CO2 remediation pilot project HRAP and uncovers the diversity of the RNA virosphere contained within it. Evidence shows that family Marnaviridae is cultured in the basin, alongside other potential microalgal infecting viruses (e.g., family Narnaviridae, family Totitiviridae, and family Yueviridae). Finally, we demonstrate that the RNA viral diversity of the HRAP is temporally dynamic across two successive culturing seasons.

Metagenomic Identification of Viral Sequences in Laboratory Reagents

Metagenomic next-generation sequencing has transformed the discovery and diagnosis of infectious disease, with the power to characterise the complete 'infectome' (bacteria, viruses, fungi, parasites) of an individual host organism. However, the identification of novel pathogens has been complicated by widespread microbial contamination in commonly used laboratory reagents. Using total RNA sequencing ("metatranscriptomics") we documented the presence of contaminant viral sequences in multiple 'blank' negative control sequencing libraries that comprise a sterile water and reagent mix. Accordingly, we identified 14 viral sequences in 7 negative control sequencing libraries. As in previous studies, several circular replication-associated protein encoding (CRESS) DNA virus-like sequences were recovered in the blank control libraries, as well as contaminating sequences from the Totiviridae, Tombusviridae and Lentiviridae families of RNA virus. These data suggest that viral contamination of common laboratory reagents is likely commonplace and can comprise a wide variety of viruses.

Novel Dicistroviruses in an Unexpected Wide Range of Invertebrates

Dicistroviruses are members of a rapidly growing family of small RNA viruses. Related sequences have been discovered in many environmental samples, indicating that our knowledge about dicistrovirus diversity and host range is still limited. In this study, we performed a systematic search against the publicly available transcriptome database, and identified large numbers of dicistrovirus-like sequences in a wide variety of eukaryotic species. The origins of these sequences were 108 invertebrates (including 77 insect species belonging to 18 orders) and 11 plants, revealing new associations between dicistroviruses and hosts. Finally, 83 transcripts corresponding to nearly-complete viral genomes were retrieved from the RNA-seq data, of which most sequences showed limited similarity to known dicistroviruses and might present previously unreported virus species. Phylogenetic analysis suggested that horizontal virus transfer has occurred between diverse hosts and has important implications for dicistrovirus evolution. The results will provide new insight into the hidden diversity of the Dicistroviridae, and help us to better understand the viral evolution, host range and the possible way of transmission.

Viral Metagenomic Analysis of Japanese Quail (Coturnix japonica) with Enteritis in the Republic of Korea

We performed viral metagenomics analysis of Japanese quail affected with enteritis to elucidate the viral etiology. Metagenomics generated 21,066,442 sequence reads via high-throughput sequencing, with a mean length of 136 nt. Enrichment in viral sequences suggested that at least three viruses were present in quail samples. Coronavirus and picornavirus were identified and are known as pathogens causing quail enteritis that match the observed morphology. Abundant reads of coronavirus from quail samples yielded four fragment sequences exhibiting six genomes of avian coronavirus. Sequence analysis showed that this quail coronavirus was related to turkey coronavirus and chicken infectious bronchitis virus. Quail picornavirus 8177 bp in size was identified and was similar to the QPV1/HUN/01 virus detected in quails without clinical symptoms in Hungary with 84.6% nucleotide and 94.6% amino acid identity. Our results are useful for understanding the genetic diversity of quail viruses. Further studies must be performed to determine whether quail coronavirus and quail picornavirus are pathogens of the digestive tract of quails.

Metatranscriptomics for the Human Microbiome and Microbial Community Functional Profiling

Shotgun metatranscriptomics (MTX) is an increasingly practical way to survey microbial community gene function and regulation at scale. This review begins by summarizing the motivations for community transcriptomics and the history of the field. We then explore the principles, best practices, and challenges of contemporary MTX workflows: beginning with laboratory methods for isolation and sequencing of community RNA, followed by informatics methods for quantifying RNA features, and finally statistical methods for detecting differential expression in a community context. In thesecond half of the review, we survey important biological findings from the MTX literature, drawing examples from the human microbiome, other (nonhuman) host-associated microbiomes, and the environment. Across these examples, MTX methods prove invaluable for probing microbe-microbe and host-microbe interactions, the dynamics of energy harvest and chemical cycling, and responses to environmental stresses. We conclude with a review of open challenges in the MTX field, including making assays and analyses more robust, accessible, and adaptable to new technologies; deciphering roles for millions of uncharacterized microbial transcripts; and solving applied problems such as biomarker discovery and development of microbial therapeutics.

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.

Role of Phylogenetic Structure in the Dynamics of Coastal Viral Assemblages

Marine microbes, including viruses, are an essential part of the marine ecosystem, forming the base of the food web and driving biogeochemical cycles. Within this system, the composition of viral assemblages changes markedly with time, and some of these changes are repeatable through time; however, the extent to which these dynamics are reflected within versus among evolutionarily related groups of viruses is largely unexplored. To examine these dynamics, changes in the composition of two groups of ecologically important viruses and communities of their potential hosts were sampled every 2 weeks for 13 months at a coastal site in British Columbia, Canada. We sequenced two marker genes for viruses-the gene encoding the major capsid protein of T4-like phages and their relatives (gp23) and the RNA-dependent RNA polymerase (RdRp) gene of marnavirus-like RNA viruses-as well as marker genes for their bacterial and eukaryotic host communities, the genes encoding 16S rRNA and 18S rRNA. There were strong lagged correlations between viral diversity and community similarity of putative hosts, implying that the viruses influenced the composition of the host communities. The results showed that for both viral assemblages, the dominant clusters of phylogenetically related viruses shifted over time, and this was correlated with environmental changes. Viral clusters contained many ephemeral taxa and few persistent taxa, but within a viral assemblage, the ephemeral and persistent taxa were closely related, implying ecological dynamics within these clusters. Furthermore, these dynamics occurred in both the RNA and DNA viral assemblages surveyed, implying that this structure is common in natural viral assemblages.IMPORTANCE Viruses are major agents of microbial mortality in marine systems, yet little is known about changes in the composition of viral assemblages in relation to those of the microbial communities that they infect. Here, we sampled coastal seawater every 2 weeks for 1 year and used high-throughput sequencing of marker genes to follow changes in the composition of two groups of ecologically important viruses, as well as the communities of bacteria and protists that serve as their respective hosts. Different subsets of genetically related viruses dominated at different times. These results demonstrate that although the genetic composition of viral assemblages is highly dynamic temporally, for the most part the shuffling of genotypes occurs within a few clusters of phylogenetically related viruses. Thus, it appears that even in temperate coastal waters with large seasonal changes, the highly dynamic shuffling of viral genotypes occurs largely within a few subsets of related individuals.

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.

RNA Viruses in Aquatic Unicellular Eukaryotes

Increasing sequence information indicates that RNA viruses constitute a major fraction of marine virus assemblages. However, only 12 RNA virus species have been described, infecting known host species of marine single-celled eukaryotes. Eight of these use diatoms as hosts, while four are resident in dinoflagellate, raphidophyte, thraustochytrid, or prasinophyte species. Most of these belong to the order Picornavirales, while two are divergent and fall into the families Alvernaviridae and Reoviridae. However, a very recent study has suggested that there is extraordinary diversity in aquatic RNA viromes, describing thousands of viruses, many of which likely use protist hosts. Thus, RNA viruses are expected to play a major ecological role for marine unicellular eukaryotic hosts. In this review, we describe in detail what has to date been discovered concerning viruses with RNA genomes that infect aquatic unicellular eukaryotes.

Virome Variation during Sea Star Wasting Disease Progression in Pisaster ochraceus (Asteroidea, Echinodermata)

Sea star wasting disease (SSWD) is a condition that has affected asteroids for over 120 years, yet mechanistic understanding of this wasting etiology remains elusive. We investigated temporal virome variation in two Pisaster ochraceus specimens that wasted in the absence of external stimuli and two specimens that did not experience SSWD for the duration of our study, and compared viromes of wasting lesion margin tissues to both artificial scar margins and grossly normal tissues over time. Global assembly of all SSWD-affected tissue libraries resulted in 24 viral genome fragments represented in >1 library. Genome fragments mostly matched densoviruses and picornaviruses with fewer matching nodaviruses, and a sobemovirus. Picornavirus-like and densovirus-like genome fragments were most similar to viral genomes recovered in metagenomic study of other marine invertebrates. Read recruitment revealed only two picornavirus-like genome fragments that recruited from only SSWD-affected specimens, but neither was unique to wasting lesions. Wasting lesion margin reads recruited to a greater number of viral genotypes (i.e., richness) than did either scar tissue and grossly normal tissue reads. Taken together, these data suggest that no single viral genome fragment was associated with SSWD. Rather, wasting lesion margins may generally support viral proliferation.

Diverged and Active Partitiviruses in Lichen

The lichen is a microbial consortium that mainly consists of fungi and either algae (Viridiplantae) or cyanobacteria. This structure also contains other bacteria, fungi, and viruses. However, RNA virus diversity associated with lichens is still unknown. Here, we analyzed RNA virus diversity in a lichen dominated by fungi and algae using dsRNA-seq technology and revealed that partitiviruses were dominant and active in the microbial consortium. The Partitiviridae sequences found in this study were classified into two genera, which have both plant- and fungi-infecting partitiviruses. This observation suggests that the lichen provides an opportunity for horizontal transfer of these partitiviruses among microbes that form the lichen consortium.

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

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

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

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

Illuminating an Ecological Blackbox: Using High Throughput Sequencing to Characterize the Plant Virome Across Scales

The ecology of plant viruses began to be explored at the end of the 19th century. Since then, major advances have revealed mechanisms of virus-host-vector interactions in various environments. These advances have been accelerated by new technlogies for virus detection and characterization, most recently including high throughput sequencing (HTS). HTS allows investigators, for the first time, to characterize all or nearly all viruses in a sample without a priori information about which viruses might be present. This powerful approach has spurred new investigation of the viral metagenome (virome). The rich virome datasets accumulated illuminate important ecological phenomena such as virus spread among host reservoirs (wild and domestic), effects of ecosystem simplification caused by human activities (and agriculture) on the biodiversity and the emergence of new viruses in crops. To be effective, however, HTS-based virome studies must successfully navigate challenges and pitfalls at each procedural step, from plant sampling to library preparation and bioinformatic analyses. This review summarizes major advances in plant virus ecology associated with technological developments, and then presents important considerations and best practices for HTS use in virome studies.

An Unconventional Flavivirus and Other RNA Viruses in the Sea Cucumber (Holothuroidea; Echinodermata) Virome

Sea cucumbers (Holothuroidea; Echinodermata) are ecologically significant constituents of benthic marine habitats. We surveilled RNA viruses inhabiting eight species (representing four families) of holothurian collected from four geographically distinct locations by viral metagenomics, including a single specimen of Apostichopus californicus affected by a hitherto undocumented wasting disease. The RNA virome comprised genome fragments of both single-stranded positive sense and double stranded RNA viruses, including those assigned to the Picornavirales, Ghabrivirales, and Amarillovirales. We discovered an unconventional flavivirus genome fragment which was most similar to a shark virus. Ghabivirales-like genome fragments were most similar to fungal totiviruses in both genome architecture and homology and had likely infected mycobiome constituents. Picornavirales, which are commonly retrieved in host-associated viral metagenomes, were similar to invertebrate transcriptome-derived picorna-like viruses. The greatest number of viral genome fragments was recovered from the wasting A. californicus library compared to the asymptomatic A. californicus library. However, reads from the asymptomatic library recruited to nearly all recovered wasting genome fragments, suggesting that they were present but not well represented in the grossly normal specimen. These results expand the known host range of flaviviruses and suggest that fungi and their viruses may play a role in holothurian ecology.

High throughput shotgun sequencing of eRNA reveals taxonomic and derived functional shifts across a benthic productivity gradient

Benthic macrofauna is regularly used in monitoring programmes, however the vast majority of benthic eukaryotic biodiversity lies mostly in microscopic organisms, such as meiofauna (invertebrates < 1 mm) and protists, that rapidly responds to environmental change. These communities have traditionally been hard to sample and handle in the laboratory, but DNA sequencing has made such work less time consuming. While DNA sequencing captures both alive and dead organisms, environmental RNA (eRNA) better targets living organisms or organisms of recent origin in the environment. Here, we assessed the biodiversity of three known bioindicator microeukaryote groups (nematodes, foraminifera, and ciliates) in sediment samples collected at seven coastal sites along an organic carbon (OC) gradient. We aimed to investigate if eRNA shotgun sequencing can be used to simultaneously detect differences in (i) biodiversity of multiple microeukaryotic communities; and (ii) functional feeding traits of nematodes. Results showed that biodiversity was lower for nematodes and foraminifera in high OC (6.2%-6.9%), when compared to low OC sediments (1.2%-2.8%). Dissimilarity in community composition increased for all three groups between Low OC and High OC, as well as the classified feeding type of nematode genera (with more nonselective deposit feeders in high OC sediment). High relative abundant genera included nematode Sabatieria and foraminifera Elphidium in high OC, and Cryptocaryon-like ciliates in low OC sediments. Considering that future sequencing technologies are likely to decrease in cost, the use of eRNA shotgun sequencing to assess biodiversity of benthic microeukaryotes could be a powerful tool in recurring monitoring programmes.

Coronaviruses in the Sea

Interest in coronaviruses because of the 2019 novel coronavirus (SARS-CoV-2) pandemic has generated concern about their occurrence and persistence in aquatic habitats. Coronaviruses are not quantitatively significant constituents of marine virioplankton. Members of the Nidovirales (to which human coronaviruses belong) infect marine mammals, teleosts and possibly invertebrates, and human coronaviruses may persist in marine plankton receiving wastewater effluent. However, virions likely experience significant particle and infectivity decay rates in surface seawater, similar to other enveloped RNA viruses.

Reconstructing ribosomal genes from large scale total RNA meta-transcriptomic data

MOTIVATION

Technological advances in meta-transcriptomics have enabled a deeper understanding of the structure and function of microbial communities. 'Total RNA' meta-transcriptomics, sequencing of total reverse transcribed RNA, provides a unique opportunity to investigate both the structure and function of active microbial communities from all three domains of life simultaneously. A major step of this approach is the reconstruction of full-length taxonomic marker genes such as the small subunit ribosomal RNA. However, current tools for this purpose are mainly targeted towards analysis of amplicon and metagenomic data and thus lack the ability to handle the massive and complex datasets typically resulting from total RNA experiments.

RESULTS

In this work, we introduce MetaRib, a new tool for reconstructing ribosomal gene sequences from total RNA meta-transcriptomic data. MetaRib is based on the popular rRNA assembly program EMIRGE, together with several improvements. We address the challenge posed by large complex datasets by integrating sub-assembly, dereplication and mapping in an iterative approach, with additional post-processing steps. We applied the method to both simulated and real-world datasets. Our results show that MetaRib can deal with larger datasets and recover more rRNA genes, which achieve around 60 times speedup and higher F1 score compared to EMIRGE in simulated datasets. In the real-world dataset, it shows similar trends but recovers more contigs compared with a previous analysis based on random sub-sampling, while enabling the comparison of individual contig abundances across samples for the first time.

AVAILABILITY AND IMPLEMENTATION

The source code of MetaRib is freely available at https://github.com/yxxue/MetaRib.

CONTACT

yaxin.xue@uib.no or Inge.Jonassen@uib.no.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Capsid Structure of a Marine Algal Virus of the Order Picornavirales

The order Picornavirales includes viruses that infect different kinds of eukaryotes and that share similar properties. The capsid proteins (CPs) of viruses in the order that infect unicellular organisms, such as algae, presumably possess certain characteristics that have changed little over the course of evolution, and thus these viruses may resemble the Picornavirales ancestor in some respects. Herein, we present the capsid structure of Chaetoceros tenuissimus RNA virus type II (CtenRNAV-II) determined using cryo-electron microscopy at a resolution of 3.1 Å, the first alga virus belonging to the family Marnaviridae of the order Picornavirales A structural comparison to related invertebrate and vertebrate viruses revealed a unique surface loop of the major CP VP1 that had not been observed previously, and further, revealed that another VP1 loop obscures the so-called canyon, which is a host-receptor binding site for many of the mammalian Picornavirales viruses. VP2 has an N-terminal tail, which has previously been reported as a primordial feature of Picornavirales viruses. The above-mentioned and other critical structural features provide new insights on three long-standing theories about Picornavirales: (i) the canyon hypothesis, (ii) the primordial VP2 domain swap, and (iii) the hypothesis that alga Picornavirales viruses could share characteristics with the Picornavirales ancestor.IMPORTANCE Identifying the acquired structural traits in virus capsids is important for elucidating what functions are essential among viruses that infect different hosts. The Picornavirales viruses infect a broad spectrum of hosts, ranging from unicellular algae to insects and mammals and include many human pathogens. Those viruses that infect unicellular protists, such as algae, are likely to have undergone fewer structural changes during the course of evolution compared to those viruses that infect multicellular eukaryotes and thus still share some characteristics with the Picornavirales ancestor. This article describes the first atomic capsid structure of an alga Marnavirus, CtenRNAV-II. A comparison to capsid structures of the related invertebrate and vertebrate viruses identified a number of structural traits that have been functionally acquired or lost during the course of evolution. These observations provide new insights on past theories on the viability and evolution of Picornavirales viruses.

High-throughput sequencing (HTS) for the analysis of viral populations

The development of High-Throughput Sequencing (HTS) technologies is having a major impact on the genomic analysis of viral populations. Current HTS platforms can capture nucleic acid variation across millions of genes for both selected amplicons and full viral genomes. HTS has already facilitated the discovery of new viruses, hinted new taxonomic classifications and provided a deeper and broader understanding of their diversity, population and genetic structure. Hence, HTS has already replaced standard Sanger sequencing in basic and applied research fields, but the next step is its implementation as a routine technology for the analysis of viruses in clinical settings. The most likely application of this implementation will be the analysis of viral genomics, because the huge population sizes, high mutation rates and very fast replacement of viral populations have demonstrated the limited information obtained with Sanger technology. In this review, we describe new technologies and provide guidelines for the high-throughput sequencing and genetic and evolutionary analyses of viral populations and metaviromes, including software applications. With the development of new HTS technologies, new and refurbished molecular and bioinformatic tools are also constantly being developed to process and integrate HTS data. These allow assembling viral genomes and inferring viral population diversity and dynamics. Finally, we also present several applications of these approaches to the analysis of viral clinical samples including transmission clusters and outbreak characterization.

Novel RNA viruses associated with Plasmodium vivax in human malaria and Leucocytozoon parasites in avian disease

Eukaryotes of the genus Plasmodium cause malaria, a parasitic disease responsible for substantial morbidity and mortality in humans. Yet, the nature and abundance of any viruses carried by these divergent eukaryotic parasites is unknown. We investigated the Plasmodium virome by performing a meta-transcriptomic analysis of blood samples taken from patients suffering from malaria and infected with P. vivax, P. falciparum or P. knowlesi. This resulted in the identification of a narnavirus-like sequence, encoding an RNA polymerase and restricted to P. vivax samples, as well as an associated viral segment of unknown function. These data, confirmed by PCR, are indicative of a novel RNA virus that we term Matryoshka RNA virus 1 (MaRNAV-1) to reflect its analogy to a "Russian doll": a virus, infecting a parasite, infecting an animal. Additional screening revealed that MaRNAV-1 was abundant in geographically diverse P. vivax derived from humans and mosquitoes, strongly supporting its association with this parasite, and not in any of the other Plasmodium samples analyzed here nor Anopheles mosquitoes in the absence of Plasmodium. Notably, related bi-segmented narnavirus-like sequences (MaRNAV-2) were retrieved from Australian birds infected with a Leucocytozoon—a genus of eukaryotic parasites that group with Plasmodium in the Apicomplexa subclass hematozoa. Together, these data support the establishment of two new phylogenetically divergent and genomically distinct viral species associated with protists, including the first virus likely infecting Plasmodium parasites. As well as broadening our understanding of the diversity and evolutionary history of the eukaryotic virosphere, the restriction to P. vivax may be of importance in understanding P. vivax-specific biology in humans and mosquitoes, and how viral co-infection might alter host responses at each stage of the P. vivax life-cycle.

While parasites are a major cause of human disease, they can themselves be infected by viruses. We asked whether three of the major malaria-causing parasites in humans—Plasmodium vivax, P. falciparum and P. knowlesi—were also infected by viruses. To this end we performed total RNA-Sequencing (“meta-transcriptomics”) on human blood samples infected with these Plasmodium species. This resulted in the discovery of an abundant bi-segmented virus—Matryoshka RNA virus 1 (MaRNAV-1)—in all P. vivax samples tested (but no other Plasmodium species) that contains a replicase segment related to those of narnaviruses, arguably the simplest type of RNA viruses discovered to date. By screening for MaRNAV-1 in a larger set of Plasmodium species we revealed a strong specificity between this virus and P. vivax, as well as the presence of a related virus—MaRNAV-2—in avian Leucocytozoon hematozoa parasites. This is the first discovery of a Plasmodium-associated virus and will assist in revealing the deep evolutionary history of RNA viruses and our understanding of Plasmodium biology and disease processes.

Next-Generation Sequencing on Insectivorous Bat Guano: An Accurate Tool to Identify Arthropod Viruses of Potential Agricultural Concern

Viruses belonging to the Dicistroviridae family have attracted a great deal of attention from scientists owing to their negative impact on agricultural economics, as well as their recent identification as potential aetiological agents of febrile illness in human patients. On the other hand, some Dicistroviruses are also studied for their potential biopesticide properties. To date, Dicistrovirus characterized in African mainland remain scarce. By using High-Throughput Sequencing technology on insectivorous bat faeces (Hipposideros Caffer) sampled in a cave used by humans to collect bat guano (bat manure) as fertilizer in Zimbabwe, we characterized the full-length sequences of three Dicistrovirus belonging to the Cripavirus and Aparavirus genus: Big Sioux River Virus-Like (BSRV-Like), Acute Bee Paralysis Virus (ABPV), and Aphid Lethal Paralysis Virus (ALPV). Phylogenetic analyses of ORF-1 and ORF-2 genes showed a complex evolutionary history between BSRV and close viruses, as well as for the Aparavirus genus. Herewith, we provide the first evidence of the presence of Dicistrovirus in Zimbabwe and highlight the need to further document the impact of such viruses on crops, as well as in beekeeping activities in Zimbabwe which represent a crucial source of income for Zimbabwean people.

Hantaviridae: Current Classification and Future Perspectives

In recent years, negative-sense RNA virus classification and taxon nomenclature have undergone considerable transformation. In 2016, the new order Bunyavirales was established, elevating the previous genus Hantavirus to family rank, thereby creating Hantaviridae. Here we summarize affirmed taxonomic modifications of this family from 2016 to 2019. Changes involve the admission of >30 new hantavirid species and the establishment of subfamilies and novel genera based on DivErsity pArtitioning by hieRarchical Clustering (DEmARC) analysis of genomic sequencing data. We outline an objective framework that can be used in future classification schemes when more hantavirids sequences will be available. Finally, we summarize current taxonomic proposals and problems in hantavirid taxonomy that will have to be addressed shortly.

Metagenomic Analysis of the Diversity of DNA Viruses in the Surface and Deep Sea of the South China Sea

A metagenomic analysis of the viral community from five surface and five deep sea water (>2000 m below the surface, mbs) samples collected from the central basin of the South China Sea and adjacent Northwest Pacific Ocean during July-August 2017 was conducted herein. We builded up a South China Sea DNA virome (SCSV) dataset of 29,967 viral Operational Taxonomic Units (vOTUs), which is comparable to the viral populations from the original Tara Ocean and Malaspina expeditions. The most abundant and widespread viral populations were from the uncultivated viruses annotated from the viral metagenomics. Only 74 and 37 vOTUs have similarity with the reported genomes from the cultivated viruses and the single-virus genomics, respectively. The community structures of deep sea viromes in the SCSV were generally different from the surface viromes. The carbon flux and nutrients (PO₄ and NOx) were related to the surface and deep sea viromes in the SCSV, respectively. In the SCSV, the annotated vOTUs could be affiliated to the cultivated viruses mainly including Pelagibacter (SAR11) phage HTVC010P, Prochlorococcus phages (P-GSP1, P-SSM4, and P-TIM68), Cyanophages (MED4-184 and MED4-117) and Mycobacterium phages (Sparky and Squirty). It indicated that phage infection to the SAR11 cluster may occur ubiquitously and has significant impacts on bathypelagic SAR11 communities in the deep sea. Meanwhile, as Prochlorococcus is prominently distributed in the euphotic ocean, the existence of their potential phages in the deep sea suggested the sedimentation mechanism might contribute to the formation of the deep sea viromes. Intriguingly, the presence of Mycobacterium phages only in the deep sea viromes, suggests inhabitance of endemic viral populations in the deep sea viromes in the SCSV. This study provided an insight of the viral community in the South China Sea and for the first time uncovered the deep sea viral diversity in the central basin of the South China Sea.

Expanding the RNA Virosphere by Unbiased Metagenomics

Although viruses comprise the most abundant genetic material in the biosphere, to date only several thousand virus species have been formally defined. Such a limited perspective on virus diversity has in part arisen because viruses were traditionally considered only as etiologic agents of overt disease in humans or economically important species and were often difficult to identify using cell culture. This view has dramatically changed with the rise of metagenomics, which is transforming virus discovery and revealing a remarkable diversity of viruses sampled from diverse cellular organisms. These newly discovered viruses help fill major gaps in the evolutionary history of viruses, revealing a near continuum of diversity among genera, families, and even orders of RNA viruses. Herein, we review some of the recent advances in our understanding of the RNA virosphere that have stemmed from metagenomics, note future directions, and highlight some of the remaining challenges to this rapidly developing field.

Marine DNA Viral Macro- and Microdiversity from Pole to Pole

Microbes drive most ecosystems and are modulated by viruses that impact their lifespan, gene flow, and metabolic outputs. However, ecosystem-level impacts of viral community diversity remain difficult to assess due to classification issues and few reference genomes. Here, we establish an ∼12-fold expanded global ocean DNA virome dataset of 195,728 viral populations, now including the Arctic Ocean, and validate that these populations form discrete genotypic clusters. Meta-community analyses revealed five ecological zones throughout the global ocean, including two distinct Arctic regions. Across the zones, local and global patterns and drivers in viral community diversity were established for both macrodiversity (inter-population diversity) and microdiversity (intra-population genetic variation). These patterns sometimes, but not always, paralleled those from macro-organisms and revealed temperate and tropical surface waters and the Arctic as biodiversity hotspots and mechanistic hypotheses to explain them. Such further understanding of ocean viruses is critical for broader inclusion in ecosystem models.

Marine RNA Virus Quasispecies Are Distributed throughout the Oceans

Very little is known about aquatic RNA virus populations and genome evolution. This is the first study that analyzes marine environmental RNA viral assemblages in an evolutionary and broad geographical context. This study contributes the largest marine RNA virus metagenomic data set to date, substantially increasing the sequencing space for RNA viruses and also providing a baseline for comparisons of marine RNA virus diversity. The new viruses discovered in this study are representative of the most abundant family of marine RNA viruses, the Marnaviridae, and expand our view of the diversity of this important group. Overall, our data and analyses provide a foundation for interpreting marine RNA virus diversity and evolution.

RNA viruses, particularly genetically diverse members of the Picornavirales, are widespread and abundant in the ocean. Gene surveys suggest that there are spatial and temporal patterns in the composition of RNA virus assemblages, but data on their diversity and genetic variability in different oceanographic settings are limited. Here, we show that specific RNA virus genomes have widespread geographic distributions and that the dominant genotypes are under purifying selection. Genomes from three previously unknown picorna-like viruses (BC-1, -2, and -3) assembled from a coastal site in British Columbia, Canada, as well as marine RNA viruses JP-A, JP-B, and Heterosigma akashiwo RNA virus exhibited different biogeographical patterns. Thus, biotic factors such as host specificity and viral life cycle, and not just abiotic processes such as dispersal, affect marine RNA virus distribution. Sequence differences relative to reference genomes imply that virus quasispecies are under purifying selection, with synonymous single-nucleotide variations dominating in genomes from geographically distinct regions resulting in conservation of amino acid sequences. Conversely, sequences from coastal South Africa that mapped to marine RNA virus JP-A exhibited more nonsynonymous mutations, probably representing amino acid changes that accumulated over a longer separation. This biogeographical analysis of marine RNA viruses demonstrates that purifying selection is occurring across oceanographic provinces. These data add to the spectrum of known marine RNA virus genomes, show the importance of dispersal and purifying selection for these viruses, and indicate that closely related RNA viruses are pathogens of eukaryotic microbes across oceans.

IMPORTANCE Very little is known about aquatic RNA virus populations and genome evolution. This is the first study that analyzes marine environmental RNA viral assemblages in an evolutionary and broad geographical context. This study contributes the largest marine RNA virus metagenomic data set to date, substantially increasing the sequencing space for RNA viruses and also providing a baseline for comparisons of marine RNA virus diversity. The new viruses discovered in this study are representative of the most abundant family of marine RNA viruses, the Marnaviridae, and expand our view of the diversity of this important group. Overall, our data and analyses provide a foundation for interpreting marine RNA virus diversity and evolution.

Viruses of Polar Aquatic Environments

The poles constitute 14% of the Earth’s biosphere: The aquatic Arctic surrounded by land in the north, and the frozen Antarctic continent surrounded by the Southern Ocean. In spite of an extremely cold climate in addition to varied topographies, the polar aquatic regions are teeming with microbial life. Even in sub-glacial regions, cellular life has adapted to these extreme environments where perhaps there are traces of early microbes on Earth. As grazing by macrofauna is limited in most of these polar regions, viruses are being recognized for their role as important agents of mortality, thereby influencing the biogeochemical cycling of nutrients that, in turn, impact community dynamics at seasonal and spatial scales. Here, we review the viral diversity in aquatic polar regions that has been discovered in the last decade, most of which has been revealed by advances in genomics-enabled technologies, and we reflect on the vast extent of the still-to-be explored polar microbial diversity and its “enigmatic virosphere”.