==== Front BMC Genomics BMC Genomics BMC Genomics 1471-2164 BioMed Central London 7298 10.1186/s12864-020-07298-z Research Article Proteome sequencing and analysis of Ophiocordyceps sinensis at different culture periods Zhang Bo zhangbo0305@zjut.edu.cn 1 Li Bo 279672107@qq.com 1 Men Xiao-Hui 1142097886@qq.com 1 Xu Zhe-Wen 570398426@qq.com 1 Wu Hui wuhui@eastprotech.com 23 Qin Xiang-Tian qinxt@eastchinapharm.com 23 Xu Feng xufeng@hdpharm.com 23 Teng Yi tengyi198777@126.com 23 Yuan Shui-Jin shuijinyuan@eastchinapharm.com 23 Jin Li-Qun jlq@ziut.edu.cn 1 http://orcid.org/0000-0003-3259-6796Liu Zhi-Qiang microliu@zjut.edu.cn 1 Zheng Yu-Guo zhengyg@zjut.edu.cn 1 1 grid.469325.f0000 0004 1761 325XKey Laboratory of Bioorganic Synthesis of Zhejiang Province, College of Biotechnology and Bioengineering, Zhejiang University of Technology, Hangzhou, 310014 China 2 HuaDong Medicine (Hangzhou) Bailing Biological Technology Co., Ltd, Hangzhou, 311220 China 3 East China Pharmaceutical Group Limited Co., Ltd, Hangzhou, 311000 China 11 12 2020 11 12 2020 2020 21 88626 10 2019 2 12 2020 © The Author(s) 2020Open AccessThis article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.Background 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. Results 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. Conclusion 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. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-020-07298-z. Keywords O. sinensis zjutProteomeDifferentially expressed proteinsActive ingredientsNational High Technology Research and Development Program of China2012AA021701Liu Zhi-Qiang Key Scientific and Technology Programs of Zhejiang Province2012C03005-1Liu Zhi-Qiang issue-copyright-statement© The Author(s) 2020 ==== Body Background Ophiocordyceps sinensis is an important traditional Chinese medicine and healthy food in China [1–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 [8], 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 [11], immunomodulatory activity [12], antioxidant activity [13] and anti-obesity effects [14]. Furthermore, the productions of active ingredients in H. sinensis by submerged fermentation were promoted and satisfying on the basis of genetic study [15–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 [18]. 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–24]. Genome sequence of O. sinensis has been reported, which revealed the pathogenic mechanism during life cycle [25]. Then, transcriptomes of H. sinensis at diffident culture periods have been sequenced and analyzed to describe metabolic pathway and infection mechanism [26]. 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 [27]. 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 [28]. In addition, there were some non-quantitative studies about the proteomics of O. sinensis [27–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 [18]. 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. Results 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. Fig. 1 Barplot of the Gene Ontology analysis. The bar chart shows the distribution of corresponding GO terms. Different colors represent different GO categories Fig. 2 Pie chart of the COG Analysis. All identified proteins are classified into 24 clusters of orthologous groups (COG) categories 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). Fig. 3 Statistics of pathway enrichment of differentially expressed proteins in the comparison groups of 6d-VS-3d, 9d-VS-3d, and 9d-VS-6d 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). Fig. 4 The metabolic pathway of D-mannitol in O. sinensis zjut. The biosynthesis pathway of D-mannitol was constructed based on the annotation of the O. sinensis zjut proteome. The precursor of D-mannitol is D-mannitol-1P which is formed from β-D-fructose-6P. The dotted line represented the protein did not annotated to the O. sinensis zjut proteome. The DEPs (HK, manA, PFK, mtlD, FBP) was shown with bold fonts. HK: Hexokinase; manA: Mannose-6-phosphate isomerase; FBP: Fructose-1,6-bisphosphatase I; pfkA: 6-phosphofructokinase; GPI: Glucose-6-phosphate isomerase; E5.1.3.15: Glucose-6-phosphate 1-epimerase; FBA: Fructose-bisphosphate aldolase; galM: Aldose 1-epimerase; pgm: Phosphoglucomutase; mtlD: Mannitol-1-phosphate 5-dehydrogenase; PFK: 6-phosphofructo-2-kinase; PFKFB2: Fructose-2,6-biphosphatase 2; PFKFB4: Fructose-2,6-biphosphatase 4 Fig. 5 The metabolic pathway of cordycepin in O. sinensis zjut. The biosynthesis pathway of cordycepin was constructed based on the annotation of the O. sinensis zjut proteome. The dotted line represented the possible protein involved in the conversion of 3′-dAMP to cordycepin. The DEPs (ADK, adk, RRM1, surE) was shown with bold fonts. purH: IMP cyclohydrolase; AMDP: AMP deaminase; adk: Adenylate kinase; purA: Adenylosuccinate synthase; purB: Adenylosuccinate lyase; ADK: Adenosine kinase; surE: 5′-nucleotidase; E3.1.3.5: 5′-nucleotidase; RRM1: Ribonucleoside-diphosphate reductase subunit M1; RRM2: Ribonucleoside-diphosphate reductase subunit M2 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. Discussion 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% [31]. 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 [34]. 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–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 [26]. 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 [26]. 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 [41]. 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 [42]. 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 [47], 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]. Conclusion 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. Methods 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 [26]. 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. Protein quantification The software Iquant was used for quantitatively analyzing the labeled peptides with isobaric tags [53]. 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. Function analysis 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: P=1−∑i=0m−1MiN−Mn−iNn 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). Statistical analysis 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 [54] with the dataset identifier PXD016899. Supplementary Information Additional file 1: Supporting information 1: 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. Abbreviations iTRAQIsobaric tags for relative and absolute quantification FDRFalse discovery rate DEPsDifferentially expressed proteins PSMsPropensity matching analysis Publisher’s Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Acknowledgements We would like to thank National High Technology Research and Development Program of China and the Key Scientific and Technology Programs of Zhejiang Province. Authors’ contributions ZQL BZ LQJ and YGZ conceived and designed the experiments, BZ, BL, FX, YT, SJY, ZWX and XHM carried out all data analyses, and wrote the manuscript. BZ, HW and XTQ revised the manuscript. All authors read and approved the final manuscript. Funding 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). Availability of data and materials All data generated or analyzed during this study are included in this article and its supplementary information files. Ethics approval and consent to participate Not applicable. Consent for publication Not applicable. Competing interests The authors declare that they have no competing interests. ==== Refs References 1. Sharma S Trade of Cordyceps sinensis from high altitudes of the Indian Himalaya: conservation and biotechnological priorities Curr Sci India 2004 86 1614 1619 2. Winkler D Yartsa Gunbu (Cordyceps sinensis ) and the fungal commodification of Tibet's rural economy Econ Bot 2008 62 291 305 10.1007/s12231-008-9038-3 3. 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 10.1089/acm.1998.4.3-289 9764768 4. 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 24063008 5. 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 1030 10.1016/j.jff.2013.04.018 32288795 6. 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 10.1016/j.jep.2004.09.029 15588646 7. 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 469 10.1016/j.resmic.2008.04.004 18555668 8. 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 10.1016/S0305-1978(00)00100-9 11336809 9. Shang HS Chou J Wu MF Shih YL Yeh MY Chung HY Evaluation of Hirsutella sinensis mycelium on food safety and anti-hepatoma activity in an animal model In Vivo 2014 28 811 817 25189893 10. 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 966 10.1016/j.ijbiomac.2015.10.075 26517958 11. Wu LY Wu MF Lu HF Liu CH Lee CH Chen YL Evaluation of Hirsutella sinensis mycelium for antifatigue effect In Vivo. 2015 29 263 267 25792655 12. He L Ji P Cheng J Wang Y Qian H Li W Structural characterization and immunostimulatory activity of a novel protein-bound polysaccharide produced by Hirsutella sinensis Food Chem 2013 141 946 953 10.1016/j.foodchem.2013.04.053 23790872 13. 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 457 10.1016/j.carbpol.2014.09.076 25498658 14. Wu TR Lin CS Chang CJ Lin TL Martel J Ko YF Gut commensal Parabacteroides goldsteinii plays a predominant role in the anti-obesity effects of polysaccharides isolated from Hirsutella sinensis Gut. 2019 68 248 262 10.1136/gutjnl-2017-315458 30007918 15. Lin S Liu ZQ Baker PJ Yi M Wu H Xu F Enhancement of cordyceps polysaccharide production via biosynthetic pathway analysis in Hirsutella sinensis Int J Biol Macromol 2016 92 872 880 10.1016/j.ijbiomac.2016.08.002 27496606 16. Liu ZQ Zhang B Lin S Baker PJ Chen MS Xue YP Enhancement of nucleoside production in Hirsutella sinensis based on biosynthetic pathway analysis Biomed Res Int 2017 2017 1 11 17. Lin S Liu ZQ Xue YP Baker PJ Wu H Xu F Biosynthetic pathway analysis for improving the cordycepin and cordycepic acid production in Hirsutella sinensis Appl Biochem Biotech 2016 179 633 649 10.1007/s12010-016-2020-0 18. Pandey A Mann M Proteomics to study genes and genomes Nature. 2000 405 837 10.1038/35015709 10866210 19. 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 89 10.1016/j.ab.2005.06.004 16168381 20. Liang Y Chen H Tang M Shen S Proteome analysis of an ectomycorrhizal fungus Boletus edulis under salt shock Mycol Res 2007 111 939 946 10.1016/j.mycres.2007.06.005 17716885 21. 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 827 10.1021/pr060474i 17269738 22. Ross PL Huang YN Marchese JN Williamson B Parker K Hattan S Multiplexed protein quantitation in Saccharomyces cerevisiae using amine-reactive isobaric tagging reagents Mol Cell Proteomics 2004 3 1154 1169 10.1074/mcp.M400129-MCP200 15385600 23. Creaney J Dick IM Leon JS Robinson BW A proteomic analysis of the malignant mesothelioma secretome using iTRAQ Cancer Genom Proteom 2017 14 103 117 10.21873/cgp.20023 24. 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 658 10.1021/pr050405o 16512681 25. Hu X Zhang YJ Xiao GH Zheng P Xia YL Zhang XY Genome survey uncovers the secrets of sex and lifestyle in caterpillar fungus Chin Sci Bull 2013 58 2846 2854 10.1007/s11434-013-5929-5 26. Liu ZQ Lin S Baker PJ Wu LF Wang XR Wu H Transcriptome sequencing and analysis of the entomopathogenic fungus Hirsutella sinensis isolated from Ophiocordyceps sinensis BMC Genomics 2015 16 106 10.1186/s12864-015-1269-y 25765329 27. Zhang X Liu Q Zhou W Li P Alolga RN Qi LW A comparative proteomic characterization and nutritional assessment of naturally- and artificially-cultivated Cordyceps sinensis J Proteome 2018 181 24 35 10.1016/j.jprot.2018.03.029 28. Dong YZ Zhang LJ Wu ZM Gao L Yao YS Tan NZ Altered proteomic polymorphisms in the caterpillar body and stroma of natural Cordyceps sinensis during maturation PLoS One 2014 9 e109083 10.1371/journal.pone.0109083 25310818 29. Zhang S Lai X Li B Wu C Wang S Chen X Application of differential proteomic analysis to authenticate Ophiocordyceps sinensis Curr Microbiol 2016 72 337 343 10.1007/s00284-015-0916-5 26660081 30. 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. 31. Maier T Güell M Serrano L Correlation of mRNA and protein in complex biological samples FEBS Lett 2009 583 3966 3973 10.1016/j.febslet.2009.10.036 19850042 32. Xi JQ Tang D Guo Y Data gene:genetic information carrier of data Comput Eng 2006 32 7 9 33. 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 964 10.1007/s10529-018-2541-4 29611067 34. Zhang Y Xu P Li W Tao Y [Advances in actinobacterial proteomics]. Chinese J Biotechnol 2014 30 1044 10.1002/biot.201400072 35. 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 10.1186/s12870-017-0977-7 28129739 36. 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 10.1016/j.jprot.2013.06.031 37. Wang Y Ma X Zhang X He X Li H Cui D 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 4289 10.1021/acs.jproteome.6b00345 27669742 38. 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 291 10.1211/jpp.61.03.0002 19222900 39. Li C Li Z Fan M Cheng W Long Y Ding T The composition of Hirsutella sinensis , anamorph of Cordyceps sinensis J Food Compos Anal 2006 19 800 805 10.1016/j.jfca.2006.04.007 40. 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 1489 10.1016/j.chembiol.2017.09.001 29056419 41. 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 315 10.1016/0005-2787(69)90181-6 4978604 42. 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 469 10.1016/S0308-8146(02)00158-9 43. Guo LX Xu XM Wu CF Lin L Zou SC Luan TG Fatty acid composition of lipids in wild Cordyceps sinensis from major habitats in China Biomed Prev Nutr 2012 2 42 50 10.1016/j.bionut.2011.07.002 44. Guo LX Xu XM Liang FR Yuan JP Peng J Wu CF 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 10.1371/journal.pone.0126095 25938484 45. 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 17008 10.1074/jbc.M114.554790 24808185 46. Kruger NJ von Schaewen A The oxidative pentose phosphate pathway: structure and organisation Curr Opin Plant Biol 2003 6 236 246 10.1016/S1369-5266(03)00039-6 12753973 47. 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 17 48. 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 238 49. 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 1411 10.1007/s11434-010-0113-7 50. Liu ZQ Dong SC Yin HH Xue YP Tang XL Zhang XJ 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 10.1016/j.biortech.2016.12.098 28092733 51. Liu ZQ Wu L Zheng L Wang WZ Zhang XJ Jin LQ 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 167 10.1016/j.biortech.2017.09.204 29040850 52. Liu ZQ Lu MM Zhang XH Cheng F Xu JM Xue YP 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 571 10.1016/j.ijbiomac.2018.05.045 29753012 53. 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 2285 10.1002/pmic.201300361 25069810 54. Ma J Chen T Wu S Yang C Bai M Shu K iProX: an integrated proteome resource Nucleic Acids Res 2019 47 D1211 D1217 10.1093/nar/gky869 30252093