Phosphate and phosphite have a differential impact on the proteome and phosphoproteome of Arabidopsis suspension cell cultures

Unveiling unique alternative splicing responses to low temperature in Zoysia japonica through ZjRTD1.0, a high-quality reference transcript dataset

Inadequate reference databases in RNA-seq analysis can hinder data utilization and interpretation. In this study, we have successfully constructed a high-quality reference transcript dataset, ZjRTD1.0, for Zoysia japonica, a widely-used turfgrass with exceptional tolerance to various abiotic stress, including low temperatures and salinity. This dataset comprises 113,089 transcripts from 57,143 genes. BUSCO analysis demonstrates exceptional completeness (92.4%) in ZjRTD1.0, with reduced proportions of fragmented (3.3%) and missing (4.3%) orthologs compared to prior datasets. ZjRTD1.0 enables more precise analyses, including transcript quantification and alternative splicing assessments using public datasets, which identified a substantial number of differentially expressed transcripts (DETs) and differential alternative splicing (DAS) events, leading to several novel findings on Z. japonica's responses to abiotic stresses. First, spliceosome gene expression influenced alternative splicing significantly under abiotic stress, with a greater impact observed during low-temperature stress. Then, a significant positive correlation was found between the number of differentially expressed genes (DEGs) encoding protein kinases and the frequency of DAS events, suggesting the role of protein phosphorylation in regulating alternative splicing. Additionally, our results suggest possible involvement of serine/arginine-rich (SR) proteins and heterogeneous nuclear ribonucleoproteins (hnRNPs) in generating inclusion/exclusion isoforms under low-temperature stress. Furthermore, our investigation revealed a significantly enhanced overlap between DEGs and differentially alternatively spliced genes (DASGs) in response to low-temperature stress, suggesting a unique co-regulatory mechanism governing transcription and splicing in the context of low-temperature response. In conclusion, we have proven that ZjRTD1.0 will serve as a reliable and useful resource for future transcriptomic analyses in Z. japonica.

Integrated proteomic analysis reveals interactions between phosphorylation and ubiquitination in rose response to Botrytis infection

As two of the most abundant post-translational modifications, phosphorylation and ubiquitination play a significant role in modulating plant-pathogen interactions and increasing evidence indicates their crosstalk in plant immunity. Rose (Rosa sp.) is one of the most important ornamental plants and can be seriously infected by Botrytis cinerea. Here, integrated proteomics analysis was performed to detect global proteome, phosphorylation, and ubiquitination changes in rose upon B. cinerea infection and investigate the possible phosphorylation and ubiquitination crosstalk. A total of 6165 proteins, 11 774 phosphorylation and 10 582 ubiquitination sites, and 77 phosphorylation and 13 ubiquitination motifs were identified. Botrytis cinerea infection resulted in 169 up-regulated and 122 down-regulated proteins, 291 up-regulated and 404 down-regulated phosphorylation sites, and 250 up-regulated and 634 down-regulated ubiquitination sites. There were 12 up-regulated PR10 proteins and half of them also showed reduced ubiquitination. A lot of kinases probably involved in plant pattern-triggered immunity signaling were up-regulated phosphoproteins. Noticeably, numerous kinases and ubiquitination-related proteins also showed a significant change in ubiquitination and phosphorylation, respectively. A cross-comparison of phosphoproteome and ubiquitylome indicated that both of two post-translational modifications of 104 proteins were dynamically regulated, and many putative pattern-triggered immunity signaling components in the plant plasma membrane were co-regulated. Moreover, five selected proteins, including four PR10 proteins and a plasma membrane aquaporin, were proven to be involved in rose resistance to B. cinerea. Our study provides insights into the molecular mechanisms underlying rose resistance to B. cinerea and also increases the database of phosphorylation and ubiquitination sites in plants.

Impaired glycosylation of GmPAP15a, a root-associated purple acid phosphatase, inhibits extracellular phytate-P utilization in soybean

