Skip to main content
  • Research article
  • Open access
  • Published:

Screening for proteins related to the biosynthesis of hispidin and its derivatives in Phellinus igniarius using iTRAQ proteomic analysis



Hispidin (HIP) and its derivatives, a class of natural fungal metabolites, possess complex chemical structures with extensive pharmacological activities. Phellinus igniarius, the most common source of HIP, can be used as both medicine and food. However, the biosynthetic pathway of HIP in P. igniarius remains unclear and we have a limited understanding of the regulatory mechanisms related to HIP. In this work, we sought to illustrate a biosynthesis system for hispidin and its derivatives at the protein level.


We found that tricetolatone (TL) is a key biosynthetic precursor in the biosynthetic pathway of hispidin and that its addition led to increased production of hispidin and various hispidin derivatives. Based on the changes in the concentrations of precursors and intermediates, key timepoints in the biosynthetic process were identified. We used isobaric tags for relative and absolute quantification (iTRAQ) to study dynamic changes of related proteins in vitro. The 270 differentially expressed proteins were determined by GO enrichment analysis to be primarily related to energy metabolism, oxidative phosphorylation, and environmental stress responses after TL supplementation. The differentially expressed proteins were related to ATP synthase, NAD binding protein, oxidoreductase, and other elements associated with electron transfer and dehydrogenation reactions during the biosynthesis of hispidin and its derivatives. Multiple reaction monitoring (MRM) technology was used to selectively verify the iTRAQ results, leading us to screen 11 proteins that were predicted to be related to the biosynthesis pathways.


These findings help to clarify the molecular mechanism of biosynthesis of hispidin and its derivatives and may serve as a foundation for future strategies to identify new hispidin derivatives.


Edible and medicinal fungi are a common source of nutrients and medicines in Asian nations and are well recognized as an important source of biologically active compounds. Phellinus igniarius (DC.Ex Fr.) Quel, a wild macrofungi, contains many bioactive compounds with properties such as antibacterial, antioxidative, antitumor, and antimutagenic activities, and has been widely used in China, Japan, and Korea for many years [1,2,3]. P. igniarius is rich in secondary metabolites with biological activities, such as polysaccharides, polyphenols, flavonoids, and other chemical species [4, 5]. Phelligridins are a class of chemical compounds from P. igniarius that possess radical-scavenging activity and a broad range of pharmacological properties, and includes chemicals such as Phelligridin D and Phelligridin G [6, 7]. Based on structural analyses, it has been postulated that hispidin is a key intermediate of phelligridins [8]. Hispidin (6-(3,4-dihydroxystyryl)-4-hydroxy-2-pyrone, HIP), is an important polyphenol found in P. igniarius. Hispidin is widely distributed in edible mushrooms, with many reports of potential pharmacologic value in the literature [9, 10]. Prior reports have suggested that HIP has activities such as anti-oxidation [11], anti-inflammatory [12], anti-cancer [13], and antiallergic [14] effects, among others. HIP is also an important luciferin precursor in some bioluminescent fungi and can be used for bioluminescent analysis of substances [15, 16]. A dose of 3 mg/g hispidin-enriched fungus mycelia has a very low level of toxicity, supporting its safety for human consumption [17]. Therefore, HIP has attracted significant attention in the fields of chemistry, pharmacy, and microbiology [18,19,20,21].

Due to the limited and unpredictable supply of wild P. igniarius, artificial fermentation has become a topic of increasing interest in recent years. Significant research has focused on optimizing fermentation conditions with the goal of increasing the yield of HIP [17]. With the development of sequencing technology, gene expression control at the transcriptional level is one method to determine composition and facilitate industrial production. Knowledge of the biosynthetic pathways of metabolites provides an important theoretical basis for such engineering approaches. To date, although the biosynthetic pathway of HIP and its derivatives in P. igniarius has been the focus of extensive research, the exact mechanisms and pathways remain incompletely understood. A previous study proposed the biosynthesis of HIP as shown in Fig. 1, where HIP is formed by the condensation of 4-hydroxy-6-methyl-2-pyrone (TL) and one molecule of 3,4-dihydroxybenzoyl-SCoA (or 3,4-dihydroxybenzaldehyde) in the P. igniarius fruit body [22]. This hypothesis has not been verified, and the biosynthetic pathway of HIP has not been categorically confirmed.

Fig. 1
figure 1

The proposed biosynthesis of HIPs

This work investigated whether TL is transformed into HIP in the fermentation broth of P. igniarius. Moreover, the extent of protein expression by the mycelium of P. igniarius was measured using iTRAQ technology following TL feeding. Proteins related to HIP synthesis were identified. The analysis explored how each protein participated following TL precursor feeding, and provides deeper insights into the mechanisms of synthesis of HIP at the proteomic level.


Metabolic analysis method determination

An HPLC method was optimized to distinguish between HIP and other derivatives in order to analyze the ethyl acetate extracts. As seen in Fig. 2, the peak at 19.7 min was attributed to HIP, 23.4 min was phelligridin D, and 9.6 min was TL. To determine that the peaks at 19.7 min and 23.4 min were indeed HIP and phelligridin D, the chemical structures of the two purified compounds A and B were assessed by NMR spectra and HRESIMS (see Fig. S1–9).

Fig. 2
figure 2

The metabolites after cultured for 7 days after TL feeding

