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

A circular single-stranded DNA mycovirus infects plants and confers broad-spectrum fungal resistance

Circular single-stranded DNA (ssDNA) viruses have been rarely found in fungi, and the evolutionary and ecological relationships among ssDNA viruses infecting fungi and other organisms remain unclear. In this study, a novel circular ssDNA virus, tentatively named Diaporthe sojae circular DNA virus 1 (DsCDV1), was identified in the phytopathogenic fungus Diaporthe sojae isolated from pear trees. DsCDV1 has a monopartite genome (3185 nt in size) encapsidated in isometric virions (21-26 nm in diameter). The genome comprises seven putative open reading frames encoding a discrete replicase (Rep) split by an intergenic region, a putative capsid protein (CP), several proteins of unknown function (P1-P4), and a long intergenic region. Notably, the two split parts of DsCDV1 Rep share high identities with the Reps of Geminiviridae and Genomoviridae, respectively, indicating an evolutionary linkage with both families. Phylogenetic analysis based on Rep or CP sequences placed DsCDV1 in a unique cluster, supporting the establishment of a new family, tentatively named Gegemycoviridae, intermediate to both families. DsCDV1 significantly attenuates fungal growth and nearly erases fungal virulence when transfected into the host fungus. Remarkably, DsCDV1 can systematically infect tobacco and pear seedlings, providing broad-spectrum resistance to fungal diseases. Subcellular localization analysis revealed that DsCDV1 P3 is systematically localized in the plasmodesmata, while its expression in trans-complementation experiments could restore systematic infection of a movement-deficient plant virus, suggesting that P3 is a movement protein. DsCDV1 exhibits unique molecular and biological traits not observed in other ssDNA viruses, serving as a link between fungal and plant ssDNA viruses and presenting an evolutionary connection between ssDNA viruses and fungi. These findings contribute to expanding our understanding of ssDNA virus diversity and evolution, offering potential biocontrol applications for managing crucial plant diseases.

New clades of viruses infecting the obligatory biotroph Bremia lactucae representing distinct evolutionary trajectory for viruses infecting oomycetes

Recent advances in high throughput sequencing (HTS) approaches allowed a broad exploration of viromes from different fungal hosts, unveiling a great diversity of mycoviruses with interesting evolutionary features. The word mycovirus historically applies also to viruses infecting oomycetes but most studies are on viruses infecting fungi, with less mycoviruses found and characterized in oomycetes, particularly in the obligatory biotrophs. We, here, describe the first virome associated to Bremia lactucae, the causal agent of lettuce downy mildew, which is an important biotrophic pathogen for lettuce production and a model system for the molecular aspects of the plant-oomycetes interactions. Among the identified viruses, we could detect (1) two new negative sense ssRNA viruses related to the yueviruses, (2) the first example of permuted RdRp in a virus infecting fungi/oomycetes, (3) a new group of bipartite dsRNA viruses showing evidence of recent bi-segmentation and concomitantly, a possible duplication event bringing a bipartite genome to tripartite, (4) a first representative of a clade of viruses with evidence of recombination between distantly related viruses, (5) a new open reading frame (ORF)an virus encoding for an RdRp with low homology to known RNA viruses, and (6) a new virus, belonging to riboviria but not conserved enough to provide a conclusive phylogenetic placement that shows evidence of a recombination event between a kitrinoviricota-like and a pisuviricota-like sequence. The results obtained show a great diversity of viruses and evolutionary mechanisms previously unreported for oomycetes-infecting viruses, supporting the existence of a large diversity of oomycetes-specific viral clades ancestral of many fungal and insect virus clades.

Discovery and Identification of Viruses Infecting Oomycetes

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

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.

Successful full-length genomic cloning and characterization of site-specific nick structures of Phytophthora endornaviruses 2 and 3 in yeast, Saccharomyces cerevisiae

Two endornaviruses, Phytophthora endornavirus 2 (PEV2) and Phytophthora endornavirus 3 (PEV3), have been discovered in pathogens targeting asparagus. In this study, we analyzed the nick structure in the RNA genomes of PEV2 and PEV3 in the host oomycetes. Northern blot hybridization using positive and negative strand-specific RNA probes targeting the 5' and 3' regions of PEV2 and PEV3 RNA genomes revealed approximately 1.0 kilobase (kb) RNA fragments located in the 5' regions of the two genomes. 3' RACE analysis determined that the size of the RNA fragments were 958 nucleotides (nt) for PEV2 and 968 nt for PEV3. We have successfully constructed full-length cDNA clones of the entire RNA genomes of PEV2 and PEV3 using a homologous recombination system in the yeast, Saccharomyces cerevisiae. These full-length cDNA sequences were ligated downstream of a constitutive expression promoter (TDH3) or a galactose-inducing promoter (GAL1) in the shuttle vector to enable the production of the full-length RNA transcripts of PEV2 and PEV3 in yeast cells. Interestingly, a 1.0 kb RNA fragment from the PEV3 positive-strand transcript was also detected with a 5'-region RNA probe, indicating that site-specific cleavage also occurred in yeast cells. Further, when PEV2 or PEV3 mRNA was overexpressed under the GAL1 promoter, yeast cell growth was suppressed. A fusion protein combining EGFP to the N-terminus of the full-length PEV2 ORF or C-terminus of the full-length PEV3 ORF was expressed, and allowed PEV2 and PEV3 ORFs to be successfully visualized in yeast cells. Expression of the fusion protein also revealed presence of heterogeneous bodies in the cells.

