Genome sequencing and analysis of fungus Hirsutella sinensis isolated from Ophiocordyceps sinensis

Combined metabolome and transcriptome analysis reveal the mechanism of water stress in Ophiocordyceps sinensis

BACKGROUND

Ophiocordyceps sinensis (O. sinensis) is the dominant bacterium in the asexual stage of Chinese cordyceps, and its growth usually suffers from water stress. Thus, simulating its ecological growth environment is crucial for artificial cultivation. This study aimed to reveal the mechanism underlying the water stress tolerance of Ophiocordyceps sinensis (O. sinensis) by combining metabolomic and transcriptome analyses to identify crucial pathways related to differentially expressed genes (DEGs) and metabolites (DEMs) involved in the response to water stress.

RESULTS

Gene coexpression analysis revealed that many genes related to 'betalain biosynthesis', 'tyrosine metabolism', 'linoleic acid metabolism', 'fructose and mannose metabolism', and 'starch and sucrose metabolism' were highly upregulated after 20d-water stress. Metabolomic analysis revealed that many metabolites regulated by these genes in these metabolic pathways were markedly decreased. On the one hand, we surmised that carbohydrate metabolism and the β-oxidation pathway worked cooperatively to generate enough acyl-CoA and then entered the TCA cycle to provide energy when exposed to water stress. On the other hand, the betalain biosynthesis and tyrosine metabolism pathway might play crucial roles in response to water stress in O. sinensis by enhancing cell osmotic potential and producing osmoregulatory substances (betaine) and antioxidant pigments (eumelanin).

CONCLUSIONS

Overall, our findings provide important information for further exploration of the mechanism underlying the water stress tolerance of O. sinensis for the industrialization of artificial cultivation of Chinese cordyceps.

Translation of Mutant Repetitive Genomic Sequences in Hirsutella sinensis and Changes in the Secondary Structures and Functional Specifications of the Encoded Proteins

