Whole genome and global gene expression analyses of the model mushroom Flammulina velutipes reveal a high capacity for lignocellulose degradation

Transcriptomic and metabolomic insights into flavor variations in wild and cultivated Agaricus bisporus

Agaricus bisporus is a widely cultivated edible fungus globally. However, the mechanisms underlying the differences in flavor and nutritional traits between wild-type (W) and cultivated-type (C) strains remain unclear, which hinders the artificial breeding of high-quality varieties. This study systematically revealed, for the first time, the molecular and metabolic basis of flavor divergence between wild and cultivated A. bisporus by integrating transcriptomics and metabolomics. A total of 43 strains (23 wild-type and 20 cultivated-type) were analyzed using high-throughput sequencing and liquid chromatography-tandem mass spectrometry (LC-MS/MS) to dissect differences in gene expression and metabolite profiles. Results showed that although total protein and amino acid contents exhibited no significant differences, transcriptomic analysis identified significant upregulation of AGABI2DRAFT_188981 and AGABI2DRAFT_191000 (genes associated with high-affinity methionine permease MUP1) in cultivated strains, suggesting their indirect regulation of flavor formation via methionine metabolism. Metabolomic analysis further revealed a marked increase in uridine levels in cultivated strains (3.2-fold higher than wild-type, p < 0.01), indicating potential medicinal value, while wild strains were enriched with flavor precursors such as fumaric acid and isoleucine (fold change ≥ 2.5). In contrast, cultivated strains accumulated metabolites like 2-hydroxybutyric acid and α-ketoglutarate (VIP > 1.5). This study pioneered the construction of a gene-metabolite correlation network, identifying a strong positive correlation between AGABI2DRAFT_191352 (6-phosphofructokinase) and 2-hydroxybutyric acid (r = 0.82), highlighting the regulatory role of glycolytic flux in flavor metabolism. These findings not only elucidate the impact of artificial cultivation on metabolic reprogramming in A. bisporus but also provide critical molecular targets for targeted breeding of strains with enhanced flavor and nutritional value, offering practical significance for advancing the edible fungi industry.

Whole-Genome Sequencing of Flammulina filiformis and Multi-Omics Analysis in Response to Low Temperature

The growth of Flammulina filiformis is strongly dependent on low-temperature cues for the initiation of primordia formation. To obtain a comprehensive understanding of the molecular mechanisms that govern the mycelial response to cold stress, de novo genome sequencing of the F. filiformis monokaryon and multi-omics data (transcriptome and metabolome) analyses of the mycelia, primordia, and fruiting bodies were conducted in the present study. Genome sequencing based on PacBio HiFi and Hi-C resulted in a 36.3 Mb genome sequence that mapped to 12 chromosomes, comprising 11,886 protein-coding genes. A total of 25 cold-responsive (COR) genes and 520 cold-adapted enzymes were identified in the genome. Multi-omics analyses showed that the pathways related to carbohydrate metabolism in the mycelia under low temperature (10 °C) were significantly enriched. Further examination of the expression profiles of carbohydrate-active enzymes (CAZymes) involved in carbohydrate metabolism revealed that out of 515 CAZyme genes in F. filiformis, 58 were specifically upregulated in mycelia under low-temperature conditions. By contrast, the expression levels of these genes in primordia and fruiting bodies reverted to those prior to low-temperature exposure. These indicate that CAZyme genes are important for the low-temperature adaptation of F. filiformis. This research contributes to the targeted breeding of F. filiformis.

Transcriptional analysis provides molecular insights into degeneration of the edible fungus Flammulina filiformis

AIMS

Flammulina filiformis, a widely cultivated edible fungus, frequently suffers from strain degeneration, leading to reduce yield and quality, and causing significant economic losses in large-scale production. This study aimed to better understand degeneration mechanisms and to develop an assay for predicting degeneration.

METHODS AND RESULTS

This study investigates strain degeneration by comparing a cultivated strain (F1) and its subcultured derivative (F7). Although both strains display similar mycelial growth, F7 shows impaired fruiting body production, indicating that degeneration occurs prior to visible changes in mycelial growth. To uncover the molecular mechanisms underlying degeneration, transcriptome analysis was carried out. A total of 352 down-regulated and 280 up-regulated differentially expressed genes (DEGs) were identified in F7 compared to F1. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses revealed significant enrichment of DEGs in functional categories related to oxidoreductase activities and peroxisome pathway. Quantitative PCR further validated the regulation of certain DEGs associated with these enriched functions. A stress tolerance assay was developed to detect degeneration in strains with unchanged mycelial growth but reduced fruiting body production.

CONCLUSIONS

Results suggested that strain degeneration in F. filiformis is closely linked to oxidative stress regulation and occurs prior to observable mycelial impairment.

Glutathione peroxidase gene regulates substrate development and prevents strain aging in Volvariella volvacea

Degradation of Volvariella volvacea is associated with the accumulation of reactive oxygen species (ROS), and glutathione peroxidase (GPX) is one of the key antioxidant enzyme. The purpose of this research is to uncover the importance of the gpx gene in the degradation and revitalization of V. volvacea. In this study, a gpx-silenced strain of the primordial strain T0 and an overexpression strain of the degraded strain T19 of V. volvacea were constructed, and their antioxidant properties, matrix degradation ability, and production traits were determined. The results showed that the expression level of gpx altered the homeostasis of the V. volvacea redox system and affected the substrate degradation ability of V. volvacea, which altered the physiological traits of the V. volvacea mycelium. Most importantly, the primordial strain T0 was unable to produce fruiting bodies due to the silence of the gpx gene. On the other hand, the gpx gene overexpression promoted the regrowth of fruiting bodies in degenerated strains of V. volvacea T19. This study provides a new biotechnological strategy to control the degeneration of V. volvacea and other edible fungi.

Chromosome-Level Genome Assembly and Annotation of the Highly Heterozygous Phallus echinovolvatus Provide New Insights into Its Genetics

Phallus echinovolvatus is a well-known edible and medicinal fungus with significant economic value. However, the available whole-genome information is lacking for this species. The chromosome-scale reference genome (Monop) and two haploid genomes (Hap1 and Hap2) of P. echinovolvatus, each assembled into 11 pseudochromosomes, were constructed using Illumina, PacBio-HiFi long-read sequencing, and Hi-C technology. The Monop had a size of 36.54 Mb, with 10,251 predicted protein-coding genes and including 433 carbohydrate-active enzyme genes, 385 cytochrome P450 enzyme genes, and 42 gene clusters related to secondary metabolite synthesis. Phylogenetic and collinearity analysis revealed a close evolutionary relationship between P. echinovolvatus and Clathrus columnatus in the core Phallales clade. Hap1 and Hap2 had sizes of 35.46 Mb and 36.11 Mb, respectively. Collinear relationships were not observed for 15.38% of the genes in the two haplotypes. Hap1 had 256 unique genes, and Hap2 had 370 unique genes. Our analysis of the P. echinovolvatus genome provides insights into the genetic basis of the mechanisms underlying the metabolic effects of bioactive substances and will aid ongoing breeding efforts and studies of genetic mechanisms.

Whole-Genome Sequence Analysis of Flammulina filiformis and Functional Validation of Gad, a Key Gene for γ-Aminobutyric Acid Synthesis

