P(3)DB: An Integrated Database for Plant Protein Phosphorylation

Deep Learning in Phosphoproteomics: Methods and Application in Cancer Drug Discovery

Protein phosphorylation is a key post-translational modification (PTM) that is a central regulatory mechanism of many cellular signaling pathways. Several protein kinases and phosphatases precisely control this biochemical process. Defects in the functions of these proteins have been implicated in many diseases, including cancer. Mass spectrometry (MS)-based analysis of biological samples provides in-depth coverage of phosphoproteome. A large amount of MS data available in public repositories has unveiled big data in the field of phosphoproteomics. To address the challenges associated with handling large data and expanding confidence in phosphorylation site prediction, the development of many computational algorithms and machine learning-based approaches have gained momentum in recent years. Together, the emergence of experimental methods with high resolution and sensitivity and data mining algorithms has provided robust analytical platforms for quantitative proteomics. In this review, we compile a comprehensive collection of bioinformatic resources used for the prediction of phosphorylation sites, and their potential therapeutic applications in the context of cancer.

Regulation of PaRBOH1-mediated ROS production in Norway spruce by Ca2+ binding and phosphorylation

Plant respiratory burst oxidase homologs (RBOHs) are plasma membrane-localized NADPH oxidases that generate superoxide anion radicals, which then dismutate to H₂O₂, into the apoplast using cytoplasmic NADPH as an electron donor. PaRBOH1 is the most highly expressed RBOH gene in developing xylem as well as in a lignin-forming cell culture of Norway spruce (Picea abies L. Karst.). Since no previous information about regulation of gymnosperm RBOHs exist, our aim was to resolve how PaRBOH1 is regulated with a focus on phosphorylation. The N-terminal part of PaRBOH1 was found to contain several putative phosphorylation sites and a four-times repeated motif with similarities to the Botrytis-induced kinase 1 target site in Arabidopsis AtRBOHD. Phosphorylation was indicated for six of the sites in in vitro kinase assays using 15 amino-acid-long peptides for each of the predicted phosphotarget site in the presence of protein extracts of developing xylem. Serine and threonine residues showing positive response in the peptide assays were individually mutated to alanine (kinase-inactive) or to aspartate (phosphomimic), and the wild type PaRBOH1 and the mutated constructs transfected to human kidney embryogenic (HEK293T) cells with a low endogenous level of extracellular ROS production. ROS-producing assays with HEK cells showed that Ca²⁺ and phosphorylation synergistically activate the enzyme and identified several serine and threonine residues that are likely to be phosphorylated including a novel phosphorylation site not characterized in other plant species. These were further investigated with a phosphoproteomic study. Results of Norway spruce, the first gymnosperm species studied in relation to RBOH regulation, show that regulation of RBOH activity is conserved among seed plants.

Dynamics of Protein Phosphorylation during Arabidopsis Seed Germination

Seed germination is critical for early plantlet development and is tightly controlled by environmental factors. Nevertheless, the signaling networks underlying germination control remain elusive. In this study, the remodeling of Arabidopsis seed phosphoproteome during imbibition was investigated using stable isotope dimethyl labeling and nanoLC-MS/MS analysis. Freshly harvested seeds were imbibed under dark or constant light to restrict or promote germination, respectively. For each light regime, phosphoproteins were extracted and identified from dry and imbibed (6 h, 16 h, and 24 h) seeds. A large repertoire of 10,244 phosphopeptides from 2546 phosphoproteins, including 110 protein kinases and key regulators of seed germination such as Delay Of Germination 1 (DOG1), was established. Most phosphoproteins were only identified in dry seeds. Early imbibition led to a similar massive downregulation in dormant and non-dormant seeds. After 24 h, 411 phosphoproteins were specifically identified in non-dormant seeds. Gene ontology analyses revealed their involvement in RNA and protein metabolism, transport, and signaling. In addition, 489 phosphopeptides were quantified, and 234 exhibited up or downregulation during imbibition. Interaction networks and motif analyses revealed their association with potential signaling modules involved in germination control. Our study provides evidence of a major role of phosphosignaling in the regulation of Arabidopsis seed germination.

Arabidopsis bioinformatics: tools and strategies

The sequencing of the Arabidopsis thaliana genome 21 years ago ushered in the genomics era for plant research. Since then, an incredible variety of bioinformatic tools permit easy access to large repositories of genomic, transcriptomic, proteomic, epigenomic and other '-omic' data. In this review, we cover some more recent tools (and highlight the 'classics') for exploring such data in order to help formulate quality, testable hypotheses, often without having to generate new experimental data. We cover tools for examining gene expression and co-expression patterns, undertaking promoter analyses and gene set enrichment analyses, and exploring protein-protein and protein-DNA interactions. We will touch on tools that integrate different data sets at the end of the article.

Novel Translational and Phosphorylation Modification Regulation Mechanisms of Tomato (Solanum lycopersicum) Fruit Ripening Revealed by Integrative Proteomics and Phosphoproteomics

The tomato is a research model for fruit-ripening, however, its fruit-ripening mechanism still needs more extensive and in-depth exploration. Here, using TMT and LC-MS, the proteome and phosphoproteome of AC⁺⁺ (wild type) and rin (ripening-inhibitor) mutant fruits were studied to investigate the translation and post-translational regulation mechanisms of tomato fruit-ripening. A total of 6141 proteins and 4011 phosphorylation sites contained quantitative information. One-hundred proteins were identified in both omics’ profiles, which were mainly found in ethylene biosynthesis and signal transduction, photosynthesis regulation, carotenoid and flavonoid biosynthesis, chlorophyll degradation, ribosomal subunit expression changes, MAPK pathway, transcription factors and kinases. The affected protein levels were correlated with their corresponding gene transcript levels, such as NAC-NOR, MADS-RIN, IMA, TAGL1, MADS-MC and TDR4. Changes in the phosphorylation levels of NAC-NOR and IMA were involved in the regulation of tomato fruit-ripening. Although photosynthesis was inhibited, there were diverse primary and secondary metabolic pathways, such as glycolysis, fatty acid metabolism, vitamin metabolism and isoprenoid biosynthesis, regulated by phosphorylation. These data constitute a map of protein—protein phosphorylation in the regulation of tomato fruit-ripening, which lays the foundation for future in-depth study of the sophisticated molecular mechanisms of fruit-ripening and provide guidance for molecular breeding.

