Characterization of hypovirus-derived small RNAs generated in the chestnut blight fungus by an inducible DCL-2-dependent pathway

Argonaute-independent, Dicer-dependent antiviral defense against RNA viruses

Antiviral RNA interference (RNAi) is conserved from yeasts to mammals. Dicer recognizes and cleaves virus-derived double-stranded RNA (dsRNA) and/or structured single-stranded RNA (ssRNA) into small-interfering RNAs, which guide effector Argonaute to homologous viral RNAs for digestion and inhibit virus replication. Thus, Argonaute is believed to be essential for antiviral RNAi. Here, we show Argonaute-independent, Dicer-dependent antiviral defense against dsRNA viruses using Cryphonectria parasitica (chestnut blight fungus), which is a model filamentous ascomycetous fungus and hosts a variety of viruses. The fungus has two dicer-like genes (dcl1 and dcl2) and four argonaute-like genes (agl1 to agl4). We prepared a suite of single to quadruple agl knockout mutants with or without dcl disruption. We tested these mutants for antiviral activities against diverse dsRNA viruses and ssRNA viruses. Although both DCL2 and AGL2 worked as antiviral players against some RNA viruses, DCL2 without argonaute was sufficient to block the replication of other RNA viruses. Overall, these results indicate the existence of a Dicer-alone defense and different degrees of susceptibility to it among RNA viruses. We discuss what determines the great difference in susceptibility to the Dicer-only defense.

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.

Mycovirus-encoded suppressors of RNA silencing: Possible allies or enemies in the use of RNAi to control fungal disease in crops

Plants, fungi, and many other eukaryotes have evolved an RNA interference (RNAi) mechanism that is key for regulating gene expression and the control of pathogens. RNAi inhibits gene expression, in a sequence-specific manner, by recognizing and deploying cognate double-stranded RNA (dsRNA) either from endogenous sources (e.g. pre-micro RNAs) or exogenous origin (e.g. viruses, dsRNA, or small interfering RNAs, siRNAs). Recent studies have demonstrated that fungal pathogens can transfer siRNAs into plant cells to suppress host immunity and aid infection, in a mechanism termed cross-kingdom RNAi. New technologies, based on RNAi are being developed for crop protection against insect pests, viruses, and more recently against fungal pathogens. One example, is host-induced gene silencing (HIGS), which is a mechanism whereby transgenic plants are modified to produce siRNAs or dsRNAs targeting key transcripts of plants, or their pathogens or pests. An alternative gene regulation strategy that also co-opts the silencing machinery is spray-induced gene silencing (SIGS), in which dsRNAs or single-stranded RNAs (ssRNAs) are applied to target genes within a pathogen or pest. Fungi also use their RNA silencing machinery against mycoviruses (fungal viruses) and mycoviruses can deploy virus-encoded suppressors of RNAi (myco-VSRs) as a counter-defence. We propose that myco-VSRs may impact new dsRNA-based management methods, resulting in unintended outcomes, including suppression of management by HIGS or SIGS. Despite a large diversity of mycoviruses being discovered using high throughput sequencing, their biology is poorly understood. In particular, the prevalence of mycoviruses and the cellular effect of their encoded VSRs are under-appreciated when considering the deployment of HIGS and SIGS strategies. This review focuses on mycoviruses, their VSR activities in fungi, and the implications for control of pathogenic fungi using RNAi.

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.

A Satellite dsRNA Attenuates the Induction of Helper Virus-Mediated Symptoms in Aspergillus flavus

Aspergillus flavus is an important fungal pathogen of animals and plants. Previously, we reported a novel partitivirus, Aspergillus flavus partitivirus 1 (AfPV1), infecting A. flavus. In this study, we obtained a small double-stranded (ds) RNA segment (734 bp), which is a satellite RNA of the helper virus, AfPV1. The presence of AfPV1 altered the colony morphology, decreased the number of conidiophores, created significantly larger vacuoles, and caused more sensitivity to osmotic, oxidative, and UV stresses in A. flavus, but the small RNA segment could attenuate the above symptoms caused by the helper virus AfPV1 in A. flavus. Moreover, AfPV1 infection reduced the pathogenicity of A. flavus in corn (Zea mays), honeycomb moth (Galleria mellonella), mice (Mus musculus), and the adhesion of conidia to host epithelial cells, and increased conidial death by macrophages. However, the small RNA segment could also attenuate the above symptoms caused by the helper virus AfPV1 in A. flavus, perhaps by reducing the genomic accumulation of the helper virus AfPV1 in A. flavus. We used this model to investigate transcriptional genes regulated by AfPV1 and the small RNA segment in A. flavus, and their role in generating different phenotypes. We found that the pathways of the genes regulated by AfPV1 in its host were similar to those of retroviral viruses. Therefore, some pathways may be of benefit to non-retroviral viral integration or endogenization into the genomes of its host. Moreover, some potential antiviral substances were also found in A. flavus using this system.

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.

Identification of an RNA Silencing Suppressor Encoded by a Symptomless Fungal Hypovirus, Cryphonectria Hypovirus 4

Simple Summary

Host antiviral defense/viral counter-defense is an interesting topic in modern virology. RNA silencing is the primary antiviral mechanism in insects, plants, and fungi, while viruses encode and utilize RNA silencing suppressors against the host defense. Hypoviruses are positive-sense single-stranded RNA viruses with phylogenetic affinity to the picorna-like supergroup, including animal poliovirus and plant potyvirus. The prototype hypovirus Cryphonectria hypovirus 1, CHV1, is one of the best-studied fungal viruses. It is known to induce hypovirulence in the chestnut blight fungus, Cryphonectria parasitica, and encode an RNA silencing suppressor. CHV4 is another hypovirus asymptomatically that infects the same host fungus. This study shows that the N-terminal portion of the CHV4 polyprotein, termed p24, is a protease that autocatalytically cleaves itself from the rest of the viral polyprotein, and functions as an antiviral RNA silencing suppressor.

Abstract

Previously, we have reported the ability of a symptomless hypovirus Cryphonectria hypovirus 4 (CHV4) of the chestnut blight fungus to facilitate stable infection by a co-infecting mycoreovirus 2 (MyRV2)—likely through the inhibitory effect of CHV4 on RNA silencing (Aulia et al., Virology, 2019). In this study, the N-terminal portion of the CHV4 polyprotein, termed p24, is identified as an autocatalytic protease capable of suppressing host antiviral RNA silencing. Using a bacterial expression system, CHV4 p24 is shown to cleave autocatalytically at the di-glycine peptide (Gly214-Gly215) of the polyprotein through its protease activity. Transgenic expression of CHV4 p24 in Cryphonectria parasitica suppresses the induction of one of the key genes of the antiviral RNA silencing, dicer-like 2, and stabilizes the infection of RNA silencing-susceptible virus MyRV2. This study shows functional similarity between CHV4 p24 and its homolog p29, encoded by the symptomatic prototype hypovirus CHV1.

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.

A Novel Mycovirus Evokes Transcriptional Rewiring in the Fungus Malassezia and Stimulates Beta Interferon Production in Macrophages

Malassezia species represent the most common fungal inhabitant of the mammalian skin microbiome and are natural skin commensal flora. However, these fungi are also associated with a variety of clinical skin disorders. Recent studies have reported associations of Malassezia with Crohn’s disease and pancreatic cancer, further implicating this fungal genus in inflammatory and neoplastic disease states. Because M. sympodialis has lost genes involved in RNA interference (RNAi), we hypothesized Malassezia could harbor dsRNA mycoviruses. Indeed, we identified a novel mycovirus of the totivirus family in several Malassezia species and characterized the MsMV1 mycovirus of M. sympodialis. We found conditions that lead to curing of the virus and analyzed isogenic virus-infected/virus-cured strains to determine MsMV1 genetic and pathogenic impacts. MsMV1 induces a strong overexpression of transcription factors and ribosomal genes, while downregulating cellular metabolism. Moreover, MsMV1 induced a significantly higher level of beta interferon expression in cultured macrophages. This study sheds light on the mechanisms of pathogenicity of Malassezia, focusing on a previously unidentified novel mycovirus.

