Comprehensive Analysis of the Membrane Phosphoproteome Regulated by Oligogalacturonides in Arabidopsis thaliana

Abstract

Early changes in the Arabidopsis thaliana membrane phosphoproteome in response to oligogalacturonides (OGs), a class of plant damage-associated molecular patterns (DAMPs), were analyzed by two complementary proteomic approaches. Differentially phosphorylated sites were determined through phosphopeptide enrichment followed by LC-MS/MS using label-free quantification; differentially phosphorylated proteins were identified by 2D-DIGE combined with phospho-specific fluorescent staining (phospho-DIGE). This large-scale phosphoproteome analysis of early OG-signaling enabled us to determine 100 regulated phosphosites using LC-MS/MS and 46 differential spots corresponding to 34 pdhosphoproteins using phospho-DIGE. Functional classification showed that the OG-responsive phosphoproteins include kinases, phosphatases and receptor-like kinases, heat shock proteins (HSPs), reactive oxygen species (ROS) scavenging enzymes, proteins related to cellular trafficking, transport, defense and signaling as well as novel candidates for a role in immunity, for which elicitor-induced phosphorylation changes have not been shown before. A comparison with previously identified elicitor-regulated phosphosites shows only a very limited overlap, uncovering the immune-related regulation of 70 phosphorylation sites and revealing novel potential players in the regulation of elicitor-dependent immunity.

Introduction

Plants have developed various mechanisms to defend themselves against biotic stresses. Inducible immune defense responses include phytoalexin accumulation, expression of pathogenesis-related proteins, production of ROS, and, in some cases, programmed cell death (Boller and Felix, 2009). Plant innate immunity is driven by the perception of danger signals mediated by recognition proteins (Chisholm et al., 2006). Pathogen-associated molecular patterns (PAMPs) are conserved molecules secreted or present on the surface of most strains of a given microbial taxonomic group that activate the so-called PAMP-triggered immunity (PTI) against a wide range of pathogens (Barrett and Heil, 2012). Plant immunity also relies on the ability to sense danger by means of endogenous molecular patterns that are present only when the tissue is infected or damaged (damage-associated molecular patterns or DAMPs). In these cases, the discrimination between an intact self and an altered self leads to the activation of the immune system (Benedetti et al., 2015). Recognition of both PAMPs and DAMPs is mediated by the so-called pattern recognition receptors (PRRs; Boller and Felix, 2009).

Oligogalacturonides (OGs) are typical plant DAMPs (Ferrari et al., 2013) produced by the action of fungal endo-polygalacturonases (PGs) on homogalacturonan (HGA), the main component of pectin (Bellincampi et al., 2014). The interaction between PGs and the plant polygalacturonase-inhibiting proteins (PGIPs) favors the formation of OGs with degree of polymerization from 10 to 15 that activate the plant innate immunity (Mattei et al., 2005; Casasoli et al., 2009; Kalunke et al., 2015). OGs induce accumulation of phytoalexins, glucanase, and chitinase, ROS production [mediated, in Arabidopsis thaliana, by the RESPIRATORY BURST OXIDASE HOMOLOG D (RBOHD)] and callose deposition (Galletti et al., 2008; Ferrari et al., 2013). Exogenous treatment with OGs protects Arabidopsis and grapevine (Vitis vinifera) leaves against infection with the necrotrophic fungus Botrytis cinerea (Aziz et al., 2004; Ferrari et al., 2007), suggesting that, when PGs are secreted by microbes at the site of infection, this elicitor is likely to be released and contributes to activate defenses responses. Because pectin is one of the most accessible targets for many microbial cell wall–degrading enzymes and among the first structures to be altered during an attempted pathogenic attack, the signaling activity of OGs is an indication that plants have evolved mechanisms to monitor HGA degradation for the early detection of tissue injury (De Lorenzo et al., 2011; Nuhse, 2012; Savatin et al., 2014b). Wall-associated kinase (WAK) receptors are potential candidates to monitor pectin integrity (De Lorenzo et al., 2011; Kohorn and Kohorn, 2012). Indeed, WAK1 has been shown to mediate the perception of OGs (Brutus et al., 2010; Gramegna et al., 2016). OGs may also regulate plant growth and development mainly through their antagonism with auxin (Savatin et al., 2011; Ferrari et al., 2013), demonstrating the potential of this molecule to modulate both developmental and defense-dependent signaling.

Perception of elicitors at the plasma membrane (PM) triggers an intracellular signaling cascade that initiates pathogen defense. Early responses induced by OGs and flg22, a peptide PAMP derived from the bacterial flagellin, largely overlap (Denoux et al., 2008). For example, both activate calcium-dependent protein kinase (CDPK; Gravino et al., 2015) and mitogen-activated protein (MAP) kinase cascades (Rasmussen et al., 2012). The MAP triple kinases indicated as Arabidopsis NUCLEUS- AND PHRAGMOPLAST-LOCALIZED KINASE1-related protein kinases (ANPs), and the MAP single kinases MPK3 and MPK6 play a role in the response to OGs and PAMPs, including the bacterial elicitors flg22 and elf18 (Galletti et al., 2011; Savatin et al., 2014a). Nonetheless, distinctive features have been described between responses induced by PAMPs and OGs. Microarray analyses show considerable differences in the late responses to these two classes of elicitors (Denoux et al., 2008), and lack of ANPs strongly reduces phosphorylation of MPK3 and MPK6 induced by OGs but enhances that induced by elf18 (Galletti et al., 2011; Savatin et al., 2014a). Moreover, while the Arabidopsis leucine-rich repeat co-receptors BAK1/SERK3 and BKK1/SERK4 are required to achieve full response to elf18 and flg22, only a subset of defense responses induced by OGs is affected by loss of these elements, pointing to a complexity in the OG signaling pathways that is unique among the characterized MAMPs and DAMPs (Gravino et al., 2016). Specific phosphorylation events might regulate the activation of distinct signaling branches leading to different downstream responses.

In previous works, several OG-regulated proteins in the Arabidopsis apoplast and nucleus have been identified by 2-D DIGE (Casasoli et al., 2007, 2008), many of them present in multiple isoforms likely due to post-translational modifications (PTMs). Here we investigated early phosphorylation events regulated by OGs in a large-scale phosphoproteomic study including membrane proteins, to facilitate the detection of low abundance proteins with a signaling role. Analysis of membrane proteins by 2-D gel electrophoresis is limited by solubility constraints, that have a minor impact on shotgun proteomics (Kleffmann et al., 2007). On the other hand, greater proteome coverage can be reached by using gel-based and gel-free methods as complementary strategies (Zhao et al., 2008; Robbins et al., 2013). For a comprehensive picture of the phosphoproteome we therefore used both LC-MS/MS and the combination of 2-D DIGE with ProQ Diamond staining, which is known to selectively stain phosphoproteins (hereon indicated as Phospho-DIGE; Bond et al., 2011; Liu et al., 2015). In the latter, phosphorylated isoforms are particularly well resolved due to the property of the acidic phosphate group(s) to lower the pI of the proteins, thereby facilitating the detection of low-abundance phospho-isoforms.

Materials and methods

Growth conditions and plant treatments

A. thaliana ecotype Columbia-0 (Col-0) was used for this study. Plants were grown on soil (Einheitserde, Germany) in a climatic chamber at 22°C and 70% relative humidity. Seedlings were grown in a growth chamber at 21°C. Sterilized seeds (20 per well) were germinated and grown in 1 mL of liquid growth medium [Murashige and Skoog (MS) medium, pH 5.7, 0.5% sucrose] in 12-well plates.

For proteomic analyses, seedlings were grown in 500-ml flasks containing 100 mL of growth medium. Flasks containing seeds (about 100 seeds per flask) were grown in a growth chamber at 21°C. After 2 weeks, the culture medium was replaced with fresh medium and seedlings were grown for an additional day. About 10 g (fresh weight) of plant material were obtained from each flask. Plants and seedlings were grown under a 16-h light/8-h dark cycle (~120 μmol m⁻² s⁻¹).

OGs with degree of polymerization of 9–16 were prepared as previously described (Pontiggia et al., 2015). The OG stock solution (10 mg/mL) was filter-sterilized before addition to the medium to a final concentration of 50 μg/ml. For proteomic studies, water- (control), and OG-treated seedlings were harvested 10 min after treatment for protein extraction.

Preparation of total protein extract and total microsomal fraction

To obtain total protein extracts for Phospho-DIGE analysis, 10 g seedlings were homogenized using mortar and pestle in liquid nitrogen in homogenization buffer (1 M NaCl, 1 mM Na₃VO₄, 1 mM Na₂MoO₄, 25 mM NaF, protease inhibitor cocktail, Sigma). Total microsomal fractions (TMF) for both Phospho-DIGE and LC-MS/MS analyses were prepared as previous reported (Alexandersson et al., 2004; Tang et al., 2008). Briefly, 10 g seedlings were grinded using a blender with three impulses of 10 s in lysis buffer containing 1 mM Na₃VO₄, 1 mM Na₂MoO₄, 25 mM NaF, protease inhibitor cocktail (Sigma). The mixture was filtered through Miracloth and centrifuged at 10,000 × g for 10 min. Supernatant was then centrifuged at 100,000 × g for 1 h and the pellet (TMF) was recovered.

H⁺-ATPase immunoblotting

Protein samples (4 μg each) were boiled for 3 min in Laemmli buffer. Following gel electrophoresis, proteins were blotted onto nitrocellulose membranes using a Trans-Blot Turbo apparatus (BioRad). Subsequently blot was first incubated with the antibody against the plasma membrane H⁺-ATPase (1:500 dilution in PBS-T, kind gift of Prof M.I. De Michelis, University of Milan) and then with a horseradish peroxidase-coupled secondary antiserum (Santa Cruz, Biotechnology, Santa Cruz, CA; 1:2000 dilution in PBS-T).

MPK3/MPK6 immunoblotting

Immunodetection of phosphorylated MAPKs MPK3 and MPK6 was performed using a polyclonal phospho-p44/42-ERK MAPK-specific antiserum (Cell Signaling Technology, Danvers, MA). Protein samples (30 μg each) were boiled for 3 min in Laemmli buffer. Following gel electrophoresis, proteins were blotted onto nitrocellulose membranes using a Trans-Blot Turbo apparatus (Biorad). Subsequently blots were first incubated with the phospho-p44/42 MAPK-specific antiserum (1:2000 dilution in TBS-T at 4°C overnight) and then with a horseradish peroxidase-coupled goat anti-rabbit secondary antiserum (Santa Cruz, Biotechnology, Santa Cruz, CA; 1:5000 dilution in PBS-T). Membrane was stripped and re-probed with anti-MAPK3 and anti-MAPK6 (Sigma) and developed as reported above.

Protein digestion and phosphopeptide enrichment for LC-MS/MS analysis

TMF from control and OG-treated samples of two independent biological replicates were resuspended in 100 μL of freshly prepared 8 M urea, 10 mM Tris-HCl pH 8.0, 5% sodium deoxycolate. Protein concentration was determined using the BCA Protein Assay Kit (Sigma). For each sample, 150 ug of proteins were subjected to reduction with 10 mM dithioreitol for 45 min at 56°C followed by alkylation of cysteines with 50 mM iodoacetamide for 25 min at room temperature in darkness. Proteolytic digestion was carried out overnight with proteomics grade Lys-C (Promega, Lys-C:protein ratio 1:50) and for additional 4 h with trypsin (Promega, trypsin:protein ratio 1:50) at room temperature. The digestion mixture was subsequently acidified with 1% (v/v) formic acid and centrifuged to remove insoluble material. Peptides were desalted using home-made microcolumns using R3 beads (Poros) packed in gel loader tips. Phosphopeptide enrichment was performed using a sequential elution from immobilized metal affinity chromatography (SIMAC) as previously described (Thingholm et al., 2008). The procedure includes an initial phosphopeptide enrichment step by IMAC, in which bound peptides are sequentially eluted with an acidic and a basic solution. Next, the IMAC flow-through fraction and the acidic elution fractions are subjected to a second phosphopeptide enrichment step using TiO₂ chromatography. The TiO₂-enriched fraction and the basic elution fractions from IMAC were lyophilized and analyzed by LC-MS/MS. For proteomic analysis, also flow-through fractions from TiO₂ chromatography were desalted on home-made R3 stage tips and analyzed by LC-MS/MS.

Proteomic and phosphoproteomic analysis by LC-MS/MS

Peptides were separated on a 15 cm PicoFrit Column (75 μm i.d., New Objective, Woburn, MA) packed with Magic C18AQ (5 μm, 200 Å, Michrom) and analyzed on a LTQ-Orbitrap Discovery (Thermo Fisher Scientific) coupled on-line with a nano-HPLC (Ultimate 3000, Thermo Fisher Scientific). Peptides were eluted using a 0–60% acetonitrile gradient in 0.1% formic acid (180 min, 300 nl/min). MS was acquired at 30,000 FWHM resolution in the FTMS (target value of 5 × 10⁵ ions) and MS/MS was carried out in the linear ion trap, with MS/MS on top 5 ions and multistage activation.

LC-MS/MS data analysis

MS data were processed using MaxQuant (Cox et al., 2014) version 1.5.2.8 with a false-discovery rate (FDR) < 0.01 at the level of proteins, peptides and modifications, using default settings. Oxidized methionine, acetylation (protein N-term) and phospho(STY) were selected as variable modifications, and carbamidomethyl cysteine as fixed modification. Proteins were identified using a target-decoy approach with a reversed database, using the Andromeda search engine against the Arabidopsis UniProt database (release-2014, 31565 entries) and LFQ quantification for phosphoproteomics and proteomics was performed by MaxQuant.

Statistical analysis was performed with Perseus (version 1.5.0.31). Hits to the reverse database and contaminants were filtered out, intensities were normalized for unequal protein amounts and log₂-transformed. For proteomic analysis quantification was performed on proteins identified with a minimum of 2 unique peptides and at least three valid values among the replicates. Significance was assessed by Student's t-test, using permutation-based FDR to control for multiple hypothesis testing (Data Sheet S2B in Supplementary Material). For the phosphoproteomic study, we first filtered to retain only class I phosphosites (localization probability > 0.75 and score difference > 5). We next considered (1) phosphosites that were exclusive for either the OG-treated samples or the control samples with two valid values in either the control or the OG-treated group; (2) phosphosites that had at least three valid values (control + OG-treated). These sites are listed in Data Sheet S1B in Supplementary Material. For phosphosites that had three valid values only, first, missing values were imputed with values representing a normal distribution around the detection limit of the mass spectrometer (downshift = 1.5; width = 0.3), to allow statistical analysis; then Significance B-values for each replicate (OG-treated vs. control) were calculated using the statistical tool in Perseus (Cox and Mann, 2008). Significance B-values represent outlier probability score weighted for the intensity. Besides the exclusive ones, phosphosites with a significance B p ≤ 0.05 in both replicates were considered as significantly regulated. In addition, a two-sample student's t-test was performed, and phosphosites that had a significance B ≤ 0.05 in one replicate and passed the t-test with p < 0.05 were also considered as significantly regulated.

2D DIGE, Pro-Q diamond staining, and image acquisition

The Clean-up kit (GE Healthcare) was used to eliminate salts and concentrate proteins in both the total protein extracts and the TMF. The protein-containing pellets were finally dissolved in a small volume (50–100 μL) of IEF buffer containing 8 M urea, 2 M thiourea, and 4% (v/v) CHAPS, centrifuged and the supernatants recovered. Protein content was measured using the BIO-RAD protein assay according to the manufacturer's instructions. Total protein extracts and TMF were labeled and analyzed using the same procedure: control and OG-treated protein samples from 3 independent biological replicates (C1, C2, C3; T1, T2, T3) were labeled with Cy3 or Cy5, according to the manufacturer's instructions for minimal labeling (Minimal labeling kit, GE Healthcare) and randomization (ETTAN DIGE System, GE Healthcare). Briefly, each sample (50 μg) was labeled with 200 pmol CyDye DIGE Fluor minimal dyes (GE Healthcare) and incubated on ice in the dark for 30 min. The internal standard (IS) was labeled with Cy2 and consisted of a pooled sample comprising equal amounts of each sample used for each replicate (C1+C2+C3+T1+T2+T3). Equal volumes of a 2X buffer [8 M urea, 2 M thiourea, 2% ASB-14 (or 4% CHAPS; v/v), 20 mg/mL DTT, and 2% IPG buffer/Pharmalytes 4-7] were added to each labeled protein samples. Samples (C+T+IS) were pooled prior to IEF, which was carried out using non-linear IPG strips (pH 4–7, 13 cm, GE Healthcare) rehydrated overnight at room temperature. IEF was performed using the Ettan IPG-phor apparatus (GE Healthcare) as follows: (1) step to 500 V (0.5 kVh); (2) step to 1000 V (0.8 kVh); (3) step to 8000 V (11.3 kVh); (4) step to 8000 V (3.0 kVh), for a total of 25 kVh at 20°C and a maximum current setting of 50 μA per strip. After IEF, the strips were equilibrated for 15 min in 100 mM Tris pH 6.8, 30% v/v glycerol, 8 M urea, 1% w/v SDS, 0.2 mg/mL bromophenol blue, 5 mg/mL DTT for reduction of disulfide bridges, and alkylated for 15 min in the same equilibration buffer containing 25 mg/mL iodoacetamide. Each strip was then loaded on top of a 12% w/v acrylamide gel for the second dimension. SDS-PAGE was carried out using the Hoefer SE 600 Ruby apparatus (GE Healthcare) at 100 V and 20 mA per gel for 15 min and then 100 V and 40 mA per gel until the bromophenol blue dye front had run off the bottom of the gel. A running buffer of 25 mM Tris pH 8.3, 192 mM glycine, and 0.1% w/v SDS was used.

Each gel was scanned using a Typhoon 9200 imager (GE Healthcare) set at the wavelengths corresponding to each CyDye [532 nm laser and 580 nm band pass (BP) 30 emission filter for Cy3; 633 nm laser and 670 nm BP30 emission filter for Cy5; 488 nm laser and 520 nm BP40 emission filter for Cy2], at high resolution (100 μm). The photo multiplier tube (PMT) voltage was adjusted for each channel between 500 and 700 V to ensure that the image intensity was within a linear range between 40,000 and 60,000 U. For phosphoprotein detection, each gel was then fixed, post-stained with Pro-Q Diamond phosphospecific fluorescent dye (Invitrogen; Bond et al., 2011; Liu et al., 2015) applying the same scanning conditions described above. Pro-Q Diamond and Cy3 have very similar spectra of fluorescence with excitation/emission maxima at 555/580 and 553/569 nm for Pro-Q Diamond and Cy3, respectively. Therefore, phosphoproteins were detected as spots with increased Cy3 fluorescence (Cy3 + Q) in comparison to original Cy3.

Phospho-DIGE statistical analysis

For data normalization and analysis, gel images acquired before and after Pro-Q Diamond post-staining were compared using the Differential In-gel Analysis (DIA) module of the DeCyder software version 6.5 (GE Healthcare) for co-detection of the three CyDye-labeled forms of each spot and calculation of the ratios between sample and internal standard abundance. Phosphoproteins were detected as spots with increased Cy3-like fluorescence due to ProQ-staining in comparison with original Cy3 emission (Stasyk et al., 2005). Statistical analysis of protein abundance changes between control and treated samples from the three independent biological replicates was performed in the Biological Variation Analysis (BVA) module for quantitative comparisons of protein abundance/phosphorylation across multiple gels. The inter-gel variability was corrected by normalization of the Cy2 internal standard spot map present in each gel. Protein spots whose intensities were significantly different among the control and treated samples were determined by paired one-way ANOVA. Proteins were considered differentially expressed when the FDR-corrected p-values of the ANOVA analysis were < 0.05.

Protein spot identification by MALDI-TOF mass spectrometry

For mass spectrometry analysis of proteins identified as differentially phosphorylated by Phospho-DIGE experiments, a preparative gel was run under the same IEF and SDS-PAGE conditions used for the DIGE gels, loading 250 μg each of control and OG treated samples mixed together. Proteins were subsequently visualized using Coomassie Brilliant Blue (CBB) R-250 stain (Sigma) according to the manufacturer's instructions. The preparative gel image was matched to the master gel image (the gel with the highest spot count in the DIGE analytical gel match-set) using DeCyder software. Matching was further improved by land marking and manually confirming potential spots of interest. By comparing the CBB-stained spot pattern with the corresponding Cy5 protein pattern, spots of interest showing differential fluorescent levels on the 2-D DIGE gels were picked manually from the preparative gel and subjected to trypsin in-gel digestion as previously described (Casasoli et al., 2007).

Protein identification was carried out using a Voyager-DE STR instrument (Applied Biosystems, Framingham, MA). Peptides were desalted using ZipTip C18 microcolumns (Millipore, Bedford, MA, USA) and spotted onto a MALDI target plate using CHCA as the matrix (10 mg/ml α-cyano-4-hydroxy-cinnamic acid in 0.1% TFA, 50% ACN). The mass spectrometer was operated in the positive ion, delayed extraction (200 ns delay time) reflector mode with an accelerating voltage of 20 kV. Each MALDI-TOF spectrum was generated by accumulating data corresponding to 200–500 laser shots. Internal mass calibration was performed by using theoretical masses of the trypsin autodigestion peaks. The mass list was then analyzed using the PeakErazor software (http://www.welcome.to/GPMAW) to eliminate contaminant peaks (keratin, trypsin added for digestion and peaks present in all mass spectra). Proteins were identified by searching the National Center for Biotechnology Information (NCBI) database using the MASCOT (www.matrixscience.com) or Aldente (xpasy.org/tools/aldente) search engines with the following criteria: cleavage by trypsin (one missed cleavage allowed), mass tolerance 30 ppm; carbamidomethyl cysteine as fixed modification, methionine oxidation, and phospho(STY) as variable modifications. Only matches with a Mascot score higher than 60 (p < 0.05; or Aldente score higher than 13.86, p < 0.05), sequence coverage higher than 10% and more than 6 peptide matches were considered significant. A number of photosynthetic proteins were identified as potentially regulated phosphoproteins. Because they corresponded to highly abundant protein spots amenable to mis-quantification, they were excluded from our dataset.

Bioinformatic analysis of phosphoproteins

The PhosPhAt (http://phosphat.mpimp-golm.mpg.de/; Durek et al., 2010) and P3DB (http://www.p3db.org/index.php; Yao et al., 2014) databases were searched to determine if phosphorylation sites had been previously reported for the identified proteins. For the analysis of significantly over-represented GO terms, differentially phosphorylated proteins were analyzed using AgriGO (http://bioinfo.cau.edu.cn/agriGO/analysis.php) for Singular Enrichment Analysis (SEA) using A. thaliana TAIR10 as selected species and Arabidopsis genome locus (TAIR10) as selected reference, Fisher as a statistical test method, Yekutieli (FDR under dependency) as z multi-test adjustment method, 0.05 as significance level, 5 as minimum number of mapping entries and complete GO as Gene Ontology type. In order to avoid redundancy at the protein level, in the case of protein groups we considered only the first isoform.

To investigate possible interactions between the OG-regulated phosphoproteins, the STRING (Search Tool for the Analysis of Interacting Genes/Proteins) algorithm was used for the creation of protein interaction networks based on published functional or informatics-predicted interactions (Szklarczyk et al., 2015). The evidence annotation in STRING was filtered out of interactions from text-mining and neighborhood, and only interactions supported by experimental evidence, co-expression and existing database information, with high-confidence score >0.7 were considered (Table S6). The SUBA3 database (http://suba3.plantenergy.uwa.edu.au/; Heazlewood et al., 2005) was used to assign the subcellular localization of differentially phosphorylated proteins.

Results

OG treatment induces early changes in the phosphoproteome of arabidopsis seedlings

In Arabidopsis, treatment with OGs induces the activation, by phosphorylation, of MPK3 and MPK6 within few minutes, a temporal kinetics similar to that observed upon treatment with PAMPs (Galletti et al., 2011). Two-week-old liquid-grown Arabidopsis ecotype Columbia 0 (Col-0) seedlings were treated with OGs (50 μg/ml) or water as a control. The 10-min time point was chosen for our analyses, because close to that used in previous studies on elicitor-induced phosphorylation (Benschop et al., 2007; Nuhse et al., 2007; Rayapuram et al., 2014). Western blot analysis using α-phospho-p44/42-ERK antibody showed the OG-induced phosphorylation of MPK3 and MPK6 (Figure 1A). TMF was obtained to facilitate the identification of less abundant membrane proteins. Enrichment of total membranes was determined by western blot using antibodies specific for the PM H⁺-ATPase (Figure 1B). In addition, total protein extracts were obtained, for Phospho-DIGE analysis only.

Figure 1

Phosphoproteomic analysis of TMF by LC-MS/MS

TMF preparations from seedlings treated with OGs or water as a control were subjected to the SIMAC phosphopeptide enrichment method. The enriched fractions were analyzed by nano-LC-MS/MS, leading to the identification of a total of 2147 unique phosphosites [Data Sheet S1C in Supplementary Material]. We could quantify 1026 phosphosites (Data Sheet S1B in Supplementary Material). A total of 99 phosphosites significantly changed phosphorylation level upon OG treatment; among them, 29 and 36 phosphosites were unique to control and OG-treated seedlings, respectively (Table 1).

Table 1

Protein nameb	TAIR IDc	Locd	Uniprot IDe	Posf	Phosphopeptide (best evidence)g	Average log2 ratio Tr/Ch	Log2 ratio Tr1/C1i	Log2 ratio Tr2/C2i	μarray OG 1 h/3 hj
TRANSPORT
Aquaporin PIP2-1; Aquaporin PIP2-2; Aquaporin PIP2-3	At3g53420; At2g37170; At2g37180	PM§	P43286; P43287; P30302	280; 278 278	SI GS(0.996)FRS(0.004)AANV	Tr	Tr	Tr	ns
ABC transporter B family member 20 (ABC-B20)	At3g55320	PM	Q9M3B9	712	GS(0.001)GVFRPQEICFDT(0.978)EES(0.054)PKAHS(0.216)PAS(0.751)EK	Tr	Tr	Tr	ns
Translocase of chloroplast 159 (TOC159)	At4g02510	PL§	O81283	686	IDGQIVT(0.003)DS(0.996)DEDVDT(0.001)EDEGEEK	Tr	Tr	Tr	ns
ABC transporter C family member 14 (ABCC14)	At3g62700	VO	Q9LZJ5	903	S(0.348)IS(0.077)IES(0.575)PRQPKS(1)PK	Tr	Tr	Tr	−1.5/ns
ABC transporter C family member 14 (ABCC14)	At3g62700	VO	Q9LZJ5	925	T(0.333)T(0.333)S(0.333)MES(0.999)PRISEINDESIK	Tr	Tr	Tr	−1.5/ns
ABC transporter G family member 36 (PEN3)	At1g59870	PM§	Q9XIE2	45	T(0.011)QS(0.989)VNDDEEAIK	5.70	4.77	6.62	2.8/ns
ABC transporter C family member 14 (ABCC14)	At3g62700	VO	Q9LZJ5	897	S(0.009)IS(0.991)IES(1)PRQPKS(1)PK	4.63	3.24	6.02	−1.5/ns
ABC transporter C family member 14 (ABCC14)	At3g62700	VO	Q9LZJ5	894	S(0.009)IS(0.991)IES(1)PRQPKS(1)PK	4.55	3.24	5.87	−1.5/ns
ATPase 1 (AHA1), plasma membrane	At2g18960	PM§	P20649	544	DS(1)NIASIPVEEIIEK	4.17	3.83	4.50	ns
Transport protein particle (TRAPP) component	At3g05000	CY	Q9CAW4	121	VSIDPSSENETQDPST(0.999)PGES(0.001)K	3.98	4.17	3.79	ns
ATPase 2 (AHA2), plasma membrane; ATPase 1 (AHA1), plasma membrane	At4g30190; At2g18960	PM§	P19456; P20649	899; 899	GS(1)YREISEIAEQAK	3.42	3.08	3.77	ns/2.3
Copper transporter 5 (CuTR5)	At5g20650	GO§	Q93VM8	68	S(0.018)S(0.018)S(0.962)GVS(0.001)APIIPK	C	C	C	ns
UDP-galactose transporter 1	At1g77610	GO	Q9C521	310	HMIS(0.001)QQT(0.006)PGT(0.821)PRT(0.821)PRT(0.35)PR	C	C	C	−1.6/ns
UDP-galactose transporter 1	At1g77610	GO	Q9C521	313	HMIS(0.001)QQT(0.006)PGT(0.821)PRT(0.821)PRT(0.35)PR	C	C	C	−1.6/ns
NSP (Nuclear shuttle protein)-interacting GTPase	At4g13350	CY	Q8W4K6	407	S(0.214)MS(0.779)APS(0.007)IQPIQGVPSGGIQSSEVKPSGR	C	C	C	na
NSP (Nuclear shuttle protein)-interacting GTPase	At4g13350	CY	Q8W4K6	219	T(0.045)S(0.045)EEGS(0.185)QS(0.725)PEQVKDIGS(0.199)AS(0.801)PPVARPVR	C	C	C	na
SIGNAL TRANSDUCTION
Transcription regulator NOT2/NOT3/NOT5 family protein	At5g18230	NU	F4JWJ6	431	NIMGVESNVQPIT(0.043)S(0.891)PIS(0.066)K	Tr	Tr	Tr	ns
Transducin family protein (WD-40 RFP)	At5g05570	MI§	F4K0R0	464	ART(0.928)PRT(0.661)PS(0.303)GES(0.107)AQWPIT(0.002)GGVPSHVDDYK	Tr	Tr	Tr	ns
Phospholipase like protein (PEARLI 4f) family	At4g38550	NU	Q9C5F6	232	EQFEDIYEQDGDVT(1)PR	Tr	Tr	Tr	6.7/ns
Phototropin-2 (PHOT2)	At5g58140	PM§	F4KDJ3	504	GEIQYFIGVQIDGSDHVEPIQNRIS(1)ER	Tr	Tr	Tr	ns
Mitogen-activated protein kinase 6 (MPK6)	At2g43790	CY§	Q39026	221	VTSESDFMT(1)EY(1)VVTR	Tr	Tr	Tr	ns
Mitogen-activated protein kinase 6 (MPK6)	At2g43790	CY§	Q39026	223	VTSESDFMT(1)EY(1)VVTR	Tr	Tr	Tr	ns
Calmodulin-binding protein(EDA39)	At4g33050	CY	O82645; F4JVX1	462; 500	VS(0.001)PANS(0.003)Y(0.003)GPIPS(0.996)PRPS(0.998)PK	Tr	Tr	Tr	24.1/2.9
Calmodulin-binding protein(EDA39)	At4g33050	CY	O82645; F4JVX1	466; 504	VS(0.001)PANS(0.003)Y(0.003)GPIPS(0.996)PRPS(0.998)PK	Tr	Tr	Tr	24.1/2.9
Phospholipase D gamma 1(PLDγ1)	At4g11850	VO	Q9T053	680	EVPVGTVS(0.001)VY(0.014)NS(0.984)PR	3.83	4.47	3.20	ns
Protein kinase family protein	At5g03320	PM	Q9LZF8	17	DEQRS(1)PKPVS(0.806)PT(0.082)S(0.027)NFS(0.085)DVNK	2.79	3.15	2.43	2.5/ns
Putative tyrosine phosphatase (PTEN2a)	At3g19420	NU	Q9LT75	509	ETENPDKDDVFS(1)DNEGDSTGPTK	2.62	2.42	2.81	1.5/ns
Cyclin-dependent protein kinase-like protein	At5g44290	NU	Q9FKV9	75	HQEIAEIGDT(1)DEDEDDDHHPPEEIK	−2.67	−2.31	−3.03	2.9/ns
Auxin response factor ARF1	At1g59750	NU	F4ID31	399	PRPPGIPS(0.084)PAT(0.06)GPS(0.163)DGVWKS(0.159)PADT(0.664)PS(0.825)S(0.825)VPIFS(0.221)PPAK	−2.74	−2.98	−2.50	ns
Auxin response factor ARF1	At1g59750	NU	F4ID31	400	PRPPGIPS(0.084)PAT(0.06)GPS(0.163)DGVWKS(0.159)PADT(0.664)PS(0.825)S(0.825)VPIFS(0.221)PPAK	−2.74	−2.98	−2.50	ns
Serine/threonine proteinkinase	At1g48210	PM	F4HWU0; Q9LUT0	353; 354	AIQPIINPPRS(0.986)APQT(0.986)PHRNPY(0.027)	C	C	C	ns
Serine/threonine proteinkinase	At1g48210	PM	F4HWU0; Q9LUT0	357; 358	AIQPIINPPRS(0.986)APQT(0.986)PHRNPY(0.027)	C	C	C	ns
MAP3K epsilon protein kinase (MKKK7)	At3g13530	PM§	Q9LJD8	788	VRS(1)GQIDPNNPIFGQNETSSISMIDQPDVIK	C	C	C	ns
Probable LRR receptor-like serine/threonine-protein kinase	At5g37450	PM	C0LGU1	268	S(0.886)IVIY(0.059)Y(0.059)IDIS(0.578)S(0.889)NKIT(0.528)GEIPKNK	C	C	C	ns
Probable LRR receptor-like serine/threonine-protein kinase	At5g37450	PM	C0LGU1	278	S(0.886)IVIY(0.059)Y(0.059)IDIS(0.578)S(0.889)NKIT(0.528)GEIPKNK	C	C	C	ns
GSK3/Shaggy-related protein kinase alpha (ASK1)	At5g26751	PM	P43288	229	GEPNIS(0.232)Y(0.768)ICSR	C	C	C	1.6/ns
Probable serine/threonine-protein kinase	At1g18390	PM	P0C5E2	643	S(0.003)GPIVAQS(0.761)PDS(0.235)VIVK	C	C	C	6.8/1.5
Phospholipase-like protein (PEARLI 4)	At2g20960	NU	Q9SKR5	314	S(0.829)KT(0.169)PEPQPT(0.001)Y(0.001)FEPSSR	C	C	C	4.0/ns
TRAFFICKING
Patellin-2 (PATL2)	At1g22530	PM	Q56ZI2	79	EIIQS(0.064)ES(0.936)FK	Tr	Tr	Tr	−2.0/−1.6
Patellin-4 (PATL4)	At1g30690	PM	Q94C59	53	S(0.126)AS(0.874)FKEESDFFADIK	Tr	Tr	Tr	ns
ADP-ribosylation factor GTPase-activating protein (ZIGA4)	At1g08680	NU	F4HXP0	236	S(0.001)DIQS(0.999)PNFQQEAEFR	Tr	Tr	Tr	ns
Kinesin-like protein KAC1	At5g10470	CY§	Q9LX99	611	TGDAIQS(0.001)QDIFS(0.999)PIPK	Tr	Tr	Tr	ns
Kinesin-like protein KAC1	At5g10470	CY§	Q9LX99	698	GEGYSAEAVAIPS(0.063)T(0.937)PNK	Tr	Tr	Tr	ns
65-kDa microtubule-associated protein 1 (MAP65-1)	At5g55230	NU; CS§	F4K3E4; Q9FLP0	615; 586	EEAASSPVSGAADHQVPAS(1)P	Tr	Tr	Tr	ns
Myosin-17 (XIK)	At5g20490	ER§	F4K5J2	990	QQAIAIS(0.957)PT(0.041)S(0.002)R	3.86	5.22	2.49	ns
Patellin-3 (PATL3)	At1g72160	CY	Q56Z59	108	SMIPQNIGS(1)FK	3.30	2.95	3.66	ns
Membrane fusion protein Use1	At3g55600	NU	Q6NKR3	104	IEDEPRS(0.999)PT(0.973)S(0.028)PQIR	2.42	2.91	1.94	1.8/ns
Syntaxin of plants 121 (SYP121)	At3g11820	PM§; GO§	Q9ZSD4	264	AS(0.213)S(0.755)FIRGGT(0.031)DQIQTAR	−3.06	−2.37	−3.74	7.5/ns
Syntaxin of plants 122 (SYP122)	At3g52400	PM§	Q9SVC2	18	MIAIFHEAFAHPPEEINS(0.983)PAS(0.017)EK	−4.08	−2.61	−5.54	11.5/ns
RESPONSE TO STRESS
Respiratory burst oxidase homolog protein D	At5g47910	PM§	Q9FIJ0	343	IIS(1)QMIS(1)QK	Tr	Tr	Tr	9.9/ns
Respiratory burst oxidase homolog protein D	At5g47910	PM§	Q9FIJ0	347	IIS(1)QMIS(1)QK	Tr	Tr	Tr	9.9/ns
Metal tolerance protein C2	At3g12100	VO	F4J8M5	37	S(0.99)FQQS(0.212)HGPRT(0.688)PAVT(0.11)EAAKPFIDR	Tr	Tr	Tr	ns
Aluminium induced protein with YGL and LRDR motifs	At5g43830	CY	Q9FG81	18	TVANSPEAIQS(0.997)PHS(0.003)SESAFAIK	5.95	6.83	5.06	ns
Aluminium induced protein with YGL and LRDR motifs	At4g27450	CY	Q93V62	219	VDS(1)EGVICGANFK	−2.83	−3.56	−2.10	ns
Protein plant cadmium resistance 8 (CdRes8)	At1g52200	PM	Q9M815	12	GRVT(0.002)T(0.002)PS(0.184)EEDS(0.812)NNGIPVQQPGT(1)PNQR	C	C	C	6.1/2.3
Putative defensin-likeprotein 135	At4g09647	EX	Q2V3K1	2	S(0.946)KT(0.849)FQPS(0.176)IVMIT(0.036)IFIIIVT(0.996)S(0.997)QR	C	C	C	ns
Putative defensin-likeprotein 135	At4g09647	EX	Q2V3K1	21	S(0.946)KT(0.849)FQPS(0.176)IVMIT(0.036)IFIIIVT(0.996)S(0.997)QR	C	C	C	ns
Putative defensin-likeprotein 135	At4g09647	EX	Q2V3K1	4	S(0.946)KT(0.849)FQPS(0.176)IVMIT(0.036)IFIIIVT(0.996)S(0.997)QR	C	C	C	ns
Putative defensin-likeprotein 135	At4g09647	EX	Q2V3K1	20	S(0.946)KT(0.849)FQPS(0.176)IVMIT(0.036)IFIIIVT(0.996)S(0.997)QR	C	C	C	ns
Heat shock protein 90-1 (HSP90-1)	At5g52640	CY; NU	P27323	219	EIS(1)DDEDEDEPK	C	C	C	5.0/4.6
PROTEIN METABOLIC PROCESS
Pumilio homolog 3; P1umilio homolog 1; Pumilio homolog 2	At2g29140; At2g29200; At2g29190	CY; NU	Q9ZW02; Q9ZW07; Q9ZW06	268; 272; 280	VPS(1)PCIT(0.999)PIGS(0.001)GR	Tr	Tr	Tr	ns
Pumilio homolog 3; Pumilio homolog 1; Pumilio homolog 2	At2g29140; At2g29200; At2g29190	CY; NU	Q9ZW02; Q9ZW07; Q9ZW06	272; 276 284	VPS(1)PCIT(0.999)PIGS(0.001)GR	Tr	Tr	Tr	ns
40S ribosomal protein S3-2	At3g53870	CY	Q9M339	212	TPIPDVVIIHS(1)PK	Tr	Tr	Tr	ns
RING/U-box domain-containing protein	At3g06330	GO	B3H497	144	IS(0.989)PVVS(0.004)T(0.004)QIS(0.002)AGVPGDPPNK	Tr	Tr	Tr	ns
Enhancer of mRNA-decapping protein 4	At3g13300	CY	Q9LTT8	98	TISYPTPPINIQS(1)PR	4.28	4.09	4.47	ns
Eukaryotic translation initiation factor 4G (EIF4G)	At3g60240	CY	A8MR97; Q76E23	1525 1529	QVIQGPSAT(0.005)VNS(0.995)PR	3.45	3.37	3.53	ns
40S ribosomal proteinS2-3	At2g41840	CY	P49688	273	AIS(0.845)T(0.149)S(0.006)KPDPVVEDQA	2.71	2.51	2.91	ns
Serine/arginine-rich splicing factor 30 (SCL30)	At3g55460	NU§	Q8L3X8	204	RS(1)Y(0.004)S(0.996)PGYEGAAAAAPDRDR	1.84	1.26	2.43	ns
Serine/arginine-rich splicing factor 30 (SCL30)	At3g55460	NU§	Q8L3X8	206	RS(1)Y(0.004)S(0.996)PGYEGAAAAAPDRDR	1.84	1.26	2.43	ns
Putative translation initiation factor eIF-2B epsilon subunit (eIF-2B)	At2g34970	CY	O64760	447	VSIIQQPT(0.049)T(0.049)EDS(0.899)DEEIEY(0.003)ADSSSGTADHISGINIQMESK	C	C	C	ns
METABOLISM
NADPH-cytochrome P450 reductase (CPR)	At4g24520	ER	F4JQY4	248	S(0.002)VAT(0.993)PYT(0.005)AVIPEYR	Tr	Tr	Tr	ns
Glutamate decarboxylase 1	At5g17330	CY; PM; GO	Q42521	8	VIS(0.064)HAVS(0.751)ES(0.232)DVS(0.935)VHS(0.008)T(0.007)FAS(0.003)R	Tr	Tr	Tr	ns
6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase (FKFBP)	At1g07110	ER§; GO§	Q9MB58	276	SVETIS(1)PFQQK	2.39	2.23	2.55	ns
Photosystem Q(B) protein (psbA)	AtCg00020	PL	P83755	2	T(1)AIIERR	−2.57	−2.42	−2.71	na
Probable pectin methyltransferase PMT13 (QUA3)	At4g00740	GO	Q93W95	280	CIIPFT(0.775)AY(0.224)NAT(0.252)Y(0.749)FIEVDRIIR	−2.93	−2.04	−3.82	−2.0/ns
Probable pectin methyltransferase PMT13 (QUA3)	At4g00740	GO	Q93W95	286	CIIPFT(0.53)AY(0.468)NAT(0.037)Y(0.966)FIEVDRIIR	−2.93	−2.04	−3.82	−2.0/ns
Cytochrome P450, family 77,subfamily A, polypeptide 6(CYP77A6)	At3g10570	ER	Q9SQY7	284	AIQKPGT(0.927)DKT(0.164)AS(0.154)S(0.774)FS(0.754)Y(0.262)IDT(0.964)IFDIK	−4.38	−4.92	−3.84	ns
Cytochrome P450, family 77,subfamily A, polypeptide 6(CYP77A6)	At3g10570	ER	Q9SQY7	282	AIQKPGT(0.927)DKT(0.164)AS(0.154)S(0.774)FS(0.754)Y(0.262)IDT(0.964)IFDIK	−4.50	−4.92	−4.08	ns
ADP-glucose pyrophosphorylase family protein	At1g74910	CY	F4HXD1	217	RVS(0.763)S(0.237)FEAIQPATR	C	C	C	ns
OTHER PROTEINS
Tetratricopeptide repeat domain protein (TSS)	At4g28080	PX	F4JKH6	1311	TSSNEISISVAGS(0.003)T(0.044)S(0.181)S(0.759)PAS(0.014)K	Tr	Tr	Tr	ns/−1.7
DEAD-box ATP-dependent RNA helicase 38 (LOS-4)	At3g53110	CY; NU	Q93ZG7	25	ADTVEKVPTVVES(0.005)S(0.004)S(0.007)S(0.013)S(0.011)T(0.011)VEAS(0.186)NS(0.762)AEK	Tr	Tr	Tr	ns
Putative bZIP protein	At3g60320	NU	Q93YU8	147	IPHIIS(0.15)ES(0.649)S(0.189)PS(0.019)S(0.08)S(0.912)PR	Tr	Tr	Tr	−3.0/ns
Nuclear/nucleolar GTPase(AtNug2)	At1g52980	NU	Q9C923	116	ERKIPMS(1)IIT(1)DNK	Tr	Tr	Tr	ns
Nuclear/nucleolar GTPase(AtNug2)	At1g52980	NU	Q9C923	119	ERKIPMS(1)IIT(1)DNK	Tr	Tr	Tr	ns
Uncharacterized protein	At1g77765	CY	Q6DYD4	94	CDVVES(0.014)DNKPERIS(0.926)PS(0.059)PK	Tr	Tr	Tr	na
Uncharacterized protein	At1g19530	NU	Q93WK6	23	S(0.865)IT(0.135)REEIDTFWK	Tr	Tr	Tr	ns
Uncharacterized membrane protein	At3g27390	PM	Q8GUM4	544	DES(0.953)IT(0.047)EPPAPVK	Tr	Tr	Tr	ns
Uncharacterized protein	At5g13260	NU	Q8VZL1	446	ISDIEIKS(1)PGGPK	3.65	3.11	4.20	ns
SAP domain-containing protein	At4g39680	NU	O65655	417	VPEAQIT(0.002)NS(0.003)AT(0.995)PT(0.017)T(0.053)T(0.931)PR	2.91	3.20	2.62	ns
SAP domain-containing protein	At4g39680	NU	O65655	421	VPEAQIT(0.44)NS(0.28)AT(0.28)PTT(0.001)T(0.999)PR	2.57	2.53	2.62	ns
Uncharacterized protein	At5g52980	CY	Q9LVV4	10	DQPDPKNGS(0.99)GIIIS(0.922)AT(0.088)EPIR	−1.99	−1.63	−2.35	ns
Uncharacterized protein	At5g52980	CY	Q9LVV4	15	DQPDPKNGS(0.99)GIIIS(0.922)AT(0.088)EPIR	−1.99	−1.63	−2.35	ns
Midasin (MDN1)	At1g67120	NU	F4HRR8	4797	CGS(1)PQKEEPGNDIEQEPETEPIEGK	C	C	C	ns
DEK domain-containing chromatin associated protein (DEK-C)	At4g26630	CY	Q9SUA1	523	SIAHS(0.004)DDES(0.996)EEEKEEEEK	C	C	C	ns
DNA topoisomerase 1 beta	At5g55310	NU	Q9FJ79	126	APS(0.002)VS(0.073)KS(0.9)DDEDS(0.024)EDDKPISAR	C	C	C	ns
BTB/POZ domain-containing protein	At5g67385	PM	Q66GP0	558	TSSSTIS(0.004)T(0.036)NPS(0.474)S(0.474)PIS(0.025)T(0.025)AS(0.068)T(0.895)GKPPIPR	C	C	C	na
Uncharacterized protein	At1g20970	PM	F4HWC3	203	NNVEEPEVEIES(0.065)DS(0.162)ET(0.773)DVEGHQGDK	C	C	C	ns
Uncharacterized protein	At2g15860	NU	Q8L4R9 ; F4IJD7	178; 190	IENS(0.005)VQQGS(0.101)S(0.894)PR	C	C	C	ns
Uncharacterized protein	At4g31880	NU	F4JTF2; Q8GUP3	558; 559	AIDEES(0.002)IHT(0.053)S(0.88)S(0.064)GDNEKPAVSSGK	C	C	C	ns
Uncharacterized protein	At5g64090	PL	Q9FMI8	8	MDFS(0.034)VKPS(0.772)GGS(0.159)PS(0.034)PSSSTSSSTPHR	C	C	C	1.6/ns

List of LC-MS/MS identified phosphosites regulated by OGs^a.

a

Phosphosites in each functional category are listed according to the average log2 ratio (Treated/Control), in descending order.

b

Full name of the identified protein. Proteins that have never been described before as phosphorylated are underlined. indicates proteins that have been identified as regulated by OGs in Kohorn et al. (2016).

c

ID of the identified protein from the TAIR database (The Arabidopsis Information Resource database. www.arabidopsis.org)

d

Subcellular localization, obtained from SUBA (the SUBcellular localization database for Arabidopsis proteins, http://suba3.plantenergy.uwa.edu.au/). PL, plastid; EX, extracellular; CY, cytosol; PM, plasma membrane; ER, endoplasmic reticulum; GO, Golgi apparatus; VO, vacuole; NU, nucleus; PX, peroxisome; MI, mitochondrion; CS, cytoskeleton. ^§indicates localization determined by GFP fusion.

e

ID of the identified protein, from the UniProtKB database (www.uniprot.org/). Proteins that have never been described before as phosphorylated at the residues reported here are ud.

f

Position of the phosphorylated amino acid on the whole protein sequence.

g

Phosphopeptide sequence with the localization probability score assigned to each site and the phosphorylated amino acid (bold and ud).

h

Average fold change (log2) between corresponding OG-treated and control phosphosites (Tr/C). Positive and negative values indicate increases and decreases in phosphorylation, respectively, upon treatment. C indicates proteins that are exclusively present in control samples and that are not detected in all replicates of OG- treated samples. Tr indicates proteins that are exclusively present in OG-treated samples and that are not detected in all replicates of control samples.

i

Fold change (log2) between corresponding OG-treated and control phosphosites (Tr/C) in each replicate.

j

Fold change of transcript levels after OG treatment at 1 and 3 h with respect to the mock-treated control (Denoux et al., 2008). Only data of genes for which fold change is significant and ≥1.5 in at least one time-point are shown. ns, not significant (P > 0.01); na, not applicable.

Proteomic analysis of the TMF did not show substantial changes in overall protein abundance (Data Sheet S2B in Supplementary Material). This is expected, as a significant de novo protein synthesis or protein degradation is unlikely to occur within 10 min of elicitor treatment (Benschop et al., 2007). Identified proteins (Data Sheet S2A in Supplementary Material) were classified according to known gene ontology using AgriGO (http://bioinfo.cau.edu.cn/agriGO/analysis.php). GO term analysis of the cellular component of the total list of proteins demonstrated a clear enrichment for PM-associated proteins, cell wall and intracellular organelles (Figure 1C).

Identification of OG-regulated phosphoproteins by phospho-DIGE

In a second approach, phospho-DIGE was exploited for detecting changes in protein phosphorylation using both TMFs and total protein extracts. Three independent biological replicates were performed for each treatment and samples were labeled according to the randomization scheme shown in Table S3, to enforce statistical analysis. About 800 fluorescent spots were detected for the total extracts, and about 500 spots for the TMFs. Statistical analysis revealed 30 spots in the total extract and 50 spots in the TMF fraction that changed their phosphorylation status (p < 0.05) after treatment with OGs (representative gel images are shown in Figures S1, S2). The spots were analyzed by MALDI-TOF-MS, leading to the identification of 46 differential phospho-isoforms, corresponding to 34 phosphoproteins (Table 2). Representative images showing spots that exhibited an increase or decrease in phosphorylation upon OG treatment, along with representative 3D images to visualize the phosphorylation changes, are shown in Figure 2. Proteomic analysis showed that, again, differences in non-phosphorylated protein abundance between OG-treated and control samples were not significant, demonstrating that most changes occurred at the PTM level.

Figure 2

Computational analysis of the differentially phosphorylated proteins

To summarize, LC-MS/MS analysis of TMF led to identification of 58 and 42 sites exhibiting an increased or decreased phosphorylation, respectively (Table 1), whereas Phospho-DIGE analysis of both TMF and total extracts revealed a total of 32 and 14 phosphoprotein spots showing an increased or decreased phosphorylation, respectively, corresponding to 34 proteins due to the presence of multispots (Table 2). Only Patellin 2 (PATL2) was found to be differentially phosphorylated in both analyses of the TMF. Two spots are shown for PATL2 (Figure 2 and Figure S2); one exhibits increased phosphorylation in response to OGs, in agreement with our detection of a phosphopeptide exclusively present in the OG-treated sample (Table 1). The other PATL2 spot shown in Figure 2 shows a slight but significant dephosphorylation (−1.15-fold). Many studies reported comparative analyses using both gel-based and gel-free methods, where the match of identified proteins is very low, and emphasized that these two strategies are indeed complementary (Kleffmann et al., 2007; Zhao et al., 2008; Robbins et al., 2013).

Previously described phosphorylation sites were searched using the PhosPhAt (http://phosphat.mpimp-golm.mpg.de/) and P3DB (http://www.p3db.org/) databases. Among all the differential phosphoproteins, sixteen [underlined in Table 1 (9 proteins) and Table 2 (7 proteins)] have never been described to be phosphorylated at the identified sites. GO term enrichment analysis with AgriGO was performed to identify significantly over-represented biological process categories in the subset (Tables 1, 2) of OG-regulated phosphoproteins (Figure 3). The GO terms referring to response to various biotic and abiotic stimuli were significantly enriched in the analyzed subsets of differential phosphoproteins, showing that the molecular processes involved in the response to OGs indeed correlate with a condition of stress. Proteins with increased phosphorylation were specifically enriched in terms referring to Signal transduction and to Cell death, while decreased phosphorylation occurred mainly in proteins associated with Response to biotic stimulus (Figure 3).

Figure 3

Connections between all the quantified phosphorylated proteins were searched by using the STRING algorithm (Figure 4). Proteins with differentially regulated phosphosites are highlighted with colors corresponding to different functional categories. Many clusters were identified, with several interesting networks standing out. The highly interconnected “transport” cluster includes several primary electrogenic proton pumps found in all eukaryotes: the PM H⁺-ATPases AHA1 and AHA2, different subunits of the vacuolar H⁺-ATPase (V-ATPase or VHA; subunits A, C, and E1; Batelli et al., 2007) as well as the subunit beta-1 of the mitochondrial ATP synthase. The V-ATPase subunit C, encoded by the single gene DE-ETIOLATED3 (DET3), is directly connected with the other V-ATPase subunit E1, but also with AHA1 and AHA2.

Figure 4

Multiple interactions are evident among elements involved in signaling (Figure 4), such as kinases, receptor kinases, and phosphatases. Many kinases appear to be connected to PP2C-type phosphatase-like protein (At2g20050), an enzyme that acts as negative modulator of protein kinase pathways involved in stress and developmental processes (Kuhn et al., 2006). Moreover, an additional smaller cluster includes several heat-shock proteins (HSP90-1, HSP70-1, HSP93-V, HSC70-1, BiP1) that are connected also to ribosomal proteins (Figure 4). Proteins involved in intracellular trafficking, namely syntaxin SYP 121 and 122, a calmodulin-binding protein (EDA39), PCaP1 (Plasma membrane-associated cation-binding protein 1), and a phospholipase like protein (PEARLI 4f) form an additional cluster.

Discussion

Previous proteomic studies have either postulated (Chivasa et al., 2006; Nuhse et al., 2007) or shown (Jones et al., 2006) that PTM-mediated regulation is important for plant immunity. Here we have studied very early protein dynamics induced by OGs, an important class of DAMPs, employing two complementary proteomic approaches, SIMAC phosphopeptide enrichment followed by LC-MS/MS and Phospho-DIGE (a gel-based approach). These analyses led to the identification of 100 regulated phosphosites and 46 differentially phosphorylated protein spots, respectively (Tables 1, 2). The functional classification of the OG-regulated phosphoproteins, based on the GO functional categories, identifies typical major PM functions. Three main categories emerged: transporter proteins, signaling proteins (receptors, kinases, phosphatases) and proteins involved in membrane trafficking.

A comparison of the OG-regulated phosphoproteome identified in this study with that regulated by flg22, which includes 127 phosphoproteins identified in three large-scale proteomics works (Nuhse et al., 2007; Benschop et al., 2007; Rayapuram et al., 2014), shows only 18 shared phosphosites, corresponding to 14 phosphoproteins (Table S4). Among these, there are MPK6, which shows OG-induced phosphorylation at the same residues that undergo flg22-induced phosphorylation (i.e., Thr221 and Tyr223), RBOHD and Penetration3 (PEN3), an ATP-binding cassette (ABC) family 36 transporter that undergoes phosphorylation also in response to H₂O₂ and methyl jasmonate (Stecker et al., 2014). ABC transporters are integral membrane proteins that transport a wide variety of substrates, such as ABA, auxin, and some plant secondary metabolites across cellular membranes (Kuromori et al., 2010).

Interestingly, twenty-two of the OG-regulated phosphoproteins identified here are reported as targets of MPK3 and MPK6 (Table S5; Lee et al., 2015). Among these, the putative tyrosine phosphatase PTEN2a shows OG-dependent increased phosphorylation. This enzyme has been shown to dephosphorylate also the 3′ phosphate group of PI3P (phosphatidylinositol 3-phosphate) in vitro and to possess strong binding affinity for PA (phosphatidic acid), an important second messenger with roles in stress and hormonal signaling (Pribat et al., 2012). Lipid signaling pathways activated downstream of the MAPKs may therefore be involved in the response to OGs. A phosphotyrosine site that disappears in response to OGs and may be a target of PTEN2a is present in the GSK/Shaggy-like kinase/ASK1 at position 229. Auto-phosphorylation of this residue is required for ASK1 trans-phosphorylation activity (de la Fuente van Bentem et al., 2008) and activation of the cytosolic glucose-6-phosphate dehydrogenase (G6PD), which is essential for maintaining the cellular redox balance (Dal Santo et al., 2012). ASK1 dephosphorylation may thus lead to reduced G6PD activity and elevated ROS levels that contribute to the elicitor-induced oxidative burst.

MPK3 and MPK6 activation promotes also the biosynthesis of indole glucosinolates and the conversion of indole-3-yl-methylglucosinolate (I3G) to 4-methoxyindole-3-yl-methylglucosinolate (4MI3G), through the phosphorylation of ETHYLENE RESPONSE FACTOR 6 (ERF6). 4MI3G is a substrate of the atypical myrosinases PEN2 and PEN3, which shows increased phosphorylation in response to OGs, and is converted to extracellular unstable anti-microbial compounds, possibly isothiocyanates (Xu et al., 2016). OGs may therefore induce production of indole glucosinolates, but this aspect has been scarcely investigated so far.

A few differentially phosphorylated proteins identified here show regulation also at the level of gene expression, with more than a 1.5-fold change in the corresponding transcript levels within 3 h (Tables 1, 2; Denoux et al., 2008), pointing to a regulation at two different levels. Most of the identified proteins, however, show mainly a PTM-based regulation, by phosphorylation/dephosphorylation, suggesting that only activation/inactivation or stabilization/degradation, but not de novo protein synthesis, is necessary to exert their physiological role.

Signal transduction proteins

Several identified proteins are likely players in signaling. Among the proteins involved in regulation through phosphorylation/dephosphorylation, we found members of the Receptor-Like Kinases (RLK) family, including members of the Leucine-Rich Repeat RLK (LRR-RLK) sub-family, cyclin-dependent protein kinases, PTEN2a, and the MAP triple kinase MKKK7 (At3g13530, also known as MAP3Kepsilon 1). MKKK7 has been described to interact with FLS2 and to be phosphorylated at Ser452/Ser854 in response to flg22, functioning as a negative regulator of flg22-triggered MPK6 activation (Mithoe et al., 2016). We find, instead, dephosphorylation of MKKK7 on Ser788 after OG perception, suggesting a distinct event in the signaling cascades initiated by the two elicitors. Another immune-related protein, the Cysteine-rich Receptor-like protein Kinase 36 (CRK36), exhibits decreased OG-dependent phosphorylation, along with up-regulation at the transcript level. CRK36 interacts with FLS2 and its overexpression leads to enhanced PTI and resistance to virulent bacteria Pseudomonas syringae pv. tomato DC3000, constitutive accumulation of callose and constitutive stomatal closure (Yeh et al., 2015). A role of CRK36 in resistance to insects has also been proposed (Barah et al., 2013). Moreover, CRK36 silencing leads to higher sensitivity to ABA and osmotic stress during the post-germination growth phase, suggesting that the protein negatively controls ABA signaling (Tanaka et al., 2012).

We also identified three phospholipases (PLDγ1 and two putative phospholipases, pEARLI4, and pEARLI4-f) as well as CBL9, a myristoylated PM-localized protein that does not possess enzymatic activity and contains four EF hands (Batistic et al., 2008). CBL9 modulates ABA responses (Kim et al., 2003) and forms complexes with CBL-interacting protein kinases (CIPKs), leading to phosphorylation of AHA2 (Fuglsang et al., 1999) and regulation of V-ATPases (Tang et al., 2012), also found among the OG-regulated phosphoproteins. CBL-CIPK complexes act as a two-module Ca²⁺⁻-decoding system where CBL phosphorylation is required for CIPK-mediated phosphorylation of target proteins (Steinhorst and Kudla, 2013) and regulation of channel activity, ion fluxes and ROS formation during environmental adaptation reactions (Kurusu et al., 2015).

It is worth noting that remorin REM1.3, identified as phosphorylated on Thr58 in response to both OGs and flg22 (Kohorn et al., 2016; Benschop et al., 2007), was also identified in our analysis as phosphorylated at the same site, but with a Significance B p-value of 0.1. A tomato remorin was the first protein described to bind and to be phosphorylated in response to OGs (Reymond et al., 1996) and potato REM1.3 has been shown to accumulate in discrete peri-haustorial domains and enhance susceptibility to Phytophthora infestans (Bozkurt et al., 2014). Arabidopsis REM1.3 is differentially recruited to detergent-insoluble membranes/detergent-resistant membranes (DIMs or DRMs), a term that indicates lipids and proteins that have been biochemically associated with sphingolipid (SL)- and sterol-enriched membranes (Tapken and Murphy, 2015).

Two proteins functioning in the transduction of the light signal, PHOT2 and Phytochrome Kinase Substrate 2 (PKS2), exhibit increased phosphorylation in response to OGs. pks2 loss-of-function mutants show altered expression patterns of auxin marker genes (Kami et al., 2014) and it is tempting to speculate that PKS2 is involved in auxin OG/antagonism, the molecular basis of which is still unknown (Savatin et al., 2011). Notably, the auxin-response factor ARF1, a repressor ARF (Chandler, 2016), showed decreased phosphorylation upon OG treatment at two novel phophosites, the adjacent residues Ser399 and Ser400 located in middle region of the protein thought to function as the repression domain (Ulmasov et al., 1999). Ser405, close to this novel dual phosphorylation site, has been described as phosphorylated in seedlings (Wang et al., 2013) but does not appear to be, regulated by OGs. Very little is known about how phosphorylation controls ARF activity. ARF2, also a repressor ARF that is phylogenetically sister to ARF1, loses its DNA-binding capability upon phosphorylation by the brassinosteroid-regulated kinase BIN2 (BRASSINOSTEROID-INSENSITIVE2), providing an example of hormonal crosstalk at the post-translational level (Vert et al., 2008). On the other hand, phosphorylation of the activator ARFs ARF7 and ARF19 increases their DNA binding ability and enhances transcription (Cho et al., 2014). In the light of the OG-auxin antagonism, OG-induced phosphorylation of ARF1 may lead to its activation as a repressor.

Trafficking-related proteins

OGs induce increased phosphorylation of cytoskeleton-related phosphoproteins such as myosin-17 (XIK), Actin7 (ACT7), and PCaP1. Myosin XIK is required for normal movements of actin filaments (Park and Nebenfuhr, 2013), turnover of actin and trafficking of vesicles carrying non-cellulosic cell wall components (Cai et al., 2014). The protein acts as a key regulator of plant antifungal immunity and contributes to cell wall-mediated penetration resistance, including deposition of callose and accumulation of PENETRATION1 (PEN1)/SYP121 and lignin-like compounds at the penetration site, when highly expressed (Yang et al., 2014). Pharmacological inhibition of myosins affects accumulation of PEN3 more than that of PEN1, suggesting that transport pathways mediating recruitment and export of the two proteins to apoplastic papillae are distinct (Underwood and Somerville, 2013). OG-induced phosphorylation of myosin XIK may regulate the motility of secretory vesicles and PEN recruitment.

ACT7 is one of the three vegetative and partially redundant isoforms of the actin family and plays a role in germination and root growth (Gilliland et al., 2003; Kijima et al., 2015). Expression of ACT7 strongly responds to auxin and wounding and is required for callus formation (Kandasamy et al., 2001); as shown in Figure 4, it is also coordinated with the expression of GAPC1 (Glyceraldehyde-3-phosphate dehydrogenase), the product of which also shows increased phosphorylation in response to OGs. How phosphorylation regulates actin function in plants is not known, except for Mimosa pudica actin, which is heavily tyrosine-phosphorylated, and dephosphorylated at different extent depending on the degree of bending of the petiols (Kameyama et al., 2000).

PCaP1 is capable of binding actin and calcium, leading to the destabilization of actin filaments, besides phosphatidylinositol phosphates (Qin et al., 2014). It is localized via N-myristoylation on the PM, from which it dissociates at increased Ca²⁺ levels (Li et al., 2011). The protein negatively regulates hypocotyl (Li et al., 2011) and pollen growth (Qin et al., 2014); moreover, it is induced by brassinosteroids and its overexpression partially rescues the morphological phenotype of the det2 mutant, defective in BR biosynthesis (Tang et al., 2008). PCap1 has never been implicated in immunity.

Components of the microtubule (MT) apparatus, namely TUBULIN ALPHA-5 (TUA5) and the 65-kDa MT-associated protein 1 (MAP65-1), a structural element involved in MT bundling and cell division, also show OG-regulated increase in phosphorylation. MTs are disrupted or rearranged in response to elicitors and stable MTs negatively affect defense (Guan et al., 2013). Phosphorylation of MAP65-1 was found exclusively in OG-treated samples, at residues located within the second MT-interacting region. Phosphorylation of this region has been shown to diminish the interaction with MTs, also in tobacco (Guan et al., 2013) and involve CDPKs and MAPKs (Smertenko et al., 2006), which are known to participate in OG signaling cascade (Ferrari et al., 2013; Gravino et al., 2015).

Vesicle trafficking appears to be another important target of the OG regulatory action. SYP122 decreases phosphorylation in response to OGs, while SYP121/PEN1 (Reichardt et al., 2011) as well as PATL2, PATL3, and PATL4 show an opposite behavior. PATL1 is essential for cell plate formation and maturation at late telophase, when dynamic vesicle trafficking and vesicle fusion occur (Peterman et al., 2004); PATL2 is related to the phosphatydylinositol transfer protein Sec14, which plays a role in the interplay between lipid metabolism and membrane trafficking, whereas PATL3 has been shown to interact with the movement protein of Alfalfa Mosaic Virus and interfere with viral movement (Peiro et al., 2014).

A phosphosite of the kinesin motor protein KAC1 (At5g10470) was found only in the presence of OGs. KAC1 is phoshorylated also in response to flg22 (Table S4) and xylanase treatment (Benschop et al., 2007). Kinesins function in vesicle transport along MTs and their phosphorylation is thought to be required for correct folding and dimerization. The C-terminal domain of KAC1 and KAC2, interact with F-actin in vitro and mediate chloroplast movement (Suetsugu et al., 2010). Both proteins belong to the kinesin-14 subfamily (Zhu and Dixit, 2012), greatly expanded in land plants and representing the largest Arabidopsis kinesin subfamily (21 members). Expansion is attributed to the lack of the MT-based motor protein dynein in plants, and some of the kinesin-14 subfamily members may substitute for dynein to perform retrograde transport functions (Muller, 2015).

ATPases and other membrane transporters

Initial pathogen recognition occurs at the PM and many of the early responses involve membrane transport processes. For example, activity of PM AHAs is dynamically regulated during plant immune responses and multiple pathogens target this family of enzymes. AHAs are responsible for creating and maintaining a negative membrane potential and a transmembrane pH gradient (acidic outside) that control multiple aspects of transport across the PM. Regulatory phosphorylation sites identified in AHAs are located either in the N- or C-terminal domain (Rudashevskaya et al., 2012). Phosphorylation of the penultimate residue, a Thr, within the auto-inhibitory C-terminal domain is critical for interaction with a 14-3-3 regulatory protein (Fuglsang et al., 1999); moreover, phosphorylation within this domain affects enzyme activity incrementally, suggesting a fine modulation (Speth et al., 2010). OGs induce phosphorylation of AHA1 and AHA2 at the conserved residue Ser899, a modification that has been reported to inhibit proton transport (Haruta et al., 2014) and to occur in response to flg22 (Nuhse et al., 2007) This modification may play a role in the extracellular alkalinization induced by elicitors, including OGs. Moreover, OGs induce phosphorylation of AHA1 at the non-conserved Ser544, located in the central part of the protein. Phosphorylation at this site has been recently described (Rudashevskaya et al., 2012) but, to our knowledge, never in response to elicitors.

Several subunits of V-ATPase (VHA-A, VHA-E1, and DET3) also change their phosphorylation status in response to OGs. V-ATPases are multi-subunit enzymes resembling the procariotic ATPases and therefore performing rotational catalysis. V-ATPases are localized on both the tonoplast and the Trans-Golgi Network (TGN) and their impairement affects morphology and function of the Golgi and TGN (Schumacher and Krebs, 2010). In both plants and animals (Robinson, 2014), V-ATPase perform uphill transport of protons and do not only energize secondary active transport, but also play a relevant role in vesicle trafficking and post-translational processing of secretory proteins, by regulating endosomal pH (Huotari and Helenius, 2011). Micro-compartments along the endocytic and secretory pathways have characteristic luminal pH values suited to their functions, and specific luminal pH and ionic concentration may be required not only for appropriate enzyme activities but also for membrane trafficking and sorting and cargo routing. Progressively decreasing pH in the endocytic pathway are thought to provide the cargo information about its location within the pathway (Huotari and Helenius, 2011; Martiniere et al., 2013).

OGs regulates phosphorylation of other plasma membrane transporters such as the acquaporin PIP2-1 and, besides PEN3, two additional ABC transporters, At3g55320 and At3g62700, belonging to subfamily 20 and 14, respectively. Ser280, Thr712 and Ser925 in PIP2-1, At3g55320, and At3g62700, respectively, showed phosphorylation only in the presence of OGs. In PIP2-1, phosphorylation at Ser280 and Ser283, has been associated to the regulation of hydraulic conductivity (Prado et al., 2013). Moreover, the latter phosphosites is necessary for correct targeting to the PM (Prak et al., 2008), showing that specific phosphorylation events control not only activity, but also targeting of transporters. Interestingly, PIP2-1 along with some members of family B and G ABC transporters not identified in this study, are markers of the so-called membrane nanodomains (NDs), ordered membrane SL- and sterol-rich microenvironments with functional assemblies of lipids and proteins that are defined through biophysical or microscopy, and not biochemical, techniques (Tapken and Murphy, 2015), and therefore are distinct from DIMs/DRMs, where PEN3 is instead found.

Response to stress

A small cluster of OG-regulated phosphoproteins is represented by heat-shock proteins, including the cytosolic heat-shock cognate 70-1(HSC70-1), the chloroplast heat shock protein 70-1 (HSP70-1), BiP1 and the cytosolic HSP90. HsP90 is a molecular chaperone responsible for the maturation and stability of a large number of signaling proteins, in particular protein kinases. We found decreased levels of the HSP90 Ser219 phosphosite, in agreement with previous reports (Reiland et al., 2009; Jones et al., 2009). The role of this modification is not known. In yeast, key phosphorylation sites on HSP90 affect the interaction with a selected subset of co-chaperones, and consequently, client proteins (Mollapour et al., 2011); a similar role is possible in plants. Arabidopsis HSP90 interacts with MPK4 and regulates its kinase activity during response to flg22 (Cui et al., 2010). HSP90 isoforms, together with the co-chaperones RAR1 and SGT1, appear to be critical also for the stability of nucleotide-binding leucine-rich repeat receptor (NLR) resistance proteins such as RPM1, MLA, RPS2, RPS4, I2, and N, and the formation of immune receptor complexes (Huang et al., 2014). In addition, they negatively regulate the accumulation of NLRs, likely to avoid autoimmunity (Huang et al., 2014). The decreased level of phosphorylation in response to OGs is likely related to the regulation of HSP90 complex assembly and substrate specificity.

Proteins involved in redox homeostasis

Many proteins involved in redox homeostasis and ROS signaling show OG-induced increase in phosphorylation, including an isocitrate dehydrogenase (CICDH; Table 2) involved in response to pathogens (Mhamdi et al., 2010), and GAPC1, a glycolytic enzyme that interacts with phospholipase D and transduces hydrogen peroxide signals during stress (Guo et al., 2012). Both proteins appear to be regulated mainly through PTM. In fact, only few proteins involved in redox regulation here identified show both phosphorylation changes and transcript up-regulation. Among these, RBOHD (Kadota et al., 2014), glutathione S-transferases (GSTF9, GST-PM24, GSTF8), the non-photosynthetic NADP-malic enzyme 2 (NADP-ME2) and flavodoxin-like quinone reductase 1 (FQR1) exhibit increased phosphorylation, whereas a monodehydroascorbate reductase (MDAR2) and the protein plant cadmium resistance 8 (CdRes8) show decreased phosphorylation. NADP-ME2 plays a role in basal defense against the hemibiotrophic fungal pathogen Colletotrichum higginsianum, and is involved in the generation of ROS (Voll et al., 2012), whereas FQR1 is a primary auxin-response gene (Laskowski et al., 2002) that belongs to quinone reductases, classified as phase II detoxification enzymes, that protect organisms from oxidative stress; interestingly, FQR1 acts as a susceptibility factor to B. cinerea (Heyno et al., 2013).

PTM-mediated regulation of oxidative-stress related proteins is in agreement with the notion that the redox state of the cell is dramatically altered in response to OGs.

Other proteins

JAL34 and JAL27 (Table 2) are jacalin lectin-like proteins similar to myrosinase-binding protein 1 (Yamada et al., 2011) and to the functionally related PYK10. The latter is a β-glucosidase localized in the endoplasmic reticulum (ER) body, a large compartment specific to the Brassicales. PYK10 forms complexes with JALs and other proteins when tissue is damaged, i.e., upon herbivore or pathogen attack. Engagement in a complex may shield active PYK10 from inhibitors, proteases and other proteins that may reduce its defensive effect (Yamada et al., 2011). JAL34 and JAL27 show also OG-induced increased abundance at a later time point (Casasoli et al., 2008). JAL34 and JAL27 as well as PCaP1, ACT7, PATL2, and MDAR2 have been described in a previous proteomic work as proteins induced by brassinosteroids (Tang et al., 2008); whether these proteins are involved in the BR-mediated inhibition of PTI is still unknown.

Conclusions

In conclusion, we have provided a large-scale study of early phosphoproteome changes in Arabidopsis following OG perception. In a recent phosphoproteomic study of Arabidopsis seedling total extracts following OG treatment, 51 phosphorylated peptides, representing 49 unique proteins, were identified (Kohorn et al., 2016). A comparison with our study shows only 9 shared phosphosites, corresponding to 9 proteins (Table 1), two of which also found in flg22-regulated phosphoproteome. Taking into consideration the other 16 shared flg22-regulated phosphosites, our work performed also on membrane proteins uncovers the immune-related regulation of 73 phosphosites. Protein interaction network analyses point to the main biological processes in which protein phosphorylation events are crucial in OG signaling, and suggest the interplay of several processes, e.g., intracellular trafficking and vesicle dynamics, cytoskeleton rearrangement, signal transduction and phospholipid signaling.

Funding

This work was supported by the European Research Council (ERC Advanced Grant233083), the Italian Ministry of Agriculture Food and Forestry Policies (MIPAAF, Project ALISAL, DM 11008/7303/10) and University of Rome Sapienza (ATENEO, 2014).

Conflict of interest statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

• 1

  AlexanderssonE.SaalbachG.LarssonC.KjellbomP. (2004). Arabidopsis plasma membrane proteomics identifies components of transport, signal transduction and membrane trafficking. Plant Cell Physiol.45, 1543–1556. 10.1093/pcp/pch209

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 2

  AzizA.HeyraudA.LambertB. (2004). Oligogalacturonide signal transduction, induction of defense-related responses and protection of grapevine against Botrytis cinerea. Planta218, 767–774. 10.1007/s00425-003-1153-x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 3

  BarahP.WingeP.KusnierczykA.TranD. H.BonesA. M. (2013). Molecular signatures in Arabidopsis thaliana in response to insect attack and bacterial infection. PLoS ONE8:e58987. 10.1371/journal.pone.0058987

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4

  BarrettL. G.HeilM. (2012). Unifying concepts and mechanisms in the specificity of plant-enemy interactions. Trends Plant Sci.17, 282–292. 10.1016/j.tplants.2012.02.009

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5

  BatelliG.VersluesP. E.AgiusF.QiuQ.FujiiH.PanS.et al. (2007). SOS2 promotes salt tolerance in part by interacting with the vacuolar H+-ATPase and upregulating its transport activity. Mol. Cell. Biol.27, 7781–7790. 10.1128/MCB.00430-07

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 6

  BatisticO.SorekN.SchultkeS.YalovskyS.KudlaJ. (2008). Dual fatty acyl modification determines the localization and plasma membrane targeting of CBL/CIPK Ca2+ signaling complexes in Arabidopsis. Plant Cell20, 1346–1362. 10.1105/tpc.108.058123

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7

  BellincampiD.CervoneF.LionettiV. (2014). Plant cell wall dynamics and wall-related susceptibility in plant-pathogen interactions. Front. Plant Sci.5:228. 10.3389/fpls.2014.00228

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 8

  BenedettiM.PontiggiaD.RaggiS.ChengZ. Y.ScaloniF.FerrariS.et al. (2015). Plant immunity triggered by engineered in vivo release of oligogalacturonides, damage-associated molecular patterns. Proc. Natl. Acad. Sci. U.S.A.112, 5533–5538. 10.1073/pnas.1504154112

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 9

  BenschopJ. J.MohammedS.O'FlahertyM.HeckA. J.SlijperM.MenkeF. L. (2007). Quantitative phosphoproteomics of early elicitor signaling in Arabidopsis. Mol. Cell. Proteomics.6, 1198–1214. 10.1074/mcp.M600429-MCP200

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 10

  BollerT.FelixG. (2009). A renaissance of elicitors: perception of microbe-associated molecular patterns and danger signals by pattern-recognition receptors. Annu. Rev. Plant Biol.60, 379–406. 10.1146/annurev.arplant.57.032905.105346

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 11

  BondA. E.RowP. E.DudleyE. (2011). Post-translation modification of proteins; methodologies and applications in plant sciences. Phytochemistry72, 975–996. 10.1016/j.phytochem.2011.01.029

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 12

  BozkurtT. O.RichardsonA.DagdasY. F.MongrandS.KamounS.RaffaeleS. (2014). The plant membrane-associated REMORIN1.3 accumulates in discrete perihaustorial domains and enhances susceptibility to Phytophthora infestans. Plant Physiol.165, 1005–1018. 10.1104/pp.114.235804

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 13

  BrutusA.SiciliaF.MaconeA.CervoneF.De LorenzoG. (2010). A domain swap approach reveals a role of the plant wall-associated kinase 1 (WAK1) as a receptor of oligogalacturonides. Proc. Natl. Acad. Sci. U.S.A.107, 9452–9457. 10.1073/pnas.1000675107

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 14

  CaiC.Henty-RidillaJ. L.SzymanskiD. B.StaigerC. J. (2014). Arabidopsis myosin XI: a motor rules the tracks. Plant Physiol.166, 1359–1370. 10.1104/pp.114.244335

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15

  CasasoliM.FedericiL.SpinelliF.Di MatteoA.VellaN.ScaloniF.et al. (2009). Integration of evolutionary and desolvation energy analysis identifies functional sites in a plant immunity protein. Proc. Natl. Acad. Sci. U.S.A.106, 7666–7671. 10.1073/pnas.0812625106

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 16

  CasasoliM.MelicianiI.CervoneF.De LorenzoG.MatteiB. (2007). Oligogalacturonide-induced changes in the nuclear proteome of Arabidopsis thaliana. Int. J. Mass Spectrom.268, 277–283. 10.1016/j.ijms.2007.07.007

  • CrossRef
  • Google Scholar

• 17

  CasasoliM.SpadoniS.LilleyK. S.CervoneF.De LorenzoG.MatteiB. (2008). Identification by 2-D DIGE of apoplastic proteins regulated by oligogalacturonides in Arabidopsis thaliana. Proteomics8, 1042–1054. 10.1002/pmic.200700523

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18

  ChandlerJ. W. (2016). Auxin response factors. Plant Cell Environ.39, 1014–1028. 10.1111/pce.12662

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 19

  ChisholmS. T.CoakerG.DayB.StaskawiczB. J. (2006). Host-microbe interactions: shaping the evolution of the plant immune response. Cell124, 803–814. 10.1016/j.cell.2006.02.008

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 20

  ChivasaS.HamiltonJ. M.PringleR. S.NdimbaB. K.SimonW. J.LindseyK.et al. (2006). Proteomic analysis of differentially expressed proteins in fungal elicitor-treated Arabidopsis cell cultures. J. Exp. Bot.57, 1553–1562. 10.1093/jxb/erj149

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 21

  ChoH.RyuH.RhoS.HillK.SmithS.AudenaertD.et al. (2014). A secreted peptide acts on BIN2-mediated phosphorylation of ARFs to potentiate auxin response during lateral root development. Nat. Cell Biol.16, 66–76. 10.1038/ncb2893

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 22

  CoxJ.HeinM. Y.LuberC. A.ParonI.NagarajN.MannM. (2014). Accurate proteome-wide label-free quantification by delayed normalization and maximal peptide ratio extraction, termed MaxLFQ. Mol. Cell. Proteomics13, 2513–2526. 10.1074/mcp.M113.031591

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 23

  CoxJ.MannM. (2008). MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat. Biotechnol.26, 1367–1372. 10.1038/nbt.1511

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 24

  CuiH. T.WangY. J.XueL.ChuJ. F.YanC. Y.FuJ. H.et al. (2010). Pseudomonas syringae effector protein avrb perturbs arabidopsis hormone signaling by activating MAP kinase 4. Cell Host Microbe7, 164–175. 10.1016/j.chom.2010.01.009

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 25

  Dal SantoS.StampflH.KrasenskyJ.KempaS.GibonY.PetutschnigE.et al. (2012). Stress-induced GSK3 regulates the redox stress response by phosphorylating glucose-6-phosphate dehydrogenase in Arabidopsis. Plant Cell24, 3380–3392. 10.1105/tpc.112.101279

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 26

  de la Fuente van BentemS.AnratherD.DohnalI.RoitingerE.CsaszarE.JooreJ.et al. (2008). Site-specific phosphorylation profiling of Arabidopsis proteins by mass spectrometry and peptide chip analysis. J. Proteome Res.7, 2458–2470. 10.1021/pr8000173

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 27

  De LorenzoG.BrutusA.SavatinD. V.SiciliaF.CervoneF. (2011). Engineering plant resistance by constructing chimeric receptors that recognize damage-associated molecular patterns (DAMPs). FEBS Lett.585, 1521–1528. 10.1016/j.febslet.2011.04.043

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 28

  DenouxC.GallettiR.MammarellaN.GopalanS.WerckD.De LorenzoG.et al. (2008). Activation of defense response pathways by OGs and Flg22 elicitors in Arabidopsis seedlings. Mol. Plant1, 423–445. 10.1093/mp/ssn019

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 29

  DurekP.SchmidtR.HeazlewoodJ. L.JonesA.MacLeanD.NagelA.et al. (2010). PhosPhAt: the Arabidopsis thaliana phosphorylation site database. An update. Nucleic Acids Res.38, D828–D834. 10.1093/nar/gkp810

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 30

  FerrariS.GallettiR.DenouxC.De LorenzoG.AusubelF. M.DewdneyJ. (2007). Resistance to Botrytis cinerea induced in Arabidopsis by elicitors is independent of salicylic acid, ethylene, or jasmonate signaling but requires PHYTOALEXIN DEFICIENT3. Plant Physiol.144, 367–379. 10.1104/pp.107.095596

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 31

  FerrariS.SavatinD. V.SiciliaF.GramegnaG.CervoneF.De LorenzoG. (2013). Oligogalacturonides: plant damage-associated molecular patterns and regulators of growth and development. Front.Plant Sci.10.3389/fpls.2013.00049

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 32

  FuglsangA. T.ViscontiS.DrummK.JahnT.StensballeA.MatteiB.et al. (1999). Binding of 14-3-3 protein to the plasma membrane H+-ATPase AHA2 involves the three C-terminal residues Tyr(946)-Thr-Val and requires phosphorylation of Thr(947). J. Biol. Chem.274, 36774–36780. 10.1074/jbc.274.51.36774

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 33

  GallettiR.DenouxC.GambettaS.DewdneyJ.AusubelF. M.De LorenzoG.et al. (2008). The AtrbohD-mediated oxidative burst elicited by oligogalacturonides in Arabidopsis is dispensable for the activation of defense responses effective against Botrytis cinerea. Plant Physiol.148, 1695–1706. 10.1104/pp.108.127845

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 34

  GallettiR.FerrariS.De LorenzoG. (2011). Arabidopsis MPK3 and MPK6 play different roles in basal and oligogalacturonide- or flagellin-induced resistance against Botrytis cinerea. Plant Physiol.157, 804–814. 10.1104/pp.111.174003

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35

  GillilandL. U.PawloskiL. C.KandasamyM. K.MeagherR. B. (2003). Arabidopsis actin gene ACT7 plays an essential role in germination and root growth. Plant J.33, 319–328. 10.1046/j.1365-313X.2003.01626.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36

  GramegnaG.ModestiV.SavatinD. V.SiciliaF.CervoneF.De LorenzoG. (2016). GRP-3 and KAPP, encoding interactors of WAK1, negatively affect defense responses induced by oligogalacturonides and local response to wounding. J. Exp. Bot.67, 1715–1729. 10.1093/jxb/erv563

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37

  GravinoM.LocciF.TundoS.CervoneF.SavatinD. V.De LorenzoG. (2016). Immune responses induced by oligogalacturonides are differentially affected by AvrPto and loss of BAK1/BKK1 and PEPR1/PEPR2. Mol. Plant Pathol. [Epub ahead of print]. 10.1111/mpp.12419

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 38

  GravinoM.SavatinD. V.MaconeA.De LorenzoG. (2015). Ethylene production in Botrytis cinerea- and oligogalacturonide-induced immunity requires calcium-dependent protein kinases. Plant J.84, 1073–1086. 10.1111/tpj.13057

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 39

  GuanX.BuchholzG.NickP. (2013). The cytoskeleton is disrupted by the bacterial effector HrpZ, but not by the bacterial PAMP flg22, in tobacco BY-2 cells. J. Exp. Bot.64, 1805–1816. 10.1093/jxb/ert042

  • CrossRef
  • Google Scholar

• 40

  GuoL.DevaiahS. P.NarasimhanR.PanX.ZhangY.ZhangW.et al. (2012). Cytosolic glyceraldehyde-3-phosphate dehydrogenases interact with phospholipase Ddelta to transduce hydrogen peroxide signals in the Arabidopsis response to stress. Plant Cell24, 2200–2212. 10.1105/tpc.111.094946

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41

  HarutaM.SabatG.SteckerK.MinkoffB. B.SussmanM. R. (2014). A peptide hormone and its receptor protein kinase regulate plant cell expansion. Science343, 408–411. 10.1126/science.1244454

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 42

  HeazlewoodJ. L.Tonti-FilippiniJ.VerboomR. E.MillarA. H. (2005). Combining experimental and predicted datasets for determination of the subcellular location of proteins in Arabidopsis. Plant Physiol.139, 598–609. 10.1104/pp.105.065532

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43

  HeynoE.AlkanN.FluhrR. (2013). A dual role for plant quinone reductases in host-fungus interaction. Physiol. Plant.149, 340–353. 10.1111/ppl.12042

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 44

  HuangS.MonaghanJ.ZhongX. H.LinL.SunT. J.DongO. X.et al. (2014). HSP90s are required for NLR immune receptor accumulation in Arabidopsis. Plant J.79, 427–439. 10.1111/tpj.12573

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 45

  HuotariJ.HeleniusA. (2011). Endosome maturation. EMBO J.30, 3481–3500. 10.1038/emboj.2011.286

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 46

  JonesA. M. E.MacLeanD.StudholmeD. J.Serna-SanzA.AndreassonE.RathjenJ. P.et al. (2009). Phosphoproteomic analysis of nuclei-enriched fractions from Arabidopsis thaliana. J. Proteomics72, 439–451. 10.1016/j.jprot.2009.02.004

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 47

  JonesA. M.ThomasV.BennettM. H.MansfieldJ.GrantM. (2006). Modifications to the Arabidopsis defense proteome occur prior to significant transcriptional change in response to inoculation with Pseudomonas syringae. Plant Physiol.142, 1603–1620. 10.1104/pp.106.086231

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 48

  KadotaY.SklenarJ.DerbyshireP.StransfeldL.AsaiS.NtoukakisV.et al. (2014). Direct regulation of the NADPH oxidase RBOHD by the PRR-associated kinase BIK1 during plant immunity. Mol. Cell54, 43–55. 10.1016/j.molcel.2014.02.021

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 49

  KalunkeR. M.TundoS.BenedettiM.CervoneF.De LorenzoG.D'OvidioR. (2015). An update on polygalacturonase-inhibiting protein (PGIP), a leucine-rich repeat protein that protects crop plants against pathogens. Front. Plant Sci.6:146. 10.3389/fpls.2015.00146

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 50

  KameyamaK.KishiY.YoshimuraM.KanzawaN.SameshimaM.TsuchiyaT. (2000). Tyrosine phosphorylation in plant bending. Nature407, 37. 10.1038/35000166

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 51

  KamiC.AllenbachL.ZourelidouM.LjungK.SchutzF.IsonoE.et al. (2014). Reduced phototropism in pks mutants may be due to altered auxin-regulated gene expression or reduced lateral auxin transport. Plant J.77, 393–403. 10.1111/tpj.12395

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 52

  KandasamyM. K.GillilandL. U.McKinneyE. C.MeagherR. B. (2001). One plant actin isovariant, ACT7, is induced by auxin and required for normal callus formation. Plant Cell13, 1541–1554. 10.1105/tpc.13.7.1541

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 53

  KijimaS. T.HiroseK.KongS. G.WadaM.UyedaT. Q. (2015). Distinct biochemical properties of Arabidopsis thaliana actin isoforms. Plant Cell Physiol.57, 46–56. 10.1093/pcp/pcv176

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 54

  KimK. N.CheongY. H.GrantJ. J.PandeyG. K.LuanS. (2003). CIPK3, a calcium sensor-associated protein kinase that regulates abscisic acid and cold signal transduction in Arabidopsis. Plant Cell15, 411–423. 10.1105/tpc.006858

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 55

  KleffmannT.von ZychlinskiA.RussenbergerD.Hirsch-HoffmannM.GehrigP.GruissemW.et al. (2007). Proteome dynamics during plastid differentiation in rice. Plant Physiol143, 912–923. 10.1104/pp.106.090738

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 56

  KohornB. D.HoonD.MinkoffB. B.SussmanM. R.KohornS. L. (2016). Rapid oligo-galacturonide induced changes in protein phosphorylation in Arabidopsis. Mol. Cell Proteomics15, 1351–1359. 10.1074/mcp.m115.055368

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 57

  KohornB. D.KohornS. L. (2012). The cell wall-associated kinases, WAKs, as pectin receptors. Front. Plant Sci.3:88. 10.3389/fpls.2012.00088

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 58

  KuhnJ. M.Boisson-DernierA.DizonM. B.MaktabiM. H.SchroederJ. I. (2006). The protein phosphatase AtPP2CA negatively regulates abscisic acid signal transduction in Arabidopsis, and effects of abh1 on AtPP2CA mRNA. Plant Physiol.140, 127–139. 10.1104/pp.105.070318

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 59

  KuromoriT.MiyajiT.YabuuchiH.ShimizuH.SugimotoE.KamiyaA.et al. (2010). ABC transporter AtABCG25 is involved in abscisic acid transport and responses. Proc. Natl. Acad. Sci. U.S.A.107, 2361–2366. 10.1073/pnas.0912516107

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 60

  KurusuT.KuchitsuK.TadaY. (2015). Plant signaling networks involving Ca(2+) and Rboh/Nox-mediated ROS production under salinity stress. Front. Plant Sci.6:427. 10.3389/fpls.2015.00427

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 61

  LaskowskiM. J.DreherK. A.GehringM. A.AbelS.GenslerA. L.SussexI. M. (2002). FQR1, a novel primary auxin-response gene, encodes a flavin mononucleotide-binding quinone reductase. Plant Physiol.128, 578–590. 10.1104/pp.010581

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 62

  LeeJ.Eschen-LippoldL.LassowskatI.BottcherC.ScheelD. (2015). Cellular reprogramming through mitogen-activated protein kinases. Front. Plant Sci.6:940. 10.3389/fpls.2015.00940

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 63

  LiJ.WangX.QinT.ZhangY.LiuX.SunJ.et al. (2011). MDP25, a novel calcium regulatory protein, mediates hypocotyl cell elongation by destabilizing cortical microtubules in Arabidopsis. Plant Cell23, 4411–4427. 10.1105/tpc.111.092684

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 64

  LiuZ.LiY.CaoH.RenD. (2015). Comparative phospho-proteomics analysis of salt-responsive phosphoproteins regulated by the MKK9-MPK6 cascade in Arabidopsis. Plant Sci.241, 138–150. 10.1016/j.plantsci.2015.10.005

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 65

  MartiniereA.BassilE.JublancE.AlconC.RegueraM.SentenacH.et al. (2013). In vivo intracellular pH measurements in tobacco and Arabidopsis reveal an unexpected pH gradient in the endomembrane system. Plant Cell25, 4028–4043. 10.1105/tpc.113.116897

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 66

  MatteiB.GallettiR.ManfrediniC.PontiggiaD.SalviG.SpadoniS.et al. (2005). Recognition and signalling in the cell wall: the case of endopolygalacturonase, PGIP and oligogalacturonides. Plant Biosyst.139, 24–27. 10.1080/112635000500059793

  • CrossRef
  • Google Scholar

• 67

  MhamdiA.MauveC.GouiaH.SaindrenanP.HodgesM.NoctorG. (2010). Cytosolic NADP-dependent isocitrate dehydrogenase contributes to redox homeostasis and the regulation of pathogen responses in Arabidopsis leaves. Plant Cell Environ.33, 1112–1123. 10.1111/j.1365-3040.2010.02133.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 68

  MithoeS. C.LudwigC.PelM. J.CucinottaM.CasartelliA.MbengueM.et al. (2016). Attenuation of pattern recognition receptor signaling is mediated by a MAP kinase kinase kinase. EMBO Rep.17, 441–454. 10.15252/embr.201540806

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 69

  MollapourM.TsutsumiS.TrumanA. W.XuW. P.VaughanC. K.BeebeK.et al. (2011). Threonine 22 phosphorylation attenuates Hsp90 interaction with cochaperones and affects its chaperone activity. Mol. Cell41, 672–681. 10.1016/j.molcel.2011.02.011

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 70

  MullerS. (2015). Plant intracellular transport: tracing functions of the retrograde kinesin. Curr. Biol.25, R808–R810. 10.1016/j.cub.2015.08.011

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 71

  NuhseT. S. (2012). Cell wall integrity signaling and innate immunity in plants. Front. Plant Sci.3:280. 10.3389/fpls.2012.00280

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 72

  NuhseT. S.BottrillA. R.JonesA. M.PeckS. C. (2007). Quantitative phosphoproteomic analysis of plasma membrane proteins reveals regulatory mechanisms of plant innate immune responses. Plant J.51, 931–940. 10.1111/j.1365-313X.2007.03192.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 73

  ParkE.NebenfuhrA. (2013). Myosin XIK of Arabidopsis thaliana accumulates at the root hair tip and is required for fast root hair growth. PLoS ONE8:e76745. 10.1371/journal.pone.0076745

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 74

  PeiroA.Izquierdo-GarciaA. C.Sanchez-NavarroJ. A.PallasV.MuletJ. M.AparicioF. (2014). Patellins 3 and 6, two members of the plant patellin family, interact with the movement protein of Alfalfa mosaic virus and interfere with viral movement. Mol. Plant Pathol.15, 881–891. 10.1111/mpp.12146

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 75

  PetermanT. K.OholY. M.McReynoldsL. J.LunaE. J. (2004). Patellin1, a novel Sec14-like protein, localizes to the cell plate and binds phosphoinositides. Plant Physiol.136, 3080–3094. 10.1104/pp.104.045369

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 76

  PontiggiaD.CiarcianelliJ.SalviG.CervoneF.De LorenzoG.MatteiB. (2015). Sensitive detection and measurement of oligogalacturonides in Arabidopsis. Front. Plant Sci.6:258. 10.3389/fpls.2015.00258

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 77

  PradoK.BoursiacY.Tournaire-RouxC.MonneuseJ. M.PostaireO.Da InesO.et al. (2013). Regulation of Arabidopsis leaf hydraulics involves light-dependent phosphorylation of aquaporins in veins. Plant Cell25, 1029–1039. 10.1105/tpc.112.108456

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 78

  PrakS.HemS.BoudetJ.ViennoisG.SommererN.RossignolM.et al. (2008). Multiple phosphorylations in the C-terminal tail of plant plasma membrane aquaporins. Mol. Cell. Proteomics7, 1019–1030. 10.1074/mcp.M700566-MCP200

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 79

  PribatA.SormaniR.Rousseau-GueutinM.JulkowskaM. M.TesterinkC.JoubesJ.et al. (2012). A novel class of PTEN protein in Arabidopsis displays unusual phosphoinositide phosphatase activity and efficiently binds phosphatidic acid. Biochem. J.441, 161–171. 10.1042/BJ20110776

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 80

  QinT.LiuX.LiJ.SunJ.SongL.MaoT. (2014). Arabidopsis microtubule-destabilizing protein 25 functions in pollen tube growth by severing actin filaments. Plant Cell26, 325–339. 10.1105/tpc.113.119768

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 81

  RasmussenM. W.RouxM.PetersenM.MundyJ. (2012). MAP kinase cascades in Arabidopsis innate immunity. Front. Plant Sci.3:169. 10.3389/fpls.2012.00169

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 82

  RayapuramN.BonhommeL.BigeardJ.HaddadouK.PrzybylskiC.HirtH.et al. (2014). Identification of novel PAMP-triggered phosphorylation and dephosphorylation events in Arabidopsis thaliana by quantitative phosphoproteomic analysis. J. Proteome Res.13, 2137–2151. 10.1021/pr401268v

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 83

  ReichardtI.SlaneD.El KasmiF.KnollC.FuchsR.MayerU.et al. (2011). Mechanisms of functional specificity among plasma-membrane syntaxins in Arabidopsis. Traffic12, 1269–1280. 10.1111/j.1600-0854.2011.01222.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 84

  ReilandS.MesserliG.BaerenfallerK.GerritsB.EndlerA.GrossmannJ.et al. (2009). Large-scale Arabidopsis phosphoproteome profiling reveals novel chloroplast kinase substrates and phosphorylation networks. Plant Physiol.150, 889–903. 10.1104/pp.109.138677

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 85

  ReymondP.KunzB.Paul-PletzerK.GrimmR.EckerskornC.FarmerE. E. (1996). Cloning of a cDNA encoding a plasma membrane-associated, uronide binding phosphoprotein with physical properties similar to viral movement proteins. Plant Cell8, 2265–2276. 10.1105/tpc.8.12.2265

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 86

  RobbinsM. L.RoyA.WangP. H.GaffoorI.SekhonR. S.de O. BuanafinaM. M.et al. (2013). Comparative proteomics analysis by DIGE and iTRAQ provides insight into the regulation of phenylpropanoids in maize. J. Proteomics.93, 254–275. 10.1016/j.jprot.2013.06.018

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 87

  RobinsonD. G. (2014). Trafficking of vacuolar sorting receptors: new data and new problems. Plant Physiol.165, 1417–1423. 10.1104/pp.114.243303

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 88

  RudashevskayaE. L.YeJ. Y.JensenO. N.FuglsangA. T.PalmgrenM. G. (2012). Phosphosite mapping of P-type plasma membrane H+-ATPase in homologous and heterologous environments. J. Biol. Chem.287, 4904–4913. 10.1074/jbc.M111.307264

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 89

  SavatinD. V.FerrariS.SiciliaF.De LorenzoG. (2011). Oligogalacturonide-auxin antagonism does not require posttranscriptional gene silencing or stabilization of auxin response repressors in Arabidopsis. Plant Physiol.157, 1163–1174. 10.1104/pp.111.184663

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 90

  SavatinD. V.GigliB. N.MartiL.FabbriC.CervoneF.De LorenzoG. (2014a). The Arabidopsis NPK1-related protein kinases ANPs are required for elicitor-induced oxidative burst and immunity. Plant Physiol.165, 1188–1202. 10.1104/pp.114.236901

  • CrossRef
  • Google Scholar

• 91

  SavatinD. V.GramegnaG.ModestiV.CervoneF. (2014b). Wounding in the plant tissue: the defense of a dangerous passage. Front. Plant Sci.5:470. 10.3389/fpls.2014.00470

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 92

  SchumacherK.KrebsM. (2010). The V-ATPase: small cargo, large effects. Curr. Opin. Plant Biol.13, 724–730. 10.1016/j.pbi.2010.07.003

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 93

  SmertenkoA. P.ChangH. Y.SonobeS.FenykS. I.WeingartnerM.BogreL.et al. (2006). Control of the AtMAP65-1 interaction with microtubules through the cell cycle. J. Cell Sci.119, 3227–3237. 10.1242/jcs.03051

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 94

  SpethC.JaspertN.MarconC.OeckingC. (2010). Regulation of the plant plasma membrane H+-ATPase by its C-terminal domain: what do we know for sure?Eur. J. Cell Biol.89, 145–151. 10.1016/j.ejcb.2009.10.015

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 95

  StasykT.MorandellS.BakryR.FeuersteinI.HuckC. W.StecherG.et al. (2005). Quantitative detection of phosphoproteins by combination of two-dimensional difference gel electrophoresis and phosphospecific fluorescent staining. Electrophoresis26, 2850–2854. 10.1002/elps.200500026

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 96

  SteckerK. E.MinkoffB. B.SussmanM. R. (2014). Phosphoproteomic analyses reveal early signaling events in the osmotic stress response. Plant Physiol.165, 1171–1187. 10.1104/pp.114.238816

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 97

  SteinhorstL.KudlaJ. (2013). Calcium and reactive oxygen species rule the waves of signaling. Plant Physiol.163, 471–485. 10.1104/pp.113.222950

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 98

  SuetsuguN.YamadaN.KagawaT.YonekuraH.UyedaT. Q. P.KadotaA.et al. (2010). Two kinesin-like proteins mediate actin-based chloroplast movement in Arabidopsis thaliana. Proc. Natl. Acad. Sci. U.S.A.107, 8860–8865. 10.1073/pnas.0912773107

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 99

  SzklarczykD.FranceschiniA.WyderS.ForslundK.HellerD.Huerta-CepasJ.et al. (2015). STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res.43, D447–D452. 10.1093/nar/gku1003

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 100

  TanakaH.OsakabeY.KatsuraS.MizunoS.MaruyamaK.KusakabeK.et al. (2012). Abiotic stress-inducible receptor-like kinases negatively control ABA signaling in Arabidopsis. Plant J.70, 599–613. 10.1111/j.1365-313X.2012.04901.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 101

  TangR. J.LiuH.YangY.YangL.GaoX. S.GarciaV. J.et al. (2012). Tonoplast calcium sensors CBL2 and CBL3 control plant growth and ion homeostasis through regulating V-ATPase activity in Arabidopsis. Cell Res.22, 1650–1665. 10.1038/cr.2012.161

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 102

  TangW. Q.DengZ. P.Oses-PrietoJ. A.SuzukiN.ZhuS. W.ZhangX.et al. (2008). Proteomics studies of brassinosteroid signal transduction using prefractionation and two-dimensional DIGE. Mol. Cell. Proteomics7, 728–738. 10.1074/mcp.M700358-MCP200

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 103

  TapkenW.MurphyA. S. (2015). Membrane nanodomains in plants: capturing form, function, and movement. J. Exp. Bot.66, 1573–1586. 10.1093/jxb/erv054

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 104

  ThingholmT. E.JensenO. N.RobinsonP. J.LarsenM. R. (2008). SIMAC (sequential elution from IMAC), a phosphoproteomics strategy for the rapid separation of monophosphorylated from multiply phosphorylated peptides. Mol. Cell. Proteomics7, 661–671. 10.1074/mcp.M700362-MCP200

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 105

  UlmasovT.HagenG.GuilfoyleT. J. (1999). Activation and repression of transcription by auxin-response factors. Proc. Natl. Acad. Sci. U.S.A.96, 5844–5849. 10.1073/pnas.96.10.5844

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 106

  UnderwoodW.SomervilleS. C. (2013). Perception of conserved pathogen elicitors at the plasma membrane leads to relocalization of the Arabidopsis PEN3 transporter. Proc. Natl. Acad. Sci. U.S.A.110, 12492–12497. 10.1073/pnas.1218701110

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 107

  VertG.WalcherC. L.ChoryJ.NemhauserJ. L. (2008). Integration of auxin and brassinosteroid pathways by Auxin Response Factor 2. Proc. Natl. Acad. Sci. U.S.A.105, 9829–9834. 10.1073/pnas.0803996105

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 108

  VizcaínoJ. A.CsordasA.del-ToroN.DianesJ. A.GrissJ.LavidasI.et al. (2016). 2016 update of the PRIDE database and its related tools. Nucleic. Acids Res. 44, D447–D456. 10.1093/nar/gkv1145

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 109

  VollL. M.ZellM. B.EngelsdorfT.SaurA.WheelerM. G.DrincovichM. F.et al. (2012). Loss of cytosolic NADP-malic enzyme 2 in Arabidopsis thaliana is associated with enhanced susceptibility to Colletotrichum higginsianum. New Phytol.195, 189–202. 10.1111/j.1469-8137.2012.04129.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 110

  WangX.BianY. Y.ChengK.GuL. F.YeM. L.ZouH. F.et al. (2013). A large-scale protein phosphorylation analysis reveals novel phosphorylation motifs and phosphoregulatory networks in Arabidopsis. J. Proteomics78, 486–498. 10.1016/j.jprot.2012.10.018

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 111

  XuJ.MengJ.MengX.ZhaoY.LiuJ.SunT.et al. (2016). Pathogen-responsive MPK3 and MPK6 reprogram the biosynthesis of indole glucosinolates and their derivatives in Arabidopsis immunity. Plant Cell28, 1144–1162. 10.1105/tpc.15.00871

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 112

  YamadaK.Hara-NishimuraI.NishimuraM. (2011). Unique defense strategy by the endoplasmic reticulum body in plants. Plant Cell Physiol.52, 2039–2049. 10.1093/pcp/pcr156

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 113

  YangL.QinL.LiuG. S.PeremyslovV. V.DoljaV. V.WeiY. D. (2014). Myosins XI modulate host cellular responses and penetration resistance to fungal pathogens. Proc. Natl. Acad. Sci. U.S.A.111, 13996–14001. 10.1073/pnas.1405292111

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 114

  YaoQ. M.GeH. Y.WuS. Q.ZhangN.ChenW.XuC. H.et al. (2014). (PDB)-D-3 3.0: from plant phosphorylation sites to protein networks. Nucleic Acids Res.42, D1206–D1213. 10.1093/nar/gkt1135

  • CrossRef
  • Google Scholar

• 115

  YehY. H.ChangY. H.HuangP. Y.HuangJ. B.ZimmerliL. (2015). Enhanced Arabidopsis pattern-triggered immunity by overexpression of cysteine-rich receptor-like kinases. Front. Plant Sci.6:322. 10.3389/fpls.2015.00322

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 116

  ZhaoZ.ZhangW.StanleyB. A.AssmannS. M. (2008). Functional proteomics of Arabidopsis thaliana guard cells uncovers new stomatal signaling pathways. Plant Cell20, 3210–3226. 10.1105/tpc.108.063263

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 117

  ZhuC.DixitR. (2012). Functions of the Arabidopsis kinesin superfamily of microtubule-based motor proteins. Protoplasma249, 887–899. 10.1007/s00709-011-0343-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar