Life history and developmental processes in the basidiomycete Coprinus cinereus

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.

Enhanced extracellular production of laccase in Coprinopsis cinerea by silencing chitinase gene

Laccase, a copper-containing polyphenol oxidase, is an important green biocatalyst. In this study, Laccase Lcc5 was homologous recombinantly expressed in Coprinopsis cinerea and a novel strategy of silencing chitinase gene expression was used to enhance recombinant Lcc5 extracellular yield. Two critical chitinase genes, ChiEn1 and ChiE2, were selected by analyzing the transcriptome data of C. cinerea FA2222, and their silent expression was performed by RNA interference (RNAi). It was found that silencing either ChiEn1 or ChiE2 reduced sporulation and growth rate, and increased cell wall sensitivity, but had no significant effect on mycelial branching. Among them, the extracellular laccase activity of the ChiE2-silenced engineered strain Cclcc5-antiChiE2-5 and the control Cclcc5-13 reached the highest values (38.2 and 25.5 U/mL, respectively) at 250 and 150 rpm agitation speeds, corresponding to productivity of 0.35 and 0.19 U/mL·h, respectively, in a 3-L fermenter culture. Moreover, since Cclcc5-antiChiE2-5 could withstand greater shear forces, its extracellular laccase activity was 2.6-fold higher than that of Cclcc5-13 when the agitation speed was all at 250 rpm. To our knowledge, this is the first report of enhanced recombinant laccase production in C. cinerea by silencing the chitinase gene. This study will pave the way for laccase industrial production and accelerate the development of a C. cinerea high-expression system. KEY POINTS: • ChiEn1 and ChiE2 are critical chitinase genes in C. cinerea FA2222 genome. • Chitinase gene silencing enhanced the tolerance of C. cinerea to shear forces. • High homologous production of Lcc5 is achieved by fermentation in a 3-L fermenter.

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.

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.

Exposure to polystyrene nanoplastics reduces bacterial and fungal biomass in microfabricated soil models

Nanoplastics have been proven to induce toxicity in diverse organisms, yet their effect on soil microbes like bacteria and fungi remains largely unexplored. In this paper, we used micro-engineered soil models to investigate the effect of polystyrene (PS) nanospheres on Pseudomonas putida and Coprinopsis cinerea. Specifically, we explored the effects of increasing concentrations of 60 nm carboxylated bovine serum albumin (BSA) coated nanospheres (0, 0.5, 2, and 10 mg/L) on these bacterial and fungal model organisms respectively, over time. We found that both microorganisms could disperse through the PS solution, but long-distance dispersal was reduced by high concentrations. Microbial biomass decreased in all treatments, in which bacteria showed a linear dose response with the strongest effect at 10 mg/L concentration, and fungi showed a non-linear response with the strongest effect at 2 mg/L concentration. At the highest nanoplastics concentration, the first colonizing fungal hyphae adsorbed most of the PS nanospheres present in their vicinity, in a process that we termed the 'vacuum cleaner effect'. As a result, the toxicity effect of the original treatment on subsequently growing fungal hyphae was reduced to a growth level indistinguishable from the control. We did not find evidence that nanoplastics are able to penetrate bacterial nor fungal cell walls. Overall, our findings provide evidence that nanoplastics can cause a direct negative effect on soil microbes and highlight the need for further studies that can explain how the microbial stress response might affect soil functions.

The first suspected disseminated Hormographiella aspergillata infection in China, diagnosed using metagenomic next-generation sequencing: a case report and literature review

Hormographiella aspergillata is a rare and emerging cause of invasive mould infections in patients with haematological malignancies, with a mortality rate of approximately 70%. Here, we present the first reported case of suspected disseminated H. aspergillata infection in China. The patient experienced a second relapse of acute myeloid leukaemia and developed neutropenia, fever, discrepant blood pressure between limbs, and cutaneous lesions limited to the left upper extremity. Since lung tissue biopsy was not feasible, metagenomic next-generation sequencing (mNGS) and panfungal polymerase chain reaction (PCR) analysis of bronchoalveolar lavage fluid and blood samples were performed, which indicated probable H. aspergillata pulmonary infection. Histopathology of cutaneous lesions revealed numerous fungal hyphae within dermal blood vessels. mNGS of a skin biopsy sample identified H. aspergillata sequences, and the fungi was subsequently recovered from fungal culture, proving cutaneous H. aspergillata infection. Despite combined antifungal therapy, the patient died owing to disease progression. Additionally, 22 previously reported cases of invasive H. aspergillata infection were reviewed in patients with haematological malignancies. Thus, mNGS is a powerful diagnostic tool for the early and effective detection of invasive H. aspergillata infections, with the advantage of sequencing all potential pathogens, and providing results within 24 h.

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.

AfLaeA, a Global Regulator of Mycelial Growth, Chlamydospore Production, Pathogenicity, Secondary Metabolism, and Energy Metabolism in the Nematode-Trapping Fungus Arthrobotrys flagrans

Arthrobotrys flagrans (Duddingtonia flagrans) is a typical nematode-trapping fungus which has been used for nematode biocontrol. The global regulator LaeA is widely distributed in filamentous fungi and plays a crucial role in secondary metabolism and development in addition to pathogenicity in fungal pathogens. In this study, the chromosome-level genome of A. flagrans CBS 565.50 was sequenced and homologous sequences of LaeA were identified in A. flagrans. A. flagrans LaeA (AfLaeA) knockout resulted in slower hyphal growth and a smoother hyphal surface. Importantly, deletion of AfLaeA resulted in the absence of chlamydospores and attenuated glycogen and lipid accumulation in hyphae. Similarly, disruption of the AfLaeA gene led to fewer traps and electron-dense bodies, lower protease activity, and a delay in capturing nematodes. The AfLaeA gene had a large effect on the secondary metabolism of A. flagrans, and both the deletion and overexpression of AfLaeA could yield new compounds, whereas some compounds were lost due to the absence of the AfLaeA. Protein-protein interactions between AfLaeA and another eight proteins were detected. Furthermore, transcriptome data analysis showed that 17.77% and 35.51% of the genes were influenced by the AfLaeA gene on days 3 and 7, respectively. AfLaeA gene deletion resulted in the higher expression level of the artA gene cluster, and multiple differentially expressed genes involved in glycogen and lipid synthesis and metabolism showed opposite expression patterns in wild-type and ΔAfLaeA strains. In summary, our results provide novel insights into the functions of AfLaeA in mycelial growth, chlamydospore production, pathogenicity, secondary metabolism, and energy metabolism in A. flagrans. IMPORTANCE The regulation of biological functions, such as the secondary metabolism, development, and pathogenicity of LaeA, has been reported in multiple fungi. But to date, no study on LaeA in nematode-trapping fungi has been reported. Moreover, it has not been investigated whether or not LaeA is involved in energy metabolism and chlamydospore formation has not been investigated. Especially in the formation mechanism of chlamydospores, several transcription factors and signaling pathways are involved in the production of chlamydospores, but the mechanism of chlamydospore formation from an epigenetic perspective has not been revealed. Concurrently, an understanding of protein-protein interactions will provide a broader perspective on the regulatory mechanism of AfLaeA in A. flagrans. This finding is critical for understanding the regulatory role of AfLaeA in the biocontrol fungus A. flagrans and establishes a foundation for developing high-efficiency nematode biocontrol agents.

