Functional convergence and divergence of mating-type genes fulfilling in Cordyceps militaris

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.

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.

Ethanolic extract from fruiting bodies of Cordyceps militaris HL8 exhibits cytotoxic activities against cancer cells, skin pathogenic yeasts, and postharvest pathogen Penicillium digitatum

Cordyceps militaris is a well-known medicinal mushroom in Asian countries. This edible fungus has been widely exploited for traditional medicine and functional food production. C. militaris is a heterothallic fungus that requires both the mating-type loci, MAT1-1 and MAT1-2, for fruiting body formation. However, recent studies also indicated two groups of C. militaris, including monokaryotic strains carrying only MAT1-1 in their genomes and heterokaryotic strains harboring both MAT1-1 and MAT1-2. These strain groups are able to produce fruiting bodies under suitable cultivating conditions. In previous work, we showed that monokaryotic strains are more stable than heterokaryotic strains in fruiting body formation through successive culturing generations. In this study, we report a high cordycepin-producing monokaryotic C. militaris strain (HL8) collected in Vietnam. This strain could form normal fruiting bodies with high biological efficiency and contain a cordycepin content of 14.43 mg/g lyophilized fruiting body biomass. The ethanol extraction of the HL8 fruiting bodies resulted in a crude extract with a cordycepin content of 69.15 mg/g. Assays of cytotoxic activity on six human cancer cell lines showed that the extract inhibited the growth of all these cell lines with the IC50 values of 6.41-11.51 µg/mL. Notably, the extract significantly reduced cell proliferation and promoted apoptosis of breast cancer cells. Furthermore, the extract also exhibited strong antifungal activity against Malassezia skin yeasts and the citrus postharvest pathogen Penicillium digitatum. Our work provides a promising monokaryotic C. militaris strain as a bioresource for medicine, cosmetics, and fruit preservation.

The MAT1 locus is required for microconidia-mediated sexual fertility in the rice blast fungus

Rice blast fungus (Pyricularia oryzae) is a heterothallic ascomycete that causes the most destructive disease in cultivated rice worldwide. This fungus reproduces sexually and asexually, and its mating type is determined by the MAT1 locus, MAT1-1 or MAT1-2. Interestingly, most rice-infecting field isolates show a loss of female fertility, but the MAT1 locus is highly conserved in female-sterile isolates. In this study, we performed a functional analysis of MAT1 using the CRISPR/Cas9 system in female- and male-fertile isolates and female-sterile (male-fertile) isolates. Consistent with a previous report, MAT1 was essential for sexual reproduction but not for asexual reproduction. Meanwhile, deletion mutants of MAT1-1-1, MAT1-1-2, and MAT1-1-3 exhibited phenotypes different from those of other previously described isolates, suggesting that the function of MAT1-1 genes and/or their target genes in sexual reproduction differs among strains or isolates. The MAT1 genes, excluding MAT1-2-6, retained their functions even in female-sterile isolates, and deletion mutants lead to loss or reduction of male fertility. Although MAT1 deletion did not affect microconidia (spermatia) production, microconidia derived from the mutants could not induce perithecia formation. These results indicated that MAT1 is required for microconidia-mediated male fertility in addition to female fertility in P. oryzae .

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.

Isolation of strains and their genome sequencing to analyze the mating system of Ophiocordyceps robertsii

The fungal genus Ophiocordyceps contains a number of insect pathogens. One of the best known of these is Ophiocordyceps sinensis, which is used in Chinese medicine and its overharvesting threatens sustainability; hence, alternative species are being sought. Ophiocordyceps robertsii, found in Australia and New Zealand, has been proposed to be a close relative to O. sinensis, but little is known about this species despite being also of historical significance. Here, O. robertsii strains were isolated into culture and high coverage draft genome sequences obtained and analyzed. This species has a large genome expansion, as also occurred in O. sinensis. The mating type locus was characterized, indicating a heterothallic arrangement whereby each strain has an idiomorphic region of two (MAT1-2-1, MAT1-2-2) or three (MAT1-1-1, MAT1-1-2, MAT1-1-3) genes flanked by the conserved APN2 and SLA2 genes. These resources provide a new opportunity for understanding the evolution of the expanded genome in the homothallic species O. sinensis, as well as capabilities to explore the pharmaceutical potential in a species endemic to Australia and New Zealand.

