Deep sequencing of mycovirus-derived small RNAs from Botrytis species

Construction of a Mycoviral Infectious Clone for Reverse Genetics in Botrytis cinerea

The most important advances in our understanding of the viral life cycle, such as genome replication, packaging, transmission, and host interactions, have been made via the development of viral infectious full-length clones. Here, we describe the detailed protocols for the construction of an infectious clone derived from Botrytis virus F (BVF), a mycoflexivirus infecting the plant pathogenic fungus Botrytis cinerea, the determination of the complete sequence of the cloned mycovirus, the preparation of fungal protoplasts, and the transfection of protoplasts using transcripts derived from the BVF infectious clone.

Determination of the Mycovirome of a Necrotrophic Fungus

Next-generation sequencing (NGS) of total RNA has allowed the detection of novel viruses infecting different hosts, such as fungi, increasing our knowledge on virus horizontal transfer events among different hosts, virus diversity, and virus evolution. Here, we describe the detailed protocols for the isolation of the plant pathogenic fungus Botrytis cinerea, from grapevine plants showing symptoms of the mold gray disease, the culture and maintenance of the isolated B. cinerea strains, the extraction of total RNA from B. cinerea strains for NGS, the bioinformatics pipeline designed and followed to detect mycoviruses in the sequenced samples, and the validation of the in silico detected mycoviruses by different approaches.

RNA interference of Aspergillus flavus in response to Aspergillus flavus partitivirus 1 infection

RNA interference (RNAi) is one of the important defense responses against viral infection, but its mechanism and impact remain unclear in mycovirus infections. In our study, reverse genetics and virus-derived small RNA sequencing were used to show the antiviral responses of RNAi components in Aspergillus flavus infected with Aspergillus flavus partitivirus 1 (AfPV1). qRT-PCR revealed that AfPV1 infection induced the expression of the RNAi components in A. flavus compared with noninfected A. flavus. Knock mutants of each RNAi component were generated, but the mutants did not exhibit any obvious phenotypic changes compared with the A. flavus parental strain. However, after AfPV1 inoculation, production of AfPV1 was significantly less than in the parental strain. Furthermore, sporulation was greater in each AfPV1-infected mutant compared with the AfPV1-infected parental A. flavus. We also investigated the sensitivity of virus-free and AfPV1-infected RNAi mutants and the parental strain to cell wall stress, osmotic stress, genotoxic stress, and oxidative stress. The mutants of DCLs and AGOs infected by AfPV1 displayed more changes than RDRP mutants in response to the first three stresses. Small RNA sequencing analysis suggested that AfPV1 infection reduced the number of unique reads of sRNA in A. flavus, although there were many vsiRNA derived from the AfPV1 genome. GO term and KEGG pathway analyses revealed that the functions of sRNA affected by AfPV1 infection were closely related to vacuole production. These results provide a better understanding of the functional role of RNAi in the impact of AfPV1 on the hypovirulence of A. flavus.

The First Virome of a Russian Vineyard

Among other pathogens, more than 80 viruses infect grapevine. The aim of this work was to study the virome diversity of grapevine viruses and mycoviruses of a vineyard using high-throughput sequencing technologies. The grapevine virome was studied in symptomatic vines of the Rkatsiteli cultivar (V. vinifera) collected at the vineyards of the Krasnodar Krai in Russia. Ribosomal-depleted total RNA and isolated small RNAs were used for library preparation and high-throughput sequencing. Six grapevine-infecting viruses and two viroids were validated by RT-PCR and analyzed phylogenetically. We identified the presence of grapevine leafroll-associated virus 3, grapevine Pinot gris virus, grapevine virus T, grapevine rupestris stem-pitting-associated virus, grapevine fleck virus, and grapevine rupestris vein feathering virus, as well as two viroids, grapevine yellow speckle viroid 1 and hop stunt viroid. We also studied the mycovirome of the vineyard and identified nine viruses with single-stranded positive-sense RNA genomes: alternaria arborescens mitovirus 1, botrytis cinerea mitovirus 1, botrytis cinerea mitovirus 2, botrytis cinerea mitovirus 3, botrytis cinerea mitovirus 4, sclerotinia sclerotiorum mitovirus 3, botrytis cinerea hypovirus 1, grapevine-associated narnavirus 1, and botrytis virus F. In addition, we identified botrytis cinerea hypovirus 1 satellite-like RNA and two single-stranded negative-sense RNA viruses. This is the first study of grapevine mycoviruses in Russia. The obtained result will contribute to the development of biocontrol strategies in the future.

Divergent RNA viruses infecting sea lice, major ectoparasites of fish

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

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.

Viral cross-class transmission results in disease of a phytopathogenic fungus

Interspecies transmission of viruses is a well-known phenomenon in animals and plants whether via contacts or vectors. In fungi, interspecies transmission between distantly related fungi is often suspected but rarely experimentally documented and may have practical implications. A newly described double-strand RNA (dsRNA) virus found asymptomatic in the phytopathogenic fungus Leptosphaeria biglobosa of cruciferous crops was successfully transmitted to an evolutionarily distant, broad-host range pathogen Botrytis cinerea. Leptosphaeria biglobosa botybirnavirus 1 (LbBV1) was characterized in L. biglobosa strain GZJS-19. Its infection in L. biglobosa was asymptomatic, as no significant differences in radial mycelial growth and pathogenicity were observed between LbBV1-infected and LbBV1-free strains. However, cross-species transmission of LbBV1 from L. biglobosa to infection in B. cinerea resulted in the hypovirulence of the recipient B. cinerea strain t-459-V. The cross-species transmission was succeeded only by inoculation of mixed spores of L. biglobosa and B. cinerea on PDA or on stems of oilseed rape with the efficiency of 4.6% and 18.8%, respectively. To investigate viral cross-species transmission between L. biglobosa and B. cinerea in nature, RNA sequencing was carried out on L. biglobosa and B. cinerea isolates obtained from Brassica samples co-infected by these two pathogens and showed that at least two mycoviruses were detected in both fungal groups. These results indicate that cross-species transmission of mycoviruses may occur frequently in nature and result in the phenotypical changes of newly invaded phytopathogenic fungi. This study also provides new insights for using asymptomatic mycoviruses as biocontrol agent.

