The blue-light receptor CmWC-1 mediates fruit body development and secondary metabolism in Cordyceps militaris

Cordyceps militaris: A novel mushroom platform for metabolic engineering

Cordyceps militaris, widely recognized as a medicinal and edible mushroom in East Asia, contains a variety of bioactive compounds, including cordycepin (COR), pentostatin (PTN) and other high-value compounds. This review explores the potential of developing C. militaris as a cell factory for the production of high-value chemicals and nutrients. This review comprehensively summarizes the fermentation advantages, metabolic networks, expression elements, and genome editing tools specific to C. militaris and discusses the challenges and barriers to further research on C. militaris across various fields, including computational biology, existing DNA elements, and genome editing approaches. This review aims to describe specific and promising opportunities for the in-depth study and development of C. militaris as a new chassis cell. Additionally, to increase the practicability of this review, examples of the construction of cell factories are provided, and promising strategies for synthetic biology development are illustrated.

Advancing Cordyceps militaris Industry: Gene Manipulation and Sustainable Biotechnological Strategies

Cordyceps militaris is considered to be of great medicinal potential due to its remarkable pharmacological effects, safety, and edible characteristics. With the completion of the genome sequence and the advancement of efficient gene-editing technologies, coupled with the identification of gene functions in Cordyceps militaris, this fungus is poised to emerge as an outstanding strain for medicinal engineering applications. This review focuses on the development and application of genomic editing techniques, including Agrobacterium tumefaciens-mediated transformation (ATMT), PEG-mediated protoplast transformation (PMT), and CRISPR/Cas9. Through the application of these techniques, researchers can engineer the biosynthetic pathways of valuable secondary metabolites to boost yields; such metabolites include cordycepin, polysaccharides, and ergothioneine. Furthermore, by identifying and modifying genes that influence the growth, disease resistance, and tolerance to environmental stress in Cordyceps militaris, it is possible to stimulate growth, enhance desirable traits, and increase resilience to unfavorable conditions. Finally, the green sustainable industrial development of C. militaris using agricultural waste to produce high-value-added products and the future research directions of C. militaris were discussed. This review will provide future directions for the large-scale production of bioactive ingredients, molecular breeding, and sustainable development of C. militaris.

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.

Construction of a New Agrobacterium tumefaciens-Mediated Transformation System based on a Dual Auxotrophic Approach in Cordyceps militaris

Cordyceps militaris is a significant edible fungus that produces a variety of bioactive compounds. We have previously established a uridine/uracil auxotrophic mutant and a corresponding Agrobacterium tumefaciens-mediated transformation (ATMT) system for genetic characterization in C. militaris using pyrG as a screening marker. In this study, we constructed an ATMT system based on a dual pyrG and hisB auxotrophic mutant of C. militaris. Using the uridine/uracil auxotrophic mutant as the background and pyrG as a selection marker, the hisB gene encoding imidazole glycerophosphate dehydratase, required for histidine biosynthesis, was knocked out by homologous recombination to construct a histidine auxotrophic C. militaris mutant. Then, pyrG in the histidine auxotrophic mutant was deleted to construct a ΔpyrG ΔhisB dual auxotrophic mutant. Further, we established an ATMT transformation system based on the dual auxotrophic C. militaris by using GFP and DsRed as reporter genes. Finally, to demonstrate the application of this dual transformation system for studies of gene function, knock out and complementation of the photoreceptor gene CmWC-1 in the dual auxotrophic C. militaris were performed. The newly constructed ATMT system with histidine and uridine/uracil auxotrophic markers provides a promising tool for genetic modifications in the medicinal fungus C. militaris.

Genomic and Transcriptome Analysis Reveals the Biosynthesis Network of Cordycepin in Cordyceps militaris

Cordycepin is the primary active compound of Cordyceps militaris. However, the definitive genetic mechanism governing cordycepin synthesis in fruiting body growth and development remains elusive, necessitating further investigation. This study consists of 64 C. militaris strains collected from northeast China. The high-yielding cordycepin strain CMS19 was selected for the analysis of cordycepin production and the genetic basis of cordycepin anabolism. First, the whole-genome sequencing of CMS19 yielded a final size of 30.96 Mb with 8 contigs and 9781 protein-coding genes. The genome component revealed the presence of four additional secondary metabolite gene clusters compared with other published genomes, suggesting the potential for the production of new natural products. The analyses of evolutionary and genetic differentiation revealed a close relationship between C. militaris and Beauveria bassiana. The population of strains distributed in northeast China exhibited the significant genetic variation. Finally, functional genes associated with cordycepin synthesis were identified using a combination of genomic and transcriptomic analyses. A large number of functional genes associated with energy and purine metabolism were significantly enriched, facilitating the reconstruction of a hypothetical cordycepin metabolic pathway. Therefore, our speculation of the cordycepin metabolism pathway involved 24 genes initiating from the glycolysis and pentose phosphate pathways, progressing through purine metabolism, and culminating in the core region of cordycepin synthesis. These findings could offer fundamental support for scientific utilizations of C. militaris germplasm resources and standardized cultivation for cordycepin production.

Light-Exposed Metabolic Responses of Cordyceps militaris through Transcriptome-Integrated Genome-Scale Modeling

The genome-scale metabolic model (GSMM) of Cordyceps militaris provides a comprehensive basis of carbon assimilation for cell growth and metabolite production. However, the model with a simple mass balance concept shows limited capability to probe the metabolic responses of C. militaris under light exposure. This study, therefore, employed the transcriptome-integrated GSMM approach to extend the investigation of C. militaris's metabolism under light conditions. Through the gene inactivity moderated by metabolism and expression (GIMME) framework, the iPS1474-tiGSMM model was furnished with the transcriptome data, thus providing a simulation that described reasonably well the metabolic responses underlying the phenotypic observation of C. militaris under the particular light conditions. The iPS1474-tiGSMM obviously showed an improved prediction of metabolic fluxes in correlation with the expressed genes involved in the cordycepin and carotenoid biosynthetic pathways under the sucrose culturing conditions. Further analysis of reporter metabolites suggested that the central carbon, purine, and fatty acid metabolisms towards carotenoid biosynthesis were the predominant metabolic processes responsible in light conditions. This finding highlights the key responsive processes enabling the acclimatization of C. militaris metabolism in varying light conditions. This study provides a valuable perspective on manipulating metabolic genes and fluxes towards the target metabolite production of C. militaris.

Transcriptome Analysis Reveals Mycelial and Fruiting Responses to Lithium Chloride in Coprinopsis cinerea