Bioinformatic Tools in Arabidopsis Research

Bioinformatic tools are now an everyday part of a plant researcher's collection of protocols. They allow almost instantaneous access to large data sets encompassing genomes, transcriptomes, proteomes, epigenomes, and other "-omes," which are now being generated with increasing speed and decreasing cost. With the appropriate queries, such tools can generate quality hypotheses, sometimes without the need for new experimental data. In this chapter, we will investigate some of the tools used for examining gene expression and coexpression patterns, performing promoter analyses and functional classification enrichment for sets of genes, and exploring protein-protein and protein-DNA interactions in Arabidopsis. We will also cover additional tools that allow integration of data from several sources for improved hypothesis generation.

Global Quantitative Proteomics Studies Revealed Tissue-Preferential Expression and Phosphorylation of Regulatory Proteins in Arabidopsis

Organogenesis in plants occurs across all stages of the life cycle. Although previous studies have identified many genes as important for either vegetative or reproductive development at the RNA level, global information on translational and post-translational levels remains limited. In this study, six Arabidopsis stages/organs were analyzed using quantitative proteomics and phosphoproteomics, identifying 2187 non-redundant proteins and evidence for 1194 phosphoproteins. Compared to the expression observed in cauline leaves, the expression of 1445, 1644, and 1377 proteins showed greater than 1.5-fold alterations in stage 1–9 flowers, stage 10–12 flowers, and open flowers, respectively. Among these, 294 phosphoproteins with 472 phosphorylation sites were newly uncovered, including 275 phosphoproteins showing differential expression patterns, providing molecular markers and possible candidates for functional studies. Proteins encoded by genes preferentially expressed in anther (15), meiocyte (4), or pollen (15) were enriched in reproductive organs, and mutants of two anther-preferentially expressed proteins, acos5 and mee48, showed obviously reduced male fertility with abnormally organized pollen exine. In addition, more phosphorylated proteins were identified in reproductive stages (1149) than in the vegetative organs (995). The floral organ-preferential phosphorylation of GRP17, CDC2/CDKA.1, and ATSK11 was confirmed with western blot analysis. Moreover, phosphorylation levels of CDPK6 and MAPK6 and their interacting proteins were elevated in reproductive tissues. Overall, our study yielded extensive data on protein expression and phosphorylation at six stages/organs and provides an important resource for future studies investigating the regulatory mechanisms governing plant development.

Comparative phosphoproteomic analysis of BR-defective mutant reveals a key role of GhSK13 in regulating cotton fiber development

Brassinosteroid (BR), a steroid phytohormone, whose signaling transduction pathways include a series of phosphorylation and dephosphorylation events, and GSK3s are the main negative regulator kinases. BRs have been shown to play vital roles in cotton fiber elongation. However, the underlying mechanism is still elusive. In this study, fibers of a BR-defective mutant Pagoda 1 (pag1), and its corresponding wild-type (ZM24) were selected for a comparative global phosphoproteome analysis at critical developmental time points: fast-growing stage (10 days after pollination (DPA)) and secondary cell wall synthesis stage (20 DPA). Based on the substrate characteristics of GSK3, 900 potential substrates were identified. Their GO and KEGG annotation results suggest that BR functions in fiber development by regulating GhSKs (GSK3s of Gossypium hirsutum L.) involved microtubule cytoskeleton organization, and pathways of glucose, sucrose and lipid metabolism. Further experimental results revealed that among the GhSK members identified, GhSK13 not only plays a role in BR signaling pathway, but also functions in developing fiber by respectively interacting with an AP2-like ethylene-responsive factor GhAP2L, a nuclear transcription factor Gh_DNF_YB19, and a homeodomain zipper member GhHDZ5. Overall, our phosphoproteomic research advances the understanding of fiber development controlled by BR signal pathways especially through GhSKs, and also offers numbers of target proteins for improving cotton fiber quality.

Bioinformatics Resources for Plant Abiotic Stress Responses: State of the Art and Opportunities in the Fast Evolving -Omics Era

Abiotic stresses are among the principal limiting factors for productivity in agriculture. In the current era of continuous climate changes, the understanding of the molecular aspects involved in abiotic stress response in plants is a priority. The rise of -omics approaches provides key strategies to promote effective research in the field, facilitating the investigations from reference models to an increasing number of species, tolerant and sensitive genotypes. Integrated multilevel approaches, based on molecular investigations at genomics, transcriptomics, proteomics and metabolomics levels, are now feasible, expanding the opportunities to clarify key molecular aspects involved in responses to abiotic stresses. To this aim, bioinformatics has become fundamental for data production, mining and integration, and necessary for extracting valuable information and for comparative efforts, paving the way to the modeling of the involved processes. We provide here an overview of bioinformatics resources for research on plant abiotic stresses, describing collections from -omics efforts in the field, ranging from raw data to complete databases or platforms, highlighting opportunities and still open challenges in abiotic stress research based on -omics technologies.

Characterization of LuWRKY36, a flax transcription factor promoting secoisolariciresinol biosynthesis in response to Fusarium oxysporum elicitors in Linum usitatissimum L. hairy roots

