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

Near-complete de novo assembly of Tricholoma bakamatsutake chromosomes revealed the structural divergence and differentiation of Tricholoma genomes

Tricholoma bakamatsutake, which is an edible ectomycorrhizal fungus associated with Fagaceae trees, may have diverged before the other species in Tricholoma section Caligata. We generated a highly contiguous whole-genome sequence for T. bakamatsutake SF-Tf05 isolated in an Oak (Quercus salicina) forest in Japan. The assembly of high-fidelity long reads, with a median read length of 12.3 kb, resulted in 13 chromosome-sized contigs comprising 142,068,211 bases with an average guanine and cytosine (GC) content of 43.94%. The 13 chromosomes were predicted to encode 11,060 genes. A contig (122,566 bases) presumably containing the whole circular mitochondrial genome was also recovered. The chromosome-wide comparison of T. bakamatsutake and Tricholoma matsutake (TMA_r1.0) indicated that the basic number of chromosomes (13) was conserved, but the structures of the corresponding chromosomes diverged, with multiple inversions and translocations. Gene conservation and cluster analyses revealed at least 3 phylogenetic clades in Tricholoma section Caligata. Specifically, all T. bakamatsutake strains belonged to the "bakamatsutake" clade, which is most proximal to the "caligatum" clade consisting of Tricholoma caligatum and Tricholoma fulvocastaneum. The constructed highly contiguous nearly telomere-to-telomere genome sequence of a T. bakamatsutake isolate will serve as a fundamental resource for future research on the evolution and differentiation of Tricholoma species.

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.

Draft Genome Sequence for the Symbiotic Pine Mushroom Tricholoma matsutake

We report the high-quality genome sequence of Tricholoma matsutake strain 2001, which was isolated from a mushroom fruiting body in South Korea. The genome has 80 contigs, a size of 162.6 Mb, and an N₅₀ value of 5,103,859 bp and will provide insight into the symbiotic association between T. matsutake and Pinus densiflora.

Gene Expression Analyses of Mutant Flammulina velutipes (Enokitake Mushroom) with Clogging Phenomenon

Regulation of proper gene expression is important for cellular and organismal survival, maintenance, and growth. Abnormal gene expression, even for a single critical gene, can thwart cellular integrity and normal physiology to cause diseases, aging, and death. Therefore, gene expression profiling serves as a powerful tool to understand the pathology of diseases and to cure them. In this study, the difference in gene expression in Flammulina velutipes was compared between the wild type (WT) mushroom and the mutant one with clogging phenomenon. Differentially expressed transcripts were screened to identify the candidate genes responsible for the mutant phenotype using the DNA microarray analysis. A total of 88 genes including 60 upregulated and 28 downregulated genes were validated using the real-time quantitative PCR analysis. In addition, proteomic differences between the WT and mutant mushroom were analyzed using two-dimensional gel electrophoresis and matrix-assisted laser desorption/ionization-time of flight (MALDI-TOF). Interestingly, the genes identified by these genomic and proteomic analyses were involved in stress response, translation, and energy/sugar metabolism, including HSP70, elongation factor 2, and pyruvate kinase. Together, our data suggest that the aberrant expression of these genes attributes to the mutant clogging phenotype. We propose that these genes can be targeted to foster normal growth in F. velutipes.

Identification of up-regulated genes in Tricholoma matsutake mycorrhiza

Many plant roots associate with fungi to form mycorrhizae; tree roots especially associate with ectomycorrhizal fungi, such as Tricholoma species. Tricholoma matsutake is an economically important fungus in Asian countries and usually inhabits forests primarily composed of Pinus densiflora (Japanese red pine). In this study, to understand the mycorrhizal association between T. matsutake and P. densiflora, genes specifically expressed in mycorrhiza compared with those expressed in mycelia and fruiting bodies were identified by RNA-seq. This revealed that genes for chromatin, proteasomes, signal transduction, pheromones, cell surface receptors, cytoskeleton, RNA processing, and transporters from T. matsutake were highly expressed in mycorrhiza. It also identified 35 mycorrhiza-induced small secreted protein (MiSSPs) that were highly expressed in mycorrhiza. Meanwhile, genes for proteases, defence-related proteins, cell-wall degradation, signal transduction, pinene synthesis, plant hormones, and transporters from P. densiflora were highly expressed in mycorrhiza. These genes may be involved in mycorrhizal formation and maintenance. A MiSSP, 1 460 819, was highly expressed in mycorrhiza, and this expression was maintained for 24 months. These results provide insight into the mycorrhizal association between T. matsutake and P. densiflora.

Draft Genome Sequence of the Reference Strain of the Korean Medicinal Mushroom Wolfiporia cocos KMCC03342