Phosphorus (P) is an essential nutrient, but easily fixed in soils. Therefore, most of soil P exists in the form of inaccessible organic phosphorus (Po), particularly phytate-P. Root-associated purple acid phosphatases (PAPs) are considered to play a crucial role in phosphate (Pi) scavenging in soils. However, evidence for regulating root-associated PAPs in utilization of extracellular phytate-P remain largely unknown in plants at both transcriptional and posttranslational levels. In this study, a Pi-starvation responsive GmPAP15a was identified in soybean (Glycine max). Overexpressing GmPAP15a led to significant increases in root-associated phytase activities, as well as total P content when phytate-P was supplied as the sole P resource in soybean hairy roots. Meanwhile, mass spectrometry (MS) analysis showed GmPAP15a was glycosylated at Asn144 and Asn502 , and its glycan structures of N-linked oligosaccharide chains exhibited microheterogeneity. Moreover, two homologues of AtPHR1, GmPHR9 and GmPHR32 were found to activate GmPAP15a transcription through luciferase activity analysis. Taken together, it is strongly suggested that GmPAP15a plays a vital role in phytate-P utilization in soybean, which might be regulated at both transcriptional and glycosylation modification levels. Our results highlight the GmPHR9/GmPHR32-GmPAP15a signalling pathway might present, and control phytate-P utilization in soybean.

In Vivo Phosphorylation of the Cytosolic Glucose-6-Phosphate Dehydrogenase Isozyme G6PD6 in Phosphate-Resupplied Arabidopsis thaliana Suspension Cells and Seedlings

Glucose-6-phosphate dehydrogenase (G6PD) catalyzes the first committed step of the oxidative pentose phosphate pathway (OPPP). Our recent phosphoproteomics study revealed that the cytosolic G6PD6 isozyme became hyperphosphorylated at Ser12, Thr13 and Ser18, 48 h following phosphate (Pi) resupply to Pi-starved (-Pi) Arabidopsis thaliana cell cultures. The aim of the present study was to assess whether G6PD6 phosphorylation also occurs in shoots or roots following Pi resupply to -Pi Arabidopsis seedlings, and to investigate its relationship with G6PD activity. Interrogation of phosphoproteomic databases indicated that N-terminal, multi-site phosphorylation of G6PD6 and its orthologs is quite prevalent. However, the functions of these phosphorylation events remain unknown. Immunoblotting with an anti-(pSer18 phosphosite-specific G6PD6) antibody confirmed that G6PD6 from Pi-resupplied, but not -Pi, Arabidopsis cell cultures or seedlings (i.e., roots) was phosphorylated at Ser18; this correlated with a significant increase in extractable G6PD activity, and biomass accumulation. Peptide kinase assays of Pi-resupplied cell culture extracts indicated that G6PD6 phosphorylation at Ser18 is catalyzed by a Ca2+-dependent protein kinase (CDPK), which correlates with the 'CDPK-like' targeting motif that flanks Ser18. Our results support the hypothesis that N-terminal phosphorylation activates G6PD6 to enhance OPPP flux and thus the production of reducing power (i.e., NADPH) and C-skeletons needed to establish the rapid resumption of growth that ensures Pi-resupply to -Pi Arabidopsis.

Glutamate decarboxylase-1 is essential for efficient acclimation of Arabidopsis thaliana to nutritional phosphorus deprivation

