Identification of microRNA-like RNAs in Cordyceps guangdongensis and their expression profile under differential developmental stages.

Cordyceps guangdongensis is a well-known fungus with high nutritional and medicinal value. The metabolite profile of C. guangdongensis is similar to that of Ophiocordyceps sinensis. In plants and animals, microRNAs play important roles in regulating gene expression at the post-transcriptional level. MicroRNA-like RNAs (milRNAs) have been documented in several macro-fungi. To comprehensively investigate the milRNAs in C. guangdongensis, three small RNA libraries from the differentially developmental stages were constructed. Twenty-six conserved milRNAs were identified, and 19 novel milRNA candidates were predicted. Among them, 20 milRNAs were differentially expressed across the developmental processes, and 12 milRNAs were verified using stem-loop quantitative real-time reverse transcription polymerase chain reaction. In addition, the potential target genes of milRNA were predicted to be involved in the development of fruiting bodies and metabolite biosynthesis. This study is the first to report the milRNAs of C. guangdongensis, and provides important insights into studies of milRNA regulation pathways in ascomycete fungi.

Mitogenome of Tolypocladium guangdongense.

Tolypocladium guangdongense is a high-value edible fungus with various medicinal and food safety properties. However, its evolutionary and genetic information is still limited. Mitochondrial genomes are potential models for molecular evolution and phylogenetic studies. In this study, we sequenced the complete mitogenome of T. guangdongense, demonstrating circular sequence of 46,102 bp, containing 14 standard protein-coding genes (PCGs), 2 ribosomal RNA subunit genes, and 28 tRNA genes. Phylogenetic analysis based on mitochondrial genes indicated that T. guangdongense was clustered into the Tolypocladium genus with high support value, based on the core PCG dataset. In addition, rps3 is also a suitable marker in the phylogenetic analysis in Hypocreales. Gene rearrangement analysis indicated that the gene order of PCGs was highly consistent in Hypocreales, and tRNA rearrangement events occurred in most species of Hypocreales; however, the rearrangement rates were not taxonomically correlated. Divergence time estimation based on the old fossil record and previous reports revealed that T. guangdongense originated approximately in the middle Cenozoic (42 Mya, 95% highest posterior density interval: 43-116) with the Tolypocladium genus differentiation. Our results provided more mitogenomic information of T. guangdongense and shed new insights into evolution of the Tolypocladium genus. KEY POINTS: • The general and unique features of T. guangdongense mitogenome are firstly reported. • Phylogenetic analysis further verified the taxonomic status of T. guangdongense. • Divergence time estimation provides more evolutionary information of T. guangdongense.

Identification of a Heat-Inducible Element of Cysteine Desulfurase Gene Promoter in Lentinula edodes.

Volatile organosulfur compounds are the main components that contribute to the unique aroma of dried Lentinula edodes. They are mainly generated during the hot-air drying process, and cysteine desulfurase is the key enzyme in this process. Temperature may be an essential factor of volatile organosulfur compound production by influencing the expression of the cysteine desulfurase gene. In this study, the promoter sequence of the cysteine desulfurase gene (pCS) was cloned and analyzed using bioinformatics tools. A series of 5'deletion fragments and site-directed mutations of pCS were constructed to identify the element that responds to heat stress. Six heat shock transcription factor (HSTF) binding sites were predicted by SCPD (The Promoter Database of Saccharomyces cerevisiae) and three of the binding sites were predicted by Yeastract (Yeast Search for Transcriptional Regulators and Consensus Tracking) in pCS. The results indicated that pCS was able to drive the expression of the EGFP (Enhanced Green Fluorescent Protein) gene in L. edodes. Moreover, the fluorescence intensity increased after heat stress. The changes in fluorescence intensity of different 5'deletion fragments showed that the heat response region was located between -500 bp and -400 bp in pCS. The site-directed mutation analysis further showed that the heat-inducible element was between -490 bp and -500 bp (TTTCTAGAAT) in pCS. Our results provide molecular insight for studying the formation of volatile organosulfur compounds in dried L. edodes.

Proteome and Transcriptome Reveal Involvement of Heat Shock Proteins and Indoleacetic Acid Metabolism Process in Lentinula Edodes Thermotolerance.