Wolfiporia cocos is a wood-decay brown rot fungus belonging to the family Polyporaceae. While the fungus grows, the sclerotium body of the strain, dubbed Bokryeong in Korean, is formed around the roots of conifer trees. The dried sclerotium has been widely used as a key component of many medicinal recipes in East Asia. Wolfiporia cocos strain KMCC03342 is the reference strain registered and maintained by the Korea Seed and Variety Service for commercial uses. Here, we present the first draft genome sequence of W. cocos KMCC03342 using a hybrid assembly technique combining both short- and long-read sequences. The genome has a total length of 55.5 Mb comprised of 343 contigs with N50 of 332 kb and 95.8% BUSCO completeness. The GC ratio was 52.2%. We predicted 14,296 protein-coding gene models based on ab initio gene prediction and evidence-based annotation procedure using RNAseq data. The annotated genome was predicted to have 19 terpene biosynthesis gene clusters, which was the same number as the previously sequenced W. cocos strain MD-104 genome but higher than Chinese W. cocos strains. The genome sequence and the predicted gene clusters allow us to study biosynthetic pathways for the active ingredients of W. cocos.

Dynamic evolution of eukaryotic mitochondrial and nuclear genomes: a case study in the gourmet pine mushroom Tricholoma matsutake

Fungi, as eukaryotic organisms, contain two genomes, the mitochondrial genome and the nuclear genome, in their cells. How the two genomes evolve and correlate to each other is debated. Herein, taking the gourmet pine mushroom Tricholoma matsutake as an example, we performed comparative mitogenomic analysis using samples collected from diverse locations and compared the evolution of the two genomes. The T. matsutake mitogenome encodes 49 genes and is rich of repetitive and non-coding DNAs. Six genes were invaded by up to 11 group I introns, with one cox1 intron cox1P372 showing presence/absence dynamics among different samples. Bioinformatic analyses suggested limited or no evidence of mitochondrial heteroplasmy. Interestingly, hundreds of mitochondrial DNA fragments were found in the nuclear genome, with several larger than 500 nt confirmed by PCR assays and read count comparisons, indicating clear evidence of transfer of mitochondrial DNA into the nuclear genome. Nuclear DNA of T. matsutake showed a higher mutation rate than mitochondrial DNA. Furthermore, we found evidence of incongruence between phylogenetic trees derived from mitogenome and nuclear DNA sequences. Together, our results reveal the dynamic genome evolution of the gourmet pine mushroom.

Spore germination and ectomycorrhizae formation of Tricholoma matsutake on pine root systems with previously established ectomycorrhizae from a dikaryotic mycelial isolate of T. matsutake

In vitro ectomycorrhizal synthesis of Tricholoma matsutake with host plants has been widely conducted to elucidate fungal symbiotic properties for future cultivation practices. Here, we report on the importance of basidiospore inocula for this fungus to provide ectomycorrhizal seedlings in vitro. Ectomycorrhizal pine seedlings synthesized in vitro with cultured mycelium of T. matsutake (isolate #45 or #84) in a 250-mL culture vessel (soil volume) were transplanted to a large 1-L culture vessel. Fresh basidiospores of this fungus were aseptically inoculated on the ectomycorrhizal root system. The ectomycorrhizal seedlings in the 1-L vessel were grown for 9 months, and some plants were further grown for 6 more months under non-aseptic conditions in 4.1-L jars. The ectomycorrhizal seedlings previously inoculated with isolate #84 in the 1-L vessel showed significant ectomycorrhizal biomass (mycorrhizal root length) after spore inoculation. The ectomycorrhizal seedlings in the 4.1-L vessel showed large shiro structures (> 10 cm in diameter). PCR amplification of intergenic spacer 1 of the rRNA gene and long terminal repeat retroelement of T. matsutake in ectomycorrhizal root tips in both the 1-L vessels and 4.1-L jars revealed the presence of amplicons of the previously inoculated culture isolate of T. matsutake and the new genet(s) that established via germination of the inoculated basidiospores. This is the first report that inoculated basidiospores of T. matsutake germinated and colonized the host root to generate ectomycorrhizae in vitro.

Biosynthesis of (R)-(-)-1-Octen-3-ol in Recombinant Saccharomyces cerevisiae with Lipoxygenase-1 and Hydroperoxide Lyase Genes from Tricholoma matsutake

Tricholoma matsutake is an ectomycorrhizal fungus, related with the host of Pinus densiflora. Most of studies on T. matsutake have focused on mycelial growth, genes and genomics, phylogenetics, symbiosis, and immune activity of this strain. T. matsutake is known for its unique fragrance in Eastern Asia. The most major component of its scent is (R)-(-)-1-octen-3-ol and is biosynthesized from the substrate linoleic acid by the sequential reaction of lipoxygenase and peroxide lyase. Here, we report for the first time the biosynthesis of (R)-(-)- 1-octen-3-ol of T. matsutake using the yeast Saccharomyces cerevisiae as a host. In this study, cDNA genes correlated with these reactions were cloned from T. matsutake, and expression studies of theses genes were carried out in the yeast Saccharomyces cerevisiae. The product of these genes expression study was carried out with Western blotting. The biosynthesis of (R)-(-)- 1-octen-3-ol of T. matsutake in recombinant Saccharomyces cerevisiae was subsequently identified with GC-MS chromatography analysis. The biosynthesis of (R)-(-)-1-octen-3-ol with S. cerevisiae represents a significant step forward.

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

[This corrects the article DOI: 10.1371/journal.pone.0227923.].