Metagenomics reshapes the concepts of RNA virus evolution by revealing extensive horizontal virus transfer

Technological advancements in deciphering RNA-RNA interactions

RNA-RNA interactions (RRIs) can dictate RNA molecules to form intricate higher-order structures and bind their RNA substrates in diverse biological processes. To elucidate the function, binding specificity, and regulatory mechanisms of various RNA molecules, especially the vast repertoire of non-coding RNAs, advanced technologies and methods that globally map RRIs are extremely valuable. In the past decades, many state-of-the-art technologies have been developed for this purpose. This review focuses on those high-throughput technologies for the global mapping of RRIs. We summarize the key concepts and the pros and cons of different technologies. In addition, we highlight the novel biological insights uncovered by these RRI mapping methods and discuss the future challenges for appreciating the crucial roles of RRIs in gene regulation across bacteria, viruses, archaea, and mammals.

Metagenomic analysis reveals ecological and functional signatures of oral phageome associated with severe early childhood caries

OBJECTIVES

Severe early childhood caries (S-ECC) is highly prevalent, affecting children's oral health. S-ECC development is closely associated with the complex oral microbial microbiome and its microorganism interactions, such as the imbalance of bacteriophages and bacteria. Till now, little is known about oral phageome on S-ECC. Therefore, this study aimed to investigate the potential role of the oral phageome in the pathogenesis of S-ECC.

METHODS

Unstimulated saliva (2 mL) was collected from 20 children with and without S-ECC for metagenomics analysis. Metagenomics sequencing and bioinformatic analysis were performed to determine the two groups' phageome diversity, taxonomic and functional annotations. Statistical analysis and visualization were performed with R and SPSS Statistics software.

RESULTS

85.7 % of the extracted viral sequences were predicted from phages, in which most phages were classified into Myoviridae, Siphoviridae, and Podoviridae. Alpha diversity decreased, and Beta diversity increased in the S-ECC phageome compared to the healthy group. The abundance of Podoviridae phages increased, and the abundance of Inoviridae, Herelleviridae, and Streptococcus phages decreased in the S-ECC group. Functional annotation revealed increased annotation on glycoside hydrolases and nucleotide metabolism, decreased glycosyl transferases, carbohydrate-binding modules, and biogenic metabolism in the S-ECC phageome.

CONCLUSIONS

Metagenomic analysis revealed reduced Streptococcus phages and significant changes in functional annotations within the S-ECC phageome. These findings suggest a potential weakening of the regulatory influence of oral bacteria, which may indicate the development of innovative prevention and treatment strategies for S-ECC. These implications deserve further investigation and hold promise for advancing our understanding and management of S-ECC.

CLINICAL SIGNIFICANCE

The findings of this study indicate that oral phageomes are associated with bacterial genomes and metabolic processes, affecting the development of S-ECC. The reduced modulatory effect of the oral phageome in counteracting S-ECC's cariogenic activity suggests a new avenue for the prevention and treatment of S-ECC.

Transmission of mycoviruses: new possibilities

Mycoviruses are viruses that infect fungi. In recent years, an increasing number of mycoviruses have been reported in a wide array of fungi. With the growing interest of scientists and society in reducing the use of agrochemicals, the debate about mycoviruses as an effective next-generation biocontrol has regained momentum. Mycoviruses can have profound effects on the host phenotype, although most viruses have neutral or no effect. We speculate that understanding multiple transmission modes of mycoviruses is central to unraveling the viral ecology and their function in regulating fungal populations. Unlike plant virus transmission via vegetative plant parts, seeds, pollen, or vectors, a widely held view is that mycoviruses are transmitted via vertical routes and only under special circumstances horizontally via hyphal contact depending on the vegetative compatibility groups (i.e., the ability of different fungal strains to undergo hyphal fusion). However, this view has been challenged over the past decades, as new possible transmission routes of mycoviruses are beginning to unravel. In this perspective, we discuss emerging studies with evidence suggesting that such novel routes of mycovirus transmission exist and are pertinent to understanding the full picture of mycovirus ecology and evolution.

Replication of single viruses across the kingdoms, Fungi, Plantae, and Animalia

It is extremely rare that a single virus crosses host barriers across multiple kingdoms. Based on phylogenetic and paleovirological analyses, it has previously been hypothesized that single members of the family Partitiviridae could cross multiple kingdoms. Partitiviridae accommodates members characterized by their simple bisegmented double-stranded RNA genome; asymptomatic infections of host organisms; the absence of an extracellular route for entry in nature; and collectively broad host range. Herein, we show the replicability of single fungal partitiviruses in three kingdoms of host organisms: Fungi, Plantae, and Animalia. Betapartitiviruses of the phytopathogenic fungusRosellinia necatrix could replicate in protoplasts of the carrot (Daucus carota), Nicotiana benthamiana and Nicotiana tabacum, in some cases reaching a level detectable by agarose gel electrophoresis. Moreover, betapartitiviruses showed more robust replication than the tested alphapartitiviruses. One of the fungal betapartitiviruses, RnPV18, could persistently and stably infect carrot plants regenerated from virion-transfected protoplasts. Both alpha- and betapartitiviruses, although with different host preference, could replicate in two insect cell lines derived from the fall armyworm Spodoptera frugiperda and the fruit fly Drosophila melanogaster. Our results indicate the replicability of single partitiviruses in members of three kingdoms and provide insights into virus adaptation, host jumping, and evolution.

Unveiling the hidden viromes across the animal tree of life: insights from a taxonomic classification pipeline applied to invertebrates of 31 metazoan phyla

Invertebrates constitute the majority of animal species on Earth, including most disease-causing agents or vectors, with more diverse viromes when compared to vertebrates. Recent advancements in high-throughput sequencing have significantly expanded our understanding of invertebrate viruses, yet this knowledge remains biased toward a few well-studied animal lineages. In this study, we analyze invertebrate DNA and RNA viromes for 31 phyla using 417 publicly available RNA-Seq data sets from diverse environments in the marine-terrestrial and marine-freshwater gradients. This study aims to (i) estimate virome compositions at the family level for the first time across the animal tree of life, including the first exploration of the virome in several phyla, (ii) quantify the diversity of invertebrate viromes and characterize the structure of invertebrate-virus infection networks, and (iii) investigate host phylum and habitat influence on virome differences. Results showed that a set of few viral families of eukaryotes, comprising Retroviridae, Flaviviridae, and several families of giant DNA viruses, were ubiquitous and highly abundant. Nevertheless, some differences emerged between phyla, revealing for instance a less diverse virome in Ctenophora compared to the other animal phyla. Compositional analysis of the viromes showed that the host phylum explained over five times more variance in composition than its habitat. Moreover, significant similarities were observed between the viromes of some phylogenetically related phyla, which could highlight the influence of co-evolution in shaping invertebrate viromes.IMPORTANCEThis study significantly enhances our understanding of the global animal virome by characterizing the viromes of previously unexamined invertebrate lineages from a large number of animal phyla. It showcases the great diversity of viromes within each phylum and investigates the role of habitat shaping animal viral communities. Furthermore, our research identifies dominant virus families in invertebrates and distinguishes phyla with analogous viromes. This study sets the road toward a deeper understanding of the virome across the animal tree of life.

Recurrent viral capture of cellular phosphodiesterases that antagonize OAS-RNase L

