Genome sequence of the insect pathogenic fungus Cordyceps militaris, a valued traditional Chinese medicine

Transcriptional variation and RNA polymorphism among different Lentinula edodes (Berk.) Pegler strains

BACKGROUND

Lentinula edodes is a commercially important mushroom known for its nutritional and therapeutic values. However, the molecular mechanisms underlying the distinct nutritional and physiological attributes of various L. edodes strains are not well understood. This study focused on three Lentinula strains (DMRO-356, DMRO-623, and DMRO-388s) with different nutritional and productivity profiles. Illumina sequencing was used to perform a whole-transcriptome analysis, conducting 100-base pair paired-end sequencing of total messenger RNA (mRNA) in duplicate, resulting in 28-48 million sequencing reads per strain. After rigorous data filtering, over 99% of high-quality reads were retained, and more than 95% were aligned to the Lentinula genome.

RESULTS

Differential gene expression analyses identified 2210 differentially expressed genes between DMRO-356 and DMRO-623, 862 between DMRO-356 and DMRO-388s, and 2212 between DMRO-623 and DMRO-388s. Significant genetic variations were found among the strains, including 7753 single nucleotide polymorphisms (SNPs) in DMRO-356 versus DMRO-623 and 4080 SNPs in DMRO-356 versus DMRO-388s. Additionally, 349 insertions/deletions (InDels) were found in DMRO-356/DMRO-623 and 218 in DMRO-356/DMRO-388 s. Non-synonymous SNPs, which alter amino acid compositions, were analyzed, showing a preference for polar over charged amino acids.

CONCLUSION

These differentially expressed genes were associated with various nutritional and developmental processes, highlighting the importance of genetic variations in shaping amino acid composition and potentially affecting protein function. This study is the first comprehensive exploration of transcriptional differences among Lentinula strains available for its cultivation, providing valuable insights to enhance mushroom quality and productivity. © 2024 Society of Chemical Industry.

Alanine supplementation enhancing cordycepin production in Cordyceps militaris via upregulation of Cns2 and Cns3 genes expression levels

Cordycepin, a key bioactive compound produced by Cordyceps militaris, faces the challenge of low productivity for commercial use. In this study, alanine supplementation in Cordyceps militaris boosted cordycepin production, peaking at 3 mg/g with 12 g/L concentration. Transcriptome analysis revealed 1711 differentially expressed genes, Pathway analysis indicates that protein processing in the endoplasmic reticulum was the most affected pathway. In addition, the transcriptome showed that adenylosuccinate lyase is essential for the synthesis of cordycepin. The modulation of four genes (Cns1-4) points to a regulatory mechanism that could increase cordycepin biosynthesis, offering a strategy to overcome low productivity for commercial applications.

Topical Infection of Cordyceps militaris in Silkworm Larvae Through the Cuticle has Lower Infectivity Compared to Beauveria bassiana and Metarhizium anisopliae

Topical infection of entomopathogenic fungi in insects occurs when the fungal conidia attach to the insect's surface (cuticle), germinate, and then form appressoria that penetrate the cuticle and enter their bodies. In this study, we inoculated silkworm larvae with three entomopathogenic fungi, Cordyceps militaris, Beauveria bassiana, and Metarhizium anisopliae, and investigated their mechanisms of infection. Attachment of the conidia of the three entomopathogenic fungi to the surface of silkworm larvae was observed under a microscope. We counted the number of conidia attached to the surface of the silkworm larvae and the number of conidia detached from the surface was counted. The number of C. militaris conidia that attached to the surface was less than that of B. bassiana and M. anisopliae; however, it germinated and formed appressoria on hydrophobic surfaces, similar to the other two strains. Mycelial growth of C. militaris was inhibited compared to that of B. bassiana in PDA medium containing 0.1% linoleic and linolenic acids. The germination of C. militaris conidia was also inhibited in PD medium containing 0.1% linoleic or linolenic acids. These results suggest that the attachment of low numbers of C. militaris conidia on the surface of silkworm larvae and presence of inhibitory linoleic or linolenic acids in the silkworm cuticles may cause low topical infectivity by C. militaris. This study improves the efficacy of topically infecting silkworms with C. militaris to produce fungal fruiting bodies for use in traditional Chinese medicine and dietary supplement production.

Discovering a novel glycosyltransferase gene CmUGT1 enhances main metabolites production of Cordyceps militaris

Cordyceps militaris is a filamentous fungus used for both medicinal and culinary purposes. It exhibits a wide range of pharmacological activities due to its valuable contents of cordycepin, polysaccharides, carotenoids, terpenoids and other metabolites. However, C. militaris strains are highly susceptible to irreversible degradation in agricultural production, which is often manifested as a prolonged color change period and a significant decrease in the production of secondary metabolites. UDP-glycosyltransferases are an important enzyme family that participates in the synthesis of terpenoids by performing the glycosylation of key residues of enzymes or molecules. However, few studies have focused on its effect on the regulation of metabolite production in C. militaris. Therefore, in this study, we performed transcriptome analysis across four different developmental stages of C. militaris to target the putative glycosyltransferase gene CmUGT1, which plays important roles in metabolite production. We further constructed and screened a CmUGT1-overexpressing strain by Agrobacterium tumefaciens-mediated infestation of C. militaris spores. The major metabolite production of the wild-type and CmUGT1-overexpressing C. militaris strains was determined after short-term shake-flask cultivation of mycelia. The results showed that the yields of carotenoids and polysaccharides in the mycelia of the CmUGT1-overexpressing strains were 3.8 and 3.4 times greater than those in the mycelia of the wild type, respectively (p < 0.01). The levels of intracellular and extracellular cordycepin produced by the overexpression strain were 4.4 and 8.0 times greater than those produced by the wild-type strain (p < 0.01). This suggests that the overexpression of CmUGT1 in C. militaris enhances the synthesis activities of the main enzymes related to metabolite production, which provides a guide for obtaining excellent recombinant strains of C. militaris.

Core microbes in Cordyceps militaris sclerotia and their nitrogen metabolism-related ecological functions

