Transcription factor genes of Schizophyllum commune involved in regulation of mushroom formation

High phenotypic and genotypic plasticity among strains of the mushroom-forming fungus Schizophyllum commune

Schizophyllum commune is a mushroom-forming fungus notable for its distinctive fruiting bodies with split gills. It is used as a model organism to study mushroom development, lignocellulose degradation and mating type loci. It is a hypervariable species with considerable genetic and phenotypic diversity between the strains. In this study, we systematically phenotyped 16 dikaryotic strains for aspects of mushroom development and 18 monokaryotic strains for lignocellulose degradation. There was considerable heterogeneity among the strains regarding these phenotypes. The majority of the strains developed mushrooms with varying morphologies, although some strains only grew vegetatively under the tested conditions. Growth on various carbon sources showed strain-specific profiles. The genomes of seven monokaryotic strains were sequenced and analyzed together with six previously published genome sequences. Moreover, the related species Schizophyllum fasciatum was sequenced. Although there was considerable genetic variation between the genome assemblies, the genes related to mushroom formation and lignocellulose degradation were well conserved. These sequenced genomes, in combination with the high phenotypic diversity, will provide a solid basis for functional genomics analyses of the strains of S. commune.

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.

Regulatory functions of homeobox domain transcription factors in fungi

Homeobox domain (HD) proteins present a crucial involvement in morphological differentiation and other functions in eukaryotes. Most HD genes encode transcription factors (TFs) that orchestrate a regulatory role in cellular and developmental decisions. In fungi, multiple studies have increased our understanding of these important HD regulators in recent years. These reports have revealed their role in fungal development, both sexual and asexual, as well as their importance in governing other biological processes in these organisms, including secondary metabolism, pathogenicity, and sensitivity to environmental stresses. Here, we provide a comprehensive review of the current knowledge on the regulatory roles of HD-TFs in fungi, with a special focus on Aspergillus species.

Snowball: a novel gene family required for developmental patterning of fruiting bodies of mushroom-forming fungi (Agaricomycetes)

The morphogenesis of sexual fruiting bodies of fungi is a complex process determined by a genetically encoded program. Fruiting bodies reached the highest complexity levels in the Agaricomycetes; yet, the underlying genetics is currently poorly known. In this work, we functionally characterized a highly conserved gene termed snb1, whose expression level increases rapidly during fruiting body initiation. According to phylogenetic analyses, orthologs of snb1 are present in almost all agaricomycetes and may represent a novel conserved gene family that plays a substantial role in fruiting body development. We disrupted snb1 using CRISPR/Cas9 in the agaricomycete model organism Coprinopsis cinerea. snb1 deletion mutants formed unique, snowball-shaped, rudimentary fruiting bodies that could not differentiate caps, stipes, and lamellae. We took advantage of this phenotype to study fruiting body differentiation using RNA-Seq analyses. This revealed differentially regulated genes and gene families that, based on wild-type RNA-Seq data, were upregulated early during development and showed tissue-specific expression, suggesting a potential role in differentiation. Taken together, the novel gene family of snb1 and the differentially expressed genes in the snb1 mutants provide valuable insights into the complex mechanisms underlying developmental patterning in the Agaricomycetes.

IMPORTANCE

Fruiting bodies of mushroom-forming fungi (Agaricomycetes) are complex multicellular structures, with a spatially and temporally integrated developmental program that is, however, currently poorly known. In this study, we present a novel, conserved gene family, Snowball (snb), termed after the unique, differentiation-less fruiting body morphology of snb1 knockout strains in the model mushroom Coprinopsis cinerea. snb is a gene of unknown function that is highly conserved among agaricomycetes and encodes a protein of unknown function. A comparative transcriptomic analysis of the early developmental stages of differentiated wild-type and non-differentiated mutant fruiting bodies revealed conserved differentially expressed genes which may be related to tissue differentiation and developmental patterning fruiting body development.

Pleurotus ostreatus as a model mushroom in genetics, cell biology, and material sciences

Pleurotus ostreatus, also known as the oyster mushroom, is a popular edible mushroom cultivated worldwide. This review aims to survey recent progress in the molecular genetics of this fungus and demonstrate its potential as a model mushroom for future research. The development of modern molecular genetic techniques and genome sequencing technologies has resulted in breakthroughs in mushroom science. With efficient transformation protocols and multiple selection markers, a powerful toolbox, including techniques such as gene knockout and genome editing, has been developed, and numerous new findings are accumulating in P. ostreatus. These include molecular mechanisms of wood component degradation, sexual development, protein secretion systems, and cell wall structure. Furthermore, these techniques enable the identification of new horizons in enzymology, biochemistry, cell biology, and material science through protein engineering, fluorescence microscopy, and molecular breeding. KEY POINTS: • Various genetic techniques are available in Pleurotus ostreatus. • P. ostreatus can be used as an alternative model mushroom in genetic analyses. • New frontiers in mushroom science are being developed using the fungus.

Anaerobic production and biosynthesis mechanism of exopolysaccharides in Schizophyllum commune 20R-7-F01

Exopolysaccharides (EPS) produced by microorganisms play a vital role in physiological and ecological processes. However, the mechanisms of EPS synthesis and release in anaerobic environments remain poorly understood. Here, we provide the first evidence of anaerobic EPS synthesis by the fungus Schizophyllum commune 20R-7-F01, isolated from coal-bearing sediments ~2.0 km below the seafloor. Under anaerobic conditions, the fungus exhibited significantly higher specific EPS production (1.57 times) than under aerobic conditions. Transcriptomic analysis revealed 2057 differentially expressed genes (DEGs) in the strain cultured anaerobically for 7 days compared to aerobically. Among these genes, 642 were significantly upregulated, while 1415 were significantly downregulated, mainly associated with carbon metabolism pathways. Genes involved in glycolysis and EPS synthesis, including hexokinase (HK), phosphoglucomutase (PGM), and (1 → 3)-β-glucan synthase (GLS), were significantly upregulated, while those related to the TCA cycle, respiratory chain, and pentose phosphate pathway were downregulated under anaerobic conditions. These findings highlight the oxygen-dependent regulation of EPS synthesis and suggest that EPS may serve as a key mechanism for fungal adaptation to anaerobic environments.

Integration of Metabolomes and Transcriptomes Provides Insights into Morphogenesis and Maturation in Morchella sextelata