The involvement of a WRKY transcription factor in the regulation of lignan biosynthesis in flax using a hairy root system is described. Secoisolariciresinol is the main flax lignan synthesized by action of LuPLR1 (pinoresinol-lariciresinol reductase 1). LuPLR1 gene promoter deletion experiments have revealed a promoter region containing W boxes potentially responsible for the response to Fusarium oxysporum. W boxes are bound by WRKY transcription factors that play a role in the response to stress. A candidate WRKY transcription factor, LuWRKY36, was isolated from both abscisic acid and Fusarium elicitor-treated flax cell cDNA libraries. This transcription factors contains two WRKY DNA-binding domains and is a homolog of AtWRKY33. Different approaches confirmed LuWRKY36 binding to a W box located in the LuPLR1 promoter occurring through a unique direct interaction mediated by its N-terminal WRKY domain. Our results propose that the positive regulator action of LuWRKY36 on the LuPLR1 gene regulation and lignan biosynthesis in response to biotic stress is positively mediated by abscisic acid and inhibited by ethylene. Additionally, we demonstrate a differential Fusarium elicitor response in susceptible and resistant flax cultivars, seen as a faster and stronger LuPLR1 gene expression response accompanied with higher secoisolariciresinol accumulation in HR of the resistant cultivar.

Genome-wide identification, structural and gene expression analysis of the bZIP transcription factor family in sweet potato wild relative Ipomoea trifida

Background

The basic leucine zipper (bZIP) transcription factor is one of the most abundant and conserved transcription factor families. In addition to being involved in growth and development, bZIP transcription factors also play an important role in plant adaption to abiotic stresses.

Results

A total of 41 bZIP genes that encode 66 proteins were identified in Ipomoea trifida. They were distributed on 14 chromosomes of Ipomoea trifida. Segmental and tandem duplication analysis showed that segmental duplication played an important role in the ItfbZIP gene amplification. ItfbZIPs were divided into ten groups (A, B, C, D, E, F, G, H, I and S groups) according to their phylogenetic relationships with Solanum lycopersicum and Arabidopsis thaliana. The regularity of the exon/intron numbers and distributions is consistent with the group classification in evolutionary tree. Prediction of the cis-acting elements found that promoter regions of ItfbZIPs harbored several stress responsive cis-acting elements. Protein three-dimensional structural analysis indicated that ItfbZIP proteins mainly consisted of α-helices and random coils. The gene expression pattern from transcriptome data and qRT-PCR analysis showed that ItfbZIP genes expressed with a tissue-specific manner and differently expressed under various abiotic stresses, suggesting that the ItfbZIPs were involved in stress response and adaption in Ipomoea trifida.

Conclusions

Genome-wide identification, gene structure, phylogeny and expression analysis of bZIP gene in Ipomoea trifida supplied a solid theoretical foundation for the functional study of bZIP gene family and further facilitated the molecular breeding of sweet potato.

Electronic supplementary material

The online version of this article (10.1186/s12863-019-0743-y) contains supplementary material, which is available to authorized users.

Global Phosphoproteomic Analysis Reveals the Defense and Response Mechanisms of Jatropha Curcas Seedling under Chilling Stress

As a promising energy plant for biodiesel, Jatropha curcas is a tropical and subtropical shrub and its growth is affected by one of major abiotic stress, chilling. Therefore, we adopt the phosphoproteomic analysis, physiological measurement and ultrastructure observation to illustrate the responsive mechanism of J. curcas seedling under chilling (4 °C) stress. After chilling for 6 h, 308 significantly changed phosphoproteins were detected. Prolonged the chilling treatment for 24 h, obvious physiological injury can be observed and a total of 332 phosphoproteins were examined to be significantly changed. After recovery (28 °C) for 24 h, 291 phosphoproteins were varied at the phosphorylation level. GO analysis showed that significantly changed phosphoproteins were mainly responsible for cellular protein modification process, transport, cellular component organization and signal transduction at the chilling and recovery periods. On the basis of protein-protein interaction network analysis, phosphorylation of several protein kinases, such as SnRK2, MEKK1, EDR1, CDPK, EIN2, EIN4, PI4K and 14-3-3 were possibly responsible for cross-talk between ABA, Ca²⁺, ethylene and phosphoinositide mediated signaling pathways. We also highlighted the phosphorylation of HOS1, APX and PIP2 might be associated with response to chilling stress in J. curcas seedling. These results will be valuable for further study from the molecular breeding perspective.

G protein subunit phosphorylation as a regulatory mechanism in heterotrimeric G protein signaling in mammals, yeast, and plants

Heterotrimeric G proteins composed of Gα, Gβ, and Gγ subunits are vital eukaryotic signaling elements that convey information from ligand-regulated G protein-coupled receptors (GPCRs) to cellular effectors. Heterotrimeric G protein-based signaling pathways are fundamental to human health [Biochimica et Biophysica Acta (2007) 1768, 994-1005] and are the target of >30% of pharmaceuticals in clinical use [Biotechnology Advances (2013) 31, 1676-1694; Nature Reviews Drug Discovery (2017) 16, 829-842]. This review focuses on phosphorylation of G protein subunits as a regulatory mechanism in mammals, budding yeast, and plants. This is a re-emerging field, as evidence for phosphoregulation of mammalian G protein subunits from biochemical studies in the early 1990s can now be complemented with contemporary phosphoproteomics and genetic approaches applied to a diversity of model systems. In addition, new evidence implicates a family of plant kinases, the receptor-like kinases, which are monophyletic with the interleukin-1 receptor-associated kinase/Pelle kinases of metazoans, as possible GPCRs that signal via subunit phosphorylation. We describe early and modern observations on G protein subunit phosphorylation and its functional consequences in these three classes of organisms, and suggest future research directions.

Evidence that the Arabidopsis thaliana 3-hydroxy-3-methylglutaryl-CoA reductase 1 is phosphorylated at Ser577 in planta