Mycoviruses infect fungi, and while most persist asymptomatically, there are examples of mycoviruses having both beneficial and detrimental effects on their host. Virus-infected Saccharomyces and Ustilago strains exhibit a killer phenotype conferring a growth advantage over uninfected strains and other competing yeast species, whereas hypovirus-infected Cryphonectria parasitica displays defects in growth, sporulation, and virulence. In this study, we identify a double-stranded RNA (dsRNA) mycovirus in five Malassezia species. Sequence analysis reveals it to be a totivirus with two dsRNA segments: a larger 4.5-kb segment with genes encoding components for viral replication and maintenance, and a smaller 1.4-kb segment encoding a novel protein. Furthermore, transcriptome sequencing (RNA-seq) of virus-infected versus virus-cured Malassezia sympodialis revealed an upregulation of dozens of ribosomal components in the cell, suggesting the virus modifies the transcriptional and translational landscapes of the cell. Given that Malassezia is the most abundant fungus on human skin, we assessed the impact of the mycovirus in a murine epicutaneous infection model. Although infection with virus-infected strains was not associated with an increased inflammatory response, we did observe enhanced skin colonization in one of two virus-infected M. sympodialis strains. Noteworthy, beta interferon expression was significantly upregulated in bone marrow-derived macrophages when challenged with virus-infected, compared to virus-cured, M. sympodialis, suggesting that the presence of the virus can induce an immunological response. Although many recent studies have illuminated how widespread mycoviruses are, there are relatively few in-depth studies about their impact on disease caused by the host fungus. We describe here a novel mycovirus in Malassezia and its possible implications in pathogenicity.

A Novel Virus Alters Gene Expression and Vacuolar Morphology in Malassezia Cells and Induces a TLR3-Mediated Inflammatory Immune Response

Malassezia is the most dominant fungal genus on the human skin surface and is associated with various skin diseases including dandruff and seborrheic dermatitis. Among Malassezia species, Malassezia restricta is the most widely observed species on the human skin. In the current study, we identified a novel dsRNA virus, named MrV40, in M. restricta and characterized the sequence and structure of the viral genome along with an independent satellite dsRNA viral segment. Moreover, expression of genes involved in ribosomal synthesis and programmed cell death was altered, indicating that virus infection affected the physiology of the fungal host cells. Our data also showed that the viral nucleic acid from MrV40 induces a TLR3-mediated inflammatory immune response in bone marrow-derived dendritic cells, indicating that a viral element likely contributes to the pathogenicity of Malassezia. This is the first study to identify and characterize a novel mycovirus in Malassezia.

Most fungal viruses have been identified in plant pathogens, whereas the presence of viral particles in human-pathogenic fungi is less well studied. In the present study, we observed extrachromosomal double-stranded RNA (dsRNA) segments in various clinical isolates of Malassezia species. Malassezia is the most dominant fungal genus on the human skin surface, and species in this group are considered etiological factors of various skin diseases including dandruff, seborrheic dermatitis, and atopic dermatitis. We identified novel dsRNA segments, and our sequencing results revealed that the virus, named MrV40, belongs to the Totiviridae family and contains an additional satellite dsRNA segment encoding a novel protein. The transcriptome of virus-infected Malassezia restricta cells was compared to that of virus-cured cells, and the results showed that transcripts involved in ribosomal biosynthesis were downregulated and those involved in energy production and programmed cell death were upregulated. Moreover, transmission electron microscopy revealed significantly larger vacuoles in virus-infected M. restricta cells, indicating that MrV40 infection dramatically altered M. restricta physiology. Our analysis also revealed that viral nucleic acid from MrV40 induced a TLR3 (Toll-like receptor 3)-mediated inflammatory immune response in bone marrow-derived dendritic cells, suggesting that a viral element contributes to the pathogenicity of Malassezia.

Analysis of an Intrinsic Mycovirus Associated With Reduced Virulence of the Human Pathogenic Fungus Aspergillus fumigatus

Aspergillus fumigatus is an airborne fungal pathogen that causes severe infections with invasive growth in immunocompromised patients. Several mycoviruses have recently been isolated from A. fumigatus strains, but there are presently no reports of mycoviral-mediated reduction or elimination of fungal pathogenicity in vertebrate models. Here, we report the biological features of a novel mycovirus, A. fumigatus chrysovirus 41362 (AfuCV41362), isolated from the hypovirulent A. fumigatus strain IFM 41362. The AfuCV41362 genome is comprised of four dsRNAs, each of which contains a single ORF (ORF1-4). ORF1 encodes a protein with sequence similarity to RNA-dependent RNA polymerases of viruses in the family Chrysoviridae, while ORF3 encodes a putative capsid protein. Viral RNAs are expressed primarily during the germination stage, and RNA-seq analysis of virus-infected A. fumigatus at the germination stage suggested that the virus suppressed expression of several pathogenicity-associated host genes, including hypoxia adaptation and nitric oxide detoxification genes. In vitro functional analysis revealed that the virus-infected strain had reduced tolerance to environmental stressors. Virus-infected A. fumigatus strain IFM 41362 had reduced virulence in vivo compared to the virus-free strain in a mouse infection model. Furthermore, introduction of the mycovirus to a natively virus-free KU A. fumigatus strain induced virus-infected phenotypes. To identify mycovirus genes responsible for the reduced virulence of A. fumigatus, each viral ORF was ectopically expressed in the virus-free KU strain. Ectopic expression of the individual ORFs only nominally reduced virulence of the host fungus in a mouse infection model. However, we found that ORF3 and ORF4 reduced tolerance to environmental stresses in in vitro analysis. Based on these results, we suggest that the AfuCV41362 mycovirus ORF3 and ORF4 reduce fungal virulence by suppressing stress tolerance together with other viral genes, rather than alone.

Exploration of the interactions between mycoviruses and Fusarium graminearum

In this review, we discuss recent studies of the interaction between Fusarium graminearum viruses (FgVs) and the fungal host, Fusarium graminearum. Comprehensive transcriptome and proteome analyses have shown changes in the expression of host genes in response to infection by diverse FgVs. Using omics data and reverse genetics, researchers have determined the effects of some fungal host proteins (including FgHex1, FgHal2, FgSwi6, and vr1) on virus accumulation, virus transmission, and host symptom development. Recent reports have revealed the functions of the RNAi component in F. graminearum and the functional redundancy of FgDICERs and FgAGOs in the antiviral defense response against different FgV infections. Studies have also documented a unique mechanism used by FgV1 to overcome the antiviral response of the fungal host.

Defective RNA of a Novel Mycovirus with High Transmissibility Detrimental to Biocontrol Properties of Trichoderma spp

Trichoderma species are a group of fungi which is widely distributed in major terrestrial ecosystems; they are also commonly used as biocontrol agents for many plant diseases. A virus, namely Trichoderma harzianum hypovirus 1 (ThHV1), was identified in T. harzianum isolate T-70, and also infected isolate T-70D, together with its defective RNA (ThHV1-S). The ThHV1 genome possessed two Open Reading Frames (ORFs), namely ORF1 and ORF2. The start codon of ORF2 overlapped with the stop codon of ORF1 in a 43 nt long region. The polypeptide encoded by ORF2 of ThHV1 shared sequence similarities with those of betahypoviruses, indicating that ThHV1 is a novel member of Hypoviridea. Isolate T-70D, carrying both ThHV1 and ThHV1-S, showed abnormal biological properties, notably a decreased mycoparasitism ability when compared with isolate T-70. Both ThHV1 and ThHV1-S could be vertically transmitted to conidia and horizontally transmitted to T. harzianum isolate T-68 and T. koningiopsis T-51. The derivative strains carrying both ThHV1 and ThHV1-S showed decreased mycoparasitism ability, whereas strains carrying ThHV1 alone were normal, indicating that ThHV1-S is closely associated with the decreased mycoparasitism ability of T. harzianum isolate T-70D. ThHV1 was widely detected in isolates of T. harzianum, T. koningiopsis and T. atroviride originating from soil of China. Therefore, viruses in fungal biocontrol agents may also be a factor associated with the stability of their application.