Flammulina filiformis is one of the widely produced edible fungi worldwide. It is rich in γ-aminobutyric acid (GABA), a non-protein amino acid with important physiological functions in humans. To investigate the functions of key genes in the GABA metabolic pathway of F. filiformis, we isolated the monokaryon Fv-HL23-1 from the factory-cultivated F. filiformis strain Fv-HL23 and then sequenced and assembled the genome using the PacBio Sequel and Illumina NovaSeq sequencing platforms. The results showed that the genome comprised 140 scaffolds with a total length of 40.96 Mb, a GC content of 49.62%, an N50 of 917,125 bp, and 14,256 protein-coding genes. Phylogenetic analysis based on the whole genome revealed a close evolutionary relationship of Fv-HL23-1 with Armillaria mellea, Lentinula edodes, and Schizophyllum commune. A total of 589 carbohydrate-active enzymes were identified in the genome of Fv-HL23-1, suggesting its strong lignocellulose degradation ability, and 108 CYP450 gene family members were identified, suggesting important functions such as resistance to stress, secondary metabolite synthesis, and growth and development. The F. filiformis proteins glutamate decarboxylase 1 (Ff-GAD1) and glutamate decarboxylase 2 (Ff-GAD2), which may be responsible for GABA synthesis, were identified by protein alignment. Molecular docking analysis showed that Ff-GAD2 may have better catalytic activity than Ff-GAD1. To verify the function of Ff-gad2, its heterologous expression in the mycelia of the mononuclear Hypsizigus marmoreus was analyzed. Compared with wild type, the GABA content of mycelia was increased by 85.40-283.90%, the growth rate was increased by 9.39 ± 2.35%, and the fresh weight was increased by 18.44 ± 7.57%. Ff-GAD2 may play a catalytic role in GABA synthesis. In addition, the expression of the full-length Ff-gad2 gene was increased by 7.96 ± 1.39 times compared with the exon expression level in H. marmoreus mycelia, suggesting that the intron may contribute to the heterologous expression of Ff-GAD2. Based on whole-genome sequencing, we analyzed the enzyme system related to the important life activities of F. filiformis, focusing on the function of Ff-GAD, a key enzyme in the GABA synthesis pathway. The results lay a foundation for elucidating the GABA metabolism pathway of edible fungi and developing targeted breeding strategies for GABA-producing edible fungi.

Pan-genome analysis reveals genomic variations during enoki mushroom domestication, with emphasis on genetic signatures of cap color and stipe length

INTRODUCTION

The domestication of edible mushrooms, including Flammulina filiformis, offers valuable insights into the genetic changes driven by artificial selection. Understanding these changes is crucial for uncovering the mechanisms behind genome evolution in domesticated mushrooms.

OBJECTIVES

This study aims to investigate the population structure, genetic diversity, and domestication-related genomic changes in F. filiformis. By comparing the genome sequences of 199 wild and cultivated strains, we aim to elucidate the impact of domestication on F. filiformis.

METHODS

We performed de novo genome assembly and gene-based pan-genome analysis on the 199 strains, which included both wild and cultivated strains. We also conducted genome-wide association studies (GWAS) using presence-absence variation (PAV) and SNP data, combined with RNA sequencing, to identify genes associated with domestication traits, such as cap color and stipe length. Gene functional confirmation was achieved through genetic transformation experiments.

RESULTS

Our analysis grouped the strains into four distinct populations, which correlated with varying intensities of artificial selection. The three cultivated populations exhibited smaller genome sizes, fewer genes, lower genetic variation, reduced gene expression diversity, and lower heterozygosity compared to the wild population. The analysis revealed the loss of genes related to the beta-lactam antibiotic catabolic process and specific MAPK pathway genes during domestication, rendering domesticated strains more susceptible to diseases. Four genes closely associated with cap color and stipe length were identified, but genetic transformation experiments confirmed the functional relevance of only two (FfB and FfD) identified through PAV-based GWAS.

CONCLUSION

This study uncovered significant genomic variations between cultivated and wild F. filiformis populations, including the loss of pathogen resistance genes during domestication. We also identified key genes linked to cap color and stipe length, demonstrating for the first time the important role of PAV variation in mushroom domestication. These insights provide a foundation for future mushroom breeding and evolutionary research.

Genome-Wide Identification and Expression Analysis of the Cys2His2 Zinc Finger Protein Gene Family in Flammulina filiformis

Zinc finger proteins (ZFPs) are essential transcription factors in eukaryotes, particularly the extensively studied C2H2 family, which is known for its involvement in various biological processes. This research provides a thorough examination and analysis of the C2H2-ZFP gene family in Flammulina filiformis. Using bioinformatics tools, 58 FfC2H2-ZFP genes spread across 11 chromosomes were identified and scrutinized in detail for their gene structures, protein characteristics, and phylogenetic relationships. The study of phylogenetics and synteny sheds light on the evolutionary relationships among C2H2-ZFPs in F. filiformis and other fungi, revealing a complex evolutionary past. The identification of conserved cis-regulatory elements in the gene promoter regions suggests intricate functionalities, particularly in the developmental and stress response pathways. By utilizing RNA-seq and qRT-PCR techniques, the expression patterns of these genes were explored across different developmental stages and tissues of F. filiformis, unveiling distinct expression profiles. Notably, significant expression variations were observed in the stipe elongation region and pilei of various sizes, indicating potential roles in fruiting body morphogenesis. This study enhances our knowledge of the C2H2-ZFP gene family in F. filiformis and lays the groundwork for future investigations into their regulatory mechanisms and applications in fungal biology and biotechnology.

Serine Rejuvenated Degenerated Volvariella volvacea by Enhancing ROS Scavenging Ability and Mitochondrial Function

Serine is a functional amino acid that effectively regulates the physiological functions of an organism. This study investigates the effects of adding exogenous serine to a culture medium to explore a feasible method for the rejuvenation of V. volvacea degenerated strains. The tissue isolation subcultured strains T6, T12, and T19 of V. volvacea were used as test strains, and the commercially cultivated strain V844 (T0) was used as a control. The results revealed that the addition of serine had no significant effect on non-degenerated strains T0 and T6, but could effectively restore the production characteristics of degenerated strains T12 and T19. Serine increased the biological efficiency of T12 and even helped the severely degenerated T19 to regrow its fruiting body. Moreover, exogenous serine up-regulated the expression of some antioxidant enzyme genes, improved antioxidase activity, reduced the accumulation of reactive oxygen species (ROS), lowered malondialdehyde (MDA) content, and restored mitochondrial membrane potential (MMP) and mitochondrial morphology. Meanwhile, serine treatment increased lignocellulase and mycelial energy levels. These findings form a theoretical basis and technical support for the rejuvenation of V. volvacea degenerated strains and other edible fungi.

What are the 100 most cited fungal genera?

