Comparative transcriptome analysis revealed genes involved in the fruiting body development of Ophiocordyceps sinensis

Comparative transcriptome analysis on candidate genes associated with fruiting body growth and development in Lyophyllum decastes

Lyophyllum decastes is a mushroom that is highly regarded for its culinary and medicinal properties. Its delectable taste and texture make it a popular choice for consumption. To gain a deeper understanding of the molecular mechanisms involved in the development of the fruiting body of L. decastes, we used RNA sequencing to conduct a comparative transcriptome analysis. The analysis encompassed various developmental stages, including the vegetative mycelium, primordial initiation, young fruiting body, medium-size fruiting body, and mature fruiting body stages. A range of 40.1 to 60.6 million clean reads were obtained, and de novo assembly generated 15,451 unigenes with an average length of 1,462.68 bp. Functional annotation of transcriptomes matched 76.84% of the unigenes to known proteins available in at least one database. The gene expression analysis revealed a significant number of differentially expressed genes (DEGs) between each stage. These genes were annotated and subjected to Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses. Highly differentially expressed unigenes were also identified, including those that encode extracellular enzymes, transcription factors, and signaling pathways. The accuracy of the RNA-Seq and DEG analyses was validated using quantitative PCR. Enzyme activity analysis experiments demonstrated that the extracellular enzymes exhibited significant differences across different developmental stages. This study provides valuable insights into the molecular mechanisms that underlie the development of the fruiting body in L. decastes.

Lessons on fruiting body morphogenesis from genomes and transcriptomes of Agaricomycetes

Fruiting bodies (sporocarps, sporophores or basidiomata) of mushroom-forming fungi (Agaricomycetes) are among the most complex structures produced by fungi. Unlike vegetative hyphae, fruiting bodies grow determinately and follow a genetically encoded developmental program that orchestrates their growth, tissue differentiation and sexual sporulation. In spite of more than a century of research, our understanding of the molecular details of fruiting body morphogenesis is still limited and a general synthesis on the genetics of this complex process is lacking. In this paper, we aim at a comprehensive identification of conserved genes related to fruiting body morphogenesis and distil novel functional hypotheses for functionally poorly characterised ones. As a result of this analysis, we report 921 conserved developmentally expressed gene families, only a few dozens of which have previously been reported to be involved in fruiting body development. Based on literature data, conserved expression patterns and functional annotations, we provide hypotheses on the potential role of these gene families in fruiting body development, yielding the most complete description of molecular processes in fruiting body morphogenesis to date. We discuss genes related to the initiation of fruiting, differentiation, growth, cell surface and cell wall, defence, transcriptional regulation as well as signal transduction. Based on these data we derive a general model of fruiting body development, which includes an early, proliferative phase that is mostly concerned with laying out the mushroom body plan (via cell division and differentiation), and a second phase of growth via cell expansion as well as meiotic events and sporulation. Altogether, our discussions cover 1 480 genes of Coprinopsis cinerea, and their orthologs in Agaricus bisporus, Cyclocybe aegerita, Armillaria ostoyae, Auriculariopsis ampla, Laccaria bicolor, Lentinula edodes, Lentinus tigrinus, Mycena kentingensis, Phanerochaete chrysosporium, Pleurotus ostreatus, and Schizophyllum commune, providing functional hypotheses for ~10 % of genes in the genomes of these species. Although experimental evidence for the role of these genes will need to be established in the future, our data provide a roadmap for guiding functional analyses of fruiting related genes in the Agaricomycetes. We anticipate that the gene compendium presented here, combined with developments in functional genomics approaches will contribute to uncovering the genetic bases of one of the most spectacular multicellular developmental processes in fungi. Citation: Nagy LG, Vonk PJ, Künzler M, Földi C, Virágh M, Ohm RA, Hennicke F, Bálint B, Csernetics Á, Hegedüs B, Hou Z, Liu XB, Nan S, M. Pareek M, Sahu N, Szathmári B, Varga T, Wu W, Yang X, Merényi Z (2023). Lessons on fruiting body morphogenesis from genomes and transcriptomes of Agaricomycetes. Studies in Mycology 104: 1-85. doi: 10.3114/sim.2022.104.01.

