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

Optimization of Fermentation Conditions of Cordyceps militaris and In Silico Analysis of Antifungal Property of Cordycepin Against Plant Pathogens

Cordyceps militaris, a medicinal fungus, has gained considerable attention owing to its potential health benefits, notably the production of bioactive compounds such as cordycepin. Cordycepin possesses significant antifungal, antibacterial, and antiviral properties. The present study focused on optimizing the fermentation conditions for C. militaris to boost the production of mycelia and cordycepin, alongside investigating its antifungal properties using in silico and in vitro approaches. The optimal conditions, yielding the highest cordycepin and mycelial biomass, were a temperature of 20°C and a pH range of 4-6, with glucose and sucrose as carbon sources and yeast extract and casein hydrolysate as nitrogen sources. Under these conditions, cordycepin production peaked at low pH (600-1000 mg/L) and with carbon and maltose (400-500 mg/L). The low temperature favored cordycepin production (400 mg/L), whereas casein hydrolysate as a nitrogen source boosted cordycepin yield (600 mg/L). The docking analysis indicated that cordycepin had the highest binding affinity for the tubulin beta chain 2 (-10.4 kcal/mol) compared to the fungicide tebuconazole (-7.9 kcal/mol for both targets). The in silico results were corroborated by in vitro studies, where the mycelial extract of C. militaris inhibited approximately 75% of fungal growth at a concentration of 6000 ppm. These findings suggest that optimizing fermentation conditions significantly enhances cordycepin production, and cordycepin shows antifungal solid activity, making it a promising agent for biocontrol in agriculture.

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.

Advances in biosynthesis and metabolic engineering strategies of cordycepin

Cordyceps militaris, also called as bei-chong-cao, is an insect-pathogenic fungus from the Ascomycota phylum and the Clavicipitaceae family. It is a valuable filamentous fungus with medicinal and edible properties that has been utilized in traditional Chinese medicine (TCM) and as a nutritious food. Cordycepin is the bioactive compound firstly isolated from C. militaris and has a variety of nutraceutical and health-promoting properties, making it widely employed in nutraceutical and pharmaceutical fields. Due to the low composition and paucity of wild resources, its availability from natural sources is limited. With the elucidation of the cordycepin biosynthetic pathway and the advent of synthetic biology, a green cordycepin biosynthesis in Saccharomyces cerevisiae and Metarhizium robertsii has been developed, indicating a potential sustainable production method of cordycepin. Given that, this review primarily focused on the metabolic engineering and heterologous biosynthesis strategies of cordycepin.

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.

Multifactorial interaction of selenium, iron, xylose, and glycine on cordycepin metabolism in Cordyceps militaris

Cordycepin, a nucleoside analog, is the main antioxidative and antimicrobial substance in Cordyceps militaris. To improve the metabolism of cordycepin, carbon sources, nitrogen sources, trace elements, and precursors were studied by single factor, Plackett-Burman, and central composite designs in C. militaris mycelial fermentation. Under the regulation of the multifactorial interactions of selenite, ferrous chloride, xylose, and glycine, cordycepin production was increased by 5.2-fold compared with the control. The gene expression of hexokinase, ATP phosphoribosyltransferase, adenylosuccinate synthetase, and cns1-3 in the glycolysis, pentose phosphate, and adenosine synthesis pathways were increased by 3.2-7.5 times due to multifactorial interactions, while the gene expression of histidine biosynthesis trifunctional protein and histidinol-phosphate aminotransferase in histidine synthesis pathway were decreased by 23.4%-56.2%. Increasing with cordycepin production, glucose uptake was accelerated, mycelia growth was inhibited, and the cell wall was damaged. Selenomethionine (SeMet), selenocysteine (SeCys), and selenium nanoparticles (SeNPs) were the major Se species in C. militaris mycelia. This study provides a new insight for promoting cordycepin production by regulating glycolysis, pentose phosphate, and histidine metabolism. KEY POINTS: • Cordycepin production in the CCDmax group was 5.2-fold than that of the control. • Glucose uptake of the CCDmax group was accelerated and cell wall was damaged. • The metabolic flux was concentrated to the cordycepin synthesis pathway.

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.

Functional insight into Cordyceps militaris sugar transporters by structure modeling, network analysis and allosteric regulation

