A member of the virus family Narnaviridae from the plant pathogenic oomycete Phytophthora infestans

Discovery and Identification of Viruses Infecting Oomycetes

This chapter describes protocols suitable for the detection and identification of RNA viruses infecting oomycetes (so-called water molds of Kingdom Heterokonta, Stramenopila), focusing on species of Phytophthora and exemplified by P. fragariae. The protocol includes laboratory procedures for oomycete cultivation and total RNA extraction from harvested mycelia, followed by instructions on suitable parameters given for sequencing companies on ribosomal RNA depletion, cDNA library preparation, and total RNA-sequencing (RNA-Seq). We also describe the bioinformatics steps needed for de novo assembly of raw reads into contigs, removal of host-associated contigs, and virus identification by database searches, as well as host validation by RT-PCR. All steps are described using an exemplar RNA-Seq library containing a yet undescribed fusagravirus hosted by a P. fragariae isolate.

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.

Bunyaviruses Affect Growth, Sporulation, and Elicitin Production in Phytophthora cactorum

Phytophthora cactorum is an important oomycetous plant pathogen with numerous host plant species, including garden strawberry (Fragaria × ananassa) and silver birch (Betula pendula). P. cactorum also hosts mycoviruses, but their phenotypic effects on the host oomycete have not been studied earlier. In the present study, we tested polyethylene glycol (PEG)-induced water stress for virus curing and created an isogenic virus-free isolate for testing viral effects in pair with the original isolate. Phytophthora cactorum bunya-like viruses 1 and 2 (PcBV1 & 2) significantly reduced hyphal growth of the P. cactorum host isolate, as well as sporangia production and size. Transcriptomic and proteomic analyses revealed an increase in the production of elicitins due to bunyavirus infection. However, the presence of bunyaviruses did not seem to alter the pathogenicity of P. cactorum. Virus transmission through anastomosis was unsuccessful in vitro.

High Diversity of Novel Viruses in the Tree Pathogen Phytophthora castaneae Revealed by High-Throughput Sequencing of Total and Small RNA

Phytophthora castaneae, an oomycete pathogen causing root and trunk rot of different tree species in Asia, was shown to harbor a rich diversity of novel viruses from different families. Four P. castaneae isolates collected from Chamaecyparis hodginsii in a semi-natural montane forest site in Vietnam were investigated for viral presence by traditional and next-generation sequencing (NGS) techniques, i.e., double-stranded RNA (dsRNA) extraction and high-throughput sequencing (HTS) of small RNAs (sRNAs) and total RNA. Genome organization, sequence similarity, and phylogenetic analyses indicated that the viruses were related to members of the order Bunyavirales and families Endornaviridae, Megabirnaviridae, Narnaviridae, Totiviridae, and the proposed family "Fusagraviridae." The study describes six novel viruses: Phytophthora castaneae RNA virus 1-5 (PcaRV1-5) and Phytophthora castaneae negative-stranded RNA virus 1 (PcaNSRV1). All six viruses were detected by sRNA sequencing, which demonstrates an active RNA interference (RNAi) system targeting viruses in P. castaneae. To our knowledge, this is the first report of viruses in P. castaneae and the whole Phytophthora major Clade 5, as well as of the activity of an RNAi mechanism targeting viral genomes among Clade 5 species. PcaRV1 is the first megabirnavirus described in oomycetes and the genus Phytophthora.

Characterization of a Novel Double-Stranded RNA Virus from Phytophthora pluvialis in New Zealand

A new dsRNA virus from the oomycete Phytophthora pluvialis has been characterized and designated as Phytophthora pluvialis RNA virus 1 (PplRV1). The genome of the PplRV1 reference genome is 6742 bp that encodes two predicted open reading frames (ORFs). ORF1 and ORF2 overlap by a 47 nt “slippery” frameshift sequence. ORF1 encodes a putative protein of unknown function. ORF2 shows high similarity to the RNA-dependent RNA polymerase (RdRp) of other dsRNA viruses. Phylogenetic analysis of the putative PplRV1 RdRp and its most closely related viruses showed PplRV1 is distinct from other known viruses (below 33% amino acid similarity), which indicates this virus may belong to a new virus family. Analyses of the geographical distribution of PplRV1 in relation to two genetically distinct classes of its host revealed two corresponding genotypes of the PplRV1 (termed a and b), which share 92.3% nt identity. The reference genome for the second genotype is 6760 bp long and a prediction of its genetic organization shows three ORFs, with ORF2 being split into two ORFs, ORF2a and ORF2b, that is conserved in seven of eleven genotype b isolates. Additionally, a quick and simple diagnostic method using qPCR has been developed, which is suitable for large scale screens to identify PplRV1 in Phytophthora.

Novel Fusari- and Toti-like Viruses, with Probable Different Origins, in the Plant Pathogenic Oomycete Globisporangiumultimum

