Genetic diversity and distribution patterns of host insects of Caterpillar Fungus Ophiocordyceps sinensis in the Qinghai-Tibet Plateau

An authentic assessment method for cordyceps sinensis

Cordyceps Sinensis, renowned for its diverse pharmacological properties and the rarity of its natural species, faces significant challenges due to rampant adulteration by counterfeit products. Consequently, there is a crucial need to reliably identify Cordyceps species to ensure their quality and efficacy. While current analytical techniques predominantly rely on LC-MS, there remains a notable deficiency and substantial demand for the development of a unified, reproducible, and fast method suitable for commercial applications. In this study, we employed a cost-effective and straightforward approach utilizing headspace GC-MS to authenticate Cordyceps sinensis. This method enables the comprehensive analysis of the chemical profile, facilitating the identification of quality and authenticity in Cordyceps samples. Through a comparative analysis of the chemical profiles of seven authentic Cordyceps samples with seven other Cordyceps samples, we propose a Quality Assessment System for Authentic Cordyceps, encompassing the following criteria: 1) the presence of 29 compounds commonly found in authentic Cordyceps within the chemical profile, and 2) the area ratio of 3-methylbutanal to 2-methylbutanal falling within the range of 2.09-3.01. This method exhibits considerable promise as a standardized, reproducible, and expeditious technique for the quality assessment and authentication of Cordyceps.

Population structure and adaptability analysis of Schizothorax o'connori based on whole-genome resequencing

BACKGROUND

Schizothorax o'connori is an endemic fish distributed in the upper and lower reaches of the Yarlung Zangbo River in China. It has experienced a fourth round of whole gene replication events and is a good model for exploring the genetic differentiation and environmental adaptability of fish in the Qinghai-Tibet Plateau. The uplift of the Qinghai-Tibet Plateau has led to changes in the river system, thereby affecting gene exchange and population differentiation between fish populations. With the release of fish whole genome data, whole genome resequencing has been widely used in genetic evolutionary analysis and screening of selected genes in fish, which can better elucidate the genetic basis and molecular environmental adaptation mechanisms of fish. Therefore, our purpose of this study was to understand the population structure and adaptive characteristics of S. o'connori using the whole-genome resequencing method.

RESULTS

The results showed that 23,602,746 SNPs were identified from seven populations, mostly distributed on chromosomes 2 and 23. There was no significant genetic differentiation between the populations, and the genetic diversity was relatively low. However, the Zangga population could be separated from the Bomi, Linzhi, and Milin populations in the cluster analysis. Based on historical dynamics analysis of the population, the size of the ancestral population of S. o'connori was affected by the late accelerated uplift of the Qinghai Tibet Plateau and the Fourth Glacial Age. The selected sites were mostly enriched in pathways related to DNA repair and energy metabolism.

CONCLUSION

Overall, the whole-genome resequencing analysis provides valuable insights into the population structure and adaptive characteristics of S. o'connori. There was no obvious genetic differentiation at the genome level between the S. o'connori populations upstream and downstream of the Yarlung Zangbo River. The current distribution pattern and genetic diversity are influenced by the late accelerated uplift of the Qinghai Tibet Plateau and the Fourth Ice Age. The selected sites of S. o'connori are enriched in the energy metabolism and DNA repair pathways to adapt to the low temperature and strong ultraviolet radiation environment at high altitude.

Conservation Priorities and Demographic History of Saussurea involucrata in the Tianshan Mountains and Altai Mountains

Rare and vulnerable endemic plants represent different evolutionary units that occur at different times, and protecting these species is a key issue in biological protection. Understanding the impact of the history of endangered plant populations on their genetic diversity helps to reveal evolutionary history and is crucial for guiding conservation efforts. Saussurea involucrata, a perennial alpine species mainly distributed in the Tianshan Mountains, is famous for its medicinal value but has become endangered due to over-exploitation. In the present study, we employed both nuclear and chloroplast DNA sequences to investigate the genetic distribution pattern and evolutionary history of S. involucrata. A total of 270 individuals covering nine S. involucrata populations were sampled for the amplification and sequencing of nrDNA Internal Transcribed Spacer (ITS) and chloroplast trnL-trnF, matK and ndhF-rpl32 sequences. Via calculation, we identified 7 nuclear and 12 plastid haplotypes. Among the nine populations, GL and BA were characterized by high haplotype diversity, whereas BG revealed the lowest haplotype diversity. Molecular dating estimations suggest that divergence among S. involucrata populations occurred around 0.75 Ma, coinciding with the uplift of Tianshan Mountains. Our results reveal that both isolation-by-distance (IBD) and isolation-by-resistance (IBR) have promoted genetic differentiation among populations of S. involucrata. The results from the ecological niche modeling analyses show a more suitable habitat for S. involucrata in the past than at present, indicating a historical distribution contraction of the species. This study provides new insight into understanding the genetic differentiation of S. involucrata, as well as the theoretical basis for conserving this species.

The phylogeny and divergence time of Ophiocordyceps sinensis and its host insects based on elongation factor 1 alpha