Comprehensive analysis of full genome sequence and Bd-milRNA/target mRNAs to discover the mechanism of hypovirulence in Botryosphaeria dothidea strains on pear infection with BdCV1 and BdPV1

Pear ring rot disease, mainly caused by Botryosphaeria dothidea, is widespread in most pear and apple-growing regions. Mycoviruses are used for biocontrol, especially in fruit tree disease. BdCV1 (Botryosphaeria dothidea chrysovirus 1) and BdPV1 (Botryosphaeria dothidea partitivirus 1) influence the biological characteristics of B. dothidea strains. BdCV1 is a potential candidate for the control of fungal disease. Therefore, it is vital to explore interactions between B. dothidea and mycovirus to clarify the pathogenic mechanisms of B. dothidea and hypovirulence of B. dothidea in pear. A high-quality full-length genome sequence of the B. dothidea LW-Hubei isolate was obtained using Single Molecule Real-Time sequencing. It has high repeat sequence with 9.3% and DNA methylation existence in the genome. The 46.34 Mb genomes contained 14,091 predicted genes, which of 13,135 were annotated. B. dothidea was predicted to express 3833 secreted proteins. In bioinformatics analysis, 351 CAZy members, 552 transporters, 128 kinases, and 1096 proteins associated with plant-host interaction (PHI) were identified. RNA-silencing components including two endoribonuclease Dicer, four argonaute (Ago) and three RNA-dependent RNA polymerase (RdRp) molecules were identified and expressed in response to mycovirus infection. Horizontal transfer of the LW-C and LW-P strains indicated that BdCV1 induced host gene silencing in LW-C to suppress BdPV1 transmission. To investigate the role of RNA-silencing in B. dothidea defense, we constructed four small RNA libraries and sequenced B. dothidea micro-like RNAs (Bd-milRNAs) produced in response to BdCV1 and BdPV1 infection. Among these, 167 conserved and 68 candidate novel Bd-milRNAs were identified, of which 161 conserved and 20 novel Bd-milRNA were differentially expressed. WEGO analysis revealed involvement of the differentially expressed Bd-milRNA-targeted genes in metabolic process, catalytic activity, cell process and response to stress or stimulus. BdCV1 had a greater effect on the phenotype, virulence, conidiomata, vertical and horizontal transmission ability, and mycelia cellular structure biological characteristics of B. dothidea strains than BdPV1 and virus-free strains. The results obtained in this study indicate that mycovirus regulates biological processes in B. dothidea through the combined interaction of antiviral defense mediated by RNA-silencing and milRNA-mediated regulation of target gene mRNA expression.

A symptomless hypovirus, CHV4, facilitates stable infection of the chestnut blight fungus by a coinfecting reovirus likely through suppression of antiviral RNA silencing

Field-collected US strain C18 of Cryphonectria parasitica, the chestnut blight fungus, was earlier reported to be infected by a double-stranded RNA virus, mycoreovirus 2 (MyRV2). Next-generation sequencing has revealed co-infection of C18 by a positive-strand RNA virus, hypovirus 4 (CHV4). The current molecular and genetic analyses showed interesting commensal interactions between the two viruses. CHV4 facilitated the stable infection and enhanced vertical transmission of MyRV2, which was readily lost during subculturing and showed reduced vertical transmission in single infections. Deletion of a key antiviral RNA silencing gene, dcl2, in isolate C18 increased stability of MyRV2 in single infections. The ability of CHV4 to facilitate stable infection with MyRV2 appears to be associated with the inhibitory effect of CHV4 on RNA silencing via compromising the induction of transcriptional up-regulation of dcl2. These results suggest that natural infection of isolate C18 by MyRV2 in the field was facilitated by CHV4 co-infection.

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.

Transcriptional and Small RNA Responses of the White Mold Fungus Sclerotinia sclerotiorum to Infection by a Virulence-Attenuating Hypovirus

Mycoviruses belonging to the family Hypoviridae cause persistent infection of many different host fungi. We previously determined that the white mold fungus, Sclerotiniasclerotiorum, infected with Sclerotinia sclerotiorum hypovirus 2-L (SsHV2-L) exhibits reduced virulence, delayed/reduced sclerotial formation, and enhanced production of aerial mycelia. To gain better insight into the cellular basis for these changes, we characterized changes in mRNA and small RNA (sRNA) accumulation in S.sclerotiorum to infection by SsHV2-L. A total of 958 mRNAs and 835 sRNA-producing loci were altered after infection by SsHV2-L, among which >100 mRNAs were predicted to encode proteins involved in the metabolism and trafficking of carbohydrates and lipids. Both S. sclerotiorum endogenous and virus-derived sRNAs were predominantly 22 nt in length suggesting one dicer-like enzyme cleaves both. Novel classes of endogenous small RNAs were predicted, including phasiRNAs and tRNA-derived small RNAs. Moreover, S. sclerotiorum phasiRNAs, which were derived from noncoding RNAs and have the potential to regulate mRNA abundance in trans, showed differential accumulation due to virus infection. tRNA fragments did not accumulate differentially after hypovirus infection. Hence, in-depth analysis showed that infection of S. sclerotiorum by a hypovirulence-inducing hypovirus produced selective, large-scale reprogramming of mRNA and sRNA production.

De novo transcriptomic assembly and mRNA expression patterns of Botryosphaeria dothidea infection with mycoviruses chrysovirus 1 (BdCV1) and partitivirus 1 (BdPV1)

BACKGROUND

Pear ring rot, caused by Botryosphaeria species, is responsible for substantial economic losses by causing severe recession of pear tree growth in China. Mycovirus-mediated hypovirulence in plant pathogenic fungi is a crucial biological method to control fungal diseases.

METHODS

We conducted a large-scale and comprehensive transcriptome analysis to identify mRNA in B. dothidea in response to mycovirus. De novo sequencing technology from four constructed libraries of LW-C (Botryosphaeria dothidea chrysovirus 1, BdCV1), LW-P (Botryosphaeria dothidea partitivirus 1, BdPV1), LW-CP (LW-1 strain infection with BdCV1 and BdPV1), and Mock (free virus) was used to investigate and compare gene expression changes in B.dothidea strains infected with mycovirus.

RESULTS

In total, 30,058 Unigenes with an average length of 2128 bp were obtained from 4 libraries of B. dothidea strains. These were annotated to specify their classified function. We demonstrate that mRNAs of B. dothidea strains in response to mycovirus are differentially expressed. In total, 5598 genes were up-regulated and 3298 were down-regulated in the LW-CP group, 4468 were up-regulated and 4291 down-regulated in the LW-C group, and 2590 were up-regulated and 2325 down-regulated in the LW-P group. RT-qPCR was used to validate the expression of 9 selected genes. The B. dothidea transcriptome was more affected by BdCV1 infection than BdPV1. We conducted GO enrichment analysis to characterize gene functions regulated by B. dothidea with mycovirus infection. These involved metabolic process, cellular process, catalytic activity, transporter activity, signaling, and other biological pathways. KEGG function analysis demonstrated that the enriched differentially expressed genes are involved in metabolism, transcription, signal transduction, and ABC transport. mRNA is therefore involved in the interaction between fungi and mycovirus. In addition, changes in differential accumulation levels of cp and RdRp of BdCV1 and BdPV1 in B. dothidea strains were evaluated, revealing that the accumulation of BdCV1 and BdPV1 is related to the phenotype and virulence of B. dothidea strain LW-1.

CONCLUSIONS

The identification and analysis of mRNAs from B. dothidea was first reported at the transcriptome level. Our analysis provides further insight into the interaction of B. dothidea strains infection with chrysovirus 1 (BdCV1) and partitivirus 1 (BdPV1) at the transcriptome level.

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.

Differential Contribution of RNA Interference Components in Response to Distinct Fusarium graminearum Virus Infections