The global diversity of fungi has been estimated between 2 to 11 million species, of which only about 155 000 have been named. Most fungi are invisible to the unaided eye, but they represent a major component of biodiversity on our planet, and play essential ecological roles, supporting life as we know it. Although approximately 20 000 fungal genera are presently recognised, the ecology of most remains undetermined. Despite all this diversity, the mycological community actively researches some fungal genera more commonly than others. This poses an interesting question: why have some fungal genera impacted mycology and related fields more than others? To address this issue, we conducted a bibliometric analysis to identify the top 100 most cited fungal genera. A thorough database search of the Web of Science, Google Scholar, and PubMed was performed to establish which genera are most cited. The most cited 10 genera are Saccharomyces, Candida, Aspergillus, Fusarium, Penicillium, Trichoderma, Botrytis, Pichia, Cryptococcus and Alternaria. Case studies are presented for the 100 most cited genera with general background, notes on their ecology and economic significance and important research advances. This paper provides a historic overview of scientific research of these genera and the prospect for further research. Citation: Bhunjun CS, Chen YJ, Phukhamsakda C, Boekhout T, Groenewald JZ, McKenzie EHC, Francisco EC, Frisvad JC, Groenewald M, Hurdeal VG, Luangsa-ard J, Perrone G, Visagie CM, Bai FY, Błaszkowski J, Braun U, de Souza FA, de Queiroz MB, Dutta AK, Gonkhom D, Goto BT, Guarnaccia V, Hagen F, Houbraken J, Lachance MA, Li JJ, Luo KY, Magurno F, Mongkolsamrit S, Robert V, Roy N, Tibpromma S, Wanasinghe DN, Wang DQ, Wei DP, Zhao CL, Aiphuk W, Ajayi-Oyetunde O, Arantes TD, Araujo JC, Begerow D, Bakhshi M, Barbosa RN, Behrens FH, Bensch K, Bezerra JDP, Bilański P, Bradley CA, Bubner B, Burgess TI, Buyck B, Čadež N, Cai L, Calaça FJS, Campbell LJ, Chaverri P, Chen YY, Chethana KWT, Coetzee B, Costa MM, Chen Q, Custódio FA, Dai YC, Damm U, de Azevedo Santiago ALCM, De Miccolis Angelini RM, Dijksterhuis J, Dissanayake AJ, Doilom M, Dong W, Alvarez-Duarte E, Fischer M, Gajanayake AJ, Gené J, Gomdola D, Gomes AAM, Hausner G, He MQ, Hou L, Iturrieta-González I, Jami F, Jankowiak R, Jayawardena RS, Kandemir H, Kiss L, Kobmoo N, Kowalski T, Landi L, Lin CG, Liu JK, Liu XB, Loizides M, Luangharn T, Maharachchikumbura SSN, Makhathini Mkhwanazi GJ, Manawasinghe IS, Marin-Felix Y, McTaggart AR, Moreau PA, Morozova OV, Mostert L, Osiewacz HD, Pem D, Phookamsak R, Pollastro S, Pordel A, Poyntner C, Phillips AJL, Phonemany M, Promputtha I, Rathnayaka AR, Rodrigues AM, Romanazzi G, Rothmann L, Salgado-Salazar C, Sandoval-Denis M, Saupe SJ, Scholler M, Scott P, Shivas RG, Silar P, Souza-Motta CM, Silva-Filho AGS, Spies CFJ, Stchigel AM, Sterflinger K, Summerbell RC, Svetasheva TY, Takamatsu S, Theelen B, Theodoro RC, Thines M, Thongklang N, Torres R, Turchetti B, van den Brule T, Wang XW, Wartchow F, Welti S, Wijesinghe SN, Wu F, Xu R, Yang ZL, Yilmaz N, Yurkov A, Zhao L, Zhao RL, Zhou N, Hyde KD, Crous PW (2024). What are the 100 most cited fungal genera? Studies in Mycology 108: 1-411. doi: 10.3114/sim.2024.108.01.

Functional Characterization of Sesquiterpene Synthases and P450 Enzymes in Flammulina velutipes for Biosynthesis of Spiro [4.5] Decane Terpene

Flammulina velutipes, a popular edible mushroom, contains sesquiterpenes with potential health benefits. We characterized 12 sesquiterpene synthases and one P450 enzyme in F. velutipes using Aspergillus oryzae as a heterologous expression system, culminating in the biosynthesis of 16 distinct sesquiterpene compounds. An enzyme encoded by the axeB gene responsible for the synthesis of the spiro [4.5] decane compound axenol was discovered, and the mechanism of spirocycle formation was elucidated through quantum mechanical calculations. Furthermore, we delineated the role of a P450 enzyme colocated with AxeB in producing the novel compound 3-oxo-axenol. Our findings highlight the diverse array of sesquiterpene skeletons and functional groups biosynthesized by these enzymes in F. velutipes and underscore the effectiveness of the A. oryzae system as a heterologous host for expressing genes in the Basidiomycota genome. These insights into the biosynthesis of bioactive compounds in F. velutipes have significant implications for functional food and drug development.

Decoding the chromosome-scale genome of the nutrient-rich Agaricus subrufescens: a resource for fungal biology and biotechnology

Agaricus subrufescens, also known as the "sun mushroom," has significant nutritional and medicinal value. However, its short shelf life due to the browning process results in post-harvest losses unless it's quickly dehydrated. This restricts its availability to consumers in the form of capsules. A genome sequence of A. subrufescens may lead to new cultivation alternatives or the application of gene editing strategies to delay the browning process. We assembled a chromosome-scale genome using a hybrid approach combining Illumina and Nanopore sequencing. The genome was assembled into 13 chromosomes and 31 unplaced scaffolds, totaling 44.5 Mb with 96.5% completeness and 47.24% GC content. 14,332 protein-coding genes were identified, with 64.6% of the genome covered by genes and 23.41% transposable elements. The mitogenome was circularized and encoded fourteen typical mitochondrial genes. Four polyphenol oxidase (PPO) genes and the Mating-type locus were identified. Phylogenomic analysis supports the placement of A. subrufescens in the Agaricomycetes clade. This is the first available genome sequence of a strain of the "sun mushroom." Results are available through a Genome Browser (https://plantgenomics.ncc.unesp.br/gen.php?id=Asub) and can support further fungal biological and genomic studies.

Chromosome-level assembly of Dictyophora rubrovolvata genome using third-generation DNA sequencing and Hi-C analysis

Dictyophora rubrovolvata, a rare edible mushroom with both nutritional and medicinal values, was regarded as the "queen of the mushroom" for its attractive appearance. Dictyophora rubrovolvata has been widely cultivated in China in recent years, and many researchers were focusing on its nutrition, culture condition, and artificial cultivation. Due to a lack of genomic information, research on bioactive substances, cross breeding, lignocellulose degradation, and molecular biology is limited. In this study, we report a chromosome-level reference genome of D. rubrovolvata using the PacBio single-molecule real-time-sequencing technique and high-throughput chromosome conformation capture (Hi-C) technologies. A total of 1.83 Gb circular consensus sequencing reads representing ∼983.34 coverage of the D. rubrovolvata genome were generated. The final genome was assembled into 136 contigs with a total length of 32.89 Mb. The scaffold and contig N50 length were 2.71 and 2.48 Mb, respectively. After chromosome-level scaffolding, 11 chromosomes with a total length of 28.24 Mb were constructed. Genome annotation further revealed that 9.86% of the genome was composed of repetitive sequences, and a total of 508 noncoding RNA (rRNA: 329, tRNA: 150, ncRNA: 29) were annotated. In addition, 9,725 protein-coding genes were predicted, among which 8,830 (90.79%) genes were predicted using homology or RNA-seq. Benchmarking Universal Single-Copy Orthologs results further revealed that there were 80.34% complete single-copy fungal orthologs. In this study, a total of 360 genes were annotated as belonging to the carbohydrate-active enzymes family. Further analysis also predicted 425 cytochromes P450 genes, which can be classified into 41 families. This highly accurate, chromosome-level reference genome of D. rubrovolvata will provide essential genomic information for understanding the molecular mechanism in its fruiting body formation during morphological development and facilitate the exploitation of medicinal compounds produced by this mushroom.

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.

High-quality genome assembly and multi-omics analysis of pigment synthesis pathway in Auricularia cornea

