FLDS: A Comprehensive dsRNA Sequencing Method for Intracellular RNA Virus Surveillance

Advancing RNA phage biology through meta-omics

Bacteriophages with RNA genomes are among the simplest biological entities on Earth. Since their discovery in the 1960s, they have been used as important models to understand the principal processes of life, including translation and the genetic code. While RNA phages were generally thought of as rare oddities in nature, meta-omics methods are rapidly changing this simplistic view by studying diverse biomes with unprecedented resolution. Metatranscriptomics dramatically expanded the number of known RNA phages from tens to tens of thousands, revealed their widespread abundance, and discovered several new families of potential RNA phages with largely unknown hosts, biology, and environmental impact. At the same time, (meta)genomic analyses of bacterial hosts are discovering an arsenal of defense systems bacteria employ to protect themselves from predation, whose functions in immunity against RNA phages we are only beginning to understand. Here, I review how meta-omics approaches are advancing the field of RNA phage biology with a focus on the discovery of new RNA phages and how bacteria might fight them.

Molecular characterization of a polymycovirus in Leptosphaeria biglobosa

Leptosphaeria biglobosa is a phytopathogenic ascomycete of Brassica napus that causes phoma stem canker/blackleg. A new double-stranded RNA (dsRNA) mycovirus from this fungus has been fully characterized. The virus genome has five dsRNA segments, ranging in length from 1,180 bp to 2,402 bp. Each dsRNA has a single open reading frame flanked by 5' and 3' untranslated regions. The proteins encoded by dsRNAs 1 and 3, an RNA-dependent RNA polymerase (RdRP) and a methyltransferase, respectively, have significant similarity to those of Plasmopara viticola lesion associated polymycovirus 1. The proline-alanine-serine-rich protein encoded by dsRNA 5 is similar to that of Erysiphe necator associated polymycovirus 1. The proteins encoded by dsRNAs 2 and 4 have significant similarity to those of a mycovirus identified in Alternaria sp. FA0703. Phylogenetic analysis based on RdRP sequences showed that this virus clusters with members of the family Polymycoviridae. Based on these observations, this virus, which we have named "Leptosphaeria biglobosa polymycovirus 1", should be classified as a member of the family Polymycoviridae. This is the first report of a polymycovirus in L. biglobosa.

Integrated analysis of protein sequence and structure redefines viral diversity and the taxonomy of the Flaviviridae

The Flaviviridae are a family of non-segmented positive-sense enveloped RNA viruses containing significant pathogens including hepatitis C virus and yellow fever virus. Recent large-scale metagenomic surveys have identified many diverse RNA viruses related to classical orthoflaviviruses and pestiviruses but quite different genome lengths and configurations, and with a hugely expanded host range that spans multiple animal phyla, including molluscs, cnidarians and stramenopiles,, and plants. Grouping of RNA-directed RNA polymerase (RdRP) hallmark gene sequences of flavivirus and 'flavi-like' viruses into four divergent clades and multiple lineages within them was congruent with helicase gene phylogeny, PPHMM profile comparisons, and comparison of RdRP protein structure predicted by AlphFold2. These results support their classification into the established order, Amarillovirales, in three families (Flaviviridae, Pestiviridae, and Hepaciviridae), and 14 genera. This taxonomic framework informed by RdRP hallmark gene evolutionary relationships provides a stable reference from which major genome re-organisational events can be understood.

Mycoviruses from Aspergillus fungi involved in fermentation of dried bonito

Fungi are exploited for fermentation of foods such as cheese, Japanese sake, and soy sauce. However, the diversity of viruses that infect fungi involved in food fermentation is poorly understood. Fermented dried bonito ("katsuobushi") is one of the most important processed marine products in Japan. Fungi involved in katsuobushi fermentation are called katsuobushi molds, and Aspergillus spp. have been reported to be dominant on the surface of katsuobushi during fermentation. Because various mycoviruses have been found in members of the genus Aspergillus, we hypothesized that katsuobushi molds are also infected with mycoviruses. Here, we describe seven novel mycoviruses belonging to six families (Chrysoviridae, Fusariviridae, Mitoviridae, Partitiviridae, Polymycoviridae, and Pseudototiviridae) from isolated katsuobushi molds (Aspergillus chevalieri and A. sulphureus) detected by fragmented and primer-ligated double-stranded RNA sequencing. Aspergillus chevalieri fusarivirus 1 has a unique bi-segmented genome, whereas other known fusariviruses have a single genomic segment. Phenotypic comparison between the parental A. chevalieri strain infected with Aspergillus chevalieri polymycovirus 1 (AchPmV1) and isogenic AchPmV1-free isolates indicated that AchPmV1 inhibits the early growth of the host. This study reveals the diversity of mycoviruses that infect katsuobushi molds, and provides insight into the effect of mycoviruses on fungi involved in fermentation.

Metagenomic Analyses of Water Samples of Two Urban Freshwaters in Berlin, Germany, Reveal New Highly Diverse Invertebrate Viruses

In an attempt to explore the RNA viromes of two German rivers, we searched the virus particle contents of one 50 L water sample each from the Teltow Canal and the Havel River for viruses assumed to infect invertebrates. More than 330 complete and partial virus genomes up to a length of 37 kb were identified, with noda-like and reo-like viruses being most abundant, followed by bunya-like and birna-like viruses. Viruses related to the Permutotetraviridae, Nidovirales, Flaviviridae, Rhabdoviridae and Chuviridae as well as the unclassified Jῑngmén virus and Negev virus groups were also present. The results indicate a broad extent of recombinant virus genomes, supporting the concept of the modularity of eukaryotic viruses. For example, novel combinations of genes encoding replicase and structural proteins with a jellyroll fold have been observed. Less than 35 viruses could be assigned to existing virus genera. These are (i) an avian deltacoronavirus which was represented by only one short contig, albeit with 98% similarity, (ii) a seadornavirus and a rotavirus, and (iii) some 30 nodaviruses. All remaining viruses are novel and too diverse for accommodation in existing genera. Many of the virus genomes exhibit ORFans encoding hypothetical proteins of up to 2000 amino acids without conserved protein domains.

A single-stranded based library preparation method for virome characterization

BACKGROUND

The gut virome is an integral component of the gut microbiome, playing a crucial role in maintaining gut health. However, accurately depicting the entire gut virome is challenging due to the inherent diversity of genome types (dsDNA, ssDNA, dsRNA, and ssRNA) and topologies (linear, circular, or fragments), with subsequently biases associated with current sequencing library preparation methods. To overcome these problems and improve reproducibility and comparability across studies, universal or standardized virome sequencing library construction methods are highly needed in the gut virome study.