BACKGROUND/AIMS: Heat stress could cause huge losses for Lentinula edodes in China and other Asian cultivation areas. Yet our understanding of mechanism how to defend to heat stress is incomplete. METHODS: Using heat-tolerant and heat-sensitive strains of L. edodes, we reported a combined proteome and transcriptome analysis of L. edodes response to 40 °C heat stress for 24 h. Meanwhile, the effect of LeDnaJ on the thermotolerance and IAA (indoleacetic acid) biosynthesis in L. edodes was analyzed via the over-expression method. RESULTS: The proteome results revealed that HSPs (heat shock proteins) such as Hsp40 (DnaJ), Hsp70, Hsp90 and key enzymes involved in tryptophan and IAA metabolism process LeTrpE, LeTrpD, LeTam-1, LeYUCCA were more highly expressed in S606 than in YS3357, demonstrating that HSPs and tryptophan as well as IAA metabolism pathway should play an important role in thermotolerance. Over-expression of LeDnaJ gene in S606 strains showed better tolerance to heat stress. It was also documented that intracellular IAA accumulation of S606 (8-fold up) was more than YS3357 (2-fold up), and exogenous IAA enhanced L. edodes tolerance to heat stress. CONCLUSION: Our data support the interest of LeTrpE, LeDnaJ, tryptophan and IAA could play a pivotal role in enhancing organism thermotolerance.

The heat shock protein 40 LeDnaJ regulates stress resistance and indole-3-acetic acid biosynthesis in Lentinula edodes.

DnaJ proteins, termed heat shock proteins based on their molecular weight, function as molecular chaperones that play critical roles in regulating organism growth and development as well as adaptation to the environment. However, little has been reported on their gene function in higher basidiomycetes. Here, the heat shock protein 40 (LeDnaJ) gene was cloned and characterized from Lentinula edodes. RNA interference was used to explore the function of LeDnaJ in response to heat stress and Trichoderma atroviride. Integration of the target gene into the L. edodes genome was confirmed by Southern blot analysis, and the silence efficiency of LeDnaJ was analyzed by qRT-PCR. The results revealed that LeDnaJ silence caused defects in mycelial growth and resistance to heat stress and T. atroviride, but increased the mycelial density compared with the wild type (WT) strain S606. Additionally, the IAA content showed a more than 10-fold increase in the WT after heat stress, but an about two-fold increase in the two LeDnaJ RNAi transfortants (LeDnaJ-i-6 and LeDnaJ-i-8). Previous study has shown that enhanced IAA (indole-3-acetic acid) content enhanced the thermotolerance of the heat-sensitive strain YS3357. In this study, it was documented that IAA amendments could partly restore the resistance to T. atroviride and thermotolerance of the two LeDnaJ RNAi transformants. Overall, LeDnaJ is nvolved in fungal growth, T. atroviride resistance, and thermotolerance by regulating the IAA biosynthesis in L. edodes.

Selection of Reference Genes for qRT-PCR Analysis in Lentinula edodes after Hot-Air Drying.

Volatile sulfur compounds gradually develop in Lentinula edodes after hot-air drying, and many genes are involved in the generation of these sulfur compounds. The expression stability of reference genes may vary in a particular experimental treatment when analyzing their expressions by quantitative real-time polymerase chain reaction (qRT-PCR). In this study, the expression profile of 17 candidate genes was assessed in L. edodes under treatment at 50 °C for 0, 1, 2, and 3 h, and the expression stability of each reference gene was analyzed by three statistical algorithms, including geNorm, NormFinder, and BestKeeper. Results indicated that the two optimal reference genes for mycelium and fruiting body were CAC and DAHP as well as CAC and NUP, respectively. Additionally, CAC and DAHP were found to be the two most stable reference genes across the mycelium and fruiting body set. Our results will provide a genetic foundation for further research on the metabolism genes of sulfur compounds in L. edodes.

Biological and Molecular Characteristics of a Novel Partitivirus Infecting the Edible Fungus Lentinula edodes.

A new partitivirus named Lentinula edodes partitivirus 1 (LePV1) was isolated from a diseased L. edodes strain with severe degeneration of the mycelium and imperfect browning in bag cultures. The nucleotide sequences of LePV1 dsRNA-1 and dsRNA-2 were determined; they were 2,382 bp and 2,245 bp in length, and each contained a single ORF encoding RNA-dependent RNA polymerase (RdRp) and coat protein (CP), respectively. The purified virus preparation contained isometric particles 34 nm in diameter encapsidating these dsRNAs. Phylogenetic analyses showed LePV1 to be a new member of Betapartitivirus, with the RdRp sequence most closely related to Grapevine partitivirus. RT-PCR analysis showed that 27 of the 56 Chinese L. edodes core collection strains carry LePV1, with the virus being more common in wild strains than cultivated strains. In addition, qPCR analysis suggested that coinfection with L. edodes mycovirus HKB (LeV-HKB) could increase replication of the RdRp gene of LePV1. This study may be essential for the development of more accurate disease diagnostics and the formulation of control strategies for viral diseases in L. edodes.