Lithium chloride (LiCl) has been used in signalling and molecular studies of animals, plants, and yeast. However, information on its roles in basidiomycetous fungi is still limited. In this study, we used RNA-Seq to study the effects of LiCl on Coprinopsis cinerea. LiCl enhanced mycelial growth and inhibited fruiting body formation in C. cinerea. RNA-Seq of the LiCl-treated C. cinerea resulted in a total of 14,128 genes. There were 1199 differentially expressed genes (DEGs) between the LiCl-treated samples and control samples in the mycelium stage (the first time point), and 1391 DEGs were detected when the control samples were forming hyphal knots while the treated samples were still in the mycelium (the second time point). Pathway enrichment analysis of the DEGs revealed a significant association between enhanced mycelium growth in the LiCl-treated C. cinerea and metabolic pathways. In addition, the DEGs involved in cellular process pathways, including "cell cycle-yeast" and "meiosis-yeast", were identified in suppressed C. cinerea fruiting body formation by LiCl under favourable environmental conditions. As LiCl can predominantly inhibit the activity of glycogen synthase kinase3 (GSK3), our findings suggest that LiCl affects the expression of genes involved in fruiting body initiation and cellular processes by inhibiting GSK3 activity which is essential for fruiting body formation.

Transcription factors: switches for regulating growth and development in macrofungi

Macrofungi (or mushrooms) act as an extraordinarily important part to human health due to their nutritional and/or medicinal value, but the detailed researches in growth and development mechanisms have yet to be explored further. Transcription factors (TFs) play indispensable roles in signal transduction and affect growth, development, and metabolism of macrofungi. In recent years, increasing research effort has been employed to probe the relationship between the development of macrofungi and TFs. Herein, the present review comprehensively summarized the functional TFs researched in macrofungi, including modulating mycelial growth, fructification, sclerotial formation, sexual reproduction, spore formation, and secondary metabolism. Meanwhile, the possible effect mechanisms of TFs on the growth and development of some macrofungi were also revealed. Specific examples of functional characterizations of TFs in macrofungi (such as Schizophyllum commune and Coprinopsis cinerea) were described to a better comprehension of regulatory effect. Future research prospects in the field of TFs of macrofungi are discussed. We illustrated the functional versatility of the TFs in macrofungi based on specific examples. A systematical realization of the interaction and possible mechanisms between TFs and macrofungi can supply possible solutions to regulate genetic characteristics, which supply novel insights into the regulation of growth, development and metabolism of macrofungi. KEY POINTS: • The functional TFs researched in macrofungi were summarized. • The possible effect mechanisms of TFs in macrofungal were described. • The multiple physiological functions of TFs in macrofungi were discussed.

Hydrophobin Gene Cmhyd4 Negatively Regulates Fruiting Body Development in Edible Fungi Cordyceps militaris

A deep understanding of the mechanism of fruiting body development is important for mushroom breeding and cultivation. Hydrophobins, small proteins exclusively secreted by fungi, have been proven to regulate the fruiting body development in many macro fungi. In this study, the hydrophobin gene Cmhyd4 was revealed to negatively regulate the fruiting body development in Cordyceps militaris, a famous edible and medicinal mushroom. Neither the overexpression nor the deletion of Cmhyd4 affected the mycelial growth rate, the hydrophobicity of the mycelia and conidia, or the conidial virulence on silkworm pupae. There was also no difference between the micromorphology of the hyphae and conidia in WT and ΔCmhyd4 strains observed by SEM. However, the ΔCmhyd4 strain showed thicker aerial mycelia in darkness and quicker growth rates under abiotic stress than the WT strain. The deletion of Cmhyd4 could promote conidia production and increase the contents of carotenoid and adenosine. The biological efficiency of the fruiting body was remarkably increased in the ΔCmhyd4 strain compared with the WT strain by improving the fruiting body density, not the height. It was indicated that Cmhyd4 played a negative role in fruiting body development. These results revealed that the diverse negative roles and regulatory effects of Cmhyd4 were totally different from those of Cmhyd1 in C. militaris and provided insights into the developmental regulatory mechanism of C. militaris and candidate genes for C. militaris strain breeding.

Role of the blue light receptor gene Icwc-1 in mycelium growth and fruiting body formation of Isaria cicadae

The Isaria cicadae, is well known highly prized medicinal mushroom with great demand in food and pharmaceutical industry. Due to its economic value and therapeutic uses, natural sources of wild I. cicadae are over-exploited and reducing continuously. Therefore, commercial cultivation in controlled environment is an utmost requirement to fulfill the consumer's demand. Due to the lack of knowledge on fruiting body (synnemata) development and regulation, commercial cultivation is currently in a difficult situation. In the growth cycle of macrofungi, such as mushrooms, light is the main factor affecting growth and development, but so far, specific effects of light on the growth and development of I. cicadae is unknown. In this study, we identified a blue light receptor white-collar-1 (Icwc-1) gene homologue with well-defined functions in morphological development in I. cicadae based on gene knockout technology and transcriptomic analysis. It was found that the Icwc-1 gene significantly affected hyphal growth and fruiting body development. This study confirms that Icwc-1 acts as an upstream regulatory gene that regulates genes associated with fruiting body formation, pigment-forming genes, and related genes for enzyme synthesis. Transcriptome data analysis also found that Icwc-1 affects many important metabolic pathways of I. cicadae, i.e., amino acid metabolism and fatty acid metabolism. The above findings will not only provide a comprehensive understanding about the molecular mechanism of light regulation in I. cicadae, but also provide new insights for future breeding program and improving this functional food production.

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.

Genome resequencing and transcriptome analysis reveal the molecular mechanism of albinism in Cordyceps militaris

Light is an important regulator of most fungal life activities and transmits signals through certain photoreceptor proteins such as phytochromes and cryptochromes. However, the light response mechanism varies across different fungi. The WCC complex composed of white collar-1 (WC-1) and white collar-2 (WC-2) is considered to be the key factor regulating fungal albinism. The photoreceptor protein Vivid (VVD) is the negative regulator of the WCC complex. In this study, we discovered an albino mutant (Alb) generated by ⁶⁰Co-γ-ray irradiation from Cordyceps militaris (C. militaris). This mutant showed albinism of the mycelia and fruiting bodies under light, and the fruiting bodies developed normally. However, this phenotype in Alb differed from that in the CmWC-1 mutant. This suggests that CmWC1 may not be mutated in Alb. A mutated polyketide synthase (CmPKS) was found through genome resequencing analysis. CmPKS was significantly induced by a light signal, and its mutation reduced melanin accumulation in C. militaris. In addition, we found that a zinc-finger domain-containing protein (CmWC-3) was induced by a light signal and interacted with CmWC-1 and CmVVD. Moreover, CmWC-2 also interacted with CmWC-1 to form the WCC complex and was inhibited by CmVVD. In addition, CmWC-3 directly bound with the CmPKS promoter, but CmWC1 did not. These results suggest that albinism and fruiting body development are two independent processes; the WCC complex of CmWC-1 with CmWC-3 regulates CmPKS expression to regulate color change, whereas CmWC-1 with CmWC-2 affects fruiting body development via the carotenoid pathway. These findings will help us to better understand the albinism mechanism of C. militaris.

