Molecular Profiling of Rice Straw Degradability Discrepancy in Stropharia rugosoannulata Core Germplasm.

The rice-S. rugosoannulata pattern is a rapidly growing agricultural practice for straw disposal and mushroom production in China. However, different S. rugosoannulata strains show a large variation in rice straw degradability. Here, we constructed a core collection of S. rugosoannulata containing 14 strains with rich genetic diversity. The molecular profiling of the lignocellulose degradability discrepancy of S. rugosoannulata strains was then explored using enzyme activity assays and transcriptome analysis. The results indicated that mycelial growth rate, lignocellulolytic enzyme activities, and rice straw degradability differed widely among the S. rugosoannulata core strains. The genes encoding lignin modifying and degrading auxiliary enzymes, oxidases, glycoside hydrolases, and detoxification proteins were differentially expressed between two representative S. rugosoannulata strains, resulting in differences in their lignocellulolytic enzyme activities and further causing differences in lignocellulose degradability. This study is useful to improve the production efficiency of S. rugosoannulata and promote the recycling of rice straw.

Transcriptome Reveals the Key Genes Related to the Metabolism of Volatile Sulfur-Containing Compounds in Lentinula edodes Mycelium.

Lentinula edodes (L. edodes) is a globally popular edible mushroom because of its characteristic sulfur-containing flavor compounds. However, the formation of the volatile sulfur-containing compounds in the mycelium of L. edodes has not been studied. We found that there were also sulfur-containing aroma compounds in the mycelium of L. edodes, and the content and composition varied at different stages of mycelial growth and development. The γ-glutamyl-transpeptidase (GGT) and cysteine sulfoxide lyase (C-S lyase) related to the generation of sulfur compounds showed the highest activities in the 15-day sample. Candidate genes for the metabolism of volatile sulfur compounds in mycelium were screened using transcriptome analysis, including encoding the GGT enzyme, C-S lyase, fatty acid oxidase, HSP20, and P450 genes. The expression patterns of Leggt3 and Leccsl3 genes were consistent with the measured activities of GGT and C-S lyase during the cultivation of mycelium and molecular dynamics simulations showed that they could stably bind to the substrate. Our findings provide insights into the formation of sulfur-containing flavor compounds in L. edodes. The mycelium of L. edodes is suggested for use as material for the production of sulfur-containing flavor compounds.

Near-gapless genome and transcriptome analyses provide insights into fruiting body development in Lentinula edodes.

Fruiting body development in macrofungi is an intensive research subject. In this study, high-quality genomes were assembled for two sexually compatible monokaryons from a heterokaryotic Lentinula edodes strain WX1, and variations in L. edodes genomes were analyzed. Specifically, differential gene expression and allele-specific expression (ASE) were analyzed using the two monokaryotic genomes and transcriptome data from four different stages of fruiting body development in WX1. Results revealed that after aeration, mycelia sensed cell wall stress, pheromones, and a decrease in CO2 concentration, leading to up-regulated expression in genes related to cell adhesion, cell wall remodeling, proteolysis, and lipid metabolism, which may promote primordium differentiation. Aquaporin genes and those related to proteolysis, mitosis, lipid, and carbohydrate metabolism may play important roles in primordium development, while genes related to tissue differentiation and sexual reproduction were active in fruiting body. Several essential genes for fruiting body development were allele-specifically expressed and the two nuclear types could synergistically regulate fruiting body development by dominantly expressing genes with different functions. ASE was probably induced by long terminal repeat-retrotransposons. Findings here contribute to the further understanding of the mechanism of fruiting body development in macrofungi.

Comparative analysis of simulated in-situ colonization and degradation by Lentinula edodes on oak wafer and corn stalk.