Glutamate decarboxylase (GAD) is a Ca2+ -calmodulin-activated, cytosolic enzyme that produces γ-aminobutyrate (GABA) as the committed step of the GABA shunt. This pathway bypasses the 2-oxoglutarate to succinate reactions of the tricarboxylic acid (TCA) cycle. GABA also accumulates during many plant stresses. We tested the hypothesis that AtGAD1 (At5G17330) facilitates Arabidopsis acclimation to Pi deprivation. Quantitative RT-PCR and immunoblotting revealed that AtGAD1 transcript and protein expression is primarily root-specific, but inducible at lower levels in shoots of Pi-deprived (-Pi) plants. Pi deprivation reduced levels of the 2-oxoglutarate dehydrogenase (2-OGDH) cofactor thiamine diphosphate (ThDP) in shoots and roots by > 50%. Growth of -Pi atgad1 T-DNA mutants was significantly attenuated relative to wild-type plants. This was accompanied by: (i) an > 60% increase in shoot and root GABA levels of -Pi wild-type, but not atgad1 plants, and (ii) markedly elevated anthocyanin and reduced free and total Pi levels in leaves of -Pi atgad1 plants. Treatment with 10 mM GABA reversed the deleterious development of -Pi atgad1 plants. Our results indicate that AtGAD1 mediates GABA shunt upregulation during Pi deprivation. This bypass is hypothesized to circumvent ThDP-limited 2-OGDH activity to facilitate TCA cycle flux and respiration by -Pi Arabidopsis.

Phosphorylation mediated regulation of RNA splicing in plants

For the past two decades, the study of alternative splicing (AS) and its involvement in plant development and stress response has grown in popularity. Only recently however, has the focus shifted to the study of how AS regulation (or lack-thereof) affects downstream mRNA and protein landscapes and how these AS regulatory events impact plant development and stress tolerance. In humans, protein phosphorylation represents one of the predominant mechanisms by which AS is regulated and thus the protein kinases governing these phosphorylation events are of interest for further study. Large-scale phosphoproteomic studies in plants have consistently found that RNA splicing-related proteins are extensively phosphorylated, however, the signaling pathways involved in AS regulation have not been resolved. In this mini-review, we summarize our current knowledge of the three major splicing-related protein kinase families in plants that are suggested to mediate AS phospho-regulation and draw comparisons to their metazoan orthologs. We also summarize and contextualize the phosphorylation events identified as occurring on splicing-related protein families to illustrate the high degree to which splicing-related proteins are modified, placing a new focus on elucidating the impacts of AS at the protein and PTM-level.

Functional assessment of AtPAP17; encoding a purple acid phosphatase involved in phosphate metabolism in Arabidopsis thaliana

PURPOSE

Purple acid phosphatases (PAPs) includ the largest classes of non-specific plant acid phosphatases. Most characterized PAPs were found to play physiological functions in phosphorus metabolism. In this study, we investigated the function of AtPAP17 gene encoding an important purple acid phosphatase in Arabidopsis thaliana.

METHODS

The full-length cDNA sequence of AtPAP17 gene under the control of CaMV-35S promoter was transferred to the A. thaliana WT plant. The generated homozygote AtPAP17-overexpressed plants were compared by the types of analyses with corresponding homozygote atpap17-mutant plant and WT in both + P (1.2 mM) and - P (0 mM) conditions.

RESULTS

In the + P condition, the highest and the lowest amount of Pi was observed in AtPAP17-overexpressed plants and atpap17-mutant plants by 111% increase and 38% decrease compared with the WT plants, respectively. Furthermore, under the same condition, APase activity of AtPAP17-overexpressed plants increased by 24% compared to the WT. Inversely, atpap17-mutant plant represented a 71% fall compared to WT plants. The comparison of fresh weight and dry weight in the studied plants showed that the highest and the lowest amount of absorbed water belonged to OE plants (with 38 and 12 mg plant^-1) and Mu plants (with 22 and 7 mg plant^-1) in + P and - P conditions, respectively.

CONCLUSION

The lack of AtPAP17 gene in the A. thaliana genome led to a remarkable reduction in the development of root biomass. Thus, AtPAP17 could have an important role in the root but not shoot developmental and structural programming. Consequently, this function enables them to absorb more water and eventually associated with more phosphate absorption.

