Vegetative development and host immune interaction of Ophiocordyceps sinensis within the hemocoel of the ghost moth larva, Thitarodes xiaojinensis

Invasion of the four kingdoms: the parasite journey across plant and non-plant hosts

Parasites have a rich and long natural history among biological entities, and it has been suggested that parasites are one of the most significant factors in the evolution of their hosts. However, it has been emphasized less frequently how co-evolution has undoubtedly also shaped the paths of parasites. It may seem safe to assume that specific differences among the array of potential hosts for particular parasites have restricted and diversified their evolutionary pathways and strategies for survival. Nevertheless, if one looks closely enough at host and parasite, one finds commonalities, both in terms of host defences and parasite strategies to out-manoeuvre them. While such analyses have been the source of numerous reviews, they are generally limited to interactions between, at most, one kingdom of parasite with two kingdoms of host (e.g. similarities in animal and plant host responses against fungi). With the aim of extending this view, we herein critically evaluate the similarities and differences across all four eukaryotic host kingdoms (plants, animals, fungi, and protists) and their parasites. In doing so, we show that hosts tend to share common strategies for defence, including both physical and behavioural barriers, and highly evolved immune responses, in particular innate immunity. Parasites have, similarly, evolved convergent strategies to counter these defences, including mechanisms of active penetration, and evading the host's innate and/or adaptive immune responses. Moreover, just as hosts have evolved behaviours to avoid parasites, many parasites have adaptations to manipulate host phenotype, physiologically, reproductively, and in terms of behaviour. Many of these strategies overlap in the host and parasite, even across wide phylogenetic expanses. That said, specific differences in host physiology and immune responses often necessitate different adaptations for parasites exploiting fundamentally different hosts. Taken together, this review facilitates hypothesis-driven investigations of parasite-host interactions that transcend the traditional kingdom-based research fields.

Neurotransmitter acetylcholine mediates the mummification of Ophiocordyceps sinensis-infected Thitarodes xiaojinensis larvae

Parasites can manipulate host behavior to facilitate parasite transmission. One such host-pathogen interaction occurs between the fungus Ophiocordyceps sinensis and the ghost moth Thitarodes xiaojinensis. O. sinensis is involved in the mummification process of infected host larvae. However, the underlying molecular and chemical mechanism for this phenomenon is unknown. We characterized the small molecules regulating host behaviors and the altered metabolites in infected and mummified host larvae. Lipid-related metabolites, such as phosphatidylcholine, were identified in infected and mummified larvae. Decreased levels of the neurotransmitter acetylcholine (ACh) and elevated choline levels were observed in the brains of both the infected and mummified larvae. The aberrant activity of acetylcholinesterase (AChE) and relative mRNA expression of ACE2 (acetylcholinesterase) may mediate the altered transformation between ACh and choline, leading to the brain dysfunction of mummified larvae. Caspofungin treatment inhibited the mummification of infected larvae and the activity of AChE. These findings indicate the importance of ACh in the mummification of host larvae after O. sinensis infection.IMPORTANCEOphiocordyceps sinensis-infected ghost moth larvae are manipulated to move to the soil surface with their heads up in death. A fruiting body then grows from the caterpillar's head, eventually producing conidia for dispersal. However, the underlying molecular and chemical mechanism has not been characterized. In this study, we describe the metabolic profile of Thitarodes xiaojinensis host larvae after O. sinensis infection. Altered metabolites, particularly lipid-related metabolites, were identified in infected and mummified larvae, suggesting that lipids are important in O. sinensis-mediated behavioral manipulation of host larvae. Decreased levels of the neurotransmitter acetylcholine were observed in both infected and mummified larvae brains. This suggests that altered or reduced acetylcholine can mediate brain dysfunction and lead to aberrant behavior. These results reveal the critical role of acetylcholine in the mummification process of infected host larvae.

Filtration effect of Cordyceps chanhua mycoderm on bacteria and its transport function on nitrogen

IMPORTANCE

During the natural growth of Cordyceps chanhua, it will form a mycoderm structure specialized from hyphae. We found that the bacterial membrane of C. chanhua not only filters environmental bacteria but also absorbs and transports nitrogen elements inside and outside the body of C. chanhua. These findings are of great significance for understanding the stable mechanism of the internal microbial community maintained by C. chanhua and how C. chanhua maintains its own nutritional balance. In addition, this study also enriched our understanding of the differences in bacterial community composition and related bacterial community functions of C. chanhua at different growth stages, which is of great value for understanding the environmental adaptation mechanism, the element distribution network, and the changing process of symbiotic microbial system after Cordyceps fungi infected the host. At the same time, it can also provide a theoretical basis for some important ecological imitation cultivation technology of Cordyceps fungi.

Histopathological observations and comparative transcriptome analysis of Ophiocordyceps sinensis infection of Hepialus xiaojinensis in the early stage

Hepialus xiaojinensis is a Lepidopteran insect and one of the hosts for the artificial cultivation of Cordyceps. Ophiocordyceps sinensis can infect and coexist with H. xiaojinensis larvae for a long time. Little studies focused on the interaction process through its early infection stage. In this research, we particularly study the interaction of infected and uninfected larvae in the 3rd (OS-3, CK-3) and 4th (OS-4, CK-4) instars. O. sinensis was distributed within the larvae and accompanied by pathological changes in some tissue structures. In response to O. sinensis infection, OS-3 enhanced the antioxidant defense ability, while OS-4 decreased. The transcriptome analysis showed that OS-3 resisted the invasion of O. sinensis by the immune and nervous systems. Correspondingly, OS-4 reduced immune response and utilized more energy for growth and development. This study provides a comprehensive resource for analyzing the mechanism of H. xiaojinensis and O. sinensis interaction.

Transcriptome and Metabolome Analyses of Thitarodes xiaojinensis in Response to Ophiocordyceps sinensis Infection

Ophiocordyceps sinensis exhibits more than 5 months of vegetative growth in Thitarodes xiaojinensis hemocoel. The peculiar development process of O. sinensis has been elucidated through morphological observation and omics technology; however, little information has been reported regarding the changes that occur in the host T. xiaojinensis. The RNA sequencing data showed that when O. sinensis blastospores were in the proliferative stage, the greatest change in the infected larval fat body was the selectively upregulated immune recognition and antimicrobial peptide genes. When O. sinensis blastospores were in the stationary stage, the immune pathways of T. xiaojinensis reverted to normal levels, which coincides with the successful settlement of O. sinensis. Pathway enrichment analysis showed a higher expression of genes involved in energy metabolism pathway in this stage. Metabolomic analyses revealed a reduction of amino acids and lipids in hemolymph, but an upregulation of lipids in the fat body of the host larvae after O. sinensis infection. We present the first transcriptome integrated with the metabolome study of T. xiaojinensis infected by O. sinensis. It will improve our understanding of the interaction mechanisms between the host and entomopathogenic fungi, and facilitate future functional studies of genes and pathways involved in these interactions.

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.

Transcriptomic Analysis Insight into the Immune Modulation during the Interaction of Ophiocordyceps sinensis and Hepialus xiaojinensis

Ophiocordyceps sinensis (Berk.) is an entomopathogenic fungus that can infect the larva of the ghost moth, Hepialus xiaojinensis, causing mummification after more than one year. This prolonged infection provides a valuable model for studying the immunological interplay between an insect host and a pathogenic fungus. A comparative transcriptome analysis of pre-infection (L) and one-year post-infection (IL) larvae was performed to investigate the immune response in the host. Here, a total of 59,668 unigenes were obtained using Illumina Sequencing in IL and L. Among the 345 identified immune-related genes, 83 out of 86 immune-related differentially expressed genes (DEGs) had a much higher expression in IL than in L. Furthermore, the immune-related DEGs were classified as pathogen recognition receptors (PRRs), signal modulators or transductors, and immune effector molecules. Serpins and protease inhibitors were found to be upregulated in the late phase of infection, suppressing the host’s immune response. Based on the above analysis, the expression levels of most immune-related genes would return to the baseline with the immune response being repressed in the late phase of infection, leading to the fungal immunological tolerance after prolonged infection. Meanwhile, the transcriptomes of IL and the mummified larva (ML) were compared to explore O. sinensis invasion. A total of 1408 novel genes were identified, with 162 of them annotated with putative functions. The gene families likely implicated in O. sinensis pathogenicity have been identified, primarily including serine carboxypeptidase, peroxidase, metalloprotease peptidase, aminopeptidases, cytochrome P450, and oxidoreductase. Furthermore, quantitative real-time PCR (qPCR) was used to assess the expression levels of some critical genes that were involved in immune response and fungal pathogenicity. The results showed that their expression levels were consistent with the transcriptomes. Taken together, our findings offered a comprehensive and precise transcriptome study to understand the immune defense in H. xiaojinensis and O. sinensis invasion, which would accelerate the large-scale artificial cultivation of this medicinal fungus.

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.