Cordyceps militaris infects insects and forms sclerotia within the insect remains, establishing insect-microbe complexes. Here, C. militaris sclerotia samples from a single location in China over a 5-year period were subjected to high-throughput DNA sequencing, and the core microbes (which were stably enriched in the sclerotia over the 5 years) were identified. Next, seven bacterial strains were isolated from the C. militaris sclerotia, their biochemical characteristics were assessed, and they were co-cultured with C. militaris to study their effects on C. militaris metabolite production and biomass. Furthermore, the effects of NH4, NO3, and peptone media on C. militaris were compared. The results showed that Rhodococcus, Phyllobacterium, Pseudomonas, Achromobacter, Ensifer, Stenotrophomonas, Sphingobacterium, Variovorax, and Acinetobacter were the core microbes. Although co-culture of C. militaris with the seven bacterial strains isolated from the sclerotia did not directly increase the cordycepin level, they all had NO3 reduction ability, and four had urea decomposition ability. Meanwhile, C. militaris in NH4 medium had an increased cordycepin level compared to C. militaris in the other two media. From this, we inferred that bacteria in the sclerotia can convert NO3 to NH4, and then cordycepin is produced using NH4, which was confirmed by RNA-seq and real-time fluorescence quantitative PCR. Thus, bacteria in the sclerotia may indirectly affect the C. militaris metabolite production by regulating nitrogen metabolism. In summary, there are stable core microbes in the C. militaris sclerotia, and they may directly and indirectly affect the growth and metabolite production of C. militaris.

IMPORTANCE

The model Cordyceps species Cordyceps militaris is rich in therapeutic compounds. It has recently been demonstrated that symbiotic microbes in sclerotia affect Cordyceps' growth, development, and secondary metabolite production. In this study, core microbes were identified based on C. militaris sclerotia samples obtained from the same site over 5 years. Additionally, bacterial strains isolated from C. militaris sclerotia were found to affect metabolite production and nitrogen utilization, based on functional tests. Moreover, based on the bacterial nitrogen metabolism capacity in the sclerotia and its influence on C. militaris metabolite production, we deduced that bacteria in the sclerotia can indirectly affect C. militaris metabolite production by regulating nitrogen metabolism. This is the first report on how bacteria in the sclerotia affect C. militaris metabolite production from the perspective of the nitrogen cycle. The results increase our understanding of microbial functions in C. militaris sclerotia.

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

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

Holistic transcriptional responses of Cordyceps militaris to different culture temperatures

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

Application and research progress of cordycepin in the treatment of tumours (Review)

Cordycepin is a nucleoside molecule found in Cordyceps sinensis and can be obtained through chemical synthesis and biotransformation. Cordycepin has been extensively studied and has been shown to have antitumour activity. This activity includes effects on the autophagy process and inhibition of the MAPK/ERK and Hedgehog pathways. Ultimately, the inhibitory effect of cordycepin on tumour cells is due to the interplay of these effects. Cordycepin was shown to enhance the therapeutic effects of radiotherapy. There is increasing evidence indicating that cordycepin plays an anticancer role in the treatment of various cancers. The present review aims to provide a clear understanding of the antitumour mechanisms of cordycepin and discuss its present application in the treatment of tumours. This information can be an important theoretical basis and provide clinical guidance for the further development of cordycepin as an antitumour drug.

An sn-2 regioselective lipase with cis-fatty acid preference from Cordyceps militaris: Biochemical characterization and insights into its regioselective mechanism

Lipase with unique regioselectivity is an attractive biocatalyst for elaborate lipid modification. However, the excavation of novel sn-2 regioselective lipases is difficult due to their scarcity in nature, with Candida antarctica lipase A (CALA) being the pronouncedly reported one. Here, we identified a novel CALA-like lipase from Cordyceps militaris (CACML7) via in silico mining. Through chiral-phase high-performance liquid chromatography, we determined that CACML7 displays sn-2 regioselectivity (>68 %) as does CALA, but exhibits distinctive chain length selectivity and bias against unsaturated fats. Notably, the curvature of the acyl-binding tunnel was expected to contribute to the 2.2-fold higher preference for cis-fatty acid (C18:1, cis-Δ9) over trans-fatty acid (C18:1, trans-Δ9) unlike trans-active CALA. Random pose docking of trioleoylglycerol (TOG) into the active site of a lid-truncated mutant of CACML7 revealed that TOG accepts a tuning fork conformation, of which the precise positioning of the reactive ester group towards the catalytic center was only favorable via sn-2 binding mode. The unique active site morphology, which we refer to as an "acyl-binding tunnel with a narrow entrance," may contribute to the sn-2 regioselectivity of CACML7. Our data provide an attractive model to better understand the mechanism underlying sn-2 regioselectivity.

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.

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.

Innovative application of CRISPR for eliminating Ustiloxin in Cordyceps militaris: Enhancing food safety and quality

Cordyceps militaris (L.) Fr. Has long been recognized as a valuable functional food consumed in numerous countries. However, biosynthetic gene clusters of this species and safety regarding mycotoxin production remain largely unexplored. In this study, a ribosomally synthesized and post-translationally modified peptide (RiPP) cluster responsible for the production of cyclopeptide mycotoxins in Cordyceps was unveiled via genome mining. Ustiloxin B and a novel, predominant and Cordyceps specific ustiloxin I were confirmed by extraction and structural analysis. The difference between Ustiloxins I and B lied in the side chain at C19, where an additional methyl substituent in Ustiloxin I resulted in an alanine moiety substitution for glycine of Ustiloxin B. The simultaneous deletion of the two adjacent core genes, CmustYb and CmustYa, using a single guide RNA designed in the intergenic region, and subsequent in-situ complementation via AMA-mediated CRISPR/Cas9 system confirmed the RiPP cluster's responsibility for ustiloxin production. The cultivation of the edited strain yielded ustiloxin-free fruiting bodies without affecting agronomic characters. PCR and genome resequencing confirmed the absence of any off-target events or foreign sequence remnants. This study marks a significant advancement in utilizing CRISPR technology to control ustiloxins in food, underscoring its broader implications for food safety and quality improvement.

Effects of mating-type ratio imbalance on the degeneration of Cordyceps militaris subculture and preventative measures

The rapid degeneration of Cordyceps militaris strains during subculture represents a bottleneck problem that affects production stability. This study explored the mechanism underlying this degeneration in three production and three wild-type strains of Cordyceps militaris, isolating single-conidium strains from each. The effects of subculturing on fructification in both original and single mating-type strains were compared. Changes in the ratio of the two mating types were analyzed in both original and degenerated strains. Based on these findings, the two mating strains were paired in different ratios to determine their effects on fruiting. The resulting five strains were heterokaryotic strains with both MAT1-1 and MAT1-2 mating-type genes. Strain jb-2 was a single mating type (MAT1-1) mutant strain that produced stable fruiting bodies but failed to produce ascospores. It was found that the loss of or imbalance in mating types was the main reason for the rapid degeneration of fruiting traits during subculture and that this occurred randomly in the MAT1-1 and MAT1-2 types. The strains differed significantly in their stability during subculture. Fruiting was stable in the single mating-type Jb-2 strain, and the eleventh-generation fruited normally. There were differences in yield between the production and wild strains after inoculation with spawn containing different proportions of mating types. The production strain was more stable when inoculated with strains with mating-type ratios of 1:9 to 9:1 without affecting the yield. However, the yield of the wild-type strain xf-1 was positively correlated with the proportion of the MAT1-2 type, while the other two strains showed no correlations. Subculturing single mating-type mycelia separately and mixing them before production effectively mitigated degeneration during subculture. For Cordyceps militaris breeding, selecting strains containing both mating types, which are insensitive to the proportion of mating-type genes, enhanced stability in subculture and reduced the risk of mating-type loss. Direct breeding of specific single-mating type strains to induce fruiting is thus an effective breeding strategy.