Ophiocordyceps sinensis Berk. is a fungal parasite that parasitizes the larvae of Hepialidae and is endemic to the Qinghai-Tibet Plateau (QTP). The phylogeny and divergence time of O. sinensis and its host insects were analyzed for 137 individuals from 48 O. sinensis populations based on the elongation factor 1 alpha (EF-1α) gene. Lower nucleotide variation, with only 7 and 16 EF-1α haplotypes, was detected in O. sinensis and its host insects, respectively. The isolated and broad distribution patterns coexisted in both O. sinensis and its host insects on the QTP. The divergence time estimates show that O. sinensis and its host insects originated later than 14.33 million years (Myr) and earlier than 23.60 Myr in the Miocene period, and the major differentiation occurred later than 4 Myr. Their origin and differentiation match well with the second and third uplifts of the QTP, respectively. The host insects from the O. sinensis populations distributed around Qinghai Lake are inferred as an ancient and relict species that has survived various geological events of the QTP. It is suitable to estimate the divergence times of both O. sinensis and its host insects from the same individuals using one gene: EF-1α. Our findings of the origin, phylogeny, and evolution of the endemic species also support the epoch of geological events on the QTP.

Out of the Qinghai-Tibetan plateau: Origin, evolution and historical biogeography of Morchella (both Elata and Esculenta clades)

Introduction

Morchella has become a research hotspot because of its wide distribution, delicious taste, and phenotypic plasticity. The Qinghai-Tibet Plateau subkingdoms (QTPs) are known as the cradle of Ice age biodiversity. However, the diversity of Morchella in the QTPs has been poorly investigated, especially in phylogenetic diversity, origin, and biogeography.

Methods

The genealogical concordance phylogenetic species recognition (GCPSR, based on Bayesian evolutionary analysis using sequences from the internal transcribed spacer (ITS), nuclear large subunit rDNA (nrLSU), translation elongation factor 1-α (EF1-α), and the largest and second largest subunits of RNA polymerase II (RPB1 and RPB2)), differentiation time estimation, and ancestral region reconstruction were used to infer Morchella's phylogenetic relationships and historical biogeography in the QTPs.

Results

Firstly, a total of 18 Morchella phylogenetic species are recognized in the QTPs, including 10 Elata clades and 8 Esculenta clades of 216 individuals Secondly, the divergences of the 18 phylogenetic species were 50.24-4.20 Mya (Eocene-Pliocene), which was closely related to the geological activities in the QTPs. Furthermore, the ancestor of Morchella probably originated in the Northern regions (Qilian Shan, Elata cade) and southwestern regions (Shangri-La, Esculenta clade) of QTPs and might have migrated from North America (Rufobrunnea clade) via Beringian Land Bridge (BLB) and Long-Distance Dispersal (LDD) expansions during the Late Cretaceous. Moreover, as the cradle of species origin and diversity, the fungi species in the QTPs have spread out and diffused to Eurasia and South Africa starting in the Paleogene Period.

Conclusion

This is the first report that Esculenta and Elata clade of Morchella originated from the QTPs because of orogenic, and rapid differentiation of fungi is strongly linked to geological uplift movement and refuge in marginal areas of the QTPs. Our findings contribute to increasing the diversity of Morchella and offer more evidence for the origin theory of the QTPs.

Interspecific Hybridization and Complete Mitochondrial Genome Analysis of Two Ghost Moth Species

Simple Summary

The Chinese cordyceps is a valuable parasitic Ophiocordyceps sinensis fungus–Thitarodes/Hepialus larva complex. In view of culturing this complex, a method for the artificial rearing of the Thitarodes/Hepialus ghost moth hosts was established. Deterioration of the host insect population and low mummification rates in infected larvae constrain effective cultivation. Hybridization of Thitarodes/Hepialus populations may overcome this problem. Thitarodes shambalaensis and Thitarodes sp. were inbred or hybridized, and the biological parameters, larval sensitivity to the fungal infection and mitochondrial genomes of the resulting populations were investigated. Hybridization of T. shambalaensis and Thitarodes sp. allowed producing a new generation. One hybrid population (T. shambalaensis females mated with Thitarodes sp. males) showed increased population growth as compared with the parental Thitarodes sp. population. The sensitivity of the inbred larval populations to four fungal isolates of O. sinensis differed. The complete mitochondrial genomes of T. shambalaensis, Thitarodes sp. and the hybrid population were 15,612 bp, 15,389 bp and 15,496 bp in length, respectively. A + T-rich regions were variable in sizes and repetitive sequences. The hybrid population was located in the same clade with T. shambalaensis, implying the maternal inheritance of mitochondrial DNA.

Abstract

The Chinese cordyceps, a parasitic Ophiocordyceps sinensis fungus–Thitarodes/Hepialus larva complex, is a valuable biological resource endemic to the Tibetan Plateau. Protection of the Plateau environment and huge market demand make it necessary to culture this complex in an artificial system. A method for the large-scale artificial rearing of the Thitarodes/Hepialus insect host has been established. However, the deterioration of the insect rearing population and low mummification of the infected larvae by the fungus constrain effective commercial cultivation. Hybridization of Thitarodes/Hepialus populations may be needed to overcome this problem. The species T. shambalaensis (GG♂ × GG♀) and an undescribed Thitarodes species (SD♂ × SD♀) were inbred or hybridized to evaluate the biological parameters, larval sensitivity to the fungal infection and mitochondrial genomes of the resulting populations. The two parental Thitarodes species exhibited significant differences in adult fresh weights and body lengths but not in pupal emergence rates. Hybridization of T. shambalaensis and Thitarodes sp. allowed producing a new generation. The SD♂ × GG♀ population showed a higher population trend index than the SD♂ × SD♀ population, implying increased population growth compared with the male parent. The sensitivity of the inbred larval populations to four fungal isolates of O. sinensis also differed. This provides possibilities to create Thitarodes/Hepialus populations with increased growth potential for the improved artificial production of the insect hosts. The mitochondrial genomes of GG♂ × GG♀, SD♂ × SD♀ and SD♂ × GG♀ were 15,612 bp, 15,389 bp and 15,496 bp in length, with an A + T content of 80.92%, 82.35% and 80.87%, respectively. The A + T-rich region contains 787 bp with two 114 bp repetitive sequences, 554 bp without repetitive sequences and 673 bp without repetitive sequences in GG♂ × GG♀, SD♂ × SD♀ and SD♂ × GG♀, respectively. The hybrid population (SD♂ × GG♀) was located in the same clade with GG♂ × GG♀, based on the phylogenetic tree constructed by 13 PCGs, implying the maternal inheritance of mitochondrial DNA.

