Genome-wide analysis of DNA methylation in the sexual stage of the insect pathogenic fungus Cordyceps militaris

Epigenetic Regulation in Insect-Microbe Interactions

Insects have evolved diverse interactions with a variety of microbes, such as pathogenic fungi, bacteria, and viruses. The immune responses of insect hosts, along with the dynamic infection process of microbes in response to the changing host environment and defenses, require rapid and fine-tuned regulation of gene expression programs. Epigenetic mechanisms, including DNA methylation, histone modifications, and noncoding RNA regulation, play important roles in regulating the expression of genes involved in insect immunity and microbial pathogenicity. This review highlights recent discoveries and insights into epigenetic regulatory mechanisms that modulate insect-microbe interactions. A deeper understanding of these regulatory mechanisms underlying insect-microbe interactions holds promise for the development of novel strategies for biological control of insect pests and mitigation of vector-borne diseases.

Epigenetic Regulation of Fungal Secondary Metabolism

Secondary metabolism is one of the important mechanisms by which fungi adapt to their living environment and promote survival and reproduction. Recent studies have shown that epigenetic regulation, such as DNA methylation, histone modifications, and non-coding RNAs, plays key roles in fungal secondary metabolism and affect fungal growth, survival, and pathogenicity. This review describes recent advances in the study of epigenetic regulation of fungal secondary metabolism. We discuss the way in which epigenetic markers respond to environmental changes and stimulate the production of biologically active compounds by fungi, and the feasibility of these new findings applied to develop new antifungal strategies and optimize secondary metabolism. In addition, we have deliberated on possible future directions of research in this field. A deeper understanding of epigenetic regulatory networks is a key focus for future research.

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.

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.

Whole-genome DNA methylome analysis of different developmental stages of the entomopathogenic fungus Beauveria bassiana NCHU-157 by nanopore sequencing

The entomopathogenic fungus (EPF), Beauveria bassiana, is an important and commonly used EPF for microbial control. However, the role of DNA methylation has not been thoroughly studied. Therefore, the whole genomic DNA methylome of one promising EPF isolate, B. bassiana NCHU-157 (Bb-NCHU-157), was investigated by Oxford Nanopore Technologies (ONT). First, the whole genome of Bb-NCHU-157 was sequenced by next-generation sequencing (NGS) and ONT. The genome of Bb-NCHU-157 contains 16 contigs with 34.19 Mb and 50% GC content, which are composed of 10,848 putative protein-coding genes. Two putative DNA methyltransferases (DNMTs) were found, including Dim-2 and C-5 cytosine-specific DNA methylases. Both DNMTs showed higher expression levels in the mycelium stage than in the conidia stage, indicating that development of DNA methylation in Bb-NCHU-157 might occur in the mycelium stage. The global methylation level of the mycelium stage (5 mC = 4.56%, CG = 3.33%, CHG = 0.74%, CHH = 0.49%) was higher than that of the conidial stage (5 mC = 2.99%, CG = 1.99%, CHG = 0.63%, CHH = 0.37%) in both the gene and transposable element (TE) regions. Furthermore, the TE regions showed higher methylation frequencies than the gene regions, especially for CHH site methylation, suggesting regulation of genomic stabilization during mycelium development. In the gene regions, high methylation frequencies were found around the transcription start site (TSS) and transcription end site (TES). Moreover, CG and CHG methylation mainly occur in the promoter and intergenic regions, while CHH methylation occurs in the TE region. Among the methylated regions, 371, 661, and 756 differentially DNA methylated regions (DMRs) were hypermethylated in the mycelium in CG, CHG, and CHH, while only 13 and 7 DMRs were hypomethylated in the mycelium in CHG, and CHH, respectively. Genes located in the DMR shared the GO terms, DNA binding (GO: 0003677), and sequence-specific DNA binding (GO: 0043565) for hypermethylation in the mycelium, suggesting that methylation might regulate gene expression from the initial process. Evaluation of the DNA methylome in Bb-NCHU-157 by ONT provided new insight into this field. These data will be further validated, and epigenetic regulation during the development of B. bassiana will be explored.

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.

