dsRNA-seq Reveals Novel RNA Virus and Virus-Like Putative Complete Genome Sequences from Hymeniacidon sp. Sponge

Double-stranded RNA sequencing reveals distinct riboviruses associated with thermoacidophilic bacteria from hot springs in Japan

Metatranscriptome sequencing expanded the known diversity of the bacterial RNA virome, suggesting that additional riboviruses infecting bacterial hosts remain to be discovered. Here we employed double-stranded RNA sequencing to recover complete genome sequences of two ribovirus groups from acidic hot springs in Japan. One group, denoted hot spring riboviruses (HsRV), consists of viruses with distinct RNA-directed RNA polymerases (RdRPs) that seem to be intermediates between typical ribovirus RdRPs and viral reverse transcriptases. This group forms a distinct phylum, Artimaviricota, or even kingdom within the realm Riboviria. We identified viruses encoding HsRV-like RdRPs in marine water, river sediments and salt marshes, indicating that this group is widespread beyond extreme ecosystems. The second group, denoted hot spring partiti-like viruses (HsPV), forms a distinct branch within the family Partitiviridae. The genome architectures of HsRV and HsPV and their identification in bacteria-dominated habitats suggest that these viruses infect thermoacidophilic bacteria.

Virus-Derived Small RNAs and microRNAs in Health and Disease

MicroRNAs (miRNAs) are short noncoding RNAs that can regulate all steps of gene expression (induction, transcription, and translation). Several virus families, primarily double-stranded DNA viruses, encode small RNAs (sRNAs), including miRNAs. These virus-derived miRNAs (v-miRNAs) help the virus evade the host's innate and adaptive immune system and maintain an environment of chronic latent infection. In this review, the functions of the sRNA-mediated virus-host interactions are highlighted, delineating their implication in chronic stress, inflammation, immunopathology, and disease. We provide insights into the latest viral RNA-based research-in silico approaches for functional characterization of v-miRNAs and other RNA types. The latest research can assist toward the identification of therapeutic targets to combat viral infections.

Distinct groups of RNA viruses associated with thermoacidophilic bacteria

Recent massive metatranscriptome mining substantially expanded the diversity of the bacterial RNA virome, suggesting that additional groups of riboviruses infecting bacterial hosts remain to be discovered. We employed full length double-stranded (ds) RNA sequencing for identification of riboviruses associated with microbial consortia dominated by bacteria and archaea in acidic hot springs in Japan. Whole sequences of two groups of multisegmented riboviruses genomes were obtained. One group, which we denoted hot spring riboviruses (HsRV), consists of unusual viruses with distinct RNA-dependent RNA polymerases (RdRPs) that seem to be intermediates between typical ribovirus RdRPs and viral reverse transcriptases. We also identified viruses encoding HsRV-like RdRPs in moderate aquatic environments, including marine water, river sediments and salt marsh, indicating that this previously overlooked ribovirus group is not restricted to the extreme ecosystem. The HsRV-like viruses are candidates for a distinct phylum or even kingdom within the viral realm Riboviria. The second group, denoted hot spring partiti-like viruses (HsPV), is a distinct branch within the family Partitiviridae. All genome segments in both these groups of viruses display the organization typical of bacterial riboviruses, where multiple open reading frames encoding individual proteins are preceded by ribosome-binding sites. Together with the identification in bacteria-dominated habitats, this genome architecture indicates that riboviruses of these distinct groups infect thermoacidophilic bacterial hosts.

NeoRdRp: A Comprehensive Dataset for Identifying RNA-dependent RNA Polymerases of Various RNA Viruses from Metatranscriptomic Data

RNA viruses are distributed throughout various environments, and most have recently been identified by metatranscriptome sequencing. However, due to the high nucleotide diversity of RNA viruses, it is still challenging to identify novel RNA viruses from metatranscriptome data. To overcome this issue, we created a dataset of RNA-dependent RNA polymerase (RdRp) domains that are essential for all RNA viruses belonging to Orthornavirae. Genes with RdRp domains from various RNA viruses were clustered based on amino acid sequence similarities. A multiple sequence alignment was generated for each cluster, and a hidden Markov model (HMM) profile was created when the number of sequences was greater than three. We further refined 426 HMM profiles by detecting RefSeq RNA virus sequences and subsequently combined the hit sequences with the RdRp domains. As a result, 1,182 HMM profiles were generated from 12,502 RdRp domain sequences, and the dataset was named NeoRdRp. The majority of NeoRdRp HMM profiles successfully detected RdRp domains, specifically in the UniProt dataset. Furthermore, we compared the NeoRdRp dataset with two previously reported methods for RNA virus detection using metatranscriptome sequencing data. Our methods successfully identified the majority of RNA viruses in the datasets; however, some RNA viruses were not detected, similar to the other two methods. NeoRdRp may be repeatedly improved by the addition of new RdRp sequences and is applicable as a system for detecting various RNA viruses from diverse metatranscriptome data.

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.

Metagenomic Assessment of DNA Viral Diversity in Freshwater Sponges, Baikalospongia bacillifera