Stage- and Rearing-Dependent Metabolomics Profiling of Ophiocordyceps sinensis and Its Pipeline Products

Cordyceps, a parasitic complex of the fungus Ophiocordyceps sinensis (Berk.) (Hypocreales: Ophiocordycipitaceae) and the ghost moth Thitarodes (Lepidoptera: Hepialidae), is a historical ethnopharmacological commodity in China. Recently, artificial cultivation of Chinese cordyceps has been established to supplement the dwindling natural resources. However, much is unknown between the natural and cultivated products in terms of nutritional aspect, which may provide essential information for quality evaluation. The current study aims to determine the metabolic profiles of 17 treatments from 3 sample groups including O. sinensis fungus, Thitarodes insect and cordyceps complex, using Gas Chromatography - Quadrupole Time-of-Flight Mass Spectrometry. A total of 98 metabolites were detected, with 90 of them varying in concentrations among groups. The tested groups could be separated, except that fungal fruiting body was clustered into the same group as Chinese cordyceps. The main distinguishing factors for the groups studied were the 24 metabolites involved in numerous different metabolic pathways. In conclusion, metabolomics of O. sinensis and its related products were determined mainly by the fruiting bodies other than culture methods. Our results suggest that artificially cultured fruiting bodies and cordyceps may share indistinguishable metabolic functions as the natural ones.

Novel formulation development from Ophiocordyceps sinensis (Berk.) for management of high-altitude maladies

Ophiocordyceps sinensis (Berk.) is a fungus closely related to medicinal mushroom, which belongs to the family Ophiocordycipitaceae. It is a well-known and rich herbal source of bioactive active constituents. The medicinal mushroom has garnered worldwide attention owing to its multifarious bioactivities. This mushroom grows on the larva of ghost moths (Hepialidae) and produces fruiting bodies, which serve as a vital natural source of medicine and supplementary diets. On account of the diverse pharmacological and bioactive constituents present in O. sinensis, it has been established as a potential antioxidant, anticancer, antibacterial, anti-proliferative, anti-inflammatory agent that has been successfully used for treating several health issues, including hypoxia-related problems encountered by mountaineers, pilgrims, tourists and soldiers occurring at high-altitude regions such as acute mountain sickness (AMS), high-altitude pulmonary edema (HAPE), high-altitude cerebral edema (HACE), frostbite, chilblains, hypothermia, etc. The most important pharmacologically active compounds present in the O. sinensis include nucleobases and its derivatives (adenosine, cordycepin, 3-deoxyadenosine, AMP, GMP, UMP, guanosine, uridine), polysaccharides (mannose, glucose, galactose, rhamnose, arabinose, xylose, galactose), proteins, peptides and steroids. This article focuses on the various research endeavors undertaken to scientifically establish the medicinal properties of O. sinensis, highlighting the various principally active compounds, their pharmacological action, drug designing and development and future perspective for various health benefits.

Uncovering fungal community composition in natural habitat of Ophiocordyceps sinensis using high-throughput sequencing and culture-dependent approaches

Background

The fungal communities inhabiting natural Ophiocordyceps sinensis play critical ecological roles in alpine meadow ecosystem, contribute to infect host insect, influence the occurrence of O. sinensis, and are repertoire of potential novel metabolites discovery. However, a comprehensive understanding of fungal communities of O. sinensis remain elusive. Therefore, the present study aimed to unravel fungal communities of natural O. sinensis using combination of high-throughput sequencing and culture-dependent approaches.

Results

A total of 280,519 high-quality sequences, belonging to 5 fungal phyla, 15 classes, 41 orders, 79 families, 112 genera, and 352 putative operational taxonomic units (OTUs) were obtained from natural O. sinensis using high-throughput sequencing. Among of which, 43 genera were identified in external mycelial cortices, Ophiocordyceps, Sebacinia and Archaeorhizomyces were predominant genera with the abundance of 95.86, 1.14, 0.85%, respectively. A total of 66 genera were identified from soil microhabitat, Inocybe, Archaeorhizomyces, unclassified Thelephoraceae, Tomentella, Thelephora, Sebacina, unclassified Ascomycota and unclassified fungi were predominant genera with an average abundance of 53.32, 8.69, 8.12, 8.12, 7.21, 4.6, 3.08 and 3.05%, respectively. The fungal communities in external mycelial cortices were significantly distinct from soil microhabitat. Meanwhile, seven types of culture media were used to isolate culturable fungi at 16 °C, resulted in 77 fungal strains identified by rDNA ITS sequence analysis, belonging to 33 genera, including Ophiocordyceps, Trichoderma, Cytospora, Truncatella, Dactylonectria, Isaria, Cephalosporium, Fusarium, Cosmospora and Paecilomyces, etc.. Among all culturable fungi, Mortierella and Trichoderma were predominant genera.

