Skip to main content

Comparative phosphoproteomic analysis of blast resistant and susceptible rice cultivars in response to salicylic acid



Salicylic acid (SA) is a significant signaling molecule that induces rice resistance against pathogen invasion. Protein phosphorylation carries out an important regulatory function in plant defense responses, while the global phosphoproteome changes in rice response to SA-mediated defense response has not been reported. In this study, a comparative phosphoproteomic profiling was conducted by two-dimensional gel electrophoresis (2-DE) and mass spectrometry (MS) analysis, with two near-isogenic rice cultivars after SA treatment.


Thirty-seven phosphoprotein spots were differentially expressed after SA treatment, twenty-nine of which were identified by MALDI-TOF/TOF MS, belonging to nine functional categories. Phosphoproteins involved in photosynthesis, antioxidative enzymes, molecular chaperones were similarly expressed in the two cultivars, suggesting SA might alleviate decreases in plant photosynthesis, regulate the antioxidant defense activities, thus improving basal resistance response in both cultivars. Meanwhile, phosphoproteins related to defense, carbohydrate metabolism, protein synthesis and degradation were differentially expressed, suggesting phosphorylation regulation mediated by SA may coordinate complex cellular activities in the two cultivars. Furthermore, the phosphorylation sites of four identified phosphoproteins were verified by NanoLC-MS/MS, and phosphorylated regulation of three enzymes (cinnamoyl-CoA reductase, phosphoglycerate mutase and ascorbate peroxidase) was validated by activity determination.


Our study suggested that phosphorylation regulation mediated by SA may contribute to the different resistance response of the two cultivars. To our knowledge, this is the first report to measure rice phosphoproteomic changes in response to SA, which provides new insights into molecular mechanisms of SA-induced rice defense.


Rice (Oryza sativa L.) is an economically important cereal crop throughout the world, providing food for over 50% of global population [1]. The ascomycetous fungus Magnaporthe oryzae causes the rice blast, one of the most devastating fungal diseases in rice production and thus poses a great threat to the world’s food security [2]. Thus far, the disease control is mainly based on using fungicides and breeding resistant cultivars. However, fungicide could fail to satisfy the requirement of environment and human health regulations, and resistant cultivars could be overcome by the quick arising/evolving of new races of M. oryzae [3, 4]. As an important signaling molecule, salicylic acid (SA) can induce plant resistance against multiple fungal, viral and bacterial pathogens [5]. SA-mediated plant defense responses are actively involved in both PTI (PAMP-triggered immunity) and ETI (effector-triggered immunity) [6]. In the previous study, we confirmed that SA can protect rice against infection by M. oryzae race ZC13 in CO39 (susceptible cultivar) and in a near isogenic line C101LAC, which carries the resistance gene Pi-1 against M. oryzae race ZC13 and thus represents a resistant cultivar [7]. Proteomic analysis further showed that SA coordinates multiple cellular activities to facilitate defense response and recovery in both rice cultivars. However, it awaits further elucidation about the detailed molecular mechanisms of SA-induced rice defense response against M. oryzae.

Phosphorylation is one of the most important post-translational protein modifications (PTMs), regulating a wide range of cellular functions in various organisms, including cell signaling, metabolism, stress responses and defense responses [8]. Phosphoproteomics can capture the dynamics and specificity of protein phosphorylation, and therefore enhance our understanding of fundamentals and complex biological processes [9]. In recent years, a large number of emerging evidences suggested that protein phosphorylation can regulate plant stress responses triggered by exogenous hormone and biotic stress [10,11,12,13]. It has been also shown that protein phosphorylation was involved in the activation of SA-induced plant resistance using traditional biochemical methodologies [14, 15]. To our knowledge, the general scope of such connections between protein phosphorylation events and SA-induced rice resistance has not been studied.

In the present study, we performed a comparative phosphoproteome to reveal the detail SA-induced mechanism in rice using two-dimensional gel electrophoresis (2DE), Pro-Q diamond phosphoprotein stain and MALDI-TOF/TOF mass spectrometer. Thirty-seven SA-responsive phosphoprotein spots were found and twenty-nine of them were identified. Phosphoproteins involved in similar or different function and expression patterns in resistant and susceptible rice cultivars were discussed. The results provided new insights about the dynamic phosphoproteomes at different time points in rice upon SA treatment, and broadened the understanding of SA-mediated rice resistance against M. oryzae infection via regulation on protein phosphorylation.


Specificity analysis of MOAC-enriched putative phosphoproteins from rice leaves

A total of 490 ± 15 μg putative phosphoproteins were enriched from 8 mg of total proteins. To test the specificity of MOAC for phosphoproteins, MOAC-enriched putative phosphoproteins were separated by 2DE and sequentially stained for phosphoproteins using Pro-Q Diamond, and for total proteins by silver stain (Fig. 1a, b). 481 ± 9 protein spots could be detected on Pro-Q Diamond-stained gels and 469 ± 12 spots on sequential silver-stained gels; of these spots, 466 were common to the two staining methods (Fig. 1c). To determine the specificity of MOAC-enriched phosphoproteins, the sequential staining images were overlapped and the protein spots were visualized in different colors using PDQuest software (Fig. 1d, e, f). The overlay image showed that most of the protein spots (over 99%) appeared in yellow, indicating that these proteins were phosphoproteins; only three protein spots appeared in red, indicating that these were nonphosphoproteins (Fig. 1f). The results demonstrated the MOAC is selective enough for detecting rice leaf phosphoproteins.

Fig. 1

2DE analysis of MOAC-enriched putative phosphoproteins from rice leaves. Original-color image of the same 2DE gel was stained by a Pro-Q diamond and b silver nitrate. c Venn diagram analysis of MOAC-enriched putative phosphoproteins in 2DE gels that overlapped between Pro-Q Diamond staining (green) and silver staining (red). False-color images of 2DE gels were visualized with different colors using PDQuest software; d Pro-Q diamond-stained protein spots were colored green and e silver-stained protein spots were colored red. f An overlay of the two images (d and e). The phosphoprotein spots were appeared yellow, and the non-phosphoprotein spots were appeared red (as shown by arrows). Asterisks indicate spots that were very abundant in the silver image but lightly stained in the Pro-Q Diamond image

Phosphoproteome changes in rice leaves upon SA treatment