True morels (Morchella, Pezizales) are a popular edible and medicinal fungus with great nutritional and economic value. The dynamics and regulatory mechanisms during the morphogenesis and maturation of morels are poorly understood. In this study, the metabolomes and transcriptomes of the mycelium (MY), primordium differentiation (PR), young fruiting body (YFB), and mature fruiting body (MFB) were comprehensively analyzed to reveal the mechanism of the morphogenesis and maturation of Morchella sextelata. A total of 748 differentially expressed metabolites (DEMs) and 5342 differentially expressed genes (DEGs) were detected, mainly enriched in the carbohydrate, amino acid, and lipid metabolism pathways, with the transition from the mycelium to the primordium being the most drastic stage at both the metabolic and transcriptional levels. The integrated metabolomics and transcriptomics highlighted significant correlations between the DEMs and DEGs, and specific amino acid and nucleotide metabolic pathways were significantly co-enriched, which may play key roles in morphological development and ascocarp maturation. A conceptual model of transcriptional and metabolic regulation was proposed during morphogenesis and maturation in M. sextelata for the first time, in which environmental factors activate the regulation of transcription factors, which then promote metabolic and transcriptional regulation from vegetative to reproductive growth. These results provide insights into the metabolic dynamics and transcriptional regulation during the morphogenesis and maturation of morels and valuable resources for future breeding enhancement and sustainable artificial cultivation.

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.

Transcription factors: switches for regulating growth and development in macrofungi

Macrofungi (or mushrooms) act as an extraordinarily important part to human health due to their nutritional and/or medicinal value, but the detailed researches in growth and development mechanisms have yet to be explored further. Transcription factors (TFs) play indispensable roles in signal transduction and affect growth, development, and metabolism of macrofungi. In recent years, increasing research effort has been employed to probe the relationship between the development of macrofungi and TFs. Herein, the present review comprehensively summarized the functional TFs researched in macrofungi, including modulating mycelial growth, fructification, sclerotial formation, sexual reproduction, spore formation, and secondary metabolism. Meanwhile, the possible effect mechanisms of TFs on the growth and development of some macrofungi were also revealed. Specific examples of functional characterizations of TFs in macrofungi (such as Schizophyllum commune and Coprinopsis cinerea) were described to a better comprehension of regulatory effect. Future research prospects in the field of TFs of macrofungi are discussed. We illustrated the functional versatility of the TFs in macrofungi based on specific examples. A systematical realization of the interaction and possible mechanisms between TFs and macrofungi can supply possible solutions to regulate genetic characteristics, which supply novel insights into the regulation of growth, development and metabolism of macrofungi. KEY POINTS: • The functional TFs researched in macrofungi were summarized. • The possible effect mechanisms of TFs in macrofungal were described. • The multiple physiological functions of TFs in macrofungi were discussed.

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.

Hydrophobin Gene Cmhyd4 Negatively Regulates Fruiting Body Development in Edible Fungi Cordyceps militaris

A deep understanding of the mechanism of fruiting body development is important for mushroom breeding and cultivation. Hydrophobins, small proteins exclusively secreted by fungi, have been proven to regulate the fruiting body development in many macro fungi. In this study, the hydrophobin gene Cmhyd4 was revealed to negatively regulate the fruiting body development in Cordyceps militaris, a famous edible and medicinal mushroom. Neither the overexpression nor the deletion of Cmhyd4 affected the mycelial growth rate, the hydrophobicity of the mycelia and conidia, or the conidial virulence on silkworm pupae. There was also no difference between the micromorphology of the hyphae and conidia in WT and ΔCmhyd4 strains observed by SEM. However, the ΔCmhyd4 strain showed thicker aerial mycelia in darkness and quicker growth rates under abiotic stress than the WT strain. The deletion of Cmhyd4 could promote conidia production and increase the contents of carotenoid and adenosine. The biological efficiency of the fruiting body was remarkably increased in the ΔCmhyd4 strain compared with the WT strain by improving the fruiting body density, not the height. It was indicated that Cmhyd4 played a negative role in fruiting body development. These results revealed that the diverse negative roles and regulatory effects of Cmhyd4 were totally different from those of Cmhyd1 in C. militaris and provided insights into the developmental regulatory mechanism of C. militaris and candidate genes for C. militaris strain breeding.

The transcription factor Ste12-like increases the mycelial abiotic stress tolerance and regulates the fruiting body development of Flammulina filiformis

Introduction

Flammulina filiformis is one of the most commercially important edible fungi worldwide, with its nutritional value and medicinal properties. It becomes a good model species to study the tolerance of abiotic stress during mycelia growth in edible mushroom cultivation. Transcription factor Ste12 has been reported to be involved in the regulation of stress tolerance and sexual reproduction in fungi.

Methods

In this study, identification and phylogenetic analysis of ste12-like was performed by bioinformatics methods. Four ste12-like overexpression transformants of F. filiformis were constructed by Agrobacterium tumefaciens-mediated transformation.

Results and Discussion

Phylogenetic analysis showed that Ste12-like contained conserved amino acid sequences. All the overexpression transformants were more tolerant to salt stress, cold stress and oxidative stress than wild-type strains. In the fruiting experiment, the number of fruiting bodies of overexpression transformants increased compared with wild-type strains, but the growth rate of stipes slowed down. It suggested that gene ste12-like was involved in the regulation of abiotic stress tolerance and fruiting body development in F. filiformis.

The Role of Chromatin and Transcriptional Control in the Formation of Sexual Fruiting Bodies in Fungi

Fungal fruiting bodies are complex, three-dimensional structures that arise from a less complex vegetative mycelium. Their formation requires the coordinated action of many genes and their gene products, and fruiting body formation is accompanied by major changes in the transcriptome. In recent years, numerous transcription factor genes as well as chromatin modifier genes that play a role in fruiting body morphogenesis were identified, and through research on several model organisms, the underlying regulatory networks that integrate chromatin structure, gene expression, and cell differentiation are becoming clearer. This review gives a summary of the current state of research on the role of transcriptional control and chromatin structure in fruiting body development. In the first part, insights from transcriptomics analyses are described, with a focus on comparative transcriptomics. In the second part, examples of more detailed functional characterizations of the role of chromatin modifiers and/or transcription factors in several model organisms (Neurospora crassa, Aspergillus nidulans, Sordaria macrospora, Coprinopsis cinerea, and Schizophyllum commune) that have led to a better understanding of regulatory networks at the level of chromatin structure and transcription are discussed.

Transcriptome profiling of transcription factors in Ganoderma lucidum in response to methyl jasmonate