Organization and Unconventional Integration of the Mating-Type Loci in Morchella Species

True morels (Morchella spp.) are a group of delicious fungi in high demand worldwide, and some species of morels have been successfully cultivated in recent years. To better understand the sexual reproductive mechanisms of these fungi, we characterized the structure of the mating-type loci from ten morel species, and seven of them were obtained using long-range PCR amplification. Among the studied species, eight were heterothallic, two were homothallic, and four types of composition were observed in the MAT loci. In three of the five black morel species, the MAT1-1-1, MAT1-1-10, and MAT1-1-11 genes were in the MAT1-1 idiomorph, and only the MAT1-2-1 gene was in the MAT1-2 idiomorph, while an integration event occurred in the other two species and resulted in the importation of the MAT1-1-11 gene into the MAT1-2 idiomorph and survival as a truncated fragment in the MAT1-1 idiomorph. However, the MAT1-1-11 gene was not available in the four yellow morels and one blushing morel species. M. rufobrunnea, a representative species of the earliest diverging branch of true morels, along with another yellow morel Mes-15, were confirmed to be homothallic, and the MAT1-1-1, MAT1-1-10, and MAT1-2-1 genes were arranged in a tandem array. Therefore, we hypothesized that homothallism should be the ancestral reproductive state in Morchella. RT-PCR analyses revealed that four mating genes could be constitutively expressed, while the MAT1-1-10 gene underwent alternative splicing to produce different splice variants.

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.

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.

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.

Identification of a key G-protein coupled receptor in mediating appressorium formation and fungal virulence against insects

Fungal G-protein coupled receptors (GPCRs) play essential roles in sensing environmental cues including host signals. The study of GPCR in mediating fungus-insect interactions is still limited. Here we report the evolution of GPCR genes encoded in the entomopathogenic Metarhizium species and found the expansion of Pth11-like GPCRs in the generalist species with a wide host range. By deletion of ten candidate genes MrGpr1-MrGpr10 selected from the six obtained subfamilies in the generalist M. robertsii, we found that each of them played a varied level of roles in mediating appressorium formation. In particular, deletion of MrGpr8 resulted in the failure of appressorium formation on different substrates and the loss of virulence during topical infection of insects but not during injection assays when compared with the wild-type (WT) strain. Further analysis revealed that disruption of MrGpr8 substantially impaired the nucleus translocation of the mitogen-activated protein kinase (MAPK) Mero-Fus3 but not the MAPK Mero-Slt2 during appressorium formation. We also found that the defect of AMrGpr8 could not be rescued with the addition of cyclic AMP for appressorium formation. Relative to the WT, differential expression of the selected genes have also been detected in AMrGpr8. The results of this study may benefit the understanding of fungus-interactions mediated by GPCRs.

Two mating-type genes MAT1-1-1 and MAT1-1-2 with significant functions in conidiation, stress response, sexual development, and pathogenicity of rice false smut fungus Villosiclava virens