3-Hydroxy-3-methylglutaryl-CoA reductase (HMGR) is an essential enzyme in the mevalonate pathway. In higher plants, mevalonate pathway involves in the production of precursor for isoprenoids biosynthesis, including essential components for cell functions. Previously, we confirmed that the Arabidopsis thaliana HMGR1S (AtHMGR1S) is phosphorylated at S577 by the combination of sucrose non-fermenting related kinase-1 (SnRK1) and geminivirus rep-interacting kinase-1 (GRIK1) in vitro. However, even in quantitative phosphoproteomics studies that were directed to find SnRK1 target substrates, AtHMGR1S phosphorylation at S577 has never been detected in planta. In this study, we expressed AtHMGR1S as a C-terminal FLAG-fusion protein in A. thaliana hmg1 mutant to confirm its phosphorylation in planta. Our results provide the first direct evidence that AtHMGR1S is phosphorylated at S577 in planta. Moreover, phosphatase inhibitors treatment to the A. thaliana seedlings induced AtHMGR1S phosphorylation at sites other than S577, suggesting the presence of a novel HMGR regulatory mechanism in planta.

Structural disorder in plant proteins: where plasticity meets sessility

Plants are sessile organisms. This intriguing nature provokes the question of how they survive despite the continual perturbations caused by their constantly changing environment. The large amount of knowledge accumulated to date demonstrates the fascinating dynamic and plastic mechanisms, which underpin the diverse strategies selected in plants in response to the fluctuating environment. This phenotypic plasticity requires an efficient integration of external cues to their growth and developmental programs that can only be achieved through the dynamic and interactive coordination of various signaling networks. Given the versatility of intrinsic structural disorder within proteins, this feature appears as one of the leading characters of such complex functional circuits, critical for plant adaptation and survival in their wild habitats. In this review, we present information of those intrinsically disordered proteins (IDPs) from plants for which their high level of predicted structural disorder has been correlated with a particular function, or where there is experimental evidence linking this structural feature with its protein function. Using examples of plant IDPs involved in the control of cell cycle, metabolism, hormonal signaling and regulation of gene expression, development and responses to stress, we demonstrate the critical importance of IDPs throughout the life of the plant.

The Kinase Activity of Calcineurin B-like Interacting Protein Kinase 26 (CIPK26) Influences Its Own Stability and that of the ABA-regulated Ubiquitin Ligase, Keep on Going (KEG)

The Really Interesting New Gene (RING)-type E3 ligase, Keep on Going (KEG) plays a critical role in Arabidopsis growth after germination and the connections between KEG and hormone signaling pathways are expanding. With regards to abscisic acid (ABA) signaling, KEG targets ABA-responsive transcription factors abscisic acid insensitive 5, ABF1 and ABF3 for ubiquitination and subsequent degradation through the 26S proteasome. Regulation of E3 ligases through self-ubiquitination is common to RING-type E3 ligases and ABA promotes KEG self-ubiquitination and degradation. ABA-mediated degradation of KEG is phosphorylation-dependent; however, upstream signaling proteins that may regulate KEG stability have not been characterized. In this report, we show that CBL-Interacting Protein Kinase (CIPK) 26 can phosphorylate KEG in vitro. Using both in vitro and in planta degradation assays we provide evidence which suggests that the kinase activity of CIPK26 promotes the degradation of KEG. Furthermore, we found that the kinase activity of CIPK26 also influences its own stability; a constitutively active version is more stable than a wild type or a kinase dead version. Our results suggest a reciprocal regulation model wherein an activated and stable CIPK26 phosphorylates KEG to promote degradation of the E3.

Bioinformatics Analysis of Protein Phosphorylation in Plant Systems Biology Using P3DB