Cordyceps guangdongensis lipid-lowering formula alleviates fat and lipid accumulation by modulating gut microbiota and short-chain fatty acids in high-fat diet mice.

Obesity has caused serious health and economic problems in the world. Cordyceps guangdongensis is a high-value macrofungus with broad application potential in the food and bio-medicine industry. This current study aimed to estimate the role of C. guangdongensis lipid-lowering compound formula (CGLC) in regulating fat and lipid accumulation, gut microbiota balance, short-chain fatty acid (SCFA) contents, and expression levels of genes involved in fat and lipid metabolism in high-fat diet (HFD) mice. The results showed that CGLC intervention markedly reduced body weights and fat accumulation in HFD mice, improved glucose tolerance and blood lipid levels, and decreased lipid droplet accumulation and fat vacuole levels in the liver. CGLC decreased the ratio of Firmicutes and Bacteroidetes and increased the relative abundances of Bacteroides (B. acidifaciens) and Bifidobacterium (B. pseudolongum). In addition, CGLC treatment significantly promoted the production of SCFAs and regulated the relative expression levels of genes involved in fat and lipid metabolism in liver. Association analysis showed that several species of Bacteroides and most of SCFAs were significantly associated with serum lipid indicators. These results suggested that CGLC is a novel candidate formulation for treating obesity and non-alcohol fatty liver by regulating gut microbiota, SCFAs, and genes involved in fat and lipid metabolism.

Biosynthetic Pathway and the Potential Role of Melatonin at Different Abiotic Stressors and Developmental Stages in Tolypocladium guangdongense.

Melatonin, a bioactive compound and an important signaling molecule produced in plants and animals, is involved in many biological processes. However, its function and synthetic pathways in fungi are poorly understood. Here, the samples from Tolypocladium guangdongense, a highly valued edible fungus with functional food properties, were collected under different experimental conditions to quantify the levels of melatonin and its intermediates. The results showed that the intracellular melatonin content was markedly improved by Congo red (CR), cold, and heat stresses; the levels of intracellular melatonin and its intermediates increased at the primordial (P) and fruiting body (FB) stages. However, the levels of most intermediates exhibited a notable decrease under CR stress. Several genes related to melatonin synthesis, excluding AADC (aromatic-L-amino-acid decarboxylase), were markedly upregulated at an early stage of CR stress but downregulated later. Compared to the mycelial stage, those genes were significantly upregulated at the P and FB stages. Additionally, exogenous melatonin promoted resistance to several abiotic stressors and P formation in T. guangdongense. This study is the first to report melatonin biosynthesis pathway in macro-fungi. Our results should help in studying the diversity of melatonin function and melatonin-synthesis pathways and provide a new viewpoint for melatonin applications in the edible-medicinal fungus.

Enhanced Expression of Thaumatin-like Protein Gene (LeTLP1) Endows Resistance to Trichoderma atroviride in Lentinula edodes.

Lentinula edodes (shiitake mushrooms) is heavily affected by the infection of Trichoderma atroviride, causing yield loss and decreases quality in shiitake mushrooms. The selection and breeding of fungal-resistant L. edodes species are an important approach to protecting L. edodes from T. atroviride infection. Herein, a highly resistant L. edodes strain (Y3334) and a susceptible strain (Y55) were obtained by using a resistance evaluation test. Transcriptome analyses and qRT-PCR detection showed that the expression level of LeTLP1 (LE01Gene05009) was strongly induced in response to T. atroviride infection in the resistant Y3334. Then, LeTLP1-silenced and LeTLP1-overexpression transformants were obtained. Overexpression of LeTLP1 resulted in resistance to T. atroviride. Compared with the parent strain Y3334, LeTLP1-silenced transformants had reduced resistance relative to T. atroviride. Additionally, the LeTLP1 protein (Y3334) exhibited significant antifungal activity against T. atroviride. These findings suggest that overexpression of LeTLP1 is a major mechanism for the resistance of L. edodes to T. atroviride. The molecular basis provides a theoretical basis for the breeding of resistant L. edodes strains and can eventually contribute to the mushroom cultivation industry and human health.

Distribution, evolution and expression of GATA-TFs provide new insights into their functions in light response and fruiting body development of Tolypocladium guangdongense.

