Comparative transcriptome analysis of dikaryotic mycelia and mature fruiting bodies in the edible mushroom Lentinula edodes

Near-gapless genome and transcriptome analyses provide insights into fruiting body development in Lentinula edodes

Fruiting body development in macrofungi is an intensive research subject. In this study, high-quality genomes were assembled for two sexually compatible monokaryons from a heterokaryotic Lentinula edodes strain WX1, and variations in L. edodes genomes were analyzed. Specifically, differential gene expression and allele-specific expression (ASE) were analyzed using the two monokaryotic genomes and transcriptome data from four different stages of fruiting body development in WX1. Results revealed that after aeration, mycelia sensed cell wall stress, pheromones, and a decrease in CO2 concentration, leading to up-regulated expression in genes related to cell adhesion, cell wall remodeling, proteolysis, and lipid metabolism, which may promote primordium differentiation. Aquaporin genes and those related to proteolysis, mitosis, lipid, and carbohydrate metabolism may play important roles in primordium development, while genes related to tissue differentiation and sexual reproduction were active in fruiting body. Several essential genes for fruiting body development were allele-specifically expressed and the two nuclear types could synergistically regulate fruiting body development by dominantly expressing genes with different functions. ASE was probably induced by long terminal repeat-retrotransposons. Findings here contribute to the further understanding of the mechanism of fruiting body development in macrofungi.

Integrative Analysis of Transcriptome and Metabolome Sheds Light on Flavonoid Biosynthesis in the Fruiting Body of Stropharia rugosoannulata

Flavonoids are a diverse family of natural compounds that are widely distributed in plants and play a critical role in plant growth, development, and stress adaptation. In recent years, the biosynthesis of flavonoids in plants has been well-researched, with the successive discovery of key genes driving this process. However, the regulation of flavonoid biosynthesis in fungi remains unclear. Stropharia rugosoannulata is an edible mushroom known for its high nutritional and pharmacological value, with flavonoids being one of its main active components. To investigate the flavonoid content of S. rugosoannulata, a study was conducted to extract and determine the total flavonoids at four stages: young mushroom (Ym), gill (Gi), maturation (Ma), and parachute-opening (Po). The findings revealed a gradual increase in total flavonoid concentration as the fruiting body developed, with significant variations observed between the Ym, Gi, and Ma stages. Subsequently, we used UPLC-MS/MS and transcriptome sequencing (RNA-seq) to quantify the flavonoids and identify regulatory genes of Ym, Gi, and Ma. In total, 53 flavonoid-related metabolites and 6726 differentially expressed genes (DEGs) were identified. Through KEGG pathway enrichment analysis, we identified 59 structural genes encoding flavonoid biosynthesis-related enzymes, most of which were up-regulated during the development of the fruiting body, consistent with the accumulation of flavonoids. This research led to the establishment of a comprehensive transcriptional metabolic regulatory network encompassing flavonoids, flavonoid synthases, and transcription factors (TFs). This represents the first systematic exploration of the molecular mechanism of flavonoids in the fruiting of fungi, offering a foundation for further research on flavonoid mechanisms and the breeding of high-quality S. rugosoannulata.

Identification and Characterization of a Predominant Hydrophobin in the Edible Mushroom Grifola frondosa

Hydrophobins (HFBs) are a group of small, secreted amphipathic proteins of fungi with multiple physiological functions and potential commercial applications. In this study, HFB genes of the edible mushroom, Grifola frondosa, were systematically identified and characterized, and their transcriptional profiles during fungal development were determined. In total, 19 typical class I HFB genes were discovered and bioinformatically analyzed. Gene expression profile examination showed that Gf.hyd9954 was particularly highly upregulated during primordia formation, suggesting its major role as the predominant HFB in the lifecycle of G. frondosa. The wettability alteration profile and the surface modification ability of recombinant rGf.hyd9954 were greater than for the Grifola HFB HGFII-his. rGf.hyd9954 was also demonstrated to form the typical class I HFB characteristic-rodlet bundles. In addition, rGf.hyd9954 was shown to possess nanoparticle characteristics and emulsification activities. This research sheds light on the regulation of fungal development and its association with the expression of HFB genes.

Hydrated lime promoted the polysaccharide content and affected the transcriptomes of Lentinula edodes during brown film formation

Brown film formation, a unique developmental stage in the life cycle of Lentinula edodes, is essential for the subsequent development of fruiting bodies in L. edodes cultivation. The pH of mushroom growth substrates are usually adjusted with hydrated lime, yet the effects of hydrated lime on cultivating L. edodes and the molecular mechanisms associated with the effects have not been studied systemically. We cultivated L. edodes on substrates supplemented with 0% (CK), 1% (T1), 3% (T2), and 5% (T3) hydrated lime (Ca (OH)2), and applied transcriptomics and qRT-PCR to study gene expression on the brown film formation stage. Hydrated lime increased polysaccharide contents in L. edodes, especially in T2, where the 5.3% polysaccharide content was approximately 1.5 times higher than in the CK. The addition of hydrated lime in the substrate promoted laccase, lignin peroxidase and manganese peroxidase activities, implying that hydrated lime improved the ability of L. edodes to decompose lignin and provide nutrition for its growth and development. Among the annotated 9,913 genes, compared to the control, 47 genes were up-regulated and 52 genes down-regulated in T1; 73 genes were up-regulated and 44 were down-regulated in T2; and 125 genes were up-regulated and 65 genes were down-regulated in T3. Differentially expressed genes (DEGs) were enriched in the amino acid metabolism, lipid metabolism and carbohydrate metabolism related pathways. The carbohydrate-active enzyme genes up-regulated in the hydrated lime treatments were mostly glycosyl hydrolase genes. The results will facilitate future optimization of L. edodes cultivation techniques and possibly shortening the production cycle.

Comparative transcriptome analysis on candidate genes associated with fruiting body growth and development in Lyophyllum decastes