To further classify the oomycete viruses that have been discovered in recent years, we investigated virus infection in the plant-parasitic oomycete Globisporangium ultimum in Japan. Double-stranded RNA detection, high-throughput sequencing, and RT-PCR revealed that the G. ultimum isolate UOP226 contained two viruses related to fusarivirus and totivirus, named Pythium ultimum RNA virus 1 (PuRV1) and Pythium ultimum RNA virus 2 (PuRV2), respectively. Phylogenetic analysis of the deduced amino acid sequence of the RNA-dependent RNA polymerase (RdRp) showed that fusari-like PuRV1 belonged to a different phylogenetic group than Plasmopara viticola lesion-associated fusari virus (PvlaFV) 1-3 from oomycete Plasmopara viticola. Codon usage bias of the PuRV1 RdRp gene was more similar to those of fungi than Globisporangium and Phytophthora, suggesting that the PuRV1 ancestor horizontally transmitted to G. ultimum ancestor from fungi. Phylogenetic analysis of the deduced amino acid sequence of the RdRp of toti-like PuRV2 showed a monophyletic group with the other toti-like oomycete viruses from Globisporangium, Phytophthora, and Pl. viticola. However, the nucleotide sequences of toti-like oomycete viruses were not so homologous, suggesting the possibility of convergent evolution of toti-like oomycete viruses.

Characterization of the Mycovirome from the Plant-Pathogenic Fungus Cercospora beticola

Cercospora leaf spot (CLS) caused by Cercospora beticola is a devastating foliar disease of sugar beet (Beta vulgaris), resulting in high yield losses worldwide. Mycoviruses are widespread fungi viruses and can be used as a potential biocontrol agent for fugal disease management. To determine the presence of mycoviruses in C. beticola, high-throughput sequencing analysis was used to determine the diversity of mycoviruses in 139 C. beticola isolates collected from major sugar beet production areas in China. The high-throughput sequencing reads were assembled and searched against the NCBI database using BLASTn and BLASTx. The results showed that the obtained 93 contigs were derived from eight novel mycoviruses, which were grouped into 3 distinct lineages, belonging to the families Hypoviridae, Narnaviridae and Botourmiaviridae, as well as some unclassified (−)ssRNA viruses in the order Bunyavirales and Mononegavirales. To the best of our knowledge, this is the first identification of highly diverse mycoviruses in C. beticola. The novel mycoviruses explored in this study will provide new viral materials to biocontrol Cercospora diseases. Future studies of these mycoviruses will aim to assess the roles of each mycovirus in biological function of C. beticola in the future.

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.

Viral diversity in Phytophthora cactorum population infecting strawberry

Eighty-eight Phytophthora cactorum strains isolated from crown or leather rot of strawberry in 1971-2019 were screened for viruses using RNA-seq and RT-PCR. Remarkably, all but one isolate were virus-infected, most of them harbouring more than one virus of different genera or species. The most common virus occurring in 94% of the isolates was the Phytophthora cactorum RNA virus 1 (PcRV1) resembling members of Totiviridae. Novel viruses related to members of Endornaviridae, named Phytophthora cactorum alphaendornaviruses 1-3 (PcAEV1-3), were found in 57% of the isolates. Four isolates hosted viruses with affinities to Bunyaviridae, named Phytophthora cactorum bunyaviruses 1-3 (PcBV1-3), and a virus resembling members of the proposed genus 'Ustivirus', named Phytophthora cactorum usti-like virus (PcUV1), was found in a single isolate. Most of the virus species were represented by several distinct strains sharing ≥81.4% aa sequence identity. We found no evidence of spatial differentiation but some temporal changes in the P. cactorum virus community were observed. Some isolates harboured two or more closely related strains of the same virus (PcAEV1 or PcRV1) sharing 86.6%-96.4% nt identity in their polymerase sequence. This was surprising as viruses with such a high similarity are typically mutually exclusive.

Multiple Viral Infections Detected in Phytophthora condilina by Total and Small RNA Sequencing

Marine oomycetes have recently been shown to be concurrently infected by (−)ssRNA viruses of the order Bunyavirales. In this work, even higher virus variability was found in a single isolate of Phytophthora condilina, a recently described member of Phytophthora phylogenetic Clade 6a, which was isolated from brackish estuarine waters in southern Portugal. Using total and small RNA-seq the full RdRp of 13 different potential novel bunya-like viruses and two complete toti-like viruses were detected. All these viruses were successfully confirmed by reverse transcription polymerase chain reaction (RT-PCR) using total RNA as template, but complementarily one of the toti-like and five of the bunya-like viruses were confirmed when dsRNA was purified for RT-PCR. In our study, total RNA-seq was by far more efficient for de novo assembling of the virus sequencing but small RNA-seq showed higher read numbers for most viruses. Two main populations of small RNAs (21 nts and 25 nts-long) were identified, which were in accordance with other Phytophthora species. To the best of our knowledge, this is the first study using small RNA sequencing to identify viruses in Phytophthora spp.