Three putative DNA methyltransferases of Verticillium dahliae differentially contribute to DNA methylation that is dispensable for growth, development and virulence

BACKGROUND

DNA methylation is an important epigenetic control mechanism that in many fungi is restricted to genomic regions containing transposable elements (TEs). Two DNA methyltransferases, Dim2 and Dnmt5, are known to perform methylation at cytosines in fungi. While most ascomycete fungi encode both Dim2 and Dnmt5, only few functional studies have been performed in species containing both.

METHODS

In this study, we report functional analysis of both Dim2 and Dnmt5 in the plant pathogenic fungus Verticillium dahliae.

RESULTS

Our results show that Dim2, but not Dnmt5 or the putative sexual-cycle-related DNA methyltransferase Rid, is responsible for the majority of DNA methylation under the tested conditions. Single or double DNA methyltransferase mutants did not show altered development, virulence, or transcription of genes or TEs. In contrast, Hp1 and Dim5 mutants that are impacted in chromatin-associated processes upstream of DNA methylation are severely affected in development and virulence and display transcriptional reprogramming in specific hypervariable genomic regions (so-called adaptive genomic regions) that contain genes associated with host colonization. As these adaptive genomic regions are largely devoid of DNA methylation and of Hp1- and Dim5-associated heterochromatin, the differential transcription is likely caused by pleiotropic effects rather than by differential DNA methylation.

CONCLUSION

Overall, our study suggests that Dim2 is the main DNA methyltransferase in V. dahliae and, in conjunction with work on other fungi, is likely the main active DNMT in ascomycetes, irrespective of Dnmt5 presence. We speculate that Dnmt5 and Rid act under specific, presently enigmatic, conditions or, alternatively, act in DNA-associated processes other than DNA methylation.

Deep sequencing of the Sanghuangporus vaninii transcriptome reveals dynamic landscapes of candidate genes involved in the biosynthesis of active compounds

The medicinal fungus Sanghuang produces diverse bioactive compounds and is widely used in Asian countries. However, little is known about the genes and pathways involved in the biosynthesis of these active compounds. Based on our previous study providing Sanghuangporus vaninii genomic information, the transcriptomes of MY (mycelium), OY (1-year-old fruiting bodies), and TY (3-year-old fruiting bodies) were determined in this study. A significant number of genes (4774) were up- or downregulated between mycelium and fruiting bodies, but only 1422 differentially expressed genes were detected between OY and TY. 138 genes encoding P450s were identified in the fungal genome and grouped into 25 P450 families; more than 64% (88) of the genes were significantly differentially expressed between the mycelium and fruiting body, suggesting that these P450s are involved in fungal sexual development. Importantly, the expression of genes involved in bioactive compound (triterpenoids, polysaccharides, and flavonoids) biosynthesis in asexual (cultured with solid and liquid media) and sexual stages was explored and combined with transcriptome and quantitative PCR analyses. More genes involved in the biosynthesis of bioactive compounds were expressed more highly in mycelium than in fruiting bodies under liquid medium culture compared with solid medium culture, which was consistent with the yields of different bioactive compounds, suggesting that liquid fermentation of S. vaninii Kangneng can be used to obtain these bioactive compounds. A comprehensive understanding of the genomic information of S. vaninii will facilitate its potential use in pharmacological and industrial applications.

Distinct Roles of Two DNA Methyltransferases from Cryphonectria parasitica in Fungal Virulence, Responses to Hypovirus Infection, and Viral Clearance