Lyophyllum decastes is a mushroom that is highly regarded for its culinary and medicinal properties. Its delectable taste and texture make it a popular choice for consumption. To gain a deeper understanding of the molecular mechanisms involved in the development of the fruiting body of L. decastes, we used RNA sequencing to conduct a comparative transcriptome analysis. The analysis encompassed various developmental stages, including the vegetative mycelium, primordial initiation, young fruiting body, medium-size fruiting body, and mature fruiting body stages. A range of 40.1 to 60.6 million clean reads were obtained, and de novo assembly generated 15,451 unigenes with an average length of 1,462.68 bp. Functional annotation of transcriptomes matched 76.84% of the unigenes to known proteins available in at least one database. The gene expression analysis revealed a significant number of differentially expressed genes (DEGs) between each stage. These genes were annotated and subjected to Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses. Highly differentially expressed unigenes were also identified, including those that encode extracellular enzymes, transcription factors, and signaling pathways. The accuracy of the RNA-Seq and DEG analyses was validated using quantitative PCR. Enzyme activity analysis experiments demonstrated that the extracellular enzymes exhibited significant differences across different developmental stages. This study provides valuable insights into the molecular mechanisms that underlie the development of the fruiting body in L. decastes.

Role and structure of the small subunit forming heterodimers with laccase-like enzymes

Unlike laccases sensu stricto, which are usually monomeric enzymes, laccase-like enzymes recently re-classified as Novel Laccases (NLACs) are characterized by the formation of heterodimers with small proteins (subunits) of unknown function. Here the NLAC from Pleurotus eryngii (PeNL) and a small protein selected from the fungal genome, that is homologous to reported POXA3 from Pleurotus ostreatus, were produced in Aspergillus oryzae separately or together. The two proteins interacted regardless of whether the small subunit was co-expressed or exogenously added to the enzyme. The stability and catalytic activity of PeNL was significantly enhanced in the presence of the small subunit. Size exclusion chromatography-multi angle light scattering (SEC-MALS) analysis confirmed that the complex PeNL-ss is a heterodimer of 77.4 kDa. The crystallographic structure of the small protein expressed in Escherichia coli was solved at 1.6 Å resolution. This is the first structure elucidated of a small subunit of a NLAC. The helix bundle structure of the small subunit accommodates well with the enzyme model structure, including interactions with specific regions of NLACs and some amino acid residues of the substrate-binding loops.

Transcriptome analysis provides insight into gamma irradiation delaying quality deterioration of postharvest Lentinula edodes during cold storage

To better determine how gamma irradiation (GI) improves abiotic stress resistance, a transcriptome analysis of postharvest L. edodes in response to 1.0 kGy GI was conducted, and further the underlying mechanism of GI in delaying quality deterioration over 20 d of cold storage was explored. The results suggested that GI was involved in multiple metabolic processes in irradiated postharvest L. edodes. In comparison with the control group, the GI group contained 430 differentially expressed genes, including 151 upregulated genes and 279 downregulated genes, which unveiled characteristic expression profiles and pathways. The genes involved in the pentose phosphate pathway were mainly upregulated and the expression level of the gene encoding deoxy-D-gluconate 3-dehydrogenase was 9.151-fold higher. In contrast, the genes related to other energy metabolism pathways were downregulated. Concurrently, GI inhibited the expression of genes associated with delta 9-fatty acid desaturase, ribosomes, and HSP20; thus, GI helped postpone the degradation of lipid components, suppress transcriptional metabolism and regulate the stress response. Additionally, the metabolic behavior of DNA repair induced by GI intensified by noticeable upregulation. These regulatory effects could play a potential and nonnegligible role in delaying the deterioration of L. edodes quality. The results provide new information on the regulatory mechanism of postharvest L. edodes when subjected to 1.0 kGy GI during cold storage.

Comparative analysis of proteomes and transcriptomes revealed the molecular mechanism of development and nutrition of Pleurotus giganteus at different fruiting body development stages

Pleurotus giganteus is a commercially cultivated high-temperature mushroom. Investigating the molecular mechanism of fruiting body development will help us to better understand the regulation of substrates and energy in this process. However, little information has been reported on the development and nutrients of the P. giganteus fruiting body. In the present study, P. giganteus is cultivated in a climate chamber, and comparative transcriptome, proteome, and nutritional analysis of P. giganteus fruiting bodies were performed. Our results revealed that Cytochrome P450 monooxygenases and hydrophobin proteins play important roles during the differentiation in the elongation stage. Later, carbon metabolism dominate the fruiting body metabolism and genes related to the carbohydrate metabolic process, glycolytic process, and gluconeogenesis were up-regulated in the mature fruiting bodies. The up-regulation of carbohydrate substrates utilization CAZymes genes and inconsistent protein expression in pileus indicated a reverse transportation of mRNA from the fruiting body to vegetative mycelia. In addition, protein concentration in the pileus is higher than that in the stem, while the stem is the major nitrogen metabolic and amino acid synthetic location. The integrated transcriptomic, proteomic, and nutritional analysis indicated a two-way transportation of substrates and mRNAs in P. giganteus. Stem synthesizes amino acids and transported them to pileus with reducing sugars, while pileus induces the expression of substrate degradation mRNA according to the needs of growth and development and transports them in the other direction.

Lessons on fruiting body morphogenesis from genomes and transcriptomes of Agaricomycetes