Insights into the structures, functions and dynamics of Cordyceps militaris (C. militaris) sugar transporters are necessary for understanding their versatile metabolic capability for fungal growth. The sequence-function relationship study of 85 C. militaris sugar transporters showed that there is a gap between phylogenetic-based subfamily classification and their functions. Beyond protein sequences, structural modeling and principal component analysis of the structural ensemble revealed the different folds of the Car and Org subfamilies. Performing channel detection and network analysis found that the Alp and Hex subfamilies can be specifically distinguished from others by the betweenness of channel residues. Signature dynamics analysis further suggested that the Hex subfamily demonstrates different dynamics, with high flexibility at the H1 region in TM11. Furthermore, the H1 region as an allosteric site was examined by network parameter calculations that guided allosteric pathways between this region and the channel cavity. Together with gene expression data of C. militaris, e.g., Hex06741 in the Hex subfamily, it was promisingly expressed when sugar utilization was altered. This work demonstrates an in silico framework for investigating C. militaris sugar transporters as an example case study of the allosteric activity of the Hex subfamily and can facilitate sugar transporter engineering design that can further optimize the preferable sugar utilization and fermentation process of C. militaris.

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.

Transcriptome and metabolome profiling unveils the mechanisms of naphthalene acetic acid in promoting cordycepin synthesis in Cordyceps militaris

Cordycepin, an important active substance in Cordyceps militaris, possesses antiviral and other beneficial activities. In addition, it has been reported to effectively promote the comprehensive treatment of COVID-19 and thus has become a research hotspot. The addition of naphthalene acetic acid (NAA) is known to significantly improve the yield of cordycepin; however, its related molecular mechanism remains unclear. We conducted a preliminary study on C. militaris with different concentrations of NAA. We found that treatment with different concentrations of NAA inhibited the growth of C. militaris, and an increase in its concentration significantly improved the cordycepin content. In addition, we conducted a transcriptome and metabolomics association analysis on C. militaris treated with NAA to understand the relevant metabolic pathway of cordycepin synthesis under NAA treatment and elucidate the relevant regulatory network of cordycepin synthesis. Weighted gene co-expression network analysis (WGCNA), transcriptome, and metabolome association analysis revealed that genes and metabolites encoding cordycepin synthesis in the purine metabolic pathway varied significantly with the concentration of NAA. Finally, we proposed a metabolic pathway by analyzing the relationship between gene-gene and gene-metabolite regulatory networks, including the interaction of cordycepin synthesis key genes; key metabolites; purine metabolism; TCA cycle; pentose phosphate pathway; alanine, aspartate, and glutamate metabolism; and histidine metabolism. In addition, we found the ABC transporter pathway to be significantly enriched. The ABC transporters are known to transport numerous amino acids, such as L-glutamate, and participate in the amino acid metabolism that affects the synthesis of cordycepin. Altogether, multiple channels work together to double the cordycepin yield, thereby providing an important reference for the molecular network relationship between the transcription and metabolism of cordycepin synthesis.

Effects of Exogenous Lanthanum Nitrate on the Active Substance Content and Antioxidant Activity of Caterpillar Medicinal Mushroom Cordyceps militaris (Ascomycetes)

Cordyceps militaris is a medicinal and edible mushroom. Researchers often add exogenous substances to the culture medium to increase the active substance content in C. militaris. However, the effect of earth elements on the active substance content in C. militaris and its antioxidant effects have not been reported. In this study, the active substance content in C. militaris treated with lanthanum nitrate was determined using high-performance liquid chromatography and ultraviolet spectrophotometry, and the effect on the antioxidant capacity of C. militaris after lanthanum nitrate spraying was further explored. The results showed that, in the experimental concentration range, the two concentrations of 10 mg/L and 50 mg/L had a significant influence on the active substance content of C. militaris. When the concentration of lanthanum nitrate was 10 mg/L, the synthesis of pentostatin and cordycepin was promoted. When the concentration of lanthanum nitrate was 50 mg/L, it significantly promoted the synthesis of cordycepin, and the ferric-reducing power and DPPH· scavenging rate of C. militaris treated at this concentration were significantly higher than those of the control group. However, lanthanum nitrate had no significant effect on ergosterol synthesis (P > 0.05). Finally, considering that the residual amount of lanthanum in C. militaris and the residual amount of lanthanum in 50 mg/L lanthanum nitrate-treated C. militaris is within the allowable daily intake of 4.2 mg for humans, the optimal concentration of lanthanum nitrate-treated C. militaris is 50 mg/L.

Enhancing Genome-Scale Model by Integrative Exometabolome and Transcriptome: Unveiling Carbon Assimilation towards Sphingolipid Biosynthetic Capability of Cordyceps militaris