Interference of Small RNAs in Fusarium graminearum through FgGMTV1 Infection

Small RNA (sRNA) plays a central role in RNA silencing in fungi. The genome of Fusarium graminearum gemytripvirus 1 (FgGMTV1) is comprised of three DNA segments: DNA-A, DNA-B, and DNA-C. DNA-A and DNA-B are associated with fungal growth and virulence reduction. To elucidate the role of RNA silencing during the interactions of fungi and viruses, the sRNA profiles of F. graminearum in association with FgGMTV1 were established, using an FgGMTV1-free library (S-S), a library for infection with the DNA-A and DNA-B segments (S-AB), and a library for infection with the DNA-A, DNA-B, and DNA-C segments (S-ABC). A large amount of virus-derived sRNA (vsiRNA) was detected in the S-AB and S-ABC libraries, accounting for 9.9% and 13.8% of the total sRNA, respectively, indicating that FgGMTV1 triggers host RNA silencing. The total numbers of sRNA reads differed among the three libraries, suggesting that FgGMTV1 infection interferes with host RNA silencing. In addition, the relative proportions of the different sRNA lengths were altered in the S-AB and S-ABC libraries. The genome distribution patterns of the mapping of vsiRNA to DNA-A and DNA-B in the S-AB and S-ABC libraries were also different. These results suggest the influence of DNA-C on host RNA silencing. Transcripts targeted by vsiRNAs were enriched in pathways that included flavin adenine dinucleotide binding, protein folding, and filamentous growth.

Similar Characteristics of siRNAs of Plant Viruses Which Replicate in Plant and Fungal Hosts

Simple Summary

RNA silencing in fungi was shown to confer antiviral defense against plant viruses. In this study, using high-throughput sequencing and bioinformatic analyses, we showed that small interfering RNAs (siRNAs) of cucumber mosaic virus and tobacco mosaic virus (TMV) which replicated in phytopathogenic fungi Rhizoctonia solani and Fusarium graminearum had similarities with viral siRNAs produced in plant hosts in regard to the size distributions, proportion of plus and minus senses, and nucleotide preference for the 5′ termini. Additionally, our results also determined that both F. graminearum DCL1 and DCL2 were involved in the production of TMV siRNAs. Thus, the fungal RNA silencing machineries have adaptive capabilities to recognize and process the genome of invading plant viruses.

Abstract

RNA silencing is a host innate antiviral mechanism which acts via the synthesis of viral-derived small interfering RNAs (vsiRNAs). We have previously reported the infection of phytopathogenic fungi by plant viruses such as cucumber mosaic virus (CMV) and tobacco mosaic virus (TMV). Furthermore, fungal RNA silencing was shown to suppress plant virus accumulation, but the characteristics of plant vsiRNAs associated with the antiviral response in this nonconventional host remain unknown. Using high-throughput sequencing, we characterized vsiRNA profiles in two plant RNA virus–fungal host pathosystems: CMV infection in phytopathogenic fungus Rhizoctonia solani and TMV infection in phytopathogenic fungus Fusarium graminearum. The relative abundances of CMV and TMV siRNAs in the respective fungal hosts were much lower than those in the respective experimental plant hosts, Nicotiana benthamiana and Nicotiana tabacum. However, CMV and TMV siRNAs in fungi had similar characteristics to those in plants, particularly in their size distributions, proportion of plus and minus senses, and nucleotide preference for the 5′ termini of vsiRNAs. The abundance of TMV siRNAs largely decreased in F. graminearum mutants with a deletion in either dicer-like 1 (dcl1) or dcl2 genes which encode key proteins for the production of siRNAs and antiviral responses. However, deletion of both dcl1 and dcl2 restored TMV siRNA accumulation in F. graminearum, indicating the production of dcl-independent siRNAs with no antiviral function in the absence of the dcl1 and dcl2 genes. Our results suggest that fungal RNA silencing recognizes and processes the invading plant RNA virus genome in a similar way as in plants.

Novel RNA viruses from the native range of Hymenoscyphus fraxineus, the causal fungal agent of ash dieback

The native Japanese population of the fungus Hymenoscyphus fraxineus, the causal agent of ash dieback in Europe, was screened for viruses using a high-throughput sequencing method. Five RNA viruses were detected in 116 fungal isolates sequenced via Illumina RNA-seq platform, with an overall virus prevalence of 11.2%. The viruses were completely sequenced by RNA ligase mediated rapid amplification of cDNA ends (RLM-RACE) followed by Sanger sequencing. The sequences appear to represent new species from three established families (Mito-, Endorna- and Partitiviridae), one recognized genus (Botybirnavirus) and a negative-sense single-stranded RNA virus in the order Bunyavirales from the proposed family "Mybuviridae". The highest prevalence was found for the mitovirus (7.8%), that had two genomic forms (linear and circular), while the other viruses were detected each in one isolate. Co-infection of a mitovirus and an endornavirus was also observed in one of the infected isolates. Here we describe the molecular characterization of the identified viruses. This study expands the diversity of viruses in H. fraxineus and provides the basis for investigating the virus-mediated control of ash dieback in Europe.

Mycovirus Diversity and Evolution Revealed/Inferred from Recent Studies

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

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.

Construction and Characterization of a Botrytis Virus F Infectious Clone

Botrytis virus F (BVF) is a positive-sense, single-stranded RNA (+ssRNA) virus within the Gammaflexiviridae family of the plant-pathogenic fungus Botrytis cinerea. In this study, the complete sequence of a BVF strain isolated from B. cinerea collected from grapevine fields in Spain was analyzed. This virus, in this work BVF-V448, has a genome of 6827 nt in length, excluding the poly(A) tail, with two open reading frames encoding an RNA dependent RNA polymerase (RdRP) and a coat protein (CP). The 5′- and 3′-terminal regions of the genome were determined by rapid amplification of cDNA ends (RACE). Furthermore, a yet undetected subgenomic RNA species in BVF-V448 was identified, indicating that the CP is expressed via 3′ coterminal subgenomic RNAs (sgRNAs). We also report the successful construction of the first BVF full-length cDNA clone and synthesized in vitro RNA transcripts using the T7 polymerase, which could efficiently transfect two different strains of B. cinerea, B05.10 and Pi258.9. The levels of growth in culture and virulence on plants of BVF-V448 transfected strains were comparable to BVF-free strains. The infectious clones generated in this work provide a useful tool for the future development of an efficient BVF foreign gene expression vector and a virus-induced gene silencing (VIGS) vector as a biological agent for the control of B. cinerea.

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.