Two Novel Endornaviruses Co-infecting a Phytophthora Pathogen of Asparagus officinalis Modulate the Developmental Stages and Fungicide Sensitivities of the Host Oomycete

Two novel endornaviruses, Phytophthora endornavirus 2 (PEV2) and Phytophthora endornavirus 3 (PEV3) were found in isolates of a Phytophthora pathogen of asparagus collected in Japan. A molecular phylogenetic analysis indicated that PEV2 and PEV3 belong to the genus Alphaendornavirus. The PEV2 and PEV3 genomes consist of 14,345 and 13,810 bp, and they contain single open reading frames of 4,640 and 4,603 codons, respectively. Their polyproteins contain the conserved domains of an RNA helicase, a UDP-glycosyltransferase, and an RNA-dependent RNA polymerase, which are conserved in other alphaendornaviruses. PEV2 is closely related to Brown algae endornavirus 2, whereas PEV3 is closely related to Phytophthora endornavirus 1 (PEV1), which infects a Phytophthora sp. specific to Douglas fir. PEV2 and PEV3 were detected at high titers in two original Phytophthora sp. isolates, and we found a sub-isolate with low titers of the viruses during subculture. We used the high- and low-titer isolates to evaluate the effects of the viruses on the growth, development, and fungicide sensitivities of the Phytophthora sp. host. The high-titer isolates produced smaller mycelial colonies and much higher numbers of zoosporangia than the low-titer isolate. These results suggest that PEV2 and PEV3 inhibited hyphal growth and stimulated zoosporangium formation. The high-titer isolates were more sensitive than the low-titer isolate to the fungicides benthiavalicarb-isopropyl, famoxadone, and chlorothalonil. In contrast, the high-titer isolates displayed lower sensitivity to the fungicide metalaxyl (an inhibitor of RNA polymerase I) when compared with the low-titer isolate. These results indicate that persistent infection with PEV2 and PEV3 may potentially affect the fungicide sensitivities of the host oomycete.

A Novel Ourmia-Like Mycovirus Confers Hypovirulence-Associated Traits on Fusarium oxysporum

Fusarium wilt caused by Fusarium oxysporum f. sp. momordicae (FoM) is an important fungal disease that affects the production of bitter gourd. Hypovirulence-associated mycoviruses have great potential and application prospects for controlling the fungal disease. In this study, a novel ourmia-like virus, named Fusarium oxysporum ourmia-like virus 1 (FoOuLV1), was isolated from FoM strain HuN8. The viral genomic RNA is 2,712 nucleotides (nt) in length and contains an open reading frame (ORF) encoding a putative RNA-dependent RNA polymerase (RdRp) using either standard or mitochondrial codes. In strain HuN8, there was also a FoOuLV1-associated RNA segment with 1,173 nt in length with no sequence homology. Phylogenetic analysis showed that FoOuLV1 is a member of the genus Magoulivirus of the family Botourmiaviridae. FoOuLV1 was found to be associated with hypovirulence in FoM. Moreover, FoOuLV1 and its hypovirulence trait can be transmitted horizontally to other FoM strains and also to other formae speciale strains of F. oxysporum. In addition, FoOuLV1 showed significant biological control effect against the bitter gourd Fusarium wilt. To our knowledge, this study reveals the first description of a hypovirulence-associated ourmia-like mycovirus, which has the potential to the biological control of Fusarium wilt.

Analysis of the virome associated to grapevine downy mildew lesions reveals new mycovirus lineages

The obligatory biotrophic oomycetes Plasmopara viticola is the causal agent of downy mildew, a destructive disease of grapevine worldwide. So far, chemical fungicides are widely employed to limit this pathogen, but their adverse effects are stimulating the quest for environmentally friendly alternative approaches. Here, we report on the search for mycoviruses that might be later developed as biocontrol agents for this pathogen. Symptomatic leaves were collected from various regions in Spain and Italy and mycelia associated to leaf lesions was harvested. Total RNA extractions were depleted of rRNA and metatranscriptomes were generated using a high-throughput sequencing approach. The virome associated to leaf lesions was then characterized through a bioinformatic pipeline relying on blast searches against current viral databases. Here, we present an inventory of 283 new RNA viruses: 222 positive strand RNA viruses, 29 negative strand RNA viruses, 27 double-stranded RNA viruses and 5 ORFan virus RdRP segments, which could not be reliably assigned to any existing group in the Riboviria. In addition to ORFan viruses, we found other surprising new evolutionary trajectories in this wide inventory of viruses. The most represented viruses in our collection are those in phylum Lenarviricota, and, among them, a group of mycovirus segments distantly related to narnaviruses, but characterized by a polymerase palm domain lacking subdomain C, with the putative GDD catalytic triad. We also provided evidence of a strict association between two RNA segments that form a new mycovirus clade of positive strand RNA in the phylum Kitrinoviricota, order Martellivirales. In the phylum Negarnaviricota, we report for the first time in the order Mononegavirales a clade of viruses that is ambisense, a feature that so far was present only in the order Bunyavirales. Furthermore, in the same phylum we detected the widespread occurrence and abundant accumulation in our libraries of a distinct mycovirus clade distantly related to the Muvirales and Goujanvirales orders, which so far include only viruses infecting invertebrates. Possible new oomycetes-specific virus clades are also described in the phylum Duplornaviricota. These data greatly expand the evolutionary history of mycoviruses adding new layers of diversity to the realm Riboviria.

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.