Application of omics technology in the research on edible fungi

Edible fungus is a large fungus distributed all over the world and used as food and medicine. But people's understanding of edible fungi is not as much as that of ordinary crops, so people have started a number of research on edible fungi in recent years. With the development of science and technology, omics technology has gradually walked into people's vision. Omics technology has high sensitivity and wide application range, which is favored by researchers. The application of omics technology to edible fungus research is a major breakthrough, which has transferred edible fungus research from artificial cultivation to basic research. Now omics technology in edible fungi has been flexibly combined with other research methods, involving multiple studies of edible fungus, such as genetic breeding, growth and development, stress resistance, and the use of special components in edible fungus as pharmaceutical additives. It is believed that in the future, the research of edible fungi will also be brought to a deeper level with the help of omics technology. This paper introduces the application progress of modern omics technology to the study on edible fungi and mentions the application prospect of edible fungi research with the constant development of omics technology, thereby providing ideas for the follow-up in-depth research on edible fungi.

The Role of Chromatin and Transcriptional Control in the Formation of Sexual Fruiting Bodies in Fungi

Fungal fruiting bodies are complex, three-dimensional structures that arise from a less complex vegetative mycelium. Their formation requires the coordinated action of many genes and their gene products, and fruiting body formation is accompanied by major changes in the transcriptome. In recent years, numerous transcription factor genes as well as chromatin modifier genes that play a role in fruiting body morphogenesis were identified, and through research on several model organisms, the underlying regulatory networks that integrate chromatin structure, gene expression, and cell differentiation are becoming clearer. This review gives a summary of the current state of research on the role of transcriptional control and chromatin structure in fruiting body development. In the first part, insights from transcriptomics analyses are described, with a focus on comparative transcriptomics. In the second part, examples of more detailed functional characterizations of the role of chromatin modifiers and/or transcription factors in several model organisms (Neurospora crassa, Aspergillus nidulans, Sordaria macrospora, Coprinopsis cinerea, and Schizophyllum commune) that have led to a better understanding of regulatory networks at the level of chromatin structure and transcription are discussed.

Cmcrf1, a Putative Zn2Cys6 Fungal Transcription Factor, Is Involved in Conidiation, Carotenoid Production, and Fruiting Body Development in Cordyceps militaris

Cordyceps militaris is a high-value medicinal and edible fungus that produces many bioactive compounds, including carotenoid, and thus, improving the carotenoid productivity of C. militaris will increase its commercial value. However, little is known about the genetic regulatory mechanism of carotenoid biosynthesis in C. militaris. To further understanding the regulatory mechanism of carotenoid biosynthesis, we performed a large-scale screen of T-DNA insertional mutant library and identified a defective mutant, denoted T111, whose colonies did not change color from white to yellow upon exposure to light. Mutation analysis confirmed that a single T-DNA insertion occurred in the gene encoding a 695-amino-acid putative fungal-specific transcription factor with a predicted Zn2Cys6 binuclear cluster DNA-binding domain found uniquely in fungi. Targeted deletion of this gene, denoted C. militaris carotenogenesis regulatory factor 1 (Cmcrf1), generated the ΔCmcrf1 mutant that exhibited drastically reduced carotenoid biosynthesis and failed to generate fruiting bodies. In addition, the ΔCmcrf1 mutant showed significantly increased conidiation and increased hypersensitivity to cell-wall-perturbing agents compared with the wild-type strain. However, the Cmcrf1 gene did not have an impact on the mycelia growth of C. militaris. These results show that Cmcrf1 is involved in carotenoid biosynthesis and is required for conidiation and fruiting body formation in C. militaris.

Efficient CRISPR/Cas9 system based on autonomously replicating plasmid with an AMA1 sequence and precisely targeted gene deletion in the edible fungus, Cordyceps militaris

Cordyceps militaris is a popular edible fungus with important economic value worldwide. In this study, an efficient CRISPR/Cas9 genome-editing system based on an autonomously replicating plasmid with an AMA1 sequence was constructed. Further, a precisely targeted gene deletion via homology-directed repair was effectively introduced in C. militaris. Gene editing was successful, with efficiencies of 55.1% and 89% for Cmwc-1 and Cmvvd, respectively. Precisely targeted gene deletion was achieved at an efficiency of 73.9% by a single guide RNA supplementation with donor DNAs. Double genes, Cmwc-1 and Cmvvd, were edited simultaneously with an efficiency of 10%. Plasmid loss was observed under non-selective culture conditions, which could permit recycling of the selectable marker and avoid the adverse effects of the CRISPR/Cas9 system on the fungus, which is beneficial for the generation of new cultivars. RNA Pol III promoters, endogenous tRNAPro of C. militaris, and chimeric AfU6-tRNAGly can be used to improve the efficiency. Polyethylene glycol-mediated protoplast transformation was markedly more efficient than Agrobacterium tumefaciens-mediated transformation of C. militaris. To our knowledge, this is the first description of genome editing and precisely targeted gene deletion in mushrooms based on AMA1 plasmids. Our findings will enable the modification of multiple genes in both functional genomics research and strain breeding.

Dynamic Genome-Wide Transcription Profiling and Direct Target Genes of CmWC-1 Reveal Hierarchical Light Signal Transduction in Cordyceps militaris

Light is necessary for primordium differentiation and fruiting body development for most edible fungi; however, light perception and signal transduction have only been well studied in model fungi. In this study, a hierarchical network of transcriptional response to light in Cordyceps militaris, one of the edible fungi, has been described on a genome-wide scale using dynamic transcriptome analysis. It was shown that light regulated the transcript of 1722 genes, making up 18% of the whole genome of C. militaris. Analysis of light-responsive genes in C. militaris identified 4 categories: immediate-early, early, late, and continuous light-responsive genes, and the gene number increased distinctly with prolonged light exposure. Light-responsive genes with distinct functional categories showed specific time-dependent regulation. The target genes of CmWC-1, the most important photoreceptor, were revealed by ChIP-seq. A total of 270 significant peaks corresponding to 427 genes were identified to be directly regulated by CmWC-1, among which 143 genes respond to light. Based on 270 ChIP-seq peaks, the binding site for CmWC-1 was identified as AAATCAGACCAC/GTGGTCTGATTT, differing from the binding site by the homolog in Neurospora crassa. Elucidating the mechanisms of light perception and signal transduction will be helpful for further research on the fruiting body development in edible fungi.

Research Progress on Cordycepin Synthesis and Methods for Enhancement of Cordycepin Production in Cordyceps militaris

