Advances in proteome-wide analysis of plant lysine acetylation

The regulation mechanism of ethephon-mediated delaying of postharvest physiological deterioration in cassava storage roots based on quantitative acetylproteomes analysis

Ethylene plays diverse roles in post-harvest processes of horticultural crops. However, its impact and regulation mechanism on the postharvest physiological deterioration (PPD) of cassava storage roots is unknown. In this study, a notable delay in PPD of cassava storage roots was observed when ethephon was utilized as an ethylene source. Physiological analyses and quantitative acetylproteomes were employed to investigate the regulation mechanism regulating cassava PPD under ethephon treatment. Ethephon was found to enhance the reactive oxygen species (ROS) scavenging system, resulting in a significant decrease in H2O2 and malondialdehyde (MDA) content. The comprehensive acetylome analysis identified 12,095 acetylation sites on 4403 proteins. Subsequent analysis demonstrated that ethephon can regulate the acetylation levels of antioxidant enzymes and members of the energy metabolism pathways. In summary, ethephon could enhance the antioxidant properties and regulate energy metabolism pathways, leading to the delayed PPD of cassava.

The conserved transcriptional regulation mechanism of ADH1 gene in Zanthoxylum armatum to waterlogging stress

Waterlogging stress negatively affects plant growth and survival. However, the ability of Zanthoxylum armatum, a valuable tree species, to tolerate and adapt to waterlogging stress remains poorly understood. Here we report how alcohol dehydrogenase 1 (ZaADH1) confers waterlogging stress tolerance in Z. armatum. ZaADH1 expression was induced after waterlogging treatment. ZaADH1 overexpression increased waterlogging stress by modulating the metabolite levels of the ADH enzyme, soluble sugar, and trehalose, promoting glycolysis and carbohydrate metabolism. The overexpression of ZaADH1 in Arabidopsis thaliana increased the total plant area and chlorophyll content, thereby increasing resistance to waterlogging stress. Physiological and overexpression transcriptome analyses in A. thaliana indicated that ZaADH1 overexpressing lines generated more carbohydrates to meet energy demands, employing a "static" strategy to increase tolerance to waterlogging stress, which confirms the conservation of the ADH1 response to waterlogging stress and represents a potential crucial measure for improving waterlogging tolerance in Z. armatum.

Advances in understanding the roles of plant HAT and HDAC in non-histone protein acetylation and deacetylation

MAIN CONCLUSION

This review focuses on HATs and HDACs that modify non-histone proteins, summarizes functional mechanisms of non-histone acetylation as well as the roles of HATs and HDACs in rice and Arabidopsis. The growth and development of plants, as well as their responses to biotic and abiotic stresses, are governed by intricate gene and protein regulatory networks, in which epigenetic modifying enzymes play a crucial role. Histone lysine acetylation levels, modulated by histone acetyltransferases (HATs) and histone deacetylases (HDACs), are well-studied in the realm of transcriptional regulation. However, the advent of advanced proteomics has unveiled that non-histone proteins also undergo acetylation, with its underlying mechanisms now being clarified. Indeed, non-histone acetylation influences protein functionality through diverse pathways, such as modulating protein stability, adjusting enzymatic activity, steering subcellular localization, influencing interactions with other post-translational modifications, and managing protein-protein and protein-DNA interactions. This review delves into the recent insights into the functional mechanisms of non-histone acetylation in plants. We also provide a summary of the roles of HATs and HDACs in rice and Arabidopsis, and explore their potential involvement in the regulation of non-histone proteins.

Lysine acetylproteome analysis reveals the lysine acetylation in developing fruit and a key acetylated protein involved in sucrose accumulation in Rosa roxburghii Tratt