Phosphodiesterases (PDEs) encoded by viruses are putatively acquired by horizontal transfer of cellular PDE ancestor genes. Viral PDEs inhibit the OAS-RNase L antiviral pathway, a key effector component of the innate immune response. Although the function of these proteins is well-characterized, the origins of these gene acquisitions are less clear. Phylogenetic analysis revealed at least five independent PDE acquisition events by ancestral viruses. We found evidence that PDE-encoding genes were horizontally transferred between coronaviruses belonging to different genera. Three clades of viruses within Nidovirales: merbecoviruses (MERS-CoV), embecoviruses (HCoV-OC43), and toroviruses encode independently acquired PDEs, and a clade of rodent alphacoronaviruses acquired an embecovirus PDE via recent horizontal transfer. Among rotaviruses, the PDE of rotavirus A was acquired independently from rotavirus B and G PDEs, which share a common ancestor. Conserved motif analysis suggests a link between all viral PDEs and a similar ancestor among the mammalian AKAP7 proteins despite low levels of sequence conservation. Additionally, we used ancestral sequence reconstruction and structural modeling to reveal that sequence and structural divergence are not well-correlated among these proteins. Specifically, merbecovirus PDEs are as structurally divergent from the ancestral protein and the solved structure of human AKAP7 PDE as they are from each other. In contrast, comparisons of rotavirus B and G PDEs reveal virtually unchanged structures despite evidence for loss of function in one, suggesting impactful changes that lie outside conserved catalytic sites. These findings highlight the complex and volatile evolutionary history of viral PDEs and provide a framework to facilitate future studies.

Long-read genome sequencing provides novel insights into the harmful algal bloom species Prymnesium parvum

Prymnesium parvum is a toxin-producing haptophyte that causes harmful algal blooms worldwide, which are often associated with massive fish-kills and subsequent economic losses. In here, we present nuclear and plastid genome assemblies using PacBio HiFi long reads and DNBseq short reads for the two P. parvum strains UTEX 2797 and CCMP 3037, representing producers of type A prymnesins. Our results show that the P. parvum strains have a moderate haptophyte genome size of 97.56 and 107.32 Mb. The genome assemblies present one of highest contiguous assembled contig sequences to date consisting of 463 and 362 contigs with a contig N50 of 596.99 kb and 968.39 kb for strain UTEX 2797 and CCMP 3037, respectively. The assembled contigs of UTEX 2797 and CCMP 3037 were anchored to 34 scaffolds, with a scaffold N50 of 5.35 Mb and 3.61 Mb, respectively, accounting for 93.2 % and 97.9 % of the total length. Each plastid genome comprises a circular contig. A total of 20,578 and 19,426 protein-coding genes were annotated for UTEX 2797 and CCMP 3037. The expanded gene family analysis showed that starch and sucrose metabolism, sulfur metabolism, energy metabolism and ABC transporters are involved in the evolution of P. parvum. Polyketide synthase (PKS) genes responsible for the production of secondary metabolites such as prymnesins displayed different expression patterns under nutrient limitation. Overlap with repeats and horizontal gene transfer may be two contributing factors to the high number of PKS genes found in this species. The two high quality P. parvum genomes will serve as valuable resources for ecological, genetic, and toxicological studies of haptophytes that can be used to monitor and potentially manage harmful blooms of ichthyotoxic P. parvum in the future.

Megataxonomy and global ecology of the virosphere

Nearly all organisms are hosts to multiple viruses that collectively appear to be the most abundant biological entities in the biosphere. With recent advances in metagenomics and metatranscriptomics, the known diversity of viruses substantially expanded. Comparative analysis of these viruses using advanced computational methods culminated in the reconstruction of the evolution of major groups of viruses and enabled the construction of a virus megataxonomy, which has been formally adopted by the International Committee on Taxonomy of Viruses. This comprehensive taxonomy consists of six virus realms, which are aspired to be monophyletic and assembled based on the conservation of hallmark proteins involved in capsid structure formation or genome replication. The viruses in different major taxa substantially differ in host range and accordingly in ecological niches. In this review article, we outline the latest developments in virus megataxonomy and the recent discoveries that will likely lead to reassessment of some major taxa, in particular, split of three of the current six realms into two or more independent realms. We then discuss the correspondence between virus taxonomy and the distribution of viruses among hosts and ecological niches, as well as the abundance of viruses versus cells in different habitats. The distribution of viruses across environments appears to be primarily determined by the host ranges, i.e. the virome is shaped by the composition of the biome in a given habitat, which itself is affected by abiotic factors.

Viruses under the Antarctic Ice Shelf are active and potentially involved in global nutrient cycles

Viruses play an important role in the marine ecosystem. However, our comprehension of viruses inhabiting the dark ocean, and in particular, under the Antarctic Ice Shelves, remains limited. Here, we mine single-cell genomic, transcriptomic, and metagenomic data to uncover the viral diversity, biogeography, activity, and their role as metabolic facilitators of microbes beneath the Ross Ice Shelf. This is the largest Antarctic ice shelf with a major impact on global carbon cycle. The viral community found in the cavity under the ice shelf mainly comprises endemic viruses adapted to polar and mesopelagic environments. The low abundance of genes related to lysogenic lifestyle (<3%) does not support a predominance of the Piggyback-the-Winner hypothesis, consistent with a low-productivity habitat. Our results indicate a viral community actively infecting key ammonium and sulfur-oxidizing chemolithoautotrophs (e.g. Nitrosopumilus spp, Thioglobus spp.), supporting a "kill-the-winner" dynamic. Based on genome analysis, these viruses carry specific auxiliary metabolic genes potentially involved in nitrogen, sulfur, and phosphorus acquisition. Altogether, the viruses under Antarctic ice shelves are putatively involved in programming the metabolism of ecologically relevant microbes that maintain primary production in these chemosynthetically-driven ecosystems, which have a major role in global nutrient cycles.

A reference genome for the long-term kleptoplast-retaining sea slug Elysia crispata morphotype clarki

Several species of sacoglossan sea slugs possess the incredible ability to sequester chloroplasts from the algae they consume. These "photosynthetic animals" incorporate stolen chloroplasts, called kleptoplasts, into the epithelial cells of tubules that extend from their digestive tracts throughout their bodies. The mechanism by which these slugs maintain functioning kleptoplasts in the absence of an algal nuclear genome is unknown. Here, we report a draft genome of the sacoglossan slug Elysia crispata morphotype clarki, a morphotype native to the Florida Keys that can retain photosynthetically active kleptoplasts for several months without feeding. We used a combination of Oxford Nanopore Technologies long reads and Illumina short reads to produce a 786-Mb assembly (N50 = 0.459 Mb) containing 68,514 predicted protein-coding genes. A phylogenetic analysis found no evidence of horizontal acquisition of genes from algae. We performed gene family and gene expression analyses to identify E. crispata genes unique to kleptoplast-containing slugs that were more highly expressed in fed versus unfed developmental life stages. Consistent with analyses in other kleptoplastic slugs, our investigation suggests that genes encoding lectin carbohydrate-binding proteins and those involved in regulation of reactive oxygen species and immunity may play a role in kleptoplast retention. Lastly, we identified four polyketide synthase genes that could potentially encode proteins producing UV- and oxidation-blocking compounds in slug cell membranes. The genome of E. crispata is a quality resource that provides potential targets for functional analyses and enables further investigation into the evolution and mechanisms of kleptoplasty in animals.

Viral Diversity in Samples of Freshwater Gastropods Benedictia baicalensis (Caenogastropoda: Benedictiidae) Revealed by Total RNA-Sequencing

Previously, the main studies were focused on viruses that cause disease in commercial and farmed shellfish and cause damage to food enterprises (for example, Ostreavirusostreidmalaco1, Aurivirus haliotidmalaco1 and Aquabirnavirus tellinae). Advances in high-throughput sequencing technologies have extended the studies to natural populations of mollusks (and other invertebrates) as unexplored niches of viral diversity and possible sources of emerging diseases. These studies have revealed a huge diversity of mostly previously unknown viruses and filled gaps in the evolutionary history of viruses. In the present study, we estimated the viral diversity in samples of the Baikal endemic gastropod Benedictia baicalensis using metatranscriptomic analysis (total RNA-sequencing); we were able to identify a wide variety of RNA-containing viruses in four samples (pools) of mollusks collected at three stations of Lake Baikal. Most of the identified viral genomes (scaffolds) had only distant similarities to known viruses or (in most cases) to metagenome-assembled viral genomes from various natural samples (mollusks, crustaceans, insects and others) mainly from freshwater ecosystems. We were able to identify viruses similar to those previously identified in mollusks (in particular to the picornaviruses Biomphalaria virus 1 and Biomphalaria virus 3 from the freshwater gastropods); it is possible that picorna-like viruses (as well as a number of other identified viruses) are pathogenic for Baikal gastropods. Our results also suggested that Baikal mollusks, like other species, may bioaccumulate or serve as a reservoir for numerous viruses that infect a variety of organisms (including vertebrates).

