Comparative transcriptomics reveals potential genes involved in the vegetative growth of Morchella importuna

Associated bacterial communities, confrontation studies, and comparative genomics reveal important interactions between Morchella with Pseudomonas spp

Members of the fungal genus Morchella are widely known for their important ecological roles and significant economic value. In this study, we used amplicon and genome sequencing to characterize bacterial communities associated with sexual fruiting bodies from wild specimens, as well as vegetative mycelium and sclerotia obtained from Morchella isolates grown in vitro. These investigations included diverse representatives from both Elata and Esculenta Morchella clades. Unique bacterial community compositions were observed across the various structures examined, both within and across individual Morchella isolates or specimens. However, specific bacterial taxa were frequently detected in association with certain structures, providing support for an associated core bacterial community. Bacteria from the genus Pseudomonas and Ralstonia constituted the core bacterial associates of Morchella mycelia and sclerotia, while other genera (e.g., Pedobacter spp., Deviosa spp., and Bradyrhizobium spp.) constituted the core bacterial community of fruiting bodies. Furthermore, the importance of Pseudomonas as a key member of the bacteriome was supported by the isolation of several Pseudomonas strains from mycelia during in vitro cultivation. Four of the six mycelial-derived Pseudomonas isolates shared 16S rDNA sequence identity with amplicon sequences recovered directly from the examined fungal structures. Distinct interaction phenotypes (antagonistic or neutral) were observed in confrontation assays between these bacteria and various Morchella isolates. Genome sequences obtained from these Pseudomonas isolates revealed intriguing differences in gene content and annotated functions, specifically with respect to toxin-antitoxin systems, cell adhesion, chitinases, and insecticidal toxins. These genetic differences correlated with the interaction phenotypes. This study provides evidence that Pseudomonas spp. are frequently associated with Morchella and these associations may greatly impact fungal physiology.

Integration of Metabolomes and Transcriptomes Provides Insights into Morphogenesis and Maturation in Morchella sextelata

True morels (Morchella, Pezizales) are a popular edible and medicinal fungus with great nutritional and economic value. The dynamics and regulatory mechanisms during the morphogenesis and maturation of morels are poorly understood. In this study, the metabolomes and transcriptomes of the mycelium (MY), primordium differentiation (PR), young fruiting body (YFB), and mature fruiting body (MFB) were comprehensively analyzed to reveal the mechanism of the morphogenesis and maturation of Morchella sextelata. A total of 748 differentially expressed metabolites (DEMs) and 5342 differentially expressed genes (DEGs) were detected, mainly enriched in the carbohydrate, amino acid, and lipid metabolism pathways, with the transition from the mycelium to the primordium being the most drastic stage at both the metabolic and transcriptional levels. The integrated metabolomics and transcriptomics highlighted significant correlations between the DEMs and DEGs, and specific amino acid and nucleotide metabolic pathways were significantly co-enriched, which may play key roles in morphological development and ascocarp maturation. A conceptual model of transcriptional and metabolic regulation was proposed during morphogenesis and maturation in M. sextelata for the first time, in which environmental factors activate the regulation of transcription factors, which then promote metabolic and transcriptional regulation from vegetative to reproductive growth. These results provide insights into the metabolic dynamics and transcriptional regulation during the morphogenesis and maturation of morels and valuable resources for future breeding enhancement and sustainable artificial cultivation.

Large-Scale Field Cultivation of Morchella and Relevance of Basic Knowledge for Its Steady Production

Morels are one of the most highly prized edible and medicinal mushrooms worldwide. Therefore, historically, there has been a large international interest in their cultivation. Numerous ecological, physiological, genetic, taxonomic, and mycochemical studies have been previously developed. At the beginning of this century, China finally achieved artificial cultivation and started a high-scale commercial development in 2012. Due to its international interest, its cultivation scale and area expanded rapidly in this country. However, along with the massive industrial scale, a number of challenges, including the maintenance of steady economic profits, arise. In order to contribute to the solution of these challenges, formal research studying selection, species recognition, strain aging, mating type structure, life cycle, nutrient metabolism, growth and development, and multi-omics has recently been boosted. This paper focuses on discussing current morel cultivation technologies, the industrial status of cultivation in China, and the relevance of basic biological research, including, e.g., the study of strain characteristics, species breeding, mating type structure, and microbial interactions. The main challenges related to the morel cultivation industry on a large scale are also analyzed. It is expected that this review will promote a steady global development of the morel industry based on permanent and robust basic scientific knowledge.

Transcriptomic and Metabolomic Profiles Provide Insights into the Red-Stipe Symptom of Morel Fruiting Bodies