RESULTS

We repurposed the ligation-based single-stranded library (SSLR) preparation method for virome studies. We demonstrate that the SSLR method exhibits exceptional efficiency in quantifying viral DNA genomes (both dsDNA and ssDNA) and outperforms existing double-stranded (Nextera) and single-stranded (xGen, MDA + Nextera) library preparation approaches in terms of minimal amplification bias, evenness of coverage, and integrity of assembling viral genomes. The SSLR method can be utilized for the simultaneous library preparation of both DNA and RNA viral genomes. Furthermore, the SSLR method showed its ability to capture highly modified phage genomes, which were often lost using other library preparation approaches.

CONCLUSION

We introduce and improve a fast, simple, and efficient ligation-based single-stranded DNA library preparation for gut virome study. This method is compatible with Illumina sequencing platforms and only requires ligation reagents within 3-h library preparation, which is similar or even better than the advanced library preparation method (xGen). We hope this method can be further optimized, validated, and widely used to make gut virome study more comparable and reproducible. Video Abstract.

Mapping glycoprotein structure reveals Flaviviridae evolutionary history

Viral glycoproteins drive membrane fusion in enveloped viruses and determine host range, tissue tropism and pathogenesis1. Despite their importance, there is a fragmentary understanding of glycoproteins within the Flaviviridae2, a large virus family that include pathogens such as hepatitis C, dengue and Zika viruses, and numerous other human, animal and emergent viruses. For many flaviviruses the glycoproteins have not yet been identified, for others, such as the hepaciviruses, the molecular mechanisms of membrane fusion remain uncharacterized3. Here we combine phylogenetic analyses with protein structure prediction to survey glycoproteins across the entire Flaviviridae. We find class II fusion systems, homologous to the Orthoflavivirus E glycoprotein in most species, including highly divergent jingmenviruses and large genome flaviviruses. However, the E1E2 glycoproteins of the hepaciviruses, pegiviruses and pestiviruses are structurally distinct, may represent a novel class of fusion mechanism, and are strictly associated with infection of vertebrate hosts. By mapping glycoprotein distribution onto the underlying phylogeny, we reveal a complex evolutionary history marked by the capture of bacterial genes and potentially inter-genus recombination. These insights, made possible through protein structure prediction, refine our understanding of viral fusion mechanisms and reveal the events that have shaped the diverse virology and ecology of the Flaviviridae.

New lineages of RNA viruses from clinical isolates of Rhizopus microsporus revealed by fragmented and primer-ligated dsRNA sequencing (FLDS) analysis

Rhizopus microsporus is a species in the order Mucorales that is known to cause mucormycosis, but it is poorly understood as a host of viruses. Here, we examined 25 clinical strains of R. microsporus for viral infection with a conventional double-stranded RNA (dsRNA) assay using agarose gel electrophoresis (AGE) and the recently established fragmented and primer-ligated dsRNA sequencing (FLDS) protocol. By AGE, five virus-infected strains were detected. Then, full-length genomic sequences of 12 novel RNA viruses were revealed by FLDS, which were related to the families Mitoviridae, Narnaviridae, and Endornaviridae, ill-defined groups of single-stranded RNA (ssRNA) viruses with similarity to the established families Virgaviridae and Phasmaviridae, and the proposed family "Ambiguiviridae." All the characterized viruses, except a potential phasmavirid with a negative-sense RNA genome, had positive-sense RNA genomes. One virus belonged to a previously established species within the family Mitoviridae, whereas the other 11 viruses represented new species or even new genera. These results show that the fungal pathogen R. microsporus harbors diverse RNA viruses and extend our understanding of the diversity of RNA viruses in the fungal order Mucorales, division Mucoromycota. Identifying RNA viruses from clinical isolates of R. microsporus may expand the repertoire of natural therapeutic agents for mucormycosis in the future.IMPORTANCEThe diversity of mycoviruses in fungal hosts in the division Mucoromycota has been underestimated, mainly within the species Rhizopus microsporus. Only five positive-sense RNA genomes had previously been discovered in this species. Because current sequencing methods poorly complete the termini of genomes, we used fragmented and primer-ligated double-stranded RNA sequencing to acquire the full-length genomes. Eleven novel mycoviruses were detected in this study, including the first negative-sense RNA genome reported in R. microsporus. Our findings extend the understanding of the viral diversity in clinical strains of Mucoromycota, may provide insights into the pathogenesis and ecology of this fungus, and may offer therapeutic options.

Metatranscriptomic analysis reveals dissimilarity in viral community activity between an ice-free and ice-covered winter in Lake Erie

Winter is a relatively under-studied season in freshwater ecology. The paucity of wintertime surveys has led to a lack of knowledge regarding microbial community activity during the winter in Lake Erie, a North American Great Lake. Viruses shape microbial communities and regulate biogeochemical cycles by acting as top-down controls, yet very few efforts have been made to examine active virus populations during the winter in Lake Erie. Furthermore, climate change-driven declines in seasonal ice cover have been shown to influence microbial community structure, but no studies have compared viral community activity between different ice cover conditions. We surveyed surface water metatranscriptomes for viral hallmark genes as a proxy for active virus populations and compared activity metrics between ice-covered and ice-free conditions from two sampled winters. Transcriptionally active viral communities were detected in both winters, spanning diverse phylogenetic clades of putative bacteriophage (Caudoviricetes), giant viruses (Nucleocytoviricota, or NCLDV), and RNA viruses (Orthornavirae). However, viral community activity metrics revealed pronounced differences between the ice-covered and ice-free winters. Viral community composition was distinct between winters and viral hallmark gene richness was reduced in the ice-covered relative to the ice-free conditions. In addition, the observed differences in viral communities correlated with microbial community activity metrics. Overall, these findings contribute to our understanding of the viral populations that are active during the winter in Lake Erie and suggest that viral community activity may be associated with ice cover extent.IMPORTANCEAs seasonal ice cover is projected to become increasingly rare on large temperate lakes, there is a need to understand how microbial communities might respond to changing ice conditions. Although it is widely recognized that viruses impact microbial community structure and function, there is little known regarding wintertime viral activity or the relationship between viral activity and ice cover extent. Our metatranscriptomic analyses indicated that viruses were transcriptionally active in the winter surface waters of Lake Erie. These findings also expanded the known diversity of viral lineages in the Great Lakes. Notably, viral community activity metrics were significantly different between the two sampled winters. The pronounced differences we observed in active viral communities between the ice-covered and ice-free samples merit further research regarding how viral communities will function in future, potentially ice-free, freshwater systems.