Multi-stage nuclear transcriptomic insights of morphogenesis and biparental role changes in Lentinula edodes

Based on six offspring with different mitochondrial (M) and parental nuclear (N) genotypes, the multi-stage morphological characteristics and nuclear transcriptomes of Lentinula edodes were compared to investigate morphogenesis mechanisms during cultivation, the key reason for cultivar resistance to genotype changes, and regulation related to biparental role changes. Six offspring had specific transcriptomic data and morphological characteristics that were mainly regulated by the two parental nuclei, followed by the cytoplasm, at different growth stages. Importing a wild N genotype easily leads to failure or instability of fruiting; however, importing wild M genotypes may improve cultivars. Major facilitator superfamily (MFS) transporter genes encoding specific metabolites in spawns may play crucial roles in fruiting body formation. Pellets from submerged cultivation and spawns from sawdust substrate cultivation showed different carbon metabolic pathways, especially in secondary metabolism, degradation of lignin, cellulose and hemicellulose, and plasma membrane transport (mainly MFS). When the stage of small young pileus (SYP) was formed on the surface of the bag, the spawns inside were mainly involved in nutrient accumulation. Just broken pileus (JBP) showed a different expression of plasma membrane transporter genes related to intracellular material transport compared to SYP and showed different ribosomal proteins and cytochrome P450 functioning in protein biosynthesis and metabolism than near spreading pileus (NSP). Biparental roles mainly regulate offspring metabolism, growth, and morphogenesis by differentially expressing specific genes during different vegetative growth stages. Additionally, some genes encoding glycine-rich RNA-binding proteins, F-box, and folliculin-interacting protein repeat-containing proteins may be related to multi-stage morphogenesis. KEY POINTS: • Replacement of nuclear genotype is not suitable for cultivar breeding of L. edodes. • Some genes show a biparental role-divergent expression at mycelial growth stage. • Transcriptomic changes of some sawdust substrate cultivation stages have been elucidated.

Chemical screening identifies the anticancer properties of Polyporous tuberaster

Most conventional anticancer drugs cause resistance to chemotherapy, which has emerged as one of the major obstacles to cancer treatment. In order to address this issue, efforts have been made to select new anticancer compounds from natural sources. The aim of this study is to identify novel anticancer compounds from mycelial culture extracts belonging to Polyporus tuberaster (P. tuberaster). Here, we found that mycelial culture extracts of P. tuberaster cultured in PDB medium (pt-PDB) effectively inhibited cancer cell growth. pt-PDB reduced the growth of cancer cells through apoptosis induction and S-phase arrest. The anticancer efficacy of pt-PDB was not to limited to one type of cancer. Furthermore, unlike traditional anticancer medications, pt-PDB did not increase the proportion of side population (SP) cells, which plays a key role in the development of chemoresistance. Taken together, we discovered a novel anticancer drug candidate that has anticancer properties without increasing the proportion of SP cells. This new drug candidate can be used for the treatment of cancer, especially chemoresistant malignancies, and will provide a breakthrough in the treatment of chemoresistant cancer.

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.

Overexpression of the Capebp2 Gene Encoding the PEBP-like Protein Promotes the Cap Redifferentiation in Cyclocybe aegerita

Phosphatidylethanolamine-binding protein (PEBP) is widely involved in various physiological behaviors, such as the transition from vegetative growth to reproductive growth in plants, tumorigenesis in the human, etc. However, few functional studies have examined pebp genes affecting the development of fungi. In this study, Capebp2 was cloned from Cyclocybe aegerita AC0007 strains based on the genome sequence and gene prediction, and the sequence alignment of CaPEBP2 with other PEBP proteins from other biological sources including plant, animal, fungi, and bacteria indicated that PEBP had low sequence similarity in fungi, whereas all protein sequences had some conserved motifs such as DPDAP and HRY. Expression analysis showed the transcription level of Capebp2 increased approximately 20-fold in fruiting bodies compared with mycelia. To uncover the function of Capebp2 in C. aegetita development, Capebp2 was cloned into a pATH vector driven by the actin promoter for obtaining overexpression transformant lines. Fruiting experiments showed the transformed strains overexpressing Capebp2 exhibited redifferentiation of the cap on their surface, including intact fruiting bodies or partial lamella during fruiting development stage, and the longitudinal section indicated that all regenerated bodies or lamella sprouted from the flesh and shared the epidermis with the mother fruiting bodies. In summary, the sequence characterization of Capebp2, expression level during different development stages, and function on fruiting body development were documented in this study, and these findings provided a reference to study the role of pebp in the development process of basidiomycetes. Importantly, gene mining of pebp, function characterization, and the regulating pathways involved need to be uncovered in further studies.

Chemical screening identifies the anticancer properties of Polyporous parvovarius

One of the biggest obstacles in cancer treatment is the development of chemoresistance. To overcome this, attempts have been made to screen novel anticancer substances derived from natural products. The purpose of this study is to find new anticancer candidates in the mycelium culture extract of mushrooms belonging to Polyporus. Here, we used a high-throughput screening to find agents capable of inhibiting cancer cell proliferation. The culture extract of Polyporus Parvovarius mycelium in DY medium (pp-DY) was effective. pp-DY inhibited cancer cell proliferation by inducing apoptosis and S-phase arrest. The anticancer property of pp-DY was not only effective against one type of cancer, but also against another type of cancer. Compound fractionation was performed, and the active ingredient exhibiting anticancer effects in pp-DY was identified as 3,4-dihydroxybenzaldehyde (Protocatechualdehyde, PCA). PCA, like pp-DY, inhibited the proliferation of cancer cells by inducing apoptosis and S-phase arrest. Furthermore, unlike conventional anticancer drugs, PCA did not increase the proportion of the side population that plays the most important role in the development of chemoresistance. Taken together, our data revealed the novel mycelium culture extract that exhibited anticancer property, and identified active ingredients that did not activate a proportion of the side population. These novel findings may have clinical applications in the treatment of cancer, particularly chemo-resistant cancer.

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.

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.