Identification and Molecular Characterization of Novel Mycoviruses in Saccharomyces and Non-Saccharomyces Yeasts of Oenological Interest

Wine yeasts can be natural hosts for dsRNA, ssRNA viruses and retrotransposon elements. In this study, high-throughput RNA sequencing combined with bioinformatic analyses unveiled the virome associated to 16 Saccharomyces cerevisiae and 8 non-Saccharomyces strains of oenological interest. Results showed the presence of six viruses and two satellite dsRNAs from four different families, two of which-Partitiviridae and Mitoviridae-were not reported before in yeasts, as well as two ORFan contigs of viral origin. According to phylogenetic analysis, four new putative mycoviruses distributed in Totivirus, Cryspovirus, and Mitovirus genera were identified. The majority of commercial S. cerevisiae strains were confirmed to be the host for helper L-A type totiviruses and satellite M dsRNAs associated with the killer phenotype, both in single and mixed infections with L-BC totiviruses, and two viral sequences belonging to a new cryspovirus putative species discovered here for the first time. Moreover, single infection by a narnavirus 20S-related sequence was also found in one S. cerevisiae strain. Considering the non-Saccharomyces yeasts, Starmerella bacillaris hosted four RNAs of viral origin-two clustering in Totivirus and Mitovirus genera, and two ORFans with putative satellite behavior. This study confirmed the infection of wine yeasts by viruses associated with useful technological characteristics and demonstrated the presence of complex mixed infections with unpredictable biological effects.

Assessment of mycoviral diversity in Pakistani fungal isolates revealed infection by 11 novel viruses of a single strain of Fusarium mangiferae isolate SP1

An extensive screening survey was conducted on Pakistani filamentous fungal isolates for the identification of viral infections. A total of 396 fungal samples were screened, of which 36 isolates were found double-stranded (ds) RNA positive with an overall frequency of 9% when analysed by a classical dsRNA isolation method. One of 36 dsRNA-positive strains, strain SP1 of a plant pathogenic fungus Fusarium mangiferae, was subjected to virome analysis. Next-generation sequencing and subsequent completion of the entire genome sequencing by a classical Sanger sequencing method showed the SP1 strain to be co-infected by 11 distinct viruses, at least seven of which should be described as new taxa at the species level according to the ICTV (International Committee on Taxonomy of Viruses) species demarcation criteria. The newly identified F. mangiferae viruses (FmVs) include two partitivirids, one betapartitivirus (FmPV1) and one gammapartitivirus (FmPV2); six mitovirids, three unuamitovirus (FmMV2, FmMV4, FmMV6), one duamitovirus (FmMV5), and two unclassified mitovirids (FmMV1, FmMV3); and three botourmiavirids, two magoulivirus (FmBOV1, FmBOV3) and one scleroulivirus (FmBOV2). The number of coinfecting viruses is among the largest ones of fungal coinfections. Their molecular features are thoroughly described here. This represents the first large virus survey in the Indian sub-continent.

Mycoviruses and their role in fungal pathogenesis

Nowadays, the focus of mycovirology research has expanded from plant pathogenic fungi and mycovirus mediated hypovirulence to include insect and human pathogenic fungi together with a range of mycovirus mediated phenotypes, such as hypervirulence, control of endophytic traits, regulation of metabolite production and drug resistance. In fungus-mycovirus-environmental interactions, the environment and both abiotic and biotic factors play crucial roles in whether and how mycovirus mediated phenotypes are manifest. Mycovirus infections result in alterations in the host transcriptome profile, via protein-protein interactions and triggering of antiviral RNA silencing in the fungus. These alterations, in combination with the environmental factors, may result in desirable phenotypic traits for the host, for us and in some cases for both.

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.

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.

Residual Effects Caused by a Past Mycovirus Infection in Fusarium circinatum

Mycoviruses are known to be difficult to cure in fungi but their spontaneous loss occurs commonly. The unexpected disappearance of mycoviruses can be explained by diverse reasons, from methodological procedures to biological events such as posttranscriptional silencing machinery. The long-term effects of a virus infection on the host organism have been well studied in the case of human viruses; however, the possible residual effect on a fungus after the degradation of a mycovirus is unknown. For that, this study analyses a possible residual effect on the transcriptome of the pathogenic fungus Fusarium circinatum after the loss of the mitovirus FcMV1. The mycovirus that previously infected the fungal isolate was not recovered after a 4-year storage period. Only 14 genes were determined as differentially expressed and were related to cell cycle regulation and amino acid metabolism. The results showed a slight acceleration in the metabolism of the host that had lost the mycovirus by the upregulation of the genes involved in essential functions for fungal development. The analysis also revealed a weak expression in the annotated genes of the RNA silencing machinery. To our knowledge, this is the first time that a potential residual effect on the host transcriptome caused by the past infection of a mycovirus is reported.

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.

The virome from a collection of endomycorrhizal fungi reveals new viral taxa with unprecedented genome organization

Mutualistic plant-associated fungi are recognized as important drivers in plant evolution, diversity, and health. The discovery that mycoviruses can take part and play important roles in symbiotic tripartite interactions has prompted us to study the viromes associated with a collection of ericoid and orchid mycorrhizal (ERM and ORM, respectively) fungi. Our study, based on high-throughput sequencing of transcriptomes (RNAseq) from fungal isolates grown in axenic cultures, revealed in both ERM and ORM fungi the presence of new mycoviruses closely related to already classified virus taxa, but also new viruses that expand the boundaries of characterized RNA virus diversity to previously undescribed evolutionary trajectories. In ERM fungi, we provide first evidence of a bipartite virus, distantly related to narnaviruses, that splits the RNA-dependent RNA polymerase (RdRP) palm domain into two distinct proteins, encoded by each of the two segments. Furthermore, in one isolate of the ORM fungus Tulasnella spp. we detected a 12 kb genomic fragment coding for an RdRP with features of bunyavirus-like RdRPs. However, this 12 kb genomic RNA has the unique features, for Bunyavirales members, of being tri-cistronic and carrying ORFs for the putative RdRP and putative nucleocapsid in ambisense orientation on the same genomic RNA. Finally, a number of ORM fungal isolates harbored a group of ambisense bicistronic viruses with a genomic size of around 5 kb, where we could identify a putative RdRP palm domain that has some features of plus strand RNA viruses; these new viruses may represent a new lineage in the Riboviria, as they could not be reliably assigned to any of the branches in the recently derived monophyletic tree that includes most viruses with an RNA genome.