Molecular and biological characterization of a novel partitivirus from Talaromyces pinophilus

Talaromyces spp. have a worldwide distribution, are ecologically diverse and have been isolated from numerous different substrates. Talaromyces spp. are considered biotechnologically important due to their ability to produce a range of enzymes and pigments. Talaromyces pinophilus, belonging to genus Talaromyces and family Trichocomaceae, is known for producing several important bioactive metabolites. Here we report the isolation and characterisation of a partitivirus from T. pinophilus which we have nominated Talaromyces pinophilus partitivirus-1 (TpPV-1). TpPV-1 possesses a genome consisting of three double stranded (ds) RNA segments i.e., dsRNAs1-3, 1824 bp, 1638 bp and 1451 bp respectively, which are encapsidated in icosahedral particles 35 nm in diameter. Both dsRNA1 and dsRNA2 contain a single open reading frame (ORF) encoding respectively a 572 amino acid (aa) protein of 65 kDa and a 504 aa protein of 50 kDa. The third segment (dsRNA3) is potentially a satellite RNA. Phylogenetic analysis revealed that the TpPV-1 belongs to the family Partitiviridae in the proposed genus Zetapartitivirus. TpPV-1 infection decreases the mycelial growth rate of the host fungus and alters pigmentation as indicated by time course experiments performed on a range of different solid media comparing virus-infected and virus-free isogenic lines. This is the first report of mycovirus infection in T. pinophilus and may provide insights into understanding the effect of the mycovirus on the production of enzymes and pigments by the host fungus.

Genomic and phylogenetic features of the Picobirnaviridae suggest microbial rather than animal hosts

The RNA virus family Picobirnaviridae has traditionally been associated with the gastrointestinal systems of terrestrial mammals and birds, with the majority of viruses detected in animal stool samples. Metatranscriptomic studies of vertebrates, invertebrates, microbial communities, and environmental samples have resulted in an enormous expansion of the genomic and phylogenetic diversity of this family. Yet picobirnaviruses remain poorly classified, with only one genus and three species formally ratified by the International Committee of Virus Taxonomy. Additionally, an inability to culture picobirnaviruses in a laboratory setting or isolate them in animal tissue samples, combined with the presence of bacterial genetic motifs in their genomes, suggests that these viruses may represent RNA bacteriophage rather than being associated with animal infection. Utilising a data set of 2,286 picobirnaviruses sourced from mammals, birds, reptiles, fish, invertebrates, microbial communities, and environmental samples, we identified seven consistent phylogenetic clusters likely representing Picobirnavirus genera that we tentatively name 'Alpha-', 'Beta-', 'Gamma-', 'Delta-', 'Epsilon-', 'Zeta-', and 'Etapicobirnavirus'. A statistical analysis of topological congruence between virus-host phylogenies revealed more frequent cross-species transmission than any other RNA virus family. In addition, bacterial ribosomal binding site motifs were more enriched in Picobirnavirus genomes than in the two groups of established RNA bacteriophage-the Leviviricetes and Cystoviridae. Overall, our findings support the hypothesis that the Picobirnaviridae have bacterial hosts and provide a lower-level taxonomic classification for this highly diverse and ubiquitous family of RNA viruses.

DsRNA sequencing revealed a previously missed terminal sequence of a +ssRNA virus that infects dinoflagellate Heterocapsa circularisquama

Heterocapsa circularisquama RNA virus (HcRNAV) is the only dinoflagellate-infecting RNA virus cultured. However, only two strains of HcRNAV have been registered with complete genome sequences (strains 34 and 109 for UA and CY types, respectively). To extend the genomic information of HcRNAV, we performed full-genome sequencing of an unsequenced strain of HcRNAV (strain A8) using the fragmented and primer-ligated double-stranded RNA (dsRNA) sequencing (FLDS) method. The complete genome of HcRNAV A8 with 4457 nucleotides (nt) was successfully determined, and sequence alignment of the major capsid protein gene suggested that A8 was a UA-type strain, consistent with its intraspecific host specificity. The complete sequence was found to be 80 nt longer at the 5' terminus than the registered sequences of HcRNAV strains (34 and 109), suggesting that FLDS is more reliable for determining the terminal sequence than conventional methods (5' Rapid Amplification of cDNA End). Our study contributes to a better understanding of dinoflagellate-infecting viruses with limited sequence data.

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.

Greetings from virologists to mycologists: A review outlining viruses that live in fungi

Viruses are genetic elements that parasitize self-replicating cells. Therefore, organisms parasitized by viruses are not limited to animals and plants but also include microorganisms. Among these, viruses that parasitize fungi are known as mycoviruses. Mycoviruses with an RNA genome persistently replicate inside fungal cells and coevolve with their host cells, similar to a cellular organelle. Within host cells, mycoviruses can modulate various fungal characteristics and activities, including pathogenicity and the production of enzymes and secondary metabolites. In this review, we provide an overview of the mycovirus research field as introduction to fungal researchers. Recognition of all genetic elements in fungi aids towards better understanding and control of fungi, and makes fungi a significant model system for studying microorganisms containing multiple genetic elements.

New clades of viruses infecting the obligatory biotroph Bremia lactucae representing distinct evolutionary trajectory for viruses infecting oomycetes

Recent advances in high throughput sequencing (HTS) approaches allowed a broad exploration of viromes from different fungal hosts, unveiling a great diversity of mycoviruses with interesting evolutionary features. The word mycovirus historically applies also to viruses infecting oomycetes but most studies are on viruses infecting fungi, with less mycoviruses found and characterized in oomycetes, particularly in the obligatory biotrophs. We, here, describe the first virome associated to Bremia lactucae, the causal agent of lettuce downy mildew, which is an important biotrophic pathogen for lettuce production and a model system for the molecular aspects of the plant-oomycetes interactions. Among the identified viruses, we could detect (1) two new negative sense ssRNA viruses related to the yueviruses, (2) the first example of permuted RdRp in a virus infecting fungi/oomycetes, (3) a new group of bipartite dsRNA viruses showing evidence of recent bi-segmentation and concomitantly, a possible duplication event bringing a bipartite genome to tripartite, (4) a first representative of a clade of viruses with evidence of recombination between distantly related viruses, (5) a new open reading frame (ORF)an virus encoding for an RdRp with low homology to known RNA viruses, and (6) a new virus, belonging to riboviria but not conserved enough to provide a conclusive phylogenetic placement that shows evidence of a recombination event between a kitrinoviricota-like and a pisuviricota-like sequence. The results obtained show a great diversity of viruses and evolutionary mechanisms previously unreported for oomycetes-infecting viruses, supporting the existence of a large diversity of oomycetes-specific viral clades ancestral of many fungal and insect virus clades.

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.