Five Post-Translational Modification Residues of CmPT2 Play Key Roles in Yeast and Rice

Chrysanthemum (Chrysanthemum morifolium Ramat.) is one of the largest cut flowers in the world. Phosphate transporter Pht1 family member CmPht1;2 protein (CmPT2) plays an important role in response to low-phosphate (LP) stress in chrysanthemum. Post-translational modification (PTM) can modulate the function of proteins in multiple ways. Here, we used yeast and rice systems to study the role of putative PTM in CmPT2 by determining the effect of mutation of key amino acid residues of putative glycosylation, phosphorylation, and myristoylation sites. We chose nine amino acid residues in the putative PTM sites and mutated them to alanine (A) (Cmphts). CmPT2 recovered the growth of yeast strain MB192 under LP conditions. However, G84A, G222A, T239A, Y242A, and N422A mutants could not grow normally under LP conditions. Analysis of phosphorus absorption kinetics showed that the Km of CmPT2 was 65.7 μM. Among the nine Cmphts, the expression of five with larger Km (124.4-397.5 μM) than CmPT2 was further evaluated in rice. Overexpression of CmPT2-OE increased plant height, effective panicle numbers, branch numbers, and yield compared with that of wild type 'Wuyunjing No. 7' (W7). Overexpression of Cmphts-OE led to decreased plant height and effective panicle numbers compared with that of the CmPT2-OE strain. The Pi content in roots of CmPT2-OE was higher than that of the W7 under both high (normal) phosphate (HP) and LP conditions. However, the Pi content in the leaves and roots was significantly lower in the N422A-OE strain than in the CmPT2-OE strain under both HP and LP conditions. Under LP conditions, the phosphorus starvation response (PSR) genes in CmPT2-OE were inhibited at the transcription level. The expression patterns of phosphorus-related genes in T239A, Y242A, and N422A-OE under LP conditions were different from those of CmPT2-OE. In conclusion, these five post-translational modification residues of CmPT2 play key roles in modulating the function of CmPT2. This work boosters our understanding of the function of phosphate transporters and provides genetic resources for improving the efficiency of phosphorus utilization in crop plants.

Dissection of the response mechanism of alfalfa under phosphite stress based on metabolomic and transcriptomic data

Phosphite, a reduced form of phosphate, inhibits the growth and even has toxic effect on plants. To learn more about the mechanism of alfalfa responses to phosphite, the morphological and physiological characteristics, and the metabolites and transcript levels were comprehensively analyzed following the exposure of alfalfa seedlings to phosphite and phosphate under greenhouse conditions. The results showed that phosphite inhibited seedling growth and photosynthesis. However, the absorption efficiency of phosphite was higher than that of phosphate in roots, which was supported by increased total phosphorus concentration of 16.29% and 52.30% on days 8 and 12. Moreover, phosphite stress affected the synthesis of lipids and carbohydrates, which were reflected in enhanced glycolipid and sulfolipid in roots and amylose in shoots. Phosphite stress resulted in a decrease in indole acetic acid (IAA) in the whole plant and zeatin in the shoots, which could enable alfalfa to adapt to the phosphite environment. Some genes involved in phosphate starvation response included SPX, phosphate response regulator2, and inorganic phosphate transporter 1-4 (PHT1;4) in roots were affected by phosphite stress. In addition, some genes that are involved in stress responses and DNA repair were induced by phosphite stress. These observations together suggest that alfalfa responds to phosphite stress by inhibiting growth, regulating the genes induced by phosphate starvation, improving oxidative protection, promoting DNA repair, and adjusting the IAA and zeatin signaling transductions. Our findings provide novel insights into the molecular response to phosphite stress in alfalfa.

Phosphorylation of S11 in PHR1 negatively controls its transcriptional activity