What are the 100 most cited fungal genera?

The global diversity of fungi has been estimated between 2 to 11 million species, of which only about 155 000 have been named. Most fungi are invisible to the unaided eye, but they represent a major component of biodiversity on our planet, and play essential ecological roles, supporting life as we know it. Although approximately 20 000 fungal genera are presently recognised, the ecology of most remains undetermined. Despite all this diversity, the mycological community actively researches some fungal genera more commonly than others. This poses an interesting question: why have some fungal genera impacted mycology and related fields more than others? To address this issue, we conducted a bibliometric analysis to identify the top 100 most cited fungal genera. A thorough database search of the Web of Science, Google Scholar, and PubMed was performed to establish which genera are most cited. The most cited 10 genera are Saccharomyces, Candida, Aspergillus, Fusarium, Penicillium, Trichoderma, Botrytis, Pichia, Cryptococcus and Alternaria. Case studies are presented for the 100 most cited genera with general background, notes on their ecology and economic significance and important research advances. This paper provides a historic overview of scientific research of these genera and the prospect for further research. Citation: Bhunjun CS, Chen YJ, Phukhamsakda C, Boekhout T, Groenewald JZ, McKenzie EHC, Francisco EC, Frisvad JC, Groenewald M, Hurdeal VG, Luangsa-ard J, Perrone G, Visagie CM, Bai FY, Błaszkowski J, Braun U, de Souza FA, de Queiroz MB, Dutta AK, Gonkhom D, Goto BT, Guarnaccia V, Hagen F, Houbraken J, Lachance MA, Li JJ, Luo KY, Magurno F, Mongkolsamrit S, Robert V, Roy N, Tibpromma S, Wanasinghe DN, Wang DQ, Wei DP, Zhao CL, Aiphuk W, Ajayi-Oyetunde O, Arantes TD, Araujo JC, Begerow D, Bakhshi M, Barbosa RN, Behrens FH, Bensch K, Bezerra JDP, Bilański P, Bradley CA, Bubner B, Burgess TI, Buyck B, Čadež N, Cai L, Calaça FJS, Campbell LJ, Chaverri P, Chen YY, Chethana KWT, Coetzee B, Costa MM, Chen Q, Custódio FA, Dai YC, Damm U, de Azevedo Santiago ALCM, De Miccolis Angelini RM, Dijksterhuis J, Dissanayake AJ, Doilom M, Dong W, Alvarez-Duarte E, Fischer M, Gajanayake AJ, Gené J, Gomdola D, Gomes AAM, Hausner G, He MQ, Hou L, Iturrieta-González I, Jami F, Jankowiak R, Jayawardena RS, Kandemir H, Kiss L, Kobmoo N, Kowalski T, Landi L, Lin CG, Liu JK, Liu XB, Loizides M, Luangharn T, Maharachchikumbura SSN, Makhathini Mkhwanazi GJ, Manawasinghe IS, Marin-Felix Y, McTaggart AR, Moreau PA, Morozova OV, Mostert L, Osiewacz HD, Pem D, Phookamsak R, Pollastro S, Pordel A, Poyntner C, Phillips AJL, Phonemany M, Promputtha I, Rathnayaka AR, Rodrigues AM, Romanazzi G, Rothmann L, Salgado-Salazar C, Sandoval-Denis M, Saupe SJ, Scholler M, Scott P, Shivas RG, Silar P, Souza-Motta CM, Silva-Filho AGS, Spies CFJ, Stchigel AM, Sterflinger K, Summerbell RC, Svetasheva TY, Takamatsu S, Theelen B, Theodoro RC, Thines M, Thongklang N, Torres R, Turchetti B, van den Brule T, Wang XW, Wartchow F, Welti S, Wijesinghe SN, Wu F, Xu R, Yang ZL, Yilmaz N, Yurkov A, Zhao L, Zhao RL, Zhou N, Hyde KD, Crous PW (2024). What are the 100 most cited fungal genera? Studies in Mycology 108: 1-411. doi: 10.3114/sim.2024.108.01.

Optimizing liquid fermentation for Wolfiporia cocos: gene expression and biosynthesis of pachymic acid and mycelial biomass

Wolfiporia cocos, a versatile fungus acclaimed for its nutritional and therapeutic benefits in Traditional Chinese Medicine, holds immense potential for pharmaceutical and industrial applications. In this study, we aimed to optimize liquid fermentation techniques and culture medium composition to maximize mycelial biomass (MB) yield, pachymic acid (PA) concentration, and overall PA production. Additionally, we investigated the molecular basis of our findings by quantifying the expression levels of genes associated with PA and MB biosynthesis using quantitative real-time polymerase chain reaction. Under the optimized fermentation conditions, significant results were achieved, with maximum MB reaching 6.68 g l-1, PA content peaking at 1.25 mg g-1, and a total PA yield of 4.76 g l-1. Notably, among the four examined genes, squalene monooxygenase, exhibited enhanced expression at 0.06 ratio under the optimized conditions. Furthermore, within the realm of carbohydrate-active enzymes, the glycoside hydrolases 16 family displayed elevated expression levels at 21 ratios, particularly during MB production. This study enhances understanding of genetic mechanism governing MB and PA production in W. cocos, highlighting the roles of squalene monooxygenase and glycoside hydrolases 16 carbohydrate-active enzymes.

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.

Biological characteristics of Cordyceps militaris single mating-type strains

Cordyceps militaris has been extensively cultivated as a model cordyceps species for commercial purposes. Nevertheless, the problems related to strain degeneration and breeding technologies remain unresolved. This study assessed the physiology and fertility traits of six C. militaris strains with distinct origins and characteristics, focusing on single mating-type strains. The results demonstrated that the three identified strains (CMDB01, CMSY01, and CMJB02) were single mating-type possessing only one mating-type gene (MAT1-1). In contrast, the other three strains (CMXF07, CMXF09, and CMMS05) were the dual mating type. The MAT1-1 strains sourced from CMDB01, CMSY01, and CMJB02 consistently produced sporocarps but failed to generate ascospores. However, when paired with MAT1-2 strains, the MAT1-1 strains with slender fruiting bodies and normal morphology were fertile. The hyphal growth rate of single mating-type strains (CMDB01, CMSY01, and CMJB02) typically surpassed that of dual mating-type strains (CMXF07, CMXF09, and CMMS05). The growth rates of MAT1-2 and MAT1-1 strains were proportional to their ratios, such that a single mating-type strain with a higher ratio exhibited an increased growth rate. As C. militaris matured, the adenosine content decreased. In summary, the C. militaris strains that consistently produce sporocarps and have a single mating type are highly promising for production and breeding.