Owing to its great market potential for food and health care, white Auricularia cornea, a rare edible fungus, has received increased attention in recent years. This study presents a high-quality genome assembly of A. cornea and multi-omics analysis of its pigment synthesis pathway. Continuous Long Reads libraries, combined with Hi-C-assisted assembly were used to assemble of white A. cornea. Based on this data, we analyzed the transcriptome and metabolome of purple and white strains during the mycelium, primordium, and fruiting body stages. Finally, we obtained the genome of A.cornea assembled from 13 clusters. Comparative and evolutionary analysis suggests that A.cornea is more closely related to Auricularia subglabra than to Auricularia heimuer. The divergence of white/purple A.cornea occurred approximately 40,000 years ago, and there were numerous inversions and translocations between homologous regions of the two genomes. Purple strain synthesized pigment via the shikimate pathway. The pigment in the fruiting body of A. cornea was γ-glutaminyl-3,4-dihydroxy-benzoate. During pigment synthesis, α-D-glucose-1P, citrate, 2-Oxoglutarate, and glutamate were four important intermediate metabolites, whereas polyphenol oxidase and other 20 enzyme genes were the key enzymes. This study sheds light on the genetic blueprint and evolutionary history of the white A.cornea genome, revealing the mechanism of pigment synthesis in A.cornea. It has important theoretical and practical implications for understanding the evolution of basidiomycetes, molecular breeding of white A.cornea, and deciphering the genetic regulations of edible fungi. Additionally, it provides valuable insights for the study of phenotypic traits in other edible fungi.

Functional Clustering of Metabolically Related Genes Is Conserved across Dikarya

Transcriptional regulation is vital for organismal survival, with many layers and mechanisms collaborating to balance gene expression. One layer of this regulation is genome organization, specifically the clustering of functionally related, co-expressed genes along the chromosomes. Spatial organization allows for position effects to stabilize RNA expression and balance transcription, which can be advantageous for a number of reasons, including reductions in stochastic influences between the gene products. The organization of co-regulated gene families into functional clusters occurs extensively in Ascomycota fungi. However, this is less characterized within the related Basidiomycota fungi despite the many uses and applications for the species within this clade. This review will provide insight into the prevalence, purpose, and significance of the clustering of functionally related genes across Dikarya, including foundational studies from Ascomycetes and the current state of our understanding throughout representative Basidiomycete species.

Phenylalanine Ammonia-Lyase: A Key Gene for Color Discrimination of Edible Mushroom Flammulina velutipes

In nature; Flammulina velutipes, also known as winter mushrooms, vary in the color of their fruiting bodies, from black, yellow, pale yellow, or beige to white. The purpose of this study was to compare the genome sequences of different colored strains of F. velutipes and to identify variations in the genes associated with fruiting body color. Comparative genomics of six F. velutipes strains revealed 70 white-strain-specific variations, including single nucleotide polymorphisms (SNPs) and insertions/deletions (indels), in the genome sequences. Among them, 36 variations were located in the open reading frames, and only one variation was identified as a mutation with a disruptive in-frame deletion (ΔGCGCAC) within the annotated gene phenylalanine ammonia-lyase 1 (Fvpal1). This mutation was found to cause a deletion, without a frameshift, of two amino acids at positions 112 and 113 (arginine and threonine, respectively) in the Fvpal1 gene of the white strain. Specific primers to detect this mutation were designed, and amplification refractory mutation system (ARMS) polymerase chain reaction (PCR) was performed to evaluate whether the mutation is color specific for the F. velutipes fruiting body. PCR analysis of a total of 95 F. velutipes strains revealed that this mutation was present only in white strains. In addition, monospores of the heterozygous mutant were isolated, and whether this mutation was related to the color of the fruiting body was evaluated by a mating assay. In the mating analysis of monospores with mutations in Fvpal1, it was found that this mutation plays an important role in determining the color of the fruiting body. Furthermore, the deletion (Δ112RT113) in Fvpal1 is located between motifs that play a key role in the catalytic function of FvPAL1. These results suggest that this mutation can be used as an effective marker for the color-specific breeding of F. velutipes, a representative edible mushroom.

Development of Multiple Nucleotide Polymorphism Molecular Markers for Enoki Mushroom (Flammulina filiformis) Cultivars Identification

The enoki mushroom (Flammulina filiformis) is one of the most important and popular edible mushrooms commercially in China. However, traditional mushroom cultivar identification is challenging due to poor accuracy, heavy workloads, and low reproducibility. To overcome this challenge, we developed a method for identifying F. filiformis strains using multiple nucleotide polymorphism sequencing (MNP-seq). This involved screening 179 universal MNP markers based on whole-genome sequencing data, constructing an MNP sequence library, and performing multiplex PCR amplification and high-sequencing. We further screened 69 core MNP markers and used them to build a neighbor-joining (NJ) phylogenetic tree of 232 cultivated and wild strains. Our analysis showed that all cultivars could be accurately separated by computing genetic similarity values and that the cultivars could be separated into 22 distinct evolutionary pedigrees. The specific value of genetic similarity can be used as the standard to distinguish F. filiformis cultivars, however, it needs to be comprehensively defined by the additional phenotype and biological characteristics of those strains in the future work.

Increasing the production of the bioactive compounds in medicinal mushrooms: an omics perspective

Macroscopic fungi, mainly higher basidiomycetes and some ascomycetes, are considered medicinal mushrooms and have long been used in different areas due to their pharmaceutically/nutritionally valuable bioactive compounds. However, the low production of these bioactive metabolites considerably limits the utilization of medicinal mushrooms both in commerce and clinical trials. As a result, many attempts, ranging from conventional methods to novel approaches, have been made to improve their production. The novel strategies include conducting omics investigations, constructing genome-scale metabolic models, and metabolic engineering. So far, genomics and the combined use of different omics studies are the most utilized omics analyses in medicinal mushroom research (both with 31% contribution), while metabolomics (with 4% contribution) is the least. This article is the first attempt for reviewing omics investigations in medicinal mushrooms with the ultimate aim of bioactive compound overproduction. In this regard, the role of these studies and systems biology in elucidating biosynthetic pathways of bioactive compounds and their contribution to metabolic engineering will be highlighted. Also, limitations of omics investigations and strategies for overcoming them will be provided in order to facilitate the overproduction of valuable bioactive metabolites in these valuable organisms.

Heterogeneous Distribution of Phospholipid Molecular Species in the Surface Culture of Flammulina velutipes: New Facts about Lipids Containing α-Linolenic Fatty Acid

Mycelial fungi grow as colonies consisting of polar growing hyphae, developing radially from spore or inoculum. Over time, the colony develops, hyphae are subject to various exogenous or endogenous stimuli, and mycelium becomes heterogeneous in growth, gene expression, biosynthesis, and secretion of proteins and metabolites. Although the biochemical and molecular mechanisms of mycelium heterogeneity have been the subject of many studies, the role of lipids in colony development and zonality is still not understood. This work was undertaken to extend our knowledge of mycelium heterogeneity and to answer the question of how different lipid molecular species are distributed in the surface colony of the basidial fungus Flammulina velutipes and how this distribution correlates with its morphology. The heterogeneity in the lipid metabolism and lipid composition of the fungal mycelium was demonstrated. According to the real-time PCR and LC-MS/MS results, the expression of genes of PC metabolism, accumulation of phospholipid classes, and degree of unsaturation of PC and PE increased in the direction from the center to the periphery of the colony. The peripheral zone of the colony was characterized by a higher value of the PC/PE ratio and a higher level of phospholipids esterified by linolenic acid. Considering that the synthesis of phospholipids in fungi occurs in different ways, we also conducted experiments with deuterium-labeled phospholipid precursors and found out that the Kennedy pathway is the predominant route for PC biosynthesis in F. velutipes. The zonal differences in gene expression and lipid composition can be explained by the participation of membrane lipids in polar growth maintenance and regulation.

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