Plant responses to phosphate starvation (-Pi) are very well characterized at the biochemical and molecular levels. The expression of thousands of genes is modified under this stress condition, depending on the action of Phosphate starvation response 1 (PHR1). Existing data indicate that neither the PHR1 transcript nor the quantity or localization of its protein increase during nutrient stress, raising the question of how its activity is regulated. Here, we present data showing that SnRK1 kinase is able to phosphorylate some phosphate starvation response proteins (PSRs), including PHR1. Based on a model of the three-dimensional structure of the catalytic subunit SnRK1α1, docking simulations predicted the binding modes of peptides from PHT1;8, PHO1 and PHR1 with SnRK1. PHR1 recombinant protein interacted in vitro with the catalytic subunits SnRK1α1 and SnRK1α2. A BiFC assay corroborated the in vivo interaction between PHR1 and SnRK1α1 in the cytoplasm and nucleus. Analysis of phosphorylated residues suggested the presence of one phosphorylated site containing the SnRK1 motif at S11, and mutation in this residue disrupted the incorporation of ³² P, suggesting that it is a major phosphorylation site. Electrophoretic mobility shift assay results indicated that the binding of PHR1 to P1BS motifs was not influenced by phosphorylation. Importantly, transient expression assays in Arabidopsis protoplasts showed a decrease in PHR1 activity in contrast with the S11A mutant, suggesting a role for Ser11 as a negative regulatory phosphorylation site. Taken together, these findings suggest that phosphorylation of PHR1 at Ser11 is a mechanism to control the PHR1-mediated adaptive response to -Pi.

Multi-omic analysis shows REVEILLE clock genes are involved in carbohydrate metabolism and proteasome function

Plants are able to sense changes in their light environments, such as the onset of day and night, as well as anticipate these changes in order to adapt and survive. Central to this ability is the plant circadian clock, a molecular circuit that precisely orchestrates plant cell processes over the course of a day. REVEILLE (RVE) proteins are recently discovered members of the plant circadian circuitry that activate the evening complex and PSEUDO-RESPONSE REGULATOR genes to maintain regular circadian oscillation. The RVE8 protein and its two homologs, RVE 4 and 6 in Arabidopsis (Arabidopsis thaliana), have been shown to limit the length of the circadian period, with rve 4 6 8 triple-knockout plants possessing an elongated period along with increased leaf surface area, biomass, cell size, and delayed flowering relative to wild-type Col-0 plants. Here, using a multi-omics approach consisting of phenomics, transcriptomics, proteomics, and metabolomics we draw new connections between RVE8-like proteins and a number of core plant cell processes. In particular, we reveal that loss of RVE8-like proteins results in altered carbohydrate, organic acid, and lipid metabolism, including a starch excess phenotype at dawn. We further demonstrate that rve 4 6 8 plants have lower levels of 20S proteasome subunits and possess significantly reduced proteasome activity, potentially explaining the increase in cell-size observed in RVE8-like mutants. Overall, this robust, multi-omic dataset provides substantial insight into the far-reaching impact RVE8-like proteins have on the diel plant cell environment.

Proteomic Analysis Dissects Molecular Mechanisms Underlying Plant Responses to Phosphorus Deficiency

Phosphorus (P) is an essential nutrient for plant growth. In recent decades, the application of phosphate (Pi) fertilizers has contributed to significant increases in crop yields all over the world. However, low efficiency of P utilization in crops leads to intensive application of Pi fertilizers, which consequently stimulates environmental pollution and exhaustion of P mineral resources. Therefore, in order to strengthen the sustainable development of agriculture, understandings of molecular mechanisms underlying P efficiency in plants are required to develop cultivars with high P utilization efficiency. Recently, a plant Pi-signaling network was established through forward and reverse genetic analysis, with the aid of the application of genomics, transcriptomics, proteomics, metabolomics, and ionomics. Among these, proteomics provides a powerful tool to investigate mechanisms underlying plant responses to Pi availability at the protein level. In this review, we summarize the recent progress of proteomic analysis in the identification of differential proteins that play roles in Pi acquisition, translocation, assimilation, and reutilization in plants. These findings could provide insights into molecular mechanisms underlying Pi acquisition and utilization efficiency, and offer new strategies in genetically engineering cultivars with high P utilization efficiency.