The dsRNA mycovirus ChNRV1 causes mild hypervirulence in the fungal phytopathogen Colletotrichum higginsianum

The genus Colletotrichum comprises a large number of filamentous fungi responsible for anthracnose diseases in many tropical and subtropical fruits and vegetables. In particular, Colletotrichum higginsianum infects Brassicaceae species, including Arabidopsis. The C. higginsianum strain IMI349063A is naturally infected with a dsRNA virus, named Colletorichum higginsianum non-segmented virus (ChNRV1). Here, we investigated the biological effect of ChNRV1 in C. higginsianum by comparing strains with and without the virus. ChNRV1 does not have an effect on C. higginsianum growth under salt and cell-wall stress conditions. However, thermal stress reduced C. higginsianum growth rate, this effect being more evident in the wild-type C. higginsianum strain containing the virus. Although ChNRV1 had no effect in conidiation, conidia were narrower when the virus is present. More importantly, ChNRV1 causes a mild increase in C. higginsianum virulence (hypervirulence) when infecting Arabidopsis plants. These findings indicated that, whereas the ChNRV1 mycovirus does not impair growth and conidiation of C. higginsianum, it confers hypervirulence to the fungal host. These findings will help in future research on the effect of mycoviral infection on pathogenic fungi in plant species of agronomical relevance.

The RNA Virome and Its Dynamics in an Invasive Fruit Fly, Bactrocera dorsalis, Imply Interactions Between Host and Viruses

The oriental fruit fly, Bactrocera dorsalis (Hendel), is an important invasive agricultural insect pest with a wide host range, and has spread around the world over the last century. This evolutionary trait may have arisen primarily from interactions between B. dorsalis and other invertebrates that share the same ecological niches. The invasive behavior of B. dorsalis also frequently exposes them to diverse species of viruses. Thereby, RNA viromes may be useful microbial markers to understand the ecological evolution of B. dorsalis as well as to investigate virus-host interactions. Here, we reported eight novel RNA viruses in B. dorsalis of a lab colony, including four positive-strand RNA viruses, two negative-strand RNA viruses, and two double-stranded RNA viruses using high-throughput sequencing technology. Analysis of the virus-derived small RNAs suggested that most of these viruses may be active and trigger the host antiviral RNAi responses. The viruses were also detected in various geographical populations of B. dorsalis, implying that there is a strong association between the viromes and host. In addition, these viruses infected specific fly tissues, predominately the central nervous system and gut. Furthermore, we explored the dynamics of the viruses when hosts were exposed to short- or long-term stressors, which showed that titers of some viruses were responsively altered in the stressed B. dorsalis. The discovery of these viruses may enrich our understanding of the species diversity of RNA viruses and also provide information on viruses in association with host adaptation in insects.

sRIS: A Small RNA Illustration System for Plant Next-Generation Sequencing Data Analysis

Small RNA (sRNA), such as microRNA (miRNA) and short interfering RNA, are well-known to control gene expression based on degradation of target mRNA in plants. A considerable amount of research has applied next-generation sequencing (NGS) to reveal the regulatory pathways of plant sRNAs. Consequently, numerous bioinformatics tools have been developed for the purpose of analyzing sRNA NGS data. However, most methods focus on the study of sRNA expression profiles or novel miRNAs predictions. The analysis of sRNA target genes is usually not integrated into their pipelines. As a result, there is still no means available for identifying the interaction mechanisms between host and virus or the synergistic effects between two viruses. For the present study, a comprehensive system, called the Small RNA Illustration System (sRIS), has been developed. This system contains two main components. The first is for sRNA overview analysis and can be used not only to identify miRNA but also to investigate virus-derived small interfering RNA. The second component is for sRNA target prediction, and it employs both bioinformatics calculations and degradome sequencing data to enhance the accuracy of target prediction. In addition, this system has been designed so that figures and tables for the outputs of each analysis can be easily retrieved and accessed, making it easier for users to quickly identify and quantify their results. sRIS is available at http://sris.itps.ncku.edu.tw/.

Movement of small RNAs in and between plants and fungi

RNA interference is a biological process whereby small RNAs inhibit gene expression through neutralizing targeted mRNA molecules. This process is conserved in eukaryotes. Here, recent work regarding the mechanisms of how small RNAs move within and between organisms is examined. Small RNAs can move locally and systemically in plants through plasmodesmata and phloem, respectively. In fungi, transportation of small RNAs may also be achieved by septal pores and vesicles. Recent evidence also supports bidirectional cross-kingdom communication of small RNAs between host plants and adapted fungal pathogens to affect the outcome of infection. We discuss several mechanisms for small RNA trafficking and describe evidence for transport through naked form, combined with RNA-binding proteins or enclosed by vesicles.

Different Genetic Sources Contribute to the Small RNA Population in the Arbuscular Mycorrhizal Fungus Gigaspora margarita

