Strand-Specific RNA-Seq Analyses of Fruiting Body Development in Coprinopsis cinerea

Full-length Chd1 of Coprinopsis cinerea is expressed after the dark period required for fruiting body maturation and impacts meiotic progression

The maturation of the fruiting body primordia in the Agaricomycete Coprinopsis cinerea is triggered by light exposure, followed by a required dark period to complete maturation. During this maturation phase, meiosis occurs within basidia arranged on the surface of the gills (lamellae) on the underside of the cap. However, the molecular events required during the dark period for fruiting body maturation remain elusive. We identified a developmental mutant that fails to mature fruiting bodies under light/dark conditions. The mutant fruiting bodies resembled those arrested by the wild-type strains cultured under continuous light. The gene responsible for this mutant phenotype encodes a chromodomain helicase DNA-binding protein 1 (Chd1) homolog, Cc.Chd1. RNA-seq revealed a low transcriptional region (LTcR) within the Cc.chd1 gene. This suggests that a short version of Cc.Chd1 (predicted 1125 aa, Cc.Chd1S) is translated from the vegetative mycelium stage until before karyogamy. In contrast, the full-length Cc.Chd1 (predicted 1441 aa, Cc.Chd1L) is translated during or after the dark period when karyogamy occurs in the basidia. Western blot analysis confirmed these types of Cc.Chd1 at the expected stages. Microscopic observations further revealed that meiotic chromosomes in basidia become arrested at prophase I in the Cc.chd1-1 mutant and wild-type strains cultured under continuous light. These findings suggest that Cc.Chd1L is required for progression from meiotic prophase I to metaphase I. Additionally, the Cc.chd1 mutant exhibits defects in light-induced secondary knot formation, suggesting a role for Cc.Chd1S in this process.

Pleurotus ostreatus mek1 is essential for meiosis and basidiospore production

Breeding strains without basidiospores is important for the mushroom industry. However, target genes for sporeless breeding remain limited. To identify a new gene essential for basidiospore production in Pleurotus ostreatus, the RNA sequencing data obtained in our previous study were re-analyzed. Among the 36 P. ostreatus genes that are exclusively expressed in the gills where basidiospores are formed and produced, candidate genes were narrowed down using transcriptome data during the fruiting stage of Coprinopsis cinerea, in which the meiotic steps progress synchronously. Three C. cinerea genes homologous to the gill-specific P. ostreatus genes were upregulated during meiosis. One of these three genes encodes a protein homologous to Saccharomyces cerevisiae Mek1p, a protein kinase important in the meiotic recombination checkpoint. Plasmids containing expression cassettes for hygromycin B-resistance screening, Cas9, and single-guide RNA targeting mek1 were introduced into the protoplasts of the dikaryotic P. ostreatus strain PC9×#64, which showed that the obtained dikaryotic gene disruptant produced no basidiospores. Microscopic analysis suggests that meiosis is suspended during telophase I. These results suggested that P. ostreatus Mek1 is essential for meiosis II progression and basidiospore production. In addition, disturbed orientation and loss of negative gravitropism during fruiting were observed.

Morphogenesis, starvation, and light responses in a mushroom-forming fungus revealed by long-read sequencing and extensive expression profiling

Mushroom-forming fungi (Agaricomycetes) are emerging as pivotal players in several fields of science and industry. Genomic data for Agaricomycetes are accumulating rapidly; however, this is not paralleled by improvements of gene annotations, which leave gene function notoriously poorly understood. We set out to improve our functional understanding of the model mushroom Coprinopsis cinerea by integrating a new, chromosome-level assembly, high-quality gene predictions, and functional information derived from broad gene-expression profiling data. The new annotation includes 5' and 3' untranslated regions (UTRs), polyadenylation sites (PASs), upstream open reading frames (uORFs), splicing isoforms, and microexons, as well as core gene sets corresponding to carbon starvation, light response, and hyphal differentiation. As a result, the genome of C. cinerea has now become the most comprehensively annotated genome among mushroom-forming fungi, which will contribute to multiple rapidly expanding fields, including research on their life history, light and stress responses, as well as multicellular development.

Whole-Genome Sequencing and Fine Map Analysis of Pholiota nameko

Pholiota nameko (T. Ito) S. Ito and S. Imai is an emerging wild mushroom species belonging to the genus Pholiota. Its unique brown-yellow appearance and significant biological activity have garnered increasing attention in recent years. However, there is a relative lack of research on the biological characteristics and genetics of P. nameko, which greatly limits the potential for an in-depth exploration of this mushroom in the research fields of molecular breeding and evolutionary biology. This study aimed to address that gap by employing Illumina and Nanopore sequencing technologies to perform whole-genome sequencing, de novo assembly, and annotation analysis of the P. nameko ZZ1 strain. Utilizing bioinformatics methods, we conducted a comprehensive analysis of the genomic characteristics of this strain and successfully identified candidate genes associated with its mating type, carbohydrate-active enzymes, virulence factors, pan-genome, and drug resistance functions. The genome of P. nameko ZZ1 is 24.58 Mb in size and comprises 33 contigs, with a contig N50 of 2.11 Mb. A hylogenetic analysis further elucidated the genetic relationship between P. nameko and other Pholiota, revealing a high degree of collinearity between P. nameko and ZZ1. In our enzyme analysis, we identified 246 enzymes in the ZZ1 genome, including 68 key carbohydrate-active enzymes (CAZymes), and predicted the presence of 11 laccases, highlighting the strain's strong potential for cellulose degradation. We conducted a pan-genomic analysis of five closely related strains of Pholiota, yielding extensive genomic information. Among these, there were 2608 core genes, accounting for 21.35% of the total genes, and 135 dispensable genes, highlighting significant genetic diversity among Pholiota and further confirming the value of pan-genomic analysis in uncovering species diversity. Notably, while we successfully identified the A-mating-type locus, composed of the homeodomain protein genes HD1 and HD2 in ZZ1, we were unable to obtain the B-mating-type locus due to technical limitations, preventing us from acquiring the pheromone receptor of the B-mating-type. We plan to supplement these data in future studies and explore the potential impact of the B-mating-type locus on the current findings. In summary, the genome data of ZZ1 presented in this study are not only valuable resources for understanding the genetic basis of this species, but also serve as a crucial foundation for subsequent genome-assisted breeding, research into cultivation technology, and the exploration of its nutritional and potential medicinal value.

Identification of Homeobox Transcription Factors in a Dimorphic Fungus Talaromyces marneffei and Protein-Protein Interaction Prediction of RfeB