BACKGROUND: Fungal GATA-type transcription factors (GATA-TFs) are a class of transcriptional regulators involved in various biological processes. However, their functions are rarely analyzed systematically, especially in edible or medicinal fungi, such as Tolypocladium guangdongense, which has various medicinal and food safety properties with a broad range of potential applications in healthcare products and the pharmaceutical industry. METHODS: GATA-TFs in T. guangdongense (TgGATAs) were identified using InterProScan. The type, distribution, and gene structure of TgGATAs were analyzed by genome-wide analyses. A phylogenetic tree was constructed to analyze their evolutionary relationships using the neighbor-joining (NJ) method. To explore the functions of GATA-TFs, conserved domains were analyzed using MEME, and cis-elements were predicted using the PlantCARE database. In addition, the expression patterns of TgGATAs under different light conditions and developmental stages were studied using qPCR. RESULTS: Seven TgGATAs were identified. They were randomly distributed on four chromosomes and contained one to four exons. Phylogenetic analysis indicated that GATA-TFs in each subgroup are highly conserved, especially for GATA1 to GATA5. Intron distribution analyses suggested that GATA1 and GATA3 possessed the most conserved gene structures. Light treatments induced the expression levels of TgGATA1 and TgGATA5-7, but the expression levels varied depending on the duration of illumination. The predicted protein structures indicate that TgGATA1 and TgGATA2 possess typical light-responsive domains and may function as photoreceptors to regulate downstream biological processes. TgGATA3 and TgGATA5 may be involved in nitrogen metabolism and siderophore biosynthesis, respectively. TgGATA6 and TgGATA7 possess unique Zn finger loop sequences, suggesting that they may have special functions. Furthermore, gene expression analysis indicated that TgGATA1 (WC1) was notably involved in mycelial color transformation, while other genes were involved in fruiting body development to some extent. These results provide valuable information to further explore the mechanisms through which TgGATAs are regulated during fruiting body development.

Transcriptome and proteome analyses reveal the regulatory networks and metabolite biosynthesis pathways during the development of Tolypocladium guangdongense.

Tolypocladium guangdongense has a similar metabolite profile to Ophiocordyceps sinensis, a highly regarded fungus used for traditional Chinese medicine with high nutritional and medicinal value. Although the genome sequence of T. guangdongense has been reported, relatively little is known about the regulatory networks for fruiting body development and about the metabolite biosynthesis pathways. In order to address this, an analysis of transcriptome and proteome at differential developmental stages of T. guangdongense was performed. In total, 9076 genes were found to be expressed and 2040 proteins were identified. There were a large number of genes that were significantly differentially expressed between the mycelial stage and the stages. Interestingly, the correlation between the transcriptomic and proteomic data was low, suggesting the importance of the post-transcriptional processes in the growth and development of T. guangdongense. Among the genes/proteins that were both differentially expressed during the developmental process, there were numerous heat shock proteins and transcription factors. In addition, there were numerous proteins involved in terpenoid, ergosterol, adenosine and polysaccharide biosynthesis that also showed significant downregulation in their expression levels during the developmental process. Furthermore, both tryptophan and tryptamine were present at higher levels in the primordium stage. However, indole-3-acetic acid (IAA) levels continuously decreased as development proceeded, and the enzymes involved in IAA biosynthesis were also clearly differentially downregulated. These data could be meaningful in studying the molecular mechanisms of fungal development, and for the industrial and medicinal application of macro-fungi.

Selection and validation of reliable reference genes for Tolypocladium guangdongense gene expression analysis under differentially developmental stages and temperature stresses.

Tolypocladium guangdongense, formerly known as Cordyceps guangdongensis, is a widely cultivated fungus of the Cordyceps s.l. species that has been investigated over the last 12 years. It has the potential to be used in a number of applications in the health and pharmaceutical industries for it has shown its high nutritional and medicinal values according to previous animal studies. qRT-PCR (quantitative reverse transcription polymerase chain reaction) is extensively used to analyze the expression pattern and molecular mechanisms of functional genes under differentially experimental conditions. The expression stability of reference genes used for normalization determines the reliability of qRT-PCR results, indicating the importance of selection and validation of reference genes before gene expression analysis. In the present study, three statistical algorithms, geNorm, NormFinder and BestKeeper, were used for analyzing the expression stability of nineteen candidate reference genes (CRGs) in T. guangdongense. Investigation were carried out under differentially experimental conditions, which included differentially developmential stages (mycelia, primordia, young and mature fruiting bodies), different carbon sources, cold and heat stresses. The results showed that histone H4 and tubulin beta chain 2 (ß-tub2) were the most and least stable genes, respectively, across all the experimental samples. Moreover, analysis of individual data sets exhibited different stability and expression profiles of reference genes. The vacuolar protein sorting gene VPS was the most stable gene expressed under the differentially developmental stages and temperature stresses, whereas H4 was the most stably expressed gene under different carbon sources. Therefore, it can be proposed that VPS and H4 are the preferred reference genes for normalization of gene expression under different experimental conditions. The results of our present study will enable more accurate evaluation of gene expression in T. guangdongense using the optimal reference gene for qRT-PCR analysis.