To study the changes of SA-induced phosphoprotein profiles, we conducted a 2DE-based phosphoproteomic analysis at 12 h and 24 h after SA treatment of C101LAC and CO39. In our previous work, we conducted SA-induced rice resistance against M. oryzae at different concentrations ranging from 0.01 mM to 1 mM SA [7]. The results showed that the optimum concentration of SA treatment to induce blast resistance of rice seedlings was at 0.1 mM. Thus, 0.1 mM SA was used in this study. At least three independent 2DE analysis was performed for each treatment, with a high level of reproducibility. Eight representative gels and other replicate gels were shown in Additional file 1: Figure S1 and Figure S2, respectively (Additional file 1). Fold changes above 1.5 in all three replicates were used as thresholds to determine the SA-responsive phosphoproteins. Additionally, differential expression pattern was shown to be similar in all three replicates, and we manually checked all the spots to ensure confidence in differentially regulated phosphoproteins. A total of 37 SA-responsive phosphoproteins were obtained from the two cultivars (Additional file 1: Figure S1). For each cultivar, 30 and 28 SA-responsive phosphoprotein spots were detected in CO39 and C101LAC, respectively; 21 of which were common in these two cultivars. A close-up view of the SA-responsive phosphoproteins on the 2DE gels was shown in (Additional file 1: Figure S3). The relative intensities of SA-responsive phosphoproteins were displayed in (Additional file 1: Fig. S4). We conclude that SA treatment resulted in reproducible and significant changes to these protein spots, which we could investigate further.

MALDI-TOF/TOF identification of SA-responsive phosphoproteins

The 37 SA-responsive phosphoprotein spots were excised from 2DE gels and further identified by MALDI-TOF/TOF MS. Of these, 29 phosphoproteins were identified with high confidence (Table 1), while the remaining 8 phosphoprotein spots (2, 10, 19, 22, 23, 26, 32 and 35) did not show a good match in the database. The spectra of protein spot 11 were provided as an example of analysis (Additional file 1: Figure S5). Generally, one protein spot in 2DE gel represented one unique protein. However, we noted that three phosphoproteins were identified in more than one spot in the same gel (Table 1, Additional file 1: Figure S1). For example, three phosphoprotein spots (11, 12 and 13) were identified as probable glutamyl endopeptidase, two (17 and 18) as glyceraldehyde-3-phosphate dehydrogenase, two (24 and 25) as alpha 1,4-glucan phosphorylase. Consistent with our results, it has been reported that protein isoforms migrated as a chain of spots, most likely due to post-translational modifications [16]. Proteins with multiple phosphorylation states could also lead to electrophoresis patterns that multiple spots are with similar molecular weight but different pI [17].

Table 1 SA-responsive rice phosphoproteins identified by MALDI-TOF/TOF MS

