Genome-wide transcriptome and proteome analysis on different developmental stages of Cordyceps militaris

Expanded Gene Regulatory Network Reveals Potential Light-Responsive Transcription Factors and Target Genes in Cordyceps militaris

Cordyceps militaris, a fungus widely used in traditional Chinese medicine and pharmacology, is recognized for its abundant bioactive compounds, including cordycepin and carotenoids. The growth, development, and metabolite production in various fungi are influenced by the complex interactions between regulatory cascades and light-signaling pathways. However, the mechanisms of gene regulation in response to light exposure in C. militaris remain largely unexplored. This study aimed to identify light-responsive genes and potential transcription factors (TFs) in C. militaris through an integrative transcriptome analysis. To achieve this, we reconstructed an expanded gene regulatory network (eGRN) comprising 507 TFs and 8662 regulated genes using both interolog-based and homolog-based methods to build the protein-protein interaction network. Aspergillus nidulans and Neurospora crassa were chosen as templates due to their relevance as fungal models and the extensive study of their light-responsive mechanisms. By utilizing the eGRN as a framework for comparing transcriptomic responses between light-exposure and dark conditions, we identified five key TFs-homeobox TF (CCM_07504), FlbC (CCM_04849), FlbB (CCM_01128), C6 zinc finger TF (CCM_05172), and mcrA (CCM_06477)-along with ten regulated genes within the light-responsive subnetwork. These TFs and regulated genes are likely crucial for the growth, development, and secondary metabolite production in C. militaris. Moreover, molecular docking analysis revealed that two novel TFs, CCM_05727 and CCM_06992, exhibit strong binding affinities and favorable docking scores with the primary light-responsive protein CmWC-1, suggesting their potential roles in light signaling pathways. This information provides an important functional interactive network for future studies on global transcriptional regulation in C. militaris and related fungi.

Holistic transcriptional responses of Cordyceps militaris to different culture temperatures

Cordyceps militaris is a medicinal entomopathogenic fungus containing valuable biometabolites for pharmaceutical applications. Its genetic inheritance and environmental factors play a crucial role in the production of biomass enriched with cordycepin. While temperature is a crucial controlled parameter for fungal cultivation, its impacts on growth and metabolite biosynthesis remains poorly characterized. This study aimed to investigate the metabolic responses and cordycepin production of C. militaris strain TBRC6039 under various temperature conditions through transcriptome analysis. Among 9599 expressed genes, 576 genes were significantly differentially expressed at culture temperatures of 15 and 25 °C. The changes in the transcriptional responses induced by these temperatures were found in several metabolisms involved in nutrient assimilation and energy source, including amino acids metabolism (e.g., glycine, serine and threonine metabolism) and lipid metabolism (e.g., biosynthesis of unsaturated fatty acids and steroid biosynthesis). At the lower temperature (15 °C), the biosynthetic pathways of lipids, specifically ergosterol and squalene, were the target for maintaining membrane function by transcriptional upregulation. Our study revealed the responsive mechanisms of C. militaris in acclimatization to temperature conditions that provide an insight on physiological manipulation for the production of metabolites by C. militaris.

A novel complementary pathway of cordycepin biosynthesis in Cordyceps militaris

We determined whether there exists a complementary pathway of cordycepin biosynthesis in wild-type Cordyceps militaris, high-cordycepin-producing strain C. militaris GYS60, and low-cordycepin-producing strain C. militaris GYS80. Differentially expressed genes were identified from the transcriptomes of the three strains. Compared with C. militaris, in GYS60 and GYS80, we identified 145 and 470 upregulated and 96 and 594 downregulated genes. Compared with GYS80, in GYS60, we identified 306 upregulated and 207 downregulated genes. Gene Ontology analysis revealed that upregulated genes were mostly involved in detoxification, antioxidant, and molecular transducer in GYS60. By Clusters of Orthologous Groups of Proteins and Kyoto Encyclopedia of Genes and Genomes analyses, eight genes were significantly upregulated: five genes related to purine metabolism, one to ATP production, one to secondary metabolite transport, and one to RNA degradation. In GYS60, cordycepin was significantly increased by upregulation of ATP production, which promoted 3',5'-cyclic AMP production. Cyclic AMP accelerated 3'-AMP accumulation, and cordycepin continued to be synthesized and exported. We verified the novel complementary pathway by adding the precursor adenosine and analyzing the expression of four key genes involved in the main pathway of cordycepin biosynthesis. Adenosine addition increased cordycepin production by 51.2% and 10.1%, respectively, in C. militaris and GYS60. Four genes in the main pathway in GYS60 were not upregulated.

Comparative genomics of fungal mutants provides a systemic view of extreme cadmium tolerance in eukaryotic microbes

Whether eukaryotic organisms can evolve for higher heavy metal resistance in laboratory conditions remains unknown. In this study, we challenged a macrofungi, Pleurotus ostreatus, in a designed microbial evolution and growth arena (MEGA)-plate with an extreme Cd gradient. Within months, the wild-type strain developed 10 mutants, exhibiting a maximum three-fold increase in Cd tolerance and slower growth rates. Genomic sequencing and re-sequencing of the wild-type and ten mutant strains generated about 51 GB data, allowing a comprehensive comparative genomics analysis. As a result, a total of 2512 common single nucleotide polymorphisms, 70 inserts and deletes, 39 copy number variations and 21 structural variations were found in the 10 mutants. The mutant genes were primarily involved in substrate transport. In combination with transcriptome analysis, we discovered that the ten mutants had a distinct Cd-resistant mechanism compared to the wild-type strain. Genes involved in oxidation-reduction, ion transmembrane transport, and metal compartment/efflux are primarily responsible for the extreme Cd tolerance in the P. ostreatus mutants. Our findings contribute to the understanding of eukaryotic Cd resistance at the genome level and establish a foundation for developing bioremediation tools utilizing highly tolerant macrofungi.

Evidence for Regulation of Cordycepin Biosynthesis by Transcription Factors Krüppel-Like Factor 4 and Retinoid X Receptor Alpha in Caterpillar Medicinal Mushroom Cordyceps militaris (Ascomycetes)