Introduction: The depolymerization of lignocellulose biomass by white-rot fungi has been an important research topic. However, few simulated in-situ analyses have been conducted to uncover the decay. Methods: In this study, the white-rot Lentinula edodes was used to colonize the wood and non-wood substrates, and then hyphal transcriptional response and substrate degradation were analyzed during the spatial-temporal colonization on different type substrates to better understand the depolymerization of lignocellulose. Results and discussion: Faster growth and thicker mat of hyphae on corn stalk were observed in comparison to oak wafer. Coincide with the higher levels of gene transcripts related to protein synthesis on corn stalk. The higher lignin oxidase activity of hyphae was detected on oak wafer, and the higher cellulase activity was observed on corn stalk containing a much higher content of soluble sugars. A large number of carbohydrate-binding module (CBM1 and CBM20)-containing enzyme genes, including lytic polysaccharide monooxygenase (AA9), cellobiohydrolase (GH6 and GH7), glucanase (GH5), xylanase (GH10 and GH11), glucoamylase (GH15), and alpha-amylase (GH13), were significantly upregulated in the back-distal hyphae colonized on corn stalk. The hyphae tended to colonize and degrade the secondary cell wall, and the deposited oxalate crystal suggested that oxalate may play an important role during lignocellulose degradation. In addition, lignin was degraded in priority in oak wafer. Of note, three lignin monomers were degraded simultaneously in oak wafer but sequentially in corn stalk. This growth Our results indicated that the white-rot degradation pattern of lignocellulose is determined by the chemical composition and structure of the colonized biomass.

Prevalence and diversity of mycoviruses occurring in Chinese Lentinula edodes germplasm resource.

Incidence and banding patterns of virus-like dsRNA elements in 215 Chinese genetically diverse Lentinula edodes strains collected from wide geographic distribution (or producing areas) were first investigated, and 17 viruses were identified including eight novel viruses. The results revealed a 63.3% incidence of dsRNA elements in the cultivated strains and a 67.2% incidence in the wild strains. A total of 10 distinguishable dsRNAs ranging from 0.6 to 12 kbp and 12 different dsRNA patterns were detected in the positive strains. The molecular information of these dsRNA elements was characterized, and the molecular information of the other 12 different viral sequences with (+) ssRNA genome was revealed in four L. edodes strains with complex dsRNA banding patterns. RT-PCR was also done to verify the five dsRNA viruses and 12 (+) ssRNA ones. The results presented may enrich our understanding of L. edodes virus diversity, and will promote further research on virus-host interactions. IMPORTANCE: Viral infections involve complicated interactions including benign, harmful or possibly beneficial to hosts. Sometimes environment could lead to a transition in lifestyles from persistent to acute, resulting in a disease phenotype. The quality of spawn, such as the vulnerability to infection of viruses, is therefore important for mushroom production. Lentinula edodes, a wood rot basidiomycete fungus, was widely cultivated in the world for its edible and medicinal properties. In this study, the profile of dsRNA elements from Chinese genetically diverse L. edodes strains collected from wide geographic distribution or producing areas was first investigated. The molecular information of the dsRNA elements was characterized. Additionally, 12 different viral sequences with (+) ssRNA genome from four L. edodes strains with complex dsRNA banding patterns were identified. The results presented here will broaden our knowledge about mushroom viruses, and promote further studies of L. edodes production and the interaction between viruses and L. edodes.

Genome Sequencing Highlights the Plant Cell Wall Degrading Capacity of Edible Mushroom Stropharia rugosoannulata.

The basidiomycetous edible mushroom Stropharia rugosoannulata has excellent nutrition, medicine, bioremediation, and biocontrol properties. S. rugosoannulata has been widely and easily cultivated using agricultural by-products showing strong lignocellulose degradation capacity. However, the unavailable high-quality genome information has hindered the research on gene function and molecular breeding of S. rugosoannulata. This study provided a high-quality genome assembly and annotation from S. rugosoannulata monokaryotic strain QGU27 based on combined Illumina-Nanopore data. The genome size was about 47.97 Mb and consisted of 20 scaffolds, with an N50 of 3.73 Mb and a GC content of 47.9%. The repetitive sequences accounted for 17.41% of the genome, mostly long terminal repeats (LTRs). A total of 15,726 coding gene sequences were putatively identified with the BUSCO score of 98.7%. There are 142 genes encoding plant cell wall degrading enzymes (PCWDEs) in the genome, and 52, 39, 30, 11, 8, and 2 genes related to lignin, cellulose, hemicellulose, pectin, chitin, and cutin degradation, respectively. Comparative genomic analysis revealed that S. rugosoannulata is superior in utilizing aldehyde-containing lignins and is possible to utilize algae during the cultivation.