Compound A was obtained as a light-yellow powder. Its quasi-molecular ion peak, observed at m/z 245.0453 (cal. 245.0450) [M - H] by HRESIMS and NMR data, showed the molecular formula as C13H10O5. The 1H-NMR spectrum of A in CH3OD (see Fig. S3) displayed resonances in pairs attributable to two analogous compounds with a content ratio of 3:2. This difference was determined to result from the cis and trans double bonds and some changes in chemical shift. The resonances were ascribed to two ABX spin coupling systems of 1,3,4-trisubstituted phenyl moiety at δH 6.95 (dd, J = 2.0, 8.0 Hz), 6.77 (d, J = 8.0 Hz), 7.03 (d, J = 2.0 Hz), and δH 6.80 (dd, J = 2.0, 8.0 Hz), 6.73 (d, J = 8.0 Hz), 6.98 (d, J = 2.0 Hz), a trans-disubstituted and a cis-disubstituted double bond at δH 7.31 (d, J = 16.0 Hz), 6.60 (d, J = 16.0 Hz), and 6.72 (d, J = 12.8 Hz), 5.98 (d, J = 12.8 Hz). The 13C-NMR spectrum of A gave two sets of 12 sp2 carbon resonances. Compared with the peak position of standards in HPLC and the data from HRESIMS, the compound was speculated to be a combination of trans-hispidin and cis-hispidin (see Fig. S1), which was also confirmed by 2D NMR. Analysis of gHSQC spectroscopic data of A furnished assignments of the proton-bearing carbon and corresponding proton resonances in the NMR spectra (Table 1). Although the signals of H3 and C3 disappeared in NMR, the visible long-range heteronuclear correlations in the HMBC of trans-hispidin for H-5/C-3 (~ 90 ppm, see Fig. S7) and the data in the literature [11, 23] supported the prior findings. Regardless of whether it is trans or cis, hispidin is a crucial intermediate of phelligridins [22].

Table 1 NMR spectroscopic data of compound Aa)

Compound B was obtained as a yellow powder. Its quasi-molecular ion peak, found at m/z 379.0466 (cal. 379.0454) [M - H] by HRESIMS, gave a molecular formula of C20H12O8. Compound B was identified as phelligridin D by comparing with standards and data from 1H-NMR spectra (see Fig. S9) [22].

The dynamic changes in concentration of TL and HIP

Compared to the negative control condition that was not given TL, the contents of HIP were increased five-fold three days after the addition of TL, and the TL content gradually decreased (Fig. 3), suggesting that TL participates in the biosynthesis of HIP. Interesting, the content of TL and hispidin, from 8 to 16 h after TL feeding, ran contrary to the general trends observed over the longer time period, suggesting that the strain was in the logarithmic phase of growth through this period. Hispidin content rose sharply over the course of 16–64 h, from < 5 mg/L at 16 h to 44.38 mg/L at 64 h. Interestingly, TL rapidly decreased during the first 128 h after addition, wherein the concentration of TL declined from 423.61 mg/L at 1 h to 50.44 mg/L at 128 h), suggesting that TL may be consumed to synthesize hispidin. After 64 h, the concentration of HIP began to decline, presumably due to the synthesis of other derivatives.

Fig. 3
figure 3

The dynamic changes in concentration of TL and HIP

iTRAQ protein analysis

Herein, the expression profiles of proteins that were produced by P. igniarius that was fed with either TL-supplemented or normal media were compared by iTRAQ analysis. A total of 339,985 spectrums were filtered by 1% FDR, which led to the identification of 5630 peptides and 1880 proteins (Fig. 4). These results indicated that iTRAQ had higher sensitivity and provided more comprehensive information than other techniques in the analysis of these fungal proteins. The molecular weights of the identified proteins were widely distributed, mostly falling between 10 and 70 kDa. There were 552 proteins with molecular weights between 20 and 30 kDa, 794 proteins between 30 and 40 kDa, 543 between 40 and 50 kDa, and 710 over 100 kDa (Fig. 5a and b). Protein coverage analysis showed that the numerical value from 0 to 100% gradually declined. There were 966 proteins with more than 11 specific peptide segments. The number of specific spectra is shown in Fig. 5c. The number distribution of specific peptide segments is shown in Fig. 5d, and suggests that the results from this analysis were generally reliable.

Fig. 4
figure 4

Overview of protein identification

Fig. 5
figure 5

Information for protein identification. a: protein quantity distribution; b: identified protein coverage; c: number of specific peptide segments; d: specific peptide number

Analysis of the difference expression of proteins (DEPs)

Screening of DEPs was carried out in order to determine which proteins were expressed with a fold change > 1.5 and a Q value < 0.05. We identified 615 DEPs in the TLPF32h group, including 385 up-regulated proteins and 230 down-regulated proteins. There were 907 DEPs in the control group compared with TLPF128h, including 249 up-regulated proteins and 458 down-regulated proteins. As seen in Fig. 6, up-regulation of proteins was more common after the addition of TL, but the number of down-regulated proteins was higher in the TLPF128h group compared to the TLPF32h group.

Fig. 6
figure 6

Three sets of differentially expressed proteins (DEPs). a GO analysis for TLPF32h-VS-control. b GO analysis for TLPF128h-VS-control. c GO analysis for TLPF128h-VS- TLPF32 h

Gene ontology (GO) is an important element of bioinformatics research that aims to unify analyses of gene expression and gene attributes of all species. GO annotation analysis showed that 297 DEPs in the TLPF32h group and control group were annotated as belonging to 35 functional groups (Fig. 7). Briefly, the DEPs were found to be mainly involved in cellular processes (38.05%) and metabolic processes (34.34%). Among the cell components, DEPs were mainly concentrated in cells (39.73%) and cell parts (39.06%). In terms of molecular function, DEPs were mainly involved in catalytic activity (55.56%) and binding activity (46.46%).

Fig. 7
figure 7

Gene ontology (GO) analysis between groups. A: the KEGG pathway analysis for TLPF32h-VS-control. B: the KEGG pathway analysis for TLPF128h- VS-control. C: the KEGG pathway analysis for TLPF128h-VS- TLPF32 h

These results indicated that TL-reactive proteins may be primarily involved in stress responses, chemical stimulation responses, primary metabolic processes, cell responses, and other related functions. We might expect that when P. igniarius experiences changing environmental stimuli, the defense system would immediately respond, increasing metabolic activity, producing defensive substances, and enhancing the activity of various enzymes to promote defense.

