The genome sequence of the commercially cultivated mushroom Agrocybe aegerita reveals a conserved repertoire of fruiting-related genes and a versatile suite of biopolymer-degrading enzymes

Lessons on fruiting body morphogenesis from genomes and transcriptomes of Agaricomycetes

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

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.

Towards Understanding the Function of Aegerolysins

Aegerolysins are remarkable proteins. They are distributed over the tree of life, being relatively widespread in bacteria and fungi, but also present in some insects, plants, protozoa, and viruses. Despite their abundance in cells of certain developmental stages and their presence in secretomes, only a few aegerolysins have been studied in detail. Their function, in particular, is intriguing. Here, we summarize previously published findings on the distribution, molecular interactions, and function of these versatile aegerolysins. They have very diverse protein sequences but a common fold. The machine learning approach of the AlphaFold2 algorithm, which incorporates physical and biological knowledge of protein structures and multisequence alignments, provides us new insights into the aegerolysins and their pore-forming partners, complemented by additional genomic support. We hypothesize that aegerolysins are involved in the mechanisms of competitive exclusion in the niche.

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

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

Genome and Comparative Transcriptome Dissection Provide Insights Into Molecular Mechanisms of Sclerotium Formation in Culinary-Medicinal Mushroom Pleurotus tuber-regium

Pleurotus tuber-regium is an edible and medicinal sclerotium-producing mushroom. The sclerotia of this mushroom also serve as food and folk medicine. Based on the description of its monokaryon genome, sequenced with Illumina and PacBio sequencing technologies, comparative transcriptomic analysis using RNA sequencing (RNA-seq) was employed to study its mechanism of sclerotium formation. The de novo assembled genome is 35.82 Mb in size with a N50 scaffold size of 4.29 Mb and encodes 12,173 putative proteins. Expression analysis demonstrated that 1,146 and 1,249 genes were upregulated and downregulated with the formation of sclerotia, respectively. The differentially expressed genes were associated with substrate decomposition, the oxidation-reduction process, cell wall synthesis, and other biological processes in P. tuber-regium. These genomic and transcriptomic resources provide useful information for the mechanism underlying sclerotium formation in P. tuber-regium.

Characterization of cheesy odor formed during fermentation of soy drink with Agrocybe aegerita

The market for plant protein-based substitutes for cheeses is growing, but the sensory properties are distinctively different from the original products. Hence, natural and vegan cheesy flavors are needed to aromatize the products. A cheesy, sweaty and parmesan-like aroma was produced by fermentation of soy drink with Agrocybe aegerita. Aroma dilution analysis revealed short-chain fatty acids (SCFAs) as main influencing cheesy odorants analyzed by gas chromatography-mass spectrometry-olfactometry. In comparison to the five cheese varieties, the SCFA profile of the fermented soy drink revealed similarities with Parmesan and Emmental cheese. Meanwhile, principal component analysis showed an approximation of the aroma profile after fermentation with A. aegerita to those of cheeses. 3-Methylbutanoic acid was synthesized from the protein fraction, while the oil fraction contributed to the formation of unbranched SCFAs like butanoic acid. Accordingly, the production of these compounds can be increased by addition of the fractions.

Peroxide-Mediated Oxygenation of Organic Compounds by Fungal Peroxygenases

Unspecific peroxygenases (UPOs), whose sequences can be found in the genomes of thousands of filamentous fungi, many yeasts and certain fungus-like protists, are fascinating biocatalysts that transfer peroxide-borne oxygen (from H2O2 or R-OOH) with high efficiency to a wide range of organic substrates, including less or unactivated carbons and heteroatoms. A twice-proline-flanked cysteine (PCP motif) typically ligates the heme that forms the heart of the active site of UPOs and enables various types of relevant oxygenation reactions (hydroxylation, epoxidation, subsequent dealkylations, deacylation, or aromatization) together with less specific one-electron oxidations (e.g., phenoxy radical formation). In consequence, the substrate portfolio of a UPO enzyme always combines prototypical monooxygenase and peroxidase activities. Here, we briefly review nearly 20 years of peroxygenase research, considering basic mechanistic, molecular, phylogenetic, and biotechnological aspects.

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.

A review of the therapeutic and biological effects of edible and wild mushrooms

Throughout history, mushrooms have occupied an inseparable part of the diet in many countries. Mushrooms are considered a rich source of phytonutrients such as polysaccharides, dietary fibers, and other micronutrients, in addition to various essential amino acids, which are building blocks of vital proteins. In general, mushrooms offer a wide range of health benefits with a large spectrum of pharmacological properties, including antidiabetic, antioxidative, antiviral, antibacterial, osteoprotective, nephroprotective, hepatoprotective, etc. Both wild edible and medicinal mushrooms possess strong therapeutic and biological activities, which are evident from their in vivo and in vitro assays. The multifunctional activities of the mushroom extracts and the targeted potential of each of the compounds in the extracts have a broad range of applications, especially in the healing and repair of various organs and cells in humans. Owing to the presence of the aforementioned properties and rich phytocomposition, mushrooms are being used in the production of nutraceuticals and pharmaceuticals. This review aims to provide a clear insight on the commercially cultivated, wild edible, and medicinal mushrooms with comprehensive information on their phytochemical constituents and properties as part of food and medicine for futuristic exploitation. Future outlook and prospective challenges associated with the cultivation and processing of these medicinal mushrooms as functional foods are also discussed.

Genome Size Versus Genome Assemblies: Are the Genomes Truly Expanded in Polyploid Fungal Symbionts?

Each day, as the amount of genomic data and bioinformatics resources grows, researchers are increasingly challenged with selecting the most appropriate approach to analyze their data. In addition, the opportunity to undertake comparative genomic analyses is growing rapidly. This is especially true for fungi due to their small genome sizes (i.e., mean 1C = 44.2 Mb). Given these opportunities and aiming to gain novel insights into the evolution of mutualisms, we focus on comparing the quality of whole genome assemblies for fungus-growing ants cultivars (Hymenoptera: Formicidae: Attini) and a free-living relative. Our analyses reveal that currently available methodologies and pipelines for analyzing whole-genome sequence data need refining. By using different genome assemblers, we show that the genome assembly size depends on what software is used. This, in turn, impacts gene number predictions, with higher gene numbers correlating positively with genome assembly size. Furthermore, the majority of fungal genome size data currently available are based on estimates derived from whole-genome assemblies generated from short-read genome data, rather than from the more accurate technique of flow cytometry. Here, we estimated the haploid genome sizes of three ant fungal symbionts by flow cytometry using the fungus Pleurotus ostreatus (Jacq.) P. Kumm. (1871) as a calibration standard. We found that published genome sizes based on genome assemblies are 2.5- to 3-fold larger than our estimates based on flow cytometry. We, therefore, recommend that flow cytometry is used to precalibrate genome assembly pipelines, to avoid incorrect estimates of genome sizes and ensure robust assemblies.

Biotransformation of [U-13C]linoleic acid suggests two independent ketonic- and aldehydic cycles within C8-oxylipin biosynthesis in Cyclocybe aegerita (V. Brig.) Vizzini

Although the typical aroma contributing compounds in fungi of the phylum Basidiomycota are known for decades, their biosynthetic pathways are still unclear. Amongst these volatiles, C8-compounds are probably the most important ones as they function, in addition to their specific perception of fungal odour, as oxylipins. Previous studies focused on C8-oxylipin production either in fruiting bodies or mycelia. However, comparisons of the C8-oxylipin biosynthesis at different developmental stages are scarce, and the biosynthesis in basidiospores was completely neglected. In this study, we addressed this gap and were able to show that the biosynthesis of C8-oxylipins differs strongly between different developmental stages. The comparison of mycelium, primordia, young fruiting bodies, mature fruiting bodies, post sporulation fruiting bodies and basidiospores revealed that the occurance of the two main C8-oxylipins octan-3-one and oct-1-en-3-ol distinguished in different stages. Whereas oct-1-en-3-ol levels peaked in the mycelium and decreased with ongoing maturation, octan-3-one levels increased during maturation. Furthermore, oct-2-en-1-ol, octan-1-ol, oct-2-enal, octan-3-ol, oct-1-en-3-one and octanal contributed to the C8-oxylipins but with drastically lower levels. Biotransformations with [U-13C]linoleic acid revealed that early developmental stages produced various [U-13C]oxylipins, whereas maturated developmental stages like post sporulation fruiting bodies and basidiospores produced predominantly [U-13C]octan-3-one. Based on the distribution of certain C8-oxylipins and biotransformations with putative precursors at different developmental stages, two distinct biosynthetic cycles were deduced with oct-2-enal (aldehydic-cycle) and oct-1-en-3-one (ketonic-cycle) as precursors.