RNA interference (RNAi) is a key regulatory pathway of gene expression in almost all eukaryotes. This mechanism relies on short non-coding RNA molecules (sRNAs) to recognize in a sequence-specific manner DNA or RNA targets leading to transcriptional or post-transcriptional gene silencing. To date, the fundamental role of sRNAs in the regulation of development, stress responses, defense against viruses and mobile elements, and cross-kingdom interactions has been extensively studied in a number of biological systems. However, the knowledge of the “RNAi world” in arbuscular mycorrhizal fungi (AMF) is still limited. AMF are obligate mutualistic endosymbionts of plants, able to provide several benefits to their partners, from improved mineral nutrition to stress tolerance. Here we described the RNAi-related genes of the AMF Gigaspora margarita and characterized, through sRNA sequencing, its complex small RNAome, considering the possible genetic sources and targets of the sRNAs. G. margarita indeed is a mosaic of different genomes since it hosts endobacteria, RNA viruses, and non-integrated DNA fragments corresponding to mitovirus sequences. Our findings show that G. margarita is equipped with a complete set of RNAi-related genes characterized by the expansion of the Argonaute-like (AGO-like) gene family that seems a common trait of AMF. With regards to sRNAs, we detected populations of sRNA reads mapping to nuclear, mitochondrial, and viral genomes that share similar features (25-nt long and 5′-end uracil read enrichments), and that clearly differ from sRNAs of endobacterial origin. Furthermore, the annotation of nuclear loci producing sRNAs suggests the occurrence of different sRNA-generating processes. In silico analyses indicate that the most abundant G. margarita sRNAs, including those of viral origin, could target transcripts in the host plant, through a hypothetical cross-kingdom RNAi.

Virome Identification and Characterization of Fusarium sacchari and F. andiyazi: Causative Agents of Pokkah Boeng Disease in Sugarcane

Fusarium sacchari and Fusarium andiyazi are two devastating sugarcane pathogens that cause pokkah boeng disease (PBD) in China. RNA_Seq was conducted to identify mycoviruses in F. sacchari and F. andiyazi isolates collected from PBD symptom-showing sugarcane plants across China. Fifteen isolates with a normal, debilitated, or abnormal phenotype in colony morphology were screened out for the existence of dsRNA from 104 Fusarium isolates. By sequencing the mixed pool of dsRNA from these Fusarium isolates, a total of 26 contigs representing complete or partial genome sequences of ten mycoviruses and their strains were identified, including one virus belonging to Hypoviridae, two mitoviruses with seven strains belonging to Narnaviridae, one virus of Chrysoviridae, and one alphavirus-like virus. RT-PCR amplification with primers specific to individual mycoviruses revealed that mitoviruses were the most prevalent and the alphavirus-like virus and chrysovirus were the least prevalent. In terms of host preference, more mitoviruses were found in F. andiyazi than in F. sacchari. Fusarium sacchari hypovirus 1 with a 13.9 kb genome and a defective genome of 12.2 kb, shares 54% identity at the amino acid level to the Wuhan insect virus 14, which is an unclassified hypovirus identified from insect meta-transcriptomics. The alphavirus-like virus, Fusarium sacchari alphavirus-like virus 1 (FsALV1), seemed to hold a distinct status amid fungal alphavirus-like viruses, with the highest identity of 27% at the amino acid level to Sclerotium rolfsii alphavirus-like virus 3 and 29% to a hepevirus, Ferret hepatitis E virus. While six of the seven mitoviruses shared 72–94% identities to known mitoviruses, Fusarium andiyazi mitovirus 2 was most similar to Alternaria brassicicola mitovirus with an identity of only 49% between the two viruses. Transmission of FsALV1 and Fusarium sacchari chrysovirus 1 (FsCV1) from F. sacharri to F. commune was observed and the characterization of the four-segment dsRNA chrysovirus was performed with aid of electron microscopy and analysis of the encapsidated RNAs. These findings provide insight into the diversity and spectrum of mycoviruses in PBD pathogens and should be useful for exploring agents to control the disease.

Extreme Diversity of Mycoviruses Present in Isolates of Rhizoctonia solani AG2-2 LP From Zoysia japonica From Brazil

Zoysia japonica, in Brazil, is commonly infected by Rhizoctonia solani (R. solani) in humid and cool weather conditions. Eight isolates of R. solani, previously identified as belonging to the AG2-2 LP anastomosis group, isolated from samples from large path symptoms, were collected from three counties in São Paulo state (Brazil) and investigated for the presence of mycoviruses. After detection of double-strand RNA (dsRNA) in all samples, RNA_Seq analysis of ribosomal RNA-depleted total RNA from in vitro cultivated mycelia was performed. Forty-seven partial or complete viral unique RNA dependent-RNA polymerase (RdRp) sequences were obtained with a high prevalence of positive sense ssRNA viruses. Sequences were sufficiently different from the first match in BLAST searches suggesting that they all qualify as possible new viral species, except for one sequence showing an almost complete match with Rhizoctonia solani dsRNA virus 2, an alphapartitivirus. Surprisingly four large contigs of putative viral RNA could not be assigned to any existing clade of viruses present in the databases, but no DNA was detected corresponding to these fragments confirming their viral replicative nature. This is the first report on the occurrence of mycoviruses in R. solani AG2-2 LP in South America.

Discovery of Two Mycoviruses by High-Throughput Sequencing and Assembly of Mycovirus-Derived Small Silencing RNAs From a Hypovirulent Strain of Sclerotinia sclerotiorum

Sclerotinia sclerotiorum, an important phytopathogenic fungus, harbors rich diversity of mycoviruses. Lately, more mycoviruses can be successfully and accurately discovered by deep sequencing, especially those that could not be detected by traditional double-stranded RNA (dsRNA) extraction. Previously, we reported that the hypovirulent S. sclerotiorum strain SZ-150 is coinfected by Sclerotinia sclerotiorum hypovirus 1 (SsHV1) and its related satellite RNA. Here, aside from SsHV1, we detected two other mycoviruses, Sclerotinia sclerotiorum botybirnavirus 3 (SsBV3/SZ-150) and Sclerotinia sclerotiorum mycotymovirus 1 (SsMTV1/SZ-150), coinfecting strain SZ-150, by deep sequencing and assembly of mycovirus-derived small RNAs and determined their full-length genomes. The genome of SsBV3/SZ-150 was found to be composed of two linear dsRNA segments, 6,212, and 5,880 bp in size, respectively. Each dsRNA segment of SsBV3/SZ-150 contains a large open reading frame (ORF) encoding RNA-dependent RNA polymerase (RdRp) and a hypothetical protein. The whole genome of SsBV3/SZ-150 shares more than 95% sequence identity with Botrytis porri botybirnavirus 1 (BpBV1) at the nucleotide (nt) or amino acid level. Thus, SsBV3/SZ-150 was assumed to be a strain of BpBV1. The genome of SsMTV1/SZ-150 consists of 6,391 nt excluding the poly(A) tail. SsMTV1/SZ-150 was predicted to contain a large ORF that encodes a putative replication-associated polyprotein (RP) with three conserved domains of viral RNA methyltransferase, viral RNA helicase, and RdRp. Phylogenetic analyses suggest that SsMTV1/SZ-150 is related, albeit distantly, to members of the family Tymoviridae. Analysis of the small RNAs derived from SsBV3/SZ-150 and SsMTV1/SZ-150 revealed that small-RNA lengths mainly range from 20 to 24 nt, with a peak at 22 nt, and the most abundant 5′-terminal nucleotide is uridine, suggesting that the Dicer 2 and Argonaute 1, two key components in the RNA inference pathway, may play important roles in the resistance to mycoviral infection in S. sclerotiorum. Neither SsBV3/SZ-150 nor SsMTV1/SZ-150 is a causal agent of hypovirulence in strain SZ-150.

