- 1Division of Crop Genetics and Plant Breeding, Crop Research Institute, Prague, Czechia
- 2Biotechnologies and Environmental Analytics Platform (BEAP), Environmental Research and Innovation (ERIN) Department, Luxembourg Institute of Science and Technology (LIST), Esch-Sur-Alzette, Luxembourg
Proteins are directly involved in plant phenotypic response to ever changing environmental conditions. The ability to produce multiple mature functional proteins, i.e., proteoforms, from a single gene sequence represents an efficient tool ensuring the diversification of protein biological functions underlying the diversity of plant phenotypic responses to environmental stresses. Basically, two major kinds of proteoforms can be distinguished: protein isoforms, i.e., alterations at protein sequence level arising from posttranscriptional modifications of a single pre-mRNA by alternative splicing or editing, and protein posttranslational modifications (PTMs), i.e., enzymatically catalyzed or spontaneous modifications of certain amino acid residues resulting in altered biological functions (or loss of biological functions, such as in non-functional proteins that raised as a product of spontaneous protein modification by reactive molecular species, RMS). Modulation of protein final sequences resulting in different protein isoforms as well as modulation of chemical properties of key amino acid residues by different PTMs (such as phosphorylation, N- and O-glycosylation, methylation, acylation, S-glutathionylation, ubiquitinylation, sumoylation, and modifications by RMS), thus, represents an efficient means to ensure the flexible modulation of protein biological functions in response to ever changing environmental conditions. The aim of this review is to provide a basic overview of the structural and functional diversity of proteoforms derived from a single gene in the context of plant evolutional adaptations underlying plant responses to the variability of environmental stresses, i.e., adverse cues mobilizing plant adaptive mechanisms to diminish their harmful effects.
Introduction
Protein versatility reflects the versatility of plant phenotypes, since proteins are directly involved in biological processes shaping plant phenotypes. Plant genome sequencing programs revealed a relatively low number of protein-coding genes (e.g., around 27,000 protein-coding genes in Arabidopsis thaliana genome) with respect to the number of different functional proteins present in plant cells, tissues, or even at whole-organism level [more than 100,000 different functional protein species and around 25 billion protein molecules are estimated to be present in Arabidopsis mesophyll cell by Heinemann et al. (2020)], thus, it is evident that complexity at the levels “gene” → “transcript” → “protein” rises due to the fact that a single gene can produce multiple functional proteins via an array of posttranscriptional and posttranslational modifications (Figure 1). It is, thus, evident that plant response to environmental cues at proteome level cannot be simply derived from transcriptome data. All proteins encoded by a single gene are termed proteoforms (Smith et al., 2013).
Figure 1. Schematic overview of posttranscriptional and posttranslational mechanisms leading to multiple functional proteoforms arising from a single gene.
Basically, two major categories of functional proteins derived from a given gene can be distinguished: “protein isoforms,” which differ in their primary protein sequence, and “protein posttranslational modifications (PTMs),” which differ in modifications of side amino acid residues ranging from a few atoms (nitrosylation, carbonylation) to whole peptides (S-glutathionylation, ubiquitinylation, sumoylation). In addition, alterations in protein/enzyme conformation in response to ligand/substrate interaction can be considered another kind of proteoforms. Differences in primary protein sequences among protein isoforms lead to differential molecular weight (MW), while conjugation of charged molecules or residues leads to altered isoelectric point (pI) values in differential PTMs. Since differential protein isoforms and PTMs often acquire differential biological functions, studying different protein isoforms and PTMs is of high importance for understanding protein functions. Methods based on 2-DE and its modifications are suitable for studies on protein isoforms and PTMs because of protein separation according to pI and MW values. Moreover, some specialized 2-DE gel-based protocols on PTMs detection were developed including either specific staining such as phosphoprotein detection by ProQ Diamond stain or PTM-specific primary antibody such as immunoblot detection of carbonylated proteins modified by dinitrophenylhydrazine (DNPH) to dinitrophenylhydrazones by a specific primary antibody [Oxyblot kit; Levine et al. (1990)], biotin-switch technique (BST) for detection of S-nitrosylated proteins (Jaffrey and Snyder, 2001) or concanavalin A lectin blot for glycoprotein detection (Hashiguchi and Komatsu, 2016) and immunoprecipitation methods for S-glutathionylated proteins detection (Butturini et al., 2019). However, the dominant methods for proteoforms detection represent MS/MS-based approaches, such as TMT labeling for phosphoprotein detection. Recently, an excellent review on MS-based approaches used for PTM identifications was published by Virág et al. (2020). Since only a relatively small fraction of a given protein undergoes PTM, the modified proteins reveal low abundances and have to be enriched in the protein samples; for example, ZrO2- or TiO2-IMAC beads are used for phosphoprotein enrichment. Approaches of immune affinity chromatography using specific anti-PTM primary antibodies are also widely used for PTM enrichment. However, the differential accumulation patterns and biological functions of different proteoforms that raised from a single gene indicate an urgent need to study plant responses to environmental cues at the proteome and PTMome levels to understand the plant phenotype under stress.
Plants as sessile organisms have developed evolutionary mechanisms to adapt to adverse environmental conditions imposing stress responses. Stress-adaptive responses involve several alterations in plant water regime (dehydration stress), redox homeostasis (enhanced risk of redox stress as a consequence of alterations in plant metabolism leading to enhanced levels of RMS) and altered physio-biochemical properties of biomolecules (e.g., altered fluidity of cellular membranes under low or high temperatures). Proteomics studies found out that these alterations are underlied by differences at proteome level including not only differential protein relative accumulation, but also differential spectrum of proteoforms.
Recently, it became evident that stress-induced alterations in protein relative abundance represent just one aspect of plant proteome stress responses, and that modifications at posttranscriptional and posttranslational levels leading to differential protein isoforms and PTMs, respectively, add another layer to the complexity and fine tuning of the plant stress response (Barua et al., 2019). The impacts of differential protein abundance patterns, structural modifications, such as isoforms and PTMs, cellular localization, and interacting partners in plant responses to environmental stresses were already discussed in recent reviews (Wang and Komatsu, 2016a; Kosová et al., 2018, 2021). Basic literature information on most studied protein PTMs in crop plants under abiotic stress was reviewed by Hashiguchi and Komatsu (2016) and Mustafa and Komatsu (2021) with a major focus on protein phosphorylation, glycosylation, acetylation, and succinylation, and by Wu et al. (2016) with a major focus on phosphorylation, ubiquitination and sumoylation, and most studied redox modifications such as protein carbonylation and S-nitrosylation. Recently, the role of PTMs in plant-pathogen interactions was reviewed by Gough and Sadanandom (2021) with a focus on phosphorylation, ubiquitination, and SUMoylation. Specialized reviews on plant protein phosphorylation under both abiotic and biotic stress signaling (Rampitsch and Bykova, 2012), rice phosphoproteomics studies (Ajadi et al., 2020), sumoylation (Benlloch and Lois, 2018), and nitrosylation (Romero-Puertas et al., 2013; Feng et al., 2019) were also published. A review focused on selected posttranscriptional and posttranslational regulations of drought and heat stress responses in plants with a focus on alternative splicing, RNA-mediated silencing, and ubiquitination and sumoylation protein modifications was published by Guerra et al. (2015).
The aim of this review is to provide a basic overview of the diversity and versatility of proteoforms derived from a single gene in order to emphasize the necessity of studying differences in protein isoforms and PTMs to understand protein biological functions. Proteoforms, i.e., multiple functional proteins derived from a single gene, highlight the necessity to study plant phenotypic response at proteome level, since it cannot be simply derived from transcriptome data. Examples of the impacts of diverse posttranscriptional and posttranslational modifications on the diversity of proteoforms arised from a single gene with respect to their biological functions are provided.
Protein Isoforms
Protein isoforms can either arise from a single gene via posttranscriptional modifications such as alternative splicing or RNA editing or they can represent products of homologous genes including both orthologous and paralogous genes. Orthologous genes represent homologs of a given gene in genomes of another plant species, while paralogous genes represent homologous genes that are products of the local multiplication of a given sequence in a single genome. As an example of gene paralogs, cluster of CBF genes at Fr2 locus on the long arm of homoeologous group 5 chromosomes in Triticeae can be given (Tondelli et al., 2011). Several important crops, such as bread wheat, potato, oilseed rape, banana, and others, which are allopolyploids, contain whole sets of homeologous genomes; thus, they contain homeologous genes. Homeologs are pairs of genes that originated by speciation and are brought back together by allopolyploidization (Glover et al., 2016). In addition to sequence homeology, the array of protein isoforms can be significantly increased by the mechanism of alternative splicing, giving rise to multiple functional proteins from a single gene. Alternative splicing, thus, represents an efficient evolutionary mechanism of enhancing phenotypic variation in response to environmental cues (Shang et al., 2017). Mechanisms of isoform differentiation at the transcript level are not included in this study, since we focus on proteins and their biological functions.
Protein isoforms can reveal differential biological functions because of differential expression patterns, differential cellular localization, differential protein interaction partners, and others. Some examples of functional diversification of isoformic proteins are provided below:
Transcriptome and proteome reprogramming under stress: non-canonical histone isoforms were found to play important roles in transcriptome reprogramming under stress. Alterations in transcript expression between canonical H2A and cold-induced H2AZ histone isoforms affecting nucleosome structure and DNA wrapping underlying transcriptome remodeling under cold acclimation were observed in A. thaliana (Kumar and Wigge, 2010) as well as in winter barley crowns (Janská et al., 2011). At proteome level, Martinez-Seidel et al. (2021) reported differential ribosomal proteins paralogs underlying alterations at ribosomal level and proteome reprogramming in barley root tip meristems subjected to cold acclimation.
Differential gene copy number determining quantitative differences in gene expression and quantitative phenotypic traits: Knox et al. (2010) and Würschum et al. (2017) revealed that the quantitative differences in acquired frost tolerance achieved by winter vs. spring Triticeae genotypes are associated with a higher copy number of CBF genes encoded by frost-tolerant winter genotypes with respect to the spring ones.
Differential expression patterns and differential interaction protein partners: Hsc70 protein is known as a cochaperone involved in protein folding. Hsc70-1 and Hsc70-3 isoforms are constitutively expressed, while Hsc70-2 isoform is pathogen-inducible and forms a complex with SGT1 protein involved in R-gene mediated resistance (Noël et al., 2007).
RubisCO activase (RCA) is crucial for maintaining the high carboxylation activity of RubisCO, since it removes pentosephosphates inhibiting the formation of carbamate critical for RubisCO carboxylation activity from RubisCO active site. Wang et al. (2020) studied RCA isoforms in relatively heat-tolerant Rhododendron hainanense. It encodes 3 RCA homologous genes named RCA1, RCA2, and RCA3, and each of them can give rise to multiple isoforms via a mechanism of alternative splicing (RCA1, RCA2 and RCA3 can potentially encode 5, 3, and 2 isoforms, respectively). While RCA2 and RCA3 isoforms are constitutively expressed, RCA1 isoforms are heat-inducible (especially the RCA1-X1 isoform) and persist in plants in recovery following a heat stress treatment; thus, they seem to provide some type of epigenetic stress memory involved in heat stress acclimation; enhanced levels of RCA1 persisting in plants reveal an improvement effect on plant photosynthesis when exposed to stress (Wang et al., 2020).
Differential cellular localization: Differential functions were described for nuclear isoforms and cytoplasmic (plastid, mitochondrial) isoforms of glycolytic enzymes. The nuclear isoform of fructose-bisphosphate aldolase (FBA) is known to act as a DNA-binding protein involved in regulation of expression of its own gene as well as other genes (Ronai et al., 1992). The nuclear isoform of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) is known to act as a tRNA-binding protein involved in tRNA export (Singh and Green, 1993). It has also been reported that FBA may be involved in the integration of intra- and extracellular signals associated with growth, development, and sugar biosynthesis (Li et al., 2012), and that non-the phosphorylating isoform of GAPDH may be involved in regulation of ROS levels (Bustos et al., 2008). The nuclear isoform of enolase encoded by Los2 locus in A. thaliana functions as a transcriptional repressor of STZ/ZAT10, which is a repressor of cold-inducible CBF pathway; the nuclear isoform of enolase, thus, acts as an indirect positive regulator of CBF pathway and CBF-regulated COR gene expression (Lee et al., 2002). Tonoplast-associated isoforms of aldolase and enolase are involved in the activation of V-ATPase in salt-treated Mesembryanthemum crystallinum plants, leading to Na+ vacuolar accumulation (Barkla et al., 2009).
Differential regulatory functions: Cytoplasmic eIF5A is known as a translation initiation factor involved in the regulation of protein biosynthesis. However, nuclear eIF5A isoforms named eIF5A-1 and eIF5A-2 were reported to be involved in the regulation of cell cycle with eIF5A-1 involved in a “switch” between cell proliferation and cell death, with eIF5A-1 leading to programmed cell death (PCD) while eIF5A-2 enhancing cell proliferation (Thompson et al., 2004).
Different cofactors and cellular localizations: Upon stress, a risk of ROS production in aerobic metabolism increases because of imbalances between the kinetics of electrontransport processes (respiratory electron transport chain, photosynthetic electrontransport chain) and the following enzymatic reactions. To diminish the harmful effects of ROS overproduction, plants express a wide array of several redox metabolism-related enzymes in order to finely tune their final levels. ROS levels are, thus, controlled by different peroxiredoxin (Prx) and thioredoxin (Trx) isoforms ensuring a fine tuning of redox processes in photosynthesis and respiratory electron transport chains in inner chloroplast and mitochondrial membranes, respectively. Superoxide dismutase (SOD) isoforms represent a well-known ROS scavenging enzyme with different cofactor ions and cellular localizations: Cu/Zn-SOD-cytoplasmic; Fe-SOD-plastid; Mn-SOD-mitochondrial localization. Several isoforms of ROS-scavenging enzymes, thus, provide fine tuning and control of oxidative stress, which is ubiquitously associated with aerobic metabolism localized to semiautonomous organelles, i.e., photosynthetic and respiratory electron transport chains localized in inner chloroplast and mitochondrial membranes, respectively (Meyer et al., 2005; Dietz, 2016).
Different isoforms form different cellular structures and reveal differential protein-protein interactions: Kijima et al. (2018) studied differential actin isoforms in A. thaliana and found that the A. thaliana genome encodes eight actin isoforms, out of which three (AtACT2, AtACT7, and AtACT8) are expressed at the vegetative stage while the remaining five (AtACT1, AtACT3, AtACT4, AtACT11, and AtACT12) are expressed at reproductive stage. They also found out that two vegetative actin isoforms, AtACT2 and AtACT7, form different types of filaments, with AtACT2 forming thinner filaments while AtACT7 forming thick bundles. They also reveal a different expression upon stress and interact with different actin-binding proteins (ABPs).
Examples of functional diversification of plant protein isoforms in plant stress responses are given in Table 1.
Table 1. Examples of cellular and functional diversification of protein isoforms in response to environmental stresses.
Protein Posttranslational Modifications
Protein posttranslational modifications (PTMs) involve all modifications of side amino acid residues occurring in nascent proteins and ranging from a few atoms, such as carbonyl or nitrosyl groups to whole peptides such as S-glutathione, ubiquitin, and SUMO. Currently, over 300 different PTMs were reported, such as phosphorylation, N- and O-glycosylation; reactive molecular species (RMS)-related modifications, such as glycation, carbonylation, nitrosylation, and others; peptide modifications, such as ubiquitylation and sumoylation, as well as acylation modifications of lysine residues such as succinylation (Xu W. et al., 2017) and crotonylation (Xu Y. et al., 2017). Protein PTMs range from single-step spontaneous non-enzymatic amino acid modifications by RMS, such as protein carbonylation or nitrosylation to large PTMs resulting from highly coordinated multi-step enzymatic processes, such as protein N-glycosylation or ubiquitination. An overview of most studied PTMs, such as target amino acids and added functional groups, is given in Table 2. Since recent articles revealed differential alterations in plant total proteome and PTMome in response to a given biological process (He et al., 2020), it is necessary to study plant responses to a given cue such as environmental stresses at PTMome level. Different PTMs result in differential modulation of protein biological functions; for example, regarding E3 RING-type ubiquitin ligase PARKIN, it was reported that nitrosylation decreases while persulfidation increases its ubiquitination capacity in mammalian neurons (Vandiver et al., 2013).
Table 2. Types of protein posttranslational modifications (PTMs), such as reactive agents, target amino acids, and resulting adducts, and methods used for PTM detection.
When studying protein PTMs, it has to be kept in mind that only a fraction of a given protein undergoes a given PTM and that the given PTM could be reversible and transient, i.e., to occur for a limited time, as it is evident for proteins in signaling cascades associated with transient phosphorylation/dephosphorylation events or for some RMS-mediated PTMs, which act as a signal indicating plant protein damage under stress. Moreover, a given amino acid residue can undergo multiple PTMs depending on the kind of stress; thus, the problem of PTM crosstalk has to be studied. In addition, only a few studies dealt with the impacts of combined stress treatments, which are common in the field conditions on plant protein PTM patterns. The results of the few studies indicate that plant response to combined stress is associated with a unique PTM pattern in comparison to the effects of single stress factors (Hu et al., 2015b).
The concept of PTM “writers,” “erasers,” and “readers” was proposed by Leutert et al. (2021). PTM mechanisms consist of “writers,” i.e., enzymes catalyzing the addition of the PTM group to a target amino acid residue (kinases, ubiquitin ligases, acylases,…); “erasers,” i.e., enzymes catalyzing reversible removal of a given PTM from a target amino acid residue (phosphatases, deubiquitinases, deacetylases,…); and “readers,” i.e., all proteins interacting with the given PTM.
Accumulation of the data from tandem MS/MS identifications on protein PTMs led to the development of specific databases and prediction tools for most studied PTMs. Plant protein PTMs are specifically mapped by the Plant PTM Viewer1, an integrative PTM resource database providing data on 19 PTM types (Willems et al., 2019), which currently contains PTMs data from five plant species, A. thaliana (165,279 PTMs in 55,988 proteins), Chlamydomonas reinhardtii (18,070 PTMs in 5,951 proteins), Oryza sativa ssp. japonica (56,606 PTMs in 19,500 proteins), Triticum aestivum (53,580 PTMs in 25,150 proteins), and Zea mays (143,869 PTMs in 37,099 proteins) (accessed 3rd November, 2021). MultiPTMs prediction tools, such as PTMscape (Li et al., 2018) and PTM-ssMP web server (Liu et al., 2018) utilizing information on site-specific profiles and enabling the identification of potential target motifs for multiple kinds of PTMs such as phosphorylation, acetylation, methylation, O-glycosylation, ubiquitination, and sumoylation, are available. Moreover, prediction tools for the most studied PTMs, such as PhosphoSitePlus (PSP) (Hornbeck et al., 2015), PHOSIDA (Gnad et al., 2011) for phosphoprotein detection, and UbiSite, as a web tool for ubiquitination prediction (Huang et al., 2016), etc., are available.
A brief overview of major high-throughput proteomic studies focused on stress-responsive PTMomes is provided in Table 3. In addition to the PTMs described below, other PTMs, such as protein C-terminal amidation or monomeric and polymeric ADP-ribosylation (MARylation and PARylation, linear or branched), are emerging to play an important role in plant innate immunity and plant-bacterial pathogen interactions. However, elucidation of their roles in plant responses to bacterial pathogens and biotic stress deserves further study (Feng et al., 2016).
Table 3. Overview of high-throughput proteomic studies focused on PTMome analysis in stress-treated plants.
Phosphorylation
Phosphorylation is the addition of a phosphate group to side amino acid chains with free hydroxyl group, i.e., serine (S), threonine (T, more frequent), and tyrosine (Y, less frequent), leading to the introduction of a negative charge and decrease in pI value of the protein. In addition to amino acids with a free hydroxyl group, phosphorylation of histidine residues resulting in phosphohistidine was also reported mostly on prokaryotes and lower eukaryotes such as yeast. However, phosphohistidine reveals stability lower than that of amino acids with phosphorylated hydroxyl groups; it, thus, represents a transient intermediate in some reactions.
Phosphorylation represents the most studied PTM at the plant proteomic level. currently, phosphoproteomics studies were published on salt-treated Arabidopsis thaliana membrane fraction (Hsu et al., 2009), osmotic stress-treated A. thaliana (Xue et al., 2013), salt-treated rice (Chitteti and Peng, 2007), drought-treated rice (Ke et al., 2009), Cd-treated rice (Zhong et al., 2017; Fang et al., 2019), chickpea under dehydration stress (Subba et al., 2013; Barua et al., 2019), wheat embryo and endosperm exposed to salinity (Luo et al., 2019), wheat seedling leaves and grains exposed to osmotic stress (Zhang et al., 2014a,b), wheat flag leaves exposed to dehydration (Luo et al., 2018), wheat embryo and endosperm under salinity (Luo et al., 2019), wheat spikes infected by Fusarium graminearum (Ding et al., 2016), einkorn wheat leaves exposed to salinity (Lv et al., 2016), barley root under salinity (Horie et al., 2011), Brachypodium distachyon exposed to drought (Yuan et al., 2016), soybean root tips exposed to flooding (Yin and Komatsu, 2015), maize leaves subjected to drought, heat, and their combination (Hu et al., 2015b), maize leaves under osmotic stress (Hu et al., 2015a), maize and drought (Bonhomme et al., 2012), common bean exposed to osmotic stress (Yang et al., 2013), sugar beet exposed to salinity (Yu et al., 2016), okra seedlings exposed to salinity (Yu et al., 2019), Picea wilsonii pollen subjected to sucrose deficit (Chen et al., 2012), and others as listed in Table 3.
Protein phosphorylation is catalyzed by protein kinases while protein dephosphorylation is catalyzed by protein phosphatases. Protein phosphorylation rates are quite low and reach only 1–2% of the entire protein abundance, especially when regarding transient protein phosphorylation in signaling pathways. Moreover, there are also differences in phosphorylation rates among different amino acid residues, indicating that some residues are quantitatively more phosphorylated while others, such as tyrosine residues, may be transiently phosphorylated only by up to 0.5%. Although tyrosine-phosphorylated proteins exist in relatively low abundance, which is 4.2% in Arabidopsis, 2.9% in rice, and 1.3% in Medicago, Y- phosphorylation plays an important role in plant responses to a changing environment (Hashiguchi and Komatsu, 2016). Due to the low phosphorylation level of target proteins, a phosphoprotein fraction has to be enriched prior to further analyses, such as using Ti4+ IMAC beads. An overview of phosphoprotein enrichment methods is provided by Fíla and Honys (2012) and Silva-Sanchez et al. (2015).
Phosphorylation was most studied in the context of protein signaling. Multiple-step protein phosphorylation cascades, such as the MAPKKK-MAPKK-MAPK kinase cascade, enable initial signal amplification during their transfer in plasmalemma-to-nucleus direction (Yamaguchi-Shinozaki and Shinozaki, 2006). In addition, calcium-dependent protein kinases are also involved in stress signaling and stress response induction. Cieśla et al. (2016) reported that Arabidopsis plants overexpressing barley calcium-dependent kinase HvCPK2a with autophosphorylation sites displayed enhanced drought sensitivity, indicating that HvCPK2a is a negative regulator of the drought stress response.
Phosphoproteomic studies revealed that phosphorylation does not only play an important role in signaling and signal transduction but also in DNA replication and transcription, stress and defense response to pathogens as well as in the regulation of ROS scavenging and cell metabolism. For example, Fang et al. (2019) reported differential phosphorylation of replication-related proteins such as DNA replication/repair enzymes, DEAD-box ATP-dependent RNA helicase 37, and transcription factors (TFs) such as MYB, zinc-finger, and F-box TFs in Cd-tolerant rice line with respect to a sensitive one. Thus, it seems that differential phosphorylation may alleviate adverse Cd impacts on replication and transcription. Luo et al. (2018) found differential phosphorylation of photosynthesis (LHCII, RubisCO activase) and energy metabolism-related (FBA, GAPDH, pyruvate decarboxylase complex components) proteins in wheat flag leaf under drought. Similarly, differential phosphorylation in photosynthesis antenna proteins and proteins involved in RNA degradation was found in salt-treated okra seedlings (Yu et al., 2019). Ding et al. (2016) studied wheat phosphoproteome as affected by F. graminearum infection and found alterations in phosphorylation of cinnamoyl-CoA reductase, isochorismate synthase, cytochrome P450, peroxidase 8, and glutathione reductase (GR) involved in lignin and phytoalexin biosynthesis as well as in ROS scavenging. Similarly, Zhang et al. (2014a) identified significant differences at phosphorylation level in proteins involved in signaling and signal transduction, RNA transcription and processing, and stress/defense pathways in two wheat cultivars exposed to osmotic stress due to PEG-6000 treatment. Phosphoproteins involved in ABA signaling include SnRK2 kinases and PP2C phosphatases, and in Ca2+ signaling include components of phosphatidylinositol-4,5-diphosphate pathway and CDPK, calmodulin2-2 and calcineurin B-like protein-interacting protein kinases (CIPKs), as well as components of MAPK signaling cascades. Their activities are often regulated by the phosphorylation status of other proteins with regulatory functions such as 14-3-3 proteins (Yu et al., 2016). Phosphorylated TFs include TFs involved in ABA-dependent pathways, such as TaABI5-1, which is phosphorylated by SnRK2 and dephosphorylated by PP6 and binds to ABRE motifs in COR/LEA gene promoters. Other ABA-dependent TFs undergoing phosphorylation are MYB TFs such as MYB1R1 and TaMYB2A. Besides ABA-dependent pathways, zinc finger CCCH domain TFs were found phosphorylated at serine residues under osmotic stress. Yin and Komatsu (2015) identified 14 differentially abundant phosphoproteins in soybean root tips exposed to flooding including nuclear-localized ABA-responsive phosphoproteins such as zinc finger/BTB domain-containing protein 47, glycine-rich protein, and rRNA-processing protein Rrp5. Similarly, Yin et al. (2014) reported a significant effect of ethylene signaling on phosphoproteome in soybean root tips under flooding, resulting in enhanced dephosphorylation of eukaryotic translation initiation factor 4G and enhanced phosphorylation of proteins involved in protein synthesis. Under high salinity, increased phosphorylation level in proteins related to stress and defense response (Cu/Zn-SOD, 2Cys Prx BAS1), carbohydrate metabolism, photosynthesis, and the transport was found in einkorn wheat (Lv et al., 2016). HvPIP2 water transport activity was modulated by phosphorylation under salinity treatment in barley roots (Horie et al., 2011). Similarly, serine phosphomotifs were identified in aquaporins PIP2,1 and PIP2, 4 in salt-treated A. thaliana (Hsu et al., 2009).
Barua et al. (2019) studied nuclear phosphoproteome in chickpea subjected to dehydration stress and identified 546 kinases such as those from MAPK kinase cascade, calcium-dependent protein kinases, SRSF protein kinase 1, and others. They also identified differential phosphorylation in E3 SUMO-protein ligase SIZ1-like isoform X1, spliceosome complex, and splicing-related proteins, indicating their altered function under stress. In addition, several regulatory proteins involved in protein degradation and regulation of flowering time and circadian clock revealed dehydration-induced dephosphorylation. A phosphomotif analysis revealed dehydration-induced enrichment in proline-directed Ser phosphorylation. Similarly, Yuan et al. (2016) identified several drought-responsive phosphoproteins involved in signal transduction (MAPK cascade, ABA-responsive SnRK2, Ca2+-responsive CDPKs), gene expression, stress and defense, photosynthesis and energy metabolism, and transmembrane transport including six phosphomotifs enriched in drought-treated B. distachyon. Hu et al. (2015b) studied maize phosphoproteome response to drought, heat, and their combination and described seven stress-enriched phosphomotifs, out of which two were common to all stresses, two were common for heat, and one was specific to the combined stress. They also found out that differential stress cues can lead to phosphorylation of differential phosphomotifs in the same proteins.
Phosphorylation as a nuclear localization signal (NLS): In dehydrins possessing a stretch of 4–10 serine residues called S-segment, phosphorylation of serine residues by SnRK2.10 kinase serves as an NLS (Close, 1997; Brini et al., 2007; Kosová et al., 2011). Enhanced phosphorylation of four dehydrin proteins without change in their relative abundance was found in common bean root tips exposed to PEG-induced osmotic stress (Yang et al., 2013).
Phosphorylation as a means of epigenetic stress memory: Ding et al. (2012) reported the importance of serine 5 (Ser5P) phosphorylation level in RNA polymerase II associated with drought stress memory-related genes in A. thaliana. The level of Ser5P in RNA polymerase II remained the same upon recovery as under drought stress in RNA polymerase II associated with drought stress memory-related COR/LEA genes RD29B and RAB18, revealing altered expression pattern under repeated stress when compared to the first stress treatment. Enhanced level of Ser5P in RNA polymerase II associated with drought stress memory-related genes persisting upon stress recovery indicates a stalled state of RNA polymerase II, i.e., its preparedness to transcription initiation upon repeating stress events that persist upon stress recovery when target transcript levels decreased. Differential phosphorylation of histone deacetylases was reported in dehydration-treated maize leaves (Bonhomme et al., 2012).
Glycosylation
Glycosylation represents the most abundant protein PTM in plants. It represents a highly spatially and temporarily coordinated enzymatic process occurring in the endoplasmic reticulum and Golgi apparatus (proteins belonging to the secretory pathway). Two kinds of protein glycosylation can be distinguished: N-glycosylation, where a relatively large oligosaccharide core is bound to an Asn residue, which is then trimmed and modified to a final form, and O-glycosylation when mostly just a single monosaccharide is linked to target Ser or Thr residues.
N-glycosylation-target amino acids represent asparagine (Asn) residues present in Asn-X-Ser/Thr motif, where X represents any amino acid except proline. In the first step, the basic tetradecyl glycan precursor GlcNAc2Man9Glc3 (Glc for glucose, Man for mannose, and GlcNAc for N-acetylglucosamine) is transferred to target asparagine residue by oligosaccharyl transferase complex (OST) in endoplasmic reticulum (ER) and the bound oligosaccharidic core is further trimmed and modified in Golgi apparatus by class I α mannosidases MNS1 and MNS2, which cleave α-1,2-mannosyl residues to generate the substrate for CGL1/GnT1, which catalyzes a GlcNAc addition to remove two additional mannose residues, to add another one GlcNAc, xylose, and fucose residues to form complex N-glycan structure (Strasser et al., 2004). Kang et al. (2008) reported a crucial role of N-glycoprotein maturation, such as complex glycan 1 (cgl1), for plasma membrane glycoproteins involved in cellulose biosynthesis and cell wall formation in salt-treated A. thaliana. Recently, a role of protein N-glycosylation in response to salinity stress was studied in A. Thaliana wild-type and two mutants, mns1mns2 and cgl1, defective in N-glycan maturation leading to the identification of salt-responsive N-glycoproteins such as class III peroxidases PRX32 and PRX34, which appeared to be involved in root growth and development under salinity stress (Liu et al., 2021).
N-glycoproteins were also studied using a concanavalin A lectin blot technique in cold-treated rice leaf sheaths (Komatsu et al., 2009) and in flooded soybean roots (Mustafa and Komatsu, 2014). In cold-treated rice leaf sheaths, 22 glycoproteins were determined by the lectin blot, out of which 12 revealed significant cold response; cold-responsive N-glycoproteins included calreticulin involved in Ca2+ signaling and several proteins involved in energy metabolism (mitochondrial F1-ATPase, 6-phosphogluconate dehydrogenase, NADP-dependent malic enzyme, emolase, UDP-glucose pyrophosphorylase) and protein folding (CPN60α, HSP90; Komatsu et al., 2009). In flooded soybean root tips, the accumulation of N-glycoproteins related to protein degradation, cell wall, and glycolysis increased, while glycosyl hydrolases, peroxidases, and other proteins related to protein glycosylation pathway, stress and defense response, and protein degradation decreased upon flooding with respect to control plants. In contrast, glycoproteins involved in glycolysis were activated (Mustafa and Komatsu, 2014). However, generally, an overall decrease in glycoprotein accumulation was found in flooded soybean root tips because of a decrease in ER-bound CNX/CRT lectin complex involved in N-glycoprotein processing (Wang and Komatsu, 2016b).
Probably the most studied protein O-glycosylation is O-GlcNAcylation, i.e., reversible addition and removal of single O-linked-β-N-acetylglucosamine (O-GlcNAc) activated as UDP-GlcNAc to target Ser or Thr residues mediated by O-GlcNAc transferases (OGTs) including two classes, SPINDLY (SPY) and SECRET AGENT (SEC), and O-GlcNAcases (OGAs).
Lectins represent an important protein group interacting with glycoproteins and, thus regulating their biological functions. One of the molecular mechanisms regulating vernalization duration in winter-type Triticeae lies in O-GlcNAc modification of Thr17 in a glycine-rich RNA-binding protein TaGRP2, which binds to VRN1 pre-mRNA first intron region, resulting in VRN1 transcript repression. During vernalization, the level of O-GlcNAc in TaGRP2 increases, leading to enhanced interaction with Jacalin lectin VER2, thus, leading to enhanced TaGRP dissociation from VRN1 pre-mRNA first intron, resulting in enhanced VRN1 transcript accumulation. O-glycosylation modification, thus, regulates VER2-TaGRP2 protein-protein interaction, resulting in altered interaction of VER2 with VRN1 first intron region and regulation of the major vernalization gene VRN1 expression (Xiao et al., 2014). O-glycosylation has also been reported for some dehydrins in cold-treated blueberry (Levi et al., 1999) and pistachio (Yakubov et al., 2005), although no specific functions for glycosylated forms of these proteins have been suggested.
Methylation
Histone methylation in arginine and lysine ε-amino groups via S-adenosylmethionine (SAM) as a universal methylating agent represents an epigenetic modification, resulting in either repression or activation of target gene expression. Methylations resulting in target gene repression include H3K9me2 (histone H3 lysine 9 dimethylation), H3K27me3 (histone H3 lysine 27 trimethylation), and H4R3me2 (histone H4 arginine 3 symmetrical dimethylation), while methylations leading to target gene transcriptional activation include H3K4me2 (histone 3 lysine 4 dimethylation) and H3K4me3 (histone H3 lysine 4 trimethylation) modifications. Histone 3 lysine K9 dimethylation (H3K9me2) followed by binding of polycomb-group (PcG) complexes results in repression of a target gene, such as FLOWERING LOCUS C (FLC), a major flowering repressor in A. thaliana, following vernalization fulfillment, thus enabling activation of flowering-inducing pathway genes SOC1 and FT (Sung and Amasino, 2005). In contrast, vernalization fulfillment in Triticeae cereals leads to a decrease in H3K27me3, while an increase in H3K4me3 at the 5′end of intron 1 region results in activation of the Triticeae VRN1 gene (Oliver et al., 2009). Similarly, the H3K4me3 of “trainable genes” revealing enhanced expression under repeated drought stress in comparison to the first drought treatment was also described as a means of drought stress memory in A. thaliana (Ding et al., 2012).
Acylation
Short Carboxylic Acids
Protein acylation, i.e., the addition of a relatively short hydrophilic negatively charged acyl group to a positively charged lysine residue, results in significant modifications of target lysine chemical properties. Lysine acylation/deacylation is enzymatically catalyzed by lysine acyltransferases (KATs)/sirtuin-class lysine deacylases (KDACs), respectively, and includes lysine acetylation, propionylation, butyrylation, β-hydroxybutyrylation, malonylation, glutarylation, succinylation, and crotonylation, which significantly modulate cellular metabolism, namely photosynthesis (Hirschey and Zhao, 2015).
Histone acylation of lysine ε-amino group represents an epigenetic modification resulting in the activation of target gene expression. Histone lysine acylation includes acetylation, propionylation, butyrylation, β-hydroxybutyrylation, malonylation, succinylation, and crotonylation employing acylCoA as an acylating agent. Crotonylation could be targeted to all histones, H1, H2A, H2B, H3, and H4. Lysine acetylation in histone H3 results in the activation gene expression such as in the case of the major vernalization gene VRN1 inducing flowering in Triticeae (Oliver et al., 2009).
Protein succinylation involves targeted succinylation of lysine residues in conserved motifs. Recently, Wang et al. (2021) studied barley succinylome in response to short-term phosphate (Pi) starvation and the following recovery using an immunoblot approach with a specific antibody raised against succinylated lysine. They identified 2,840 succinylation sites across 884 proteins, of which 11 conserved succinylation-related protein motifs were determined. Pi starvation enriched protein succinylation of ribosomal proteins (40S ribosomal protein S6, 60S ribosomal proteins L4-1, L6), proteins involved in glycolysis, protein translation (eukaryotic translation initiation factor 5A-1 EIF5A, elongation factors 1-α and Ts), and RNA degradation, while Pi recovery enriched lysine succinylome in TCA cycle, glycolysis (GAPDH), and oxidative phosphorylation (ATP synthase subunits β and δ of FoF1-ATPase) pathways.
Long Carboxylic Acids (Fatty Acids)
Addition of long hydrophobic fatty acid chains, such as palmitic acid (C16:0) or myristic acid (C14:0) (palmitoylation, myristoylation), to membrane proteins provides a hydrophobic “anchor” attaching the protein to hydrophobic lipid bilayer. N-myristoylation of lysine residues in target proteins affects protein-membrane interactions (McLaughlin and Aderem, 1995), thus modulating the function of membrane-associated proteins involved in stress signaling, such as SOS3 function (Ishitani et al., 2000) and SnRK1 kinase pathway (Pierre et al., 2007). Moreover, N-myristoylation also modulates the function of proteins involved in ubiquitin-dependent protein degradation, namely, F-box proteins and parts of 26S proteasome as well as MYB transcription factors (Boisson et al., 2003; Podell and Gribskov, 2004).
Peptide Conjugation
Protein S-Glutathionylation
S-glutathione is a three amino acid peptide (γ-glutamyl-cysteinyl-glycine) ubiquitously found in cell cytosols where it acts as a cellular redox state regulator because of reversible GSH ↔ GSSG transitions. Recently, it was found that glutathione conjugation to cysteine residues in cellular proteins acts as a protection shield against ROS, thus preventing protein oxidative modifications (Diaz-Vivancos et al., 2015). Martinez-Seidel et al. (2021) found enhanced levels of enzymes involved in glutathione biosynthesis (hydroxyacylglutathione hydrolase-like) and S-glutathione conjugation (five members of glutathione S-transferase family) in barley root tips subjected to cold acclimation, indicating a protective role of protein S-glutathionylation against oxidative stress associated with most environmental stresses.
Protein Ubiquitination
Ubiquitin is a small peptide (8.6 kDa, 76 amino acids) attached to ε-amino group in lysine residues of a target protein by enzyme-regulated series of reactions called ubiquitination. Ubiquitination is a three-step process that includes ubiquitin activation, conjugation, and ligation of ubiquitin to target lysine residues via a diglycine motif in ubiquitin C-terminus catalyzed by ubiquitin activase (E1), ubiquitin conjugase (E2), and ubiquitin ligase (E3) enzymes, respectively. Ubiquitin also underlies automodification in each of its seven lysine residues (K6, K11, K27, K29, K33, K48, and K63), out of which K48-linked polyubiquitination is the predominant type involved in protein targeting to proteasome-mediated degradation (Hershko et al., 2000). Ubiquitination as a mark of 26S proteasome-targeted protein degradation plays an important role in plant stress acclimation, since this process is associated with an enhanced rate of both protein biosynthesis and protein degradation (Sharma et al., 2016). Bueso et al. (2014) reported an interaction of RING-type E3 ubiquitin ligase RSL1 with PYL4 and PYR1 ABA receptors resulting in their targeted degradation, thus implying that the E3 ubiquitin ligase affects ABA signaling in A. thaliana. In A. thaliana, RING finger E3 ubiquitin ligases AtAIRP4 and SDIR1 were reported to positively modulate salt stress response via ABA signaling (Yang et al., 2015; Zhang et al., 2015).
Protein Sumoylation
SUMO (a small ubiquitin-like modifier; 12 kDa, 100 amino acids) are small proteins structurally similar to ubiquitin which are covalently attached via their C-terminal diglycine motif to lysine ε-amino group in the target proteins via an enzymatic cascade of SUMO-conjugating enzymes (SCE) analogous to that involved in ubiquitination, i.e., SUMO-activating enzyme (E1, SAE), SUMO-conjugating enzyme (E2, SCE) nad SUMO-ligase (E3). However, sumoylation does not lead to target protein degradation (Benlloch and Lois, 2018). Ghimire et al. (2020) studied alterations in potato (Solanum tuberosum) seven StSUMO genes, and nine StSCE enzymes in response to salt and PEG treatments at transcript levels by qRT-PCR and found an upregulation in StSCE1/5/6/7 under both salt and PEG treatments, while StSUMO1/2/4 were upregulated only under salt stress, whereas StSUMO2/4 were downregulated under PEG treatment.
Protein Posttranslational Modifications by Reactive Molecular Species
Environmental stress leads to an enhanced imbalance in key processes of plant energy metabolism, especially photosynthesis and aerobic respiration, resulting in enhanced production of RMS. According to reactive atoms, major types of RMS include reactive oxygen species (ROS), reactive nitrogen species (RNS), reactive carbonyl species (RCS), reactive sulfur species (RSS), reactive halogen species (RXS), and reactive metal species. RMS-mediated PTMs, thus, represent products of spontaneous reactions, which, in case of amino acids in protein active sites, can potentially lead to loss of target protein biological activity. RMS-mediated PTMs have recently been reviewed by Mock and Dietz (2016). ROS, RNS, and RSS in legumes subjected to abiotic stresses were discussed by Matamoros and Becana (2021).
Reactive Oxygen Species
Reactive oxygen species (ROS) include superoxide anion radical, hydroxyl radical, hydrogen peroxide, singlet oxygen. Major ROS in plant cells arise from imbalance at acceptor sides of photosynthetic electron transport chain and respiratory electron transport located in inner membranes of chloroplast and mitochondria, respectively. Photosynthetic electron transport chain reveals significant alterations in absorbed light energy leading to electron excitation and the following transport in response to sunlight whose intensity can significantly change within fractions of a second. Excess electrons at the acceptor side of the photosynthetic electron transport chain react with oxygen (the so-called Mehler reaction/Fenton reaction), resulting in superoxide anion formation. ROS also react with a wide array of target amino acids such as sulfur-containing cysteinyl and methionyl residues (Cys, Met), His, Trp, and Tyr, resulting in the formation of peroxide or dioxetane intermediates. Oxidized Tyr and Trp convert to indole derivatives. Histidyl endoperoxides can form His-His or His-Lys crosslinks. Peroxidized Met residues convert to methionine sulfoxide (Mock and Dietz, 2016).
Reactive Carbonyls: Protein Glycation
Protein glycation represents a spontaneous non-enzymatic modification of side free amino groups in target amino acids lysine, arginine, proline, and threonine by active carbonyl groups found in some metabolic intermediates (active dicarbonyls such as glyoxal and methylglyoxal, MG) and reducing sugars, such as glucose, through the so-called Maillard reaction. The reaction is spontaneous and can lead to loss of protein biological activity, since target lysine and arginine residues prone to glycation are often found in protein active sites and domains underlying protein biological activity. Moreover, the risk of protein glycation increases under stress because of enhanced levels of active dicarbonyls glyoxal, MG, and 3-deoxyglucosone as major glycating agents. In contrast, sucrose, as a main sugar transport form in plants, is a non-reducing sugar (no active carbonyl group); thus, it does not reveal a glycating activity (Suzuki et al., 2010; Rabbani et al., 2020).
Major glycation adducts represent Nε-(1-deoxy-D-fructose-1-yl) lysine (FL), MG-derived hydroimidazolone (MG-H1; arginine-derived adduct from MG), Nε-carboxymethyl-lysine (CML), glucosepane (formed by degradation of FL residues with a proximate arginine residue), and pentosidine [for more details, see review by Rabbani et al. (2020)].
Detection of carbonylated proteins in 2-DE gels is usually performed using an OxyBlot kit, which lies in the derivatization of carbonylated amino acids by dinitrophenylhydrazine (DNPH) to form dinitrophenylhydrazones, which can be detected by a specific primary antibody against dinitrophenylhydrazones. Enhanced protein carbonylation detected by DNPH immunoblots was found in cadmium-treated pea (Romero-Puertas et al., 2002).
Protein carbonylation belongs to the earliest studied plant protein PTMs. In A. thaliana, studies on carbonylated proteins during the vegetative stage (Johansson et al., 2004) and during seed germination (Job et al., 2005) using an Oxyblot kit for carbonylated protein detection and led to the identification of a wide range of carbonylated proteins, such as HSP70, ATP synthase subunits, and RubisCO LSU, and publication of several glycolytic enzymes. Most stress factors (except flooding causing hypoxia) enhance ROS levels, resulting in enhanced protein carbonylation. Pandey and Dubey (2021) studied the combined effect of PEG-induced dehydration and aluminum toxicity on protein carbonylation in rice seedlings belonging to two contrasting cultivars, Malvyia-36 (sensitive) and Vandana (tolerant), and found differential carbonylation in several photosynthesis-related proteins such as RubisCO LSU, PS I reaction center subunit II and IV proteins, and OEE3 protein, as well as in calmodulin involved in Ca2+ signaling.
Protein Nitrosylation by Reactive Nitrogen Species
Major reactive nitrogen species (RNS) include nitric oxide (NO), nitrogen dioxide (NOO.), and peroxynitrite (ONOO–), which are generated by the reaction of NO with superoxide radicals (O2–) and oxygen, respectively. Nitric oxide (NO) is synthesized by cytosolic nitrate reductase as a side product and reacts with thiol groups in Cys residues, resulting in the formation of S-nitrosothiols, while peroxynitrite and nitrogen dioxide are major NO donors for Tyr nitration. The major buffer for NO in plant cells is glutathione (GSH), resulting in the formation of S-nitrosoglutathione (GSNO). Recently, the molecular mechanisms involved in protein S-nitrosylation and current knowledge of the biological roles of S-nitrosoproteins in plants were reviewed by Feng et al. (2019).
Abat and Deswal (2009) studied NO-induced S-nitrosoproteome in Brassica juncea treated with various environmental stresses such as cold, drought, heat, and salinity, and reported S-nitrosylation (SNO)-related inactivation of RubisCO carboxylase activity under cold stress. RubisCO depletion led to the identification of S-nitrosylation in further proteins such as thioredoxin, SOD, GAPDH, FBA, and sedoheptulose bisphosphatase, indicating an important role of S-nitrosylation in the regulation of Calvin cycle enzymes in cold-treated B. juncea seedlings (Sehrawat et al., 2013). In addition, further study on S-nitrosoproteins in cold-treated B. juncea seedlings (Sehrawat and Deswal, 2014) identified novel targets such as storage proteins (12S globulin cruciferin), stress- and defense-related proteins (lactoylglutathione lyase, germin-like protein), and redox metabolism-related proteins (DHAR, Cu/Zn-SOD). Camejo et al. (2013) reported enhanced S-nitrosylation of several mitochondrial proteins in salt-treated pea, namely ATP synthase β subunit, Hsp 90, and peroxiredoxin IIF (PRXIIF), leading to a decrease in PRXIIF enzymatic activity.
Tanou et al. (2009, 2012) studied the effects of hydrogen peroxide (H2O2) and nitric oxide (NO) as an ROS and an RNS, respectively, on citrus (Citrus aurantium) plants exposed to salinity and identified 40 carbonylated and 49 S-nitrosylated proteins, indicating an overlap between H2O2 and NO-induced signaling pathways. Ziogas et al. (2013) studied the impacts of six abiotic stresses namely, cold, heat, drought, salinity, continuous light, and darkness on citrus plants with respect to their nitrosoproteome and RNS levels and found differential impacts of the different stresses on RNS and related enzymes activities. Cold led to enhanced S-nitrosoglutathione reductase (GSNOR) and nitrate reductase (NR) gene expression and enzymatic activity. Peroxynitrite scavenging activity was elicited by continuous light, darkness, or drought but it was suppressed by salinity. Nitration activity was enhanced by salinity and suppressed by continuous light or darkness. Protein S-nitrosylation levels were enhanced by heat, cold, and drought but were suppressed by darkness and salinity.
The impacts of heat stress on tyrosine nitration (Y-nitration) were studied in sunflower by Chaki et al. (2011, 2013) who found that tyrosine nitration leads to the inhibition of ferredoxin-NADP reductase and β-carbonic anhydrase (β-CA) enzymatic activities, respectively, implicating adverse impacts of heat on photosynthesis via enhanced RNS.
Protein Modification by Reactive Sulfur Species
Major reactive sulfur species (RSS) include thiyl radicals, sulfenic acid, and disulfide monooxide.
Protein persulfidation, also known as protein sulfhydration, represents a modification of thiol groups (-SH) in target cysteine residues to -SSH by H2S, which is produced from cysteine in stomatal guard cells by the L-cysteine desulfhydrase 1 (DES1) enzyme, which is involved in ABA-dependent signaling pathway regulating stomatal closure. Enhanced levels of persulfidated proteins were found as being associated with autophagy in response to severe environmental stresses such as drought and nutrient starvation (e.g., nitrogen starvation). Protein persulfidation can be detected by tag-switch label-free LC-MS/MS approach. Persulfidation was reported for H2S-dependent PTM of SnRK2.6 kinase in ABA-dependent signaling pathway and DES1 itself. Persulfidation was reported for several proteins involved in ABA-dependent signaling pathways such as ABA receptors PYR/PYL and their downstream elements such as SnRK2.1 protein kinase and PP2C-7 protein phosphatase as well as proteins involved in jasmonic acid (JA) and brassinosteroid (BR) signaling, BRI1 or BSK1. Several proteins involved in signaling related to autophagy, such as SNF1-related protein kinase catalytic subunits alpha KIN10 and KIN11 and receptor-like kinases FERONIA and THESEUS1, were found persulfidated under N starvation in A. thaliana (Jurado-Flores et al., 2021). Protein persulfidation plays a role in the processes of autophagy, ubiquitin-dependent protein degradation, and ABA-dependent signaling.
Non-Covalent Protein Modifications
In addition to covalent PTMs, non-covalent protein modifications such as protein/enzyme cofactors and/or substrates resulting in different protein/enzyme conformations represent a further level leading to novel functional proteoforms. Proteins/enzymes with multiple cofactors/substrates, probably the most studied plant enzyme with multiple substrates and multiple catalytic functions, is RubisCO, i.e., ribulose-1,5-bisphosphate carboxylase/oxygenase, which can bind either CO2 or O2 as competing substrate and, thus, catalyze either carboxylation or oxygenation of the substrate ribulose-1,5-bisphosphate. RubisCO needs RubisCO activase (RCA) to activate the RubisCO active site for CO2 binding.
Enzymes can modify their affinity to a substrate based on ambient substrate concentration. In oligomeric enzymes with multiple catalytic sites, allosteric effect occurs, i.e., transition of the catalytic site from a low-affinity state to a high-affinity state is transmitted from one subunit to the other. Allosteric regulation, thus, means the spatial effect of the subunit on the others in oligomeric enzymes. This type of regulation of the affinity of the binding site is not possible in monomeric enzymes. However, it was also found that in monomeric enzymes, such as glucokinase, increasing glucose concentration as a substrate leads to an increased shift of the enzyme from low-affinity to high-affinity state analogously to the allosteric regulation in oligomeric enzymes (Whittington et al., 2015). As substrate concentration increases, the probability of encountering other substrate molecules after turnover also increases, thus preventing the enzyme from relaxing back to the low-affinity state. This time-dependent shift in the ratio of low-affinity to high-affinity enzyme proteoforms with respect to increased substrate concentration is called allokairy (Hilser et al., 2015). Mnemonical enzymes are monomeric enzymes that can retain high-affinity conformation after product release when substrate levels increase in an ambient environment (Ricard et al., 1974). The shifts between enzyme low-affinity and high-affinity states, such as allokairy in monomeric human glucokinase, play an important role in disease development (Whittington et al., 2015). However, to the best of our knowledge, no analogous research has been published for plant enzymes.
Conclusion
Proteins are directly involved in plant responses to environmental cues. Mechanisms of posttranscriptional and posttranslational modifications of protein primary sequence and encoded amino acid residues, respectively, resulting in multiple functional proteins (proteoforms) derived from a single gene. Mechanisms of posttranscriptional and posttranslational modifications increase phenotypic variability with respect to genetic information encoded by the genome. Duplications of local genomic regions, as well as whole genome duplications, and the presence of multiple genomes in polyploid species lead to paralogous and orthologous genes, respectively, which can undergo functional diversifications represented by differential expression patterns. From a biochemical point of view, the twenty gene-encoded amino acids represent a very diverse set of chemical species such as non-polar (hydrophobic) amino acids with aliphatic (alanine, leucine, isoleucine, valine) or aromatic (phenylalanine) side chains as well as both acidic (aspartic and glutamic acid) and basic (lysine, arginine, asparagine, glutamine) groups. The coded amino acids can, thus, undergo several modifications resulting in altered properties. Some PTMs represent a result of a series of spatially and temporally coordinated enzymatic reactions (phosphorylation, N-glycosylation, succinylation, crotonylation, ubiquitination, sumoylation), while others occur spontaneously (PTMs arising from reactions with RMS such as oxidation, glycation, and S-nitrosylation). The level of spontaneous PTMs rises under environmental stress because of imbalances in cellular metabolism, leading to enhanced formation of reactive molecular species (RMS). RMS-mediated PTMs can result in loss of target protein biological activity, but they also serve as a signal-inducing stress response (Chaki et al., 2011, 2013; Camejo et al., 2013; Mock and Dietz, 2016).
It has to be summarized that PTMs represent a modification of only a relatively small fraction of the given protein; moreover, the modification is often reversible and/or transient, and the resulting PTM, thus, acts as a transient signal of the given stress cue, which is the case of phosphorylation and RMS-induced PTMs. Several proteins can undergo multiple PTMs under the given stress treatment, e.g., calreticulin was reported to undergo both phosphorylation and N-glycosylation in cold-treated rice leaf sheaths (Komatsu et al., 2009). Moreover, the same protein can undergo the same kind of PTM in different target motifs under different stress treatments; for example, differential phosphomotifs were identified in the same phosphoproteins in maize phosphoproteome under drought, heat, and their combination (Hu et al., 2015b). All these factors determine the protein’s biological activity.
The major roles of differential protein isoforms and PTMs in plant cells subjected to stress could be summarized as follows (Table 4):
Table 4. Overview on major protein biological functions affected by different protein isoforms or PTMs in response to environmental stress conditions.
Signaling and signal transduction: Stress-induced signal is transduced and amplified via reversible phosphorylation in protein kinase cascade. The reversibility of phosphorylation enables efficient regulation of signal spread and switching off of it when target cellular mechanisms are activated (Yamaguchi-Shinozaki and Shinozaki, 2006).
Transcriptome/proteome reprogramming under stress: Expression of differential histone isoforms and PTMs underlying modular reprogramming of plant transcriptome under stress enables expression of whole arrays of newly induced transcripts (Kumar and Wigge, 2010; Janská et al., 2011); ribosomal isoforms underlying proteome reprogramming under stress (Martinez-Seidel et al., 2021).
PTMs of histones and transcription-associated proteins determine the preparedness of stress memory-related genes to prompt transcription under repeated stress events resulting in enhanced stress tolerance. Stress memory mechanisms include both target gene repressions, such as FLC repression in A. thaliana by H3K9me2 and H3K27me3 followed by polycomb group proteins binding (Sung and Amasino, 2005), and target genes activation such as H3K4me3 and H3 acetylation in association with VRN1 gene activation in Triticeae following vernalization fulfillment (Oliver et al., 2009).
Differential isoforms and PTMs underlie differential enzymatic kinetics/activity and protein-protein interactions under stress treatments. These modifications are crucial for the final protein biological function. For example, the level of O-GlcNAc modification at Thr17 in TaGRP2 glycine-rich RNA binding protein determines its interaction with VER2 Jacalin lectin protein, resulting in TaGRP2 dissociation from VRN1 pre-mRNA first intron and induction of VRN1 transcript accumulation (Xiao et al., 2014). Phosphorylation of penultimate Thr in plasma membrane H+-ATPase leads to binding to phosphorylated C-terminus in 14-3-3 protein involved in the regulation of stomatal aperture in guard cells of A. thaliana rosette leaves (Hayashi et al., 2011). In salinity signaling, a physical interaction of myristoylated Ca2+-binding protein SOS3 with a Ser/Thr protein kinase SOS2 leads to SOS3-SOS2 complex formation and SOS2-mediated phosphorylation of SOS1 plasma membrane-localized Na+/H+ antiporter, resulting in its activation (Sanchez-Barrena et al., 2007). In addition, it was found that SOS2 activity is modulated by the phosphorylation status of 14-3-3 protein. Phosphorylated 14-3-3 protein inhibits SOS2 kinase activity, while 14-3-3 dephosphorylation upon salt stress leads to SOS2 activation (Yu et al., 2016).
Differential peroxiredoxin (Prx) and thioredoxin (Trx) isoforms finely tune redox reactions at the end of electrontransport chains in inner membranes of chloroplast and mitochondria. RubisCO activase (RCA) isoforms were found to finely tune RubisCO activity in heat-treated Rhododendron (Wang et al., 2020).
Besides modulation of enzyme activity, protein isoforms and PTMs are also involved in the modulation of structural protein properties such as chaperone Hsc70 isoforms determining other protein conformations and interactions (Noël et al., 2007), and actin isoforms underlying altered thickness of actin filaments (Kijima et al., 2018); altered patterns of N-glycosylation in cell wall proteins were found in salt-treated A. thaliana (Liu et al., 2021).
It can be concluded that protein isoforms and PTMs are involved in a broad scale of biological processes ranging from cell signaling, transcriptome reprogramming, protein expression, and protein/enzyme function/activity to stress memory mechanisms. Proteoforms represent most probably a universal means of how plants can modulate their biological processes in response to environmental cues; however, the current state of knowledge is far from being complete. Some examples of cellular processes affected by protein isoforms and PTMs in plant responses to environmental stresses are given in Figure 2.
Figure 2. Some examples of cellular processes affected by protein isoforms (blue) and protein posttranslational modifications (PTMs) (red) involved in plant responses to environmental stresses. ER, endoplasmic reticulum; GA, Golgi apparatus; R, receptor protein(s); RCA, RubisCO activase.
Mechanisms of protein posttranscriptional and posttranslational modifications, thus, represent an efficient means of enhancing biological variability in response to variability in an ambient environment. It is, thus, becoming evident that just the sole identification of protein primary sequence encoded by a given gene does not fully define protein final structure and biological activity. In plants, the study of proteoforms is in its beginnings; however, from comparisons of the diversity of some plant proteins, such as Arabidopsis actin isoforms, with their human homologs, i.e., non-muscle actins, it is evident that plant isoforms reveal greater sequence variability than human isoforms that may reflect the greater diversity of environmental factors that plants, as sessile organisms, have to face during their life cycle (Kijima et al., 2018). When interpreting proteomic results, we have to consider differential proteoforms that can be derived from the identified protein and can reveal differential biological functions. The study of proteoforms, thus, indicates that proteome represents a specific level of an organism that cannot be derived from the simple transformation of the transcriptome. The study of proteoforms also represents a key to understanding the versatility of plant responses to environmental cues such as abiotic and biotic stresses.
Author Contributions
KK outlined the idea and prepared the draft manuscript text. PV prepared the figures and tables. PV, MK, IP, and JR actively searched for relevant literature, suggested comments on the manuscript draft, and read and approved the final version of the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by Czech Ministry of Agriculture institutional project MZe CR RO0418 and projects by Czech Ministry of Agriculture QK1710302 and QK1910269, and the Luxembourg National Research Fund (FNR), project SMARTWALL, C15/SR/10240550.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors would like to thank Tereza Nešporová, a Ph.D. student, for her search on plant protein PTM databases and prediction tools.
Footnotes
References
Abat, J. K., and Deswal, R. (2009). Differential modulation of S-nitrosoproteome of Brassica juncea by low temperature: change in S-nitrosylation of Rubisco is responsible for the inactivation of its carboxylase activity. Proteomics 9, 4368–4380. doi: 10.1002/pmic.200800985
Ajadi, A. A., Cisse, A., Ahmad, S., Wang, Y., Shu, Y., Li, S., et al. (2020). Protein phosphorylation and phosphoproteome: an overview. Rice Sci. 27, 184–200. doi: 10.1016/j.rsci.2020.04.003
Barkla, B. J., Vera-Estrella, R., Hernández-Coronado, M., and Pantoja, O. (2009). Quantitative proteomics of the tonoplast reveals a role for glycolytic enzymes in salt tolerance. Plant Cell 21, 4044–4058. doi: 10.1105/tpc.109.069211
Barua, P., Lande, N. V., Subba, P., Gayen, D., Pinto, S., Prasad, T. S. K., et al. (2019). Dehydration-responsive nuclear proteome landscape off chickpea (Cicer arietinum L.) reveals phosphorylation-mediated regulation of stress response. Plant Cell Environ. 42, 230–244. doi: 10.1111/pce.13334
Benlloch, R., and Lois, L. M. (2018). Sumoylation in plants: mechanistic insights and its role in drought stress. J. Exp. Bot. 69, 4539–4554. doi: 10.1093/jxb/ery233
Boisson, B., Giglione, C., and Meinnel, T. (2003). Unexpected protein families including cell defense components feature in the N-myristoylome of a higher eukaryote. J. Biol. Chem. 278, 43418–43429. doi: 10.1074/jbc.M307321200
Bonhomme, L., Valot, B., Tardieu, F., and Zivy, M. (2012). Phosphoproteome dynamics upon changes in plant water status reveal early events associated with rapid growth adjustment in maize leaves. Mol. Cell. Proteomics 11, 957–972. doi: 10.1074/mcp.M111.015867
Brini, F., Hanin, M., Lumbreras, V., Irar, S., Pages, M., and Masmoudi, K. (2007). Functional characterization of DHN-5, a dehydrin showing a differential phosphorylation pattern in two Tunisian durum wheat (Triticum durum Desf.) varieties with marked differences in salt and drought tolerance. Plant Sci. 172, 20–28.
Bueso, E., Rodriguez, L., Lorenzo-Orts, L., Gonzalez-Guzman, M., Sayas, E., Munoz-Bertomeu, J., et al. (2014). The single-subunit RING-type E3 ubiquitin ligase RSL1 targets PYL4 and PYR1 ABA receptors in plasma membrane to modulate abscisic acid signaling. Plant J. 80, 1057–1071. doi: 10.1111/tpj.12708
Bustos, D. M., Bustamante, C. A., and Iglesias, A. A. (2008). Involvement of nonphosphorylating glyceraldehyde-3-phosphate dehydrogenase in response to oxidative stress. J. Plant Physiol. 165, 456–461. doi: 10.1016/j.jplph.2007.06.005
Butturini, E., Boriero, D., Carcerrei de Prati, A., and Mariotto, S. (2019). Immunoprecipitation methods to identify S-glutathionylation in target proteins. MethodsX 6, 1992–1998. doi: 10.1016/j.mex.2019.09.001
Camejo, D., Romero-Puertas, M. D., Rodríguez-Serrano, M., Sandalio, L. M., Lázaro, J. J., Jiménez, A., et al. (2013). Salinity-induced changes in S-nitrosylation of pea mitochondrial proteins. J. Proteomics 79, 87–99. doi: 10.1016/j.jprot.2012.12.003
Carlson, S. M., and Gozani, O. (2014). Emerging technologies to map the protein methylome. J. Mol. Biol. 426, 3350–3362. doi: 10.1016/j.jmb.2014.04.024
Cieśla, A., Mitula, F., Misztal, L., Fedorowicz-Stronska, O., Janicka, S., Tajdel-Zielinska, M., et al. (2016). A role for barley calcium-dependent protein kinase CPK2a in the response to drought. Front. Plant Sci. 7:1550. doi: 10.3389/fpls.2016.01550
Close, T. J. (1997). Dehydrins: a commonalty in the response of plants to dehydration and low temperature. Physiol. Plant. 100, 291–296. doi: 10.1111/j.1399-3054.1997.tp04785.x
Diaz-Vivancos, P., De Simone, A., Kiddle, G., and Foyer, C. H. (2015). Glutathione – linking cell proliferation to oxidative stress. Free Radic. Biol. Med. 89, 1154–1164. doi: 10.1016/j.freeradbiomed.2015.09.023
Dietz, K. J. (2016). Thiol-based peroxidases and ascorbate peroxidases: why plantsrely on multiple peroxidase systems in the photosynthesizing chloroplast? Mol. Cells 39, 20–25. doi: 10.14348/molcells.2016.2324
Ding, L., Yang, R., Yang, G., Cao, J., Li, P., and Zhou, Y. (2016). Identification of putative phosphoproteins in wheat spikes induced by Fusarium graminearum. Planta 243, 719–731. doi: 10.1007/s00425-015-2441-y
Ding, Y., Fromm, M., and Avramova, Z. (2012). Multiple exposures to drought “train” transcriptional responses in Arabidopsis. Nat. Commun. 3:740. doi: 10.1038/ncomms1732
Fang, Y., Deng, X., Lu, X., Zheng, J., Jiang, H., Rao, Y., et al. (2019). Differential phosphoproteome study of the response to cadmium stress in rice. Ecotoxicol. Environ. Saf. 180, 780–788. doi: 10.1016/j.ecoenv.2019.05.068
Feng, B., Liu, C., Shan, L., and He, P. (2016). Protein ADP-ribosylation takes control in plant-bacterium interactions. PLoS Pathog. 12:e1005941. doi: 10.1371/journal.ppat.1005941
Feng, J., Chen, L., and Zuo, J. (2019). Protein S-nitrosylation in plants: current progresses and challenges. J. Integr. Plant Biol. 61, 1206–1223. doi: 10.1111/jipb.12780
Fíla, J., and Honys, D. (2012). Enrichment techniques employed in phosphoproteomics. Amino Acids 43, 1025–1047. doi: 10.1007/s00726-011-1111-z
Ghimire, S., Tang, X., Zhang, N., Liu, W., Qi, X., Fu, X., et al. (2020). Genomic analysis of the SUMO-conjugating enzyme and genes under abiotic stress in potato (Solanum tuberosum L.). Int. J. Genomics 2020:9703638. doi: 10.1155/2020/9703638
Glover, N. M., Redestig, H., and Dessimoz, C. (2016). Homoeologs: what are they and how do we infer them? Trends Plant Sci. 21, 609–621. doi: 10.1016/j.tplants.2016.02.005
Gnad, F., Gunawardena, J., and Mann, M. (2011). PHOSIDA 2011: the posttranslational modification database. Nucleic Acids Res. 39, D253–D260. doi: 10.1093/nar/gkq1159
Gough, C., and Sadanandom, A. (2021). Understanding and exploiting post-translational modifications for plant disease resistance. Biomolecules 11:1122. doi: 10.3390/biom11081122
Guerra, D., Crosatti, C., Khoshro, H. H., Mastrangelo, A. M., Mica, E., and Mazzucotelli, E. (2015). Post-transcriptional and post-translational regulations of drought and heat response in plants: a spider’s web of mechanisms. Front. Plant Sci. 6:57. doi: 10.3389/fpls.2015.00057
Hashiguchi, A., and Komatsu, S. (2016). Impact of post-translational modifications of crop proteins under abiotic stress. Proteomes 4:42. doi: 10.3390/proteomes4040042
Hayashi, M., Inoue, S., Takahashi, K., and Kinoshita, T. (2011). Immunochemical detection of blue light-induced phosphorylation of the plasma membrane H+-ATPase in stomatal guard cells. Plant Cell Physiol. 52, 1238–1248. doi: 10.1093/pcp/pcr072
He, D., Li, M., Damaris, R. N., Bu, C., Xue, J., and Yang, P. (2020). Quantitative ubiquitylomics approach for characterizing the dynamic change and extensive modulation of ubiquitylation in rice seed germination. Plant J. 101, 1430–1447. doi: 10.1111/tpj.14593
Heinemann, B., Künzler, P., Eubel, H., Braun, H. P., and Hildebrandt, T. M. (2020). Estimating the number of protein molecules in a plant cell: protein and amino acid homeostasis during drought. Plant Physiol. 185, 385–404. doi: 10.1093/plphys/kiaa050
Hershko, A., Ciechanover, A., and Varshavsky, A. (2000). The ubiquitin system. Nat. Med. 6, 1073–1081. doi: 10.1038/80384
Hilser, V. J., Anderson, J. A., and Motlagh, H. N. (2015). Allostery vs. “allokairy“. Proc. Natl. Acad. Sci. U.S.A. 112, 11430–11431. doi: 10.1073/pnas.1515239112
Hirschey, M. D., and Zhao, Y. (2015). Metabolic regulation by lysine malonylation, succinylation, and glutarylation. Mol. Cell. Proteomics 14, 2308–2315. doi: 10.1074/mcp.R114.046664
Horie, T., Kaneko, T., Sugimoto, G., Sasano, S., Panda, S. K., Shibasaka, M., et al. (2011). Mechanisms of water transport mediated by PIP aquaporins and their regulation via phosphorylation events under salinity stress in barley roots. Plant Cell Physiol. 52, 663–675. doi: 10.1093/pcp/pcr027
Hornbeck, P. V., Zhang, B., Murray, B., Kornhauser, J. M., Latham, V., and Skrzypek, E. (2015). PhosphoSitePlus, 2014: mutations, PTMs and recalibrations. Nucleic Acids Res. 43, D512–D520. doi: 10.1093/nar/gku1267
Hsu, J. L., Wang, L. Y., Wang, S. Y., Lin, C. H., Ho, K. C., Shi, F. K., et al. (2009). Functional phosphoproteomic profiling of phosphorylation sites in membrane fractions of salts-stressed Arabidopsis thaliana. Proteome Sci. 7:42. doi: 10.1186/1477-5956-7-42
Hu, X., Wu, L., Zhao, F., Zhang, D., Li, N., Zhu, G., et al. (2015b). Phosphoproteomic analysis of the response of maize leaves to drought, heat and their combination stress. Front. Plant Sci. 6:298. doi: 10.3389/fpls.2015.00298
Hu, X., Li, N., Wu, L., Li, C., Zhang, L., Liu, T., et al. (2015a). Quantitative iTRAQ-based proteomic analysis of phosphoproteins and ABA-regulated phosphoproteins in maize leaves under osmotic stress. Sci. Rep. 5:15626. doi: 10.1038/srep15626
Huang, C. H., Su, M. G., Kao, H. J., Jhong, J. H., Weng, S. L., and Lee, T. Y. (2016). UbiSite: incorporating two-layered machine learning method with substrate motifs to predict ubiquitin-conjugation site on lysines. BMC Syst. Biol. 10:S6.
Chaki, M., Carreras, A., López-Jaramillo, J., Begara-Morales, J. C., Sánchez-Calvo, B., Valderrama, R., et al. (2013). Tyrosine nitration provokes inhibition of sunflower carbonic anhydrase (β-CA) activity under high temperature stress. Nitric Oxide 29, 30–33. doi: 10.1016/j.niox.2012.12.003
Chaki, M., Valderrama, R., Fernandez-Ocana, A. M., Carreras, A., Gomez-Rodriguez, M. V., Lopez-Jaramillo, J., et al. (2011). High temperature triggers the metabolism of S.nitrosothiols in sunflower mediating a process of nitrosative stress which provokes the inhibition of ferredoxin-NADP reductase by tyrosine nitration. Plant Cell Environ. 34, 1803–1818. doi: 10.1111/j.1365-3040.2011.02376.x
Chen, Y., Liu, P., Hoehenwarter, W., and Lin, J. (2012). Proteomic and phosphoproteomic analysis of Picea wilsonii pollen development under nutrient limitation. J. Proteome Res. 11, 4180–4190. doi: 10.1021/pr300295m
Chitteti, B. R., and Peng, Z. (2007). Proteome and phosphoproteome differential expression under salinity stress in rice (Oryza sativa) roots. J. Proteome Res. 6, 1718–1727. doi: 10.1021/pr060678z
Ishitani, M., Liu, J., Halfter, U., Kim, C. S., Shi, W., and Zhu, J. K. (2000). SOS3 function in plant salt tolerance requires N-myristoylation and calcium binding. Plant Cell 12, 1667–1677. doi: 10.1105/tpc.12.9.1667
Jaffrey, S. R., and Snyder, S. H. (2001). The biotin switch method for the detection of S-nitrosylated proteins. Sci. STKE 86:pl1. doi: 10.1126/stke.2001.86.pl1
Janská, A., Aprile, A., Zámeèník, J., Cattivelli, L., and Ovesná, J. (2011). Transcriptional responses of winter barley to cold indicate nucleosome remodelling as a specific feature of crown tissues. Funct. Integr. Genomics 11, 307–325. doi: 10.1007/s10142-011-0213-8
Job, C., Rajjou, L., Lovigny, Y., Belghazi, M., and Job, D. (2005). Patterns of protein oxidation in Arabidopsis seeds and during germination. Plant Physiol. 138, 790–802. doi: 10.1104/pp.105.062778
Johansson, E., Olsson, O., and Nyström, T. (2004). Progression and specificity of protein oxidation in the life cycle of Arabidopsis thaliana. J. Biol. Chem. 279, 22204–22208. doi: 10.1074/jbc.M402652200
Jurado-Flores, A., Romero, L. C., and Gotor, C. (2021). Label-free quantitative proteomic analysis of nitrogen stravation in Arabidopsis root reveals new aspects of H2S signaling by protein persulfidation. Antioxidants 10:508. doi: 10.3390/antiox10040508
Kang, J. S., Frank, J., Kang, C. H., Kajiura, H., Vikram, M., Ueda, A., et al. (2008). Salt tolerance of Arabidopsis thaliana requires maturation of N-glycosylated proteins in the Golgi apparatus. Proc. Natl. Acad. Sci. U.S.A. 105, 5933–5938. doi: 10.1073/pnas.0800237105
Ke, Y., Han, G., He, H., and Li, J. (2009). Differential regulation of proteins and phosphoproteins in rice under drought stress. Biochem. Biophys. Res. Commun. 379, 133–138. doi: 10.1016/j.bbrc.2008.12.067
Kijima, S. T., Staiger, C. J., Katoh, K., Nagasaki, A., Ito, K., and Uyeda, T. Q. P. (2018). Arabidopsis vegetative actin isoforms, AtACT2 and AtACT7, generate distinct filament arrays in living plant cells. Sci. Rep. 8:4381. doi: 10.1038/s41598-018-22707-w
Knox, A. K., Dhillon, T., Cheng, H., Tondelli, A., Pecchioni, N., and Stockinger, E. J. (2010). CBF gene copy number variation at frost resistance-2 is associated with levels of freezing tolerance in temperate-climate cereals. Theor. Appl. Genet. 121, 21–35. doi: 10.1007/s00122-010-1288-7
Komatsu, S., Yamada, E., and Furukawa, K. (2009). Cold stress changes the concanavalin A-positive glycosylation pattern of proteins expressed in the basal parts of rice leaf sheaths. Amino Acids 36, 115–123. doi: 10.1007/s00726-008-0039-4
Kosová, K., Prášil, I. T., and Vítámvás, P. (2011). “Role of dehydrins in plant stress response,” in Handbook of Plant and Crop Stress, 3rd Edn, ed. M. Pessarakli (Boca Raton, FL: CRC Press), 239–285.
Kosová, K., Vítámvás, P., and Prášil, I. T. (2021). “Physiological adaptations in temperate crops to environmental constraints during the growing season,” in Handbook of Plant and Crop Physiology, 4th Edn, ed. M. Pessarakli (Boca Raton, FL: Taylor and Francis).
Kosová, K., Vítámvás, P., Urban, M. O., Prášil, I. T., and Renaut, J. (2018). Plant abiotic stress proteomics: the major factors determining alterations in cellular proteome. Front. Plant Sci. 9:122. doi: 10.3389/fpls.2018.00122
Kumar, S. V., and Wigge, P. A. (2010). H2A.Z-containing nucleosomes mediate the thermosensory response in Arabidopsis. Cell 140, 136–147. doi: 10.1016/j.cell.2009.11.006
Lee, H., Guo, Y., Ohta, M., Xiong, L. M., Stevenson, B., and Zhu, J. K. (2002). LOS2, a genetic locus required for cold-responsive gene transcription encodes a bi-functional enolase. EMBO J. 21, 2692–2702. doi: 10.1093/emboj/21.11.2692
Leutert, M., Entwisle, S. W., and Villén, J. (2021). Decoding post-translational modification crosstalk with proteomics. Mol. Cell. Proteomics 20, 100129. doi: 10.1016/j.mcpro.2021.100129
Levi, A., Panta, G. R., Parmentier, C. M., Muthalif, M. M., Arora, R., Shanker, S., et al. (1999). Complementary DNA cloning, sequencing and expression of an unusual dehydrin from blueberry floral buds. Physiol. Plant. 107, 98–109. doi: 10.1034/j.1399-3054.1999.100114.x
Levine, R. L., Garland, D., Oliver, C. N., Amici, A., Climent, I., Lenz, A. G., et al. (1990). Determination of carbonyl content in oxidatively modified proteins. Methods Enzymol. 186, 464–478. doi: 10.1016/0076-6879(90)86141-h
Li, G. X. H., Vogel, C., and Choi, H. (2018). PTMscape: an open source tool to predict generic post-translational modifications and map modification crosstalk in protein domains and biological processes. Mol. Omics 14, 197–209. doi: 10.1039/c8mo00027a
Li, G., Zhang, Z. S., Gao, H. Y., Liu, P., Dong, S. T., Zhang, J. W., et al. (2012). Effects of nitrogen on photosynthetic characteristics of leaves from two different stay-green corn (Zea mays L.) varieties at the grain-filling stage. Can. J. Plant Sci. 92, 671–680. doi: 10.4141/CJPS2012-039
Liu, C., Niu, G., Li, X., Zhang, H., Chen, H., Lan, P., et al. (2021). Comparative label-free quantitative proteomics analysis revealed essential roles of N-glycans in salt tolerance by modulating protein abundance in Arabidopsis. Front. Plant Sci. 12:646425. doi: 10.3389/fpls.2021.646425
Liu, Y., Wang, M., Xi, J., Luo, F., and Li, A. (2018). PTM-ssMP: A web server for predicting different types of post-translational modification sites using novel site-specific modification profile. Int. J. Biol. Sci. 14, 946–956. doi: 10.7150/ijbs.24121
Luo, F., Deng, X., Liu, Y., and Yan, Y. (2018). Identification of phosphorylation proteins in response to water deficit during wheat flag leaf and grain development. Bot. Stud. 59:28. doi: 10.1186/s40529-018-0245-7
Luo, X., Han, C., Deng, X., Zhu, D., Liu, Y., and Yan, Y. (2019). Identification of phosphorylated proteins in response to salt stress in wheat embryo and endosperm during seed germination. Cereal Res. Commun. 47, 53–66. doi: 10.1556/0806.46.02018.061
Lv, D. W., Zhu, G. R., Zhu, D., Bian, Y. W., Liang, X. N., Cheng, Z. W., et al. (2016). Proteomic and phosphoproteomic analysis reveals the response and defense mechanism in leaves of diploid wheat Triticum monococcum under salt stress and recovery. J. Proteomics 143, 93–105. doi: 10.1016/jprot.2016.04.013
Martinez-Seidel, F., Suwanchaikasem, P., Nie, S., Leeming, M. G., Pereira Firmino, A. A., Williamson, N. A., et al. (2021). Membrane-enriched proteomics link ribosome accumulation and proteome reprogramming with cold acclimation in barley root meristems. Front. Plant Sci. 12:656683. doi: 10.3389/fpls.2021.656683
Matamoros, M. A., and Becana, M. (2021). Molecular responses of legumes to abiotic stress: post-translational modifications of proteins and redox signaling. J. Exp. Bot. 72, 5876–5892. doi: 10.1093/jxb/erab008
McLaughlin, S., and Aderem, A. (1995). The myristoyl-electrostatic switch: a modulator of reversible protein-membrane interactions. Trends Biochem. Sci. 20, 272–276. doi: 10.1016/s0968-0004(00)89042-8
Meyer, Y., Reichheld, J. P., and Vignols, F. (2005). Thioredoxins in Arabidopsis andother plants. Photosynth. Res. 86, 419–433. doi: 10.1007/s11120-005-5220-y
Mock, H. P., and Dietz, K. J. (2016). Redox proteomics for the assessment of redox-related posttranslational regulation in plants. Biochim. Biophys. Acta 1864, 967–973. doi: 10.1016/j.bbapap.2016.01.005
Mustafa, G., and Komatsu, S. (2014). Quantitative proteomics reveals the effect of protein glycosylation in soybean root under flooding stress. Front. Plant Sci. 5:627. doi: 10.3389/fpls.2021.656683
Mustafa, G., and Komatsu, S. (2021). Plant proteomic research for improvement of food crops under stresses: a review. Mol. Omics doi: 10.1039/D1MOD00151E
Noël, L. D., Cagna, G., Stuttmann, J., Wirthmüller, L., Betsuyaku, S., Witte, C. P., et al. (2007). Interaction between SGT1 and cytosolic/nuclear HSC70 chaperones regulates Arabidopsis immune responses. Plant Cell 19, 4061–4076. doi: 10.1105/tpc.107.051896
Oliver, S. N., Finnegan, E. J., Dennis, E. S., Peacock, W. J., and Trevaskis, B. (2009). Vernalization induced flowering in cereals is associated with changes in histone methylation at the VERNALIZATION1 gene. Proc. Natl. Acad. Sci. U.S.A. 106, 8386–8391. doi: 10.1073/pnas.0903566106
Pandey, P., and Dubey, R. S. (2021). Cultivar specific oxidative modification of proteins, proteolytic activity and alterations in proteomes of rice seedlings to simultaneous water deficit and aluminum toxicity. J. Plant Growth Regul. doi: 10.1007/s00344-021-10363-w
Pierre, M., Traverso, J. A., Boisson, B., Domenichini, S., Bouchez, D., Giglione, C., et al. (2007). N-myristoylation regulates the SnRK1 pathway in Arabidopsis. Plant Cell 19, 2804–2821. doi: 10.1105/tpc.107.051870
Podell, S., and Gribskov, M. (2004). Predicting N-terminal myristoylation sites in plant proteins. BMC Genomics 5:37. doi: 10.1186/1471-2164-5-37
Rabbani, N., Al-Motawa, M., and Thornalley, P. J. (2020). Protein glycation in plants - an under-researched field with much still to discover. Int. J. Mol. Sci. 21:3942. doi: 10.3390/ijms21113942
Rampitsch, C., and Bykova, N. (2012). The beginnings of crop phosphoproteomics: exploring early warning systems of stress. Front. Plant Sci. 3:144. doi: 10.3389/fpls.2012.00144
Ricard, J., Meunier, J. C., and Buc, J. (1974). Regulatory behavior of monomeric enzymes. 1. The mnemonical enzyme concept. Eur. J. Biotechnol. 49, 195–208. doi: 10.1111/j.1432-1033.1974.tb03825.x
Romero-Puertas, M. C., Palma, J. M., Gómez, M., del Río, L. A., and Sandalio, L. M. (2002). Cadmium causes the oxidative modification of proteins in pea plants. Plant Cell Environ. 25, 677–686. doi: 10.1046/j.1365-3040.2002.00850.x
Romero-Puertas, M. C., Rodriguez-Serrano, M., and Sandalio, L. M. (2013). Protein S-nitrosylation in plants under abiotic stress: an overview. Front. Plant Sci. 4:373. doi: 10.3389/fpls.2013.00373
Ronai, Z., Robinson, R., Rutberg, S., Lazarus, P., and Sardana, M. (1992). Aldolase-DNA interactions in a SEWA cell system. Biochim. Biophys. Acta 1130, 20–28. doi: 10.1016/0167-4781(92)90456-A
Sanchez-Barrena, M. J., Fujii, H., Angulo, I., Martinez-Ripoll, M., Zhu, J. K., and Albert, A. (2007). The structure of the C-terminal domain of the protein kinase atsos2 bound to the calcium sensor atsos3. Mol. Cell 26, 427–435. doi: 10.1016/j.molcel.2007.04.013
Sarge, K. D., and Park-Sarge, O. K. (2009). Detection of proteins sumolyated in vivo and in vitro. Methods Mol. Biol. 590, 265–277. doi: 10.1007/978-1-603-27-378-7_17
Sehrawat, A., Abat, J. K., and Deswal, R. (2013). RuBisCO depletion improved proteome coverage of cold responsive S-nitrosylated targets in Brassica juncea. Front. Plant Sci. 4:342. doi: 10.3389/fpls.2013.00342
Sehrawat, A., and Deswal, R. (2014). S-nitrosylation analysis in Brassica juncea apoplast highlights the importance of nitric oxide in cold-stress signaling. J. Proteome Res. 13, 2599–2619. doi: 10.1021/pr500082u
Shang, X., Cao, Y., and Ma, L. (2017). Alternative splicing in plant genes: a means of regulating the environmental fitness of plants. Int. J. Mol. Sci. 18:432. doi: 10.3390/ijms18020432
Sharma, B., Joshi, D., Yadav, P. K., Gupta, A. K., and Bhatt, T. K. (2016). Role of ubiquitin-mediated degradation system in plant biology. Front. Plant Sci. 7:806. doi: 10.3389/fpls.2016.00806
Silva-Sanchez, C., Li, H., and Chen, S. (2015). Recent advances and challenges in plant phosphoproteomics. Proteomics 15, 1127–1141. doi: 10.1002/pmic.201400410
Singh, R., and Green, M. R. (1993). Sequence-specific binding of transfer-RNA by glyceraldehyde-3-phosphate dehydrogenase. Science 259, 365–368. doi: 10.1126/science.8420004
Smith, L. M., Kelleher, N. L., and The Consortium for Top Down Proteomics (2013). Proteoform: a single term describing protein complexity. Nat. Methods 10, 186–187. doi: 10.1038/nmth.2369
Strasser, R., Altmann, F., Mach, L., Glossl, J., and Steinkellner, H. (2004). Generation of Arabidopsis thaliana plants with complex N-glycans lacking beta-1,2-linked xylose and core alpha-1,3-linked fucose. FEBS Lett. 561, 132–136. doi: 10.1016/S0014-5793(04)00150-4
Subba, P., Barua, P., Kumar, R., Datta, A., Soni, K. K., Chakraborty, S., et al. (2013). Phosphoproteomic dynamics of chickpea (Cicer arietinum L.) reveals shared and distinct components of dehydration response. J. Proteome Res. 12, 5025–5047. doi: 10.1021/pr400628j
Sung, S., and Amasino, R. M. (2005). Remembering winter: toward a molecular understanding of vernalization. Annu. Rev. Plant Biol. 56, 491–508. doi: 10.1146/annurev.arplant.56.032604.144307
Suzuki, Y. J., Carini, M., and Butterfield, D. A. (2010). Protein carbonylation. Antioxid. Redox Signal. 12, 323–325. doi: 10.1089/ars.2009.2887
Tanou, G., Filippou, P., Belghazi, M., Job, D., Diamantidis, G., Fotopoulos, V., et al. (2012). Oxidative and nitrosative-based signalling and associated post-translational modifications orchestrate the acclimation of citrus plants to salinity stress. Plant J. 72, 585–599. doi: 10.1111/j.1365-313X.2012.05100.x
Tanou, G., Job, C., Rajjou, L., Arc, E., Belghazi, M., Diamantidis, G., et al. (2009). Proteomics reveals the overlapping roles of hydrogen peroxide and nitric oxide in the acclimation of citrus plants to salinity. Plant J. 60, 795–804. doi: 10.1111/j.1365-313X.2009.04000.x
Thompson, J. E., Hopkins, M. T., Taylor, C., and Wang, T. W. (2004). Regulation of senescence by eukaryotic translation initiation factor 5A: implications for plant growth and development. Trends Plant Sci. 9, 174–179. doi: 10.1016/j.tplants.2004.02.008
Tondelli, A., Francia, E., Barabaschi, D., Pasquariello, M., and Pecchioni, N. (2011). Inside the CBF locus in Poaceae. Plant Sci. 180, 39–45. doi: 10.1016/j.plantsci.2010.08.012
Vandiver, M. S., Paul, B. D., Xu, R., Karuppagounder, S., Rao, F., Snowman, A. M., et al. (2013). Sulfhydration mediates neuroprotective actions of parkin. Nat. Commun. 4:1626. doi: 10.1038/ncomms2623
Virág, D., Dalmadi-Kiss, B., Vékey, K., Drahos, L., Klebovich, I., Antal, I., et al. (2020). Current trends in the analysis of post-translational modifications. Chromatographia 83, 1–10. doi: 10.1007/s10337-0419-03796-9
Wang, J., Ma, Z., Li, C., Ren, P., Yao, L., Li, B., et al. (2021). Dynamic responses of barley root quantitative succinyl-proteome to short-term phosphate starvation and recovery. Front. Plant Sci. 12:649147. doi: 10.3389/fpls.2021.649147
Wang, X., and Komatsu, S. (2016a). Plant subcellular proteomics: application for exploring optimal cell function in soybean. J. Proteomics 143, 45–56. doi: 10.1016/j.jprot.2016.01.011
Wang, X., and Komatsu, S. (2016b). Gel-free/label-free proteomic analysis of endoplasmic reticulum proteins in soybean root tips under flooding and drought stresses. J. Proteome Res. 15, 2211–2227. doi: 10.1021/acs.jproteome.6b00190
Wang, X., Li, L., Liu, B., Zhou, H., Elmongy, M. S., and Xia, Y. (2020). Combined proteome and transcriptome analysis of heat-primed azalea reveals new insights into plant heat acclimation memory. Front. Plant Sci. 11:1278. doi: 10.3389/fpls.2020.01278
Whittington, A. C., Larion, M., Bowler, J. M., Ramsey, K. M., Brüschweiler, R., and Miller, B. G. (2015). Dual allosteric activation mechanisms in monomeric human glucokinase. Proc. Natl. Acad. Sci. U.S.A. 112, 11553–11558. doi: 10.1073/pnas.1506664112
Willems, P., Horne, A., Van Parys, T., Goormachtig, S., De Smet, I., Botzki, A., et al. (2019). The plant PTM viewer, a central resource for exploring plant protein modifications. Plant J. 99, 752–762.
Wu, X., Gong, F., Cao, D., Hu, X., and Wang, W. (2016). Advances in crop proteomics: PTMs of proteins under abiotic stress. Proteomics 16, 847–865. doi: 10.1002/pmic.201500301
Würschum, T., Longin, C. F., Hahn, V., Tucker, M. R., and Leiser, W. L. (2017). Copy number variations of CBF genes at the Fr-A2 locus are essential components of winetr hardiness in wheat. Plant J. 89, 764–773. doi: 10.1111/tpj.13424
Xiao, J., Xu, S., Li, C., Xu, Y., Xing, L., Niu, Y., et al. (2014). O-GlcNAc-mediated interaction between VER2 and TaGRP2 elicits TaVRN1 mRNA accumulation during vernalization in winter wheat. Nat. Commun. 5:4572. doi: 10.1038/ncomms5572
Xu, W., Wan, J., Zhan, J., Li, X., He, H., Shi, Z., et al. (2017). Global profiling of crotonylation on non-histone proteins. Cell Res. 27, 946–949. doi: 10.1038/cr.2017.60
Xu, Y., Shen, C., Ma, J., Chen, W., Mao, J., Zhou, Y., et al. (2017). Quantitative succinyl-proteome profiling of Camellia sinensis cv. ‘Anji Baicha’ during periodic albinism. Sci. Rep. 7:187. doi: 10.1038/s41598-017-02128-x
Xue, L., Wang, P., Wang, L., Renzi, E., Radivojac, P., Tang, H., et al. (2013). Quantitative measurement of phosphoproteome response to osmotic stress in Arabidopsis based on library-assisted extracted ion chromatogram (LAXIC). Mol. Cell. Proteomics 12, 2354–2369. doi: 10.1074/mcp.O113.027284
Yakubov, B., Barazani, O., Shachack, A., Rowland, L. J., Shoseyov, O., and Golan-Goldhirsh, A. (2005). Cloning and expression of a dehydrin-like protein from Pistacia vera L. Trees 19, 224–230. doi: 10.1007/s00468-004-0385-0
Yamaguchi-Shinozaki, K., and Shinozaki, K. (2006). Transcriptional regulatory networks in cellular responses and tolerance to dehydration and cold stresses. Annu. Rev. Plant Biol. 57, 781–803. doi: 10.1146/annurev.arplant.57.032905.105444
Yang, I., Liu, Q., Liu, Z., Yang, H., Wang, J., Li, X., et al. (2015). Arabidopsis C3HC4-RING finger E3 ubiquitin ligase AtAIRP4 positively regulates stress-responsive abscisic acid signalling. J. Integr. Plant Biol. 58, 67–80. doi: 10.1111/jipb.12364
Yang, Z. B., Eticha, D., Führs, H., Heintz, D., Ajoub, D., Van Dorsselaer, A., et al. (2013). Proteomic and phosphoproteomic analysis of polyethylene glycol-induced osmotic stress in root tips of common bean (Phaseolus vulgaris L.). J. Exp. Bot. 64, 5569–5586. doi: 10.1093/jxb/ert328
Yin, X., and Komatsu, S. (2015). Quantitative proteomics of nuclear phosphoproteins in the root tip of soybean during the initial stages of flooding stress. J. Proteomics 119, 183–195. doi: 10.1016/j.prot.2015.02.004
Yin, X., Sakata, K., and Komatsu, S. (2014). Phosphoproteomics reveals the effect of ethylene in soybean root under flooding stress. J. Proteome Res. 13, 5618–5634. doi: 10.1021/pr500621c
Yu, B., Li, J., Koh, J., Dufresne, C., Yang, N., Qi, S., et al. (2016). Quantitative proteomics and phosphoproteomics of sugar beet monosomic addition line M14 in response to salt stress. J. Proteomics 143, 286–297. doi: 10.1016/j.jprot.2016.04.011
Yu, C., Wu, Q., Sun, C., Tang, M., Sun, J., and Zhan, Y. (2019). The phosphoproteomic response of okra (Abelmoschus esculentus L.) seedlings to salt stress. Int. J. Mol. Sci. 20:1262. doi: 10.3390/ijms20061262
Yuan, L. L., Zhang, M., Yan, X., Bian, Y. W., Zhen, S. M., and Yan, Y. M. (2016). Dynamic phosphoproteome analysis of seedling leaves in Brachypodium distachyon L. reveals central phosphorylated proteins involved in the drought stress response. Sci. Rep. 6:35280. doi: 10.1038/srep35280
Zhang, H., Cui, F., Wu, Y., Lou, I., Tian, M., Ning, Y., et al. (2015). The Ring finger ubiquitin E3 ligase SDIR1 targets SDIR1-interacting protein 1 for degradation to modulate the salt stress response and ABA signaling in Arabidopsis. Plant Cell 27, 214–227. doi: 10.1105/tpc.114.134163
Zhang, M., Lv, D., Ge, P., Bian, Y., Chen, G., Zhu, G., et al. (2014a). Phosphoproteome analysis reveals new drought response and defense mechanisms of seedling leaves in bread wheat (Triticum aestivum L.). J. Proteomics 109, 290–308. doi: 10.1016/j.jprot.2014.07.010
Zhang, M., Ma, C. Y., Lv, D. W., Zhen, S. M., Li, X. H., and Yan, Y. M. (2014b). Comparative phosphoproteome analysis of the developing grains in bread wheat (Triticum aestivum L.) under well-watered and water-deficit conditions. J. Proteome Res. 13, 4281–4297. doi: 10.1021/pr500400t
Zhong, M., Li, S. F., Huang, F. L., Qiu, J. H., Zhang, J., Sheng, Z. H., et al. (2017). The phosphoproteomic response of rice seedlings to cadmium stress. Int. J. Mol. Sci. 18:2055. doi: 10.3390/ijms18102055
Keywords: protein isoforms, protein posttranslational modifications (PTMs), protein-protein interactions, biological functions, environmental stresses, crops
Citation: Kosová K, Vítámvás P, Prášil IT, Klíma M and Renaut J (2021) Plant Proteoforms Under Environmental Stress: Functional Proteins Arising From a Single Gene. Front. Plant Sci. 12:793113. doi: 10.3389/fpls.2021.793113
Received: 11 October 2021; Accepted: 16 November 2021;
Published: 14 December 2021.
Edited by:
Sebastien Christian Carpentier, Bioversity International, BelgiumReviewed by:
Dominique Job, UMR 5240 Microbiologie, Adaptation et Pathogenie (MAP), FranceXiaojian Yin, China Pharmaceutical University, China
Copyright © 2021 Kosová, Vítámvás, Prášil, Klíma and Renaut. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Klára Kosová, a29zb3ZhQHZ1cnYuY3o=