The cultivation of true morels (Morchella spp., Morchellaceae, Ascomycota) has rapidly expanded in recent years, especially in China. Red stipe is a symptom wherein the stipe of morel fruiting bodies becomes red-gray, resulting in the gradual death of the affected fruiting bodies. The impact of red-stipe symptom occurrence on the development and nutritional quality of morel fruiting bodies remains unclear. Herein, morel ascocarps with the red-stipe symptom (R) and normal (N), artificially cultivated in the Fujian Province of China, were selected for the transcriptome and metabolome analysis to study the physiological and biochemical responses of morel fruiting bodies to the red-stipe symptom. Transcriptome data revealed several differentially expressed genes between the R and N groups significantly enriched in the tyrosine, riboflavin, and glycerophospholipid metabolism pathways. Similarly, the differentially accumulated metabolites were mainly assigned to metabolic pathways, including tyrosine, the biosynthesis of plant secondary metabolites, and the biosynthesis of amino acids. Moreover, the transcriptome and metabolome data combination revealed that tyrosine metabolism was the most enriched pathway, which was followed by ATP-binding cassette (ABC) transport, alanine, aspartate, and glutamate metabolism. Overall, the integration of transcriptomic and metabolomic data of M. sextelata affected by red-stipe symptoms identified several important genes, metabolites, and pathways. These findings further improve our understanding of the mechanisms underlying the red-stipe symptom development of M. sextelata and provide new insights into how to optimize its cultivation methods.

Transcriptomics combined with metabolomics unveiled the key genes and metabolites of mycelium growth in Morchella importuna

Morels (Morchella) are one of the most popular edible fungi in the world, especially known for their rich nutrition and delicious taste. Earlier research indicates that the production of fruiting bodies can be affected by the growth of mycelium. To investigate the molecular mechanisms underlying mycelium growth in Morchella importuna, we performed transcriptome analysis and metabolomics analysis of three growth stages of the hypha of M. importuna. As a result, 24 differentially expressed genes, such as transketolase (tktA), glucose-6-phosphate dehydrogenase (G6PDH), fructose-diphosphate aldolase (Fba), and ribose-5-phosphate isomerase (rpiA), as well as 15 differentially accumulated metabolites, including succinate and oxaloacetate, were identified and considered as the key genes and metabolites to mycelium growth in M. importuna. In addition, guanosine 3',5'-cyclic monophosphate (cGMP), guanosine-5'-monophosphate (GMP), and several small peptides were found to differentially accumulate in different growth stages. Furthermore, five pathways, namely, starch and sucrose metabolism, pentose and glucuronate interconversions, fructose and mannose metabolism, tyrosine metabolism, and purine nucleotides, enriched by most DEGs, existed in the three compared groups and were also recognized as important pathways for the development of mycelium in morels. The comprehensive transcriptomics and metabolomics data generated in our study provided valuable information for understanding the mycelium growth of M. importuna, and these data also unveiled the key genes, metabolites, and pathways involved in mycelium growth. This research provides a great theoretical basis for the stable production and breeding of morels.

Out of the Qinghai-Tibetan plateau: Origin, evolution and historical biogeography of Morchella (both Elata and Esculenta clades)

Introduction

Morchella has become a research hotspot because of its wide distribution, delicious taste, and phenotypic plasticity. The Qinghai-Tibet Plateau subkingdoms (QTPs) are known as the cradle of Ice age biodiversity. However, the diversity of Morchella in the QTPs has been poorly investigated, especially in phylogenetic diversity, origin, and biogeography.

Methods

The genealogical concordance phylogenetic species recognition (GCPSR, based on Bayesian evolutionary analysis using sequences from the internal transcribed spacer (ITS), nuclear large subunit rDNA (nrLSU), translation elongation factor 1-α (EF1-α), and the largest and second largest subunits of RNA polymerase II (RPB1 and RPB2)), differentiation time estimation, and ancestral region reconstruction were used to infer Morchella's phylogenetic relationships and historical biogeography in the QTPs.

Results

Firstly, a total of 18 Morchella phylogenetic species are recognized in the QTPs, including 10 Elata clades and 8 Esculenta clades of 216 individuals Secondly, the divergences of the 18 phylogenetic species were 50.24-4.20 Mya (Eocene-Pliocene), which was closely related to the geological activities in the QTPs. Furthermore, the ancestor of Morchella probably originated in the Northern regions (Qilian Shan, Elata cade) and southwestern regions (Shangri-La, Esculenta clade) of QTPs and might have migrated from North America (Rufobrunnea clade) via Beringian Land Bridge (BLB) and Long-Distance Dispersal (LDD) expansions during the Late Cretaceous. Moreover, as the cradle of species origin and diversity, the fungi species in the QTPs have spread out and diffused to Eurasia and South Africa starting in the Paleogene Period.