Rice false smut caused by Villosiclava virens is one of the destructive diseases on panicles of rice. Sexual development of V. virens, controlled by mating-type locus, plays an important role in the prevalence of rice false smut and genetic diversity of the pathogen. However, how the mating-type genes mediate sexual development of the V. virens remains largely unknown. In this study, we characterized the two mating-type genes, MAT1-1-1 and MAT1-1-2, in V. virens. MAT1-1-1 knockout mutant showed defects in hyphal growth, conidia morphogenesis, sexual development, and increase in the tolerance to salt and osmotic stress. Targeted deletion of MAT1-1-2 not only impaired the sclerotia formation and pathogenicity of V. virens, but also reduced the production of conidia. The MAT1-1-2 mutant showed increases in tolerance to salt and hydrogen peroxide stress, but decreases in tolerance to osmotic stress. Yeast two-hybrid assay showed that MAT1-1-1 interacted with MAT1-1-2, indicating that those proteins might form a complex to regulate sexual development. In addition, MAT1-1-1 localized in the nucleus, and MAT1-1-2 localized in the cytoplasm. Collectively, our results demonstrate that MAT1-1-1 and MAT1-1-2 play important roles in the conidiation, stress response, sexual development, and pathogenicity of V. virens, thus providing new insights into the function of mating-type gene.

Bioactive Metabolites and Potential Mycotoxins Produced by Cordyceps Fungi: A Review of Safety

Ascomycete Cordyceps fungi such as C. militaris, C. cicadae, and C. guangdongensis have been mass produced on artificial media either as food supplements or health additives while the byproducts of culture substrates are largely used as animal feed. The safety concerns associated with the daily consumption of Cordyceps fungi or related products are still being debated. On the one hand, the known compounds from these fungi such as adenosine analogs cordycepin and pentostatin have demonstrated different beneficial or pharmaceutical activities but also dose-dependent cytotoxicities, neurological toxicities and or toxicological effects in humans and animals. On the other hand, the possibility of mycotoxin production by Cordyceps fungi has not been completely ruled out. In contrast to a few metabolites identified, an array of biosynthetic gene clusters (BGCs) are encoded in each genome of these fungi with the potential to produce a plethora of as yet unknown secondary metabolites. Conservation analysis of BGCs suggests that mycotoxin analogs of PR-toxin and trichothecenes might be produced by Cordyceps fungi. Future elucidation of the compounds produced by these functionally unknown BGCs, and in-depth assessments of metabolite bioactivity and chemical safety, will not only facilitate the safe use of Cordyceps fungi as human food or alternative medicine, but will also benefit the use of mass production byproducts as animal feed. To corroborate the long record of use as a traditional medicine, future efforts will also benefit the exploration of Cordyceps fungi for pharmaceutical purposes.

Diverse function and regulation of CmSnf1 in entomopathogenic fungus Cordyceps militaris

SNF1/AMPK protein kinases play important roles in fungal development and activation of catabolite-repressed genes. In this study, we characterized the role of SNF1 ortholog in Cordyceps militaris (CmSnf1). The vegetative growth of a CmSnf1 deletion mutant was (ΔCmSnf1) reduced by 42.2% with arabinose as a sole carbon source. Most strikingly, the ΔCmSnf1 produced only a few conidia and exhibited delayed conidial germination. We found that CmSnf1 was necessary for mycelium to penetrate the insect cuticle to form the fruiting body on silkworm pupae, consistent with the down-regulation of chitinase- and protease-encoding genes in ΔCmSnf1. However, cordycepin content increased by more than 7 times in culture supernatants. Correspondingly, the relative expression levels of cordycepin gene cluster members were also elevated. In particular, the expression of cns4 associated with cordycepin transfer was up-regulated >10-fold. Furthermore, transcriptional analysis showed that CmSnf1 regulated the expression of genes involved in cell autophagy and oxidative stress tolerance. We speculated that under environmental stress, both the ATG and SNF1 pathways might collaborate to sustain adverse environments. Our study provides an initial framework to probe the diverse function and regulation of CmSnf1 in C. militaris, which will shed more light on the direction of molecular improvement of medicinal fungi.

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.

Cloning and analysis of Ophiocordyceps xuefengensis mating type (MAT) loci