Sequence and phylogenetic analysis of a novel alphaendornavirus, the first virus described from the oomycete plant pathogen Phytophthora heveae

Here, we report the discovery and complete genome sequence of a novel virus, designated as "Phytophthora heveae alphaendornavirus 1" (PhAEV1), from a single isolate of the plant pathogenic oomycete Phytophthora heveae (kingdom Stramenipila) isolated from a tropical evergreen lowland rainforest in northern Vietnam. PhAEV1 was detected by both cellulose affinity chromatography of dsRNA and high-throughput sequencing of total RNA, and its presence and sequence were confirmed by RT-PCR and Sanger sequencing. The PhAEV1 genome, 12,820 nucleotides (nt) in length, was predicted to encode a single large polyprotein with the catalytic core domain of viral (superfamily 1) RNA helicase (HEL, amino acid [aa] positions 1,287-1,531), glycosyltransferase (GT, aa positions ca. 2,800-3,125), and RNA-directed RNA polymerase (RdRp, aa positions 3,875-4,112). PhAEV1 is the most similar to Phytophthora cactorum alphaendornavirus 3, sharing 39.4% and 39.1% nt and aa sequence identity, respectively. In addition to the first 5'-terminal AUG codon, three additional in-frame methionine codons were found in close proximity (nt 14-16, 96-98, and 176-178), suggesting potential additional translation initiation sites. Conserved RdRp motifs (A-E) similar to those detected in related endornaviruses were identified in PhAEV1, as well as in several previously described alphaendornaviruses from other Phytophthora species in which these motifs had not been identified previously. Phylogenetic analysis showed that PhAEV1 clusters with members of the genus Alphaendornavirus in the family Endornaviridae and is basal to two other alphaendornaviruses described from another oomycete, Phytophthora cactorum. PhAEV1 is the first virus reported in P. heveae.

Identification and Molecular Characterization of a Novel Carlavirus Infecting Chrysanthemum morifolium in China

Chrysanthemum (Chrysanthemum morifolium) is an important ornamental and medicinal plant suffering from many viruses and viroids worldwide. In this study, a new carlavirus, tentatively named Chinese isolate of Carya illinoinensis carlavirus 1 (CiCV1-CN), was identified from chrysanthemum plants in Zhejiang Province, China. The genome sequence of CiCV1-CN was 8795 nucleotides (nt) in length, with a 68-nt 5'-untranslated region (UTR) and a 76-nt 3'-UTR, which contained six predicted open reading frames (ORFs) that encode six corresponding proteins of various sizes. Phylogenetic analyses based on full-length genome and coat protein sequences revealed that CiCV1-CN is in an evolutionary branch with chrysanthemum virus R (CVR) in the Carlavirus genus. Pairwise sequence identity analysis showed that, except for CiCV1, CiCV1-CN has the highest whole-genome sequence identity of 71.3% to CVR-X6. At the amino acid level, the highest identities of predicted proteins encoded by the ORF1, ORF2, ORF3, ORF4, ORF5, and ORF6 of CiCV1-CN were 77.1% in the CVR-X21 ORF1, 80.3% in the CVR-X13 ORF2, 74.8% in the CVR-X21 ORF3, 60.9% in the CVR-BJ ORF4, 90.2% in the CVR-X6 and CVR-TX ORF5s, and 79.4% in the CVR-X21 ORF6. Furthermore, we also found a transient expression of the cysteine-rich protein (CRP) encoded by the ORF6 of CiCV1-CN in Nicotiana benthamiana plants using a potato virus X-based vector, which can result in a downward leaf curl and hypersensitive cell death over the time course. These results demonstrated that CiCV1-CN is a pathogenic virus and C. morifolium is a natural host of CiCV1.

Discovery and Characterization of Putative Glycoprotein-Encoding Mycoviruses in the Bunyavirales