Expression profiling of Cordyceps DnaJ protein family in Tolypocladium guangdongense during developmental and temperature stress processes.

DnaJ is an important molecular chaperone, with significant roles in growth, development, and stress resistance. Studies on the DnaJ gene family in macro-fungi such as Cordyceps spp. s.l. is scare. In this study, 22, 20, and 24 putative DnaJ genes were identified in Tolypocladium guangdongense, Ophiocordyceps sinensis, and C. militaris, respectively. They were classified into four groups based on the presence of the J, zinc finger, and C-terminal domains. We mainly studied the T. guangdongense DnaJ genes being located in the endoplasmic reticulum, cytoplasm, mitochondrion, and nucleus. Phylogenetic analysis revealed gene duplications during the evolutionary process. Multiple cis-elements and transcription factor binding sites were observed in the promoter, suggesting their involvement in the response to multiple stresses. qRT-PCR analysis showed that 63.63% and 45.45% of T. guangdongense DnaJ genes were differentially expressed under cold and heat stress, respectively, indicating their involvement in the response to temperature stress. Many T. guangdongense DnaJ genes in the primordium and fruiting body exhibited differential expression, in comparison to those in the mycelium, suggesting a regulatory role in its growth and development process. These findings will facilitate further functional analysis, and provide information on the classification and conservative functions of DnaJ proteins in macro-fungi.

Hsp40 Protein LeDnaJ07 Enhances the Thermotolerance of Lentinula edodes and Regulates IAA Biosynthesis by Interacting LetrpE.

Our previous study found that LeDnaJ07 RNAi decreased Lentinula edodes resistance to heat stress and Trichoderma atroviride infection. In this study, the structure and function of the LeDnaJ07 gene was analyzed by gene cloning and overexpression in L. edodes stress-sensitive strain YS55 via the Agrobacterium-mediated transformation method. Transformants were confirmed by qRT-PCR, fluorescence observation and Southern blotting. Overexpression of LeDnaJ07 in YS55 not only enhanced L. edodes mycelial resistance to heat stress but also facilitated mycelial growth. In the presence of heat stress, the intracellular IAA content showed a significant increase in the two LeDnaJ07 overexpression strains but only a slight change in the YS55 wild type strain. Moreover, the interaction between LeDnaJ07 and LetrpE was demonstrated via Y2H and BiFC assays. These results suggested that LeDnaJ07 may be involved in regulating IAA biosynthesis and the resistance of L. edodes to heat stresses via interacting with LetrpE.

Selection and Validation of Reference Genes for qRT-PCR in Lentinula edodes under Different Experimental Conditions.

Lentinula edodes is the most consumed mushroom in Asia due to its nutritional and medicinal values, and the optimal reference gene is crucial for normalization of its gene expression analysis. Here, the expression stability of 18 candidate reference genes (CRGs) in L. edodes was analyzed by three statistical algorithms (geNorm, NormFinder and BestKeeper) under different stresses (heat, cadmium excess and Trichoderma atroviride infection), different substrates (straw, sawdust and corn stalk) and different development stages (mycelia, primordia and fruit bodies). Among the 18 CRGs, 28S, Actin and α-tub exhibited the highest expression stability in L. edodes under all conditions, while GPD, SPRYP and MSF showed the least stable expression. The best reference gene in different conditions was different. The pairwise variation values showed that two genes would be sufficient for accurate normalization under different conditions of L. edodes. This study will contribute to more accurate estimation of the gene relative expression levels under different conditions using the optimal reference gene in qRT-PCR (quantitative reverse transcription polymerase chain reaction) analysis.

The DnaJ Gene Family in Shiitake Culinary-Medicinal Mushroom, Lentinus edodes (Agaricomycetes): Comprehensive Identification, Characterization, and Expression Profiles under Different Conditions.

