Polysaccharide Fraction Extracted from Endophytic Fungus Trichoderma atroviride D16 Has an Influence on the Proteomics Profile of the Salvia miltiorrhiza Hairy Roots

Trichoderma atroviride develops a symbiont relationship with Salvia miltiorrhiza and this association involves a number of signaling pathways and proteomic responses between both partners. In our previous study, we have reported that polysaccharide fraction (PSF) of T. atroviride could promote tanshinones accumulation in S. miltiorrhiza hairy roots. Consequently, the present data elucidates the broad proteomics changes under treatment of PSF. Furthermore, we reported several previously undescribed and unexpected responses, containing gene expression patterns consistent with biochemical stresses and metabolic patterns inside the host. In summary, the PSF-induced tanshinones accumulation in S. miltiorrhiza hairy roots may be closely related to Ca2+ triggering, peroxide reaction, protein phosphorylation, and jasmonic acid (JA) signal transduction, leading to an increase in leucine-rich repeat (LRR) protein synthesis. This results in the changes in basic metabolic flux of sugars, amino acids, and protein synthesis, along with signal defense reactions. The results reported here increase our understanding of the interaction between T. atroviride and S. miltiorrhiza and specifically confirm the proteomic responses underlying the activities of PSF.


Introduction
Salvia miltiorrhiza is a Chinese traditional medical herb widely used for preventing and treating disorders of liver, vascular, menstrual, and blood circulation systems. Diterpenoid pigments are the main bioactive constituents of S. miltiorrhiza roots, and exert anti-inflammatory properties via significant inhibition of production of NO, IL-1β, and TNF-α [1,2]. In particular, tanshinones, the predominant active constituents in the roots of S. miltiorrhiza, possess various promising bioactivities, such as anti-inflammation, anticoagulation, and liver protection, etc. [2,3]. However, the current production of tanshinones is not meeting the current medicinal market needs. Hairy root can yield a

The Culture and Treatment of Hairy Root
Simple hairy roots were separated into different 250 mL conical flasks equipped with 100 mL half-strength B5 liquid medium at 25 • C, 135 rpm in the dark. Every Erlenmeyer flask was loaded with 0.1 g hairy root and was cultured for 3 weeks and the liquid medium was changed weekly. The blank and the treatment group were cultured in half B5 liquid medium, except that the treatment group also contained 60 mg/L PSF. The hairy root samples were taken out for further experiments on the third and sixth days after treatment commenced.
The samples were handled with liquid nitrogen to remove the moisture and then deposited in −80 • C refrigerator. The blank group samples taken on day 3 were numbered A1 and A2, while treatment group samples were numbered B1 and B2. Similarly, the blank group samples taken on day 6 were numbered C1 and C2, while treatment group samples were numbered D1 and D2. However, hairy root samples were disposed by the same method on the twelfth day after treatment for RNA extraction.

The Preparation of PSF
The hyphae of T. atroviride D16 was transferred into the newly-made Potato Dextrose Agar (PDA) medium for 4 days and then moved into 250 mL conical flasks containing 100 mL half B5 liquid medium and cultured for 3 days. Finally, the cultured D16 was added into 5 L conical flask for amplified culturing between 7 to 10 days. Subsequently, the mycelia were collected with vacuum suction filtration, washed three times with distilled water and homogenized by placing in a blender for 10 min, followed by addition of three volumes of distilled water equivalent to the wet weight of mycelia. The water solution of mycelia was then heated to a temperature of 121 • C under high pressure for 60 min and was filtered as fungal elicitor. The filtrate was decompressed and concentrated to an appropriate volume at 75 • C and mixed with 4 times volume of 95% ethanol for 2 days at room temperature. After lyophilization, the precipitate was further subjected to deproteinization with Sevag reagent (chloroform:n-butanol, 4:1, v/v), and small molecule impurities were removed as well as proteins less than 2000 kDa. Thereafter, the solution was freeze-dried under vacuum into powder as PSF.