C. militaris is an insect-born fungus that belongs to Ascomycota and Cordyceps. It has a variety of biological activities that can be applied in medicine, health-care products, cosmeceuticals and other fields. Cordycepin (COR) is one of the major bioactive components identified from C. militaris. Thus, C. militaris and COR have attracted extensive attention. In this study, chemical synthetic methods and the biosynthesis pathway of COR were reviewed. As commercially COR was mainly isolated from C. militaris fermentation, the optimizations for liquid and solid fermentation and genetic modifications of C. militaris to increase COR content were also summarized. Moreover, the research progress of genetic modifications of C. militaris and methods for separation and purification COR were introduced. Finally, the existing problems and future research direction of C. militaris were discussed. This study provides a reference for the production of COR in the future.

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.

Transcriptome Analysis Reveals Candidate Genes Involved in Light-Induced Primordium Differentiation in Pleurotus eryngii

Pleurotus eryngii, a highly valued edible fungus, is one of the major commercially cultivated mushrooms in China. The development of P. eryngii, especially during the stage of primordium differentiation, is easily affected by light. However, the molecular mechanism underlying the response of primordium differentiation to light remains unknown. In the present study, primordium expression profiles under blue-light stimulation, red-light stimulation, and exposure to darkness were compared using high-throughput sequencing. A total of 16,321 differentially expressed genes (DEGs) were identified from three comparisons. GO enrichment analysis showed that a large number of DEGs were related to light stimulation and amino acid biosynthesis. KEGG analyses demonstrated that the MAPK signaling pathway, oxidative phosphorylation pathway, and RNA transport were most active during primordium differentiation. Furthermore, it was predicted that the blue-light photoreceptor WC-1 and Deoxyribodipyrimidine photolyase PHR play important roles in the primordium differentiation of P. eryngii. Taken together, the results of this study provide a speculative mechanism that light induces primordium differentiation and a foundation for further research on fruiting body development in P. eryngii.

Screening and Functional Verification of Selectable Marker Genes for Cordyceps militaris

The selectable marker genes are necessary resistance genes for gene knockout, gene complementation, and gene overexpression in filamentous fungi. Moreover, the more sensitive the filamentous fungi are to antibiotics, the more helpful it is to screen the target transformants. In order to obtain the antibiotic (or herbicide) which can effectively inhibit the growth of Cordyceps militaris and verify the function of the corresponding resistance gene in C. militaris, the sensitivity of C. militaris to hygromycin and glufosinate ammonium was compared to determine the resistance gene that was more suitable for the screening of C. militaris transformants. The binary vector of the selectable marker gene was constructed by combining the double-joint PCR (DJ-PCR) method and the homologous recombination method, and the function of the selectable marker gene in C. militaris was verified by the Agrobacterium tumefaciens-mediated transformation method. The results showed that C. militaris was more sensitive to glufosinate ammonium than hygromycin. The growth of C. militaris could be completely inhibited by 250 μg/mL glufosinate ammonium. The expression cassette of the glufosinate ammonium resistance gene (bar gene) was successfully constructed by DJ-PCR. The binary vector pCAMBIA0390-Bar was successfully constructed by homologous recombination. The bar gene of the vector pCAMBIA0390-Bar was successfully integrated into the C. militaris genome and could be highly expressed in the transformants of C. militaris. This study will promote the identification of C. militaris gene function and reveal the biosynthetic pathways of bioactive components in C. militaris.

Biomass and Cordycepin Production by the Medicinal Mushroom Cordyceps militaris-A Review of Various Aspects and Recent Trends towards the Exploitation of a Valuable Fungus

Cordyceps militaris is an entomopathogenic ascomycete with similar pharmacological importance to that of the wild caterpillar fungus Ophiocordyceps sinensis. C. militaris has attracted significant research and commercial interest due to its content in bioactive compounds beneficial to human health and the relative ease of cultivation under laboratory conditions. However, room for improvement exists in the commercial-scale cultivation of C. militaris and concerns issues principally related to appropriate strain selection, genetic degeneration of cultures, and substrate optimization. In particular, culture degeneration-usually expressed by abnormal fruit body formation and reduced sporulation-results in important economic losses and is holding back investors and potential growers (mainly in Western countries) from further developing this highly promising sector. In the present review, the main factors that influence the generation of biomass and metabolites (with emphasis on cordycepin biosynthesis) by C. militaris are presented and evaluated in conjunction with the use of a wide range of supplements or additives towards the enhancement of fungal productivity in large-scale cultivation processes. Moreover, physiological and genetic factors that increase or reduce the manifestation of strain degeneration in C. militaris are outlined. Finally, methodologies for developing protocols to be used in C. militaris functional biology studies are discussed.

Transcriptome sequencing and global analysis of blue light-responsive genes provide clues for high carotenoid yields in Blakeslea trispora

Blakeslea trispora has great potential uses in industrial production because of the excellent capability of producing a large quantity of carotenoids. However, the mechanisms of light-induced carotenoid biosynthesis even the structural and regulatory genes in pathways remain unclear. In this paper, we reported the first transcriptome study in B. trispora in which we have carried out global survey of expression changes of genes participated in blue light response. We verified that the yield of β-carotene increased 3-fold when transferred from darkness to blue light for 24 h and the enhancement of transcription levels of carRA and carB presented a positive correlation with the increase in carotenoid production. RNA-seq analysis revealed that 1124 genes were upregulated and 740 genes were downregulated respectively after blue light exposure. Annotation through GO, KEGG, Swissprot, and COG databases showed 11119 unigenes compared well with known gene sequences, 5514 unigenes were classified into Gene Ontology, and 4675 unigenes were involved in distinct pathways. Among the blue light-responsive genes, 4 genes (carG1, carG3, carRA and carB) identified to function in carotenoid metabolic pathways were dominantly upregulated. We also discovered that 142 TF genes belonging to 45 different superfamilies showed significant differential expression (p≤ 0.05), 62 of which were obviously repressed by blue light. The detailed profile of transcription data will not only allow us to conduct further functional genomics study in B. trispora, but also enhance our understanding of potential metabolic pathway and regulatory network involved in light-regulated carotenoid synthesis.

Hydrophobin CmHYD1 Is Involved in Conidiation, Infection and Primordium Formation, and Regulated by GATA Transcription Factor CmAreA in Edible Fungus, Cordyceps militaris