Plant virology: an RNA treasure trove

Key principles pertaining to RNA biology not infrequently have their origins in plant virology. Examples have arisen from studies on viral RNA-intrinsic properties and the infection process from gene expression, replication, movement, and defense evasion to biotechnological applications. Since RNA is at the core of the central dogma in molecular biology, how plant virology assisted in the reinforcement or adaptations of this concept, while at other instances shook up elements of the doctrine, is discussed. Moreover, despite the negative effects of viral diseases in agriculture worldwide, plant viruses can be considered a scientific treasure trove. Today they remain tools of discovery for biotechnology, studying evolution, cell biology, and host-microbe interactions.

Alternative probe hybridization buffers for target RNA depletion and viral sequence recovery in NGS for poultry samples

Non-targeted next generation sequencing (NGS) is widely applied to identify the diversity of pathogens in field samples. However, abundance of host RNA (especially rRNA) and other environmental nucleic acids can reduce the abundance of pathogen specific reads of interest, reduce depth of coverage and increase surveillance costs. We presently deplete chicken- and selected bacterial-specific rRNAs in poultry field RNA samples with complementary DNA probes in a commercially available probe hybridization buffer followed by digestion of the RNA:DNA hybrids with RNase H. Because the current buffer is an expensive special order reagent of proprietary composition, we tested in-house and other commercially available buffers and identified a viable alternative that yields equivalent host rRNA depletion and viral-specific reads in poultry samples as the current special order reagent but at a reduced cost.

Molecular Characterization of Two Totiviruses from the Commensal Yeast Geotrichum candidum

Mycoviruses can infect many of the major taxa of fungi including yeasts. Mycoviruses in the yeast fungus Geotrichum candidum are not well studied with only three G. candidum-associated viral species characterized to date, all of which belong to the Totiviridae genus Totivirus. In this study, we report the molecular characteristics of another two totiviruses co-infecting isolate Gc6 of G. candidum. The two totiviruses were tentatively named Geotrichum candidum totivirus 2 isolate Gc6 (GcTV2-Gc6) and Geotrichum candidum totivirus 4 isolate Gc6 (GcTV4-Gc6). Both viruses have the typical genome organization of totiviruses comprising two ORFs encoding capsid protein (CP) and RNA-dependent RNA polymerase (RdRp) at the N and C termini, respectively. The genomes of GcTV2-Gc6 and GcTV4-Gc6 are 4592 and 4530 bp long, respectively. Both viruses contain the-frameshifting elements and their proteins could be expressed as a single fusion protein. GcTV2-Gc6 is closely related to a totivirus isolated from the same host whereas GcTV4-Gc6 is related to insect-associated totiviruses. The phylogenetic analysis indicated that GcTV2-Gc6 and GcTV4-Gc6 belong to two different sister clades, I-A and I-B, respectively. It is interesting that all viruses identified from G. candidum belong to the genus Totivirus; however, this might be due to the lack of research reporting the characterization of mycoviruses from this fungal host. It is possible that the RNA interference (RNAi) mechanism cannot actively suppress totivirus accumulation in G. candidum Gc6.

Recurrent Viral Capture of Cellular Phosphodiesterases that Antagonize OAS-RNase L

Phosphodiesterases (PDEs) encoded by viruses are putatively acquired by horizontal transfer of cellular PDE ancestor genes. Viral PDEs inhibit the OAS-RNase L antiviral pathway, a key effector component of the innate immune response. Although the function of these proteins is well-characterized, the origins of these gene acquisitions is less clear. Phylogenetic analysis revealed at least five independent PDE acquisition events by ancestral viruses. We found evidence that PDE-encoding genes were horizontally transferred between coronavirus genera. Three clades of viruses within Nidovirales: merbecoviruses (MERS-CoV), embecoviruses (OC43), and toroviruses encode independently acquired PDEs, and a clade of rodent alphacoronaviruses acquired an embecovirus PDE via recent horizontal transfer. Among rotaviruses, the PDE of Rotavirus A was acquired independently from Rotavirus B and G PDEs, which share a common ancestor. Conserved motif analysis suggests a link between all viral PDEs and a similar ancestor among the mammalian AKAP7 proteins despite low levels of sequence conservation. Additionally, we used ancestral sequence reconstruction and structural modeling to reveal that sequence and structural divergence are not well-correlated among these proteins. Specifically, merbecovirus PDEs are as structurally divergent from the ancestral protein and the solved structure of human AKAP7 PDE as they are from each other. In contrast, comparisons of Rotavirus B and G PDEs reveal virtually unchanged structures despite evidence for loss of function in one, suggesting impactful changes that lie outside conserved catalytic sites. These findings highlight the complex and volatile evolutionary history of viral PDEs and provide a framework to facilitate future studies.

Virome analysis of an ectomycorrhizal fungus Suillus luteus revealing potential evolutionary implications

Suillus luteus is a widespread edible ectomycorrhizal fungus that holds significant importance in both ecological and economic value. Mycoviruses are ubiquitous infectious agents hosted in different fungi, with some known to exert beneficial or detrimental effects on their hosts. However, mycoviruses hosted in ectomycorrhizal fungi remain poorly studied. To address this gap in knowledge, we employed next-generation sequencing (NGS) to investigate the virome of S. luteus. Using BLASTp analysis and phylogenetic tree construction, we identified 33 mycovirus species, with over half of them belonging to the phylum Lenarviricota, and 29 of these viruses were novel. These mycoviruses were further grouped into 11 lineages, with the discovery of a new negative-sense single-stranded RNA viral family in the order Bunyavirales. In addition, our findings suggest the occurrence of cross-species transmission (CST) between the fungus and ticks, shedding light on potential evolutionary events that have shaped the viral community in different hosts. This study is not only the first study to characterize mycoviruses in S. luteus but highlights the enormous diversity of mycoviruses and their implications for virus evolution.

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.

Two novel flavi-like viruses shed light on the plant-infecting koshoviruses

The family Flaviviridae is composed of viruses with a positive-sense single-stranded RNA genome and includes viruses that are important veterinary and human pathogens. Most members of the family are arthropod- and vertebrate-infecting viruses, but more recently, divergent flavi-like viruses have been identified in marine invertebrate and vertebrate hosts. The striking discovery of gentian Kobu-sho-associated virus (GKaV), along with a recent report of a related virus from carrot, has expanded the known host range of flavi-like viruses to plants, suggesting they could be grouped in a proposed genus tentatively named "Koshovirus". Here, we report the identification and characterization of two novel RNA viruses that show a genetic and evolutionary relationship to the previously identified "koshoviruses". Their genome sequences were obtained from transcriptomic datasets of the flowering plants Coptis teeta and Sonchus asper. These two new viruses, which we have named "coptis flavi-like virus 1" (CopFLV1) and "sonchus flavi-like virus 1" (SonFLV1), are members of novel species characterized by the longest monopartite RNA genome observed so far among plant-associated RNA viruses, which is ca. 24 kb in size. Structural and functional annotations of the polyproteins of all koshoviruses resulted in the detection not only of the expected helicase and RNA-dependent RNA polymerase but also of several additional divergent domains, including AlkB oxygenase, trypsin-like serine protease, methyltransferase, and envelope E1 flavi-like domains. Phylogenetic analysis showed that CopFLV1, SonFLV1, GKaV, and the carrot flavi-like virus were grouped together in a monophyletic clade, strongly supporting the recent proposal for creation of the genus "Koshovirus" for the group of related plant-infecting flavi-like viruses.

Plant virus movement proteins originated from jelly-roll capsid proteins