Multiple repetitive sequences of authentic genes commonly exist in fungal genomes. AT-biased genotypes of Ophiocordyceps sinensis have been hypothesized as repetitive pseudogenes in the genome of Hirsutella sinensis (GC-biased Genotype #1 of O. sinensis) and are generated through repeat-induced point mutation (RIP), which is charactered by cytosine-to-thymine and guanine-to-adenine transitions, concurrent epigenetic methylation, and dysfunctionality. This multilocus study examined repetitive sequences in the H. sinensis genome and transcriptome using a bioinformatic approach and revealed that 8.2% of the authentic genes had repetitive copies, including various allelic insertions/deletions, transversions, and transitions. The transcripts for the repetitive sequences, regardless of the decreases, increases, or bidirectional changes in the AT content, were identified in the H. sinensis transcriptome, resulting in changes in the secondary protein structure and functional specification. Multiple repetitive internal transcribed spacer (ITS) copies containing multiple insertion/deletion and transversion alleles in the genome of H. sinensis were GC-biased and were theoretically not generated through RIP mutagenesis. The repetitive ITS copies were genetically and phylogenetically distinct from the AT-biased O. sinensis genotypes that possess multiple transition alleles. The sequences of Genotypes #2-17 of O. sinensis, both GC- and AT-biased, were absent from the H. sinensis genome, belong to the interindividual fungi, and differentially occur in different compartments of the natural Cordyceps sinensis insect-fungi complex, which contains >90 fungal species from >37 genera. Metatranscriptomic analyses of natural C. sinensis revealed the transcriptional silencing of 5.8S genes in all C. sinensis-colonizing fungi in natural settings, including H. sinensis and other genotypes of O. sinensis. Thus, AT-biased genotypes of O. sinensis might have evolved through advanced evolutionary mechanisms, not through RIP mutagenesis, in parallel with GC-biased Genotype #1 of H. sinensis from a common genetic ancestor over the long course of evolution.

Establishment of an Agrobacterium tumefaciens-Mediated Transformation System for Hirsutella sinensis

Ophiocordyceps sinensis (Berk.) is a complex is formed by Hepialidae larvae and Hirsutella sinensis. Infestation by H. sinensis, interaction with host larvae, and fruiting body development are three crucial processes affecting the formation of O. sinensis. However, research on the molecular mechanism of O. sinensis formation has been hindered by the lack of effective genetic transformation protocols. Therefore, Agrobacterium tumefaciens-mediated transformation (ATMT) was adopted to genetically transform two H. sinensis strains and optimize the transformation conditions. The results revealed that the most suitable Agrobacterium strain for H. sinensis transformation was AGL1, and that the surfactant Triton X-100 could also induce ATMT, although less effectively than acetosyringone (AS). In addition, the endogenous promoters of H. sinensis genes had a stronger ability to drive the expression of the target gene than did the exogenous promoter. The optimal transformation conditions were as follows: AS and hygromycin B concentrations of 100 μM and 50 μg/mL, respectively; A. tumefaciens OD600 of 0.4; cocultivation at 18 °C for 24 h; and H. sinensis used within three passages. The results lay a foundation for the functional study of key regulatory genes involved in the formation of O. 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.

The transcriptome of the entomopathogenic fungus Culicinomyces clavisporus contains an ortholog of the insecticidal ribotoxin Hirsutellin

The entomopathogenic fungus Culicinomyces clavisporus is known to infect and kill mosquito larvae and therefore has been seen as a potential biological control agent against disease vector mosquitoes. Whereas most fungal entomopathogens infect hosts by penetrating the external cuticle, C. clavisporus initiates infection through ingestion (per os). This unique infection strategy suggests that the C. clavisporus genome may be mined for novel pathogenicity factors with potential for vector control. To this end, an Isoseq-based transcriptome analysis was initiated, and resulted in a total of 3,512,145 sequences, with an average length of 1,732 bp. Transcripts assembly and annotation suggested that the C. clavisporus transcriptome lacked the cuticle-degrading proteins that have been associated with other entomopathogenic fungi, supporting the per os pathogenicity process. Furthermore, mining of the sequence data unexpectedly revealed C. clavisporus transcripts homologous to the Hirsutellin toxin. Comparative sequence analyses indicated that the C. clavisporus Hirsutellin predicted protein has retained the canonical molecular features that have been associated with the ribotoxic and insecticidal properties of the original toxin isolated from Hirsutella thompsonii. The identification of an Hirsutellin ortholog in C. clavisporus was supported by phylogenetic analyses demonstrating that Culicinomyces and Hirsutella were closely related genera in the Ophiocordycipitaceae family. Validation of the mosquitocidal activity of this novel C. clavisporus protein has yet to be performed but may help position Hirsutellin orthologs as prime candidates for the development of alternative biocontrol approaches complementing the current toolbox of vector mosquito management strategies.

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

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

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.

Unraveling the Secrets of a Double-Life Fungus by Genomics: Ophiocordyceps australis CCMB661 Displays Molecular Machinery for Both Parasitic and Endophytic Lifestyles

Ophiocordyceps australis (Ascomycota, Hypocreales, Ophiocordycipitaceae) is a classic entomopathogenic fungus that parasitizes ants (Hymenoptera, Ponerinae, Ponerini). Nonetheless, according to our results, this fungal species also exhibits a complete set of genes coding for plant cell wall degrading Carbohydrate-Active enZymes (CAZymes), enabling a full endophytic stage and, consequently, its dual ability to both parasitize insects and live inside plant tissue. The main objective of our study was the sequencing and full characterization of the genome of the fungal strain of O. australis (CCMB661) and its predicted secretome. The assembled genome had a total length of 30.31 Mb, N50 of 92.624 bp, GC content of 46.36%, and 8,043 protein-coding genes, 175 of which encoded CAZymes. In addition, the primary genes encoding proteins and critical enzymes during the infection process and those responsible for the host-pathogen interaction have been identified, including proteases (Pr1, Pr4), aminopeptidases, chitinases (Cht2), adhesins, lectins, lipases, and behavioral manipulators, such as enterotoxins, Protein Tyrosine Phosphatases (PTPs), and Glycoside Hydrolases (GHs). Our findings indicate that the presence of genes coding for Mad2 and GHs in O. australis may facilitate the infection process in plants, suggesting interkingdom colonization. Furthermore, our study elucidated the pathogenicity mechanisms for this Ophiocordyceps species, which still is scarcely studied.

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.

The complete mitochondrial genome of Ahamus yushuensis Chu et Wang 1985 (Lepidoptera: Hepialidae) and phylogenetic analysis

The complete mitochondrial genome (mitogenome) of Ahamus yushuensis was determined in this study. This mitogenome is 15,336 bp and encodes 37 mitochondrial genes, including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), and two ribosomal RNA genes (rrnL and rrnS). The A. yushuensis mitogenome has an A + T content of 82.2% and presents a positive AT-skew (0.052) and a negative GC-skew (−0.236). Twelve PCGs start with a typical ATN codon, whereas a single PCG uses CGA (coxI) as the initial codon. The maximum likelihood phylogenetic analysis based on the concatenated nucleotide sequences of 13 PCGs strongly supported the monophyletic relationship of A. yushuensis to the clade of Thitarodes damxungensis and A. yunnanensis.

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.