Identification of Diverse Mycoviruses through Metatranscriptomics Characterization of the Viromes of Five Major Fungal Plant Pathogens

New Insights on the Integrated Management of Plant Diseases by RNA Strategies: Mycoviruses and RNA Interference

RNA-based strategies for plant disease management offer an attractive alternative to agrochemicals that negatively impact human and ecosystem health and lead to pathogen resistance. There has been recent interest in using mycoviruses for fungal disease control after it was discovered that some cause hypovirulence in fungal pathogens, which refers to a decline in the ability of a pathogen to cause disease. Cryphonectria parasitica, the causal agent of chestnut blight, has set an ideal model of management through the release of hypovirulent strains. However, mycovirus-based management of plant diseases is still restricted by limited approaches to search for viruses causing hypovirulence and the lack of protocols allowing effective and systemic virus infection in pathogens. RNA interference (RNAi), the eukaryotic cell system that recognizes RNA sequences and specifically degrades them, represents a promising. RNA-based disease management method. The natural occurrence of cross-kingdom RNAi provides a basis for host-induced gene silencing, while the ability of most pathogens to uptake exogenous small RNAs enables the use of spray-induced gene silencing techniques. This review describes the mechanisms behind and the potential of two RNA-based strategies, mycoviruses and RNAi, for plant disease management. Successful applications are discussed, as well as the research gaps and limitations that remain to be addressed.

Metatranscriptomic Analysis Reveals Rich Mycoviral Diversity in Three Major Fungal Pathogens of Rice

In recent years, three major fungal diseases of rice, i.e., rice blast, rice false smut, and rice-sheath blight, have caused serious worldwide rice-yield reductions and are threatening global food security. Mycoviruses are ubiquitous in almost all major groups of filamentous fungi, oomycetes, and yeasts. To reveal the mycoviral diversity in three major fungal pathogens of rice, we performed a metatranscriptomic analysis of 343 strains, representing the three major fungal pathogens of rice, Pyricularia oryzae, Ustilaginoidea virens, and Rhizoctonia solani, sampled in southern China. The analysis identified 682 contigs representing the partial or complete genomes of 68 mycoviruses, with 42 described for the first time. These mycoviruses showed affinity with eight distinct lineages: Botourmiaviridae, Partitiviridae, Totiviridae, Chrysoviridae, Hypoviridae, Mitoviridae, Narnaviridae, and Polymycoviridae. More than half (36/68, 52.9%) of the viral sequences were predicted to be members of the families Narnaviridae and Botourmiaviridae. The members of the family Polymycoviridae were also identified for the first time in the three major fungal pathogens of rice. These findings are of great significance for understanding the diversity, origin, and evolution of, as well as the relationship between, genome structures and functions of mycoviruses in three major fungal pathogens of rice.

A novel fungal negative-stranded RNA virus related to mymonaviruses in Auricularia heimuer

Here, we report the characterization of a novel (-)ssRNA mycovirus isolated from Auricularia heimuer CCMJ1222, using a combination of RNA-seq, reverse transcription polymerase chain reaction, 5' and 3' rapid amplification of cDNA ends, and Sanger sequencing. Based on database searches, sequence alignment, and phylogenetic analysis, we designated the virus as "Auricularia heimuer negative-stranded RNA virus 1" (AhNsRV1). This virus has a monopartite RNA genome related to mymonaviruses (order Mononegavirales). The AhNsRV1 genome consists of 11,441 nucleotides and contains six open reading frames (ORFs). The largest ORF encodes a putative RNA-dependent RNA polymerase; the other ORFs encode hypothetical proteins with no conserved domains or known function. AhNsRV1 is the first (-)ssRNA virus and the third virus known to infect A. heimuer.

Fusarium Graminearum Virus-1 Strain FgV1-SD4 Infection Eliminates Mycotoxin Deoxynivalenol Synthesis by Fusarium graminearum in FHB

Deoxynivalenol (DON) toxin production during the infection of F. graminearum in small grain crops is one of the most harmful virulence factors associated with economic losses. Metatranscriptome sequencing and RT-qPCR traced back that the only mycovirus infecting an F. graminearum isolate, designated as Fg-4-2, was a novel strain of Fusarium graminearum virus 1 (FgV1), designated as FgV1-SD4. The isolate Fg-4-2 showed significantly reduced virulence against wheat compared to the virus-free culture, designated as isolate Fg-4-1, which was obtained by deep freezing and single conidial germination. Notably, no DON accumulation was detected in the harvested wheat seeds infected by Fg-4-2, whereas ~18 ppm DON was detected in seeds infected by Fg-4-1. Comparison of the genome sequence of FgV1-SD4 with other identified strains of FgV1, i.e., FgV1-DK21 and FgV1-ch, indicates mutations on ORF-2 and the 3′-UTR in the genome that might be associated with hypovirulence. This mycovirus strain alone and specific genetic components of FgV1-SD4 can be further optimized to be developed as a biocontrol agent to reduce Fusarium head blight and to lower the DON accumulation levels in small grain crops due to this fungal disease.

Mycovirus Hunting Revealed the Presence of Diverse Viruses in a Single Isolate of the Phytopathogenic Fungus Diplodia seriata From Pakistan

Diplodia seriata in the family Botryosphaeriaceae is a cosmopolitan phytopathogenic fungus and is responsible for causing cankers, fruit rot and leaf spots on economically important plants. In this study, we characterized the virome of a single Pakistani strain (L3) of D. seriata. Several viral-like contig sequences were obtained via a previously conducted next-generation sequencing analysis. Multiple infection of the L3 strain by eight RNA mycoviruses was confirmed through RT-PCR using total RNA samples extracted from this strain; the entire genomes were determined via Sanger sequencing of RT-PCR and RACE clones. A BLAST search and phylogenetic analyses indicated that these eight mycoviruses belong to seven different viral families. Four identified mycoviruses belong to double-stranded RNA viral families, including Polymycoviridae, Chrysoviridae, Totiviridae and Partitiviridae, and the remaining four identified mycoviruses belong to single-stranded RNA viral families, i.e., Botourmiaviridae, and two previously proposed families "Ambiguiviridae" and "Splipalmiviridae". Of the eight, five mycoviruses appear to represent new virus species. A morphological comparison of L3 and partially cured strain L3ht1 suggested that one or more of the three viruses belonging to Polymycoviridae, "Splipalmiviridae" and "Ambiguiviridae" are involved in the irregular colony phenotype of L3. To our knowledge, this is the first report of diverse virome characterization from D. seriata.

Resolving deep evolutionary relationships within the RNA virus phylum Lenarviricota

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

Diversity of mycoviruses in edible fungi

Mycoviruses (fungal viruses) are widespread in all major taxonomic groups of fungi. Although most mycovirus infections are latent, some mycoviruses, such as La France isometric virus, mushroom virus X, and oyster mushroom spherical virus, can cause severe diseases in edible fungi and lead to significant production losses. Recently, deep sequencing has been employed as a powerful research tool to identify new mycoviruses and to enhance our understanding of virus diversity and evolution. An increasing number of novel mycoviruses that can infect edible fungi have been reported, including double-stranded (ds) RNA, positive-sense ( +)ssRNA, and negative-sense (-)ssRNA viruses. To date, approximately 60 mycoviruses have been reported in edible fungi. In this review, we summarize the recent advances in the diversity and evolution of mycoviruses that can infect edible fungi. We also discuss mycovirus transmission, co-infections, and genetic variations, as well as the methods used to detect and control of mycoviruses in edible fungi, and provide insights for future research on mushroom viral diseases.

The Mycovirome in a Worldwide Collection of the Brown Rot Fungus Monilinia fructicola

The fungus Monilinia fructicola is responsible for brown rot on stone and pome fruit and causes heavy yield losses both pre- and post-harvest. Several mycoviruses are known to infect fungal plant pathogens. In this study, a metagenomic approach was applied to obtain a comprehensive characterization of the mycovirome in a worldwide collection of 58 M. fructicola strains. Deep sequencing of double-stranded (ds)RNA extracts revealed a great abundance and variety of mycoviruses. A total of 32 phylogenetically distinct positive-sense (+) single-stranded (ss)RNA viruses were identified. They included twelve mitoviruses, one in the proposed family Splipalmiviridae, and twelve botourmiaviruses (phylum Lenarviricota), eleven of which were novel viral species; two hypoviruses, three in the proposed family Fusariviridae, and one barnavirus (phylum Pisuviricota); as well as one novel beny-like virus (phylum Kitrinoviricota), the first one identified in Ascomycetes. A partial sequence of a new putative ssDNA mycovirus related to viruses within the Parvoviridae family was detected in a M. fructicola isolate from Serbia. The availability of genomic sequences of mycoviruses will serve as a solid basis for further research aimed at deepening the knowledge on virus–host and virus–virus interactions and to explore their potential as biocontrol agents against brown rot disease.

A Novel Heptasegmented Positive-Sense Single-Stranded RNA Virus from the Phytopathogenic Fungus Colletotrichum fructicola

In this study, a novel positive-sense single-stranded RNA (+ssRNA) mycovirus, tentatively named Colletotrichum fructicola RNA virus 1 (CfRV1), was identified in the phytopathogenic fungus Colletotrichum fructicola. CfRV1 has seven genomic components, encoding seven proteins from open reading frames (ORFs) flanked by highly conserved untranslated regions (UTRs). Proteins encoded by ORFs 1, 2, 3, 5, and 6 are more similar to the putative RNA-dependent RNA polymerase (RdRp), hypothetical protein (P2), methyltransferase, and two hypothetical proteins of Hadaka virus 1 (HadV1), a capsidless 10- or 11-segmented +ssRNA virus, while proteins encoded by ORFs 4 and 7 showed no detectable similarity to any known proteins. Notably, proteins encoded by ORFs 1 to 3 also share considerably high similarity with the corresponding proteins of polymycoviruses. Phylogenetic analysis conducted based on the amino acid sequence of CfRV1 RdRp and related viruses placed CfRV1 and HadV1 together in the same clade, close to polymycoviruses and astroviruses. CfRV1-infected C. fructicola strains demonstrate a moderately attenuated growth rate and virulence compared to uninfected isolates. CfRV1 is capsidless and potentially encapsulated in vesicles inside fungal cells, as revealed by transmission electron microscopy. CfRV1 and HadV1 are +ssRNA mycoviruses closely related to polymycoviruses and astroviruses, represent a new linkage between +ssRNA viruses and the intermediate double-stranded RNA (dsRNA) polymycoviruses, and expand our understanding of virus diversity, taxonomy, evolution, and biological traits. IMPORTANCE A scenario proposing that dsRNA viruses evolved from +ssRNA viruses is still considered controversial due to intergroup knowledge gaps in virus diversity. Recently, polymycoviruses and hadakaviruses were found as intermediate dsRNA and +ssRNA stages, respectively, between +ssRNA and dsRNA viruses. Here, we identified a novel +ssRNA mycovirus, Colletotrichum fructicola RNA virus 1 (CfRV1), isolated from Colletotrichum fructicola in China. CfRV1 is phylogenetically related to the 10- or 11-segmented Hadaka virus 1 (HadV1) but consists of only seven genomic segments encoding two novel proteins. CfRV1 is naked and may be encapsulated in vesicles inside fungal cells, representing a potential novel lifestyle for multisegmented RNA viruses. CfRV1 and HadV1 are intermediate +ssRNA mycoviruses in the linkage between +ssRNA viruses and the intermediate dsRNA polymycoviruses and expand our understanding of virus diversity, taxonomy, and evolution.

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.

Two Novel Rhabdoviruses Related to Hypervirulence in a Phytopathogenic Fungus

Rhabdoviruses are ubiquitous and diverse viruses that propagate owing to bidirectional interactions with their vertebrate, arthropod, and plant hosts, and some of them could pose global health or agricultural threats. However, rhabdoviruses have rarely been reported in fungi. Here, two newly identified fungal rhabdoviruses, Rhizoctonia solani rhabdovirus 1 (RsRhV1) and RsRhV2, were discovered and molecularly characterized from the phytopathogenic fungus Rhizoctonia solani. The genomic organizations of RsRhV1 and RsRhV2 are 11,716 and 11,496 nucleotides (nt) in length, respectively, and consist of five open reading frames (ORFs) (ORFs I to V). ORF I, ORF IV, and ORF V encode the viral nucleocapsid (N), glycoprotein (G), and RNA polymerase (L), respectively. The putative protein encoded by ORF III has a lower level of identity with the matrix protein of rhabdoviruses. ORF II encodes a hypothetical protein with unknown function. Phylogenetic trees based on multiple alignments of N, L, and G proteins revealed that RsRhV1 and RsRhV2 are new members of the family Rhabdoviridae, but they form an independent evolutionary branch significantly distinct from other known nonfungal rhabdoviruses, suggesting that they represent a novel viral evolutionary lineage within Rhabdoviridae. Compared to strains lacking rhabdoviruses, strains harboring RsRhV2 and RsRhV1 showed hypervirulence, suggesting that RsRhV1 and RsRhV2 might be associated with the virulence of R. solani. Taken together, this study enriches our understanding of the diversity and host range of rhabdoviruses. IMPORTANCE Mycoviruses have been attracting an increasing amount of attention due to their impact on important medical, agricultural, and industrial fungi. Rhabdoviruses are prevalent across a wide spectrum of hosts, from plants to invertebrates and vertebrates. This study molecularly characterized two novel rhabdoviruses from four Rhizoctonia solani strains, based on their genomic structures, transcription strategy, phylogenetic relationships, and biological impact on their host. Our study makes a significant contribution to the literature because it not only enriches the mycovirus database but also expands the known host range of rhabdoviruses. It also offers insight into the evolutionary linkage between animal viruses and mycoviruses and the transmission of viruses from one host to another. Our study will also help expand the contemporary knowledge of the classification of rhabdoviruses, as well as providing a new model to study rhabdovirus-host interactions, which will benefit the agriculture and medical areas of human welfare.

A novel hexa-segmented dsRNA mycovirus confers hypovirulence in the phytopathogenic fungus Diaporthe pseudophoenicicola

A novel hexa-segmented double-stranded RNA (dsRNA) mycovirus was isolated and characterized from the filamentous phytopathogenic fungus Diaporthe pseudophoenicicola, and was named Diaporthe pseudophoenicicola chrysovirus 1 (DpCV1). The full-length cDNAs of dsRNA1-6 were 3335, 3030, 3039, 2980, 963, and 780 bp, respectively. Sequence analysis indicated the presence of nine open reading frames (ORFs) in the DpCV1 genome. ORF1 in dsRNA1 putatively encoded the RNA-dependent RNA polymerase (RdRp), and ORF3 in dsRNA2 encoded a capsid protein (CP). The seven remaining ORFs, ORF2 in dsRNA2, ORF4 in dsRNA3, ORF6, 7 in dsRNA4, ORF8 in dsRNA5, and ORF9 in dsRNA6, encoded proteins with unknown functions. Phylogenetic analysis revealed that DpCV1 is closely related to members of the cluster I group within the family Chrysoviridae but formed a separate clade. Importantly, all the six segments of DpCV1 were cured successfully through single spore isolation to obtain the isogenic virus-free strains. DpCV1 can confer hypovirulence to the fungal host of Diaporthe pseudophoenicicola. Compared with the virus-free strain, WC02 harboring the DpCV1 is more sensitive to fungicide prochloraz. Furthermore, the cell wall of DpCV1 infected strain was loose and enlarged. This is the first report of a hexa-segmented tentative chrysovirus in D. pseudophoenicicola. This article is protected by copyright. All rights reserved.

Genome Characterization and Phylogenetic Analysis of a Novel Endornavirus That Infects Fungal Pathogen Sclerotinia sclerotiorum

Endornaviruses are capsidless linear (+) ssRNA viruses in the family Endornaviridae. In this study, Scelrotinia sclerotiorum endornavirus 11 (SsEV11), a novel endornavirus infecting hypovirulent Sclerotinia sclerotiorum strain XY79, was identified and cloned using virome sequencing analysis and rapid amplification of cDNA ends (RACE) techniques. The full-length genome of SsEV11 is 11906 nt in length with a large ORF, which encodes a large polyprotein of 3928 amino acid residues, containing a viral methyltransferase domain, a cysteine-rich region, a putative DEADc, a viral helicase domain, and an RNA-dependent RNA polymerase (RdRp) 2 domain. The 5’ and 3’ untranslated regions (UTR) are 31 nt and 90 nt, respectively. According to the BLAST result of the nucleotide sequence, SsEV11 shows the highest identity (45%) with Sclerotinia minor endornavirus 1 (SmEV1). Phylogenetic analysis based on amino acid sequence of RdRp demonstrated that SsEV11 clusters to endornavirus and has a close relationship with Betaendornavirus. Phylogenetic analysis based on the sequence of endornaviral RdRp domain indicated that there were three large clusters in the phylogenetic tree. Combining the results of alignment analysis, Cluster I at least has five subclusters including typical members of Alphaendornavirus and many unclassified endornaviruses that isolated from fungi, oomycetes, algae, and insects; Cluster II also has five subclusters including typical members of Betaendornavirus, SsEV11, and other unclassified viruses that infected fungi; Cluster III includes many endorna-like viruses that infect nematodes, mites, and insects. Viruses in Cluster I and Cluster II are close to each other and relatively distant to those in Cluster III. Our study characterized a novel betaendornavirus, SsEV11, infected fungal pathogen S. sclerotiorum, and suggested that notable phylogenetic diverse exists in endornaviruses. In addition, at least, one novel genus, Gammaendornavirus, should be established to accommodate those endorna-like viruses in Cluster III.

Complete nucleotide sequence of a novel alphapartitivirus from Rhizoctonia solani AG-4 HG III isolate SM03

In this report, the full genome sequence of a novel mycovirus, designated as "Rhizoctonia solani partitivirus SM03" (RsPV-SM03), was determined in Rhizoctonia solani AG-4 HG III isolate SM03. RsPV-SM03 genome consists of two dsRNAs (dsRNA-1 and dsRNA-2), each of them contains one single open reading frame (ORF). ORF1 of dsRNA-1 encodes a putative RNA-dependent RNA polymerase (RdRp), while ORF2 of dsRNA-2 encodes a putative viral coat protein (CP). Phylogenetic analysis indicated that the RdRp and CP of RsPV-SM03 are closely related to those of other members of the genus Alphapartitivirus, family Partitiviridae, suggesting that RsPV-SM03 represents a novel species in the genus Alphapartitivirus.

Molecular characterization of a novel penoulivirus from the phytopathogenic fungus Colletotrichum camelliae

Colletotrichum camelliae is a widespread filamentous phytopathogenic fungus. In this study, a novel mycovirus designated as "Colletotrichum camelliae botourmiavirus 1" (CcBV1) was isolated from strain ZJQT11 of C. camelliae, and its complete genome sequence was determined. CcBV1 has a genome of 2,506 nucleotides and contains a large open reading frame (ORF) that encodes an RNA-dependent RNA polymerase (RdRp) with 672 amino acids and a predicted molecular mass of 75.23 kDa. A BLASTp search showed that RdRp encoded by CcBV1 is closely related to that of Pyricularia oryzae ourmia-like virus 1 with 73.22% identity. Phylogenetic analysis indicated that CcBV1 clustered in the penoulivirus clade within the family Botourmiaviridae. To the best of our knowledge, this is the first report of a penoulivirus in C. camelliae.

Four Novel Mycoviruses from the Hypovirulent Botrytis cinerea SZ-2-3y Isolate from Paris polyphylla: Molecular Characterisation and Mitoviral Sequence Transboundary Entry into Plants

A hypovirulent SZ-2-3y strain isolated from diseased Paris polyphylla was identified as Botrytis cinerea. Interestingly, SZ-2-3y was coinfected with a mitovirus, two botouliviruses, and a 3074 nt fusarivirus, designated Botrytis cinerea fusarivirus 8 (BcFV8); it shares an 87.2% sequence identity with the previously identified Botrytis cinerea fusarivirus 6 (BcFV6). The full-length 2945 nt genome sequence of the mitovirus, termed Botrytis cinerea mitovirus 10 (BcMV10), shares a 54% sequence identity with Fusarium boothii mitovirus 1 (FbMV1), and clusters with fungus mitoviruses, plant mitoviruses and plant mitochondria; hence BcMV10 is a new Mitoviridae member. The full-length 2759 nt and 2812 nt genome sequences of the other two botouliviruses, named Botrytis cinerea botoulivirus 18 and 19 (BcBoV18 and 19), share a 40% amino acid sequence identity with RNA-dependent RNA polymerase protein (RdRp), and these are new members of the Botoulivirus genus of Botourmiaviridae. Horizontal transmission analysis showed that BcBoV18, BcBoV19 and BcFV8 are not related to hypovirulence, suggesting that BcMV10 may induce hypovirulence. Intriguingly, a partial BcMV10 sequence was detected in cucumber plants inoculated with SZ-2-3y mycelium or pXT1/BcMV10 agrobacterium. In conclusion, we identified a hypovirulent SZ-2-3y fungal strain from P. polyphylla, coinfected with four novel mycoviruses that could serve as potential biocontrol agents. Our findings provide evidence of cross-kingdom mycoviral sequence transmission.