Transcriptional response of mushrooms to artificial sun exposure

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

Genome reduction and relaxed selection is associated with the transition to symbiosis in the basidiomycete genus Podaxis

Insights into the genomic consequences of symbiosis for basidiomycete fungi associated with social insects remain sparse. Capitalizing on viability of spores from centuries-old herbarium specimens of free-living, facultative, and specialist termite-associated Podaxis fungi, we obtained genomes of 10 specimens, including two type species described by Linnaeus >240 years ago. We document that the transition to termite association was accompanied by significant reductions in genome size and gene content, accelerated evolution in protein-coding genes, and reduced functional capacities for oxidative stress responses and lignin degradation. Functional testing confirmed that termite specialists perform worse under oxidative stress, while all lineages retained some capacity to cleave lignin. Mitochondrial genomes of termite associates were significantly larger; possibly driven by smaller population sizes or reduced competition, supported by apparent loss of certain biosynthetic gene clusters. Our findings point to relaxed selection that mirrors genome traits observed among obligate endosymbiotic bacteria of many insects.

The Pacific Tree-Parasitic Fungus Cyclocybe parasitica Exhibits Monokaryotic Fruiting, Showing Phenotypes Known from Bracket Fungi and from Cyclocybe aegerita

Cyclocybe parasitica is a wood-destroying parasitic edible mushroom growing on diverse broad-leafed trees in New Zealand and other Pacific areas. Recent molecular systematics of European Cyclocybe aegerita, a newly delimited Asian phylum and of related species, corroborated the distinction of the chiefly saprobic cultivated edible mushroom C. aegerita from C. parasitica. Here, we show that C. parasitica exhibits a morpho-physiological trait characteristic to its European cousin, i.e., monokaryotic fruiting sensu stricto (basidiome formation without mating). Monokaryotic fruiting structures formed by C. parasitica ICMP 11668-derived monokaryons were categorized into four phenotypes. One of them displays ulcer-like structures previously reported from bracket fungi. Histology of dikaryotic and monokaryotic C. parasitica fruiting structures revealed anatomical commonalities and differences between them, and towards monokaryotic fruiting structures of C. aegerita. Mating experiments with C. parasitica strains representative of each fruiting phenotype identified compatible sibling monokaryons. Given reports on hypothetically monokaryotic basidiome field populations of ‘C. aegerita sensu lato’, it seems worthwhile to prospectively investigate whether monokaryotic fruiting s.str. occurs in nature. Sampling from such populations including karyotyping, comparative -omics, and competition assays may help to answer this question and provide evidence whether this trait may confer competitive advantages to a species capable of it.

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.

Ageritin from Pioppino Mushroom: The Prototype of Ribotoxin-Like Proteins, a Novel Family of Specific Ribonucleases in Edible Mushrooms

Ageritin is a specific ribonuclease, extracted from the edible mushroom Cyclocybe aegerita (synonym Agrocybe aegerita), which cleaves a single phosphodiester bond located within the universally conserved alpha-sarcin loop (SRL) of 23-28S rRNAs. This cleavage leads to the inhibition of protein biosynthesis, followed by cellular death through apoptosis. The structural and enzymatic properties show that Ageritin is the prototype of a novel specific ribonucleases family named 'ribotoxin-like proteins', recently found in fruiting bodies of other edible basidiomycetes mushrooms (e.g., Ostreatin from Pleurotus ostreatus, Edulitins from Boletus edulis, and Gambositin from Calocybe gambosa). Although the putative role of this toxin, present in high amount in fruiting body (>2.5 mg per 100 g) of C. aegerita, is unknown, its antifungal and insecticidal actions strongly support a role in defense mechanisms. Thus, in this review, we focus on structural, biological, antipathogenic, and enzymatic characteristics of this ribotoxin-like protein. We also highlight its biological relevance and potential biotechnological applications in agriculture as a bio-pesticide and in biomedicine as a therapeutic and diagnostic agent.

Bioinformatics-aided identification, characterization and applications of mushroom linalool synthases