The mechanisms of RNA interference (RNAi) as a defense response against viruses remain unclear in many plant-pathogenic fungi. In this study, we used reverse genetics and virus-derived small RNA profiling to investigate the contributions of RNAi components to the antiviral response against Fusarium graminearum viruses 1 to 3 (FgV1, -2, and -3). Real-time reverse transcription-quantitative PCR (qRT-PCR) indicated that infection of Fusarium graminearum by FgV1, -2, or -3 differentially induces the gene expression of RNAi components in F. graminearum Transcripts of the DICER-2 and AGO-1 genes of F. graminearum (FgDICER-2 and FgAGO-1) accumulated at lower levels following FgV1 infection than following FgV2 or FgV3 infection. We constructed gene disruption and overexpression mutants for each of the Argonaute and dicer genes and for two RNA-dependent RNA polymerase (RdRP) genes and generated virus-infected strains of each mutant. Interestingly, mycelial growth was significantly faster for the FgV1-infected FgAGO-1 overexpression mutant than for the FgV1-infected wild type, while neither FgV2 nor FgV3 infection altered the colony morphology of the gene deletion and overexpression mutants. FgV1 RNA accumulation was significantly decreased in the FgAGO-1 overexpression mutant. Furthermore, the levels of induction of FgAGO-1, FgDICER-2, and some of the FgRdRP genes caused by FgV2 and FgV3 infection were similar to those caused by hairpin RNA-induced gene silencing. Using small RNA sequencing analysis, we documented different patterns of virus-derived small interfering RNA (vsiRNA) production in strains infected with FgV1, -2, and -3. Our results suggest that the Argonaute protein encoded by FgAGO-1 is required for RNAi in F. graminearum, that FgAGO-1 induction differs in response to FgV1, -2, and -3, and that FgAGO-1 might contribute to the accumulation of vsiRNAs in FgV1-infected F. graminearumIMPORTANCE To increase our understanding of how RNAi components in Fusarium graminearum react to mycovirus infections, we characterized the role(s) of RNAi components involved in the antiviral defense response against Fusarium graminearum viruses (FgVs). We observed differences in the levels of induction of RNA silencing-related genes, including FgDICER-2 and FgAGO-1, in response to infection by three different FgVs. FgAGO-1 can efficiently induce a robust RNAi response against FgV1 infection, but FgDICER genes might be relatively redundant to FgAGO-1 with respect to antiviral defense. However, the contribution of this gene in the response to the other FgV infections might be small. Compared to previous studies of Cryphonectria parasitica, which showed dicer-like protein 2 and Argonaute-like protein 2 to be important in antiviral RNA silencing, our results showed that F. graminearum developed a more complex and robust RNA silencing system against mycoviruses and that FgDICER-1 and FgDICER-2 and FgAGO-1 and FgAGO-2 had redundant roles in antiviral RNA silencing.

Mycoviruses as Triggers and Targets of RNA Silencing in White Mold Fungus Sclerotinia sclerotiorum

This study aimed to demonstrate the existence of antiviral RNA silencing mechanisms in Sclerotinia sclerotiorum by infecting wild-type and RNA-silencing-deficient strains of the fungus with an RNA virus and a DNA virus. Key silencing-related genes were disrupted to dissect the RNA silencing pathway. Specifically, dicer genes (dcl-1, dcl-2, and both dcl-1/dcl-2) were displaced by selective marker(s). Disruption mutants were then compared for changes in phenotype, virulence, and susceptibility to virus infections. Wild-type and mutant strains were transfected with a single-stranded RNA virus, SsHV2-L, and copies of a single-stranded DNA mycovirus, SsHADV-1, as a synthetic virus constructed in this study. Disruption of dcl-1 or dcl-2 resulted in no changes in phenotype compared to wild-type S. sclerotiorum; however, the double dicer mutant strain exhibited significantly slower growth. Furthermore, the Δdcl-1/dcl-2 double mutant, which was slow growing without virus infection, exhibited much more severe debilitation following virus infections including phenotypic changes such as slower growth, reduced pigmentation, and delayed sclerotial formation. These phenotypic changes were absent in the single mutants, Δdcl-1 and Δdcl-2. Complementation of a single dicer in the double disruption mutant reversed viral susceptibility to the wild-type state. Virus-derived small RNAs were accumulated from virus-infected wild-type strains with strand bias towards the negative sense. The findings of these studies indicate that S. sclerotiorum has robust RNA silencing mechanisms that process both DNA and RNA mycoviruses and that, when both dicers are silenced, invasive nucleic acids can greatly debilitate the virulence of this fungus.

Trans-Kingdom RNA Silencing in Plant-Fungal Pathogen Interactions

Fungal pathogens represent a major group of plant invaders that are the causative agents of many notorious plant diseases. Large quantities of RNAs, especially small RNAs involved in gene silencing, have been found to transmit bidirectionally between fungal pathogens and their hosts. Although host-induced gene silencing (HIGS) technology has been developed and applied to protect crops from fungal infections, the mechanisms of RNA transmission, especially small RNAs regulating trans-kingdom RNA silencing in plant immunity, are largely unknown. In this review, we summarize and discuss recent important findings regarding trans-kingdom sRNAs and RNA silencing in plant-fungal pathogen interactions compared with the well-known RNAi mechanisms in plants and fungi. We focus on the interactions between plant and fungal pathogens with broad hosts, represented by the vascular pathogen Verticillium dahliae and non-vascular pathogen Botrytis cinerea, and discuss the known instances of natural RNAi transmission between fungal pathogens and host plants. Given that HIGS has been developed and recently applied in controlling Verticillium wilt diseases, we propose an ideal research system exploiting plant vasculature-Verticillium interaction to further study trans-kingdom RNA silencing.

Development of an RNA interference (RNAi) gene knockdown protocol in the anaerobic gut fungus Pecoramyces ruminantium strain C1A

Members of the anaerobic gut fungi (AGF) reside in rumen, hindgut, and feces of ruminant and non-ruminant herbivorous mammals and reptilian herbivores. No protocols for gene insertion, deletion, silencing, or mutation are currently available for the AGF, rendering gene-targeted molecular biological manipulations unfeasible. Here, we developed and optimized an RNA interference (RNAi)-based protocol for targeted gene silencing in the anaerobic gut fungus Pecoramyces ruminantium strain C1A. Analysis of the C1A genome identified genes encoding enzymes required for RNA silencing in fungi (Dicer, Argonaute, Neurospora crassa QDE-3 homolog DNA helicase, Argonaute-interacting protein, and Neurospora crassa QIP homolog exonuclease); and the competency of C1A germinating spores for RNA uptake was confirmed using fluorescently labeled small interfering RNAs (siRNA). Addition of chemically-synthesized siRNAs targeting D-lactate dehydrogenase (ldhD) gene to C1A germinating spores resulted in marked target gene silencing; as evident by significantly lower ldhD transcriptional levels, a marked reduction in the D-LDH specific enzymatic activity in intracellular protein extracts, and a reduction in D-lactate levels accumulating in the culture supernatant. Comparative transcriptomic analysis of untreated versus siRNA-treated cultures identified a few off-target siRNA-mediated gene silencing effects. As well, significant differential up-regulation of the gene encoding NAD-dependent 2-hydroxyacid dehydrogenase (Pfam00389) in siRNA-treated C1A cultures was observed, which could possibly compensate for loss of D-LDH as an electron sink mechanism in C1A. The results demonstrate the feasibility of RNAi in anaerobic fungi, and opens the door for gene silencing-based studies in this fungal clade.

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.

Deep sequencing of mycovirus-derived small RNAs from Botrytis species