Non-detection of mycoviruses in amphibian chytrid fungus (Batrachochytrium dendrobatidis) from Australia

Mycoviruses may influence the pathogenicity of disease-causing fungi. Although mycoviruses have been found in some chytrid fungi, limited testing has not detected them in Batrachochytrium dendrobatidis (Bd), the cause of the devastating amphibian disease, chytridiomycosis. Here we conducted a survey for mycovirus presence in 38 Bd isolates from Australia (n = 31), Brazil (n = 5) and South Korea (n = 2) with a combination of modern high-throughput sequencing and conventional dsRNA cellulose chromatography. Mycoviruses were not detected in any isolates. This result was unexpected, given the long evolutionary history of Bd, as well as the high prevalence of mycoviruses in related fungal species. Given our widespread sampling in Australia and the limited number of Bd introductions, we suggest that mycoviruses are uncommon or absent from Australian Bd. Testing more isolates from regions where Bd originated, as well as regions with high diversity or low fungal virulence may identify mycoviruses that could aid in disease control.

Metagenomics Analysis of the Wheat Virome Identifies Novel Plant and Fungal-Associated Viral Sequences

Wheat viruses including wheat streak mosaic virus, Triticum mosaic virus, and barley yellow dwarf virus cost substantial losses in crop yields every year. Although there have been extensive studies conducted on these known wheat viruses, currently, there is limited knowledge about all components of the wheat (Triticum aestivum L.) virome. Here, we determined the composition of the wheat virome through total RNA deep sequencing of field-collected leaf samples. Sequences were de novo assembled after removing the host reads, and BLASTx searches were conducted. In addition to the documented wheat viruses, novel plant and fungal-associated viral sequences were identified. We obtained the full genome sequence of the first umbra-like associated RNA virus tentatively named wheat umbra-like virus in cereals. Moreover, a novel bi-segmented putative virus tentatively named wheat-associated vipovirus sharing low but significant similarity with both plant and fungal-associated viruses was identified. Additionally, a new putative fungal-associated tobamo-like virus and novel putative Mitovirus were discovered in wheat samples. The discovery and characterization of novel viral sequences associated with wheat is important to determine if these putative viruses may pose a threat to the wheat industry or have the potential to be used as new biological control agents for wheat pathogens either as wild-type or recombinant viruses.

Complete genome sequence of a novel mitovirus from binucleate Rhizoctonia AG-K strain FAS2909W

The full-length AU-rich (63.14%) 2,794-nucleotide sequence of Rhizoctonia mitovirus K1 (RMV-K1) isolated from the binucleate AG-K strain FAS2909W was determined. The positive strand of RMV-K1 contains a large open reading frame (ORF) when the fungal mitochondrial genetic code is used. This ORF was predicted to encode an RdRp protein exhibiting the highest sequence identity (41.77%) to Rhizoctonia solani mitovirus 30. Phylogenetic analysis showed that RMV-K1 is a novel member of the genus Mitovirus, family Mitoviridae.

Metagenomic Identification of Viral Sequences in Laboratory Reagents

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

Cyanolichen microbiome contains novel viruses that encode genes to promote microbial metabolism

Lichen thalli are formed through the symbiotic association of a filamentous fungus and photosynthetic green alga and/or cyanobacterium. Recent studies have revealed lichens also host highly diverse communities of secondary fungal and bacterial symbionts, yet few studies have examined the viral component within these complex symbioses. Here, we describe viral biodiversity and functions in cyanolichens collected from across North America and Europe. As current machine-learning viral-detection tools are not trained on complex eukaryotic metagenomes, we first developed efficient methods to remove eukaryotic reads prior to viral detection and a custom pipeline to validate viral contigs predicted with three machine-learning methods. Our resulting high-quality viral data illustrate that every cyanolichen thallus contains diverse viruses that are distinct from viruses in other terrestrial ecosystems. In addition to cyanobacteria, predicted viral hosts include other lichen-associated bacterial lineages and algae, although a large fraction of viral contigs had no host prediction. Functional annotation of cyanolichen viral sequences predicts numerous viral-encoded auxiliary metabolic genes (AMGs) involved in amino acid, nucleotide, and carbohydrate metabolism, including AMGs for secondary metabolism (antibiotics and antimicrobials) and fatty acid biosynthesis. Overall, the diversity of cyanolichen AMGs suggests that viruses may alter microbial interactions within these complex symbiotic assemblages.

A New Double-Stranded RNA Mycovirus in Cryphonectria naterciae Is Able to Cross the Species Barrier and Is Deleterious to a New Host

Cryphonectria is a fungal genus associated with economically significant disease of trees. Herein we characterized a novel double-stranded RNA virus from the fungal species Cryphonectria naterciae, a species unexplored as a virus host. De novo assembly of RNA-seq data and Sanger sequencing of RACE (rapid amplification of cDNA ends) clones gave the complete, non-segmented genome (10,164 bp) of the virus termed Cryphonectria naterciae fusagravirus (CnFGV1) that was phylogenetically placed within the previously proposed viral family Fusagraviridae. Of 31 field-collected strains of C. naterciae, 40% tested CnFGV1-positive. Cocultivation resulted in within-species transmission of CnFGV1 to virus-free strains of C. naterciae. Comparison of the mycelium phenotype and the growth rate of CnFGV1-infected and virus-free isogenic strains revealed frequent sectoring and growth reduction in C. naterciae upon virus infection. Co-culturing also led to cross-species transmission of CnFGV1 to Cryphonectria carpinicola and Cryphonectria radicalis, but not to Cryphonectria parasitica. The virus-infected C. naterciae and the experimentally infected Cryphonectria spp. readily transmitted CnFGV1 through asexual spores to the next generation. CnFGV1 strongly reduced conidiation and in some cases vegetative growth of C. carpinicola, which is involved in the European hornbeam disease. This study is the first report of a fusagravirus in the family Cryphonectriaceae and lays the groundwork for assessing a hypovirulence effect of CnFGV1 against the hornbeam decline in Europe.

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.

Mycoviruses infecting Colletotrichum spp.: A comprehensive review

Colletotrichum is one of the most economically important fungal genera, which affects a wide range of hosts, specifically tropical and subtropical crops. Thus far, there have been several records of mycovirus infection in Colletotrichum spp., primarily by viruses of the Partitiviridae family. There have also been records of infections by mycoviruses of the Chrysoviridae family. Mycoviruses are (+)ssRNA and dsRNA genome viruses, which may or may not be enveloped. To date, no mycovirus with a DNA genome has been isolated from Colletotrichum spp. Typically, mycoviruses cause latent infections, although hypo- and hypervirulence have also been reported in Colletotrichum spp. In addition to its effects on pathogenic behavior, mycovirus infection can lead to important physiological changes, such as altered morphological characteristics, reduced vegetative growth, and suppressed conidia production. Therefore, research on mycoviruses infecting phytopathogenic fungi can help develop alternative methods to chemical control, which can cause irreversible damage to humans and the environment. From an agricultural perspective, mycoviruses can contribute to sustainable agriculture as biological control agents via changes in fungal physiology, ultimately resulting in the total loss of or reduction in the virulence of these pathogens.