Talaromyces marneffei is a thermally dimorphic fungus that can cause life-threatening systemic mycoses, particularly in immunocompromised individuals. Fungal homeobox transcription factors control various developmental processes, including the regulation of sexual reproduction, morphology, metabolism, and virulence. However, the function of homeobox proteins in T. marneffei has not been fully explored. Here, we searched the T. marneffei genome for the total homeobox transcription factors and predicted their biological relevance by performing gene expression analysis in different cell types, including conidia, mycelia, yeasts, and during phase transition. RfeB is selected for further computational analysis since (i) its transcripts were differentially expressed in different phases of T. marneffei, and (ii) this protein contains the highly conserved protein-protein interaction region (IR), which could be important for pathobiology and have therapeutic application. To assess the structure-function of the IR region, in silico alanine substitutions were performed at three-conserved IR residues (Asp276, Glu279, and Gln282) of RfeB, generating a triple RfeB mutated protein. Using 3D modeling and molecular dynamics simulations, we compared the protein complex formation of wild-type and mutated RfeB proteins with the putative partner candidate TmSwi5. Our results demonstrated that the mutated RfeB protein exhibited increased free binding energy, elevated protein compactness, and a reduced number of atomic contacts, suggesting disrupted protein stability and interaction. Notably, our model revealed that the IR residues primarily stabilized the RfeB binding sites located in the central region (CR). This computational approach for protein mutagenesis could provide a foundation for future experimental studies on the functional characterization of RfeB and other homeodomain-containing proteins in T. marneffei.

Characterization of linoleate dioxygenases in basidiomycetes and the functional role of CcLdo1 in regulating fruiting body development in Coprinopsis cinerea

Coprinopsis cinerea, a model fungus, is utilized for investigating the developmental mechanisms of basidiomycetes. The development of basidiomycetes is a highly organized process that requires coordination among genetic, environmental, and physiological factors. Oxylipins, a class of widely distributed signaling molecules, play crucial roles in fungal biology. Among oxylipins, the sexual pheromone-inducing factors (psi factors) have been identified as key regulators of the balance between asexual and sexual spore development in Ascomycetes. Linoleate dioxygenases are enzymes involved in the biosynthesis of psi factors, yet their specific physiological functions in basidiomycete development remain unclear. In this study, linoleate dioxygenases in basidiomycetes were identified and characterized. Phylogenetic analysis revealed that linoleate dioxygenases from Basidiomycota formed a distinct clade, with linoleate dioxygenases from Agaricomycetes segregating into three groups and those from Ustilaginomycetes forming a separate group. Both basidiomycete and ascomycete linoleate dioxygenases shared two characteristic domains: the N-terminal of linoleate dioxygenase domain and the C-terminal of cytochrome P450 domain. While the linoleate dioxygenase domains exhibited similarity between basidiomycetes and ascomycetes, the cytochrome P450 domains displayed high diversity in key sites. Furthermore, the gene encoding the linoleate dioxygenase Ccldo1 in C. cinerea was knocked out, resulting in a significant increase in fruiting body formation without affecting asexual conidia production. This observation suggests that secondary metabolites synthesized by CcLdo1 negatively regulate the sexual reproduction process in C. cinerea while not influencing the asexual reproductive process. This study represents the first identification of a gene involved in secondary metabolite synthesis that regulates basidiocarp development in a basidiomycete.

Identification of two genes essential for basidiospore formation during the postmeiotic stages in Pleurotus ostreatus

A sporeless strain is an important breeding target in the mushroom industry. However, basidiospore production in the oyster mushroom Pleurotus ostreatus has been shown to be impaired by single-gene mutations in only two meiosis-related genes, mer3 and msh4. This study proposed a strategy for identifying the genes essential for basidiospore formation after meiotic division to determine new targets for molecular breeding. RNA-seq analysis was performed to identify P. ostreatus genes that are specifically expressed in the gill tissue of fruiting bodies, where basidiospore formation occurs. Transcriptome data during fruiting development of Coprinopsis cinerea, in which the meiotic steps progress synchronously, were then used to identify genes that are active in the postmeiotic stages. Based on these comparative analyses, five P. ostreatus genes were identified. Plasmids containing expression cassettes for hygromycin B-resistance screening, Cas9, and single-guide RNA targeting each gene were introduced into the protoplasts of dikaryotic strain, PC9×#64, to generate dikaryotic gene disruptants. Among the obtained transformants, three dikaryotic pcl1 disruptants and two cro6c disruptants did not produce basidiospores. Microscopic analyses indicated that spore formation was arrested at particular stages in these gene disruptants. These results indicate that these two genes are essential for mature spore formation in this fungus.

The exp2 gene, which encodes a protein with two zinc finger domains, regulates cap expansion and autolysis in Coprinopsis cinerea

Cap expansion in agaricoid mushroom species is an important event for sexual reproduction because meiosis occurs in basidia under the cap, and basidiospores can be released by opening the cap. However, molecular mechanisms underlying cap expansion in basidiomycetes remain poorly understood. We aimed to elucidate the molecular mechanisms of cap expansion in basidiomycetes by analyzing the unique cap-expansionless UV mutant #13 (exp2-1) in Coprinopsis cinerea. Linkage analysis and consequent genome sequence analysis revealed that the gene responsible for the mutant phenotypes encodes a putative transcription factor with two C2H2 zinc finger motifs. The mutant that was genome-edited to lack exp2 exhibited an expansionless phenotype. Some of the genes encoding cell wall degradation-related enzymes showed decreased expression during cap expansion and autolysis in the exp2 UV and genome-edited mutant. The exp2 gene is widely conserved in Agaricomycetes, suggesting that Exp2 homologs regulate fruiting body maturation in Agaricomycetes, especially cap expansion in Agaricoid-type mushroom-forming fungi. Therefore, exp2 homologs could be a target for mushroom breeding to maintain shape after harvest for some cultivating mushrooms, presenting a promising avenue for further research in breeding techniques.

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.

Transcriptome Analysis Reveals Mycelial and Fruiting Responses to Lithium Chloride in Coprinopsis cinerea