Conclusion

This is the first report that Esculenta and Elata clade of Morchella originated from the QTPs because of orogenic, and rapid differentiation of fungi is strongly linked to geological uplift movement and refuge in marginal areas of the QTPs. Our findings contribute to increasing the diversity of Morchella and offer more evidence for the origin theory of the QTPs.

Epidemic Identification of Fungal Diseases in Morchella Cultivation across China

True morels (Morchella, Pezizales) are world-renowned edible mushrooms (ascocarps) that are widely demanded in international markets. Morchella has been successfully artificially cultivated since 2012 in China and is rapidly becoming a new edible mushroom industry occupying up to 16,466 hectares in the 2021-2022 season. However, nearly 25% of the total cultivation area has annually suffered from fungal diseases. While a variety of morel pathogenic fungi have been reported their epidemic characteristics are unknown, particularly in regional or national scales. In this paper, ITS amplicon sequencing and microscopic examination were concurrently performed on the morel ascocarp lesions from 32 sites in 18 provinces across China. Results showed that Diploöspora longispora (75.48%), Clonostachys solani (5.04%), Mortierella gamsii (0.83%), Mortierella amoeboidea (0.37%) and Penicillium kongii (0.15%) were the putative pathogenic fungi. The long, oval, septate conidia of D. longispora was observed on all ascocarps. Oval asexual spores and sporogenic structures, such as those of Clonostachys, were also detected in C. solani infected samples with high ITS read abundance. Seven isolates of D. longispora were isolated from seven selected ascocarps lesions. The microscopic characteristics of pure cultures of these isolates were consistent with the morphological characteristics of ascocarps lesions. Diploöspora longispora had the highest amplification abundance in 93.75% of the samples, while C. solani had the highest amplification abundance in six biological samples (6.25%) of the remaining two sampling sites. The results demonstrate that D. longispora is a major culprit of morel fungal diseases. Other low-abundance non-host fungi appear to be saprophytic fungi infecting after D. longispora. This study provides data supporting the morphological and molecular identification and prevention of fungal diseases of morel ascocarps.

Genome and Comparative Transcriptome Dissection Provide Insights Into Molecular Mechanisms of Sclerotium Formation in Culinary-Medicinal Mushroom Pleurotus tuber-regium

Pleurotus tuber-regium is an edible and medicinal sclerotium-producing mushroom. The sclerotia of this mushroom also serve as food and folk medicine. Based on the description of its monokaryon genome, sequenced with Illumina and PacBio sequencing technologies, comparative transcriptomic analysis using RNA sequencing (RNA-seq) was employed to study its mechanism of sclerotium formation. The de novo assembled genome is 35.82 Mb in size with a N50 scaffold size of 4.29 Mb and encodes 12,173 putative proteins. Expression analysis demonstrated that 1,146 and 1,249 genes were upregulated and downregulated with the formation of sclerotia, respectively. The differentially expressed genes were associated with substrate decomposition, the oxidation-reduction process, cell wall synthesis, and other biological processes in P. tuber-regium. These genomic and transcriptomic resources provide useful information for the mechanism underlying sclerotium formation in P. tuber-regium.

Transcriptome analysis of genes associated with autolysis of Coprinus comatus

Coprinus comatus, widely known as "Jituigu", is an important commodity and food in China. The yield of C. comatus, however, is substantially reduced by the autolysis of the fruiting bodies after harvest. To gain insight into the molecular mechanism underlying this autolysis, we divided the growth of C. comatus fruiting bodies into four stages: infant stage (I), mature stage (M), discolored stage (D), and autolysis stage (A). We then subjected these stages to de novo transcriptomic analysis using high-throughput Illumina sequencing. A total of 12,946 unigenes were annotated and analyzed with the Gene Ontology (GO), Clusters of Orthologous Groups of proteins (COG), and Kyoto Encyclopedia of Genes and Genomes (KEGG). We analyzed the differentially expressed genes (DEGs) between stages I and M, M and D, and D and A. Because the changes from M to D are thought to be related to autolysis, we focused on the DEGs between these two stages. We found that the pathways related to metabolic activity began to vary in the transition from M to D, including pathways named as autophagy-yeast, peroxisome, and starch and sucrose metabolism. This study also speculates the possible process of the autolysis of Coprinus comatus. In addition, 20 genes of interest were analyzed by quantitative real-time PCR to verify their expression profiles at the four developmental stages. This study, which is the first to describe the transcriptome of C. comatus, provides a foundation for future studies concerning the molecular basis of the autolysis of its fruiting bodies.