A review on polysaccharide biosynthesis in Cordyceps militaris

Cordyceps militaris (C. militaris) is an edible parasitic fungus with medicinal properties. Its bioactive polysaccharides are structurally diverse and exhibit various metabolic and biological activities, including antitumor, hypoglycemic, antioxidant, hypolipidemic, anti-inflammatory, immunostimulatory, and anti-atherosclerotic effects. These properties make C. militaris-derived polysaccharides a promising candidate for future development. Recent advancements in microbial fermentation technology have enabled successful laboratory cultivation and extraction of these polysaccharides. These polysaccharides are structurally diverse and exhibit various biological activities, such as immunostimulatory, antioxidant, antitumor, hypolipidemic, and anti-atherosclerotic effects. This review aims to summarize the structure and production mechanisms of polysaccharides from C. militaris, covering extraction methods, key genes and pathways involved in biosynthesis, and fermentation factors that influence yield and activity. Furthermore, the future potential and challenges of utilizing polysaccharides in the development of health foods and pharmaceuticals are addressed. This review serves as a valuable reference in the fields of food and medicine, and provides a theoretical foundation for the study of polysaccharides.

Rapid analysis of insecticidal metabolites from the entomopathogenic fungus Beauveria bassiana 331R using UPLC-Q-Orbitrap MS

Beauveria bassiana, a representative entomopathogenic fungus, is increasingly being utilized as an eco-friendly pest management alternative to chemical insecticides. This fungus produces a range of insecticidal secondary metabolites that act as antimicrobial and immunosuppressive agents. However, detailed qualitative and quantitative analysis related to these compounds remains scarce, we developed a method for the rapid analysis of these metabolites. Eight secondary metabolites (bassianin, bassianolide, beauvericin, beauveriolide I, enniatin A, A1, and B, and tenellin) were efficiently extracted when B. bassiana-infected Tenebrio molitor larvae were ground in 70% EtOH extraction solvent and subsequently subjected to ultrasonic treatment for 30 min. The eight metabolites were rapidly and simultaneously analyzed using ultra-performance liquid chromatography-quadrupole-Orbitrap mass spectrometry (UPLC-Q-Orbitrap MS). Bassianolide (20.6-51.1 µg/g) and beauvericin (63.6-109.8 µg/g) were identified as the main metabolites in B. basssiana-infected larvae, indicating that they are likely major toxins of B. bassiana. Validation of the method exhibited recovery rates in the range of 80-115% and precision in the range of 0.1-8.0%, indicating no significant interference from compounds in the matrix. We developed a method to rapidly analyze eight insecticidal metabolites using UPLC-Q-Orbitrap MS. This can be extensively utilized for detecting and producing insecticidal fungal secondary metabolites.

Identification of the chaA and fwA Spore Color Genes of Aspergillus nidulans

Wild-type Aspergillus nidulans asexual spores (conidia) are green due to a pigment that protects the spores against ultraviolet light. The pigment is produced by a biosynthetic pathway, the genes of which are dispersed in the genome. The backbone molecule of the pigment is a polyketide synthesized by a polyketide synthase encoded by the wA gene. If wA is not functional, the conidia are white. The polyketide is modified by a laccase encoded by the yA gene and inactivation of yA in an otherwise wild-type background results in yellow spores. Additional spore color mutations have been isolated and mapped to a locus genetically, but the genes that correspond to these loci have not been determined. Spore color markers have been useful historically, and they remain valuable in the molecular genetics era. One can determine if a transforming fragment has been successfully integrated at the wA or yA locus by simply looking at the color of transformant conidia. The genes of the potentially useful color loci chaA (chartreuse conidia) and fwA (fawn conidia) have not been identified previously. We chose a set of candidate genes for each locus by comparing the assembled genome with the genetic map. By systematically deleting these candidate genes, we identified a cytochrome P450 gene (AN10028) corresponding to chaA. Deletions of this gene result in chartreuse conidia and chartreuse mutations can be complemented in trans by a functional copy of this gene. With fwA, we found that the existing fawn mutation, fwA1, is a deletion of 2241 base pairs that inactivates three genes. By deleting each of these genes, we determined that fwA is AN1088, an EthD domain protein. Deletion of AN1088 results in fawn conidia as expected. Neither deletion of chaA nor fwA restricts growth and both should be valuable target loci for transformations. Combinations of deletions have allowed us to investigate the epistasis relationships of wA, yA, chaA and fwA.

Molecular mechanisms underlying phenotypic degeneration in Cordyceps militaris: insights from transcriptome reanalysis and osmotic stress studies

Phenotypic degeneration in Cordyceps militaris poses a significant concern for producers, yet the mechanisms underlying this phenomenon remain elusive. To address this concern, we isolated two strains that differ in their abilities to form fruiting bodies. Our observations revealed that the degenerated strain lost the capacity to develop fruiting bodies, exhibited limited radial expansion, increased spore density, and elevated intracellular glycerol levels. Transcriptome reanalysis uncovered dysregulation of genes involved in the MAPK signaling pathway in the degenerate strain. Our RT-qPCR results demonstrated reduced expression of sexual development genes, along with upregulation of genes involved in asexual sporulation, glycerol synthesis, and MAPK regulation, when compared to the wild-type strain. Additionally, we discovered that osmotic stress reduced radial growth but increased conidia sporulation and glycerol accumulation in all strains. Furthermore, hyperosmotic stress inhibited fruiting body formation in all neutralized strains. These findings indicate dysregulation of the MAPK signaling pathway, the possibility of the activation of the high-osmolarity glycerol and spore formation modules, as well as the downregulation of the pheromone response and filamentous growth cascades in the degenerate strain. Overall, our study sheds light on the mechanisms underlying Cordyceps militaris degeneration and identifies potential targets for improving cultivation practices.

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.

Integrated Comparative Transcriptome and Weighted Gene Co-Expression Network Analysis Provide Valuable Insights into the Mechanisms of Pinhead Initiation in Chinese Caterpillar Mushroom Ophiocordyceps sinensis (Ascomycota)