Ganoderma lucidum is a traditional Chinese medicine and its major active ingredients are ganoderma triterpenoids (GTs). To screen for transcription factors (TFs) that involved in the biosynthetic pathway of GTs in G. lucidum, the chemical composition in mycelia, primordium and fruiting body were analyzed, and the transcriptomes of mycelia induced by methyl jasmonate (MeJA) were analyzed. In addition, the expression level data of MeJA-responsive TFs in mycelia, primordia and fruiting body were downloaded from the database, and the correlation analysis was carried out between their expression profiles and the content of total triterpenoids. The results showed that a total of 89 components were identified, and the content of total triterpenoids was the highest in primordium, followed by fruiting body and mycelia. There were 103 differentially expressed TFs that response to MeJA-induction including 95 upregulated and 8 downregulated genes. These TFs were classified into 22 families including C2H2 (15), TFII-related (12), HTH (9), fungal (8), bZIP (6), HMG (5), DADS (2), etc. Correlation analysis showed that the expression level of GL23559 (MADS), GL26472 (HTH), and GL31187 (HMG) showed a positive correlation with the GTs content, respectively. While the expression level of GL25628 (fungal) and GL26980 (PHD) showed a negative correlation with the GTs content, respectively. Furthermore, the over expression of the Glmhr1 gene (GL25628) in Pichia pastoris GS115 indicated that it might be a negative regulator of GT biosynthesis through decreasing the production of lanosterol. This study provided useful information for a better understanding of the regulation of TFs involved in GT biosynthesis and fungal growth in G. lucidum.

Comparative transcriptome and WGCNA reveal key genes involved in lignocellulose degradation in Sarcomyxa edulis

The developmental transcriptomes of Sarcomyxa edulis were assessed to explore the molecular mechanisms underlying lignocellulose degradation. Six stages were analyzed, spanning the entire developmental process: growth of mycelium until occupying half the bag (B1), mycelium under low-temperature stimulation after occupying the entire bag (B2), appearance of mycelium in primordia (B3), primordia (B4), mycelium at the harvest stage (B5), and mature fruiting body (B6). Samples from all six developmental stages were used for transcriptome sequencing, with three biological replicates for all experiments. A co-expression network of weighted genes associated with extracellular enzyme physiological traits was constructed using weighted gene co-expression network analysis (WGCNA). We obtained 19 gene co-expression modules significantly associated with lignocellulose degradation. In addition, 12 key genes and 8 kinds of TF families involved in lignocellulose degradation pathways were discovered from the four modules that exhibited the highest correlation with the target traits. These results provide new insights that advance our understanding of the molecular genetic mechanisms of lignocellulose degradation in S. edulis to facilitate its utilization by the edible mushroom industry.

Identification of Long Non-Coding RNAs and Their Target Genes from Mycelium and Primordium in Model Mushroom Schizophyllum commune

Schizophyllum commune has emerged as the most promising model mushroom to study developmental stages (mycelium, primordium), which are two primary processes of fruit body development. Long non-coding RNA (lncRNA) has been proved to participate in fruit development and sex differentiation in fungi. However, potential lncRNAs have not been identified in S. commune from mycelium to primordium developmental stages. In this study, lncRNA-seq was performed in S. commune and 61.56 Gb clean data were generated from mycelium and primordium developmental stages. Furthermore, 191 lncRNAs had been obtained and a total of 49 lncRNAs were classified as differently expressed lncRNAs. Additionally, 26 up-regulated differently expressed lncRNAs and 23 down-regulated between mycelium and primordia libraries were detected. Further, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that differentially expressed lncRNAs target genes from the MAPK pathway, phosphatidylinositol signal, ubiquitin-mediated proteolysis, autophagy, and cell cycle. This study provides a new resource for further research on the relationship between lncRNA and two developmental stages (mycelium, primordium) in S. commune.

Whole-Genome Sequencing and Transcriptome Analysis of Ganoderma lucidum Strain Yw-1-5 Provides New Insights into the Enhanced Effect of Tween80 on Exopolysaccharide Production

Ganoderma lucidum is an important medicinal mushroom widely cultured in Asian countries. Exopolysaccharides are bioactive compounds of G. lucidum with health benefits. Limited exopolysaccharide content hinders its extraction from G. lucidum. The addition of Tween80 had an enhanced effect on G. lucidum exopolysaccharide production in submerged fermentation. However, the mechanism of this effect remains unclear. In this study, we report on a high-quality assembly of G. lucidum strain yw-1-5 to lay the foundation for further transcriptome analysis. The genome sequence was 58.16 Mb and consisted of 58 scaffolds with an N50 of 4.78 Mb. A total of 13,957 protein-coding genes were annotated and Hi-C data mapped to 12 pseudo-chromosomes. Genes encoding glycosyltransferases and glycoside hydrolases were also obtained. Furthermore, RNA-seq was performed in a Tween80-treated group and control group for revealing the enhanced effect of Tween80 on exopolysaccharide production. In total, 655 genes were identified as differentially expressed, including 341 up-regulated and 314 down-regulated. Further analysis of differentially expressed genes showed that groups of MAPK, amino sugar and nucleotide sugar metabolism, autophagy, ubiquitin-mediated proteolysis, peroxisome, starch and sucrose metabolism, TCA cycle, glycolysis/gluconeogenesis KEGG pathway, glycosyltransferases and glycoside hydrolases played important roles in the enhanced effect of Tween80 on exopolysaccharide production. This work provides a valuable resource for facilitating our understanding of the synthesis of polysaccharides and accelerating the breeding of new strains with a high content of exopolysaccharides.

De Novo Assembly Transcriptome Analysis Reveals the Preliminary Molecular Mechanism of Primordium Formation in Pleurotus tuoliensis

Primordium formation is extremely important for yield of Pleurotus tuoliensis. However, the molecular mechanism underlying primordium formation is largely unknown. This study investigated the transcriptional properties during primordium formation of P. tuoliensis by comparing transcriptome. Clean reads were assembled into 57,075 transcripts and 6874 unigenes. A total of 1397 differentially expressed genes were identified (26 DEGs altered in all stages). GO and KEGG enrichment analysis showed that these DEGs were involved in “oxidoreductase activity”, “glycolysis/gluconeogenesis”, “MAPK signaling pathways”, and “ribosomes”. Our results support further understanding of the transcriptional changes and molecular processes underlying primordium formation and differentiation of P. tuoliensis.

Gene targeting of dikaryotic Pleurotus ostreatus nuclei using the CRISPR/Cas9 system

Clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated protein 9 (Cas9)-assisted gene targeting is a promising method used in molecular breeding. We recently reported the successful introduction of this method in the monokaryotic Pleurotus ostreatus (oyster mushroom), PC9. However, considering their application in mushroom breeding, dikaryotic strains (with targeted gene mutations in both nuclei) need to be generated. This is laborious and time-consuming because a classical crossing technique is used. Herein, we report a technique that targets both nuclei of dikaryotic P. ostreatus, PC9×#64 in a transformation experiment using plasmid-based CRISPR/Cas9, with the aim of developing a method for efficient and rapid molecular breeding. As an example, we targeted strains with low basidiospore production ability through the meiosis-related genes mer3 or msh4. Four different plasmids containing expression cassettes for Cas9 and two different gRNAs targeting mer3 or msh4 were constructed and separately introduced into PC9×#64. Eight of the 38 dikaryotic transformants analyzed produced no basidiospores. Genomic PCR suggested that msh4 or mer3 mutations were introduced into both nuclei of seven out of eight strains. Thus, in this study, we demonstrated simultaneous gene targeting using our CRISPR/Cas9 system, which may be useful for the molecular breeding of cultivated agaricomycetes.

Transcriptome and metabolome analyses reveal transcription factors regulating ganoderic acid biosynthesis in Ganoderma lucidum development

Ganoderma lucidum is an important medicinal fungus in Asian countries. Ganoderic acid (GA) is the major variety of bioactive and medicative components in G. lucidum. Biosynthesis of secondary metabolites is usually associated with cell differentiation and development. However, the mechanism underlying these phenomena remain unclear. Transcription factors play an essential regulatory role in the signal transduction pathway, owing to the fact that they represent the major link between signal transduction and expression of target genes. In the present study, we performed transcriptome and metabolome analyses to identify transcription factors involved in GA biosynthesis during development of G. lucidum. Transcriptome data revealed differentially expressed genes between mycelia and primordia, as well as between mycelia and the fruiting body. Results from gene ontology enrichment analysis and metabolome analyses suggested that GAs and flavonoids biosynthetic process significantly changed during fungal development. The analysis of predicted occurrences of DNA-binding domains revealed a set of 53 potential transcription factor families in G. lucidum. Notably, we found homeobox transcription factor and velvet family protein played important role in GA biosynthesis. Combined with previous studies, we provided a model diagram of transcription factors involved in GA biosynthesis during fruiting body formation. Collectively, these results are expected to enhance our understanding into the mechanisms underlying secondary metabolite biosynthesis and development in fungi.

The Transcription Factor Roc1 Is a Key Regulator of Cellulose Degradation in the Wood-Decaying Mushroom Schizophyllum commune

Wood-degrading fungi in the phylum Basidiomycota play a crucial role in nutrient recycling by breaking down all components of wood. Fungi have evolved transcriptional networks that regulate expression of wood-degrading enzymes, allowing them to prioritize one nutrient source over another.

Sexual Crossing, Chromosome-Level Genome Sequences, and Comparative Genomic Analyses for the Medicinal Mushroom Taiwanofungus Camphoratus (Syn. Antrodia Cinnamomea, Antrodia Camphorata)

Taiwanofungus camphoratus mushrooms are a complementary and alternative medicine for hangovers, cancer, hypertension, obesity, diabetes, and inflammation. Though Taiwanofungus camphoratus has attracted considerable biotechnological and pharmacological attention, neither classical genetic nor genomic approaches have been properly established for it. We isolated four sexually competent monokaryons from two T. camphoratus dikaryons used for the commercial cultivation of orange-red (HC1) and milky-white (SN1) mushrooms, respectively. We also sequenced, annotated, and comparatively analyzed high-quality and chromosome-level genome sequences of these four monokaryons. These genomic resources represent a valuable basis for understanding the biology, evolution, and secondary metabolite biosynthesis of this economically important mushrooms. We demonstrate that T. camphoratus has a tetrapolar mating system and that HC1 and SN1 represent two intraspecies isolates displaying karyotypic variation. Compared with several edible mushroom model organisms, T. camphoratus underwent a significant contraction in the gene family and individual gene numbers, most notably for plant, fungal, and bacterial cell-wall-degrading enzymes, explaining why T. camphoratus mushrooms are rare in natural environments, are difficult and time-consuming to artificially cultivate, and are susceptible to fungal and bacterial infections. Our results lay the foundation for an in-depth T. camphoratus study, including precise genetic manipulation, improvements to mushroom fruiting, and synthetic biology applications for producing natural medicinal products. IMPORTANCETaiwanofungus camphoratus (Tc) is a basidiomycete fungus that causes brown heart rot of the aromatic tree Cinnamomum kanehirae. The Tc fruiting bodies have been used to treat hangovers, abdominal pain, diarrhea, hypertension, and other diseases first by aboriginal Taiwanese and later by people in many countries. To establish classical genetic and genomic approaches for this economically important medicinal mushroom, we first isolated and characterized four sexually competent monokaryons from two dikaryons wildly used for commercial production of Tc mushrooms. We applied PacBio single molecule, real-time sequencing technology to determine the near-completed genome sequences of four monokaryons. These telomere-to-telomere and gapless haploid genome sequences reveal all genomic variants needed to be studied and discovered, including centromeres, telomeres, retrotransposons, mating type loci, biosynthetic, and metabolic gene clusters. Substantial interspecies diversities are also discovered between Tc and several other mushroom model organisms, including Agrocybe aegerita, Coprinopsis cinerea, and Schizophyllum commune, and Ganoderma lucidum.

Gene age shapes the transcriptional landscape of sexual morphogenesis in mushroom forming fungi (Agaricomycetes)

Multicellularity has been one of the most important innovations in the history of life. The role of gene regulatory changes in driving transitions to multicellularity is being increasingly recognized; however, factors influencing gene expression patterns are poorly known in many clades. Here, we compared the developmental transcriptomes of complex multicellular fruiting bodies of eight Agaricomycetes and Cryptococcus neoformans, a closely related human pathogen with a simple morphology. In-depth analysis in Pleurotus ostreatus revealed that allele-specific expression, natural antisense transcripts, and developmental gene expression, but not RNA editing or a ‘developmental hourglass,’ act in concert to shape its transcriptome during fruiting body development. We found that transcriptional patterns of genes strongly depend on their evolutionary ages. Young genes showed more developmental and allele-specific expression variation, possibly because of weaker evolutionary constraint, suggestive of nonadaptive expression variance in fruiting bodies. These results prompted us to define a set of conserved genes specifically regulated only during complex morphogenesis by excluding young genes and accounting for deeply conserved ones shared with species showing simple sexual development. Analysis of the resulting gene set revealed evolutionary and functional associations with complex multicellularity, which allowed us to speculate they are involved in complex multicellular morphogenesis of mushroom fruiting bodies.

Genome-wide identification and expression analyses of C2H2 zinc finger transcription factors in Pleurotus ostreatus

The C2H2-type zinc finger proteins (C2H2-ZFPs) regulate various developmental processes and abiotic stress responses in eukaryotes. Yet, a comprehensive analysis of these transcription factors which could be used to find candidate genes related to the control the development and abiotic stress tolerance has not been performed in Pleurotus ostreatus. To fill this knowledge gap, 18 C2H2-ZFs were identified in the P. ostreatus genome. Phylogenetic analysis indicated that these proteins have dissimilar amino acid sequences. In addition, these proteins had variable protein characteristics, gene intron-exon structures, and motif compositions. The expression patterns of PoC2H2-ZFs in mycelia, primordia, and young and mature fruiting bodies were investigated using qRT-PCR. The expression of some PoC2H2-ZFs is regulated by auxin and cytokinin. Moreover, members of PoC2H2-ZFs expression levels are changed dramatically under heat and cold stress, suggesting that these genes may participate in abiotic stress responses. These findings could be used to study the role of P. ostreatus-derived C2H2-ZFs in development and stress tolerance.

Evolutionary Morphogenesis of Sexual Fruiting Bodies in Basidiomycota: Toward a New Evo-Devo Synthesis

The development of sexual fruiting bodies is one of the most complex morphogenetic processes in fungi. Mycologists have long been fascinated by the morphological and developmental diversity of fruiting bodies; however, evolutionary developmental biology of fungi still lags significantly behind that of animals or plants. Here, we summarize the current state of knowledge on fruiting bodies of mushroom-forming Basidiomycota, focusing on phylogenetic and developmental biology. Phylogenetic approaches have revealed a complex history of morphological transformations and convergence in fruiting body morphologies. Frequent transformations and convergence is characteristic of fruiting bodies in contrast to animals or plants, where main body plans are highly conserved. At the same time, insights into the genetic bases of fruiting body development have been achieved using forward and reverse genetic approaches in selected model systems. Phylogenetic and developmental studies of fruiting bodies have each yielded major advances, but they have produced largely disjunct bodies of knowledge. An integrative approach, combining phylogenetic, developmental, and functional biology, is needed to achieve a true fungal evolutionary developmental biology (evo-devo) synthesis for fungal fruiting bodies.

Unveiling the Role Displayed by Penicillium digitatum PdMut3 Transcription Factor in Pathogen-Fruit Interaction

Zn2Cys6 transcription factors are unique to fungi and are involved in different regulatory functions. In this study, we have identified the Penicillium digitatumPdMut3 gene, which encodes a putative Zn (II) 2Cys6 DNA-binding protein. Elimination of PdMut3 in Pd1 strain caused increased virulence during citrus infection. The transcription of the PdMut3 gene showed a higher expression rate during fungal growth and less transcription during fruit infection. Furthermore, the deletion of the gene in the wild-type isolate of P. digitatum did not produce any modification of the sensitivity to different fungicides, indicating that the gene is not associated with resistance to fungicides. In contrast, PdMut3 null mutants showed a reduction in growth in minimal media, which was associated with severe alterations in conidiophore development and morphological alterations of the hyphae. Mutants showed greater sensitivity to compounds that interfere with the cell wall and an invasive growth block. Thus, PdMut3 might have an indirect role in fungi virulence through metabolism and peroxisomes development.

Effect of oxygen concentrations and branched-chain amino acids on the growth and development of sub-seafloor fungus, Schizophyllum commune 20R-7-F01

Fungi have been reported to be the dominant eukaryotic group in anoxic sub-seafloor sediments, but how fungi subsist in the anoxic sub-marine sedimental environment is rarely understood. Our previous study demonstrated that the fungus, Schizophyllum commune 20R-7-F01 isolated from a ~2 km sediment below the seafloor, can grow and produce primordia in the complete absence of oxygen with enhanced production of branched-chain amino acids (BCAAs), but the primordia cannot be developed into fruit bodies without oxygen. Here, we present the individual and synergistic effects of oxygen and BCAAs on the fruit-body development of this strain. It was found that the fungus required a minimum oxygen concentration of 0.5% pO₂ to generate primordia and 1% pO₂ to convert primordia into mature fruit body. However, if BCAAs (20 mM) were added to the medium, the primordium could be developed into fruit body at a lower oxygen concentration up to 0.5% pO₂ where genes fst4 and c2h2 playing an important role in compensating oxygen deficiency. Moreover, under hypoxic conditions, the fungus showed an increase in mitochondrial number and initiation of auto-phagocytosis. These findings suggest that the fruit-body formation of S. commune may have multiple mechanisms, including energy and amino acid metabolism in response to oxygen concentrations.

Effects of exogenous ascorbic acid on the mycelia growth and primordia formation of Pleurotus ostreatus

Primordia formation is the first and most critical step in the development of fruiting bodies of edible fungi. In this study, the effects of exogenous ascorbic acid (ASA) on the Pleurotus ostreatus mycelia growth and primordia formation were researched and the results showed that the growth rate of P. ostreatus mycelia was accelerated and the time of primordia formation was advanced. The protein content and ascorbate oxidase (AAO) activity analysis showed that with the increase of ASA concentration, the protein content of mycelia first decreased and then increased, and in a certain concentration range, exogenous ASA could significantly promote the activity of AAO. Further expression analysis of the development regulating genes (Pofst3 and Pofst4) as well as blue light receptor coding genes (PoWC-1 and PoWC-2) showed the expression levels of those four genes all changed after the exogenous ASA addition, which indicated that the expression changes of PoWC-1 and PoWC-2, two key genes in the light morphogenesis, might affect the expression levels of development regulating genes Pofst3 and Pofst4, so as to lead to the formation of primordia in advance.

What Role Might Non-Mating Receptors Play in Schizophyllum commune?

The B mating-type locus of the tetrapolar basidiomycete Schizophyllum commune encodes pheromones and pheromone receptors in multiple allelic specificities. This work adds substantial new evidence into the organization of the B mating-type loci of distantly related S. commune strains showing a high level of synteny in gene order and neighboring genes. Four pheromone receptor-like genes were found in the genome of S. commune with brl1, brl2 and brl3 located at the B mating-type locus, whereas brl4 is located separately. Expression analysis of brl genes in different developmental stages indicates a function in filamentous growth and mating. Based on the extensive sequence analysis and functional characterization of brl-overexpression mutants, a function of Brl1 in mating is proposed, while Brl3, Brl4 and Brl2 (to a lower extent) have a role in vegetative growth, possible determination of growth direction. The brl3 and brl4 overexpression mutants had a dikaryon-like, irregular and feathery phenotype, and they avoided the formation of same-clone colonies on solid medium, which points towards enhanced detection of self-signals. These data are supported by localization of Brl fusion proteins in tips, at septa and in not-yet-fused clamps of a dikaryon, confirming their importance for growth and development in S. commune.

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.

H3K4me2 ChIP-Seq reveals the epigenetic landscape during mushroom formation and novel developmental regulators of Schizophyllum commune

Mushroom formation represents the most complex multicellular development in fungi. In the model mushroom Schizophyllum commune, comparative genomics and transcriptomics have previously resulted in a regulatory model of mushroom development. However, little is known about the role of epigenetic regulation. We used chromatin immunoprecipitation sequencing (ChIP-Seq) to determine the distribution of dimethylation of lysine 4 on histone H3 (H3K4me2), a mark for transcriptionally active genes, during monokaryotic and dikaryotic development. We identified a total of 6032 and 5889 sites during monokaryotic and dikaryotic development, respectively. The sites were strongly enriched near translation initiation sites of genes. Although the overall epigenetic landscape was similar between both conditions, we identified 837 sites of differential enrichment during monokaryotic or dikaryotic development, associated with 965 genes. Six transcription factor genes were enriched in H3K4me2 during dikaryotic development, indicating that these are epigenetically regulated during development. Deletion of two of these genes (fst1 and zfc7) resulted in arrested development of fruiting bodies, resulting in immature mushrooms. Together these results indicate that H3K4me2 ChIP-Seq is a powerful new tool to map the restructuring of the epigenetic landscape during mushroom development. Moreover, it can be used to identify novel developmental regulators.

The transcription factor FvHmg1 negatively regulates fruiting body development in Winter Mushroom Flammulina velutipes

Fruiting body formation in Agaricomycetes represents the most complex and unclear process in the fungi. Mating type pathways (matA and matB) and transcription factors are important regulators in the process. Here, we report a new High-mobility-group (HMG) box domain protein FvHmg1 that acts as a negative transcription regulator in fruiting body development in Winter Mushroom Flammulina velutipes. However, the expression of Fvhmg1 in dikaryon and primordial stages was significantly lower than that of monokaryon. The Fvhmg1-RNAi mutants had a better ability of fruiting than wild type strain. Overall expression of Fvhmg1 was controlled under compatible matA and matB genes where compatible matA genes could increase its expression level, while compatible matB genes had the opposite effect. It means when two monokaryons with compatible matA and matB genes were crossed, the negatively transcription factor FvHmg1 was inhibited, and normal fully fruiting body could formation and develop. The relationship between FvHmg1 and mating type pathway would advance to understand of sexual reproduction and fruiting body development in edible mushrooms.

Molecular Mechanism by Which the GATA Transcription Factor CcNsdD2 Regulates the Developmental Fate of Coprinopsis cinerea under Dark or Light Conditions

Coprinopsis cinerea has seven homologs of the Aspergillus nidulans transcription factor NsdD. Of these, CcNsdD1 and CcNsdD2 from C. cinerea show the best identities of 62 and 50% to A. nidulans NsdD, respectively. After 4 days of constant darkness cultivation, CcnsdD2, but not CcnsdD1, was upregulated on the first day of light/dark cultivation to induce fruiting bodies, and overexpression of CcnsdD2, but not CcnsdD1, produced more fruiting bodies under a light/dark rhythm. Although single knockdown of CcnsdD2 did not affect fruiting body production due to upregulation of its homolog CcnsdD1, the double-knockdown CcNsdD1/NsdD2-RNAi transformant showed defects in fruiting body formation under a light/dark rhythm. Knockdown of CcnsdD1/nsdD2 led to the differentiation of primary hyphal knots into sclerotia rather than secondary hyphal knots under a light/dark rhythm, similar to the differentiation of primary hyphal knots into sclerotia of the wild-type strain under darkness. The CcNsdD2-overexpressing transformant produced more primary hyphal knots, secondary hyphal knots, and fruiting bodies under a light/dark rhythm but only more primary hyphal knots and sclerotia under darkness. RNA-seq revealed that some genes reported previously to be involved in formation of hyphal knots and primordia, cyclopropane-fatty-acyl-phospholipid synthases cfs1-3, galectins cgl1-3, and hydrophobins hyd1-3 were downregulated in the CcNsdD1/NsdD2-RNAi transformant compared to the mock transformant. ChIP-seq and electrophoretic mobility shift assay demonstrated that CcNsdD2 bound to promoter regulatory sequences containing a GATC motif in cfs1, cfs2, cgl1, and hyd1. A molecular mechanism by which CcNsdD2 regulates the developmental fate of C. cinerea under dark or light conditions is proposed. IMPORTANCE The model mushroom Coprinopsis cinerea exhibits remarkable photomorphogenesis during fruiting body development. This study reports that the C. cinerea transcription factor CcNsdD2 promotes primary hyphal knot formation by upregulating cfs1, cfs2, cgl1, and hyd1. Although the induction of CcnsdD2 is not under direct control of light and photoreceptors, the CcNsdD2-mediated developmental fates of the primary hyphal knots depend on the following light/dark rhythm cultivation or dark cultivation after full growth of mycelia in the constant dark cultivation. This study provides new insight into the molecular mechanism by which CcNsdD2 regulates the developmental fate of C. cinerea under dark or light conditions. In addition, the result that overexpression of CcnsdD2 induced more secondary hyphal knots, primordia, and fruiting bodies under light/dark rhythm cultivation conditions has potential applied value in the edible mushroom industry.

The anaerobic survival mechanism of Schizophyllum commune 20R-7-F01, isolated from deep sediment 2 km below the seafloor

Fungi dominated the eukaryotic group in the anaerobic sedimentary environment below the ocean floor where they play an essential ecological role. However, the adaptive mechanism of fungi to these anaerobic environments is still unclear. Here, we reported the anaerobic adaptive mechanism of Schizophyllum commune 20R-7-F01, isolated from deep coal-bearing sediment down to ~2 km below the seafloor, through biochemical, metabolomic and transcriptome analyses. The fungus grows well, but the morphology changes obviously and the fruit body develops incompletely under complete hypoxia. Compared with aerobic conditions, the fungus has enhanced branched-chain amino acid biosynthesis and ethanol fermentation under anaerobic conditions, and genes related to these metabolisms have been significantly up-regulated. Additionally, the fungus shows novel strategies for synthesizing ethanol by utilizing both glycolysis and ethanol fermentation pathways. These findings suggest that the subseafloor fungi may adopt multiple mechanisms to cope with lack of oxygen.

Transcriptome data reveal conserved patterns of fruiting body development and response to heat stress in the mushroom-forming fungus Flammulina filiformis