In KEGG pathway enrichment analysis, a total of total 82 DEPs were matched with the KEGG pathway database using Blast_v2.2.26 software. Compared with the control group (Fig. 8), 265 DEPs were labeled as 52 KEGG pathways, and metabolic pathways (ko01100, 77 DEPs) were the primary pathways that were enriched. Indeed, the DEPs were clearly concentrated in metabolic pathways and oxidative phosphorylation, and most of these were up-regulated. Thus, it appears that TL supplementation accelerates P. igniarius metabolism and produces oxidative stress. The other DEPs that were identified were mainly related to secondary metabolites and ribosomal protein biosynthesis in the TLPF128h group. Ribosomal proteins represent a common type of RNA-binding protein, which mainly participate in protein translation, post-translational modification, or protein folding. Because we found that the ribosomal protein pathway was down-regulated, we hypothesize that the synthesis of new protein would be slowed. Compared with the TLPF32h group, metabolic processes and ribosomal proteins were the differentially regulated main pathways in the TLPF128h group. Most of the differentially expressed proteins (DEPs) in the 128 h group were down-regulated relative to 32 h, suggesting that the metabolic rate decreased over time.

Fig. 8
figure 8

The KEGG pathway analysis between groups

Selective verification of iTRAQ data by mass spectrometry multi-reaction monitoring technology (MRM)

In light of the above protein analysis results, 15 candidate DEPs related to the TL reaction were selected to establish an MRM method. Of the 15 target proteins, 11 had MS/MS spectra and unique peptide(s). Therefore, those 11 DEPs were detected and quantified by MRM (Table 2). The results showed that the expression levels of DEP and iTRAQ were almost in direct agreement with each other.

Table 2 Quantity analysis of the DEPs


Genomic mining and functional determination are commonly used to study the mechanisms of fungal biosynthesis and molecular composition. The synthetic genes and regulatory mechanisms of some fungal products have been elucidated in detail, such as monacolin K [24] and griseofulvin. Hispidin and its derivatives are active constituents of edible fungi [25]. Proteomic analysis is an important element of mechanistic studies because of the uncertainties involved with genome sequencing and biosynthetic pathways of compounds. This method was utilized in early studies of the biosynthesis of aflatoxin, and the biosynthesis mechanism was confirmed by feeding the intermediate versicolorin A and using LC-MS/MS to analyze the resulting proteins [26].

In this study, the proteins involved in the biosynthesis of secondary metabolites were identified by iTRAQ method for the first time, and a large number of DEPs were detected in P. igniarius after supplementation with the precursor molecule TL. In addition, there were 270 DEPs that were separated into groups including energy metabolism, oxidative phosphorylation, and environmental stress responses. These proteins may provide a new perspective for the study of the biosynthesis of hispidin and its derivatives. Comparative analyses showed that DEPs found after TL feeding were mainly involved in cell processes, metabolic processes, located in cells or cell parts, catalytic activities, and binding functions.

Comparing the proteomic changes after TL supplementation

After comparing the dynamic changes in concentration of HIP, it appeared that TL supplementation led to significant improvement of HIP production. KEGG pathway enrichment analysis of three groups indicated that DEPs of TLPF32h vs control group were most closely associated with the metabolic pathway (Fig. 8). Almost all proteins had high expression levels that participated in oxidative phosphorylation. The proton pump F-ATPase is an ATP synthase, which plays an important role in energy transduction in most species. Further, the F-ATPase helps maintain the pH of the cytoplasm by secreting protons, and F-ATPase inhibitors reduce the acid resistance of bacteria [27]. In this study, a large number of F-ATPase were up-regulated after TL feeding. The enhancement of ATP synthase, consistent with increased expression of proteins related to oxidative phosphorylation, could increase the ATP supply to the biosynthesis of HIP. Consistently, proteomic analysis showed that NADH dehydrogenase related energy production were also up-regulated in P. igniarius. NADH dehydrogenase is the major element of respiratory chain complex I, Electrons mainly move from NADH to the mitochondrial respiratory chain. NADH dehydrogenase is a small protein (< 10 kDa) and relatively simple in structure [28]. Its sub-structures and functional model have been well established. In this study, several complex I subunits were up-regulated to produce more protons, which promote ATP synthesis. That is to say, in P. igniarius, increased synthesis of ATP may play a central role in the differential expression seen after supplementation with TL. This ATP-related enzyme was down-regulated on KEGG pathway enrichment analysis of TLPF128h vs the TLPF32h group. The inexistence of riboflavin, thiamine, glutathione, cysteine and methionine, and selenocompound metabolism at TLPF128h vs TLPF32h group showed mycelium metabolism basically returned to normal. It is worth noting that proteins related to ribosomes were significantly down-regulated in TLPF128h. This suggest that synthesis of protein was decreased at that time.

Environmental stress may lead to the accumulation of HIPs

The regulatory mechanisms of macrofungi are largely unclear because of the complexity of the physiological process and genomes. In past years, the development of proteomics has provided insights and a systematic understanding of changes in filamentous fungi and yeast at the molecular level in response to environmental changes, such as the response of Flammulina velutipes to cold and light stress [29]. Among the 270 DEPs, some were stress response proteins produced by organisms in response to environmental changes. For example, cysteine proteases respond to biological and abiotic stressors in the external environment. They not only participate in the regulation of transcription factors, but also process proteins [30]. Glutathione S-transferases, which are mainly found in the cytoplasm, are multifunctional supergene family proteins. Fungal glutathione S-transferases have a variety of structures and functions, regulate the growth and development of organisms, and protect the enzyme system from environmental stress [31]. Heat shock proteins (HSPs) were perhaps the most significantly changed DEPs. In emergency circumstances, HSPs help to maintain proper protein folding, promote the transmembrane transport of protein molecules, and promote the recovery or degradation of damaged proteins. In order to promote the hydrolysis of ATP, HSPs activate protein kinase C, thereby reducing the damage caused by oxygen free radicals. Recent research has shown that pyruvate accumulation is a protective mechanism against stress conditions in fungi [32]. In this study, our data showed that pyruvate dehydrogenase and dihydrolipoyllysine-residue acetyltransferase content were increased, probably because of pyruvate accumulation (Fig. 9). These proteins play an important role in the defense response and were up-regulated in TLPF32h. Our previous studies showed that there was trace hispidin and hispidin derivatives in the fermentation broth of P. igniarius, found at concentrations too low to be separated and extracted. However, the compounds can be extracted from the fruiting bodies of wild fungi. The dramatic increase in stress proteins indicated that TL feeding caused a stress reaction in P. igniarius mycelium. The content of HIP tended to increase over the same timeframe. Therefore, we speculate that the metabolic pathway of hispidin in P. igniarius may be related to the accumulation of secondary metabolites in stressful environments.