Conclusions

The significantly differences and overlap in fungal community structure between two approaches highlight that the integration of high-throughput sequencing and culture-dependent approaches would generate more information. Our result reveal a comprehensive understanding of fungal community structure of natural O. sinensis, provide new insight into O. sinensis associated fungi, and support that microbiota of natural O. sinensis is an untapped source for novel bioactive metabolites discovery.

Haplotype Diversity of NADPH-Cytochrome P450 Reductase Gene of Ophiocordyceps sinensis and the Effect on Fungal Infection in Host Insects

Ophiocordyceps sinensis Berk. is a fungal parasite that parasitizes the larvae of Hepialidae and is used as a traditional Chinese medicine. However, it is not clear how O. sinensis infects its host. The encoding gene haplotype diversity and predicted function of the nicotinamide adenine dinucleotide phosphate (NADPH)-cytochrome P450 reductase (CPR) related to the fungal pathogenicity was analyzed for 219 individuals from 47 O. sinensis populations. Two NADPH CPR genes of O. sinensis were detected and their dominant haplotypes were widely distributed throughout the entire distribution range in Western China. Only 5.43% of all O. sinensis individuals possessed the specific private haplotypes of NADPH CPR-1 and CPR-2 genes. Bioinformatic analyses predicted that the phosphorylation sites, motifs, and domains of NADPH CPR of O. sinensis were different between those encoding by the dominant and private gene haplotypes. The one-to-one match fungus–host correspondence of the same individual suggested that the widely distributed O. sinensis with the dominant NADPH CPR gene haplotypes may strongly infect almost all host insects through a random infection by oral or respiratory pores. Conversely, O. sinensis with the specific private NADPH CPR gene haplotypes is likely to infect only a few corresponding host insects by breaching the cuticle, due to the changed NADPH CPR structure and function.

Analysis of Volatile Components in Different Ophiocordyceps sinensis and Insect Host Products

The artificial production of Ophiocordyceps sinensis mycelia and fruiting bodies and the Chinese cordyceps has been established. However, the volatile components from these O. sinensis products are not fully identified. An efficient, convenient, and widely used approach based on headspace solid-phase microextraction (HS-SPME) combined with comprehensive two-dimensional gas chromatography and quadrupole time-of-flight mass spectrometry (GC×GC-QTOFMS) was developed for the extraction and the analysis of volatile compounds from three categories of 16 products, including O. sinensis fungus, Thitarodes hosts of O. sinensis, and the Chinese cordyceps. A total of 120 volatile components including 36 alkanes, 25 terpenes, 17 aromatic hydrocarbons, 10 ketones, 5 olefines, 5 alcohols, 3 phenols, and 19 other compounds were identified. The contents of these components varied greatly among the products but alkanes, especially 2,5,6-trimethyldecane, 2,3-dimethylundecane and 2,2,4,4-tetramethyloctane, are the dominant compounds in general. Three categories of volatile compounds were confirmed by partial least squares-discriminant analysis (PLS-DA). This study provided an ideal method for characterizing and distinguishing different O. sinensis and insect hosts-based products.

Thitarodes shambalaensis sp. nov. (Lepidoptera, Hepialidae): a new host of the caterpillar fungus Ophiocordyceps sinensis supported by genome-wide SNP data

A new species of ghost moth, Thitarodesshambalaensissp. nov., is described from Yanzigou glacier, Mt. Gongga, Sichuan, China. The species is a host of the economically important caterpillar fungus Ophiocordycepssinensis. Establishment of this new species is supported by morphology and genetic differentiation measured in a CO1 phylogeny and in genome-wide SNP coverage. A summary tree from 538 sequences of different genetic markers from Thitarodes (including sequences extracted from caterpillar fungus sclerotium samples) support the genus Thitarodes as a monophyletic group, and indicate that Thitarodes is the host genus for O.sinensis. Sampling efforts so far have centered on half of the known phylogenetic diversity of Thitarodes, with some species-level clusters (separated by < 2.5% genetic distance) containing 17 described species. Fifteen clusters are known from either a single “orphan taxon” or a single sequence from a caterpillar fungus sclerotium sample. We provide suggestions for building a more robust phylogeny of the genus Thitarodes and highlight some of the conservation threats that species from this genus face due to unprecedented habitat exploitation.

Internal and External Microbial Community of the Thitarodes Moth, the Host of Ophiocordyceps sinensis