Protein phosphorylation is one of the most pervasive protein post-translational modification events in plant cells. It is involved in many plant biological processes, such as plant growth, organ development, and plant immunology, by regulating or switching signaling and metabolic pathways. High-throughput experimental methods like mass spectrometry can easily characterize hundreds to thousands of phosphorylation events in a single experiment. With the increasing volume of the data sets, Plant Protein Phosphorylation DataBase (P3DB, http://p3db.org ) provides a comprehensive, systematic, and interactive online platform to deposit, query, analyze, and visualize these phosphorylation events in many plant species. It stores the protein phosphorylation sites in the context of identified mass spectra, phosphopeptides, and phosphoproteins contributed from various plant proteome studies. In addition, P3DB associates these plant phosphorylation sites to protein physicochemical information in the protein charts and tertiary structures, while various protein annotations from hierarchical kinase phosphatase families, protein domains, and gene ontology are also added into the database. P3DB not only provides rich information, but also interconnects and provides visualization of the data in networks, in systems biology context. Currently, P3DB includes the KiC (Kinase Client) assay network, the protein-protein interaction network, the kinase-substrate network, the phosphatase-substrate network, and the protein domain co-occurrence network. All of these are available to query for and visualize existing phosphorylation events. Although P3DB only hosts experimentally identified phosphorylation data, it provides a plant phosphorylation prediction model for any unknown queries on the fly. P3DB is an entry point to the plant phosphorylation community to deposit and visualize any customized data sets within this systems biology framework. Nowadays, P3DB has become one of the major bioinformatics platforms of protein phosphorylation in plant biology.

Dynamic Phosphoproteome Analysis of Seedling Leaves in Brachypodium distachyon L. Reveals Central Phosphorylated Proteins Involved in the Drought Stress Response

Drought stress is a major abiotic stress affecting plant growth and development. In this study, we performed the first dynamic phosphoproteome analysis of Brachypodium distachyon L. seedling leaves under drought stress for different times. A total of 4924 phosphopeptides, contained 6362 phosphosites belonging to 2748 phosphoproteins. Rigorous standards were imposed to screen 484 phosphorylation sites, representing 442 unique phosphoproteins. Comparative analyses revealed significant changes in phosphorylation levels at 0, 6, and 24 h under drought stress. The most phosphorylated proteins and the highest phosphorylation level occurred at 6 h. Venn analysis showed that the up-regulated phosphopeptides at 6 h were almost two-fold those at 24 h. Motif-X analysis identified the six motifs: [sP], [Rxxs], [LxRxxs], [sxD], [sF], and [TP], among which [LxRxxs] was also previously identified in B. distachyon. Results from molecular function and protein-protein interaction analyses suggested that phosphoproteins mainly participate in signal transduction, gene expression, drought response and defense, photosynthesis and energy metabolism, and material transmembrane transport. These phosphoproteins, which showed significant changes in phosphorylation levels, play important roles in signal transduction and material transmembrane transport in response to drought conditions. Our results provide new insights into the molecular mechanism of this plant’s abiotic stress response through phosphorylation modification.

Quantitative phosphoproteomics reveals the role of the AMPK plant ortholog SnRK1 as a metabolic master regulator under energy deprivation

Since years, research on SnRK1, the major cellular energy sensor in plants, has tried to define its role in energy signalling. However, these attempts were notoriously hampered by the lethality of a complete knockout of SnRK1. Therefore, we generated an inducible amiRNA::SnRK1α2 in a snrk1α1 knock out background (snrk1α1/α2) to abolish SnRK1 activity to understand major systemic functions of SnRK1 signalling under energy deprivation triggered by extended night treatment. We analysed the in vivo phosphoproteome, proteome and metabolome and found that activation of SnRK1 is essential for repression of high energy demanding cell processes such as protein synthesis. The most abundant effect was the constitutively high phosphorylation of ribosomal protein S6 (RPS6) in the snrk1α1/α2 mutant. RPS6 is a major target of TOR signalling and its phosphorylation correlates with translation. Further evidence for an antagonistic SnRK1 and TOR crosstalk comparable to the animal system was demonstrated by the in vivo interaction of SnRK1α1 and RAPTOR1B in the cytosol and by phosphorylation of RAPTOR1B by SnRK1α1 in kinase assays. Moreover, changed levels of phosphorylation states of several chloroplastic proteins in the snrk1α1/α2 mutant indicated an unexpected link to regulation of photosynthesis, the main energy source in plants.

DAPPLE 2: a Tool for the Homology-Based Prediction of Post-Translational Modification Sites

The post-translational modification of proteins is critical for regulating their function. Although many post-translational modification sites have been experimentally determined, particularly in certain model organisms, experimental knowledge of these sites is severely lacking for many species. Thus, it is important to be able to predict sites of post-translational modification in such species. Previously, we described DAPPLE, a tool that facilitates the homology-based prediction of one particular post-translational modification, phosphorylation, in an organism of interest using known phosphorylation sites from other organisms. Here, we describe DAPPLE 2, which expands and improves upon DAPPLE in three major ways. First, it predicts sites for many post-translational modifications (20 different types) using data from several sources (15 online databases). Second, it has the ability to make predictions approximately 2-7 times faster than DAPPLE depending on the database size and the organism of interest. Third, it simplifies and accelerates the process of selecting predicted sites of interest by categorizing them based on gene ontology terms, keywords, and signaling pathways. We show that DAPPLE 2 can successfully predict known human post-translational modification sites using, as input, known sites from species that are either closely (e.g., mouse) or distantly (e.g., yeast) related to humans. DAPPLE 2 can be accessed at http://saphire.usask.ca/saphire/dapple2 .

An autophosphorylation site database for leucine-rich repeat receptor-like kinases in Arabidopsis thaliana

Leucine-rich repeat receptor-like kinases (LRR RLKs) form a large family of plant signaling proteins consisting of an extracellular domain connected by a single-pass transmembrane sequence to a cytoplasmic kinase domain. Autophosphorylation on specific Ser and/or Thr residues in the cytoplasmic domain is often critical for the activation of several LRR RLK family members with proven functional roles in plant growth regulation, morphogenesis, disease resistance, and stress responses. While identification and functional characterization of in vivo phosphorylation sites is ultimately required for a full understanding of LRR RLK biology and function, bacterial expression of recombinant LRR RLK cytoplasmic catalytic domains for identification of in vitro autophosphorylation sites provides a useful resource for further targeted identification and functional analysis of in vivo sites. In this study we employed high-throughput cloning and a variety of mass spectrometry approaches to generate an autophosphorylation site database representative of more than 30% of the approximately 223 LRR RLKs in Arabidopsis thaliana. We used His-tagged constructs of complete cytoplasmic domains to identify a total of 592 phosphorylation events across 73 LRR RLKs, with 497 sites uniquely assigned to specific Ser (268 sites) or Thr (229 sites) residues in 68 LRR RLKs. Multiple autophosphorylation sites per LRR RLK were the norm, with an average of seven sites per cytoplasmic domain, while some proteins showed more than 20 unique autophosphorylation sites. The database was used to analyze trends in the localization of phosphorylation sites across cytoplasmic kinase subdomains and to derive a statistically significant sequence motif for phospho-Ser autophosphorylation.

Case study: using sequence homology to identify putative phosphorylation sites in an evolutionarily distant species (honeybee)

The majority of scientific resources are devoted to studying a relatively small number of model species, meaning that the ability to translate knowledge across species is of considerable importance. Obtaining species-specific knowledge enables targeted investigations of the biology and pathobiology of a particular species, and facilitates comparative analyses. Phosphorylation is the most widespread posttranslational modification in eukaryotes, and although many phosphorylation sites have been experimentally identified for some species, little or no data are available for others. Using the honeybee as a test organism, this case study illustrates the process of using protein sequence homology to identify putative phosphorylation sites in a species of interest using experimentally determined sites from other species. A number of issues associated with this process are examined and discussed. Several databases of experimentally determined phosphorylation sites exist; however, it can be difficult for the nonspecialist to ascertain how their contents compare. Thus, this case study assesses the content and comparability of several phosphorylation site databases. Additional issues examined include the efficacy of homology-based phosphorylation site prediction, the impact of the level of evolutionary relatedness between species in making these predictions, the ability to translate knowledge of phosphorylation sites across large evolutionary distances and the criteria that should be used in selecting probable phosphorylation sites in the species of interest. Although focusing on phosphorylation, the issues discussed here also apply to the homology-based cross-species prediction of other posttranslational modifications, as well as to sequence motifs in general.

The pea seedling mitochondrial Nε-lysine acetylome

Posttranslational lysine acetylation is believed to occur in all taxa and to affect thousands of proteins. In contrast to the hundreds of mitochondrial proteins reported to be lysine-acetylated in non-plant species, only a handful have been reported from the plant taxa previously examined. To investigate whether this reflects a biologically significant difference or merely a peculiarity of the samples thus far examined, we immunoenriched and analyzed acetylated peptides from highly purified pea seedling mitochondria using mass spectrometry. Our results indicate that a multitude of mitochondrial proteins, involved in a variety of processes, are acetylated in pea seedlings.

The global phosphoproteome of Chlamydomonas reinhardtii reveals complex organellar phosphorylation in the flagella and thylakoid membrane

Chlamydomonas reinhardtii is the most intensively-studied and well-developed model for investigation of a wide-range of microalgal processes ranging from basic development through understanding triacylglycerol production. Although proteomic technologies permit interrogation of these processes at the protein level and efforts to date indicate phosphorylation-based regulation of proteins in C. reinhardtii is essential for its underlying biology, characterization of the C. reinhardtii phosphoproteome has been limited. Herein, we report the richest exploration of the C. reinhardtii proteome to date. Complementary enrichment strategies were used to detect 4588 phosphoproteins distributed among every cellular component in C. reinhardtii. Additionally, we report 18,160 unique phosphopeptides at <1% false discovery rate, which comprise 15,862 unique phosphosites - 98% of which are novel. Given that an estimated 30% of proteins in a eukaryotic cell are subject to phosphorylation, we report the majority of the phosphoproteome (23%) of C. reinhardtii. Proteins in key biological pathways were phosphorylated, including photosynthesis, pigment production, carbon assimilation, glycolysis, and protein and carbohydrate metabolism, and it is noteworthy that hyperphosphorylation was observed in flagellar proteins. This rich data set is available via ProteomeXchange (ID: PXD000783) and will significantly enhance understanding of a range of regulatory mechanisms controlling a variety of cellular process and will serve as a critical resource for the microalgal community.

Meta-Analysis of Arabidopsis thaliana Phospho-Proteomics Data Reveals Compartmentalization of Phosphorylation Motifs

Protein (de)phosphorylation plays an important role in plants. To provide a robust foundation for subcellular phosphorylation signaling network analysis and kinase-substrate relationships, we performed a meta-analysis of 27 published and unpublished in-house mass spectrometry-based phospho-proteome data sets for Arabidopsis thaliana covering a range of processes, (non)photosynthetic tissue types, and cell cultures. This resulted in an assembly of 60,366 phospho-peptides matching to 8141 nonredundant proteins. Filtering the data for quality and consistency generated a set of medium and a set of high confidence phospho-proteins and their assigned phospho-sites. The relation between single and multiphosphorylated peptides is discussed. The distribution of p-proteins across cellular functions and subcellular compartments was determined and showed overrepresentation of protein kinases. Extensive differences in frequency of pY were found between individual studies due to proteomics and mass spectrometry workflows. Interestingly, pY was underrepresented in peroxisomes but overrepresented in mitochondria. Using motif-finding algorithms motif-x and MMFPh at high stringency, we identified compartmentalization of phosphorylation motifs likely reflecting localized kinase activity. The filtering of the data assembly improved signal/noise ratio for such motifs. Identified motifs were linked to kinases through (bioinformatic) enrichment analysis. This study also provides insight into the challenges/pitfalls of using large-scale phospho-proteomic data sets to nonexperts.

Large-scale phosphoproteome analysis in seedling leaves of Brachypodium distachyon L

Background

Protein phosphorylation is one of the most important post-translational modifications involved in the regulation of plant growth and development as well as diverse stress response. As a member of the Poaceae, Brachypodium distachyon L. is a new model plant for wheat and barley as well as several potential biofuel grasses such as switchgrass. Vegetative growth is vital for biomass accumulation of plants, but knowledge regarding the role of protein phosphorylation modification during vegetative growth, especially in biofuel plants, is far from comprehensive.

Results

In this study, we carried out the first large-scale phosphoproteome analysis of seedling leaves in Brachypodium accession Bd21 using TiO₂ microcolumns combined with liquid chromatography-tandem mass spectrometry (LC-MS/MS) and MaxQuant software. A total of 1470 phosphorylation sites in 950 phosphoproteins were identified, and these phosphoproteins were implicated in various molecular functions and basic cellular processes by gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses. Among the 950 phosphoproteins identified, 127 contained 3 to 8 phosphorylation sites. Conservation analysis showed that 93.4% of the 950 phosphoproteins had phosphorylation orthologs in other plant species. Motif-X analysis of the phosphorylation sites identified 13 significantly enriched phosphorylation motifs, of which 3 were novel phosphorylation motifs. Meanwhile, there were 91 phosphoproteins with both multiple phosphorylation sites and multiple phosphorylation motifs. In addition, we identified 58 phosphorylated transcription factors across 21 families and found out 6 significantly over-represented transcription factor families (C3H, Trihelix, CAMTA, TALE, MYB_related and CPP). Eighty-four protein kinases (PKs), 8 protein phosphatases (PPs) and 6 CESAs were recognized as phosphoproteins.

Conclusions

Through a large-scale bioinformatics analysis of the phosphorylation data in seedling leaves, a complicated PKs- and PPs- centered network related to rapid vegetative growth was deciphered in B. distachyon. We revealed a MAPK cascade network that might play the crucial roles during the phosphorylation signal transduction in leaf growth and development. The phosphoproteins and phosphosites identified from our study expanded our knowledge of protein phosphorylation modification in plants, especially in monocots.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-375) contains supplementary material, which is available to authorized users.