Lysine acetylation is a common post-translational modification of proteins in plants. Rosa roxburghii Tratt. is an economically important fruit tree known for its high nutritional value. However, the characteristics of acetylome-related proteins during fruit development in this crop remain unknown. This study aimed to explore the global acetylproteome of R. roxburghii fruit to identify key lysine-acetylated proteins associated with its quality traits. A total of 4280 acetylated proteins were identified, among them, 981 proteins exhibited differential acetylation (DA) while 19 proteins showed increased acetylation level consistently on individual sites. Functional classification revealed that these DA proteins were primarily associated with central metabolic pathways, carbohydrate metabolism, terpenoids and polyketides metabolism, lipid metabolism, and amino acid metabolism, highlighting the importance of lysine acetylation in fruit quality formation. Notably, the most significant up-regulated acetylation occurred in sucrose synthase (SuS1), a key enzyme in sucrose biosynthesis. Enzyme assays, RNA-seq and proteome analysis indicated that SuS activity, which was independent of its transcriptome and proteome level, may be enhanced by up-acetylation, ultimately increasing sucrose accumulation. Thus, these findings offer a better understanding of the global acetylproteome of R. roxburghii fruit, while also uncover a novel mechanism of acetylated SuS-mediated in sucrose metabolism in plant. SIGNIFICANCE: Rosa roxburghii Tratt. is an important horticultural crop whose commercial value is closely linked to its fruit quality. Acetylation modification is a post-translational mechanism observed in plants, which regulates the physiological functions and metabolic fluxes involved in various biological processes. The regulatory mechanism of lysine acetylation in the fruit quality formation in perennial woody plants has not been fully elucidated, while most of the research has primarily focused on annual crops. Therefore, this study, for the first time, uses Rosaceae fruits as the research material to elucidate the regulatory role of lysine-acetylated proteins in fruit development, identify key metabolic processes influencing fruit quality formation, and provide valuable insights for cultivation strategies.

Sexually dimorphic acetyl-CoA biosynthesis and utilization in response to drought and exogenous acetic acid

Female willows exhibit greater drought tolerance and benefit more from exogenous acetic acid (AA)-improved drought tolerance than males. However, the potential mechanisms driving these sex-specific responses remain unclear. To comprehensively investigate the sexually dimorphic responsive mechanisms of willows to drought and exogenous AA, here, we performed physiological, proteomic, Lys-acetylproteomic, and transgenic analyses in female and male Salix myrtillacea exposed to drought and AA-applicated drought treatments, focusing on protein abundance and lysine acetylation (LysAc) changes. Drought-tolerant females suffered less drought-induced photosynthetic and oxidative damage, did not activate AA and acetyl-CoA biosynthesis, TCA cycle, fatty acid metabolism, and jasmonic acid signaling as strongly as drought-sensitive males. Exogenous AA caused overaccumulation of endogenous AA and inhibition of acetyl-CoA biosynthesis and utilization in males. However, exogenous AA greatly enhanced acetyl-CoA biosynthesis and utilization and further enhanced drought performance of females, possibly determining that AA improved drought tolerance more in females than in males. Interestingly, overexpression of acetyl-CoA synthetase (ACS) could reprogram fatty acids, increase LysAc levels, and improve drought tolerance, highlighting the involvement of ACS-derived acetyl-CoA in drought responses. In addition, drought and exogenous AA induced sexually dimorphic LysAc associated with histones, transcription factors, and metabolic enzymes in willows. Especially, exogenous AA may greatly improve the photosynthetic capacity of S. myrtillacea males by decreasing LysAc levels and increasing the abundances of photosynthetic proteins. While hyperacetylation in glycolysis, TCA cycle, and fatty acid biosynthesis potentially possibly serve as negative feedback to acclimate acetyl-CoA biosynthesis and utilization in drought-stressed males and AA-applicated females. Thus, acetyl-CoA biosynthesis and utilization determine the sexually dimorphic responses of S. myrtillacea to drought and exogenous AA.

Research Progress on the Mechanism and Function of Histone Acetylation Regulating the Interaction between Pathogenic Fungi and Plant Hosts