Construction of nucleus-directed fluorescent reporter systems and its application to verification of heterokaryon formation in Morchella importuna.

Introduction: Morchella importuna (M. importuna) is a rare fungus with high nutrition value and distinct flavor. Despite the successful artificial cultivation, its genetic characteristics and biological processes such as life cycle, reproductive system, and trophic mode remain poorly understood. Methods: Considering this, we constructed pEH2B and pMH2B vectors by fusing M. importuna endogenous histone protein H2B with fluorescent proteins eGFP or mCherry, respectively. Based on the constructed pEH2B and pMH2B vectors, nuclear fluorescence localization was performed via Agrobacterium tumefaciens-mediated transformation (ATMT). These two vectors were both driven by two endogenous promoters glyceraldehyde 3-phosphate dehydrogenase (GPD) and ubiquitin (UBI). The vector-based reporter systems were tested by the paired culture of two genetically modified strains pEH2B-labeled M04M24 (24e, MAT1-1-1) and pMH2B-abeled M04M26 (26m, MAT1-2-1). Results: The fluorescence observation and molecular identification results indicated the successful hyphal fusion and heterokaryon formation. We found that the expression of the reporter genes was stable, and it did not interfere with the growth of the fungus. Discussion: Our constructed nucleus-directed fluorescent systems in M. importuna can be used for monitoring the dynamic development and reproductive processes in living cells and also for monitoring the interaction between morels and plant roots. Therefore, morels exhibit the potential to be a candidate organism used for the research on basic biology and genetics of ascomycetes.

Comparative Transcriptome Analysis Provides Insights Into the Mechanism by Which 2,4-Dichlorophenoxyacetic Acid Improves Thermotolerance in Lentinula edodes.

As the widest cultivated edible mushroom worldwide, Lentinula edodes suffers serious yield and quality losses from heat stress during growth and development, and in our previous study, exogenous 2,4-Dichlorophenoxyacetic acid (2,4-D) was found to improve the thermotolerance of L. edodes strain YS3357, but the molecular mechanism remains unclear. Here, we explored the potential protective mechanism of exogenous 2,4-D against heat stress by transcriptome analysis. 2,4-D possible improve the thermotolerance of L. edodes through regulating antioxidant genes, transcription factors, energy-provision system, membrane fluidity, and cell wall remodeling. Furthermore, 2,4-D was also found to regulate the saturation levels of fatty acids and ATP content in L. edodes mycelium under heat stress. This study proposed a regulatory network of 2,4-D in regulating L. edodes response to heat stress, providing a theoretical basis for improving L. edodes thermotolerance, and facilitating the understanding of the molecular mechanism of exogenous hormones in alleviating abiotic stress damage to macrofungi.

Population genomics provides insights into the genetic basis of adaptive evolution in the mushroom-forming fungus Lentinula edodes.