Cordyceps militaris, Chinese traditional medicinal fungus, has many bioactive properties. Cordycepin (3'-deoxyadenosine) is a major bioactive component of C. militaris. Various methods can significantly elevate cordycepin production, which suggests a diverse set of metabolic regulatory mechanisms. Thus, we aimed to identify transcription factors that regulate cordycepin biosynthesis pathways. Transcriptome analysis of wild-type C. militaris, C. militaris GYS60, a cordycepin high-producing strain, and C. militaris GYS80, a low-producing strain, were used to measure expression and function of genes related to cordycepin biosynthesis. The transcriptome expression data were confirmed by quantitative real-time polymerase chain reaction. We identified 155 relevant transcription factors in 19 families that included Fork head/winged helix factors, other C4 zinc finger-type factors, C2H2 zinc finger factors, tryptophan cluster factors, nuclear receptors with C4 zinc fingers, homeodomain factors, and Rel homology region factors. Energy generation and amino acid conversion pathways were activated in GYS60 so that abundance of cordycepin precursors was increased. Genes and transcription factors for rate-limiting enzymes in these pathways were identified. Overexpression of two key transcription factors, Kruppel-like factor 4 (Klf4) and Retinoid X receptor alpha (Rxra), promoted high cordycepin production in GYS60. In GYS60, Klf4 and Rxra were responsible for upregulation of genes in cordycepin biosynthesis, namely an oxidoreductase, 3',5'-cyclic AMP phosphodiesterase, a transferase, and adenylate cyclase. Upregulation of these genes increased 3'-AMP content, thereby elevating cordycepin synthesis.

Comparative proteomics reveals the mechanism of cyclosporine production and mycelial growth in Tolypocladium inflatum affected by different carbon sources

Cyclosporine A (CsA) is a secondary cyclopeptide metabolite produced by Tolypocladium inflatum that is widely used clinically as an immunosuppressant. CsA production and mycelial growth differed when T. inflatum was cultured in different carbon source media. During early fermentation, CsA was preferred to be produced in fructose medium, while the mycelium preferred to accumulate in sucrose medium. On the sixth day, the difference was most pronounced. In this study, high-throughput comparative proteomics methods were applied to analyze differences in protein expression of mycelial samples on day 6, revealing the proteins and mechanisms that positively regulate CsA production related to carbon metabolism. The differences included small molecule acid metabolism, lipid metabolism, organic catabolism, exocrine secretion, CsA substrate Bmt synthesis, and transcriptional regulation processes. The proteins involved in the regulation of mycelial growth related to carbon metabolism were also revealed and were associated with waste reoxidation processes or coenzyme metabolism, small molecule synthesis or metabolism, the stress response, genetic information or epigenetic changes, cell component assembly, cell wall integrity, membrane metabolism, vesicle transport, intramembrane localization, and the regulation of filamentous growth. This study provides a reliable reference for CsA production from high-efficiency fermentation. This study provides key information for obtaining more CsA high-yielding strains through metabolic engineering strategies.

High-level production of cordycepin by the xylose-utilising Cordyceps militaris strain 147 in an optimised medium

Cordycepin is an important active metabolite of Cordyceps militaris. Xylose, an attractive feedstock for producing chemicals through microbial fermentation, cannot be effectively utilised by many reported C. militaris strains. Herein, a xylose-utilising C. militaris strain 147 produced the highest level of cordycepin (3.03 g/L) in xylose culture. Xylose, alanine, and ammonium citrate were determined as the main affecting factors on the cordycepin production using a Plackett-Burman design. The combination of these factors was optimised using response surface methodology, and the maximal 6.54 g/L of cordycepin was produced by the fungus in the optimal medium. Transcriptome analysis revealed that xylose utilisation upregulated the transcriptional levels of genes participating in purine and energy metabolisms in the fungus, which may facilitate the formation of precursors for cordycepin biosynthesis. This investigation provides new insights into the efficient production of cordycepin and is conducive to the valorisation of biomass rich in xylose.

Comparative analysis of proteomes and transcriptomes revealed the molecular mechanism of development and nutrition of Pleurotus giganteus at different fruiting body development stages

Pleurotus giganteus is a commercially cultivated high-temperature mushroom. Investigating the molecular mechanism of fruiting body development will help us to better understand the regulation of substrates and energy in this process. However, little information has been reported on the development and nutrients of the P. giganteus fruiting body. In the present study, P. giganteus is cultivated in a climate chamber, and comparative transcriptome, proteome, and nutritional analysis of P. giganteus fruiting bodies were performed. Our results revealed that Cytochrome P450 monooxygenases and hydrophobin proteins play important roles during the differentiation in the elongation stage. Later, carbon metabolism dominate the fruiting body metabolism and genes related to the carbohydrate metabolic process, glycolytic process, and gluconeogenesis were up-regulated in the mature fruiting bodies. The up-regulation of carbohydrate substrates utilization CAZymes genes and inconsistent protein expression in pileus indicated a reverse transportation of mRNA from the fruiting body to vegetative mycelia. In addition, protein concentration in the pileus is higher than that in the stem, while the stem is the major nitrogen metabolic and amino acid synthetic location. The integrated transcriptomic, proteomic, and nutritional analysis indicated a two-way transportation of substrates and mRNAs in P. giganteus. Stem synthesizes amino acids and transported them to pileus with reducing sugars, while pileus induces the expression of substrate degradation mRNA according to the needs of growth and development and transports them in the other direction.

Functional genomics and systems biology of Cordyceps species for biotechnological applications

The use of Cordyceps species for the manufacture of natural products has been established; however, the tremendous advances observed in recent years in genetic engineering and molecular biology have revolutionized the optimization of Cordyceps as cell factories and drastically expanded the biotechnological potential of these fungi. Here, we present a review of systems and synthetic biology studies of Cordyceps and their implications for fungal biology and industrial applications. We summarize the current status of synthetic biology for enhancing targeted metabolites in Cordyceps species, such as cordycepin, adenosine, polysaccharide, and pentostatin. Progress in the systems and synthetic biology of Cordyceps provides a strategy for comprehensively comprehensive controlling efficient cell factories of natural bioproducts and novel synthetic biology toolbox for targeted engineering.

Comparative proteome profiling of fusarium chlamydosporum and elucidation of pigment biosynthetic pathway under nitrogen stress

The main aim of this study was to understand the protein expression of F. chlamydosporum in two different medium composition in varying concentrations of nitrogen. The interesting phenomenon of producing diverse pigments by a single strain in different concentrations of Nitrogen made us further to explore the difference in the protein expression of the fungus when grown in these two media. For this, we had adopted non-gel-based method of protein separation by LC-MS/MS analysis followed by label free identification of proteins by SWATH analysis. The molecular and biological functions of each protein and their Gene Ontology annotations were analyzed by UniProt KB and KEGG pathway; the secondary metabolite pathways and the carbohydrate metabolic pathways were analyzed by DAVID bioinformatics tool. The positively regulated proteins biologically functioned for the secondary metabolite production in optimized medium were Diphosphomevalonate decarboxylase (terpenoid backbone biosynthesis), Phytoene synthase (carotenoid biosynthesis), 6,7-dimethyl-8-ribityllumazine synthase (riboflavin biosynthesis). The main characteristic change observed was that proteins related to carotenoid biosynthesis and terpenoid synthesis were not regulated in nitrogen limited medium. Except for the protein 6,7-dimethyl-8-ribityllumazine synthase, all the enzymes related to fatty acid biosynthesis and polyketide chain elongation were up regulated. Apart from the proteins related to secondary metabolite production, two novel proteins were found to be up regulated in Nitrogen Limited Medium; C-fem protein responsible for fungal pathogenesis and DAO domain containing protein which functions as a neuromodulator and catalyzes the synthesis of dopamine. SIGNIFICANCE: This particular strain of F. chlamydosporum of immense genetic and biochemical diversity represents an interesting example of a microorganism which can produce a variety of bioactive compounds and this can be exploited in various industries. The production of carotenoids and polyketides by this fungus when grown in the same media with different concentrations of Nitrogen has been published by us following which we analyzed the proteome sequence of the fungus in varying Nutrient conditions. Following the proteome analysis and expression, we could derive the pathway leading to the biosynthesis of varying secondary metabolites by the fungus which has not been published or studied so far.

Increasing the production of the bioactive compounds in medicinal mushrooms: an omics perspective

Macroscopic fungi, mainly higher basidiomycetes and some ascomycetes, are considered medicinal mushrooms and have long been used in different areas due to their pharmaceutically/nutritionally valuable bioactive compounds. However, the low production of these bioactive metabolites considerably limits the utilization of medicinal mushrooms both in commerce and clinical trials. As a result, many attempts, ranging from conventional methods to novel approaches, have been made to improve their production. The novel strategies include conducting omics investigations, constructing genome-scale metabolic models, and metabolic engineering. So far, genomics and the combined use of different omics studies are the most utilized omics analyses in medicinal mushroom research (both with 31% contribution), while metabolomics (with 4% contribution) is the least. This article is the first attempt for reviewing omics investigations in medicinal mushrooms with the ultimate aim of bioactive compound overproduction. In this regard, the role of these studies and systems biology in elucidating biosynthetic pathways of bioactive compounds and their contribution to metabolic engineering will be highlighted. Also, limitations of omics investigations and strategies for overcoming them will be provided in order to facilitate the overproduction of valuable bioactive metabolites in these valuable organisms.

Construction of Cordycepin High-Production Strain and Optimization of Culture Conditions

This study aimed to increase cordycepin production by over-expressing bio-synthetic enzyme genes, including the adenylosuccinate synthase, adenylosuccinate lyase, and 5'-nucleotidase genes. Research data showed that the extracellular and intracellular cordycepin concent of 24 recombinant strains were higher than those of C. militaris WT, indicating that over-expression of key enzyme genes increased cordycepin production. Among them, the CM-adss-5 strain had highest cordycepin production, and the extracellular and intracellular cordycepin concent were 1119.75 ± 1.61 and 65.56 ± 0.97 mg/L, which were 1.26 and 2.61 times that of C. militaris WT. This study also optimized the culture conditions of CM-adss-5 strain through single factor experiments to obtain the best culture conditions. The best culture condition was 25 °C constant temperature, 180-rpm shaking culture, fermentation period 12 days, inoculate amount 5%, initial pH 6, seed age 108 h, and liquid volume 110/250 mL. Then, the extracellular and intracellular cordycepin content of CM-adss-5 strain reached 2581.96 ± 21.07 and 164.08 ± 1.44 mg/L, which were higher by 130.6% and 150.3%, respectively. Therefore, our research provides a way to efficiently produce cordycepin for the development of cordycepin and its downstream products.

Ttc39c is a potential target for the treatment of lung cancer

BACKGROUND

The novel TTC gene, tetratricopeptide repeat domain 39 C (Ttc39c), mainly mediates the interaction between proteins. It is involved in the progression of various tumors. In this study, we determined the effect of Ttc39c on lung adenocarcinoma and found that it might be used as a potential intervention target.

METHODS

We performed a difference analysis of Ttc39c samples from the TCGA database. Transwell experiments were conducted to determine the ability of cell metastasis. Celigo and MTT assays were performed to determine the effect of Ttc39c gene subtraction on cell proliferation. FACS was performed to determine the effect of Ttc39c gene subtraction on apoptosis. Clone-formation experiments were conducted to determine the effect of Ttc39c gene subtraction on cloning ability. Transcriptomics, proteomics, and metabolomics were used to elucidate the enrichment pathway of the Ttc39c gene in the progression of lung adenocarcinoma.

RESULTS

The expression of Ttc39c increased significantly in lung adenocarcinoma. The proliferation, metastasis, and cloning ability of human lung cancer cells were inhibited, while the apoptosis of cells increased significantly after the depletion of Ttc39c. Our results based on the transcriptomics, proteomics, and metabolomics analyses indicated that Ttc39c might be involved in the progression of lung adenocarcinoma (LUAD) mainly through the metabolic pathway and the p53 pathway.

CONCLUSION

To summarize, Ttc39c strongly regulates the proliferation and metastasis of lung adenocarcinoma cells. The main pathways involved in Ttc39c in lung adenocarcinoma include the energy metabolism and p53 pathways.

Genomic Comparative Analysis of Cordyceps pseudotenuipes with Other Species from Cordyceps

The whole genome of Cordyceps pseudotenuipes was sequenced, annotated, and compared with three related species to characterize the genome. The antibiotics and Secondary Metabolites Analysis Shell (antiSMASH) and local BLAST analysis were used to explore the secondary metabolites (SMs) and biosynthesis gene clusters (BGCs) of the genus Cordyceps. The genome-wide basic characteristics of C. pseudotenuipes, C. tenuipes, C. cicadae, and C. militaris revealed unequal genome size, with C. cicadae as the largest (34.11 Mb), followed by C. militaris (32.27 Mb). However, the total gene lengths of C. pseudotenuipes and C. tenuipes were similar (30.1 Mb and 30.06 Mb). The GC contents of C. pseudotenuipes, C. tenuipes, C. cicadae, and C. militaris genomes differed slightly (51.40% to 54.11%). AntiSMASH and local BLAST analysis showed that C. pseudotenuipes, C. tenuipes, C. cicadae, and C. militaris had 31, 28, 31, and 29 putative SM BGCs, respectively. The SM BGCs contained different quantities of polyketide synthetase (PKS), nonribosomal peptide synthetase (NRPS), terpene, hybrid PKS + NRPS, and hybrid NRPS + Other. Moreover, C. pseudotenuipes, C. tenuipes, C. cicadae, and C. militaris had BGCs for the synthesis of dimethylcoprogen. C. pseudotenuipes, C. tenuipes, and C. cicadae had BGCs for the synthesis of leucinostatin A/B, neosartorin, dimethylcoprogen, wortmanamide A/B, and beauvericin. In addition, the SM BGCs unique to C. pseudotenuipes were clavaric acid, communesin, and deoxynivalenol. Synteny analysis indicated that the scaffolds where the SM BGC was located were divided into more than 70 collinear blocks, and there might be rearrangements. Altogether, these findings improved our understanding of the molecular biology of the genus Cordyceps and will facilitate the discovery of new biologically active SMs from the genus Cordyceps using heterologous expression and gene knockdown methods.

Comprehensive Transcriptomic Analysis of Cordyceps militaris Cultivated on Germinated Soybeans

The ascomycete fungus Cordyceps militaris infects lepidopteran larvae and pupae and forms characteristic fruiting bodies. Owing to its immune-enhancing effects, the fungus has been used as a medicine. For industrial application, this fungus can be grown on geminated soybeans as an alternative protein source. In our study, we performed a comprehensive transcriptomic analysis to identify core gene sets during C. militaris cultivation on germinated soybeans. RNA-Seq technology was applied to the fungal cultures at seven-time points (2, 4, and 7-day and 2, 3, 5, 7-week old cultures) to investigate the global transcriptomic change. We conducted a time-series analysis using a two-step regression strategy and chose 1460 significant genes and assigned them into five clusters. Characterization of each cluster based on Gene Ontology and Kyoto Encyclopedia of Genes and Genomes databases revealed that transcription profiles changed after two weeks of incubation. Gene mapping of cordycepin biosynthesis and isoflavone modification pathways also confirmed that gene expression in the early stage of GSC cultivation is important for these metabolic pathways. Our transcriptomic analysis and selected genes provided a comprehensive molecular basis for the cultivation of C. militaris on germinated soybeans.

Transcriptomic Responses of Cordyceps militaris to Salt Treatment During Cordycepins Production

Cordycepin is a major bioactive compound found in Cordyceps militaris (C. militaris) that exhibits a broad spectrum of biological activities. Hence, it is potentially a bioactive ingredient of pharmaceutical and cosmetic products. However, overexploitation and low productivity of natural C. militaris is a barrier to commercialization, which leads to insufficient supply to meet its existing market demands. In this study, a preliminary study of distinct concentrations of salt treatments toward C. militaris was conducted. Although the growth of C. militaris was inhibited by different salt treatments, the cordycepin production increased significantly accompanied by the increment of salt concentration. Among them, the content of cordycepin in the 7% salt-treated group was five-fold higher than that of the control group. Further transcriptome analysis of samples with four salt concentrations, coupled with Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment, several differentially expressed genes (DEGs) were found. Finally, dynamic changes of the expression patterns of four genes involved in the cordycepin biosynthesis pathway were observed by the quantitative real-time PCR. Taken together, our study provides a global transcriptome characterization of the salt treatment adaptation process in C. militaris and facilitates the construction of industrial strains with a high cordycepin production and salt tolerance.

Comparative Proteomic Analysis within the Developmental Stages of the Mushroom White Hypsizygus marmoreus

(1) Background: The white Hypsizygus marmoreus is a popular edible mushroom in East Asia markets. Research on the systematic investigation of the protein expression changes in the cultivation process of this mushroom are few. (2) Methods: Label-free LC-MS/MS quantitative proteomics analysis technique was adopted to obtain the protein expression profiles of six groups of samples collected in different growth stages. A total of 3468 proteins were identified. The UpSetR plot analysis, Pearson correlation coefficient (PCC) analysis, and principal component (PC) analysis were performed to reveal the correlation among the six groups of samples. The differentially expressed proteins (DEPs) were organised by One-way ANOVA test and divided into four clusters. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were performed to divide the DEPs into different metabolic processes and pathways in each cluster. (3) Results: The DEPs in cluster 1 are of the highest abundance in the mycelium and are mainly involved in protein biosynthesis, biosynthesis of cofactors, lipid metabolism, spliceosome, cell cycle regulation, and MAPK signaling pathway. The DEPs in cluster 2 are enriched in the stem and are mainly associated with protein biosynthesis, biosynthesis of cofactors, carbon, and energy metabolism. The DEPs in cluster 3 are highly expressed in the primordia and unmatured fruiting bodies and are related to amino acids metabolism, carbon and carbohydrate metabolism, protein biosynthesis and processing, biosynthesis of cofactors, cell cycle regulation, MAPK signaling pathway, ubiquitin-mediated proteolysis, and proteasome. The DEPs in cluster 4 are of the highest abundance in the cap and are mainly associated with spliceosome, endocytosis, nucleocytoplasmic transport, protein processing, oxidative phosphorylation, biosynthesis of cofactors, amino acids metabolism, and lipid metabolism. (4) Conclusions: This research reports the proteome analysis of different developmental stages during the cultivation of the commercially relevant edible fungi the white H. marmoreus. In the mycelium stage, most of the DEPs are associated with cell proliferation, signal response, and mycelium growth. In the primordia and unmatured fruiting bodies stage, the DEPs are mainly involved in biomass increase, cell proliferation, signal response, and differentiation. In the mature fruiting body stage, the DEPs in the stem are largely associated with cell elongation and increase in biomass, and most of the DEPs in the cap are mainly related to pileus expansion. Several carbohydrate-active enzymes, transcription factors, heat shock proteins, and some DEPs involved in MAPK and cAMP signaling pathways were determined. These proteins might play vital roles in metabolic processes and activities. This research can add value to the understanding of mechanisms concerning mushroom development during commercial production.

Comparative transcriptome analysis reveals the differential response to cadmium stress of two Pleurotus fungi: Pleurotus cornucopiae and Pleurotus ostreatus

Pleurotus has great potential for heavy metal mycoremediation. Using comparative transcriptome analysis, the response of Pleurotus ostreatus and Pleurotus cornucopiae under Cd contamination was evaluated. P. ostreatus and P. cornucopia accumulated 0.34 and 0.46 mg/g Cd in mycelium, respectively. Cd removal elevated with its concentration elevation, which reached 56.47% and 54.60% for P. ostreatus and P. cornucopia with Cd at 20 mg/L. Low-level Cd (≤ 1 mg/L) had no significant influence on either fungus, while varied response was observed under high-level Cd. 705 differentially expressed genes (DEGs) were identified in P. cornucopia at Cd1 and Cd20, whereas 12,551 DEGs in P. ostreatus. Differentially regulated functional categories and pathways were also identified. ATP-binding cassette transporters were involved in Cd transport in P. cornucopia, whereas the endocytosis and phagosome pathways were more enhanced in P. ostreatus. 26 enzymes including peroxisomal enzymes catalase and superoxide dismutase were upregulated in P. ostreatus, whereas only cytosolic catalase was overexpressed in P. cornucopia, suggesting their different Cd detoxification pathways. Also, the mitogen-activated protein kinase signaling pathway involved in Cd resistance in both species instead of glutathione metabolism, although more active in P. ostreatus. These findings provided new insight into the molecular mechanism of mycoremediation and accumulator screening.

Deep sequencing of the Sanghuangporus vaninii transcriptome reveals dynamic landscapes of candidate genes involved in the biosynthesis of active compounds

The medicinal fungus Sanghuang produces diverse bioactive compounds and is widely used in Asian countries. However, little is known about the genes and pathways involved in the biosynthesis of these active compounds. Based on our previous study providing Sanghuangporus vaninii genomic information, the transcriptomes of MY (mycelium), OY (1-year-old fruiting bodies), and TY (3-year-old fruiting bodies) were determined in this study. A significant number of genes (4774) were up- or downregulated between mycelium and fruiting bodies, but only 1422 differentially expressed genes were detected between OY and TY. 138 genes encoding P450s were identified in the fungal genome and grouped into 25 P450 families; more than 64% (88) of the genes were significantly differentially expressed between the mycelium and fruiting body, suggesting that these P450s are involved in fungal sexual development. Importantly, the expression of genes involved in bioactive compound (triterpenoids, polysaccharides, and flavonoids) biosynthesis in asexual (cultured with solid and liquid media) and sexual stages was explored and combined with transcriptome and quantitative PCR analyses. More genes involved in the biosynthesis of bioactive compounds were expressed more highly in mycelium than in fruiting bodies under liquid medium culture compared with solid medium culture, which was consistent with the yields of different bioactive compounds, suggesting that liquid fermentation of S. vaninii Kangneng can be used to obtain these bioactive compounds. A comprehensive understanding of the genomic information of S. vaninii will facilitate its potential use in pharmacological and industrial applications.

New insight into the rapid growth of the Mikania micrantha stem based on DIA proteomic and RNA-Seq analysis

Mikania micrantha is one of the world's most invasive plants, which causes severe damage to natural ecosystems and agroforestry systems due to its rapid stem growth. This work investigated the proteomic and transcriptomic profiles of M. micrantha in different stem tissues (pre-internode, post-internode, and internode), as well as in adventitious roots and primary roots with the final goal of elucidating differentially expressed genes and proteins responsible for the rapid growth of stem. The objective was approached by using DIA-based proteomic and RNA-Seq technologies. More than seven giga-transcriptome clean reads were sequenced, and 5196 protein species were identified. Differentially expressed genes identified in all stem tissues were significantly enriched in photosynthesis and carbon fixation, suggesting that the stem possesses a strong photosynthetic capacity in order to maintain the energy supply for this species. Analysis of differentially expressed proteins showed that proteins related to photosystem I/II and the cytochrome b6/f complex, such as D1, D2, and cp43, were also highly accumulated in the adventitious roots, corroborating the transcriptome analysis results. These results provided basic proteomic and transcriptional expression information about the M. micrantha stem and adventitious root, thereby improving our understanding of the molecular mechanism underlying rapid growth in this species. SIGNIFICANCE: This is the first study to investigate the proteomic and transcriptomic profiles of Mikania micrantha, a highly invasive plant, in different stem tissues (pre-internode, post-internode, and internode), as well as in adventitious and primary roots, using the latest DIA-based (data-independent acquisition mode) proteomic and RNA-Seq technologies. A comprehensive study was carried out, and differentially expressed genes and differentially expressed proteins identified in the pre-internode, post-internode, and internode tissues were significantly enriched during photosynthesis and carbon fixation, suggesting that the M. micrantha stem possesses a strong photosynthetic capacity that allows the plant to maintain a high energy supply. Enriched plant hormone signal transduction pathway analysis revealed an interaction between auxin and other phytohormones involved in adventitious root development. The study provided basic data on the molecular mechanism of M. micrantha vegetative propagation and the rapid growth of its stem. The novel scientific content of this study successfully builds upon the limited information currently available on the subject, therefore warranting publication.

Proteome sequencing and analysis of Ophiocordyceps sinensis at different culture periods

Background

Ophiocordyceps sinensis is an important traditional Chinese medicine for its comprehensive active ingredients, such as cordycepin, cordycepic acid, and Cordyceps polysaccharide. O. sinensis zjut, a special strain isolated from O. sinensis, has similar pharmacological functions to wild O. sinensis. Currently, O. sinensis with artificial cultivation has been widely studied, but systematic fundamental research at protein levels has not been determined.

Results

Proteomes of O. sinensis zjut at different culture periods (growth period, 3rd day; pre-stable period, 6th day; and stable period, 9th day) were relatively quantified by relative isotope markers and absolute quantitative technology. In total, 4005 proteins were obtained and further annotated with Gene Ontology, Kyoto Encyclopedia of Genes and Genomes database. Based on the result of the annotations, metabolic pathways of active ingredients, amino acids and fatty acid were constructed, and the related enzymes were exhibited. Subsequently, comparative proteomics of O. sinensis zjut identified the differentially expressed proteins (DEPs) by growth in different culture periods, to find the important proteins involved in metabolic pathways of active ingredients. 605 DEPs between 6d-VS-3d, 1188 DEPs between 9d-VS-3d, and 428 DEPs between 9d-VS-6d were obtained, respectively.

Conclusion

This work provided scientific basis to study protein profile and comparison of protein expression levels of O. sinensis zjut, and it will be helpful for metabolic engineering works to active ingredients for exploration, application and improvement of this fungus.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07298-z.

Transcriptome and proteome analyses reveal the regulatory networks and metabolite biosynthesis pathways during the development of Tolypocladium guangdongense

Tolypocladium guangdongense has a similar metabolite profile to Ophiocordyceps sinensis, a highly regarded fungus used for traditional Chinese medicine with high nutritional and medicinal value. Although the genome sequence of T. guangdongense has been reported, relatively little is known about the regulatory networks for fruiting body development and about the metabolite biosynthesis pathways. In order to address this, an analysis of transcriptome and proteome at differential developmental stages of T. guangdongense was performed. In total, 9076 genes were found to be expressed and 2040 proteins were identified. There were a large number of genes that were significantly differentially expressed between the mycelial stage and the stages. Interestingly, the correlation between the transcriptomic and proteomic data was low, suggesting the importance of the post-transcriptional processes in the growth and development of T. guangdongense. Among the genes/proteins that were both differentially expressed during the developmental process, there were numerous heat shock proteins and transcription factors. In addition, there were numerous proteins involved in terpenoid, ergosterol, adenosine and polysaccharide biosynthesis that also showed significant downregulation in their expression levels during the developmental process. Furthermore, both tryptophan and tryptamine were present at higher levels in the primordium stage. However, indole-3-acetic acid (IAA) levels continuously decreased as development proceeded, and the enzymes involved in IAA biosynthesis were also clearly differentially downregulated. These data could be meaningful in studying the molecular mechanisms of fungal development, and for the industrial and medicinal application of macro-fungi.

Transcriptome Analysis Reveals the Flexibility of Cordycepin Network in Cordyceps militaris Activated by L-Alanine Addition

Cordycepin, isolated from the traditional medicinal fungus Cordyceps militaris, has gained much attention due to its various clinical functions. Previous reports of L-alanine addition could significantly improve cordycepin production, but the molecular mechanism remains unknown. In this study, transcriptome analysis of C. militaris with doubled cordycepin production induced by L-alanine addition provides an insight into the flexibility of the cordycepin network. The biopathways of energy generation and amino acid conversion were activated so that cordycepin substrate generation was consequently improved. Specific genes of rate-limiting enzymes in these pathways, as well as related transcription factors, were figured out. Two key Zn2Cys6-type transcription factors CmTf1 and CmTf2 were verified to play the roles of doubling the cordycepin production by overexpression of their coding genes in C. militaris wild type. These results provide a complete map of the cordycepin network in C. militaris with a distinct understanding of the flexibility of joints, giving a better foundation for increasing cordycepin yield and strain breeding in the future.

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

Ophiocordyceps sinensis is a highly valued fungus that has been used as traditional Asian medicine. This fungus is one of the most important sources of income for the nomadic populations of the Tibetan Plateau. With global warming and excessive collection, the wild O. sinensis resources declined dramatically. The cultivation of O. sinensis hasn’t been fully operational due to the unclear genetic basis of the fruiting body development. Here, our study conducted pairwise comparisons between transcriptomes acquired from different growth stages of O. sinensis including asexual mycelium (CM), developing fruiting body (DF) and mature fruiting body (FB). All RNA-Seq reads were aligned to the genome of O. sinensis CO18 prior to comparative analyses. Cluster analysis showed that the expression profiles of FB and DF were highly similar compared to CM. Alternative splicing analysis (AS) revealed that the stage-specific splicing genes may have important functions in the development of fruiting body. Functional enrichment analyses showed that differentially expressed genes (DEGs) were enriched in protein synthesis and baseline metabolism during fruiting body development, indicating that more protein and energy might be required for fruiting body development. In addition, some fruiting body development-associated genes impacted by ecological factors were up-regulated in FB samples, such as the nucleoside diphosphate kinase gene (ndk), β subunit of the fatty acid synthase gene (cel-2) and the superoxide dismutase gene (sod). Moreover, the expression levels of several cytoskeletons genes were significantly altered during all these growth stages, suggesting that these genes play crucial roles in both vegetative growth and the fruiting body development. Quantitative PCR (qPCR) was used to validate the gene expression profile and the results supported the accuracy of the RNA-Seq and DEGs analysis. Our study offers a novel perspective to understand the underlying growth stage-specific molecular differences and the biology of O. sinensis fruiting body development.

Combining transcriptomics and metabolomics to reveal the underlying molecular mechanism of ergosterol biosynthesis during the fruiting process of Flammulina velutipes

Background

Flammulina velutipes has been recognized as a useful basidiomycete with nutritional and medicinal values. Ergosterol, one of the main sterols of F. velutipes is an important precursor of novel anticancer and anti-HIV drugs. Therefore, many studies have focused on the biosynthesis of ergosterol and have attempted to upregulate its content in multiple organisms. Great progress has been made in understanding the regulation of ergosterol biosynthesis in Saccharomyces cerevisiae. However, this molecular mechanism in F. velutipes remains largely uncharacterized.

Results

In this study, nine cDNA libraries, prepared from mycelia, young fruiting bodies and mature fruiting bodies of F. velutipes (three replicate sets for each stage), were sequenced using the Illumina HiSeq™ 4000 platform, resulting in at least 6.63 Gb of clean reads from each library. We studied the changes in genes and metabolites in the ergosterol biosynthesis pathway of F. velutipes during the development of fruiting bodies. A total of 13 genes (6 upregulated and 7 downregulated) were differentially expressed during the development from mycelia to young fruiting bodies (T1), while only 1 gene (1 downregulated) was differentially expressed during the development from young fruiting bodies to mature fruiting bodies (T2). A total of 7 metabolites (3 increased and 4 reduced) were found to have changed in content during T1, and 4 metabolites (4 increased) were found to be different during T2. A conjoint analysis of the genome-wide connection network revealed that the metabolites that were more likely to be regulated were primarily in the post-squalene pathway.

Conclusions

This study provides useful information for understanding the regulation of ergosterol biosynthesis and the regulatory relationship between metabolites and genes in the ergosterol biosynthesis pathway during the development of fruiting bodies in F. velutipes.

Neutral Lipid Content in Lipid Droplets: Potential Biomarker of Cordycepin Accumulation in Cordycepin-Producing Fungi

To clarify the relationship between neutral lipid content and cordycepin accumulation in Cordyceps militaris, mutants were generated from mixed spores of two C. militaris strains with varying cordycepin-producing capacities. Fifteen stable mutants producing from 0.001 to 2.363 mg/mL cordycepin were finally selected. The relative fluorescence intensities of the 15 mutants, two C. militaris strains and an Aspergillus nidulans strain at different concentrations of lyophilized mycelium powder were then investigated using the Nile red method. The mutant CM1-1-1 with the highest relative fluorescence intensity among the eighteen strains was selected for optimizing the Nile red method. Relative fluorescence intensity was linearly correlated with cordycepin concentration in liquid broth (R² = 0.9514) and in lyophilized mycelium powder (R² = 0.9378) for the 18 cordycepin-producing strains under identical culture conditions and with cordycepin concentration in liquid broth (R² = 0.9727) and in lyophilized mycelium powder (R² = 0.9613) for CM1-1-1 under eight different sets of conditions. In addition, the cordycepin content in lyophilized mycelium powder measured by the Nile red method was linearly correlated with that determined by an HPLC method (R² = 0.9627). In conclusion, neutral lipids in lipid droplets are required during cordycepin accumulation; these neutral lipids are potential biomarkers of cordycepin biosynthesis.

Therapeutic Potential and Biological Applications of Cordycepin and Metabolic Mechanisms in Cordycepin-Producing Fungi

Cordycepin (3′-deoxyadenosine), a cytotoxic nucleoside analogue found in Cordyceps militaris, has attracted much attention due to its therapeutic potential and biological value. Cordycepin interacts with multiple medicinal targets associated with cancer, tumor, inflammation, oxidant, polyadenylation of mRNA, etc. The investigation of the medicinal drug actions supports the discovery of novel targets and the development of new drugs to enhance the therapeutic potency and reduce toxicity. Cordycepin may be of great value owing to its medicinal potential as an external drug, such as in cosmeceutical, traumatic, antalgic and muscle strain applications. In addition, the biological application of cordycepin, for example, as a ligand, has been used to uncover molecular structures. Notably, studies that investigated the metabolic mechanisms of cordycepin-producing fungi have yielded significant information related to the biosynthesis of high levels of cordycepin. Here, we summarized the medicinal targets, biological applications, cytotoxicity, delivery carriers, stability, and pros/cons of cordycepin in clinical applications, as well as described the metabolic mechanisms of cordycepin in cordycepin-producing fungi. We posit that new approaches, including single-cell analysis, have the potential to enhance medicinal potency and unravel all facets of metabolic mechanisms of cordycepin in Cordyceps militaris.

Alternative metabolic routes in channeling xylose to cordycepin production of Cordyceps militaris identified by comparative transcriptome analysis

The responsive mechanism of C. militaris TBRC7358 on xylose utilization was investigated by comparative analysis of transcriptomes, growth kinetics and cordycepin productions. The result showed that the culture grown on xylose exhibited high production yield of cordycepin on dry biomass. Comparing xylose to other carbon sources, a set of significantly up-regulated genes in xylose were enriched in pentose and glucuronate interconversion, and cordycepin biosynthesis. After validating up-regulated genes using quantitative real-time PCR, interestingly, putative alternative 3'-AMP-associated metabolic route on cordycepin biosynthesis was identified. Through reporter metabolites analysis of C. militaris, significant metabolites (e.g., AMP, glycine and L-glutamate) were identified guiding involvement of growth and cordycepin production. These findings suggested that there was a cooperative mechanism in transcriptional control of the supplying precursors pool directed towards the cordycepin biosynthesis through main and putative alternative metabolic routes for leverage of cell growth and cordycepin production on xylose of C. militaris strain TBRC7358.

An Integrated Transcriptome and Proteome Analysis Reveals Putative Regulators of Adventitious Root Formation in Taxodium 'Zhongshanshan'

Adventitious root (AR) formation from cuttings is the primary manner for the commercial vegetative propagation of trees. Cuttings is also the main method for the vegetative reproduction of Taxodium ‘Zhongshanshan’, while knowledge of the molecular mechanisms regulating the processes is limited. Here, we used mRNA sequencing and an isobaric tag for relative and absolute quantitation-based quantitative proteomic (iTRAQ) analysis to measure changes in gene and protein expression levels during AR formation in Taxodium ‘Zhongshanshan’. Three comparison groups were established to represent the three developmental stages in the AR formation process. At the transcript level, 4743 genes showed an expression difference in the comparison groups as detected by RNA sequencing. At the protein level, 4005 proteins differed in their relative abundance levels, as indicated by the quantitative proteomic analysis. A comparison of the transcriptome and proteome data revealed regulatory aspects of metabolism during AR formation and development. In summary, hormonal signal transduction is different at different developmental stages during AR formation. Other factors related to carbohydrate and energy metabolism and protein degradation and some transcription factor activity levels, were also correlated with AR formation. Studying the identified genes and proteins will provide further insights into the molecular mechanisms controlling AR formation.

Study of the whole genome, methylome and transcriptome of Cordyceps militaris

The complete genome of Cordyceps militaris was sequenced using single-molecule real-time (SMRT) sequencing technology at a coverage over 300×. The genome size was 32.57 Mb, and 14 contigs ranging from 0.35 to 4.58 Mb with an N50 of 2.86 Mb were assembled, including 4 contigs with telomeric sequences on both ends and an additional 8 contigs with telomeric sequences on either the 5′ or 3′ end. A methylome database of the genome was constructed using SMRT and m4C and m6A methylated nucleotides, and many unknown modification types were identified. The major m6A methylation motif is GA and GGAG, and the major m4C methylation motif is GC or CG/GC. In the C. militaris genome DNA, there were four types of methylated nucleotides that we confirmed using high-resolution LCMS-IT-TOF. Using PacBio Iso-Seq, a total of 31,133 complete cDNA sequences were obtained in the fruiting body. The conserved domains of the nontranscribed regions of the genome include TATA boxes, which are the initial regions of genome replication. There were 406 structural variants between the HN and CM01 strains, and there were 1,114 structural variants between the HN and ATCC strains.

Comparative transcriptome analysis reveals relationship of three major domesticated varieties of Auricularia auricula-judae

Auricularia auricula-judae is an edible mushroom and a traditional medicine in China as well as the fourth largest cultivated mushroom species in the world. Here for the first time, we present comparative transcriptome analyses of the fruiting bodies of three morphologically distinguishable A. auricula-judae cultivated varieties (Wujin, smooth; Banjin, partially wrinkled; and Quanjin, fully wrinkled) collected from Jilin Province, China. Biological triplicates were performed to determine the expression levels of 13,937 unigenes. Among them, only 13 unigenes were annotated to A. auricula-judae, highlighting the lack of publicly available reference sequences for this economically important species. Principal component analysis (PCA) determined that the gene expression profile of Quanjin was unique when compared to those of Banjin and Wujin. Such relationships were further supported by analyses of annotated and unannotated unigenes, differentially expressed unigenes, gene ontology functions, and the family of peroxidase genes. Using the KEGG database, significant alternations in biological pathways were detected among the three cultivars. This work contributes a large set of A. auricula-judae sequences to public database, establishes the relationships among major cultivars, and provides molecular guidance for breeding and cultivation.

Identification of cordycepin biosynthesis-related genes through de novo transcriptome assembly and analysis in Cordyceps cicadae

Cordyceps cicadae (Chanhua) is a parasitic fungus that grows on Cicada flammata larvae and is used to relieve exhaustion and treat numerous diseases, in part through its active constituent, cordycepin. We used de novo Illumina HiSeq 4000 sequencing to obtain transcriptomes of C. cicadae mycelium, fruiting body, and sclerotium, and identify differentially expressed genes. In the mycelium versus sclerotium libraries, 1576 upregulated and 2300 downregulated genes were identified. In the mycelium versus fruiting body and fruiting body versus sclerotium body libraries, 1604 and 1474 upregulated and 1365 and 1320 downregulated genes, respectively, were identified. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses identified 19 genes differentially expressed in mycelium versus fruiting body as related to the purine pathway, along with 28 and 16 genes differentially expressed in the mycelium versus sclerotium and fruiting body versus sclerotium groups, respectively. Gene expression of six key enzymes was validated by quantitative polymerase chain reaction. Specifically, 5'-nucleotidase (c62060g1) and adenosine deaminase (c35629g1) in purine nucleotide metabolism, which are involved in cordycepin biosynthesis, were significantly upregulated in the sclerotium group. These findings improved our understanding of genes involved in the biosynthesis of cordycepin and other characteristic secondary metabolites in C. cicadae.

Uncovering global metabolic response to cordycepin production in Cordyceps militaris through transcriptome and genome-scale network-driven analysis

The cellular metabolic adaptations of Cordyceps militaris have been progressively studied. In particular, the cordycepin pathway is of interest in medicinal applications. Even though the metabolic pathways for cordycepin production are known to be related to different carbon sources, the regulatory mechanisms at a systems level are poorly characterized. To explore the regulatory mechanisms, this study therefore aimed to investigate the global metabolic response to cordycepin production in C. militaris through transcriptome analysis and genome-scale network-driven analysis. Here, transcriptome analysis of 16,805 expressed genes in C. militaris strain TBRC6039 grown on different carbon sources was performed. Of these genes, 2,883 were significantly differentially expressed genes, uncovering sucrose- and glucose-mediated changes in the transcriptional regulation of central carbon metabolism in C. militaris, which was shown using the CmSNF1 mechanism as an example. After applying genome-scale metabolic network-driven analysis, reporter metabolites and key metabolic subnetworks involving adenosine, cordycepin and methionine were proposed through the up-regulation of cordycepin biosynthetic genes. Our findings suggest that the transcriptional regulation of these pathways is a ubiquitous feature in response to specific culture conditions during cordycepin overproduction.

The pharmacological properties of Ophiocordyceps xuefengensis revealed by transcriptome analysis

ETHNOPHARMACOLOGICAL RELEVANCE

The Yao ethnic group in Xuefeng Mountains area have used Xuefeng cordyceps, the caterpillar-fungus complex of Ophiocordyceps xuefengensis, for treating a variety of diseases for long. Just like some other cordyceps, O. xuefengensis, which is identified as the sister taxon of O. sinensis in 2013, also seems to have broad pharmacological properties, not only enhancing human immunity, anti-bacteria, anti-virus, but also anti-tumor. However, investigation of the medicinal fugal species O. xuefengensis can be found only in few literature records since its pharmacological and therapeutic use is mainly in traditional Yao communities by local healers.

AIM OF THE STUDY

The aim of this study is to collect samples of Xuefeng cordyceps and isolate the strain of O. xuefengensis, to determine bioactive components and evaluate the anti-tumor activity, to obtain the gene expression profile of O. xuefengensis and reveal its pharmacological properties by de novo transcriptome analysis. Accordingly, we attempt to provide information and give a comprehensive understanding of this mysterious medicinal fugal species from traditional Yao communities of China.

MATERIAL AND METHODS

Bioactive components were determined with HPLC-DAD-Q-TOF-MS technology; in vitro anti-tumor activity against 6 cell lines was evaluated using standard MTT assay; transcriptome analysis was done by de novo sequencing; unique genes were functionally profiled basing on Gene Ontology Database and the targeted genes were examined by blast.

RESULTS

Trace cordycepin, an anti-tumor agent, was detected in O. xuefengensis water extract. To some extent, the raw water extract of O. xuefengensis showed in vitro anti-tumor activity, against A549, HepG2, MCF-7, PC-3 and Raji cell lines. A total of 94,858 transcripts and 49,001 unique genes were obtained, amongst, 43.4% unique genes were matched with those of O. sinensis. Not all supposed genes related to cordycepin biosynthetic pathways were found by transcriptome analysis.

CONCLUSION

According to the gene expression profile, O. xuefengensis is very close to medicinal fungus O. sinensis. Raw water extract of O. xuefengensis, to a certain degree, could inhibit the growth of tumor cells, indicating that this fungus could be a new resource for the exploration of anti-tumor drug.

Comparative transcriptome analysis of dikaryotic mycelia and mature fruiting bodies in the edible mushroom Lentinula edodes

Lentinula edodes is a popular cultivated edible mushroom with high nutritional and medicinal value. To understand the regulation of gene expression in the dikaryotic mycelium and mature fruiting body in the commercially important Korean L. edodes strain, we first performed comparative transcriptomic analysis, using Illumina HiSeq platform. De novo assembly of these sequences revealed 11,675 representative transcripts in two different stages of L. edodes. A total of 9,092 unigenes were annotated and subjected to Gene Ontology, EuKaryotic Orthologous Groups, and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. Gene expression analysis revealed that 2,080 genes were differentially expressed, with 1,503 and 577 upregulated in the mycelium and a mature fruiting body, respectively. Analysis of 18 KEGG categories indicated that fruiting body-specific transcripts were significantly enriched in ‘replication and repair’ and ‘transcription’ pathways, which are important for premeiotic replication, karyogamy, and meiosis during maturation. We also searched for fruiting body-specific proteins such as aspartic protease, gamma-glutamyl transpeptidase, and cyclohexanone monooxygenase, which are involved in fruiting body maturation and isolation of functional substances. These transcriptomes will be useful in elucidating the molecular mechanisms of mature fruiting body development and beneficial properties, and contribute to the characterization of novel genes in L. edodes.

Whole Genome Sequence of an Edible and Potential Medicinal Fungus, Cordyceps guangdongensis

Cordyceps guangdongensis is an edible fungus which was approved as a novel food by the Chinese Ministry of Public Health in 2013. It also has a broad prospect of application in pharmaceutical industries, with many medicinal activities. In this study, the whole genome of C. guangdongensis GD15, a single spore isolate from a wild strain, was sequenced and assembled with Illumina and PacBio sequencing technology. The generated genome is 29.05 Mb in size, comprising nine scaffolds with an average GC content of 57.01%. It is predicted to contain a total of 9150 protein-coding genes. Sequence identification and comparative analysis indicated that the assembled scaffolds contained two complete chromosomes and four single-end chromosomes, showing a high level assembly. Gene annotation revealed a diversity of transposons that could contribute to the genome size and evolution. Besides, approximately 15.57% and 12.01% genes involved in metabolic processes were annotated by KEGG and COG respectively. Genes belonging to CAZymes accounted for 3.15% of the total genes. In addition, 435 transcription factors, involved in various biological processes, were identified. Among the identified transcription factors, the fungal transcription regulatory proteins (18.39%) and fungal-specific transcription factors (19.77%) represented the two largest classes of transcription factors. This genomic resource provided a new insight into better understanding the relevance of phenotypic characters and genetic mechanisms in C. guangdongensis.

Expression of a recombinant hybrid antimicrobial peptide magainin II-cecropin B in the mycelium of the medicinal fungus Cordyceps militaris and its validation in mice

Background

Antibiotic residues can cause antibiotic resistance in livestock and their food safety-related issues have increased the consumer demand for products lacking these residues. Hence, developing safe and effective antibiotic alternatives is important to the animal feed industry. With their strong antibacterial actions, antimicrobial peptides have potential as antibiotic alternatives.

Results

We investigated the antibacterial and immunomodulatory activities and the mechanisms of action of an antimicrobial peptide. The hybrid antimicrobial peptide magainin II-cecropin B (Mag II-CB) gene was transformed into the medicinal Cordyceps militaris fungus. Recombinant Mag II-CB exhibited broad-spectrum antibacterial activity in vitro and its antibacterial and immunomodulatory functions were evaluated in BALB/c mice infected with Escherichia coli (ATCC 25922). Histologically, Mag II-CB ameliorated E. coli-related intestinal damage and maintained the integrity of the intestinal mucosal barrier by up-regulating tight junction proteins (zonula occludens-1, claudin-1 and occludin). The intestinal microbial flora was positively modulated in the Mag II-CB-treated mice infected with E. coli. Mag II-CB treatment also supported immune functioning in the mice by regulating their plasma immunoglobulin and ileum secreted immunoglobulin A levels, by attenuating their pro-inflammatory cytokine levels, and by elevating their anti-inflammatory cytokines levels. Moreover, directly feeding the infected mice with the C. militaris mycelium producing Mag II-CB further proofed the antibacterial and immunomodulatory functions of recombinant hybrid antimicrobial peptide.

Conclusion

Our findings suggest that both purified recombinant AMPs and C. militaris mycelium producing AMPs display antibacterial and immunomodulatory activities in mice. And C. militaris producing AMPs has the potential to become a substitute to antibiotics as a feed additive for livestock in future.

Electronic supplementary material

The online version of this article (10.1186/s12934-018-0865-3) contains supplementary material, which is available to authorized users.

Insight into cordycepin biosynthesis of Cordyceps militaris: Comparison between a liquid surface culture and a submerged culture through transcriptomic analysis

Cordyceps militaris produces cordycepin, which is known to be a bioactive compound. Currently, cordycepin hyperproduction of C. militaris was carried out in a liquid surface culture because of its low productivity in a submerged culture, however the reason was not known. In this study, 4.92 g/L of cordycepin was produced at the 15^th day of C. militaris NBRC 103752 liquid surface culture, but only 1 mg/L was produced in the submerged culture. RNA-Seq was used to clarify the gene expression profiles of the cordycepin biosynthetic pathways of the submerged culture and the liquid surface culture. From this analysis, 1036 genes were shown to be upregulated and 557 genes were downregulated in the liquid surface culture compared with the submerged culture. Specifically, adenylosuccinate synthetase and phosphoribosylaminoimidazole-succinocarboxamide (SAICAR) synthase in purine nucleotide metabolism were significantly upregulated in the liquid surface culture. Thick mycelia formation in the liquid surface culture was found to induce the expression of hypoxia-related genes (GABA shunt, glutamate synthetase precursor, and succinate-semialdehyde dehydrogenase). Cytochrome P450 oxidoreductases containing heme were also found to be significantly enriched, suggesting that a hypoxic condition might be created in the liquid surface culture. These results suggest that hypoxic conditions are more suitable for cordycepin production in the liquid surface culture compared with the submerged culture. Our analysis paves the way for unraveling the cordycepin biosynthesis pathway and for improving cordycepin production in C. militaris.

Transcriptome-wide analysis reveals the progress of Cordyceps militaris subculture degeneration

BACKGROUND

The entomopathogenic mushroom Cordyceps militaris is an important medicinal and food resource owing to its various medicinal components and pharmacological effects. However, the high frequency of strain degeneration during subculture seriously restricts the large-scale production of C. militaris, and the mechanism underlying strain degeneration remains unclear. In this study, we artificially cultured C. militaris for six generations and compared changes during fruiting body growth. The transcriptome of six generations of C. militaris strains were sequenced with the Illumine Hiseq4000.

RESULTS

The subcultured C. militaris strains degenerated beginning at the third generation, with incomplete fruiting body growth beginning at the fourth generation. Over 9,015 unigenes and 731 new genes were identified. In addition, 35,323 alternative splicing (AS) events were detected in all samples, and more AS events occurred in the second, fourth and sixth generations. Compared with the first generation, the third generation (degenerated strain) included 2,498 differentially expressed genes (DEGs) including 1,729 up-regulated and 769 down-regulated genes. This number was higher than the number of DEGs in the second (1,892 DEGs), fourth (2,006 DEGs), fifth (2,273 DEGs) and sixth (2,188 DEGs) generations. Validation of RNA-seq by qRT-PCR showed that the expression patterns of 51 DEGs were in accordance with the transcriptome data.

CONCLUSION

Our results suggest that the mechanism of C. militaris strain degeneration is associated with gene involved in toxin biosynthesis, energy metabolism, and DNA methylation and chromosome remodeling.

Transcriptome analysis of Paecilomyces hepiali at different growth stages and culture additives to reveal putative genes in cordycepin biosynthesis

Paecilomyces hepialid (PH) is an endoparasitic fungus of Cordyceps sinensis (CS) and has become a substitute for CS due to their similar pharmacological activities. Because the market demand for CS is difficult to satisfy, and cordycepin, the effective compound of CS, is difficult to industrially produce, we produced 5 samples of PH by culturing for different durations and adding different additives to the culture broth, and detected their cordycepin content with UPLC ESI MS/MS. Then we grouped these cultures into five transcriptome comparisons containing 3 time variable groups and 2 additive variable groups. We used next-generation (NG) sequencing methods to acquire transcriptomic information and investigated the response of gene expression to the additives and the influence of different growth stages. This work will contribute to a better understanding of purine metabolism in PH, and possibly in other Cordyceps species. It will provide a useful resource to further advance transcriptomics studies in Cordyceps species.

A Lentinus edodes polysaccharide induces mitochondrial-mediated apoptosis in human cervical carcinoma HeLa cells

In this study, a homogeneous polysaccharide (LEP1) with an average molecular weight of 53kDa was successfully purified from the fruiting bodies of Lentinus edodes and its anticancer efficacy on human cervical carcinoma HeLa cells in vitro and associated possible molecular mechanism were also evaluated. MTT assay showed that LEP1 exhibited a dose-dependent inhibitory effect on the proliferation of HeLa cells and caused apoptotic death. Our present findings provided the first evidence that LEP1 induced the apoptosis of HeLa cells via a mitochondria dependent pathway, as indicated by an increase in Bax/Bcl-2 ratio, a loss of mitochondrial membrane potential (Δym), the release of cytochrome c from the mitochondria to the cytosol, activation of caspase-9 and caspase-3, and cleavage of poly (ADP-ribose) polymerase (PARP) in HeLa cells. These combined results unequivocally indicated that the involvement of mitochondria-mediated signaling pathway in LEP1-induced apoptosis and strongly provided experimental evidence for the use of LEP1 as a potential therapeutic agent in the prevention and treatment of human cervical carcinoma.

The transcriptional landscape of basidiosporogenesis in mature Pisolithus microcarpus basidiocarp

Background

Pisolithus microcarpus (Cooke & Massee) G. Cunn is a gasteromycete that produces closed basidiocarps in symbiosis with eucalypts and acacias. The fungus produces a complex basidiocarp composed of peridioles at different developmental stages and an upper layer of basidiospores free of the hyphae and ready for wind dispersal upon the rupture of the basidiocarp pellis. During basidiosporogenesis, a process that takes place inside the basidiocarp peridioles, a conspicuous reserve of fatty acids is present throughout development. While several previous studies have described basidiosporogenesis inside peridioles, very little is known about gene expression changes that may occur during this part of the fungal life cycle. The objective of this work was to analyze gene transcription during peridiole and basidiospore development, while focusing specifically on cell cycle progression and lipid metabolism.

Results

Throughout different developmental stages of the peridioles we analyzed, 737 genes were regulated between adjacent compartments (>5 fold, FDR-corrected p-value < 0.05) corresponding to 3.49% of the genes present in the P. microcarpus genome. We identified three clusters among the regulated genes which showed differential expression between the peridiole developmental stages and the basidiospores. During peridiole development, transcripts for proteins involved in cellular processes, signaling, and information storage were detected, notably those for coding transcription factors, DNA polymerase subunits, DNA repair proteins, and genes involved in chromatin structure. For both internal embedded basidiospores (hereto referred to as “Internal spores”, IS) and external free basidiospores (hereto referred to as “Free spores”, FS), upregulated transcripts were found to involve primary metabolism, particularly fatty acid metabolism (FA). High expression of transcripts related to β-oxidation and the glyoxylate shunt indicated that fatty acids served as a major carbon source for basidiosporogenesis.

Conclusion

Our results show that basidiocarp formation in P. microcarpus involves a complex array of genes that are regulated throughout peridiole development. We identified waves of transcripts with coordinated regulation and identified transcription factors which may play a role in this regulation. This is the first work to describe gene expression patterns during basidiocarp formation in an ectomycorrhizal gasteromycete fungus and sheds light on genes that may play important roles in the developmental process.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3545-5) contains supplementary material, which is available to authorized users.