Histone acetylation is a crucial epigenetic modification, one that holds the key to regulating gene expression by meticulously modulating the conformation of chromatin. Most histone acetylation enzymes (HATs) and deacetylation enzymes (HDACs) in fungi were originally discovered in yeast. The functions and mechanisms of HATs and HDACs in yeast that have been documented offer us an excellent entry point for gaining insights into these two types of enzymes. In the interaction between plants and pathogenic fungi, histone acetylation assumes a critical role, governing fungal pathogenicity and plant immunity. This review paper delves deep into the recent advancements in understanding how histone acetylation shapes the interaction between plants and fungi. It explores how this epigenetic modification influences the intricate balance of power between these two kingdoms of life, highlighting the intricate network of interactions and the subtle shifts in these interactions that can lead to either mutual coexistence or hostile confrontation.

Acetylation of transcription factor BpTCP20 by acetyltransferase BpPDCE23 modulates salt tolerance in birch

Teosinte branched 1/Cycloidea/Proliferating cell factor (TCP) transcription factors function in abiotic stress responses. However, how TCPs confer salt tolerance is unclear. Here, we characterized a TCP transcription factor, BpTCP20, that responds to salt stress in birch (Betula platyphylla Suk). Plants overexpressing BpTCP20 displayed increased salt tolerance, and Bptcp20 knockout mutants displayed reduced salt tolerance relative to the wild-type (WT) birch. BpTCP20 conferred salt tolerance by mediating stomatal closure and reducing reactive oxygen species (ROS) accumulation. Chromatin immunoprecipitation sequencing showed that BpTCP20 binds to NeuroD1, T-box, and two unknown elements (termed TBS1 and TBS2) to regulate target genes. In birch, salt stress led to acetylation of BpTCP20 acetylation at lysine 259. A mutated BpTCP20 variant (abolished for acetylation, termed BpTCP20259) was overexpressed in birch, which led to decreased salt tolerance compared with plants overexpressing BpTCP20. However, BpTCP20259-overexpressing plants still displayed increased salt tolerance relative to untransformed WT plants. BpTCP20259 showed reduced binding to the promoters of target genes and decreased target gene activation, leading to decreased salt tolerance. In addition, we identified dihydrolipoyllysine-residue acetyltransferase component of pyruvate dehydrogenase complex (BpPDCE23), an acetyltransferase that interacts with and acetylates BpTCP20 to enhance its binding to DNA motifs. Together, these results suggest that BpTCP20 is a transcriptional regulator of salt tolerance, whose activity is modulated by BpPDCE23-mediated acetylation.

Comprehensive Proteome and Acetylome Analysis of Needle Senescence in Larix gmelinii

Leaf senescence is essential for the growth and development of deciduous trees in the next season. Larix gmelinii, a deciduous coniferous tree, exhibits its most distinctive feature by turning yellow in the autumn and eventually shedding its leaves, resulting in significant changes in its appearance during the fall. Lysine acetylation plays an important role in diverse cellular processes; however, limited knowledge is available regarding acetylations in the needle senescence of L. gmelinii. In this study, the proteomics and acetylated modification omics of two phenotypic leaves, yellow and green (senescent and non-senescent) needles, were analyzed before autumn defoliation. In total, 5022 proteins and 4469 unique acetylation sites in 2414 lysine acylated proteins were identified, and this resulted in the discovery of 1335 differentially expressed proteins (DEPs) and 605 differentially expressed acetylated proteins (DAPs) in yellow versus green needles. There are significant differences between the proteome and acetylome; only 269 proteins were found to be DEP and DAP, of which 136 proteins were consistently expressed in both the DEP and DAP, 91 proteins were upregulated, and 45 proteins were down-regulated. The DEPs participate in the metabolism of starch and sucrose, while the DAPs are involved in glycolysis and the tricarboxylic acid cycle. Among them, DEPs underwent significant changes in glycolysis and citric acid cycling. Most of the enzymes involved in glycolysis and the citrate cycle were acetylated. DAPs were down-regulated in glycolysis and up-regulated in the citrate cycle. In all, the results of this study reveal the important role of lysine acetylation in the senescence of L. gmelinii needles and provide a new perspective for understanding the molecular mechanism of leaf senescence and tree seasonal growth.