Medium optimization for high mycelial soluble protein content of Ophiocordyceps sinensis using response surface methodology

Ophiocordyceps sinensis is widely utilized due to its pharmaceutical value. Mycelial protein forms a key active component of O. sinensis and determines the medicinal potential of fungus. Here, we describe the development of an optimized fermentation medium to obtain more mycelial soluble protein from O. sinensis using response surface methodology (RSM) and investigate the increased mycelial protein content using transcriptomics. The maximum mycelial protein content of 2.11% was obtained using a medium consisting of 20% beef broth, 0.10% peptone, 2% glucose, 0.15% yeast extract, 0.20% KH₂PO₄, and 0.02% MgSO₄. Transcriptome analysis identified 790 differentially expressed genes (DEGs), including 592 up-regulated genes and 198 down-regulated genes, optimisation resulted in more up-regulated genes. The main DEGs were enriched in metabolic pathways, ABC transporters, starch and sucrose metabolism, tyrosine metabolism, and glutathione metabolism. In addition, some DEGs associated with mycelial protein enhancement such as tyrosinase (TYR), glutathione S-transferase (GST), glutamine synthetase (glnA), and β-glucosidase may contribute to increased mycelial protein content. Real-time quantitative PCR (RT-qPCR) was used to confirm gene expression and the results support the accuracy of RNA-Seq and DEG analysis. This study provides an optimized fermentation method for enhancing the mycelial protein content of O. sinensis and a reference for the effective development of O. sinensis protein.

Comparative transcriptomics reveals unique pine wood decay strategies in the Sparassis latifolia

Sparassis latifolia is a valuable edible mushroom, growing on fresh pine wood sawdust substrate. However, the mechanistic bases are poorly understood. The gene expression profiles of S. latifolia were analyzed from submerged cultures with fresh pine wood sawdust substrate for different time (0 h, 1 h, 6 h, 1 day, 5 days, and 10 days, respectively). The total number of differentially expressed genes (DEGs) identified under pine sawdust inducing was 2,659 compared to 0 h (CK). And 1,073, 520, 385, 424, and 257 DEGs were identified at the five time points, respectively. There were 34 genes in common at all inoculated time points, including FAD/NAD(P)-binding domain-containing protein, glucose methanol choline (GMC) oxidoreductase, flavin-containing monooxygenase, and taurine catabolism dioxygenase. Weighted gene co-expression analysis (WGCNA) was then used to compare the molecular characteristics among the groups and identified that the blue module had the highest correlation with the time induced by pine wood sawdust. There were 102 DEGs out of 125 genes in the blue model, which were most enriched in nitronate monooxygenase activity, dioxygenase activity, and oxidation-reduction process GO terms (p < 0.05), and peroxisome in KEGG pathway. This may provide clues into mechanisms that S. latifolia can grow on fresh pine wood sawdust substrate.

Transcriptome Analysis Explored the Differential Genes’ Expression During the Development of the Stropharia rugosoannulata Fruiting Body

Stropharia rugosoannulata (S. rugosoannulata) is a fungus with great edible and nutritional values; however, the development mechanism of its fruiting body has not been studied. Thus, this study aimed to analyze the differentially expressed genes (DEGs) in four stages; primordia stage (Sra1), young mushroom stage (Sra2), picking stage (Sra3), and opening umbrella stage (Sra4). Therefore, total RNA was extracted for further RNA-sequencing analysis. In three pairwise comparison groups (PCGs), Sra1 vs. Sra2, Sra2 vs. Sra3, and Sra3 vs. Sra4, a total of 3,112 DEGs were identified among the three PCGs. A GO analysis of the DEGs showed that there were 21 terms significantly enriched in Sra1 vs. Sra2 PCG. There was no significantly enriched GO term in the other two PCGs. Furthermore, KEGG pathway analysis showed that these DEGs were mainly enriched in glucose and amino acid metabolisms. Moreover we found that intron retention (IR) and the alternative 3′ splice site (A3SS) accounted for more than 80%. The development of the S. rugosoannulata fruiting body mainly involved glucose and amino acid metabolisms. IR and A3SS were the two main types of ASE, which played an important role in the development and maturation of the S. rugosoannulata fruiting body.