Direct Metatranscriptomic Survey of the Sunflower Microbiome and Virome

Sunflowers (Helianthus annuus L.) are susceptible to multiple diseases in field production. In this study, we collected diseased sunflower leaves in fields located in South Dakota, USA, for virome investigation. The leaves showed visible symptoms on the foliage, indicating phomopsis and rust infections. To identify the viruses potentially associated with the disease diagnosed, symptomatic leaves were obtained from diseased plants. Total RNA was extracted corresponding to each disease diagnosed to generate libraries for paired-end high throughput sequencing. Short sequencing reads were assembled de novo and the contigs with similarities to viruses were identified by aligning against a custom protein database. We report the discovery of two novel mitoviruses, four novel partitiviruses, one novel victorivirus, and nine novel totiviruses based on similarities to RNA-dependent RNA polymerases and capsid proteins. Contigs similar to bean yellow mosaic virus and Sclerotinia sclerotiorum hypovirulence-associated DNA virus were also detected. To the best of our knowledge, this is the first report of direct metatranscriptomics discovery of viruses associated with fungal infections of sunflowers bypassing culturing. These newly discovered viruses represent a natural genetic resource from which we can further develop potential biopesticide to control sunflower diseases.

Nine viruses from eight lineages exhibiting new evolutionary modes that co-infect a hypovirulent phytopathogenic fungus

Mycoviruses are an important component of the virosphere, but our current knowledge of their genome organization diversity and evolution remains rudimentary. In this study, the mycovirus composition in a hypovirulent strain of Sclerotinia sclerotiorum was molecularly characterized. Nine mycoviruses were identified and assigned into eight potential families. Of them, six were close relatives of known mycoviruses, while the other three had unique genome organizations and evolutionary positions. A deltaflexivirus with a tripartite genome has evolved via arrangement and horizontal gene transfer events, which could be an evolutionary connection from unsegmented to segmented RNA viruses. Two mycoviruses had acquired a second helicase gene by two different evolutionary mechanisms. A rhabdovirus representing an independent viral evolutionary branch was the first to be confirmed to occur naturally in fungi. The major hypovirulence-associated factor, an endornavirus, was finally corroborated. Our study expands the diversity of mycoviruses and potential virocontrol agents, and also provides new insights into virus evolutionary modes including virus genome segmentation.

Identification of mycoviruses in phytopathogenic fungi is necessary for understanding the origin of viruses and developing virocontrol strategies to protect plants. Nine mycoviruses with RNA genomes were identified in a hypovirulent strain of Sclerotinia sclerotiorum and were classified into eight potential viral families, suggesting that the composition of mycoviral communities was complex in this single fungal strain. They included four previously characterized mycoviruses and three distant relatives of known mycoviruses, as well as the first reports of a deltaflexivirus with a tripartite genome, and a fungal rhabdovirus. In addition, we found an endornavirus associated with hypovirulence in a phytopathogenic fungus. Our study makes a significant contribution because it not only expands the diversity-related knowledge of mycoviruses and potential virocontrol agents, but also provides new insights into mycovirus evolution.

Towards the Forest Virome: High-Throughput Sequencing Drastically Expands Our Understanding on Virosphere in Temperate Forest Ecosystems

Thanks to the development of HTS technologies, a vast amount of genetic information on the virosphere of temperate forests has been gained in the last seven years. To estimate the qualitative/quantitative impact of HTS on forest virology, we have summarized viruses affecting major tree/shrub species and their fungal associates, including fungal plant pathogens, mutualists and saprotrophs. The contribution of HTS methods is extremely significant for forest virology. Reviewed data on viral presence in holobionts allowed us a first attempt to address the role of virome in holobionts. Forest health is dependent on the variability of microorganisms interacting with the host tree/holobiont; symbiotic microbiota and pathogens engage in a permanent interplay, which influences the host. Through virus–virus interplays synergistic or antagonistic relations may evolve, which may drastically affect the health of the holobiont. Novel insights of these interplays may allow practical applications for forest plant protection based on endophytes and mycovirus biocontrol agents. The current analysis is conceived in light of the prospect that novel viruses may initiate an emergent infectious disease and that measures for the avoidance of future outbreaks in forests should be considered.

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.

The Sisal Virome: Uncovering the Viral Diversity of Agave Varieties Reveals New and Organ-Specific Viruses

Sisal is a common name for different plant varieties in the genus Agave (especially Agave sisalana) used for high-quality natural leaf fiber extraction. Despite the economic value of these plants, we still lack information about the diversity of viruses (virome) in non-tequilana species from the genus Agave. In this work, by associating RNA and DNA deep sequencing we were able to identify 25 putative viral species infecting A. sisalana, A. fourcroydes, and Agave hybrid 11648, including one strain of Cowpea Mild Mottle Virus (CPMMV) and 24 elements likely representing new viruses. Phylogenetic analysis indicated they belong to at least six viral families: Alphaflexiviridae, Betaflexiviridae, Botourmiaviridae, Closteroviridae, Partitiviridae, Virgaviridae, and three distinct unclassified groups. We observed higher viral taxa richness in roots when compared to leaves and stems. Furthermore, leaves and stems are very similar diversity-wise, with a lower number of taxa and dominance of a single viral species. Finally, approximately 50% of the identified viruses were found in all Agave organs investigated, which suggests that they likely produce a systemic infection. This is the first metatranscriptomics study focused on viral identification in species from the genus Agave. Despite having analyzed symptomless individuals, we identified several viruses supposedly infecting Agave species, including organ-specific and systemic species. Surprisingly, some of these putative viruses are probably infecting microorganisms composing the plant microbiota. Altogether, our results reinforce the importance of unbiased strategies for the identification and monitoring of viruses in plant species, including those with asymptomatic phenotypes.

The full-length genome sequence of a novel mitovirus from Botryosphaeria dothidea, the causal agent of pear ring rot disease

Here, we describe a novel mycovirus, tentatively designated as "Botryosphaeria dothidea mitovirus 3" (BdMV3), isolated from Botryosphaeria dothidea strain FJ, which causes pear ring rot disease in Fujian Province, China. The complete genome nucleotide sequence of BdMV3 is 2538 nt in length and contains a single 2070-nt open reading frame (ORF) encoding a putative RNA-dependent RNA polymerase (RdRp) of 689 amino acids (aa) using the fungal mitochondrial genetic code. BLASTp analysis revealed that the RdRp of BdMV3 shares 28.91%-69.36% sequence identity (query sequence coverage more than 90%) with those of members of the genus Mitovirus, with the highest sequence identity of 69.36% and 68.79% to the corresponding RdRp aa sequences of Rhizoctonia solani mitovirus 10 and Macrophomina phaseolina mitovirus 4, respectively. Phylogenetic analysis based on RdRp aa sequences indicated that BdMV3 is a new member of the genus Mitovirus in the family Mitoviridae.