RNA silencing is an ancient regulatory mechanism operating in all eukaryotic cells. In fungi, it was first discovered in Neurospora crassa, although its potential as a defence mechanism against mycoviruses was first reported in Cryphonectria parasitica and, later, in several fungal species. There is little evidence of the antiviral potential of RNA silencing in the phytopathogenic species of the fungal genus Botrytis. Moreover, little is known about the RNA silencing components in these fungi, although the analysis of public genome databases identified two Dicer-like genes in B. cinerea, as in most of the ascomycetes sequenced to date. In this work, we used deep sequencing to study the virus-derived small RNA (vsiRNA) populations from different mycoviruses infecting field isolates of Botrytis spp. The mycoviruses under study belong to different genera and species, and have different types of genome [double-stranded RNA (dsRNA), (+)single-stranded RNA (ssRNA) and (-)ssRNA]. In general, vsiRNAs derived from mycoviruses are mostly of 21, 20 and 22 nucleotides in length, possess sense or antisense orientation, either in a similar ratio or with a predominance of sense polarity depending on the virus species, have predominantly U at their 5' end, and are unevenly distributed along the viral genome, showing conspicuous hotspots of vsiRNA accumulation. These characteristics reveal striking similarities with vsiRNAs produced by plant viruses, suggesting similar pathways of viral targeting in plants and fungi. We have shown that the fungal RNA silencing machinery acts against the mycoviruses used in this work in a similar manner independent of their viral or fungal origin.

The RNAi Universe in Fungi: A Varied Landscape of Small RNAs and Biological Functions

RNA interference (RNAi) is a conserved eukaryotic mechanism that uses small RNA molecules to suppress gene expression through sequence-specific messenger RNA degradation, translational repression, or transcriptional inhibition. In filamentous fungi, the protective function of RNAi in the maintenance of genome integrity is well known. However, knowledge of the regulatory role of RNAi in fungi has had to wait until the recent identification of different endogenous small RNA classes, which are generated by distinct RNAi pathways. In addition, RNAi research on new fungal models has uncovered the role of small RNAs and RNAi pathways in the regulation of diverse biological functions. In this review, we give an up-to-date overview of the different classes of small RNAs and RNAi pathways in fungi and their roles in the defense of genome integrity and regulation of fungal physiology and development, as well as in the interaction of fungi with biotic and abiotic environments.

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.

Genome-wide exonic small interference RNA-mediated gene silencing regulates sexual reproduction in the homothallic fungus Fusarium graminearum

Various ascomycete fungi possess sex-specific molecular mechanisms, such as repeat-induced point mutations, meiotic silencing by unpaired DNA, and unusual adenosine-to-inosine RNA editing, for genome defense or gene regulation. Using a combined analysis of functional genetics and deep sequencing of small noncoding RNA (sRNA), mRNA, and the degradome, we found that the sex-specifically induced exonic small interference RNA (ex-siRNA)-mediated RNA interference (RNAi) mechanism has an important role in fine-tuning the transcriptome during ascospore formation in the head blight fungus Fusarium graminearum. Approximately one-third of the total sRNAs were produced from the gene region, and sRNAs with an antisense direction or 5'-U were involved in post-transcriptional gene regulation by reducing the stability of the corresponding gene transcripts. Although both Dicers and Argonautes partially share their functions, the sex-specific RNAi pathway is primarily mediated by FgDicer1 and FgAgo2, while the constitutively expressed RNAi components FgDicer2 and FgAgo1 are responsible for hairpin-induced RNAi. Based on our results, we concluded that F. graminearum primarily utilizes ex-siRNA-mediated RNAi for ascosporogenesis but not for genome defenses and other developmental stages. Each fungal species appears to have evolved RNAi-based gene regulation for specific developmental stages or stress responses. This study provides new insights into the regulatory role of sRNAs in fungi and other lower eukaryotes.

Small RNA-Based Antiviral Defense in the Phytopathogenic Fungus Colletotrichum higginsianum

Even though the fungal kingdom contains more than 3 million species, little is known about the biological roles of RNA silencing in fungi. The Colletotrichum genus comprises fungal species that are pathogenic for a wide range of crop species worldwide. To investigate the role of RNA silencing in the ascomycete fungus Colletotrichum higginsianum, knock-out mutants affecting genes for three RNA-dependent RNA polymerase (RDR), two Dicer-like (DCL), and two Argonaute (AGO) proteins were generated by targeted gene replacement. No effects were observed on vegetative growth for any mutant strain when grown on complex or minimal media. However, Δdcl1, Δdcl1Δdcl2 double mutant, and Δago1 strains showed severe defects in conidiation and conidia morphology. Total RNA transcripts and small RNA populations were analyzed in parental and mutant strains. The greatest effects on both RNA populations was observed in the Δdcl1, Δdcl1Δdcl2, and Δago1 strains, in which a previously uncharacterized dsRNA mycovirus [termed Colletotrichum higginsianum non-segmented dsRNA virus 1 (ChNRV1)] was derepressed. Phylogenetic analyses clearly showed a close relationship between ChNRV1 and members of the segmented Partitiviridae family, despite the non-segmented nature of the genome. Immunoprecipitation of small RNAs associated with AGO1 showed abundant loading of 5’U-containing viral siRNA. C. higginsianum parental and Δdcl1 mutant strains cured of ChNRV1 revealed that the conidiation and spore morphology defects were primarily caused by ChNRV1. Based on these results, RNA silencing involving ChDCL1 and ChAGO1 in C. higginsianum is proposed to function as an antiviral mechanism.

Colletotrichum sp. comprises a diverse group of fungal pathogens that attack over 3000 plant species worldwide. Understanding the underlying mechanisms that govern fungal development and pathogenicity may enable more effective and sustainable approaches to crop disease management and control. In most organisms, RNA silencing is an important mechanism to control endogenous and exogenous RNA. RNA silencing utilizes small regulatory molecules (small RNAs) produced by proteins called Dicer (DCL), and exercise their function though effector proteins named Argonaute (AGO). Here, we investigated the role of RNA silencing machinery in the fungus Colletotrichum higginsianum, by generating deletions in genes encoding RNA silencing components. Severe defects were observed in both conidiation and conidia morphology in the Δdcl1, Δdcl1Δdcl2, and Δago1 strains. Analysis of transcripts and small RNAs revealed an uncharacterized dsRNA virus persistently infecting C. higginsianum. The virus was shown (1) to be de-repressed in the Δdcl1, Δdcl1Δdcl2 and Δago1 strains, and (2) to cause the conidiation and spore mutant phenotypes. Our results indicate that C. higginsianum employs RNA silencing as an antiviral mechanism to suppress viruses and their debilitating effects.

Differential Inductions of RNA Silencing among Encapsidated Double-Stranded RNA Mycoviruses in the White Root Rot Fungus Rosellinia necatrix

RNA silencing acts as a defense mechanism against virus infection in a wide variety of organisms. Here, we investigated inductions of RNA silencing against encapsidated double-stranded RNA (dsRNA) fungal viruses (mycoviruses), including a partitivirus (RnPV1), a quadrivirus (RnQV1), a victorivirus (RnVV1), a mycoreovirus (RnMyRV3), and a megabirnavirus (RnMBV1) in the phytopathogenic fungus Rosellinia necatrix Expression profiling of RNA silencing-related genes revealed that a dicer-like gene, an Argonaute-like gene, and two RNA-dependent RNA polymerase genes were upregulated by RnMyRV3 or RnMBV1 infection but not by other virus infections or by constitutive expression of dsRNA in R. necatrix Massive analysis of viral small RNAs (vsRNAs) from the five mycoviruses showed that 19- to 22-nucleotide (nt) vsRNAs were predominant; however, their ability to form duplexes with 3' overhangs and the 5' nucleotide preferences of vsRNAs differed among the five mycoviruses. The abundances of 19- to 22-nt vsRNAs from RnPV1, RnQV1, RnVV1, RnMyRV3, and RnMBV1 were 6.8%, 1.2%, 0.3%, 13.0%, and 24.9%, respectively. Importantly, the vsRNA abundances and accumulation levels of viral RNA were not always correlated, and the origins of the vsRNAs were distinguishable among the five mycoviruses. These data corroborated diverse interactions between encapsidated dsRNA mycoviruses and RNA silencing. Moreover, a green fluorescent protein (GFP)-based sensor assay in R. necatrix revealed that RnMBV1 infection induced silencing of the target sensor gene (GFP gene and the partial RnMBV1 sequence), suggesting that vsRNAs from RnMBV1 activated the RNA-induced silencing complex. Overall, this study provides insights into RNA silencing against encapsidated dsRNA mycoviruses.