Towards revealing the functions of all genes in plants

The great recent progress made in identifying the molecular parts lists of organisms revealed the paucity of our understanding of what most of the parts do. In this review, we introduce computational and statistical approaches and omics data used for inferring gene function in plants, with an emphasis on network-based inference. We also discuss caveats associated with network-based function predictions such as performance assessment, annotation propagation, the guilt-by-association concept, and the meaning of hubs. Finally, we note the current limitations and possible future directions such as the need for gold standard data from several species, unified access to data and tools, quantitative comparison of data and tool quality, and high-throughput experimental validation platforms for systematic gene function elucidation in plants.

Phosphoproteomic studies in Arabidopsis and tobacco male gametophytes

Mature pollen represents an extremely resistant quiescent structure surrounded by a tough cell wall. After its hydration on stigma papillary cells, pollen tube growth starts rapidly. Massive metabolic changes are likely to be accompanied by changes in protein phosphorylation. Protein phosphorylation belongs among the most rapid post-translational modifications. To date, only Arabidopsis thaliana and tobacco (Nicotiana tabacum) mature pollen have been subjected to phosphoproteomic studies in order to identify the phosphoproteins present. In the present mini-review, Arabidopsis and tobacco datasets were compared with each other. The representation of the O-phosphorylated amino acids was compared between these two datasets, and the putative pollen-specific or pollen-abundant phosphopeptides were highlighted. Finally, the phosphorylation sites common for both Arabidopsis and tobacco phosphoproteins are listed as well as the phosphorylation motifs identified.