Lithium chloride (LiCl) has been used in signalling and molecular studies of animals, plants, and yeast. However, information on its roles in basidiomycetous fungi is still limited. In this study, we used RNA-Seq to study the effects of LiCl on Coprinopsis cinerea. LiCl enhanced mycelial growth and inhibited fruiting body formation in C. cinerea. RNA-Seq of the LiCl-treated C. cinerea resulted in a total of 14,128 genes. There were 1199 differentially expressed genes (DEGs) between the LiCl-treated samples and control samples in the mycelium stage (the first time point), and 1391 DEGs were detected when the control samples were forming hyphal knots while the treated samples were still in the mycelium (the second time point). Pathway enrichment analysis of the DEGs revealed a significant association between enhanced mycelium growth in the LiCl-treated C. cinerea and metabolic pathways. In addition, the DEGs involved in cellular process pathways, including "cell cycle-yeast" and "meiosis-yeast", were identified in suppressed C. cinerea fruiting body formation by LiCl under favourable environmental conditions. As LiCl can predominantly inhibit the activity of glycogen synthase kinase3 (GSK3), our findings suggest that LiCl affects the expression of genes involved in fruiting body initiation and cellular processes by inhibiting GSK3 activity which is essential for fruiting body formation.

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.

The Genome-Wide Characterization of Alternative Splicing and RNA Editing in the Development of Coprinopsis cinerea

Coprinopsis cinerea is one of the model species used in fungal developmental studies. This mushroom-forming Basidiomycetes fungus has several developmental destinies in response to changing environments, with dynamic developmental regulations of the organism. Although the gene expression in C. cinerea development has already been profiled broadly, previous studies have only focused on a specific stage or process of fungal development. A comprehensive perspective across different developmental paths is lacking, and a global view on the dynamic transcriptional regulations in the life cycle and the developmental paths is far from complete. In addition, knowledge on co- and post-transcriptional modifications in this fungus remains rare. In this study, we investigated the transcriptional changes and modifications in C. cinerea during the processes of spore germination, vegetative growth, oidiation, sclerotia formation, and fruiting body formation by inducing different developmental paths of the organism and profiling the transcriptomes using the high-throughput sequencing method. Transition in the identity and abundance of expressed genes drive the physiological and morphological alterations of the organism, including metabolism and multicellularity construction. Moreover, stage- and tissue-specific alternative splicing and RNA editing took place and functioned in C. cinerea. These modifications were negatively correlated to the conservation features of genes and could provide extra plasticity to the transcriptome during fungal development. We suggest that C. cinerea applies different molecular strategies in its developmental regulation, including shifts in expressed gene sets, diversifications of genetic information, and reversible diversifications of RNA molecules. Such features would increase the fungal adaptability in the rapidly changing environment, especially in the transition of developmental programs and the maintenance and balance of genetic and transcriptomic divergence. The multi-layer regulatory network of gene expression serves as the molecular basis of the functioning of developmental regulation.

Revisiting the AA14 family of lytic polysaccharide monooxygenases and their catalytic activity

Lytic polysaccharide monooxygenases (LPMOs) belonging to the AA14 family are believed to contribute to the enzymatic degradation of lignocellulosic biomass by specifically acting on xylan in recalcitrant cellulose-xylan complexes. Functional characterization of an AA14 LPMO from Trichoderma reesei, TrAA14A, and a re-evaluation of the properties of the previously described AA14 from Pycnoporus coccineus, PcoAA14A, showed that these proteins have oxidase and peroxidase activities that are common for LPMOs. However, we were not able to detect activity on cellulose-associated xylan or any other tested polysaccharide substrate, meaning that the substrate of these enzymes remains unknown. Next to raising questions regarding the true nature of AA14 LPMOs, the present data illustrate possible pitfalls in the functional characterization of these intriguing enzymes.

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.

Live-cell imaging of septins and cell polarity proteins in the growing dikaryotic vegetative hypha of the model mushroom Coprinopsis cinerea

The developmental biology underlying the morphogenesis of mushrooms remains poorly understood despite the essential role of fungi in the terrestrial environment and global carbon cycle. The mushroom Coprinopsis cinerea is a leading model system for the molecular and cellular basis of fungal morphogenesis. The dikaryotic vegetative hyphae of this fungus grow by tip growth with clamp cell formation, conjugate nuclear division, septation, subapical peg formation, and fusion of the clamp cell to the peg. Studying these processes provides many opportunities to gain insights into fungal cell morphogenesis. Here, we report the dynamics of five septins, as well as the regulators CcCla4, CcSpa2, and F-actin, visualized by tagging with fluorescent proteins, EGFP, PA-GFP or mCherry, in the growing dikaryotic vegetative hyphae. We also observed the nuclei using tagged Sumo proteins and histone H1. The five septins colocalized at the hyphal tip in the shape of a dome with a hole (DwH). CcSpa2-EGFP signals were observed in the hole, while CcCla4 signals were observed as the fluctuating dome at the hyphal tip. Before septation, CcCla4-EGFP was also occasionally recruited transiently around the future septum site. Fluorescent protein-tagged septins and F-actin together formed a contractile ring at the septum site. These distinct specialized growth machineries at different sites of dikaryotic vegetative hyphae provide a foundation to explore the differentiation program of various types of cells required for fruiting body formation.

Transcriptome Profiling Reveals Candidate Genes Related to Stipe Gradient Elongation of Flammulina filiformis

Stipe gradient elongation is an important and remarkable feature in the development of most mushroom fruiting bodies. However, its molecular mechanism has rarely been described. Here, the decreasing trend of stipe elongation and increasing trend of cell length in a gradient from the top to the base of the stipe were determined in a model basidiomycete mushroom: Flammulina filiformis. According to RNA-seq results, 1409 differentially expressed genes (DEGs) were identified among elongation region (ER), transition region (TR), and stable region (SR) samples, including 26 transcription factors (TFs). Based on Short Time-series Expression Miner (STEM) clustering of DEGs, clusters 1 and 3, with obvious expression trends that were consistent with or in contrast to the elongation rate, were screened. The cluster 1 DEGs were mainly involved in the GO cellular component category and KEGG genetic information processing class; however, the cluster 3 DEGs were mainly involved in metabolic processes. Furthermore, qRT-PCR confirmed that key genes of the long-chain fatty acid synthesis pathway were involved in stipe gradient elongation and regulated by NADPH oxidase-derived ROS signaling molecules. These findings provide an essential basis for understanding the molecular mechanism of stipe gradient elongation.

Structural and Functional Analysis of Peptides Derived from KEX2-Processed Repeat Proteins in Agaricomycetes Using Reverse Genetics and Peptidomics