Ophiocordyceps sinensis is a widely known medicinal entomogenous fungus, which parasitizes the soil-borne larva of Thitarodes (Hepialidae, Lepidoptera) distributed in the Qinghai–Tibetan Plateau and its adjacent areas. Previous research has involved artificial cultivation of Chinese cordyceps (the fungus-caterpillar complex), but it is difficult to achieve large-scale cultivation because the coupling relation between the crucial microbes and their hosts is not quite clear. To clarify the influence of the internal microbial community on the occurrence of Chinese cordyceps, in this study, the unfertilized eggs of Thitarodes of different sampling sites were chosen to analyze the bacterial and fungal communities via 16S rRNA and ITS sequencing for the first time. The results showed that for bacteria, 348 genera (dominant genera include Wolbachia, Spiroplasma, Carnobacterium, Sphingobium, and Acinetobacter) belonging to 26 phyla (dominant phyla include Proteobacteria, Firmicutes, Tenericutes, Actinobacteria, Acidobacteria, and Bacteroidetes), 58 classes, 84 orders, and 120 families were identified from 1294 operational taxonomic units (OTUs). The dominant bacterial genus (Spiroplasma) may be an important bacterial factor promoting the occurrence of Chinese cordyceps. For fungi, 289 genera, mainly including Aureobasidium, Candida, and Cryptococcus, were identified, and they belonged to 5 phyla (Ascomycota, Basidiomycota, Chytridiomycota, Glomeromycota, and Zygomycota), 26 classes, 82 orders, and 165 families. Eight bacterial OTUs and 12 fungal OTUs were shared among all of the detected samples and were considered as core species. Among them, Wolbachia, Spiroplasma, Carnobacterium, Aureobasidium, and Phoma may play important roles in helping the host larva to digest foods, adapt to extreme environments, or resist pathogens. On the other hand, the external (soil) microbial community was synchronously and comparatively analyzed. Comparative analysis revealed that external microbial factors might play a more significant role in the occurrence of Chinese cordyceps, owing to the significant differences revealed by α-diversity and β-diversity analyses among different groups. In summary, the results of this study may contribute to the large-scale cultivation of Chinese cordyceps.

Diversity and Co-Occurrence Patterns of Soil Bacterial and Fungal Communities of Chinese Cordyceps Habitats at Shergyla Mountain, Tibet: Implications for the Occurrence

Chinese Cordyceps is a well-known medicinal larva-fungus symbiote distributed in the Qinghai-Tibetan Plateau and adjacent areas. Previous studies have involved its artificial cultivation but commercial cultivation is difficult to perform because the crucial factors triggering the occurrence of Chinese Cordyceps are not quite clear. The occurrence of Chinese Cordyceps is greatly affected by the soil environment, including the soil’s physicochemical and microecological properties. In this study, the effects of these soil properties on the occurrence of Chinese Cordyceps were investigated. The results show that the physicochemical properties, including easily oxidizable organic carbon (EOC), soil organic carbon (SOC), humic acid carbon (HAC), humin carbon (HMC), and pH, might be negatively related to the occurrence of Chinese Cordyceps, and soil water content (SWC) might be positively related. Several soil physicochemical parameters (pH, SOC, HMC, HAC, available potassium (APO), available phosphorus (APH), microbial biomass carbon (MBC), and the ratio of NH₄⁺ to NO₃⁻ (NH₄⁺/NO₃⁻)) and microbial properties interact and mix together, which might affect the occurrence of Chinese Cordyceps. Soil microbial community structure was also a possible factor, and a low level of bacterial and fungal diversity was suitable for the occurrence of Chinese Cordyceps. The intra-kingdom network revealed that a closer correlation of the bacterial community might help the occurrence of Chinese Cordyceps, while a closer correlation of the fungal community might suppress it. The inter-kingdom network revealed that the occurrence rate of Chinese Cordyceps might be negatively correlated with the stability of the correlation state of the soil habitat. In conclusion, this study shows that soil physicochemical properties and microbial communities could be greatly related with the occurrence of Chinese Cordyceps. In addition, soil physicochemical properties, the level of bacterial and fungal diversity, and correlations of bacterial and fungal communities should be controlled to a certain level to increase the production of Chinese Cordyceps in artificial cultivation.

Complete mitochondrial genome of two Thitarodes species (Lepidoptera, Hepialidae), the host moths of Ophiocordyceps sinensis and phylogenetic implications

Thitarodes (Lepidoptera, Hepialidae) is the only genus that hosts to the Ophiocordyceps sinensis, a traditional Chinese medicine considered as a powerful medicinal supplement. In this study, the complete mitochondrial genomes (mitogenomes) of two species, T. damxungensis and T. pui, have been sequenced, which are 15,928 bp and 15,362 bp in size respectively, and both contain 13 protein-coding genes (PCGs), 2 rRNAs, 22 tRNAs and an AT-rich region. Like other hepialoids, the gene arrangement of the mitogenomes of T. damxungensis and T. pui is identical to the ancestral arrangement but differs from those of other lepidopteran species on account of the different arrangements of trnM, trnI, and trnQ. The size of AT-rich region is 545 bp in T. damxungensis and 1030 bp in T. pui. Tandem repetition in the AT-rich region is responsible for the length difference of the A + T-rich region in both species. In Hepialidae, the phylogenetic study based on the dataset of the sequences that combined the protein-coding genes and RNA genes suggested that the species T. yunnanensis should still belong to the genus Thitarodes rather than Ahamns, which is different from the results based on the traditional phylogeny.

Phylogeographic structures of the host insects of Ophiocordyceps sinensis