The initiation and formation of the "pinhead" is the key node in growth process of Ophiocordyceps sinensis (Chinese Cordyceps). The research on the mechanism of changes in this growth stage is the basis for realizing the industrialization of its artificial cultivation. Clarifying the mechanisms of pinhead initiation is essential for its further application. Here, we performed a comprehensive transcriptome analysis of pinhead initiation process in O. sinensis. Comparative transcriptome analysis revealed remarkable variation in gene expression and enriched pathways at different pinhead initiation stages. Gene co-expression network analysis by WGCNA identified 4 modules highly relevant to different pinhead initiation stages, and 23 hub genes. The biological function analysis and hub gene annotation of these identified modules demonstrated that transmembrane transport and nucleotide excision repair were the topmost enriched in pre-pinhead initiation stage, carbohydrate metabolism and protein glycosylation were specially enriched in pinhead initiation stage, nucleotide binding and DNA metabolic process were over-represented after pinhead stage. These key regulators are mainly involved in carbohydrate metabolism, synthesis of proteins and nucleic acids. This work excavated the candidate pathways and hub genes related to the pinhead initiation stage, which will serve as a reference for realizing the industrialization of artificial cultivation in O. sinensis.

Two new diketopiperazines from the Cordyceps fungus Samsoniella sp. XY4

Two new diketopiperazines, namely samsoniellain A (1) and samsoniellain B (2), together with two known compounds (3, 4) were isolated from Cordyceps fungus Samsoniella sp. XY4. The planar structures of 1 and 2 were determined by HRESIMS, 1D and 2D NMR spectroscopy. The absolute configurations of 1 and 2 were determined by comparison of quantum chemical TDDFT calculated and experimental ECD spectra. Results of antimicrobial activity indicated that compound 2 showed weak bacteriostatic activities against S. typhimurium χ 8956, H. influenza ATCC 10211, MRSA 2024 with the MIC values of 128, 256, and 256 μg ml-1, respectively. This is the first report about secondary metabolites of Samsoniella sp.

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.

Safe-Harbor-Targeted CRISPR/Cas9 System and Cmhyd1 Overexpression Enhances Disease Resistance in Cordyceps militaris

Fungal disease of mushroomCordyceps militaris (CM) caused byCalcarisporium cordycipiticola (CC) is destructive to fruiting body cultivation, resulting in significant economic loss and potential food safety risks. CRISPR/Cas9 genome editing has proven to be a powerful tool for crop improvement but seldom succeeded in mushrooms. Here, the first genomic safe-harbor site, CmSH1 locus, was identified in the CM genome. A safe-harbor-targeted CRISPR/Cas9 system based on an autonomously replicating plasmid was designed to facilitate alien gene integration at the CmSH1 locus. Cmhyd1, one of the hydrophobin genes, was confirmed as a defensive factor against CC infection, and Cmhyd1 overexpression by this system showed enhancement of disease resistance with negligible effect on the agronomic traits of CM. No off-target events and residues of plasmid sequence were tested by PCR and genome resequencing. This study provided the first safe harbor site for genetic manipulations, a safe harbor-targeted CRISPR/Cas9 system, and the first disease-resistant gene-editing breeding system in mushrooms.

Effects of MAT1-2 Spore Ratios on Fruiting Body Formation and Degeneration in the Heterothallic Fungus Cordyceps militaris

The medicinal mushroom Cordyceps militaris is widely exploited in traditional medicine and nutraceuticals in Asian countries. However, fruiting body production in C. militaris is facing degeneration through cultivation batches, and the molecular mechanism of this phenomenon remains unclear. This study showed that fruiting body formation in three different C. militaris strains, namely G12, B12, and HQ1, severely declined after three successive culturing generations using the spore isolation method. PCR analyses revealed that these strains exist as heterokaryons and possess both the mating-type loci, MAT1-1 and MAT1-2. Further, monokaryotic isolates carrying MAT1-1 or MAT1-2 were successfully separated from the fruiting bodies of all three heterokaryotic strains. A spore combination of the MAT1-1 monokaryotic isolate and the MAT1-2 monokaryotic isolate promoted fruiting body formation, while the single monokaryotic isolates could not do that themselves. Notably, we found that changes in ratios of the MAT1-2 spores strongly influenced fruiting body formation in these strains. When the ratios of the MAT1-2 spores increased to more than 15 times compared to the MAT1-1 spores, the fruiting body formation decreased sharply. In contrast, when MAT1-1 spores were increased proportionally, fruiting body formation was only slightly reduced. Our study also proposes a new solution to mitigate the degeneration in the heterokaryotic C. militaris strains caused by successive culturing generations.

A bacterial-like Pictet-Spenglerase drives the evolution of fungi to produce β-carboline glycosides together with separate genes

Diverse β-carboline (βC) alkaloids are produced by microbes, plants, and animals with myriad bioactivities and drug potentials. However, the biosynthetic mechanism of βCs remains largely elusive, especially regarding the hydroxyl and glucosyl modifications of βCs. Here, we report the presence of the bacterial-like Pictet-Spenglerase gene Fcs1 in the entomopathogenic Beauveria fungi that can catalyze the biosynthesis of the βC skeleton. The overexpression of Fcs1 in Beauveria bassiana led to the identification of six βC methyl glycosides, termed bassicarbosides (BCSs) A-F. We verified that the cytochrome P450 (CYP) genes adjacent to Fcs1 cannot oxidize βCs. Alternatively, the separated CYP684B2 family gene Fcs2 was identified to catalyze βC hydroxylation together with its cofactor gene Fcs3. The functional homologue of Fcs2 is only present in the Fcs1-containing fungi and highly similar to the Fcs1-connected yet nonfunctional CYP. Both evolved quicker than those from fungi without Fcs1 homologues. Finally, the paired methyl/glucosyl transferase genes were verified to mediate the production of BCSs from hydroxy-βCs. All these functionally verified genes are located on different chromosomes of Beauveria, which is in contrast to the typical content-clustered feature of fungal biosynthetic gene clusters (BGCs). We also found that the production of BCSs selectively contributed to fungal infection of different insect species. Our findings shed light on the biosynthetic mechanism of βC glycosides, including the identification of a βC hydroxylase. The results of this study also propose an evolving process of fungal BGC formation following the horizontal transfer of a bacterial gene to fungi.

Systematic analyses with genomic and metabolomic insights reveal a new species, Ophiocordyceps indica sp. nov. from treeline area of Indian Western Himalayan region