IMPORTANCE

Encapsidated dsRNA fungal viruses (mycoviruses) are believed to replicate inside their virions; therefore, there is a question of whether they induce RNA silencing. Here, we investigated inductions of RNA silencing against encapsidated dsRNA mycoviruses (a partitivirus, a quadrivirus, a victorivirus, a mycoreovirus, and a megabirnavirus) in Rosellinia necatrix We revealed upregulation of RNA silencing-related genes in R. necatrix infected with a mycoreovirus or a megabirnavirus but not with other viruses, which was consistent with the relatively high abundances of vsRNAs from the two mycoviruses. We also showed common and different molecular features and origins of the vsRNAs from the five mycoviruses. Furthermore, we demonstrated the activation of RNA-induced silencing complex by mycoviruses in R. necatrix Taken together, our data provide insights into an RNA silencing pathway against encapsidated dsRNA mycoviruses which is differentially induced among encapsidated dsRNA mycoviruses; that is, diverse replication strategies exist among encapsidated dsRNA mycoviruses.

Generation of a high resolution map of sRNAs from Fusarium graminearum and analysis of responses to viral infection

Previously, we characterized F. graminearum hypovirus 1 (FgHV1) and F. graminearum hypovirus 2 (FgHV2), which are the only two hypoviruses in F. graminearum that are closely related to Cryphonectria hypovirus 1 (CHV1) and Cryphonectria hypovirus 2 (CHV2) in the Hypoviridae family. In this study, we preliminarily elucidated the RNA silencing mechanism of the F. graminearum/hypovirus system from a small RNA (sRNA) perspective by using HiSeq deep sequencing. The length distributions of F. graminearum sRNA were altered by hypoviral infection. Potential microRNA-like (milRNA) candidates were differentially expressed between the hypovirus-free and hypovirus-infected library types. Extensive virus-derived small interfering RNAs (vsiRNAs) were also principally defined. The 1,831,081 and 3,254,758 total reads generated from the FgHV1 and FgHV2 genomes in F. graminearum yielded the first high-resolution sRNA maps of fungal viruses. In addition, extensive bioinformatics searches identified a large number of transcripts that are potentially targeted by vsiRNAs, several of which were effectively down-regulated. In particular, the RNA silencing-related genes FgDicer1 and FgRdRp5 were predicted targets of FgHV1- and FgHV2-derived siRNAs, possibly revealing a novel anti-RNA silencing strategy employed by mycoviruses.

Virus-derived small RNAs: molecular footprints of host-pathogen interactions

Viruses are obligatory intracellular parasites that require the host machinery to replicate. During their replication cycle, viral RNA intermediates can be recognized and degraded by different antiviral mechanisms that include RNA decay, RNA interference, and RNase L pathways. As a consequence of viral RNA degradation, infected cells can accumulate virus‐derived small RNAs at high levels compared to cellular molecules. These small RNAs are imprinted with molecular characteristics that reflect their origin. First, small RNAs can be used to reconstruct viral sequences and identify the virus from which they originated. Second, other molecular features of small RNAs such as size, polarity, and base preferences depend on the type of viral substrate and host mechanism of degradation. Thus, the pattern of small RNAs generated in infected cells can be used as a molecular footprint to identify and characterize viruses independent on sequence homology searches against known references. Hence, sequencing of small RNAs obtained from infected cells enables virus discovery and characterization using both sequence‐dependent strategies and novel pattern‐based approaches. Recent studies are helping unlock the full application of small RNA sequencing for virus discovery and characterization. WIREs RNA 2016, 7:824–837. doi: 10.1002/wrna.1361

This article is categorized under: 1

RNA Processing > Processing of Small RNAs

2

RNA Turnover and Surveillance > Turnover/Surveillance Mechanisms

3

RNA Methods > RNA Analyses In Vitro and In Silico

Multiple approaches for the detection and characterization of viral and plasmid symbionts from a collection of marine fungi

The number of reported mycoviruses is increasing exponentially due to the current ability to detect mycoviruses using next-generation sequencing (NGS) approaches, with a large number of viral genomes built in-silico using data from fungal transcriptome projects. We decided to screen a collection of fungi originating from a specific marine environment (associated with the seagrass Posidonia oceanica) for the presence of mycoviruses: our findings reveal a wealth of diversity among these symbionts and this complexity will require further studies to address their specific role in this ecological niche. In specific, we identified twelve new virus species belonging to nine distinct lineages: they are members of megabirnavirus, totivirus, chrysovirus, partitivirus and five still undefined clades. We showed evidence of an endogenized virus ORF, and evidence of accumulation of dsRNA from metaviridae retroviral elements. We applied different techniques for detecting the presence of mycoviruses including (i) dsRNA extraction and cDNA cloning, (ii) small and total RNA sequencing through NGS techniques, (iii) rolling circle amplification (RCA) and total DNA extraction analyses, (iv) virus purifications and electron microscopy. We tried also to critically evaluate the intrinsic value and limitations of each of these techniques. Based on the samples we could compare directly, RNAseq analysis is superior to sRNA for de novo assembly of mycoviruses. To our knowledge this is the first report on the virome of fungi isolated from marine environment. The GenBank/eMBL/DDBJ accession numbers of the sequences reported in this paper are: KT601099-KT601110; KT601114-KT601120; KT592305; KT950836-KT950841.

Taxonomy, biogeography and importance of Heterobasidion viruses

The genus Heterobasidion consists of several species of necrotrophic and saprotrophic fungi, and includes some of the most detrimental organisms in boreal conifer forests. These fungi host a widespread and diverse mycovirus community composed of more than 16 species of Partitiviridae, a species of Narnaviridae and one taxonomically unassigned virus related to the Curvularia thermal tolerance virus. These viruses are able to cross species borders, co-infect single host strains and cause phenotypic changes in their hosts. The abundance of viruses increases over time in Heterobasidion infection centers, and they are targeted by fungal RNA interference. Long-term field studies are essential for obtaining a comprehensive view of virus effects in the nature.

A novel betapartitivirus RnPV6 from Rosellinia necatrix tolerates host RNA silencing but is interfered by its defective RNAs

The family Partitiviridae comprises of five genera with bi-segmented dsRNA genomes that accommodate members infecting plants, fungi or protists. All partitiviruses with only a few exceptions cause asymptomatic infections. We report the characterization of a novel betapartitivirus termed Rosellinia necatrix partitivirus 6 (RnPV6) from a field isolate of a plant pathogenic fungus, white root rot fungus. RnPV6 has typical partitivirus features: dsRNA1 and dsRNA2 are 2462 and 2499bps in length encoding RNA-dependent RNA polymerase and capsid protein. Purified particles are spherical with a diameter of 30nm. Taking advantage of infectivity as virions, RnPV6 was introduced into a model filamentous fungal host, chestnut blight fungus to investigate virus/host interactions. Unlike other partitiviruses tested previously, RnPV6 induced profound phenotypic alterations with symptoms characterized by a reduced growth rate and enhanced pigmentation and was tolerant to host RNA silencing. In addition, a variety of defective RNAs derived from dsRNA1 appear after virion transfection. These sub-viral RNAs were shown to interfere with RnPV6 replication, at least for that of cognate segment dsRNA1. Presence of these sub-viral elements resulted in reduced symptom expression by RnPV6, suggesting their nature as defective-interfering RNAs. The features of RnPV6 are similar to but distinct from those of a previously reported alphapartitivirus, Rosellinia necatrix partitivirus 2 that is susceptible to RNA silencing.

Molecular characterization of a novel hypovirus from the plant pathogenic fungus Fusarium graminearum