Numerous, diverse plant viruses encode movement proteins (MPs) that aid the virus movement through plasmodesmata, the plant intercellular channels. MPs are essential for virus spread and propagation in distal tissues, and several unrelated MPs have been identified. The 30K superfamily of MPs (named after the molecular mass of tobacco mosaic virus (TMV) MP, the classical model of plant virology) is the largest and most diverse MP variety, represented in 16 virus families, but its evolutionary origin remained obscure. Here, we show that the core structural domain of the 30K MPs is homologous to the jelly-roll domain of the capsid proteins (CPs) of small RNA and DNA viruses, in particular, those infecting plants. The closest similarity was observed between the 30K MPs and the CPs of the viruses in the families Bromoviridae and Geminiviridae. We hypothesize that the MPs evolved via duplication or horizontal acquisition of the CP gene in a virus that infected an ancestor of vascular plants, followed by neofunctionalization of one of the paralogous CPs, potentially through the acquisition of unique N- and C-terminal regions. During the subsequent coevolution of viruses with diversifying vascular plants, the 30K MP genes underwent explosive horizontal spread among emergent RNA and DNA viruses, likely permitting viruses of insects and fungi that coinfected plants to expand their host ranges, molding the contemporary plant virome.

On the nature of picobirnaviruses

The picobirnaviruses (Picobirnaviridae, Picobirnavirus, PBVs) are currently thought to be animal viruses, as they are usually found in animal stool samples. However, no animal model or cell culture for their propagation has yet been found. In 2018, a hypothetical assumption about PBVs belonging to prokaryotic viruses was put forward and experimentally substantiated. This hypothesis is based on the presence of Shine-Dalgarno sequences in the genome of all PBVs before three reading frames (ORF) at the ribosomal binding site, with which the prokaryotic genome is saturated, while in the eukaryotic genome such regions occur with low frequency. The genome saturation with the Shine-Dalgarno sequences, as well as the preservation of this saturation in the progeny, according to scientists, allows us to attribute PBVs to prokaryotic viruses. On the other hand, there is a possibility that PBVs belong to viruses of eukaryotic hosts - fungi or invertebrates, since PBV-like sequences similar to the genome of fungal viruses from the families of mitoviruses and partitiviruses have been identified. In this regard, the idea arose that, in terms of reproduction mode, PBVs resemble fungal viruses. The divergence of views on the true PBV host(s) has sparked discussions among scientists and required further research to elucidate their nature. The review highlights the results of the search for a PBV host. The reasons for the occurrence of atypical sequences among the PBV genome sequences that use an alternative mitochondrial code of lower eukaryotes (fungi and invertebrates) for the translation of viral RNA-dependent RNA polymerase (RdRp) instead of the standard genetic code are analyzed. The purpose of the review was to collect arguments in support of the hypothesis about the phage nature of PBVs and to find the most realistic explanation of the reasons for identifying non-standard genomic sequences for PBVs. Based on the hypothesis about the genealogical relationship of PBVs with RNA viruses from other families with similar segmented genomes, such as Reoviridae, Cystoviridae, Totiviridae and Partitiviridae, virologists support the assumption of a decisive role in the origin of atypical PBV-like reassortment strains between PBVs and viruses of the listed families. The collected arguments given in this review indicate a high probability of a phage nature of PBVs. The data presented in the review show that the belonging of PBV-like progeny to prokaryotic or eukaryotic viruses is determined not only by its genome saturation level with a prokaryotic motif, standard or mitochondrial genetic code. The primary structure of the gene encoding the viral capsid protein responsible for the presence or absence of specific proteolytic properties of the virus that determine its ability for independent horizontal transmission into new cells may also be a decisive factor.

Divergent RNA viruses infecting sea lice, major ectoparasites of fish

Sea lice, the major ectoparasites of fish, have significant economic impacts on wild and farmed finfish, and have been implicated in the decline of wild salmon populations. As blood-feeding arthropods, sea lice may also be reservoirs for viruses infecting fish. However, except for two groups of negative-strand RNA viruses within the order Mononegavirales, nothing is known about viruses of sea lice. Here, we used transcriptomic data from three key species of sea lice (Lepeophtheirus salmonis, Caligus clemensi, and Caligus rogercresseyi) to identify 32 previously unknown RNA viruses. The viruses encompassed all the existing phyla of RNA viruses, with many placed in deeply branching lineages that likely represent new families and genera. Importantly, the presence of canonical virus-derived small interfering RNAs (viRNAs) indicates that most of these viruses infect sea lice, even though in some cases their closest classified relatives are only known to infect plants or fungi. We also identified both viRNAs and PIWI-interacting RNAs (piRNAs) from sequences of a bunya-like and two qin-like viruses in C. rogercresseyi. Our analyses showed that most of the viruses found in C. rogercresseyi occurred in multiple life stages, spanning from planktonic to parasitic stages. Phylogenetic analysis revealed that many of the viruses infecting sea lice were closely related to those that infect a wide array of eukaryotes with which arthropods associate, including fungi and parasitic tapeworms, implying that over evolutionary time there has been cross-phylum and cross-kingdom switching of viruses between arthropods and other eukaryotes. Overall, this study greatly expands our view of virus diversity in crustaceans, identifies viruses that infect and replicate in sea lice, and provides evidence that over evolutionary time, viruses have switched between arthropods and eukaryotic hosts in other phyla and kingdoms.

The virome of the last eukaryotic common ancestor and eukaryogenesis

All extant eukaryotes descend from the last eukaryotic common ancestor (LECA), which is thought to have featured complex cellular organization. To gain insight into LECA biology and eukaryogenesis-the origin of the eukaryotic cell, which remains poorly understood-we reconstructed the LECA virus repertoire. We compiled an inventory of eukaryotic hosts of all major virus taxa and reconstructed the LECA virome by inferring the origins of these groups of viruses. The origin of the LECA virome can be traced back to a small set of bacterial-not archaeal-viruses. This provenance of the LECA virome is probably due to the bacterial origin of eukaryotic membranes, which is most compatible with two endosymbiosis events in a syntrophic model of eukaryogenesis. In the first endosymbiosis, a bacterial host engulfed an Asgard archaeon, preventing archaeal viruses from entry owing to a lack of archaeal virus receptors on the external membranes.

Identification of RNA Virus-Derived RdRp Sequences in Publicly Available Transcriptomic Data Sets

RNA viruses are abundant and highly diverse and infect all or most eukaryotic organisms. However, only a tiny fraction of the number and diversity of RNA virus species have been catalogued. To cost-effectively expand the diversity of known RNA virus sequences, we mined publicly available transcriptomic data sets. We developed 77 family-level Hidden Markov Model profiles for the viral RNA-dependent RNA polymerase (RdRp)-the only universal "hallmark" gene of RNA viruses. By using these to search the National Center for Biotechnology Information Transcriptome Shotgun Assembly database, we identified 5,867 contigs encoding RNA virus RdRps or fragments thereof and analyzed their diversity, taxonomic classification, phylogeny, and host associations. Our study expands the known diversity of RNA viruses, and the 77 curated RdRp Profile Hidden Markov Models provide a useful resource for the virus discovery community.

Arboviruses and symbiotic viruses cooperatively hijack insect sperm-specific proteins for paternal transmission

Arboviruses and symbiotic viruses can be paternally transmitted by male insects to their offspring for long-term viral persistence in nature, but the mechanism remains largely unknown. Here, we identify the sperm-specific serpin protein HongrES1 of leafhopper Recilia dorsalis as a mediator of paternal transmission of the reovirus Rice gall dwarf virus (RGDV) and a previously undescribed symbiotic virus of the Virgaviridae family, Recilia dorsalis filamentous virus (RdFV). We show that HongrES1 mediates the direct binding of virions to leafhopper sperm surfaces and subsequent paternal transmission via interaction with both viral capsid proteins. Direct interaction of viral capsid proteins mediates simultaneously invasion of two viruses into male reproductive organs. Moreover, arbovirus activates HongrES1 expression to suppress the conversion of prophenoloxidase to active phenoloxidase, potentially producing a mild antiviral melanization defense. Paternal virus transmission scarcely affects offspring fitness. These findings provide insights into how different viruses cooperatively hijack insect sperm-specific proteins for paternal transmission without disturbing sperm functions.

Host–Pathogen Interactions Influencing Zoonotic Spillover Potential and Transmission in Humans