Organic Nitrogen Supplementation Increases Vegetative and Reproductive Biomass in a Versatile White Rot Fungus

The Black Poplar Mushroom Cyclocybe aegerita (syn. Agrocybe aegerita) is a white-rot fungus that naturally fruits from woody substrates, including buried wood. It is known for its substrate versatility and is equipped with a respective carbohydrate-active enzyme repertoire being intermediate between typical white-rot fungi and plant litter decomposers. Given relative nitrogen scarcity in wood, mobilization of nitrogen from surrounding litter is known as a way to meet nitrogen requirements for cellular homeostasis and reproduction of wood decay fungi. However, the effect of added nitrogen on vegetative and reproductive biomass has not yet been studied in a uniform minimalistic laboratory setup. For C. aegerita, such a growth and fruiting setup has been developed. In the present study, this white-rot fungus has been grown with and without additional β-adenosine, an organic nitrogen source present in plant litter. Elevated β-adenosine levels increased aerial mycelium weight by 30% (1 × β-adenosine) and 55% (10 × β-adenosine), reproductive biomass by 75% (1 × β-adenosine) and by 100% (10 × β-adenosine), number of primordia by 127% (10 × β-adenosine) and accelerated primordium formation by 1.6 days (10 × β-adenosine), compared to the control treatment. These findings imply that C. aegerita invests additional organic nitrogen resources into direct vegetative and reproductive biomass build-up at the same time. Colonization of niches with accessory nitrogen sources, like buried wood, which is near the plant litter layer, may thus provide an evolutionary fitness advantage. Globally anthropogenically altered nitrogen dynamics may affect hyphal-driven processes as well as fruit body-driven food webs.

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.

Mitochondrial Effects on the Physiological Characteristics of Lentinula edodes

In the mating of filamentous basidiomycetes, dikaryotic mycelia are generated through the reciprocal movement of nuclei to a monokaryotic cytoplasm where a nucleus of compatible mating type resides, resulting in the establishment of two different dikaryotic strains having the same nuclei but different mitochondria. To better understand the role of mitochondria in mushrooms, we created four sets of dikaryotic strains of Lentinula edodes, including B2 × E13 (B2 side) and B2 × E13 (E13 side), B5 × E13 (B5 side) and B5 × E13 (E13 side), E8 × H3 (E8 side) and E8 × H3 (H3 side), and K3 × H3 (K3 side) and K3 × H3 (H3 side). The karyotypes and mitochondrial types of the dikaryotic strains were successfully identified by the A mating type markers and the mitochondrial variable length tandem repeat markers, respectively. Comparative analyses of the dikaryotic strains on the mycelial growth, substrate browning, fruiting characteristics, and mitochondrial gene expression revealed that certain mitochondria are more effective in the mycelial growth and the production of fruiting body, possibly through the activated energy metabolism. Our findings indicate that mitochondria affect the physiology of dikaryotic strains having the same nuclear information and therefore a selection strategy aimed at mitochondrial function is needed in the development of new mushroom strain.

Coprinopsis cinerea dioxygenase is an oxygenase forming 10(S)-hydroperoxide of linoleic acid, essential for mushroom alcohol, 1-octen-3-ol synthesis

1-Octen-3-ol is a volatile oxylipin found ubiquitously in Basidiomycota and Ascomycota. The biosynthetic pathway forming 1-octen-3-ol from linoleic acid via the linoleic acid 10(S)-hydroperoxide was characterized 40 years ago in mushrooms, yet the enzymes involved are not identified. The dioxygenase 1 and 2 genes (Ccdox1 and Ccdox2) in the mushroom Coprinopsis cinerea contain an N-terminal cyclooxygenase-like heme peroxidase domain and a C-terminal cytochrome P450-related domain. Herein, we show that recombinant CcDOX1 is responsible for dioxygenation of linoleic acid to form the 10(S)-hydroperoxide, the first step in 1-octen-3-ol synthesis, whereas CcDOX2 conceivably forms linoleate 8-hydroperoxide. We demonstrate that knockout of the Ccdox1 gene suppressed 1-octen-3-ol synthesis, although added linoleic acid 10(S)-hydroperoxide was still efficiently converted. The P450-related domain of CcDOX1 lacks the characteristic Cys heme ligand and the evidence indicates that a second uncharacterized enzyme converts the 10(S)-hydroperoxide to 1-octen-3-ol. Additionally, we determined the gene knockout strain (ΔCcdox1) was less attractive to fruit fly larvae, while the feeding behavior of fungus gnats on ΔCcdox1 mycelia showed little difference from that on the mycelia of the wild-type strain. The proliferation of fungivorous nematodes on ΔCcdox1 mycelia was similar to or slightly worse than that on wild-type mycelia. Thus, 1-octen-3-ol seems to be an attractive compound involved in emitter-receiver ecological communication in mushrooms.

A novel gene, Le-Dd10, is involved in fruiting body formation of Lentinula edodes

The cDNA library prepared from Lentinula edodes, Hokken 600 (H600), primordia was screened using cDNA expressed specifically in Dictyostelium discoideum prestalk as a probe. Twenty-one clones, Le-Dd1 ~ 21, were isolated from the L. edodes primordia cDNA library. Functional analysis of each gene was carried out by transformation into protoplast cells from L. edodes Mori 252 (M252) mycelia with the overexpression vector pLG-RasF1 of each gene because M252 protoplast cells were transformed with an 11-fold higher efficiency than H600 cells. Transformants with the overexpression vector of Le-Dd10 formed a fruiting body at almost the same time as H600, a positive control, although M252, a negative control, did not form a fruiting body under culture conditions. This suggested that Le-Dd10 is involved in the formation of fruiting bodies. Single-strand conformation polymorphism analysis revealed that Le-Dd10 is located on No. 4 linkage group of L. edodes. The properties of Le-Dd10 products were investigated by Western blotting analysis using polyclonal antibodies against GST:Le-Dd10 fusion proteins. As a result, 56-kDa, 27-kDa, and 14-kDa protein bands appeared in primordial and fruiting body stages, although the expected molecular weight of the Le-Dd10 product was 50 kDa.