Analysis of Hyperosmotic Tolerance Mechanisms in Gracilariopsis lemaneiformis Based on Weighted Co-Expression Network Analysis

We conducted transcriptome sequencing on salt-tolerant mutants X5 and X3, and a control (Ctr) strain of Gracilariopsis lemaneiformis after treatment with artificial seawater at varying salinities (30‱, 45‱, and 60‱) for 3 weeks. Differentially expressed genes were identified and a weighted co-expression network analysis was conducted. The blue, red, and tan modules were most closely associated with salinity, while the black, cyan, light cyan, and yellow modules showed a close correlation with strain attributes. KEGG enrichment of genes from the aforementioned modules revealed that the key enrichment pathways for salinity attributes included the proteasome and carbon fixation in photosynthesis, whereas the key pathways for strain attributes consisted of lipid metabolism, oxidative phosphorylation, soluble N-ethylmaleimide-sensitive factor-activating protein receptor (SNARE) interactions in vesicular transport, and porphyrin and chlorophyll metabolism. Gene expression for the proteasome and carbon fixation in photosynthesis was higher in all strains at 60‱. In addition, gene expression in the proteasome pathway was higher in the X5-60 than Ctr-60 and X3-60. Based on the above data and relevant literature, we speculated that mutant X5 likely copes with high salt stress by upregulating genes related to lysosome and carbon fixation in photosynthesis. The proteasome may be reset to adjust the organism's proteome composition to adapt to high-salt environments, while carbon fixation may aid in maintaining material and energy metabolism for normal life activities by enhancing carbon dioxide uptake via photosynthesis. The differences between the X5-30 and Ctr-30 expression of genes involved in the synthesis of secondary metabolites, oxidative phosphorylation, and SNARE interactions in vesicular transport suggested that the X5-30 may differ from Ctr-30 in lipid metabolism, energy metabolism, and vesicular transport. Finally, among the key pathways with good correlation with salinity and strain traits, the key genes with significant correlation with salinity and strain traits were identified by correlation analysis.

Multi-omics explores the potential regulatory role of acetylation modification in flavonoid biosynthesis of Ginkgo biloba

Flavonoids are crucial medicinal active ingredients in Ginkgo biloba L. However, the effect of protein post-translational modifications on flavonoid biosynthesis remains poorly explored. Lysine acetylation, a reversible post-translational modification, plays a crucial role in metabolic regulation. This study aims to investigate the potential role of acetylation in G. biloba flavonoid biosynthesis. Through comprehensive analysis of transcriptomes, metabolomes, proteomes and acetylated proteins in different tissues, a total of 11,788 lysine acetylation sites were identified on 4324 acetylated proteins, including 89 acetylation sites on 23 proteins. Additionally, 128 types of differentially accumulated flavonoids were identified among tissues, and a dataset of differentially expressed genes related to the flavonoid biosynthesis pathway was constructed. Twelve (CHI, C3H1, ANR, DFR, CCoAOMT1, F3H1, F3H2, CCoAOMT2, C3H2, HCT, F3'5'H and FG2) acetylated proteins that might be involved in flavonoid biosynthesis were identified. Specifically, we found that the modification levels of CCoAOMT1 and F3'5'H sites correlated with the catalytic production of homoeriodictyol and dihydromyricetin, respectively. Inhibitors of lysine deacetylase (trichostatin A) impacted total flavonoid content in different tissues and increased flavonoid levels in G. biloba roots. Treatment with trichostatin A revealed that expression levels of GbF3'5'H and GbCCoAOMT1 in stems and leaves aligned with total flavonoid content variations, while in roots, expression levels of GbC3H2 and GbFG2 corresponded to total flavonoid content changes. Collectively, these findings reveal for the first time the important role of acetylation in flavonoid biosynthesis.

Proteomic Analysis of Lysine Acetylation and Succinylation to Investigate the Pathogenicity of Virulent Pseudomonas syringae pv. tomato DC3000 and Avirulent Line Pseudomonas syringae pv. tomato DC3000 avrRpm1 on Arabidopsis thaliana