A fungus-insect complex, known as DongChong XiaCao, is formed from the infection of the hepialid larvae by the fungus Ophiocordyceps sinensis, which is endemic to the Qinghai-Tibetan Plateau (QTP). Due to previously limited sample collection size, the data about the diversity and structure of the host insect was insufficient and lacked details. The purpose of this study was aimed to discuss the diversity and phylogeography of the host insects of O. sinensis with a large-scale sampling. The mitochondrial cytochrome oxidase I gene (cox1) was sequenced and analyzed among 710 samples representing 88 geographic locations. 205 haplotypes of cox1 were identified from all the 710 samples and 4 phylogenetic clades with 12 subclades were identified. Instead of the single latitude-based divergence suggested previously, three distribution patterns were deduced to correspond to the phylogeographic structures, including but not limited to the co-existence of a wide and specific local phylogeographic distribution structures. Two separate genetic diversity and differentiation centers, namely the northwestern Yunnan and the southeastern Tibet were identified. Dating analyses from three calibrations supported that the divergence of the 4 clades occurred in the Oligocene-Miocene period (30.54-13.66 million years ago) (Ma), which were connected with the second and third geological movements of the QTP (17-25, 8-13 Ma). Our results provide a more detailed understanding of the divergence and distribution patterns of the host insects of O. sinensis.

Characterization of Humic Substances in the Soils of Ophiocordyceps sinensis Habitats in the Sejila Mountain, Tibet: Implication for the Food Source of Thitarodes Larvae

Humic substances in soil are considered to be an alternative food to the tender plant roots for Thitarodes larvae in the habitats of Ophiocordyceps sinensis in the Qinghai-Tibetan Plateau. However, there is no report involving the evaluation of their potential as a food source from the composition and structure of habitat soils. In this work, the composition and structure of humic substances in habitat soils from the Sejila Mountain, Tibet were characterized by diverse techniques for evaluating the nutritional value and possibility of humus as the food source for Thitarodes larvae. Fourier transform infrared spectroscopy revealed that humic acid may possess superior ability to provide the molecular segments for biosynthesizing lipids more than other humic fractions. Combining with the analysis of solid-state 13C nuclear magnetic resonance spectrum, the fractions of hydrophobic fulvic acid and hydrophilic fulvic acid are further considered as a potential food source for Thitarodes larvae. Overall, humic substances in habitat soils are rich in the molecular segments for biosynthesizing lipids and other important nutrients, which may provide the energy and material sources for maintaining the survival of Thitarodes larvae in the absence of tender plant roots, particularly in the annual cold winter. Combining with the evidence of physico-chemical parameters of habitat soils and stable carbon isotopic composition of major tender plant roots in the Sejila Mountain, the composition and structure of humic substances in habitat soils may provide a novel idea for the eco-friendly and semi-wild cultivation of Thitarodes larvae with low cost.

Discriminatory Power Evaluation of Nuclear Ribosomal RNA Barcoding Sequences Through Ophiocordyceps sinensis Related Samples

Since the cost of Ophiocordyceps sinensis has increased dramatically and the counterfeits may have adverse effect to health, a rapid and precise species-level DNA barcoding identification system could be a potent approach and significantly enhance the regulatory capacity. The discrimination power of three subunits sequences from nuclear ribosomal RNA gene cluster were determined by Simpson’s index of discrimination using 43 wild O. sinensis fruiting bodies, pure cultures, commercial mycelium fermented powder and counterfeits. The internal transcribed spacer (ITS) sequences showed the highest variance and discrimination power among 43 samples, as determined by Simpson’s index of discrimination (D = 0.972), followed by large subunit (LSU; D = 0.963) and small subunit (SSU; D = 0.921). ITS-2 sequences showed the highest discrimination power for 43 samples among ITS-1, ITS-2, and 5.8S region of ITS sequences. All O. sinensis samples were grouped into a unique ITS sequence cluster under 95% similarity and two O. sinensis samples and six non-O. sinensis samples showed false claims. Our data showed that the ITS region could provide accurate species identification for O. sinensis samples, especially when macroscopic and microscopic method could not be applied in the highly processed commercial products. Since the authentication of O. sinensis related products is essential to ensure its safety and efficacy, identification of O. sinensis through ITS sequence comparison or unique PCR amplification of the species specific target, such as the ITS region, should be considered in the next revision of Chinese pharmacopeia.

Discovery of differential sequences for improving breeding and yield of cultivated Ophiocordyceps sinensis through ITS sequencing and phylogenetic analysis

To accelerate the breeding process of cultivated Ophiocordyceps sinensis and increase its yield, it is important to identify molecular fingerprint of dominant O. sinensis. In the present study, we collected 3 batches of industrially cultivated O. sinensis product with higher yield than the others and compared their internal transcribed spacer (ITS) sequences with the wild and the reported. The ITS sequence was obtained by bidirectional sequencing and analyzed with molecular systematics as a DNA barcode for rapid and accurate identification of wild and cultivated O. sinensis collected. The ITS sequences of O. sinensis with detailed collection loci on NCBI were downloaded to construct a phylogenetic tree together with the sequences obtained from the present study by using neighbor-joining method based on their evolution relationship. The information on collection loci was analyzed with ArcGIS 10.2 to demonstrate the geographic distribution of these samples and thus to determine the origin of the dominant samples. The results showed that all wild and cultivated samples were identified as O. sinensis and all sequences were divided into seven phylogenetic groups in the tree. Those groups were precisely distributed on the map and the process of their system evolution was clearly presented. The three cultivated samples were clustered into two dominant groups, showing the correlation between the industrially cultivated samples and the dominant wild samples, which can provide references for its optimized breeding in the future.

Comparative Genomics Reveals the Core Gene Toolbox for the Fungus-Insect Symbiosis