Marine Oomycetes of the Genus Halophytophthora Harbor Viruses Related to Bunyaviruses

We investigated the incidence of RNA viruses in a collection of Halophytophthora spp. from estuarine ecosystems in southern Portugal. The first approach to detect the presence of viruses was based on the occurrence of dsRNA, typically considered as a viral molecule in plants and fungi. Two dsRNA-banding patterns (∼7 and 9 kb) were observed in seven of 73 Halophytophthora isolates tested (9.6%). Consequently, two dsRNA-hosting isolates were chosen to perform stranded RNA sequencing for de novo virus sequence assembly. A total of eight putative novel virus species with genomic affinities to members of the order Bunyavirales were detected and their full-length RdRp gene characterized by RACE. Based on the direct partial amplification of their RdRp gene by RT-PCR multiple viral infections occur in both isolates selected. Likewise, the screening of those viruses in the whole collection of Halophytophthora isolates showed that their occurrence is limited to one single Halophytophthora species. To our knowledge, this is the first report demonstrating the presence of negative (−) ssRNA viruses in marine oomycetes.

Narnaviruses: novel players in fungal-bacterial symbioses

Rhizopus microsporus is an early-diverging fungal species with importance in ecology, agriculture, food production, and public health. Pathogenic strains of R. microsporus harbor an intracellular bacterial symbiont, Mycetohabitans (formerly named Burkholderia). This vertically transmitted bacterial symbiont is responsible for the production of toxins crucial to the pathogenicity of Rhizopus and remarkably also for fungal reproduction. Here we show that R. microsporus can live not only in symbiosis with bacteria but also with two viral members of the genus Narnavirus. Our experiments revealed that both viruses replicated similarly in the growth conditions we tested. Viral copies were affected by the developmental stage of the fungus, the substrate, and the presence or absence of Mycetohabitans. Absolute quantification of narnaviruses in isolated asexual sporangiospores and sexual zygospores indicates their vertical transmission. By curing R. microsporus of its viral and bacterial symbionts and reinfecting bacteria to reestablish symbiosis, we demonstrate that these viruses affect fungal biology. Narnaviruses decrease asexual reproduction, but together with Mycetohabitans, are required for sexual reproductive success. This fungal-bacterial-viral system represents an outstanding model to investigate three-way microbial symbioses and their evolution.

A case for a negative-strand coding sequence in a group of positive-sense RNA viruses

Positive-sense single-stranded RNA viruses form the largest and most diverse group of eukaryote-infecting viruses. Their genomes comprise one or more segments of coding-sense RNA that function directly as messenger RNAs upon release into the cytoplasm of infected cells. Positive-sense RNA viruses are generally accepted to encode proteins solely on the positive strand. However, we previously identified a surprisingly long (∼1,000-codon) open reading frame (ORF) on the negative strand of some members of the family Narnaviridae which, together with RNA bacteriophages of the family Leviviridae, form a sister group to all other positive-sense RNA viruses. Here, we completed the genomes of three mosquito-associated narnaviruses, all of which have the long reverse-frame ORF. We systematically identified narnaviral sequences in public data sets from a wide range of sources, including arthropod, fungal, and plant transcriptomic data sets. Long reverse-frame ORFs are widespread in one clade of narnaviruses, where they frequently occupy >95 per cent of the genome. The reverse-frame ORFs correspond to a specific avoidance of CUA, UUA, and UCA codons (i.e. stop codon reverse complements) in the forward-frame RNA-dependent RNA polymerase ORF. However, absence of these codons cannot be explained by other factors such as inability to decode these codons or GC3 bias. Together with other analyses, we provide the strongest evidence yet of coding capacity on the negative strand of a positive-sense RNA virus. As these ORFs comprise some of the longest known overlapping genes, their study may be of broad relevance to understanding overlapping gene evolution and de novo origin of genes.

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

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

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

An exploration of ambigrammatic sequences in narnaviruses