Pseudomonas syringae pv. tomato DC3000 (Pst DC3000) is able to infect many economically important crops and thus causes substantial losses in the global agricultural economy. Pst DC3000 can be divided into virulent lines and avirulent lines. For instance, the pathogen effector avrRPM1 of avirulent line Pst-avrRpm1 (Pst DC3000 avrRpm1) can be recognized and detoxified by the plant. To further compare the pathogenicity mechanisms of virulent and avirulent Pst DC3000, a comprehensive analysis of the acetylome and succinylome in Arabidopsis thaliana was conducted following infection with virulent line Pst DC3000 and avirulent line Pst-avrRpm1. In this study, a total of 1625 acetylated proteins encompassing 3423 distinct acetylation sites were successfully identified. Additionally, 229 succinylated proteins with 527 unique succinylation sites were detected. A comparison of these modification profiles between plants infected with Pst DC3000 and Pst-avrRpm1 revealed significant differences. Specifically, modification sites demonstrated inconsistencies, with a variance of up to 10% compared to the control group. Moreover, lysine acetylation (Kac) and lysine succinylation (Ksu) displayed distinct preferences in their modification patterns. Lysine acetylation is observed to exhibit a tendency towards up-regulation in Arabidopsis infected with Pst-avrRpm1. Conversely, the disparity in the number of Ksu up-regulated and down-regulated sites was not as pronounced. Motif enrichment analysis disclosed that acetylation modification sequences are relatively conserved, and regions rich in polar acidic/basic and non-polar hydrophobic amino acids are hotspots for acetylation modifications. Functional enrichment analysis indicated that the differentially modified proteins are primarily enriched in the photosynthesis pathway, particularly in relation to light-capturing proteins. In conclusion, this study provides an insightful profile of the lysine acetylome and succinylome in A. thaliana infected with virulent and avirulent lines of Pst DC3000. Our findings revealed the potential impact of these post-translational modifications (PTMs) on the physiological functions of the host plant during pathogen infection. This study offers valuable insights into the complex interactions between plant pathogens and their hosts, laying the groundwork for future research on disease resistance and pathogenesis mechanisms.

Deciphering the structure, function, and mechanism of lysine acetyltransferase cGNAT2 in cyanobacteria

Lysine acetylation is a conserved regulatory posttranslational protein modification that is performed by lysine acetyltransferases (KATs). By catalyzing the transfer of acetyl groups to substrate proteins, KATs play critical regulatory roles in all domains of life; however, no KATs have yet been identified in cyanobacteria. Here, we tested all predicted KATs in the cyanobacterium Synechococcus sp. PCC 7002 (Syn7002) and demonstrated that A1596, which we named cyanobacterial Gcn5-related N-acetyltransferase (cGNAT2), can catalyze lysine acetylation in vivo and in vitro. Eight amino acid residues were identified as the key residues in the putative active site of cGNAT2, as indicated by structural simulation and site-directed mutagenesis. The loss of cGNAT2 altered both growth and photosynthetic electron transport in Syn7002. In addition, quantitative analysis of the lysine acetylome identified 548 endogenous substrates of cGNAT2 in Syn7002. We further demonstrated that cGNAT2 can acetylate NAD(P)H dehydrogenase J (NdhJ) in vivo and in vitro, with the inability to acetylate K89 residues, thus decreasing NdhJ activity and affecting both growth and electron transport in Syn7002. In summary, this study identified a KAT in cyanobacteria and revealed that cGNAT2 regulates growth and photosynthesis in Syn7002 through an acetylation-mediated mechanism.

HOOKLESS1 acetylates AUTOPHAGY-RELATED PROTEIN18a to promote autophagy during nutrient starvation in Arabidopsis