Bioactivities of fungal peptides are of interest for basic research and therapeutic drug development. Some of these peptides are derived from "KEX2-processed repeat proteins" (KEPs), a recently defined class of precursor proteins that contain multiple peptide cores flanked by KEX2 protease cleavage sites. Genome mining has revealed that KEPs are widespread in the fungal kingdom. Their functions are largely unknown. Here, we present the first in-depth structural and functional analysis of KEPs in a basidiomycete. We bioinformatically identified KEP-encoding genes in the genome of the model agaricomycete Coprinopsis cinerea and established a detection protocol for the derived peptides by overexpressing the C. cinerea KEPs in the yeast Pichia pastoris. Using this protocol, which includes peptide extraction and mass spectrometry with data analysis using the search engine Mascot, we confirmed the presence of several KEP-derived peptides in C. cinerea, as well as in the edible mushrooms Lentinula edodes, Pleurotus ostreatus, and Pleurotus eryngii. While CRISPR-mediated knockout of C. cinerea kep genes did not result in any detectable phenotype, knockout of kex genes caused defects in mycelial growth and fruiting body formation. These results suggest that KEP-derived peptides may play a role in the interaction of C. cinerea with the biotic environment and that the KEP-processing KEX proteases target a variety of substrates in agaricomycetes, including some important for mycelial growth and differentiation. IMPORTANCE Two recent bioinformatics studies have demonstrated that KEX2-processed repeat proteins are widespread in the fungal kingdom. However, despite the prevalence of KEPs in fungal genomes, only few KEP-derived peptides have been detected and studied so far. Here, we present a protocol for the extraction and structural characterization of KEP-derived peptides from fungal culture supernatants and tissues. The protocol was successfully used to detect several linear and minimally modified KEP-derived peptides in the agaricomycetes C. cinerea, L. edodes, P. ostreatus, and P. eryngii. Our study establishes a new protocol for the targeted search of KEP-derived peptides in fungi, which will hopefully lead to the discovery of more of these interesting fungal peptides and allow a further characterization of KEPs.

Preassembled Cas9 Ribonucleoprotein-Mediated Gene Deletion Identifies the Carbon Catabolite Repressor and Its Target Genes in Coprinopsis cinerea

Cre1 is an important transcription factor that regulates carbon catabolite repression (CCR) and is widely conserved across fungi. The cre1 gene has been extensively studied in several Ascomycota species, whereas its role in gene expression regulation in the Basidiomycota species remains poorly understood. Here, we identified and investigated the role of cre1 in Coprinopsis cinerea, a basidiomycete model mushroom that can efficiently degrade lignocellulosic plant wastes. We used a rapid and efficient gene deletion approach based on PCR-amplified split-marker DNA cassettes together with in vitro assembled Cas9-guide RNA ribonucleoproteins (Cas9 RNPs) to generate C. cinerea cre1 gene deletion strains. Gene expression profiling of two independent C. cinerea cre1 mutants showed significant deregulation of carbohydrate metabolism, plant cell wall degrading enzymes (PCWDEs), plasma membrane transporter-related and several transcription factor-encoding genes, among others. Our results support the notion that, like reports in the ascomycetes, Cre1 of C. cinerea orchestrates CCR through a combined regulation of diverse genes, including PCWDEs, transcription factors that positively regulate PCWDEs, and membrane transporters which could import simple sugars that can induce the expression of PWCDEs. Somewhat paradoxically, though in accordance with other Agaricomycetes, genes related to lignin degradation were mostly downregulated in cre1 mutants, indicating they fall under different regulation than other PCWDEs. The gene deletion approach and the data presented here will expand our knowledge of CCR in the Basidiomycota and provide functional hypotheses on genes related to plant biomass degradation. IMPORTANCE Mushroom-forming fungi include some of the most efficient lignocellulosic plant biomass degraders. They degrade dead plant materials by a battery of lignin-, cellulose-, hemicellulose-, and pectin-degrading enzymes, the encoding genes of which are under tight transcriptional control. One of the highest-level regulations of these metabolic enzymes is known as carbon catabolite repression, which is orchestrated by the transcription factor Cre1, and ensures that costly lignocellulose-degrading enzyme genes are expressed only when simple carbon sources (e.g., glucose) are not available. Here, we identified the Cre1 ortholog in a litter decomposer Agaricomycete, Coprinopsis cinerea, knocked it out, and characterized transcriptional changes in the mutants. We identified several dozen lignocellulolytic enzyme genes as well as membrane transporters and other transcription factors as putative target genes of C. cinerea cre1. These results extend knowledge on carbon catabolite repression to litter decomposer Basidiomycota.

Identification of Key Regulatory Pathways of Basidiocarp Formation in Pleurotus spp. Using Modeling, Simulation and System Biology Studies

Pleurotus (Oyster mushroom) is an important cultivated edible mushroom across the world. It has several therapeutic effects as it contains various useful bio-molecules. The cultivation and crop management of these basidiomycete fungi depends on many extrinsic and intrinsic factors such as substrate composition, growing environment, enzymatic properties, and the genetic makeup, etc. Moreover, for efficient crop production, a comprehensive understanding of the fundamental properties viz. intrinsic–extrinsic factors and genotype-environment interaction analysis is required. The present study explores the basidiocarp formation biology in Pleurotus mushroom using an in silico response to the environmental factors and involvement of the major regulatory genes. The predictive model developed in this study indicates involvement of the key regulatory pathways in the pinhead to fruit body development process. Notably, the major regulatory pathways involved in the conversion of mycelium aggregation to pinhead formation and White Collar protein (PoWC1) binding flavin-chromophore (FAD) to activate respiratory enzymes. Overall, cell differentiation and higher expression of respiratory enzymes are the two important steps for basidiocarp formation. PoWC1 and pofst genes were participate in the structural changes process. Besides this, the PoWC1 gene is also involved in the respiratory requirement, while the OLYA6 gene is the triggering point of fruiting. The findings of the present study could be utilized to understand the detailed mechanism associated with the basidiocarp formation and to cultivate mushrooms at a sustainable level.

De novo transcriptome assembly and comprehensive assessment provide insight into fruiting body formation of Sparassis latifolia

The genes associated with fruiting body formation of Sparasis latifolia are valuable for improving mushroom breeding. To investigate this process, 4.8 × 108 RNA-Seq reads were acquired from three stages: hyphal knot (SM), primordium (SP), and primordium differentiation (SPD). The de novo assembly generated a total of 48,549 unigenes, of which 71.53% (34,728) unigenes could be annotated by at least one of the KEGG (Kyoto Encyclopedia of Genes and Genomes), GO (Gene Ontology), and KOG (Eukaryotic Orthologous Group) databases. KEGG and KOG analyses respectively mapped 32,765 unigenes to 202 pathways and 19,408 unigenes to 25 categories. KEGG pathway enrichment analysis of DEGs (differentially expressed genes) indicated primordium initiation was significantly related to 66 pathways, such as "Ribosome", "metabolism of xenobiotics by cytochrome P450", and "glutathione metabolism" (among others). The MAPK and mTOR signal transduction pathways underwent significant adjustments during the SM to SP transition. Further, our research revealed the PI3K-Akt signaling pathway related to cell proliferation could play crucial functions during the development of SP and SPD. These findings provide crucial candidate genes and pathways related to primordium differentiation and development in S. latifolia, and advances our knowledge about mushroom morphogenesis.