Fruiting bodies (sporocarps, sporophores or basidiomata) of mushroom-forming fungi (Agaricomycetes) are among the most complex structures produced by fungi. Unlike vegetative hyphae, fruiting bodies grow determinately and follow a genetically encoded developmental program that orchestrates their growth, tissue differentiation and sexual sporulation. In spite of more than a century of research, our understanding of the molecular details of fruiting body morphogenesis is still limited and a general synthesis on the genetics of this complex process is lacking. In this paper, we aim at a comprehensive identification of conserved genes related to fruiting body morphogenesis and distil novel functional hypotheses for functionally poorly characterised ones. As a result of this analysis, we report 921 conserved developmentally expressed gene families, only a few dozens of which have previously been reported to be involved in fruiting body development. Based on literature data, conserved expression patterns and functional annotations, we provide hypotheses on the potential role of these gene families in fruiting body development, yielding the most complete description of molecular processes in fruiting body morphogenesis to date. We discuss genes related to the initiation of fruiting, differentiation, growth, cell surface and cell wall, defence, transcriptional regulation as well as signal transduction. Based on these data we derive a general model of fruiting body development, which includes an early, proliferative phase that is mostly concerned with laying out the mushroom body plan (via cell division and differentiation), and a second phase of growth via cell expansion as well as meiotic events and sporulation. Altogether, our discussions cover 1 480 genes of Coprinopsis cinerea, and their orthologs in Agaricus bisporus, Cyclocybe aegerita, Armillaria ostoyae, Auriculariopsis ampla, Laccaria bicolor, Lentinula edodes, Lentinus tigrinus, Mycena kentingensis, Phanerochaete chrysosporium, Pleurotus ostreatus, and Schizophyllum commune, providing functional hypotheses for ~10 % of genes in the genomes of these species. Although experimental evidence for the role of these genes will need to be established in the future, our data provide a roadmap for guiding functional analyses of fruiting related genes in the Agaricomycetes. We anticipate that the gene compendium presented here, combined with developments in functional genomics approaches will contribute to uncovering the genetic bases of one of the most spectacular multicellular developmental processes in fungi. Citation: Nagy LG, Vonk PJ, Künzler M, Földi C, Virágh M, Ohm RA, Hennicke F, Bálint B, Csernetics Á, Hegedüs B, Hou Z, Liu XB, Nan S, M. Pareek M, Sahu N, Szathmári B, Varga T, Wu W, Yang X, Merényi Z (2023). Lessons on fruiting body morphogenesis from genomes and transcriptomes of Agaricomycetes. Studies in Mycology 104: 1-85. doi: 10.3114/sim.2022.104.01.

HOG1 Mitogen-Activated Protein Kinase Pathway–Related Autophagy Induced by H2O2 in Lentinula edodes Mycelia

Mycelial ageing is associated with ROS and autophagy in Lentinula edodes. However, the underlying cellular and molecular mechanisms between ROS and autophagy remain obscure. This study induced autophagy in L. edodes mycelia through exogenous H2O2 treatment. Results showed that 100 μM H2O2 treatment for 24 h significantly inhibited mycelial growth. H2O2 caused the depolarisation of MMP and accumulation of TUNEL-positive nuclei, which was similar to the ageing phenotype of L. edodes mycelia. Transcriptome analysis showed that differentially expressed genes were enriched in the mitophagic, autophagic, and MAPK pathways. LeAtg8 and LeHog1 were selected as hub genes. RNA and protein levels of LeATG8 increased in the H2O2-treated mycelia. Using fluorescent labelling, we observed for the first time the classic ring structure of autophagosomes in a mushroom, while 3D imaging suggested that these autophagosomes surrounded the nuclei to degrade them at specific growth stages. Phospho-LeHOG1 protein can translocate from the cytoplasm to the nucleus to regulate mycelial cells, resisting ROS-induced oxidative stress. Furthermore, LeATG8 expression was suppressed when LeHOG1 phosphorylation was inhibited. These results suggest that the LeATG8-dependent autophagy in L. edodes mycelial is closely associated with the activity or even phosphorylation of LeHOG1.

Application of omics technology in the research on edible fungi

Edible fungus is a large fungus distributed all over the world and used as food and medicine. But people's understanding of edible fungi is not as much as that of ordinary crops, so people have started a number of research on edible fungi in recent years. With the development of science and technology, omics technology has gradually walked into people's vision. Omics technology has high sensitivity and wide application range, which is favored by researchers. The application of omics technology to edible fungus research is a major breakthrough, which has transferred edible fungus research from artificial cultivation to basic research. Now omics technology in edible fungi has been flexibly combined with other research methods, involving multiple studies of edible fungus, such as genetic breeding, growth and development, stress resistance, and the use of special components in edible fungus as pharmaceutical additives. It is believed that in the future, the research of edible fungi will also be brought to a deeper level with the help of omics technology. This paper introduces the application progress of modern omics technology to the study on edible fungi and mentions the application prospect of edible fungi research with the constant development of omics technology, thereby providing ideas for the follow-up in-depth research on edible fungi.

Secrets of Flavonoid Synthesis in Mushroom Cells

Flavonoids are chemical compounds that occur widely across the plant kingdom. They are considered valuable food additives with pro-health properties, and their sources have also been identified in other kingdoms. Especially interesting is the ability of edible mushrooms to synthesize flavonoids. Mushrooms are usually defined as a group of fungal species capable of producing macroscopic fruiting bodies, and there are many articles considering the content of flavonoids in this group of fungi. Whereas the synthesis of flavonoids was revealed in mycelial cells, the ability of mushroom fruiting bodies to produce flavonoids does not seem to be clearly resolved. This article, as an overview of the latest key scientific findings on flavonoids in mushrooms, outlines and organizes the current state of knowledge on the ability of mushroom fruiting bodies to synthesize this important group of compounds for vital processes. Putting the puzzle of the current state of knowledge on flavonoid biosynthesis in mushroom cells together, we propose a universal scheme of studies to unambiguously decide whether the fruiting bodies of individual mushrooms are capable of synthesizing flavonoids.

Comparative transcriptome analysis revealed candidate genes involved in fruiting body development and sporulation in Ganoderma lucidum

Ganoderma lucidum is an edible mushroom highly regarded in the traditional Chinese medicine. To better understand the molecular mechanisms underlying fruiting body development in G. lucidum, transcriptome analysis based on RNA sequencing was carried out on different developmental stages: mycelium (G1); primordium (G2); young fruiting body (G3); mature fruiting body (G4); fruiting body in post-sporulation stage (G5). In total, 26,137 unigenes with an average length of 1078 bp were de novo assembled. Functional annotation of transcriptomes matched 72.49% of the unigenes to known proteins available in at least one database. Differentially expressed genes (DEGs) were identified between the evaluated stages: 3135 DEGs in G1 versus G2; 120 in G2 versus G3; 3919 in G3 versus G4; and 1012 in G4 versus G5. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of DEGs identified in G1 versus G2 revealed that, in addition to global and overview maps, enriched pathways were related to amino acid metabolism and carbohydrate metabolism. In contrast, DEGs identified in G2 versus G3 were mainly assigned to the category of metabolism of amino acids and their derivatives, comprising mostly upregulated unigenes. In addition, highly expressed unigenes associated with the transition between different developmental stages were identified, including those encoding hydrophobins, cytochrome P450s, extracellular proteases, and several transcription factors. Meanwhile, highly expressed unigenes related to meiosis such as DMC1, MSH4, HOP1, and Mek1 were also analyzed. Our study provides important insights into the molecular mechanisms underlying fruiting body development and sporulation in G. lucidum.

Mycelium vs. Fruiting Bodies of Edible Fungi-A Comparison of Metabolites

Edible mushrooms are widely appreciated for their appealing flavours, low caloric values and high content of presumably health-protecting metabolites. Their long history of safe use together with the looming worldwide food crisis have revived the idea of generating meat analogues and protein isolates by the controlled fermentation of mycelia of these edible fungi as a dietary option. The occurrence of proteins, polysaccharides, smaller metabolites, metal ions and toxins in mycelia and fruiting bodies is compared among the three most popular species, Agaricus bisporus (button mushroom), Pleurotus ostreatus (oyster mushroom), Lentinus edodes (shiitake) and some closely related species. Large effects of substrate chemistry, strain, developmental stage and ecological interactions result in a wide variation of the concentrations of some metabolites in both mycelial cells and fruiting bodies. This is obviously a result of the high adaptation abilities required to survive in natural habitats. Fungal bioprocesses are decoupled from agricultural production and can be operated anytime, anywhere, and on any scale according to demand. It is concluded that fungal biomass, if produced under food-grade conditions and on an industrial scale, could provide a safe and nutritious meat substitute and protein isolates with a high biological value for future vegan foods.

The Molecular Mechanism of Yellow Mushroom (Floccularia luteovirens) Response to Strong Ultraviolet Radiation on the Qinghai-Tibet Plateau

The Qinghai-Tibet Plateau (QTP) is the highest plateau in the world, and its ultraviolet (UV) radiation is much greater than that of other regions in the world. Yellow mushroom (Floccularia luteovirens) is a unique and widely distributed edible fungus on the QTP. However, the molecular mechanism of F. luteovirens’s response to strong UV radiation remains unclear. Herein, we reported the 205 environmental adaptation and information processing genes from genome of F. luteovirens. In addition, we assembled the RNA sequence of UV-affected F. luteovirens at different growth stages. The results showed that in response to strong UV radiation, a total of 11,871 significantly different genes were identified, of which 4,444 genes in the vegetative mycelium (VM) stage were significantly different from the young fruiting bodies (YFB) stage, and only 2,431 genes in the YFB stage were significantly different from fruiting bodies (FB) stage. A total of 225 differentially expressed genes (DEGs) were found to be involved in environmental signal transduction, biochemical reaction preparation and stress response pathway, pigment metabolism pathway, and growth cycle regulation, so as to sense UV radiation, promote repair damage, regulate intracellular homeostasis, and reduce oxidative damage of UV radiation. On the basis of these results, a molecular regulation model was proposed for the response of F. luteovirens to strong UV radiation. These results revealed the molecular mechanism of adaptation of F. luteovirens adapting to strong UV radiation, and provided novel insights into mechanisms of fungi adapting to extreme environmental conditions on the QTP; the production the riboflavin pigment of the endemic fungi (Yellow mushroom) in the QTP was one of the response to extreme environment of the strong UV radiation.

Transcriptome Analysis Reveals Candidate Genes Involved in Light-Induced Primordium Differentiation in Pleurotus eryngii

Pleurotus eryngii, a highly valued edible fungus, is one of the major commercially cultivated mushrooms in China. The development of P. eryngii, especially during the stage of primordium differentiation, is easily affected by light. However, the molecular mechanism underlying the response of primordium differentiation to light remains unknown. In the present study, primordium expression profiles under blue-light stimulation, red-light stimulation, and exposure to darkness were compared using high-throughput sequencing. A total of 16,321 differentially expressed genes (DEGs) were identified from three comparisons. GO enrichment analysis showed that a large number of DEGs were related to light stimulation and amino acid biosynthesis. KEGG analyses demonstrated that the MAPK signaling pathway, oxidative phosphorylation pathway, and RNA transport were most active during primordium differentiation. Furthermore, it was predicted that the blue-light photoreceptor WC-1 and Deoxyribodipyrimidine photolyase PHR play important roles in the primordium differentiation of P. eryngii. Taken together, the results of this study provide a speculative mechanism that light induces primordium differentiation and a foundation for further research on fruiting body development in P. eryngii.

Transcriptomics Analysis of Primordium Formation in Pleurotus eryngii

Primordium formation is an important stage preceding the growth and development of the Pleurotus eryngii fruiting body. However, the molecular mechanisms underlying primordium formation remain unclear. In the present study, comparative transcriptomics was performed between mature mycelia and primordium to analyze the transcriptional properties during primordium formation in P. eryngii. A total of 19,655 differentially expressed genes (10,718 upregulated genes and 8937 downregulated genes) were identified. These differentially expressed genes were involved in cell wall degradation, carbohydrate hydrolysis, light perception, and cAMP signal transduction. These results aid further understanding of the transcriptional changes and the molecular processes underlying primordium formation and differentiation, which may lay the foundation for improving the cultivation and quality control of P. eryngii.

Whole-Genome and Transcriptome Sequencing of Phlebopus portentosus Reveals Its Associated Ectomycorrhizal Niche and Conserved Pathways Involved in Fruiting Body Development

Phlebopus portentosus (Berk. and Broome) Boedijin, a widely consumed mushroom in China and Thailand, is the first species in the order Boletaceae to have been industrially cultivated on a large scale. However, to date, the lignocellulose degradation system and molecular basis of fruiting body development in P. portentosus have remained cryptic. In the present study, genome and transcriptome sequencing of P. portentosus was performed during the mycelium (S), primordium (P), and fruiting body (F) stages. A genome of 32.74 Mb with a 48.92% GC content across 62 scaffolds was obtained. A total of 9,464 putative genes were predicted from the genome, of which the number of genes related to plant cell wall-degrading enzymes was much lower than that of some saprophytic mushrooms with specific ectomycorrhizal niches. Principal component analysis of RNA-Seq data revealed that the gene expression profiles at all three stages were different. The low expression of plant cell wall-degrading genes also confirmed the limited ability to degrade lignocellulose. The expression profiles also revealed that some conserved and specific pathways were enriched in the different developmental stages of P. portentosus. Starch and sucrose metabolic pathways were enriched in the mycelium stage, while DNA replication, the proteasome and MAPK signaling pathways may be associated with maturation. These results provide a new perspective for understanding the key pathways and hub genes involved in P. portentosus development.

Comparative transcriptome analysis reveals the differential response to cadmium stress of two Pleurotus fungi: Pleurotus cornucopiae and Pleurotus ostreatus

Pleurotus has great potential for heavy metal mycoremediation. Using comparative transcriptome analysis, the response of Pleurotus ostreatus and Pleurotus cornucopiae under Cd contamination was evaluated. P. ostreatus and P. cornucopia accumulated 0.34 and 0.46 mg/g Cd in mycelium, respectively. Cd removal elevated with its concentration elevation, which reached 56.47% and 54.60% for P. ostreatus and P. cornucopia with Cd at 20 mg/L. Low-level Cd (≤ 1 mg/L) had no significant influence on either fungus, while varied response was observed under high-level Cd. 705 differentially expressed genes (DEGs) were identified in P. cornucopia at Cd1 and Cd20, whereas 12,551 DEGs in P. ostreatus. Differentially regulated functional categories and pathways were also identified. ATP-binding cassette transporters were involved in Cd transport in P. cornucopia, whereas the endocytosis and phagosome pathways were more enhanced in P. ostreatus. 26 enzymes including peroxisomal enzymes catalase and superoxide dismutase were upregulated in P. ostreatus, whereas only cytosolic catalase was overexpressed in P. cornucopia, suggesting their different Cd detoxification pathways. Also, the mitogen-activated protein kinase signaling pathway involved in Cd resistance in both species instead of glutathione metabolism, although more active in P. ostreatus. These findings provided new insight into the molecular mechanism of mycoremediation and accumulator screening.

Transcriptional response of mushrooms to artificial sun exposure

Climate change causes increased tree mortality leading to canopy loss and thus sun-exposed forest floors. Sun exposure creates extreme temperatures and radiation, with potentially more drastic effects on forest organisms than the current increase in mean temperature. Such conditions might potentially negatively affect the maturation of mushrooms of forest fungi. A failure of reaching maturation would mean no sexual spore release and, thus, entail a loss of genetic diversity. However, we currently have a limited understanding of the quality and quantity of mushroom-specific molecular responses caused by sun exposure. Thus, to understand the short-term responses toward enhanced sun exposure, we exposed mushrooms of the wood-inhabiting forest species Lentinula edodes, while still attached to their mycelium and substrate, to artificial solar light (ca. 30°C and 100,000 lux) for 5, 30, and 60 min. We found significant differentially expressed genes at 30 and 60 min. Eukaryotic Orthologous Groups (KOG) class enrichment pointed to defense mechanisms. The 20 most significant differentially expressed genes showed the expression of heat-shock proteins, an important family of proteins under heat stress. Although preliminary, our results suggest mushroom-specific molecular responses to tolerate enhanced sun exposure as expected under climate change. Whether mushroom-specific molecular responses are able to maintain fungal fitness under opening forest canopies remains to be tested.

Comparative Transcriptome and Endophytic Bacterial Community Analysis of Morchella conica SH

The precious rare edible fungus Morchella conica is popular worldwide for its rich nutrition, savory flavor, and varieties of bioactive components. Due to its high commercial, nutritional, and medicinal value, it has always been a hot spot. However, the molecular mechanism and endophytic bacterial communities in M. conica were poorly understood. In this study, we sequenced, assembled, and analyzed the genome of M. conica SH. Transcriptome analysis reveals significant differences between the mycelia and fruiting body. As shown in this study, 1,329 and 2,796 genes were specifically expressed in the mycelia and fruiting body, respectively. The Gene Ontology (GO) enrichment showed that RNA polymerase II transcription activity-related genes were enriched in the mycelium-specific gene cluster, and nucleotide binding-related genes were enriched in the fruiting body-specific gene cluster. Further analysis of differentially expressed genes in different development stages resulted in finding two groups with distinct expression patterns. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment displays that glycan degradation and ABC transporters were enriched in the group 1 with low expressed level in the mycelia, while taurine and hypotaurine metabolismand tyrosine metabolism-related genes were significantly enriched in the group 2 with high expressed level in mycelia. Moreover, a dynamic shift of bacterial communities in the developing fruiting body was detected by 16S rRNA sequencing, and co-expression analysis suggested that bacterial communities might play an important role in regulating gene expression. Taken together, our study provided a better understanding of the molecular biology of M. conica SH and direction for future research on artificial cultivation.

Mycopharmaceuticals and Nutraceuticals: Promising Agents to Improve Human Well-Being and Life Quality

Fungi, especially edible mushrooms, are considered as high-quality food with nutritive and functional values. They are of considerable interest and have been used in the synthesis of nutraceutical supplements due to their medicinal properties and economic significance. Specific fungal groups, including predominantly filamentous endophytic fungi from Ascomycete phylum and several Basidiomycetes, produce secondary metabolites (SMs) with bioactive properties that are involved in the antimicrobial and antioxidant activities. These beneficial fungi, while high in protein and important fat contents, are also a great source of several minerals and vitamins, in particular B vitamins that play important roles in carbohydrate and fat metabolism and the maintenance of the nervous system. This review article will summarize and discuss the abilities of fungi to produce antioxidant, anticancer, antiobesity, and antidiabetic molecules while also reviewing the evidence from the last decade on the importance of research in fungi related products with direct and indirect impact on human health.

Transcriptome of different fruiting stages in the cultivated mushroom Cyclocybe aegerita suggests a complex regulation of fruiting and reveals enzymes putatively involved in fungal oxylipin biosynthesis

BACKGROUND

Cyclocybe aegerita (syn. Agrocybe aegerita) is a commercially cultivated mushroom. Its archetypal agaric morphology and its ability to undergo its whole life cycle under laboratory conditions makes this fungus a well-suited model for studying fruiting body (basidiome, basidiocarp) development. To elucidate the so far barely understood biosynthesis of fungal volatiles, alterations in the transcriptome during different developmental stages of C. aegerita were analyzed and combined with changes in the volatile profile during its different fruiting stages.

RESULTS

A transcriptomic study at seven points in time during fruiting body development of C. aegerita with seven mycelial and five fruiting body stages was conducted. Differential gene expression was observed for genes involved in fungal fruiting body formation showing interesting transcriptional patterns and correlations of these fruiting-related genes with the developmental stages. Combining transcriptome and volatilome data, enzymes putatively involved in the biosynthesis of C8 oxylipins in C. aegerita including lipoxygenases (LOXs), dioxygenases (DOXs), hydroperoxide lyases (HPLs), alcohol dehydrogenases (ADHs) and ene-reductases could be identified. Furthermore, we were able to localize the mycelium as the main source for sesquiterpenes predominant during sporulation in the headspace of C. aegerita cultures. In contrast, changes in the C8 profile detected in late stages of development are probably due to the activity of enzymes located in the fruiting bodies.

CONCLUSIONS

In this study, the combination of volatilome and transcriptome data of C. aegerita revealed interesting candidates both for functional genetics-based analysis of fruiting-related genes and for prospective enzyme characterization studies to further elucidate the so far barely understood biosynthesis of fungal C8 oxylipins.

Analyses of transcriptomes and the first complete genome of Leucocalocybe mongolica provide new insights into phylogenetic relationships and conservation

In this study, we report a de novo assembly of the first high-quality genome for a wild mushroom species Leucocalocybe mongolica (LM). We performed high-throughput transcriptome sequencing to analyze the genetic basis for the life history of LM. Our results show that the genome size of LM is 46.0 Mb, including 26 contigs with a contig N50 size of 3.6 Mb. In total, we predicted 11,599 protein-coding genes, of which 65.7% (7630) could be aligned with high confidence to annotated homologous genes in other species. We performed phylogenetic analyses using genes form 3269 single-copy gene families and showed support for distinguishing LM from the genus Tricholoma (L.) P.Kumm., in which it is sometimes circumscribed. We believe that one reason for limited wild occurrences of LM may be the loss of key metabolic genes, especially carbohydrate-active enzymes (CAZymes), based on comparisons with other closely related species. The results of our transcriptome analyses between vegetative (mycelia) and reproductive (fruiting bodies) organs indicated that changes in gene expression among some key CAZyme genes may help to determine the switch from asexual to sexual reproduction. Taken together, our genomic and transcriptome data for LM comprise a valuable resource for both understanding the evolutionary and life history of this species.

Genome-wide analysis and prediction of genes involved in the biosynthesis of polysaccharides and bioactive secondary metabolites in high-temperature-tolerant wild Flammulina filiformis

Background

Flammulina filiformis (previously known as Asian F. velutipes) is a popular commercial edible mushroom. Many bioactive compounds with medicinal effects, such as polysaccharides and sesquiterpenoids, have been isolated and identified from F. filiformis, but their biosynthesis and regulation at the molecular level remains unclear. In this study, we sequenced the genome of the wild strain F. filiformis Liu355, predicted its biosynthetic gene clusters (BGCs) and profiled the expression of these genes in wild and cultivar strains and in different developmental stages of the wild F. filiformis strain by a comparative transcriptomic analysis.

Results

We found that the genome of the F. filiformis was 35.01 Mb in length and harbored 10,396 gene models. Thirteen putative terpenoid gene clusters were predicted and 12 sesquiterpene synthase genes belonging to four different groups and two type I polyketide synthase gene clusters were identified in the F. filiformis genome. The number of genes related to terpenoid biosynthesis was higher in the wild strain (119 genes) than in the cultivar strain (81 genes). Most terpenoid biosynthesis genes were upregulated in the primordium and fruiting body of the wild strain, while the polyketide synthase genes were generally upregulated in the mycelium of the wild strain. Moreover, genes encoding UDP-glucose pyrophosphorylase and UDP-glucose dehydrogenase, which are involved in polysaccharide biosynthesis, had relatively high transcript levels both in the mycelium and fruiting body of the wild F. filiformis strain.

Conclusions

F. filiformis is enriched in a number of gene clusters involved in the biosynthesis of polysaccharides and terpenoid bioactive compounds and these genes usually display differential expression between wild and cultivar strains, even in different developmental stages. This study expands our knowledge of the biology of F. filiformis and provides valuable data for elucidating the regulation of secondary metabolites in this unique F. filiformis strain.

Untargeted Metabolite Profiling of Antimicrobial Compounds in the Brown Film of Lentinula edodes Mycelium via LC-MS/MS Analysis

The brown film (BF) of Lentinula edodes mycelium has been reported to exert biological activities during mushroom cultivation; however, to date, there is limited information on its chemical composition. In this study, untargeted metabolomics analysis was performed via liquid chromatography-mass spectrometry (LC-MS), and the results were used to screen the antimicrobial compounds. A total of 236 differential metabolites were found among the BF stages compared with the white hyphal stage. Among them, five important antimicrobial metabolites related to antimicrobial activities, namely, 6-deoxyerythronolide B, tanikolide, hydroxyanthraquinone, benzylideneacetone, and 9-OxooTrE, were present at high levels in the BF samples. The score plots of the principal component analysis indicated that the samples from four time points could be classified into two groups. This study provided a comprehensive profile of the antimicrobial compounds produced during BF formation and partly clarified the antibacterial and antifungal mechanism of the BF of L. edodes mycelium.

Transcriptome analysis of the edible mushroom Lentinula edodes in response to blue light

Lentinula edodes is one of the most popular edible mushrooms worldwide and contains important medicinal components such as lentinan, ergosterol, and eritadenine. Mushroom metabolism is regulated by the mycelia and fruit body using light; however, in mushrooms, the underlying molecular mechanisms controlling this process as well as light-induced gene expression remain unclear. Therefore, in this study, we compared morphological changes and gene expression in the fruit bodies of L. edodes cultivated under blue light and continuous darkness. Our results showed that blue light primarily induced pileus growth (diameter and thickness) compared to dark cultivation. Alternatively, stipe length development was promoted by dark cultivation. We also performed RNAseq on L. edodes under the blue light/dark cultivation conditions. A total of 12,051 genes were used for aligning the Illumina raw reads and 762 genes that showed fold change cut-offs of >|2| and significance p-values of <0.05 were selected under blue light condition. Among the genes which showed two-fold changed genes, 221 were upregulated and 541 were downregulated. In order to identify blue light induced candidate genes, differentially expressed genes (DEGs) were selected according to 4-fold changes and validated by RT-PCR. We identified 8 upregulated genes under blue light condition, such as DDR48-heat shock protein, Fasciclin-domain-containing protein and carbohydrate esterase family 4 protein, FAD NAD-binding domain-containing protein that are involved in morphological development of primordium and embryonic muscle development, cell adhesion and affect the structure of cellulosic and non-cellulosic cell walls of fruit body development, and photoreceptor of blue light signaling for fruit body and pigment development, respectively. This study provides valuable insights into the molecular mechanisms underlying the role of blue light in mushroom growth and development and can thus contribute to breeding programs to improve mushroom cultivation.

The PoV mycovirus affects extracellular enzyme expression and fruiting body yield in the oyster mushroom, Pleurotus ostreatus

Isogenic virus-cured and virus-infected fungal strains were previously obtained and compared to investigate mycoviral diseases and, specifically, the influence of viral infection on the vegetative growth of Pleurotus ostreatus. The present study demonstrated that infection with mycovirus PoV-ASI2792 (PoV) caused phenotypic and physiological changes in fungal cells and mycelia. The microscopically determined growth rate of the virus-infected strain was lower than that of the virus-cured strain, due to the conglomerate phenomenon during the mycelial growth process. An exploration of the viral effects of PoV on fruiting bodies yield showed significantly lower than that on virus-cured P. ostreatus. A colorimetric assay of polyphenol oxidase activity in the strains showed very weak activity in the virus-infected strain. To estimate the activity levels of enzymes related to the growth and fruiting body formation, the relative expression levels of genes encoding various extracellular enzymes such as Carbohydrate-Active Enzymes (CAZymes) were measured by quantitative RT-PCR. The expression levels of the assayed genes were significantly lower in virus-infected than in virus-cured P. ostreatus. Together, these results indicate that PoV infection affects the spawn growth and fruiting body formation of P. ostreatus via decreased expression and activity of some extracellular enzymes including lignocellulolytic enzymes.

Integration of ATAC-Seq and RNA-Seq Identifies Key Genes in Light-Induced Primordia Formation of Sparassis latifolia

Light is an essential environmental factor for Sparassis latifolia primordia formation, but the molecular mechanism is still unclear. In this study, differential expression profiling of light-induced primordia formation (LIPF) was established by integrating the assay for transposase accessible chromatin by sequencing (ATAC-seq) and RNA-seq technology. The integrated results from the ATAC-seq and RNA-seq showed 13 down-regulated genes and 17 up-regulated genes in both the L vs. D and P vs. D groups, for both methods. According to the gene ontology (GO) annotation of these differentially expressed genes (DEGs), the top three biological process categories were cysteine biosynthetic process via cystathionine, vitamin B6 catabolic, and glycine metabolic; the top three molecular function categories were 5-methyltetrahydropteroyltriglutamate-homocysteine S-methyltransferase activity, glycine binding, and pyridoxal phosphate binding; cellular component categories were significantly enriched in the glycine cleavage complex. The KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis revealed that these genes were associated with vitamin B6 metabolism; selenocompound metabolism; cysteine and methionine metabolism; glycine, serine, and threonine metabolism; and glyoxylate and dicarboxylate metabolism pathways. The expression of most of the DEGs was validated by qRT-PCR. To the best of our knowledge, this study is the first integrative analysis of ATAC-seq and RNA-seq for macro-fungi. These results provided a new perspective on the understanding of key pathways and hub genes in LIPF in S. latifolia. It will be helpful in understanding the primary environmental response, and provides new information to the existing models of primordia formation in edible and medicinal fungi.

Comparative genomics reveals unique wood-decay strategies and fruiting body development in the Schizophyllaceae

Agaricomycetes are fruiting body-forming fungi that produce some of the most efficient enzyme systems to degrade wood. Despite decades-long interest in their biology, the evolution and functional diversity of both wood-decay and fruiting body formation are incompletely known. We performed comparative genomic and transcriptomic analyses of wood-decay and fruiting body development in Auriculariopsis ampla and Schizophyllum commune (Schizophyllaceae), species with secondarily simplified morphologies, an enigmatic wood-decay strategy and weak pathogenicity to woody plants. The plant cell wall-degrading enzyme repertoires of Schizophyllaceae are transitional between those of white rot species and less efficient wood-degraders such as brown rot or mycorrhizal fungi. Rich repertoires of suberinase and tannase genes were found in both species, with tannases restricted to Agaricomycetes that preferentially colonize bark-covered wood, suggesting potential complementation of their weaker wood-decaying abilities and adaptations to wood colonization through the bark. Fruiting body transcriptomes revealed a high rate of divergence in developmental gene expression, but also several genes with conserved expression patterns, including novel transcription factors and small-secreted proteins, some of the latter which might represent fruiting body effectors. Taken together, our analyses highlighted novel aspects of wood-decay and fruiting body development in an important family of mushroom-forming fungi.

Effects of Medium Composition and Genetic Background on Agrobacterium-Mediated Transformation Efficiency of Lentinula edodes

The establishment of genetic transformation method is crucial for the functional genomics research in filamentous fungi. Although the transformation method has been developed in several types of fungi, a highly efficient and convenient transformation system is desperately needed in Lentinula edodes. Present work established the Agrobacterium-mediated transformation (ATMT) of basidiomycete L. edodes in both monokaryon and dikaryon mycelia by using constructed binary plasmid pCAMBIA-1300-GFP. Then, the transformation efficiency of ATMT was evaluated by using different mediums for recipient incubation and different varieties of L. edodes. The results showed that in dikaryon strain W1, the positive hygromycin-resistant transformants was observed in all medium with the positive frequency of selected transformants that ranged from 0 to 30%. While in the monokaryon strain W1-26, only the millet medium group obtained positive transformants with a positive frequency of 75.48%. Moreover, three dikaryotic wild strains (YS55, YS3334, and YS3357) and two dikaryotic cultivated strains (W1 and S606) showed the highest transformation efficiency, with 32.96% of the germination frequency, and 85.12% of positive frequency for hygromycin-resistant transformants. This work demonstrated that Agrobacterium-mediated transformation was successfully performed in L. edodes, and the genotype of recipients as well as the medium for mycelial incubation were suggested to play key roles in determining the transformation efficiency. These findings may provide new avenues for the genetic modification of edible mushroom and may extend the cognition of DNA-mediated transformation in filamentous fungi.

Transcriptome Changes during Major Developmental Transitions Accompanied with Little Alteration of DNA Methylome in Two Pleurotus Species

Pleurotus tuoliensis (Pt) and P. eryngii var. eryngii (Pe) are important edible mushrooms. The epigenetic and gene expression signatures characterizing major developmental transitions in these two mushrooms remain largely unknown. Here, we report global analyses of DNA methylation and gene expression in both mushrooms across three major developmental transitions, from mycelium to primordium and to fruit body, by whole-genome bisulfite sequencing (WGBS) and RNA-seq-based transcriptome profiling. Our results revealed that in both Pt and Pe the landscapes of methylome are largely stable irrespective of genomic features, e.g., in both protein-coding genes and transposable elements (TEs), across the developmental transitions. The repressive impact of DNA methylation on expression of a small subset of genes is likely due to TE-associated effects rather than their own developmental dynamics. Global expression of gene orthologs was also broadly conserved between Pt and Pe, but discernible interspecific differences exist especially at the fruit body formation stage, and which are primarily due to differences in trans-acting factors. The methylome and transcriptome repertories we established for the two mushroom species may facilitate further studies of the epigenetic and transcriptional regulatory mechanisms underpinning gene during development in Pleurotus and related genera.

High-throughput targeted gene deletion in the model mushroom Schizophyllum commune using pre-assembled Cas9 ribonucleoproteins

Efficient gene deletion methods are essential for the high-throughput study of gene function. Compared to most ascomycete model systems, gene deletion is more laborious in mushroom-forming basidiomycetes due to the relatively low incidence of homologous recombination (HR) and relatively high incidence of non-homologous end-joining (NHEJ). Here, we describe the use of pre-assembled Cas9-sgRNA ribonucleoproteins (RNPs) to efficiently delete the homeodomain transcription factor gene hom2 in the mushroom-forming basidiomycete Schizophyllum commune by replacing it with a selectable marker. All components (Cas9 protein, sgRNA, and repair template with selectable marker) were supplied to wild type protoplasts by PEG-mediated transformation, abolishing the need to optimize the expression of cas9 and sgRNAs. A Δku80 background further increased the efficiency of gene deletion. A repair template with homology arms of 250 bp was sufficient to efficiently induce homologous recombination. This is the first report of the use of pre-assembled Cas9 RNPs in a mushroom-forming basidiomycete and this approach may also improve the genetic accessibility of non-model species.

Alternative metabolic routes in channeling xylose to cordycepin production of Cordyceps militaris identified by comparative transcriptome analysis

The responsive mechanism of C. militaris TBRC7358 on xylose utilization was investigated by comparative analysis of transcriptomes, growth kinetics and cordycepin productions. The result showed that the culture grown on xylose exhibited high production yield of cordycepin on dry biomass. Comparing xylose to other carbon sources, a set of significantly up-regulated genes in xylose were enriched in pentose and glucuronate interconversion, and cordycepin biosynthesis. After validating up-regulated genes using quantitative real-time PCR, interestingly, putative alternative 3'-AMP-associated metabolic route on cordycepin biosynthesis was identified. Through reporter metabolites analysis of C. militaris, significant metabolites (e.g., AMP, glycine and L-glutamate) were identified guiding involvement of growth and cordycepin production. These findings suggested that there was a cooperative mechanism in transcriptional control of the supplying precursors pool directed towards the cordycepin biosynthesis through main and putative alternative metabolic routes for leverage of cell growth and cordycepin production on xylose of C. militaris strain TBRC7358.

A Novel Rapid Fungal Promoter Analysis System Using the Phosphopantetheinyl Transferase Gene, npgA, in Aspergillus nidulans

To develop a convenient promoter analysis system for fungi, a null-pigment mutant (NPG) of Aspergillus nidulans was used with the 4'-phosphopantetheinyl transferase (PPTase) gene, npgA, which restores the normal pigmentation in A. nidulans, as a new reporter gene. The functional organization of serially deleted promoter regions of the A. nidulans trpC gene and the Cryphonectria parasitica crp gene in filamentous fungi was representatively investigated to establish a novel fungal promoter assay system that depends on color complementation of the NPG mutant with the PPTase npgA gene. Several promoter regions of the trpC and crp genes were fused to the npgA gene containing the 1,034-bp open reading frame and the 966-bp 3' downstream region from the TAA, and the constructed fusions were introduced into the NPG mutant in A. nidulans to evaluate color recovery due to the transcriptional activity of the sequence elements. Serial deletion of the trpC and crp promoter regions in this PPTase reporter assay system reaffirmed results in previous reports by using the fungal transformation step without a laborious verification process. This approach suggests a more rapid and convenient system than conventional analyses for fungal gene expression studies.