Proteome and phosphoproteome characterization reveals new response and defense mechanisms of Brachypodium distachyon leaves under salt stress

Salinity is a major abiotic stress affecting plant growth and development. Understanding the molecular mechanisms of salt response and defense in plants will help in efforts to improve the salt tolerance of crops. Brachypodium distachyon is a new model plant for wheat, barley, and several potential biofuel grasses. In the current study, proteome and phosphoproteome changes induced by salt stress were the focus. The Bd21 leaves were initially treated with salt in concentrations ranging from 80 to 320 mm and then underwent a recovery process prior to proteome analysis. A total of 80 differentially expressed protein spots corresponding to 60 unique proteins were identified. The sample treated with a median salt level of 240 mm and the control were selected for phosphopeptide purification using TiO2 microcolumns and LC-MS/MS for phosphoproteome analysis to identify the phosphorylation sites and phosphoproteins. A total of 1509 phosphoproteins and 2839 phosphorylation sites were identified. Among them, 468 phosphoproteins containing 496 phosphorylation sites demonstrated significant changes at the phosphorylation level. Nine phosphorylation motifs were extracted from the 496 phosphorylation sites. Of the 60 unique differentially expressed proteins, 14 were also identified as phosphoproteins. Many proteins and phosphoproteins, as well as potential signal pathways associated with salt response and defense, were found, including three 14-3-3s (GF14A, GF14B, and 14-3-3A) for signal transduction and several ABA signal-associated proteins such as ABF2, TRAB1, and SAPK8. Finally, a schematic salt response and defense mechanism in B. distachyon was proposed.

Genomics and bioinformatics resources for translational science in Rosaceae

Recent technological advances in biology promise unprecedented opportunities for rapid and sustainable advancement of crop quality. Following this trend, the Rosaceae research community continues to generate large amounts of genomic, genetic and breeding data. These include annotated whole genome sequences, transcriptome and expression data, proteomic and metabolomic data, genotypic and phenotypic data, and genetic and physical maps. Analysis, storage, integration and dissemination of these data using bioinformatics tools and databases are essential to provide utility of the data for basic, translational and applied research. This review discusses the currently available genomics and bioinformatics resources for the Rosaceae family.

PAPE (Prefractionation-Assisted Phosphoprotein Enrichment): A Novel Approach for Phosphoproteomic Analysis of Green Tissues from Plants

Phosphorylation is an important post-translational protein modification with regulatory roles in diverse cellular signaling pathways. Despite recent advances in mass spectrometry, the detection of phosphoproteins involved in signaling is still challenging, as protein phosphorylation is typically transient and/or occurs at low levels. In green plant tissues, the presence of highly abundant proteins, such as the subunits of the RuBisCO complex, further complicates phosphoprotein analysis. Here, we describe a simple, but powerful, method, which we named prefractionation-assisted phosphoprotein enrichment (PAPE), to increase the yield of phosphoproteins from Arabidopsis thaliana leaf material. The first step, a prefractionation via ammonium sulfate precipitation, not only depleted RuBisCO almost completely, but, serendipitously, also served as an efficient phosphoprotein enrichment step. When coupled with a subsequent metal oxide affinity chromatography (MOAC) step, the phosphoprotein content was highly enriched. The reproducibility and efficiency of phosphoprotein enrichment was verified by phospho-specific staining and, further, by mass spectrometry, where it could be shown that the final PAPE fraction contained a significant number of known and additionally novel (potential) phosphoproteins. Hence, this facile two-step procedure is a good prerequisite to probe the phosphoproteome and gain deeper insight into plant phosphorylation-based signaling events.

P³DB 3.0: From plant phosphorylation sites to protein networks