Enzymes empower chemical industries and are the keystone for metabolic engineering. For example, linalool synthases are indispensable for the biosynthesis of linalool, an important fragrance used in 60–80% cosmetic and personal care products. However, plant linalool synthases have low activities while expressed in microbes. Aided by bioinformatics analysis, four linalool/nerolidol synthases (LNSs) from various Agaricomycetes were accurately predicted and validated experimentally. Furthermore, we discovered a linalool synthase (Ap.LS) with exceptionally high levels of selectivity and activity from Agrocybe pediades, ideal for linalool bioproduction. It effectively converted glucose into enantiopure (R)-linalool in Escherichia coli, 44-fold and 287-fold more efficient than its bacterial and plant counterparts, respectively. Phylogenetic analysis indicated the divergent evolution paths for plant, bacterial and fungal linalool synthases. More critically, structural comparison provided catalytic insights into Ap.LS superior specificity and activity, and mutational experiments validated the key residues responsible for the specificity.

Zhang et al. identified four linalool/nerolidol synthases from fungi using bioinformatics and phylogenetic analysis and validated their functions with in vitro and in vivo methods. One of them is a rare and highly specific monoterpene synthase and responsible for impressive titres of the commercially sought-after fragrance (R)-linalool.

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

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

Multilocus phylogeny- and fruiting feature-assisted delimitation of European Cyclocybe aegerita from a new Asian species complex and related species

Cyclocybe aegerita (synonym: Agrocybe aegerita) is a widely cultivated edible and reportedly almost cosmopolitan mushroom species that serves as a model fungus for basidiome formation and as producer of useful natural products and enzymes. Focusing on strains from different continents, here, we present a phylogenetic analysis of this species and some adjacent taxa that employs four phylogenetic markers. In addition, we tested the strains’ capability to fructify on agar media. Our analysis reveals that “C. aegerita sensu lato” splits up into the following two well-supported monophyletic geographic lineages: a European clade and an Asian clade. The European one is closely associated with the Chinese species Cyclocybe salicaceicola. In contrast, the Asian lineage, which we preliminarily designate as Cyclocybe chaxingu agg., may comprise several species (species complex) and clusters with the Pacific species Cyclocybe parasitica (New Zealand). In addition, fruiting properties differ across C. aegerita and its Asian and Pacific relatives; however, strains from the Asian clade and C. parasitica tend to form larger basidiomes with relatively big caps and long stipes and strains from the European clade exhibit a more variable fruiting productivity with the tendency to form more basidiomes, with smaller caps and shorter stipes. Moreover, some strains showed individual fruiting patterns, such as the preference to fruit where they were exposed to injuring stimuli. In conclusion, the delimitation of the newly delimited Asian species complex from our multilocus phylogeny of “C. aegerita sensu lato”, which is supported by phenotypic data, depicts an exemplary case of biogeographic diversity within a previously thought homogeneous species of near worldwide distribution.

Gene Organization, Expression, and Localization of Ribotoxin-Like Protein Ageritin in Fruiting Body and Mycelium of Agrocybe aegerita

The edible mushroom Agrocybe aegerita produces a ribotoxin-like protein known as Ageritin. In this work, the gene encoding Ageritin was characterized by sequence analysis. It contains several typical features of fungal genes such as three short introns (60, 55 and 69 bp) located at the 5' region of the coding sequence and typical splice junctions. This sequence codes for a precursor of 156 amino acids (~17-kDa) containing an additional N-terminal peptide of 21 amino acid residues, absent in the purified toxin (135 amino acid residues; ~15-kDa). The presence of 17-kDa and 15-kDa forms was investigated by Western blot in specific parts of fruiting body and in mycelia of A. aegerita. Data show that the 15-kDa Ageritin is the only form retrieved in the fruiting body and the principal form in mycelium. The immunolocalization by confocal laser scanning microscopy and transmission electron microscopy proves that Ageritin has vacuolar localization in hyphae. Coupling these data with a bioinformatics approach, we suggest that the N-terminal peptide of Ageritin (not found in the purified toxin) is a new signal peptide in fungi involved in intracellular routing from endoplasmic reticulum to vacuole, necessary for self-defense of A. aegerita ribosomes from Ageritin toxicity.

Degradative Ability of Mushrooms Cultivated on Corn Silage Digestate

