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.

Fungi infestation as a disease has serious impacts on the cultivation of Morchella species. To investigate the effects of fungi infestation on the microbial diversity and community structure of soil when cultivating Morchella sextelata, we sampled soil samples of Morchella cultivars in the Qinghai–Xizang Platea and used metagenome sequencing technology to identify the disease fungi and analyze the differences in microbial diversity and structure between disease-infested and healthy soils. The disease fungi identified were Tricharina gilva and Peziza lohjaoensis, and the microbial diversity of T. gilva-infected soil was higher than that of healthy soil, while the diversity of P. lohjaoensis-infected soil was lower. Interestingly, whether infected with T. gilva or P. lohjaoensis, the soil microbial community was changed, and the dominant phyla and genera were different in different soil samples. When infected with P. lohjaoensis, the dominant phyla with relatively high abundances included Proteobacteria, Bacteroidetes, and Ascomycota, with average relative abundances of 44%, 18%, and 15%, respectively, and the dominant genera with high relative abundances encompassed Pseudomonadaceae, Terfezia, and Pedobacter, with average relative abundances of 8%, 9%, and 5%, respectively. Following infection with T. gilva, the dominant phyla with higher relative abundances were Proteobacteria, Acidobacteria, and Bacteroidetes, with average relative abundances of 46%, 15%, and 12%, respectively, and the dominant genera with high relative abundances included Hydrogenophaga, Sphingomonas, and Polaromonas, with average relative abundances of 9%, 3%, and 2%, respectively. Additionally, we found that lipid-metabolism-related genes were less abundant in the soil infected with P. lohjaoensis than in the other soil samples, and glycoside hydrolase diversity was lower in the soil infected with T. gilva than in other healthy soils. The results showed that the effects of different disease fungi on soil microbial communities and functional genes were different, which provided a theoretical basis for the sustainable cultivation of Morchella.

ATP-binding cassette (ABC) transporters constitute a superfamily of transmembrane transport proteins that are ubiquitously present across prokaryotic and eukaryotic organisms, exhibiting extensive substrate diversity and significant functional versatility. Yet, systematic genome-wide characterization of ABC transporters in Morchellasextelata remains largely unexplored. Here, a total of 33 ABC transporter genes were identified in M. sextelata and classified into eight subfamilies: ABCA (1), ABCB (7), ABCC (5), ABCD (6), ABCE (1), ABCF (4), ABCG (7), and ABCI (2), based on the features of domain composition and topology. These genes were unevenly distributed across 15 scaffolds, encoding proteins of 276–1626 aa (30.54–180.90 kDa) with pI values of 5.40–9.62 (16 acidic, 17 basic). Subcellular localization predictions indicated 24 MsABCs target the cell membrane, with others localizing to mitochondria, cytoplasm, or nucleus. Motif analysis identified 10 conserved motifs; all members harbor 1–2 copies of Motif 1, corresponding to the ATP-binding/Walker B region of the nucleotide-binding domain. Phylogenetically, ABC transporters form 7 major clades that align tightly with subfamily grouping. Promoter analysis revealed 67 distinct cis-element types grouped into 7 functional categories, with MYB (myeloblastosis) transcription factor binding motifs being the most frequent. Expression profiling showed 20 MsABCs were induced by 10 mg/L potassium selenite; during fruiting body morphogenesis, 19 MsABCs exhibited induced expression at late developmental stages, while nine genes maintained constitutively high expression across all stages. Thirteen MsABCs displayed higher transcript abundance in mycelium, four MsABCs were highly expressed in the stipe, and 14 MsABCs were highly expressed in the pileus. Collectively, our findings elucidate the genetic architecture of the M. sextelata ABC transporter family and delineate their expression dynamics under selenite supplementation, during fruiting body development, and across tissues, thereby establishing a framework for future functional characterization and utilization.