Current research and future directions for realizing the ideal One-Health approach: A summary of key-informant interviews in Japan and a literature review

The COVID-19 pandemic has highlighted the importance of the One Health (OH) approach, which considers the health of humans, animals, and the environment in preventing future pandemics. A wide range of sustainable interdisciplinary collaborations are required to truly fulfill the purpose of the OH approach. It is well-recognized, however, that such collaborations are challenging. In this study, we undertook key-informant interviews with a panel of stakeholders from Japan to identify their perceived needs and challenges related to OH research. This panel included scientists, government officials, journalists, and industry stakeholders. By combining a thematic analysis of these interviews and a literature review, we summarized two key themes pertinent to the effective implementation of OH research: types of required research and systems to support that research. As a technological issue, interviewees suggested the importance of research and development of methodologies that can promote the integration and collaboration of research fields that are currently fragmented. An example of such a methodology would allow researchers to obtain high-resolution metadata (e.g. ecological and wildlife data) with high throughput and then maximize the use of the obtained metadata in research, such as in environmental DNA analysis, database construction, or the use of computational algorithms to find novel viral genomes. In terms of systems surrounding OH research, some interviewees stressed the importance of creating a sustainable research system, such as one that has continuous budget support and allows researchers to pursue their academic careers and interests. These perceptions and challenges held by Japanese stakeholders may be common to others around the world. We hope this review will encourage more researchers and others to work together to create a resilient society against future pandemics.

Association of virome dynamics with mosquito species and environmental factors

BACKGROUND

The pathogenic viruses transmitted by mosquitoes cause a variety of animal and human diseases and public health concerns. Virome surveillance is important for the discovery, and control of mosquito-borne pathogenic viruses, as well as early warning systems. Virome composition in mosquitoes is affected by mosquito species, food source, and geographic region. However, the complex associations of virome composition remain largely unknown.

RESULTS

Here, we profiled the high-depth RNA viromes of 15 species of field-caught adult mosquitoes, especially from Culex, Aedes, Anopheles, and Armigeres in Hainan Island from 2018 to 2020. We detected 57 known and 39 novel viruses belonging to 15 families. We established the associations of the RNA viruses with mosquito species and their foods, indicating the importance of feeding acquisition of RNA viruses in determining virome composition. A large fraction of RNA viruses were persistent in the same mosquito species across the 3 years and different locations, showing the species-specific stability of viromes in Hainan Island. In contrast, the virome compositions of single mosquito species in different geographic regions worldwide are visibly distinct. This is consistent with the differences in food sources of mosquitoes distributed broadly across continents.

CONCLUSIONS

Thus, species-specific viromes in a relatively small area are limited by viral interspecific competition and food sources, whereas the viromes of mosquito species in large geographic regions may be governed by ecological interactions between mosquitoes and local environmental factors. Video Abstract.

Exploring the Mycovirus Universe: Identification, Diversity, and Biotechnological Applications

Viruses that infect fungi are known as mycoviruses and are characterized by the lack of an extracellular phase. In recent years, the advances on nucleic acids sequencing technologies have led to a considerable increase in the number of fungi-infecting viral species described in the literature, with a special interest in assessing potential applications as fungal biocontrol agents. In the present study, we performed a comprehensive review using Scopus, Web of Science, and PubMed databases to mine mycoviruses data to explore their molecular features and their use in biotechnology. Our results showed the existence of 267 mycovirus species, of which 189 are recognized by the International Committee on Taxonomy of Viruses (ICTV). The majority of the mycoviruses identified have a dsRNA genome (38.6%), whereas the Botourmiaviridae (ssRNA+) alone represents 14% of all mycoviruses diversity. Regarding fungal hosts, members from the Sclerotinicaeae appeared as the most common species described to be infected by mycoviruses, with 16 different viral families identified so far. It is noteworthy that such results are directly associated with the high number of studies and strategies used to investigate the presence of viruses in members of the Sclerotinicaeae family. The knowledge about replication strategy and possible impact on fungi biology is available for only a small fraction of the mycoviruses studied, which is the main limitation for considering these elements potential targets for biotechnological applications. Altogether, our investigation allowed us to summarize the general characteristics of mycoviruses and their hosts, the consequences, and the implications of this knowledge on mycovirus–fungi interactions, providing an important source of information for future studies.

The First Identification of a Narnavirus in Bigyra, a Marine Protist

Current information on the diversity and evolution of eukaryotic RNA viruses is biased towards host lineages, such as animals, plants, and fungi. Although protists represent the majority of eukaryotic diversity, our understanding of the protist RNA virosphere is still limited. To reveal untapped RNA viral diversity, we screened RNA viruses from 30 marine protist isolates and identified a novel RNA virus named Haloplacidia narnavirus 1 (HpNV1). A phylogenetic ana-lysis revealed that HpNV1 is a new member of the family Narnaviridae. The present study filled a gap in the distribution of narnaviruses and implies their wide distribution in Stramenopiles.

Mycoviruses as a part of the global virome: Diversity, evolutionary links and lifestyle

Knowledge of mycovirus diversity, evolution, horizontal gene transfer and shared ancestry with viruses infecting distantly related hosts, such as plants and arthropods, has increased vastly during the last few years due to advances in the high throughput sequencing methodologies. This also has enabled the discovery of novel mycoviruses with previously unknown genome types, mainly new positive and negative single-stranded RNA mycoviruses ((+) ssRNA and (-) ssRNA) and single-stranded DNA mycoviruses (ssDNA), and has increased our knowledge of double-stranded RNA mycoviruses (dsRNA), which in the past were thought to be the most common viruses infecting fungi. Fungi and oomycetes (Stramenopila) share similar lifestyles and also have similar viromes. Hypothesis about the origin and cross-kingdom transmission events of viruses have been raised and are supported by phylogenetic analysis and by the discovery of natural exchange of viruses between different hosts during virus-fungus coinfection in planta. In this review we make a compilation of the current information on the genome organization, diversity and taxonomy of mycoviruses, discussing their possible origins. Our focus is in recent findings suggesting the expansion of the host range of many viral taxa previously considered to be exclusively fungal, but we also address factors affecting virus transmissibility and coexistence in single fungal or oomycete isolates, as well as the development of synthetic mycoviruses and their use in investigating mycovirus replication cycles and pathogenicity.