The current management practice of digestate from biogas plants involves its use for land application as a fertilizer. Nevertheless, the inadequate handling of digestate may cause environmental risks due to losses of ammonia, methane and nitrous oxide. Therefore, the key goals of digestate management are to maximize its value by developing new digestate products, reducing its dependency on soil application and the consequent air pollution. The high nitrogen and lignin content in solid digestate make it a suitable substrate for edible and medicinal mushroom cultivation. To this aim, the mycelial growth rate and degradation capacity of the lignocellulosic component from corn silage digestate, undigested wheat straw and their mixture were investigated on Cyclocybe aegerita, Coprinus comatus, Morchella importuna, Pleurotus cornucopiae and Pleurotus ostreatus. The structural modification of the substrates was performed by using attenuated total reflectance-Fourier transform infrared (ATR-FTIR) spectroscopy. Preliminary in vitro results demonstrated the ability of P. ostreatus, P. cornucopiae and M. importuna to grow and decay hemicellulose and lignin of digestate. Cultivation trials were carried out on C. aegerita, P. cornucopiae and P. ostreatus. Pleurotus ostreatus showed the highest biological efficiency and fruiting body production in the presence of the digestate; moreover, P. ostreatus and P. cornucopiae were able to degrade the lignin. These results provide attractive perspectives both for more sustainable digestate management and for the improvement of mushroom cultivation efficiency.

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

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

Genome Assembly and Pathway Analysis of Edible Mushroom Agrocybe cylindracea

Agrocybe cylindracea, an edible mushroom, is widely cultivated for its abundance of nutrients and flavor, and many of its metabolites are reported to have beneficial roles, such as medicinal effects on tumors and chronical illnesses. However, the lack of genomic information has hindered further molecular studies on this fungus. Here, we present a genome assembly of A. cylindracea together with comparative genomics and pathway analyses of Agaricales species. The draft, generated from both next-generation sequencing (NGS) and single-molecule real-time (SMRT) sequencing platforms to overcome high genetic heterozygosity, is composed of a 56.5 Mb sequence and 15,384 predicted genes. This mushroom possesses a complex reproductive system, including tetrapolar heterothallic and secondary homothallic mechanisms, and harbors several hydrolases and peptidases for gradual and effective degradation of various carbon sources. Our pathway analysis reveals complex processes involved in the biosynthesis of polysaccharides and other active substances, including B vitamins, unsaturated fatty acids, and N-acetylglucosamine. RNA-seq data show that A. cylindracea stipes tend to synthesize carbohydrate for carbon sequestration and energy storage, whereas pilei are more active in carbon utilization and unsaturated fatty acid biosynthesis. These results reflect diverse functions of the two anatomical structures of the fruiting body. Our comprehensive genomic and transcriptomic data, as well as preliminary comparative analyses, provide insights into the molecular details of the medicinal effects in terms of active compounds and nutrient components.

Agrocybe aegerita Serves As a Gateway for Identifying Sesquiterpene Biosynthetic Enzymes in Higher Fungi

Terpenoids constitute a structurally diverse group of natural products with wide applications in the pharmaceutical, nutritional, flavor and fragrance industries. Fungi are known to produce a large variety of terpenoids, yet fungal terpene synthases remain largely unexploited. Here, we report the sesquiterpene network and gene clusters of the black poplar mushroom Agrocybe aegerita. Among 11 putative sesquiterpene synthases (STSs) identified in its genome, nine are functional, including two novel synthases producing viridiflorol and viridiflorene. On this basis, an additional 1133 STS homologues from higher fungi have been curated and used for a sequence similarity network to probe isofunctional STS groups. With the focus on two STS groups, one producing viridiflorene/viridiflorol and one Δ⁶-protoilludene, the isofunctionality was probed and verified. Three new Δ⁶-protoilludene synthases and two new viridflorene/viridiflorol synthases from five different fungi were correctly predicted. The study herein serves as a fundamental predictive framework for the discovery of fungal STSs and biosynthesis of novel terpenoids. Furthermore, it becomes clear that fungal STS function differs between the phyla Ascomycota and Basidiomycota with the latter phylum being more dominant in the overall number and variability. This study aims to encourage the scientific community to further work on fungal STS and the products, biological functions, and potential applications of this vast source of natural products.

Selection and validation of reliable reference genes for Tolypocladium guangdongense gene expression analysis under differentially developmental stages and temperature stresses