Nuclear-specific gene expression in heterokaryons of the filamentous ascomycete Neurospora tetrasperma

Heterokaryosis is a system in which genetically distinct nuclei coexist within the same cytoplasm. While heterokaryosis dominates the life cycle of many fungal species, the transcriptomic changes associated with the transition from homokaryosis to heterokaryosis is not well understood. Here, we analyse gene expression profiles of homokaryons and heterokaryons from three phylogenetically and reproductively isolated lineages of the filamentous ascomycete Neurospora tetrasperma. We show that heterokaryons are transcriptionally distinct from homokaryons in the sexual stage of development, but not in the vegetative stage, suggesting that the phenotypic switch to fertility in heterokaryons is associated with major changes in gene expression. Heterokaryon expression is predominantly defined by additive effects of its two nuclear components. Furthermore, allele-specific expression analysis of heterokaryons with varying nuclear ratios show patterns of expression ratios strongly dependent on nuclear ratios in the vegetative stage. By contrast, in the sexual stage, strong deviations of expression ratios indicate a co-regulation of nuclear gene expression in all three lineages. Taken together, our results show two levels of expression control: additive effects suggest a nuclear level of expression, whereas co-regulation of gene expression indicate a heterokaryon level of control.

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

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

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 × 10⁸ 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.

Systematic analysis of changes across different developmental stages of the mushroom Sarcomyxa edulis

Sarcomyxa edulis is a widely harvested mushroom of Northeastern Asia. Its development can be divided into six stages: 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). Differentially expressed gene (DEG) analysis and weighted gene coexpression network analysis (WGCNA) are important bioinformatic methods for screening key genes. To explore the growth and development mechanisms of the mushroom S. edulis and clarify its genetic background, DEG and WGCNA analyses were combined to screen key genes at different developmental stages. From A1 to A6, respectively, 459, 97, 885, 169, 277, and 712 key genes were identified. Then the Gene Ontology (GO) terms and KEGG pathways of key genes were analyzed, and GO and KEGG analyses were performed on all genes across different periods using GSEA. In summary, the genes in A1 were mainly involved in amino sugar and nucleotide sugar metabolism, structural molecule activity, and oxidative phosphorylation. At the A2 stage, genes were mainly involved in peptidase activity, peroxidase activity, oxidoreductase activity, antioxidant activity, biosynthesis of secondary metabolites, and glycolysis and gluconeogenesis. A3 genes were involved in gene expression, RNA metabolism, spliceosome, RNA transport, and ribosome biogenesis. A4 genes were mainly involved in the biosynthesis of secondary metabolites, proteasome complex, cellular protein complex assembly, actin filament-based processes, oxidative phosphorylation, and carbon metabolism. The A5 stage genes were involved in the carbohydrate metabolic process, polysaccharide metabolic process, and the biosynthesis of secondary metabolites, leucine, isoleucine, and ABC transporters. Finally, A6 genes were mainly involved in the cell cycle, meiosis of yeast, MAPK signaling pathway, cellular response to DNA damage stimulus, DNA metabolic process, DNA replication, and DNA repair. The combination of multiple analyses provides us with an in-depth understanding of the network that regulates mushroom development.

MAPK CcSakA of the HOG Pathway Is Involved in Stipe Elongation during Fruiting Body Development in Coprinopsis cinerea

Mitogen-activated protein kinase (MAPK) pathways, such as the high-osmolarity glycerol mitogen-activated protein kinase (HOG) pathway, are evolutionarily conserved signaling modules responsible for transmitting environmental stress signals in eukaryotic organisms. Here, we identified the MAPK homologue in the HOG pathway of Coprinopsis cinerea, which was named CcSakA. Furthermore, during the development of the fruiting body, CcSakA was phosphorylated in the fast elongating apical part of the stipe, which meant that CcSakA was activated in the apical elongating stipe region of the fruiting body. The knockdown of CcSakA resulted in a shorter stipe of the fruiting body compared to the control strain, and the expression of phosphomimicking mutant CcSakA led to a longer stipe of the fruiting body compared to the control strain. The chitinase CcChiE1, which plays a key role during stipe elongation, was downregulated in the CcSakA knockdown strains and upregulated in the CcSakA phosphomimicking mutant strains. The results indicated that CcSakA participated in the elongation of stipes in the fruiting body development of C. cinerea by regulating the expression of CcChiE1. Analysis of the H₂O₂ concentration in different parts of the stipe showed that the oxidative stress in the elongating part of the stipe was higher than those in the non-elongating part. The results indicated that CcSakA of the HOG pathway may be activated by oxidative stress. Our results demonstrated that the HOG pathway transmits stress signals and regulates the expression of CcChiE1 during fruiting body development in C. cinerea.

Are Mushrooms Parametric?

Designing with biological materials as a burgeoning approach in the architecture field requires the development of new design strategies and fabrication methods. In this paper, we question if designers can use a parametric design approach while working with living materials. The research uses fungi as a biomaterial probe to experiment with the parametric behavior of living systems. Running design experiments using fungi helps to understand the extent to which biological systems can be considered parametric and, if so, what kind of parametric systems they are. Answering these questions provides a method to work with complex biological systems and may lead to new approaches of fabricating materials by tuning the environmental parameters of biological growth.

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.

A Velvet Transcription Factor Specifically Activates Mating through a Novel Mating-Responsive Protein in the Human Fungal Pathogen Cryptococcus deneoformans

Sexual reproduction facilitates infection by the production of both a lineage advantage and infectious sexual spores in the ubiquitous human fungal pathogen Cryptococcus deneoformans. However, the regulatory determinants specific for initiating mating remain poorly understood. Here, we identified a velvet family regulator, Cva1, that strongly promotes sexual reproduction in C. deneoformans. This regulation was determined to be specific, based on a comprehensive phenotypic analysis of cva1Δ under 26 distinct in vitro and in vivo growth conditions. We further revealed that Cva1 plays a critical role in the initiation of early mating events, including sexual cell-cell fusion, but is not important for the late sexual development stages or meiosis. Thus, Cva1 specifically contributes to mating activation. Importantly, a novel mating-responsive protein, Cfs1, serves as the key target of Cva1 during mating, since its absence nearly blocks cell-cell fusion in C. deneoformans and its sister species C. neoformans. Together, our findings provide insight into how C. deneoformans ensures the regulatory specificity of mating.

IMPORTANCE The human fungal pathogen C. deneoformans is a model organism for studying fungal sexual reproduction, which is considered to be important to infection. However, the specific regulatory determinants for activation of sexual reproduction remain poorly understood. In this study, by combining transcriptomic and comprehensive phenotypic analysis, we identified a velvet family regulator Cva1 that specifically and critically elicits early mating events, including sexual cell-cell fusion. Significantly, Cva1 induces mating through the novel mating-responsive protein Cfs1, which is essential for cell-cell fusion in C. deneoformans and its sister species C. neoformans. Considering that Cva1 and Cfs1 are highly conserved in species belonging to Cryptococcaeceae, they may play conserved and specific roles in the initiation of sexual reproduction in this important fungal clade, which includes multiple human fungal pathogens.

Identification of volatile producing enzymes in higher fungi: Combining analytical and bioinformatic methods

Filamentous fungi harbor the genetic potential for the biosynthesis of several secondary metabolites including various volatile organic compounds (VOCs). Nonetheless, under standard laboratory conditions, many of these VOCs are not formed. Furthermore, little is known about enzymes involved in the production of fungal VOCs. To tap these interesting topics, we developed an approach to identify enzymes putatively involved in the fungal VOC biosynthesis. In this chapter, we highlight different fungal cultivation methods and techniques for the extraction of VOCs, including a method that allows the noninvasive analysis of VOCs. In addition using terpene synthases as an example, it is depicted how enzymes putatively involved in VOC synthesis can be identified by means of bioinformatic approaches. Transcriptomic data of chosen genes combined with volatilome data obtained during different developmental stages is demonstrated as a powerful tool to identify enzymes putatively involved in fungal VOC biosynthesis. Especially with regard to subsequent enzyme characterization, this procedure is a target-oriented way to save time and efforts by considering only the most important enzymes.

Injection into and extraction from single fungal cells

The direct delivery of molecules and the sampling of endogenous compounds into and from living cells provide powerful means to modulate and study cellular functions. Intracellular injection and extraction remain challenging for fungal cells that possess a cell wall. The most common methods for intracellular delivery into fungi rely on the initial degradation of the cell wall to generate protoplasts, a step that represents a major bottleneck in terms of time, efficiency, standardization, and cell viability. Here, we show that fluidic force microscopy enables the injection of solutions and cytoplasmic fluid extraction into and out of individual fungal cells, including unicellular model yeasts and multicellular filamentous fungi. The approach is strain- and cargo-independent and opens new opportunities for manipulating and analyzing fungi. We also perturb individual hyphal compartments within intact mycelial networks to study the cellular response at the single cell level.

Guillaume-Gentil et al. describe a method that employs a modified AFM tip for selectively sampling from and injecting into individual fungal cells of differing morphology. The authors describe extensive modifications on their system previously used for mammalian cells to overcome many of the challenges associated with working on single fungal cells.

Large-scale fungal strain sequencing unravels the molecular diversity in mating loci maintained by long-term balancing selection

Balancing selection, an evolutionary force that retains genetic diversity, has been detected in multiple genes and organisms, such as the sexual mating loci in fungi. However, to quantify the strength of balancing selection and define the mating-related genes require a large number of strains. In tetrapolar basidiomycete fungi, sexual type is determined by two unlinked loci, MATA and MATB. Genes in both loci define mating type identity, control successful mating and completion of the life cycle. These loci are usually highly diverse. Previous studies have speculated, based on culture crosses, that species of the non-model genus Trichaptum (Hymenochaetales, Basidiomycota) possess a tetrapolar mating system, with multiple alleles. Here, we sequenced a hundred and eighty strains of three Trichaptum species. We characterized the chromosomal location of MATA and MATB, the molecular structure of MAT regions and their allelic richness. The sequencing effort was sufficient to molecularly characterize multiple MAT alleles segregating before the speciation event of Trichaptum species. Analyses suggested that long-term balancing selection has generated trans-species polymorphisms. Mating sequences were classified in different allelic classes based on an amino acid identity (AAI) threshold supported by phylogenetics. 17,550 mating types were predicted based on the allelic classes. In vitro crosses allowed us to support the degree of allelic divergence needed for successful mating. Even with the high amount of divergence, key amino acids in functional domains are conserved. We conclude that the genetic diversity of mating loci in Trichaptum is due to long-term balancing selection, with limited recombination and duplication activity. The large number of sequenced strains highlighted the importance of sequencing multiple individuals from different species to detect the mating-related genes, the mechanisms generating diversity and the evolutionary forces maintaining them.

Response of Fruit Body Assemblage Color Lightness to Macroclimate and Vegetation Cover

Understanding how species relate mechanistically to their environment via traits is a central goal in ecology. Many macroecological rules were found for macroorganisms, however, whether they can explain microorganismal macroecological patterns still requires investigation. Further, whether macroecological rules are also applicable in microclimates is largely unexplored. Here we use fruit body-forming fungi to understand both aspects better. A recent study showed first evidence for the thermal-melanism hypothesis (Bogert’s rule) in fruit body-forming fungi and relied on a continental spatial scale with large grid size. At large spatial extent and grid sizes, other factors like dispersal limitation or local microclimatic variability might influence observed patterns besides the rule of interest. Therefore, we test fungal assemblage fruit body color lightness along a local elevational gradient (mean annual temperature gradient of 7°C) while considering the vegetation cover as a proxy for local variability in microclimate. Using multivariate linear modeling, we found that fungal fruiting assemblages are significantly darker at lower mean annual temperatures supporting the thermal-melanism hypothesis. Further, we found a non-significant trend of assemblage color lightness with vegetation cover. Our results support Bogert’s rule for microorganisms with macroclimate, which was also found for macroorganisms.

Pattern formation features might explain homoplasy: fertile surfaces in higher fungi as an example

Fungi show a high degree of morphological convergence. Regarded for a long time as an obstacle for phylogenetic studies, homoplasy has also been proposed as a source of information about underlying morphogenetic patterning mechanisms. The "local-activation and long-range inhibition principle" (LALIP), underlying the famous reaction-diffusion model proposed by Alan Turing in 1952, appears to be one of the universal phenomena that can explain the ontogenetic origin of seriate patterns in living organisms. Reproductive structures of fungi in the class Agaricomycetes show a highly periodic structure resulting in, for example, poroid, odontoid, lamellate or labyrinthic hymenophores. In this paper, we claim that self-organized patterns might underlie the basic ontogenetic processes of these structures. Simulations based on LALIP-driven models and covering a wide range of parameters show an absolute mutual correspondence with the morphospace explored by extant agaricomycetes. This could not only explain geometric particularities but could also account for the limited possibilities displayed by hymenial configurations, thus making homoplasy a direct consequence of the limited morphospace resulting from the proposed patterning dynamics.

Genetic structure and evolutionary diversity of mating-type (MAT) loci in Hypsizygus marmoreus

The mating compatibility in fungi is generally governed by genes located within a single or two unlinked mating type (MAT) loci. Hypsizygus marmoreus is an edible mushroom in the order Agaricales with a tetrapolar system, which contains two unlinked MAT loci-homeodomain (HD) transcription factor genes and pheromone/pheromone receptor genes (P/R). In this study, we analyzed the genetic structure and diversity of MAT loci in tetrapolar system of H. marmoreus through sequencing of 54 heterokaryon and 8 homokaryon strains. Although within the HD loci, the gene order was conserved, the gene contents were variable, and the HD loci haplotypes were further classified into four types. By analyzing the structure, phylogeny, and the HD transmissibility based on the progeny of these four HD mating-type loci types, we found that they were heritable and tightly linked at the HD loci. The P/R loci genes were found to comprise three pheromone receptors, three pheromones, and two pheromone receptor-like genes. Intra- and inter-specific phylogenetic analyses of pheromone receptors revealed that the STE3 genes were divided into three groups, and we thus theorize that they diverged before speciation. Comparative analysis of the MAT regions among 73 Basidiomycete species indicated that the diversity of HD and P/R loci in Agaricales and Boletales may contribute to mating compatibility. The number of HD genes were not correlated with the tetrapolar or bipolar systems. In H. marmoreus, the expression levels of these genes at HD and P/R loci of compatible strains were found higher than in those of homonuclear/homokaryotic strains, indicating that these mating genes acted as switches for mating processes. Further collinear analysis of HD loci in interspecific species found that HD loci contains conserved recombination hotspots showing major rearrangements in Coprinopsis cinerea and Schizophyllum commune, suggesting different mechanisms for evolution of physically linked MAT loci in these groups. It seems likely that gene rearrangements are common in Agaricales fungi around HD loci. Together, our study provides insights into the genomic basis of mating compatibility in H. marmoreus.

A Transcriptomic Atlas of the Ectomycorrhizal Fungus Laccaria bicolor

Trees are able to colonize, establish and survive in a wide range of soils through associations with ectomycorrhizal (EcM) fungi. Proper functioning of EcM fungi implies the differentiation of structures within the fungal colony. A symbiotic structure is dedicated to nutrient exchange and the extramatricular mycelium explores soil for nutrients. Eventually, basidiocarps develop to assure last stages of sexual reproduction. The aim of this study is to understand how an EcM fungus uses its gene set to support functional differentiation and development of specialized morphological structures. We examined the transcriptomes of Laccaria bicolor under a series of experimental setups, including the growth with Populus tremula x alba at different developmental stages, basidiocarps and free-living mycelium, under various conditions of N, P and C supply. In particular, N supply induced global transcriptional changes, whereas responses to P supply seemed to be independent from it. Symbiosis development with poplar is characterized by transcriptional waves. Basidiocarp development shares transcriptional signatures with other basidiomycetes. Overlaps in transcriptional responses of L. bicolor hyphae to a host plant and N/C supply next to co-regulation of genes in basidiocarps and mature mycorrhiza were detected. Few genes are induced in a single condition only, but functional and morphological differentiation rather involves fine tuning of larger gene sets. Overall, this transcriptomic atlas builds a reference to study the function and stability of EcM symbiosis in distinct conditions using L. bicolor as a model and indicates both similarities and differences with other ectomycorrhizal fungi, allowing researchers to distinguish conserved processes such as basidiocarp development from nutrient homeostasis.

Discovery and Functional Study of a Novel Genomic Locus Homologous to Bα-Mating-Type Sublocus of Lentinula edodes

The interaction of mating pheromone and pheromone receptor from the B mating-type locus is the first step in the activation of the mushroom mating signal transduction pathway. The B mating-type locus of Lentinula edodes is composed of Bα and Bβ subloci, each of which contains genes for mating pheromone and pheromone receptor. Allelic variations in both subloci generate multiple B mating-types through which L. edodes maintains genetic diversity. In addition to the B mating-type locus, our genomic sequence analysis revealed the presence of a novel chromosomal locus 43.3 kb away from the B mating-type locus, containing genes for a pair of mating pheromones (PHBN1 and PHBN2) and a pheromone receptor (RCBN). The new locus (Bα-N) was homologous to the Bα sublocus, but unlike the multiallelic Bα sublocus, it was highly conserved across the wild and cultivated strains. The interactions of RcbN with various mating pheromones from the B and Bα-N mating-type loci were investigated using yeast model that replaced endogenous yeast mating pheromone receptor STE2 with RCBN. The yeast mating signal transduction pathway was only activated in the presence of PHBN1 or PHBN2 in the RcbN producing yeast, indicating that RcbN interacts with self-pheromones (PHBN1 and PHBN2), not with pheromones from the B mating-type locus. The biological function of the Bα-N locus was suggested to control the expression of A mating-type genes, as evidenced by the increased expression of two A-genes HD1 and HD2 upon the treatment of synthetic PHBN1 and PHBN2 peptides to the monokaryotic strain of L. edodes.

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.

The extracellular β-glucosidase BGL2 has two variants with different molecular sizes and hydrolytic activities in the stipe or pilei of Coprinopsis cinerea