A novel mycovirus, termed Fusarium graminearum Hypovirus 2 (FgHV2/JS16), isolated from a plant pathogenic fungus, Fusarium graminearum strain JS16, was molecularly and biologically characterized. The genome of FgHV2/JS16 is 12,800 nucleotides (nts) long, excluding the poly (A) tail. This genome has only one large putative open reading frame, which encodes a polyprotein containing three normal functional domains, papain-like protease, RNA-dependent RNA polymerase, RNA helicase, and a novel domain with homologous bacterial SMC (structural maintenance of chromosomes) chromosome segregation proteins. A defective RNA segment that is 4553-nts long, excluding the poly (A) tail, was also detected in strain JS16. The polyprotein shared significant aa identities with Cryphonectria hypovirus 1 (CHV1) (16.8%) and CHV2 (16.2%). Phylogenetic analyses based on multiple alignments of the polyprotein clearly divided the members of Hypoviridae into two major groups, suggesting that FgHV2/JS16 was a novel hypovirus of a newly proposed genus-Alphahypovirus-composed of the members of Group 1, including CHV1, CHV2, FgHV1 and Sclerotinia sclerotiorum hypovirus 2. FgHV2/JS16 was shown to be associated with hypovirulence phenotypes according to comparisons of the biological properties shared between FgHV2/JS16-infected and FgHV2/JS16-free isogenic strains. Furthermore, we demonstrated that FgHV2/JS16 infection activated the RNA interference pathway in Fusarium graminearum by relative quantitative real time RT-PCR.

50-plus years of fungal viruses

Mycoviruses are widespread in all major taxa of fungi. They are transmitted intracellularly during cell division, sporogenesis, and/or cell-to-cell fusion (hyphal anastomosis), and thus their life cycles generally lack an extracellular phase. Their natural host ranges are limited to individuals within the same or closely related vegetative compatibility groups, although recent advances have established expanded experimental host ranges for some mycoviruses. Most known mycoviruses have dsRNA genomes packaged in isometric particles, but an increasing number of positive- or negative-strand ssRNA and ssDNA viruses have been isolated and characterized. Although many mycoviruses do not have marked effects on their hosts, those that reduce the virulence of their phytopathogenic fungal hosts are of considerable interest for development of novel biocontrol strategies. Mycoviruses that infect endophytic fungi and those that encode killer toxins are also of special interest. Structural analyses of mycoviruses have promoted better understanding of virus assembly, function, and evolution.

Discovery of replicating circular RNAs by RNA-seq and computational algorithms

Replicating circular RNAs are independent plant pathogens known as viroids, or act to modulate the pathogenesis of plant and animal viruses as their satellite RNAs. The rate of discovery of these subviral pathogens was low over the past 40 years because the classical approaches are technical demanding and time-consuming. We previously described an approach for homology-independent discovery of replicating circular RNAs by analysing the total small RNA populations from samples of diseased tissues with a computational program known as progressive filtering of overlapping small RNAs (PFOR). However, PFOR written in PERL language is extremely slow and is unable to discover those subviral pathogens that do not trigger in vivo accumulation of extensively overlapping small RNAs. Moreover, PFOR is yet to identify a new viroid capable of initiating independent infection. Here we report the development of PFOR2 that adopted parallel programming in the C++ language and was 3 to 8 times faster than PFOR. A new computational program was further developed and incorporated into PFOR2 to allow the identification of circular RNAs by deep sequencing of long RNAs instead of small RNAs. PFOR2 analysis of the small RNA libraries from grapevine and apple plants led to the discovery of Grapevine latent viroid (GLVd) and Apple hammerhead viroid-like RNA (AHVd-like RNA), respectively. GLVd was proposed as a new species in the genus Apscaviroid, because it contained the typical structural elements found in this group of viroids and initiated independent infection in grapevine seedlings. AHVd-like RNA encoded a biologically active hammerhead ribozyme in both polarities, and was not specifically associated with any of the viruses found in apple plants. We propose that these computational algorithms have the potential to discover novel circular RNAs in plants, invertebrates and vertebrates regardless of whether they replicate and/or induce the in vivo accumulation of small RNAs.

Diagnosis and discovery of fungal viruses using deep sequencing of small RNAs

Analysis of virus-derived small RNAs with high-throughput sequencing has been successful for detecting novel viruses in plants and invertebrates. However, the applicability of this method has not been demonstrated in fungi, although fungi were among the first organisms reported to utilize RNA silencing. Here, we used virus-infected isolates of the fungal species complex Heterobasidion annosum sensu lato as a model system to test whether mycovirus genome segments can be detected with small RNA deep sequencing. Species of the genus Heterobasidion are some of the most devastating forest pathogens in boreal forests. These fungi cause wood decay and are commonly infected with species of the family Partitiviridae and the unassigned virus species Heterobasidion RNA virus 6. Small RNA deep sequencing allowed the simultaneous detection of all eight double-stranded RNA virus strains known to be present in the tested samples and one putative mitovirus species (family Narnaviridae) with a single-stranded RNA genome, designated here as Heterobasidion mitovirus 1. Prior to this study, no members of the family Narnaviridae had been described as infecting species of Heterobasidion. Quantification of viral double- and single-stranded RNA with quantitative PCR indicated that co-infecting viral species and viruses with segmented genomes can be detected with small RNA deep sequencing despite vast differences in the amount of RNA. This is the first study demonstrating the usefulness of this method for detecting fungal viruses. Moreover, the results suggest that viral genomes are processed into small RNAs by different species of Heterobasidion.

Characterizing the roles of Cryphonectria parasitica RNA-dependent RNA polymerase-like genes in antiviral defense, viral recombination and transposon transcript accumulation

An inducible RNA-silencing pathway, involving a single Dicer protein, DCL2, and a single Argonaute protein, AGL2, was recently shown to serve as an effective antiviral defense response in the chestnut blight fungus Cryphonectria parasitica. Eukaryotic RNA-dependent RNA polymerases (RdRPs) are frequently involved in transcriptional and posttranscriptional gene silencing and antiviral defense. We report here the identification and characterization of four RdRP genes (rdr1–4) in the C. parasitica genome. Sequence relationships with other eukaryotic RdRPs indicated that RDR1 and RDR2 were closely related to QDE-1, an RdRP involved in RNA silencing (“quelling”) in Neurospora crassa, whereas RDR3 was more closely related to the meiotic silencing gene SAD-1 in N. crassa. The RdRP domain of RDR4, related to N. crassa RRP-3 of unknown function, was truncated and showed evidence of alternative splicing. Similar to reports for dcl2 and agl2, the expression levels for rdr3 and rdr4 increased after hypovirus CHV-1/EP713 infection, while expression levels of rdr1 and rdr2 were unchanged. The virus-responsive induction patterns for rdr3 and rdr4 were altered in the Δdcl2 and Δagl2 strains, suggesting some level of interaction between rdr3 and rdr4 and the dcl2/agl2 silencing pathway. Single rdr gene knockouts Δrdr1–4, double knockouts Δrdr1/2, Δrdr2/3, Δrdr1/3, and a triple knockout, Δrdr1/2/3, were generated and evaluated for effects on fungal phenotype, the antiviral defense response, viral RNA recombination activity and transposon expression. None of the single or multiple rdr knockout strains displayed any phenotypic differences from the parental strains with or without viral infection or any significant changes in viral RNA accumulation or recombination activity or transposon RNA accumulation, indicating no detectable contribution by the C. parasitica rdr genes to these processes.

Mutagenesis of the catalytic and cleavage site residues of the hypovirus papain-like proteases p29 and p48 reveals alternative processing and contributions to optimal viral RNA accumulation