Introduction: Mushroom-forming fungi comprise diverse species that develop complex multicellular structures. In cultivated species, both ecological adaptation and artificial selection have driven genome evolution. However, little is known about the connections among genotype, phenotype and adaptation in mushroom-forming fungi. Objectives: This study aimed to (1) uncover the population structure and demographic history of Lentinula edodes, (2) dissect the genetic basis of adaptive evolution in L. edodes, and (3) determine if genes related to fruiting body development are involved in adaptive evolution. Methods: We analyzed genomes and fruiting body-related traits (FBRTs) in 133 L. edodes strains and conducted RNA-seq analysis of fruiting body development in the YS69 strain. Combined methods of genomic scan for divergence, genome-wide association studies (GWAS), and RNA-seq were used to dissect the genetic basis of adaptive evolution. Results: We detected three distinct subgroups of L. edodes via single nucleotide polymorphisms, which showed robust phenotypic and temperature response differentiation and correlation with geographical distribution. Demographic history inference suggests that the subgroups diverged 36,871 generations ago. Moreover, L. edodes cultivars in China may have originated from the vicinity of Northeast China. A total of 942 genes were found to be related to genetic divergence by genomic scan, and 719 genes were identified to be candidates underlying FBRTs by GWAS. Integrating results of genomic scan and GWAS, 80 genes were detected to be related to phenotypic differentiation. A total of 364 genes related to fruiting body development were involved in genetic divergence and phenotypic differentiation. Conclusion: Adaptation to the local environment, especially temperature, triggered genetic divergence and phenotypic differentiation of L. edodes. A general model for genetic divergence and phenotypic differentiation during adaptive evolution in L. edodes, which involves in signal perception and transduction, transcriptional regulation, and fruiting body morphogenesis, was also integrated here.

Curing two predominant viruses occurring in Lentinula edodes by chemotherapy and mycelial fragmentation methods.

Previous research has established that Lentinula edodes mycovirus HKB (LeV-HKB) and L. edodes partitivirus 1(LePV1) are major mycoviruses identified in L.edodes germplasm. In this paper, two different methods for curing these two dsRNA mycoviruses, ribavirin treatment and mycelial fragmentation, were evaluated for the first time. Mycelial fragmentation was found to resulted in LeV-HKB- and LePV1-cured fungal strains, whereas ribavirin treatment could eliminate LeV-HKB only. Although no LePV1-cured strain was obtained via ribavirin treatment by the end of the experiment, the relative LePV1 concentration in the eighth successive subcultures was lower than that of the untreated control. The culture features of several virus-cured strains had faster mycelial growth rate and higher colony density than the infected ones. It was also suggested that LeV-HKB infection may affect the pigmentation in plate- and bag-cultivated mycelia of L. edodes strain L135.

Chromosome-Wide Characterization of Intragenic Crossover in Shiitake Mushroom, Lentinula edodes.

Meiotic crossover plays a critical role in generating genetic variations and is a central component of breeding. However, our understanding of crossover in mushroom-forming fungi is limited. Here, in Lentinula edodes, we characterized the chromosome-wide intragenic crossovers, by utilizing the single-nucleotide polymorphisms (SNPs) datasets of an F1 haploid progeny. A total of 884 intragenic crossovers were identified in 110 single-spore isolates, the majority of which were closer to transcript start sites. About 71.5% of the intragenic crossovers were clustered into 65 crossover hotspots. A 10 bp motif (GCTCTCGAAA) was significantly enriched in the hotspot regions. Crossover frequencies around mating-type A (MAT-A) loci were enhanced and formed a hotspot in L. edodes. Genome-wide quantitative trait loci (QTLs) mapping identified sixteen crossover-QTLs, contributing 8.5-29.1% of variations. Most of the detected crossover-QTLs were co-located with crossover hotspots. Both cis- and trans-QTLs contributed to the nonuniformity of crossover along chromosomes. On chr2, we identified a QTL hotspot that regulated local, global crossover variation and crossover hotspot in L. edodes. These findings and observations provide a comprehensive view of the crossover landscape in L. edodes, and advance our understandings of conservation and diversity of meiotic recombination in mushroom-forming fungi.

Enhanced Expression of Thaumatin-like Protein Gene (LeTLP1) Endows Resistance to Trichoderma atroviride in Lentinula edodes.