Emerging infectious diseases of zoonotic origin are an ever-increasing public health risk and economic burden. The factors that determine if and when an animal virus is able to spill over into the human population with sufficient success to achieve ongoing transmission in humans are complex and dynamic. We are currently unable to fully predict which pathogens may appear in humans, where and with what impact. In this review, we highlight current knowledge of the key host–pathogen interactions known to influence zoonotic spillover potential and transmission in humans, with a particular focus on two important human viruses of zoonotic origin, the Nipah virus and the Ebola virus. Namely, key factors determining spillover potential include cellular and tissue tropism, as well as the virulence and pathogenic characteristics of the pathogen and the capacity of the pathogen to adapt and evolve within a novel host environment. We also detail our emerging understanding of the importance of steric hindrance of host cell factors by viral proteins using a “flytrap”-type mechanism of protein amyloidogenesis that could be crucial in developing future antiviral therapies against emerging pathogens. Finally, we discuss strategies to prepare for and to reduce the frequency of zoonotic spillover occurrences in order to minimize the risk of new outbreaks.

The Virome of Bats Inhabiting Brazilian Biomes: Knowledge Gaps and Biases towards Zoonotic Viruses

Bats host a large variety of viruses, including some that may infect other vertebrates and humans. Research on bat-borne viruses attracted significant attention in recent years mainly due to epizootics caused by viruses having bats as hosts. The characterization of the viral communities of bats was then prioritized, but despite increasing efforts, there are large disparities in the geographical ranges covered and the methodologies employed around the world. As a result, large gaps remain in our current understanding of bat viromes and their role in disease emergence. This is particularly true for megadiverse regions in Latin America. This review aims to summarize the current understanding about bat viruses that inhabit Brazilian biomes, one of the most bat species-rich and diverse regions of the globe. Taking into account all known bat-associated viral families studied in Brazilian biomes, we found that almost half of all bat species (86/181 species) were not investigated for viruses at all. Moreover, only a small fraction of viral lineages or families have been studied more in depth, usually employing targeted methods with limited power to characterize a broad virus diversity. Additionally, these studies relied on limited spatiotemporal sampling and small sample sizes. Therefore, our current understanding of bat viral communities in the Brazilian biomes is limited and biased at different levels, limiting zoonotic risk assessments of bat-borne viruses. Considering these limitations, we propose strategies to bridge the existing gaps in the near future. IMPORTANCE Bat-borne viruses have attracted much attention due to zoonotic outbreaks with large consequences to humans. Because of that, virus characterization in bats has been prioritized in tropical regions of the globe. However, bat-virus research in Latin America and particularly in Brazil, which are among the most bat species-rich regions of the world, are highly biased toward zoonotic viruses and known bat reservoir species. These results have direct implication for virus studies in general but also for new zoonotic virus and spillover events characterization. The limited knowledge we currently have about the virome of Brazilian bats drastically limits any broad assessment of zoonotic viruses they carry and calls for coordinated and large-scale studies to fill this crucial knowledge gap.

Putative Mitoviruses without In-Frame UGA(W) Codons: Evolutionary Implications

Mitoviruses are small vertically transmitted RNA viruses found in fungi, plants and animals. Taxonomically, a total of 105 species and 4 genera have been formally recognized by ICTV, and recently, 18 new putative species have been included in a new proposed genus. Transcriptomic and metatranscriptomic studies are a major source of countless new virus-like sequences that are continually being added to open databases and these may be good sources for identifying new putative mitoviruses. The search for mitovirus-like sequences in the NCBI databases resulted in the discovery of more than one hundred new putative mitoviruses, with important implications for taxonomy and also for the evolutionary scenario. Here, we propose the inclusion of four new putative members to the genus Kvaramitovirus, and the existence of a new large basally divergent lineage composed of 144 members that lack internal UGA codons (subfamily "Arkeomitovirinae"), a feature not shared by the vast majority of mitoviruses. Finally, a taxonomic categorization proposal and a detailed description of the evolutionary history of mitoviruses were carried out. This in silico study supports the hypothesis of the existence of a basally divergent lineage that could have had an impact on the early evolutionary history of mitoviruses.

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.

The global virome: How much diversity and how many independent origins?

Viruses are considered to be the most abundant biological entities on earth. They also display striking genetic diversity as emphatically demonstrated by the recent advances of metagenomics and metatranscriptomics. But what are the limits of this diversity, that is, how many virus species in the earth virome? By combining the available estimates of the number of prokaryote species with those of the virome size, we obtain back-of-the-envelope estimates of the total number of distinct virus species, which come out astronomically large, from about 10⁷ to about 10⁹ . The route of virus origins apparently involved non-viral replicators capturing and exapting various cellular proteins to become virus capsid subunits. How many times in the history of life has this happened? In other words, how many realms of viruses, the highest rank taxa that are supposed to be monophyletic, comprise the global virome? We argue that viruses emerged on a number (even if far from astronomical) independent occasions, so the number of realms will considerably increase from the current 6, by splitting some of the current realms, giving the realm status to some of the currently unclassified groups of viruses and discovery of new distinct groups.

Transcriptome Mining Expands Knowledge of RNA Viruses across the Plant Kingdom

Our current understanding of plant viruses stems largely from those affecting economically important plants. Yet plant species in cultivation represent a small and biased subset of the plant kingdom. Here, we describe virus diversity and abundance in 1,079 transcriptomes from species across the breadth of the plant kingdom (Archaeplastida) by analyzing open-source data from the 1000 Plant Transcriptomes Initiative (1KP). We identified 104 potentially novel viruses, of which 40% were single-stranded positive-sense RNA viruses across eight orders, including members of the Hepelivirales, Tymovirales, Cryppavirales, Martellivirales, and Picornavirales. One-third of the newly described viruses were double-stranded RNA viruses from the orders Durnavirales and Ghabrivirales. The remaining were negative-sense RNA viruses from the Rhabdoviridae, Aspiviridae, Yueviridae, and Phenuiviridae and the newly proposed Viridisbunyaviridae. Our analysis considerably expands the known host range of 13 virus families to include lower plants (e.g., Benyviridae and Secoviridae) and 4 virus families to include alga hosts (e.g., Tymoviridae and Chrysoviridae). More broadly, however, a cophylogeny analysis revealed that the evolutionary history of these families is largely driven by cross-species transmission events. The discovery of the first 30-kDa movement protein in a nonvascular plant suggests that the acquisition of plant virus movement proteins occurred prior to the emergence of the plant vascular system. Together, these data highlight that numerous RNA virus families are associated with older evolutionary plant lineages than previously thought and that the apparent scarcity of RNA viruses found in lower plants likely reflects a lack of investigation rather than their absence. IMPORTANCE Our knowledge of plant viruses is mainly limited to those infecting economically important host species. In particular, we know little about those viruses infecting basal plant lineages such as the ferns, lycophytes, bryophytes, and charophytes. To expand this understanding, we conducted a broad-scale viral survey of species across the breadth of the plant kingdom. We found that basal plants harbor a wide diversity of RNA viruses, including some that are sufficiently divergent to likely compose a new virus family. The basal plant virome revealed offers key insights into the evolutionary history of core plant virus gene modules and genome segments. More broadly, this work emphasizes that the scarcity of viruses found in these species to date most likely reflects the limited research in this area.

Virome Analysis of Normal and Growth Retardation Disease-Affected Macrobrachium rosenbergii