Modern genomics has shed light on many entomopathogenic fungi and expanded our knowledge widely; however, little is known about the genomic features of the insect-commensal fungi. Harpellales are obligate commensals living in the digestive tracts of disease-bearing insects (black flies, midges, and mosquitoes). In this study, we produced and annotated whole-genome sequences of nine Harpellales taxa and conducted the first comparative analyses to infer the genomic diversity within the members of the Harpellales. The genomes of the insect gut fungi feature low (26% to 37%) GC content and large genome size variations (25 to 102 Mb). Further comparisons with insect-pathogenic fungi (from both Ascomycota and Zoopagomycota), as well as with free-living relatives (as negative controls), helped to identify a gene toolbox that is essential to the fungus-insect symbiosis. The results not only narrow the genomic scope of fungus-insect interactions from several thousands to eight core players but also distinguish host invasion strategies employed by insect pathogens and commensals. The genomic content suggests that insect commensal fungi rely mostly on adhesion protein anchors that target digestive system, while entomopathogenic fungi have higher numbers of transmembrane helices, signal peptides, and pathogen-host interaction (PHI) genes across the whole genome and enrich genes as well as functional domains to inactivate the host inflammation system and suppress the host defense. Phylogenomic analyses have revealed that genome sizes of Harpellales fungi vary among lineages with an integer-multiple pattern, which implies that ancient genome duplications may have occurred within the gut of insects.

Insect guts harbor various microbes that are important for host digestion, immune response, and disease dispersal in certain cases. Bacteria, which are among the primary endosymbionts, have been studied extensively. However, fungi, which are also frequently encountered, are poorly known with respect to their biology within the insect guts. To understand the genomic features and related biology, we produced the whole-genome sequences of nine gut commensal fungi from disease-bearing insects (black flies, midges, and mosquitoes). The results show that insect gut fungi tend to have low GC content across their genomes. By comparing these commensals with entomopathogenic and free-living fungi that have available genome sequences, we found a universal core gene toolbox that is unique and thus potentially important for the insect-fungus symbiosis. This comparative work also uncovered different host invasion strategies employed by insect pathogens and commensals, as well as a model system to study ancient fungal genome duplication within the gut of insects.

Stable Carbon Isotope Composition of the Lipids in Natural Ophiocordyceps sinensis from Major Habitats in China and Its Substitutes

Ophiocordyceps sinensis is one rare medicinal fungus produced in the Qinghai-Tibetan Plateau. Its quality and price varies hugely with different habitat, and its numerous substitutes have sprung up in functional food markets. This paper aims to discriminate the geographic origin of wild O. sinensis and its substitutes via element analyzer-isotope ratio mass spectrometry and gas chromatography-isotope ratio mass spectrometry. The δ13C values of major fatty acids in the lipids of O. sinensis are characterized unanimously by the variation relation C18:0 < C18:2 ≈ C16:0 < C18:1, while their fluctuation intervals are notably different between those of neutral and polar lipids. The comparative analysis of the δ13C ratios of major fatty acids in lipids of O. sinensis suggests that the δ13C patterns may be sensitive potential indicators to discriminate its geographical origin. The δ13C values of individual major fatty acids of lipids from the cultivated stromata of Cordyceps militaris (SCM), the fermented mycelia of Hirsurella sinensis (FMH) and Paecilomyces epiali (FMP) range from -31.2‰ to -29.7‰, -16.9‰ to -14.3‰, and -26.5‰ to -23.9‰, respectively. Their δ13C pattern of individual major fatty acids may be used as a potential indicator to discriminate the products of natural O. sinensis and its substitutes.

The Mechanisms of Pharmacological Activities of Ophiocordyceps sinensis Fungi

The entomopathogenic fungus Ophiocordyceps sinensis, formerly known as Cordyceps sinensis, has long been used as a traditional Chinese medicine for the treatment of many illnesses. In recent years its usage has increased dramatically because of the improvement of people's living standard and the emphasis on health. Such demands have resulted in over-harvesting of this fungus in the wild. Fortunately, scientists have demonstrated that artificially cultured and fermented mycelial products of O. sinensis have similar pharmacological activities to wild O. sinensis. The availability of laboratory cultures will likely to further expand its usage for the treatment of various illnesses. In this review, we summarize recent results on the pharmacological activities of the components of O. sinensis and their putative mechanisms of actions. Copyright © 2016 John Wiley & Sons, Ltd.

Isolation and Classification of Fungal Whitefly Entomopathogens from Soils of Qinghai-Tibet Plateau and Gansu Corridor in China

Qinghai-Tibet Plateau and Gansu Corridor of China with distinct geographic and climatic conditions are remote and less disturbed by humans, in which are likely to find some new strains of fungal entomopathogens against B-biotype whiteflies that is a very important invading pest worldwide. In this research, nineteen strains among six species of entomogenous fungi were isolated from the soil samples collected from 32 locations in Qinghai-Tibet Plateau and Gansu Corridor. From the data of isolation rates, it was indicated that the good biodiversity of entomogenous fungi was found in the soil covered good vegetations. On the contrary, no strains were isolated from the desert areas. In addition, the dominant species, Isaria fumosorosea and Metarhizium anisopliae var. anisopliae in the Qinghai-Tibet Plateau are different from the strains of other places based on ITS genetic homology analysis. It was verified that the Qinghai-Tibet Plateau area was less disturbed by human, and the fungi in this place exchanged less compared with other regional species. All of these strains showed the pathogenicity against the B-biotype whitefly with the mortality of more than 30%. However, a few strains of Paecilomyces lilacinus, Lecanicillium psalliotae, Aspergillus ustus, I. fumosorosea and M. anisopliae var. anisopliae had better virulence with LC₅₀s of 0.36–26.44×10⁶ spores/mL on post-treatment day 6–7. Especially, the L. psalliotae strain LpTS01 was the greatest virulence with LC₅₀ of 0.36×10⁶spores/mL and LT₅₀ of 4.23d. Our research thus presents some new insights to discover new entomopathogenic fungal strains used for B-biotype whitefly biocontrol.