Lentinula edodes (shiitake mushrooms) is heavily affected by the infection of Trichoderma atroviride, causing yield loss and decreases quality in shiitake mushrooms. The selection and breeding of fungal-resistant L. edodes species are an important approach to protecting L. edodes from T. atroviride infection. Herein, a highly resistant L. edodes strain (Y3334) and a susceptible strain (Y55) were obtained by using a resistance evaluation test. Transcriptome analyses and qRT-PCR detection showed that the expression level of LeTLP1 (LE01Gene05009) was strongly induced in response to T. atroviride infection in the resistant Y3334. Then, LeTLP1-silenced and LeTLP1-overexpression transformants were obtained. Overexpression of LeTLP1 resulted in resistance to T. atroviride. Compared with the parent strain Y3334, LeTLP1-silenced transformants had reduced resistance relative to T. atroviride. Additionally, the LeTLP1 protein (Y3334) exhibited significant antifungal activity against T. atroviride. These findings suggest that overexpression of LeTLP1 is a major mechanism for the resistance of L. edodes to T. atroviride. The molecular basis provides a theoretical basis for the breeding of resistant L. edodes strains and can eventually contribute to the mushroom cultivation industry and human health.

RNA-Seq-based high-resolution linkage map reveals the genetic architecture of fruiting body development in shiitake mushroom, Lentinula edodes.

Fruiting body development (FBD) of mushroom-forming fungi has attracted tremendous interest. However, the genetic and molecular basis of FBD is poorly known. Here, using Lentinula edodes (shiitake) as a model, we deciphered the genetic architecture underlying fruiting body-related traits (FBRTs) by combined genomic, genetic and phenotypic data. Using RNA-Seq of fruiting bodies from 110 dikaryons in a bi-parental mapping population, we constructed an ultra-high-density genetic map of L. edodes (Lemap2.0) with a total length of 810.14 cM, which covered 81.7% of the shiitake genome. A total of 94 scaffolds of the shiitake genome were aligned to Lemap2.0 and re-anchored into nine pseudo-chromosomes. Then via quantitative trait locus (QTL) analysis, we disclosed an outline of the genetic architecture of FBD in shiitake. Twenty-nine QTLs and three main genomic regions associated with FBD of shiitake were identified. Using meta-QTL analysis, seven pleiotropic QTLs for multiple traits were detected, which contributed to the correlations of FBRTs. In the mapped QTLs, the expression of 246 genes were found to significantly correlate with the phenotypic traits. Thirty-three of them were involved in FBD and could represent candidate genes controlling the shape and size of fruiting bodies. Collectively, our findings have advanced our understanding of the genetic regulation of FBD in shiitake and mushroom-forming fungi at large.

Biological Characterization and Antagonist Screening of Cladosporium anthropophilum, a Novel Pathogen Causing Stipe Black Rot on Commercial Medicinal Mushroom, Flammulina filiformis (Agaricomycetes).

We report a new fungal disease, named stipe black rot, in a cultivation factory of Flammulina filiformis (an edible mushroom cultivated worldwide) in China. The pathogen was identified as Cladosporium anthropophilum by phylogenetic analysis and morphology characterization. C. anthropophilum was characterized to mainly infect the stipe bottom and cause stipe blackening and rot, with its optimal mycelial growth conditions consisting of 25°C, pH 7, and carbon and nitrogen sources of soluble starch and sodium nitrate, respectively. Furthermore, inhibitory evaluation showed that hydrogen peroxide silver disinfectant (HPSD) can efficiently inhibit the mycelial growth of C. anthropophilum, followed by the aqueous extracts of garlic and onion. This study identified C. anthropophilum as the pathogen for the new F. filiformis black rot disease and HPSD as an effective antagonist against the pathogen, which facilitates the understanding of fungal diseases and their control in edible mushrooms.

Mycoviral diversity and characteristics of a negative-stranded RNA virus LeNSRV1 in the edible mushroom Lentinula edodes.