Biological and Molecular Characterization of Chenopodium quinoa Mitovirus 1 Reveals a Distinct Small RNA Response Compared to Those of Cytoplasmic RNA Viruses

Indirect evidence of mitochondrial viruses in plants comes from discovery of genomic fragments integrated into the nuclear and mitochondrial DNA of a number of plant species. Here, we report the existence of replicating mitochondrial virus in plants: from transcriptome sequencing (RNA-seq) data of infected Chenopodium quinoa, a plant species commonly used as a test plant in virus host range experiments, among other virus contigs, we could assemble a 2.7-kb contig that had highest similarity to mitoviruses found in plant genomes. Northern blot analyses confirmed the existence of plus- and minus-strand RNA corresponding to the mitovirus genome. No DNA corresponding to the genomic RNA was detected, excluding the endogenization of such virus. We have tested a number of C. quinoa accessions, and the virus was present in a number of commercial varieties but absent from a large collection of Bolivian and Peruvian accessions. The virus could not be transmitted mechanically or by grafting, but it is transmitted vertically through seeds at a 100% rate. Small RNA analysis of a C. quinoa line carrying the mitovirus and infected by alfalfa mosaic virus showed that the typical antiviral silencing response active against cytoplasmic viruses (21- to 22-nucleotide [nt] vsRNA peaks) is not active against CqMV1, since in this specific case the longest accumulating vsRNA length is 16 nt, which is the same as that corresponding to RNA from mitochondrial genes. This is evidence of a distinct viral RNA degradation mechanism active inside mitochondria that also may have an antiviral effect.IMPORTANCE This paper reports the first biological characterization of a bona fide plant mitovirus in an important crop, Chenopodium quinoa, providing data supporting that mitoviruses have the typical features of cryptic (persistent) plant viruses. We, for the first time, demonstrate that plant mitoviruses are associated with mitochondria in plants. In contrast to fungal mitoviruses, plant mitoviruses are not substantially affected by the antiviral silencing pathway, and the most abundant mitovirus small RNA length is 16 nt.

Investigation of Host Range of and Host Defense against a Mitochondrially Replicating Mitovirus

Mitoviruses (genus Mitovirus, family Narnaviridae) are mitochondrially replicating viruses that have the simplest positive-sense RNA genomes of 2.2 to 4.4 kb with a single open reading frame (ORF) encoding an RNA-dependent RNA polymerase. Cryphonectria parasitica mitovirus 1 (CpMV1) from U.S. strain NB631 of the chestnut blight fungus, Cryphonectria parasitica, was the first virus identified as a mitochondrially replicating virus. Despite subsequent discovery of many other mitoviruses from diverse fungi, no great advances in understanding mitovirus biology have emerged, partly because of the lack of inoculation methods. Here we developed a protoplast fusion-based protocol for horizontal transmission of CpMV1 that entailed fusion of recipient and donor protoplasts, hyphal anastomosis, and single-conidium isolation. This method allowed expansion of the host range to many other C. parasitica strains. Species within and outside the family Cryphonectriaceae, Cryphonectria radicalis and Valsa ceratosperma, also supported the replication of CpMV1 at a level comparable to that in the natural host. No stable maintenance of CpMV1 was observed in Helminthosporium victoriae PCR-based haplotyping of virus-infected fungal strains confirmed the recipient mitochondrial genetic background. Phenotypic comparison between CpMV1-free and -infected isogenic strains revealed no overt effects of the virus. Taking advantage of the infectivity to the standard strain C. parasitica EP155, accumulation levels were compared among antiviral RNA silencing-proficient and -deficient strains in the EP155 background. Comparable accumulation levels were observed among these strains, suggesting the avoidance of antiviral RNA silencing by CpMV1, which is consistent with its mitochondrial replication. Collectively, the results of study provide a foundation to further explore the biology of mitoviruses.IMPORTANCE Capsidless mitoviruses, which are ubiquitously detected in filamentous fungi, have the simplest RNA genomes of 2.2 to 4.4 kb, encoding only RNA-dependent RNA polymerase. Despite their simple genomes, detailed biological characterization of mitoviruses has been hampered by their mitochondrial location within the cell, posing challenges to their experimental introduction and study. Here we developed a protoplast fusion-based protocol for horizontal transfer of the prototype mitovirus, Cryphonectria parasitica mitovirus 1 (CpMV1), which was isolated from strain NB631 of the chestnut blight fungus (Cryphonectria parasitica), a model filamentous fungus for studying virus-host interactions. The host range of CpMV1 has been expanded to many different strains of C. parasitica and different fungal species within and outside the Cryphonectriaceae. Comparison of CpMV1 accumulation among various RNA silencing-deficient and -competent strains showed clearly that the virus was unaffected by RNA silencing. This study provides a solid foundation for further exploration of mitovirus-host interactions.

Small RNA-Seq to Characterize Viruses Responsible of Lettuce Big Vein Disease in Spain