Introduction

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

Methods

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

Results and Discussion

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

Morphological and transcriptomic analysis revealing morphological variations and genetic clues in one Lentinula edodes abnormal browning strain

Strain abnormal browning is a common problem during cultivation of Lentinula edodes. In this study, the L. edodes strain mycelia isolated from Le-WB and cultured on MYG (Le-WP) isolated from an abnormal browning bag was compared with its normal control mycelia isolated from Le-BB and cultured on MYG (Le-BP). The aerial hyphae of Le-WP were white, and the hyphal growth was significantly reduced. Morphological observation of Le-WP under scanning electron microscope (SEM) and transmission electron microscopy (TEM) revealed abnormal organelle structures. Through transcriptomic analysis, more differentially expressed genes (DEGs) were expressed in the metabolic process and catalytic activity in Le-WP than Le-BP. Two Kyoto encyclopedia of genes and genomes (KEGG) pathways named pentose and glucorunate interconversions, and starch and sucrose metabolism were found to be enriched in Le-WP. The gene expression profiles involved in these two pathways were further analyzed and 12 key genes were selected to be verified by quantitative real-time PCR (qRT-PCR), and the results showed that most of these genes were upregulated in Le-WP. Additionally, the content of 1,3-beta-glucan in Le-WP was also significantly higher than in other samples. This research suggests that abnormal strains may be related to the abnormal synthesis of 1,3-beta-glucan, and it needs further research. This research exhibits possible morphological and genetic clues of Le-WP and lays the foundation for understanding the degeneration of L. edodes strains.

Expanding the valorization of waste mushroom substrates in agricultural production: progress and challenges

Waste mushroom substrate (WMS) generated in large quantities from mushroom production process has caused severe environmental pollution. As a sustainable resource, the valorization of WMS in the agricultural field has attracted attention due to the abundant active components. A comprehensive review of valorization of WMS in agricultural production is meaningful to promote the further utilization of this resource. This paper provided an overview of the valorization in sustainable agricultural production using WMS, including animal and crop farming improvement, and agricultural environmental restoration. Moreover, the limitations and the possible development directions of WMS in agricultural production were discussed. Different sustainable cycle models for WMS in agricultural production were proposed. The aim of this review is to provide a feasible solution for the favorable treatment of WMS and improvement of agricultural production quality.

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.

Transcription factor FfMYB15 regulates the expression of cellulase gene FfCEL6B during mycelial growth of Flammulina filiformis

Background

Cellulose degradation can determine mycelial growth rate and affect yield during the growth of Flammulina filiformis. The degradation of cellulose requires the joint action of a variety of cellulases, and some cellulase-related genes have been detected in mushrooms. However, little is known about the transcriptional regulatory mechanisms of cellulose degradation.

Results

In this study, FfMYB15 that may regulate the expression of cellulase gene FfCEL6B in F. filiformis was identified. RNA interference (RNAi) showed that FfCEL6B positively regulated mycelial growth. Gene expression analyses indicated that the expression patterns of FfCEL6B and FfMYB15 in mycelia cultured on the 0.9% cellulose medium for different times were similar with a correlation coefficient of 0.953. Subcellular localization and transcriptional activity analyses implied that FfMYB15 was located in the nucleus and was a transcriptional activator. Electrophoretic mobility shift assay (EMSA) and dual-luciferase assays demonstrated that FfMYB15 could bind and activate FfCEL6B promoter by recognizing MYB cis-acting element.

Conclusions

This study indicated that FfCEL6B played an active role in mycelial growth of F. filiformis and was regulated by FfMYB15.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-022-01932-z.

Comparative analyses of Flammulina filiformis mitochondrial genomes reveal high length polymorphism in intergenic regions and multiple intron gain/loss in cox1

The golden-needle mushroom Flammulina filiformis is one of the bulk mushroom products in the world. This study obtained complete mitogenomes of 44 wild isolates collected from nine provinces and two artificially bred cultivars of F. filiformis, together with three Flammulina rossica isolates and one Flammulina fennae isolate for comparison. The mitogenome of F. filiformis ranged from 83,540 bp to 90,938 bp, consisting of 14 conserved protein-coding genes (PCGs), two rRNA genes, and 25 tRNA genes. To the best of our knowledge, it contained the highest proportion of intergenic regions compared to the other known Basidiomycota mitogenomes. Introns and intergenic regions were two major contributing factors to the total size of the F. filiformis mitogenome. The conserved PCG cox3 is located in an intron of another conserved PCG, nad5. This is a unique phenomenon in all known fungal mitogenomes. Gain/loss of introns was observed in cox1, nad5, and rnl. Length polymorphism was widely observed in intergenic regions. Accordingly, primers were designed as useful markers for rapid identification of F. filiformis isolates with differentiated mitogenomes. Our findings provide a basis for further studies related to variety identification and population genetics of this economically important mushroom.

Successive mycelial subculturing decreased lignocellulase activity and increased ROS accumulation in Volvariella volvacea

Strain degradation is a common problem in many artificially-cultivated edible mushrooms. As a fungus with poor tolerance to low-temperature, Volvariella volvacea cannot delay its degradation by long-term low temperature storage like other fungi, so its degradation is particularly severe, which hinders industrial applications. Periodic mycelial subculture is a common storage method for V. volvacea, but excessive subculturing can also lead to strain degeneration. After 20 months of continuous subculturing every 3 days, V. volvacea strains S1–S20 were obtained, and their characteristics throughout the subculture process were analyzed. With increasing number of subculture, the growth rate, mycelial biomass, the number of fruiting bodies and biological efficiency gradually decreased while the production cycle and the time to primordium formation was lengthened. Strains S13–S20, obtained after 13–20 months of mycelial subculturing, also lacked the ability to produce fruiting bodies during cultivation experiments. Determination of reactive oxygen species (ROS) content as well as enzyme activity showed that decreased lignocellulase activity, along with excessive accumulation of ROS, was concomitant with the subculture-associated degeneration of V. volvacea. Reverse transcription polymerase chain reaction (RT-PCR) was eventually used to analyze the gene expression for lignocellulase and antioxidant enzymes in subcultured V. volvacea strains, with the results found to be consistent with prior observations regarding enzyme activities. These findings could form the basis of further studies on the degeneration mechanism of V. volvacea and other fungi.

Whole-Genome Sequencing and Comparative Genomics Analysis of the Wild Edible Mushroom (Gomphus purpuraceus) Provide Insights into Its Potential Food Application and Artificial Domestication

Gomphus purpuraceus (Iwade) Yokoyama is a species of wild fungi that grows in southwest China, considered an edible and medicinal fungus with potential commercial prospects. However, the detailed mechanisms related to the development of mycelium and the formation of the fruiting body are unclear. To obtain a comprehensive overview of genetic features, whole-genome and comparative genomics analyses of G. purpuraceus were performed. High-quality DNA was extracted from the mycelium, which was isolated from a fresh fruiting body of G. purpuraceus. The DNA sample was subjected to sequencing using Illumina and Oxford Nanopore sequencing platforms. A genome assembly totaling 40.15 Mb in 50 contigs with an N50 length of 2.06 Mb was generated, and 8705 putative predicted genes were found. Subsequently, phylogenetic analysis revealed a close evolutionary relationship between G. purpuraceus and Gomphus bonarii. Moreover, a total of 403 carbohydrate-active enzymes (CAZymes) were identified in G. purpuraceus, which included 147 glycoside hydrolases (GHs), 85 glycosyl transferases (GTs), 8 polysaccharide lyases (PLs), 76 carbohydrate esterases (CEs), 57 auxiliary activities (AAs) and 30 carbohydrate-binding modules (CBMs). Compared with the other 13 fungi (Laccaria bicolor, Russula virescens, Boletus edulis, etc.), the number and distribution of CAZymes in G. purpuraceus were similar to other mycorrhizal fungi. Furthermore, the optimization of culture medium for G. purpuraceus showed the efficient utilization of disaccharides such as sucrose and maltose. The genome of G. purpuraceus provides new insights into its niche, food applications and potential artificial domestication.