The giant freshwater prawn, Macrobrachium rosenbergii, is an important aquaculture species in China. Growth retardation disease (GRD) is a common contagious disease in M. rosenbergii, resulting in slow growth and precocious puberty in prawns, and has caused growing economic losses in the M. rosenbergii industry. To investigate the viral diversity of M. rosenbergii and identify potentially high-risk viruses linked to GRD, virome analysis of the GRD-affected and normal M. rosenbergii was carried out using next-generation sequencing (NGS). A total of 327 contigs (>500 bp) were related to viral sequences belonging to 23 families/orders and a group of unclassified viruses. The majority of the viral contigs in M. rosenbergii belonged to the order Picornavirales, with the Solinviviridae family being the most abundant in both the diseased and normal groups. Furthermore, 16 RNA viral sequences with nearly complete genomes were characterized and phylogenetically analyzed, belonging to the families Solinviviridae, Flaviviridae, Polycipiviridae, Marnaviridae, and Dicistroviridae as well as three new clades of the order Picornavirales. Notably, the cross-species transmission of a picorna-like virus was observed between M. rosenbergii and plants. The "core virome" seemed to be present in the diseased and normal prawns. Still, a clear difference in viral abundance was observed between the two groups. These results showed that the broad diversity of viruses is present in M. rosenbergii and that the association between viruses and disease of M. rosenbergii needs to be further investigated. IMPORTANCE Growth retardation disease (GRD) has seriously affected the development and economic growth of the M. rosenbergii aquaculture industry. Our virome analysis showed that diverse viral sequences were present in M. rosenbergii, significantly expanding our knowledge of viral diversity in M. rosenbergii. Some differences in viral composition were noted between the diseased and normal prawns, indicating that some viruses become more abundant in occurrences or outbreaks of diseases. In the future, more research will be needed to determine which viruses pose a risk for M. rosenbergii. Our study provides important baseline information contributing to disease surveillance and risk assessment in M. rosenbergii aquaculture.

Over two decades of research on the marine RNA virosphere

RNA viruses (realm: Riboviria), including RNA phages and eukaryote-infecting RNA viruses, are essential components of marine ecosystems. A large number of marine RNA viruses have been discovered in the last two decades because of the rapid development of next-generation sequencing (NGS) technology. Indeed, the combination of NGS and state-of-the-art meta-omics methods (viromics, the study of all viruses in a specific environment) has led to a fundamental understanding of the taxonomy and genetic diversity of RNA viruses in the sea, suggesting the complex ecological roles played by RNA viruses in this complex ecosystem. Furthermore, comparisons of viromes in the context of highly variable marine niches reveal the biogeographic patterns and ecological impact of marine RNA viruses, whose role in global ecology is becoming increasingly clearer. In this review, we summarize the characteristics of the global marine RNA virosphere and outline the taxonomic hierarchy of RNA viruses with a specific focus on their ancient evolutionary history. We also review the development of methodology and the major progress resulting from its applications in RNA viromics. The aim of this review is not only to provide an in-depth understanding of multifaceted aspects of marine RNA viruses, but to offer future perspectives on developing a better methodology for discovery, and exploring the evolutionary origin and major ecological significance of marine RNA virosphere.

Human Clinical Isolates of Pathogenic Fungi Are Host to Diverse Mycoviruses

Fungi host viruses from many families, and next-generation sequencing can be used to discover previously unknown genomes. Some fungus-infecting viruses (mycoviruses) confer hypovirulence on their pathogenic hosts, raising the possibility of therapeutic application in the treatment of fungal diseases. Though all fungi probably host mycoviruses, many human pathogens have none documented, implying the mycoviral catalogue remains at an early stage. Here, we carried out virus discovery on 61 cultures of pathogenic fungi covering 27 genera and at least 56 species. Using next-generation sequencing of total nucleic acids, we found no DNA viruses but did find a surprising RNA virus diversity of 11 genomes from six classified families and two unclassified lineages, including eight genomes likely representing new species. Among these was the first jivivirus detected in a fungal host (Aspergillus lentulus). We separately utilized rolling circle amplification and next-generation sequencing to identify ssDNA viruses specifically. We identified 13 new cressdnaviruses across all libraries, but unlike the RNA viruses, they could not be confirmed by PCR in either the original unamplified samples or freshly amplified nucleic acids. Their distributions among sequencing libraries and inconsistent detection suggest low-level contamination of reagents. This highlights both the importance of validation assays and the risks of viral host prediction on the basis of highly amplified sequencing libraries. Meanwhile, the detected RNA viruses provide a basis for experimentation to characterize possible hypovirulent effects, and hint at a wealth of uncharted viral diversity currently frozen in biobanks. IMPORTANCE Fungal pathogens of humans are a growing global health burden. Viruses of fungi may represent future therapeutic tools, but for many fungal pathogens there are no known viruses. Our study examined the viral content of diverse human-pathogenic fungi in a clinical biobank, identifying numerous viral genomes, including one lineage previously not known to infect fungi.

Expansion of the global RNA virome reveals diverse clades of bacteriophages

High-throughput RNA sequencing offers broad opportunities to explore the Earth RNA virome. Mining 5,150 diverse metatranscriptomes uncovered >2.5 million RNA virus contigs. Analysis of >330,000 RNA-dependent RNA polymerases (RdRPs) shows that this expansion corresponds to a 5-fold increase of the known RNA virus diversity. Gene content analysis revealed multiple protein domains previously not found in RNA viruses and implicated in virus-host interactions. Extended RdRP phylogeny supports the monophyly of the five established phyla and reveals two putative additional bacteriophage phyla and numerous putative additional classes and orders. The dramatically expanded phylum Lenarviricota, consisting of bacterial and related eukaryotic viruses, now accounts for a third of the RNA virome. Identification of CRISPR spacer matches and bacteriolytic proteins suggests that subsets of picobirnaviruses and partitiviruses, previously associated with eukaryotes, infect prokaryotic hosts.

From Contagium vivum fluidum to Riboviria: A Tobacco Mosaic Virus-Centric History of Virus Taxonomy

Viruses were discovered as agents of disease in the late 19th century, but it was not until the 1930s that the nature of these agents was elucidated. Nevertheless, as soon as viral diseases started to be recognized and cataloged, there were attempts to classify and name viruses. Although these early attempts failed to be adopted by the nascent virology community, they are evidence of the human compulsion to try to organize the natural world into well-defined categories. Different classification schemes were proposed during the 20th century, but again none were widely embraced by virologists. In 1966, with the creation of the International Committee on Nomenclature of Viruses (eventually renamed as the International Committee on Taxonomy of Viruses), a more organized effort led to an official taxonomy in which viruses were classified into families and genera. At present, a much better understanding of the evolutionary relationships among viruses has led to the establishment of a 15-rank taxonomy based primarily on these evolutionary relationships. This review of virus taxonomy will be centered on the tobacco mosaic virus (TMV), the agent of the disease studied by Dmitry Ivanovsky and the first virus to be recognized as such, which was often historically at the center of major advancements in virology during the 20th century.

Aquatic Flaviviruses

Flaviviruses are positive-sense single-stranded RNA viruses, including some well-known human pathogens such as Zika, dengue, and yellow fever viruses, which are primarily associated with mosquito and tick vectors. The vast majority of flavivirus research has focused on terrestrial environments; however, recent findings indicate that a range of flaviviruses are also present in aquatic environments, both marine and freshwater. These flaviviruses are found in various hosts, including fish, crustaceans, molluscs, and echinoderms. Although the effects of aquatic flaviviruses on the hosts they infect are not all known, some have been detected in farmed species and may have detrimental effects on the aquaculture industry. Exploration of the evolutionary history through the discovery of the Wenzhou shark flavivirus in both a shark and crab host is of particular interest since the potential dual-host nature of this virus may indicate that the invertebrate-vertebrate relationship seen in other flaviviruses may have a more profound evolutionary root than previously expected. Potential endogenous viral elements and the range of novel aquatic flaviviruses discovered thus shed light on virus origins and evolutionary history and may indicate that, like terrestrial life, the origins of flaviviruses may lie in aquatic environments.

Poxviruses capture host genes by LINE-1 retrotransposition

Horizontal gene transfer (HGT) provides a major source of genetic variation. Many viruses, including poxviruses, encode genes with crucial functions directly gained by gene transfer from hosts. The mechanism of transfer to poxvirus genomes is unknown. Using genome analysis and experimental screens of infected cells, we discovered a central role for Long Interspersed Nuclear Element-1 retrotransposition in HGT to virus genomes. The process recapitulates processed pseudogene generation, but with host messenger RNA directed into virus genomes. Intriguingly, hallmark features of retrotransposition appear to favor virus adaption through rapid duplication of captured host genes on arrival. Our study reveals a previously unrecognized conduit of genetic traffic with fundamental implications for the evolution of many virus classes and their hosts.

