In vivo phosphoproteome characterization reveals key starch granule-binding phosphoproteins involved in wheat water-deficit response
BMC Plant Biology volume 17, Article number: 168 (2017)
Drought stress during grain development causes significant yield loss in cereal production. The phosphorylated modification of starch granule-binding proteins (SGBPs) is an important mechanism regulating wheat starch biosynthesis. In this study, we performed the first proteomics and phosphoproteomics analyses of SGBPs in elite Chinese bread wheat (Triticum aestivum L.) cultivar Jingdong 17 under well-watered and water-stress conditions.
Water stress treatment caused significant reductions in spike grain numbers and weight, total starch and amylopectin content, and grain yield. Two-dimensional gel electrophoresis revealed that the quantity of SGBPs was reduced significantly by water-deficit treatment. Phosphoproteome characterization of SGBPs under water-deficit treatment demonstrated a reduced level of phosphorylation of main starch synthesis enzymes, particularly for granule-bound starch synthase (GBSS I), starch synthase II-a (SS II-a), and starch synthase III (SS III). Specifically, the Ser34 site of the GBSSI protein, the Tyr358 site of SS II-a, and the Ser837 site of SS III-a exhibited significant less phosphorylation under water-deficit treatment than well-watered treatment. Furthermore, the expression levels of several key genes related with starch biosynthesis detected by qRT-PCR were decreased significantly at 15 days post-anthesis under water-deficit treatment. Immunolocalization showed a clear movement of GBSS I from the periphery to the interior of starch granules during grain development, under both water-deficit and well-watered conditions.
Our results demonstrated that the reduction in gene expression or transcription level, protein expression and phosphorylation levels of starch biosynthesis related enzymes under water-deficit conditions is responsible for the significant decrease in total starch content and grain yield.
Wheat (Triticum aestivum L.) grain yield and quality are affected by many environmental factors. Among the external stresses, water-deficit is one of the most important stress factor, as it causes significant yield loss in cereal production [1, 2]. The reductions of grain yield are due mainly to reduced starch accumulation, with starch representing 65–75% of grain dry weight [3, 4]. Thus, water stress affects starch biosynthesis and the composition of grain starch [1, 5].
Starch consists of the glucose polymers amylose and amylopectin. Amylose, as a relatively linear molecule composing of (1–4)-linked units of d-glucopyranosyl, is produced by GBSS I . Amylopectin consists mainly of long chains of (1–4)-linked d-glucopyranosyl units with occasional branching (1–6) linkages, resulting in tandem linked clusters (~9–10 nm long each) [7, 8]. Four key enzymes, ADP-glucose pyrophosphorylase (AGPase), SS, starch-branching enzyme (SBE), and starch-debranching enzyme (DBE) are involved in amylopectin synthesis .
Based on its diameter, starch granules can be classified into three types: A-granules with 10–50 μm in diameter representing up to 70% of the volume and 10% of the total number of starch granules [10, 11]; B-granules with diameter 5–9 μm and accounting for about 30% of the volume and 90% of the total number of granules, and C-granules with less than 5 μm in diameter. In general, C-granules are difficult to isolate and quantify, leading to misclassification as B-granules [12, 13]. For B-granules, the amylose content is 5% less than in A-granules .
Protein phosphorylation is one of the most common post-translational modifications (PTMs) in vivo, which regulates different biological processes, such as transcription and translation, proliferation, cellular signaling, communication and differentiation . Phosphorylation modifications in eukaryotes mostly occur at the residues of serine (Ser), threonine (Thr), and tyrosine (Tyr).
To date, the phosphoproteins in starch are detected mainly by three methods: Pro-Q Diamond staining, phosphorylation isotope labeling by γ-32P-ATP in vitro, and liquid chromatography-mass spectrometry (LC-MS/MS) technology [16,17,18,19]. For example, GBSS I can be stained with a phosphoprotein-specific dye in maize (Zea mays L.) . These phosphorylated enzymes (proteins) binding starch granules labeled by γ-32P-ATP in vitro can improve amylase activity and amylose synthesis [16, 20, 21].
In a recent study in our laboratory, we identified phosphorylation sites of GBSS I and SS I by LC-MS/MS which play important roles in starch biosynthesis [18, 22]. In bread wheat, large-scale phosphoproteomic analyses of seeds and seedling leaves have been performed using phosphopeptide enrichment and LS-MS/MS [15, 23,24,25,26,27]. More recently, we performed the first comparative proteomics analysis of A- and B-type starch granule-associated proteins in wheat and Aegilops crassa through starch granule-binding protein (SGBP) extraction and two-dimensional gel electrophoresis (2-DE) maps . However, because the enrichment of starch granule-binding proteins with high purity in sufficient quantity from starch granules is difficult, characterization of the phosphoproteome of SGBPs in vivo during wheat grain development and under water-deficit has been limited.
In the present study, we reported the first phosphoproteome analysis of SGBPs in vivo in developing wheat grains subjected to water-deficit treatment. We determined the dynamic transcriptional expression profiles of the key genes associated with starch biosynthesis during grain development and under water-deficit treatment. Our results revealed the changes of SGBP phosphorylation level in vivo in response to drought stress, as well as starch biosynthesis, which provides new insights into the molecular mechanisms of the response to drought stress.
Plant materials, planting, and sampling
The elite Chinese winter wheat cultivar Jingdong 17 with good drought resistance was used as material and planted in winter at the experimental station of China Agricultural University during the 2014–2015 growing season from October to June of next year, with three biological replicates (each plot was 30 m2). The soil is glue test bed for low light loamy salinization soil, and the cultivation and management practices used were the same as the local wheat management practices. Two treatments were imposed at the jointing and flowering stages: well-watered with regular water irrigation (CK) and water-deficit without irrigation (DS). The average annual amount of sunshine at the experimental location is 2690 h, the average annual temperature is 12.5 °C, and average rainfall during the wheat growing season is approximately 150 mm. To harvest grain materials at different developmental stages, individual flowers were marked with colored tape at post-anthesis. Developing grains from the middle of spikes were harvested at 8, 10, 15, 20, 25, and 30 days post-anthesis (DPA) from 9:00–11:00 am, and stored at −80°C.
Measurement of soil water content, agronomic traits, and quality
The measurement of soil water content during different development stages, including the double-ridge, jointing, flowering, and maturity stages, followed the previously described soil drilling method of Guan et al. .
The physiological traits of flag leaves and relative water content at five time points after flowering under water-deficit and well-watered conditions were measured according to Larbi and Mekliche . The main phenotypic and agronomic traits of mature plants were recorded, including the number of spikelets, infertile spikelets, and kernels per spike. We also recorded data on grain weight per spike, weight of 1000 kernels, and grain yield. Statistical analyses were conducted with SPSS statistics software (version 17.0). Total starch and amylose content were measured with assay kits (Megazyme Int. Ireland, Ireland), according to the manufacturer’s protocol and Chen et al. . Farinograph parameters were obtained by analyzing 10-g samples using a Farinograph-E instrument (C.W. Brabender Instruments, Inc.), based on the American Association of Cereal Chemists Approved Method (2000) 54 21.
Purification of starch granules, ultrastructural observation by SEM, and extraction of SGBPs
Starch granules from three biological replicates at various grain development stages were separated and purified according to Chen et al. . Scanning electron microscopy (SEM) examination of starch granules was performed according to Chen et al.  with a Hitachi Model S-4700 scanning electron microscope. SGBPs were extracted and purified according to the previous reports. [19, 22, 30].
2-DE and protein identification by MALDI-TOF/TOF-MS
Two-dimensional gel electrophoresis (2-DE), image analysis, and matrix-assisted laser desorption/ionization time-of-flight/time-of-flight mass spectrometry (MALDI-TOF/TOF-MS) identification of extracted SGBPs in three biological replicates were performed according to Hao et al. . For the linear gradient IPG strip, each sample includes 50 μg SGPBs in 360 μL rehydration buffer. For each biological sample, the 2-DE experiments were repeated three times for error control. The gels were stained with CBB (R-250/G-250 = 4:1) and analyzed the expression quantity using ImageMaster 2D Platinum software version 7.0 (Amersham Biosciences).
The differentially expressed protein spots were excised and analyzed according to the recent reports [22, 30]. Statistical significance was determined according to the Protein Score C.I. % and Total Ion Score C.I. %, both above 95% and a significance threshold of p < 0.05 for the MS/MS .
Phosphopeptide enrichment using TiO2 microcolumns
The enrichment procedure for phosphopeptides from three biological replicates was performed as reported by Wu et al.  and Zhang et al. . SGBPs extracted at 20 DPA were reduced directly with dithiothreitol (DTT), alkylated with iodoacetamide, and, subsequently, digested with endoproteinase Lys-C and trypsin. TiO2 beads (GL Sciences, Tokyo, Japan) were incubated in 400 μL of loading buffer, containing 65% acetonitrile (ACN) and 2% trifluoroacetic acid (TFA), saturated with glutamic acid. A total of 3 mg of tryptic peptides from each biological replicate was dissolved in 600 μL of loading buffer, and then incubated with the appropriate amount of TiO2 beads. After washing with 600 μL of wash buffer (65% ACN, 0.1% TFA), the phosphopeptides were eluted twice with 300 μL of elution buffer (500 mM NH4OH, 60% ACN). The eluates were dried and reconstituted in 0.1% formic acid (FA) in H2O for MS analysis.
Phosphopeptide identification, phosphorylation site localization, and Bioinformatic analysis
The enriched phosphopeptides of each biological replicates were separated on a self-packed C18 reversed-phase column that was directly connected with a nanoelectrospray ion source in an LTQ-Orbitrap XL mass spectrometer (Thermo Fisher Scientific, America). The inner diameter of reversed-phase column is 75 μm and 150 mm length (Column Technology Innovation CTI, Fremont, CA). The pump flow was split to achieve a flow rate of 1 μL/min for sample loading and 300 nL/min for MS analysis. The mobile phases consisted of A (0.1% FA) and B (0.1% FA/80% ACN). A five-step linear gradient was applied (from 5% to 30% B in 105 min, 35% to 90% B in 16 min, 90% B in 4 min, 90% to 2% B for 0.5 min, and 2% B for 14.5 min). The spray voltage was set to 2.0 kV, and the temperature of the heated capillary was set to 240 °C. Every MS scan was acquired at a resolution of 60,000 (at 400 m/z) with the lock-mass option enabled, followed by top 10 data-dependent MS/MS scans performed using collision induced dissociation (CID). We set 500 as the threshold for precursor ion selection, and the mass window for precursor ion selection was 2.0 Da. The dynamic exclusion duration was 120 s, the repeat count and duration was 1 and 30s, respectively. The analyzer for the MS scans was Orbitrap, and the analyzer for the MS/MS scans was LTQ (37% relative collision energy). Three biological replicates were performed independently from sample collection to phosphopeptide identification. Raw files were processed using MaxQuant (version 22.214.171.124) software, according to Cox et al. , and then searched against the wheat database (77,037 entries).
Phosphorylation residue localization was evaluated based upon the post-translational modification (PTM) scores. Potential phosphorylation residues were then grouped into three categories based upon Schwartz et al.  and Sugiyama et al. .
Phosphopeptides that met the following conditions were considered as having undergone a significant change in phosphorylation level, based on the method of Lv et al. : (1) phosphopeptide detected in all three biological replicates; (2) phosphopeptides with credible Student’s t-test (p < 0.05), (3) phosphorylation localization probability ≥0.75, and (4) phosphorylation site score differences ≥5.
The motif-X algorithm (http://motif-x.med.harvard.edu/) was used to extract significant enriched phosphorylation motif from phosphopeptides . A Phyre2 online server was used to predict the 3D structures of the proteins  and the 3D structures and the phosphorylated sites were displayed using Swiss-PdbViewer (SPDBV version 4.1) software .
Antibody development, western blotting, and immunolabeling
Polyclonal antibodies were produced in rabbits derived from the N-terminal sequences of T. aestivum SBE I, and T. aestivum SBE IIa . Monoclonal antibody anti-wheat GBSS I, SS I, and SS II antisera were prepared as described by Li et al. . After staining, the expected protein bands were collected and digested with trypsin, and then identified by MALDI-TOF/TOF-MS as described above. The identified proteins were further confirmed by Western blotting according to Li et al. .
Immunolabeling experiments were performed with samples of 15 DPA grains, using the recently described method by Chen et al. . The primary antibodies of GBSS I were diluted to 1:3000. In the control experiment, the primary antibody was omitted to test for nonspecific secondary antibody binding [42, 43].
RNA extraction, cDNA synthesis, and real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR)
Developing grains from the central parts of the spikes were harvested at six developmental stages of Jingdong 17 (8, 10, 15, 20, 25, and 30 DPA). Total RNA of individual samples was extracted with TRIzol Reagent (Invitrogen) according to Li et al. . cDNA were synthesised based on Chen et al. . Transcriptional expression patterns of the selected genes were detected using qRT-PCR according to Chen et al.  with minor modifications. The genes and the candidate reference genes were listed in Additional file 1: Table S1.
The strategy and experimental design of this study are presented in Additional file 2: Figure S1. Developing Jingdong 17 grains were harvested at 25 DPA from plants grown under water-deficit and well-watered conditions. Starch granules were then separated and purified. SGBPs were extracted and separated by 2-DE, and identified by MALDI-TOF/TOF-MS. Subcellular localization of key SGBPs was determined using immunogold labeling. Subsequently, SGBPs were subjected to large-scale phosphoproteome analysis with TiO2 microcolumns enrichment and LC-MS/MS, and we studied phosphorylated SGBPs that regulate starch biosynthesis and drought stress responses. We determined the dynamic transcriptional expression profiles of key genes related to starch biosynthesis during grain development under the two conditions.
Dynamic changes in soil water and leaf water content
Grain developmental stages and soil water content are shown in Fig. 1. The seeds of plants subjected to water-deficit treatment ripened faster than the seeds under well-watered conditions (Fig. 1a). Soil water content increased gradually with depth of the soil layers, from 0 to 200 cm. In the 0–100 cm soil layer, the water content ranged from 7% to 8%, 9% to 14%, under water-deficit condition and well-watered condition, respectively. Soil water content was approximately 20%–30% lower under water-deficit than under well-watered condition as shown in Fig. 1b. Thus, according to Guttieri et al. , the effect of the water-deficit corresponded to medium-to-severe drought stress. The medium-to-severe drought stress resulted in significantly reduced relative water content of the flag leaves (Fig. 1c).
Effects of water-deficit treatment on yield traits and starch composition
Changes in main agronomic traits and starch contents in response to control and water-deficit treatments are shown in Table 1. These results indicated that some important agronomic traits related to yield were affected significantly by water-deficit. In particular, the 1000-kernel weight and grain yield were reduced significantly under water-deficit condition. Total starch and amylopectin contents were also reduced significantly. These results indicated that the water-deficit treatment could cause significant yield losses and premature grain development, at the key jointing and flowering stages.
We determined changes in the main farinograph parameters in response to the control and water-deficit treatments, (Table 2). As shown in Table 2, important farinograph parameters, including peak viscosity, minimum viscosity, attenuation, final viscosity, and retrograde values, were affected significantly by water deficit. But the parameters of attenuation, peak time and pasting temperature did not change significantly under water deficit in this study. And our results indicated that drought stress affected starch quality.
Comparative proteome analysis of SGBPs under water-deficit conditions
Starch granules were separated and purified, and then observed using SEM. Microscope images showed that the granules had no protein or cell debris contamination (Additional file 2: Fig. S2). After removing the surface-associated proteins by washing three times, SDS-PAGE analysis further confirmed the complete removal of all contamination (Additional file 2: Fig. S3). Moreover, SDS extraction and sonication did not change the visual appearance of starch granules observed by SEM.
A 2-DE separation was performed using 60 μg of SGBPs extracted from developing grains at 25 DPA. The SGBPs displayed similar proteome profiles under the two treatments (Fig. 2a). In total, we identified 33 protein spots containing SS I, SS II-a, and GBSS I, of which 12 spots were identified as GBSS I. A total of five and eight spots were identified as SS I and SS II-a, respectively. The remaining spots were identified as GBSS I partial (Additional file 1: Tables S2 and S3). Further analysis showed that the GBSS I, SS I, and SS II-a contents accounted for 71 and 71.9% of total SGBPs for the well-watered and water-deficit treatments, respectively. However, the respective proportions showed greater differences. For the well-watered treatment, the waxy protein (GBSS I) accounted for approximately 48.9% of the total SGBP content, while, under drought stress, it increased slightly to 52.4%. SS II-a accounted for 11.3% under the well-watered treatment, and its content increased to 12.5% in response to water-deficit. SS I accounted for 10.8% in well-watered plots, but it decreased significantly to 7.0% for the water-deficit treatment (Fig. 2b).
Phosphoproteome characterization of SGBPs under well-watered and water-deficit conditions
SGBPs (5 mg) were extracted from purified starch granules at 25 DPA (three biological replicates) and the phosphopeptides were enriched and analyzed by LC-MS/MS. In total, we identified 76 phosphopeptides containing 80 phosphorylated sites, representing 66 phosphoproteins (Fig. 3a–c). Under well-watered conditions, 67 phosphopeptides containing 71 phosphorylated sites and 57 phosphoproteins were identified (Additional file 1: Table S4). A total of 69 phosphopeptides with 73 phosphorylated sites were identified in 59 phosphoproteins under water-deficit conditions (Additional file 1: Table S4). Among the 76 phosphopeptides, 74 included one phosphorylated site, two had two phosphorylated sites, and one had more than three phosphorylated sites (Fig. 3d). In addition, among the 80 phosphorylated sites, 65 serines, 11 threonines, and 4 tyrosines were phosphorylated (Fig. 3e). To improve the accuracy of the dataset, only phosphorylation sites with p ≥ 0.75 were used in subsequent analyses. All mass spectrometry proteomics data obtained in this study (Additional file 1: Table S4) were deposited in the ProteomeXchange Consortium with the dataset identifier PXD002805 .
The identified phosphoproteins were mainly divided into seven categories: starch synthesis enzymes, protein kinase/phosphatase, ATPase, ubiquitin protein, transporter proteins, transcription/translation factors, and others (Additional file 1: Table S4). Starch synthesis enzymes included GBSS I, SS II-a, and SS III. The phosphorylation site distributions and 3-D structures of GBSS I and SS II-a are shown in Fig. 4. Three phosphorylated sites were detected under the water-deficit and well-watered conditions: Ser34, Ser450, and Thr323 in GBSS I (Fig. 4a, b) and Ser228, Ser776, and Tyr358 in SS II-a (Fig. 4c). It is worth noting that the degrees of phosphorylation at the Ser34 site in waxy protein, Tyr358 site in SS II-a, and the Ser837 site in SS III-a were reduced significantly in response to water-deficit conditions (Table 3). We also found that the degrees of phosphorylation of two kinase proteins (phosphatidylinositol 4-kinase type 2-beta and hydroxysteroid 11-beta-dehydrogenase) were reduced significantly in response to drought stress. The extent of phosphorylation at two sites in globulin 3 and one site in photosynthesis-related proteins PSII D2 increased significantly under the water-deficit treatment (Table 3).
Starch granule-binding proteins verification and western blotting
The SGBPs from 25 DPA were fractionated by SDS-PAGE and silver-stained (Fig. 5a). The isolated GBSS I and SS II-a proteins were identified by MALDI-TOF/TOF-MS (Additional file 1: Table S5). Moreover, monoclonal antibodies against GBSS I and SS II-a demonstrated their high specificity (Fig. 5b, c).
Subcellular localization of GBSS I
The subcellular localization of GBSS I was determined by immunolocalization. As is shown in Fig. 6, the hybridization signal of GBSS I was present mainly around the starch granules, under both conditions at 15 DPA, while most signals were located within starch granules at 25 DPA (Fig. 6).
Dynamic transcriptional expression profiling of starch synthesis-related genes during grain development
We determined the dynamic expression profiles of 16 starch synthesis-related genes, during six grain developmental stages in Jingdong 17 under both conditions by qRT-PCR analysis (Fig. 7). In general, all starch synthesis-related genes displayed a similar expression trend during grain development under both conditions. Most of the genes, such as GBSS I, SS I, SS III, SBE I, SBE II-a, SBE II-b, Pullulanase (PUL), isoamylase I (ISA I), isoamylase II (ISA II), adenosine 5 diphosphate glucose pyrophosphorylase small subunit II (AGPS II) and adenosine 5 diphosphate glucose pyrophosphorylase large subunit I (AGPL I), exhibited an up-downregulated expression trend and a higher expression level at the mid stages of grain development. And some genes (SS II-a, PGM, AGPS I) exhibited a downregulated expression trend. However, some genes exhibited expression differences between well-watered and water-deficit conditions. For example, SS III, SBE II-b and AGPLI genes were downregulated significantly at 15 DPA under water-deficit condition (Fig. 7).
Previous studies showed that drought stress severely affects wheat growth and development, especially at the grain filling stage, which generally results in a significant grain yield reduction. Photosynthesis is weakened in response to drought stress, and the accumulation of dry matter is reduced. Earlier studies showed that the grain filling rate and spike grain weight are reduced when drought stress is imposed at grain filling [45,46,47]. Continuous drought results in significant loss of total dry matter and wheat grain yield . Similarly, our study showed that water-deficit contributed to a marked decline in wheat leaf water content, spike grain number, 1000-kernel weight, and grain yield (Table 1). Wheat endosperm consists of starch, which serves as a major contributor to grain yield. And our results showed that total starch and amylopectin contents were reduced significantly under water-deficit conditions (Table 1). And the previous studies had shown that drought stress could decrease the number of endosperm cells and influence the formation of amyloids at early stages of seed growth, which subsequently causes a decline in endosperm starch accumulation [49, 50]. The reduction in endosperm starch accumulation is the primary cause of the grain weight reduction under drought conditions . Our study indicated that the significant reduction in amylopectins, which account for 70% of the starch content, could contribute to a marked reduction in total starch content and grain yield (Table 1).
As is known, starch is composed of amylopectins and amyloses, of which amylose synthesis is controlled by GBSS I, and amylopectins are synthesized by concerted reactions catalyzed by four enzyme classes: AGPase, SS, SBE, and DBE . In this study, our proteomic results revealed that the amount of the protein SSI was compromised significantly under water-deficit conditions compared to well-watered conditions. While the amylose content increased sightly after water-deficit treatment. GBSS I expression was increased during all stages, with the exception of 15 DPA after drought stress treatment according to the results of qRT-PCR shown in Fig. 7. And the quantity of the GBSS I protein was also increased, but not significantly for the water-deficit treatment group (Fig. 2). Meanwhile, immunolocalization results suggest that GBSS I had a clear movement from the periphery to the interior of starch granules during grain development (Fig. 6).
The synthesis of amylopectins are mainly controlled and catalyzed by starch synthase (SS). And the SS group consists of four isoforms (SS-I, SS-II, SS-III, and SS-IV), which are localized at the granule surface predominantly [52,53,54,55]. In this study, we found that the content of SS I was reduced significantly in response to the water-deficit treatment (Fig. 2), which finally resulted in significant reduction of amylopectin and total starch content (Table 1). Studies involving Arabidopsis thaliana L. and Oryza sativa L. demonstrated that SSI is required for the elongation of short A-chains within amylopectin [56, 57]. The function of SS-II is elongation of amylopectin chains of DP 6–10 to produce intermediate-length chains of DP 12–25 . Similar to the GBSS I gene, the transcriptional expression changes of several amylopectin synthesis-related genes such as SS III, SBE II-b and AGPL I were reduced significantly at 15 DPA (Fig. 7), suggesting that the effects of drought stress on starch biosynthesis occur at the transcriptional level, protein translation and posttranslational phosphorylation modification levels.
Phosphorylation could improve the activity of the SGBPs, thereby increasing the rate of starch synthesis and increasing total starch content [20, 21]. In a previous study, we used LC-MS/MS to identify phosphorylation sites present in GBSS I and SS I [18, 19]. In this study, we performed large-scale phosphoproteomic analysis of starch granule-binding proteins in Jingdong 17 in response to well-watered and water-deficit treatments. And our results indicated that the phosphorylation levels of key starch synthesis-related proteins were reduced significantly by the water-deficit treatment. Four types of enzyme related to starch synthesis, including GBSS I, SS II-a, SS III, and AGPS, were phosphorylated, and the phosphorylation sites were highly conserved [18, 19]. However, there were some differences in the degree of phosphorylation. For example, the Ser34 site of GBSSI, the Tyr358 site of SS II-a, and the Ser837 site of SS III-a had lower degrees of phosphorylation under drought conditions than under the well-watered treatment (Table 3). Based on previous study, the phosphorylation sites in starch synthase might play important roles in recognizing and attracting glucan substrates . A study on the 14–3-3 protein showed that phosphorylation could increase the affinity between the interacting proteins . For the GBSSI, reduced phosphorylation at Ser34 might affect enzyme activity, influence amylase synthesis, and reduce amylose content, for the reason that phosphorylation improves amylase activity and increases amylose synthesis . Similarly, for amylopectin, reduced phosphorylation at Tyr358 of SS II-a could also affect enzyme activity, influence amylopectin synthesis, and reduce amylopectin content. The obvious reduction of amylopectin contents would inevitably reduce total starch content (Table 1). Thus, our results showed that SGBP phosphorylation plays a critical role in starch biosynthesis.
Based on our results and previous reports, we propose putative pathways of amylose and amylopectin biosynthesis in response to water-deficit conditions (Fig. 8). The expression of starch biosynthesis-related proteins, such as GBSS I and SS II-a showed tiny increase, SSI was decreased significantly when subjected to the water-deficit treatment (Fig. 2). The phosphorylation levels of three key starch biosynthesis-related proteins were also reduced markedly by the water-deficit treatment as shown in Table 3, which could affect the activity of starch biosynthesis enzymes and their interaction, and result in reduced total starch content and grain yield. This also showed that the transcription levels of these genes are generally consistent with their phosphorylation levels but a little different with protein levels.
Water-deficit treatment significantly reduced spike grain numbers and weight, total starch and amylopectin content, and grain yield. Proteome analyses found that the expression quantity of SS II-a was reduced markedly when subjected to water deficit stress. Phosphoproteome analysis of SGBPs found that the phosphorylation levels of the main starch synthesis enzymes such as GBSS I, SS II-a, and SS III were decreased significantly under water-deficit. Furthermore, several key genes related with starch biosynthesis such as SS III, SBE II-b and AGPL I were downregulated significantly at 15 DPA. Our results demonstrated that the decrease in transcription expression, protein expression and phosphorylation levels of starch biosynthesis-related enzymes under water-deficit condition are responsible for the reduction in total starch content and grain yield.
- AGPL I:
Adenosine 5 diphosphate glucose pyrophosphorylase large subunit I;
- AGPS II:
Adenosine 5 diphosphate glucose pyrophosphorylase small subunit II
Collision induced dissociation
- GBSS I:
Granule-bound starch synthase
Liquid chromatography coupled to tandem mass spectrometry
Matrix-assisted laser desorption/ionisation time-of-flight tandem mass spectrometry
National Center for Biotechnology Information
Plant protein phosphorylation data base
Real-time quantitative reverse transcription polymerase chain reaction
Scanning electronic microscope
Starch granule-binding proteins
Ahmadi A, Baker DA. The effect of water stress on the activities of key regulatory enzymes of the sucrose to starch pathway in wheat. Plant Growth Regul. 2001;35(1):81–91. doi:10.1023/A:1013827600528.
Bartels D, Sunkar R. Drought and salt tolerance in plants. Crit Rev Plant Sci. 2005;24(1):23–58. doi.org/10.1080/07352680590910410
Duffus CM. Control of starch biosynthesis in developing cereal grains. Biochem Soc Tran. 1992;20(1):13–8. doi:10.1042/bst0200013.
Hurkman WJ, McCue KF, Altenbach SB, Korn A, Tanaka CK, Kothari KM, et al. Effect of temperature on expression of genes encoding enzymes for starch biosynthesis in developing wheat endosperm. Plant Sci. 2003;164(5):873–81. doi.org/10.1016/S0168-9452(03)00076-1
Jenner CF, Rathjen AJ. Factors regulating the accumulation of starch in ripening wheat grains. Functional Plant Bio. 1975;2(3):311–22. doi.org/10.1071/ PP9750311
Vrinten P, Nakamura T. Wheat granule-bound starch synthase I and II are encoded by separate genes that are expressed in different tissues. Plant Physiol. 2000;122(1):255–63.
Nakarnura T, Yamamori M, Hirano H, Hidaka S. Identification of three Wx proteins in wheat (Triticum aestivum L.). Biochem Genet. 1993;31(1):75–86. doi:10.1007/BF02399821.
Miura H, Sugawara A. Dosage effects of the three Wx genes on amylase synthesis in wheat endosperm. Theor Appl Genet. 1996;93(7):1066–70. doi:10.1007/s001220050336.
Whistler RL, Daniel JR. Molecular structure of starch. In: Whistler RL, Bemiller JN, Paschall EF, editors. Starch: chemistry and technology. Ed 2 ed. New York: Academic; 1984. p. 153–82.
Hughes CE, Briarty LG. Stereological analysis of contribution made to mature wheat endosperm starch by large and small granules. Starch-Stärke. 1976;28(10):336–7. doi:10.1002/star.19760281003.
Langeveld SMJ, Wijk RV, Stuurman N, Kijne JW, Pater SD. B-Granule granule containing protrusions and interconnections between amyloplasts in developing wheat endosperm revealed by transmission electron microscopy and GFP expression. J Exp Bot. 2000;51(349):1357–61. doi:10.1093/jexbot/51.349.1357.
Bechtel DB, Wilson JD. Amyloplast formation and starch granule development in hard red winter wheat. Cereal Chem. 2003;80(2):175–83. doi:10.1094/CCHEM.2003.80.2.175.
Wilson JD, Bechtel DB, Todd TC, Seib PA. Measurement of wheat starch granule size distribution using image analysis and laser diffraction technology. Cereal Chem. 2006;83(3):259–68. doi:10.1094/CC-83-0259.
Ao ZH, Jane JL. Characterization and modeling of the A- and B-granule starches of wheat, triticale, and barley. Carbohydr Polym. 2007;67(1):46–55. doi:10.1016/j.carbpol.2006.04.013.
Thingholm TE, Jensen ON, Larsen MR. Analytical strategies for phosphoproteomics. Proteomics. 2009;9(6):1451–68. doi:10.1002/pmic. 200800454.
Tetlow IJ, Beisel KG, Cameron S, Makhmoudova A, Liu F, Bresolin NS, et al. Analysis of protein complexes in wheat amyloplasts reveals functional interactions among starch biosynthetic enzymes. Plant Physiol. 2008;146(4):1878–91. doi:10.1104/pp.108.116244.
Grimaud F, Rogniaux H, James MG, Myers AM, Planchot V. Proteome and phosphoproteome analysis of starch granule-associated proteins from normal maize and mutants affected in starch biosynthesis. J Exp Bot. 2008;59(12):3395–406. doi:10.1093/jxb/ern198.
Chen GX, Zhu JT, Zhou JW, Subburaj S, Zhang M, Han CX, et al. Dynamic development of starch granules and the regulation of starch biosynthesis in Brachypodium Distachyon: comparison with common wheat and Aegilops peregrina. BMC Plant Biol. 2014;14(1):198–213. doi:10.1186/s12870- 014-0198-2.
Cao H, Yan X, Chen GX, Zhou JW, Li XH, Ma WJ, et al. Comparative proteome analysis of A- and B-type starch granule-associated proteins in bread wheat (Triticum aestivum L.) and Aegilops crassa. J. Proteomics. 2015;112:95–112. doi.org/10.1016/j.jprot.2014.08.002
Tetlow IJ, Morell MK, Emes MJ. Recent developments in understanding the regulation of starch metabolism in higher plants. J Exp Bot. 2004;55(406):2131–45. doi:10.1093/jxb/erh248.
Tetlow IJ, Wait R, ZX L, Akkasaeng R, Bowsher CG, Esposito S, et al. Protein phosphorylation in amyloplasts regulates starch branching enzyme activity and protein-protein interactions. Plant Cell. 2004;16(3):694–708. doi:10.1105/ tpc.017400.
Chen GX, Zhou JW, Liu YL, XB L, Han CX, Zhang WY, et al. Biosynthesis and regulation of wheat amylose and amylopectin from proteomic and phosphoproteomic characterization of granule binding proteins. Sci Rep. 2016;6:33111. doi:10.1038/srep33111.
Lv DW, Subburaj S, Ma C, Yan X, Li X, Appels R, et al. Proteome and phosphoproteome characterization reveals new response and defense mechanisms of Brachypodium distachyon leaves under salt stress. Mol Cell Proteomics. 2014a;13(2):632–52. doi:10.1074/mcp.M113.030171.
Lv DW, Ge P, Zhang M, Cheng ZW, Li XH, Yan YM. Integrative network analysis of the signaling cascades in seedling leaves of bread wheat by large-scale phosphoproteomic profiling. J Proteome Res. 2014b;13(5):2381–95. doi:10.1021/pr401184v.
Lv DW, Li X, Zhang M, AQ G, Zhen SM, Wang C, et al. Large-scale phosphoproteome analysis in seedling leaves of Brachypodium distachyon L. BMC Genomics. 2014c;15(1):375. doi:10.1186/1471-2164-15-375.
Zhang M, Ma CY, Lv DW, Zhen SM, Li XH, Yan YM. Comparative phosphoproteome analysis of the developing grains in bread wheat (Triticum aestivum L.) under well-watered and water deficit conditions. J Proteome Res. 2014;13(10):4281–97. doi:10.1021/pr500400t.
Zhang M, Lv DW, Ge P, Bian YW, Chen GX, Zhu GR, et al. Phosphoproteome analysis reveals new drought response and defense mechanisms of seedling leaves in bread wheat (Triticum aestivum L.). J Proteome. 2014;109:290–308. doi:10.1016/j.jprot.2014.07.010.
Guan JH, Zhang YP, Jiang A. Effect of different irrigation on water consumption characteristics and yield of spring wheat. Chin Agric Sci Bull. 2009;25(8):272–6.
Larbi A, Mekliche A. Relative water content (RWC) and leaf senescence as screening tool for drought tolerance in wheat. Options Méditerranéennes: Seŕie A Seminaires Méditerranéens. 2004;60:193–206.
Bancel E, Rogniaux H, Debiton C, Chambon C, Branlard G. Extraction and proteome analysis of starch granule-associated proteins in mature wheat kernel (Triticum aestivum L.). J Proteome Res. 2010;9(6):3299–310. doi:10.1021/pr9010525.
Hao PC, Zhu JT, AQ G, Lv DW, Ge P, Chen GX, et al. An integrative proteome analysis of different seedling organs in tolerant and sensitive wheat cultivars under drought stress and recovery. Proteomics. 2015;15(9):1544–63. doi:10.1002/pmic.201400179.
Ge P, Ma C, Wang S, Gao L, Li X, Guo G, et al. Comparative proteomic analysis of grain development in two spring wheat varieties under drought stress. Anal Bioanal Chem. 2012;402(3):1297–313. doi:10.1007/s00216-011-5532-z.
Wu J, Shakey Q, Liu W, Schuller A, Follettie MT. Global profiling of phosphopeptides by titania affinity enrichment. J Proteome Res. 2007;6(12):4684–9. doi:10.1021/pr070481m.
Cox J, Mann M. MaxQuant enables high peptide identification rates, individualized ppb-range mass accuracies and proteome-wide protein quantification. Nat Biotechnol. 2008;26(12):1367–72. doi:10.1038/nbt.1511.
Schwartz D, Gygi SP. An iterative statistical approach to the identification of protein phosphorylation motifs from large-scale datasets. Nat Biotechnol. 2005;23(11):1391–8. doi:10.1038/nbt1146.
Sugiyama N, Masuda T, Shinoda K, Nakamura A, Tomita M, Ishihama Y. Phosphopeptide enrichment by aliphatic hydroxyl acid-modified metal oxide chromatography for nano-LC-MS/MS in proteomics applications. Mol Cell Proteomics. 2007;6(6):1103–9. doi:10.1074/mcp.T600060-MCP200.
Kelley LA, Sternberg MJ. Protein structure prediction on the web: a case study using the Phyre server. Nat Protoc. 2009;4(3):363–71. doi:10.1038/ nprot.2009.2.
Guex N, Peitsch MC. SWISS-MODEL and the Swiss-Pdb viewer: an environment for comparative protein modeling. Electrophoresis. 1997;18(15):2714–23. doi:10.1002/elps.1150181505.
Vizcaíno JA, Cote RG, Csordas A, Dianes JA, Fabregat A, Foster JM, et al. The proteomics identifications (PRIDE) database and associated tools: status in 2013. Nucleic Acids Res. 2013;41:1063–9. doi:10.1093/nar/gks1262.
Rahman S, Regina A, Li Z, Mukai Y, Yamamoto M, Kosar-Hashemi B, et al. Comparison of starch-branching enzyme genes reveals evolutionary relationships among isoforms. Characterization of a gene for starch-branching enzyme IIa from wheat D genome donor Aegilops tauschii. Plant Physiol. 2001;125(3):1314–24. doi:10.1104/pp.125.3.1314.
Li XH, Ma WJ, Gao LY, Zhang YZ, Wang AL, Ji KM, et al. A novel chimeric LMW-GS gene from the wild relatives of wheat Ae. kotschyi and Ae. juvenalis: evolution at the Glu-3 loci. Genetics. 2008;180(1):93–101. doi:10.1534/ genetics.108.092403.
Danscher G. Localization of gold in biological tissue. A photochemical method for light and electron microscopy Histochem J. 1981;71(1):81–8. doi:10.1007/BF00592572.
Steinbrecht RA, Ozaki M, Ziegelberger G. Immunocytochemical localization of pheromone-binding protein in moth antennae. Cell Tissue Res. 1992;270(2):287–302. doi:10.1007/BF00328015.
Guttieri MJ, Stark JC, O’Brien K, Souza E. Relative sensitivity of spring wheat grain yield and quality parameters to moisture deficit. Crop Sci. 2001;41(2):327–35. doi:10.2135/cropsci2001.412327x.
Jamieson PD, Martin RJ, Francis GS. Drought influences on grain yield of barley, wheat, and maize. New Zealand J Crop Hort Sci. 1995;23(1):55–66. doi.org/10.1080/01140671.1995.9513868
Zhang J, Sui X, Li B, Su B, Li J, Zhou D. An improved water-use efficiency for winter wheat grown under reduced irrigation. Field Crops Res. 1998;59(2):91–8. doi.org/10.1016/S0378-4290(98)00104-X
Zhang JJ, Dell B, Conocono E, Waters I, Setter T, Appels R. Water deficits in wheat: fructan exohydrolase (1-FEH) mRNA expression and relationship to soluble carbohydrate concentrations in two varieties. New Phytol. 2009;181(4):843–50. doi:10.1111/j.1469-8137.2008.02713.x.
Hajheidari M, Eivazi A, Buchanan BB, Wong JH, Majidi I, Salekdeh GH. Proteomics uncovers a role for redox in drought tolerance in wheat. J Proteome Res. 2007;6(4):1451–60. doi:10.1021/pr060570j.
Nicolas ME, Lambers H, Simpson RJ, Dalling MJ. Effect of post-anthesis drought on cell division and starch accumulation in developing wheat grains. Ann Bot. 1985;55(3):433–44. doi:10.1093/oxfordjournals.aob.a086922.
Saini HS, Westgate ME. Reproductive development in grain crops during drought. Adv Agron. 1999;68(1):59–96. doi.org/10.1016/S0065-2113(08)60843-3
Nicolas ME, Lambers H, Simpson RJ, Dalling MJ. Effect of drought on metabolism and partitioning of carbon in two wheat varieties differing in drought-tolerance. Ann Bot. 1985b;55(5):727–47. doi:10.1093/oxfordjournals.aob.a086951
Ball S, Morell M. From bacterial glycogen to starch: understanding the biogenesis of the plant starch granule. Ann Rev Plant Biol. 2003;54(1):207–33. doi:10.1146/annurev.arplant.54.031902.134927.
Gao M, Wanat J, Stinard PS, James MG, Myers AM. Characterization of dull1, a maize gene coding for a novel starch synthase. Plant Cell. 1998;10(3):399–412. doi:10.1105/tpc.10.3.399.
Jane J, Chen YY, Lee LF, McPherson AE, Wong KS, Radosavljevic M, et al. Effects of amylopectin branch chain length and amylose content on the gelatizination and pasting properties of starch. Cereal Chem. 1999;76(5):629–37. doi:10.1094/CCHEM.19126.96.36.1999.
Roldan I, Wattebled F, Lucas MM, Delvallé D, Planchot V, Jiménez S, et al. The phenotype of soluble starch synthase IV defective mutants of Arabidopsis thaliana suggests a novel function of elongation enzymes in the control of starch granule formation. Plant J. 2007;49(3):492–504. doi:10.1111/j.1365-313X. 2006.02968.x.
Delvalle D, Dumez S, Wattebled F. Soluble starch synthase I: a major determinant for the synthesis of amylopectin in Arabidopsis thaliana leaves. Plant J. 2005;43(3):398–412. doi:10.1111/j.1365-313X.2005.02462.x.
Fujita N, Yoshida M, Asakura N, Ohdan T, Miyao A, Hirochika H, et al. Function and characterization of starch synthase I using mutants in rice. Plant Physiol. 2006;140(3):1070–84. doi:10.1104/pp.105.071845.
Zhang XL, Colleoni C, Ratushna V, Sirghle-Colleoni M, James MG, Myers AM. Molecular characterization demonstrates that the Zea mays gene sugary2 codes for the starch synthase isoform SSIIa. Plant Mol Biol. 2004;54(6):865–79. doi:10.1007/s11103-004-0312-1.
Riou P, Kjær S, Garg R, Purkiss A, George R, Cain RJ, et al. 14-3-3 proteins interact with a hybrid prenyl-phosphorylation motif to inhibit G proteins. Cell. 2013;153(3):640–53. doi:10.1016/j.cell.2013.03.044.
The English in this document has been checked by at least two professional editors of text-check, both native speakers of English. For a certificate, please see: http://www.textcheck.com/certificate/index/KKXgt2.
This research was financially supported by the grants from the Ministry of Science and Technology of China (2016YFD0100502, 2016ZX08009003–004) and the National Natural Science Foundation of China (31271703). These funding bodies did not play any roles in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Availability of data and materials
All relevant data are within this article and its additional files.
Ethics approval and consent to participate
The source of our grain seeds were from the Pro Zhimin Wang of China Agricultural University. The experimental research on plants, including collection of plant material, was complied with institutional, national, or international guidelines. And field studies were conducted in accordance with local legislation. And also comply with the Convention on the Trade in Endangered Species of Wild Fauna and Flora. This material has been deposited in a publicly available herbarium.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Efficiency and R2 values (coefficient of determination) of primer pairs. These values were determined using standard curves. Table S2. Granule-binding proteins identified by MALDI-TOF/TOF-MS. Table S3. Peptide information for the identified starch granule-associated proteins with at least twofold differences. Table S4. Total phosphorylated sites and phosphoproteins identified in this study. Table S5. Identification of proteins related to starch synthesis by matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF MS). (ZIP 124 kb)
Experimental workflow used in this study. DPA: days post anthesis. 1, 2, and 3 represent the three biological replicates. Figure S2. Scanning electron microscopy images for evaluation of the purity of the starch granules from the Chinese wheat cultivar Jingdong 17. Figure S3. SDS-PAGE of proteins in supernatant after successive washing steps, one, two and three times, respectively (1–3). (ZIP 1339 kb)
About this article
Cite this article
Chen, GX., Zhen, SM., Liu, YL. et al. In vivo phosphoproteome characterization reveals key starch granule-binding phosphoproteins involved in wheat water-deficit response. BMC Plant Biol 17, 168 (2017). https://doi.org/10.1186/s12870-017-1118-z
- Triticum aestivum L.
- Starch granule-binding proteins