Two variants of extracellular β-glucosidase (BGL2) were purified from the stipe and pilei of Coprinopsis cinerea. In the stipe, BGL2 was a monomeric protein with an apparent molecular mass of approximately 220 kDa, representing a mature full-length peptide of BGL2. However, in the pilei, the apparent molecular mass of BGL2 was only approximately 120 kDa, consisting of the 60 kDa N-terminal fragment and 55 kDa C-terminal fragment. The hydrolytic activities of BGL2 purified from the pilei were higher than those of BGL2 purified from the stipe. No mRNA splice variants of bgl2 were detected. Therefore, the different variants of BGL2 in the stipe and pilei were not formed by differential RNA splicing. Furthermore, in vitro experiments showed that full-length BGL2 could be cleaved by endogenous proteases from pilei or commercial trypsin at a similar site to form an oligomeric protein consisting of the N-terminal fragment and C-terminal fragment similar to BGL2 from pilei. The hydrolytic activity of BGL2 increased after cleavage by those proteases in vitro. We conclude that the 120 kDa variant of BGL2 in the pilei of C. cinerea is formed by posttranslational proteolytic cleavage. Posttranslational proteolytic cleavage is an efficient way to regulate the activity of BGL2 to adapt to the needs of different physiological functions in the elongation stipe and expansion pilei of C. cinerea.

Comparison of the Nutritional Properties and Transcriptome Profiling Between the Two Different Harvesting Periods of Auricularia polytricha

Auricularia polytricha (A. polytricha), regarded as an edible and medical mushroom, has attracted toward the research interests because of the high nutrition and bioactivity. The nutritional and medical properties of A. polytricha have been well-studied; however, research about the difference of the nutritional properties and transcriptome profiling between the two different harvesting periods of A. polytricha was limited. In this study, the nutritional properties and transcriptome profiling were compared between the two different harvesting periods of A. polytricha: AP_S1 (the stage for the first harvesting period) and AP_S2 (the stage for the third harvesting period). This study showed that AP_S1 had the more growth advantages than AP_S2 including biomass, auricle area and thickness, protein and calcium contents, and most species of the amino acid contents, which contributed to the higher sensory evaluation and acceptability of AP_S1. Transcriptome profiling showed that a total of 30,298 unigenes were successfully annotated in the two different harvesting periods of A. polytricha. At a threshold of two-fold change, 1,415 and 3,213 unigenes were up- and downregulated, respectively. All the differentially expressed genes (DEGs) analysis showed that the some synthesis and metabolic processes were strengthened in AP_S1, especially the synthesis and metabolism of the amino acids and protein. The enhanced energy metabolism pathways could provide more energy for AP_S1 to synthesize the nutritional substance. Moreover, the expressions of 10 selected DEGs involved in the amino acid and protein synthesis pathways and energy metabolism pathways were higher in AP_S1 compared to AP_S2, consistent with Illumina analysis. To the best of our knowledge, this is the first study that compares the nutritional properties and transcriptome profiling between the two different harvesting periods of A. polytricha and the results can present insights into the growth and genetic characteristics of A. polytricha.

Evolutionary Genomics of Sex-Related Chromosomes at the Base of the Green Lineage

Although sex is now accepted as a ubiquitous and ancestral feature of eukaryotes, direct observation of sex is still lacking in most unicellular eukaryotic lineages. Evidence of sex is frequently indirect and inferred from the identification of genes involved in meiosis from whole genome data and/or the detection of recombination signatures from genetic diversity in natural populations. In haploid unicellular eukaryotes, sex-related chromosomes are named mating-type (MTs) chromosomes and generally carry large genomic regions where recombination is suppressed. These regions have been characterized in Fungi and Chlorophyta and determine gamete compatibility and fusion. Two candidate MT+ and MT− alleles, spanning 450–650 kb, have recently been described in Ostreococcus tauri, a marine phytoplanktonic alga from the Mamiellophyceae class, an early diverging branch in the green lineage. Here, we investigate the architecture and evolution of these candidate MT+ and MT− alleles. We analyzed the phylogenetic profile and GC content of MT gene families in eight different genomes whose divergence has been previously estimated at up to 640 Myr, and found evidence that the divergence of the two MT alleles predates speciation in the Ostreococcus genus. Phylogenetic profiles of MT trans-specific polymorphisms in gametologs disclosed candidate MTs in two additional species, and possibly a third. These Mamiellales MT candidates are likely to be the oldest mating-type loci described to date, which makes them fascinating models to investigate the evolutionary mechanisms of haploid sex determination in eukaryotes.

Nuclear behavior in Coprinus comatus from Japan and estimation of the lifecycle

In order to elucidate the lifecycle of Coprinus comatus, we examined the number of nuclei in basidiospores, hyphal cells and oidia. Basidiospores isolated from the fruiting bodies of four Japanese strains were binucleate. In both primary and secondary mycelia, most of the cells were binucleate. In addition, oidia and oidiophores were observed for the first time in this mushroom and most of the oidia were binucleate. Based on these results, the lifecycle of C. comatus was inferred to be as follows. A homokaryotic binucleate basidiospore germinates and produces homokaryotic binucleate hyphae. After mating between compatible homokaryotic binucleate hyphae, a heterokaryotic binucleate secondary mycelium is produced. If environmental conditions are suitable for fruiting, homokaryotic binucleate basidiospores in the fruiting body are produced.

Intracellular pathways for lignin catabolism in white-rot fungi

Lignin is a biopolymer found in plant cell walls that accounts for 30% of the organic carbon in the biosphere. White-rot fungi (WRF) are considered the most efficient organisms at degrading lignin in nature. While lignin depolymerization by WRF has been extensively studied, the possibility that WRF are able to utilize lignin as a carbon source is still a matter of controversy. Here, we employ ¹³C-isotope labeling, systems biology approaches, and in vitro enzyme assays to demonstrate that two WRF, Trametes versicolor and Gelatoporia subvermispora, funnel carbon from lignin-derived aromatic compounds into central carbon metabolism via intracellular catabolic pathways. These results provide insights into global carbon cycling in soil ecosystems and furthermore establish a foundation for employing WRF in simultaneous lignin depolymerization and bioconversion to bioproducts-a key step toward enabling a sustainable bioeconomy.

Asexual reproduction and vegetative growth of Bionectria ochroleuca in response to temperature and photoperiod