The entomopathogenic fungus Ophiocordyceps xuefengensis, a recently described species and identified as the sister taxon of Ophiocordyceps sinensis, is a desirable alternative to O. sinensis. The mating systems of fungi play a vitally important role in the regulation of sexual reproduction and evolution, but the mating type loci of O. xuefengensis were completely unknown. In this study, the mating systems of O. xuefengensis were analyzed. The conserved α-box region of the MAT1-1-1 and HMG-box of MAT1-2-1 were successfully obtained by PCR amplification. The distribution of both mating types in different tissues of wild and cultivated O. xuefengensis growth was detected and analyzed. The results showed that the asci always harbored both mating types, whereas the sclerotium, the stipe and each isolated strain of wild O. xuefengensis always had only one idiomorph, either MAT1-1 or MAT1-2, which confirmed that O. xuefengensis is heterothallic. The MAT1-1 locus of O. xuefengensis harbors MAT1-1-1, MAT1-1-2 and MAT1-1-3, and MAT1-2 contains the MAT1-2-1 gene. Southern blot analysis showed the MAT-1-1-1 and MAT-1-2-1 genes were single-copy in O. xuefengensis. These results will help to understand its life cycle and support artificial cultivation of O. xuefengensis.

Genome-wide analysis of DNA methylation in subcultured Cordyceps militaris

The entomopathogenic mushroom Cordyceps militaris is a storehouse of various medicinal compounds and pharmacological effects. However, the high frequency of strain degeneration during subculture and preservation severely limits the large-scale production of C. militaris. DNA methylation is an important epigenomic modification involved in gene regulation. In this study, we used bisulfite sequencing for DNA methylation profiling of wild-type and mutant C. militaris. The differentially methylated regions (DMRs) of the two types were analyzed using Gene Ontology (GO) clustering and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. DNA methylation levels of the wild-type and mutant-type C. militaris were 0.48% and 0.56%, respectively. Methylation appeared at CHH dinucleotides in 58.62% and 58.20% of all methylated cytosine sites in the wild and mutant types, respectively. In all, 188 DMRs were identified from the wild and mutant types. Most of the DMRs ranged from 200 to 350 bp in length. KEGG pathways of the expression of DMR-related genes, which are involved in pyruvate metabolism, glycerophospholipid metabolism, DNA replication, and N-glycan biosynthesis. This contributes to the knowledge and understanding of the possible mechanisms of C. militaris strain degeneration.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Omics data reveal the unusual asexual-fruiting nature and secondary metabolic potentials of the medicinal fungus Cordyceps cicadae

Background

Ascomycete Cordyceps species have been using as valued traditional Chinese medicines. Particularly, the fruiting bodies of Cordyceps cicadae (syn. Isaria cicadae) have long been utilized for the treatment of chronic kidney disease. However, the genetics and bioactive chemicals in this fungus have been largely unexplored.

Results

In this study, we performed comprehensive omics analyses of C. cicadae, and found that, in contrast to other Cordyceps fungi, C. cicadae produces asexual fruiting bodies with the production of conidial spores instead of the meiotic ascospores. Genome sequencing and comparative genomic analysis indicate that the protein families encoded by C. cicadae are typical of entomopathogenic fungi, including the expansion of proteases and chitinases for targeting insect hosts. Interestingly, we found that the MAT1-2 mating-type locus of the sequenced strain contains an abnormally truncated MAT1-1-1 gene. Gene deletions revealed that asexual fruiting of C. cicadae is independent of the MAT locus control. RNA-seq transcriptome data also indicate that, compared to growth in a liquid culture, the putative genes involved in mating and meiosis processes were not up-regulated during fungal fruiting, further supporting asexual reproduction in this fungus. The genome of C. cicadae encodes an array of conservative and divergent gene clusters for secondary metabolisms. Based on our analysis, the production of known carcinogenic metabolites by this fungus could be potentially precluded. However, the confirmed production of oosporein raises health concerns about the frequent consumption of fungal fruiting bodies.

Conclusions

The results of this study expand our knowledge of fungal genetics that asexual fruiting can occur independent of the MAT locus control. The obtained genomic and metabolomic data will benefit future investigations of this fungus for medicinal uses.

Electronic supplementary material

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