Immulectin-2 from the ghost moth, Thitarodes xiaojinensis (Lepidoptera: Hepialidae), modulates cellular and humoral responses against fungal infection

C type-lectins constitute a large family of pattern recognition receptors, playing important roles in insect immune defenses. Thitarodes xiaojinensis larvae showed distinct immune features after Ophiocordyceps sinensis, Cordyceps militaris, or Beauveria bassiana infection. Based on transcriptome and immunoblot analysis, we found that immulectin-2 (IML2) was induced after T. xiaojinensis larvae were infected by C. militaris or B. bassiana but maintained at a low level after larvae injected with O. sinensis or Ringer's buffer. Recombinant IML2 (rIML2) could promote melanization, encapsulation, phagocytosis, and hemocyte aggregation in vitro. RNA interference with IML2 induced a significant reduction in the transcript levels of various antimicrobial peptides. Importantly, we found that the abundance of O. sinensis blastospores coated with rIML2 dramatically decreased in the host hemolymph. Overall, this study demonstrated that T. xiaojinensis IML2 modulates cellular and humoral responses to entomopathogenic fungi, broadening our view of the immune interaction between O. sinensis and its host.

Distinct Responses of Thitarodes xiaojinensis β-1,3-Glucan Recognition Protein-1 and Immulectin-8 to Ophiocordyceps sinensis and Cordyceps militaris Infection

Melanization and encapsulation are prominent defense responses against microbes detected by pattern recognition receptors of their host insects. In the ghost moth Thitarodes xiaojinensis, an activated immune system can melanize and encapsulate the fungus Cordyceps militaris However, these responses were hardly detected in the host hemolymph postinfection of another fungus Ophiocordyceps sinensis The immune interaction between O. sinensis and the host remains largely unknown, which hinders the artificial cultivation of Chinese cordyceps. We found that T. xiaojinensis β-1,3-glucan recognition protein-1 (βGRP1) was needed for prophenoloxidase activation induced by C. militaris Failure of βGRP1 to recognize O. sinensis is a primary reason for the lack of melanization in the infected host. Lyticase or snailase treatment combined with binding and immunofluorescence detection showed the existence of a protective layer preventing the fungus from βGRP1 recognition. Coimmunoprecipitation and mass spectrometry analysis indicated that βGRP1 interacted with immulectin-8 (IML8) via binding to C. militaris IML8 promotes encapsulation. This study suggests the roles of T. xiaojinensis βGRP1 and IML8 in modulating immune responses against C. militaris Most importantly, the data indicate that O. sinensis may evade melanization by preventing βGRP1 recognition.

iTRAQ-based comparative proteome analyses of different growth stages revealing the regulatory role of reactive oxygen species in the fruiting body development of Ophiocordyceps sinensis

In this study, using an isobaric tags for relative and absolute quantitation (iTRAQ ) approach coupled with LC-MS / MS and bioinformatics, the proteomes were analyzed for the crucial three stages covering the fruiting body development of Ophiocordyceps sinensis, including sclerotium (ST), primordium (PR) and mature fruiting body (MF), with a focus on fruiting body development-related proteins and the potential mechanisms of the development. A total of 1,875 proteins were identified. Principal Component Analysis (PCA) demonstrated that the protein patterns between PR and MF were more similar than ST. Differentially accumulated proteins (DAPs) analysis showed that there were 510, 173 and 514 DAPs in the comparisons of ST vs. PR, PR vs. MF and ST vs. MF, respectively. A total of 62 shared DAPs were identified and primarily enriched in proteins related to ‘carbon transport and mechanism’, ‘the response to oxidative stress’, ‘antioxidative activity’ and ‘translation’. KEGG and GO databases showed that the DAPs were enriched in terms of ‘primary metabolisms (amino acid/fatty acid/energy metabolism)’, ‘the response to oxidative stress’ and ‘peroxidase’. Furthermore, 34 DAPs involved in reactive oxygen species (ROS) metabolism were identified and clustered across the three stages using hierarchical clustering implemented in hCluster R package . It was suggested that their roles and the underlying mechanisms may be stage-specific. ROS may play a role in fungal pathogenicity in ST, the fruit-body initiation in PR, sexual reproduction and highland adaptation in MF. Crucial ROS-related proteins were identified, such as superoxide dismutase (SOD, T5A6F1), Nor-1 (T5AFX3), electron transport protein (T5AHD1), histidine phosphotransferase (HPt, T5A9Z5) and Glutathione peroxidase (T5A9V1). Besides, the accumulation of ROS at the three stages were assayed using 2,7-dichlorofuorescin diacetate (DCFH-DA) stanning. A much stronger ROS accumulation was detected at the stage MF, compared to the stages of PR and ST. Sections of ST and fruit-body part of MF were stained by DCFH-DA and observed under the fluorescencemicroscope, showing ROS was distributed within the conidiospore and ascus. Besides, SOD activity increased across the three stages, while CAT activity has a strong increasement in MF compared to the stages of ST and PR. It was suggested that ROS may act in gradient-dependent manner to regulate the fruiting body development. The coding region sequences of six DAPs were analyzed at mRNA level by quantitative real-time PCR (qRT-PCR). The results support the result of DAPs analysis and the proteome sequencing data. Our findings offer the perspective of proteome to understand the biology of fruiting body development and highland adaptation in O. sinensis, which would inform the big industry of this valuable fungus.

