Molecular characterization of a novel mycovirus in the cultivated mushroom, Lentinula edodes
Virology Journal volume 9, Article number: 60 (2012)
In the 1970s, mycoviruses were identified that infected the edible mushroom Lentinula edodes (shiitake), but they were not regarded as causal agents for mushroom diseases. None of their genes has been sequenced. In this study, the dsRNA genome of a mycovirus recently found in a shiitake commercial strain was sequenced and its molecular structure was characterized.
A cDNA library was constructed from a dsRNA purified from the fruiting body of L. edodes. The virus was tentatively named L. edodes mycovirus HKB (LeV). Based on the deduced RNA-dependent RNA polymerase (RdRp) sequence, phylogenetic analysis of LeV was conducted. Because no virion particles associated with the dsRNA were observed by electron microscopic observation, atomic force microscopy (AFM) observation was chosen for achieving molecular imaging of the virus.
The 11,282-bp genome of LeV was obtained. The genome encoded two open reading frames (ORFs). ORF1 coded for a hypothetical protein and ORF2 for a putative RdRp, respectively. In addition, a region coding for a NUDIX domain was present in ORF1. There was a 62-bp intergenic region between ORF1 and RdRp. Similarity with coat protein of mycoviruses was not found within the whole sequence. Based on phylogenetic analysis of the putative RdRp sequence, LeV grouped into a clade with dsRNA found in the basidiomycetes Phlebiopsis gigantea and Helicobasidium mompa. The clade was placed apart from the Totiviridae and Chrysoviridae families. As suggested from the genome sequence, AFM revealed that the structure of LeV was linear unencapsidated dsRNA.
The results suggest that LeV represents a novel family of mycoviruses, found thus far only among the basidiomycetes.
In the 1970s, viruses that infect the cultivated mushroom Lentinula edodes, or shiitake, were extensively studied in Japan [1–3], and three morphologically distinct viruses were detected by electron microscopy [1, 3]. However, unlike La France disease of the white button mushroom [4, 5], mycoviruses have not been associated with shiitake diseases because these mycoviruses have commonly been found in healthy fruiting bodies [1, 3]. In the USA, dsRNAs have also been observed in shiitake strains, but these appeared to be latent . In the 1970s, shiitake cultivation was performed by inoculating mycelium spawn on oak logs; however, this labor-intensive method was gradually replaced by indoor cultivation using sawdust substrate supplemented with rice bran. Currently in Japan, about 75,016 t (82% of all shiitake) are produced indoors annually using bag cultures with a sawdust-based substrate . In a bag culture, the shiitake mycelium is fully grown in the substrate until brown pigment is produced outside and the substrate becomes stiff. Complete browning of the exterior surface of the substrate is an important marker that normal fruiting bodies will develop in the following stage of cultivation. However, abnormal symptoms are occasionally observed in bag cultures, such as the growth of white or fluffy mycelia on the surface of substrate, inadequate or imperfect substrate browning , and malformations of the fruiting body. These symptoms often result in serious economic losses. Whether or not these abnormalities are associated with mycovirus is unknown. As a first step toward addressing that question, 46 shiitake isolates belonging to 11 commercial strains were examined for mycovirus infections . As a result, dsRNA was found in two isolates; one showed imperfect browning and the other was asymptomatic. Agarose gel analysis showed that the isolate with imperfect browning contained several dsRNAs, but the asymptomatic isolate contained only a single dsRNA. In this study, the single dsRNA was tentatively designated as Lentinula edodes mycovirus HKB (LeV) and was sequenced.
Three fruiting bodies were disrupted in 60 ml of 0.1 M phosphate buffer, and the virus fraction was precipitated with 10% PEG 8000 and 0.15 M NaCl, as described previously . The PEG precipitate was suspended in TES (10 mM Tris-HCl, 1 mM EDTA, 0.15 M NaCl, pH 7.0), and total RNA was isolated using the QIAmp Viral RNA Mini Kit (Qiagen) according to the manufacturer's instructions. Then dsRNA was isolated from the viral RNA by DNase I (Promega) digestion for 30 min at 37°C, followed by S1 nuclease (TaKaRa) digestion. The resulting dsRNA was concentrated in nuclease-free water by filtration through Ultrafree 0.5 100 K centrifugal filters (Millipore) several times to remove degraded nucleic acids and salts. Finally, dsRNA was purified with the RNeasy Mini Elute Cleanup Kit (Qiagen).
cDNA library construction and sequencing of dsRNA
The purified dsRNA served as a template for cDNA synthesis by random priming with the PrimeScript 1st strand cDNA Synthesis Kit (TaKaRa) according to the standard protocol, except that the denaturing condition was changed from 65°C, 5 min to 98°C, 1 min. The resulting cDNAs were electrophoresed in an agarose gel, and cDNAs sized 1.5-2.0 kb were extracted. The cDNAs were blunt-ended with T4 DNA polymerase, ligated into pUC118/HincII/BAP (TaKaRa), and transformed into E. coli DH10B cells by electroporation (Gene Pulser, Bio-Rad). The resulting cloned DNA was sequenced with the BigDye Terminator v.3.1 Cycle Sequencing Kit (Applied Biosystems) with M13 forward and reverse primers. A contiguous region (contig) was assembled with Sequencher™ 4.6 (Gene Codes Corp).
AFM microscopic observation
The PEG precipitate was suspended in 500 μl of TE and filtered through a 0.2-μm filter (Millipore), and 1 μl of the filtrate was diluted with 50 μl TE + 10 mM MgCl2. The sample was observed using atomic force microscopy (AFM) as previously described . A total of 10 μl of the RNA solution was dropped onto freshly cleaved muscovite mica (1 × 1 cm), which after standing for several minutes was washed with distilled water. The sample was dried under a stream of nitrogen. Observations were performed on a MFP-3D (Asylum Research) in the tapping mode in air, using a silicon cantilever OMCL-AC240TS (Olympus). Fields of 2 μm × 2 μm were scanned at a frequency of < 2 Hz. To confirm that the AFM image observed was RNA, AFM imaging was also performed with PEG precipitate after RNaseA digestion. The length of the RNA molecule was measured using software developed by the Research Institute of Biomolecule Metrology Co., Ltd. (Japanese patent P2000-230823A).
An unrooted polygenetic tree was constructed with sequences retrieved by a PSI-BLAST  search (Table 1) using the neighbor-joining method with a bootstrap of 1000. The sequence identity was below 40%; thus, we used the Multiple Alignment using Fast Fourier Transform (MAFFT) program L9INS-I http://mafft.cbrc.jp/alignment/server/ for making the multiple alignments and constructing the phylogenetic tree. Sequences with bootstrap values above 70% were included in the tree and visualized on the web http://www.genome.jp/tools/mafft/
Results and discussion
A cDNA library was constructed from dsRNA purified from the fruiting body of Lentinula edodes strain HKB. An 11,282-bp cDNA contig containing two open reading frames (ORFs) was obtained (deposited under Accession No.AB429556) (Figure 1). ORF1 encoded a 218,428-Da protein composed of 1,975 amino acids, containing the conserved motif of a NUDIX domain (between 319aa and 451aa) (pfam00293, NUDIX, 1.56e-04) . Viral coding NUDIX domain has been found only in Poxviruses . LeV, as well as Phlebiopsis gigantea mycovirus dsRNA 1 is the first mycovirus described, that codes for NUDIX domain. ORF2 encoded a 162,240-Da protein composed of 1426 amino acids, containing the conserved motif of RNA-dependent RNA polymerase (RdRp) (pfam02123, RDRP_4, 4.48e-22). The deduced amino acid sequence of ORF1 showed significant similarity (3e-83) to hypothetical protein PgV-1_gp1 of Phlebiopsis gigantea mycovirus dsRNA 1 . The deduced amino acid sequence of ORF2 showed similarity with members of Totivirus  and Chrysovirus . The family Totiviridae includes all the viruses of fungi and protozoa that have virions and a single-component dsRNA genome. The family Chrysovirus includes dsRNA viruses that infect fungi or plants and have four genome components in isometric virions. Because LeV contained a single component dsRNA genome, it could be classified as a member of Totivirus . But unlike Totivirus, there was no coding region for a coat protein within the LeV genome.
AFM microscopic observation
No virion or vesicle particles were detected by electron microscopy observation of the PEG precipitate of LeV. Purification of the virus by sucrose density-gradient ultracentrifugation was unsuccessful. When the PEG precipitate was directly digested with RNase, the viral dsRNA could no longer be detected by agarose gel electrophoresis (data not shown). These experimental data suggested that the virus did not form a virion and therefore was vulnerable to excess purification. Because the method is nondestructive, AFM was chosen for further observation . The PEG precipitate was suspended in TE and directly observed under AFM after being filtered through a 0.2-mm filter (Millipore). Numerous linear particles were observed (Figure 2A). The structure (Figure 2B) was very similar to the AFM images of dsRNA observed by Abels et al. and Vilfan et al. [26, 27]. To confirm that these structures were RNA, PEG precipitate after RNase digestion was also observed by AFM (Figure 3B). After the RNase digestion, the linear materials were no longer present. Thus, LeV was revealed as a linear dsRNA. The length of the dsRNA was 3,539 nm (Figure 3A), and its molecular weight was estimated (based on 1 μm = ~3 kb dsRNA ) to be 11.8 kb. The estimated size was very close to that of the sequenced dsRNA (11,282 bp).
In the current classification of viruses, unencapsidated mycoviruses are placed in either the Hypoviridae [28, 29] or Endornaviridae [30, 31] families. However, the putative RdRp sequence of LeV showed no similarity to the RdRp sequences of Hypovirus or Endornavirus. An unrooted polygenetic tree was constructed with sequences retrieved by a PSI-BLAST  search (Table 1) using the neighbor-joining method with a bootstrap of 1000. Sequences with bootstrap values above 70% were included in the tree (Figure 4). This analysis grouped LeV with Phlebiopsis gigantea mycovirus dsRNA1 (E value; 4e-150) (CAJ34333) and Helicobasidium mompa V670 L2-dsRNA virus (translated in this study from AB275288) (E value; 7e-33) in a clade independent of Totiviridae and Chrysoviridae (Figure 4). Both P. gigantea and H. mompa are basidiomycetes. No virion particles associated with P. gigantea mycovirus dsRNA have been described . Currently, whether H. mompa dsRNA forms virions or not is unknown. If H. mompa dsRNA was also unencapsidated, this would suggest that there is a novel family of mycoviruses so far found only in the basidiomycetes, having a monopartite dsRNA genome but do not package into virions.
Hypothesis of unencapsidated virus in the basidiomycetes
L. edodes is a white-rot basidiomycete and P. gigantea is also a basidiomycete that causes white rot in conifer logs and stumps [32, 33]. Often, P. gigantea is isolated from bark beetle, as is Lentinula boryana, a fungus belonging to the same family as L. edodes . Virion structure is necessary for a virus to infect and exit from the host cell and gain a greater probability of propagation. Because bark beetles feed on basidiomycete fungi , the evolution of unencapsidated virus in white-rot basidiomycetes may have been achieved through their association with wood-feeding insects. Because their host can be efficiently transferred to new environments by bark beetles, the coat protein would no longer be necessary for virus propagation. Additionally, the conserved amino acids of putative RdRp molecules of unencapsidated viruses isolated from plant-feeding insects, such as the alfalfa hopper and beet leafhopper , are significantly similar to LeV (Figure 4) (Table 1). This fact also supports the hypothesis that unencapsidated species of mycovirus present in the basidiomycetes might be closely associated with insects.
dsRNA found in a commercial strain of Lentinula edodes (shiitake) (designated as L. edodes mycovirus HKB; LeV) was sequenced. The 11.8-kb genome contained two ORFs. ORF1 coded for a hypothetical protein containing a NUDIX domain. Previously, viral NUDIX domain was found only in Poxviruses. LeV is the only mycovirus that codes NUDIX domain. ORF2 coded for RdRp with high similarity with Totivirus and Chrysovirus. The genome did not code for a coat protein, and AFM observation revealed LeV to be unencapsidated linear dsRNA. The results show that LeV is a novel mycovirus with a monopartite dsRNA genome that is closely related to Totivirus, but that it does not form a virion particle and might represent a new class of mycovirus.
Ushiyama R, Hashioka Y: Viruses associated with Hymenomycetes I Filamentous virus-like particles in the cells of a fruit body of shiitake, Lentinus edode (Berk) SING. Rept Tottori Mycol Inst. 1973, 10: 797-805.
Ushiyama N, Nakai Y, Ikegami M: Evidence for double-stranded RNA from polyhedral virus-like particles in Lentinus edodes (BERK) Sing. Virology. 1977, 77: 880-883. 10.1016/0042-6822(77)90512-8.
Mori K, Kuida K, Hosokawa D, Takehara M: Virus-like particle in several mushrooms. Mushroom Sci: Proceedings of the Tenth Internationl Congress on the Science and Cultivation of Edible Fungi. 1978, 10: 773-787.
Sinden JW, Hauser E: Report of two new mushroom diseases. Mushroom Sci. 1950, 1: 96-
Fletcher JT, Gaze RH: La France disease. In Mushroom Pest and Disease Control. Edited by: Fletcher JT, Gaze RH. 2008, London: Manson Publishing Ltd, 114-118.
Rytter JL, Royse DJ, Romaine CP: Incidence and diversity of double-stranded RNA in Lentinula edodes. Mycologia. 1991, 83: 506-510. 10.2307/3760362.
Statics of Forestry Agency, Ministry of Agriculture, Forestry and Fisheries. (in Japanese), [http://www.rinya.maff.go.jp/j/tokuyou/tokusan/1.html]
Ohta C, Taguchi T, Takahashi S, Ohtsuka K, Eda K, Ayusawa S, Magae Y: Detection of double stranded RNA elements in cultivated Lentinula edodes (in Japanese). Mushroom Sci Biotechnol. 2008, 16: 155-158.
Magae Y, Sunagawa M: Characterization of a mycovirus associated with the brown discoloration of edible mushroom, Flammulina velutipes. Virol J. 2010, 7: 342-
Umemura K, Nogami F, Okada T, Kuroda R: AFM characterization of single strand-specific endonuclease activity on linear DNA. Nucleic Acid Res. 2000, 28: E39-10.1093/nar/28.9.e39.
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acid Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.
Katoh K, Kuma K, Toh H, Miyata M: MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic Acid Res. 2005, 33: 511-518. 10.1093/nar/gki198.
Kozlakidis Z, Hacker CV, Bradley D, Jamal A, Phoon X, Webber J, Brasier CM, Buck KW, Coutts RH: Molecular characterisation of two novel double-stranded RNA elements from Phlebiopsis gigantea. Virus Genes. 2009, 39: 132-136. 10.1007/s11262-009-0364-z.
Kim JM, Kim JA, Park JA, Park SM, Cha BJ, Yang MS, Kim DH: Molecular diversity of chrysoviruses in Korean isolates of a new fungal species, Cryphonectria nitschkei. J Microbiol. 2009, 47: 441-447. 10.1007/s12275-009-0206-7.
Caston JR, Ghabrial SA, Jiang D, Rivas G, Alfonso C, Roca R, Luque D, Carrascosa JL: Three-dimensional structure of Penicillium chrysogenus virus: a double-stranded RNA virus with a genuine T = 1 capsid. J Mol Biol. 2003, 331: 417-431. 10.1016/S0022-2836(03)00695-8.
Rwahnih AM, Daubert S, Urbez-Torres JR, Cordero F, Rowhani A: Deep sequencing evidence from single grapevine plants reveals a virome dominated by mycoviruses. Arch Virol. 2011, 156: 397-403. 10.1007/s00705-010-0869-8.
Jamal A, Bignell EM, Coutts RH: Complete nucleotide sequences of four dsRNAs associated with a new chrysovirus infecting Aspergillus fumigatus. Virus Res. 2010, 153: 64-70. 10.1016/j.virusres.2010.07.008.
Spear A, Sisterson MS, Yokomi R, Stenger DC: Plant-feeding insects harbor double-stranded RNA viruses encoding a novel proline-alanine rich protein and a polymerase distantly related to that of fungal viruses. Virology. 2010, 404: 304-311. 10.1016/j.virol.2010.05.015.
Icho T, Wickner RB: The double-stranded RNA genome of yeast virus L-A encodes its own putative RNA polymerase by fusing two open reading frames. J Biol Chem. 1989, 264: 6716-6723.
Preisig O, Wingfield BD, Wingfield MJ: Coinfection of a fungal pathogen by two distinct double-stranded RNA viruses. Virology. 1998, 252: 399-406. 10.1006/viro.1998.9480.
Parrish S, Resch W, Moss B: Vaccinia virus D10 protein has mRNA decapping activity, providing a mechanism for control of host and viral gene expression. Proc Natl Acad Sci USA. 2007, 104: 2139-2144.
Parrish S, Moss B: Characterization of a second vaccinia virus mRNA-decapping enzyme conserved in poxviruses. J Virol. 2007, 81: 12973-12978. 10.1128/JVI.01668-07.
Wickner RB, Wang CC: Patterson JL: Totiviridae. In Virus Taxonomy 8th Report of the ICTV. Edited by: Fauquet CM, Mayo MA, Maniloff J, Desselberger U, Ball LA. 2005, Amsterdam: Elsevier, 571-580.
Ghabrial SA, Buck KW, Hillman BI, Martelli GP, Milne RG: Chrysoviridae. In Virus Taxonomy 8th Report of the ICTV. Edited by: Fauquet CM, Mayo MA, Maniloff J, Desselberger U, Ball LA. 2005, Amsterdam: Elsevier, 581-590.
Kuznetsov YG, McPherson A: Atomic force microscopy in imaging of viruses and virus-infected cells. Microbiol Mol Biol Rev. 2011, 75: 268-285. 10.1128/MMBR.00041-10.
Abels JA, Moreno-Herrero F, van der Heijden T, Dekker C, Dekker NH: Single-molecule measurements of the persistence length of double-stranded RNA. Biophys J. 2005, 88: 2737-2744. 10.1529/biophysj.104.052811.
Vilfan ID, Kamping W, van den Hout M, Candelli A: An RNA toolbox for single-molecule force spectroscopy studies. Nucleic Acid Res. 2007, 35: 6625-6639. 10.1093/nar/gkm585.
Nuss DL: Hypovirulence: mycoviruses at the fungal-plant interface. Nat Rev Microbiol. 2005, 8: 632-642.
Nuss DL, Hillman BI, Rigling D, Suzuki N: Hypoviridae. In Virus Taxonomy 8th Report of the ICTV. Edited by: Fauquet CM, Mayo MA, Maniloff J, Desselberger U, Ball LA. 2005, Amsterdam: Elsevier, 597-601.
Gibbs MJ, Pfeiffer P, Fukuhara T: Endornavirus. In Virus Taxonomy 8th Report of the ICTV. Edited by: Fauquet CM, Mayo MA, Maniloff J, Desselberger U, Ball LA. 2005, 605: 603-
Fukuhara T, Moriyama H: Endornaviruses. In Encyclopedia of Virology. Edited by: Mahy BWJ, Van Regenmortel MHV. 2008, Oxford: Elsevier, 109-116.
Behrendt CJ, Blanchette RA: Biological processing of pine logs for pulp and paper production with Phlebiopsis gigantea. Appl Environ Microbiol. 1997, 63: 1995-2000.
Cram M: Phlebiopsis (=Phanerochaete, =Peniophora, =Phlebia, =Peniophora) gigantea (Basidiomycetes: Corticiaceae). Biological Control: A Guide to Natural Enemies in North America. Edited by: Shelton AM. 2007, [http://www.biocontrol.entomology.cornell.edu/pathogens/phlebiopsis.html]
Hsiau PTW, Harrington TC: Phylogenetics and adaptations of basidiomycetous fungi fed upon by bark beetles (Coleoptera; Scolytidae). Symbiosis. 2003, 34: 111-131.
I would like to thank China Ohta, Katsumasa Eda, and Sumio Ayusawa of Hokken Co., Ltd. for providing fruiting bodies of strain HKB and all the collaboration during the project granted by the Agriculture, Forestry and Fisheries Research Council of Japan. I would also like to thank Dr. Junji Magae for valuable suggestions. Construction of cDNA libraries and sequencing was performed by TaKaRa Bio Inc. AFM was performed at the Research Institute of Biomolecule Metrology Co., Ltd., (RIBM) Enokido, Tsukuba, Ibaraki, 305-0853, Japan.
The author declares that they have no competing interests.
About this article
Cite this article
Magae, Y. Molecular characterization of a novel mycovirus in the cultivated mushroom, Lentinula edodes. Virol J 9, 60 (2012). https://doi.org/10.1186/1743-422X-9-60