Two DNA methyltransferase (DNMTase) genes from Cryphonectria parasitica have been previously identified as CpDmt1 and CpDmt2, which are orthologous to rid and dim-2 of Neurospora crassa, respectively. While global changes in DNA methylation have been associated with fungal sectorization and CpDmt1 but not CpDmt2 has been implicated in the sporadic sectorization, the present study continues to investigate the biological functions of both DNMTase genes. Transcription of both DNMTases is regulated in response to infection with the Cryphonectria hypovirus 1 (CHV1-EP713). CpDmt1 is upregulated and CpDmt2 is downregulated by CHV1 infection. Conidium production and response to heat stress are affected only by mutation of CpDmt1, not by CpDmt2 mutation. Significant changes in virulence are observed in opposite directions; i.e., the CpDmt1-null mutant is hypervirulent, while the CpDmt2-null mutant is hypovirulent. Compared to the CHV1-infected wild type, CHV1-transferred single and double mutants show severe growth retardation: the colony size is less than 10% that of the parental virus-free null mutants, and their titers of transferred CHV1 are higher than that of the wild type, implying that no defect in viral replication occurs. However, as cultivation proceeds, spontaneous viral clearance is observed in hypovirus-infected colonies of the null mutants, which has never been reported in this fungus-virus interaction. This study demonstrates that both DNMTases are significant factors in fungal development and virulence. Each fungal DNMTase affects fungal biology in both common and separate ways. In addition, both genes are essential to the antiviral responses, including viral clearance which depends on their mutations.

Comparative analysis of genome-wide DNA methylation in Neurospora

DNA methylation is an epigenetic mark that plays an important role in genetic regulation in eukaryotes. Major progress has been made in dissecting the molecular pathways that regulate DNA methylation. Yet, little is known about DNA methylation variation over evolutionary time. Here we present an investigation of the variation of DNA methylation and transposable element (TE) content in species of the filamentous ascomycetes Neurospora. We generated genome-wide DNA methylation data at single-base resolution, together with genomic TE content and gene expression data, of 10 individuals representing five closely related Neurospora species. We found that the methylation levels were low (ranging from 1.3% to 2.5%) and varied among the genomes in a species-specific way. Furthermore, we found that the TEs over 400 bp long were targeted by DNA methylation, and in all genomes, high methylation correlated with low GC, confirming a conserved link between DNA methylation and Repeat Induced Point (RIP) mutations in this group of fungi. Both TE content and DNA methylation pattern showed phylogenetic signal, and the species with the highest TE load (N. crassa) also exhibited the highest methylation level per TE. Our results suggest that DNA methylation is an evolvable trait and indicate that the genomes of Neurospora are shaped by an evolutionary arms race between TEs and host defence.

DNA methyltransferase implicated in the recovery of conidiation, through successive plant passages, in phenotypically degenerated Metarhizium

Metarhizium robertsii is a fungus with two lifestyles; it is a plant root symbiont and an insect pathogen. A spontaneously phenotypically degenerated strain of M. robertsii strain ARSEF 2575 (M. robertsii lc-2575; lc = low conidiation) showed a reduction in conidiation and fungal virulence after successive subculturing on agar medium. In order to recover conidiation, we experimentally passaged M. robertsii lc-2575 through plant (soldier bean and switchgrass) root or insect (Galleria mellonella) larvae. After five passages, the resultant strains had significantly increased conidial yields on agar and increased virulence in insect bioassays. Concomitantly, DNA methyltransferase, MrDIM-2 expression was downregulated in BR5 (a strain after 5 bean root passages) and isolates after switchgrass and insect passages. Bisulfite sequencing showed little difference in overall genomic DNA methylation levels (~ 0.37%) between M. robertsii lc-2575 and BR5. However, a finer comparison of the different methylated regions (DMRs) showed that DMRs of BR5 were more abundant in the intergenic regions (69.32%) compared with that of M. robertsii lc-2575 (33.33%). The addition of DNA methyltransferase inhibitor, 5-azacytidine, to agar supported the role of DNA methyltransferases and resulted in an increase in conidiation of M. robertsii lc-2575. Differential gene expression was observed in selected DMRs in BR5 when compared with M. robertsii lc-2575. Here we implicated epigenetic regulation in the recovery of conidiation through the effects of DNA methyltransferase and that plant passage could be used as a method to recover fungal conidiation and virulence in a phenotypically degenerated M. robertsii. KEY POINTS: • Passage of Metarhizium through plant root or insect results in increased conidiation. • DNA methyltransferase is downregulated after host passage. • Bisulfite sequencing identified potentially methylated genes involved in conidiation.