Comprehensive Genetic Analysis of Monokaryon and Dikaryon Populations Provides Insight Into Cross-Breeding of Flammulina filiformis

Flammulina filiformis, as one of the most popular edible fungi in East Asia, is produced in an industrialized and standardized way. However, its monotonous variety and product convergence have seriously restricted the development of the industry. In this study, 11 cultivated strains and 13 wild strains of F. filiformis were collected from multiple regions of China and Japan and were performed genome sequencing. Together with genome data of six strains previously released, in total 23 dikaryons (formed by two monokaryons mating, can making fruiting body), 35 monokaryons (formed by protoplast-regenerating of dikaryon and isolating) were used for genetic diversity and population structure analysis based on the high-throughput genotyping. Firstly, a set of SNP markers with intrapopulation polymorphism including 849,987 bi-allelic SNPs were developed and basically covered all of 11 chromosomes with a high distribution density of 24.16 SNP markers per kb. The cultivated dikaryotic strains were divided into three subgroups, and their breeding history was made inferences, which is consistent with the available pedigree records. The wild dikaryotic strains were divided into two subgroups and showed varied contributions of genetic components with high genetic diversity. All the investigated dikaryons have a symmetric distribution pattern with their two constituent monokaryons in principal component analysis. Finally, we summarized the pedigree relationship diagram of F. filiformis main strains including six modules, and the genotypes of hybrids can be directly phased by the known parental allele according to it. This study provides a method to distinguish two sets of monokaryon haplotypes, and several valuable genetic resources of wild F. filiformis, and an effective strategy for guiding F. filiformis breeding based on the population structure and pedigree relationship in future.

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.

Reactive Oxygen Species Distribution Involved in Stipe Gradient Elongation in the Mushroom Flammulina filiformis

The mushroom stipe raises the pileus above the substrate into a suitable position for dispersing spores. The stipe elongates at different speeds along its length, with the rate of elongation decreasing in a gradient from the top to the base. However, the molecular mechanisms underlying stipe gradient elongation are largely unknown. Here, we used the model basidiomycete mushroom Flammulina filiformis to investigate the mechanism of mushroom stipe elongation and the role of reactive oxygen species (ROS) signaling in this process. Our results show that O₂^- and H₂O₂ exhibit opposite gradient distributions in the stipe, with higher O₂^- levels in the elongation region (ER), and higher H₂O₂ levels in the stable region (SR). Moreover, NADPH-oxidase-encoding genes are up-regulated in the ER, have a function in producing O₂^-, and positively regulate stipe elongation. Genes encoding manganese superoxide dismutase (MnSOD) are up-regulated in the SR, have a function in producing H₂O_2, and negatively regulate stipe elongation. Altogether, our data demonstrate that ROS (O₂^-/H₂O₂) redistribution mediated by NADPH oxidase and MnSODs is linked to the gradient elongation of the F. filiformis stipe.

FFGA1 Protein Is Essential for Regulating Vegetative Growth, Cell Wall Integrity, and Protection against Stress in Flammunina filiformis

Flammulina filiformis is a popular mushroom which has been regarded as a potential model fungus for mycelium growth, fruiting body development, and stress response studies. Based on a genome-wide search, four genes encoding heterotrimeric G protein α subunits were identified in F. filiformis. The data of conserved domain analysis showed that these genes contain only one subgroup I of Gα subunit (Gαi), similar to many other fungi. To explore the function of Gαi, FfGa1 over-expression (OE) and RNA interference (RNAi) strains were generated using the Agrobacterium tumefaciens-mediated transformation (ATMT) approach. RNAi strains showed remarkably reduced growth on PDA medium and sensitivity to cell wall-perturbing agents, with maximum growth inhibition, but showed better growth in response to hypertonic stress-causing agents, while OE strains exhibited more resistance to thermal stress and mycoparasite Trichoderma as compared to the wild-type and RNAi strains. Taken together, our results indicated that FfGa1 positively regulates hyphal extension, and is crucial for the maintenance of cell wall integrity and protection against biotic and abiotic (hypertonic and thermal) stress.

The Fvclp1 gene regulates mycelial growth and fruiting body development in edible mushroom Flammulina velutipes

Fruiting body development in Agaricomycetes represents the most complex and unclear process in the fungi. Mating type pathways (A and B) and transcription factors are important regulators in the sexual development of mushrooms. It is known that clampless1 (clp1) is an additional gene that participate under the homeodomain (HD) genes in the matA pathway and clp1 inactivation blocks clamps formation in Coprinopsis cinerea. In this study we identified and analyzed a homologous Fvclp1 gene in the edible mushroom Flammulina velutipes. The coding sequence of the Fvclp1 was 1011 bp without intron interruption, encoding a protein of 336 amino acids. To exhibit the role of Fvclp1 in clamp development and fruiting body formation, knockdown and overexpression mutants were prepared. No significant difference was observed in the monokaryotic hyphal morphology of overexpression and knockdown transformants. In the dikaryotic hyphae from the compatible crossings between the wild-type L22 strain and Fvclp1 knockdown or overexpression mutants, clamp connections developed. However, knockdown mutants could generate fewer fruiting bodies than the wild-type strain. On the contrary, reduced mycelial growth rate but improved fruiting ability was observed in the dikaryotic Fvclp1 overexpression mutants as compared to the wild-type strain. These results indicate that Fvclp1 is necessary and actively involved in fruiting body development in F. velutipes. Overall, these findings suggest that further studies on the function of Fvclp1 would advance our understanding of sexual reproduction and fruiting body development in edible mushrooms.

The Genomic and Transcriptomic Analyses of Floccularia luteovirens, a Rare Edible Fungus in the Qinghai-Tibet Plateau, Provide Insights into the Taxonomy Placement and Fruiting Body Formation

Floccularia luteovirens is a famous and precious edible mushroom (Huang Mogu) on the Qinghai–Tibet plateau that has a unique flavor and remarkable medical functions. Herein, we report a reference-grade 27 Mb genome of F. luteovirens containing 7068 protein-coding genes. The genome component and gene functions were predicted. Genome ontology enrichment and pathway analyses indicated the potential production capacity for terpenoids, polyketides and polysaccharides. Moreover, 16 putative gene clusters and 145 genes coding for secondary metabolites were obtained, including guadinomine and melleolides. In addition, phylogenetic and comparative genomic analyses shed light on the precise classification of F. luteovirens suggesting that it belongs to the genus Floccularia instead of Armillaria. RNA-sequencing and comparative transcriptomic analysis revealed differentially expressed genes during four developmental stages of F. luteovirens, that of which helps to identify important genes regulating fruiting body formation for strain modification. This study will provide insight into artificial cultivation and increase the production of useful metabolites.

Genome of Ganoderma Species Provides Insights Into the Evolution, Conifers Substrate Utilization, and Terpene Synthesis for Ganoderma tsugae