Functional Analysis of Ribonucleotide Reductase from Cordyceps militaris Expressed in Escherichia coli

Cordyceps militaris produces cordycepin (3'-deoxyadenosine), which has various activities, including anti-oxidant, anti-tumoral, anti-viral, and anti-inflammatory. Ribonucleotide reductase (RNR) seems to be a candidate to produce cordycepin in C. militaris because RNR catalyzes the reduction of nucleotides to 2'-deoxynucleotides, whose structures are similar to that of cordycepin. However, the role of RNR has not been confirmed yet. In this study, complementary DNAs (cDNAs) of C. militaris RNR (CmRNR) large and small subunits (CmR1 and CmR2) were cloned from C. militaris NBRC9787 to investigate the function of CmRNR for its cordycepin production. C. militaris NBRC9787 began to produce cordycepin when grown in a liquid surface culture in medium composed of glucose and yeast extract for 15 days. CmR1 cDNA and CmR2 cDNA were obtained from its genomic DNA and from total RNA extracted from its mycelia after cultivation for 21 days, respectively. Recombinant CmR1 and CmR2 were expressed individually in Escherichia coli and purified. Purified recombinant CmR1 and CmR2 showed RNR activity toward adenosine diphosphate (ADP) only when two subunits were mixed but only show the reduction of ADP to 2'-deoxyADP. These results indicate that the pathway from ADP to 3'deoxyADP via CmRNR does not exist in C. militaris and cordycepin production in C. militaris may be mediated by other enzymes.

Comparative mitochondrial genomics toward exploring molecular markers in the medicinal fungus Cordyceps militaris

Cordyceps militaris is a fungus used for developing health food, but knowledge about its intraspecific differentiation is limited due to lack of efficient markers. Herein, we assembled the mitochondrial genomes of eight C. militaris strains and performed a comparative mitochondrial genomic analysis together with three previously reported mitochondrial genomes of the fungus. Sizes of the 11 mitochondrial genomes varied from 26.5 to 33.9 kb mainly due to variable intron contents (from two to eight introns per strain). Nucleotide variability varied according to different regions with non-coding regions showing higher variation frequency than coding regions. Recombination events were identified between some locus pairs but seemed not to contribute greatly to genetic variations of the fungus. Based on nucleotide diversity fluctuations across the alignment of all mitochondrial genomes, molecular markers with the potential to be used for future typing studies were determined.