Comparative analysis of genome-encoded viral sequences reveals the evolutionary history of flavivirids (family Flaviviridae)

Flavivirids (family Flaviviridae) are a group of positive-strand ribonucleic acid (RNA) viruses that pose serious risks to human and animal health on a global scale. Here, we use flavivirid-derived deoxyribonucleic acid (DNA) sequences, identified in animal genomes, to reconstruct the long-term evolutionary history of family Flaviviridae. We demonstrate that flavivirids are >100 million years old and show that this timing can be combined with dates inferred from co-phyletic analysis to produce a cohesive overview of their evolution, distribution, and diversity wherein the main flavivirid subgroups originate in early animals and broadly co-diverge with major animal phyla. In addition, we reveal evidence that the 'classical flaviviruses' of vertebrates, most of which are transmitted via blood-feeding arthropod vectors, originally evolved in haematophagous arachnids and later acquired the capacity to be transmitted by insects. Our findings imply that the biological properties of flavivirids have been acquired gradually over the course of animal evolution. Thus, broad-scale comparative analysis will likely reveal fundamental insights into their biology. We therefore published our results via an open, extensible, database (Flavivirid-GLUE), which we constructed to facilitate the wider utilisation of genomic data and evolution-related domain knowledge in flavivirid research.

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.

Primordial Capsid and Spooled ssDNA Genome Structures Unravel Ancestral Events of Eukaryotic Viruses

Marine algae viruses are important for controlling microorganism communities in the marine ecosystem and played fundamental roles during the early events of viral evolution. Here, we have focused on one major group of marine algae viruses, the single-stranded DNA (ssDNA) viruses from the Bacilladnaviridae family. We present the capsid structure of the bacilladnavirus Chaetoceros tenuissimus DNA virus type II (CtenDNAV-II), determined at 2.4-Å resolution. A structure-based phylogenetic analysis supported the previous theory that bacilladnaviruses have acquired their capsid protein via horizontal gene transfer from a ssRNA virus. The capsid protein contains the widespread virus jelly-roll fold but has additional unique features; a third β-sheet and a long C-terminal tail. Furthermore, a low-resolution reconstruction of the CtenDNAV-II genome revealed a partially spooled structure, an arrangement previously only described for dsRNA and dsDNA viruses. Together, these results exemplify the importance of genetic recombination for the emergence and evolution of ssDNA viruses and provide important insights into the underlying mechanisms that dictate genome organization. IMPORTANCE Single-stranded DNA (ssDNA) viruses are an extremely widespread group of viruses that infect diverse hosts from all three domains of life, consequently having great economic, medical, and ecological importance. In particular, bacilladnaviruses are highly abundant in marine sediments and greatly influence the dynamic appearance and disappearance of certain algae species. Despite the importance of ssDNA viruses and the last couple of years' advancements in cryo-electron microscopy, structural information on the genomes of ssDNA viruses remains limited. This paper describes two important achievements: (i) the first atomic structure of a bacilladnavirus capsid, which revealed that the capsid protein gene presumably was acquired from a ssRNA virus in early evolutionary events; and (ii) the structural organization of a ssDNA genome, which retains a spooled arrangement that previously only been observed for double-stranded viruses.

Novel RNA viruses associated with avian haemosporidian parasites

Avian haemosporidian parasites can cause malaria-like symptoms in their hosts and have been implicated in the demise of some bird species. The newly described Matryoshka RNA viruses (MaRNAV1 and MaRNAV2) infect haemosporidian parasites that in turn infect their vertebrate hosts. MaRNAV2 was the first RNA virus discovered associated with parasites of the genus Leucocytozoon. By analyzing metatranscriptomes from the NCBI SRA database with local sequence alignment tools, we detected two novel RNA viruses; we describe them as MaRNAV3 associated with Leucocytozoon and MaRNAV4 associated with Parahaemoproteus. MaRNAV3 had ~59% amino acid identity to the RNA-dependent RNA-polymerase (RdRp) of MaRNAV1 and ~63% amino acid identity to MaRNAV2. MaRNAV4 had ~44% amino acid identity to MaRNAV1 and ~47% amino acid identity to MaRNAV2. These findings lead us to hypothesize that MaRNAV_like viruses are widespread and tightly associated with the order Haemosporida since they have been described in human Plasmodium vivax, and now two genera of avian haemosporidians.

RNA enigma: "From origin of life to novel Coronavirus- COVID-19"

Recent events of the viral catastrophe have shown the rapidity of spread of new disease through emergence of virulent strains. Proper control measures can be developed only through understanding the evolution of virulence in RNA viruses. To understand the evolution of this novel Coronavirus, COVID-19, it is imperative to delineate the evolution of RNA, its transformation into first life forms, the steady and continuous evolution and emergence through modification in their genome and nevertheless the natural selection. This review will throw light on these aspects to understand the possible origin of COVID-19 to control and eradicate this viral outbreak.

A World of Viruses Nested within Parasites: Unraveling Viral Diversity within Parasitic Flatworms (Platyhelminthes)

Because parasites have an inextricable relationship with their host, they have the potential to serve as viral reservoirs or facilitate virus host shifts. And yet, little is known about viruses infecting parasitic hosts except for blood-feeding arthropods that are well-known vectors of zoonotic viruses. Herein, we uncovered viruses of flatworms (phylum Platyhelminthes, group Neodermata) that specialize in parasitizing vertebrates and their ancestral free-living relatives. We discovered 115 novel viral sequences, including 1 in Macrostomorpha, 5 in Polycladida, 44 in Tricladida, 1 in Monogenea, 15 in Cestoda, and 49 in Trematoda, through data mining. The majority of newly identified viruses constitute novel families or genera. Phylogenetic analyses show that the virome of flatworms changed dramatically during the transition of neodermatans to a parasitic lifestyle. Most Neodermata viruses seem to codiversify with their host, with the exception of rhabdoviruses, which may switch hosts more often, based on phylogenetic relationships. Neodermata rhabdoviruses also have a position ancestral to vertebrate-associated rhabdo viruses, including lyssaviruses, suggesting that vertebrate-associated rhabdoviruses emerged from a flatworm rhabdovirus in a parasitized host. This study reveals an extensive diversity of viruses in Platyhelminthes and highlights the need to evaluate the role of viral infection in flatworm-associated diseases.

IMPORTANCE Little is known about the diversity of parasite-associated viruses and how these viruses may impact parasite fitness, parasite-host interactions, and virus evolution. The discovery of over a hundred viruses associated with a range of free-living and parasitic flatworms, including parasites of economic and clinical relevance, allowed us to compare the viromes of flatworms with contrasting lifestyles. The results suggest that flatworms acquired novel viruses after their transition to a parasitic lifestyle and highlight the possibility that they acquired viruses from their hosts and vice versa. An interesting example is the discovery of flatworm rhabdoviruses that have a position ancestral to rabies viruses and other vertebrate-associated rhabdoviruses, demonstrating that flatworm-associated viruses have emerged in a vertebrate host at least once in history. Therefore, parasitic flatworms may play a role in virus diversity and emergence. The roles that parasite-infecting viruses play in parasite-associated diseases remain to be investigated.

Genomic evolution of the Coronaviridae family

The current outbreak of coronavirus disease-2019 (COVID-19) caused by SARS-CoV-2 poses unparalleled challenges to global public health. SARS-CoV-2 is a Betacoronavirus, one of four genera belonging to the Coronaviridae subfamily Orthocoronavirinae. Coronaviridae, in turn, are members of the order Nidovirales, a group of enveloped, positive-stranded RNA viruses. Here we present a systematic phylogenetic and evolutionary study based on protein domain architecture, encompassing the entire proteomes of all Orthocoronavirinae, as well as other Nidovirales. This analysis has revealed that the genomic evolution of Nidovirales is associated with extensive gains and losses of protein domains. In Orthocoronavirinae, the sections of the genomes that show the largest divergence in protein domains are found in the proteins encoded in the amino-terminal end of the polyprotein (PP1ab), the spike protein (S), and many of the accessory proteins. The diversity among the accessory proteins is particularly striking, as each subgenus possesses a set of accessory proteins that is almost entirely specific to that subgenus. The only notable exception to this is ORF3b, which is present and orthologous over all Alphacoronaviruses. In contrast, the membrane protein (M), envelope small membrane protein (E), nucleoprotein (N), as well as proteins encoded in the central and carboxy-terminal end of PP1ab (such as the 3C-like protease, RNA-dependent RNA polymerase, and Helicase) show stable domain architectures across all Orthocoronavirinae. This comprehensive analysis of the Coronaviridae domain architecture has important implication for efforts to develop broadly cross-protective coronavirus vaccines.