The positive-stranded RNA genome of the prototypic virulence-attenuating hypovirus CHV-1/EP713 contains two open reading frames (ORF), each encoding an autocatalytic papain-like leader protease. Protease p29, derived from the N-terminal portion of ORF A, functions as a suppressor of RNA silencing, while protease p48, derived from the N-terminal portion of ORF B, is required for viral RNA replication. The catalytic and cleavage site residues required for autoproteolytic processing have been functionally mapped in vitro for both proteases but not confirmed in the infected fungal host. We report here the mutagenesis of the CHV-1/EP713 infectious cDNA clone to define the requirements for p29 and p48 cleavage and the role of autoproteolysis in the context of hypovirus replication. Mutation of the catalytic cysteine and histidine residues for either p29 or p48 was tolerated but reduced viral RNA accumulation to ca. 20 to 50% of the wild-type level. Mutation of the p29 catalytic residues caused an accumulation of unprocessed ORF A product p69. Surprisingly, the release of p48 from the ORF B-encoded polyprotein was not prevented by mutation of the p48 catalytic and cleavage site residues and was independent of p29. The results show that, while dispensable for hypovirus replication, the autocatalytic processing of the leader proteases p29 and p48 contributes to optimal virus RNA accumulation. The role of the predicted catalytic residues in autoproteolytic processing of p29 was confirmed in the infected host, while p48 was found to also undergo alternative processing independent of the encoded papain-like protease activities. Importance: Hypoviruses are positive-strand RNA mycoviruses that attenuate virulence of their pathogenic fungal hosts. The prototypic hypovirus CHV-1/EP713, which infects the chestnut bight fungus Cryphonetria parasitica, encodes two papain-like autocatalytic leader proteases, p29 and p48, that also have important functions in suppressing the RNA silencing antiviral defense response and in viral RNA replication, respectively. The mutational analyses of the CHV-1/EP713 infectious cDNA clone, reported here, define the requirements for p29 and p48 cleavage and the functional importance of autoproteolysis in the context of hypovirus replication and exposed an alternative p48 processing pathway independent of the encoded papain-like protease activities. These findings provide additional insights into hypovirus gene expression, replication, and evolution and inform ongoing efforts to engineer hypoviruses for interrogating and modulating fungal virulence.

A comparison of transcriptional patterns and mycological phenotypes following infection of Fusarium graminearum by four mycoviruses

Many fungi-infecting viruses, which are termed mycoviruses, have been identified, and most do not cause any visible symptoms. Some mycoviruses, however, can attenuate the virulence of the infected fungi, a phenomenon referred to as hypovirulence. To study fungus responses to virus infection, we established a model system composed of Fusarium graminearum and four mycoviruses including FgV1 (Fusarium graminearum virus 1), FgV2, FgV3, and FgV4. FgV1 and FgV2 infections caused several phenotypic alterations in F. graminearum including abnormal colony morphology, defects in perithecium development, and reductions in growth rate, conidiation, and virulence. In contrast, FgV3 and FgV4 infections did not cause any phenotypic change. An RNA-Seq-based analysis of the host transcriptome identified four unique Fusarium transcriptomes, one for each of the four mycoviruses. Unexpectedly, the fungal host transcriptome was more affected by FgV1 and FgV4 infections than by FgV2 and FgV3 infections. Gene ontology (GO) enrichment analysis revealed that FgV1 and FgV3 infections resulted in down-regulation of host genes required for cellular transport systems. FgV4 infection reduced the expression of genes involved in RNA processing and ribosome assembly. We also found 12 genes that were differentially expressed in response to all four mycovirus infections. Unfortunately, functions of most of these genes are still unknown. Taken together, our analysis provides further detailed insights into the interactions between mycoviruses and F. graminearum.

[Cryphonectria parasitica as a host of fungal viruses: a tool useful to unravel the mycovirus world]

There appear to be over a million of fungal species including those that have been unidentified and unreported, where a variety of viruses make a world as well. Studies on a very small number of them conducted during the last two decades demonstrated the infectivity of fungal viruses that had previously been assumed to be inheritable, indigenus and non-infectious. Also, great technical advances were achieved. The chest blight fungus (Cryphonectria parasitica), a phytopathogenic ascomycetous fungus, has emerged as a model filamentous fungus for fungal virology. The genome sequence with annotations, albeit not thorough, many useful research tools, and gene manipulation technologies are available for this fungus. Importantly, C. parasitica can support replication of homologous viruses naturally infecting it, in addition to heterologous viruses infecting another plant pathogenic fungus, Rosellinia necatrix taxonomically belonging to a different order. In this article, I overview general properties of fungal viruses and advantages of the chestnut blight fungus as a mycovirus host. Furthermore, I introduce two recent studies carried out using this fungal host:''Defective interfering RNA and RNA silencing that regulate the replication of a partitivirus'' and'' RNA silencing and RNA recombination''.

Hypovirus molecular biology: from Koch's postulates to host self-recognition genes that restrict virus transmission

The idea that viruses can be used to control fungal diseases has been a driving force in mycovirus research since the earliest days. Viruses in the family Hypoviridae associated with reduced virulence (hypovirulence) of the chestnut blight fungus, Cryphonectria parasitica, have held a prominent place in this research. This has been due in part to the severity of the chestnut blight epidemics in North America and Europe and early reports of hypovirulence-mediated mitigation of disease in European forests and successful application for control of chestnut blight in chestnut orchards. A more recent contributing factor has been the development of a hypovirus/C. parasitica experimental system that has overcome many of the challenges associated with mycovirus research, stemming primarily from the exclusive intracellular lifestyle shared by all mycoviruses. This chapter will focus on hypovirus molecular biology with an emphasis on the development of the hypovirus/C. parasitica experimental system and its contributions to fundamental and practical advances in mycovirology and the broader understanding of virus-host interactions and fungal pathogenesis.

A mycoreovirus suppresses RNA silencing in the white root rot fungus, Rosellinia necatrix

RNA silencing is a fundamental antiviral response in eukaryotic organisms. We investigated the counterdefense strategy of a fungal virus (mycovirus) against RNA silencing in the white root rot fungus, Rosellinia necatrix. We generated an R. necatrix strain that constitutively induced RNA silencing of the exogenous green fluorescent protein (GFP) gene, and infected it with each of four unrelated mycoviruses, including a partitivirus, a mycoreovirus, a megabirnavirus, and a quadrivirus. Infection with a mycoreovirus (R. necatrix mycoreovirus 3; RnMyRV3) suppressed RNA silencing of GFP, while the other mycoviruses did not. RnMyRV3 reduced accumulation of GFP-small interfering (si) RNAs and increased accumulation of GFP-double-stranded (ds) RNA; suggesting that the virus interferes with the dicing of dsRNA. Moreover, an agroinfiltration assay in planta revealed that the S10 gene of RnMyRV3 has RNA silencing suppressor activity. These data corroborate the counterdefense strategy of RnMyRV3 against host RNA silencing.

Functional diversity of RNAi-associated sRNAs in fungi

Yeast and filamentous fungi have been essential model systems for unveiling the secrets of RNA interference (RNAi). Research on these organisms has contributed to identifying general mechanisms and conserved eukaryotic RNAi machinery that can be found from fungi to mammals. The development of deep sequencing technologies has brought on the last wave of studies on RNAi in fungi, which has been focused on the identification of new types of functional small RNAs (sRNAs). These studies have discovered an unexpected diversity of sRNA, biogenesis pathways and new functions that are the focus of this review.

Small RNA populations for two unrelated viruses exhibit different biases in strand polarity and proximity to terminal sequences in the insect host Homalodisca vitripennis

Next generation sequence analyses were used to assess virus-derived small RNA (vsRNA) profiles for Homalodisca coagulata virus-1 (HoCV-1), family Dicistroviridae, and Homalodisca vitripennis reovirus (HoVRV), family Reoviridae, from virus-infected H. vitripennis, the glassy-winged sharpshooter. The vsRNA reads were mapped against the monopartite genome of HoCV-1 and all 12 genome segments of HoVRV, and 21nt vsRNAs were most common. However, strikingly contrasting patterns for the HoCV-1 and HoVRV genomic RNAs were observed. The majority of HoCV-1 vsRNAs mapped to the genomic positive-strand RNA and, although minor hotspots were observed, vsRNAs mapped across the entire genomic RNA. In contrast, HoVRV vsRNAs mapped to both positive and negative-sense strands for all genome segments, but different genomic segments showed distinct hotspots. The HoVRV vsRNAs were more common for 5' and 3' regions of HoVRV regions of all segments. These data suggest that taxonomically different viruses in the same host offer different targets for RNA-antiviral defense.