Arabidopsis PAP17 is a dual-localized purple acid phosphatase up-regulated during phosphate deprivation, senescence, and oxidative stress

A 35 kDa monomeric purple acid phosphatase (APase) was purified from cell wall extracts of Pi starved (-Pi) Arabidopsis thaliana suspension cells and identified as AtPAP17 (At3g17790) by mass spectrometry and N-terminal microsequencing. AtPAP17 was de novo synthesized and dual-localized to the secretome and/or intracellular fraction of -Pi or salt-stressed plants, or senescing leaves. Transiently expressed AtPAP17-green fluorescent protein localized to lytic vacuoles of the Arabidopsis suspension cells. No significant biochemical or phenotypical changes associated with AtPAP17 loss of function were observed in an atpap17 mutant during Pi deprivation, leaf senescence, or salinity stress. Nevertheless, AtPAP17 is hypothesized to contribute to Pi metabolism owing to its marked up-regulation during Pi starvation and leaf senescence, broad APase substrate selectivity and pH activity profile, and rapid repression and turnover following Pi resupply to -Pi plants. While AtPAP17 also catalyzed the peroxidation of luminol, which was optimal at pH 9.2, it exhibited a low Vmax and affinity for hydrogen peroxide relative to horseradish peroxidase. These results, coupled with absence of a phenotype in the salt-stressed or -Pi atpap17 mutant, do not support proposals that the peroxidase activity of AtPAP17 contributes to the detoxification of reactive oxygen species during stresses that trigger AtPAP17 up-regulation.

ptxD/Phi as alternative selectable marker system for genetic transformation for bio-safety concerns: a review

Antibiotic and herbicide resistance genes are the most common marker genes for plant transformation to improve crop yield and food quality. However, there is public concern about the use of resistance marker genes in food crops due to the risk of potential gene flow from transgenic plants to compatible weedy relatives, leading to the possible development of “superweeds” and antibiotic resistance. Several selectable marker genes such as aph, nptII, aaC3, aadA, pat, bar, epsp and gat, which have been synthesized to generate transgenic plants by genetic transformation, have shown some limitations. These marker genes, which confer antibiotic or herbicide resistance and are introduced into crops along with economically valuable genes, have three main problems: selective agents have negative effects on plant cell proliferation and differentiation, uncertainty about the environmental effects of many selectable marker genes, and difficulty in performing recurrent transformations with the same selectable marker to pyramid desired genes. Recently, a simple, novel, and affordable method was presented for plant cells to convert non-metabolizable phosphite (Phi) to an important phosphate (Pi) for developing cells by gene expression encoding a phosphite oxidoreductase (PTXD) enzyme. The ptxD gene, in combination with a selection medium containing Phi as the sole phosphorus (P) source, can serve as an effective and efficient system for selecting transformed cells. The selection system adds nutrients to transgenic plants without potential risks to the environment. The ptxD/Phi system has been shown to be a promising transgenic selection system with several advantages in cost and safety compared to other antibiotic-based selection systems. In this review, we have summarized the development of selection markers for genetic transformation and the potential use of the ptxD/Phi scheme as an alternative selection marker system to minimize the future use of antibiotic and herbicide marker genes.

Proteomics and plant biology: contributions to date and a look towards the next decade

INTRODUCTION

This review presents the view of the author, that is opinionable and even speculative, on the field of proteomics, its application to plant biology knowledge, and translation to biotechnology. Written in a more academic than scientific style, it is based on past original and review articles by the author´s group, and those published by leading scientists in the last two years.

AREAS COVERED