Whole-Genome and Transcriptome Sequencing of Phlebopus portentosus Reveals Its Associated Ectomycorrhizal Niche and Conserved Pathways Involved in Fruiting Body Development

Phlebopus portentosus (Berk. and Broome) Boedijin, a widely consumed mushroom in China and Thailand, is the first species in the order Boletaceae to have been industrially cultivated on a large scale. However, to date, the lignocellulose degradation system and molecular basis of fruiting body development in P. portentosus have remained cryptic. In the present study, genome and transcriptome sequencing of P. portentosus was performed during the mycelium (S), primordium (P), and fruiting body (F) stages. A genome of 32.74 Mb with a 48.92% GC content across 62 scaffolds was obtained. A total of 9,464 putative genes were predicted from the genome, of which the number of genes related to plant cell wall-degrading enzymes was much lower than that of some saprophytic mushrooms with specific ectomycorrhizal niches. Principal component analysis of RNA-Seq data revealed that the gene expression profiles at all three stages were different. The low expression of plant cell wall-degrading genes also confirmed the limited ability to degrade lignocellulose. The expression profiles also revealed that some conserved and specific pathways were enriched in the different developmental stages of P. portentosus. Starch and sucrose metabolic pathways were enriched in the mycelium stage, while DNA replication, the proteasome and MAPK signaling pathways may be associated with maturation. These results provide a new perspective for understanding the key pathways and hub genes involved in P. portentosus development.

Comparative Transcriptome and Endophytic Bacterial Community Analysis of Morchella conica SH

The precious rare edible fungus Morchella conica is popular worldwide for its rich nutrition, savory flavor, and varieties of bioactive components. Due to its high commercial, nutritional, and medicinal value, it has always been a hot spot. However, the molecular mechanism and endophytic bacterial communities in M. conica were poorly understood. In this study, we sequenced, assembled, and analyzed the genome of M. conica SH. Transcriptome analysis reveals significant differences between the mycelia and fruiting body. As shown in this study, 1,329 and 2,796 genes were specifically expressed in the mycelia and fruiting body, respectively. The Gene Ontology (GO) enrichment showed that RNA polymerase II transcription activity-related genes were enriched in the mycelium-specific gene cluster, and nucleotide binding-related genes were enriched in the fruiting body-specific gene cluster. Further analysis of differentially expressed genes in different development stages resulted in finding two groups with distinct expression patterns. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment displays that glycan degradation and ABC transporters were enriched in the group 1 with low expressed level in the mycelia, while taurine and hypotaurine metabolismand tyrosine metabolism-related genes were significantly enriched in the group 2 with high expressed level in mycelia. Moreover, a dynamic shift of bacterial communities in the developing fruiting body was detected by 16S rRNA sequencing, and co-expression analysis suggested that bacterial communities might play an important role in regulating gene expression. Taken together, our study provided a better understanding of the molecular biology of M. conica SH and direction for future research on artificial cultivation.

Comparative transcriptome analysis of cells from different areas reveals ROS responsive mechanism at sclerotial initiation stage in Morchella importuna

Morels are some of the most highly prized edible and medicinal mushrooms, with great economic and scientific value. Outdoor cultivation has been achieved and expanded on a large scale in China in recent years. Sclerotial formation is one of the most important phases during the morel life cycle, and previous reports indicated that reactive oxygen species (ROS) play an important role. However, ROS response mechanisms at sclerotial initiation (SI) stage are poorly understood. In this study, comparative transcriptome analyses were performed with sclerotial and hyphal cells at different areas in the same plate at SI stage. Gene expression was significantly different at SI stage between sclerotial formation and mycelia growth areas. GO and KEGG analyses indicated more vigorous metabolic characteristics in the hyphae area, while transcription process, DNA repair, and protein processing were enriched in sclerotial cells. Gene expression related to H₂O₂ production was high in the hyphae area, while expression of H₂O₂-scavenging genes was high in sclerotial cells, leading to a higher H₂O₂ concentration in the hyphal region than in the sclerotium. Minor differences were observed in gene expression of H₂O₂-induced signaling pathway in sclerotial and hyphal cells; however, expression levels of the target genes of transcription factor MSN2, important in the H₂O₂-induced signaling pathways, were significantly different. MSN2 enhanced stress response regulation in sclerotia by regulating these target genes. Small molecular HSPs were also found upregulated in sclerotial cells. This study indicated that sclerotial cells are more resistant to ROS stress than hyphal cells through transcriptional regulation of related genes.

