Figures
Abstract
Satellites associated with plant or animal viruses have been largely detected and characterized, while those from mycoviruses together with their roles remain far less determined. Three dsRNA segments (dsRNA 1 to 3 termed according to their decreasing sizes) were identified in a strain of phytopathogenic fungus Pestalotiopsis fici AH1-1 isolated from a tea leaf. The complete sequences of dsRNAs 1 to 3, with the sizes of 10316, 5511, and 631 bp, were determined by random cloning together with a RACE protocol. Sequence analyses support that dsRNA1 is a genome of a novel hypovirus belonging to genus Alphahypovirus of the family Hypoviridae, tentatively named Pestalotiopsis fici hypovirus 1 (PfHV1); dsRNA2 is a defective RNA (D-RNA) generating from dsRNA1 with septal deletions; and dsRNA3 is the satellite component of PfHV1 since it could be co-precipitated with other dsRNA components in the same sucrose fraction by ultra-centrifuge, suggesting that it is encapsulated together with PfHV1 genomic dsRNAs. Moreover, dsRNA3 shares an identical stretch (170 bp) with dsRNAs 1 and 2 at their 5′ termini and the remaining are heterogenous, which is distinct from a typical satellite that generally has very little or no sequence similarity with helper viruses. More importantly, dsRNA3 lacks a substantial open reading frame (ORF) and a poly (A) tail, which is unlike the known satellite RNAs of hypoviruses, as well as unlike those in association with Totiviridae and Partitiviridae since the latters are encapsidated in coat proteins. As up-regulated expression of RNA3, dsRNA1 was significantly down-regulated, suggesting that dsRNA3 negatively regulates the expression of dsRNA1, whereas dsRNAs 1 to 3 have no obvious impact on the biological traits of the host fungus including morphologies and virulence. This study indicates that PfHV1 dsRNA3 is a special type of satellite-like nucleic acid that has substantial sequence homology with the host viral genome without encapsidation in a coat protein, which broadens the definition of fungal satellite.
Author summary
A novel hypovirus belonging to the genus Alphahypovirus of the family Hypoviridae was identified in a phytopathogenetic P. fici isolated from a tea leaf and tentatively named Pestalotiopsis fici hypovirus 1 (PfHV1), which represents the first hypovirus from Pestalotiopsis fungi. Three RNA components (RNA 1 to 3 termed according to their decreasing sizes, dsRNA 1 to 3 termed for their replicative form are corresponding with RNA 1 to 3) associated with PfHV1 including the genomic RNA, defective RNA and satellite-like RNA (SatL-RNA) were determined, and their molecular and biological traits were characterized. The SatL-RNA negatively regulates the accumulation of PfHV1 RNA1, and represents a novel class of satellite nucleic acids since it has substantial sequence homology with the host viral genome and is not encapsidated in the coat protein of its helper virus. This study contributes to our better understanding of fungal satellites, as well as the molecular and biological traits of hypoviruses.
Citation: Han Z, Liu J, Kong L, He Y, Wu H, Xu W (2023) A special satellite-like RNA of a novel hypovirus from Pestalotiopsis fici broadens the definition of fungal satellite. PLoS Pathog 19(6): e1010889. https://doi.org/10.1371/journal.ppat.1010889
Editor: Aiming Wang, Agriculture and Agri-Food Canada, CANADA
Received: September 21, 2022; Accepted: May 23, 2023; Published: June 7, 2023
Copyright: © 2023 Han et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: PfHV1 dsRNAs 1 to 3 were deposited in GenBank with accession numbers OP441373-OP441375, and all relevant data are within the manuscript and its Supporting Information files.
Funding: This work was financially supported by the National Key R&D Program of China (No. 2021YFD1000400), National Natural Science Foundation of China (No. 31872014 and 32172475), and Enshi Tujia and Miao Autonomous Prefecture Bureau of Science and Technology (No. XYJ2021000083) to W.X. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors declare no competing interests.
Introduction
Hypoviridae is a family of capsidless viruses with positive-sense (+), single-stranded (ss) RNA genomes of 9.1–12.7 kb that possess either a single large ORF or another small ORF [1]. Several years ago, only one genus, Hypovirus, is included in this family, which contains Cryphonectria hypovirus 1–4 (CHV1-4) and among them, CHV1 is associated with biologic control of the filamentous fungus that causes chestnut blight [2]. With the number increase in recent years, the genus Hypovirus has been rapidly expanded into eight genera, namely Alpha-, Beta-, Gamma-, Delta-, Epsilon-, Zeta-, Eta-, and Thetahypovirus [3–5]. Hypoviruses have been heavily detected in ascomycetous and basidiomycetous filamentous fungi, and some of these were able to alter fungal host phenotypes and attenuate the host virulence, e.g., Alternaria alternata hypovirus 1 (AaHV1) [6], Botrytis cinerea hypovirus 1 (BcHV1) [7], while others were not, e.g., Fusarium graminearum hypovirus 1 (FgHV1) [8]. Defective RNAs (D-RNAs) were generally detected in association with hypoviruses, e.g., CHV1-EP713 [9,10], Sclerotinia sclerotiorum hypovirus 1 (SsHV1) [11], Fusarium graminearum Hypovirus 2 (FgHV2) [12], Botrytis cinerea hypovirus 1 (BcHV1) [7], AaHV1 [6], which were a type of generations of internally deleted genomic RNAs.
Satellites are subviral agents which lack genes that could encode functions needed for replication, and depend on the co-infection of a host cell with a helper virus, which consist of satellite viruses and satellite nucleic acids, discriminating by the presence of a structural protein [13]. Satellite nucleic acids include ssDNAs and ssRNAs, associated with plant or animal viruses, and double stranded RNAs (dsRNAs), are mainly associated with mycoviruses, namely the families Totiviridae and Partitiviridae. Satellites do not constitute a homogeneous taxonomic group, and have no taxonomic correlation with their helper viruses, and have arisen independently a number of times during virus evolution [14]. Since satellites were mostly characterized as ssRNA satellites that use ssRNA plant viruses as helpers, and those associated with mycoviruses remain rarely characterized.
Pestalotiopsis spp. belong to the family Amphisphaeriaceae that are distributed widely throughout tropical and temperate regions. They are common phytopathogens reducing production and causing enormous losses in either horticultural plants or forest trees [15]. For example, in tea plants, P. theae and P. longiseta are responsible for tea gray blight disease, which leads to large brown spots with the formation of apparent concentric rings in the late period of disease on mature and old leaves of tea plants, and even leads to death of tender branches of infected tea seedlings; P. theae, P. camelliae, and P. clavispora were also identified in association with small brown-black spots on tea tender leaves [16,17]. However, endophytic Pestalotiopsis species can confer fitness to the host plants, and are also thought to be a rich source for bioprospecting compared to other fungal genera since they produce novel compounds with medicinal, agricultural and industrial applications [17]. Therefore, Pestalotiopsis fungi live in host plants with a complex lifestyle, being parasitic, commensal, or mutualistic. Mycoviruses might have participated in the transition of these lifestyles, since a chrysovirus, the first mycovirus isolated from Pestalotiopsis fungi, can convert its host fungus from a phytopathogenetic one to a nonpathogenetic endophyte on tea leaves, conferring high resistance to the host plants against the virulent P. theae strains [15]. It is fascinating to characterize more mycoviruses to understand the virus taxonomy, evolution, molecular and biological traits related to this important fungal species. Specifically, characterization of a hypovirus in this fungal species remains a great interesting since hypoviruses represent the classic and successful examples for biological control of plant diseases [18].
In this study, a novel special satellite RNA associated with a novel hypovirus from P. fici with three components (RNAs 1 to 3) was identified and characterized, which represent a novel class of satellite nucleic acids, illustrating some novel molecular and biological traits of a fungal satellite.
Materials and methods
Fungal strains and cultures
P. fici strain AH1-1 was isolated from a tea leaf showing the typical symptoms of tea grey blight disease collected in Liu’an city, Anhui province, China, and was identified based on morphologies and multi-loci sequences as previously described [19], and strain AH1-1V-, AH1-1-14 and AH1-1-24 are subisolates generated from single conidia of strain AH1-1. Strain AH1-1V- is a cured subisolate of AH1-1 by elimination of dsRNA components, and conferred hygromycin (hyg)-resistant ability by transformation of a hyg gene to generate strain AH1-1hygV-. P. theae TP-2-2W, JWX-3-1, and CJB-4-1 were collected from tea leaves in Yichang city, Hubei province, China, and were identified on their morphologies and internally transcribed spacer (ITS) sequence. P. theae strain LI41 was isolated from a healthy tea leaf collected in Xuan’en county, En’shi prefecture, Hubei province, China, and LI41-1P1 is a subisolate of P. theae LI41, both of which were infected by Pestalotiopsis theae chrysovirus-1 (PtCV1) [15]. All strains were grown at 25°C in the dark for 3–5 days on potato dextrose agar (PDA; 20% diced potatoes, 2% glucose, and 1.5% agar) unless otherwise stated.
dsRNA extraction, purification, and enzymatic treatments
For dsRNA extraction, fungal mycelial plugs were inoculated onto sterilized cellophane disks on PDA plates at 25°C in the dark for 4 to 5 days, and the mycelia were collected and subjected to dsRNA extraction as previously described [20]. The resulting nucleic acids were treated with 2 U DNase I (New England Biolabs), 10 U S1 nuclease (Thermo Scientific) at 37°C for 1 h, or those without treatments were fractionated by electrophoresis on 1.2% agarose gels with Tris-acetate-EDTA (TAE) buffer and visualized by staining with ethidium bromide. Each dsRNA band was excised, purified using a gel extraction kit (Qiagen), dissolved in DEPC-treated water and stored at −80°C until use.
The dsRNA preparation, BdRV1 dsRNAs (extracted from mycelia infected by BdRV1) [21], in vitro transcripts of complementary DNAs (cDNAs) of peach latent mosaic viroid (PLMVd) [16] with its linearized pGEM-T plasmids harboring the cDNAs, and DNA controls generated from polymerase chain reaction (PCR) products were treated with S1 nuclease or/and DNase I as indicated above, and with 200 ng/mL RNase A (Thermo Scientific) in 2×saline sodium citrate (SSC) (300 mM NaCl, 30 mM sodium citrate, pH 7.0) or 0.1×SSC as previously described [20].
Cloning, sequencing, and sequence analysis
The sequences of the dsRNAs were determined by cloning and sequencing amplicons generated by reverse transcription and polymerase chain reaction (RT-PCR) using the random primers 05RACE-3RT and 05RACE-3 (S1 Table) as previously described [19]. The 5′- and 3′-terminal sequences of the dsRNAs were obtained by cloning and sequencing the RT-PCR amplicons generated using a standard RNA ligase mediated rapid amplification of cDNA ends (RLM-RACE) protocol, including the usage of PC3-T7loop and PC2 (S1 Table). The oligonucleotide primers used for RLM-RACE were designed based on sequence information obtained from the randomly primed amplicons [22]. At least three independent clones of each amplicon were sequenced in both directions, by Sangon Biotech Co., Ltd, Shanghai, China. Sequence similarity searches were performed using the BLASTn program for nucleic acids or BLASTp for putative proteins against the National Center for Biotechnology Information (NCBI) databases. Multiple alignments of nucleic and amino acid (aa) sequences were conducted using MEGA7 and MAFFT online (MAFFT alignment and NJ / UPGMA phylogeny (cbrc.jp)). The phylogenetic tree for RNA-dependent RNA polymerase (RdRp) sequences was constructed using MEGA 11 with the Maximum Likelihood method [23]. Open reading frames (ORFs) were deduced using an ORF finder (https://www.ncbi.nlm.nih.gov/orffinder/). Conserved domains were predicted using Phyre2 (https://www.sbg.bio.ic.ac.uk/phyre2). Internal ribosomal entry site (IRES) elements were predicted by aligning with IRESite database online (IRESite: The database of experimentally verified IRES structures) [24].
RT-PCR detection, RNA blot analysis, co-precipitate, and RNase A treatment
RT-PCR amplification was performed using a specific primer pair derived from the dsRNA 1 sequence (PfHV1-1F: 5′-TTCGATTTCAACGCCAGGTC-3′; PfHV1-1R: 5′-GCCGGGTCTATCGTCTTTTC-3′) (S1 Table), generating a 307-bp fragment with an annealing temperature of 56°C in a PCR Thermal Cycler (Model PTC-100, MJ Research).
Digoxigenin (DIG)-labeled riboprobes 1 (binding to positions from 1588 to 1934 nt of dsRNA1), 2 and 3 (116–355 and 177–485 nt of dsRNA3, respectively) were synthesized in vitro transcription based on the corresponding cDNAs inserted into the pGEM-T easy vector with T7 RNA polymerase (Takara, Beijing, China) as previously described [25]. Nucleic acids were spotted (for dot blotting) or electro-transferred (for Northern blot) to positively-charged nylon membranes (Roche Diagnostics), and hybridized with DIG-labeled riboprobes as previously described [25].
Co-precipitate analysis of dsRNAs was conducted by sucrose gradient centrifugation according to the methods used for the purification of viral particles previously described [15]. The resulted aliquots of each fraction (100 μL) were subjected to dsRNA extraction to monitor for the presence of viral dsRNAs.
RNase A treatment was conducted as previously described [26]. Briefly, fungal mycelia were cultured on cellophane-PDA for 3 days at room temperature and powdered in the presence of liquid nitrogen. The mycelial powder was homogenized in 2 volumes (vol/wt) of 0.1 M sodium phosphate (pH 7.2) and centrifugated at 10000 r/min for 10 min at 4°C. The supernatant was treated with RNase A (10 μg/mL) (Thermo scientific, Shanghai) for 30 min at 37°C. The solution was extracted with phenol and chloroform to inactivate RNase A, followed by chloroform extraction.
Engineering of RNA3-over-expressed strain and quantitative real-time PCR (qRT-PCR) analysis
To generate RNA3-over-expressed strains, RNA3 cDNAs were inserted into the multi-cloning sites of vector pCAMBgfp after linearization by Pst I [27]. The constructed vector was then transformed to the strain AH1-1-24 via PEG4000 mediated protoplasts transformation as previously described [15].
The total RNAs from mycelia were extracted by Trizol agent and subjected to digestion by 5×gDNA digesture (Foregene, Chengdu, China) to remove the genomic DNAs, followed by cDNA synthesis using 2× reverse transcriptase (Foregene, Chengdu, China) and amplification using 2×TSINGKE master qPCR mix (Tsingke, Beijing, China) using specific primers designed for quantitative analysis in CFX96 Real-Time PCR Detection System (Bio-Rad, USA). Amplification was conducted for 40 cycles of 95°C for 10 s, 56°C for 10 s and 72°C for 15 s after initial heating at 95°C for 1 min. A partial fragment of the Actin gene was used to normalize the RNA samples for each qRT-PCR, and each treatment was conducted with three technical replicates.
Contact cultures of P. fici isolates
Horizontal transmission of the virus was conducted as previously described [20]. Briefly, both donor and recipient strains were cultured together on the same Petri dishes at 25°C for 7 days, and allowed to physically contact each other. Following contact, mycelial agar plugs were excised from the contact area of two strains and sub-cultured onto fresh PDA plates. Seven independent donor-recipient pairs were assessed, and two mycelial agar plugs were selected from each pair for further analysis, resulting in a total of 14 isolates.
Growth rate, morphology, and virulence assays
Fungal growth rates and morphologies were assessed as previously described [15]. Three biological replicates for each strain were monitored and the resulting data were subjected to statistical analysis as described below. Fungal virulence was determined following inoculation of detached tea leaves (Camellia sinensis, C. sinensis vars. Echa no.1 or Fuyun no.6) as previously described [15]. At 4 day past inoculation (dpi), lesions that developed on the inoculated leaves were measured. Eight biological replicates for each strain were monitored and the resulting data were subjected to statistical analysis as described below.
Statistical analyses
The biological data were statistical analysis using SPSS Statistics 17.0 with one-way analysis of variance (ANOVA) and Tukey post-hoc tests (Win Wrap Basic; http://www.winwrap.com). The mean values for the biological replicates are presented as column charts with error bars representing standard error of mean (SEM). The graphs were produced in Excel (Microsoft) and GraphPad Prism 7 (GraphPad software). P-values < 0.05 were considered to indicate statistical significance, while P-values < 0.01, extremely significant.
Results
A complex of dsRNA segments were extracted from P. fici strain AH1-1
Nucleic acid preparations enriched in dsRNA were obtained from P. fici strains AH1-1 and TP-2-2W isolated from tea leaves collected in China, and were subjected to digestion with DNase I and S1 nuclease and then agarose gel electrophoresis. The results showed that three dsRNAs (nominated 1–3 according to their decreasing sizes) were detected in preparations of strain AH1-1 but not in a control strain TP-2-2W (Fig 1B). Of these components, dsRNA1 migrated together with the fungal genomic DNAs as electrophoresis on an agarose gel, and became apparent after treatment by DNase I (Fig 1B). The sequences of the full-length cDNAs of dsRNAs 1–3 were determined by assembling partial-length cDNAs that were amplified from the purified dsRNAs using RT-PCR and RLM-RACE protocols. The corresponding sequences were deposited in GenBank with accession numbers OP441373-OP441375.
(A) Colonies of Pestalotiopsis strain AH1-1 and TP-2-2W grown on PDA medium for 6 days. (B) Electrophoresis analysis of nucleic acids extracted from strain AH1-1 (left panel) and TP-2-2W (right) without treatment (lane 2) and treated with S1 nuclease (3) and DNase I (4), respectively, on 1% agarose gel. (C) Determined the dsRNA nature by enzymatic treatment with S1 nuclease, DNase I and RNase A. Botryosphaeria dothidea RNA virus 1 (BdRV 1), peach latent mosaic viroid (PLMVd), and PCR products were involved as dsRNA, ssRNA and DNA controls, respectively. M, DNA size marker.
The dsRNA nature of the three observed bands was further assessed by treatments with DNase I, S1 nuclease or RNase A (in 2× and 0.1×SSC), together with an ssRNA control (in vitro dimeric transcripts of peach latent mosaic viroid (PLMVd) and dsRNAs from a dsRNA mycovirus (Botryosphaeria dothidea RNA virus 1, BdRV 1). The RNAs extracted from strain AH1-1 together with BdRV 1 dsRNAs were digested by RNase A in 0.1×SSC, but they resisted digestion by DNase I, S1 nuclease, and RNase A in 2×SSC. In sharp contrast, PLMVd transcripts were completely degraded by S1 nuclease and by RNase A under both ionic conditions, but resisted digestion by DNase I; the genomic DNAs (extracted together with BdRV 1) and the cDNA plasmid of PLMVd (used for in vitro transcription) were completely degraded by DNase I while resistant to RNase A and S1 nuclease (Fig 1C). These data strongly support that the nucleic acids extracted from strain AH1-1 were indeed dsRNAs instead of DNAs or ssRNAs.
dsRNA 1 composes the genome of a novel hypovirus
The complete nucleotide sequence of dsRNA1 is 10316 base pair (bp) in length, consisting of a single putative ORF beginning at AUG (nt positions 513–515) and terminating at UAG (9600–9602), coding for a polyprotein of 3029 aa with an approximate molecular mass of 345.9 kDa.
The deduced polyprotein aa sequence contains conserved domains of RNA dependent RNA polymerase (RdRp), viral RNA Helicase (Hel), and a most likely papain-like protease (Pro) (Fig 2B). Within the RdRp domain, nine conserved motifs (Ia-VII), characteristic for an RNA virus [28], were detected (Fig 3A), and seven conserved motifs (I-VI) relatively conserved in (+) ssRNA were detected within the Hel domain (Fig 3C); while in the possible Pro domain, only cysteine and glycine were detected whereas the core residue histidine was absent (Fig 3B). The 5′-untranslated region (UTR) of dsRNA1 is 512 nts in size, contains six AUG codons (S1C Fig), and the 3′-UTR is 714 nts including a 6 nt-sized adenine tail (poly A) (Fig 2B). Each UTR region forms a compact and complex stem-loop structure as predicted in RNAfold (S1A and S1B Fig).
(A) An ML phylogenetic tree constructed based on the polyprotein sequences encoded by PfHV1 ORF1 and those of known hypoviruses. (B) Genomic organization of the ORF1-encoding polyprotein (grey color). The conserved motifs for protease (Pro), RNA dependent RNA polyprotein (RdRp) and helicase (Hel) domain blocks on the polyprotein are marked in dark red, red and purple respectively, with the lengths corresponding to their aa sizes. The numbers under the line indicate the start and end positions of genome, 5′- and 3′-untranslated regions (UTRs), and the conserved domains.
(A to C) The conserved domains for RdRp, Pro, Hel. The conserved motifs or residues are highlighted.
BLASTp searches revealed that the polypeptide has the highest identity (55.13%) with the polyprotein (accession no. YP_009342443.1, coverage 96%, E value = 0) encoded by Wuhan insect virus 14 (WIV14), as well as high identities (50.25–55.10%, coverage 55–92%, E value = 0) to those of Fusarium sacchari hypovirus 1 (FsHV1; no. QIQ28422.1;) [29], Apis hypovirus 1–3 (no. UCR92523.1, UCR92524.1, UCR92525.1; direct submission), and Alternaria alternata hypovirus 1 (AaHV1; no. QFR36339.1; S2 Table)[6]. Phylogenetic analysis of the polyprotein sequence with those of all the known members of the family Hypoviridae illustrated that the sequence clustered with FsHV1 and AaHV1 in the genus Alphahypovirus (Fig 2A, [3,6]. These results suggest that dsRNA1 is a genomic component of a novel hypovirus, and it is tentatively named Pestalotiopsis fici hypovirus 1 (PfHV1).
dsRNA2 is a D-RNA of PfHV1
The complete nucleotide sequence of dsRNA2 is 5511 bp in length excluding the poly (A) tail, consisting of a single putative ORF (nt positions 510–4043), coding for a polyprotein of 1177 aa with an approximate molecular mass of 133.5 kDa. Alignment with dsRNA1 indicated that the total sequence of dsRNA2 is separately identical with dsRNA1 in five regions, i.e., nt position 4 to 3939, 5988 to 6225, 6220 to 6508, 6512 to 6770, and 6774 to 7562. Namely, dsRNA2 is a D-RNA generated from dsRNA1 by deletion three regions from nts 1 to 3, 3940 to 5987 and 7563 to 10316 (Fig 4A). Moreover, an inversion was observed in dsRNA2 sequence between nt position 4175 to 4433 and 4434 to 4722 as compared with dsRNA1 (Fig 4A).
(A) Schematic diagram of the corresponding positions of dsRNA1, D-RNA, and SatL-RNA sequences. The purple and black blocks indicate the identical and deleted regions in dsRNA1 and D-RNA, respectively, and yellow and pink lines indicate the inversed regions in them. The red and black blocks indicate identical and exogenous regions in SatL-RNA and dsRNA1 or and D-RNA, respectively. (B) Northern blotting analysis of dsRNA1, D-RNA and SatL-RNA using riboprobe 1 to 3, whose corresponding positions in these dsRNAs are marked in green, blue, and grey, respectively. (C) P. fici strain AH1-1V- free of PfHV1, P. fici strain AH1-1 infected by PfHV1, and P. theae strain LI41-1P1 infected by PtCV1 (encapsidated virion, a chrysovirus) were used. Mycelial homogenates were incubated at 37°C or on ice for 30 min in the presence or absence of RNase A. Subsequently, they were electrophoresed on a 1.5% agarose gel after RNA extraction with chloroform. (D) Co-precipitation analysis of dsRNA1, D-RNA and SatL-RNA by ultra-centrifugation in stepwise sucrose gradients (from 10% to 35% sucrose with 5% increments). (E) Relative expression analysis of dsRNA1 and SatL-RNA using RT-qPCR. Data are means±SEM (n = 3). Different letter indicates a significant difference at p < 0.05 (one-way ANOVA).
To further confirm that dsRNA2 as a D-RNA generated from dsRNA1, a digoxigenin (DIG)-labeled riboprobe (termed riboprobe 1) binding to the nt 1588 to 1934 of dsRNA1 was synthesized in vitro transcription and hybridization with the nucleic acid preparations of P. fici strain AH1-1. As expected, two hybridized signals were observed corresponding to the electrophoresis position of dsRNAs 1 and 2, further supporting the D-RNA nature of dsRNA2.
To further check whether PfHV1 dsRNAs are encapsidated, nucleic acid preparations of P. fici strain AH1-1 without protein-denaturing treatment process were subjected to RNase A treatment at 37°C or on ice, and P. theae strain LI41-1P1 infected by PtCV1 (an encapsidated dsRNA mycovirus) was involved in parallel as a control. The results showed that PfHV1 dsRNAs, similar to PtCV1, are resistant to RNase A treatment, suggesting that they are encapsulated (Fig 4C).
dsRNA3 is a satellite-like RNA (SatL-RNA) of PfHV1
The complete nucleotide sequence of dsRNA3 is 631 bp in length without a poly (A) tail, and it has a small putative ORF (nt positions 248–400) on its positive strand coding for a protein of 50 aa with an approximate molecular mass of 5.75 kDa. dsRNA3 contains nine putative IRES elements in nt positions 1–247 (e-values 2.8–31) as aligned with those deposited in IRESite database (S3 Table). Alignment of dsRNA3 with dsRNAs 1 and 2 indicated that their first 170 bp were 100% identical, excluding 3 nts deletion at 5′-termini of dsRNA 2, while the remaining sequence had no detectable similarity with both dsRNAs. To further confirm that dsRNA3 is identical with dsRNAs 1 and 2 at the 5′-termini while heterogenetic in the remaining, two DIG-labeled riboprobes (termed riboprobe 2 and 3) binding to the identical and heterogenetic regions of dsRNA1, respectively, were synthesized in vitro transcription and hybridization with the nucleic acid preparations of P. fici strain AH1-1. As riboprobe 2 was utilized, three hybridized signals were observed corresponding to the electrophoresis position of dsRNAs 1 to 3; while as riboprobe 3 was utilized, only one hybridization signal was observed corresponding to the electrophoresis position of dsRNA3 (Fig 4B). Thus, the hybridization analysis supports the conclusion from the sequence alignment.
To further confirm that dsNRA3 is a SatL-RNA hosted by PfHV1, the presumed viral vesicles were purified by centrifugation in stepwise sucrose gradients (10–35% with 5% sucrose increments), and subjected to nucleic acids extraction. Agarose gel electrophoresis of the resulted nucleic acids showed that dsRNA3 together with the typical pattern of PfHV1 dsRNAs 1 and 2 was most recovered from the 15% fractions (Fig 4D). These data suggest that dsRNAs 1 to 3 are encapsulated together, and it supports that dsNRA3 is a SatL-RNA of PfHV1.
dsRNA3 negatively regulates the replication of PfHV1
To investigate whether the dsRNA3 affects the replication of PfHV1, the full-length of dsRNA3 cDNAs was engineered into an expression vector containing a hygromycin resistance marker, and transformed into the protoplasts prepared from the dsRNA1-only infected strain AH1-1-24 of P. fici. With the resistance selection and qPCR identification, three derivative subisolates, termed AH1-1-24SatL-OE1 to AH1-1-24SatL-OE3, transformed by RNA3-expression vector were obtained. Expression analysis of dsRNA1 together with RNA3 using qRT-PCR in the three resulting subisolates revealed that the expression levels of RNA3 were significantly up-regulated up to 1.96 to 9.72 folds in these transformed subisolates (Fig 4E), whereas dsRNA1 were significantly down-regulated, up to 30.03% to 70.63% (Fig 4E). These results suggest that the overexpression of the satellite-like RNA inhibits the replication of PfHV1. Whereas, no obvious accumulation of SatL-RNA was observed in these over-expressed strains as checked by electrophoresis analysis of their dsRNAs on agarose gel (S2A Fig).
PfHV1 dsRNAs are transmitted vertically and horizontally in different efficiency
To investigate potential vertical transmission of PfHV1 in P. fici strain AH1-1, individual conidia were isolated from mycelia and cultured on PDA. A total of 51 single subisolates were randomly selected and analyzed for the presence of PfHV1 dsRNAs. The experiments showed all (100%) of the subisolates were infected with PfHV1 dsRNA1 as detected by agarose gel analysis of their nucleic acids after treated by DNase I (Fig 5A). This was further confirmed by dot blotting analysis of 25 randomly selected subisolates using riboprobe 1 and RT-PCR analysis of 8 subisolates using primers PfHV1-1F/1R for dsRNA 1, respectively (Fig 5B). It is worth noting that the RT-PCR products of PfHV1 dsRNA1 showed fewer contents for all the derived subisolates as compared with the parental strain AH1-1 (Fig 5B) whereas dsRNAs 2 or 3 were not detected in all these conidium-generated subisolates, suggesting that they are not vertically transmitted via conidia.
(A) Electrophoresis analysis of dsRNAs extracted from the conidium-generated subisolates of strain AH1-1 on 1% agarose gel after being treated by DNase I. (B) Dot blotting analysis with riboprobe 1 (I) and RT-PCR detection with primer pair PfHV1-1F/-1R (II, S1 Table) of PfHV1 in some randomly selected conidium-generated subisolates. (C) Electrophoresis analysis on 1% agarose gel of dsRNA-enriched fractions (without treatment by DNase I) and RT-PCR amplicons with primer pair PfHV1-1F/-1R of PfHV1 in the subisolates (AH1-1hygV+-1 to -4) of AH1-1hygV- (as a recipient strain) contact culture with AH1-1 (donor). Since the fungal genomic bands were kept in the same electrophoresis positions of dsRNA1, the PfHV1 presentence was detected by RT-PCR correspondingly. (Note: D-RNA and dsRNA3 turned faint for unknown reason but they were in the strains identified by both RT-PCR and dot blotting).
To investigate horizontal transmission of PfHV1, PfHV1-infected AH1-1 (donor) and AH1-1hygV- (receptor) were dual-cultured, and four AH1-1hygV- mycelium discs next to the contact area were randomly selected for RT-PCR analysis of the presence of PfHV1 dsRNAs. This experiment revealed that all the selected subisolates (nominated AH1-1hygV+-1 to -4) were infected by PfHV1 dsRNAs 1 and 2, while not by dsRNA3 (Fig 5C). Whereas as dual-culture of AH1-1 and PfHV1-free Pestalotiopsis sp. CJB-1, none of 12 isolates from the later were infected by PfHV1 dsRNAs 1 to 3.
These results suggest that PfHV1 dsRNA1 is efficiently transmitted vertically and horizontally in the host strain, while not easily to other Pestalotiopsis strains. DsRNA2 is efficiently transmitted horizontally but not vertically in the same host strain, whereas dsRNA3 could not be transmitted horizontally and vertically.
PfHV1 confers no obvious effects on the biological traits of the host fungus
To check whether PfHV1 confers some biological effects on the host fungus, the morphologies and growth rates were accessed for strain AH1-1 (PfHV1+) and its cured subisolate AH1-1V- (PfHV1-) after culture on PDA in parallel at 4 dpi, no obvious difference was observed in these aspects (Fig 6A and 6B). Moreover, they induced similar sized lesions as being inoculated on wounded tea leaves (C. sinensis var. Fuyun no.6) (Fig 6C). As compared with other strains free of PfHV1 including Pestalotiopsis sp. strains (AH1-1-14, TP-2-2W, JWX-3-1, and CJB-4-1), they also showed similar morphologies, growth rates, and virulence (Fig 6D and 6E). These results suggest that PfHV1 confers no obvious effects on the morphologies and virulence of P. fici.
(A and B) Colonies and growth rates of AH1-1V- (PfHV1-) and AH1-1 (PfHV1+) strains cultured on PDA medium for 6 days, respectively. The growth rates were measured at 2, 4, and 6 day, respectively. (C) The lesion lengths on tea leaves (Camellia sinensis var. Fuyun no.6) induced by inoculation with the mycelium discs of AH1-1V- and AH1-1 strains at 4 day post inoculation (dpi). (D) Colonies and growth rates of AH1-1 and other Pestalotiopsis spp. strains on PDA medium cultured for 5 days. (E) The symptoms and lesion lengths caused by AH1-1 and other Pestalotiopsis spp. strains on tea leaves (C. sinensis var. Fuyun no.6), respectively. Data are means ± SEM (n = 6 or 3 for virulence and growth rate assessment, respectively). Different letters indicate a significant difference at p < 0.05 (one-way ANOVA).
Discussion
In this study, three dsRNA components were detected in P. fici, and their full-length sequences were determined and characterized, suggesting that they belong to a novel hypovirus, tentatively named PfHV1, under the genus Alphahypovirus [3,4,30]. Here, the dsRNA nature of PfHV1 genome was confirmed with enzymatic treatments, suggesting that PfHV1 genomic components accumulate in dsRNA nature in vivo in abundant. The predominant dsRNA form found in infected mycelia appears to be replicative intermediate or replicative form RNA, while the accumulation of plus strand RNA is low due to the RNA silencing antiviral defense response of the fungal host [31]. Other than CHV 1 to CHV4 which have been well studied in Cryphonectria parasitica [1], hypoviruses or hypo-like viruses have been largely discovered from other phytopathogenic fungi, including Sclerotinia sclerotiorum [4,11,30,32], Valsa ceratosperma [3], Fusarium graminearum [8,12], Phomopsis longicolla [33], Macrophomina phaseolina [34], Botrytis cinerea [7], and Rosellinia necatrix [35], and even in non-fungal eukaryotes (invertebrates) [36] revealed by large-scale meta-transcriptomic analysis. To our knowledge, this is the first report of hypovirus isolated from Pestalotiopsis fungus, and the second mycovirus besides of a chrysovirus infecting this fungal genus.
Genomic analysis of PfHV1 dsRNA1 suggests that it contains one ORF coding a polyprotein, which is somewhat different from other members in “Alphahypovirus”, since the group genomes typically harbor two ORFs (A and B), with two papain-like cysteine Pro encoded by ORF A and the leading regions of ORF B, respectively, as exemplified by CHV1 [28,37,38], while it is similar to the genomic organization of AaHV1 and those of genera Betahypovirus and Gammahypovirus. Moreover, a complete papain-like cysteine Pro motif was not detected in the polyprotein of PfHV1, whereas alignment with the polyproteins of other hypoviruses uncovered the presence of two conserved cysteine Pro core residues (cysteine, and glycine in Fig 3B;) [22] but lacked the histidine at its N-terminal region (29–132 aa region), conserved among hypo- and hypo-like viruses [6], suggesting that the cysteine Pro might be encoded by PfHV1 dsRNA1 to process the polyprotein, but it possess an obvious diversity and evolution at this region of PfHV1 genome. Moreover, the conserved domains of ORFs coding for the functional proteins including RdRp and Hel were detected in the polyprotein of PfHV1. Of these, RdRp is important for RNA viral replication and always conserved in polyproteins encoded by hypo- and hypo-like viruses [6,29]. While several aa were divergent in the conserved motifs as compared with the known members, e.g., in RdRp motif IV and VI. More importantly, a conserved SDD tripeptide was observed in RdRp motif IV, and it is different from most other known (+) ssRNA viruses, in which the consensus tripeptide is GDD [39,40].
dsRNA 2 is determined as a defective component of dsRNA 1, which is very stable since it could be horizontally transmitted, as well as easily observed in other natural Pestalotiopsis spp. strains, like the defective RNA 2 associated with CHV3 [41]. However, dsRNA 2 could not be vertically transmitted together with dsRNA1, suggesting that is not a normal and necessary component of PfHV1 but rather a by-product of virus replication. It is worth noting that it is composed of two inversed sequences of considerately large size, suggesting that a deletion-junction process had happened after the transcription of dsRNA1, in which the defective RNAs is maintenance of an ORF through a deletion-junction to maintain the translation of two different genes [42,43]. The coding capacity of the defective RNA and ability to be translated may affect defective RNA accumulation [44,45], since we have observed the PfHV1 dsRNA1 turned into an obviously lower titer after vertical transmission in the fungal progeny subisolates without the presence of the dsRNA2 components. It is unknown what viral factors are involved in generation and maintenance of the defective dsRNA, which is likely cleaved by dicer-like 2 (DCL2) and argonaute-like 2 (AGL2) proteins responsible for the antiviral RNA silencing pathways in a filamentous fungus [10,46].
dsRNA3 segment is determined to be a SatL-RNA of PfHV1 due to the following reasons: 1) it was co-precipitated with dsRNAs 1 and 2 in the same sucrose fraction of 15% by ultra-centrifuge, suggesting that it is encapsulated together with PfHV1 genomic dsRNAs; 2) it codes no structural proteins or RdRp proteins, suggesting its replication depends on the helper virus; 3) no dsRNA3 component was solely observed in any subisolates absent of dsRNA1 after vertical or horizontal transmission. Until now, three satellite or SatL-RNAs have been observed in association with hypoviruses, i.e., dsRNA4 (937 bp in size) and dsRNA3 (dimeric dsRNA4 linking by a poly (A) of CHV3 in C. parasitica and S-dsRNA (3.6 kb) of SsHV1 in S. sclerotiorum, both of which encode authentic proteins and contain a poly (A) tail [11,41]. Of these, CHV3 dsRNA-4 was confirmed to encode a small polypeptide of 9.4 kDa with an in vitro translation experiment [41,47]; S-dsRNA (3.6 kb) of SsHV1 was proposed to encode a protein of 639 aa with an approximate molecular mass of 71 kDa [11]. PfHV1 dsRNA3 shares no detectable identities with the both satellite (or satellite-like) RNAs, contains no a poly (A) tail, and most likely encodes no proteins, suggesting that dsRNA3 is a completely different SatL-RNA as compared the ones of CHV3 and SsHV1, and belongs to a new class of satellite nucleic acids. Generally, satellites are distinct from their helper virus with a nucleotide sequence that is substantially different from that of their helper virus, generally having very little or no sequence similarity with helper viruses, encapsidated by the capsid protein (CP) of helper viruses [48]. Here, dsRNA3 should be encapsulated in the replication vesicles, composed of host-derived lipids [49], instead of CP by the helper virus PfHV1 since the latter has no structural proteins; moreover, dsRNA3 shares a substantially identical sequence (170 bp in size) with the helper virus genome, suggest that it partially originated from its helper virus and is distinct from a typical satellite nucleic acid. Therefore, PfHV1 dsRNA3 represents a novel class of satellite nucleic acids unreported before, unlike the dsRNA satellites in association with Totiviridae and Partitiviridae. We propose these satellite nucleic acids are classified into two classes including those coding dsRNAs with a poly (A) tail (like CHV3 dsRNAs 3 and 4) and those no-coding dsRNAs without a poly (A) tail (like PfHV1 dsRNA3). Additionally, we suggest to broaden the definition of fungal satellites to accommodate these satellite RNAs by including a special type of satellite nucleic acid that has substantial sequence homology with the host viral genome without encapsidation in a coat protein.
Satellite RNAs have been described in association with several mycoviruses in filamentous fungi and yeast [50–52], and it has been proved that the amplification used viral RdRp [13,53]. Of them, the satellite RNAs associated with Totiviridae encode a secreted preprotoxin that is lethal to sensitive cells (virus-free or containing helper virus only), and imparts self-protection against the secreted toxin and confers ecological advantage by killing competing virus- or satellite-free fungi; the satellite RNA associated with SsHV1 most likely enhances the hypovirulence traits of S. sclerotiorum. However, the biological functions of other satellite RNAs remain largely unknown. Here, as compared with the dsRNA3-infected and -free strains of P. fici isolates, i.e., strain AH1-1 subisolates derived from horizontal transmission and those from vertical transmission, no obvious morphologies, growth rates or virulence were observed, suggesting that the SatL-RNA has no impact on these aspects. However, as the satellite-like RNA was overexpressed, it inhibits the replication of PfHV1, suggesting that dsRNA3 negatively regulates the replication of PfHV1. These results suggest that the SatL-RNAs harbored by hypoviruses have a similar biological trait to the satellite RNAs associated with animal and plant viruses. E.g., amplification of satellites by viral RdRp may down-regulate synthesis of viral RNAs and expression of their products in the research of Trichomonas vaginalis virus [53]; cucumber mosaic virus (CMV) satellite RNA reduced the yield of accumulated helper virus [54,55]. Based on the in vitro competition assay, the reduction was proposed to be due to the competition of satellite RNA for limited amount of viral replicase with CMV gRNAs [56]. Furthermore, a three-way branched secondary structure was identified in the satellite RNA, and was indispensable for the helper virus inhibition [57]. Similarly, a conserved apical hairpin stem-loop structure was identified in the 5′-UTR of bamboo mosaic virus (BaMV) satellite RNA as a potent determinant of the down-regulation of helper virus replication [58]. Since the SatL-RNA encodes no functional proteins, we conclude that it inhibits the replication of helper virus most likely with similar mechanism to the CMV and BaMV satellite RNAs.
In summary, to gain insight into mycoviruses in P. fici, belonging to an important fungus genus related to agriculture, industry and medicine, a novel hypovirus from P. fici together with its D-RNA and a SatL-RNA, were identified and characterized. To our knowledge, PfHV1 represents the first hypovirus and the second mycovirus isolated from Pestalotiopsis spp., and its SatL-RNA represents a novel class of satellite nucleic acids, which should contribute to our better understanding of the taxonomy, evolution, molecular and biological traits of mycoviruses, as well as the related satellites.
Supporting information
S2 Table. Blastp searches of dsRNA1-coding proteins with those deposited in NCBI.
https://doi.org/10.1371/journal.ppat.1010889.s002
(XLSX)
S3 Table. Predicted IRES elements in SatL-RNA UTR.
https://doi.org/10.1371/journal.ppat.1010889.s003
(XLSX)
S1 Fig. Predicted secondary structures of both UTR of PfHV1.
(A and B) The secondary structures of 5′- and 3′-UTR of PfHV1 were predicted in RNAfold, respectively. (C) Six AUGs in the 5′-UTR of PfHV1 were found.
https://doi.org/10.1371/journal.ppat.1010889.s004
(TIF)
S2 Fig. Electrophoresis of dsRNAs extracted from strains AH1-1-24 and its subisolates with SatL-RNA being overpressed.
https://doi.org/10.1371/journal.ppat.1010889.s005
(TIF)
S3 Fig. The prediction of secondary structure of SatL-RNA.
Two black curves refer to the nucleotides that make up a potential pseudoknot.
https://doi.org/10.1371/journal.ppat.1010889.s006
(TIF)
References
- 1. Suzuki N, Ghabrial SA, Kim K-H, Pearson M, Marzano S-YL, Yaegashi H, et al. ICTV Virus Taxonomy Profile: Hypoviridae. J Gen Virol. 2018;99: 615–616. pmid:29589826
- 2. Hillman BI, Suzuki N. Viruses of the Chestnut Blight Fungus, Cryphonectria parasitica. Adv Virus Res. 2004;63: 423–472. pmid:15530566
- 3. Yaegashi H, Kanematsu S, Ito T. Molecular characterization of a new hypovirus infecting a phytopathogenic fungus, Valsa ceratosperma. Virus Res. 2012;165: 143–150. pmid:22366520
- 4. Khalifa ME, Pearson MN. Characterisation of a novel hypovirus from Sclerotinia sclerotiorum potentially representing a new genus within the Hypoviridae. Virology. 2014;464–465: 441–449. pmid:25108682
- 5. Turina M, Chiba S, Velasco L, Ayllon M, Suzuki N, Lee-Marzano S, et al. Create 35 new species and eight new genera, abolish one genus (Hypovirus) and rename four existing species (Durnavirales: Hypoviridae). ICTV. 2020;2021.006F. Available from: https://ictv.global/files/proposals/approved?fid=4894#block-teamplus-page-title
- 6. Li H, Bian RL, Liu Q, Yang L, Pang TX, Salaipeth L, et al. Identification of a novel hypovirulence-inducing hypovirus from Alternaria alternata. Front Microbiol. 2019;10: 1076. pmid:31156589
- 7. Hao FM, Ding T, Wu MD, Zhang J, Yang L, Chen WD, et al. Two Novel Hypovirulence-Associated Mycoviruses in the Phytopathogenic Fungus Botrytis cinerea: Molecular Characterization and Suppression of Infection Cushion Formation. Viruses. 2018;10: E254. pmid:29757259
- 8. Wang SC, Kondo H, Liu L, Guo LH, Qiu DW. A novel virus in the family Hypoviridae from the plant pathogenic fungus Fusarium graminearum. Virus Res. 2013;174: 69–77. pmid:23499998
- 9. Shapira R, Choi GH, Hillman BI, Nuss DL. The contribution of defective RNAs to the complexity of viral-encoded double-stranded RNA populations present in hypovirulent strains of the chestnut blight fungus Cryphonectria parasitica. EMBO J. 1991;10: 741–746. pmid:2009855
- 10. Zhang XM, Nuss DL. A host dicer is required for defective viral RNA production and recombinant virus vector RNA instability for a positive sense RNA virus. Proc Natl Acad Sci U S A. 2008;105: 16749–16754. pmid:18922782
- 11. Xie JT, Xiao XQ, Fu YP, Liu HQ, Cheng JS, Ghabrial SA, et al. A novel mycovirus closely related to hypoviruses that infects the plant pathogenic fungus Sclerotinia sclerotiorum. Virology. 2011;418: 49–56. pmid:21813149
- 12. Li PF, Zhang HL, Chen XG, Qiu DW, Guo LH. Molecular characterization of a novel hypovirus from the plant pathogenic fungus Fusarium graminearum. Virology. 2015;481: 151–160. pmid:25781585
- 13. Gnanasekaran P, Chakraborty S. Biology of viral satellites and their role in pathogenesis. Curr Opin Virol. 2018;33: 96–105. pmid:30144641
- 14.
Briddon RW, Ghabrial S, Lin NS, Palukaitis P, Scholthof K-BG, Vetten HJ. The subviral agents. In: Andrew MQK, Michael JA, Eric BC, Elliot JL, editors. Virus taxonomy: Ninth report of the International Committee on Taxonomy of Viruses. Boston: Academic Press; 2012. pp. 1209–1219.
- 15. Zhou LL, Li XP, Kotta-Loizou I, Dong KL, Li SF, Ni DJ, et al. A mycovirus modulates the endophytic and pathogenic traits of a plant associated fungus. ISME J. 2021;15: 1893–1906. pmid:33531623
- 16. Wang LP, He Y, Kang YP, Hong N, Farooq ABU, Wang GP, et al. Virulence determination and molecular features of peach latent mosaic viroid isolates derived from phenotypically different peach leaves: A nucleotide polymorphism in L11 contributes to symptom alteration. Virus Res. 2013;177: 171–178. pmid:23973915
- 17. Maharachchikumbura SSN, Hyde KD, Groenewald JZ, Xu J, Crous PW. Pestalotiopsis revisited. Stud Mycol. 2014;79: 121–186. pmid:25492988
- 18. Nuss DL. Hypovirulence: Mycoviruses at the fungal–plant interface. Nat Rev Microbiol. 2005;3: 632–642. pmid:16064055
- 19. Wang ZH, Zhao ZX, Hong N, Ni DJ, Cai L, Xu WX, et al. Characterization of Causal Agents of a Novel Disease Inducing Brown-Black Spots on Tender Tea Leaves in China. Plant Dis. 2017;101: 1802–1811. pmid:30676920
- 20. Jia HX, Dong KL, Zhou LL, Wang GP, Hong N, Jiang DH, et al. A dsRNA virus with filamentous viral particles. Nat Commun. 2017;8: 168. pmid:28761042
- 21. Zhai LF, Xiang J, Zhang MX, Fu M, Yang ZK, Hong N, et al. Characterization of a novel double-stranded RNA mycovirus conferring hypovirulence from the phytopathogenic fungus Botryosphaeria dothidea. Virology. 2016;493: 75–85. pmid:27015523
- 22. Liu HQ, Fu YP, Jiang DH, Li GQ, Xie J, Peng YL, et al. A novel mycovirus that is related to the human pathogen hepatitis E virus and rubi-like viruses. J Virol. 2009;83: 1981–1991. pmid:19073734
- 23. Tamura K, Stecher G, Kumar S. MEGA11: Molecular evolutionary genetics analysis version 11. Mol Biol Evol. 2021;38: 3022–3027. pmid:33892491
- 24. Chiba S, Jamal A, Suzuki N. First evidence for internal ribosomal entry sites in diverse fungal virus genomes. mBio. 2018;9: e02350–17. pmid:29559577
- 25.
Pallás V, Sánchez-Navarro JA, Kinard GR, Di Serio F. Chapter 35—Molecular hybridization techniques for detecting and studying viroids. In: Hadidi A, Flores R, Randles JW, Palukaitis P, editors. Viroids and Satellites. Boston: Academic Press; 2017. pp. 369–379. https://doi.org/10.1016/B978-0-12-801498-1.00035–8
- 26. Sato Y, Shamsi W, Jamal A, Bhatti MF, Kondo H, Suzuki N. Hadaka virus 1: a capsidless eleven-segmented positive-sense single-stranded RNA virus from a phytopathogenic fungus, Fusarium oxysporum. mBio. 2020;11: e00450–20. pmid:32457242
- 27. Sesma A, Osbourn AE. The rice leaf blast pathogen undergoes developmental processes typical of root-infecting fungi. Nature. 2004;431: 582–586. pmid:15457264
- 28. Koonin EV, Choi GH, Nuss DL, Shapira R, Carrington JC. Evidence for common ancestry of a chestnut blight hypovirulence-associated double-stranded RNA and a group of positive-strand RNA plant viruses. Proc Natl Acad Sci U S A. 1991;88: 10647–10651. pmid:1961731
- 29. Yao ZT, Zou CW, Peng N, Zhu Y, Bao YX, Zhou QJ, et al. Virome identification and characterization of Fusarium sacchari and F. andiyazi: Causative agents of pokkah boeng disease in sugarcane. Front Microbiol. 2020;11: 240. pmid:32140150
- 30. Hu ZJ, Wu SS, Cheng JS, Fu YP, Jiang DH, Xie JT. Molecular characterization of two positive-strand RNA viruses co-infecting a hypovirulent strain of Sclerotinia sclerotiorum. Virology. 2014;464–465: 450–459. pmid:25104554
- 31.
Dawe AL, Nuss DL. Chapter five—Hypovirus molecular biology: From Koch’s postulates to host self-recognition genes that restrict virus transmission. In: Ghabrial SA, editor. Advances in Virus Research. Boston: Academic Press; 2013. pp. 109–147. https://doi.org/10.1016/B978-0-12-394315-6.00005–2
- 32. Marzano SYL, Hobbs HA, Nelson BD, Hartman GL, Eastburn DM, McCoppin NK, et al. Transfection of Sclerotinia sclerotiorum with in vitro transcripts of a naturally occurring interspecific recombinant of Sclerotinia sclerotiorum hypovirus 2 significantly reduces virulence of the fungus. J Virol. 2015;89: 5060–5071. pmid:25694604
- 33. Koloniuk I, El-Habbak MH, Petrzik K, Ghabrial SA. Complete genome sequence of a novel hypovirus infecting Phomopsis longicolla. Arch Virol. 2014;159: 1861–1863. pmid:24473711
- 34. Marzano SYL, Domier LL. Novel mycoviruses discovered from metatranscriptomics survey of soybean phyllosphere phytobiomes. Virus Res. 2016;213: 332–342. pmid:26547008
- 35. Arjona-Lopez JM, Telengech P, Jamal A, Hisano S, Kondo H, Yelin MD, et al. Novel, diverse RNA viruses from Mediterranean isolates of the phytopathogenic fungus, Rosellinia necatrix: Insights into evolutionary biology of fungal viruses. Environ Microbiol. 2018;20: 1464–1483. pmid:29411500
- 36. Shi M, Lin XD, Tian JH, Chen LJ, Chen X, Li CX, et al. Redefining the invertebrate RNA virosphere. Nature. 2016;540: 539–543. pmid:27880757
- 37. Choi GH, Pawlyk DM, Nuss DL. The autocatalytic protease p29 encoded by a hypovirulence-associated virus of the chestnut blight fungus resembles the potyvirus-encoded protease HC-Pro. Virology. 1991;183: 747–752. pmid:1853573
- 38. Shapira R, Nuss DL. Gene expression by a hypovirulence-associated virus of the chestnut blight fungus involves two papain-like protease activities. Essential residues and cleavage site requirements for p48 autoproteolysis. J Biol Chem. 1991;266: 19419–19425. pmid:1918054
- 39. Poch O, Blumberg BM, Bougueleret L, Tordo N. Sequence comparison of five polymerases (L proteins) of unsegmented negative-strand RNA viruses: theoretical assignment of functional domains. J Gen Virol. 1990;71 (Pt5): 1153–1162. pmid:2161049
- 40. Xu X, Liu YQ, Weiss S, Arnold E, Sarafianos SG, Ding JP. Molecular model of SARS coronavirus polymerase: implications for biochemical functions and drug design. Nucleic Acids Res. 2003;31: 7117–7130. pmid:14654687
- 41. Hillman BI, Foglia R, Yuan W. Satellite and defective RNAs of Cryphonectria hypovirus 3-Grand Haven 2, a virus species in the family Hypoviridae with a single open reading frame. Virology. 2000;276: 181–189. pmid:11022006
- 42. Romero J, Huang Q, Pogany J, Bujarski JJ. Characterization of defective interfering RNA components that increase symptom severity of broad bean mottle virus infections. Virology. 1993;194: 576–584. pmid:8503176
- 43. White KA, Bancroft JB, Mackie GA. Defective RNAs of clover yellow mosaic virus encode nonstructural/coat protein fusion products. Virology. 1991;183: 479–486. pmid:1830181
- 44. Pogany J, Romero J, Huang Q, Sgro JY, Shang H, Bujarski JJ. De novo generation of defective interfering-like RNAs in broad bean mottle bromovirus. Virology. 1995;212: 574–586. pmid:7571427
- 45. White KA, Bancroft JB, Mackie GA. Coding capacity determines in vivo accumulation of a defective RNA of clover yellow mosaic virus. J Virol. 1992;66: 3069–3076. pmid:1560537
- 46. Chiba S, Suzuki N. Highly activated RNA silencing via strong induction of dicer by one virus can interfere with the replication of an unrelated virus. Proc Natl Acad Sci U S A. 2015;112: E4911–4918. pmid:26283371
- 47. Yuan W, Hillman BI. In vitro translational analysis of genomic, defective, and satellite RNAs of Cryphonectria hypovirus 3-GH2. Virology. 2001;281: 117–123. pmid:11222102
- 48.
Maloy OC, Murray TD. Encyclopedia of Plant Pathology. 2001, Wiley J. ed. Hoboken: John Wiley & Sons; 2001.
- 49. Newhouse JR, MacDonald WL, Hoch HC. Virus-like particles in hyphae and conidia of European hypovirulent (dsRNA-containing) strains of Cryphonectria parasitica. Can J Bot. 1990;68: 90–101.
- 50. Deng F, Boland GJ. A satellite RNA of Ophiostoma novo-ulmi mitovirus 3a in hypovirulent isolates of Sclerotinia homoeocarpa. Phytopathology. 2004;94: 917–923. pmid:18943067
- 51. Wickner RB. Double-stranded RNA viruses of Saccharomyces cerevisiae. Microbiol Rev. 1996;60: 250–265. pmid:8852903
- 52. Romanos MA, Buck KW, Rawlinson CJ. A satellite double-stranded RNA in a virus from Gaeumannomyces graminis. J Gen Virol. 1981;57(DEC): 375–385.
- 53. Khoshnan A, Alderete JF. Characterization of double-stranded RNA satellites associated with the Trichomonas vaginalis virus. J Virol. 1995;69: 6892–6897. pmid:7474105
- 54. Gal-On A, Kaplan I, Palukaitis P. Differential effects of satellite RNA on the accumulation of cucumber mosaic virus RNAs and their encoded proteins in tobacco vs zucchini squash with two strains of CMV helper virus. Virology. 1995;208: 58–66. pmid:11831731
- 55. Liao QS, Zhu LP, Du ZY, Zeng R, Peng JL, Chen JS. Satellite RNA-mediated reduction of cucumber mosaic virus genomic RNAs accumulation in Nicotiana tabacum. Acta Biochim Biophys Sin (Shanghai). 2007;39: 217–223. pmid:17342261
- 56. Wu G, Kaper JM. Competition of viral and satellite RNAs of cucumber mosaic virus for replication in vitro by viral RNA-dependent RNA polymerase. Res Virol. 1995;146: 61–67. pmid:7538689
- 57. He L, Wang Q, Gu ZH, Liao QS, Palukaitis P, Du ZY. A conserved RNA structure is essential for a satellite RNA-mediated inhibition of helper virus accumulation. Nucleic Acids Res. 2019;47: 8255–8271. pmid:31269212
- 58. Hsu YH, Chen HC, Cheng J, Annamali P, Lin BY, Wu CT, et al. Crucial role of the 5′ conserved structure of bamboo mosaic virus satellite RNA in downregulation of helper viral RNA replication. J Virol. 2006;80: 2566–2574. pmid:16474162