Proteome sequencing and analysis of Ophiocordyceps sinensis at different culture periods
BMC Genomics volume 21, Article number: 886 (2020)
Ophiocordyceps sinensis is an important traditional Chinese medicine for its comprehensive active ingredients, such as cordycepin, cordycepic acid, and Cordyceps polysaccharide. O. sinensis zjut, a special strain isolated from O. sinensis, has similar pharmacological functions to wild O. sinensis. Currently, O. sinensis with artificial cultivation has been widely studied, but systematic fundamental research at protein levels has not been determined.
Proteomes of O. sinensis zjut at different culture periods (growth period, 3rd day; pre-stable period, 6th day; and stable period, 9th day) were relatively quantified by relative isotope markers and absolute quantitative technology. In total, 4005 proteins were obtained and further annotated with Gene Ontology, Kyoto Encyclopedia of Genes and Genomes database. Based on the result of the annotations, metabolic pathways of active ingredients, amino acids and fatty acid were constructed, and the related enzymes were exhibited. Subsequently, comparative proteomics of O. sinensis zjut identified the differentially expressed proteins (DEPs) by growth in different culture periods, to find the important proteins involved in metabolic pathways of active ingredients. 605 DEPs between 6d-VS-3d, 1188 DEPs between 9d-VS-3d, and 428 DEPs between 9d-VS-6d were obtained, respectively.
This work provided scientific basis to study protein profile and comparison of protein expression levels of O. sinensis zjut, and it will be helpful for metabolic engineering works to active ingredients for exploration, application and improvement of this fungus.
Ophiocordyceps sinensis is an important traditional Chinese medicine and healthy food in China [1,2,3]. Previous studies have revealed that active ingredients of O. sinensis, including D-mannitol, cordycepin, purine nucleotides and polysaccharide, have various pharmacological functions, such as adaptogenic activity, immunomodulatory effects, antioxidant activity and anti-cancer [4, 5]. However, due to high demand and insufficient supply for wild O. sinensis, the medicinal value development has been seriously limited [6, 7]. Recently, several strains that isolated from wild O. sinensis have similar pharmacological functions to wild O. sinensis , representing an useful alternative for production. Moreover, new potential drugs and active compounds derived from natural sources have been screened as effective disease treatment from Hirsutella sinensis, which is the potential anamorph of O. sinensis. The antitumor activity of H. sinensis mycelium was found as seen in studies on human tissue, including prostate (PC3), breast (MCF7), hepatocellular (HepG2, Hep3B) and colorectal (HT-29) [9, 10]. Previous studies have revealed the active ingredients of artificially-cultivated O. sinensis had extensive medicinal value for human health, such as anti-fatigue activity , immunomodulatory activity , antioxidant activity  and anti-obesity effects . Furthermore, the productions of active ingredients in H. sinensis by submerged fermentation were promoted and satisfying on the basis of genetic study [15,16,17].
Proteomics is the large-scale study of proteins, and it can be divided into the areas of large-scale identification of proteins and their post-translational modifications, comparison of protein expression levels, and protein-protein interactions . In previous studies, biochemical methods have been widely used to study proteome, including polyacrylamide gel electrophoresis (PAGE), two-dimensional electrophoresis (2-DE) liquid chromatography, surface-enhanced laser desorption ionization-time of flight-mass spectrometry (SELDI-TOF-MS), and matrix-assisted laser desorption ionization-time of flight-mass spectrometry (MALDI-TOF-MS) [19, 20]. Recently, with the ability to carry out relative (or absolute) quantification in up to eight phenotypes, isobaric tags for relative and absolute quantification (iTRAQ) have caught the attention of proteomics community [21, 22]. iTRAQ-based proteomics has been used in studying secretome, plasma membrane proteome, and intracellular proteome [22,23,24].
Genome sequence of O. sinensis has been reported, which revealed the pathogenic mechanism during life cycle . Then, transcriptomes of H. sinensis at diffident culture periods have been sequenced and analyzed to describe metabolic pathway and infection mechanism . Meanwhile, the comparative proteomic has shown a snapshot proteome profile and revealed the similarity of the proteins and metabolites composition between naturally- and artificially-cultivated of O. sinensis . Dong et al studied the dynamic polymorphic alterations among differentially expressed proteins of multiple intrinsic fungi in the caterpillar body and stroma of natural O. sinensis during maturation, and the results revealed there were the apparent proteomic polymorphism dissimilarity of these organisms to support the integrated micro-ecosystem hypothesis for natural O. sinensis . In addition, there were some non-quantitative studies about the proteomics of O. sinensis [27,28,29,30]. In fact, without the complement of proteomics, only genome and transcriptome sequences are not sufficient to elucidate biological functions. Moreover, there is no strict linear relationship between genes and the corresponding proteins . Therefore, proteome of the O. sinensis zjut was used to the large-scale study of gene functions directly at the protein level and tried to illuminate the synthesis mechanism of active ingredients in this study.
In this study, proteomes of the O. sinensis zjut at different culture periods were relatively quantified by iTRAQ together with two-dimensional liquid chromatography tandem mass spectrometry (LC-MS/MS). Large-scale identification and different expression analysis of the proteins were performed, and metabolic pathways were constructed, especially active ingredients.
Protein profiling and iTRAQ quantification
In order to obtain the overview of the O. sinensis zjut proteome, the protein samples were prepared from the mycelium at different culture periods (3d, 6d and 9d) and relatively quantified by iTRAQ. A total of 371,999 spectra were generated and 4005 proteins included 22,202 peptides were identified (with 1% FDR) (Supplementary information 3: Table S2). GO analysis of total proteins was based on biological process, cellular component and molecular function, and all proteins were classified into 41 functional groups (Fig. 1). Biological processes were associated with the following brief pathways: metabolic process (29.27%), cellular process (25.99%) and single-organism process (17.72%); cellular components were assigned to the following cellular compartments: cell (24.98%), cell part (24.98%) and organelle (15.73%); The most highly enriched molecular functions were binding (48.34%), followed by catalytic activity (40.73%). Moreover, the proteins were further classified into 24 functional categories using COG classifications (Fig. 2). The largest group category was found to have general functions only (18.68%), followed by translation, ribosomal structure and biogenesis (9.00%), and posttranslational modification, protein turnover, chaperones (7.72%). Only a small fraction of the protein was functionally related to the categories of cell motility (0.22%) and nuclear structure (0.06%). About 2.95% of the identified proteins (94 proteins) were related to secondary metabolites biosynthesis, transport and catabolism.
In this study, 6d-VS-3d, 9d-VS-3d, and 9d-VS-6d were set as comparison groups. Corresponding proteins, exhibiting a greater than 1.2-fold change and Q-value less than 0.05, were defined as differentially expressed proteins. Ultimately, 605 DEPs (340 up-regulated proteins and 265 down-regulated proteins) between 6d-VS-3d, 1188 DEPs (545 up-regulated proteins and 643 down-regulated proteins) between 9d-VS-3d, and 428 DEPs (215 up-regulated proteins and 213 down-regulated proteins) between 9d-VS-6d were obtained (Supplementary information 3: Table S3).
GO and pathway enrichment analysis of DEPs
The GO enrichment analysis of DEPs was shown in Supplementary information 4 (Figure S2, Figure S3, Figure S4). The DEPs were further analyzed by KEGG. As shown in Fig. 3, among the comparison group of 6d-VS-3d, the DEPs were mainly enriched in “ribosome” (39 members), followed by “tryptophan metabolism” (12 members), “arginine and proline metabolism” (10 members), “tyrosine metabolism” (9 members), and “phenylalanine metabolism” (8 members). For the comparison group of 9d-VS-6d, the DEPs were mainly enriched in “carbon metabolism” (19 members) and” phenylalanine metabolism” (6 members). Among the comparison group of 9d-VS-3d, the DEPs were mainly enriched in “ribosome” (43 members), followed by “tyrosine metabolism” (15 members), “phenylalanine metabolism” (12 members), “tryptophan metabolism” (16 members), “arginine and proline metabolism” (17 members), “alanine, aspartate and glutamate metabolism” (17 members), “one carbon pool by folate” (10 members), and “glyoxylate and dicarboxylate metabolism” (15 members).
Metabolic pathways of D-mannitol, cordycepin and purine nucleotides
Based on the glycolytic pathway and fructose-mannose pathway, the biosynthetic pathway of D-mannitol was successfully constructed (Fig. 4), and the average content of D-mannitol was 10.12% (3rd day 7.16%, 6th day 8.63%, and 9th day 14.58%, Supplementary information 1: Table S1). According to the constructed pathway, the initial precursors of D-mannitol are various saccharide, such as glucose, fructose and mannose. A total of 18 proteins, including 7 DEPs (two up-regulated HK, one up-regulated FBP, one up-regulated mtlD, two down-regulated manA and one down-regulated FBP), were located in the D-mannitol metabolic pathway (Supplementary information 5: Table S4). After several catalytic steps, β-D-fructose-6P is formed, and then converts to D-mannitol-1P which is the immediate precursor of D-mannitol. Unfortunately, we did not find mannitol-1-phosphate phosphatase which catalyzes D-mannitol-1-P to D-mannitol. The biosynthetic pathway of cordycepin, which might originate from histidine and end up with 3′-deoxyadenosine, was constructed (Fig. 5), and the average content of cordycepin was 0.308 mg/g (3rd day 0.327 mg/g, 6th day 0.140 mg/g, and 9th day 0.457 mg/g, Supplementary information 1: Table S1). In the constructed pathway, cordycepin is synthesis following the biosynthesis of adenosine. AMP, the precursor of adenosine, can convert to ADP by the catalysis of adk, and then RRM1 or RRM2 catalyzes ADP to form 3′-dADP. 3′-dADP dephosphorylates and forms 3′-dAMP which may be the immediate precursor of cordycepin. A total of 13 proteins, including 4 DEPs (one up-regulated ADK, one down-regulated adk, one down-regulated surE and one down-regulated RRM1) participated in the biosynthetic pathway of cordycepin (Supplementary information 5: Table S5). In addition, several types of purine nucleosides including inosine, adenosine, guanosine and xanthosine could be synthesized from IMP. The possible metabolic pathways of purine nucleosides and the corresponding proteins were shown in Supplementary information 5 (Figure S5, Table S6).
Metabolic pathway of monosaccharides
Polysaccharides from O. sinensis zjut were made up of mannose, galactose and glucose, and its total content was 3.848% (Supplementary information 1: Table S1). According to the metabolic pathways, there are 18 proteins involved in the biosynthetic pathway of D-mannose (Supplementary information 6: Figure S6, Table S7), 14 proteins involved in D-glucose biosynthesis (Supplementary information 6: Figure S7, Table S8), and 14 proteins involved in D-galactose biosynthesis (Supplementary information 6: Figure S8, Table S9). The metabolic pathways, including corresponding proteins and mRNA involved in these pathways, were shown in Supplementary information 6.
Amino acid metabolic pathways
Sixteen types of amino acids in O. sinensis zjut were detected by amino acid analyzer (Supplementary information 7: Table S10). The amino acid with the highest content in O. sinensis zjut was histidine (40.93 mg/g), followed by arginine (36.14 mg/g). The amino acids with the lowest contents was tyrosine (6.23 mg/g), followed by isoleucine (7.25 mg/g), and phenylalanine (8.08 mg/g). The proteins involved in histidine, arginine, phenylalanine and tyrosine biosynthesis were successfully obtained (Supplementary information 8: Table S11, Table S12, Table S13). The synthetic pathways of histidine (Supplementary information 8: Figure S9), arginine (Supplementary information 8: Figure S10), phenylalanine and tyrosine (Supplementary information 8: Figure S11) were successfully constructed, respectively.
Noteworthy proteins in O. sinensis zjut
Based on the analysis of KEGG database, the metabolic pathway of glycolysis/gluconeogenesis and citrate cycle (Supplementary information 10: Figure S12, Figure S13,) were constructed, and the proteins involved in these energy pathways were obtained (Supplementary information 9 and 10: Table S14, Table S15, Table S16). Moreover, 9 proteins (two superoxide dismutase, three catalase, a thioredoxin reductase and three peroxiredoxin) with antioxidant activity were obtained (Supplementary information 11: Table S17), which might play an important role in oxidation resistance.
Recently, the genetic information of Cordyceps was gradual enrichment [25, 26], and the productions of active compounds in H. sinensis by submerged fermentation were promoted and satisfying on the basis of genetic study [15, 16]. With the studies focusing on genetic information, there were little proteomics reports for artificially-cultivated O. sinensis. As we know, gene expression regulation was a complex multilevel process, and the correlation between mRNA and protein abundances was only approximately 27 to 40% . In fact, genes were the carriers of genetic information [32, 33], and proteins were the executors of the physiological functions and direct manifestations of life activity . Therefore, it is urgent and important to enhance proteome analysis of artificially-cultivated O. sinensis.
iTRAQ, as a reliable quantitative approach, has widely used in the field of crop proteomics, allowing simultaneous identification and quantification of proteins from multiple samples with high coverage [35,36,37]. In this study, proteomes of O. sinensis zjut at different culture periods were investigated by iTRAQ, and this systematically large-scale study of proteins provided a novel way to study the gene functions and expression difference directly at protein levels in submerged fermentation for O. sinensis. Subsequent studies, including function annotation and metabolic pathways construction, were not only the complement and verification to the genetic analysis with genome and transcriptome but also the necessary to fundamental research. As an analytic system, the genome, transcriptome and proteome sequencing analysis could provide a scientific basis for medicinal mechanism to carry out great exploration, application and improvement of submerged fermentation for O. sinensis.
Active ingredients of submerged fermentation for O. sinensis, such as mannitol, cordycepin and polysaccharide, are used to treat weakness after sickness, lung and kidney-associated diseases, and sexual dysfunction [6, 38, 39]. Recently, the metabolic pathways of active ingredients were constructed at the transcriptional level in H. sinensis . In this study, the metabolic pathways were also predicted directly at protein level, which played significant roles in further study. For metabolic pathway of mannitol (Fig. 4), mannitol-1-phosphate phosphatase that catalyzes D-mannitol-1-P to D-mannitol was not annotated in proteome, and it was more confidence to the result of the transcriptome study for other phosphatases replaced the function of mannitol-1-phosphate phosphatase or mannitol-1-phosphate phosphatase sequence of O. sinensis zjut, which was un-annotated to protein database since its low homology with currently reported mannitol-1-phosphate phosphatases from other organisms . With the proteome analysis of O. sinensis zjut, the new metabolic pathway of cordycepin (Fig. 5) was different with the conjecture in the transcriptome study. Following the biosynthesis of adenosine, the proteomics result indicated SurE (E3.1.35) (5′-nucleotidase) catalyzes the conversion of 3′-dAMP (3′-adenine deoxynucleotide) to cordycepin rather than N-glycosylation lyase catalyzes the glycosylation exchange between adenosine and cordycepose to generate cordycepin in genetic studies [26, 40]. The constructed metabolic pathways could be the guidance in metabolic engineering study of this organism in future.
Based on the proteomics at different culture periods, DEPs were obtained to be analyzed. As shown in Fig. 4 and Table S4, two up-regulated HKs may boost the conversion of α-D-glucose to α-D-glucose-6P as well as the conversion of D-fructose to β-D-fructose-6P, and the up-regulated FBP may boost the conversion of β-D-fructose-1,6P2 to β-D-fructose-6P, leading to the accumulation of β-D-fructose-6P. Two down-regulated manA, one down-regulated PFK and one up-regulated mtlD may convert more β-D-fructose-6P to D-mannit-1P, which might indicate the growth rate of D-mannitol was rise with the culture periods of the O. sinensis zjut. There were three down-regulated proteins that might be the key enzymes in the metabolic pathway of cordycepin (Fig. 5 and Table S5). Comparing with the 3rd day, adk expression level was down-regulated 0.86-fold on the 6th day and 0.77-fold on the 9th day, RRM1 expression level was down-regulated 0.87-fold on the 6th day and 0.8-fold on the 9th day, and surE expression level was down-regulated 0.65-fold on the 6th day and 0.57-fold on the 9th day. These down-regulated proteins catalyze AMP to ADP, ADP to 3′-dADP, 3′-dADP to 3′-dAMP, 3′-dAMP to cordycepin, respectively. The growth rate of cordycepin was decreased with the culture periods of the O. sinensis zjut, since adenosine or IMP could be an intermediate for other compounds such as inosine, xanthosine and guanosine . As a secondary metabolic pathway, the low content of cordycepin on 6th day (3rd day 0.327 mg/g, 6th day 0.140 mg/g, and 9th day 0.457 mg/g) might be caused by periodical release. The DEPs involved in mannitol and cordycepin biosynthesis play important functions, and studies focus on these enzymes would pave a theoretical foundation to regulate and promote the active ingredients productions of submerged fermentation by O. sinensis.
Previous study reported that sixteen to eighteen types of amino acids existed in O. sinensis . In this study, sixteen types of amino acids were obtained in O. sinensis zjut including 7 essential amino acids and 13 pharmaceutical amino acids. As shown in Supplementary information 7: Table S10, the contents of most individual amino acids were discrepant between O. sinensis zjut and natural O. sinensis. The two principal amino acids and their levels in O. sinensis zjut were 13.36% histidine and 11.80% arginine, compared with 3.08% histidine and 5.54% arginine in natural O. sinensis. Furthermore, the content of valine, which was the most different amino acid, was 7.98% in O. sinensis zjut, compared with 1.03% in O. sinensis. The variation in O. sinensis zjut and natural O. sinensis might be associated with geography, climate and nutrition differences.
In this study, we found some noteworthy pathways and proteins in O. sinensis zjut. Fatty acids, which might protect organism from cold, are the major energy sources for O. sinensis zjut [43, 44]. As the major sources for metabolic intermediates [45, 46], two energy pathways (glycolysis/gluconeogenesis pathway and the citrate cycle) might also provide energy for O. sinensis zjut and protect it from cold. Recently, superoxide dismutase was extracted from O. sinensis , to study for powerful anti-inflammatory and pharmacological activities [30, 48]. Two superoxide dismutases, which were responsible for converting superoxide to O2 and H2O2, were found in O. sinensis zjut. Among the proteinases, some special functional proteins have been found, including two secreted aspartic proteinases, which are important hydrolases and potent antifungal agents [49, 50]. Moreover, the proteins with special function could be cloned and then expressed in Escherichia coli to study the enzymatic property for exploring its industry application [51, 52].
In this study, iTRAQ technology together with two-dimensional liquid chromatography tandem mass spectrometry were first used to relatively quantify the proteins of O. sinensis zjut at different culture periods. In total, 4005 proteins were obtained and further annotated, and the metabolic pathways of D-mannitol, cordycepin, monosaccharides, two high content (histidine and arginine) and two low content amino acids (phenylalanine and tyrosine), and fatty acid were constructed at the protein level. The important proteins involved in metabolic pathways of active ingredients were found by comparative proteomics of O. sinensis zjut at different culture periods. In addition, two energy pathways (glycolysis/gluconeogenesis pathway and the citrate cycle) were described. Proteins with antioxidant and antifungal activity were also obtained, which could be cloned, expressed and studied in the future. This work provided scientific basis to study protein profile and comparison of protein expression levels of submerged fermentation for O. sinensis, and it will be helpful for metabolic engineering works to active ingredients for exploration, application and improvement of this strain.
Strains, medium and culture conditions
Wild O. sinensis samples were collected from the slopes of 4000–4300 m above sea level on the Qinghai-Tibet plateau in Yushu city, Qinghai Province. The isolation and identification processes as well as preservation information of O. sinensis zjut were referred to a previous study . The medium (pH 7.0) for O. sinensis zjut cultivation contained glucose (1.5%), corn powder (1%), silkworm chrysalis meal (2%), MgSO4 (0.05%), and KH2PO4 (0.05%). The strain was grown on the medium using 200 L submerged stirred fermenter at 16 °C. The single colonies, growth curve and active compounds contents of O. sinensis zjut were performed (Supplementary information 1: Figure S1, Table S1). The samples for proteome analysis were harvested on the 3rd, 6th and 9th days.
iTRAQ quantitative proteome
The quantitative proteome was carried out by iTRAQ according to the standard experimental protocol. The details of the iTRAQ procedure were shown in Supplementary information 2.
The software Iquant was used for quantitatively analyzing the labeled peptides with isobaric tags . To assess the confidence of peptides, the PSMs (propensity matching analysis) were pre-filtered at a PSM-level FDR of 1%. Identified peptide sequences were assembled into a set of confident proteins with the parsimony principle.
Enrichment analysis of the DEPs was conducted according to the information from Kyoto Encyclopedia of Genes and Genome (KEGG) pathway and Gene Ontology (GO) databases, respectively, using the following formula:
N is the number of all proteins that can be annotated to GO or KEGG, n is the number of DEPs in N, M is the number of proteins that can be annotated to a certain GO term or KEGG pathway, and m is the number of DEPs in M.
Blast2GO tool (23) was used to categorize the identified proteins according to their functions, biological process and cellular component. The metabolic pathway analysis was conducted according to the KEGG Pathway Database (http://www.genome.jp/kegg/pathway.html).
Unless otherwise noted, all the experiments in this study were performed in triplicate. An analysis of variance (ANOVA) was performed using the SAS program version 8.1 (SAS Institute Inc., Cary, NC, USA). The least significant difference (LSD) was computed at p < 0.05. All the figures in this study were processed using the Origin software version 8.0 (OriginLab Corp., Northampton, MA, USA).
Availability of data and materials
The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium (http://proteomecentral.proteomexchange.org/) via the iProX partner repository  with the dataset identifier PXD016899.
Availability of data and materials
All data generated or analyzed during this study are included in this article and its supplementary information files.
Isobaric tags for relative and absolute quantification
False discovery rate
Differentially expressed proteins
Propensity matching analysis
Sharma S. Trade of Cordyceps sinensis from high altitudes of the Indian Himalaya: conservation and biotechnological priorities. Curr Sci India. 2004;86:1614–9.
Winkler D. Yartsa Gunbu (Cordyceps sinensis) and the fungal commodification of Tibet's rural economy. Econ Bot. 2008;62:291–305.
Zhu JS, Halpern GM, Jones K. The scientific rediscovery of an ancient Chinese herbal medicine: Cordyceps sinensis: part I. J Altern Complement Med. 1998;4:289–303.
Singh M, Tulsawani R, Koganti P, Chauhan A, Manickam M, Misra K. Cordyceps sinensis increases hypoxia tolerance by inducing heme oxygenase-1 and metallothionein via Nrf2 activation in human lung epithelial cells. Biomed Res Int. 2013;2013:569206.
Shashidhar MG, Giridhar P, Sankar KU, Manohar B. Bioactive principles from Cordyceps sinensis: a potent food supplement - a review. J Funct Foods. 2013;5:1013–30.
Buenz EJ, Bauer BA, Osmundson TW, Motley TJ. The traditional Chinese medicine Cordyceps sinensis and its effects on apoptotic homeostasis. J Ethnopharmacol. 2005;96:19–29.
Zhang Y, Liu X, Wang M. Cloning, expression, and characterization of two novel cuticle-degrading serine proteases from the entomopathogenic fungus Cordyceps sinensis. Res Microbiol. 2008;159:462–9.
Chen YQ, Wang N, Qu LH, Li TH, Zhang WM. Determination of the anamorph of Cordyceps sinensis inferred from the analysis of the ribosomal DNA internal transcribed spacers and 5.8S rDNA. Biochem Syst Ecol. 2001;29:597–607.
Shang HS, Chou J, Wu MF, Shih YL, Yeh MY, Chung HY, et al. Evaluation of Hirsutella sinensis mycelium on food safety and anti-hepatoma activity in an animal model. In Vivo. 2014;28:811–7.
Zhu ZY, Liu XC, Fang XN, Sun HQ, Yang XY, Zhang YM. Structural characterization and anti-tumor activity of polysaccharide produced by Hirsutella sinensis. Int J Biol Macromol. 2016;82:959–66.
Wu LY, Wu MF, Lu HF, Liu CH, Lee CH, Chen YL, et al. Evaluation of Hirsutella sinensis mycelium for antifatigue effect. In Vivo. 2015;29:263–7.
He L, Ji P, Cheng J, Wang Y, Qian H, Li W, et al. Structural characterization and immunostimulatory activity of a novel protein-bound polysaccharide produced by Hirsutella sinensis. Food Chem. 2013;141:946–53.
Meng L, Sun S, Li R, Shen Z, Wang P, Jiang X. Antioxidant activity of polysaccharides produced by Hirsutella sp. and relation with their chemical characteristics. Carbohydr Polym. 2015;117:452–7.
Wu TR, Lin CS, Chang CJ, Lin TL, Martel J, Ko YF, et al. Gut commensal Parabacteroides goldsteinii plays a predominant role in the anti-obesity effects of polysaccharides isolated from Hirsutella sinensis. Gut. 2019;68:248–62.
Lin S, Liu ZQ, Baker PJ, Yi M, Wu H, Xu F, et al. Enhancement of cordyceps polysaccharide production via biosynthetic pathway analysis in Hirsutella sinensis. Int J Biol Macromol. 2016;92:872–80.
Liu ZQ, Zhang B, Lin S, Baker PJ, Chen MS, Xue YP, et al. Enhancement of nucleoside production in Hirsutella sinensis based on biosynthetic pathway analysis. Biomed Res Int. 2017;2017:1–11.
Lin S, Liu ZQ, Xue YP, Baker PJ, Wu H, Xu F, et al. Biosynthetic pathway analysis for improving the cordycepin and cordycepic acid production in Hirsutella sinensis. Appl Biochem Biotech. 2016;179:633–49.
Pandey A, Mann M. Proteomics to study genes and genomes. Nature. 2000;405:837.
Kumarathasan P, Mohottalage S, Goegan P, Vincent R. An optimized protein in-gel digest method for reliable proteome characterization by MALDI-TOF-MS analysis. Anal Biochem. 2005;346:85–9.
Liang Y, Chen H, Tang M, Shen S. Proteome analysis of an ectomycorrhizal fungus Boletus edulis under salt shock. Mycol Res. 2007;111:939–46.
Gan CS, Chong PK, Pham TK, Wright PC. Technical, experimental, and biological variations in isobaric tags for relative and absolute quantitation (iTRAQ). J Proteome Res. 2007;6:821–7.
Ross PL, Huang YN, Marchese JN, Williamson B, Parker K, Hattan S, et al. Multiplexed protein quantitation in Saccharomyces cerevisiae using amine-reactive isobaric tagging reagents. Mol Cell Proteomics. 2004;3:1154–69.
Creaney J, Dick IM, Leon JS, Robinson BW. A proteomic analysis of the malignant mesothelioma secretome using iTRAQ. Cancer Genom Proteom. 2017;14:103–17.
Wu WW, Wang G, Baek SJ, Shen RF. Comparative study of three proteomic quantitative methods, DIGE, cICAT, and iTRAQ, using 2D gel- or LC-MALDI TOF/TOF. J Proteome Res. 2006;5:651–8.
Hu X, Zhang YJ, Xiao GH, Zheng P, Xia YL, Zhang XY, et al. Genome survey uncovers the secrets of sex and lifestyle in caterpillar fungus. Chin Sci Bull. 2013;58:2846–54.
Liu ZQ, Lin S, Baker PJ, Wu LF, Wang XR, Wu H, et al. Transcriptome sequencing and analysis of the entomopathogenic fungus Hirsutella sinensis isolated from Ophiocordyceps sinensis. BMC Genomics. 2015;16:106.
Zhang X, Liu Q, Zhou W, Li P, Alolga RN, Qi LW, et al. A comparative proteomic characterization and nutritional assessment of naturally- and artificially-cultivated Cordyceps sinensis. J Proteome. 2018;181:24–35.
Dong YZ, Zhang LJ, Wu ZM, Gao L, Yao YS, Tan NZ, et al. Altered proteomic polymorphisms in the caterpillar body and stroma of natural Cordyceps sinensis during maturation. PLoS One. 2014;9:e109083.
Zhang S, Lai X, Li B, Wu C, Wang S, Chen X, et al. Application of differential proteomic analysis to authenticate Ophiocordyceps sinensis. Curr Microbiol. 2016;72:337–43.
Yin YL, Yu GJ, Chen YJ, Jiang S, Wang M, Jin YX, et al. Genome-wide transcriptome and proteome analysis on different developmental stages of Cordyceps militaris. PLoS One. 2012;7:e51853.
Maier T, Güell M, Serrano L. Correlation of mRNA and protein in complex biological samples. FEBS Lett. 2009;583:3966–73.
Xi JQ, Tang D, Guo Y. Data gene:genetic information carrier of data. Comput Eng. 2006;32:7–9.
Cheng F, Xiang C, Zhang XJ, Liu ZQ, Zheng YG. ReToAd: simple method for the rapid replacement of promoters to improve protein production. Biotechnol Lett. 2018;40:957–64.
Zhang Y, Xu P, Li W, Tao Y. [Advances in actinobacterial proteomics]. Chinese. J Biotechnol. 2014;30:1044.
Guo X, Xu J, Cui X, Chen H, Qi H. iTRAQ-based protein profiling and fruit quality changes at different development stages of Oriental Melon. BMC Plant Biol. 2017;17:28.
Lin XQ, Liang SL, Han SY, Zheng SP, Ye YR, Lin Y. Quantitative iTRAQ LC-MS/MS proteomics reveals the cellular response to heterologous protein overexpression and the regulation of HAC1 in Pichia pastoris. J Proteome. 2013;91:58–72.
Wang Y, Ma X, Zhang X, He X, Li H, Cui D, et al. ITRAQ-based proteomic analysis of the metabolic mechanisms behind lipid accumulation and degradation during peanut seed development and postgermination. J Proteome Res. 2016;15:4277–89.
Zhou X, Gong Z, Su Y, Lin J, Tang K. Cordyceps fungi: natural products, pharmacological functions and developmental products. J Pharm Pharmacol. 2009;61:279–91.
Li C, Li Z, Fan M, Cheng W, Long Y, Ding T, et al. The composition of Hirsutella sinensis, anamorph of Cordyceps sinensis. J Food Compos Anal. 2006;19:800–5.
Xia YL, Luo FF, Shang YF, Chen PL, Lu YZ, Wang CS. Fungal cordycepin biosynthesis is coupled with the production of the safeguard molecule pentostatin. Cell Chem Biol. 2017;24:1479–89.
Chassy BM, Suhadolnik RJ. Nucleoside antibiotics. IV. Metabolic fate of adenosine and cordycepin by Cordyceps militaris during cordycepin biosynthesis. Biochim Biophys Acta. 1969;182:307–15.
Hsu TH, Shiao LH, Hsieh CY, Chang DM. A comparison of the chemical composition and bioactive ingredients of the Chinese medicinal mushroom DongChongXiaCao, its counterfeit and mimic, and fermented mycelium of Cordyceps sinensis. Food Chem. 2002;78:463–9.
Guo LX, Xu XM, Wu CF, Lin L, Zou SC, Luan TG, et al. Fatty acid composition of lipids in wild Cordyceps sinensis from major habitats in China. Biomed Prev Nutr. 2012;2:42–50.
Guo LX, Xu XM, Liang FR, Yuan JP, Peng J, Wu CF, et al. Morphological observations and fatty acid composition of indoor-cultivated Cordyceps sinensis at a high-altitude laboratory on Sejila mountain, Tibet. Plos One. 2015;10:e0126095.
Schlessinger A, Sun NN, Colas C, Pajor AM. Determinants of substrate and cation transport in the human Na+/dicarboxylate cotransporter NaDC3. J Biol Chem. 2014;289:16998–7008.
Kruger NJ, von Schaewen A. The oxidative pentose phosphate pathway: structure and organisation. Curr Opin Plant Biol. 2003;6:236–46.
Kong X, Wang H, Xiao Z, Chu J, Yin X, Yang X. Experimental study on the antiaging effect of Cordyceps sinensis cultivated by artificial fermentation. J Capital Univ Med Sci. 1997;18:15–7.
Sung-Jun J, Tae-Hee L, Dae-Hoon C, Yeong-Hwan H. Effect of light conditions on production of cordycepin of Cordyceps militaris. Korean J Microbiology. 2005;41:236–8.
Zhang SP, Feng H, Li XY, Jin YS, Dong W. Genome research profile of two Cordyceps sinensis cDNA libraries. Chin Sci Bull. 2010;55:1403–11.
Liu ZQ, Dong SC, Yin HH, Xue YP, Tang XL, Zhang XJ, et al. Enzymatic synthesis of an ezetimibe intermediate using carbonyl reductase coupled with glucose dehydrogenase in an aqueous-organic solvent system. Bioresour Technol. 2017;229:26–32.
Liu ZQ, Wu L, Zheng L, Wang WZ, Zhang XJ, Jin LQ, et al. Biosynthesis of tert-butyl (3R,5S)-6-chloro-3,5-dihydroxyhexanoate by carbonyl reductase from Rhodosporidium toruloides in mono and biphasic media. Bioresour Technol. 2018;249:161–7.
Liu ZQ, Lu MM, Zhang XH, Cheng F, Xu JM, Xue YP, et al. Significant improvement of the nitrilase activity by semi-rational protein engineering and its application in the production of iminodiacetic acid. Int J Biol Macromol. 2018;116:563–71.
Wen B, Zhou R, Feng Q, Wang QH, Wang J, Liu SQ. IQuant: an automated pipeline for quantitative proteomics based upon isobaric tags. Proteomics. 2014;14:2280–5.
Ma J, Chen T, Wu S, Yang C, Bai M, Shu K, et al. iProX: an integrated proteome resource. Nucleic Acids Res. 2019;47:D1211–7.
We would like to thank National High Technology Research and Development Program of China and the Key Scientific and Technology Programs of Zhejiang Province.
This work has been partly supported by the National High Technology Research and Development Program of China (No.2012AA021701) and the Key Scientific and Technology Programs of Zhejiang Province (No.2012C03005–1).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
: The single colonies, growth curve and active compounds contents in O. sinensis zjut; Supporting information 2 The details of iTRAQ quantitative proteomics; Supporting information 3: The identified results of the samples and DEPs in O. sinensis zjut; Supporting information 4: GO enrichment analysis of DEPs in O. sinensis zjut; Supporting information 5: The proteins involved in metabolic pathway of D-mannital, cordycepin and purine nucleotides in O. sinensis zjut, and purine nucleotides metabolic pathway; Supporting information 6: The metabolic pathways of D-mannose, D-galactose and D-glucose as well as the proteins involved in these pathways; Supporting information 7: The contents of sixteen types of amino acids in O. sinensis zjut; Supporting information 8: The metabolic pathways of the histidine, arginine, phenylalanine and tyrosine as well as the proteins involved in these pathways; Supporting information 9: The proteins involved in the biosynthetic pathway of fatty acids; Supporting information 10: The metabolic pathways of glycolysis/gluconeogenesis and citrate cycle as well as the proteins involved in these pathways; Supporting information 11: The proteins with anti-oxidant activity in O. sinensis zjut.
About this article
Cite this article
Zhang, B., Li, B., Men, XH. et al. Proteome sequencing and analysis of Ophiocordyceps sinensis at different culture periods. BMC Genomics 21, 886 (2020). https://doi.org/10.1186/s12864-020-07298-z
- O. sinensis zjut
- Differentially expressed proteins
- Active ingredients