Based on biological annotations from UniProtKB database (, the 29 identified SA-responsive phosphoproteins were functionally classified into 9 groups: photosynthesis, defense, antioxidative enzymes, protein synthesis and degradation, molecular chaperones, amino acid metabolism, carbohydrate metabolism, energy metabolism and other metabolism (Additional file 1: Figure S6). Among them, phosphoproteins involved in carbohydrate metabolism, protein synthesis and degradation were the most abundant, both representing 24.14% (7/29) of the phosphoproteins identified, respectively.

Phosphorylation site identification

Next, we tried to determine the phosphorylation patterns of total proteins by using NanoLC-MS/MS analysis, following the workflow as depicted in Fig. 2a. Totally, 1815 phosphosites were identified, which come from 1537 phosphopeptides (Additional file 2: Table S1) on 839 phosphoproteins (Additional file 3: Table S2). Among these 1815 phosphorylated residues, there were 1539 phosphoserine (pS), 246 phosphothreonine (pT), and 30 phosphotyrosine (pY), corresponding to 84.79, 13.55, and 1.65% respectively of all phosphorylated residues (Fig. 2b). We noticed that 257 phosphopeptides were multiply phosphorylated, and 1279 phosphopeptides were singly phosphorylated (Fig. 2c). Then, the 29 identified SA-responsive phosphoproteins were scanned in the NanoLC-MS/MS analysis data. Phosphorylation sites within 4 identified SA-responsive phosphoproteins were verified (Additional file 4: Table S3). The MS/MS spectrum of a representative phosphorylated peptide (ELLS*YEYDGDEVPIVAGSALK), Elongation factor Tu (spot 30) was shown in (Additional file 1: Figure S7), as an example.

Fig. 2

NanoLC-MS/MS identification for phosphorylation sites of total proteins from rice leaves. a Overview of the experimental design for the phosphorylation patterns of total proteins. The total proteins were digested with trypsin. Phosphopeptides are enriched from the pooled peptide mixture with titanium dioxide (TiO2) beads and subsequently analyzed with nanoLC−MS/MS. b The distribution of peptides having one, two, three, and four and more phosphorylation sites. c The distribution of phosphorylated residues. pS, phosphoserine; pT, phosphothreonine; pY, phosphotyrosine

Enzyme activities in rice leaves induced by SA treatment

To further validate possible regulation of enzyme activity by protein phosphorylation, we selected three important enzymes from the identified phosphoproteins, and assessed their activities. A significant increase in cinnamoyl-CoA reductase (CCR) activity was detected in C101LAC at 24 h after SA treatment (Fig. 3a). CCR showed an increased phosphorylation level in C101LAC at 24 h after SA treatment, suggesting that the enzyme activity of CCR was regulated by phosphorylation. Consistent with the phosphoproteomic results, a significant decrease in APX activity was noted in both cultivars 24 h after SA treatment. A significant decrease in PGAM activity was detected only in C101LAC 12 h after SA treatment, but little change was found in CO39 (Fig. 3b, c). Phosphoproteomic analysis also showed that phosphorylated PGAM and phosphorylated APX were down-regulated only in C101LAC at 12 h after SA treatment, but little change was found in CO39 (Table 1). These results above strongly suggesting that the enzyme activity of PGAM and APX was regulated by phosphorylation.

Fig. 3

Quantitative analysis of a CCR, b PGAM, and c APX activities in the leaves of rice seedlings after SA treatment. Bars indicate ± standard error of the mean. Different small letters in each group indicate significant differences at P ≤ 0.05

ROS accumulation in rice leaves induced by SA treatment

SA and ROS production in rice defense response to biotic and abiotic stress has been well documented [18, 19]. In this study, the phosphoproteomic analysis revealed several phosphoproteins involved in biogenesis of ROS. In agreement with the phosphoproteomic results, a significant increase in O2.-, H2O2 and MDA contents was noted in both rice cultivars after SA treatment, compared with the corresponding controls (Fig. 4). Relative to the control, SA resulted in an increase of H2O2 contents by 1.77 and 2.72 times in CO39, 2.47 and 3.15 times in CO101LAC at 12 h and 24 h, respectively (Fig. 4a). An increase of O2.- content was also noted by 1.94 and 3.75 times in CO39, 2.59 and 4.46 times in CO101LAC at 12 h and 24 h after SA treatment, respectively (Fig. 4b). MDA contents was significantly increased by 1.83 and 2.61 times in CO39, 2.61 and 2.79 times in CO101LAC 12 h and 24 h post SA treatment, respectively, compared with the controls (Fig. 4c). Taken together, SA treatment significantly increased the contents of O2.-, H2O2 and MDA in both resistant (C101LAC) and susceptible (CO39) cultivars, but the fold changes in C101LAC were significantly higher than that in CO39.

Fig. 4

ROS production in rice leaves induced by SA. Water was used as a control (CK). a O2.- production. b H2O2 content. c MDA contents. Bars indicate ± standard error of the mean. Different capital letters in each group indicate significant differences at P ≤ 0.01

Transcriptional expression analysis of SA-responsive phosphoproteins

Six genes encoding the identified phosphoproteins (Additional file 4: Table S4) were selected for expression analysis via qRT-PCR. The gene encoding putative chaperonin 60 beta was significantly down-regulated in CO39 after SA treatment, but up-regulated in C101LAC (Fig. 5a). Elongation factor Tu significantly decreased in CO39 at mRNA levels, did not significantly vary in C101LAC (Fig. 5b). The expression changes of the above two phosphoproteins at mRNA levels were consistent with the 2DE results both in CO39 and C101LAC. Two phosphoproteins (eukaryotic initiation factor 4A-1 and phosphoribulokinase) showed a significant decrease in mRNA levels but a significant increase in phosphorylation levels in C101LAC; however, they decreased both in mRNA levels and phosphorylation levels in CO39 (Fig. 5c, d). Interestingly, glyceraldehyde-3-phosphate dehydrogenase significantly decreased in mRNA levels, but increased in phosphorylation levels in CO39; however, the expression decreased both in mRNA levels and phosphorylation levels in C101LAC (Fig. 5e). L-ascorbate peroxidase 1 significantly increased in mRNA levels both in CO39 and C101LAC; while its phosphorylation level showed a typical decrease in C101LAC and no significant variation in CO39 (Fig. 5f). In total, these data showed a lack of correlation between transcriptional regulation and post-tranlational regulation (by protein phosphorylation) in rice plants treated with SA. Consistent with this results, several reported phosphoproteomics studies had also revealed different changes in the mRNA levels and their corresponding proteins levels [20, 21]. The results highlighted the importance of SA-induced resistance mechanism at multiple molecular levels in rice plant.

Fig. 5

Transcript analysis by qRT-PCR of six differentially expressed genes after SA treatment. a Putative chaperonin 60 beta, b Elongation factor Tu, c Eukaryotic initiation factor 4A-1, d Phosphoribulokinase, e Glyceraldehyde-3-phosphate dehydrogenase, and f L-ascorbate peroxidase 1. Bars indicate ± standard error of the mean. Different small letters in each group indicate significant differences at P ≤ 0.05


SA is an important signaling molecule that plays key roles in the regulation of plant defense against pathogens. In our previous work, we explored the molecular mechanisms of SA-mediated protection of rice against M. oryzae infection by proteomic profiling analysis with blast-resistant vs. -susceptible rice cultivars after SA treatment [7]. Protein phosphorylation, as one of the most common and best characterized post-translational modifications, is a key process that regulates a large number of biological processed in plants, including cell signaling, metabolism, hormone and stress responses [22]. However, it is very difficult to detect phosphoproteins by the “proteomic” approach, due to the low abundance of phosphoproteins [23]. For further understanding the molecular mechanisms, we performed a comparative phosphoproteome analysis of two near rice isogenic lines after SA treatment. In total, we identified 29 phosphoproteins in response to SA treatment, belonging to 9 functional categories. Among the 29 phosphoproteins, 17 protein spots were common in both resistant (C101LAC) and susceptible (CO39) cultivars, suggesting that some physiological processes could be commonly influenced by SA. However, the other 12 phosphoproteins showed different changes between CO39 and C101LAC. The expression patterns of these phosphoproteins could help to understand the different molecular mechanism of SA-induced resistance against M. oryzae in the two rice cultivars.

Photosynthesis has been well-known to be highly sensitive to various exogenous stimulus, including SA, JA and some plant hormone [24, 25]. In this study, total 5 phosphoproteins related to photosynthesis were differentially regulated after SA treatment, in which chloroplast 23 kDa polypeptide of photosystem II (PSII, spots 34) was firstly reported to be phosphorylated which was validated by NanoLC-MS/MS. It has been previously demonstrated that photosynthesis could be regulated through phosphorylation of photosynthesis-related proteins. For example, phosphorylation of rubisco activase (RCA), known as an ancillary photosynthetic protein essential for Rubisco activity, was reduced following abiotic stress, resulting in the malfunctioning of photosynthesis [26]. Transketolase (TKL), a key enzyme linking the non-oxidative pentose phosphate pathway and the Calvin cycle, is presumed to participate in the functional regulation of numerous metabolic pathways by phosphorylation [27]. Previous studies have shown that dephosphorylation leads to a dramatic decrease in the TKL activity, which declines the rates of photosynthesis [28]. In this study, most of the phosphoproteins, including RCA (spots 5), Rubisco large subunit (Spots 20) and TKL (spots 27) were up-regulated in both rice cultivars at 12 h and 24 h post SA treatment. Consistent with our results, it has previously been reported that SA could alleviate decreases in plant photosynthesis against pathogen infection to produce various metabolites and energy, thus improving plant defense response [7]. However, we also noted that HAD-superfamily hydrolase (subfamily IA, variant 3 containing protein, spot 1) were down-regulated in C101LAC 12 h post SA treatment. It is suggested that SA regulates plant photosynthesis via phosphorylation of photosynthesis-related proteins and the effect of SA on the photosynthetic machinery is complex in rice plants.

We identified a putative cinnamoyl-CoA reductase (CCR, spot 37) as SA-responsive phosphoprotein. CCR catalyzes the first step of the monolignol pathway for lignin biosynthesis and therefore plays an essential role in defense-related processes in rice [29]. In our previous study, CCR was firstly reported to be dephosphorylated in susceptible rice plants after M. oryzae infection, which resulted in a decrease in enzyme activities [30]. In this study, CCR was notably up-regulated only in C101LAC at 24 h after SA treatment (Table 1). Correspondingly, CCR activity was significantly increased in the resistant cultivar at 24 h after SA treatment, but not in the susceptible cultivar (Fig. 3a). CCR, as an effector of small GTPase Rac in rice defense signaling, is considered to be activated by OsRac1, which further leads to efficient production of monolignols, deposition of lignin, increased ROS accumulation [31]. Taken together, we speculate that posphorylation of CCR may increase enzyme activities, accelerate lignin synthesis and ROS accumulation, thus further activate rice defense responses against infection of pathogens.

Phosphorylation of antioxidant enzymes has been reported to be involved in antioxidant defense [32]. In this study, two phosphoproteins were identified as GDP-mannose 3,5-epimerase 1 (OsGME 1, spot 15) and L-ascorbate peroxidase 1 (OsAPx01, spot 33), which are known as antioxidant enzymes. Phosphorylated OsGME 1 was down-regulated only in the susceptible cultivar at 12 h after SA treatment, and phosphorylated OsAPx01 was down-regulated in both cultivars. Previous studies have showed that dephosphorylation decreases the activities for APX and OsGME 1 upon biotic or abiotic stress in several plants [26, 33]. To validate phosphorylated regulation of enzyme activity, APX activity was performed (Fig. 3b). The results showed that APX activity was significantly decreased in both cultivars at 24 h after SA treatment, which were consistent with the phosphoproteomic analysis. It has been suggested that down-regulation of scavenging/antioxidant systems and/or activities may contribute to increase ROS accumulation [34]. Earlier study reported that plant cells presumably regulated ROS levels by coordinating activities of ROS-generating enzymes such as SOD and ROS-degrading enzymes such as APX, POX, and CAT [35]. A detailed comparison of the influence of SA on H2O2 production, SOD activities and H2O2-degrading enzymes showed that tobacco leaves treated with SA may have enhanced H2O2 largely by activating enzymes capable of generating H2O2, and by inactivating enzymes that are capable of degrading H2O2 [36]. In this study, the findings that phosphorylation of APX and OsGME 1 was down-regulated in rice suggests a significant decrease of the activities of ROS-scavenging antioxidant enzymes, further accelerates oxidative burst in both cultivars. Consistent with this notion, significant ROS accumulation was observed in both rice cultivars after SA treatment (Fig. 4), which indicated that SA could regulate antioxidant enzymes’ activity by phosphorylation, then enhance ROS accumulation in rice plants as a defense response. However, SA and ROS interactions are complicated and awaits further investigation in future to fully elucidate.

One phosphoprotein was identified as 60 kDa chaperonin beta (Cpn60β, spot 23), which is a molecular chaperone involved in protein destination and assembly. Molecular chaperones play regulatory roles in prevention of stress injury, immune response, and cell death in plants, which might help plant enhance defense system against pathogen infection [37]. Phosphorylation is a common regulatory mode for the function of different classes of molecular chaperones, which markedly improves the capacity to enhance its binding to unfolded proteins, facilitate rapid degradation of certain abnormal proteins, and protect against oxidative stress injury [38, 39]. In this study, phosphorylated Cpn60β was up-regulated at 24 h after SA treatment in both cultivars (Table 1). A previously proteomic study showed that 60 kDa chaperonin was significantly up-regulated in resistant tomato cultivar compared with susceptible tomato cultivar after bacterial infection, which might provide enhanced defense system against Pseudomonas solanacearum [40]. Taken together, we speculate that SA-induced phosphorylation of Cpn60β may improve the activities of Cpn60β, further contributing towards rice defense against pathogen attack.

Seven phosphoproteins were involved in carbohydrate metabolism, including phosphoglycerate mutase (PGAM, spot 9), glyceraldehyde-3-phosphate dehydrogenase (GAPDH, spot 17 and 18), alpha 1,4-glucan phosphorylase L isozyme (α-GP, spots 24 and 25), phosphoribulokinase (PRK, spot 31 and 32). PGAM catalyzes the conversion of 3-phosphoglycerate to 2-phosphoglycerate, which is a crucial step in glycolysis. Previous phosphoproteomics studies showed that PGAM can be phosphorylated at multiple tyrosine sites, which enhances PGAM activity and upregulates glycolysis [41, 42]. In this study, phosphorylated PGAM was down-regulated only in C101LAC 12 h after SA treatment. To validate regulation of enzyme activity by SA-induced protein phosphorylation, PGAM activity was assessed (Fig. 3b). Agreeable with the phosphoproteomic analysis, the results showed that a significant decrease in PGAM activity was noted in C101LAC 12 h after SA treatment, but little change was found in CO39. α-GP catalyzes the reversible phosphorolysis of glucan chains and releases glucose-1-phosphate for starch resynthesis. Phosphorylated α-GP is active, and dephosphorylated α-GP is inactive [43]. GAPDH catalyzes a critical step in glycolysis; GAPDH activity was significantly decreased if phosphorylated [44]. In this study, phosphorylated α-GP was up-regulated both in CO39 and C101LAC 12 h after SA treatment, suggesting that SA might increase α-GP activity. However, phosphorylated GAPDH was down-regulated in C101LAC, but up-regulated in CO39 12 h after SA treatment. The results above suggest that accumulation of carbohydrates was differently regulated in two different rice cultivars by SA induction, indicating that carbohydrate metabolism pathways may be regulated by reverse phosphorylation of important enzymes. However, little is known about the potential role for phosphorylated regulation in these enzymes of carbohydrate metabolism and further investigation is needed to decipher their role in rice.


In this study, we performed a comparative phosphoproteomic analysis to investigate the molecular mechanisms of SA-induced defense response in different rice cultivars. A total of 29 SA-responsive phosphoproteins were successfully identified by MAIDL-TOF/TOF analysis. Phosphoproteins involved in photosynthesis, antioxidative enzymes, molecular chaperones showed similar changes in the two cultivars, while phosphoproteins related to protein synthesis and degradation, defense, amino acid metabolism, and carbohydrate metabolism were differentially expressed. Furthermore, phosphorylation within four identified phosphoproteins was validated by NanoLC-MS/MS analysis, and phosphorylated regulation of three important enzymes (CCR, PGAM and APX) was verified by activity determination. To best of our knowledge, it is the first report to measure rice phosphoproteomic changes induced by SA, which may broaden our understanding of SA-responsive mechanisms in rice.



The following chemical reagents were used in this study: SA (Sigma-Aldrich, St. Louis, MO, USA); Immobiline™ DryStrip pH 4–7 NL, 18 cm and IPG buffer pH 4–7 (GE Healthcare, Uppsala, Sweden); Pro-Q Diamond phosphoprotein gel stain (Molecular Probes, Eugene, OR, USA); Glycine, 1.5 mol/L Tris-HCl buffer pH 8.8, 30% acrylamide/bis solution (37.5:1), and overlay agrose (Bio-Rad, Hercules, CA, USA).

Plant materials

Two rice near isogenic lines (Oryza sativa indica) were obtained from the International Rice Research Institute, including C101LAC carrying the resistance gene Pi-1 against M. oryzae, and background line CO39 carrying no known resistance gene. Rice seedlings were sprayed with 0.1 mM SA solution (containing 0.02% v/v Tween 20) at the four-leaf stage [7]. The fourth leaves were harvested at 12 and 24 h after SA treatment. Spraying with sterilized water containing 0.02% v/v Tween 20 served as blank control. The leaves were sampled by freezing in liquid nitrogen, and stored at − 80 °C before assessment.

Phosphoproteome enrichment, 2DE and gel analysis

Total proteins were extracted essentially from 5 g of rice leaf samples by using a PEG-mediated prefractionation method [45]. Enrichment of phosphorylated proteins follows the well-established Al (OH)3-MOAC method [30]. The protein content was determined using the coomassie blue dye-binding method, with BSA as the standard [46]. Three replicates with different pools of leaf samples were performed, and all the procedures were carried out at 4 °C.

2DE follows the established protocol [47]. Phosphoproteins were visualized with Pro-Q Diamond fluorescent gel stain according to the methods described previously [48], before imaging with a Typhoon Trio Variable Mode Imager (GE Healthcare, Uppsala Sweden). PDQuest software (Version 8.0, Bio-Rad, Hercules, CA, USA) was used for quantitative analysis with gel spots, including spot detection, measurement, matching and calculation. Each phosphoprotein sample was analyzed by 2DE for at least three times. The protein spots showing ≥1.5-fold increase or decrease in all three biological repeats were selected as putative differentially regulated phosphoproteins.

Phosphoprotein identification by MALDI-TOF/TOF MS

The differentially expressed phosphoprotein spots were manually excised from the gels, before in-gel digestion [7]. The peptides were subsequently analyzed using the ABI 4800 Proteomics Analyzer MALDI-TOF/TOF (Applied Biosystems, Foster City, CA). Database search was performed in the Oryza sativa database (Uniprot, v.2016.08.24) using the MASCOT search engine 2.2 (Matrix Science, Ltd.) with GPS-Explorer Software 3.6 (Applied Biosystems). The parameter settings were as following: peptide mass tolerance: 100 ppm; fragment tolerance: ±0.3 Da; protein score C.I.%: ≥95%; total ion score C.I.%: ≥95% and significance threshold: p < 0.05. Besides, to eliminate the redundancy of proteins that appeared in the database under different names and accession numbers, the single-protein member belonging to the species of O. sativa or others with the highest protein score (top rank) was singled out from the multi-protein family.

Identification of phosphorylation sites by NanoLC-MS/MS

Total proteins extraction from rice leaves was performed using a PEG-mediated prefractionation method [30], and protein was digested with trypsin following the FASP method [49]. The trypsin-digested peptide mixture was then loaded onto aliquot of titanium dioxide (TiO2) beads (5 μm Titansphere, GL Sciences, Japan), which were then collected by centrifugation after washing twice with 30 mg/mL DHB (2,5-dihydroxybenzoic acid) buffer. The beads were further washed for twice with 60% ACN/0.1% TFA and 0.1% TFA respectively, before elution with a 60% ACN/4% ammonium solution. NanoLC-MS/MS was performed with a Q Exactive MS (Thermo Finnigan) equipped with Easy nLC1000 (ThermoFisher, San Jose, CA). The peptide mixture was seperated on a C18-reversed phase column with a flow rate of 250 nL/min over 240 min. Peptides were analyzed by MS/MS in positive ion mode, and the MS/MS spectra search was performed against the Uniprot_Oryza database (v.2018.02.27), using Mascot 2.2 engine. Proteome Discoverer 1.3 (Thermo Electron, San Jose, CA) was used for identification of phosphorylation peptides, with the threshold setting as pRS score above 50 indicating a good PSM (Peptide Spectrum Matches) and pRS probabilities above 75% indicating a truly phosphorylated site.

Determination of enzyme activities and reactive oxygen species

The fourth leaves were sampled at 12 and 24 h after SA treatment. Enzymatic activity of APX (ascorbate peroxidase), CCR (cinnamoyl-CoA reductase) activity, and phosphoglycerate mutase (PGAM) was assayed respectively following the established protocols [50,51,52]. Hydrogen peroxide (H2O2) follows Brennan and Frenkel’s method [53]. The rate of superoxide (O2.−) production was measured based on nitroblue tetrazolium (NBT) reduction [54]. Determination of malondialdehyde (MDA) content follows the previously described method [55], using the following formula to calculate MDA content (C): C (μmol/L) = 6.45 (A532- A600)-0.56 A450.

Quantitative real-time PCR (qRT-PCR) analysis

Total RNA isolation from rice leaves was performed using Eastep@ Super Total RNA Extraction Kit (Promega, Shanghai, China). The first-strand cDNA was prepared from 2 μg of normalized total RNA, with FastKing RT Kit (with gDNase) (Tiangen Biotech, Beijing, China). Gene-specific primers used for selected gene transcription assessment were designed using the Primer 5.0 software and were listed in (Additional file 4: Table S4). The tubulin gene (Uniprot Accession No. Q58G87) was used as reference. qRT-PCR was conducted on a CFX Coxnnect™ Real-Time System (Bio-Rad, Hercules, CA, USA) with the iTaq™ SYBR® Green Supermix (Bio-Rad, Hercules, CA, USA) according to the manufacturer’s protocol. Three independent biological replicates were performed for each gene. Relative transcript levels for each gene were calculated by 2Ct method [56].

Statistical analysis

Means ± standard error (SE) was derived from three biological replicates. The data were analyzed using the one-way analyses of variance (ANOVA) and the significant differences were determined as p ≤ 0.05, by the Duncan’s test using SPSS software (version 19.0, SPSS Inc., Chicago, IL, USA).

Availability of data and materials

The data generated or analyzed during this study are included in this published article and its supplementary information files.



Two-dimensional gel electrophoresis


Ascorbate peroxidase


Cinnamoyl-CoA reductase


Isoelectric focusing


Immobiline pH gradien


Liquid chromatography


Metal oxide affinity chromatography


Mass spectrometer


Tandem mass spectrometry


Polyethylene glycol


Phosphoglycerate mutase

Pro-Q DPS:

Pro-Q diamond phosphoprotein stain


Post-translational modifications


Quantitative real-time PCR


Reactive oxygen species


Salicylic acid


  1. 1.

    Sasaya T, Nakazono-Nagaoka E, Saika H, Aoki H, Hiraguri A, Netsu O, Uehara-Ichiki T, Onuki M, Toki S, Saito K, Yatou O. Transgenic strategies to confer resistance against viruses in rice plants. Front in Microbiology. 2014;4:409.

    Article  Google Scholar 

  2. 2.

    Jones K, Kim DW, Park JS, Khang CH. Live-cell fluorescence imaging to investigate the dynamics of plant cell death during infection by the rice blast fungus Magnaporthe oryzae. BMC Plant Biol. 2016;16:69.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Ribot C, Hirsch J, Balzergue S, Tharreau D, Notteghem JL, Lebrun MH, Morel JB. Susceptibility of rice to the blast fungus, Magnaporthe grisea. J Plant Physiol. 2008;165(1):114–24.

    Article  CAS  Google Scholar 

  4. 4.

    Abed-Ashtiani F, Kadir JB, Selamat AB, Hanif BM, Nasehi A. Effect of foliar and root application of silicon against rice blast fungus in MR219 rice variety. Plant Pathol J. 2012;28(2):164–71.

    Article  CAS  Google Scholar 

  5. 5.

    RivasSan VM, Plasencia J. Salicylic acid beyond defence: its role in plant growth and development. J Exp Bot. 2011;62(10):3321–38.

    Article  CAS  Google Scholar 

  6. 6.

    Dempsey DA, Klessig DF. How does the multifaceted plant hormone salicylic acid combat disease in plants and are similar mechanisms utilized in humans? BMC Biol. 2017;15:23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Li YF, Zhang ZH, Nie YF, Zhang LH, Wang ZZ. Proteomic analysis of salicylic acid-induced resistance to Magnaporthe oryzae in susceptible and resistant rice. Proteomics. 2012;12(14):2340–54.

    Article  CAS  Google Scholar 

  8. 8.

    Lozano-Duran R, Robatzek S. 14–3-3 proteins in plant-pathogen interactions. Mol Plant-Microbe Interact. 2015;28(5):511–8.

    Article  CAS  Google Scholar 

  9. 9.

    Silva-Sanchez C, Li HY, Chen SX. Recent advances and challenges in plant phosphoproteomics. Proteomics. 2015;15(5–6):1127–41.

    Article  CAS  Google Scholar 

  10. 10.

    Wu LJ, Hu XL. Wang SX, Tian L, Pang YJ, Han ZP, Wu LC. Chen YH Quantitative analysis of changes in the phosphoproteome of maize induced by the plant hormone salicylic acid Scientific Reports. 2015;5:18155.

    CAS  Google Scholar 

  11. 11.

    Gupta R, Min CW, Meng Q, Agrawal GK, Rakwal R, Kim ST. Comparative phosphoproteome analysis upon ethylene and abscisic acid treatment in Glycine max leaves. Plant Physiol Biochem. 2018;130:173–80.

    Article  CAS  Google Scholar 

  12. 12.

    Hou YX, Qiu JH, Tong XH, Wei XJ, Nallamilli BR, Wu WH, Huang SW, Zhang JA. Comprehensive quantitative phosphoproteome analysis of rice in response to bacterial blight. BMC Plant Biol. 2015;15:163.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. 13.

    Qiu JH, Hou YX, Wang YF, Li ZY, Zhao J, Tong XH. Lin HY. Wei XJ, Ao H, Zhang JA. Comprehensive proteomic survey of ABA-induced protein phosphorylation in rice (Oryza sativa L.). International Journal of Molecular Sciences. 2017;18(1):60.

  14. 14.

    Gu YQ, Yang C, Thara VK, Zhou J, Martin GB. Pti4 is induced by ethylene and salicylic acid, and its product is phosphorylated by the Pto kinase. Plant Cell. 2000;12(5):771–86.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Ueno Y, Yoshida R, Kishikaboshi M, Matsushita A, Jiang CJ, Goto S, Takahashi A, Hirochika H, Takatsuji H. MAP kinases phosphorylate rice WRKY45. Plant Signal Behav. 2013;8(6):e24510.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Zhang Z, Chen J, Lin S, Li Z, Cheng R, Fang C, Chen H, Lin W. Proteomic and phosphoproteomic determination of ABA's effects on grain-filling of Oryza sativa L. inferior spikelets. Plant Science. 2012:185–186–259–73.

  17. 17.

    Naryzhny SN, Zgoda VG, Maynskova MA, Novikova SE, Ronzhina NL, Vakhrushev IV, Khryapova EV, Lisitsa AV, Tikhonova OV, Ponomarenko EA, Archakov AI. Combination of virtual and experimental 2DE together with ESI LC-MS/MS gives a clearer view about proteomes of human cells and plasma. Electrophoresis. 2016;37(2):302–9.

    Article  CAS  Google Scholar 

  18. 18.

    Kim Y, Mun BG, Khan AL, Waqas M, Kim HH, Shahzad R, Imran M, Yun BW, Lee IJ. Regulation of reactive oxygen and nitrogen species by salicylic acid in rice plants under salinity stress conditions. PLoS One. 2018;13(3):e0192650.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Sheteiwy MS, An J, Yin M, Jia X, Guan Y, He F, Hu J. Cold plasma treatment and exogenous salicylic acid priming enhances salinity tolerance of Oryza sativa seedlings. Protoplasma. 2019;256(1):79–99.

    Article  CAS  Google Scholar 

  20. 20.

    Pi Z, Zhao ML, Peng XJ, Shen SH. Phosphoproteomic analysis of Paper mulberry reveals phosphorylation functions in chilling tolerance. J Proteome Res. 2017;16(5):1944–61.

    Article  CAS  Google Scholar 

  21. 21.

    Pan DZ, Wang LX, Tan FL, Lu S, Lv XJ, Zaynab M, Cheng CL, Abubakar YS, Chen SP, Chen W. Phosphoproteomics unveils stable energy supply as key to flooding tolerance in Kandelia candel. J Proteome. 2018;176:1–12.

    Article  CAS  Google Scholar 

  22. 22.

    Kline-Jonakin KG, Barrett-Wilt GA, Sussman MR. Quantitative plant phosphoproteomics. Curr Opin Plant Biol. 2011;14:507–11.

    Article  CAS  Google Scholar 

  23. 23.

    Wolschin F, Wienkoop S, Weckwerth W. Enrichment of phosphorylated proteins and peptides from complex mixtures using metal oxide/hyroxide affinity chromatography (MOAC). Proteomics. 2005;5:4389–97.

    Article  CAS  Google Scholar 

  24. 24.

    Moya JL, Ros R, Picazo I. Heavy metal-hormone interactions in rice plants: effects on growth, net photosynthesis, and carbohydrate distribution. J Plant Growth Regul. 1995;14(2):61–7.

    Article  CAS  Google Scholar 

  25. 25.

    Nazar R, Umar S, Khan NA. Exogenous salicylic acid improves photosynthesis and growth through increase in ascorbate-glutathione metabolism and S assimilation in mustard under salt stress. Plant Signal Behav. 2015;10(3):e1003751.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Margaria P, Abbà S, Palmano S. Novel aspects of grapevine response to phytoplasma infection investigated by a proteomic and phosphoproteomic approach with data integration into functional networks. BMC Genomics. 2013;14:38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Rocha AG, Mehlmer N, Staez S, Mair A, Parvin N, Chigri F, Teige M, Vothknecht UC. Phosphorylation of Arabidopsis transketolase at Ser428 provides a potential paradigm for the metabolic control of chloroplast carbon metabolism. The Biochemical Journal. 2014;458(2):313–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Suzuki Y, Kondo E, Makino A. Effects of co-overexpression of the genes of Rubisco and transketolase on photosynthesis in rice. Photosynth Res. 2017;131(3):281–9.

    Article  CAS  Google Scholar 

  29. 29.

    Park HL, Bhoo SH, Kwon M, Lee SW, Cho MH. Biochemical and expression analyses of the rice cinnamoyl-CoA reductase gene family. Front Plant Sci. 2017;8:2099.

    Article  PubMed  PubMed Central  Google Scholar 

  30. 30.

    Li YF, Ye ZJ, Nie YF, Zhang J, Wang GL, Wang ZZ. Comparative phosphoproteome analysis of Magnaporthe oryzae-responsive proteins in susceptible and resistant rice cultivars. J Proteome. 2015;115:66–80.

    Article  CAS  Google Scholar 

  31. 31.

    Kawasaki T, Koita H, Nakatsubo T, Hasegawa K, Wakabayashi K, Takahashi H, Umemura K, Umezawa T, Shimamoto K. Cinnamoyl-CoA reductase, a key enzyme in lignin biosynthesis, is an effector of small GTPase Rac in defense signaling in rice. Proc Natl Acad Sci U S A. 2006;103(1):230–5.

    Article  CAS  Google Scholar 

  32. 32.

    Gou JY, Li K, Wu KT, Wang XD, Lin HQ, Cantu D, Uauy C, Dobon-Alonso A, Midorikawa T, Inoue K, Sanchez J, Fu D, Blechl A, Wallington E, Fahima T, Meeta M, Epstein L, Dubcovsky J. Wheat stripe rust resistance protein WKS1 reduces the ability of the thylakoid-associated ascorbate peroxidase to detoxify reactive oxygen species. Plant Cell. 2015;27(6):1755–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Han C, Wang K, Yang PF. Gel-based comparative phosphoproteomic analysis on rice embryo during germination. Plant & Cell Physiology. 2014;55(8):1376–94.

    Article  CAS  Google Scholar 

  34. 34.

    Xia XJ, Zhou YH, Shi K, Zhou J, Foyer CH, Yu JQ. Interplay between reactive oxygen species and hormones in the control of plant development and stress tolerance. J Exp Bot. 2015;66(10):2839–56.

    Article  CAS  Google Scholar 

  35. 35.

    Herrera-Vásquez A, Salinas P, Holuigue L. Salicylic acid and reactive oxygen species interplay in the transcriptional control of defense genes expression. Front Plant Sci. 2015;6:171.

    Article  PubMed  PubMed Central  Google Scholar 

  36. 36.

    Chen Z, Silva H, Klessig DF. Active oxygen species in the induction of plant systemic acquired resistance by salicylic acid. Science. 1993;262(5141):1883–6.

    Article  CAS  Google Scholar 

  37. 37.

    Ishikawa A, Tanaka H, Nakai M, Asahi T. Deletion of a chaperonin 60β gene leads to cell death in the Arabidopsis lesion initiation 1 mutant. Plant Cell Physiol. 2003;44(3):255–61.

    Article  CAS  Google Scholar 

  38. 38.

    Sherman MY, Goldberg AL. Heat shock in Escherichia coli alters the protein-binding properties of the chaperonin groEL by inducing its phosphorylation. Nature. 1992;357(6374):167–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Canova MJ, Kremer L, Molle V. The mycobacterium tuberculosis GroEL1 chaperone is a substrate of Ser/Thr protein kinases. J Bacteriol. 2009;191(8):2876–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. 40.

    Afroz A, Khan MR, Ahsan N, Komatsu S. Comparative proteomic analysis of bacterial wilt susceptible and resistant tomato cultivars. Peptides. 2009;30(9):1600–7.

    Article  CAS  Google Scholar 

  41. 41.

    Hitosugi T, Zhou L, Fan J, Elf S, Zhang L, Xie JX, Wang Y, Gu TL, Aleckovic M, LeRoy G, Kang YB, Kang HB, Seo JH, Shan CL, Jin P, Gong WM, Lonial S, Arellano ML, Khoury HJ, Chen GZ, Shin DM, Khuri FR, Boggon TJ, Kang S, He C, Chen J. Tyr26 phosphorylation of PGAM1 provides a metabolic advantage to tumours by stabilizing the active conformation. Nat Commun. 2013;4:1790.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. 42.

    Wang Y, Cai WS, Chen LN, Wang GY. Molecular dynamics simulation reveals how phosphorylation of tyrosine 26 of phosphoglycerate mutase 1 upregulates glycolysis and promotes tumor growth. Oncotarget. 2017;8(7):12093–107.

    PubMed  PubMed Central  Google Scholar 

  43. 43.

    Tan AW, Nuttall FQ. Characteristics of the dephosphoryated form of phosphorylase purified from rat liver and measurement of its activity in crude liver preparations. Biochim Biophys Acta. 1975;410(1):45–60.

    Article  CAS  Google Scholar 

  44. 44.

    Piattoni CV, Ferrero DML, Dellaferrera I, Vegetti A, Iglesias AÁ. Cytosolic glyceraldehyde-3-phosphate dehydrogenase is phosphorylated during seed development. Front Plant Sci. 2017;8:522.

    Article  PubMed  PubMed Central  Google Scholar 

  45. 45.

    Kim ST, Cho KS, Yu S, Kim SG, Hong JC, Han CD, Bae DW, Nam MH, Kang KY. Proteomic analysis of differentially expressed proteins induced by rice blast fungus and elicitor in suspension-cultured rice cells. Proteomics. 2003;3(12):2368–78.

    Article  CAS  Google Scholar 

  46. 46.

    Bradford MM. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Analytical Biochemistry. 1976;72:248–54.

    Article  CAS  Google Scholar 

  47. 47.

    Li YF, Nie YF, Zhang ZH, Ye ZJ, Zou XT, Zhang LH, Wang ZZ. Comparative proteomic analysis of methyl jasmonate-induced defense responses in different rice cultivars. Proteomics. 2014;14(9):1088–101.

    Article  CAS  Google Scholar 

  48. 48.

    Agrawal GK, Thelen JJ. Development of a simplified, economical polyacrylamide gel staining protocol for phosphoproteins. Proteomics. 2005;5(18):4684–8.

    Article  CAS  Google Scholar 

  49. 49.

    Wiśniewski JR, Zougman A, Nagaraj N, Mann M. Universal sample preparation method for proteome analysis. Nat Methods. 2009;6(5):359–62.

    Article  CAS  Google Scholar 

  50. 50.

    Sofo A, Scopa A, Nuzzaci M, Vitti A. Ascorbate peroxidase and catalase activities and their genetic regulation in plants subjected to drought and salinity stresses. Int J Mol Sci. 2015;16(6):13561–78.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. 51.

    Lüderitz T, Grisebach H. Enzymic synthesis of lignin precursors. Comparison of cinnamoyl-COA reductase and cinnamyl alcohol: NADP+ dehydrogenase from spruce (Picea abies L.) and soybean (Glycine max L.). Eur J Biochem. 1981;119(1):115–24.

    Article  Google Scholar 

  52. 52.

    Bourgis F, Botha FC, Mani S, Hiten FN, Rigden DJ, Verbruggen N. Characterization and functional investigation of an Arabidopsis cDNA encoding a homologue to the d-PGMase superfamily. J Exp Bot. 2005;56(414):1129–42.

    Article  CAS  Google Scholar 

  53. 53.

    Brennan T, Frenkel C. Involvement of hydrogen peroxide in the regulation of senescence in pear. Plant Physiol. 1977;59(3):411–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. 54.

    Shah J, Kachroo PK, Nandi A, Klessig DF. A recessive mutation in the Arabidopsis SSI2 gene confers SA- and NPR1-independent expression of PR genes and resistance against bacterial and oomycete pathogens. Plant J. 2001;25(5):563–74.

    Article  CAS  Google Scholar 

  55. 55.

    Li H. Principles and techniques of plant physiological biochemical experimental. Beijing: Higher Education Press; 2000. p. 196–7.

    Google Scholar 

  56. 56.

    Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method. Methods. 2001;25(4):402–8.

    Article  CAS  Google Scholar 

Download references


We thank Professor Jianchi Chen (United States Department of Agriculture-Agricultural Research Service, San Joaquin Valley Agricultural Sciences Center) and Professor Yizhen Deng (Guangdong Province Key Laboratory of Microbial Signals and Disease Control, South China Agricultural University) for critical reviewing the manuscript.


This work was supported by the National Natural Science Foundation of China (31671968), Natural Science Foundation of Guangdong Province (2015A030313406), Guangdong Science and Technology Program (2016A020210099 and 2016A020210098), and Guangzhou Science and Technology Program (201804010119). The funding body had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information




RS, SQ, TZ, ZW, HL, YL and YN conceived the study, participated in the design of the study and discussed the results critically. TZ and HL performed qRT-PCR analysis. RS, SQ, TZ performed the phosphoproteomic analysis and biochemical analysis. ZW, HL, YL and YN was responsible for NanoLC-MS/MS analysis. RS, YL and YN wrote and revised the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Yunfeng Li or Yanfang Nie.

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: Figure S1.

Representative 2DE patterns of phosphoproteins from rice leaves treated with MQ water (as the control) and SA. Figure S2 All additional 2DE gels of rice phosphoproteins were shown as replicate gels. Figure S3. Close-up views of the regions of 2DE gels showing all SA-responsive phosphoprotein spots in two rice cultivars. Figure S4. Quantitative analysis of the SA-responsive phosphoproteins in rice leaves. Figure S5. Identification of spot 11 by MALDI-TOF/TOF MS. Figure S6. The functional category distribution of the 29 SA-responsive phosphoproteins. Figure S7. The MS/MS spectra of representative phosphorylated peptides of ELLS*YEYDGDEVPIVAGSALK, corresponding to Elongation factor Tu (Q6ZI53).

Additional file 2: Table S1.

Phosphopeptides and their phosphorylation sites (marked by lowercase letters of the amino acid residue) identified by NanoLC-MS/MS analysis.

Additional file 3: Table S2.

Rice phosphoproteins identified by NanoLC-MS/MS analysis.

Additional file 4: Table S3.

Mapping of SA-responsive phosphoproteins with the NanoLC-MS/MS data. Table S4. Gene-specific primers designed for qRT-PCR.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Sun, R., Qin, S., Zhang, T. et al. Comparative phosphoproteomic analysis of blast resistant and susceptible rice cultivars in response to salicylic acid. BMC Plant Biol 19, 454 (2019).

Download citation


  • Salicylic acid
  • Rice
  • Phosphoproteome
  • Two-dimensional gel electrophoresis
  • Protein phosphorylation