Ganoderma tsugae is an endemic medicinal mushroom in Northeast China, providing important source of pharmaceutical product. Comparing with other Ganoderma species, wild G. tsugae can utilize coniferous wood. However, functional genes related to medicinal component synthesis and the genetic mechanism of conifer substrate utilization is still obscure. Here, we assembled a high-quality G. tsugae genome with 18 contigs and 98.5% BUSCO genes and performed the comparative genomics with other Ganoderma species. G. tsugae diverged from their common ancestor of G. lingzhi and G. sinense about 21 million years ago. Genes in G. tsugae-specific and G. tsugae-expanded gene families, such as salh, phea, cyp53a1, and cyp102a, and positively selected genes, such as glpk and amie, were functionally enriched in plant-pathogen interaction, benzoate degradation, and fanconi anemia pathway. Those functional genes might contribute to conifer substrate utilization of G. tsugae. Meanwhile, gene families in the terpene synthesis were identified and genome-wide SNP variants were detected in population. Finally, the study provided valuable genomic resources and offered useful hints for the functional gene mapping and investigation of key gene contributing to conifer cultivation substrate utilization and medicinal component biosynthesis.

Screening candidate genes related to volatile synthesis in shiitake mushrooms and construction of regulatory networks to effectively improve mushroom aroma

BACKGROUND

Metabolite formation is a physiological stress response during the growth and development of shiitake mushrooms (Lentinula edodes). The characteristic flavor metabolites are important quality components in shiitake mushrooms. To investigate the formation mechanisms of characteristic flavor metabolites, transcriptome analyses were performed on shiitake mushrooms harvested at different growth stages.

RESULTS

In total, 30 genes related to the synthesis of characteristic volatiles of mushrooms were identified via screening. Through KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis of the selected genes and correlation analyses of gene expressions, the main volatile synthesis pathways were determined as histidine metabolism, glutathione metabolism and biosynthesis of unsaturated fatty acids. Gene cluster and correlation analyses were performed to clarify the combined effects of different genes in the enzymatic reactions. Further, a correlation network of candidate genes was built based on the gene expression levels.

CONCLUSION

The activities of flavor synthases and the content of characteristic flavor metabolites were analyzed; the enzyme activity changes and metabolic product distribution sites were clarified. A synthesis and regulation network was constructed for the candidate genes and characteristic volatiles, and information was obtained for 16 hub genes. Moreover, it was essential to identify and characterize the key genes and synthases involved in the synthesis of the characteristic volatiles of mushrooms. This information provides us with a better understanding of the biosynthesis and regulation of the volatiles, which will lay the foundation for improving the quality of shiitake mushrooms. © 2021 Society of Chemical Industry.

Analysis of the Genome Sequence of Strain GiC-126 of Gloeostereum incarnatum with Genetic Linkage Map

Gloeostereum incarnatum has edible and medicinal value and was first cultivated and domesticated in China. We sequenced the G. incarnatum monokaryotic strain GiC-126 on an Illumina HiSeq X Ten system and obtained a 34.52-Mb genome assembly sequence that encoded 16,895 predicted genes. We combined the GiC-126 genome with the published genome of G. incarnatum strain CCMJ2665 to construct a genetic linkage map (GiC-126 genome) that had 10 linkage groups (LGs), and the 15 assembly sequences of CCMJ2665 were integrated into 8 LGs. We identified 1912 simple sequence repeat (SSR) loci and detected 700 genes containing 768 SSRs in the genome; 65 and 100 of them were annotated with gene ontology (GO) terms and KEGG pathways, respectively. Carbohydrate-active enzymes (CAZymes) were identified in 20 fungal genomes and annotated; among them, 144 CAZymes were annotated in the GiC-126 genome. The A mating-type locus (MAT-A) of G. incarnatum was located on scaffold885 at 38.9 cM of LG1 and was flanked by two homeodomain (HD1) genes, mip and beta-fg. Fourteen segregation distortion markers were detected in the genetic linkage map, all of which were skewed toward the parent GiC-126. They formed three segregation distortion regions (SDR1-SDR3), and 22 predictive genes were found in scaffold1920 where three segregation distortion markers were located in SDR1. In this study, we corrected and updated the genomic information of G. incarnatum. Our results will provide a theoretical basis for fine gene mapping, functional gene cloning, and genetic breeding the follow-up of G. incarnatum.

Symbiotic nitrogen fixation in the reproductive structures of a basidiomycete fungus

Nitrogen (N) fixation is a driving force for the formation of symbiotic associations between N₂-fixing bacteria and eukaryotes.¹ Limited examples of these associations are known in fungi, and none with sexual structures of non-lichenized species.^2-6 The basidiomycete Guyanagaster necrorhizus is a sequestrate fungus endemic to the Guiana Shield.⁷ Like the root rot-causing species in its sister genera Armillaria and Desarmillaria, G. necrorhizus sporocarps fruit from roots of decaying trees (Figures 1A-1C),⁸ and genome sequencing is consistent with observations that G. necrorhizus is a white-rotting decomposer. This species also represents the first documentation of an arthropod-dispersed sequestrate fungus. Numerous species of distantly related wood-feeding termites, which scavenge for N-rich food, feed on the mature spore-bearing tissue, or gleba, of G. necrorhizus. During feeding, mature spores adhere to termites for subsequent dispersal.⁹ Using chemical assays, isotope analysis, and high-throughput sequencing, we show that the sporocarps harbor actively N₂-fixing Enterobacteriaceae species and that the N content within fungal tissue increases with maturation. Untargeted proteomic profiling suggests that ATP generation in the gleba is accomplished via fermentation. The use of fermentation-an anaerobic process-indicates that the sporocarp environment is anoxic, likely an adaptation to protect the oxygen-sensitive nitrogenase enzyme. Sporocarps also have a thick outer covering, possibly to limit oxygen diffusion. The enriched N content within mature sporocarps may offer a dietary inducement for termites in exchange for spore dispersal. These results show that the flexible metabolic capacity of fungi may facilitate N₂-fixing associations, as well as higher-level organismal associations.

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

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

Comparative Analysis of Carbohydrate Active Enzymes in the Flammulina velutipes var. lupinicola Genome

The purpose of this study was to determine the genome sequence of Flammulina velutipes var. lupinicola based on next-generation sequencing (NGS) and to identify the genes encoding carbohydrate-active enzymes (CAZymes) in the genome. The optimal assembly (71 kmer) based on ABySS de novo assembly revealed a total length of 33,223,357 bp (49.53% GC content). A total of 15,337 gene structures were identified in the F. velutipes var. lupinicola genome using ab initio gene prediction method with Funannotate pipeline. Analysis of the orthologs revealed that 11,966 (96.6%) out of the 15,337 predicted genes belonged to the orthogroups and 170 genes were specific for F. velutipes var. lupinicola. CAZymes are divided into six classes: auxiliary activities (AAs), glycosyltransferases (GTs), carbohydrate esterases (CEs), polysaccharide lyases (PLs), glycoside hydrolases (GHs), and carbohydrate-binding modules (CBMs). A total of 551 genes encoding CAZymes were identified in the F. velutipes var. lupinicola genome by analyzing the dbCAN meta server database (HMMER, Hotpep, and DIAMOND searches), which consisted of 54-95 AAs, 145-188 GHs, 55-73 GTs, 6-19 PLs, 13-59 CEs, and 7-67 CBMs. CAZymes can be widely used to produce bio-based products (food, paper, textiles, animal feed, and biofuels). Therefore, information about the CAZyme repertoire of the F. velutipes var. lupinicola genome will help in understanding the lignocellulosic machinery and in-depth studies will provide opportunities for using this fungus for biotechnological and industrial applications.