Discovery, Genomic Sequence Characterization and Phylogenetic Analysis of Novel RNA Viruses in the Turfgrass Pathogenic Colletotrichum spp. in Japan

Turfgrass used in various areas of the golf course has been found to present anthracnose disease, which is caused by Colletotrichum spp. To obtain potential biological agents, we identified four novel RNA viruses and obtained full-length viral genomes from turfgrass pathogenic Colletotrichum spp. in Japan. We characterized two novel dsRNA partitiviruses: Colletotrichum associated partitivirus 1 (CaPV1) and Colletotrichum associated partitivirus 2 (CaPV2), as well as two negative single-stranded (ss) RNA viruses: Colletotrichum associated negative-stranded RNA virus 1 (CaNSRV1) and Colletotrichum associated negative-stranded RNA virus 2 (CaNSRV2). Using specific RT-PCR assays, we confirmed the presence of CaPV1, CaPV2 and CaNSRV1 in dsRNAs from original and sub-isolates of Colletotrichum sp. MBCT-264, as well as CaNSRV2 in dsRNAs from original and sub-isolates of Colletotrichum sp. MBCT-288. This is the first time mycoviruses have been discovered in turfgrass pathogenic Colletotrichum spp. in Japan. CaPV1 and CaPV2 are new members of the newly proposed genus "Zetapartitivirus" and genus Alphapartitivirus, respectively, in the family Partitiviridae, according to genomic characterization and phylogenetic analysis. Negative sense ssRNA viruses CaNSRV1 and CaNSRV2, on the other hand, are new members of the family Phenuiviridae and the proposed family "Mycoaspirividae", respectively. These findings reveal previously unknown RNA virus diversity and evolution in turfgrass pathogenic Colletotrichum spp.

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.

Molecular Characterization of Novel Mycoviruses in Seven Umbelopsis Strains

The presence of viruses is less explored in Mucoromycota as compared to other fungal groups such as Ascomycota and Basidiomycota. Recently, more and more mycoviruses are identified from the early-diverging lineages of fungi. We have determined the genome of 11 novel dsRNA viruses in seven different Umbelopsis strains with next-generation sequencing (NGS). The identified viruses were named Umbelopsis ramanniana virus 5 (UrV5), 6a (UrV6a); 6b (UrV6b); 7 (UrV7); 8a (UrV8a); 8b (UrV8b); Umbelopsis gibberispora virus 1 (UgV1); 2 (UgV2) and Umbelopsis dimorpha virus 1a (UdV1a), 1b (UdV1b) and 2 (UdV2). All the newly identified viruses contain two open reading frames (ORFs), which putatively encode the coat protein (CP) and the RNA-dependent RNA polymerase (RdRp), respectively. Based on the phylogeny inferred from the RdRp sequences, eight viruses (UrV7, UrV8a, UrV8b, UgV1, UgV2, UdV1a, UdV1b and UdV2) belong to the genus Totivirus, while UrV5, UrV6a and UrV6b are placed into a yet unclassified but well-defined Totiviridae-related group. In UrV5, UgV1, UgV2, UrV8b, UdV1a, UdV2 and UdV1b, ORF2 is predicted to be translated as a fusion protein via a rare +1 (or -2) ribosomal frameshift, which is not characteristic to most members of the Totivirus genus. Virus particles 31 to 32 nm in diameter could be detected in the examined fungal strains by transmission electron microscopy. Through the identification and characterization of new viruses of Mucoromycota fungi, we can gain insight into the diversity of mycoviruses, as well as into their phylogeny and genome organization.

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.

Mycovirus Diversity and Evolution Revealed/Inferred from Recent Studies

High-throughput virome analyses with various fungi, from cultured or uncultured sources, have led to the discovery of diverse viruses with unique genome structures and even neo-lifestyles. Examples in the former category include splipalmiviruses and ambiviruses. Splipalmiviruses, related to yeast narnaviruses, have multiple positive-sense (+) single-stranded (ss) RNA genomic segments that separately encode the RNA-dependent RNA polymerase motifs, the hallmark of RNA viruses (members of the kingdom Orthornavirae). Ambiviruses appear to have an undivided ssRNA genome of 3∼5 kb with two large open reading frames (ORFs) separated by intergenic regions. Another narna-like virus group has two fully overlapping ORFs on both strands of a genomic segment that span more than 90% of the genome size. New virus lifestyles exhibited by mycoviruses include the yado-kari/yado-nushi nature characterized by the partnership between the (+)ssRNA yadokarivirus and an unrelated dsRNA virus (donor of the capsid for the former) and the hadaka nature of capsidless 10-11 segmented (+)ssRNA accessible by RNase in infected mycelial homogenates. Furthermore, dsRNA polymycoviruses with phylogenetic affinity to (+)ssRNA animal caliciviruses have been shown to be infectious as dsRNA-protein complexes or deproteinized naked dsRNA. Many previous phylogenetic gaps have been filled by recently discovered fungal and other viruses, which haveprovided interesting evolutionary insights. Phylogenetic analyses and the discovery of natural and experimental cross-kingdom infections suggest that horizontal virus transfer may have occurred and continue to occur between fungi and other kingdoms.

Freshwater macrophytes harbor viruses representing all five major phyla of the RNA viral kingdom Orthornavirae

Research on aquatic plant viruses is lagging behind that of their terrestrial counterparts. To address this knowledge gap, here we identified viruses associated with freshwater macrophytes, a taxonomically diverse group of aquatic phototrophs that are visible with the naked eye. We surveyed pooled macrophyte samples collected at four spring sites in Florida, USA through next generation sequencing of RNA extracted from purified viral particles. Sequencing efforts resulted in the detection of 156 freshwater macrophyte associated (FMA) viral contigs, 37 of which approximate complete genomes or segments. FMA viral contigs represent putative members from all five major phyla of the RNA viral kingdom Orthornavirae. Similar to viral types found in land plants, viral sequences identified in macrophytes were dominated by positive-sense RNA viruses. Over half of the FMA viral contigs were most similar to viruses reported from diverse hosts in aquatic environments, including phototrophs, invertebrates, and fungi. The detection of FMA viruses from orders dominated by plant viruses, namely Patatavirales and Tymovirales, indicate that members of these orders may thrive in aquatic hosts. PCR assays confirmed the presence of putative FMA plant viruses in asymptomatic vascular plants, indicating that viruses with persistent lifestyles are widespread in macrophytes. The detection of potato virus Y and oat blue dwarf virus in submerged macrophytes suggests that terrestrial plant viruses infect underwater plants and highlights a potential terrestrial-freshwater plant virus continuum. Defining the virome of unexplored macrophytes will improve our understanding of virus evolution in terrestrial and aquatic primary producers and reveal the potential ecological impacts of viral infection in macrophytes.

