A New High-Quality Draft Genome Assembly of the Chinese Cordyceps Ophiocordyceps sinensis

Mutations and Differential Transcription of Mating-Type and Pheromone Receptor Genes in Hirsutella sinensis and the Natural Cordyceps sinensis Insect-Fungi Complex

Sexual reproduction in ascomycetes is controlled by the mating-type (MAT) locus. (Pseudo)homothallic reproduction has been hypothesized on the basis of genetic data from Hirsutella sinensis (Genotype #1 of Ophiocordyceps sinensis). However, the differential occurrence and differential transcription of mating-type genes in the MAT1-1 and MAT1-2 idiomorphs were found in the genome and transcriptome assemblies of H. sinensis, and the introns of the MAT1-2-1 transcript were alternatively spliced with an unspliced intron I that contains stop codons. These findings reveal that O. sinensis reproduction is controlled at the genetic, transcriptional, and coupled transcriptional-translational levels. This study revealed that mutant mating proteins could potentially have various secondary structures. Differential occurrence and transcription of the a-/α-pheromone receptor genes were also found in H. sinensis. The data were inconsistent with self-fertilization under (pseudo)homothallism but suggest the self-sterility of H. sinensis and the requirement of mating partners to achieve O. sinensis sexual outcrossing under heterothallism or hybridization. Although consistent occurrence and transcription of the mating-type genes of both the MAT1-1 and MAT1-2 idiomorphs have been reported in natural and cultivated Cordyceps sinensis insect-fungi complexes, the mutant MAT1-1-1 and α-pheromone receptor transcripts in natural C. sinensis result in N-terminal or middle-truncated proteins with significantly altered overall hydrophobicity and secondary structures of the proteins, suggesting heterogeneous fungal source(s) of the proteins and hybridization reproduction because of the co-occurrence of multiple genomically independent genotypes of O. sinensis and >90 fungal species in natural C. sinensis.

Selection and validation of reference genes for RT-qPCR in ophiocordyceps sinensis under different experimental conditions

The Chinese caterpillar mushroom, Ophiocordyceps sinensis (O. sinensis), is a rarely medicinal fungus in traditional chinese herbal medicine due to its unique medicinal values, and the expression stability of reference genes is essential to normalize its gene expression analysis. In this study, BestKeeper, NormFinder and geNorm, three authoritative statistical arithmetics, were applied to evaluate the expression stability of sixteen candidate reference genes (CRGs) in O. sinensis under different stress [low temperature (4°C), light treatment (300 lx), NaCl (3.8%)] and different development stages (mycelia, primordia and fruit bodies) and formation of morphologic mycelium (aeriasubstrate, hyphae knot mycelium). The paired variation values indicated that two genes could be enough to accurate standardization exposed to different conditions of O.sinensis. Among these sixteen CRGs, 18S ribosomal RNA (18S rRNA) and beta-Tubulin (β-TUB) showed the topmost expression stability in O.sinensis exposed to all conditions, while glutathione hydrolase proenzym (GGT) and Phosphoglucose isomerase (PGI) showed the least expression stability. The optimal reference gene in different conditions was various. β-TUB and Ubiquitin (UBQ) were identified as the two most stable genes in different primordia developmental stage, while phosphoglucomutase (PGM) with elongation factor 1-alpha (EF1-α) and 18S rRNA with UBQ were the most stably expressed for differentially morphologic mycelium stages and different stresses, respectively. These results will contribute to more accurate evaluation of the gene relative expression levels in O.sinensis under different conditions using the optimal reference gene in real-time quantitative PCR (RT-qPCR) analysis.

Characterization of the Biosynthetic Gene Cluster for the Ribosomally Synthesized Cyclic Peptide Epichloëcyclins in Epichloë festucae

The various grass-induced epichloëcyclins of the Epichloë spp. are ribosomally synthesized and post-translationally modified peptides (RiPPs), produced as small, secreted cyclopeptides from a single gene, gigA. Here, four clustered and coregulated genes (gigA, gigB, gigC, and kexB) with predicted roles in epichloëcyclin production in Epichloë festucae were evaluated through gene disruption. Subsequent chemical analysis indicates that GigB is a DUF3328 domain-containing protein associated with cyclization of epichloëcyclins; GigC is a methyltransferase enzyme responsible for N-methylation of desmethylepichloëcyclins; and KexB is a subtilisin-like enzyme, partly responsible for the propeptide cleavage of epichloëcyclin intermediates. Symbiotic effects on the host phenotype were not observed for gigA, gigC, or kexB mutants, although ΔgigB infection correlated with increased host tiller height and biomass, while only ΔkexB exhibited an effect on endophyte morphology. Disrupting epichloëcyclin biosynthesis showed negligible influence on the biosynthesis of E. festucae-associated alkaloids. Epichloëcyclins may perform other secondary metabolism functions in Epichloë and other fungi.

Transcriptome and metabonomics combined analysis revealed the energy supply mechanism involved in fruiting body initiation in Chinese cordyceps

Chinese cordyceps was one of most valuable traditional Chinese medicine fungi. To elucidate the molecular mechanisms related to energy supply mechanism involved in the initiation and formation of primordium in Chinese cordyceps, we performed the integrated metabolomic and transcriptomic analyses of it at pre-primordium period, primordium germination period and after-primordium period, respectively. Transcriptome analysis showed that many genes related to 'starch and sucrose metabolism', 'fructose and mannose metabolism', 'linoleic acid metabolism', 'fatty acids degradation' and 'glycerophospholipid metabolism' were highly up-regulated at primordium germination period. Metabolomic analysis showed many metabolites regulated by these genes in these metabolism pathways were also markedly accumulated at this period. Consequently, we inferred that carbohydrate metabolism and β-oxidation pathway of palmitic acid and linoleic acid worked cooperatively to generate enough acyl-CoA, and then entered TCA cycle to provide energy for fruiting body initiation. Overall, our finding provided important information for further exploring the energy metabolic mechanisms of realizing the industrialization of Chinese cordyceps artificial cultivation.

Altered GC- and AT-biased genotypes of Ophiocordyceps sinensis in the stromal fertile portions and ascospores of natural Cordyceps sinensis

OBJECTIVE

To examine multiple genotypes of Ophiocordyceps sinensis in a semi-quantitative manner in the stromal fertile portion (SFP) densely covered with numerous ascocarps and ascospores of natural Cordyceps sinensis and to outline the dynamic alterations of the coexisting O. sinensis genotypes in different developmental phases.

METHODS

Mature Cordyceps sinensis specimens were harvested and continuously cultivated in our laboratory (altitude 2,254 m). The SFPs (with ascocarps) and fully and semi-ejected ascospores were collected for histological and molecular examinations. Biochip-based single nucleotide polymorphism (SNP) MALDI-TOF mass spectrometry (MS) was used to genotype multiple O. sinensis mutants in the SFPs and ascospores.

RESULTS

Microscopic analysis revealed distinct morphologies of the SFPs (with ascocarps) before and after ascospore ejection and SFP of developmental failure, which, along with the fully and semi-ejected ascospores, were subjected to SNP MS genotyping analysis. Mass spectra showed the coexistence of GC- and AT-biased genotypes of O. sinensis that were genetically and phylogenetically distinct in the SFPs before and after ejection and of developmental failure and in fully and semi-ejected ascospores. The intensity ratios of MS peaks were dynamically altered in the SFPs and the fully and semi-ejected ascospores. Mass spectra also showed transversion mutation alleles of unknown upstream and downstream sequences with altered intensities in the SFPs and ascospores. Genotype #5 of AT-biased Cluster-A maintained a high intensity in all SFPs and ascospores. An MS peak with a high intensity containing AT-biased Genotypes #6 and #15 in pre-ejection SFPs was significantly attenuated after ascospore ejection. The abundance of Genotypes #5‒6 and #16 of AT-biased Cluster-A was differentially altered in the fully and semi-ejected ascospores that were collected from the same Cordyceps sinensis specimens.

CONCLUSION

Multiple O. sinensis genotypes coexisted in different combinations with altered abundances in the SFPs prior to and after ejection, the SFP of developmental failure, and the two types of ascospores of Cordyceps sinensis, demonstrating their genomic independence. Metagenomic fungal members present in different combinations and with dynamic alterations play symbiotic roles in different compartments of natural Cordyceps sinensis.

Differential coexistence of multiple genotypes of Ophiocordyceps sinensis in the stromata, ascocarps and ascospores of natural Cordyceps sinensis

OBJECTIVE

To examine the differential occurrence of Ophiocordyceps sinensis genotypes in the stroma, stromal fertile portion (SFP) densely covered with numerous ascocarps, and ascospores of natural Cordyceps sinensis.

METHODS

Immature and mature C. sinensis specimens were harvested. Mature C. sinensis specimens were continuously cultivated in our laboratory (altitude 2,200 m). The SFPs (with ascocarps) and ascospores of C. sinensis were collected for microscopic and molecular analyses using species-/genotype-specific primers. Sequences of mutant genotypes of O. sinensis were aligned with that of Genotype #1 Hirsutella sinensis and compared phylogenetically using a Bayesian majority-rule method.

RESULTS

Fully and semiejected ascospores were collected from the same specimens. The semiejected ascospores tightly adhered to the surface of the asci as observed by the naked eye and under optical and confocal microscopies. The multicellular heterokaryotic ascospores showed uneven staining of nuclei. The immature and mature stromata, SFPs (with ascocarps) and ascospores were found to differentially contain several GC- and AT-biased genotypes of O. sinensis, Samsoniella hepiali, and an AB067719-type fungus. The genotypes within AT-biased Cluster-A in the Bayesian tree occurred in all compartments of C. sinensis, but those within AT-biased Cluster-B were present in immature and mature stromata and SPFs but absent in the ascospores. Genotype #13 of O. sinensis was present in semi-ejected ascospores and Genotype #14 in fully ejected ascospores. GC-biased Genotypes #13-14 featured large DNA segment substitutions and genetic material recombination between the genomes of the parental fungi (H. sinensis and the AB067719-type fungus). These ascosporic offspring genotypes combined with varying abundances of S. hepiali in the 2 types of ascospores participated in the control of the development, maturation and ejection of the ascospores.

CONCLUSION

Multiple genotypes of O. sinensis coexist differentially in the stromata, SFPs and 2 types of C. sinensis ascospores, along with S. hepiali and the AB067719-type fungus. The fungal components in different combinations and their dynamic alterations in the compartments of C. sinensis during maturation play symbiotic roles in the lifecycle of natural C. sinensis.

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.

Integrated metabolomics and transcriptomics reveal metabolites difference between wild and cultivated Ophiocordyceps sinensis

Ophiocordyceps sinensis is a traditional medicinal fungus endemic to the alpine and high-altitude areas of the Qinghai-Tibet plateau. The scarcity of the wild resource has led to increased attention to artificially cultivated O. sinensis. However, little is known about the metabolic differences and the regulatory mechanisms between cultivated and wild O. sinensis. This study exploited untargeted metabolomics and transcriptomics to uncover the differences in accumulated metabolites and expressed genes between wild and cultivated O. sinensis. Metabolomics results revealed that 368 differentially accumulated metabolites were mainly enriched in biosynthesis of amino acids, biosynthesis of plant secondary metabolites and purine nucleotide metabolism. Cultivated O. sinensis contained more amino acids and derivatives, carbohydrates and derivatives, and phenolic acids than wild O. sinensis, whereas the contents of most nucleosides and nucleotides in wild O. sinensis were significantly higher than in cultivated O. sinensis. Transcriptome analysis indicated that 4430 annotated differentially expressed genes were identified between two types. Integrated metabolomics and transcriptomics analyses suggested that IMPDH, AK, ADSS, guaA and GUK genes might be related to the synthesis of purine nucleotides and nucleosides. Our findings will provide a new insight into the molecular basis of metabolic variations of this medicinal fungus.

The genome of a hadal sea cucumber reveals novel adaptive strategies to deep-sea environments

How organisms cope with coldness and high pressure in the hadal zone remains poorly understood. Here, we sequenced and assembled the genome of hadal sea cucumber Paelopatides sp. Yap with high quality and explored its potential mechanisms for deep-sea adaptation. First, the expansion of ACOX1 for rate-limiting enzyme in the DHA synthesis pathway, increased DHA content in the phospholipid bilayer, and positive selection of EPT1 may maintain cell membrane fluidity. Second, three genes for translation initiation factors and two for ribosomal proteins underwent expansion, and three ribosomal protein genes were positively selected, which may ameliorate the protein synthesis inhibition or ribosome dissociation in the hadal zone. Third, expansion and positive selection of genes associated with stalled replication fork recovery and DNA repair suggest improvements in DNA protection. This is the first genome sequence of a hadal invertebrate. Our results provide insights into the genetic adaptations used by invertebrate in deep oceans.

Medium optimization for high mycelial soluble protein content of Ophiocordyceps sinensis using response surface methodology

Ophiocordyceps sinensis is widely utilized due to its pharmaceutical value. Mycelial protein forms a key active component of O. sinensis and determines the medicinal potential of fungus. Here, we describe the development of an optimized fermentation medium to obtain more mycelial soluble protein from O. sinensis using response surface methodology (RSM) and investigate the increased mycelial protein content using transcriptomics. The maximum mycelial protein content of 2.11% was obtained using a medium consisting of 20% beef broth, 0.10% peptone, 2% glucose, 0.15% yeast extract, 0.20% KH₂PO₄, and 0.02% MgSO₄. Transcriptome analysis identified 790 differentially expressed genes (DEGs), including 592 up-regulated genes and 198 down-regulated genes, optimisation resulted in more up-regulated genes. The main DEGs were enriched in metabolic pathways, ABC transporters, starch and sucrose metabolism, tyrosine metabolism, and glutathione metabolism. In addition, some DEGs associated with mycelial protein enhancement such as tyrosinase (TYR), glutathione S-transferase (GST), glutamine synthetase (glnA), and β-glucosidase may contribute to increased mycelial protein content. Real-time quantitative PCR (RT-qPCR) was used to confirm gene expression and the results support the accuracy of RNA-Seq and DEG analysis. This study provides an optimized fermentation method for enhancing the mycelial protein content of O. sinensis and a reference for the effective development of O. sinensis protein.

The invasion process of the entomopathogenic fungus Ophiocordyceps sinensis into the larvae of ghost moths (Thitarodes xiaojinensis) using a GFP-labeled strain

Chinese cordyceps is a well-known and valuable traditional Chinese medicine that forms after Ophiocordyceps sinensis parasitizes ghost moth larvae. The low natural infection rate of O. sinensis limits large-scale artificial cultivation of Chinese cordyceps, and the invasion process is unclear. To investigate the temporal and spatial regulation when O. sinensis enters ghost moths, we constructed an O. sinensis transformant that stably expresses green fluorescent protein (GFP). Inoculating Thitarodes xiaojinensis larvae with a high concentration of GFP-labeled O. sinensis, we observed that O. sinensis conidia could adhere to the host cuticle within 2 days, germinate penetration pegs within 4 days, and produce blastospores in the host hemocoel within 6 days. The reconstructed three-dimensional (3D) structures of the invasion sites showed that penetration pegs germinated directly from O. sinensis conidia at the joining site with the larval cuticle. Differentiated appressoria or hyphae along the host epicuticle are not required for O. sinensis to invade ghost moths. Overall, the specific invasion process of O. sinensis into its host is clarified, and we provided a new perspective on the invasion process of entomopathogenic fungi.

Genomic Sequence and Transcriptome Analysis of the Medicinal Fungus Keithomyces neogunnii

The filamentous fungus Keithomyces neogunnii can infect the larvae of Lepidoptera (Hepialus sp.) and form an insect-fungi complex, which is utilized as an important traditional Chinese medicine. As a valuable medicinal fungus, K. neogunnii produces diverse bioactive substances (e.g., polysaccharide, vitamins, cordycepic acid, and adenosine) under cultivation conditions. Herein, we report the first high-quality genome of the K. neogunnii single-spore isolate Cg7.2a using single-molecule real-time sequencing technology in combination with Illumina sequencing. The assembled genome was 32.6 Mb in size, containing 8,641 predicted genes and having a GC content of 52.16%. RNA sequencing analysis revealed the maximum number of differentially expressed genes in the fungus during the stroma formation stage compared with those during the mycelium stage. These data are valuable to enhance our understanding of the biology, development, evolution, and physiological metabolism of K. neogunnii.

Lipid-Binding Aegerolysin from Biocontrol Fungus Beauveria bassiana

Fungi are the most common pathogens of insects and thus important regulators of their populations. Lipid-binding aegerolysin proteins, which are commonly found in the fungal kingdom, may be involved in several biologically relevant processes including attack and defense against other organisms. Aegerolysins act alone or together with membrane-attack-complex/perforin (MACPF)-like proteins to form transmembrane pores that lead to cell lysis. We performed an in-depth bioinformatics analysis of aegerolysins in entomopathogenic fungi and selected a candidate aegerolysin, beauveriolysin A (BlyA) from Beauveria bassiana. BlyA was expressed as a recombinant protein in Escherichia coli, and purified to further determine its functional and structural properties, including lipid-binding ability. Aegerolysins were found to be encoded in genomes of entomopathogenic fungi, such as Beauveria, Cordyceps, Metarhizium and Ophiocordyceps. Detailed bioinformatics analysis revealed that they are linked to MACPF-like genes in most genomes. We also show that BlyA interacts with an insect-specific membrane lipid. These results were placed in the context of other fungal and bacterial aegerolysins and their partner proteins. We believe that aegerolysins play a role in promoting the entomopathogenic and antagonistic activity of B. bassiana, which is an active ingredient of bioinsecticides.

An Overview of Genomics, Phylogenomics and Proteomics Approaches in Ascomycota

Fungi are among the most successful eukaryotes on Earth: they have evolved strategies to survive in the most diverse environments and stressful conditions and have been selected and exploited for multiple aims by humans. The characteristic features intrinsic of Fungi have required evolutionary changes and adaptations at deep molecular levels. Omics approaches, nowadays including genomics, metagenomics, phylogenomics, transcriptomics, metabolomics, and proteomics have enormously advanced the way to understand fungal diversity at diverse taxonomic levels, under changeable conditions and in still under-investigated environments. These approaches can be applied both on environmental communities and on individual organisms, either in nature or in axenic culture and have led the traditional morphology-based fungal systematic to increasingly implement molecular-based approaches. The advent of next-generation sequencing technologies was key to boost advances in fungal genomics and proteomics research. Much effort has also been directed towards the development of methodologies for optimal genomic DNA and protein extraction and separation. To date, the amount of proteomics investigations in Ascomycetes exceeds those carried out in any other fungal group. This is primarily due to the preponderance of their involvement in plant and animal diseases and multiple industrial applications, and therefore the need to understand the biological basis of the infectious process to develop mechanisms for biologic control, as well as to detect key proteins with roles in stress survival. Here we chose to present an overview as much comprehensive as possible of the major advances, mainly of the past decade, in the fields of genomics (including phylogenomics) and proteomics of Ascomycota, focusing particularly on those reporting on opportunistic pathogenic, extremophilic, polyextremotolerant and lichenized fungi. We also present a review of the mostly used genome sequencing technologies and methods for DNA sequence and protein analyses applied so far for fungi.