High Throughput Identification of the Potential Antioxidant Peptides in Ophiocordyceps sinensis

Ophiocordyceps sinensis, an ascomycete caterpillar fungus, has been used as a Traditional Chinese Medicine owing to its bioactive properties. However, until now the bio-active peptides have not been identified in this fungus. Here, the raw RNA sequences of three crucial growth stages of the artificially cultivated O. sinensis and the wild-grown mature fruit-body were aligned to the genome of O. sinensis. Both homology-based prediction and de novo-based prediction methods were used to identify 8541 putative antioxidant peptides (pAOPs). The expression profiles of the cultivated mature fruiting body were similar to those found in the wild specimens. The differential expression of 1008 pAOPs matched genes had the highest difference between ST and MF, suggesting that the pAOPs were primarily induced and play important roles in the process of the fruit-body maturation. Gene ontology analysis showed that most of pAOPs matched genes were enriched in terms of ‘cell redox homeostasis’, ‘response to oxidative stresses’, ‘catalase activity’, and ‘ integral component of cell membrane’. A total of 1655 pAOPs was identified in our protein-seqs, and some crucial pAOPs were selected, including catalase, peroxiredoxin, and SOD [Cu–Zn]. Our findings offer the first identification of the active peptide ingredients in O. sinensis, facilitating the discovery of anti-infectious bio-activity and the understanding of the roles of AOPs in fungal pathogenicity and the high-altitude adaptation in this medicinal fungus.

Comparative analysis of fatty acid metabolism based on transcriptome sequencing of wild and cultivated Ophiocordyceps sinensis

Background

Ophiocordyceps sinensis is a species endemic to the alpine and high-altitude areas of the Qinghai-Tibet plateau. Although O. sinensis has been cultivated since the past few years, whether cultivated O. sinensis can completely replace wild O. sinensis remains to be determined.

Methods

To explore the differences of O. sinensis grown in varied environments, we conducted morphological and transcriptomic comparisons between wild and cultivated samples who with the same genetic background.

Results

The results of morphological anatomy showed that there were significant differences between wild and cultivated O. sinensis, which were caused by different growth environments. Then, a total of 9,360 transcripts were identified using Illumina paired-end sequencing. Differential expression analysis revealed that 73.89% differentially expressed genes (DEGs) were upregulated in O. sinensis grown under natural conditions compared with that grown under artificial conditions. Functional enrichment analysis showed that some key DEGs related to fatty acid metabolism, including acyl-CoA dehydrogenase, enoyl-CoA hydratase, 3-ketoacyl-CoA thiolase, and acetyl-CoA acetyltransferase, were upregulated in wild O. sinensis. Furthermore, gas chromatography-mass spectrometry results confirmed that the fatty acid content of wild O. sinensis was significantly higher than that of cultivated O. sinensis and that unsaturated fatty acids accounted for a larger proportion.

Conclusion

These results provide a theoretical insight to the molecular regulation mechanism that causes differences between wild and cultivated O. sinensis and improving artificial breeding.

mRNA-seq and miRNA-seq profiling analyses reveal molecular mechanisms regulating induction of fruiting body in Ophiocordyceps sinensis