Transcriptome and Differentially Expressed Gene Profiles in Mycelium, Primordium and Fruiting Body Development in Stropharia rugosoannulata

Stropharia rugosoannulata uses straw as a growth substrate during artificial cultivation and has been widely promoted in China. However, its fruiting body formation and development processes have not been elucidated. In this study, the developmental transcriptomes were analyzed at three stages: the mycelium (G-S), primordium (P-S) and fruiting body (M-F) stages. A total of 9690 differentially expressed genes (DEGs) were identified in the different developmental stages. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses showed that these DEGs were involved mainly in hydrolase activity, structural molecule activity and oxidoreductase activity as well as xenobiotic biodegradation and metabolism and energy metabolism pathways. We further found that the higher expression of most carbohydrate enzyme (i.e., GH, CE, CBM, AA and PL) genes in the hyphal (i.e., G-S) stage was related mainly to substrate degradation, while the upregulation of glycosyltransferase (GT) gene expression in the P-S and M-F stages may be related to cell wall synthesis. In addition, we found that CO₂-sensing-related genes (i.e., CA-2, CA-3, PKA-1 and PKA-2) were upregulated in the P-S and M-F stages, heat shock protein genes (HSP60 and HSP90) were significantly downregulated in the P-S stage and upregulated in the M-F stage and the transcription factors (i.e., steA, MYB, nosA, HAP1, and GATA-4/5/6) involved in growth and development were significantly upregulated in the P-S stage. These results suggest that environmental factors (i.e., CO₂ and temperature) and transcription factors may play a key role in primordium formation. In short, this study provides new insights into the study of stimulating primordia formation affecting the development of fruiting bodies of S. rugosoannulata.

Structure-function relationship of a novel fucoside-binding fruiting body lectin from Coprinopsis cinerea exhibiting nematotoxic activity

Lectins are non-immunoglobulin-type proteins that bind to specific carbohydrate epitopes and play important roles in intra- and inter-organismic interactions. Here, we describe a novel fucose-specific lectin, termed CML1, which we identified from fruiting body extracts of Coprinopsis cinerea. For further characterization, the coding sequence for CML1 was cloned and heterologously expressed in Escherichia coli. Feeding of CML1-producing bacteria inhibited larval development of the bacterivorous nematode Caenorhabditis tropicalis, but not of C. elegans. The crystal structure of the recombinant protein in its apo-form and in complex with H type I or Lewis A blood group antigens was determined by X-ray crystallography. The protein folds as a sandwich of 2 antiparallel β-sheets and forms hexamers resulting from a trimer of dimers. The hexameric arrangement was confirmed by small-angle X-ray scattering (SAXS). One carbohydrate-binding site per protomer was found at the dimer interface with both protomers contributing to ligand binding, resulting in a hexavalent lectin. In terms of lectin activity of recombinant CML1, substitution of the carbohydrate-interacting residues His54, Asn55, Trp94, and Arg114 by Ala abolished carbohydrate-binding and nematotoxicity. Although no similarities to any characterized lectin were found, sequence alignments identified many non-characterized agaricomycete proteins. These results suggest that CML1 is the founding member of a novel family of fucoside-binding lectins involved in the defense of agaricomycete fruiting bodies against predation by fungivorous nematodes.

Cocaprins, β-Trefoil Fold Inhibitors of Cysteine and Aspartic Proteases from Coprinopsis cinerea

We introduce a new family of fungal protease inhibitors with β-trefoil fold from the mushroom Coprinopsis cinerea, named cocaprins, which inhibit both cysteine and aspartic proteases. Two cocaprin-encoding genes are differentially expressed in fungal tissues. One is highly transcribed in vegetative mycelium and the other in the stipes of mature fruiting bodies. Cocaprins are small proteins (15 kDa) with acidic isoelectric points that form dimers. The three-dimensional structure of cocaprin 1 showed similarity to fungal β-trefoil lectins. Cocaprins inhibit plant C1 family cysteine proteases with Ki in the micromolar range, but do not inhibit the C13 family protease legumain, which distinguishes them from mycocypins. Cocaprins also inhibit the aspartic protease pepsin with Ki in the low micromolar range. Mutagenesis revealed that the β2-β3 loop is involved in the inhibition of cysteine proteases and that the inhibitory reactive sites for aspartic and cysteine proteases are located at different positions on the protein. Their biological function is thought to be the regulation of endogenous proteolytic activities or in defense against fungal antagonists. Cocaprins are the first characterized aspartic protease inhibitors with β-trefoil fold from fungi, and demonstrate the incredible plasticity of loop functionalization in fungal proteins with β-trefoil fold.

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.

Pan-genome Analyses of 3 Strains of Inonotus obliquus and Prediction of Polysaccharide and Terpenoid Genes

Inonotus obliquus is a rare, edible and medicinal fungus that is widely used as a remedy for various diseases. Its main bioactive substances are polysaccharides and terpenoids. In this study, we characterized and investigated the pan-genome of three strains of I. obliquus. The genome sizes of JL01, HE, and NBRC8681 were 32.04, 29.04, and 31.78 Mb, respectively. There were 6 543 core gene families and 6 197 accessory gene families among the three strains, with 14 polysaccharide-related core gene families and seven accessory gene families. For terpenoids, there were 13 core gene families and 17 accessory gene families. Pan-genome sequencing of I. obliquus has improved our understanding of biological characteristics related to the biosynthesis of polysaccharides and terpenoids at the molecular level, which in turn will enable us to increase the production of polysaccharides and terpenoids by this mushroom.

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.

Microbes as a production host to produce natural activecompounds from mushrooms: biosynthetic pathway elucidationand metabolic engineering

Mushrooms are abundant in bioactive natural compounds. Due to strict growth conditions and long fermentation-time, microbe as a production host is an alternative and sustainable approach for the production of natural compounds. This review focuses on the biosynthetic pathways of mushroom originated natural compounds and microbes as the production host for the production of the above natural compounds.

Comparative transcriptome analysis of abnormal cap and healthy fruiting bodies of the edible mushroom Lentinula edodes

Lentinula edodes, a commercially important mushroom, is cultivated worldwide. Artificially cultivated L. edodes often present with abnormal symptoms in the fruiting body, which affect their commercial value and reduce production efficiency. In this study, we carried out a comparative transcriptome analysis of normal fruiting body pileus (LeNP), normal margin in abnormal fruiting body pileus (LeAPNM), and abnormal margin in abnormal fruiting body pileus (LeAPAM). Metabolic pathways such as those involved in transmembrane transport, ribosome production, tryptophan metabolism, arginine and proline metabolism, and the metabolism of other amino acids were significantly enriched in LeAPAM. F-box, short-chain dehydrogenases/reductases, the major facilitator superfamily, and the FMN_red superfamily are related to malformation in L. edodes. Genes encoding heat shock proteins, G protein, and β-1,3-glucanase in the GH5 family showed different expression patterns, suggesting that these genes are involved in the development of L. edodes fruiting bodies. In particular, CAZymes, which are involved in the development of cell walls in L. edodes, were highly expressed in LeAPAM. According to TEM observation, the cell wall of LeAPAM samples showed significant thickening compared to the other samples. These results suggested that cell wall anabolism in LeAPAM samples was more active than that in normal fruiting bodies, enhancing the environmental adaptability of the fungus. This study provides preliminary data for future research aimed at solving the phenomenon of abnormal fruiting bodies of L. edodes.

Distribution of methionine sulfoxide reductases in fungi and conservation of the free-methionine-R-sulfoxide reductase in multicellular eukaryotes

Methionine, either as a free amino acid or included in proteins, can be oxidized into methionine sulfoxide (MetO), which exists as R and S diastereomers. Almost all characterized organisms possess thiol-oxidoreductases named methionine sulfoxide reductase (Msr) enzymes to reduce MetO back to Met. MsrA and MsrB reduce the S and R diastereomers of MetO, respectively, with strict stereospecificity and are found in almost all organisms. Another type of thiol-oxidoreductase, the free-methionine-R-sulfoxide reductase (fRMsr), identified so far in prokaryotes and a few unicellular eukaryotes, reduces the R MetO diastereomer of the free amino acid. Moreover, some bacteria possess molybdenum-containing enzymes that reduce MetO, either in the free or protein-bound forms. All these Msrs play important roles in the protection of organisms against oxidative stress. Fungi are heterotrophic eukaryotes that colonize all niches on Earth and play fundamental functions, in organic matter recycling, as symbionts, or as pathogens of numerous organisms. However, our knowledge on fungal Msrs is still limited. Here, we performed a survey of msr genes in almost 700 genomes across the fungal kingdom. We show that most fungi possess one gene coding for each type of methionine sulfoxide reductase: MsrA, MsrB, and fRMsr. However, several fungi living in anaerobic environments or as obligate intracellular parasites were devoid of msr genes. Sequence inspection and phylogenetic analyses allowed us to identify non-canonical sequences with potentially novel enzymatic properties. Finaly, we identified several ocurences of msr horizontal gene transfer from bacteria to fungi.

Cloning and analysis of the photoreceptor genes Oxwc-1 and Oxwc-2 from ethnopharmacological fungus Ophiocordyceps xuefengensis

Ophiocordyceps xuefengensis is an ethnopharmacological fungus with broad pharmacological properties. Light is a critical environmental factor for the stromata formation and development of many fungi. In this study, photomorphogenesis and blue light receptor genes were studied using a strain of O. xuefengensis. Light represses vegetative growth, but conidia linked to stromata can be observed under both light and dark conditions. Light and dark conditions had little effect on the accumulation of polysaccharides and adenosine. The genes Oxwc-1 and Oxwc-2 encoding photoreceptors of O. xuefengensis were cloned and predicted to possess polypeptides of 937 and 525 amino acids, respectively. A phylogenetic analysis based on fungal WC-1/2 supported OxWC-1 and OxWC-2 were photoreceptor. The expression of both the Oxwc-1 and Oxwc-2 genes reached a maximum after receiving light stimulation for 15 min, which might relate to the inhibition of stromata growth.

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.

Selection and Evaluation of Potential Reference Genes for Quantitative Real-Time PCR in Agaricus blazei Based on Transcriptome Sequencing Data

Quantitative real-time PCR (qRT-PCR) is widely used to detect gene expression due to its high sensitivity, high throughput, and convenience. The accurate choice of reference genes is required for normalization of gene expression in qRT-PCR analysis. In order to identify the optimal candidates for gene expression analysis using qRT-PCR in Agaricus blazei, we studied the potential reference genes in this economically important edible fungus. In this study, transcriptome datasets were used as source for identification of candidate reference genes. And 27 potential reference genes including 21 newly stable genes, three classical housekeeping genes, and homologous genes of three ideal reference genes in Volvariella volvacea, were screened based on transcriptome datasets of A. blazei and previous studies. The expression stability of these genes was investigated by qRT-PCR analysis and further evaluated by four software packages, geNorm, NormFinder, BestKeeper, and RefFinder. Among these candidates, α-TUB (Tubulin alpha) and Cox5a (COX5A subunit VA of cytochrome c oxidase) were revealed as the most stable in fruit body, and suitable for 5 different developmental stages. α-TUB and ATP3 (ATP3 gamma subunit of the F1 sector of mitochondrial F1F0 ATP synthase) showed the most stable expression in stipe tissues and, Uqcrc (core subunit of the ubiquinol-cytochrome c reductase complex) and PUP3 (20S proteasome subunit beta 3) performed well in pileus tissues during the process of A. blazei development, while GAPDH (glyceraldehyde-3-phosphate dehydrogenase) was among the least stable genes in all sample sets. Finally, the Ableln3 (homology of eln3 gene of Coprinus cinereus) was adopted to validate the reliability of these stable and unstable reference genes, indicating that the use of unsuitable reference genes as internal controls could change the target gene’s expression pattern. This study can provide guidance for choosing reference genes for analyzing the expression pattern of target genes and facilitate the functional genomic investigation on fruit body formation and development, as well as stipe elongation and pileus expansion in A. blazei.

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.

Identification of hydrophobin genes and their physiological functions related to growth and development in Pleurotus ostreatus