Changes in transcriptomic and metabolomic profiles of morphotypes of Ophiocordyceps sinensis within the hemocoel of its host larvae, Thitarodes xiaojinensis

Background

Ophiocordyceps sinensis (Berk.) is a well-known entomopathogenic and medicinal fungus. It parasitizes and mummifies the underground ghost moth larvae to produce a fruiting body named Chinese cordyceps. Specific for the fungus, O. sinensis experiences a biotrophic vegetative growth period spanning over 5 months. During this vegetative growth, it appears successively in the host hemocoel in three/four morphotypes, namely, the yeast-like blastospores (subdivided into proliferative (BP) and stationary phase (BS)), prehyphae (PreHy) and the hyphae (Hy). This peculiar morphogenesis has been elucidated through morphological and ultrastructural observations, but its molecular basis remains cryptic. In this study, transcriptome and metabolome profiling of BP, BS, PreHy and Hy stages were performed to characterize the key genes, metabolites, and signaling pathways that regulated the vegetative development of O. sinensis in Thitarodes xiaojinensis larva.

Results

The molecular events and metabolic pathways that regulated different intracellular processes at various stages were examined. Cluster analyses of differentially expressed genes across the four stages revealed the stage specifically enriched pathways. Analysis of metabolome profiles showed that carbon metabolism and several amino acids biosynthesis were significantly perturbed during the tested development stages of O. sinensis in the host hemocoel. Genes homologous to Saccharomyces cerevisiae MAPK cascade were significantly up-regulated during the transition from blastospore to hypha. The up-regulation of Sho1, a regulator protein, suggested nutrient starvation act a role in activation of MAPK pathway and filamentous growth. In addition, up-regulation of several fatty acid synthesis genes and their corresponding products accumulation in the samples of BS might explain more lipid droplets were observed in BS than in BP. Coupled with the up-regulation of fatty acid degradation during PreHy and Hy stages, it is presumed that lipid accumulation and mobilization play important roles in filamentous development.

Conclusions

This is the first report comprehensively describing developmental transcriptomics and metabolomics of O. sinensis in vivo. Our findings provide new perspectives into the key pathways and hub genes involved in morphological changes of fungus developed in the hemocoel of its host, and are expected to guide future studies on morphogenesis and morphotype changes of entomopathogenic fungi in vivo.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07209-2.

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

Ophiocordyceps sinensis (Berk.) is an entomopathogenic fungus endemic to the Qinghai-Tibet Plateau. It parasitizes and mummifies the underground ghost moth larvae, then produces a fruiting body. The fungus-insect complex, called Chinese cordyceps or "DongChongXiaCao," is not only a valuable traditional Chinese medicine, but also a major source of income for numerous Himalayan residents. Here, taking advantage of rapid advances in single-molecule sequencing, we assembled a highly contiguous genome assembly of O. sinensis. The assembly of 23 contigs was ∼110.8 Mb with a N50 length of 18.2 Mb. We used RNA-seq and homologous protein sequences to identify 8,916 protein-coding genes in the IOZ07 assembly. Moreover, 63 secondary metabolite gene clusters were identified in the improved assembly. The improved assembly and genome features described in this study will further inform the evolutionary study and resource utilization of Chinese cordyceps.