Although segmented negative-sense RNA viruses (SNSRVs) have been frequently discovered in various fungi, most SNSRVs reported only the large segments. In this study, we investigated the diversity of the mycoviruses in the phytopathogenic fungus Fusarium asiaticum using the metatranscriptomic technique. We identified 17 fungal single-stranded RNA (ssRNA) viruses including nine viruses within Mitoviridae, one each in Narnaviridae, Botourmiaviridae, Hypoviridae, Fusariviridae, and Narliviridae, two in Mymonaviridae, and one trisegmented virus temporarily named Fusarium asiaticum mycobunyavirus 1 (FaMBV1). The FaMBV1 genome comprises three RNA segments, large (L), medium (M), and small (S) with 6,468, 2,639, and 1,420 nucleotides, respectively. These L, M, and S segments putatively encode the L protein, glycoprotein, and nucleocapsid, respectively. Phylogenetic analysis based on the L protein showed that FaMBV1 is phylogenetically clustered with Alternaria tenuissima negative-stranded RNA virus 2 (AtNSRV2) and Sclerotinia sclerotiorum negative-stranded RNA virus 5 (SsNSRV5) but distantly related to the members of the family Phenuiviridae. FaMBV1 could be vertically transmitted by asexual spores with lower efficiency (16.7%, 2/42). Comparison between FaMBV1-free and -infected fungal strains revealed that FaMBV1 has little effect on hyphal growth, pathogenicity, and conidium production, and its M segment is dispensable for viral replication and lost during subculture and asexual conidiation. The M and S segments of AtNSRV2 and SsNSRV5 were found using bioinformatics methods, indicating that the two fungal NSRVs harbor trisegmented genomes. Our results provide a new example of the existence and evolution of the segmented negative-sense RNA viruses in fungi. IMPORTANCE Fungal segmented negative-sense RNA viruses (SNSRVs) have been frequently found. Only the large segment encoding RNA-dependent RNA polymerase (RdRp) has been reported in most fungal SNSRVs, except for a few fungal SNSRVs reported to encode nucleocapsids, nonstructural proteins, or movement proteins. Virome analysis of the Fusarium spp. that cause Fusarium head blight discovered a novel virus, Fusarium asiaticum mycobunyavirus 1 (FaMBV1), representing a novel lineage of the family Phenuiviridae. FaMBV1 harbors a trisegmented genome that putatively encodes RdRp, glycoproteins, and nucleocapsids. The putative glycoprotein was first described in fungal SNSRVs and shared homology with glycoprotein of animal phenuivirus but was dispensable for its replication in F. asiaticum. Two other trisegmented fungal SNSRVs that also encode glycoproteins were discovered, implying that three-segment bunyavirus infections may be common in fungi. These findings provide new insights into the ecology and evolution of SNSRVs, particularly those infecting fungi.

Mycoviruses as a part of the global virome: Diversity, evolutionary links and lifestyle

Knowledge of mycovirus diversity, evolution, horizontal gene transfer and shared ancestry with viruses infecting distantly related hosts, such as plants and arthropods, has increased vastly during the last few years due to advances in the high throughput sequencing methodologies. This also has enabled the discovery of novel mycoviruses with previously unknown genome types, mainly new positive and negative single-stranded RNA mycoviruses ((+) ssRNA and (-) ssRNA) and single-stranded DNA mycoviruses (ssDNA), and has increased our knowledge of double-stranded RNA mycoviruses (dsRNA), which in the past were thought to be the most common viruses infecting fungi. Fungi and oomycetes (Stramenopila) share similar lifestyles and also have similar viromes. Hypothesis about the origin and cross-kingdom transmission events of viruses have been raised and are supported by phylogenetic analysis and by the discovery of natural exchange of viruses between different hosts during virus-fungus coinfection in planta. In this review we make a compilation of the current information on the genome organization, diversity and taxonomy of mycoviruses, discussing their possible origins. Our focus is in recent findings suggesting the expansion of the host range of many viral taxa previously considered to be exclusively fungal, but we also address factors affecting virus transmissibility and coexistence in single fungal or oomycete isolates, as well as the development of synthetic mycoviruses and their use in investigating mycovirus replication cycles and pathogenicity.

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

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

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.

Natural Populations from the Phytophthora palustris Complex Show a High Diversity and Abundance of ssRNA and dsRNA Viruses

We explored the virome of the “Phytophthora palustris complex”, a group of aquatic specialists geographically limited to Southeast and East Asia, the native origin of many destructive invasive forest Phytophthora spp. Based on high-throughput sequencing (RNAseq) of 112 isolates of “P. palustris” collected from rivers, mangroves, and ponds, and natural forests in subtropical and tropical areas in Indonesia, Taiwan, and Japan, 52 putative viruses were identified, which, to varying degrees, were phylogenetically related to the families Botybirnaviridae, Narnaviridae, Tombusviridae, and Totiviridae, and the order Bunyavirales. The prevalence of all viruses in their hosts was investigated and confirmed by RT-PCR. The rich virus composition, high abundance, and distribution discovered in our study indicate that viruses are naturally infecting taxa from the “P. palustris complex” in their natural niche, and that they are predominant members of the host cellular environment. Certain Indonesian localities are the viruses’ hotspots and particular “P. palustris” isolates show complex multiviral infections. This study defines the first bi-segmented bunya-like virus together with the first tombus-like and botybirna-like viruses in the genus Phytophthora and provides insights into the spread and evolution of RNA viruses in the natural populations of an oomycete species.

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.