Tolypocladium guangdongense, formerly known as Cordyceps guangdongensis, is a widely cultivated fungus of the Cordyceps s.l. species that has been investigated over the last 12 years. It has the potential to be used in a number of applications in the health and pharmaceutical industries for it has shown its high nutritional and medicinal values according to previous animal studies. qRT-PCR (quantitative reverse transcription polymerase chain reaction) is extensively used to analyze the expression pattern and molecular mechanisms of functional genes under differentially experimental conditions. The expression stability of reference genes used for normalization determines the reliability of qRT-PCR results, indicating the importance of selection and validation of reference genes before gene expression analysis. In the present study, three statistical algorithms, geNorm, NormFinder and BestKeeper, were used for analyzing the expression stability of nineteen candidate reference genes (CRGs) in T. guangdongense. Investigation were carried out under differentially experimental conditions, which included differentially developmential stages (mycelia, primordia, young and mature fruiting bodies), different carbon sources, cold and heat stresses. The results showed that histone H4 and tubulin beta chain 2 (β-tub2) were the most and least stable genes, respectively, across all the experimental samples. Moreover, analysis of individual data sets exhibited different stability and expression profiles of reference genes. The vacuolar protein sorting gene VPS was the most stable gene expressed under the differentially developmental stages and temperature stresses, whereas H4 was the most stably expressed gene under different carbon sources. Therefore, it can be proposed that VPS and H4 are the preferred reference genes for normalization of gene expression under different experimental conditions. The results of our present study will enable more accurate evaluation of gene expression in T. guangdongense using the optimal reference gene for qRT-PCR analysis.

Fungal Adaptation to the Advanced Stages of Wood Decomposition: Insights from the Steccherinum ochraceum

Steccherinum ochraceum is a white rot basidiomycete with wide ecological amplitude. It occurs in different regions of Russia and throughout the world, occupying different climatic zones. S. ochraceum colonizes stumps, trunks, and branches of various deciduous (seldom coniferous) trees. As a secondary colonizing fungus, S. ochraceum is mainly observed at the late decay stages. Here, we present the de novo assembly and annotation of the genome of S. ochraceum, LE-BIN 3174. This is the 8th published genome of fungus from the residual polyporoid clade and the first from the Steccherinaceae family. The obtained genome provides a first glimpse into the genetic and enzymatic mechanisms governing adaptation of S. ochraceum to an ecological niche of pre-degraded wood. It is proposed that increased number of carbohydrate-active enzymes (CAZymes) belonging to the AA superfamily and decreased number of CAZymes belonging to the GH superfamily reflects substrate preferences of S. ochraceum. This proposition is further substantiated by the results of the biochemical plate tests and exoproteomic study, which demonstrates that S. ochraceum assumes the intermediate position between typical primary colonizing fungi and litter decomposers or humus saprotrophs. Phylogenetic analysis of S. ochraceum laccase and class II peroxidase genes revealed the distinct evolutional origin of these genes in the Steccherinaceae family.

Heterologous Production and Functional Characterization of Ageritin, a Novel Type of Ribotoxin Highly Expressed during Fruiting of the Edible Mushroom Agrocybe aegerita

Fungi produce various defense proteins against antagonists, including ribotoxins. These toxins cleave a single phosphodiester bond within the universally conserved sarcin-ricin loop of ribosomes and inhibit protein biosynthesis. Here, we report on the structure and function of ageritin, a previously reported ribotoxin from the edible mushroom Agrocybe aegerita The amino acid sequence of ageritin was derived from cDNA isolated from the dikaryon A. aegerita AAE-3 and lacks, according to in silico prediction, a signal peptide for classical secretion, predicting a cytoplasmic localization of the protein. The calculated molecular weight of the protein is slightly higher than the one reported for native ageritin. The A. aegerita ageritin-encoding gene, AaeAGT1, is highly induced during fruiting, and toxicity assays with AaeAGT1 heterologously expressed in Escherichia coli showed a strong toxicity against Aedes aegypti larvae yet not against nematodes. The activity of recombinant A. aegerita ageritin toward rabbit ribosomes was confirmed in vitro Mutagenesis studies revealed a correlation between in vivo and in vitro activities, indicating that entomotoxicity is mediated by ribonucleolytic cleavage. The strong larvicidal activity of ageritin makes this protein a promising candidate for novel biopesticide development.IMPORTANCE Our results suggest a pronounced organismal specificity of a protein toxin with a very conserved intracellular molecular target. The molecular details of the toxin-target interaction will provide important insight into the mechanism of action of protein toxins and the ribosome. This insight might be exploited to develop novel bioinsecticides.

A new lipoxygenase from the agaric fungus Agrocybe aegerita: Biochemical characterization and kinetic properties

