Morchella is a kind of precious edible, medicinal fungi with a series of important effects, including anti-tumor and anti-oxidation effects. Based on the data of 18 environmental variables and the distribution sites of wild Morchella species, this study used a maximum entropy (MaxEnt) model to predict the changes in the geographic distribution of Morchella species in different historical periods (the Last Glacial Maximum (LGM), Mid Holocene (MH), current, 2050s and 2070s). The results revealed that the area under the curve (AUC) values of the receiver operating characteristic curves of different periods were all relatively high (>0.83), indicating that the results of the maximum entropy model are good. Species distribution modeling showed that the major factors influencing the geographical distribution of Morchella species were the precipitation of the driest quarter (Bio17), elevation, the mean temperature of the coldest quarter (Bio11) and the annual mean temperature (Bio1). The simulation of geographic distribution suggested that the current suitable habitat of Morchella was mainly located in Yunnan, Sichuan, Gansu, Shaanxi, Xinjiang Uygur Autonomous Region (XUAR) and other provinces in China. Compared with current times, the suitable area in Northwest and Northeast China decreased in the LGM and MH periods. As for the future periods, the suitable habitats all increased under the different scenarios compared with those in contemporary times, showing a trend of expansion to Northeast and Northwest China. These results could provide a theoretical basis for the protection, rational exploitation and utilization of wild Morchella resources under scenarios of climate change.

True morels (Morchella) are a well-known edible fungi, with economically and medicinally important values. However, molecular identification and species taxonomy of the genus Morchella have long been controversial, due to numerous intermediate morphologies among species. In this study, we determined the identification efficiency of DNA barcoding and species classification of 260 individuals from 45 Morchella species, on the basis of multiple nuclear DNA markers. DNA barcoding analysis showed that the individual DNA fragment has a lower resolution of species identification than that of combined multiple DNA markers. ITS showed the highest level of species discrimination among the individual genetic markers. Interestingly, the combined DNA markers significantly increased the resolution of species identification. A combination of four DNA genes (EF1-α, RPB1, RPB2 and ITS) showed a higher species delimitation than that any combination of two or three markers. Phylogenetic analysis suggested that the species in genus Morchella could have been divided into two large genetic clades, the Elata Clade and Esculenta Clade lineages. The two lineages divided approximately 133.11 Mya [95% HPD interval: 82.77–197.95] in the early Cretaceous period. However, some phylogenetic species of Morchella showed inconsistent evolutionary relationships with the traditional morphological classifications, which may have resulted from incomplete lineage sorting and/or introgressive hybridization among species. These findings demonstrate that the interspecific gene introgression may have affected the species identification of true morels, and that the combined DNA markers significantly improve the resolution of species discrimination.

Morels (Morchella, Ascomycota) are an extremely desired group of edible mushrooms with worldwide distribution. Morchella eohespera is a typical black morel species, belonging to the Elata clade of Morchella species. The biological and genetic studies of this mushroom are rare, largely hindering the studies of molecular breeding and evolutionary aspects. In this study, we performed de novo sequencing and assembly of the M. eohespera strain m200 genome using the third-generation nanopore sequencing platform. The whole-genome size of M. eohespera was 53.81 Mb with a contig N50 of 1.93 Mb, and the GC content was 47.70%. A total of 9,189 protein-coding genes were annotated. Molecular dating showed that M. eohespera differentiated from its relative M. conica at ~19.03 Mya (million years ago) in Burdigalian. Evolutionary analysis showed that 657 gene families were contracted and 244 gene families expanded in M. eohespera versus the related morel species. The non-coding RNA prediction results showed that there were 336 tRNAs, 76 rRNAs, and 45 snRNAs in the M. eohespera genome. Interestingly, there was a high degree of repetition (20.93%) in the M. eohespera genome, and the sizes of long interspersed nuclear elements, short interspersed nuclear elements, and long terminal repeats were 0.83 Mb, 0.009 Mb, and 4.56 Mb, respectively. Additionally, selection pressure analysis identified that a total of 492 genes in the M. eohespera genome have undergone signatures of positive selection. The results of this study provide new insights into the genome evolution of M. eohespera and lay the foundation for in-depth research into the molecular biology of the genus Morchella in the future.

Morchella eohespera Beug, Voitk & O’Donnell is a typical black morel species. In this study, using the Nanopore sequencing platform, we characterized its whole mitochondrial (mt) genome sequence. Mt genome of M. eohespera is composed of circular DNA molecules of 243,963 bp, which encoded 102 protein-coding genes (PCGs), two ribosomal RNA genes (rRNA), and 31 transfer RNA (tRNA) genes. The base composition of M. eohespera mitogenome is as follows: A (30.40%), T (29.30%), G (20.8%), and C (19.5%). The phylogenetic analysis suggested that M. eohespera was closely related to the congeneric M. importuna.