Starting with a general definition and references to historical milestones, it covers sections devoted to the different platforms employed, the plant biology discourse in the protein language, challenges and future prospects, ending with the author opinion.

EXPERT OPINION

In 25 years, five proteomics platform generations have appeared. We are now moving from proteomics to Systems Biology. While feasible with model organisms, proteomics of orphan species remains challenging. Proteomics, even in its simplest approach, sheds light on plant biological processes, central dogma, and molecular bases of phenotypes of interest, and it can be translated to areas such as food traceability and allergen detection. Proteomics should be validated and optimized to each experimental system, objectives, and hypothesis. It has limitations, artifacts, and biases. We should not blindly accept proteomics data and just create a list of proteins, networks, and avoid speculative biological interpretations. From the hundred to thousand proteins identified and quantified, it is important to obtain a focus and validate some of them, otherwise it is merely. We are starting to have the protein pieces, so let, from now, build the proteomics and biological puzzle.

Recent insights into the metabolic adaptations of phosphorus-deprived plants

Inorganic phosphate (Pi) is an essential macronutrient required for many fundamental processes in plants, including photosynthesis and respiration, as well as nucleic acid, protein, and membrane phospholipid synthesis. The huge use of Pi-containing fertilizers in agriculture demonstrates that the soluble Pi levels of most soils are suboptimal for crop growth. This review explores recent advances concerning the understanding of adaptive metabolic processes that plants have evolved to alleviate the negative impact of nutritional Pi deficiency. Plant Pi starvation responses arise from complex signaling pathways that integrate altered gene expression with post-transcriptional and post-translational mechanisms. The resultant remodeling of the transcriptome, proteome, and metabolome enhances the efficiency of root Pi acquisition from the soil, as well as the use of assimilated Pi throughout the plant. We emphasize how the up-regulation of high-affinity Pi transporters and intra- and extracellular Pi scavenging and recycling enzymes, organic acid anion efflux, membrane remodeling, and the remarkable flexibility of plant metabolism and bioenergetics contribute to the survival of Pi-deficient plants. This research field is enabling the development of a broad range of innovative and promising strategies for engineering phosphorus-efficient crops. Such cultivars are urgently needed to reduce inputs of unsustainable and non-renewable Pi fertilizers for maximum agronomic benefit and long-term global food security and ecosystem preservation.

Global Profiling of Phosphorylation Reveals the Barley Roots Response to Phosphorus Starvation and Resupply

Phosphorus (P) deficiency is a major threat to the crop production, and for understanding the response mechanism of plant roots, P stress may facilitate the development of crops with increased tolerance. Phosphorylation plays a critical role in the regulation of proteins for plant responses to biotic and abiotic stress; however, its functions in P starvation/resupply are largely unknown for barley (Hordeum vulgare) growth. Here, we performed a global review of phosphorylation in barley roots treated by P starvation/resupply. We identified 7,710 phosphorylation sites on 3,373 proteins, of which 76 types of conserved motifs were extracted from 10,428 phosphorylated peptides. Most phosphorylated proteins were located in the nucleus (36%) and chloroplast (32%). Compared with the control, 186 and 131 phosphorylated proteins under P starvation condition and 156 and 111 phosphorylated proteins under P resupply condition showed significant differences at 6 and 48 h, respectively. These proteins mainly participated in carbohydrate metabolism, phytohormones, signal transduction, cell wall stress, and oxidases stress. Moreover, the pathways of the ribosome, RNA binding, protein transport, and metal binding were significantly enriched under P starvation, and only two pathways of ribosome and RNA binding were greatly enriched under Pi resupply according to the protein-protein interaction analysis. The results suggested that the phosphorylation proteins might play important roles in the metabolic processes of barley roots in response to Pi deficiency/resupply. The data not only provide unique access to phosphorylation reprogramming of plant roots under deficiency/resupply but also demonstrate the close cooperation between these phosphorylation proteins and key metabolic functions.