Mushroom-forming fungi are complex multicellular organisms that form the basis of a large industry, yet, our understanding of the mechanisms of mushroom development and its responses to various stresses remains limited. The winter mushroom (Flammulina filiformis) is cultivated at a large commercial scale in East Asia and is a species with a preference for low temperatures. This study investigated fruiting body development in F. filiformis by comparing transcriptomes of 4 developmental stages, and compared the developmental genes to a 200-genome dataset to identify conserved genes involved in fruiting body development, and examined the response of heat sensitive and -resistant strains to heat stress. Our data revealed widely conserved genes involved in primordium development of F. filiformis, many of which originated before the emergence of the Agaricomycetes, indicating co-option for complex multicellularity during evolution. We also revealed several notable fruiting-specific genes, including the genes with conserved stipe-specific expression patterns and the others which related to sexual development, water absorption, basidium formation and sporulation, among others. Comparative analysis revealed that heat stress induced more genes in the heat resistant strain (M1) than in the heat sensitive one (XR). Of particular importance are the hsp70, hsp90 and fes1 genes, which may facilitate the adjustment to heat stress in the early stages of fruiting body development. These data highlighted novel genes involved in complex multicellular development in fungi and aid further studies on gene function and efforts to improve the productivity and heat tolerance in mushroom-forming fungi.

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.

Critical Factors Involved in Primordia Building in Agaricus bisporus: A Review

The button mushroom Agaricus bisporus is an economically important crop worldwide. Many aspects of its cultivation are well known, except for the precise biological triggers for its fructification. By and large, for most basidiomycete species, nutrient availability, light and a drop in temperature are critical factors for fructification. A. bisporus deviates from this pattern in the sense that it does not require light for fructification. Furthermore its fructification seems to be inhibited by a self-generated factor which needs to be removed by microorganisms in order to initiate fruiting. This review explores what is known about the morphogenesis of fruiting initiation in A. bisporus, the microflora, the self-inhibitors for fruiting initiation and transcription factors involved. This information is subsequently contrasted with an overall model of the regulatory system involved in the initiation of the formation of primordia in basidiomycetes. The comparison reveals a number of the blank spots in our understanding of the fruiting process in A. bisporus.

Comparative analysis of transcriptomes revealed the molecular mechanism of development of Tricholoma matsutake at different stages of fruiting bodies

The purpose of the study is to investigate the molecular mechanisms of development of Tricholoma matsutake fruiting body at the primordial stage (TM-1), the intermediate stage (TM-2) and the mature stage (TM-3) using RNA-Seq sequencing technology. The analysis of gene expression level revealed that the Spn2 and Eef1a1 gene were the key genes in the primordial stage of T. matsutake by regulating cytokinesis, protein synthesis, and cell growth. And the Ubc, Atp6, Cytb, and Pth2 gene were the key genes in the mature stage of T. matsutake by regulating energy metabolism and protein synthesis. Differential expression genes (DEGs) analysis results showed that Cdc28, Rad53, Dun1, Pho85 and Pho81 were the key DEGs regulating cell cycle genes of T. matsutake from primordial stage to intermediate stage. And APC, Cyr1, Cdc45, Spo11 and Rec8 genes were the key DEGs for the meiosis and sporogenesis of T. matsutake from the intermediate stage to the mature stage.

Genome-wide characterization of the Zn(II)2Cys6 zinc cluster-encoding gene family in Pleurotus ostreatus and expression analyses of this family during developmental stages and under heat stress

Pleurotus ostreatus is one of the most widely cultivated mushrooms in China. The regulatory mechanisms of fruiting body formation and the response to heat stress in P. ostreatus are main research focuses. The Zn(II)₂Cys₆ family is one of the largest families of transcriptional factors and plays important roles in multiple biological processes in fungi. In this study, we identified 66 zinc cluster proteins in P. ostreatus (PoZCPs) through a genome-wide search. The PoZCPs were classified into 15 types according to their zinc cluster domain. Physical and chemical property analyses showed a huge diversity among the PoZCPs. Phylogenetic analysis of PoZCPs classified these proteins into six groups and conserved motif combinations and similar gene structures were observed in each group. The expression profiles of these PoZCP genes during different developmental stages and under heat stress were further investigated by RNA-sequencing (RNA-seq), revealing diverse expression patterns. A total of 13 PoZCPs that may participate in development or the heat stress response were selected for validation of their expression levels through real-time quantitative PCR (RT-qPCR) analysis, and some developmental stage-specific and heat stress-responsive candidates were identified. The findings contribute to our understanding of the roles and regulatory mechanisms of ZCPs in P. ostreatus.

A putative transcription factor LFC1 negatively regulates development and yield of winter mushroom

Basidioma is the fruiting body of mushroom species. The deep understanding on the mechanism of basidioma development is valuable for mushroom breeding and cultivation. From winter mushroom (Flammulina velutipes), one of the top five industrially cultivated mushrooms, a novel putative Zn(II)2Cys6 transcription factor LFC1 with negative regulatory function in basidioma development was identified. The transcript level of lfc1 was dramatically decreased during basidioma development. Neither overexpression nor knockdown of lfc1 affected hyphal vegetative growth. However, knockdown of lfc1 could promote basidioma development and shorten cultivation time by 2 days, while overexpression of lfc1 delayed the optimal harvest time by 3 days. In the lfc1 knockdown strain, in which the lfc1 expression was reduced by 72%, mushroom yield and biological efficiency could be increased at least by 24%. Knockdown of lfc1 did not affect the shape of caps but significantly increased basidioma length and number, while its overexpression did not affect basidioma length but dramatically reduced basidioma number. In addition, rather than producing basidiomata with round caps as in wild type, the caps of basidiomata in the lfc1 overexpression mutants were significantly larger and the cap edge was wrinkled. RNA-seq analysis revealed that 455 genes had opposite transcriptional responses to lfc1 overexpression and knockdown. Some of them were previously reported as genes involved in basidioma development, including 3 hydrophobin encoding genes, 2 lectin encoding genes, FVFD16, an Eln2 ortholog encoding gene, and 3 genes encoding membrane components. As LFC1 homologs are widely present in mushroom species, lfc1 can be useful in mushroom breeding.Key Points• A novel transcription factor LFC1 negatively regulates fruiting in winter mushroom• LFC1 regulated transcription of more than 400 genes.• Reduction of LFC1 expression could shorten cultivation time and increase yield.• lfc1 could be a potentially useful reference gene for mushroom breeding.

Unusual genome expansion and transcription suppression in ectomycorrhizal Tricholoma matsutake by insertions of transposable elements