Monosaccharides Composition of PSF
A total of 5.00 mg D16 PSF was hydrolyzed with trifluoroacetic acid (TFA, 2 M) for 2 h at 120 • C hermetically, the TFA was then evaporated under vacuum, and the residue was subsequently washed three times with methanol. The residue was then mixed with 30 mg NaBH 4 and 1mL H 2 O and left overnight for reduction followed by the addition of methanol and acetic acid at a ratio of 5:1, the solvent was evaporated under reduced pressure, followed by the addition and evaporation of methanol three times. The sample was then left to dry at 105 • C for 10 min followed by acetylation of the sample by the addition of 3 mL acetic anhydride (AC 2 O), for 60 min at 105 • C, after which the reaction was terminated by the addition of 1 mL water followed by vortex mixing. The acetylation products were extracted with 2 mL chloroform, which was then dried by the addition of anhydrous sodium sulfate prior to GC-MS analysis. The standard products (including rhamnose, arabinose, xylose, fructose, mannose, glucose, and galactose, which were purchased from the National Institute for Food and Drug Control, Beijing, China) were subjected to the same procedure as the sample without hydrolyzation.
GC-MS analysis conditions: The size of TR-5MS (Thermo) chromatographic column was 60 m × 0.25 mm × 0.25 mm; initial conditions: The programmed temperature was 140 to 198 • C at 2 • C/min and kept for 4 min. The temperature was then increased to 214 • C at 4 • C/min and 217 • C at 1 • C/min for 4 min. Finally, the temperature was raised to 250 • C at 3 • C/min for 5 min; injection port temperature was 250 • C; the carrier gas was He; the flow rate was 1 mL/min; ion source 250 • C; m/z = 40-500; EI = 70 ev.
The methanol extract of hairy roots was analyzed and the peaks identified and contrasted in comparison with the available standards.

RNA Isolation and Real-Time Quantitative PCR Analysis
TRIzol reagent (Invitrogen Co., Carlsbad, CA, USA) was used to extract the total RNA of hairy root and 1 µL RNA was diluted to 100 µL with RNase ddH 2 O. According to the RNA concentration ratio of 260 and 280 nm UV absorbance value, OD260/280 should be between 1.9 to 2.1 and the total RNA concentration was calculated according to the formula total RNA (µg/mL) = OD260 × 100 (dilution multiple) × 40 µg/mL. The reverse transcription reaction was conducted at 37 • C for 15 min, 85 • C for 5 s termination reaction at ABI 9700 PCR. Finally, 90 µL RNase-free ddH 2 O was added to 100 µL for storing at −20 • C for Real-time PCR. Primers with the following sequences, LRRK F (5 -TGTGGTAGCTTTGTGGGGTT-3 ) and LRRK R

Protein Extraction of Salvia Miltorrhiza Hariy Root
A total of 0.5000 g of S. miltiorrhiza hairy root sample was grounded into powder in liquid nitrogen and transferred into 10% trichloroacetic acid and precooling acetone solution containing 65 mM Dithiothreitol (DTT). Followed precipitation at −20 • C for 6 h, the sample was centrifuged at 10,000 rpm, at 4 • C for 45 min, and the sediment was re-suspended by the addition of PH standard reagent (STD buffer) according to the volume ratio of 10:1. The sample was then vortexed and blended and placed in a boiling water bath for 5 min and then subjected to ultrasound 10 times (ultrasound for 10 s, interval for 15 s). Finally, the sample was heated in a water bath for 5 min and subjected to centrifuging for 15 min at 10,000 rpm. The supernatant containing the total protein extract of S. miltiorrhiza hairy root was then stored at −80 • C until assayed for protein content by the bicinchonininc acid (BCA) and Bradford method.

