The change of gravity vector induces short-term phosphoproteomic alterations in Arabidopsis

The Quantitative Biotinylproteomics Studies Reveal a WInd-Related Kinase 1 (Raf-Like Kinase 36) Functioning as an Early Signaling Component in Wind-Induced Thigmomorphogenesis and Gravitropism

Wind is one of the most prevalent environmental forces entraining plants to develop various mechano-responses, collectively called thigmomorphogenesis. Largely unknown is how plants transduce these versatile wind force signals downstream to nuclear events and to the development of thigmomorphogenic phenotype or anemotropic response. To identify molecular components at the early steps of the wind force signaling, two mechanical signaling-related phosphoproteins, identified from our previous phosphoproteomic study of Arabidopsis touch response, mitogen-activated protein kinase kinase 1 (MKK1) and 2 (MKK2), were selected for performing in planta TurboID (ID)-based quantitative proximity-labeling (PL) proteomics. This quantitative biotinylproteomics was separately performed on MKK1-ID and MKK2-ID transgenic plants, respectively, using the genetically engineered TurboID biotin ligase expression transgenics as a universal control. This unique PTM proteomics successfully identified 11 and 71 MKK1 and MKK2 putative interactors, respectively. Biotin occupancy ratio (BOR) was found to be an alternative parameter to measure the extent of proximity and specificity between the proximal target proteins and the bait fusion protein. Bioinformatics analysis of these biotinylprotein data also found that TurboID biotin ligase favorably labels the loop region of target proteins. A WInd-Related Kinase 1 (WIRK1), previously known as rapidly accelerated fibrosarcoma (Raf)-like kinase 36 (RAF36), was found to be a putative common interactor for both MKK1 and MKK2 and preferentially interacts with MKK2. Further molecular biology studies of the Arabidopsis RAF36 kinase found that it plays a role in wind regulation of the touch-responsive TCH3 and CML38 gene expression and the phosphorylation of a touch-regulated PATL3 phosphoprotein. Measurement of leaf morphology and shoot gravitropic response of wirk1 (raf36) mutant revealed that the WIRK1 gene is involved in both wind-triggered rosette thigmomorphogenesis and gravitropism of Arabidopsis stems, suggesting that the WIRK1 (RAF36) protein probably functioning upstream of both MKK1 and MKK2 and that it may serve as the crosstalk point among multiple mechano-signal transduction pathways mediating both wind mechano-response and gravitropism.

SILIA-Based 4C Quantitative PTM Proteomics

To absolutely and relatively quantitate the alteration of a posttranslationally modified (PTM) proteome in response to a specific internal or external signal, a ¹⁵N-stable isotope labeling in Arabidopsis (SILIA) protocol has been integrated into the 4C quantitative PTM proteomics, named as SILIA-based 4C quantitative PTM proteomics (S4Quap). The isotope metabolic labeling produces both forward (F) and reciprocal (R) mixings of either ¹⁴N/¹⁵N-coded tissues or the ¹⁴N/¹⁵N-coded total cellular proteins. Plant protein is isolated using a urea-based extraction buffer (UEB). The presence of 8 M urea, 2% polyvinylpolypyrrolidone (PVPP), and 5 mM ascorbic acid allows to instantly denature protein, remove the phenolic compounds, and curb the oxidation by free radicals once plant cells are broken. The total cellular proteins are routinely processed into peptides by trypsin. The PTM peptide yield of affinity enrichment and preparation is 0.1-0.2% in general. Ion exchange chromatographic fractionation prepares the PTM peptides for LC-MS/MS analysis. The collected mass spectrograms are subjected to a target-decoy sequence analysis using various search engines. The computational programs are subsequently applied to analyze the ratios of the extracted ion chromatogram (XIC) of the ¹⁴N/¹⁵N isotope-coded PTM peptide ions and to perform the statistical evaluation of the quantitation results. The Student t-test values of ratios of quantifiable ¹⁴N/¹⁵N-coded PTM peptides are normally corrected using a Benjamini-Hochberg (BH) multiple hypothesis test to select the significantly regulated PTM peptide groups (BH-FDR < 5%). Consequently, the highly selected prospect candidate(s) of PTM proteins are confirmed and validated using biochemical, molecular, cellular, and transgenic plant analysis.

qPTMplants: an integrative database of quantitative post-translational modifications in plants

As a crucial molecular mechanism, post-translational modifications (PTMs) play critical roles in a wide range of biological processes in plants. Recent advances in mass spectrometry-based proteomic technologies have greatly accelerated the profiling and quantification of plant PTM events. Although several databases have been constructed to store plant PTM data, a resource including more plant species and more PTM types with quantitative dynamics still remains to be developed. In this paper, we present an integrative database of quantitative PTMs in plants named qPTMplants (http://qptmplants.omicsbio.info), which hosts 1 242 365 experimentally identified PTM events for 429 821 nonredundant sites on 123 551 proteins under 583 conditions for 23 PTM types in 43 plant species from 293 published studies, with 620 509 quantification events for 136 700 PTM sites on 55 361 proteins under 354 conditions. Moreover, the experimental details, such as conditions, samples, instruments and methods, were manually curated, while a variety of annotations, including the sequence and structural characteristics, were integrated into qPTMplants. Then, various search and browse functions were implemented to access the qPTMplants data in a user-friendly manner. Overall, we anticipate that the qPTMplants database will be a valuable resource for further research on PTMs in plants.

Isotopically Dimethyl Labeling-Based Quantitative Proteomic Analysis of Phosphoproteomes of Soybean Cultivars

Isotopically dimethyl labeling was applied in a quantitative post-translational modification (PTM) proteomic study of phosphoproteomic changes in the drought responses of two contrasting soybean cultivars. A total of 9457 phosphopeptides were identified subsequently, corresponding to 4571 phosphoprotein groups and 3889 leading phosphoproteins, which contained nine kinase families consisting of 279 kinases. These phosphoproteins contained a total of 8087 phosphosites, 6106 of which were newly identified and constituted 54% of the current soybean phosphosite repository. These phosphosites were converted into the highly conserved kinase docking sites by bioinformatics analysis, which predicted six kinase families that matched with those newly found nine kinase families. The overly post-translationally modified proteins (OPP) occupies 2.1% of these leading phosphoproteins. Most of these OPPs are photoreceptors, mRNA-, histone-, and phospholipid-binding proteins, as well as protein kinase/phosphatases. The subgroup population distribution of phosphoproteins over the number of phosphosites of phosphoproteins follows the exponential decay law, Y = 4.13e-0.098X - 0.04. Out of 218 significantly regulated unique phosphopeptide groups, 188 phosphoproteins were regulated by the drought-tolerant cultivar under the water loss condition. These significantly regulated phosphoproteins (SRP) are mainly enriched in the biological functions of water transport and deprivation, methionine metabolic processes, photosynthesis/light reaction, and response to cadmium ion, osmotic stress, and ABA response. Seventeen and 15 SRPs are protein kinases/phosphatases and transcription factors, respectively. Bioinformatics analysis again revealed that three members of the calcium dependent protein kinase family (CAMK family), GmSRK2I, GmCIPK25, and GmAKINβ1 kinases, constitute a phosphor-relay-mediated signal transduction network, regulating ion channel activities and many nuclear events in this drought-tolerant cultivar, which presumably contributes to the development of the soybean drought tolerance under water deprivation process.