Genome sequencing of Tricholoma matsutake revealed its unusually large size as 189.0 Mbp, which is a consequence of extraordinarily high transposable element (TE) content. We identified that 702 genes were surrounded by TEs, and 83.2% of these genes were not transcribed at any developmental stage. This observation indicated that the insertion of TEs alters the transcription of the genes neighboring these TEs. Repeat-induced point mutation, such as C to T hypermutation with a bias over "CpG" dinucleotides, was also recognized in this genome, representing a typical defense mechanism against TEs during evolution. Many transcription factor genes were activated in both the primordia and fruiting body stages, which indicates that many regulatory processes are shared during the developmental stages. Small secreted protein genes (<300 aa) were dominantly transcribed in the hyphae, where symbiotic interactions occur with the hosts. Comparative analysis with 37 Agaricomycetes genomes revealed that IstB-like domains (PF01695) were conserved across taxonomically diverse mycorrhizal genomes, where the T. matsutake genome contained four copies of this domain. Three of the IstB-like genes were overexpressed in the hyphae. Similar to other ectomycorrhizal genomes, the CAZyme gene set was reduced in T. matsutake, including losses in the glycoside hydrolase genes. The T. matsutake genome sequence provides insight into the causes and consequences of genome size inflation.

Comparative Proteomic Analysis of Pleurotus ostreatus Reveals Great Metabolic Differences in the Cap and Stipe Development and the Potential Role of Ca2+ in the Primordium Differentiation

Pleurotus ostreatus is a widely cultivated edible fungus around the world. At present, studies on the developmental process of the fruiting body are limited. In our study, we compared the differentially expressed proteins (DEPs) in the stipe and cap of the fruiting body by high-throughput proteomics. GO and pathway analysis revealed the great differences in the metabolic levels, including sucrose and starch metabolism, and sphingolipid signaling and metabolism, and the differences of 16 important DEPs were validated further by qPCR analysis in expression level. In order to control the cap and stipe development, several chemical inducers were applied to the primordium of the fruiting body according to the pathway enrichment results. We found that CaCl₂ can affect the primordium differentiation through inhibiting the stipe development. EGTA (ethyleneglycol bis (β-aminoethyl ether)-N,N,N′,N′-tetraacetic acid) treatment confirmed the inhibitory role of Ca²⁺ in the stipe development. Our study not only shows great metabolic differences during the cap and stipe development but also reveals the underlying mechanism directing the primordium differentiation in the early development of the fruiting body for the first time. Most importantly, we provide a reliable application strategy for the cultivation and improvement of the Pleurotus ostreatus, which can be an example and reference for a more edible fungus.

A Single Transcription Factor (PDD1) Determines Development and Yield of Winter Mushroom (Flammulina velutipes)

Most of the edible mushrooms cannot be cultivated or have low yield under industrial conditions, partially due to the lack of knowledge on how basidioma (fruiting body) development is regulated. From winter mushroom (Flammulina velutipes), one of the most popular industrially cultivated mushrooms, a transcription factor, PDD1, with a high-mobility group (HMG)-box domain was identified based on its increased transcription during basidioma development. pdd1 knockdown by RNA interference affected vegetative growth and dramatically impaired basidioma development. A strain with an 89.9% reduction in the level of pdd1 transcription failed to produce primordia, while overexpression of pdd1 promoted basidioma development. When the transcriptional level of pdd1 was increased to 5 times the base level, the mushroom cultivation time was shortened by 9.8% and the yield was increased by at least 33%. RNA sequencing (RNA-seq) analysis revealed that pdd1 knockdown downregulated 331 genes and upregulated 463 genes. PDD1 positively regulated several genes related to fruiting, including 6 pheromone receptor-encoding genes, 3 jacalin-related lectin-encoding genes, FVFD16, and 2 FVFD16 homolog-encoding genes. PDD1 is a novel transcription factor with regulatory function in basidioma development found in industrially cultivated mushrooms. Since its orthologs are widely present in fungal species of the Basidiomycota phylum, PDD1 might have important application prospects in mushroom breeding.IMPORTANCE Mushrooms are sources of food and medicine and provide abundant nutrients and bioactive compounds. However, most of the edible mushrooms cannot be cultivated commercially due to the limited understanding of basidioma development. From winter mushroom (Flammulina velutipes; also known as Enokitake), one of the most commonly cultivated mushrooms, we identified a novel transcription factor, PDD1, positively regulating basidioma development. PDD1 increases expression during basidioma development. Artificially increasing its expression promoted basidioma formation and dramatically increased mushroom yield, while reducing its expression dramatically impaired its development. In its PDD1 overexpression mutants, mushroom number, height, yield, and biological efficiency were significantly increased. PDD1 regulates the expression of some genes that are important in or related to basidioma development. PDD1 is the first identified transcription factor with defined functions in mushroom development among commercially cultivated mushroom species, and it might be useful in mushroom breeding.

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.

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.

A hydrophobin gene, Hyd9, plays an important role in the formation of aerial hyphae and primordia in Flammulina filiformis

Flammulina filiformis is an edible fungus that is largely cultivated and widely consumed around the world. The quantity and quality of the primordia, which gives rise to the fruiting body, affects its production efficiency. Hydrophobins are involved in the formation of the fruiting body of macrofungi. However, functional verification of the hydrophobin genes is limited to date. In this study, we used gene silencing and overexpression analyses to investigate the function of one F. filiformis hydrophobin gene (Hyd9) during the development of the fruiting body. The Hyd9-silenced transformants exhibited sparse aerial hyphae, resulting in fewer primordia and fruiting bodies. In contrast, the Hyd9 overexpression strain displayed denser aerial hyphae and more primordia. The phenotypes of these transgenic lines strongly suggested that Hyd9 plays an important role in the formation of aerial hyphal knots (the primary stage of primordia) and primordia in F. filiformis. These results will be beneficial for developing more efficient methods to induce primordia formation in F. filiformis and other commercially valuable mushrooms.

Transcriptional profiling of Auricularia cornea in selenium accumulation

Auricularia cornea is a widely cultivated edible fungus with substantial nutritive value. This study aimed to enrich the multifunctional bionutrient element selenium in A. cornea to improve its quality and explore the accumulation of selenium in the fungus using high-throughput RNA-Seq technology. In general, the treatment group with a 100 µg/g supply of selenium outperformed the other treatment groups in terms of high yield, rich crude polysaccharides and a high total selenium concentration. Additional evidences demonstrated the budding and mature phases were two typical growth stages of A. cornea and were important for the accumulation of selenium. Therefore, the budding and mature phase tissues of A. cornea in the treatment group with a 100 µg/g supply of selenium were used for transcriptome analysis and compared to those of a control group that lacked additional selenium. A total of 2.56 × 105 unigenes from A. cornea transcriptome were assembled and annotated to five frequently used databases including NR, GO, KEGG, eggNOG and SwissProt. GO and KEGG pathway analysis revealed that genes involved in metabolic process and translation were up-expressed at the budding stage in response to selenium supplementation, including amino acid metabolism, lipid metabolism, ribosome. In addition, the differential gene expression patterns of A. cornea suggested that the up-expressed genes were more likely to be detected at the budding stage than at the mature stage. These results provide insights into the transcriptional response of A. cornea to selenium accumulation.