Characterization of a novel botoulivirus isolated from the phytopathogenic fungus Sclerotinia sclerotiorum

Sclerotinia sclerotiorum ourmiavirus 17 (SsOV17) was isolated from the hypovirulent strain GF3 of Sclerotinia sclerotiorum. The genome of SsOV17 is 2,802 nt in length and contains a single long open reading frame (ORF) flanked by a short structured 5'-untranslated region (5'-UTR) (28 nt) and a long 3'-UTR (788 nt), respectively. The ORF encodes a protein with 663 amino acids and a predicted molecular mass of 75.0 kDa. A BLASTp search indicated that the protein encoded by SsOV17 is closely related to the putative RNA-dependent RNA polymerase (RdRp) of Sclerotinia sclerotiorum ourmiavirus 13 (71% identity). A multiple sequence alignment indicated that eight conserved amino acid motifs were present in the RdRp conserved region of SsOV17. Phylogenetic analysis demonstrated that SsOV17 clustered with members of the genus Botoulivirus.

A Novel Partitivirus That Confer Hypovirulence to the Plant Pathogenic Fungus Colletotrichum liriopes

Here, we report a novel double-stranded RNA virus designated Colletotrichum liriopes partitivirus 1 (ClPV1) from the plant pathogenic fungus C. liriopes. ClPV1 genome has two double stranded RNAs (dsRNAs), named as dsRNA 1 and dsRNA 2, which in the lengths of 1,807 and 1,706 bp, respectively. The dsRNA 1 and dsRNA 2 encoded proteins showing significant amino acid (aa) sequence identity to the RNA-dependent RNA polymerase (RdRp) and coat protein (CP) of partitiviruses, respectively. Phylogenetic analysis using the aa sequences of RdRp and CP indicated that ClPV1 was grouped to members of the putative Epsilonpartitivirus genus in the Partitiviridae family. Spherical viral particles in approximately 35 nm in diameter and packaging the ClPV1 genome were isolated. Virus elimination and virus transfection with purified viral particles, and biological comparison revealed that ClPV1 could reduce the virulence and conidia production of C. liriopes. To the best of our knowledge, this is the first report of mycovirus in C. liriopes fungus.

The unique genome organization of two novel fusariviruses hosted by the true morel mushroom Morchella esculenta

Two putative mycoviruses belonging to the proposed family "Fusariviridae" were identified in Morchella esculenta by sequencing of double-stranded RNAs extracted from the morel mushroom. These viruses were tentatively named "Morchella esculenta fusarivirus 1″ (MeFV1) and "Morchella esculenta fusarivirus 2″ (MeFV2). Including the poly(A) tail the complete genomes of MeFV1 and MeFV2 are composed of 9096 and 9011 nucleotides (nt) respectively. Both genomes contain four non-overlapping open reading frames (ORFs) in which the largest and the smallest ORFs are ORF2 and ORF3 for both genomes respectively. The ORF1 of MeFV1 and MeFV2 are preceded by the 5' untranslated regions (UTRs) of 27 and 37 nt respectively and encode 341 and 339 aa long proteins that do not exhibit significant similarity to any of the protein sequences present in GenBank database. The 1502 and 1511 aa long proteins encoded by ORF2 of MeFV1 and MeFV2 share 84.42% sequence identity to each other and are 58.54% and 58.57% identical to the RNA-dependent RNA polymerase (RdRp) of Morchella importuna fusarivirus 1 (MiFV1) respectively. Interestingly, a Promethin/LDAF1 protein domain that is associated with the endoplasmic reticulum (ER) and lipid droplet (LD) membranes was identified at the N terminal regions of MeFV1 and MeFV2 RdRps, implying that the replication of these viruses is linked to the lipid membranes. The ORF3 and ORF4 of MeFV1 and MeFV2 encode proteins (268 and 333 aa long, and 645 and 647 aa long respectively) that only share significant sequence similarities with the proteins encoded by the ORF2 and ORF3 of MiFV1 respectively. The 3' UTRs of MeFV1 and MeFV2 are 162 and 159 nt long respectively and both of them have 51 nt long terminal poly(A) traits. To our knowledge, MeFV1 and MeFV2 are the first fusariviruses identified in M. esculenta and this is the first study reporting on the presence of Promethin/LDAF1 domain in viral RdRps.

Molecular characterization of a novel partitivirus and a fusarivirus coinfecting the fungus Nigrospora sphaerica

In this report, we describe the molecular characterization of two novel mycoviruses coinfecting the plant pathogenic fungus Nigrospora sphaerica, which were designated "Nigrospora sphaerica fusarivirus 1" (NsFV1) and "Nigrospora sphaerica partitivirus 1" (NsPV1). NsFV1 has an undivided genome measuring 6,147 nt, excluding the polyA tail, and was predicted to contain two nonoverlapping open reading frames (ORF1 and 2). The larger ORF1 encodes a polyprotein containing a conserved RNA-dependent RNA polymerase (RdRp) and a helicase domain that has functions related to RNA replication, and the smaller ORF2 encodes a putative protein with an unknown function. NsPV1 consists of two genome segments, which measure 1,796 bp and 1,455 bp in length. Each of the two dsRNAs has a single ORF, and they are predicted to encode proteins with homology to viral RdRps and coat proteins of members of the family Partitiviridae. Phylogenetic analysis indicated that NsFV1 is a member of the recently proposed family "Fusariviridae", while NsPV1 was determined to belong to the genus Gammapartitivirus in the family Partitiviridae. To the best of our knowledge, this report is the first to describe mycoviruses infecting N. sphaerica.

Novel Mycoviruses Discovered in the Mycovirome of a Necrotrophic Fungus

Botrytis cinerea is one of the most important plant-pathogenic fungus. Products based on microorganisms can be used in biocontrol strategies alternative to chemical control, and mycoviruses have been explored as putative biological agents in such approaches. Here, we have explored the mycovirome of B. cinerea isolates from grapevine of Italy and Spain to increase the knowledge about mycoviral diversity and evolution, and to search for new widely distributed mycoviruses that could be active ingredients in biological products to control this hazardous fungus. A total of 248 B. cinerea field isolates were used for our metatranscriptomic study. Ninety-two mycoviruses were identified: 62 new mycoviral species constituting putative novel viral genera and families. Of these mycoviruses, 57 had a positive-sense single-stranded RNA (ssRNA) genome, 19 contained a double-stranded RNA (dsRNA) genome, 15 had a negative-sense ssRNA genome, and 1 contained a single-stranded DNA (ssDNA) genome. In general, ssRNA mycoviruses were widely distributed in all sampled regions, the ssDNA mycovirus was more frequently found in Spain, and dsRNA mycoviruses were scattered in some pools of both countries. Some of the identified mycoviruses belong to clades that have never been found associated with Botrytis species: Botrytis-infecting narnaviruses; alpha-like, umbra-like, and tymo-like ssRNA+ mycoviruses; trisegmented ssRNA- mycovirus; bisegmented and tetrasegmented dsRNA mycoviruses; and finally, an ssDNA mycovirus. Among the results obtained in this massive mycovirus screening, the discovery of novel bisegmented viruses, phylogenetically related to narnaviruses, is remarkable.IMPORTANCE The results obtained here have expanded our knowledge of mycoviral diversity, horizontal transfers, and putative cross-kingdom events. To date, this study presents the most extensive and wide diversity collection of mycoviruses infecting the necrotrophic fungus B. cinerea The collection included all types of mycoviruses, with dsRNA, ssRNA+, ssRNA-, and ssDNA genomes, most of which were discovered here, and some of which were previously reported as infecting B. cinerea or other plant-pathogenic fungi. Some of these mycoviruses are reported for the first time here associated with B. cinerea, as a trisegmented ssRNA- mycovirus and as an ssDNA mycovirus, but even more remarkablly, we also describe here four novel bisegmented viruses (binarnaviruses) not previously described in nature. The present findings significantly contribute to general knowledge in virology and more particularly in the field of mycovirology.

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

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