Bioinformatics and RT-PCR analysis of RNA from four Lentinula edodes samples identified 22 different virus-like contigs comprising 15 novel and 3 previously reported viruses. We further investigated the Lentinula edodes negative-stranded RNA virus 1 (LeNSRV1) isolated from a symptomatic sample, whose virion is a filamentous particle with a diameter of ~15 nm and a length of ~1200 nm. RT-PCR analysis detected LeNSRV1 in 10 of the 56 Chinese L. edodes core collection strains and 6 of the 22 monokaryotic strains from the L. edodes strain HNZMD. Genetic variation analysis showed that the sequences encoding the nucleocapsid protein (ORF2) from all the aforementioned LeNSRV1 positive strains are very conservative. The results presented here may enrich our understanding of L. edodes virus diversity and the characteristics of LeNSRV1, and will promote further research on virus-host interaction in L. edodes.

The development of an efficient RNAi system based on Agrobacterium-mediated transformation approach for studying functional genomics in medical fungus Wolfiporia cocos.

Genetic transformation methods reported for Wolfiporia cocos are limited. In this study, we describe an efficient RNA interference (RNAi) system based on Agrobacterium-mediated transformation approach in W. cocos for the first time. Actively growing mycelial plugs were used as recipients for transformation using endogenous orotidine-5'-phosphate decarboxylase gene (URA3) as both a selective marker and a silencing gene, under the control of the dual promoters of Legpd and Leactin from Lentinula edodes and the single promoter of Wcgpd from W. cocos, respectively. The results showed that both the two kinds of promoters effectively drive the expression of URA3 gene, and the URA3-silenced transformants could be selected on CYM medium containing 5'-fluoroorotic acid. In addition, silencing URA3 gene has no effect on the growth of W. cocos hyphae. The incomplete silencing of the URA3 locus was also observed in this study. This study will promote further study on the mechanism of substrate degradation, sclerotial formation, and biosynthesis network of pharmacological compounds in W. cocos.

Characterization of Two Mitochondrial Genomes and Gene Expression Analysis Reveal Clues for Variations, Evolution, and Large-Sclerotium Formation in Medical Fungus Wolfiporia cocos.

Wolfiporia cocos, a precious mushroom with a long history as an edible food and Asian traditional medicine, remains unclear in the genetic mechanism underlying the formation of large sclerotia. Here, two complete circular mitogenomes (BL16, 135,686 bp and MD-104 SS10, 124,842 bp, respectively) were presented in detail first. The salient features in the mitogenomes of W. cocos include an intron in the tRNA (trnQ-UUG2), and an obvious gene rearrangement identified between the two mitogenomes from the widely geographically separated W. cocos strains. Genome comparison and phylogenetic analyses reveal some variations and evolutional characteristics in W. cocos. Whether the mitochondrion is functional in W. cocos sclerotium development was investigated by analyzing the mitogenome synteny of 10 sclerotium-forming fungi and mitochondrial gene expression patterns in different W. cocos sclerotium-developmental stages. Three common homologous genes identified across ten sclerotium-forming fungi were also found to exhibit significant differential expression levels during W. cocos sclerotium development. Most of the mitogenomic genes are not expressed in the mycelial stage but highly expressed in the sclerotium initial or developmental stage. These results indicate that some of mitochondrial genes may play a role in the development of sclerotium in W. cocos, which needs to be further elucidated in future studies. This study will stimulate new ideas on cytoplasmic inheritance of W. cocos and facilitate the research on the role of mitochondria in large sclerotium formation.

Hsp40 Protein LeDnaJ07 Enhances the Thermotolerance of Lentinula edodes and Regulates IAA Biosynthesis by Interacting LetrpE.