Hydrophobins are a family of small proteins exclusively secreted by fungi, and play a variety of roles in the life cycle. Cmhyd1, one of the hydrophobin class II members in Cordyceps militaris, has been shown to have a high transcript level during fruiting body development. Here, deletion of Cmhyd1 results in reduction in aerial mycelia, conidiation, hydrophobicity and infection ability, and complete inhibition of pigmentation and primordium differentiation. Cmhyd1 plays roles in conidiation and cuticle-bypassing infection by regulating the transcripts of frequency clock protein, Cmfrq, and velvet protein, Cmvosa, as well as primordium formation via the mitogen-activated protein kinase signaling pathway. Cmhyd1 also participates in stress response, including tolerance of mycelia to osmotic and oxidative stresses, and conidia to high or low temperatures. CmAreA, a transcription factor of nitrogen regulatory, is recruited to the promoter of Cmhyd1 and activates the transcription of Cmhyd1 with coactivator CmOTam using electrophoretic mobility shift assays and transient luciferase expression in tobacco. Furthermore, CmHYD1 is proved to regulate the transcription of Cmarea at different developmental stages via a positive feedback loop. These results reveal the diverse roles and regulation of Cmhyd1 in C. militaris, and provide insights into the developmental regulatory mechanism of mushrooms.

Structure prediction and function characterization of WC-2 proteins in Blakeslea trispora

Blakeslea trispora is known for its potential to produce an excess of carotenoids in mixed cultures of strains of opposite sex. The biosynthesis of β-carotene in B. trispora is activated not only by sex hormone trisporic acid but also by light, especially blue light. In fungi, the most intensively investigated blue-light reception proteins are WC-1 and WC-2, and the two proteins form a transcription factor complex which is called WCC by their PAS domains. Notably, multiple genes similar to wc-1 and wc-2 have been identified and characterized in Phycomyces, Mucor, and Rhizopus. Here we report that there are four members of wc-2-like gene family in B. trispora genome: Btwc-2a, Btwc-2b, Btwc-2c, and Btwc-2d. When the mycelia were exposed to blue light, their transcription levels are regulated differentially. Except for BtWC-2b, which only has a PAS domain, the other three proteins contain both a PAS domain and a ZnF domain. BtWC-2a interacts with either BtWC-1a or BtWC-1c to form different photoreceptor complexes in yeast two-hybrid assays, which is the unique situation not yet described in other fungi. In addition, the protein-protein docking analysis by the predicted 3D structures showed that the two complexes are structurally different. These results suggested that WC proteins of B. trispora are still involved in light regulation by forming WCC and the regulation mechanism of the photobiology appears to be more complex.

A light life together: photosensing in the plant microbiota

Bacteria and fungi of the plant microbiota can be phytopathogens, parasites or symbionts that establish mutually advantageous relationships with plants. They are often rich in photoreceptors for UVA-Visible light, and in many cases, they exhibit light regulation of growth patterns, infectivity or virulence, reproductive traits, and production of pigments and of metabolites. In addition to the light-driven effects, often demonstrated via the generation of photoreceptor gene knock-outs, microbial photoreceptors can exert effects also in the dark. Interestingly, some fungi switch their attitude towards plants in dependence of illumination or dark conditions in as much as they may be symbiotic or pathogenic. This review summarizes the current knowledge about the roles of light and photoreceptors in plant-associated bacteria and fungi aiming at the identification of common traits and general working ideas. Still, reports on light-driven infection of plants are often restricted to the description of macroscopically observable phenomena, whereas detailed information on the molecular level, e.g., protein-protein interaction during signal transduction or induction mechanisms of infectivity/virulence initiation remains sparse. As it becomes apparent from still only few molecular studies, photoreceptors, often from the red- and the blue light sensitive groups interact and mutually modulate their individual effects. The topic is of great relevance, even in economic terms, referring to plant-pathogen or plant-symbionts interactions, considering the increasing usage of artificial illumination in greenhouses, the possible light-regulation of the synthesis of plant-growth stimulating substances or herbicides by certain symbionts, and the biocontrol of pests by selected fungi and bacteria in a sustainable agriculture.

Cysteine-Rich Hydrophobin Gene Family: Genome Wide Analysis, Phylogeny and Transcript Profiling in Cordyceps militaris

Hydrophobins are a family of small secreted proteins found exclusively in fungi, and they play various roles in the life cycle. In the present study, genome wide analysis and transcript profiling of the hydrophobin family in Cordyceps militaris, a well-known edible and medicinal mushroom, were studied. The distribution of hydrophobins in ascomycetes with different lifestyles showed that pathogenic fungi had significantly more hydrophobins than saprotrophic fungi, and class II members accounted for the majority. Phylogenetic analysis of hydrophobin proteins from the species of Cordyceps s.l. indicated that there was more variability among the class II members than class I. Only a few hydrophobin-encoding genes evolved by duplication in Cordyceps s.l., which was inconsistent with the important role of gene duplication in basidiomycetes. Different transcript patterns of four hydrophobin-encoding genes during the life cycle indicated the possible different functions for each. The transcripts of Cmhyd2, 3 and 4 can respond to light and were related with the photoreceptors. CmQHYD, with four hydrophobin II domains, was first found in C. militaris, and multi-domain hydrophobins were only distributed in the species of Cordycipitaceae and Clavicipitaceae. These results could be helpful for further function research of hydrophobins and could provide valuable information for the evolution of hydrophobins.

Distribution, evolution and expression of GATA-TFs provide new insights into their functions in light response and fruiting body development of Tolypocladium guangdongense

Background

Fungal GATA-type transcription factors (GATA-TFs) are a class of transcriptional regulators involved in various biological processes. However, their functions are rarely analyzed systematically, especially in edible or medicinal fungi, such as Tolypocladium guangdongense, which has various medicinal and food safety properties with a broad range of potential applications in healthcare products and the pharmaceutical industry.

Methods

GATA-TFs in T. guangdongense (TgGATAs) were identified using InterProScan. The type, distribution, and gene structure of TgGATAs were analyzed by genome-wide analyses. A phylogenetic tree was constructed to analyze their evolutionary relationships using the neighbor-joining (NJ) method. To explore the functions of GATA-TFs, conserved domains were analyzed using MEME, and cis-elements were predicted using the PlantCARE database. In addition, the expression patterns of TgGATAs under different light conditions and developmental stages were studied using qPCR.

Results

Seven TgGATAs were identified. They were randomly distributed on four chromosomes and contained one to four exons. Phylogenetic analysis indicated that GATA-TFs in each subgroup are highly conserved, especially for GATA1 to GATA5. Intron distribution analyses suggested that GATA1 and GATA3 possessed the most conserved gene structures. Light treatments induced the expression levels of TgGATA1 and TgGATA5-7, but the expression levels varied depending on the duration of illumination. The predicted protein structures indicate that TgGATA1 and TgGATA2 possess typical light-responsive domains and may function as photoreceptors to regulate downstream biological processes. TgGATA3 and TgGATA5 may be involved in nitrogen metabolism and siderophore biosynthesis, respectively. TgGATA6 and TgGATA7 possess unique Zn finger loop sequences, suggesting that they may have special functions. Furthermore, gene expression analysis indicated that TgGATA1 (WC1) was notably involved in mycelial color transformation, while other genes were involved in fruiting body development to some extent. These results provide valuable information to further explore the mechanisms through which TgGATAs are regulated during fruiting body development.

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.

Analysis of reference genes stability and histidine kinase expression under cold stress in Cordyceps militaris

The development of fungal fruiting bodies from a hyphal thallus is inducible under low temperature (cold stress). The molecular mechanism has been subject to surprisingly few studies. Analysis of gene expression level has become an important means to study gene function and its regulation mechanism. But identification of reference genes (RGs) stability under cold stress have not been reported in famous medicinal mushroom-forming fungi Cordyceps militaris. Herein, 12 candidate RGs had been systematically validated under cold stress in C. militaris. Three different algorithms, geNorm, NormFinder and BestKeeper were applied to evaluate the expression stability of the RGs. Our results showed that UBC and UBQ were the most stable RGs for cold treatments in short and long periods, respectively. 2 RGs (UBC and PP2A) and 3 RGs (UBQ, TUB and CYP) were the suitable RGs for cold treatments in short and long periods, respectively. Moreover, target genes, two-component-system histidine kinase genes, were selected to validate the most and least stable RGs under cold treatment, which indicated that use of unstable expressed genes as RGs leads to biased results. Our results provide a good starting point for accurate reverse transcriptase quantitative polymerase chain reaction normalization by using UBC and UBQ in C. militaris under cold stress and better support for understanding the mechanism of response to cold stress and fruiting body formation in C. militaris and other mushroom-forming fungi in future research.

MAT1-1-3, a Mating Type Gene in the Villosiclava virens, Is Required for Fruiting Bodies and Sclerotia Formation, Asexual Development and Pathogenicity

Villosiclava virens is the prevalent causative pathogen of rice false smut, a destructive rice disease. Mating-type genes play a vital role in the evolution of mating systems in fungi. Some fungi have lost MAT1-1-3, one of the mating-type genes, during evolution, whereas others still retain MAT1-1-3. However, how MAT1-1-3 regulates the sexual development of heterothallic V. virens remains unknown. Here, we generated the MAT1-1-3 mutants, which exhibited defects in vegetative growth, stress response, pathogenicity, sclerotia formation and fruiting body maturation. An artificial outcrossing inoculation assay showed that the Δmat1-1-3 mutant was unable to produce sclerotia. Unexpectedly, the Δmat1-1-3 mutant could form immature fruiting bodies without mating on potato sucrose agar medium (PSA) compared with the wild-type strain, most likely by activating the truncated MAT1-2-1 transcription to regulate the sexual development. Moreover, RNA-seq data showed that knockout of MAT1-1-3 results in misregulation of a subset of genes involved in sexual development, MAPK signaling, cell wall integrity, autophagy, epigenetic modification, and transcriptional regulation. Collectively, this study reveals that MAT1-1-3 is required for asexual and sexual development, and pathogenicity of V. virens, thereby provides new insights into the function of mating-type genes in the fungi life cycle and infection process.

Genome-wide characterization of the Zn(II)2Cys6 zinc cluster-encoding gene family in Pleurotus ostreatus and expression analyses of this family during developmental stages and under heat stress

Pleurotus ostreatus is one of the most widely cultivated mushrooms in China. The regulatory mechanisms of fruiting body formation and the response to heat stress in P. ostreatus are main research focuses. The Zn(II)₂Cys₆ family is one of the largest families of transcriptional factors and plays important roles in multiple biological processes in fungi. In this study, we identified 66 zinc cluster proteins in P. ostreatus (PoZCPs) through a genome-wide search. The PoZCPs were classified into 15 types according to their zinc cluster domain. Physical and chemical property analyses showed a huge diversity among the PoZCPs. Phylogenetic analysis of PoZCPs classified these proteins into six groups and conserved motif combinations and similar gene structures were observed in each group. The expression profiles of these PoZCP genes during different developmental stages and under heat stress were further investigated by RNA-sequencing (RNA-seq), revealing diverse expression patterns. A total of 13 PoZCPs that may participate in development or the heat stress response were selected for validation of their expression levels through real-time quantitative PCR (RT-qPCR) analysis, and some developmental stage-specific and heat stress-responsive candidates were identified. The findings contribute to our understanding of the roles and regulatory mechanisms of ZCPs in P. ostreatus.

A WD40 Protein Encoding Gene Fvcpc2 Positively Regulates Mushroom Development and Yield in Flammulina velutipes

Ascomycota and Basidiomycota are two closely related phyla and fungi in two phyla share some common morphological developmental process during fruiting body formation. In Neurospora crassa, the Gβ-like protein CPC-2 with a seven-WD40 repeat domain was previously reported. By transforming CPC-2 ortholog encoding genes, from 7 different fungal species across Ascomycota and Basidiomycota, into the cpc-2 deletion mutant of N. crassa, we demonstrate that all tested CPC-2 ortholog genes were able to complement the defects of the cpc-2 deletion mutant in sexual development, indicating that CPC-2 proteins from Ascomycota and Basidiomycota have the similar cellular function. Using Flammulina velutipes as a model system for mushroom species, the CPC-2 ortholog FvCPC2 was characterized. Fvcpc2 increased transcription during fruiting body development. Knockdown of Fvcpc2 by RNAi completely impaired fruiting body formation. In three Fvcpc2 knockdown mutants, transcriptional levels of genes encoding adenylate cyclase and protein kinase A catalytic subunit were significantly lower and colony growth became slower than wild type. The addition of cAMP or the PKA-activator 8-Bromo-cAMP into the medium restored the Fvcpc2 knockdown mutants to the wild-type colony growth phenotype, suggesting that the involvement of cAMP production in the regulatory mechanisms of FvCPC2. Knockdown of Fvcpc2 also weakened transcriptional responses to sexual development induction by some genes related to fruiting body development, including 4 jacalin-related lectin encoding genes, 4 hydrophobin encoding genes, and 3 functionally-unknown genes, suggesting the participation of these genes in the mechanisms by which FvCPC2 regulates fruiting body development. All three Fvcpc2 overexpression strains displayed increased mushroom yield and shortened cultivation time compared to wild type, suggesting that Fvcpc2 can be a promising reference gene for Winter Mushroom breeding. Since the orthologs of FvCPC2 were highly conserved and specifically expressed during fruiting body development in different edible mushrooms, genes encoding FvCPC2 orthologs in other mushroom species also have potential application in breeding.

Cmfhp Gene Mediates Fruiting Body Development and Carotenoid Production in Cordyceps militaris