The emerging lettuce big-vein disease (LBVD) is causing losses in lettuce production ranging from 30 to 70% worldwide. Several studies have associated this disease with Mirafiori lettuce big-vein virus (MiLBVV) alone or in mixed infection with lettuce big-vein associated virus (LBVaV). We used Illumina small RNA sequencing (sRNA-seq) to identify viruses present in symptomatic lettuce plants from commercial fields in Southern Spain. Data analysis using the VirusDetect tool showed the consistent presence of MiLBVV and LBVaV in diseased plants. Populations of MiLBVV and LBVaV viral small RNAs (sRNAs) were characterized, showing features essentially similar to those of other viruses, with the peculiarity of an uneven asymmetric distribution of MiLBVV virus-derived small RNAs (vsRNAs) for the different polarities of genomic RNA4 vs. RNAs1 to 3. Sanger sequencing of coat protein genes was used to study MiLBVV and LBVaV phylogenetic relationships and population genetics. The Spanish MiLBVV population was composed of isolates from three well-differentiated lineages and reflected almost all of the diversity reported for the MiLBVV species, whereas the LBVaV population showed very little genetic differentiation at the regional scale but lineage differentiation at a global geographical scale. Universal primers were used to detect and quantify the accumulation of MiLBVV and LBVaV in field samples; both symptomatic and asymptomatic plants from affected fields carried equal viral loads, with LBVaV accumulating at higher levels than MiLBVV.

Diverse, Novel Mycoviruses From the Virome of a Hypovirulent Sclerotium rolfsii Strain

Sclerotium rolfsii, which causes southern blight in a wide variety of crops, is a devastating plant pathogen worldwide. Mycoviruses that induce hypovirulence in phytopathogenic fungi are potential biological control resources against fungal plant diseases. However, in S. rolfsii, mycoviruses are rarely reported. In a previous study, we found a hypovirulent strain carrying a diverse pattern of dsRNAs. Here, we utilized the RNA_Seq technique to detect viral sequences. Deep sequencing, RT-PCR and Sanger sequencing validation analyses revealed that this strain harbors various new viral species that show affinity to the distinctly established and proposed families Benyviridae, Endornaviridae, Fusariviridae, Hypoviridae, and Fusagraviridae. Moreover, some viral sequences that could not be assigned to any of the existing families or orders were also identified and showed similarities to the Alphavirus, Ourmiavirus, phlegivirus-like and Curvularia thermal tolerance virus-like groups. In addition, we also conducted deep sequencing analysis of small RNAs in the virus-infecting fugal strain. The results indicated that the Dicer-mediated gene silencing mechanism was present in S. rolfsii. This is the first report of viral diversity in a single S. rolfsii fungal strain, and the results presented herein might provide insight into the taxonomy and evolution of mycoviruses and be useful for the exploration of mycoviruses as biocontrol agents.

Small RNA-Omics for Plant Virus Identification, Virome Reconstruction, and Antiviral Defense Characterization

RNA interference (RNAi)-based antiviral defense generates small interfering RNAs that represent the entire genome sequences of both RNA and DNA viruses as well as viroids and viral satellites. Therefore, deep sequencing and bioinformatics analysis of small RNA population (small RNA-ome) allows not only for universal virus detection and genome reconstruction but also for complete virome reconstruction in mixed infections. Viral infections (like other stress factors) can also perturb the RNAi and gene silencing pathways regulating endogenous gene expression and repressing transposons and host genome-integrated endogenous viral elements which can potentially be released from the genome and contribute to disease. This review describes the application of small RNA-omics for virus detection, virome reconstruction and antiviral defense characterization in cultivated and non-cultivated plants. Reviewing available evidence from a large and ever growing number of studies of naturally or experimentally infected hosts revealed that all families of land plant viruses, their satellites and viroids spawn characteristic small RNAs which can be assembled into contigs of sufficient length for virus, satellite or viroid identification and for exhaustive reconstruction of complex viromes. Moreover, the small RNA size, polarity and hotspot profiles reflect virome interactions with the plant RNAi machinery and allow to distinguish between silent endogenous viral elements and their replicating episomal counterparts. Models for the biogenesis and functions of small interfering RNAs derived from all types of RNA and DNA viruses, satellites and viroids as well as endogenous viral elements are presented and discussed.

Molecular characterization of a novel double-stranded RNA mycovirus of Trichoderma asperellum strain JLM45-3

In this study, we describe a novel mycovirus isolated from Trichoderma asperellum, which was designated Trichoderma asperellum dsRNA Virus 1 (TaRV1). The sequence analysis revealed that TaRV1 has two discontinuous open reading frames (ORF), ORF1 and ORF2. A hypothetical protein and an RNA-dependent RNA polymerase are encoded by ORF1 and ORF2, respectively. Phylogenetic analysis based on RdRp sequences clearly places TaRV1 in a taxonomically unassigned dsRNA mycovirus group.

Genome characterization of a divergent isolate of the mycovirus Botrytis virus F from a grapevine metagenome

As part of a grapevine metagenome study, total RNA extracted from grapevine phloem scrapings was analyzed by Illumina sequencing. For one 420A rootstock sample, reads mapping against a reference database and BLAST annotation of contigs identified the presence of a divergent isolate of Botrytis virus F (BVF). The full genome sequence of this isolate (IVC-5-77) was determined (6,828 nucleotides [nt], excluding the poly(A) tail) and shown to be collinear with that of the BVF reference isolate, with the two open reading frames encoding a replication-associated protein (REP) and a coat protein (CP). The IVC-5-77 isolate, however, is very divergent, showing only 81.3-81.6% nucleotide sequence identity to the two other sequenced BVF isolates. The internal non-coding region was also found to be highly variable between BVF isolates. Analysis of the RNASeq reads demonstrated that close to 20% of them belong to Botrytis cinerea, the putative host of the IVC-5-77 isolate. These results extend our knowledge of the diversity and variability of BVF and demonstrate its detectability, together with that of its B. cinerea host, in total RNA RNASeq data from grapevine phloem scrapings.

Metagenomic sequencing suggests a diversity of RNA interference-like responses to viruses across multicellular eukaryotes