Acetylation is an important posttranslational modification (PTM) that regulates almost all core processes of autophagy in yeast and mammals. However, the role of protein acetylation in plant autophagy and the underlying regulatory mechanisms remain unclear. Here, we show the essential role of the putative acetyltransferase HOOKLESS1 (HLS1) in acetylation of the autophagy-related protein ATG18a, a key autophagy component that regulates autophagosome formation in Arabidopsis (Arabidopsis thaliana). Loss of HLS1 function suppressed starvation-induced autophagy and increased plant susceptibility to nutrient deprivation. We discovered that HLS1 physically interacts with and directly acetylates ATG18a both in vitro and in vivo. In contrast, mutating putative active sites in HLS1 inhibited ATG18a acetylation and suppressed autophagy upon nutrient deprivation. Accordingly, overexpression of ATG18a mutant variants with lower acetylation levels inhibited the binding activity of ATG18a to PtdIns(3)P and autophagosome formation under starvation conditions. Moreover, HLS1-modulated autophagy was uncoupled from its function in hook development. Taken together, these findings shed light on a key regulator of autophagy and further elucidate the importance of PTMs in modulating autophagy in plants.

The many faces of lysine acylation in proteins: Phytohormones as unexplored substrates

Protein post-translational modification (PTM) is a ubiquitous process that occurs in most proteins. Lysine residues containing an ε-amino group are recognized as hotspots for the addition of different chemical groups. Lysine acetylation, extensively studied in histones, serves as an epigenetic hallmark capable of promoting changes in chromatin structure and availability. Acyl groups derived from molecules involved in carbohydrate and lipid metabolisms, such as lactate, succinate and hydroxybutyrate, were identified as lysine modifications of histones and other proteins. Lysine-acyltransferases do not exhibit significant substrate specificity concerning acyl donors. Furthermore, plant hormones harboring acyl groups often form conjugates with free amino acids to regulate their activity and function during plant physiological processes and responses, a process mediated by GH3 enzymes. Besides forming low-molecular weight conjugates, auxins have been shown to covalently modify proteins in bean seeds. Aside from auxins, other phytohormones with acyl groups are unexplored potential substrates for post-translational acylation of proteins. Using MS data searches, we revealed various proteins with lysine residues linked to auxin, abscisic acid, gibberellic acid, jasmonic acid, and salicylic acid. These findings raise compelling questions about the ability of plant hormones harboring carboxyl groups to serve as new candidates for protein acylation and acting in protein PTM and modulation.

A comprehensive dynamic immune acetylproteomics profiling induced by Puccinia polysora in maize

Lysine-ε-acetylation (Kac) is a reversible post-translational modification that plays important roles during plant-pathogen interactions. Some pathogens can deliver secreted effectors encoding acetyltransferases or deacetylases into host cell to directly modify acetylation of host proteins. However, the function of these acetylated host proteins in plant-pathogen defense remains to be determined. Employing high-resolution tandem mass spectrometry, we analyzed protein abundance and lysine acetylation changes in maize infected with Puccinia polysora (P. polysora) at 0 h, 12 h, 24 h, 48 h and 72 h. A total of 7412 Kac sites from 4697 proteins were identified, and 1732 Kac sites from 1006 proteins were quantified. Analyzed the features of lysine acetylation, we found that Kac is ubiquitous in cellular compartments and preferentially targets lysine residues in the -F/W/Y-X-X-K (ac)-N/S/T/P/Y/G- motif of the protein, this Kac motif contained proteins enriched in basic metabolism and defense-associated pathways during fungal infection. Further analysis of acetylproteomics data indicated that maize regulates cellular processes in response to P. polysora infection by altering Kac levels of histones and non-histones. In addition, acetylation of pathogen defense-related proteins presented converse patterns in signaling transduction, defense response, cell wall fortification, ROS scavenging, redox reaction and proteostasis. Our results provide informative resources for studying protein acetylation in plant-pathogen interactions, not only greatly extending the understanding on the roles of acetylation in vivo, but also providing a comprehensive dynamic pattern of Kac modifications in the process of plant immune response.

Proteome-wide analyses reveal diverse functions of protein acetylation and succinylation modifications in fast growing stolons of bermudagrass (Cynodon dactylon L.)

BACKGROUND