Cordyceps militaris fruiting bodies contain a variety of bioactive components that are beneficial to the human body. However, the low yield of fruiting bodies and the low carotenoid content in C. militaris have seriously hindered the development of the C. militaris industry. To elucidate the developmental mechanism of the fruiting bodies of C. militaris and the biosynthesis mechanism of carotenoids, the function of the flavohemoprotein-like Cmfhp gene of C. militaris was identified for the first time. The Cmfhp gene was knocked out by the split-marker method, and the targeted gene deletion mutant ΔCmfhp was obtained. An increased nitric oxide (NO) content, no fruiting body production, decreased carotenoid content, and reduced conidial production were found in the mutant ΔCmfhp. These characteristics were restored when the Cmfhp gene expression cassette was complemented into the ΔCmfhp strain by the Agrobacterium tumefaciens-mediated transformation method. Nonetheless, the Cmfhp gene had no significant effect on the mycelial growth rate of C. militaris. These results indicated that the Cmfhp gene regulated the biosynthesis of NO and carotenoids, the development of fruiting bodies, and the formation of conidia. These findings potentially pave the way to reveal the developmental mechanism of fruiting bodies and the biosynthesis mechanism of carotenoids in C. militaris.

Integration of ATAC-Seq and RNA-Seq Identifies Key Genes in Light-Induced Primordia Formation of Sparassis latifolia

Light is an essential environmental factor for Sparassis latifolia primordia formation, but the molecular mechanism is still unclear. In this study, differential expression profiling of light-induced primordia formation (LIPF) was established by integrating the assay for transposase accessible chromatin by sequencing (ATAC-seq) and RNA-seq technology. The integrated results from the ATAC-seq and RNA-seq showed 13 down-regulated genes and 17 up-regulated genes in both the L vs. D and P vs. D groups, for both methods. According to the gene ontology (GO) annotation of these differentially expressed genes (DEGs), the top three biological process categories were cysteine biosynthetic process via cystathionine, vitamin B6 catabolic, and glycine metabolic; the top three molecular function categories were 5-methyltetrahydropteroyltriglutamate-homocysteine S-methyltransferase activity, glycine binding, and pyridoxal phosphate binding; cellular component categories were significantly enriched in the glycine cleavage complex. The KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis revealed that these genes were associated with vitamin B6 metabolism; selenocompound metabolism; cysteine and methionine metabolism; glycine, serine, and threonine metabolism; and glyoxylate and dicarboxylate metabolism pathways. The expression of most of the DEGs was validated by qRT-PCR. To the best of our knowledge, this study is the first integrative analysis of ATAC-seq and RNA-seq for macro-fungi. These results provided a new perspective on the understanding of key pathways and hub genes in LIPF in S. latifolia. It will be helpful in understanding the primary environmental response, and provides new information to the existing models of primordia formation in edible and medicinal fungi.

CmVVD is involved in fruiting body development and carotenoid production and the transcriptional linkage among three blue-light receptors in edible fungus Cordyceps militaris

Fruiting body development and carotenoid production are light-induced in Cordyceps militaris. Our previous studies have shown that two blue-light receptors, CmWC-1 and CmCRY-DASH, regulate fruiting body development and secondary metabolism. However, the photosensory system of C. militaris remains unclear. Here, gene deletion of Cmvvd, coding for another blue-light receptor, resulted in reduced conidiation level and significant promotion of carotenoid content. Cmvvd transcription levels at fruiting body stages were higher than at other stages, and fruiting bodies could not develop normally in ΔCmvvd strains, indicating that Cmvvd might play an important role in fruiting body development. Rhythm loops were not affected in ΔCmvvd strains but were regulated by Cmwc-1, and the expression of the rhythm regulator gene Cmfrq was dependent on CmWC-1. Chromatin immunoprecipitation assay confirmed that Cmvvd is the direct target of CmWC-1 in this fungus. Our results also revealed interdependent transcriptional relationships between Cmwc-1 and Cmvvd, and between Cmwc-1 and Cmcry-DASH. Cmcry-DASH expression was affected by Cmvvd, and the function-loss of Cmcry-DASH might be compensated by the high transcription of Cmvvd. This is the first report of the transcriptional linkage among the three blue-light receptors in edible fungi and will be helpful for studies of multicellular development in this fungus.

Identification of two Pleurotus ostreatus blue light receptor genes (PoWC-1 and PoWC-2) and in vivo confirmation of complex PoWC-12 formation through yeast two hybrid system

Blue light is necessary for initiation of mushroom formation of Pleurotus ostreatus. In this study, we isolated homologues of the blue light receptor genes wc-1 and wc-2 from P. ostreatus, PoWC-1 and PoWC-2. The PoWC-1 contained three typical PAS domains and one PAS domain exhibited significant similarity to the LOV domain of known blue light receptors. The PoWC-2 had one typical PAS domain and one ZnF domain. The qRT-PCR analysis showed that PoWC-1 and PoWC-2 expression increased in a short time, and the final level tended to be stable along with the light illumination. The PoWC-1 and PoWC-2 expression levels of the primordium period was higher than that of mature fruiting-body period; and in the pileus were the highest, followed by the stipe and the gills. The expression of PoWC-1 and PoWC-2 in pre-primordial mycelia was induced by light exposure. In vivo analysis through yeast two-hybrid experiment disclosed that PoWC-1 and PoWC-2 could form heterologous complex to activate the reporter genes and the complex perform the transcription factor function requiring the addition of FAD.

Transcriptome Analysis of Cordyceps militaris Reveals Genes Associated With Carotenoid Synthesis and Identification of the Function of the Cmtns Gene

Cordyceps militaris, a valuable edible and medicinal fungus, has attracted increasing attention because of its various bioactive ingredients. However, the biosynthetic pathway of C. militaris carotenoids is still unknown due to lack of transcriptome information. To uncover genes related to the biosynthesis of C. militaris carotenoids, the transcriptomes of mycelia CM10_D cultured under dark conditions and mycelia CM10_L cultured under light exposure conditions were sequenced. Compared with mycelia CM10_D, 866 up-regulated genes and 856 down-regulated genes were found in mycelia CM10_L. Gene ontology (GO) analysis of differentially expressed genes (DEGs) indicated that DEGs were mainly classified into the “metabolic process,” “membrane,” and “catalytic activity” terms. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of DEGs suggested that DEGs were mainly enriched in “metabolic pathways,” “MAPK signaling pathway-yeast,” and “biosynthesis of secondary metabolites.” In addition, the carotenoid content of the Cmtns gene deletion mutant (ΔCmtns) was significantly lower than that of the wild-type C. militaris CM10, while the carotenoid content of the complementary strain (ΔCmtns-c) of the Cmtns gene was not significantly different from that of C. militaris CM10, suggesting that the Cmtns gene significantly affected the biosynthesis of carotenoids in C. militaris. These results potentially pave the way for revealing the biosynthetic pathway of carotenoids and improving carotenoids production in C. militaris.

The Blue-Light Photoreceptor Sfwc-1 Gene Regulates the Phototropic Response and Fruiting-Body Development in the Homothallic Ascomycete Sordaria fimicola

Sordaria fimicola, a coprophilous ascomycete, is a homothallic fungus that can undergo sexual differentiation with cellular and morphological changes followed by multicellular tissue development to complete its sexual cycle. In this study, we identified and characterized the blue-light photoreceptor gene in S. fimicola The S. fimicola white collar-1 photoreceptor (SfWC-1) contains light-oxygen-voltage-sensing (LOV), Per-Arnt-Sim (PAS), and other conserved domains and is homologous to the WC-1 blue-light photoreceptor of Neurospora crassa The LOV domain of Sfwc-1 was deleted by homologous recombination using Agrobacterium-mediated protoplast transformation. The Sfwc-1(Δlov) mutant showed normal vegetative growth but produced less carotenoid pigment under illumination. The mutant showed delayed and less-pronounced fruiting-body formation, was defective in phototropism of the perithecial beaks, and lacked the fruiting-body zonation pattern compared with the wild type under the illumination condition. Gene expression analyses supported the light-induced functions of the Sfwc-1 gene in the physiology and developmental process of perithecial formation in S. fimicola Moreover, green fluorescent protein (GFP)-tagged SfWC-1 fluorescence signals were transiently strong upon light induction and prominently located inside the nuclei of living hyphae. Our studies focused on the putative blue-light photoreceptor in a model ascomycete and contribute to a better understanding of the photoregulatory functions and networks mediated by the evolutionarily conserved blue-light photoreceptors across diverse fungal phyla.IMPORTANCESordaria sp. has been a model for study of fruiting-body differentiation in fungi. Several environmental factors, including light, affect cellular and morphological changes during multicellular tissue development. Here, we created a light-oxygen-voltage-sensing (LOV) domain-deleted Sfwc-1 mutant to study blue-light photoresponses in Sordaria fimicola Phototropism and rhythmic zonation of perithecia were defective in the Sfwc-1(Δlov) mutant. Moreover, fruiting-body development in the mutant was reduced and also significantly delayed. Gene expression analysis and subcellular localization study further revealed the light-induced differential gene expression and cellular responses upon light stimulation in S. fimicola.

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.

Developmental transcriptomics of Chinese cordyceps reveals gene regulatory network and expression profiles of sexual development-related genes

BACKGROUND

Chinese cordyceps, also known as Chinese caterpillar fungus (Ophiocordyceps sinensis, syn. Cordyceps sinensis), is of particular interest for its cryptic life cycle and economic and ecological importance. The large-scale artificial cultivation was succeeded recently after several decades of efforts and attempts. However, the induction of primordium, sexual development of O. sinensis and the molecular basis of its lifestyle still remain cryptic.

RESULTS

The developmental transcriptomes were analyzed for six stages covering the whole developmental process, including hyphae (HY), sclerotium (ST), primordium (PR), young fruiting body (YF), developed fruiting body (DF) and mature fruiting body (MF), with a focus on the expression of sexual development-related genes. Principal component analysis revealed that the gene expression profiles at the stages of primordium formation and fruiting body development are more similar than those of the undifferentiated HY stage. The PR and MF stages grouped together, suggesting that primordium differentiation and sexual maturation have similar expression patterns. Many more DEGs were identified between the ST and HY stages, covering 47.5% of the O. sinensis genome, followed by the comparisons between the ST and PR stages. Using pairwise comparisons and weighted gene coexpression network analysis, modules of coexpressed genes and candidate hub genes for each developmental stage were identified. The four mating type loci genes expressed during primordium differentiation and sexual maturation; however, spatiotemporal specificity of gene expression indicated that they also expressed during the anamorphic HY stage. The four mating type genes were not coordinately expressed, suggesting they may have divergent roles. The expression of the four mating type genes was highest in the fertile part and lowest in the sclerotium of the MF stage, indicating that there is tissue specificity. Half of genes related to mating signaling showed as the highest expression in the ST stage, indicating fruiting was initiated in the ST stage.

CONCLUSIONS

These results provide a new perspective to understanding of the key pathways and hub genes, and sexual development-related gene profile in the development of Chinese cordyceps. It will be helpful for underlying sexual reproduction, and add new information to existing models of fruiting body development in edible fungi.

Genome-Wide Analysis of the Zn(II)₂Cys₆ Zinc Cluster-Encoding Gene Family in Tolypocladium guangdongense and Its Light-Induced Expression

The Zn(II)₂Cys₆ zinc cluster gene family is a subclass of zinc-finger proteins, which are transcriptional regulators involved in a wide variety of biological processes in fungi. We performed genome-wide identification and characterization of Zn(II)₂Cys₆ zinc-cluster gene (C6 zinc gene) family in Tolypocladiumguangdongense, Cordycepsmilitaris and Ophiocordycepssinensis. Based on the structures of the C6 zinc domains, these proteins were observed to be evolutionarily conserved in ascomycete fungi. We focused on T.guangdongense, a medicinal fungus, and identified 139 C6 zinc genes which could be divided into three groups. Among them, 49.6% belonged to the fungal specific transcriptional factors, and 16% had a DUF3468 domain. Homologous and phylogenetic analysis indicated that 29 C6 zinc genes were possibly involved in the metabolic process, while five C6 zinc genes were supposed to be involved in asexual or sexual development. Gene expression analysis revealed that 54 C6 zinc genes were differentially expressed under light, including two genes that possibly influenced the development, and seven genes that possibly influenced the metabolic processes. This indicated that light may affect the development and metabolic processes, at least partially, through the regulation of C6 zinc genes in T.guangdongense. Our results provide comprehensive data for further analyzing the functions of the C6 zinc genes.

milR4 and milR16 Mediated Fruiting Body Development in the Medicinal Fungus Cordyceps militaris

Cordyceps militaris readily performs sexual reproduction, thus providing a remarkably rich model for understanding the processes involved in sexual development. It could regulate expression of human genes by diet-derived miRNA-like RNAs (milRNAs). However, the study of miRNAs in C. militaris has been limited. In the present study, genes encoding Dicers, Argonautes, and RNA-dependent RNA polymerases were identified. Illumina deep sequencing was performed to characterize the milRNAs in C. militaris at asexual and sexual development stages. Total 38 milRNAs were identified and five milRNAs were validated by northern blot and qRT-PCR, out of which, 19 were specific for sexual development. Importantly, the fungi could not form fruiting bodies after disruption of milR4, while the perithecium was formed in advance after over-expression of milR4. Abnormal pale yellow fruiting body primordium, covered with abnormal primordium, was formed in the strain with miR16 disruption. Although no milR4 or milR16 target genes were identified, differential expression of many different genes involved in mycelium growth and sexual development (mating process, mating signaling, and fruiting body development) among these mutants were found. Overall, milRNAs play vital roles in sexual development in C. militaris.