Hydrophobins are small secreted proteins with important physiological functions and potential applications. Here, Pleurotus ostreatus hydrophobin genes were systematically analyzed: they were characterized, classified, and their expression profiles and gene functions were explored. In total, 40 P. ostreatus hydrophobin genes were found and showed genetic diversity, of which 15 were newly identified. The hydrophobin protein sequences were diverse but all contained eight cysteine residues with a conserved spacing pattern, and 33 of them were class I hydrophobins. The expression profile analyses showed that Vmh3 and Hydph20 were abundant in monokaryotic and dikaryotic mycelia, whereas Hydph17, Po.hyd16, Hydph8 were specifically expressed in monokaryotic mycelia and Po.hyd10 were specific in dikaryotic mycelia. Furthermore, Vmh3, Hydph20, Po.hyd7, and Po.hyd10 were abundant when dikaryotic mycelia cultivated on PDA, which are different from on substrate (Vmh2, Vmh3, Hydph7, Po.hyd3, Po.hyd7, Po.hyd9); Hydph12, POH1, and Po.hyd4 can be induced by natural light and cold stimulation during development from mycelia to primordia; Vmh3, FBH1, Hydph12, Po.hyd1-Po.hyd5, and Po.hyd8 were highly expressed in primordia and young fruiting bodies; Hydph12, Po.hyd1, Po.hyd4, and Po.hyd5 were specifically expressed in pilei. In addition, RNAi transformants of FBH1 exhibited slower growth rates and had fewer primordia and fruiting bodies, which suggests FBH1 affects the growth rate and primordia formation of P. ostreatus. Therefore, P. ostreatus hydrophobin genes belong to a large family and are temporally and spatially expressed to meet the developmental needs of mushroom.

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.

Chromosome-level de novo assembly of Coprinopsis cinerea A43mut B43mut pab1-1 #326 and genetic variant identification of mutants using Nanopore MinION sequencing

The homokaryotic Coprinopsis cinerea strain A43mut B43mut pab1-1 #326 is a widely used experimental model for developmental studies in mushroom-forming fungi. It can grow on defined artificial media and complete the whole lifecycle within two weeks. The mutations in mating type factors A and B result in the special feature of clamp formation and fruiting without mating. This feature allows investigations and manipulations with a homokaryotic genetic background. Current genome assembly of strain #326 was based on short-read sequencing data and was highly fragmented, leading to the bias in gene annotation and downstream analyses. Here, we report a chromosome-level genome assembly of strain #326. Oxford Nanopore Technology (ONT) MinION sequencing was used to get long reads. Illumina short reads was used to polish the sequences. A combined assembly yield 13 chromosomes and a mitochondrial genome as individual scaffolds. The assembly has 15,250 annotated genes with a high synteny with the C. cinerea strain Okayama-7 #130. This assembly has great improvement on contiguity and annotations. It is a suitable reference for further genomic studies, especially for the genetic, genomic and transcriptomic analyses in ONT long reads. Single nucleotide variants and structural variants in six mutagenized and cisplatin-screened mutants could be identified and validated. A 66 bp deletion in Ras GTPase-activating protein (RasGAP) was found in all mutants. To make a better use of ONT sequencing platform, we modified a high-molecular-weight genomic DNA isolation protocol based on magnetic beads for filamentous fungi. This study showed the use of MinION to construct a fungal reference genome and to perform downstream studies in an individual laboratory. An experimental workflow was proposed, from DNA isolation and whole genome sequencing, to genome assembly and variant calling. Our results provided solutions and parameters for fungal genomic analysis on MinION sequencing platform.

Comparative genomics of pyrophilous fungi reveals a link between fire events and developmental genes

Forest fires generate a large amount of carbon that remains resident on the site as dead and partially 'pyrolysed' (i.e. burnt) material that has long residency times and constitutes a significant pool in fire-prone ecosystems. In addition, fire-induced hydrophobic soil layers, caused by condensation of pyrolysed waxes and lipids, increase post-fire erosion and can lead to long-term productivity losses. A small set of pyrophilous fungi dominate post-fire soils and are likely to be involved with the degradation of all these compounds, yet almost nothing is currently known about what these fungi do or the metabolic processes they employ. In this study, we sequenced and analysed genomes from fungi isolated after Rim fire near Yosemite National Park in 2013 and showed the enrichment/expansion of CAZymes and families known to be involved in fruiting body initiation when compared to other basidiomycete fungi. We found gene families potentially involved in the degradation of the hydrophobic layer and pyrolysed organic matter, such as hydrophobic surface binding proteins, laccases (AA1_1), xylanases (GH10, GH11), fatty acid desaturases and tannases. Thus, pyrophilous fungi are important actors to restate the soil's functional capabilities.

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.

Selection markers for transformation of the sequenced reference monokaryon Okayama 7/#130 and homokaryon AmutBmut of Coprinopsis cinerea

BACKGROUND

Two reference strains have been sequenced from the mushroom Coprinopsis cinerea, monokaryon Okayama 7/#130 (OK130) and the self-compatible homokaryon AmutBmut. An adenine-auxotrophy in OK130 (ade8-1) and a para-aminobenzoic acid (PABA)-auxotrophy in AmutBmut (pab1-1) offer selection markers for transformations. Of these two strains, homokaryon AmutBmut had been transformed before to PABA-prototrophy and with the bacterial hygromycin resistance marker hph, respectively.

RESULTS

Gene ade8 encodes a bifunctional enzyme with an N-terminal glycinamide ribonucleotide synthase (GARS) and a C-terminal aminoimidazole ribonucleotide synthase (AIRS) domain required for steps 2 and 5 in the de novo biosynthesis of purines, respectively. In OK130, a missense mutation in ade8-1 rendered residue N231 for ribose recognition by the A loop of the GARS domain into D231. The new ade8 + vector pCcAde8 complements the auxotrophy of OK130 in transformations. Transformation rates with pCcAde8 in single-vector and co-transformations with ade8 +-selection were similarly high, unlike for trp1 + plasmids which exhibit suicidal feedback-effects in single-vector transformations with complementation of tryptophan synthase defects. As various other plasmids, unselected pCcAde8 helped in co-transformations of trp1 strains with a trp1 +-selection vector to overcome suicidal effects by transferred trp1 +. Co-transformation rates of pCcAde8 in OK130 under adenine selection with nuclear integration of unselected DNA were as high as 80% of clones. Co-transformation rates of expressed genes reached 26-42% for various laccase genes and up to 67% with lcc9 silencing vectors. The bacterial gene hph can also be used as another, albeit less efficient, selection marker for OK130 transformants, but with similarly high co-transformation rates. We further show that the pab1-1 defect in AmutBmut is due to a missense mutation which changed the conserved PIKGT motif for chorismate binding in the C-terminal PabB domain to PIEGT in the mutated 4-amino-4-deoxychorismate synthase.

CONCLUSIONS