Fig. 9
figure 9

Representative metabolic pathway maps of differentially expressed proteins involved in pyruvate metabolic in KEGG. The enzyme codes were noted as follows: pyruvate dehydrogenase (EC, dihydrolipoyllysine-residue acetyltransferase (EC, dihydrolipoamide dehydrogenase (EC Red for upregulated proteins

Biosynthesis pathway of HIPs

The pathway for biosynthesis of HIP was inferred to be related to the cinnamate pathway for the metabolism of phenylalanine as far back as 1973 [33]. Recent research has revealed that HIP is a precursor of luciferin 3-hydroxyhispidin [34]. It has been found that the luminescence of luminous fungal fruiting body extracts is enhanced by the addition of hispidin biosynthetic components, namely caffeic acid, ATP, and malonyl-CoA [35]. Another theory, that HIP is condensated by TL and 3,4-dihydroxybenzaldehyde based on structural analysis of compounds isolated from the fungus fruiting body [22]. The biogenesis of many HIP derivatives has led to speculation that it is generated by the condensation of HIP, a process that could be catalyzed by peroxidase [36]. The mechanism for synthesis of HIP and other complex derivatives, however, is not specific.

TL is a common, small molecule, polyketone precursor compound that participates in the biosynthesis of HIPs. Expression of the protein GME4094_g, a polyketide synthase, was down-regulated in TLPF32h and not up-regulated in TLPF128h (Table 3). In the synthesis of HIP, dehydrogenase is necessary. Dehydrogenases GME1855_g, GME4242_g, GME4316_g, and GME9582_g were essential in the HIP synthesis pathway and were increased during the HIP production period. A series of enzymes related to the biosynthesis of HIPs is listed in Table 3. In addition, some oxygenases (Table 4) appeared to play a role in the synthesis of HIP’s structural derivatives (such as phelligridin D, Fig. 10). In this process, phenolic compounds are oxidized to benzoquinones that stimulate free radical transfer and then form a variety of compounds with novel chemical skeletons. With further elucidation of the secondary metabolites and regulation of its biosynthetic pathway, it will be increasingly possible to control the production of hispidin via genetic engineering.

Table 3 The DEPs related to the biosynthesis of hispidin
Table 4 The DEPs related to the biosynthesis of phelligridins
Fig. 10
figure 10

The proposed biosynthesis mechanism of phelligridins


In this study, a new synthetic pathway of hispidin was analyzed and validated by the precursor feeding method. The DEPs of P. igniarius after TL feeding were analyzed using the iTRAQ technique, and DEPs involved in the biosynthesis of hispidin and its derivatives were subsequently screened. The synthesis of hispidin and its derivatives appears to be related to a large number of oxidoreductases and dehydrogenases as well as to stress response systems. According to the results of proteomic analyses and previous data, it appears that TL affects the biosynthesis of hispidin through several routes, including by acting on reactive oxygen species scavenging, signaling pathways, secondary metabolite synthesis, and other reactions. This work sought to provide further mechanistic clarification for the biosynthetic pathways important for the production of hispidin and its derivatives through analysis of whole gene sequences, functional gene annotations, and prediction of secondary metabolite synthesis gene clustering, combined with analysis of DEPs from iTRAQ results. Further research should focus on the biosynthetic gene cluster and synthetase of hispidin.



P. igniarius (CGMCC 5.95) was purchased from the China Microbial Preservation Center (Beijing, China) and was cultured in modified Martin medium (MMM), including 5 g/L peptone, 1 g/L K2HPO3, 2 g/L yeast extract, 20 g/L glucose, 20 g/L agar, pH 6.2–6.5. The culture medium was sterilized by autoclave at 121 °C for 30 min. The strain was stored at 4 °C.

HPLC detection of metabolites by precursor feeding

P. igniarius was cultured in MMM medium for seven days at 28 °C. An agar patch (1 cm × 1 cm) was added into a 250 mL conical bottle containing 150 mL liquid medium (the above MMM without agar). The culture conditions were 28 °C, 180 rpm. Then, 15 days later, the TL, which was dissolved in double distilled water, was added into the fermentation liquor at a final concentration of ~ 0.1 mg/mL. The co-culture was then allowed to grow for another 9 days under the same conditions. Samples (2.0 mL) were taken at 1, 3, 5, 7, and 9 days after the addition of TL, and were extracted by EtOAc (3 × 2.0 mL). The organic layer was removed by volatilization and the residue was dissolved in 2 mL methanol for HPLC detection. HPLC was conducted using Agilent 1220 Infinity II equipment with Agilent Eclipse XDB-C18 (4.6 mm × 250 mm, 5 μm) as the separation column. Two-phase gradient elution methods were used in this work including 0.2% formic acid (an aqueous solution) as phase A, and MeOH as phase B. The MeOH fraction in phase B was gradually increased from 40% (V/V) to 80% (V/V) over 25 min. The spectra were recorded by a DAD detector at 380 nm. The flow velocity was set to 1 mL/min, and the injection volume was 10 μL.

Separation and identification of metabolites

A 5000 mL scale-up of P. igniarius was cultured according to the above-mentioned fermentation conditions. After the 0.1 mg/mL TL was added, the co-culture was grown for another 3 days. The filtrate of the fermentation broth was concentrated to 500 mL and partitioned with EtOAc (3 × 500 mL). The EtOAc extract was then evaporated under reduced pressure to yield 290 mg of residue, which was subjected to Sephadex LH-20 column chromatography, eluted with petrol-CHCl3-MeOH (5:4:1) to produce five fractions (A-E). Fraction B was purified through reverse-phase preparative HPLC using a mobile phase of MeOH-H2O (45:55) to afford compound A (25.5 mg). Fraction D was separated by preparative RP-HPLC using MeOH-H2O (62:38) to afford compound B (7.0 mg). The structures of the target compounds were identified by 1D and 2D NMR as well as HRESIMS analysis. 1D- and 2D-NMR spectra were obtained at 400 MHz for 1H and 100 MHz for 13C, respectively, on Bluker 400 MHz spectrometers in methanol-d4 with solvent peaks used as references. HRESIMS data were measured using an Agilent 1290 Infinity II Accurate Mass Q-ToF LC/MS spectrometer.

Determination of HIP yield by external standard method

The concentrations of HIP and TL were determined by HPLC with an external standard at the wavelength of 254 nm, and a 20–80% methanol gradient was used for elution. The concentration points for the HIP standard curve were 1.0, 0.5, 0.25, 0.125, 0.0625, and 0.03125 mg/mL, and the concentrations for the TL standard curve were 8, 4, 2, 1, 0.5, and 0.25 mg/mL. The linear standard equation was obtained by using the linear least square method. As in the previous methods description, P. igniarius was cultured in MMM liquid medium for 15 days and then fed TL at a final concentration of 0.1 mg/mL. Fermentation broth samples were taken at 0 h, 1 h, 2 h, 4 h, 8 h, 32 h, 64 h, and 128 h after the addition of TL into the fermentation system, and then the samples were extracted with ethyl acetate and analyzed by HPLC.

Protein extraction and enzymatic hydrolysis

A 5 mm magnetic bead and 25 μL of lysis buffer were added to 5 mg mycelium samples. The final concentrations of PMSF and EDTA were 1 mM and 2 mM, respectively. The mycelium samples were allowed to rest for 5 min after the eddy oscillation, and the final concentration of DTT was 10 mM. The supernatant was collected by centrifugation at 4 °C for 20 min after 2 min oscillation with a tissue abrasive apparatus. The supernatant was treated with 10 mM DTT for 1 h in a water bath at 56 °C. After coming back down to room temperature, IAM was added at a final concentration of 55 mM, and then the sample was kept in the dark for 45 min. The supernatant was precooled with 5× volume of acetone, precipitated at − 20 °C for 2 h, and was centrifuged at 9000×g, 4 °C for 20 min. The addition of acetone was repeated three times, followed each time by centrifugation and discarding the supernatant until the supernatant was colorless. Then, 25 μL of lysis buffer was added into the precipitation mixture. After 5 min of ultrasonication in an ice bath, the supernatant was centrifuged at 9000×g, 4 °C for 20 min.

iTRAQ marker and peptide segment separation

First, 100 μg protein solution was extracted from each sample, and the trypsin enzyme was added in the ratio of protein:enzyme of 40:1. The enzyme was hydrolyzed for 4 h at 37 °C, at which point trypsin was again added (the same amount), and the mixture was continuously hydrolyzed for another 8 h. The resulting enzymatically generated peptides were desalted using a Strata X column and then vacuum-dried.

Eight groups of iTRAQ labeling reagents (113, 114, 115, 116, 117, 118, 119, 121) were selected. Then, 50 μL isopropanol were added to each tube at 25 °C, and the mixture was centrifuged after swirl oscillation. The supernatants were transferred to another clean sample tube along with the enzymatic peptides in 0.5 M TEAB. The different peptide fragments were labeled with their respective peptide’s labels.

The above samples were separated by liquid phase chromatography with a LC-20AB system (Shimadzu, Japan), and the separation column used was a Gemini C18 column (5 μm, 4.6 × 250 mm). The above dried samples were re-dissolved with phase A (5% ACN, pH 9.8), and the flow rate gradient was set as: 5% mobile phase B (95% ACN, pH 9.8) for 10 min, 5–35% mobile phase B for 40 min, 35–95% mobile phase B for 3 min, and 5% mobile phase B for 10 min, at rate of 1.0 × 10− 3 L/min. The detection wavelength used was 214 nm. The fractions were collected once per minute, and the sample components were merged by chromatographic elution peaks to yield a total of 20 fractions. Those 20 fractions were frozen and concentrated each to the same volume of 1.0 mL.

LC-MS/MS analysis

After centrifugation for 10 min, the supernatant samples were separated using LC-20 AD nanoflow liquid chromatography. A trap column was used to concentrate protein and remove salts. It was connected with a self-assembled C18 column (75 μm inner diameter, 3.6 μm column diameter, 15 cm column length) in series. The flow rate was set to 300 nL/min. Separation was carried out using the following gradients:

(1) 0–8 min, 5% mobile phase B (98% acetonitrile aqueous solution containing 0.1% formic acid (FA));

(2) 8–43 min, 8–35% mobile phase B gradient;

(3) 43–48 min, 35–60% mobile phase B gradient;

(4) 48–50 min, 60–80% mobile phase B gradient;

(5) 50–55 min, equivalent 80% mobile phase B;

(6) 55–65 min, equivalent 5% mobile phase B.

After ionization from a nano ESI source, the peptide segments were analyzed by high resolution liquid chromatography-mass spectrometry (HR-LC-MS) TripleTOF 5600. Using Proteome Discoverer, a Thermo Scientific tool, the original mass spectrometry file was converted into a MGF format file containing the information of secondary mass spectrometry (MS/MS) spectra, in which “BEGIN IONS” and “END IONS” were the starting and ending positions of each spectrum. The UniProt protein database was used to identify the proteins. The MGF file and protein database were searched to obtain the final protein identification results using the Mascot 2.3.02 identification software.

iTRAQ data analysis

iTRAQ data were quantified by IQuant software. The spectra and the list of peptide segments were filtered by 1% FDR (PSM-level FDR < 0.01) to identify significant results. According to the parsimony principle, the peptide segments were used to assemble proteins and produce a series of proteomes. In order to control the false positive rate of proteins, the process was again filtered at the protein level with 1% FDR (Protein-level FDR < 0.01). Quant’s workflow included the following steps: protein filtering, purity correction of report group labels, normalization of quantitative values, complementation of missing values, calculation of quantitative values of proteins, statistical analysis, and display of final results.

Bioinformatic analysis

In this study, three sets of experiments including Control, TLPF32h, TLPF128h, and three sets of comparisons of TLPF32h /control, TLPF128h/control, TLPF128h/TLPF32h were conducted. The differentially expressed proteins (DEPs) were screened on the basis of a fold change > 1.5 and P < 0.05.

In the GO enrichment analysis of DEPs, the GO entries with significant enrichment were identified by hypergeometric test, compared with all identified proteins as background. The principle of Pathway enrichment analysis was similar. The hypergeometric test formula was as follows:

$$ \mathrm{P}=1\hbox{-} \sum \limits_{\mathrm{i}=0}^{\mathrm{m}-1}\frac{\left({}_{\mathrm{i}}^{\mathrm{M}}\right)\left({}_{n-i}^{N-M}\right)}{\left({}_n^N\right)} $$

N: The number of GO entries matched in all identified proteins; n: the number of GO entries matched in DEPs; M: one of the GO entries matched in all identified proteins; m: one of the GO entries matched in DEPs. DEPs were determined to be significantly enriched in the GO entry if the P value of the hypergeometric test was less than 0.05.

Validation of MRM technology

According to the results of iTRAQ data, differentially expressed proteins were screened and selected for MRM validation. Each sample was treated with 100 μg protein solution. Trypsin enzyme was added at a ratio of protein:enzyme of 40:1, with a total amount of 2.5 μg of trypsin enzyme, and enzymatic hydrolysis was carried out at 37 °C for 4 h. Trypsin was added one more time at the same ratio, and the enzymatic hydrolysis was continued at 37 °C for another 8 h. The enzymatic peptides were desalted using a Strata X column and then vacuum-dried.

The samples were scanned with a HPLC-TripleTOF 5600 mass spectrometer, and the resulting data were searched by Mascot v2.3 using the Fungi protein database (11,243 sequences) added to the internal standard peptide sequence. The DAT files were imported into Skyline software to establish the atlas library (credibility > 0.95). Skyline software was used to select target peptide segments under the following conditions:

(1) The peptide segments had matched second-order ions.

(2) The length of the peptide segments were between 5 and 40 amino acids.

(3) The peptide segment was the only one for the target protein.

(4) The cysteine in the peptide segment was modified to carbamidomethyl.

(5) There was no variable modification of the peptide segment.

(6) No methionine was found in the peptide segment.

(7) There was no missing cut of the peptide segment.

Transition selection set was determined as follows:

(1) The fragment ions were B and Y ions;

(2) The parent ion charges were 2,3,4;

(3) The charges of fragment ions were 1,2;

(4) Debris ion charge-mass ratio < 1250 (four-stage rod scanning range);

(5) The number of transitions was 6.

Validation of MRM method

Skyline set up the MRM method and output it to the QTRAP 5500 mass spectrometer for MRM scanning verification. The success of MRM mass spectrometry for target proteins depended on the validation results, which must conform to:

(1) There were co-elution peaks in transitions with different peptide segments.

(2) The chromatographic peak area intensity of transitions with different peptide segments was correlated with the fragmentation intensity of spectral data.

(3) The retention time of the peptide segments in MR and full-spectrum scans was good.

MRM mass spectrometry detection

The peptide fragments were separated by liquid phase and entered into the QTRAP 5500 tandem mass spectrometer. The ion source was Nanospray Illsource. In data acquisition, the instrument parameters are set as follows: spray voltage was 2400 V, spray gas was 23. With MRM scanning mode, the resolution of Q1 and Q3 was set to Unit mode.

Data analysis

Each transition signal of the target protein was normalized to the signal of beta-galactosidase. After normalized intensity, a linear mixed model integrated with the MS stats tool was used to quantify the target protein in the sample. This model gave the ratio of protein in the comparison group and the adjusted p-value. The corrected P value reflected the false positive rate of the original statistical test (Benjamin and Hochberg). If the difference in final target protein concentration was at least 1.5-fold with P < 0.05 (false positive < 0.05), the protein was considered to be significantly differentially expressed.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.





Adenosine triphosphate


Diode array detector


Difference expression of proteins




Ethylene diamine tetraacetic acid


False discovery rate


Gene ontology




Heteronuclear multiple bond correlation


High performance liquid chromatography


High resolution electrospray ionization mass spectroscopy


Heat shock protein


Heteronuclear single quantum correlation




Kyoto encyclopedia of genes and genomes


Modified Martin medium


Multiple reaction monitoring


Nuclear Magnetic Resonance


Phenylmethanesulfonyl fluoride


Tetraethylammonium bromide


Tricetolatone precursor feeding


  1. Wang Y, Shang XY, Wang SJ, Mo SY, Li S, Yang YC, Ye F, Shi JG, He L. Structures, biogenesis, and biological activities of pyrano [4,3-c]isochromen-4-one derivatives from the fungus Phellinus igniarius. J Nat Prod. 2007;70(2):296–9.

    Article  Google Scholar 

  2. Zhang H, Ma H, Liu W, Pei J, Wang Z, Zhou H, Yan J. Ultrasound enhanced production and antioxidant activity of polysaccharides from mycelial fermentation of Phellinus igniarius. Carbohydr Polym. 2014;113:380–7.

    Article  CAS  Google Scholar 

  3. Konno S, Chu K, Feuer N, Phillips J, Choudhury M. Potent anticancer effects of bioactive mushroom extracts (Phellinus linteus) on a variety of human Cancer cells. J Clin Med Res. 2015;7(2):76–82.

    Article  Google Scholar 

  4. Gao W, Wang W, Sun W, Wang M, Zhang N, Yu S. Antitumor and immunomodulating activities of six Phellinus igniarius polysaccharides of different origins. Exp Ther Med. 2017;14(5):4627–32.

    CAS  PubMed  PubMed Central  Google Scholar 

  5. Lee MS, Hwang BS, Lee IK, Seo GS, Yun BS. Chemical constituents of the culture broth of Phellinus linteus and their antioxidant activity. Mycobiology. 2015;43(1):43–8.

    Article  Google Scholar 

  6. Kim JE, Takanche JS, Yun BS, Yi HK. Anti-inflammatory character of Phelligridin D modulates periodontal regeneration in lipopolysaccharide-induced human periodontal ligament cells. J Periodontal Res. 2018;53(5):816–24.

    Article  CAS  Google Scholar 

  7. Kim JY, Kim DW, Hwang BS, Woo EE, Lee YJ, Jeong KW, Lee IK, Yun BS. Neuraminidase inhibitors from the fruiting body of Phellinus igniarius. Mycobiology. 2016;44(2):117–20.

    Article  Google Scholar 

  8. Wang Y, Mo SY, Wang SJ, Li S, Yang YC, Shi JG. A unique highly oxygenated pyrano [4,3-c][2] benzopyran-1,6-dione derivative with antioxidant and cytotoxic activities from the fungus Phellinus igniarius. Org Lett. 2005;7(9):1675–8.

    Article  CAS  Google Scholar 

  9. Zan LF, Bao HY, Bau T, Li Y. A new antioxidant Pyrano [4,3-c][2] benzopyran-1,6-dione derivative from the medicinal mushroom Fomitiporia ellipsoidea. Nat Prod Commun. 2015;10(2):315–6.

    PubMed  Google Scholar 

  10. Hwang BS, Lee IK, Choi HJ, Yun BS. Anti-influenza activities of polyphenols from the medicinal mushroom Phellinus baumii. Bioorg Med Chem Lett. 2015;25(16):3256–60.

    Article  CAS  Google Scholar 

  11. Park IH, Chung SK, Lee KB, Yoo YC, Kim SK, Kim GS, Song KS. An antioxidant hispidin from the mycelial cultures of Phellinus linteus. Arch Pharm Res. 2004;27(6):615–8.

    Article  CAS  Google Scholar 

  12. Shao HJ, Jeong JB, Kim KJ, Lee SH. Anti-inflammatory activity of mushroom-derived hispidin through blocking of NF-B activation. J Sci Food Agr. 2015;95(12):2482–6.

    Article  CAS  Google Scholar 

  13. Lim JH, Lee YM, Park SR, Kim DH, Lim BO. Anticancer activity of Hispidin via reactive oxygen species-mediated apoptosis in Colon Cancer cells. Anticancer Res. 2014;34(8):4087–93.

    CAS  PubMed  Google Scholar 

  14. Tamrakar S, Fukami K, Parajuli GP, Shimizu K. Antiallergic activity of the wild mushrooms of Nepal and the pure compound Hispidin. J Med Food. 2019;22(2):225–7.

    Article  CAS  Google Scholar 

  15. Purtov KV, Petushkov VN, Rodionova NS, Gitelson JI. Why does the bioluminescent fungus Armillaria mellea have luminous mycelium but nonluminous fruiting body? Dokl Biochem Biophys. 2017;474(1):217–9.

    Article  CAS  Google Scholar 

  16. Puzyr AP, Medvedeva SE, Burov AE, Zernov YP, Bondar VS. Detection of Hispidin by a luminescent system from Basidiomycete Armillaria borealis. Dokl Biochem Biophys. 2018;480(1):173–6.

    Article  CAS  Google Scholar 

  17. Li IC, Chen CC, Sheu S, Huang IH, Chen CC. Optimized production and safety evaluation of hispidin-enriched Sanghuangporus sanghuang mycelia. Food Sci Nutr. 2020;8(4):1864–73.

    Article  CAS  Google Scholar 

  18. Grundemann C, Arnhold M, Meier S, Backer C, Garcia-Kaufer M, Grunewald F, Steinborn C, Klemd AM, Wille R, Huber R, et al. Effects of Inonotus hispidus extracts and compounds on human Immunocompetent cells. Planta Med. 2016;82(15):1359–67.

    Article  Google Scholar 

  19. Kim HG, Yoon DH, Lee WH, Han SK, Shrestha B, Kim CH, Lim MH, Chang W, Lim S, Choi S, et al. Phellinus linteus inhibits inflammatory mediators by suppressing redox-based NF-kappaB and MAPKs activation in lipopolysaccharide-induced RAW 264.7 macrophage. J Ethnopharmacol. 2007;114(3):307–15.

    Article  CAS  Google Scholar 

  20. Lv LX, Zhou ZX, Zhou Z, Zhang LJ, Yan R, Zhao Z, Yang LY, Bian XY, Jiang HY, Li YD, et al. Hispidin induces autophagic and necrotic death in SGC-7901 gastric cancer cells through lysosomal membrane permeabilization by inhibiting tubulin polymerization. Oncotarget. 2017;8(16):26992–7006.

    Article  Google Scholar 

  21. Kim DE, Kim B, Shin HS, Kwon HJ, Park ES. The protective effect of hispidin against hydrogen peroxide-induced apoptosis in H9c2 cardiomyoblast cells through Akt/GSK-3 beta and ERK1/2 signaling pathway. Exp Cell Res. 2014;327(2):264–75.

    Article  CAS  Google Scholar 

  22. Mo S, Wang S, Zhou G, Yang Y, Li Y, Chen X, Shi J. Phelligridins C-F: cytotoxic pyrano [4,3-c][2] benzopyran-1,6-dione and furo [3,2-c]pyran-4-one derivatives from the fungus Phellinus igniarius. J Nat Prod. 2004;67(5):823–8.

    Article  CAS  Google Scholar 

  23. Zan LF, Qin JC, Zhang YM, Yao YH, Bao HY, Li X. Antioxidant Hispidin derivatives from medicinal mushroom Inonotus hispidus. Chem Pharm Bull. 2011;59(6):770–2.

    Article  CAS  Google Scholar 

  24. Zhang C, Liang J, Yang L, Sun B, Wang C. De novo RNA sequencing and Transcriptome analysis of Monascus purpureus and analysis of key genes involved in Monacolin K biosynthesis. PLoS One. 2017;12(1):e0170149.

    Article  Google Scholar 

  25. Han JJ, Bao L, He LW, Zhang XQ, Yang XL, Li SJ, Yao YJ, Liu HW. Phaeolschidins A-E, five hispidin derivatives with antioxidant activity from the fruiting body of Phaeolus schweinitzii collected in the Tibetan plateau. J Nat Prod. 2013;76(8):1448–53.

    Article  CAS  Google Scholar 

  26. Linz JE, Chanda A, Hong SY, Whitten DA, Wilkerson C, Roze LV. Proteomic and biochemical evidence support a role for transport vesicles and endosomes in stress response and secondary metabolism in Aspergillus parasiticus. J Proteome Res. 2012;11(2):767–75.

    Article  CAS  Google Scholar 

  27. Sekiya M, Izumisawa S, Iwamoto-Kihara A, Fan Y, Shimoyama Y, Sasaki M, Nakanishi-Matsui M. Proton-pumping F-ATPase plays an important role in Streptococcus mutans under acidic conditions. Arch Biochem Biophys. 2019;666:46–51.

    Article  CAS  Google Scholar 

  28. Yagi T, Matsuno-Yagi A. The proton-translocating NADH-quinone oxidoreductase in the respiratory chain: the secret unlocked. Biochemistry. 2003;42(8):2266–74.

    Article  CAS  Google Scholar 

  29. Liu JY, Chang MC, Meng JL, Feng CP, Wang Y. A comparative proteome approach reveals metabolic changes associated with Flammulina velutipes mycelia in response to cold and light stress. J Agr Food Chem. 2018;66(14):3716–25.

    Article  CAS  Google Scholar 

  30. Wang W, Zhao P, Zhou XM, Xiong HX, Sun MX. Genome-wide identification and characterization of cystatin family genes in rice (Oryza sativa L.). Plant Cell Rep. 2015;34(9):1579–92.

    Article  CAS  Google Scholar 

  31. Shen M, Zhao DK, Qiao Q, Liu L, Wang JL, Cao GH, Li T, Zhao ZW. Identification of glutathione S-transferase (GST) genes from a dark septate endophytic fungus (Exophiala pisciphila) and their expression patterns under varied metals stress. PLoS One. 2015;10(4):e0123418.

    Article  Google Scholar 

  32. Zhang X, Leger RJS, Fang WG. Stress-induced pyruvate accumulation contributes to cross protection in a fungus. Environ Microbiol. 2018;20(3):1158–69.

    Article  CAS  Google Scholar 

  33. Nambudiri AM, Vance CP, Towers GH. Effect of light on enzymes of phenylpropanoid metabolism and hispidin biosynthesis in Polyporus hispidus. Biochem J. 1973;134(4):891–7.

    Article  CAS  Google Scholar 

  34. Purtov KV, Petushkov VN, Baranov MS, Mineev KS, Rodionova NS, Kaskova ZM, Tsarkova AS, Petunin AI, Bondar VS, Rodicheva EK, et al. The chemical basis of fungal bioluminescence. Angew Chem Int Edit. 2015;54(28):8124–8.

    Article  CAS  Google Scholar 

  35. Oba Y, Suzuki Y, Martins GNR, Carvalho RP, Pereira TA, Waldenmaier HE, Kanie S, Naito M, Oliveira AG, Dorr FA, et al. Identification of hispidin as a bioluminescent active compound and its recycling biosynthesis in the luminous fungal fruiting body. Photochem Photobiol Sci. 2017;16(9):1435–40.

    Article  CAS  Google Scholar 

  36. Lee IK, Yun BS. Highly oxygenated and unsaturated metabolites providing a diversity of hispidin class antioxidants in the medicinal mushrooms Inonotus and Phellinus. Bioorg Med Chem. 2007;15(10):3309–14.

    Article  CAS  Google Scholar 

Download references


Not applicable.


The project was financially supported by the National Natural Science Foundation of China (NNSFC; Project No. 81960643 and 81560567) and Key R&D projects in Ningxia province (Project No. 2020BFG02006). The funds helped purchase reagent and sample detection, while played no role in the design of the study, interpretation of data and in writing the manuscript.

Author information

Authors and Affiliations



WXL and JZB responsible for the overall arrangements for the study. GJJ and LXX made substantial contributions to the experimental operation and data acquisition. GJJ drafted the manuscript and JHY revised it critically for important intellectual content. LYJ creates charts and illustrations. LC, ZL, HY and BCC analyzed and interpreted the data. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Xiuli Wu.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figs. S1–9.

MS, 1D and 2D NMR spectra of compounds A and B. This material is available free of charge via the Internet at

Rights and permissions

Open Access This 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 The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Guo, J., Liu, X., Li, Y. et al. Screening for proteins related to the biosynthesis of hispidin and its derivatives in Phellinus igniarius using iTRAQ proteomic analysis. BMC Microbiol 21, 81 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: