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.
Keywords: Analysis; Bisulfite sequencing; Cordyceps militaris; DNA methylation; Degenerative.