Sponges (type Porifera) are multicellular organisms that give shelter to a variety of microorganisms: fungi, algae, archaea, bacteria, and viruses. The studies concerning the composition of viral communities in sponges have appeared rather recently, and the diversity and role of viruses in sponge holobionts remain largely undisclosed. In this study, we assessed the diversity of DNA viruses in the associated community of the Baikal endemic sponge, Baikalospongia bacillifera, using a metagenomic approach, and compared the virome data from samples of sponges and Baikal water (control sample). Significant differences in terms of taxonomy, putative host range of identified scaffolds, and functional annotation of predicted viral proteins were revealed in viromes of sponge B. bacillifera and the Baikal water. This is the evidence in favor of specificity of viral communities in sponges. The diversity shift of viral communities in a diseased specimen, in comparison with a visually healthy sponge, probably reflects the changes in the composition of microbial communities in affected sponges. We identified many viral genes encoding the proteins with metabolic functions; therefore, viruses in Baikal sponges regulate the number and diversity of their associated community, and also take a part in the vital activity of the holobiont, and this is especially significant in the case of damage (or disease) of these organisms in unfavorable conditions. When comparing the Baikal viromes with similar datasets of marine sponge (Ianthella basta), in addition to significant differences in the taxonomic and functional composition of viral communities, we revealed common scaffolds/virotypes in the cross-assembly of reads, which may indicate the presence of some closely related sponge-specific viruses in marine and freshwater sponges.

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.

Detection and Characterization of RNA Viruses in Red Macroalgae (Bangiaceae) and Their Food Product (Nori Sheets)

Persistent RNA viruses, which have been suggested to form symbiotic relationships with their hosts, have been reported to occur in eukaryotes, such as plants, fungi, and algae. Based on empirical findings, these viruses may also be present in commercially cultivated macroalgae. Accordingly, the present study aimed to screen red macroalgae (family Bangiaceae conchocelis and Neopyropia yezoensis thallus) and processed nori sheets (N. yezoensis) for persistent RNA viruses using fragmented and primer-ligated dsRNA sequencing (FLDS) and targeted reverse transcription PCR (RT-PCR). A Totiviridae-related virus was detected in the conchocelis of Neoporphyra haitanensis, which is widely cultivated in China, while two Mitoviridae-related viruses were found in several conchocelis samples and all N. yezoensis-derived samples (thallus and nori sheets). Mitoviridae-related viruses in N. yezoensis are widespread among cultivated species and not expected to inhibit host growth. Mitoviridae-related viruses were also detected in several phylogenetically distant species in the family Bangiaceae, which suggests that these viruses persisted and coexist in the family Bangiaceae over a long period of time. The present study is the first to report persistent RNA viruses in nori sheets and their raw materials.

RNA Viral Metagenome Analysis of Subnanogram dsRNA Using Fragmented and Primer Ligated dsRNA Sequencing (FLDS)

Fragmented and primer ligated dsRNA sequencing (FLDS) is a sequencing method applicable to long double-stranded RNA (dsRNA) that enables the complete genome sequencing of both double- and single-stranded RNA viruses. However, the application of this method on a low amount of dsRNA has been hindered by adaptor dimer formation during cDNA amplification and sequence library preparation. We herein developed FLDS ver. 3 by optimizing the terminal modification of an oligonucleotide adaptor and the conditions of adaptor ligation. We also examined the concentration of Mg²⁺ in the PCR reaction for cDNA amplification and the purification method of amplified cDNA. Fine sequence reads were successfully obtained from metagenomic shotgun sequencing libraries constructed from 10 and 100 pg dsRNA, and these libraries exhibited weaker detection sensitivity for low-abundance dsRNAs (viral genomes and genome segments) than that constructed from 1‍ ‍ng of dsRNA. We also report the utility of capillary electrophoresis for dsRNA quantification. The FLDS ver. 3 package expands the frontiers of our knowledge in RNA virus diversity and evolution.

Viral RNA Genomes Identified from Marine Macroalgae and a Diatom

Protists provide insights into the diversity and function of RNA viruses in marine systems. Among them, marine macroalgae are good targets for RNA virome analyses because they have a sufficient biomass in nature. However, RNA viruses in macroalgae have not yet been examined in detail, and only partial genome sequences have been reported for the majority of RNA viruses. Therefore, to obtain further insights into the distribution and diversity of RNA viruses associated with marine protists, we herein examined RNA viruses in macroalgae and a diatom. We report the putative complete genome sequences of six novel RNA viruses from two marine macroalgae and one diatom holobiont. Four viruses were not classified into established viral genera or families. Furthermore, a virus classified into Totiviridae showed a genome structure that has not yet been reported in this family. These results suggest that a number of distinct RNA viruses are widespread in a broad range of protists.

Discovery of divided RdRp sequences and a hitherto unknown genomic complexity in fungal viruses