Oxylipins are metabolites with a variety of biological functions. However, the biosynthetic pathway is widely unknown. It is considered that the first step is the oxygenation of polyunsaturated fatty acids like linoleic acid. Therefore, a lipoxygenase (LOX) from the edible basidiomycete Agrocybe aegerita was investigated. The AaeLOX4 was heterologously expressed in E. coli and purified via affinity chromatography and gel filtration. Biochemical properties and kinetic parameters of the purified AaeLOX4 were determined with linoleic acid and linolenic acid as substrates. The obtained K_m, v_max and k_cat values for linoleic acid were 295.5 μM, 16.5 μM · min^-1 · mg^-1 and 103.9 s^-1, respectively. For linolenic acid K_m, v_max and k_cat values of 634.2 μM, 19.5 μM · min^-1 · mg^-1 and 18.3 s^-1 were calculated. Maximum activities were observed at pH 7.5 and 25 °C. The main product of linoleic acid conversion was identified with normal-phase HPLC. This analysis revealed an explicit production of 13-hydroperoxy-9,11-octadecadienoic acid (13-HPOD). The experimental regio specificity is underpinned by the amino acid residues W384, F450, R594 and V635 considered relevant for regio specificity in LOX. In conclusion, HPLC-analysis and alignments revealed that AaeLOX4 is a 13-LOX.

New terpenoids from the fermentation broth of the edible mushroom Cyclocybe aegerita

The strophariaceous basidiomycete Cyclocybe aegerita (synonyms Agrocybe aegerita and A. cylindracea) is one of the most praised cultivated edible mushrooms and is being cultivated at large scale for food production. Furthermore, the fungus serves as a model organism to study fruiting body formation and the production of secondary metabolites during the life cycle of Basidiomycota. By studying the secondary metabolite profiles of C. aegerita, we found several terpenoids in submerged cultures. Aside from the main metabolite, bovistol (1), two new bovistol derivatives B and C (2, 3) and pasteurestin C as a new protoilludane (4) were isolated by preparative HPLC. Their structures were elucidated by mass spectrometry and NMR spectroscopy. The relative configurations of 2–4 were assigned by ROESY correlations, and ³J_H,H coupling constants in the case of 4. Applying quantitative PCR for gene expression validation, we linked the production of bovistol and its derivatives to the respective biosynthesis gene clusters.

Genome and secretome of Chondrostereum purpureum correspond to saprotrophic and phytopathogenic life styles

The basidiomycete Chondrostereum purpureum (Silverleaf fungus) is a saprotroph and plant pathogen commercially used for combatting forest "weed" trees in vegetation management. However, little is known about its lignocellulose-degrading capabilities and the enzymatic machinery that is responsible for the degradative potential, and it is not yet clear to which group of wood-rot fungi it actually belongs. Here, we sequenced and analyzed the draft genome of C. purpureum (41.2 Mbp) and performed a quantitative proteomic approach during growth in submerged and solid-state cultures based on soybean meal suspension or containing beech wood supplemented with phenol-rich olive mill residues, respectively. The fungus harbors characteristic lignocellulolytic hydrolases (GH6 and GH7) and oxidoreductases (e.g. laccase, heme peroxidases). High abundance of some of these genes (e.g. 45 laccases, nine GH7) can be explained by gene expansion, e.g. identified for the laccase orthogroup ORTHOMCL11 that exhibits a total of 18 lineage-specific duplications. Other expanded genes families encode for proteins more related to a pathogenic lifestyle (e.g. protease and cytochrome P450s). The fungus responds to the presence of complex growth substrates (lignocellulose, phenolic residues) by the secretion of most of these lignocellulolytic and lignin-modifying enzymes (e.g. alcohol and aryl alcohol oxidases, laccases, GH6, GH7). Based on the genetic and enzymatic constitution, we consider the 'marasmioid' fungus C. purpureum as a 'phytopathogenic' white-rot fungus (WRF) that possesses a complex extracellular enzyme machinery to accomplish efficient lignocellulose degradation during both saprotrophic and phytopathogenic life phases.

Genome Sequencing Illustrates the Genetic Basis of the Pharmacological Properties of Gloeostereum incarnatum