Epigenetic manipulation of filamentous fungi for biotechnological applications: a systematic review

The study of the epigenetic regulation of gene function has reached pivotal importance in life sciences in the last decades. The mechanisms and effects of processes such as DNA methylation, histone posttranslational modifications and non-coding RNAs, as well as their impact on chromatin structure and dynamics, are clearly involved in physiology homeostasis in plants, animals and microorganisms. In the fungal kingdom, studies on the model yeasts Saccharomyces cerevisiae and Schizosaccharomyces pombe contributed enormously to the elucidation of the eukaryote epigenetic landscape. Epigenetic regulation plays a central role in the expression of virulence attributes of human pathogens such as Candida albicans. In this article, we review the most recent studies on the effects of drugs capable of altering epigenetic states and on the impact of chromatin structure-related genes deletion in filamentous fungi. Emphasis is given on plant and insect pathogens, endophytes, secondary metabolites and cellulases/xylanases producing species.

The Pattern and Function of DNA Methylation in Fungal Plant Pathogens

To successfully infect plants and trigger disease, fungal plant pathogens use various strategies that are dependent on characteristics of their biology and genomes. Although pathogenic fungi are different from animals and plants in the genomic heritability, sequence feature, and epigenetic modification, an increasing number of phytopathogenic fungi have been demonstrated to share DNA methyltransferases (MTases) responsible for DNA methylation with animals and plants. Fungal plant pathogens predominantly possess four types of DNA MTase homologs, including DIM-2, DNMT1, DNMT5, and RID. Numerous studies have indicated that DNA methylation in phytopathogenic fungi mainly distributes in transposable elements (TEs), gene promoter regions, and the repetitive DNA sequences. As an important and heritable epigenetic modification, DNA methylation is associated with silencing of gene expression and transposon, and it is responsible for a wide range of biological phenomena in fungi. This review highlights the relevant reports and insights into the important roles of DNA methylation in the modulation of development, pathogenicity, and secondary metabolism of fungal plant pathogens. Recent evidences prove that there are massive links between DNA and histone methylation in fungi, and they commonly regulate fungal development and mycotoxin biosynthesis.

Revealing mitogenome-wide DNA methylation and RNA editing of three Ascomycotina fungi using SMRT sequencing

Beauveria bassiana, Cordyceps militaris and Ophiocordyceps sinensis (Ascomycotina) are traditional Chinese medicines. Here, mitogenomes of these three Ascomycotina fungi were sequenced and de-novo assembled using single-molecule real-time sequencing. The results showed that their complete mitogenomes were 31,258, 31,854 and 157,584 bp, respectively, with sequencing depth approximately 278,760×, 326,283× and 69,385×. Types of repeat sequences were mainly (AA)n, (AAT)n, (TA)n and (TATT)n. DNA methylation motifs were revealed in DNA modifications of these three fungi. We discovered new models of RNA editing through analysis of transcriptomes from B. bassiana and C. militaris. These data lay a solid foundation for further genetic and biological studies about these three fungi, especially for elucidating the mitogenome evolution and exploring the regulatory mechanism of adapting environment.

Transcriptome Analysis Reveals the Molecular Mechanisms Underlying Adenosine Biosynthesis in Anamorph Strain of Caterpillar Fungus