RNA interference (RNAi)-related pathways target viruses and transposable element (TE) transcripts in plants, fungi, and ecdysozoans (nematodes and arthropods), giving protection against infection and transmission. In each case, this produces abundant TE and virus-derived 20-30nt small RNAs, which provide a characteristic signature of RNAi-mediated defence. The broad phylogenetic distribution of the Argonaute and Dicer-family genes that mediate these pathways suggests that defensive RNAi is ancient, and probably shared by most animal (metazoan) phyla. Indeed, while vertebrates had been thought an exception, it has recently been argued that mammals also possess an antiviral RNAi pathway, although its immunological relevance is currently uncertain and the viral small RNAs (viRNAs) are not easily detectable. Here we use a metagenomic approach to test for the presence of viRNAs in five species from divergent animal phyla (Porifera, Cnidaria, Echinodermata, Mollusca, and Annelida), and in a brown alga-which represents an independent origin of multicellularity from plants, fungi, and animals. We use metagenomic RNA sequencing to identify around 80 virus-like contigs in these lineages, and small RNA sequencing to identify viRNAs derived from those viruses. We identified 21U small RNAs derived from an RNA virus in the brown alga, reminiscent of plant and fungal viRNAs, despite the deep divergence between these lineages. However, contrary to our expectations, we were unable to identify canonical (i.e. Drosophila- or nematode-like) viRNAs in any of the animals, despite the widespread presence of abundant micro-RNAs, and somatic transposon-derived piwi-interacting RNAs. We did identify a distinctive group of small RNAs derived from RNA viruses in the mollusc. However, unlike ecdysozoan viRNAs, these had a piRNA-like length distribution but lacked key signatures of piRNA biogenesis. We also identified primary piRNAs derived from putatively endogenous copies of DNA viruses in the cnidarian and the echinoderm, and an endogenous RNA virus in the mollusc. The absence of canonical virus-derived small RNAs from our samples may suggest that the majority of animal phyla lack an antiviral RNAi response. Alternatively, these phyla could possess an antiviral RNAi response resembling that reported for vertebrates, with cryptic viRNAs not detectable through simple metagenomic sequencing of wild-type individuals. In either case, our findings show that the antiviral RNAi responses of arthropods and nematodes, which are highly divergent from each other and from that of plants and fungi, are also highly diverged from the most likely ancestral metazoan state.

Two Novel Hypovirulence-Associated Mycoviruses in the Phytopathogenic Fungus Botrytis cinerea: Molecular Characterization and Suppression of Infection Cushion Formation

Botrytis cinerea is a necrotrophic fungus causing disease on many important agricultural crops. Two novel mycoviruses, namely Botrytis cinerea hypovirus 1 (BcHV1) and Botrytis cinerea fusarivirus 1 (BcFV1), were fully sequenced. The genome of BcHV1 is 10,214 nt long excluding a poly-A tail and possesses one large open reading frame (ORF) encoding a polyprotein possessing several conserved domains including RNA-dependent RNA polymerase (RdRp), showing homology to hypovirus-encoded polyproteins. Phylogenetic analysis indicated that BcHV1 may belong to the proposed genus Betahypovirus in the viral family Hypoviridae. The genome of BcFV1 is 8411 nt in length excluding the poly A tail and theoretically processes two major ORFs, namely ORF1 and ORF2. The larger ORF1 encoded polypeptide contains protein domains of an RdRp and a viral helicase, whereas the function of smaller ORF2 remains unknown. The BcFV1 was phylogenetically clustered with other fusariviruses forming an independent branch, indicating BcFV1 was a member in Fusariviridae. Both BcHV1 and BcFV1 were capable of being transmitted horizontally through hyphal anastomosis. Infection by BcHV1 alone caused attenuated virulence without affecting mycelial growth, significantly inhibited infection cushion (IC) formation, and altered expression of several IC-formation-associated genes. However, wound inoculation could fully rescue the virulence phenotype of the BcHV1 infected isolate. These results indicate the BcHV1-associated hypovirulence is caused by the viral influence on IC-formation-associated pathways.

Characterization of small RNAs originating from mitoviruses infecting the conifer pathogen Fusarium circinatum

Deep sequencing of small RNAs has proved effective in the diagnosis of mycovirus infections. In this study, the presence of mycoviruses in ten isolates of the phytopathogenic fungus Fusarium circinatum was investigated by high-throughput sequencing (HTS) of small RNAs. The contigs resulting from de novo assembly of the reads were aligned to viral genome sequences. The presence of each mycovirus detected in the isolates was confirmed by RT-PCR analysis with four previously described primer pairs and seven new pairs designed on the basis of sequencing data. The findings demonstrate the potential use of HTS for reconstructing previously identified mitoviruses infecting F. circinatum.

Profile and functional analysis of small RNAs derived from Aspergillus fumigatus infected with double-stranded RNA mycoviruses

Background

Mycoviruses are viruses that naturally infect and replicate in fungi. Aspergillus fumigatus, an opportunistic pathogen causing fungal lung diseases in humans and animals, was recently shown to harbour several different types of mycoviruses. A well-characterised defence against virus infection is RNA silencing. The A. fumigatus genome encodes essential components of the RNA silencing machinery, including Dicer, Argonaute and RNA-dependent RNA polymerase (RdRP) homologues. Active silencing of double-stranded (ds)RNA and the generation of small RNAs (sRNAs) has been shown for several mycoviruses and it is anticipated that a similar mechanism will be activated in A. fumigatus isolates infected with mycoviruses.

Results

To investigate the existence and nature of A. fumigatus sRNAs, sRNA-seq libraries of virus-free and virus-infected isolates were created using Scriptminer adapters and compared. Three dsRNA viruses were investigated: Aspergillus fumigatus partitivirus-1 (AfuPV-1, PV), Aspergillus fumigatus chrysovirus (AfuCV, CV) and Aspergillus fumigatus tetramycovirus-1 (AfuTmV-1, NK) which were selected because they induce phenotypic changes such as coloration and sectoring. The dsRNAs of all three viruses, which included two conventionally encapsidated ones PV and CV and one unencapsidated example NK, were silenced and yielded characteristic vsiRNAs together with co-incidental silencing of host fungal genes which shared sequence homology with the viral genomes.

Conclusions

Virus-derived sRNAs were detected and characterised in the presence of virus infection. Differentially expressed A. fumigatus microRNA-like (miRNA-like) sRNAs and small interfering RNAs (siRNAs) were detected and validated. Host sRNA loci which were differentially expressed as a result of virus infection were also identified. To our knowledge, this is the first study reporting the sRNA profiles of A. fumigatus isolates.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3773-8) contains supplementary material, which is available to authorized users.