The DnaJ proteins, also called heat-shock protein 40 based on their molecular weight, play significant roles in organism growth and development and resistance to abiotic and biotic stresses. However, studies on the DnaJ gene family in Lentinus edodes (= Lentinula edodes) are less well known. In this study, 29 putative L. edode DnaJ genes (LeDnaJ01 to LeDnaJ29) were identified using bioinformatics analysis and were classified into four groups according to the presence of the J protein and zinc finger as well as C-terminal domain. Multiple cis elements related to the phytohormone and stresses were found in the promotor region of the LeDnaJ genes. In addition, qRT-PCR analysis revealed that 79.31% of LeDnaJ genes were induced by cadmium, 55.17% were induced by Trichoderma atroviride, and 37.93% were induced by heat stress, indicating that the LeDnaJ proteins might participate in the response of L. edodes to the multiple stresses. Meanwhile, qRT-PCR analysis also revealed that all LeDnaJs are expressed in at least one development stage, indicating that they could be involved in the process of L. edodes growth and development and the response to the abiotic and biotic stresses. Taken together, these results advance the functional analysis of DnaJ genes in Basidiomycetes.

First report of a cross-kingdom pathogenic bacterium, Achromobacter xylosoxidans isolated from stipe-rot Coprinus comatus.

Coprinus comatus is an edible mushroom widely cultivated in China as a delicious food. Various diseases have occurred on C. comatus with the cultivated area increasing. In this study, the pathogenic bacterium JTG-B1, identified as Achromobacter xylosoxidans by 16S rDNA and nrdA gene sequencing, was isolated from edible mushroom Coprinus comatus with serious rot disease on its stipe. A. xylosoxidans has been confirmed as an important opportunistic human pathogenic bacterium and has been isolated from respiratory samples from cystic fibrosis. It is widely distributed in the environment. Here, we first report that fungi can also serve as a host for A. xylosoxidans. We confirmed that it can cross-kingdom infect between animals (mice) and fungi (C. comatus). The results of pathogenicity tests, physiological, biochemical and genotyping analysis of A. xylosoxidans from different hosts suggested that different strain of A. xylosoxidans may have pathogenicity differentiation. A. xylosoxidans not only is pathogenic to C. comatus but also may threaten human health.

Diversity and effect of Trichoderma spp. associated with green mold disease on Lentinula edodes in China.

Lentinula edodes, one of the most important edible mushrooms in China, is affected heavily by the infection of green mold that overgrows mushroom mycelia. We collected the diseased samples from main L. edodes cultivation regions in China to characterize the pathogen and to study the effect of Trichoderma spp. on L. edodes species. We identified six Trichoderma species, that is, T. harzianum, T. atroviride, T. viride, T. pleuroticola, T. longibrachiatum, and T. oblongisporum based on the internal transcribed spacer or tef1-α sequences and morphology characteristics. In confrontation cultures on Petri plates or in tubes, and in L. edodes cultures in a medium containing Trichoderma metabolites, L. edodes mycelia were not only distorted and swollen, but also inhibited by Trichoderma isolates. It is not possible that adjusting pH value or temperature is used for controlling L. edodes green disease, because the growth of most of Trichoderma isolates and L. edodes shared similar pH and temperature conditions.

Genome Sequence of the Edible Cultivated Mushroom Lentinula edodes (Shiitake) Reveals Insights into Lignocellulose Degradation.

Lentinula edodes, one of the most popular, edible mushroom species with a high content of proteins and polysaccharides as well as unique aroma, is widely cultivated in many Asian countries, especially in China, Japan and Korea. As a white rot fungus with lignocellulose degradation ability, L. edodes has the potential for application in the utilization of agriculture straw resources. Here, we report its 41.8-Mb genome, encoding 14,889 predicted genes. Through a phylogenetic analysis with model species of fungi, the evolutionary divergence time of L. edodes and Gymnopus luxurians was estimated to be 39 MYA. The carbohydrate-active enzyme genes in L. edodes were compared with those of the other 25 fungal species, and 101 lignocellulolytic enzymes were identified in L. edodes, similar to other white rot fungi. Transcriptome analysis showed that the expression of genes encoding two cellulases and 16 transcription factor was up-regulated when mycelia were cultivated for 120 minutes in cellulose medium versus glucose medium. Our results will foster a better understanding of the molecular mechanism of lignocellulose degradation and provide the basis for partial replacement of wood sawdust with agricultural wastes in L. edodes cultivation.