Our previous study found that LeDnaJ07 RNAi decreased Lentinula edodes resistance to heat stress and Trichoderma atroviride infection. In this study, the structure and function of the LeDnaJ07 gene was analyzed by gene cloning and overexpression in L. edodes stress-sensitive strain YS55 via the Agrobacterium-mediated transformation method. Transformants were confirmed by qRT-PCR, fluorescence observation and Southern blotting. Overexpression of LeDnaJ07 in YS55 not only enhanced L. edodes mycelial resistance to heat stress but also facilitated mycelial growth. In the presence of heat stress, the intracellular IAA content showed a significant increase in the two LeDnaJ07 overexpression strains but only a slight change in the YS55 wild type strain. Moreover, the interaction between LeDnaJ07 and LetrpE was demonstrated via Y2H and BiFC assays. These results suggested that LeDnaJ07 may be involved in regulating IAA biosynthesis and the resistance of L. edodes to heat stresses via interacting with LetrpE.

Bioconversion of rice straw agro-residues by Lentinula edodes and evaluation of non-volatile taste compounds in mushrooms.

Rice straw was substituted for sawdust at five different ratios of 0, 20%, 40%, 60%, and 80% (Control, RS20, RS40, RS60 and RS80, respectively) to obtain five kinds of Lentinula edodes. The effects of adding cropped rice straw to substrate formulas on the proximate composition and non-volatile taste compounds in mushrooms were investigated. The control group had the highest level of MY and BE among the five formulations. The protein levels in mushrooms decreased with the addition of rice straw and the ash levels increased. We found that trehalose, mannitol, and arabitol were the main soluble sugars in the five kinds of mushrooms. The contents of total free amino acids varied from 16.29 to 24.59 mg/g and the highest level of free amino acids was found in mushrooms cultivated from RS20 and RS40. Moreover, the addition of rice straw improved the contents of monosodium glutamate (MSG)-like amino acids in mushrooms. The 5'-Nucleotide levels ranged from 1.66 to 4.48 mg/g and equivalent umami concentration (EUC) value increased with the addition of rice straw. Our results suggest that rice straw is a potential substitute for sawdust to cultivate L. edodes with more non-volatile taste compounds.

Subchromosome-Scale Nuclear and Complete Mitochondrial Genome Characteristics of Morchella crassipes.

Morchella crassipes (Vent.) Pers., a typical yellow morel species with high economic value, is mainly distributed in the low altitude plains of Eurasia. However, rare research has been performed on its genomics and polarity, thus limiting its research and development. Here, we reported a fine physical map of the nuclear genome at the subchromosomal-scale and the complete mitochondrial genome of M. crassipes. The complete size of the nuclear genome was 56.7 Mb, and 23 scaffolds were assembled, with eight of them being complete chromosomes. A total of 11,565 encoding proteins were predicted. The divergence time analysis showed that M. crassipes representing yellow morels differentiated with black morels at ~33.98 Mya (million years), with 150 gene families contracted and expanded in M. crassipes versus the two black morels (M. snyderi and M. importuna). Furthermore, 409 CAZYme genes were annotated in M. crassipes, containing almost all plant cell wall degrading enzymes compared with the mycorrhizal fungi (truffles). Genomic annotation of mating type loci and amplification of the mating genes in the monospore population was conducted, the results indicated that M. crassipes is a heterothallic fungus. Additionally, a complete circular mitochondrial genome of M. crassipes was assembled, the size reached as large as 531,195 bp. It can be observed that the strikingly large size was the biggest up till now, coupled with 14 core conserved mitochondrial protein-coding genes, two rRNAs, 31 tRNAs, 51 introns, and 412 ncORFs. The total length of intron sequences accounted for 53.67% of the mitochondrial genome, with 19 introns having a length over 5 kb. Particularly, 221 of 412 ncORFs were distributed within 51 introns, and the total length of the ncORFs sequence accounted for 40.83% of the mitochondrial genome, and 297 ncORFs had expression activity in the mycelium stage, suggesting their potential functions in M. crassipes. Meanwhile, there was a high degree of repetition (51.31%) in the mitochondria of M. crassipes. Thus, the large number of introns, ncORFs and internal repeat sequences may contribute jointly to the largest fungal mitochondrial genome to date. The fine physical maps of nuclear genome and mitochondrial genome obtained in this study will open a new door for better understanding of the mysterious species of M. crassipes.