Circulative Transmission of Cileviruses in Brevipalpus Mites May Involve the Paracellular Movement of Virions

Plant viruses transmitted by mites of the genus Brevipalpus are members of the genera Cilevirus, family Kitaviridae, or Dichorhavirus, family Rhabdoviridae. They produce non-systemic infections that typically display necrotic and/or chlorotic lesions around the inoculation loci. The cilevirus citrus leprosis virus C (CiLV-C) causes citrus leprosis, rated as one of the most destructive diseases affecting this crop in the Americas. CiLV-C is vectored in a persistent manner by the flat mite Brevipalpus yothersi. Upon the ingestion of viral particles with the content of the infected plant cell, virions must pass through the midgut epithelium and the anterior podocephalic gland of the mites. Following the duct from this gland, virions reach the salivary canal before their inoculation into a new plant cell through the stylet canal. It is still unclear whether CiLV-C multiplies in mite cells and what mechanisms contribute to its movement through mite tissues. In this study, based on direct observation of histological sections from viruliferous mites using the transmission electron microscope, we posit the hypothesis of the paracellular movement of CiLV-C in mites which may involve the manipulation of septate junctions. We detail the presence of viral particles aligned in the intercellular spaces between cells and the gastrovascular system of Brevipalpus mites. Accordingly, we propose putative genes that could control either active or passive paracellular circulation of viral particles inside the mites.

Six Novel Mycoviruses Containing Positive Single-Stranded RNA and Double-Stranded RNA Genomes Co-Infect a Single Strain of the Rhizoctoniasolani AG-3 PT

Six novel mycoviruses that collectively represent the mycovirome of Rhizoctonia solani anastomosis group (AG)-3 PT strain ZJ-2H, which causes potato black scurf, were identified through metatranscriptome sequencing and putatively designated as Rhizoctonia solani fusarivirus 4 [RsFV4, positive single-stranded RNA (+ssRNA)], Rhizoctonia solani fusarivirus 5 (RsFV5, +ssRNA), Rhizoctonia solani mitovirus 40 (RsMV40, +ssRNA), Rhizoctonia solani partitivirus 10 [RsPV10, double-stranded RNA (dsRNA)], Rhizoctonia solani partitivirus 11 (RsPV11, dsRNA), and Rhizoctonia solani RNA virus 11 (RsRV11, dsRNA). Whole genome sequences of RsFV4, RsMV40, RsPV10, RsPV11, and RsRV11, as well as a partial genome sequence of RsFV5, were obtained. The 3'- and 5'- untranslated regions of the five mycoviruses with complete genome sequences were folded into stable stem-loop or panhandle secondary structures. RsFV4 and RsFV5 are most closely related to Rhizoctonia solani fusarivirus 1 (RsFV1), however, the first open reading frame (ORF) of RsFV4 and RsFV5 encode a hypothetical protein that differs from the first ORF of RsFV1, which encodes a helicase. We confirmed that RsPV10 and RsPV11 assemble into the spherical virus particles (approximately 30 nm in diameter) that were extracted from strain ZJ-2H. This is the first report that +ssRNA and dsRNA viruses co-infect a single strain of R. solani AG-3 PT.

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.

Intrinsic variation in the vertically transmitted core virome of the mosquito Aedes aegypti

Virome studies among metazoans have revealed the ubiquity of RNA viruses in animals, contributing to a fundamental re-thinking of the relationships between organisms and their microbiota. Mosquito viromes, often scrutinized due to their public health relevance, may also provide insight into broadly applicable concepts, such as a "core virome," a set of viruses consistently associated with a host species or population that may fundamentally impact its basic biology. A subset of mosquito-associated viruses (MAVs) could comprise such a core, and MAVs can be categorized as (i) arboviruses, which alternate between mosquito and vertebrate hosts, (ii) insect-specific viruses, which cannot replicate in vertebrate cells, and (iii) viruses with unknown specificity. MAVs have been widely characterized in the disease vector Aedes aegypti, and the occurrence of a core virome in this species has been proposed but remains unclear. Using a wild population previously surveyed for MAVs and a common laboratory strain, we investigated viromes in reproductive tissue via metagenomic RNA sequencing. Virome composition varied across samples, but four groups comprised >97% of virus sequences: a novel partiti-like virus (Partitiviridae), a toti-like virus (Totiviridae), unclassified Riboviria, and four orthomyxo-like viruses (Orthormyxoviridae). Whole or partial genomes for the partiti-like virus, toti-like virus, and one orthomyxo-like virus were assembled and analyzed phylogenetically. Multigenerational maintenance of these MAVs was confirmed by RT-PCR, indicating vertical transmission as a mechanism for persistence. This study provides fundamental information regarding MAV ecology and variability in A. aegypti and the potential for vertically maintained core viromes at the population level.

Maf/ham1-like pyrophosphatases of non-canonical nucleotides are host-specific partners of viral RNA-dependent RNA polymerases

Cassava brown streak disease (CBSD), dubbed the “Ebola of plants”, is a serious threat to food security in Africa caused by two viruses of the family Potyviridae: cassava brown streak virus (CBSV) and Ugandan (U)CBSV. Intriguingly, U/CBSV, along with another member of this family and one secoviridae, are the only known RNA viruses encoding a protein of the Maf/ham1-like family, a group of widespread pyrophosphatase of non-canonical nucleotides (ITPase) expressed by all living organisms. Despite the socio-economic impact of CDSD, the relevance and role of this atypical viral factor has not been yet established. Here, using an infectious cDNA clone and reverse genetics, we demonstrate that UCBSV requires the ITPase activity for infectivity in cassava, but not in the model plant Nicotiana benthamiana. HPLC-MS/MS experiments showed that, quite likely, this host-specific constraint is due to an unexpected high concentration of non-canonical nucleotides in cassava. Finally, protein analyses and experimental evolution of mutant viruses indicated that keeping a fraction of the yielded UCBSV ITPase covalently bound to the viral RNA-dependent RNA polymerase (RdRP) optimizes viral fitness, and this seems to be a feature shared by the other members of the Potyviridae family expressing Maf/ham1-like proteins. All in all, our work (i) reveals that the over-accumulation of non-canonical nucleotides in the host might have a key role in antiviral defense, and (ii) provides the first example of an RdRP-ITPase partnership, reinforcing the idea that RNA viruses are incredibly versatile at adaptation to different host setups.

Cassava is one the most important staple food around the world in term of caloric intake. The cassava brown streak disease, caused by cassava brown streak virus (CBSV) and Ugandan (U)CBSV–Ipomovirus genus, Potyviridae family-, produces massive losses in cassava production. Curiously, these two viruses, unlike the vast majority of members of the family, encode a Maf1/ham1-like pyrophosphatase (HAM1) of non-canonical nucleotides with unknown relevance and function in viruses. This study aims to fill this gap in our knowledge by using reverse genetics, biochemistry, metabolomics and directed virus evolution. Hence, we found that HAM1 is required for UCBSV to infect cassava, where its pyrophosphatase activity resulted critical, but not to propagate in the model plant Nicotiana benthamiana. In addition, we demonstrated that HAM1 works in partnership with the viral RdRP during infection. Unexpected high levels of ITP/XTP non-canonical nucleotides found in cassava, and the known flexibility of RNA viruses to incorporate additional factors when required, supports the idea that the high concentration of ITP/XTP worked as a selection pressure to promote the acquisition of HAM1 into the virus in order to promote a successful infection.