Caterpillar fungus is a well-known fungal Chinese medicine. To reveal molecular changes during early and late stages of adenosine biosynthesis, transcriptome analysis was performed with the anamorph strain of caterpillar fungus. A total of 2,764 differentially expressed genes (DEGs) were identified (p ≤ 0.05, |log₂ Ratio| ≥ 1), of which 1,737 were up-regulated and 1,027 were down-regulated. Gene expression profiling on 4–10 d revealed a distinct shift in expression of the purine metabolism pathway. Differential expression of 17 selected DEGs which involved in purine metabolism (map00230) were validated by qPCR, and the expression trends were consistent with the RNA-Seq results. Subsequently, the predicted adenosine biosynthesis pathway combined with qPCR and gene expression data of RNA-Seq indicated that the increased adenosine accumulation is a result of down-regulation of ndk, ADK, and APRT genes combined with up-regulation of AK gene. This study will be valuable for understanding the molecular mechanisms of the adenosine biosynthesis in caterpillar fungus.

Transcriptome Changes during Major Developmental Transitions Accompanied with Little Alteration of DNA Methylome in Two Pleurotus Species

Pleurotus tuoliensis (Pt) and P. eryngii var. eryngii (Pe) are important edible mushrooms. The epigenetic and gene expression signatures characterizing major developmental transitions in these two mushrooms remain largely unknown. Here, we report global analyses of DNA methylation and gene expression in both mushrooms across three major developmental transitions, from mycelium to primordium and to fruit body, by whole-genome bisulfite sequencing (WGBS) and RNA-seq-based transcriptome profiling. Our results revealed that in both Pt and Pe the landscapes of methylome are largely stable irrespective of genomic features, e.g., in both protein-coding genes and transposable elements (TEs), across the developmental transitions. The repressive impact of DNA methylation on expression of a small subset of genes is likely due to TE-associated effects rather than their own developmental dynamics. Global expression of gene orthologs was also broadly conserved between Pt and Pe, but discernible interspecific differences exist especially at the fruit body formation stage, and which are primarily due to differences in trans-acting factors. The methylome and transcriptome repertories we established for the two mushroom species may facilitate further studies of the epigenetic and transcriptional regulatory mechanisms underpinning gene during development in Pleurotus and related genera.

Diversity of cytosine methylation across the fungal tree of life

The generation of thousands of fungal genomes is leading to a better understanding of genes and genomic organization within the kingdom. However, the epigenome, which includes DNA and chromatin modifications, remains poorly investigated in fungi. Large comparative studies in animals and plants have deepened our understanding of epigenomic variation, particularly of the modified base 5-methylcytosine (5mC), but taxonomic sampling of disparate groups is needed to develop unifying explanations for 5mC variation. Here we utilize the largest phylogenetic resolution of 5mC methyltransferases (5mC MTases) and genome evolution to better understand levels and patterns of 5mC across fungi. We show that extant 5mC MTase genotypes are descendent from ancestral maintenance and de novo genotypes, whereas the 5mC MTases DIM-2 and RID are more recently derived, and that 5mC levels are correlated with 5mC MTase genotype and transposon content. Our survey also revealed that fungi lack canonical gene body methylation, which distinguishes fungal epigenomes from certain insect and plant species. However, some fungal species possess independently derived clusters of contiguous 5mC encompassing many genes. In some cases, DNA repair pathways and the N6-methyladenine (6mA) DNA modification negatively coevolved with 5mC pathways, which additionally contributed to interspecific epigenomic variation across fungi.

Study of the whole genome, methylome and transcriptome of Cordyceps militaris

The complete genome of Cordyceps militaris was sequenced using single-molecule real-time (SMRT) sequencing technology at a coverage over 300×. The genome size was 32.57 Mb, and 14 contigs ranging from 0.35 to 4.58 Mb with an N50 of 2.86 Mb were assembled, including 4 contigs with telomeric sequences on both ends and an additional 8 contigs with telomeric sequences on either the 5′ or 3′ end. A methylome database of the genome was constructed using SMRT and m4C and m6A methylated nucleotides, and many unknown modification types were identified. The major m6A methylation motif is GA and GGAG, and the major m4C methylation motif is GC or CG/GC. In the C. militaris genome DNA, there were four types of methylated nucleotides that we confirmed using high-resolution LCMS-IT-TOF. Using PacBio Iso-Seq, a total of 31,133 complete cDNA sequences were obtained in the fruiting body. The conserved domains of the nontranscribed regions of the genome include TATA boxes, which are the initial regions of genome replication. There were 406 structural variants between the HN and CM01 strains, and there were 1,114 structural variants between the HN and ATCC strains.

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

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

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.