Ophiocordyceps is a species-rich genus in the order Hypocreales (Sordariomycetes, Ascomycota) depicting a fascinating relationship between microbes and insects. In the present study, a new species, Ophiocordyceps indica sp. nov., is discovered infecting lepidopteran larvae from tree line locations (2,202-2,653 m AMSL) of the Kullu District, Himachal Pradesh, Indian Western Himalayan region, using combinations of morphological and molecular phylogenetic analyses. A phylogeny for Ophiocordyceps based on a combined multigene (nrSSU, nrLSU, tef-1α, and RPB1) dataset is provided, and its taxonomic status within Ophiocordycipitaceae is briefly discussed. Its genome size (~59 Mb) revealed 94% genetic similarity with O. sinensis; however, it differs from other extant Ophiocordyceps species based on morphological characteristics, molecular phylogenetic relationships, and genetic distance. O. indica is identified as the second homothallic species in the family Ophiocordycipitaceae, after O. sinensis. The presence of targeted marker components, viz. nucleosides (2,303.25 μg/g), amino acids (6.15%), mannitol (10.13%), and biological activity data, suggests it to be a new potential source of nutraceutical importance. Data generated around this economically important species will expand our understanding regarding the diversity of Ophiocordyceps-like taxa from new locations, thus providing new research avenues.

Genomic virulence features of Beauveria bassiana as a biocontrol agent for the mountain pine beetle population

BACKGROUND

The mountain pine beetle, Dendroctonus ponderosae, is an irruptive bark beetle that causes extensive mortality to many pine species within the forests of western North America. Driven by climate change and wildfire suppression, a recent mountain pine beetle (MPB) outbreak has spread across more than 18 million hectares, including areas to the east of the Rocky Mountains that comprise populations and species of pines not previously affected. Despite its impacts, there are few tactics available to control MPB populations. Beauveria bassiana is an entomopathogenic fungus used as a biological agent in agriculture and forestry and has potential as a management tactic for the mountain pine beetle population. This work investigates the phenotypic and genomic variation between B. bassiana strains to identify optimal strains against a specific insect.

RESULTS

Using comparative genome and transcriptome analyses of eight B. bassiana isolates, we have identified the genetic basis of virulence, which includes oosporein production. Genes unique to the more virulent strains included functions in biosynthesis of mycotoxins, membrane transporters, and transcription factors. Significant differential expression of genes related to virulence, transmembrane transport, and stress response was identified between the different strains, as well as up to nine-fold upregulation of genes involved in the biosynthesis of oosporein. Differential correlation analysis revealed transcription factors that may be involved in regulating oosporein production.

CONCLUSION

This study provides a foundation for the selection and/or engineering of the most effective strain of B. bassiana for the biological control of mountain pine beetle and other insect pests populations.

The Sordariomycetes: an expanding resource with Big Data for mining in evolutionary genomics and transcriptomics

Advances in genomics and transcriptomics accompanying the rapid accumulation of omics data have provided new tools that have transformed and expanded the traditional concepts of model fungi. Evolutionary genomics and transcriptomics have flourished with the use of classical and newer fungal models that facilitate the study of diverse topics encompassing fungal biology and development. Technological advances have also created the opportunity to obtain and mine large datasets. One such continuously growing dataset is that of the Sordariomycetes, which exhibit a richness of species, ecological diversity, economic importance, and a profound research history on amenable models. Currently, 3,574 species of this class have been sequenced, comprising nearly one-third of the available ascomycete genomes. Among these genomes, multiple representatives of the model genera Fusarium, Neurospora, and Trichoderma are present. In this review, we examine recently published studies and data on the Sordariomycetes that have contributed novel insights to the field of fungal evolution via integrative analyses of the genetic, pathogenic, and other biological characteristics of the fungi. Some of these studies applied ancestral state analysis of gene expression among divergent lineages to infer regulatory network models, identify key genetic elements in fungal sexual development, and investigate the regulation of conidial germination and secondary metabolism. Such multispecies investigations address challenges in the study of fungal evolutionary genomics derived from studies that are often based on limited model genomes and that primarily focus on the aspects of biology driven by knowledge drawn from a few model species. Rapidly accumulating information and expanding capabilities for systems biological analysis of Big Data are setting the stage for the expansion of the concept of model systems from unitary taxonomic species/genera to inclusive clusters of well-studied models that can facilitate both the in-depth study of specific lineages and also investigation of trait diversity across lineages. The Sordariomycetes class, in particular, offers abundant omics data and a large and active global research community. As such, the Sordariomycetes can form a core omics clade, providing a blueprint for the expansion of our knowledge of evolution at the genomic scale in the exciting era of Big Data and artificial intelligence, and serving as a reference for the future analysis of different taxonomic levels within the fungal kingdom.

Genomic Characterization and Functional Description of Beauveria bassiana Isolates from Latin America

Beauveria bassiana is an entomopathogenic fungus used in agriculture as a biological controller worldwide. Despite being a well-studied organism, there are no genomic studies of B. bassiana isolates from Central American and Caribbean countries. This work characterized the functional potential of eight Neotropical isolates and provided an overview of their genomic characteristics, targeting genes associated with pathogenicity, the production of secondary metabolites, and the identification of CAZYmes as tools for future biotechnological applications. In addition, a comparison between these isolates and reference genomes was performed. Differences were observed according to geographical location and the lineages of the B. bassiana complex to which each isolate belonged.

Entomopathogenicity of Ascomycete Fungus Cordyceps militaris on the Cotton Bollworm, Helicoverpa zea (Boddie) (Lepidoptera: Noctuidae)

This study investigated the exposure of the cotton bollworm, Helicoverpa zea (Boddie) (Lepidoptera: Noctuidae) to a novel pathogenic fungal agent historically associated with human medicinal value, a commercial strain of Cordyceps militaris ((L.) Fr.) Vuill. (Hypocreales). A series of comparative studies were conducted to evaluate the efficacy of two different exposure methods using four concentrations (n × 10⁹, n × 10⁸, n × 10⁷, n × 10⁶) of C. militaris, where n × 10⁹ provided a concentration of approximately 420 ± 37 spores per mm² with 398 ± 28 viable spores. Survival of cotton bollworms of all stages was not affected by C. militaris at any concentration 1 d post-exposure. The greatest reduction in survival and highest sporulation rates were observed primarily on or after 7 d post-exposure for early instars (first and second). Significant declines in the survival of early instars were observed for all concentrations at 7 d, and 95% mortality by 10 d, with the exception of the fifth instars that experienced a less severe reduction in survival (35%) when exposed to any concentrations used in the study. Survival of late instars (third to fifth) ranged from 44% to 68% on day 10, while adult survival was near 99% across the duration of the experiment. The relatively narrow range observed for both the lethal concentration and sporulation of second, third, and fifth instar cotton bollworms exposed to the C. militaris strain may demonstrate potential field application for control of larval populations of cotton bollworms.

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.

Complete genome analysis of pathogenic Metschnikowia bicuspidata strain MQ2101 isolated from diseased ridgetail white prawn, Exopalaemon carinicauda

BACKGROUND

Metschnikowia bicuspidata is a pathogenic yesst that can cause disease in many different economic aquatic animal species. In recent years, there was a new disease outbreak in ridgetail white prawn (Exopalaemon carinicauda) in coastal areas of Jiangsu Province China that was referred to as zombie disease by local farmers. The pathogen was first isolated and identified as M. bicuspidata. Although the pathogenicity and pathogenesis of this pathogen in other animals have been reported in some previous studies, research on its molecular mechanisms is still very limited. Therefore, a genome-wide study is necessary to better understand the physiological and pathogenic mechanisms of M. bicuspidata.

RESULT

In this study, we obtained a pathogenic strain, MQ2101, of M. bicuspidata from diseased E. carinicauda and sequenced its whole genome. The size of the whole genome was 15.98 Mb, and it was assembled into 5 scaffolds. The genome contained 3934 coding genes, among which 3899 genes with biological functions were annotated in multiple underlying databases. In KOG database, 2627 genes were annotated, which were categorized into 25 classes including general function prediction only, posttranslational modification, protein turnover, chaperones, and signal transduction mechanisms. In KEGG database, 2493 genes were annotated, which were categorized into five classes, including cellular processes, environmental information processing, genetic information processing, metabolism and organismal systems. In GO database, 2893 genes were annotated, which were mainly classified in cell, cell part, cellular processes and metabolic processes. There were 1055 genes annotated in the PHI database, accounting for 26.81% of the total genome, among which 5 genes were directly related to pathogenicity (identity ≥ 50%), including hsp90, PacC, and PHO84. There were also some genes related to the activity of the yeast itself that could be targeted by antiyeast drugs. Analysis based on the DFVF database showed that strain MQ2101 contained 235 potential virulence genes. BLAST searches in the CAZy database showed that strain MQ2101 may have a more complex carbohydrate metabolism system than other yeasts of the same family. In addition, two gene clusters and 168 putative secretory proteins were predicted in strain MQ2101, and functional analysis showed that some of the secretory proteins may be directly involved in the pathogenesis of the strain. Gene family analysis with five other yeasts revealed that strain MQ2101 has 245 unique gene families, including 274 genes involved in pathogenicity that could serve as potential targets.

CONCLUSION

Genome-wide analysis elucidated the pathogenicity-associated genes of M. bicuspidate while also revealing a complex metabolic mechanism and providing putative targets of action for the development of antiyeast drugs for this pathogen. The obtained whole-genome sequencing data provide an important theoretical basis for transcriptomic, proteomic and metabolic studies of M. bicuspidata and lay a foundation for defining its specific mechanism of host infestation.

Functional genomics and systems biology of Cordyceps species for biotechnological applications

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

Dual Transcriptomics Reveals Interspecific Interactions between the Mycoparasite Calcarisporium cordycipiticola and Its Host Cordyceps militaris

Calcarisporium cordycipiticola is a mycoparasite of the edible fungus Cordyceps militaris, and mycoparasitism causes devastating diseases of mushrooms. In this study, dual-transcriptomic analysis was performed to reveal interspecific interactions between the mycoparasite C. cordycipiticola and its host C. militaris. At 4 and 8 days postinfection (dpi), 2,959 and 2,077 differentially expressed genes (DEGs) of C. cordycipiticola and 914 and 1,548 DEGs of C. militaris were identified compared with the mycelial stage, respectively, indicating that C. cordycipiticola responded more quickly than C. militaris. Lectins of the pathogen may play a role in the recognition of fungal prey. Both Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that primary metabolism was vigorous for the pathogen to colonize the host and that the pathogen's attack substantially altered C. militaris' primary metabolism. C. cordycipiticola upregulated some carbohydrate-active enzyme (CAZyme) genes, including CBM18, GH18, GH16, and GH76, for degrading the host cell wall and defending against host immunity. C. militaris produced excessive reactive oxygen species (ROS) to respond to the infection. The GO term "heme binding" was the only shared term enriched at both stages at 4 and 8 dpi, indicating that iron was important for both the pathogen and the host. The uptake of iron by pathogens through multiple pathways promoted colonization and removed high ROS levels produced by the host. The transcription levels of Cmhsp78, Cmhsp70, and Cmhyd1 in C. militaris responded quickly, and these genes have potential as candidates for the breeding of resistant varieties. This study provides clues for understanding the interactions between a mycoparasite and its mushroom host and will be helpful for the breeding of resistant varieties and disease prevention and control for this edible fungus. IMPORTANCE White mildew disease caused by Calcarisporium cordycipiticola is devastating for the fruiting body cultivation of Cordyceps militaris, a popular and highly valued edible fungus. Here, the pathogenic mechanisms of C. cordycipiticola, the responses of C. militaris to the infection, and the interaction of these two phylogenetically close species were revealed by time course dual-transcriptome profiles. In general, the host C. militaris responds more slowly than the pathogen C. cordycipiticola. For the first time, we found that iron was important for both the mycoparasite and the host. C. cordycipiticola takes up iron by multiple pathways to promote colonization and remove high ROS levels produced by the host. The rapidly responding genes Cmhsp70, Cmhsp78, and Cmhyd1 in C. militaris may have the potential as candidate genes for the breeding of resistant varieties. This study expands our understanding of the mycoparasitic interactions of two species from sister families and will be helpful for the breeding of and disease prevention and control in mushrooms.

Increased Cordycepin Production in Yarrowia lipolytica Using Combinatorial Metabolic Engineering Strategies

As the first nucleoside antibiotic discovered in fungi, cordycepin, with its various biological activities, has wide applications. At present, cordycepin is mainly obtained from the natural fruiting bodies of Cordyceps militaris. However, due to long production periods, low yields, and low extraction efficiency, harvesting cordycepin from natural C. militaris is not ideal, making it difficult to meet market demands. In this study, an engineered Yarrowia lipolytica YlCor-18 strain, constructed by combining metabolic engineering strategies, achieved efficient de novo cordycepin production from glucose. First, the cordycepin biosynthetic pathway derived from C. militaris was introduced into Y. lipolytica. Furthermore, metabolic engineering strategies including promoter, protein, adenosine triphosphate, and precursor engineering were combined to enhance the synthetic ability of engineered strains of cordycepin. Fermentation conditions were also optimized, after which, the production titer and yields of cordycepin in the engineered strain YlCor-18 under fed-batch fermentation were improved to 4362.54 mg/L and 213.85 mg/g, respectively, after 168 h. This study demonstrates the potential of Y. lipolytica as a cell factory for cordycepin synthesis, which will serve as the model for the green biomanufacturing of other nucleoside antibiotics using artificial cell factories.

Essential Role of WetA, but No Role of VosA, in Asexual Development, Conidial Maturation and Insect Pathogenicity of Metarhizium robertsii