In the past few years, the Plant Protein Phosphorylation Database (P³DB, http://p3db.org) has become one of the most significant in vivo data resources for studying plant phosphoproteomics. We have substantially updated P³DB with respect to format, new datasets and analytic tools. In the P³DB 3.0, there are altogether 47 923 phosphosites in 16 477 phosphoproteins curated across nine plant organisms from 32 studies, which have met our multiple quality standards for acquisition of in vivo phosphorylation site data. Centralized by these phosphorylation data, multiple related data and annotations are provided, including protein–protein interaction (PPI), gene ontology, protein tertiary structures, orthologous sequences, kinase/phosphatase classification and Kinase Client Assay (KiC Assay) data—all of which provides context for the phosphorylation event. In addition, P³DB 3.0 incorporates multiple network viewers for the above features, such as PPI network, kinase-substrate network, phosphatase-substrate network, and domain co-occurrence network to help study phosphorylation from a systems point of view. Furthermore, the new P³DB reflects a community-based design through which users can share datasets and automate data depository processes for publication purposes. Each of these new features supports the goal of making P³DB a comprehensive, systematic and interactive platform for phosphoproteomics research.

MASCP gator: an overview of the Arabidopsis proteomic aggregation portal

A key challenge in the area of bioinformatics in the coming decades is the ability to manage the wealth of information that is being generated from the variety of high throughput methodologies currently being undertaken in laboratories across the world. While these approaches have made available large volumes of data to the research community, less attention has been given to the problem of how to intuitively present the data to enable greater biological insights. Recently, an attempt was made to tackle this problem in the area of Arabidopsis proteomics. The model plant has been the target of countless proteomics surveys producing an exhaustive array of data and online repositories. The MASCP Gator is an aggregation portal for proteomic data currently being produced by the community and unites a large collection of specialized resources to a single portal (http://gator.masc-proteomics.org/). Here we describe the latest additions, upgrades and features to this resource further expanding its role into protein modifications and genome sequence variations.

Construction of protein phosphorylation networks by data mining, text mining and ontology integration: analysis of the spindle checkpoint

Knowledge representation of the role of phosphorylation is essential for the meaningful understanding of many biological processes. However, such a representation is challenging because proteins can exist in numerous phosphorylated forms with each one having its own characteristic protein–protein interactions (PPIs), functions and subcellular localization. In this article, we evaluate the current state of phosphorylation event curation and then present a bioinformatics framework for the annotation and representation of phosphorylated proteins and construction of phosphorylation networks that addresses some of the gaps in current curation efforts. The integrated approach involves (i) text mining guided by RLIMS-P, a tool that identifies phosphorylation-related information in scientific literature; (ii) data mining from curated PPI databases; (iii) protein form and complex representation using the Protein Ontology (PRO); (iv) functional annotation using the Gene Ontology (GO); and (v) network visualization and analysis with Cytoscape. We use this framework to study the spindle checkpoint, the process that monitors the assembly of the mitotic spindle and blocks cell cycle progression at metaphase until all chromosomes have made bipolar spindle attachments. The phosphorylation networks we construct, centered on the human checkpoint kinase BUB1B (BubR1) and its yeast counterpart MAD3, offer a unique view of the spindle checkpoint that emphasizes biologically relevant phosphorylated forms, phosphorylation-state–specific PPIs and kinase–substrate relationships. Our approach for constructing protein phosphorylation networks can be applied to any biological process that is affected by phosphorylation.

Database URL:http://www.yeastgenome.org/

Advancements in the analysis of the Arabidopsis plasma membrane proteome

The plasma membrane (PM) regulates diverse processes essential to plant growth, development, and survival in an ever-changing environment. In addition to maintaining normal cellular homeostasis and plant nutrient status, PM proteins perceive and respond to a myriad of environmental cues. Here we review recent advances in the analysis of the plant PM proteome with a focus on the model plant Arabidopsis thaliana. Due to membrane heterogeneity, hydrophobicity, and low relative abundance, analysis of the PM proteome has been a special challenge. Various experimental techniques to enrich PM proteins and different protein and peptide separation strategies have facilitated the identification of thousands of integral and membrane-associated proteins. Numerous classes of proteins are present at the PM with diverse biological functions. PM microdomains have attracted much attention. However, it still remains a challenge to characterize these cell membrane compartments. Dynamic changes in the PM proteome in response to different biotic and abiotic stimuli are highlighted. Future prospects for PM proteomics research are also discussed.

Management and dissemination of MS proteomic data with PROTICdb: example of a quantitative comparison between methods of protein extraction

High throughput MS-based proteomic experiments generate large volumes of complex data and necessitate bioinformatics tools to facilitate their handling. Needs include means to archive data, to disseminate them to the scientific communities, and to organize and annotate them to facilitate their interpretation. We present here an evolution of PROTICdb, a database software that now handles MS data, including quantification. PROTICdb has been developed to be as independent as possible from tools used to produce the data. Biological samples and proteomics data are described using ontology terms. A Taverna workflow is embedded, thus permitting to automatically retrieve information related to identified proteins by querying external databases. Stored data can be displayed graphically and a "Query Builder" allows users to make sophisticated queries without knowledge on the underlying database structure. All resources can be accessed programmatically using a Java client API or RESTful web services, allowing the integration of PROTICdb in any portal. An example of application is presented, where proteins extracted from a maize leaf sample by four different methods were compared using a label-free shotgun method. Data are available at http://moulon.inra.fr/protic/public. PROTICdb thus provides means for data storage, enrichment, and dissemination of proteomics data.

Biochemistry, proteomics, and phosphoproteomics of plant mitochondria from non-photosynthetic cells

Mitochondria fulfill some basic roles in all plant cells. They supply the cell with energy in the form of ATP and reducing equivalents [NAD(P)H] and they provide the cell with intermediates for a range of biosynthetic pathways. In addition to this, mitochondria contribute to a number of specialized functions depending on the tissue and cell type, as well as environmental conditions. We will here review the biochemistry and proteomics of mitochondria from non-green cells and organs, which differ from those of photosynthetic organs in a number of respects. We will briefly cover purification of mitochondria and general biochemical properties such as oxidative phosphorylation. We will then mention a few adaptive properties in response to water stress, seed maturation and germination, and the ability to function under hypoxic conditions. The discussion will mainly focus on Arabidopsis cell cultures, etiolated germinating rice seedlings and potato tubers as model plants. It will cover the general proteome as well as the posttranslational modification protein phosphorylation. To date 64 phosphorylated mitochondrial proteins with a total of 103 phosphorylation sites have been identified.