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.

Morchella is a popular edible fungus worldwide due to its rich nutrition and unique flavor. Many research efforts were made on the domestication and cultivation of Morchella all over the world. In recent years, the cultivation of Morchella was successfully commercialized in China. However, the biology is not well understood, which restricts the further development of the morel fungus cultivation industry. In this paper, we performed de novo sequencing and assembly of the genomes of two monospores with a different mating type (M04M24 and M04M26) isolated from the commercially cultivated strain M04. Gene annotation and comparative genome analysis were performed to study differences in CAZyme (Carbohydrate-active enzyme) enzyme content, transcription factors, duplicated sequences, structure of mating type sites, and differences at the gene and functional levels between the two monospore strains of M. importuna. Results showed that the de novo assembled haploid M04M24 and M04M26 genomes were 48.98 and 51.07 Mb, respectively. A complete fine physical map of M. importuna was obtained from genome coverage and gene completeness evaluation. A total of 10,852 and 10,902 common genes and 667 and 868 endemic genes were identified from the two monospore strains, respectively. The Gene Ontology (GO) and KAAS (KEGG Automatic Annotation Serve) enrichment analyses showed that the endemic genes performed different functions. The two monospore strains had 99.22% collinearity with each other, accompanied with certain position and rearrangement events. Analysis of complete mating-type loci revealed that the two monospore M. importuna strains contained an independent mating-type structure and remained conserved in sequence and location. The phylogenetic and divergence time of M. importuna was analyzed at the whole-genome level for the first time. The bifurcation time of morel and tuber was estimated to be 201.14 million years ago (Mya); the two monospore strains with a different mating type represented the evolution of different nuclei, and the single copy homologous genes between them were also different due to a genetic differentiation distance about 0.65 Mya. Compared with truffles, M. importuna had an extension of 28 clusters of orthologous genes (COGs) and a contraction of two COGs. The two different polar nuclei with different degrees of contraction and expansion suggested that they might have undergone different evolutionary processes. The different mating-type structures, together with the functional clustering and enrichment analysis results of the endemic genes of the two different polar nuclei, imply that M. importuna might be a heterothallic fungus and the interaction between the endemic genes may be necessary for its complete life history. Studies on the genome of M. importuna facilitate a better understanding of morel biology and evolution.

The reliability of qRT-PCR results depend on the stability of reference genes used for normalization, suggesting the necessity of identification of reference genes before gene expression analysis. Morels are edible mushrooms well-known across the world and highly prized by many culinary kitchens. Here, several candidate genes were selected and designed according to the Morchella importuna transcriptome data. The stability of the candidate genes was evaluated with geNorm and NormFinder under three different experimental conditions, and several genes with excellent stability were selected. The extensive adaptability of the selected genes was tested in ten Morchella species. Results from the three experimental conditions revealed that ACT1 and INTF7 were the most prominent genes in Morchella, CYC3 was the most stable gene in different development stages, INTF4/AEF3 were the top-ranked genes across carbon sources, while INTF3/CYC3 pair showed the robust stability for temperature stress treatment. We suggest using ACT1, AEF3, CYC3, INTF3, INTF4 and INTF7 as reference genes for gene expression analysis studies for any of the 10 Morchella strains tested in this study. The stability and practicality of the gene, vacuolar protein sorting (INTF3), vacuolar ATP synthase (INTF4) and14-3-3 protein (INTF7) involving the basic biological processes were validated for the first time as the candidate reference genes for quantitative PCR. Furthermore, the stability of the reference genes was found to vary under the three different experimental conditions, indicating the importance of identifying specific reference genes for particular conditions.