The Mass Spectrum Analysis of Protein
The protein samples were enzymolyzed and the peptides solution was qualified by measuring at OD280. The sample were then analyzed by LC-MS preliminarily and performed on AKTA Purifier 100 (GE Healthcare, Chicago, IL, USA) for SCX classification via strong cation chromatographic column. Finally, the samples were separated by the Easy nLC liquid phase system of the nanoliter velocity. Q-Exactive mass spectrometer (Thermo Finnigan, Silicon Valley, CA, USA) was used to analyze the mass spectra of the sample solution after capillary HPLC separation. The analysis time was 120 min, testing method was positive ion mode, and the mother ion scanning ranged from 300 to 1800 m/z. The level of mass spectrum resolution was 70,000 (m/z 200), AGC target was 3e6, first level of Maximum IT was 10 ms, number of scan ranges was 1, and dynamic exclusion was 40 s. The ion-charge ratio of peptide and the peptide fragments was collected according to the following methods: Full scan acquiring 10 secondary debris each map; secondary excitation types: heated capillary dissociation (HCD); separate window: 2 m/z; secondary mass spectrum resolution: 17,500 (200) m/z; microscans: 1; secondary maximum IT: 60 ms; standard collision energy: 30 eV; under fill the wire: 0.1%.

Data Analysis
The raw data of mass spectrometry was analyzed with software Mascot 2.2 and the identified proteins were blasted with UniProt database for the Gene Ontology (GO) classification annotation. The peak strength value of the reported ion was analyzed quantitatively with Proteome Discoverer1.4 software and the differentially expressed genes of Fc (fold change, Fc (B/A) = B/A, Fc (D/C) = D/C) were screened according to the standard of Fc ≥ 1.50 or Fc ≤ 0.66. The differentially expressed proteins identified in the Arabidopsis database were compared to the database (UniProt, SWISSPROT and TREMBL) by using the software BLASTP (2.2.23+) and the comparison results were extracted. These comments were corresponded to the GO number by the gene association file (gene association. Goa uniprot) and the number of the corresponding small items were counted in the three large catalogs of the biological process, cellular component, and molecular function. According to the bi-directional best hits (BBH) analysis, the corresponding number was mapped to KEGG pathway maps and the total protein of Arabidopsis was calculated by using the hyper-geometric distribution of p-value. The false discovery rate (FDR) was then used to correct it and z-score was calculated at the same time. The pathway of significant enrichment was the pathway of FDR is less than 0.05 and z-score > 0. INTACT and String was used to construct differentially expressed protein interaction networks: (a) Searching INTACT to build network (http://www.ebi.ac.uk/intact/); (b) using the software String to build the network (http://string.embl.de).
All experiments, including both control and different treatments of hairy root cultures, HPLC analysis, and semi-quantitative real-time PCR, were performed in triplicate. The results are presented as their mean values and standard deviations (SD). The error bars in the figures represent the standard deviation in biological triplicates. The statistical significance of the differences in root growth and the accumulation of phenolic acids and tanshinones were analyzed by one-way analysis of variance (ANOVA) with SPSS software (version 18.0, Chicago, IL, USA), the heatmap was drawn using Excel software (version: 2010, Microsoft Co., Redmond, WA, USA), and the term significant was used to denote the differences for which p-value is < 0.05. The statistical significance of differences in gene transcripts was analyzed by one-sample t-test.

Monosacharide Composition and Its Effects on Accumulation of Tanshinones of the D16 PSF
D16 PSF was extracted from T. atroviride D16 and it had greatly induced the accumulation of tanshinones in S. miltiorrhiza hairy roots, as reported in the previous study [6]. To further determine the effects of hydrolytic constitutes on S. miltiorrhiza hairy root, the D16 PSF was hydrolyzed to determine its monosaccharides compositions, and the results revealed that D16 PSF was composited by mannose, glucose, and galactose at a ratio of 5:10:2 ( Figure 1). The standard samples of mannose, glucose, and galactose were made into the solution with the concentration of 60 mg/L, respectively, and then put into the S. miltiorrhiza hairy root culture medium in triplet for 14 days. The results showed the content of tanshinones in S. miltiorrhiza hairy root after 14 days' culture varied in the presence of PSF and the monosaccharides (Table 1). Compared with the control group, PSF (60 mg/L) stimulated the tanshinones accumulation at the greatest extent. Different components of PSF had various effects on promoting synthesis of tanshinones; however, the singular result of each D16 PSF constitute could not catch up with the influence of PSF on the accumulation of tanshinones. Under the treatment of PSF (60 mg/L), the content of dihydrotanshinone I was similar to the results at the same concentration of mannose and galactose respectively, while the equally-increased effects on tanshinone I were among the different treatments of PSF (60 mg/L), glucose, and mannose. Under the treatment of PSF (60 mg/L), the content of cryptotanshinones was 3.08-fold higher than the control group, while less increased effects occurred followed by the treatment of glucose. Galactose made the cryptotanshinones content decreased, while mannose had no effect on the accumulation of cryptotanshinones.

Differential Proteomics of S. miltiorrhiza in Response to D16 PSF
The differences of protein content in S. miltiorrhiza hairy roots treated with D16 PSF were studied by proteomics technique of iTRAQ. Based on the database of Arabidopsis thaliana, 3685 feature peptides (unique peptide) belonged to 1476 proteins, which were identified and analyzed for GO classification annotation. As shown in Figure 2, most of proteins can be annotated into the three ontologies of biological process, cellular component, and molecular function. In the biological process annotation, the number of proteins in the metabolic process (1162) and response to stimulus (584) was the largest, and 187 proteins were not annotated. In cellular components annotation, positioning to the cytoplasm (1036), the cell membrane (602), and the cytosol (517) account for the greatest number of proteins, while 211 proteins had not been annotated. In the molecular function annotation, a large number of proteins were divided into catalytic activity (894), nucleotide binding (507), protein binding (417), as well as metal ion binding (315) under the entry and 146 proteins had not been annotated. Among the 1476 proteins, 89 differentially expressed proteins were picked out for quantitative analysis (Figure 3). The classification of 89 different proteins was identified and analyzed in the Arabidopsis thaliana database. In addition, the number of proteins was corresponded to each category, as shown in Figure 4. The 89 proteins corresponded to 39 KEGG pathways, including amino acid metabolism, sugar metabolism, energy metabolism, and secondary metabolites synthesis. The first 10 KEGG pathways were annotated as shown in Table 2. The 89 different proteins were analyzed by adopting the INTACT and String interaction networks and the results are shown in Figure 5.   Table 2. The 10 KEGG pathways annotated by major differential proteins.

Quantitative Analysis of the 89 Differential Proteins
For the quantitative analysis, 89 differential proteins were mainly involved in signal transduction, redox, amino acid synthesis and metabolism, protein synthesis and metabolism, the biosynthesis of terpenoids, carbohydrate synthesis, metabolism and transport, and other physiological and biochemical process (Table 3, Figure 6). The 89 differential proteins in signal transduction involves the leucine-rich repeat (LRR) protein, Ca 2+ -related protein, peroxidase-related protein, protein phosphorylation protein, and the synthesis of jasmonic acid (JA). We found that two LRR-related proteins significantly increased leucine-rich repeat protein kinase family protein and leucine-rich repeat receptor-like protein kinase PEPR2 under the action of T. atroviride D16PSF on the third day and sixth day in this study, which revealed that the LRR protein plays an important role in promoting S. miltiorrhiza hairy root growth and tanshinone synthesis under the action of T. atroviride D16 PSF. The relative expression of LRR and LRR-PEPR on the fourteenth day of S. miltiorrhiza after culturing in the presence of D16 PSF was also detected. The results revealed that under the sustained action of D16 PSF, a trend showing an increase in LRR and LRR-PEPR is observed (Figure 7). In our study, we also detected that Ca 2+ -related protein expression changed significantly. The cation-transporting ATPase and calcium-binding protein increased significantly at day three and calmodulin-binding protein (AT3g52870/F8J2_40) decreased significantly on the sixth day and the relative expression of cation-transporting ATPase still increased greatly on the fourteenth day ( Figure 7). 4-coumarate-CoA ligase-like 10 increased significantly responding to Ca 2+ on the sixth day and nucleoporin related to the level of intracellular Ca 2+ decreased on the third day, but increased on the sixth day significantly [18]. These significant changes of Ca 2+ -related proteins showed calcium plays a very important role under the action of T. atroviride D16 PSF in S. miltiorrhiza hairy root. At the same time, we also found superoxide dismutase, peroxidase, and peroxidase 50 associated with peroxide was lowered significantly. 12-oxophytodienoate reductase is a very important key enzyme in the jasmonic acid (JA) biosynthetic pathway, and two proteins, 12-oxophytodienoate reductase-like protein and 12-oxophytodienoate reductase (At1g04380), related with jasmonic acid were found to be significantly higher in this study [19]. Therefore, the effect of T. atroviride D16 PSF on S. miltiorrhiza hairy root may also be related to peroxide, protein phosphorylation, and jasmonic acid. The 89 differential proteins are mainly involved at the levels of DNA, RNA, and protein translation in protein synthesis. At the level of DNA, the protein-DNA replication licensing factor MCM3 appeared significantly decreased after the first significant rise. DNA and protein damage repair/toleration protein expression was significantly lower on the sixth day. At the level of RNA, transcription factor ABA-inducible bHLH TYPE and transcription factor Pur-alpha 1 was raised significantly on the third day; the pre-mRNA splicing factor SF2-like protein and nuclear protein was also significantly increased. At the level of protein translation, 40 S ribosomal protein was raised significantly on the third day and the rest of the proteins related to the protein degradation were significantly lower such as chaperone DNAJ domain containing protein, the ubiquitin-like 5 and E3 ubiquitin protein-protein ligase XBAT33. These changes showed T. atroviride D16 PSF played a positive regulatory role in protein synthesis of S. miltiorrhiza hairy root at the protein expression level.
In terms of the differential proteins in secondary metabolism, they were mainly involved in the cytochrome P450 enzymes, redox enzyme, and enzymes related to the terpenoid biosynthetic pathway, and most of these proteins increased significantly under the influence of T. atroviride D16 PSF. At the same time, significant changes have taken place in some metabolic activities, like synthesis and metabolism of sugar for providing energy, the amino acid synthesis, and metabolism to provide raw materials for protein synthesis and transporters which have transportation function.
With the deep analysis of these differential expression proteins, it can be presumed that D16 PSF regulates the synthesis and metabolism of saccharides and amino acids, the transcription of genes and the translation of protein through the signal transduction pathways involved in leucine-rich repeat proteins, calcium, peroxides, protein phosphorylation, and jasmonic acid (JA) synthesis. This then results in changes of protein expression involved in secondary metabolism and thus induces metabolic profile changes in S. miltiorrhiza hairy roots.  Transcription factor Pur-alpha 1 1.91 -Note: B/A and D/C are 3 and 6 days, respectively, dealing with the ratio of the amount of protein group and the blank group, 1.5 or higher to increase protein, 0.66 or less for cut protein, "-" indicates no difference. GASA5, GA-Stimulated in Arabidopsis; LRR, Leucine-rich repeat.

Discussion
It is known that protein is the carrier of gene function and the executor of the life activities. The high-throughput and large-scale proteomics to study time changes of the expression of all the proteins in the cell or tissue can reveal the physiological and biochemical processes in detail [20,21]. At present, proteomics research mainly is based on the separation of the gel and the liquid phase. The methods based on gel separation include two-dimensional electrophoresis proteomics research and fluorescent difference gel electrophoresis (DIGE) proteomics research, mainly through the gel electrophoresis separation of mixed protein samples. After getting various single protein points, we can screen differential proteins according to the dyed color and identify proteins through the MALDI-MS finally. Methods based on the liquid phase include label-free proteomics research, marked iTRAQ/SILAC proteomics research, and SRM/MRM proteomics research, mainly by means of liquid phase to separate mixed protein samples and through the mass spectrometer connected directly with liquid phase for qualitative and quantitative analysis [22]. Proteomics technology has been widely used in plant growth and development, plant adaptation mechanism to biological and abiotic stress, and the interaction mechanism of microbes [23,24]. iTRAQ (isobaric tags for relative and absolute quantitation) technology is a kind of proteomics technology through four or eight kinds of isotope labels specifically to mark amino groups of polypeptides, and then to identify and separate proteins by LC-MS and analyze relative and absolute content of different protein samples according to the isotope intensity [25,26]. This experiment adopted the iTRAQ technology to tag eight S. miltiorrhiza hairy root samples of control group and T. atroviride D16 PSF on the third and sixth day.
Leucine-rich repeat (LRR) protein kinase family proteins are the largest known type of transmembrane receptor protein kinases in plants, which have the function of regulation in plant growth and development, hormone signal transduction, and biological and non-biological stress response [27,28]. LRR proteins play a critical role in the process of Piriformospora indica promoting growth of Arabidopsis thaliana [29,30]. As studied, Pep1 is a 23-amino acid peptide that enhances resistance to a root pathogen, Pythium irregulare. Pep1 and its homologs (Pep2 to Pep7) are endogenous amplifiers of innate immunity of Arabidopsis thaliana that induce the transcription of defense-related genes and bind to PEPR1, a plasma membrane leucine-rich repeat (LRR) receptor kinase [31]. Our present results revealed that T. atroviride D16 PSF could gradually enhance the LRR protein kinase family protein and LRR receptor-like protein kinase PEPR2 in hairy roots of S. miltiorrhiza. The LRR receptor-like protein kinase PEPR2 might strengthen the ability of resisting T. atroviride D16 PSF, whereas the expression of probable LRR receptor-like serine/threonine-protein kinase At4g20940 decreased. CML42 are calcium-binding proteins that are thought to function as plant signal transduction elements, and were up-regulated and induced by Spodoptera littoralis in Arabidopsis thaliana followed by Ca 2+ and phytohormone elevation [32,33]. In addition, CML42 enhanced the content of aliphatic glucosinolate and hyperactivated transcript accumulation of the jasmonic acid (JA)-responsive genes through the negative regulation of the jasmonate receptor Coronatine Insensitive1 (COI 1). As proteomics analysis revealed, this was expressed 1.7-fold higher in S. miltiorrhiza hairy root on the third day, which indicated T. atroviride D16 PSF may regulate Ca 2+ and JA elevation. It has been reported that CML42 was not only involved in abiotic stress responses and insect herbivory defense, but also related to trichomes branching and endophytic fungi stimulation as a Ca 2+ sensor, indicating CML42 is an important calcium-binding protein in plant growth and defense processes. Nucleoporin is a component of the nuclear pore complex, which has strongly attracted the attention of its involvement in hormonal and pathogen/symbiotic signaling [34,35]. Ca 2+ signal transduction is important in the interaction between microorganisms and plants [36,37]. In the present study, as a key nucleoprotein, the Nup133/Nup155-like protein had a sharp rise on the third day as well as a substantial decrease in S. miltiorrhiza hairy root under the treatment of T. atroviride D16 PSF. The present results revealed that Nup133 are essential for mRNA export, and the Nup133/Nup155-like protein may regulate the symbiotic signaling transduction and S. miltiorrhiza response via mRNA export.
Cation-transporting ATPase is an enzyme protein widely found in biological membranes. Its structure and function is complex and plays an important regulatory role in the biological activities of cells [38]. PDE1 is an encoded P-type ATPase, which is required for the maintenance of phospholipid asymmetry in biological membranes, as membrane asymmetry may play a critical role in the development of infection hyphae by phytopathogenic fungi. The cation-transporting ATPase was up-regulated and induced by T. atroviride D16 PSF on the third day, and may provide a chance for interacting with S. miltiorrhiza hairy root. The 4-coumarate-Co-A ligase (4CL)-like proteins belong to the adenosine monophosphate (AMP)-binding domain-containing proteins family and widely exist in various plant species [39]. AMP-binding domain-containing 4CLs are critical enzymes in the phenylpropanoid metabolism pathway and drive the carbon flow from primary metabolism to different branches of secondary metabolism in plants [40]. Along with the tanshinones accumulation, the 4-coumarate-Co-A ligase (4CL)-like 10 protein were raised on the sixth day. In support, the 4-coumarate-Co-A ligase (4CL) like protein was up-regulated by Magnaporthe oryzae infection, which may be a defense-related AMP-binding protein (AMPBP) that is involved in the regulation of the defense response through salicylic acid (SA) and/or jasmonic acid (JA)/ ethylene (ET) signaling pathways. Phosphoprotein phosphatases (PPP) are present in all eukaryotic organisms, which is an ancient and important regulatory enzyme. The protein phosphatase 2Cs (PP2Cs) from various organisms have been implicated to act as negative modulators of protein kinase pathways involved in diverse environmental stress responses and developmental processes [41]. The Ser/Thr-specific phosphatases are metal-dependent enzymes divided into two major families: The PPP family, which includes protein phosphatases 1, 2A, and 2B (PP1, PP2A, PP2B/calcineurin); and the PPM family, which includes PP2C, which is a highly ABA-induced protein in guard cells once found in Arabidopsis and induced by endophyte fungi T. atroviride D16 PSF as well in this study. In addition, dual specificity protein phosphatases were significantly higher on the third day, referring to serine/threonine protein phosphatase and probable protein phosphatase related to protein phosphorylation [42]. T. atroviride D16 PSF induced high expression of protein phosphatase 2C and serine/threonine-protein phosphatase PP2A-2 through ABA signaling in S. miltiorrhiza hairy root. 12-oxophytodienoate reductase-like protein 1 is an enzyme in the jasmonic acid (JA) biosynthesis pathway [43], and the expression of 12-oxophytodienoate reductase-like protein 1 was gradually enhanced under the treatment of T. atroviride D16 PSF involved in the JA biosynthetic pathway. Superoxide dismutases (SODs) are involved in plant adaptive responses to biotic and abiotic stresses although the upstream signaling process which modulates their expression is not clearly understood [44]. Surprisingly, the level of SOD and peroxides were decreased at the treatment of T. atroviride D16 PSF, while plant overexpressing antioxidant enzymes had higher tolerance to external stress [45]. Phospholipid hydroperoxide glutathione peroxidase (PHGPX) is the principal enzymatic defense against oxidative destruction of biomembranes [46], and its major function appears to be the scavenging of phospholipid hydroperoxides, thereby protecting cell membranes from peroxidative damage. Gene expression analysis has shown an increase in the levels of PHGPX mRNA in several plant species undergoing different biotic and abiotic stresses, such as pathogen infections, high salt concentrations, exposure to heavy metals, mechanical stimulation, aluminum toxicity, seed germination, salt and osmotic stress, oxidative stress, and chilling stress [47]. PHGPXs may play dual roles as a redox transducer in addition to acting as a H 2 O 2 scavenger under stress, thus PHGPX proteins may have different functions in plant cells, with some isoforms functioning in the signal transduction pathway, while others are involved in catalyzing the reduction of harmful products formed by ROS. The effect of various signaling molecules (salicylic acid, JA, ABA, ethylene) and certain protein phosphatase inhibitors (cantharidin and endothall) on the expression of the PHGPX gene in rice seedling leaves has demonstrated an up-regulation of the mRNA levels. These data suggest the role of this enzyme in the induction of defense mechanisms in plant cells subjected to oxidative damage, as a result of exposure to various environmental stresses, is by reducing the phospholipid hydroperoxides formed in the biomembranes.

Conclusions
Collectively, the analysis indicated that 89 differentially abundant proteins were involved mainly in protein synthesis, protein folding and degradation, biotic stress defense, photosynthesis, RNA process, signal transduction, and other functions. When induced by T. atroviride D16 PSF, S. miltiorrhiza hairy roots generally respond through elucidating the synthesis of tanshinones. These results provided valuable information for T. atroviride D16 PSF inducing tanshinones accumulation of S. miltiorrhiza. Characterization of PSF and proteomics provided an important bioinformatic resource for investigating mechanisms in inducing tanshinones accumulation.