Narnaviruses have been described as positive-sense RNA viruses with a remarkably simple genome of ~3 kb, encoding only a highly conserved RNA-dependent RNA polymerase (RdRp). Many narnaviruses, however, are 'ambigrammatic' and harbour an additional uninterrupted open reading frame (ORF) covering almost the entire length of the reverse complement strand. No function has been described for this ORF, yet the absence of stops is conserved across diverse narnaviruses, and in every case the codons in the reverse ORF and the RdRp are aligned. The >3 kb ORF overlap on opposite strands, unprecedented among RNA viruses, motivates an exploration of the constraints imposed or alleviated by the codon alignment. Here, we show that only when the codon frames are aligned can all stop codons be eliminated from the reverse strand by synonymous single-nucleotide substitutions in the RdRp gene, suggesting a mechanism for de novo gene creation within a strongly conserved amino-acid sequence. It will be fascinating to explore what implications this coding strategy has for other aspects of narnavirus biology. Beyond narnaviruses, our rapidly expanding catalogue of viral diversity may yet reveal additional examples of this broadly-extensible principle for ambigrammatic-sequence development.

Viruses of fungi and oomycetes in the soil environment

Soils support a myriad of organisms hosting highly diverse viromes. In this minireview, we focus on viruses hosted by true fungi and oomycetes (members of Stamenopila, Chromalveolata) inhabiting bulk soil, rhizosphere and litter layer, and representing different ecological guilds, including fungal saprotrophs, mycorrhizal fungi, mutualistic endophytes and pathogens. Viruses infecting fungi and oomycetes are characterized by persistent intracellular nonlytic lifestyles and transmission via spores and/or hyphal contacts. Almost all fungal and oomycete viruses have genomes composed of single-stranded or double-stranded RNA, and recent studies have revealed numerous novel viruses representing yet unclassified family-level groups. Depending on the virus–host combination, infections can be asymptomatic, beneficial or detrimental to the host. Thus, mycovirus infections may contribute to the multiplex interactions of hosts, therefore likely affecting the dynamics of fungal communities required for the functioning of soil ecosystems. However, the effects of fungal and oomycete viruses on soil ecological processes are still mostly unknown. Interestingly, new metagenomics data suggest an extensive level of horizontal virus transfer between plants, fungi and insects.

Characterization of a Novel Ourmia-Like Mycovirus Infecting Magnaporthe oryzae and Implications for Viral Diversity and Evolution

Here, the molecular characterization of a novel mycovirus that was isolated from a phytopathogenic fungus Magnaporthe oryzae and designed as Magnaporthe oryzae ourmia-like virus 4 (MOLV4) is reported. MOLV4 has a genome that is 2497 bp long and possesses a single open reading frame (ORF), which encodes the product RNA-dependent RNA polymerase (RdRp). Sequence similarities were found between the MOLV4 encoded RdRp and the counterparts of a few previously reported ourmia-like mycoviruses. Virus-curing and biological comparison indicate that the virus has no or mild effects on the morphology and mycelium growth rate of the host fungus. Phylogenetic analysis using the RdRp aa sequences was performed. The results show that MOLV4 is clustered with the ourmia-like mycoviruses, forming a clade closely related to ourmiaviruses but distinct from narnaviruses. In addition, database searches revealed that several MOLV4-related sequences are present in the transcriptome shotgun assembly (TSA) library, expressed sequence tag database (ESTdb), whole-genome shotgun (WGS) library, and genomic survey sequences (GSS) libraries of a few other species of eukaryote organisms. Our results show that MOLV4, together with other similar ourmia-like mycoviruses, might represent a virus clade that links the plant ourmiaviruses and fungal narnaviruses and has a wide range of hosts.

RNA Viruses in Blechomonas (Trypanosomatidae) and Evolution of Leishmaniavirus

In this work, we analyzed viral prevalence in trypanosomatid parasites (Blechomonas spp.) infecting Siphonaptera and discovered nine species of viruses from three different groups (leishbunyaviruses, narnaviruses, and leishmaniaviruses). Most of the flagellate isolates bore two or three viral types (mixed infections). Although no new viral groups were documented in Blechomonas spp., our findings are important for the comprehension of viral evolution. The discovery of bunyaviruses in blechomonads was anticipated, since these viruses have envelopes facilitating their interspecific transmission and have already been found in various trypanosomatids and metatranscriptomes with trypanosomatid signatures. In this work, we also provided evidence that even representatives of the family Narnaviridae are capable of host switching and evidently have accomplished switches multiple times in the course of their evolution. The most unexpected finding was the presence of leishmaniaviruses, a group previously solely confined to the human pathogens Leishmania spp. From phylogenetic inferences and analyses of the life cycles of Leishmania and Blechomonas, we concluded that a common ancestor of leishmaniaviruses most likely infected Leishmania first and was acquired by Blechomonas by horizontal transfer. Our findings demonstrate that evolution of leishmaniaviruses is more complex than previously thought and includes occasional host switching.IMPORTANCE Flagellates belonging to the genus Leishmania are important human parasites. Some strains of different Leishmania species harbor viruses (leishmaniaviruses), which facilitate metastatic spread of the parasites, thus aggravating the disease. Up until now, these viruses were known to be hosted only by Leishmania Here, we analyzed viral distribution in Blechomonas, a related group of flagellates parasitizing fleas, and revealed that they also bear leishmaniaviruses. Our findings shed light on the entangled evolution of these viruses. In addition, we documented that Blechomonas can be also infected by leishbunyaviruses and narnaviruses, viral groups known from other insects' flagellates.

A Novel Hypovirus Species From Xylariaceae Fungi Infecting Avocado

The white rot root disease caused by Rosellinia necatrix is a major concern for avocado cultivation in Spain. Healthy escapes of avocado trees surrounded by diseased trees prompted us to hypothesize the presence of hypovirulent R. necatrix due to mycovirus infections. Recently, we reported the presence of another fungal species, Entoleuca sp., belonging to the Xylariaceae, that was also found in healthy avocado trees and frequently co-infecting the same roots than R. necatrix. We investigated the presence of mycoviruses that might explain the hypovirulence. For that, we performed deep sequencing of dsRNAs from two isolates of Entoleuca sp. that revealed the simultaneous infection of several mycoviruses, not described previously. In this work, we report a new member of the Hypoviridae, tentatively named Entoleuca hypovirus 1 (EnHV1). The complete genome sequence was obtained for two EnHV1 strains, which lengths resulted to be 14,958 and 14,984 nt, respectively, excluding the poly(A) tails. The genome shows two ORFs separated by a 32-nt inter-ORF, and both 5′- and 3′-UTRs longer than any other hypovirus reported to date. The analysis of virus-derived siRNA populations obtained from Entoleuca sp. demonstrated antiviral silencing activity in this fungus. We screened a collection of Entoleuca sp. and R. necatrix isolates and found that EnHV1 was present in both fungal species. A genetic population analysis of EnHV1 strains revealed the presence of two main clades, each of them including members from both Entoleuca sp. and R. necatrix, which suggests intra- and interspecific virus transmission in the field. Several attempts failed to cure Entoleuca sp. from EnHV1. However, all Entoleuca sp. isolates collected from avocado, whether harboring the virus or not, showed hypovirulence. Conversely, all R. necatrix isolates were pathogenic to that crop, regardless of being infected by EnHV1.

Viral discovery and diversity in trypanosomatid protozoa with a focus on relatives of the human parasite Leishmania

Knowledge of viral diversity is expanding greatly, but many lineages remain underexplored. We surveyed RNA viruses in 52 cultured monoxenous relatives of the human parasite Leishmania (Crithidia and Leptomonas), as well as plant-infecting PhytomonasLeptomonas pyrrhocoris was a hotbed for viral discovery, carrying a virus (Leptomonas pyrrhocoris ostravirus 1) with a highly divergent RNA-dependent RNA polymerase missed by conventional BLAST searches, an emergent clade of tombus-like viruses, and an example of viral endogenization. A deep-branching clade of trypanosomatid narnaviruses was found, notable as Leptomonas seymouri bearing Narna-like virus 1 (LepseyNLV1) have been reported in cultures recovered from patients with visceral leishmaniasis. A deep-branching trypanosomatid viral lineage showing strong affinities to bunyaviruses was termed "Leishbunyavirus" (LBV) and judged sufficiently distinct to warrant assignment within a proposed family termed "Leishbunyaviridae" Numerous relatives of trypanosomatid viruses were found in insect metatranscriptomic surveys, which likely arise from trypanosomatid microbiota. Despite extensive sampling we found no relatives of the totivirus Leishmaniavirus (LRV1/2), implying that it was acquired at about the same time the Leishmania became able to parasitize vertebrates. As viruses were found in over a quarter of isolates tested, many more are likely to be found in the >600 unsurveyed trypanosomatid species. Viral loss was occasionally observed in culture, providing potentially isogenic virus-free lines enabling studies probing the biological role of trypanosomatid viruses. These data shed important insights on the emergence of viruses within an important trypanosomatid clade relevant to human disease.

Phomopsis longicolla RNA virus 1 - Novel virus at the edge of myco- and plant viruses

The complete nucleotide sequence of a new RNA mycovirus in the KY isolate of Phomopsis longicolla Hobbs 1985 and its protoplasts subcultures p5, p9, and ME711 was discovered. The virus, provisionally named Phomopsis longicolla RNA virus 1 (PlRV1), was localized in mitochondria and was determined to have a genome 2822 nucleotides long. A single open reading frame could be translated in silico by both standard and mitochondrial genetic codes into a product featuring conservative domains for an RNA-dependent RNA polymerase (RdRp). The RdRp of PlRV1 has no counterpart among mycoviruses, but it is about 30% identical with the RdRp of plant ourmiaviruses. Recently, new mycoviruses related to plant ourmiaviruses and forming one clade with PlRV1 have been discovered. This separate clade could represent the crucial link between plant and fungal viruses.

Studies on the Virome of the Entomopathogenic Fungus Beauveria bassiana Reveal Novel dsRNA Elements and Mild Hypervirulence

The entomopathogenic fungus Beauveria bassiana has a wide host range and is used as a biocontrol agent against arthropod pests. Mycoviruses have been described in phytopathogenic fungi while in entomopathogenic fungi their presence has been reported only rarely. Here we show that 21.3% of a collection of B. bassiana isolates sourced from worldwide locations, harbor dsRNA elements. Molecular characterization of these elements revealed the prevalence of mycoviruses belonging to the Partitiviridae and Totiviridae families, the smallest reported virus to date, belonging to the family Narnaviridae, and viruses unassigned to a family or genus. Of particular importance is the discovery of members of a newly proposed family Polymycoviridae in B. bassiana. Polymycoviruses, previously designated as tetramycoviruses, consist of four non-conventionally encapsidated capped dsRNAs. The presence of additional non-homologous genomic segments in B. bassiana polymycoviruses and other fungi illustrates the unprecedented dynamic nature of the viral genome. Finally, a comparison of virus-free and virus-infected isogenic lines derived from an exemplar B. bassiana isolate revealed a mild hypervirulent effect of mycoviruses on the growth of their host isolate and on its pathogenicity against the greater wax moth Galleria mellonella, highlighting for the first time the potential of mycoviruses as enhancers of biocontrol agents.

A Narnavirus in the Trypanosomatid Protist Plant Pathogen Phytomonas serpens

We describe here a new RNA virus (PserNV1) from the plant protist parasite Phytomonas serpens (family Trypanosomatidae, Kinetoplastida, supergroup Excavata). The properties of PserNV1 permit assignment to the genus Narnavirus (Narnaviridae), the first reported from a host other than fungi or oomycetes.

Molecular characterization of Botrytis ourmia-like virus, a mycovirus close to the plant pathogenic genus Ourmiavirus

The molecular characterization of a novel single-stranded RNA virus, obtained by next generation sequencing using Illumina platform, in a field grapevine isolate of the plant pathogenic fungus Botrytis, is reported in this work. The sequence comparison of this virus against the NCBI database showed a strong identity with RNA dependent RNA polymerases (RdRps) of plant pathogenic viruses belonging to the genus Ourmiavirus, therefore, this novel virus was named Botrytis ourmia-like virus (BOLV). BOLV has one open reading frame of 2169 nucleotides, which encodes a protein of 722 amino acids showing conserved domains of plant RNA viruses RdRps such as the most conserved GDD active domain. Our analyses showed that BOLV is phylogenetically closer to the fungal Narnavirus and the plant Ourmiavirus than to Mitovirus of the family Narnaviridae. Hence, we proposed that BOLV might represent the link between fungal viruses of the family Narnaviridae and the plant ourmiaviruses.

Widespread mitovirus sequences in plant genomes

The exploration of the evolution of RNA viruses has been aided recently by the discovery of copies of fragments or complete genomes of non-retroviral RNA viruses (Non-retroviral Endogenous RNA Viral Elements, or NERVEs) in many eukaryotic nuclear genomes. Among the most prominent NERVEs are partial copies of the RNA dependent RNA polymerase (RdRP) of the mitoviruses in plant mitochondrial genomes. Mitoviruses are in the family Narnaviridae, which are the simplest viruses, encoding only a single protein (the RdRP) in their unencapsidated viral plus strand. Narnaviruses are known only in fungi, and the origin of plant mitochondrial mitovirus NERVEs appears to be horizontal transfer from plant pathogenic fungi. At least one mitochondrial mitovirus NERVE, but not its nuclear copy, is expressed.

Virus world as an evolutionary network of viruses and capsidless selfish elements

Viruses were defined as one of the two principal types of organisms in the biosphere, namely, as capsid-encoding organisms in contrast to ribosome-encoding organisms, i.e., all cellular life forms. Structurally similar, apparently homologous capsids are present in a huge variety of icosahedral viruses that infect bacteria, archaea, and eukaryotes. These findings prompted the concept of the capsid as the virus "self" that defines the identity of deep, ancient viral lineages. However, several other widespread viral "hallmark genes" encode key components of the viral replication apparatus (such as polymerases and helicases) and combine with different capsid proteins, given the inherently modular character of viral evolution. Furthermore, diverse, widespread, capsidless selfish genetic elements, such as plasmids and various types of transposons, share hallmark genes with viruses. Viruses appear to have evolved from capsidless selfish elements, and vice versa, on multiple occasions during evolution. At the earliest, precellular stage of life's evolution, capsidless genetic parasites most likely emerged first and subsequently gave rise to different classes of viruses. In this review, we develop the concept of a greater virus world which forms an evolutionary network that is held together by shared conserved genes and includes both bona fide capsid-encoding viruses and different classes of capsidless replicons. Theoretical studies indicate that selfish replicons (genetic parasites) inevitably emerge in any sufficiently complex evolving ensemble of replicators. Therefore, the key signature of the greater virus world is not the presence of a capsid but rather genetic, informational parasitism itself, i.e., various degrees of reliance on the information processing systems of the host.