Armillaria root rot fungi host single-stranded RNA viruses

Species of Armillaria are distributed globally and include some of the most important pathogens of forest and ornamental trees. Some of them form large long-living clones that are considered as one of the largest organisms on earth and are capable of long-range spore-mediated transfer as well as vegetative spread by drought-resistant hyphal cords called rhizomorphs. However, the virus community infecting these species has remained unknown. In this study we used dsRNA screening and high-throughput sequencing to search for possible virus infections in a collection of Armillaria isolates representing three different species: Armillaria mellea from South Africa, A. borealis from Finland and Russia (Siberia) and A. cepistipes from Finland. Our analysis revealed the presence of both negative-sense RNA viruses and positive-sense RNA viruses, while no dsRNA viruses were detected. The viruses included putative new members of virus families Mymonaviridae, Botourmiaviridae and Virgaviridae and members of a recently discovered virus group tentatively named "ambiviruses" with ambisense bicistronic genomic organization. We demonstrated that Armillaria isolates can be cured of viruses by thermal treatment, which enables the examination of virus effects on host growth and phenotype using isogenic virus-infected and virus-free strains.

Interannual dynamics, diversity and evolution of the virome in Sclerotinia sclerotiorum from a single crop field

Mycovirus diversity is generally analyzed from isolates of fungal culture isolates at a single point in time as a snapshot. The stability of mycovirus composition within the same geographical location over time remains unclear. Not knowing how the population fluctuates in the field can be a source of unpredictability in the successful application of virocontrol. To better understand the changes over time, we monitored the interannual dynamics and abundance of mycoviruses infecting Sclerotinia sclerotiorum at a rapeseed-growing field for three years. We found that the virome in S. sclerotiorum harbors unique mycovirus compositions each year. In total, sixty-eight mycoviruses were identified, among which twenty-four were detected in all three successive years. These twenty-four mycoviruses can be classified as the members of the core virome in this S. sclerotiorum population, which show persistence and relatively high transmissibility under field conditions. Nearly two-thirds of the mycoviruses have positive-sense, single-stranded RNA genomes and were found consistently across all three years. Moreover, twenty-eight mycoviruses are newly described, including four novel, multi-segmented narnaviruses, and four unique bunyaviruses. Overall, the newly discovered mycoviruses in this study belong to as many as twenty families, into which eight were first identified in S. sclerotiorum, demonstrating evolutionarily diverse viromes. Our findings not only shed light on the annual variation of mycovirus diversity but also provide important virus evolutionary clues.

Genome structure and diversity of novel endornaviruses from wheat sharp eyespot pathogen Rhizoctonia cerealis

Rhizoctonia cerealis (teleomorph Ceratobasidium cereale) is a soil-borne plant pathogenic fungus that can cause sharp eyespot in wheat or yellow patch in grasses. In this study, 21 new endornavirus genomes were obtained from five R. cerealis strains through the high-throughput sequencing of viral double-stranded RNA. Eighteen viruses were identified as Alphaendornavirus, and three viruses were identified as new species of Betaendornavirus on the basis of the phylogenetic analysis of the deduced amino acid sequences of RNA-dependent RNA polymerase. Notably, 12 of the new alphaendornaviruses could encode two open reading frames (ORFs), which were a rare feature of Endornaviridae. The amino acid sequences encoded by ORF2 from different endornaviruses had very low identity, and their functions and evolution origins remained unclear. Different endornavirus species with remarkably different genome structures could be found in the same R. cerealis strain. This study indicated that endornaviruses are common in R. cerealis and display wide diversity. Betaendornaviruses were found in R. cerealis, and a new species was proposed. This study is the first to report that the endornaviruses from R. cerealis can encode two ORFs and enhances our understanding of the viruses in the Endornaviridae family.

Multiple mycoviruses identified in Pestalotiopsis spp. from Chinese bayberry

Background

Chinese bayberry (Myrica rubra) is a subtropical fruit crop widely grown in southern China. Twig dieback is a disease of Chinese bayberry caused by Pestalotiopsis spp. and results in great economic losses to Chinese bayberry production. A virus survey was conducted in the population of Pestalotiopsis spp. infecting M. rubra in China. We explored the viral diversity in Pestalotiopsis spp., which may provide resources for further development as biocontrol agents of twig dieback.

Methods

Strains of Pestalotiopsis spp. were isolated from diseased twigs of M. rubra, and cultured on potato dextrose agar for RNA extraction. The total RNA of each strain was extracted, mixed, and used for RNA sequencing. The resulting sequences were deduplicated, annotated, and then used for phylogenetic analysis.

Results

Seven novel viruses were characterized from 59 isolates of M. rubra collected from 14 localities in China. Based on the phylogenetic analysis, these viruses were classified into five viral families/orders, Botourmiaviridae, Mitoviridae, Partitiviridae, Tymovirales and Bunyavirales, and one virus, Pestalotiopsis negative-stranded RNA virus 1, which likely belongs to a new viral family.

Conclusions

Metatranscriptomics analysis showed the presence of various mycoviruses in Pestalotiopsis spp. isolated from M. rubra in China. The genomes of eight putative viruses were identified, seven of which were nearly full-length. Some of these viruses of Pestalotiopsis spp. may have the potential for the biological control of twig dieback of M. rubra.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12985-021-01513-3.

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 gammapartitivirus from the entomopathogenic fungus Metarhizium brunneum

Here, we report a novel partitivirus infecting Metarhizium brunneum, which was designated "Metarhizium brunneum partitivirus 2" (MbPV2). The complete genome of MbPV2 consists of two segments, dsRNA1 and dsRNA2, with each dsRNA possessing a single open reading frame (ORF). dsRNA1 (1,775 bp) encodes a conserved RNA-dependent RNA polymerase (RdRp) with the highest sequence similarity to Plasmopara viticola associated partitivirus 1 (PvAPV1), while dsRNA2 (1,568 bp) encodes a coat protein (CP) with the highest sequence similarity to Colletotrichum partitivirus 1 (CtParV1). Phylogenetic analysis based on RdRp sequences showed that MbPV2 is a new member of the genus Gammapartitivirus, family Partitiviridae. This is the first report of a gammapartitivirus that infects the entomopathogenic fungus Metarhizium brunneum.