Gloeostereum incarnatum is a precious edible mushroom that is widely grown in Asia and known for its useful medicinal properties. Here, we present a high-quality genome of G. incarnatum using the single-molecule real-time (SMRT) sequencing platform. The G. incarnatum genome, which is the first complete genome to be sequenced in the family Cyphellaceae, was 38.67 Mbp, with an N50 of 3.5 Mbp, encoding 15,251 proteins. Based on our phylogenetic analysis, the Cyphellaceae diverged ~174 million years ago. Several genes and gene clusters associated with lignocellulose degradation, secondary metabolites, and polysaccharide biosynthesis were identified in G. incarnatum, and compared with other medicinal mushrooms. In particular, we identified two terpenoid-associated gene clusters, each containing a gene encoding a sesterterpenoid synthase adjacent to a gene encoding a cytochrome P450 enzyme. These clusters might participate in the biosynthesis of incarnal, a known bioactive sesterterpenoid produced by G. incarnatum. Through a transcriptomic analysis comparing the G. incarnatum mycelium and fruiting body, we also demonstrated that the genes associated with terpenoid biosynthesis were generally upregulated in the mycelium, while those associated with polysaccharide biosynthesis were generally upregulated in the fruiting body. This study provides insights into the genetic basis of the medicinal properties of G. incarnatum, laying a framework for future characterization of bioactive proteins and pharmaceutical uses of this fungus.

Exploring molecular tools for transformation and gene expression in the cultivated edible mushroom Agrocybe aegerita

Agrocybe aegerita is a cultivated edible mushroom in numerous countries, which also serves as a model basidiomycete to study fruiting body formation. Aiming to create an easily expandable customised molecular toolset for transformation and constitutive gene of interest expression, we first created a homologous dominant marker for transformant selection. Progeny monokaryons of the genome-sequenced dikaryon A. aegerita AAE-3 used here were identified as sensitive to the systemic fungicide carboxin. We cloned the wild-type gene encoding the iron-sulphur protein subunit of succinate dehydrogenase AaeSdi1 including its up- and downstream regions, and introduced a single-point mutation (His237 to Leu) to make it confer carboxin resistance. PEG-mediated transformation of protoplasts derived from either oidia or vegetative monokaryotic mycelium with the resulting carboxin resistance marker (Cbx^R) plasmid pSDI1E3 yielded carboxin-resistant transformants in both cases. Plasmid DNA linearised within the selection marker resulted in transformants with ectopic multiple insertions of plasmid DNA in a head-to-tail repeat-like fashion. When circular plasmid was used, ectopic single integration into the fungal genome was favoured, but also gene conversion at the homologous locus was seen in 1 out of 11 analysed transformants. Employing Cbx^R as selection marker, two versions of a reporter gene construct were assembled via Golden Gate cloning which allows easy recombination of its modules. These consisted of an eGFP expression cassette controlled by the native promoter P_AaeGPDII and the heterologous terminator T_nos, once with and once without an intron in front of the eGFP start codon. After protoplast transformation with either construct as circular plasmid DNA, GFP fluorescence was detected with either transformants, indicating that expression of eGFP is intron-independent in A. aegerita. This paves the way for functional genetics approaches to A. aegerita, e.g., via constitutive expression of fruiting-related genes.

Genomic discovery of the hypsin gene and biosynthetic pathways for terpenoids in Hypsizygus marmoreus

Background

Hypsizygus marmoreus (Beech mushroom) is a popular ingredient in Asian cuisine. The medicinal effects of its bioactive compounds such as hypsin and hypsiziprenol have been reported, but the genetic basis or biosynthesis of these components is unknown.

Results

In this study, we sequenced a reference strain of H. marmoreus (Haemi 51,987–8). We evaluated various assembly strategies, and as a result the Allpaths and PBJelly produced the best assembly. The resulting genome was 42.7 Mbp in length and annotated with 16,627 gene models. A putative gene (Hypma_04324) encoding the antifungal and antiproliferative hypsin protein with 75% sequence identity with the previously known N-terminal sequence was identified. Carbohydrate active enzyme analysis displayed the typical feature of white-rot fungi where auxiliary activity and carbohydrate-binding modules were enriched. The genome annotation revealed four terpene synthase genes responsible for terpenoid biosynthesis. From the gene tree analysis, we identified that terpene synthase genes can be classified into six clades. Four terpene synthase genes of H. marmoreus belonged to four different groups that implies they may be involved in the synthesis of different structures of terpenes. A terpene synthase gene cluster was well-conserved in Agaricomycetes genomes, which contained known biosynthesis and regulatory genes.

Conclusions

Genome sequence analysis of this mushroom led to the discovery of the hypsin gene. Comparative genome analysis revealed the conserved gene cluster for terpenoid biosynthesis in the genome. These discoveries will further our understanding of the biosynthesis of medicinal bioactive molecules in this edible mushroom.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5159-y) contains supplementary material, which is available to authorized users.