Growth and reproduction are two essential life-history traits for fungi. Understanding life-history strategies provides insight into the environmental adaption of species. Here, we investigated the colonial morphology, vegetative growth, and asexual reproduction of the ascomycete fungus Bionectria ochroleuca in response to a variety of environmental conditions. We demonstrated that the increased temperature from 15 to 25°C induced mycelial growth and conidiation in B. ochroleuca. We also found that the optimal temperatures for mycelial growth and conidial formation in this fungus species were 25 and 30°C, respectively. However, as the temperature increased from 25 to 30°C, mycelial growth was suppressed, but the total number of conidia was significantly increased. The shift in light-dark cycles dramatically changed the morphological features of the colonies and affected both vegetative growth and asexual reproduction. Under incubation environments of alternating light and dark (16:8 and 8:16 light:dark cycles), conidiophores and conidia in the colonies formed dense-sparse rings and displayed synchronous wave structures. When the light duration was prolonged in the sequence of 0, 8, 16, and 24 hr per day, mycelial growth was suppressed, but conidiation was promoted. Together, our results indicate that temperature and light period may trigger a trade-off between vegetative growth and asexual reproduction in B. ochroleuca.

Changes in soil bacteria functional ecology associated with Morchella rufobrunnea fruiting in a natural habitat

Morchella rufobrunnea is a saprobic edible mushroom, found in a range of ecological niches, indicating nutritional adjustment to different habitats and possible interaction with soil prokaryotic microbiome (SPM). Using the 16S rRNA gene, we examined the SPM of M. rufobrunnea that appeared in a natural habitat in Northern Israel. Three sample types were included: bare soil without mushroom, soil beneath young mushroom initials and soil beneath the mature fruiting body. Morchella rufobrunnea developmental stage was significantly associated with changes in bacterial populations (PERMANOVA, p < 0.0005). Indicator analysis with point-biserial correlation coefficient found 180 operational taxonomic units (OTU) uniquely associated with distinct stages of development. The Functional Annotation of Prokaryotic Taxonomy (FAPROTAX) database helped to infer ecological roles for indicator OTU. The functional ecological progression begins with establishment of a photoautotrophic N-fixing bacterial mat on bare soil. Pioneer heterotrophs including oligotrophs, acidifying nutrient mobilizers and nitrifiers are congruent with appearance of young M. rufobrunnea initials. Under the mature fruiting body, the population changed to saprobes, organic-N degraders, denitrifiers, insect endosymbionts and fungal antagonists. Based on this work, M. rufobrunnea may be able to influence SPM and change the soil nutritional profile.

Functional analyses of Pleurotus ostreatus pcc1 and clp1 using CRISPR/Cas9

Understanding the molecular mechanisms controlling dikaryon formation in Agaricomycetes, which is basically controlled by A and B mating-type loci, contributes to improving mushroom cultivation and breeding. In Coprinopsis cinerea, various mutations in the SRY-type high mobility group protein-encoding gene, pcc1, were shown to activate the A-regulated pathway to induce pseudoclamp (clamp cells without clamp connection) and fruiting body formation in monokaryons. The formation of clamp cells was blocked in AmutBmut strain 326 with clp1-1 mutation in C. cinerea. However, considering the diverse mechanisms of sexual development among Agaricomycetes, it remains unclear whether similar phenotypes are also observed in clp1 or pcc1 mutants in cultivated mushrooms. Therefore, phenotypic analyses of Pleurotus ostreatus pcc1 or clp1 (Popcc1 or Poclp1) mutants generated using CRISPR/Cas9 were performed in this study. Plasmids with Cas9 expression cassette and different single guide RNAs targeting Popcc1 or Poclp1 were individually introduced into a monokaryotic P. ostreatus strain PC9 to obtain the mutants. Unlike in C. cinerea, the pseudoclamp cell was not observed in monokaryotic Popcc1 mutants, but it was observed after crossing two compatible strains with Popcc1 mutations. In Poclp1 mutants, dikaryosis was impaired as clamp cells were not observed after crossing, suggesting that Poclp1 functions may be essential for clamp cell formation, like in C. cinerea. These results provided a clue with respect to conserved and diverse mechanisms underlying sexual development in Agaricomycetes (at least between C. cinerea and P. ostreatus).

Transcriptome Dynamics Underlying Chlamydospore Formation in Trichoderma virens GV29-8

Trichoderma spp. are widely used biocontrol agents which are antagonistic to a variety of plant pathogens. Chlamydospores are a type of propagules produced by many fungi that have thick walls and are highly resistant to adverse environmental conditions. Chlamydospore preparations of Trichoderma spp. can withstand various storage conditions, have a longer shelf life than conidial preparations and have better application potential. However, large-scale production of chlamydospores has proven difficult. To understand the molecular mechanisms governing chlamydospore formation (CF) in Trichoderma fungi, we performed a comprehensive analysis of transcriptome dynamics during CF across 8 different developmental time points, which were divided into 4 stages according to PCA analysis: the mycelium growth stage (S1), early and middle stage of CF (S2), flourishing stage of CF (S3), and late stage of CF and mycelia initial autolysis (S4). 2864, 3206, and 3630 DEGs were screened from S2 vs S1, S3 vs S2, and S4 vs S3, respectively. We then identified the pathways and genes that play important roles in each stage of CF by GO, KEGG, STC and WGCNA analysis. The results showed that DEGs in the S2 vs S1 were mainly enriched in organonitrogen compound metabolism, those in S3 vs S2 were mainly involved in secondary metabolite, cell cycle, and N-glycan biosynthesis, and DEGs in S4 vs S3 were mainly involved in lipid, glycogen, and chitin metabolic processes. We speculated that mycelial assimilation and absorption of exogenous nitrogen in the early growth stage (S1), resulted in subsequent nitrogen deficiency (S2). At the same time, secondary metabolites and active oxygen free radicals released during mycelial growth produced an adverse growth environment. The resulting nitrogen-deficient and toxin enriched medium may stimulate cell differentiation by initiating cell cycle regulation to induce morphological transformation of mycelia into chlamydospores. High expression of genes relating to glycogen, lipid, mannan, and chitin synthetic metabolic pathways during the flourishing (S3) and late stages (S4) of CF may be conducive to energy storage and cell wall construction in chlamydospores. For further verifying the functions of the amino sugar and nucleotide sugar metabolism (tre00520) pathway in the CF of T. virens GV29-8 strain, the chitin synthase gene (TRIVIDRAFT_90152), one key gene of the pathway, was deleted and resulted in the dysplasia of mycelia and an incapability to form normal chlamydospores, which illustrated the pathway affecting the CF of T. virens GV29-8 strain. Our results provide a new perspective for understanding the genetics of biochemical pathways involved in CF of Trichoderma spp.

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.