Cordyceps militaris is an industrially important fungus, which is often used in Asia as traditional medicine. There has been a published genome-scale metabolic model (GSMM) of C. militaris useful for predicting its growth behaviors; however, lipid metabolism, which plays a vital role in cellular functions, remains incomplete in the GSMM of C. militaris. A comprehensive study on C. militaris was thus performed by enhancing GSMM through integrative analysis of metabolic footprint and transcriptome data. Through the enhanced GSMM of C. militaris (called iPC1469), it contained 1469 genes, 1904 metabolic reactions and 1229 metabolites. After model evaluation, in silico growth simulation results agreed well with the experimental data of the fungal growths on different carbon sources. Beyond the model-driven integrative data analysis, interestingly, we found key metabolic responses in alteration of lipid metabolism in C. militaris upon different carbon sources. The sphingoid bases (e.g., sphinganine, sphingosine, and phytosphingosine) and ceramide were statistically significant accumulated in the xylose culture when compared with other cultures; this study suggests that the sphingolipid biosynthetic capability in C. militaris was dependent on the carbon source assimilated for cell growth; this finding provides a comprehensive basis for the sphingolipid biosynthesis in C. militaris that can help to further redesign its metabolic control for medicinal and functional food applications.

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.

Construction of Light-Responsive Gene Regulatory Network for Growth, Development and Secondary Metabolite Production in Cordyceps militaris

Cordyceps militaris is an edible fungus that produces many beneficial compounds, including cordycepin and carotenoid. In many fungi, growth, development and secondary metabolite production are controlled by crosstalk between light-signaling pathways and other regulatory cascades. However, little is known about the gene regulation upon light exposure in C. militaris. This study aims to construct a gene regulatory network (GRN) that responds to light in C. militaris. First, a genome-scale GRN was built based on transcription factor (TF)-target gene interactions predicted from the Regulatory Sequence Analysis Tools (RSAT). Then, a light-responsive GRN was extracted by integrating the transcriptomic data onto the genome-scale GRN. The light-responsive network contains 2689 genes and 6837 interactions. From the network, five TFs, Snf21 (CCM_04586), an AT-hook DNA-binding motif TF (CCM_08536), a homeobox TF (CCM_07504), a forkhead box protein L2 (CCM_02646) and a heat shock factor Hsf1 (CCM_05142), were identified as key regulators that co-regulate a large group of growth and developmental genes. The identified regulatory network and expression profiles from our analysis suggested how light may induce the growth and development of C. militaris into a sexual cycle. The light-mediated regulation also couples fungal development with cordycepin and carotenoid production. This study leads to an enhanced understanding of the light-responsive regulation of growth, development and secondary metabolite production in the fungi.

Dynamic content changes of cordycepin and adenosine and transcriptome in Cordyceps kyushuensis Kob at different fermentation stages

20% (w/w) Astragali radix was added to the rice medium to cultivate C. kyushuensis Kob. The fermentation product was collected at mycelium stage, coloring stage, stromata-forming initial stage and fruiting body stage of C. kyushuensis Kob. The dynamic content changes of cordycepin and adenosine were detected at different fermentation stages. In the rice medium with Astragalus radix, both cordycepin and adenosine reached the highest content value on the 30th day of fermentation, 17.31 mg/g and 0.94 mg/g, respectively, which were 8.6 times and 2.0 times of that in rice medium at the same stage. At the same time, transcriptomics technology was used to analyze C. kyushuensis Kob during these four periods.

Metabolic Responses of Carotenoid and Cordycepin Biosynthetic Pathways in Cordyceps militaris under Light-Programming Exposure through Genome-Wide Transcriptional Analysis

Cordyceps militaris is currently exploited for commercial production of specialty products as its biomass constituents are enriched in bioactive compounds, such as cordycepin. The rational process development is important for economically feasible production of high quality bioproducts. Light is an abiotic factor affecting the cultivation process of this entomopathogenic fungus, particularly in its carotenoid formation. To uncover the cell response to light exposure, this study aimed to systematically investigate the metabolic responses of C. militaris strain TBRC6039 using integrative genome-wide transcriptome and genome-scale metabolic network (GSMN)-driven analysis. The genome-wide transcriptome analysis showed 8747 expressed genes in the glucose and sucrose cultures grown under light-programming and dark conditions. Of them, 689 differentially expressed genes were significant in response to the light-programming exposure. Through integration with the GSMN-driven analysis using the improved network (iRT1467), the reporter metabolites, e.g., adenosine-5'-monophosphate (AMP) and 2-oxoglutarate, were identified when cultivated under the carotenoid-producing condition controlled by light-programming exposure, linking to up-regulations of the metabolic genes involved in glyoxalase system, as well as cordycepin and carotenoid biosynthesis. These results indicated that C. militaris had a metabolic control in acclimatization to light exposure through transcriptional co-regulation, which supported the cell growth and cordycepin production in addition to the accumulation of carotenoid as a photo-protective bio-pigment. This study provides a perspective in manipulating the metabolic fluxes towards the target metabolites through either genetic or physiological approaches.

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.

Synthesis of cordycepin: Current scenario and future perspectives