Occurrence and genetic diversity of the Plasmopara halstedii virus in sunflower downy mildew populations of the world

Plasmopara halstedii virus (PhV) is a ss(+)RNA virus that exclusively occurs in the sunflower downy mildew pathogen Plasmopara halstedii, a biotrophic oomycete of severe economic impact. The virus origin and its genomic variability are unknown. A PCR-based screening of 128 samples of P. halstedii from five continents and up to 40 y old was conducted. PhV RNA was found in over 90 % of the isolates with no correlation to geographic origin or pathotype of its host. Sequence analyses of the two open reading frames (ORFs) revealed only 18 single nucleotide polymorphisms (SNPs) in 3873 nucleotides. The SNPs had no recognizable effect on the two encoded virus proteins. In 398 nucleotides of the untranslated regions (UTRs) of the RNA 2 strand eight additional SNPs and one short deletion was found. Modelling experiments revealed no effects of these variations on the secondary structure of the RNA. The results showed the presence of PhV in P. halstedii isolates of global origin and the existence of the virus since more than 40 y. The virus genome revealed a surprisingly low variation in both coding and noncoding parts. No sequence differences were correlated with host pathotype or geographic populations of the oomycete.

Novel virus discovery and genome reconstruction from field RNA samples reveals highly divergent viruses in dipteran hosts

We investigated whether small RNA (sRNA) sequenced from field-collected mosquitoes and chironomids (Diptera) can be used as a proxy signature of viral prevalence within a range of species and viral groups, using sRNAs sequenced from wild-caught specimens, to inform total RNA deep sequencing of samples of particular interest. Using this strategy, we sequenced from adult Anopheles maculipennis s.l. mosquitoes the apparently nearly complete genome of one previously undescribed virus related to chronic bee paralysis virus, and, from a pool of Ochlerotatus caspius and Oc. detritus mosquitoes, a nearly complete entomobirnavirus genome. We also reconstructed long sequences (1503-6557 nt) related to at least nine other viruses. Crucially, several of the sequences detected were reconstructed from host organisms highly divergent from those in which related viruses have been previously isolated or discovered. It is clear that viral transmission and maintenance cycles in nature are likely to be significantly more complex and taxonomically diverse than previously expected.

Phytophthora viruses

Phytophthora sp. is a genus in the oomycetes, which are similar to filamentous fungi in morphology and habitat, but phylogenetically more closely related to brown algae and diatoms and fall in the kingdom Stramenopila. In the past few years, several viruses have been characterized in Phytophthora species, including four viruses from Phytophthora infestans, the late blight pathogen, and an endornavirus from an unnamed Phytophthora species from Douglas fir. Studies on Phytophthora viruses have revealed several interesting systems. Phytophthora infestans RNA virus 1 (PiRV-1) and PiRV-2 are likely the first members of two new virus families; studies on PiRV-3 support the establishment of a new virus genus that is not affiliated with established virus families; PiRV-4 is a member of Narnaviridae, most likely in the genus Narnavirus; and Phytophthora endornavirus 1 (PEV1) was the first nonplant endornavirus at the time of reporting. Viral capsids have not been found in any of the above-mentioned viruses. PiRV-1 demonstrated a unique genome organization that requires further examination, and PiRV-2 may have played a role in late blight resurgence in 1980s-1990s.

The family narnaviridae: simplest of RNA viruses

Members of the virus family Narnaviridae contain the simplest genomes of any RNA virus, ranging from 2.3 to 3.6 kb and encoding only a single polypeptide that has an RNA-dependent RNA polymerase domain. The family is subdivided into two genera based on subcellular location: members of the genus Narnavirus have been found in the yeast Saccharomyces cerevisiae and in the oomycete Phytophthora infestans and are confined to the cytosol, while members of the genus Mitovirus have been found only in filamentous fungi and are found in mitochondria. None identified thus far encodes a capsid protein; like several other RNA viruses of lower eukaryotes, their genomes are confined within lipid vesicles. As more family members are discovered, their importance as genetic elements is becoming evident. The unique association of the genus Mitovirus with mitochondria renders them potentially valuable tools to study biology of lower eukaryotes.