RNA viromes from terrestrial sites across China expand environmental viral diversity

Environmental RNA viruses are ubiquitous and diverse, and probably have important ecological and biogeochemical impacts. Understanding the global diversity of RNA viruses is limited by sampling biases, dependence on cell culture and PCR for virus discovery, and a focus on viruses pathogenic to humans or economically important animals and plants. To address this knowledge gap, we generated metatranscriptomic sequence data from 32 diverse environments in 16 provinces and regions of China. We identified 6,624 putatively novel virus operational taxonomic units from soil, sediment and faecal samples, greatly expanding known diversity of the RNA virosphere. These newly identified viruses included positive-sense, negative-sense and double-strand RNA viruses from at least 62 families. Sediments and animal faeces were rich sources of viruses. Virome compositions were affected by local environmental factors, including organic content and eukaryote species abundance. Notably, environmental factors had a greater impact on the abundance and diversity of plant, fungal and bacterial viruses than of animal viromes. Our data confirm that RNA viruses are an integral part of both terrestrial and aquatic ecosystems.

Resolving deep evolutionary relationships within the RNA virus phylum Lenarviricota

The RNA virus phylum Lenarviricota is composed of the fungi-associated families Narnaviridae and Mitoviridae, the RNA bacteriophage Leviviridae, and the plant and fungi-associated Botourmiaviridae. Members of the Lenarviricota are abundant in most environments and boast remarkable phylogenetic and genomic diversity. As this phylum includes both RNA bacteriophage and fungi- and plant-associated species, the Lenarviricota likely mark a major evolutionary transition between those RNA viruses associated with prokaryotes and eukaryotes. Despite the remarkable expansion of this phylum following metagenomic studies, the phylogenetic relationships among the families within the Lenarviricota remain uncertain. Utilising a large data set of relevant viral sequences, we performed phylogenetic and genomic analyses to resolve the complex evolutionary history within this phylum and identify patterns in the evolution of virus genome organisation. Despite limitations reflecting very high levels of sequence diversity, our phylogenetic analyses suggest that the Leviviridae comprise the basal lineage within the Lenarviricota. Our phylogenetic results also support the construction of a new virus family-the Narliviridae-comprising a set of diverse and phylogenetically distinct species, including a number of uniquely encapsidated viruses. We propose a taxonomic restructuring within the Lenarviricota to better reflect the phylogenetic relationships documented here, with the Botourmiaviridae and Narliviridae combined into the order Ourlivirales, the Narnaviridae remaining in the order Wolframvirales, and these orders combined into the single class, the Amabiliviricetes. In sum, this study provides insights into the complex evolutionary relationships among the diverse families that make up the Lenarviricota.

Increasing the Sensitivity of Leishmania RNA Virus 2 (LRV2) Detection with a Modification in cDNA Synthesis

OBJECTIVE

Leishmania RNA virus was detected the first time in the New World Leishmania species. Recent studies were also showed the presence of Leishmania RNA virus 2 (LRV2) in Old Word Leishmania species including Turkish L. major and L. tropica isolates. This study aimed to increase the sensitivity of qPCR with a modification in the denaturation step of cDNA preparation protocol.

METHODS

In this study, LRV2+ three L. major, two L. tropica strains and L. major control strain (MHOM/SU/73/5-ASKH) were included. Total RNA isolation was done using different numbers of Leishmania promastigotes (10⁸, 10⁵ and 10³). Before cDNA synthesis, samples were denatured at 95 °C for 2 min, as a modification of the kit procedure. qPCR was undertaken using 0.5 mM primers (LRV F-HR/LRV R-HR) diluted in SYBR Green Master mix.

RESULTS

We observed lower Ct values in amplicons with the modified version than with the classical kit protocol for cDNA synthesis, in all of the strains used in the study. The addition of pre-denaturation step at 95 °C showed lower Ct values meaning the sensitivity increased. Different parasite dilutions showed similar results.

CONCLUSION

It is important to increase the sensitivity especially with the aim for detecting LRV in clinical samples obtained from patients probably have less number of parasites. The presence and burden of the virus can help to understand the relationship between the clinical findings and the pathogenicity of the parasite which may lead to changes in the course of treatment.

The Intrapopulation Genetic Diversity of RNA Virus May Influence the Sensitivity of Chlorine Disinfection

RNA virus populations are not clonal; rather, they comprise a mutant swarm in which sequences are slightly different from the master sequence. Genetic diversity within a population (intrapopulation genetic diversity) is critical for RNA viruses to survive under environmental stresses. Disinfection has become an important practice in the control of pathogenic viruses; however, the impact of intrapopulation genetic diversity on the sensitivity of disinfection, defined as -log10 (postdisinfected infectious titer/predisinfected titer), has not been elucidated. In this study, we serially passaged populations of rhesus rotavirus. We demonstrated that populations with reduced chlorine sensitivity emerged at random and independently of chlorine exposure. Sequencing analysis revealed that compared with sensitive populations, less-sensitive ones had higher non-synonymous genetic diversity of the outer capsid protein gene, suggesting that changes in the amino acid sequences of the outer capsid protein were the main factors influencing chlorine sensitivity. No common mutations were found among less-sensitive populations, indicating that rather than specific mutations, the diversity of the outer capsid protein itself was associated with the disinfection sensitivity and that the disinfection sensitivity changed stochastically. Simulation results suggest that the disinfection sensitivity of a genetically diverse population is destabilized if cooperative viral clusters including multiple sequences are formed. These results advocate that any prevention measures leading to low intrapopulation genetic diversity are important to prevent the spread and evolution of pathogenic RNA viruses in society.

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.

Characterization of Two Novel Toti-Like Viruses Co-infecting the Atlantic Blue Crab, Callinectes sapidus, in Its Northern Range of the United States