Cordyceps genus, such as C. militaris and C. kyushuensis, is a source of a rare traditional Chinese medicine that has been used for the treatment of numerous chronic and malignant diseases. Cordycepin, 3'-deoxyadenosine, is a major active compound found in most Cordyceps. Cordycepin exhibits a variety of biological activities, including anti-tumor, immunomodulation, antioxidant, and anti-aging, among others, which could be applied in health products, medicine, cosmeceutical etc. fields. This review focuses on the synthesis methods for cordycepin. The current methods for cordycepin synthesis involve chemical synthesis, microbial fermentation, in vitro synthesis and biosynthesis; however, some defects are unavoidable and the production is still far from the demand of cordycepin. For the future study of cordycepin synthesis, based on the illumination of cordycepin biosynthesis pathway, genetical engineering of the Cordyceps strain or introducing microbes by virtue of synthetic biology will be the great potential strategies for cordycepin synthesis. This review will aid the future synthesis of the valuable cordycepin.

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.

Probing Carbon Utilization of Cordyceps militaris by Sugar Transportome and Protein Structural Analysis

Beyond comparative genomics, we identified 85 sugar transporter genes in Cordyceps militaris, clustering into nine subfamilies as sequence- and phylogenetic-based functional classification, presuming the versatile capability of the fungal growths on a range of sugars. Further analysis of the global gene expression patterns of C. militaris showed 123 genes were significantly expressed across the sucrose, glucose, and xylose cultures. The sugar transporters specific for pentose were then identified by gene-set enrichment analysis. Of them, the putative pentose transporter, CCM_06358 gene, was highest expressed in the xylose culture, and its functional role in xylose transport was discovered by the analysis of conserved structural motifs. In addition, a battery of molecular modeling methods, including homology modeling, transport pathway analysis, residue interaction network combined with molecular mechanics Poisson–Boltzmann surface area simulation (MM-PBSA), was implemented for probing the structure and function of the selected pentose transporter (CCM_06358) as a representative of sugar transportome in C. militaris. Considering the network bottlenecks and structural organizations, we further identified key amino acids (Phe38 and Trp441) and their interactions with other residues, contributing the xylose transport function, as verified by binding free energy calculation. The strategy used herein generated remarkably valuable biological information, which is applicable for the study of sugar transportome and the structure engineering of targeted transporter proteins that might link to the production of bioactive compounds derived from xylose metabolism, such as cordycepin.

Optimizing cultivation of Cordyceps militaris for fast growth and cordycepin overproduction using rational design of synthetic media

Cordyceps militaris is an entomopathogenic fungus which is often used in Asia as a traditional medicine developed from age-old wisdom. Presently, cordycepin from C. militaris is a great interest in medicinal applications. However, cellular growth of C. militaris and the association with cordycepin production remain poorly understood. To explore the metabolism of C. militaris as potential cell factories in medical and biotechnology applications, this study developed a high-quality genome-scale metabolic model of C. militaris, iNR1329, based on its genomic content and physiological data. The model included a total of 1329 genes, 1821 biochemical reactions, and 1171 metabolites among 4 different cellular compartments. Its in silico growth simulation results agreed well with experimental data on different carbon sources. iNR1329 was further used for optimizing the growth and cordycepin overproduction using a novel approach, POPCORN, for rational design of synthetic media. In addition to the high-quality GEM iNR1329, the presented POPCORN approach was successfully used to rationally design an optimal synthetic medium with C:N ratio of 8:1 for enhancing 3.5-fold increase in cordycepin production. This study thus provides a novel insight into C. militaris physiology and highlights a potential GEM-driven method for synthetic media design and metabolic engineering application. The iNR1329 and the POPCORN approach are available at the GitHub repository: https://github.com/sysbiomics/Cordyceps_militaris-GEM.

Improved Cordycepin Production by Cordyceps militaris KYL05 Using Casein Hydrolysate in Submerged Conditions

Cordycepin, a beneficial bioactive product specifically found in Cordyceps, has received attention in various bioindustrial applications such as in pharmaceuticals, functional foods, and cosmetics, due to its significant functions. However, low productivity of cordycepin is a barrier to commercialization. In this study, Cordyceps militaris was mutated by UV irradiation to improve the cordycepin production. The highest producer KYL05 strain was finally selected and its cordycepin production was increased about 1.5-fold compared to wild type. In addition, the effects of culture conditions were fundamentally investigated. Optimal conditions were as follows: pH 6, temperature of 25 °C, shaking speed of 150 rpm, and culture time of 6 days. Effects of medium component on cordycepin production were also investigated by using various carbon and nitrogen sources. It was found that glucose and casein hydrolysate (CH) were most effective as carbon and nitrogen sources in cordycepin production (2.3-fold improvement) with maximum cordycepin production of about 445 mg/L. In particular, production was significantly affected by CH. These results should be of value in improving the efficiency of mass production of cordycepin.

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.