Bermudagrass (Cynodon dactylon L.) is an important warm-season turfgrass species with well-developed stolons, which lay the foundation for the fast propagation of bermudagrass plants through asexual clonal growth. However, the growth and development of bermudagrass stolons are still poorly understood at the molecular level.

RESULTS

In this study, we comprehensively analyzed the acetylation and succinylation modifications of proteins in fast-growing stolons of the bermudagrass cultivar Yangjiang. A total of 4657 lysine acetylation sites on 1914 proteins and 226 lysine succinylation sites on 128 proteins were successfully identified using liquid chromatography coupled to tandem mass spectrometry, respectively. Furthermore, 78 proteins and 81 lysine sites were found to be both acetylated and succinylated. Functional enrichment analysis revealed that acetylated proteins regulate diverse reactions of carbohydrate metabolism and protein turnover, whereas succinylated proteins mainly regulate the citrate cycle. These results partly explained the different growth disturbances of bermudagrass stolons under treatment with sodium butyrate and sodium malonate, which interfere with protein acetylation and succinylation, respectively. Moreover, 140 acetylated proteins and 42 succinylated proteins were further characterized having similarly modified orthologs in other grass species. Site-specific mutations combined with enzymatic activity assays indicated that the conserved acetylation of catalase and succinylation of malate dehydrogenase both inhibited their activities, further implying important regulatory roles of the two modifications.

CONCLUSION

In summary, our study implied that lysine acetylation and succinylation of proteins possibly play important regulatory roles in the fast growth of bermudagrass stolons. The results not only provide new insights into clonal growth of bermudagrass but also offer a rich resource for functional analyses of protein lysine acetylation and succinylation in plants.

Global proteome analyses of phosphorylation and succinylation of barley root proteins in response to phosphate starvation and recovery

Phosphate (Pi) stress is an important environmental factor that limits plant growth and development. Of various posttranslational modifications (PTMs), protein phosphorylation and succinylation are the two most important PTMs that regulate multiple biological processes in response to Pi stress. However, these PTMs have been investigated individually but their interactions with proteins in response to Pi stress remain poorly understood. In this study, to elucidate the underlying mechanisms of protein phosphorylation and succinylation in response to Pi stress, we performed a global analysis of the barley root phosphorylome and succinylome in Pi starvation and recovery stages, respectively. A total of 3,634 and 884 unique phosphorylated and succinylated proteins, respectively, corresponding to 11,538 and 2,840 phospho- and succinyl-sites, were identified; of these, 275 proteins were found to be simultaneously phosphorylated and succinylated. Gene Set Enrichment Analysis was performed with a Kyoto Encyclopedia of Genes and Genomes pathway database revealing pathways that significantly enriched in the phosphorylome and succinylome. Such pathways, were dynamically regulated by Pi starvation and recovery treatments, and could be partitioned into distinct metabolic processes. In particular, phosphorylated proteins related to purine, the mitogen-activated protein kinase (MAPK) signaling pathway, pyrimidine, and ATP-binding cassette (ABC) transporters were upregulated in both Pi deprivation and recovery stages. Succinylated proteins, significantly upregulated by both Pi starvation and recovery, were enriched in nitrogen metabolism and phenylpropanoid biosynthesis. Meanwhile, succinylated proteins that were significantly downregulated by both Pi starvation and recovery were enriched in lysine degradation and tryptophan metabolism. This highlighted the importance of these metabolic pathways in regulating Pi homeostasis. Furthermore, protein-protein interaction network analyses showed that the response of central metabolic pathways to Pi starvation and recovery was significantly modulated by phosphorylation or succinylation, both individually and together. In addition, we discovered relevant proteins involved in MAPK signaling and phenylpropanoid biosynthetic pathways existing in interactions between phosphorylated and succinylated proteins in response to Pi recovery. The current study not only provides a comprehensive analysis of phosphorylated and succinylated proteins in plant responses to Pi starvation and recovery, but also reveals detailed interactions between phosphorylated and succinylated proteins in barley roots.