The advancement of high throughput sequencing has greatly facilitated the exploration of viruses that infect marine hosts. For example, a number of putative virus genomes belonging to the Totiviridae family have been described in crustacean hosts. However, there has been no characterization of the most newly discovered putative viruses beyond description of their genomes. In this study, two novel double-stranded RNA (dsRNA) virus genomes were discovered in the Atlantic blue crab (Callinectes sapidus) and further investigated. Sequencing of both virus genomes revealed that they each encode RNA dependent RNA polymerase proteins (RdRps) with similarities to toti-like viruses. The viruses were tentatively named Callinectes sapidus toti-like virus 1 (CsTLV1) and Callinectes sapidus toti-like virus 2 (CsTLV2). Both genomes have typical elements required for −1 ribosomal frameshifting, which may induce the expression of an encoded ORF1–ORF2 (gag-pol) fusion protein. Phylogenetic analyses of CsTLV1 and CsTLV2 RdRp amino acid sequences suggested that they are members of two new genera in the family Totiviridae. The CsTLV1 and CsTLV2 genomes were detected in muscle, gill, and hepatopancreas of blue crabs by real-time reverse transcription quantitative PCR (RT-qPCR). The presence of ~40 nm totivirus-like viral particles in all three tissues was verified by transmission electron microscopy, and pathology associated with CsTLV1 and CsTLV2 infections were observed by histology. PCR assays showed the prevalence and geographic range of these viruses, to be restricted to the northeast United States sites sampled. The two virus genomes co-occurred in almost all cases, with the CsTLV2 genome being found on its own in 8.5% cases, and the CsTLV1 genome not yet found on its own. To our knowledge, this is the first report of toti-like viruses in C. sapidus. The information reported here provides the knowledge and tools to investigate transmission and potential pathogenicity of these viruses.

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.

Current challenges to virus discovery by meta-transcriptomics

Meta-transcriptomic next-generation sequencing has transformed virus discovery, dramatically expanding our knowledge of the known virosphere. Nevertheless, the use of meta-transcriptomics for virus discovery faces important challenges. As this technology becomes more widely adopted, the proportion of viral sequences in public databases with incorrect (e.g. mis-assignment of host) or limited information (e.g. lacking taxonomic classification) is likely to grow, limiting their utility in bioinformatic pipelines for virus discovery. In addition, we currently lack the bioinformatic tools that can accurately identify viruses showing little or no sequence similarity to database viruses or those that represent likely reagent contaminants. Herein, we outline some of the challenges to effective meta-transcriptomic virus discovery as well as their potential solutions.

Revealing RNA virus diversity and evolution in unicellular algae transcriptomes

Remarkably little is known about the diversity and evolution of RNA viruses in unicellular eukaryotes. We screened a total of 570 transcriptomes from the Marine Microbial Eukaryote Transcriptome Sequencing Project that encompasses a wide diversity of microbial eukaryotes, including most major photosynthetic lineages (i.e. the microalgae). From this, we identified thirty new and divergent RNA virus species, occupying a range of phylogenetic positions within the overall diversity of RNA viruses. Approximately one-third of the newly described viruses comprised single-stranded positive-sense RNA viruses from the order Lenarviricota associated with fungi, plants, and protists, while another third were related to the order Ghabrivirales, including members of the protist and fungi-associated Totiviridae. Other viral species showed sequence similarity to positive-sense RNA viruses from the algae-associated Marnaviridae, the double-stranded RNA (ds-RNA) Partitiviridae, as well as tentative evidence for one negative-sense RNA virus related to the Qinviridae. Importantly, we were able to identify divergent RNA viruses from distant host taxa, revealing the ancestry of these viral families and greatly extending our knowledge of the RNA viromes of microalgal cultures. Both the limited number of viruses detected per sample and the low sequence identity to known RNA viruses imply that additional microalgal viruses exist that could not be detected at the current sequencing depth or were too divergent to be identified using sequence similarity. Together, these results highlight the need for further investigation of algal-associated RNA viruses as well as the development of new tools to identify RNA viruses that exhibit very high levels of sequence divergence.

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.

De novo Sequencing of Novel Mycoviruses From Fusarium sambucinum: An Attempt on Direct RNA Sequencing of Viral dsRNAs

An increasing number of viruses are continuously being found in a wide range of organisms, including fungi. Recent studies have revealed a wide viral diversity in microbes and a potential importance of these viruses in the natural environment. Although virus exploration has been accelerated by short-read, high-throughput sequencing (HTS), and viral de novo sequencing is still challenging because of several biological/molecular features such as micro-diversity and secondary structure of RNA genomes. This study conducted de novo sequencing of multiple double-stranded (ds) RNA (dsRNA) elements that were obtained from fungal viruses infecting two Fusarium sambucinum strains, FA1837 and FA2242, using conventional HTS and long-read direct RNA sequencing (DRS). De novo assembly of the read data from both technologies generated near-entire genomic sequence of the viruses, and the sequence homology search and phylogenetic analysis suggested that these represented novel species of the Hypoviridae, Totiviridae, and Mitoviridae families. However, the DRS-based consensus sequences contained numerous indel errors that differed from the HTS consensus sequences, and these errors hampered accurate open reading frame (ORF) prediction. Although with its present performance, the use of DRS is premature to determine viral genome sequences, the DRS-mediated sequencing shows great potential as a user-friendly platform for a one-shot, whole-genome sequencing of RNA viruses due to its long-reading ability and relative structure-tolerant nature.

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.

Complete genome sequence of Sikte (Sitke) waterborne virus, a member of the genus Tombusvirus

The 4704-nt genome sequence of Sikte waterborne virus (SWV), determined by fragmented and primer ligated dsRNA sequencing and by direct Sanger sequencing, is linear, nonsegmented and has the five ORFs of other tombusviruses. The 5' and 3' untranslated regions (UTRs) are 150 and 335 nt long, respectively. Phylogenetic analysis of the coat protein revealed that SWV is related to CymRSV and PNSV, but that of the SWV replicase protein, the p92 readthrough protein, indicated a close relationship to CNV. These phylogenetic analyses suggest the occurrence of recombination events in SWV, as reported previously for other tombusviruses.

Late Male-Killing Viruses in Homona magnanima Identified as Osugoroshi Viruses, Novel Members of Partitiviridae

Late male-killing, a male-specific death after hatching, is a unique phenotype found in Homona magnanima, oriental tea tortrix. The male-killing agent was suspected to be an RNA virus, but details were unknown. We herein successfully isolated and identified the putative male-killing virus as Osugoroshi viruses (OGVs). The three RNA-dependent RNA polymerase genes detected were phylogenetically related to Partitiviridae, a group of segmented double-stranded RNA viruses. Purified dsRNA from a late male-killing strain of H. magnanima revealed 24 segments, in addition to the RdRps, with consensus terminal sequences. These segments included the previously found male-killing agents MK1068 (herein OGV-related RNA16) and MK1241 (OGV-related RNA7) RNAs. Ultramicroscopic observation of purified virions, which induced late male-killing in the progeny of injected moths, showed sizes typical of Partitiviridae. Mathematical modeling showed the importance of late male-killing in facilitating horizontal transmission of OGVs in an H. magnanima population. This study is the first report on the isolation of partiti-like virus from insects, and one thought to be associated with late male-killing, although the viral genomic contents and combinations in each virus are still unknown.

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.