Ophiocordyceps sinensis has been a source of valuable materials in traditional Asian medicine for over two thousand years. With recent global warming and overharvest, however, the availability of these wild fungi has decreased dramatically. While fruiting body of O. sinensis has been artificially cultivated, the molecular mechanisms that govern the induction of fruiting body at the transcriptional and post-transcriptional levels are unclear. In this study, we carried out both mRNA and small RNA sequencing to identify crucial genes and miRNA-like RNAs (milRNAs) involved in the development of fruiting body. A total of 2875 differentially expressed genes (DEGs), and 71 differentially expressed milRNAs (DEMs) were identified among the mycoparasite complex, the sclerotium (ST) and the fruiting body stage. Functional enrichment and Gene Set Enrichment Analysis indicated that the ST had increased oxidative stress and energy metabolism and that mitogen-activated protein kinase signaling might induce the formation of fruiting body. Integrated analysis of DEGs and DEMs revealed that n_os_milR16, n_os_milR21, n_os_milR34, and n_os_milR90 could be candidate milRNAs that regulate the induction of fruiting body. This study provides transcriptome-wide insight into the molecular basis of fruiting body formation in O. Sinensis and identifies potential candidate genes for improving induction rate.

iTRAQ-based comparative proteome analyses of different growth stages revealing the regulatory role of reactive oxygen species in the fruiting body development of Ophiocordyceps sinensis

In this study, using an isobaric tags for relative and absolute quantitation (iTRAQ ) approach coupled with LC-MS / MS and bioinformatics, the proteomes were analyzed for the crucial three stages covering the fruiting body development of Ophiocordyceps sinensis, including sclerotium (ST), primordium (PR) and mature fruiting body (MF), with a focus on fruiting body development-related proteins and the potential mechanisms of the development. A total of 1,875 proteins were identified. Principal Component Analysis (PCA) demonstrated that the protein patterns between PR and MF were more similar than ST. Differentially accumulated proteins (DAPs) analysis showed that there were 510, 173 and 514 DAPs in the comparisons of ST vs. PR, PR vs. MF and ST vs. MF, respectively. A total of 62 shared DAPs were identified and primarily enriched in proteins related to ‘carbon transport and mechanism’, ‘the response to oxidative stress’, ‘antioxidative activity’ and ‘translation’. KEGG and GO databases showed that the DAPs were enriched in terms of ‘primary metabolisms (amino acid/fatty acid/energy metabolism)’, ‘the response to oxidative stress’ and ‘peroxidase’. Furthermore, 34 DAPs involved in reactive oxygen species (ROS) metabolism were identified and clustered across the three stages using hierarchical clustering implemented in hCluster R package . It was suggested that their roles and the underlying mechanisms may be stage-specific. ROS may play a role in fungal pathogenicity in ST, the fruit-body initiation in PR, sexual reproduction and highland adaptation in MF. Crucial ROS-related proteins were identified, such as superoxide dismutase (SOD, T5A6F1), Nor-1 (T5AFX3), electron transport protein (T5AHD1), histidine phosphotransferase (HPt, T5A9Z5) and Glutathione peroxidase (T5A9V1). Besides, the accumulation of ROS at the three stages were assayed using 2,7-dichlorofuorescin diacetate (DCFH-DA) stanning. A much stronger ROS accumulation was detected at the stage MF, compared to the stages of PR and ST. Sections of ST and fruit-body part of MF were stained by DCFH-DA and observed under the fluorescencemicroscope, showing ROS was distributed within the conidiospore and ascus. Besides, SOD activity increased across the three stages, while CAT activity has a strong increasement in MF compared to the stages of ST and PR. It was suggested that ROS may act in gradient-dependent manner to regulate the fruiting body development. The coding region sequences of six DAPs were analyzed at mRNA level by quantitative real-time PCR (qRT-PCR). The results support the result of DAPs analysis and the proteome sequencing data. Our findings offer the perspective of proteome to understand the biology of fruiting body development and highland adaptation in O. sinensis, which would inform the big industry of this valuable fungus.