Insights from the genome of Ophiocordyceps polyrhachis-furcata to pathogenicity and host specificity in insect fungi

Background

Ophiocordyceps unilateralis is an outstanding insect fungus for its biology to manipulate host ants’ behavior and for its extreme host-specificity. Through the sequencing and annotation of Ophiocordyceps polyrhachis-furcata, a species in the O. unilateralis species complex specific to the ant Polyrhachis furcata, comparative analyses on genes involved in pathogenicity and virulence between this fungus and other fungi were undertaken in order to gain insights into its biology and the emergence of host specificity.

Results

O. polyrhachis-furcata possesses various genes implicated in pathogenicity and virulence common with other fungi. Overall, this fungus possesses protein-coding genes similar to those found on other insect fungi with available genomic resources (Beauveria bassiana, Metarhizium robertsii (formerly classified as M. anisopliae s.l.), Metarhizium acridum, Cordyceps militaris, Ophiocordyceps sinensis). Comparative analyses in regard of the host ranges of insect fungi showed a tendency toward contractions of various gene families for narrow host-range species, including cuticle-degrading genes (proteases, carbohydrate esterases) and some families of pathogen-host interaction (PHI) genes. For many families of genes, O. polyrhachis-furcata had the least number of genes found; some genes commonly found in other insect fungi are even absent (e.g. Class 1 hydrophobin). However, there are expansions of genes involved in 1) the production of bacterial-like toxins in O. polyrhachis-furcata, compared with other entomopathogenic fungi, and 2) retrotransposable elements.

Conclusions

The gain and loss of gene families helps us understand how fungal pathogenicity in insect hosts evolved. The loss of various genes involved throughout the pathogenesis for O. unilateralis would result in a reduced capacity to exploit larger ranges of hosts and therefore in the different level of host specificity, while the expansions of other gene families suggest an adaptation to particular environments with unexpected strategies like oral toxicity, through the production of bacterial-like toxins, or sophisticated mechanisms underlying pathogenicity through retrotransposons.

Electronic supplementary material

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

Genomics and Evolution in Traditional Medicinal Plants: Road to a Healthier Life

Medicinal plants have long been utilized in traditional medicine and ethnomedicine worldwide. This review presents a glimpse of the current status of and future trends in medicinal plant genomics, evolution, and phylogeny. These dynamic fields are at the intersection of phytochemistry and plant biology and are concerned with the evolution mechanisms and systematics of medicinal plant genomes, origin and evolution of the plant genotype and metabolic phenotype, interaction between medicinal plant genomes and their environment, the correlation between genomic diversity and metabolite diversity, and so on. Use of the emerging high-end genomic technologies can be expanded from crop plants to traditional medicinal plants, in order to expedite medicinal plant breeding and transform them into living factories of medicinal compounds. The utility of molecular phylogeny and phylogenomics in predicting chemodiversity and bioprospecting is also highlighted within the context of natural-product-based drug discovery and development. Representative case studies of medicinal plant genome, phylogeny, and evolution are summarized to exemplify the expansion of knowledge pedigree and the paradigm shift to the omics-based approaches, which update our awareness about plant genome evolution and enable the molecular breeding of medicinal plants and the sustainable utilization of plant pharmaceutical resources.

Morphological Observations and Fatty Acid Composition of Indoor-Cultivated Cordyceps sinensis at a High-Altitude Laboratory on Sejila Mountain, Tibet

Cordyceps sinensis, a caterpillar entomopathogenic fungus-host larva complex, is a rare medicinal herb found in the Qinghai-Tibetan Plateau and its surrounding high-altitude areas. The alternation of generations in the life cycle, whatever the fungus or its host insect, requires special growth conditions. However, it is difficult to simulate the growth conditions of C. sinensis, which hinders its artificial cultivation. In this work, the life cycle from the host larva to C. sinensis was observed in an indoor-cultivation laboratory at 4,200 m a.s.l. on Sejila Mountain, Tibet. Comparative examinations between indoor-cultivated and wild C. sinensis demonstrated that the indoor-cultivated C. sinensis preferred to germinate multiple long, slim stromata at diverse positions on dead larvae, including but not limited to their heads. Their fatty acid composition shows a significant difference in the levels of polyunsaturated fatty acids (PUFAs). In indoor-cultivated C. sinensis, PUFAs constituted 24.59% and 49.43%, respectively, of neutral and polar lipids; meanwhile, in wild C. sinensis, PUFAs represented 34.34% and 61.25% of neutral and polar lipids, respectively. These observations and fatty acid data suggest that environmental factors, particularly temperature, soil pressure and light intensity, strongly affect the growth of C. sinensis. Our new findings may provide important information for improving techniques for the large-scale artificial cultivation of C. sinensis.