A novel previously undescribed fusarivirus from the phytopathogenic fungus Setosphaeria turcica

A putative mycovirus belonging to the proposed family "Fusariviridae" was discovered in Setosphaeria turcica by sequencing a double-stranded RNA extracted from this phytopathogenic fungus. The virus was tentatively named "Setosphaeria turcica fusarivirus 1" (StFV1). StFV1 has a genome comprising 6685 nucleotides. The genome contains three open reading frames (ORF). The largest ORF, ORF1, is preceded by an untranslated region (UTR) of 16 nucleotides and separated from ORF2 by an intergenic region of 63 nucleotides. The smallest ORF, ORF3, overlaps ORF2 by 16 nucleotides and is followed by a 3'-UTR of 82 nucleotides. The protein encoded by ORF1 is 71.8%, 67.4% and 68.1% identical to the RNA-dependent RNA polymerases (RdRps) of Pleospora typhicola fusarivirus 1 (PtFV1), Plasmopara viticola lesion-associated fusarivirus 1 (PvlaFV1), and Plasmopara viticola lesion-associated fusarivirus 3 (PvlaFV3), respectively, but has less than 47% amino acid sequence identity to the RdRps of other fusariviruses. To our knowledge, this is the first fusarivirus discovered in S. turcica and the first virus to be identified in this fungus using conventional cloning methods.

Divergent RNA viruses in Macrophomina phaseolina exhibit potential as virocontrol agents

Macrophomina phaseolina is an important necrotrophic phytopathogenic fungus and cause extensive damage in many oilseed crops. Twelve M.phaseolina isolates with diverse biological phenotypes were selected for a high-throughput sequencing-based metatranscriptomic and bioinformatics analysis to identify viruses infecting M.phaseolina. The analysis identified 40 partial or nearly complete viral genome segments, 31 of which were novel viruses. Among these viral sequences, 43% of the viral genomes were double-stranded RNA (dsRNA), 47% were positive single-stranded RNA (ssRNA+), and the remaining 10% were negative sense-stranded RNA (ssRNA−). The 40 viruses showed affinity to 13 distinct viral lineages, including Bunyavirales (four viruses), Totiviridae (three viruses), Chrysoviridae (five viruses), Partitiviridae (four viruses), Hypoviridae (one virus), Endornaviridae (two viruses), Tombusviridae (three viruses), Narnaviridae (one virus), Potyviridae (one virus), Bromoviridae (one virus), Virgaviridae (six viruses), ‘Fusagraviridae’ (five viruses), and Ourmiavirus (four viruses). Two viruses are closely related to two families, Potyviridae and Bromoviridae, which previously contained no mycovirus species. Moreover, nine novel viruses associated with M.phaseolina were identified in the family Totiviridae, Endornaviridae, and Partitiviridae. Coinfection with multiple viruses is prevalent in M.phaseolina, with each isolate harboring different numbers of viruses, ranging from three to eighteen. Furthermore, the effects of the viruses on the fungal host were analyzed according to the biological characteristics of each isolate. The results suggested that M.phaseolina hypovirus 2, M.phaseolina fusagravirus virus 1-5 (MpFV1-5), M.phaseolina endornavirus 1-2 (MpEV1-2), M.phaseolina ourmia-like virus 1-3 (MpOLV1-3), M.phaseolina mitovirus 4 (MpMV4), and M.phaseolina mycobunyavirus 1-4 (MpMBV1-4) were only detected in hypovirulent isolates. Those viruses associated with hypovirulence might be used as biological control agents as an environmentally friendly alternative to chemical fungicides. These findings considerably expand our understanding of mycoviruses in M.phaseolina and unvailed the presence of a huge difference among viruses in isolates from different hosts in distant geographical regions. Together, the present study provides new knowledge about viral evolution and fungus-virus coevolution.

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

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

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.

Mycoviral diversity and characteristics of a negative-stranded RNA virus LeNSRV1 in the edible mushroom Lentinula edodes

Bioinformatics and RT-PCR analysis of RNA from four Lentinula edodes samples identified 22 different virus-like contigs comprising 15 novel and 3 previously reported viruses. We further investigated the Lentinula edodes negative-stranded RNA virus 1 (LeNSRV1) isolated from a symptomatic sample, whose virion is a filamentous particle with a diameter of ~15 nm and a length of ~1200 nm. RT-PCR analysis detected LeNSRV1 in 10 of the 56 Chinese L. edodes core collection strains and 6 of the 22 monokaryotic strains from the L. edodes strain HNZMD. Genetic variation analysis showed that the sequences encoding the nucleocapsid protein (ORF2) from all the aforementioned LeNSRV1 positive strains are very conservative. The results presented here may enrich our understanding of L. edodes virus diversity and the characteristics of LeNSRV1, and will promote further research on virus-host interaction in L. edodes.

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.

Virome characterization of Cryphonectria parasitica isolates from Azerbaijan unveiled a new mymonavirus and a putative new RNA virus unrelated to described viral sequences

Cryphonectria parasitica, the causal agent of chestnut blight is controlled in Europe through natural spread of Cryphonectria hypovirus 1 (CHV1), a mycovirus able to induce hypovirulence to the host. In recent years C. parasitica was reported infecting Azerbaijani population of chestnut, but the presence of CHV1 still needs to be confirmed. Aim of this work was to investigate fifty-five C. parasitica isolates collected in Azerbaijan to describe the associated viruses. Our work found i) the first negative-sense ssRNA virus known to infect C. parasitica naturally for which we propose the name Cryphonectria parasitica sclerotimonavirus 1 (CpSV1) and ii) an RNA sequence showing peculiar features suggesting a viral nature for which we propose the name Cryphonectria parasitica ambivirus 1 (CpaV1). The discovery of CpaV1 expands our knowledge of the RNA virosphere suggesting the existence of a new lineage that cannot presently be reliably associated to the monophyletic Riboviria.

Establishment of Neurospora crassa as a model organism for fungal virology

The filamentous fungus Neurospora crassa is used as a model organism for genetics, developmental biology and molecular biology. Remarkably, it is not known to host or to be susceptible to infection with any viruses. Here, we identify diverse RNA viruses in N. crassa and other Neurospora species, and show that N. crassa supports the replication of these viruses as well as some viruses from other fungi. Several encapsidated double-stranded RNA viruses and capsid-less positive-sense single-stranded RNA viruses can be experimentally introduced into N. crassa protoplasts or spheroplasts. This allowed us to examine viral replication and RNAi-mediated antiviral responses in this organism. We show that viral infection upregulates the transcription of RNAi components, and that Dicer proteins (DCL-1, DCL-2) and an Argonaute (QDE-2) participate in suppression of viral replication. Our study thus establishes N. crassa as a model system for the study of host-virus interactions.

The fungus Neurospora crassa is a model organism for the study of various biological processes, but it is not known to be infected by any viruses. Here, Honda et al. identify RNA viruses that infect N. crassa and examine viral replication and RNAi-mediated antiviral responses, thus establishing this fungus as a model for the study of host-virus interactions.