Untargeted Metabolomics Reveals Alterations in the Primary Metabolites and Potential Pathways in the Vegetative Growth of Morchella sextelata

Morchella sextelata, one of the true morels, has recently been artificially cultivated with stable production in China. Analysis of the variations in primary metabolites during the vegetative stages of M. sextelata is critical for understanding the metabolic process. In this study, three developmental stages were categorized based on morphological and developmental changes, including the young mushroom period, marketable mature period, and physiological maturity period. Untargeted metabolomics-based mass spectrometry was used to analyze the change of metabolites during the growth stages of M. sextelata. The result showed that the metabolites’ content at the different growth stages were significantly different. The relative contents of linoleic acid, mannitol, oleamide, and betaine were higher at each growth stage. Flavor substances were significantly metabolizable during commodity maturity, while amino acids, organic acids, and lipids were significantly metabolizing at physiological maturity. Pathway analysis of the most significant changes involved Pyrimidine metabolism, Vitamin B6 metabolism, Arginine biosynthesis, Lysine biosynthesis, and Lysine degradation. The results can provide a theoretical basis for further clarifying the metabolic regulation mechanism and lay the foundation for optimizing the cultivation process of M. sextelata.

A comprehensive review on Morchella importuna: cultivation aspects, phytochemistry, and other significant applications

Morchella importuna is one of the most highly priced edible mushrooms and a rich source of bioactive substances with numerous beneficial medicinal properties. It has been artificially cultivated in the last few years but due to the unclear mechanism of its fruiting body formation, the stable production has not been achieved yet. This review aims to provide the detailed information about the current research status of the indoor and outdoor cultivation and significant applications of M. importuna worldwide with a vision to know more about its potential therapeutic index and edible values. It will provide the basis of better understanding about the need to develop biotechnological processes for morel farming under controlled conditions.

The mitochondrial genome of Morchella importuna (272.2 kb) is the largest among fungi and contains numerous introns, mitochondrial non-conserved open reading frames and repetitive sequences

The complete mitochondrial genome of Morchella importuna, the famous edible and medicinal mushroom, was assembled as a 272,238 bp single circular dsDNA. As the largest mitogenome among fungi, it exhibits several distinct characteristics. The mitogenome of M. importuna encoded 14 core conserved mitochondrial protein-coding genes and 151 mitochondrial non-conserved open reading frames (ncORFs) were predicted, of which 61 were annotated as homing endonuclease genes, and 108 were confirmed to be expressed during the vegetative growth stages of M. importuna. In addition, 34 introns were identified in seven core genes (cob, cox1, cox2, cox3, nad1, nad4 and nad5) and two rRNA genes (rrnS and rrnL) with a length from 383 bp to 7453 bp, and eight large introns with a length range of 2340 bp to 7453 bp contained multiple intronic mtORFs. Moreover, 34 group I (IA, IB, IC1, IC2, ID and derived group I introns) and four group II intron domains were identified for the 34 introns, including five hybrid ones. Furthermore, the M. importuna mitogenome showed the presence of about 18.7% mitogenomic interspersed repeats. These and the aforementioned ncORFs and introns, contributed to the enlarged size of the mitogenome.

Comparative transcriptome analysis reveals potential fruiting body formation mechanisms in Morchella importuna

Morchella importuna has been artificially cultivated, but stable production remains difficult because its mechanisms of fruiting body formation are unclear. To investigate the fruiting body formation mechanisms, we sequenced the transcriptomes of Morchella importuna at the mycelial and young fruiting body stages. Among the 12,561 differentially expressed genes (DEGs), 9215 were upregulated, and 3346 were downregulated. DEG enrichment analysis showed that these genes were enriched in the "generation of precursor metabolites and energy", "carbohydrate catabolic process", and "oxidoreductase activity" Gene Ontology (GO) functional categories. Enzyme activity assay results indicated that the activity levels of CAZymes (carbohydrate-active enzymes), oxidoreductases (SOD (superoxide dismutase), CAT (catalase)) and mitochondrial complex (complex I, II, III) proteins were significantly increased from the mycelial stage to the young fruiting body stage. In addition, the genes encoding CAZymes, mitochondrial proteins, oxidoreductases and heat shock proteins had higher expression levels in the young fruiting body stage than in the mycelial stage, and the qRT-PCR results showed similar trends to the RNA-Seq results. In summary, these results suggest that carbohydrate catabolism and energy metabolism are significantly enhanced in the young fruiting body stage and that growth environment temperature changes affect the formation of fruiting bodies.