ade8-1 and pab1-1 auxotrophic defects in C. cinerea reference strains OK130 and AmutBmut for complementation in transformation are described. pCcAde8 is a new transformation vector useful for selection in single and co-transformations of the sequenced monokaryon OK130 which was transformed for the first time. The bacterial gene hph can also be used as an additional selection marker in OK130, making in combination with ade8 + successive rounds of transformation possible.

Unmatched Level of Molecular Convergence among Deeply Divergent Complex Multicellular Fungi

Convergent evolution is pervasive in nature, but it is poorly understood how various constraints and natural selection limit the diversity of evolvable phenotypes. Here, we analyze the transcriptome across fruiting body development to understand the independent evolution of complex multicellularity in the two largest clades of fungi-the Agarico- and Pezizomycotina. Despite >650 My of divergence between these clades, we find that very similar sets of genes have convergently been co-opted for complex multicellularity, followed by expansions of their gene families by duplications. Over 82% of shared multicellularity-related gene families were expanding in both clades, indicating a high prevalence of convergence also at the gene family level. This convergence is coupled with a rich inferred repertoire of multicellularity-related genes in the most recent common ancestor of the Agarico- and Pezizomycotina, consistent with the hypothesis that the coding capacity of ancestral fungal genomes might have promoted the repeated evolution of complex multicellularity. We interpret this repertoire as an indication of evolutionary predisposition of fungal ancestors for evolving complex multicellular fruiting bodies. Our work suggests that evolutionary convergence may happen not only when organisms are closely related or are under similar selection pressures, but also when ancestral genomic repertoires render certain evolutionary trajectories more likely than others, even across large phylogenetic distances.

Carbon metabolism and transcriptome in developmental paths differentiation of a homokaryotic Coprinopsis cinerea strain

The balance and interplay between sexual and asexual reproduction is one of the most intriguing mysteries in the study of fungi. The choice of developmental strategy reflects the ability of fungi to adapt to the changing environment. However, the evolution of developmental paths and the metabolic regulation during differentiation and morphogenesis are poorly understood. Here, an analysis was performed of carbohydrate metabolism and gene expression regulation during the early differentiation process from the vegetative mycelium, to the differentiated structures, fruiting body, oidia and sclerotia, of a homokaryotic fruiting Coprinopsis cinerea strain A43mutB43mut pab1-1 #326. Changes during morphogenesis and the evolution of developmental strategies were followed. Conversion between glucose and glycogen and between glucose and beta-glucan were the main carbon flows in the differentiation processes. Genes related to carbohydrate transport and metabolism were significantly differentially expressed among paths. Sclerotia displayed a set of specifically up-regulated genes that were enriched in the carbon metabolism and energy production and conversion processes. Evolutionary transcriptomic analysis of four developmental paths showed that all transcriptomes were under the purifying selection, and the more stressful the environment, the younger the transcriptome age. Oidiation has the lowest value of transcriptome age index (TAI) and transcriptome divergence index (TDI), while the fruiting process has the highest of both indexes. These findings provide new insights into the regulations of carbon metabolism and gene expressions during the early stages of fungal developmental paths differentiation, and improve our understanding of the evolutionary process of life history and reproductive strategy in fungi.

Phylogenetic origins and family classification of typhuloid fungi, with emphasis on Ceratellopsis, Macrotyphula and Typhula (Basidiomycota)

Typhuloid fungi are a very poorly known group of tiny clavarioid homobasidiomycetes. The phylogenetic position and family classification of the genera targeted here, Ceratellopsis, Macrotyphula, Pterula sensu lato and Typhula, are controversial and based on unresolved phylogenies. Our six-gene phylogeny with an expanded taxon sampling shows that typhuloid fungi evolved at least twice in the Agaricales (Pleurotineae, Clavariineae) and once in the Hymenochaetales. Macrotyphula, Pterulicium and Typhula are nested within the Pleurotineae. The type of Typhula (1818) and Sclerotium (1790), T. phacorrhiza and S. complanatum (synonym T. phacorrhiza), are encompassed in the Macrotyphula clade that is distantly related to a monophyletic group formed by species usually assigned to Typhula. Thus, the correct name for Macrotyphula (1972) and Typhula is Sclerotium and all Typhula species but those in the T. phacorrhiza group need to be transferred to Pistillaria (1821). To avoid undesirable nomenclatural changes, we suggest to conserve Typhula with T. incarnata as type. Clavariaceae is supported as a separate, early diverging lineage within Agaricales, with Hygrophoraceae as a successive sister taxon to the rest of the Agaricales. Ceratellopsis s. auct. is polyphyletic because C. acuminata nests in Clavariaceae and C. sagittiformis in the Hymenochaetales. Ceratellopsis is found to be an earlier name for Pterulicium, because the type, C. queletii, represents Pterulicium gracile (synonym Pterula gracilis), deeply nested in the Pterulicium clade. To avoid re-combining a large number of names in Ceratellopsis we suggest to conserve it with C. acuminata as type. The new genus Bryopistillaria is created to include C. sagittiformis. The families Sarcomyxaceae and Phyllotopsidaceae, and the suborder Clavariineae, are described as new. Six new combinations are proposed and 15 names typified.

Genome-wide mRNA and miRNA analysis in the early stages of germ tube outgrowth in Coprinopsis cinerea

Coprinopsis cinerea is a model mushroom-forming basidiomycete which produces basidiospores during sexual reproduction. This fungus is widely used to study fruiting body formation and development. Molecular mechanisms controlling its growth from vegetative mycelium to multicellular mature fruiting body have been studied extensively. However, little is known about the underlying biological processes during germ tube outgrowth or the transition from basidiospores to multinucleate hyphae. To better understand sexual spore germination in fungi, here we examined the time-dependent cellular events at resting, germinating and fully germinated basidiospores of C. cinerea by genome-wide transcriptional and post-transcriptional analyses and by carbohydrate composition analysis. Our results revealed a high demand of protein degradation, and biosynthesis of various compounds at the early stage of basidiospore gemination and dynamic changes of carbohydrate metabolism throughout the germination process. Seven microRNA-like RNAs (milRNAs) were identified in the resting basidiospores of C. cinerea, six of which were basidiospore-specific. Glycogen and trehalose were shown to be the carbon sources supporting the initiation of germ tube outgrowth. One basidiospore-specific milRNA, cci-milR-37, was found to be a potential regulator of glycogen metabolic pathways related to vegetative hyphal growth. Our results demonstrated the mRNA and miRNA-mediated regulation on energy production, protein and carbohydrate metabolisms at the early developmental stages of germ tube to form totipotent hyphae. To our knowledge, this is the first study to show the roles of miRNAs in mushroom basidiospore germination and out-growth.

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.

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.