Comparative Methylome Analysis Reveals Perturbation of Host Epigenome in Chestnut Blight Fungus by a Hypovirus

In eukaryotic genomes, DNA methylation is an important type of epigenetic modification that plays crucial roles in many biological processes. To investigate the impact of a hypovirus infection on the methylome of Cryphonectria parasitica, the chestnut blight fungus, whole-genome bisulfite sequencing (WGBS) was employed to generate single-base resolution methylomes of the fungus with/without hypovirus infection. The results showed that hypovirus infection alters methylation in all three contexts (CG, CHG, and CHH), especially in gene promoters. A total of 600 differentially methylated regions (DMRs) were identified, of which 144 could be annotated to functional genes. RNA-seq analysis revealed that DNA methylation in promoter is negatively correlated with gene expression. Among DMRs, four genes were shown to be involved in conidiation, orange pigment production, and virulence. Taken together, our DNA methylomes analysis provide valuable insights into the understanding of the relationship between DNA methylation and hypovirus infection, as well as phenotypic traits in C. parasitica.

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.

DNA methyltransferases contribute to the fungal development, stress tolerance and virulence of the entomopathogenic fungus Metarhizium robertsii

DNA methylation is an important epigenetic mark in mammals, plants, and fungi and depends on multiple genetic pathways involving de novo and maintenance DNA methyltransferases (DNMTases). Metarhizium robertsii, a model system for investigating insect-fungus interactions, has been used as an environmentally friendly alternative to chemical insecticides. However, little is known concerning the molecular basis for DNA methylation. Here, we report on the roles of two DNMTases (MrRID and MrDIM-2) by characterizing ΔMrRID, ΔMrDIM-2, and ΔRID/ΔDIM-2 mutants. The results showed that approximately 71, 10, and 8% of ^mC sites remained in the ΔMrRID, ΔMrDIM-2, and ΔRID/ΔDIM-2 strains, respectively, compared with the wild-type (WT) strain. Further analysis showed that MrRID regulates the specificity of DNA methylation and MrDIM-2 is responsible for most DNA methylation, implying an interaction or cooperation between MrRID and MrDIM-2 for DNA methylation. Moreover, the ΔMrDIM-2 and ΔRID/ΔDIM-2 strains showed more defects in radial growth and conidial production compared to the WT. Under ultraviolet (UV) irradiation or heat stress, an obvious reduction in spore viability was observed for all the mutant strains compared to the WT. The spore median lethal times (LT₅₀s) for the ΔMrDIM-2 and ΔRID/ΔDIM-2 strains in the greater wax moth, Galleria mellonella, were decreased by 47.7 and 65.9%, respectively, which showed that MrDIM-2 is required for full fungal virulence. Our data advances the understanding of the function of DNMTase in entomopathogenic fungi, which should contribute to future epigenetic investigations in fungi.

Comparative mitochondrial genomics toward exploring molecular markers in the medicinal fungus Cordyceps militaris

Cordyceps militaris is a fungus used for developing health food, but knowledge about its intraspecific differentiation is limited due to lack of efficient markers. Herein, we assembled the mitochondrial genomes of eight C. militaris strains and performed a comparative mitochondrial genomic analysis together with three previously reported mitochondrial genomes of the fungus. Sizes of the 11 mitochondrial genomes varied from 26.5 to 33.9 kb mainly due to variable intron contents (from two to eight introns per strain). Nucleotide variability varied according to different regions with non-coding regions showing higher variation frequency than coding regions. Recombination events were identified between some locus pairs but seemed not to contribute greatly to genetic variations of the fungus. Based on nucleotide diversity fluctuations across the alignment of all mitochondrial genomes, molecular markers with the potential to be used for future typing studies were determined.