Conidial maturation, which is crucial for conidial quality, is controlled by the asexual development activator WetA and the downstream, velvety protein VosA in Aspergillus. Their orthologs have proved functional in conidial quality control of Beauveria bassiana, as seen in Aspergillus, but are functionally unexplored, in Metarhizium robertsii, another hypocrealean insect pathogen. Here, WetA and VosA prove essential and nonessential for M. robertsii's life cycle, respectively. Disruption of wetA increased hyphal sensitivity to oxidative stress and Congo red-induced cell wall stress, but had little impact on radial growth. The ΔwetA mutant was severely compromised in conidiation capacity and conidial quality, which was featured by slower germination, decreased UV resistance, reduced hydrophobicity, and deformed hydrophobin rodlet bundles that were assembled onto conidial coat. The mutant's virulence was greatly attenuated via normal infection due to a blockage of infection-required cellular processes. All examined phenotypes were unaffected for the ΔvosA mutant. Intriguingly, mannitol was much less accumulated in the 7- and 15-day-old cultures of ΔwetA and ΔvosA than of control strains, while accumulated trehalose was not detectable at all, revealing little a link of intracellular polyol accumulation to conidial maturation. Transcriptomic analysis revealed differential regulation of 160 genes (up/down ratio: 72:88) in ΔwetA. These genes were mostly involved in cellular component, biological process, and molecular function but rarely associated with asexual development. Conclusively, WetA plays a relatively conserved role in M. robertsii's spore surface structure, and also a differentiated role in some other cellular processes associated with conidial maturation. VosA is functionally redundant in M. robertsii unlike its ortholog in B. bassiana. IMPORTANCE WetA and VosA regulate conidiation and conidial maturation required for the life cycle of Beauveria bassiana, like they do in Aspergillus, but remain functionally unexplored in Metarhizium robertsii, another hypocrealean pathogen considered to have evolved insect pathogenicity ~130 million years later than B. bassiana. This study reveals a similar role of WetA ortholog in asexual development, conidial maturation, and insect pathogenicity, and also its distinctive role in mediating some other conidial maturation-related cellular events, but has functional redundancy of VosA in M. robertsii. The maturation process vital for conidial quality proves dependent on a role of WetA in spore wall assembly but is independent of its role in intracellular polyol accumulation. Transcriptomic analysis reveals a link of WetA to 160 genes involved in cellular component, biological process, and molecular function. Our study unveils that M. robertsii WetA or VosA is functionally differential or different from those learned in B. bassiana and other ascomycetes.

Chitinase Is Involved in the Fruiting Body Development of Medicinal Fungus Cordyceps militaris

Cordyceps militaris is a famous traditional edible and medicinal fungus in Asia, and its fruiting body has rich medicinal value. The molecular mechanism of fruiting body development is still not well understood in C. militaris. In this study, phylogenetically analysis and protein domains prediction of the 14 putative chitinases were performed. The transcription level and enzyme activity of chitinase were significant increased during fruiting body development of C. militaris. Then, two chitinase genes (Chi1 and Chi4) were selected to construct gene silencing strain by RNA interference. When Chi1 and Chi4 genes were knockdown, the differentiation of the primordium was blocked, and the number of fruiting body was significantly decreased approximately by 50% compared to wild-type (WT) strain. The length of the single mature fruiting body was shortened by 27% and 38% in Chi1- and Chi4-silenced strains, respectively. In addition, the chitin content and cell wall thickness were significantly increased in Chi1- and Chi4-silenced strains. These results provide new insights into the biological functions of chitinase in fruiting body development of C. militaris.

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.

Efficient disruption of CmHk1 using CRISPR/Cas9 ribonucleoprotein delivery in Cordyceps militaris

Cordyceps militaris, an entomopathogenic ascomycete, produces edible medicinal mushrooms known to have medicinal and therapeutic functions. To develop the genetic transformation system in C. militaris, green fluorescent protein (GFP) mutants of C. militaris were generated by PEG-mediated protoplast transformation. The CRISPR/Cas9 ribonucleoprotein (RNP) targeting the class III histidine kinase of C. militaris (CmHk1) was then delivered into protoplasts of C. militaris through the transformation system. Mutations induced by the RNP in selected mutants were detected: 1 nt deletion (6 mutants), 3 nt deletion with substitution of 1 nt (1 mutant), insertion of 85 nts (1 mutant), 41 nts (2 mutants), and 35 nts (5 mutants). An in vitro sensitivity assay of the mutants indicated that knockout of CmHk1 reduced sensitivity to two fungicides, iprodione and fludioxonil, but increased sensitivity to osmotic stresses compared to the wild type. Summing up, the CRISPR/Cas9 RNP delivery system was successfully developed, and our results revealed that CmHk1 was involved in the fungicide resistance and osmotic stress in 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.

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

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

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.

De Novo Assembly Transcriptome Analysis Reveals the Preliminary Molecular Mechanism of Primordium Formation in Pleurotus tuoliensis

Primordium formation is extremely important for yield of Pleurotus tuoliensis. However, the molecular mechanism underlying primordium formation is largely unknown. This study investigated the transcriptional properties during primordium formation of P. tuoliensis by comparing transcriptome. Clean reads were assembled into 57,075 transcripts and 6874 unigenes. A total of 1397 differentially expressed genes were identified (26 DEGs altered in all stages). GO and KEGG enrichment analysis showed that these DEGs were involved in “oxidoreductase activity”, “glycolysis/gluconeogenesis”, “MAPK signaling pathways”, and “ribosomes”. Our results support further understanding of the transcriptional changes and molecular processes underlying primordium formation and differentiation of P. tuoliensis.

Multiplex gene precise editing and large DNA fragment deletion by the CRISPR-Cas9-TRAMA system in edible mushroom Cordyceps militaris

The medicinal mushroom Cordyceps militaris contains abundant valuable bioactive ingredients that have attracted a great deal of attention in the pharmaceutical and cosmetic industries. However, the development of this valuable mushroom faces the obstacle of lacking powerful genomic engineering tools. Here, by excavating the endogenous tRNA-processed element, introducing the extrachromosomal plasmid and alongside with homologous template, we develop a marker-free CRISPR-Cas9-TRAMA genomic editing system to achieve the multiplex gene precise editing and large synthetic cluster deletion in C. militaris. We further operated editing in the synthetases of cordycepin and ergothioneine to demonstrate the application of Cas9-TRAMA system in protein modification, promoter strength evaluation and 10 kb metabolic synthetic cluster deletion. The Cas9-TRAMA system provides a scalable method for excavating the valuable metabolic resource of medicinal mushrooms and constructing a mystical cellular pathway to elucidate the complex cell behaviours of the edible mushroom.