By identifying variations in viral RNA genomes, cutting-edge metagenome technology has potential to reshape current concepts about the evolution of RNA viruses. This technology, however, cannot process low-homology genomic regions properly, leaving the true diversity of RNA viruses unappreciated. To overcome this technological limitation, we applied an advanced method, Fragmented and Primer-Ligated Double-stranded (ds) RNA Sequencing (FLDS), to screen RNA viruses from 155 fungal isolates, which allowed us to obtain complete viral genomes in a homology-independent manner. We created a high-quality catalog of 19 RNA viruses (12 viral species) that infect Aspergillus isolates. Among them, nine viruses were not detectable by the conventional methodology involving agarose gel electrophoresis of dsRNA, a hallmark of RNA virus infections. Segmented genome structures were determined in 42 per cent of the viruses. Some RNA viruses had novel genome architectures; one contained a dual methyltransferase domain and another had a separated RNA-dependent RNA polymerase (RdRp) gene. A virus from a different fungal taxon (Pyricularia) had an RdRp sequence that was separated on different segments, suggesting that a divided RdRp is widely present among fungal viruses, despite the belief that all RNA viruses encode RdRp as a single gene. These findings illustrate the previously hidden diversity and evolution of RNA viruses, and prompt reconsideration of the structural plasticity of RdRp.

Genome Organization of a New Double-Stranded RNA LA Helper Virus From Wine Torulaspora delbrueckii Killer Yeast as Compared With Its Saccharomyces Counterparts

Wine killer yeasts such as killer strains of Torulaspora delbrueckii and Saccharomyces cerevisiae contain helper large-size (4.6 kb) dsRNA viruses (V-LA) required for the stable maintenance and replication of killer medium-size dsRNA viruses (V-M) which bear the genes that encode for the killer toxin. The genome of the new V-LA dsRNA from the T. delbrueckii Kbarr1 killer yeast (TdV-LAbarr1) was characterized by high-throughput sequencing (HTS). The canonical genome of TdV-LAbarr1 shares a high sequence identity and similar genome organization with its Saccharomyces counterparts. It contains all the known conserved motifs predicted to be necessary for virus translation, packaging, and replication. Similarly, the Gag-Pol amino-acid sequence of this virus contains all the features required for cap-snatching and RNA polymerase activity, as well as the expected regional variables previously found in other LA viruses. Sequence comparison showed that two main clusters (99.2–100% and 96.3–98.8% identity) include most LA viruses from Saccharomyces, with TdV-LAbarr1 being the most distant from all these viruses (61.5–62.5% identity). Viral co-evolution and cross transmission between different yeast species are discussed based on this sequence comparison. Additional 5′ and 3′ sequences were found in the TdV-LAbarr1 genome as well as in some newly sequenced V-LA genomes from S. cerevisiae. A stretch involving the 5′ extra sequence of TdV-LAbarr1 is identical to a homologous stretch close to the 5′ end of the canonical sequence of the same virus (self-identity). Our modeling suggests that these stretches can form single-strand stem loops, whose unpaired nucleotides could anneal to create an intramolecular kissing complex. Similar stem loops are also found in the 3′ extra sequence of the same virus as well as in the extra sequences of some LA viruses from S. cerevisiae. A possible origin of these extra sequences as well as their function in obviating ssRNA degradation and allowing RNA transcription and replication are discussed.

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.

Viral ecogenomics across the Porifera

BACKGROUND

Viruses directly affect the most important biological processes in the ocean via their regulation of prokaryotic and eukaryotic populations. Marine sponges form stable symbiotic partnerships with a wide diversity of microorganisms and this high symbiont complexity makes them an ideal model for studying viral ecology. Here, we used morphological and molecular approaches to illuminate the diversity and function of viruses inhabiting nine sponge species from the Great Barrier Reef and seven from the Red Sea.

RESULTS

Viromic sequencing revealed host-specific and site-specific patterns in the viral assemblages, with all sponge species dominated by the bacteriophage order Caudovirales but also containing variable representation from the nucleocytoplasmic large DNA virus families Mimiviridae, Marseilleviridae, Phycodnaviridae, Ascoviridae, Iridoviridae, Asfarviridae and Poxviridae. Whilst core viral functions related to replication, infection and structure were largely consistent across the sponge viromes, functional profiles varied significantly between species and sites largely due to differential representation of putative auxiliary metabolic genes (AMGs) and accessory genes, including those associated with herbicide resistance, heavy metal resistance and nylon degradation. Furthermore, putative AMGs varied with the composition and abundance of the sponge-associated microbiome. For instance, genes associated with antimicrobial activity were enriched in low microbial abundance sponges, genes associated with nitrogen metabolism were enriched in high microbial abundance sponges and genes related to cellulose biosynthesis were enriched in species that host photosynthetic symbionts.

CONCLUSIONS

Our results highlight the diverse functional roles that viruses can play in marine sponges and are consistent with our current understanding of sponge ecology. Differential representation of putative viral AMGs and accessory genes across sponge species illustrate the diverse suite of beneficial roles viruses can play in the functional ecology of these complex reef holobionts. Video Abstract.