Phenotypic and Molecular Biological Analysis of Polymycovirus AfuPmV-1M From Aspergillus fumigatus: Reduced Fungal Virulence in a Mouse Infection Model

The filamentous fungal pathogen Aspergillus fumigatus is one of the most common causal agents of invasive fungal infection in humans; the infection is associated with an alarmingly high mortality rate. In this study, we investigated whether a mycovirus, named AfuPmV-1M, can reduce the virulence of A. fumigatus in a mouse infection model. AfuPmV-1M has high sequence similarity to AfuPmV-1, one of the polymycovirus that is a capsidless four-segment double-stranded RNA (dsRNA) virus, previously isolated from the genome reference strain of A. fumigatus, Af293. However, we found the isolate had an additional fifth dsRNA segment, referred to as open reading frame 5 (ORF5), which has not been reported in AfuPmV-1. We then established isogenic lines of virus-infected and virus-free A. fumigatus strains. Mycovirus infection had apparent influences on fungal phenotypes, with the virus-infected strain producing a reduced mycelial mass and reduced conidial number in comparison with these features of the virus-free strain. Also, resting conidia of the infected strain showed reduced adherence to pulmonary epithelial cells and reduced tolerance to macrophage phagocytosis. In an immunosuppressed mouse infection model, the virus-infected strain showed reduced mortality in comparison with mortality due to the virus-free strain. RNA sequencing and high-performance liquid chromatography (HPLC) analysis showed that the virus suppressed the expression of genes for gliotoxin synthesis and its production at the mycelial stage. Conversely, the virus enhanced gene expression and biosynthesis of fumagillin. Viral RNA expression was enhanced during conidial maturation, conidial germination, and the mycelial stage. We presume that the RNA or translation products of the virus affected fungal phenotypes, including spore formation and toxin synthesis. To identify the mycovirus genes responsible for attenuation of fungal virulence, each viral ORF was ectopically expressed in the virus-free KU strain. We found that the expression of ORF2 and ORF5 reduced fungal virulence in the mouse model. In addition, ORF3 affected the stress tolerance of host A. fumigatus in culture. We hypothesize that the respective viral genes work cooperatively to suppress the pathogenicity of the fungal host.

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.

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

Invertebrates are a source of previously unknown RNA viruses that fill gaps in the viral phylogenetic tree. Although limited information is currently available on RNA viral diversity in the marine sponge, a primordial multicellular animal that belongs to the phylum Porifera, the marine sponge is one of the well-studied holobiont systems. In the present study, we elucidated the putative complete genome sequences of five novel RNA viruses from Hymeniacidon sponge using a combination of double-stranded RNA sequencing, called fragmented and primer ligated dsRNA sequencing, and a conventional transcriptome method targeting single-stranded RNA. We identified highly diverged RNA-dependent RNA polymerase sequences, including a potential novel RNA viral lineage, in the sponge and three viruses presumed to infect sponge cells.

Comparative Analysis of RNA Virome Composition in Rabbits and Associated Ectoparasites

Ectoparasites play an important role in the transmission of many vertebrate-infecting viruses, including Zika and dengue viruses. Although it is becoming increasingly clear that invertebrate species harbor substantial virus diversity, it is unclear how many of the viruses carried by invertebrates have the potential to infect vertebrate species. We used the European rabbit (Oryctolagus cuniculus) as a model species to compare virome compositions in a vertebrate host and known associated ectoparasite mechanical vectors, in this case, fleas and blowflies. In particular, we aimed to infer the extent of viral transfer between these distinct types of host. Our analysis revealed that despite extensive viral diversity in both rabbits and associated ectoparasites, and the close interaction of these vertebrate and invertebrate species, biological viral transmission from ectoparasites to vertebrate species is rare. We did, however, find evidence to support the idea of a role of blowflies in transmitting viruses without active replication in the insect.

Ectoparasites play an important role in virus transmission among vertebrates. Little, however, is known about the nature of those viruses that pass between invertebrates and vertebrates. In Australia, flies and fleas support the mechanical transmission of two viral biological controls against wild rabbits—rabbit hemorrhagic disease virus (RHDV) and myxoma virus. We compared virome compositions in rabbits and these ectoparasites, sequencing total RNA from multiple tissues and gut contents of wild rabbits, fleas collected from these rabbits, and flies trapped sympatrically. Meta-transcriptomic analyses identified 50 novel viruses from multiple RNA virus families. Rabbits and their ectoparasites were characterized by markedly different viromes, with virus abundance greatest in flies. Although viral contigs from six virus families/groups were found in both rabbits and ectoparasites, they clustered in distinct host-dependent lineages. A novel calicivirus and a picornavirus detected in rabbit cecal content were vertebrate specific; the newly detected calicivirus was distinct from known rabbit caliciviruses, while the picornavirus clustered with sapeloviruses. Several picobirnaviruses were also identified that fell in diverse phylogenetic positions, compatible with the idea that they are associated with bacteria. Further comparative analysis revealed that the remaining viruses found in rabbits, and all those from ectoparasites, were likely associated with invertebrates, plants, and coinfecting endosymbionts. While no full genomes of vertebrate-associated viruses were detected in ectoparasites, small numbers of reads from rabbit astrovirus, RHDV, and other lagoviruses were present in flies. This supports a role for flies in the mechanical transmission of RHDV, while their involvement in astrovirus transmission merits additional exploration.

IMPORTANCE Ectoparasites play an important role in the transmission of many vertebrate-infecting viruses, including Zika and dengue viruses. Although it is becoming increasingly clear that invertebrate species harbor substantial virus diversity, it is unclear how many of the viruses carried by invertebrates have the potential to infect vertebrate species. We used the European rabbit (Oryctolagus cuniculus) as a model species to compare virome compositions in a vertebrate host and known associated ectoparasite mechanical vectors, in this case, fleas and blowflies. In particular, we aimed to infer the extent of viral transfer between these distinct types of host. Our analysis revealed that despite extensive viral diversity in both rabbits and associated ectoparasites, and the close interaction of these vertebrate and invertebrate species, biological viral transmission from ectoparasites to vertebrate species is rare. We did, however, find evidence to support the idea of a role of blowflies in transmitting viruses without active replication in the insect.