Mycena genomes resolve the evolution of fungal bioluminescence

Mushroom-forming fungi in the order Agaricales represent an independent origin of bioluminescence in the tree of life; yet the diversity, evolutionary history, and timing of the origin of fungal luciferases remain elusive. We sequenced the genomes and transcriptomes of five bonnet mushroom species (Mycena spp.), a diverse lineage comprising the majority of bioluminescent fungi. Two species with haploid genome assemblies ∼150 Mb are among the largest in Agaricales, and we found that a variety of repeats between Mycena species were differentially mediated by DNA methylation. We show that bioluminescence evolved in the last common ancestor of mycenoid and the marasmioid clade of Agaricales and was maintained through at least 160 million years of evolution. Analyses of synteny across genomes of bioluminescent species resolved how the luciferase cluster was derived by duplication and translocation, frequently rearranged and lost in most Mycena species, but conserved in the Armillaria lineage. Luciferase cluster members were coexpressed across developmental stages, with the highest expression in fruiting body caps and stipes, suggesting fruiting-related adaptive functions. Our results contribute to understanding a de novo origin of bioluminescence and the corresponding gene cluster in a diverse group of enigmatic fungal species.

Immunomodulatory Effects of Edible and Medicinal Mushrooms and Their Bioactive Immunoregulatory Products

Mushrooms have been valued as food and health supplements by humans for centuries. They are rich in dietary fiber, essential amino acids, minerals, and many bioactive compounds, especially those related to human immune system functions. Mushrooms contain diverse immunoregulatory compounds such as terpenes and terpenoids, lectins, fungal immunomodulatory proteins (FIPs) and polysaccharides. The distributions of these compounds differ among mushroom species and their potent immune modulation activities vary depending on their core structures and fraction composition chemical modifications. Here we review the current status of clinical studies on immunomodulatory activities of mushrooms and mushroom products. The potential mechanisms for their activities both in vitro and in vivo were summarized. We describe the approaches that have been used in the development and application of bioactive compounds extracted from mushrooms. These developments have led to the commercialization of a large number of mushroom products. Finally, we discuss the problems in pharmacological applications of mushrooms and mushroom products and highlight a few areas that should be improved before immunomodulatory compounds from mushrooms can be widely used as therapeutic agents.

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

Background

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

Results

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

Conclusions

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

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.

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.

Characterization and expression pattern analysis of pheromone receptor-like genes in Winter Mushroom Flammulina filiformis

Pheromone receptor-like genes (PRLGs) belong to the G protein-coupled receptors (GPCRs) family that interacts with biotic and abiotic stimulants and transmits signals to intracellular downstream pathways in eukaryotic cells. In this study, we investigated the structure and expressions patterns of PRLGs in Winter Mushroom Flammulina filiformis. Based on the alignment analysis, the structure of PRLGs was found conserved in F. filiformis strains expect few single-nucleotide polymorphism (SNP) sites. Six PRLGs were found at five different unlinked loci, scattered in the genomes of F. filiformis strains. These genes contain 2-5 introns; however, the introns were not found in the same relative positions regarding the encoded protein sequences in tested strains of F. filiformis. Three conserved motifs were identified in peptides structures of PRLGs, however, FfSte3.s6 contained only two types, suggests its difference in evolution and function. We have further analyzed the expression patterns of each PRLGs in different developmental stages of the fruiting body in F. filiformis by quantitative real-time polymerase chain reaction (qRT-PCR). The results exhibited expression variation of PRLGs at different developmental stages of the F. filiformis. Especially, FfSte3.s1 and FfSte3.s2 exhibited maximum expression level in mycelia stage. Other PRLGs exhibited high expression level in fruiting body stages. This study suggests that PRLGs could be vital genes involving in fruiting body development in F. filiformis. However, further studies could be performed to reveal their specific functional pathways in the fruiting body development.

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.

A WD40 Protein Encoding Gene Fvcpc2 Positively Regulates Mushroom Development and Yield in Flammulina velutipes

Ascomycota and Basidiomycota are two closely related phyla and fungi in two phyla share some common morphological developmental process during fruiting body formation. In Neurospora crassa, the Gβ-like protein CPC-2 with a seven-WD40 repeat domain was previously reported. By transforming CPC-2 ortholog encoding genes, from 7 different fungal species across Ascomycota and Basidiomycota, into the cpc-2 deletion mutant of N. crassa, we demonstrate that all tested CPC-2 ortholog genes were able to complement the defects of the cpc-2 deletion mutant in sexual development, indicating that CPC-2 proteins from Ascomycota and Basidiomycota have the similar cellular function. Using Flammulina velutipes as a model system for mushroom species, the CPC-2 ortholog FvCPC2 was characterized. Fvcpc2 increased transcription during fruiting body development. Knockdown of Fvcpc2 by RNAi completely impaired fruiting body formation. In three Fvcpc2 knockdown mutants, transcriptional levels of genes encoding adenylate cyclase and protein kinase A catalytic subunit were significantly lower and colony growth became slower than wild type. The addition of cAMP or the PKA-activator 8-Bromo-cAMP into the medium restored the Fvcpc2 knockdown mutants to the wild-type colony growth phenotype, suggesting that the involvement of cAMP production in the regulatory mechanisms of FvCPC2. Knockdown of Fvcpc2 also weakened transcriptional responses to sexual development induction by some genes related to fruiting body development, including 4 jacalin-related lectin encoding genes, 4 hydrophobin encoding genes, and 3 functionally-unknown genes, suggesting the participation of these genes in the mechanisms by which FvCPC2 regulates fruiting body development. All three Fvcpc2 overexpression strains displayed increased mushroom yield and shortened cultivation time compared to wild type, suggesting that Fvcpc2 can be a promising reference gene for Winter Mushroom breeding. Since the orthologs of FvCPC2 were highly conserved and specifically expressed during fruiting body development in different edible mushrooms, genes encoding FvCPC2 orthologs in other mushroom species also have potential application in breeding.

Genome Sequence Analysis of Auricularia heimuer Combined with Genetic Linkage Map

Auricularia heimuer is one of the most popular edible fungi in China. In this study, the whole genome of A. heimuer was sequenced on the Illumina HiSeq X system and compared with other mushrooms genomes. As a wood-rotting fungus, a total of 509 carbohydrate-active enzymes (CAZymes) were annotated in order to explore its potential capabilities on wood degradation. The glycoside hydrolases (GH) family genes in the A. heimuer genome were more abundant than the genes in the other 11 mushrooms genomes. The A. heimuer genome contained 102 genes encoding class III, IV, and V ethanol dehydrogenases. Evolutionary analysis based on 562 orthologous single-copy genes from 15 mushrooms showed that Auricularia formed an early independent branch of Agaricomycetes. The mating-type locus of A. heimuer was located on linkage group 8 by genetic linkage analysis. By combining the genome sequence analysis with the genetic linkage map, the mating-type locus of A. heimuer was located on scaffold45 and consisted of two subunits, α and β. Each subunit consisted of a pair of homeodomain mating-type protein genes HD1 and HD2. The mapping revealed conserved synteny at the whole mating-type loci and mirror symmetry relations near the β subunit between A. heimuer and Exidia glandulosa. This study proposed the potential for the bioethanol production by consolidated bioprocessing of A. heimuer. It will promote understanding of the lignocellulose degradation system and facilitate more efficient conversion of the agricultural wastes used for mushroom cultivation. It also will advance the research on the fruiting body development and evolution of A. heimuer.