Physcomitrella patens DNA methyltransferase 2is required for recovery from salt and osmotic stress

Potential Role of DNA Methylation as a Driver of Plastic Responses to the Environment Across Cells, Organisms, and Populations

There is great interest in exploring epigenetic modifications as drivers of adaptive organismal responses to environmental change. Extending this hypothesis to populations, epigenetically driven plasticity could influence phenotypic changes across environments. The canonical model posits that epigenetic modifications alter gene regulation and subsequently impact phenotypes. We first discuss origins of epigenetic variation in nature, which may arise from genetic variation, spontaneous epimutations, epigenetic drift, or variation in epigenetic capacitors. We then review and synthesize literature addressing three facets of the aforementioned model: (i) causal effects of epigenetic modifications on phenotypic plasticity at the organismal level, (ii) divergence of epigenetic patterns in natural populations distributed across environmental gradients, and (iii) the relationship between environmentally induced epigenetic changes and gene expression at the molecular level. We focus on DNA methylation, the most extensively studied epigenetic modification. We find support for environmentally associated epigenetic structure in populations and selection on stable epigenetic variants, and that inhibition of epigenetic enzymes frequently bears causal effects on plasticity. However, there are pervasive confounding issues in the literature. Effects of chromatin-modifying enzymes on phenotype may be independent of epigenetic marks, alternatively resulting from functions and protein interactions extrinsic of epigenetics. Associations between environmentally induced changes in DNA methylation and expression are strong in plants and mammals but notably absent in invertebrates and nonmammalian vertebrates. Given these challenges, we describe emerging approaches to better investigate how epigenetic modifications affect gene regulation, phenotypic plasticity, and divergence among populations.

Arabidopsis T-DNA mutants affected in TRDMT1/DNMT2 show differential protein synthesis and compromised stress tolerance

TRDMT1/DNMT2 belongs to the conserved family of nucleic acid methyltransferases. Unlike the animal systems, studies on TRDMT1/DNMT2 in land plants have been limited. We show that TRDMT1/DNMT2 is strongly conserved in the green lineage. Studies in mosses have previously shown that TRDMT1/DNMT2 plays a crucial role in modulating molecular networks involved in stress perception and signalling and in transcription/stability of specific tRNAs under stress. To gain deeper insight into its biological roles in a flowering plant, we examined more closely the previously reported Arabidopsis SALK_136635C line deficient in TRDMT1/DNMT2 function [Goll MG et al. (2006) Science 311, 395-398]. RNAs derived from Arabidopsis Dnmt2-deficient plants lacked m5 C38 in tRNAAsp . In this study, by transient expression assays we show that Arabidopsis TRDMT1/DNMT2 is distributed in the nucleus, cytoplasm and RNA-processing bodies, suggesting a role for TRDMT1/DNMT2 in RNA metabolic processes possibly by shuttling between cellular compartments. Bright-field and high-resolution SEM and qPCR analysis reveal roles of TRDMT1/DNMT2 in proper growth and developmental progression. Quantitative proteome analysis by LC-MS/MS coupled with qPCR shows AtTRDMT1/AtDNMT2 function to be crucial for protein synthesis and cellular homeostasis via housekeeping roles and proteins with poly-Asp stretches and RNA pol II activity on selected genes are affected in attrdmt1/atdnmt2. This shift in metabolic pathways primes the mutant plants to become increasingly sensitive to oxidative and osmotic stress. Taken together, our study sheds light on the mechanistic role of TRDMT1/DNMT2 in a flowering plant.

Restricted tRNA methylation by intermolecular disulfide bonds in DNMT2/TRDMT1

Reportedly, DNMT2/TRDMT1 mainly methylates tRNAs at C38 and prevents them from the cleavage under stress. It also plays an essential role in the survival and physiological homeostasis of organisms. Nevertheless, DNMT2/TRDMT1 exhibits much weaker tRNA methylation activity in vitro than other tRNA methyltransferases, TrmD and Trm5. Here, we explored the restricted tRNA methylation mechanism by DNMT2/TRDMT1. In the current study, the optimized buffer C at 37 °C was the best condition for DNMT2/TRDMT1 activation. Of note, Dithiothreitol (DTT) was an indispensable component for this enzyme catalysis. Moreover, reductants took similar effects on the conformation change and oligomeric formation of DNMT2/TRDMT1. Ultimately, LC-MS/MS result revealed that C292-C292 and C292-C287 were predominant intermolecular disulfide bonds in recombinant DNMT2/TRDMT1. Notably, DNMT2/TRDMT1 existed primarily as dimers via intermolecular disulfide bonds C79-C24, C292-C292, and C222-C24 in HEK293T cells. GSSG stress enhanced tRNA methylation level in the early stage of stress, whereas the DNMT2/TRDMT1 activity might be unfavorable along with this enzyme accumulation in the nucleus. Excitingly, GSH stress downregulated the DNMT2/TRDMT1 expression and promoted tRNA methylation in cells, probably through breaking intermolecular disulfide bonds in this enzyme. Thus, our findings demonstrated restricted tRNA methylation by disulfide bonds in DNMT2/TRDMT1, and will provide important implications for redox stress related-diseases.

Role of TRDMT1/DNMT2 in stress adaptation and its influence on transcriptome and proteome dynamics under osmotic stress in Physcomitrium patens

Early land plants such as the moss Physcomitrium patens lack several morphological traits that offer protection to tracheophytes from environmental stresses. These plants instead have evolved several physiological and biochemical mechanisms that facilitate them to adapt to terrestrial stresses such as drought. We have previously shown that loss-of-function mutants of tRNA (cytosine(38)-C(5))-methyltransferase TRDMT1/DNMT2 in P. patens are highly sensitive to oxidative and osmotic stress. To gain insight into the role of PpTRDMT1/PpDNMT2 in modulating genetic networks under osmotic stress, genome-wide transcriptome and proteome studies were undertaken in wild-type and ppdnmt2 plants. Transcriptome analysis revealed 375 genes to be differentially expressed in the ppdnmt2 under stress compared to the WT. Most of these genes are affiliated with carbohydrate metabolic pathways, photosynthesis, cell wall biogenesis, pathways related to isotropic and polarised cell growth and transcription factors among others. Histochemical staining showed elevated levels of reactive oxygen species in ppdnmt2 while transmission electron microscopy revealed no distinct defects in the ultrastructure of chloroplasts. Immunoprecipitation using PpDNMT2-specific antibody coupled with mass spectrometry revealed core proteins of the glycolytic pathway, antioxidant enzymes, proteins of amino acid biosynthetic pathways and photosynthesis-related proteins among others to co-purify with PpTRDMT1/PpDNMT2 under osmotic stress. Yeast two-hybrid assays, protein deletion and α-galactosidase assays showed the cytosol glycolytic protein glyceraldehyde 3-phosphate dehydrogenase to bind to the catalytic motifs in PpTRDMT1/PpDNMT2. Results presented in this study allow us to better understand genetic networks linking enzymes of energy metabolism, epigenetic processes and RNA pol II-mediated transcription during osmotic stress tolerance in P. patens.

Elucidation of transfer RNAs as stress regulating agents and the experimental strategies to conceive the functional role of tRNA-derived fragments in plants

In plants, the transfer RNAs (tRNAs) exhibit their profound influence in orchestrating diverse physiological activities like cell growth, development, and response to several surrounding stimuli. The tRNAs, which were known to restrict their function solely in deciphering the codons, are now emerging as frontline defenders in stress biology. The plants that are constantly confronted with a huge panoply of stresses rely on tRNA-mediated stress regulation by altering the tRNA abundance, curbing the transport of tRNAs, fragmenting the mature tRNAs during stress. Among them, the studies on the generation of transfer RNA-derived fragments (tRFs) and their biological implication in stress response have attained huge interest. In plants, the tRFs hold stable expression patterns and regulate biological functions under diverse environmental conditions. In this review, we discuss the fate of plant tRNAs upon stress and thereafter how the tRFs are metamorphosed into sharp ammunition to wrestle with stress. We also address the various methods developed to date for uncovering the role of tRFs and their function in plants.

Determinants of DNMT2/TRDMT1 preference for substrates tRNA and DNA during the evolution

RNA methyltransferase DNMT2/TRDMT1 is the most conserved member of the DNMT family from bacteria to plants and mammals. In previous studies, we found some determinants for tRNA recognition of DNMT2/TRDMT1, but the preference mechanism of this enzyme for substrates tRNA and DNA remains to be explored. In the present study, CFT-containing target recognition domain (TRD) and target recognition extension domain (TRED) in DNMT2/TRDMT1 play a crucial role in the substrate DNA and RNA selection during the evolution. Moreover, the classical substrate tRNA for DNMT2/TRDMT1 had a characteristic sequence CUXXCAC in the anticodon loop. Position 35 was occupied by U, making cytosine-38 (C38) twist into the loop, whereas C, G or A was located at position 35, keeping the C38-flipping state. Hence, the substrate preference could be modulated by the easily flipped state of target cytosine in tRNA, as well as TRD and TRED. Additionally, DNMT2/TRDMT1 cancer mutant activity was collectively mediated by five enzymatic characteristics, which might impact gene expressions. Importantly, G155C, G155V and G155S mutations reduced enzymatic activities and showed significant associations with diseases using seven prediction methods. Altogether, these findings will assist in illustrating the substrate preference mechanism of DNMT2/TRDMT1 and provide a promising therapeutic strategy for cancer.

RNA methyltransferases in plants: Breakthroughs in function and evolution

Each day it is becoming increasingly difficult not to notice the completely new, fast growing, extremely intricate and challenging world of epitranscriptomics as the understanding of RNA methylation is expanding at a hasty rate. Writers (methyltransferases), erasers (demethylases) and readers (RNA-binding proteins) are responsible for adding, removing and recognising methyl groups on RNA, respectively. Several methyltransferases identified in plants are now being investigated and recent studies have shown a connection between RNA-methyltransferases (RNA-MTases) and stress and development processes. However, compared to their animal and bacteria counterparts, the understanding of RNA methyltransferases is still incipient, particularly those located in organelles. Comparative and systematic analyses allowed the tracing of the evolution of these enzymes suggesting the existence of several methyltransferases yet to be characterised. This review outlines the functions of plant nuclear and organellar RNA-MTases in plant development and stress responses and the comparative and evolutionary discoveries made on RNA-MTases across kingdoms.

The DEAD-box RNA helicase eIF4A1 interacts with the SWI2/SNF2-related chromatin remodelling ATPase DDM1 in the moss Physcomitrella

eIF4A is a DEAD box containing RNA helicase that plays crucial roles in regulating translation initiation, growth and abiotic stress tolerance in plants. It also functions as an ATP-dependent RNA binding protein to curb granule formation by limiting RNA-RNA interactions that promote RNA condensation and formation of ribonucleoprotein particles in vivo. Helicase activity of eIF4A is known to be dictated by its binding partners. Proteins interacting with eIF4A have been identified across land plants. In monocots a close link between eIF4A regulated processes and DNA methylation in epigenetic regulation of plant development is inferred from interaction between OseIF4A and the de novo methyltransferase OsDRM2 and loss-of-function studies of these genes in Oryza sativa and Brachypodium distachyon. In the moss Physcomitrella patens, eIF4A1 encoded by Pp3c6_1080V3.1 interacts with the heterogeneous nuclear ribonucleoprotein (hnRNP) PpLIF2L1, homolog of which in Arabidopsis regulates transcription of stress-responsive genes. In this study, using different protein-protein interaction methods, targeted gene knockout strategy and quantitative expression analysis we show genetic interaction between PpeIF4A1 and the putative nucleosome remodeler protein PpDDM1 and between PpDDM1 and PpLIF2L1 in vivo. Stress-induced co-expression of PpeIF4A1, PpDDM1 and PpLIF2L1, their roles in salt stress tolerance and differences in subnuclear distribution of PpLIF2L1 in ppeif4a1 cells in comparison to wild type suggest existence of a regulatory network comprising of RNA helicases, chromatin remodelling proteins and hnRNP active in stress-responsive biological processes in P. patens.

Small DNA Methylation, Big Player in Plant Abiotic Stress Responses and Memory

DNA methylation is a conserved epigenetic mark that plays important roles in maintaining genome stability and regulating gene expression. As sessile organisms, plants have evolved sophisticated regulatory systems to endure or respond to diverse adverse abiotic environmental challenges, i.e., abiotic stresses, such as extreme temperatures (cold and heat), drought and salinity. Plant stress responses are often accompanied by changes in chromatin modifications at diverse responsive loci, such as 5-methylcytosine (5mC) and N 6-methyladenine (6mA) DNA methylation. Some abiotic stress responses are memorized for several hours or days through mitotic cell divisions and quickly reset to baseline levels after normal conditions are restored, which is referred to as somatic memory. In some cases, stress-induced chromatin marks are meiotically heritable and can impart the memory of stress exposure from parent plants to at least the next stress-free offspring generation through the mechanisms of transgenerational epigenetic inheritance, which may offer the descendants the potential to be adaptive for better fitness. In this review, we briefly summarize recent achievements regarding the establishment, maintenance and reset of DNA methylation, and highlight the diverse roles of DNA methylation in plant responses to abiotic stresses. Further, we discuss the potential role of DNA methylation in abiotic stress-induced somatic memory and transgenerational inheritance. Future research directions are proposed to develop stress-tolerant engineered crops to reduce the negative effects of abiotic stresses.

Comprehensive analysis of the Ppatg3 mutant reveals that autophagy plays important roles in gametophore senescence in Physcomitrella patens

BACKGROUND

Autophagy is an evolutionarily conserved system for the degradation of intracellular components in eukaryotic organisms. Autophagy plays essential roles in preventing premature senescence and extending the longevity of vascular plants. However, the mechanisms and physiological roles of autophagy in preventing senescence in basal land plants are still obscure.

RESULTS

Here, we investigated the functional roles of the autophagy-related gene PpATG3 from Physcomitrella patens and demonstrated that its deletion prevents autophagy. In addition, Ppatg3 mutant showed premature gametophore senescence and reduced protonema formation compared to wild-type (WT) plants under normal growth conditions. The abundance of nitrogen (N) but not carbon (C) differed significantly between Ppatg3 mutant and WT plants, as did relative fatty acid levels. In vivo protein localization indicated that PpATG3 localizes to the cytoplasm, and in vitro Y2H assays confirmed that PpATG3 interacts with PpATG7 and PpATG12. Plastoglobuli (PGs) accumulated in Ppatg3, indicating that the process that degrades damaged chloroplasts in senescent gametophore cells was impaired in this mutant. RNA-Seq uncovered a detailed, comprehensive set of regulatory pathways that were affected by the autophagy mutation.

CONCLUSIONS

The autophagy-related gene PpATG3 is essential for autophagosome formation in P. patens. Our findings provide evidence that autophagy functions in N utilization, fatty acid metabolism and damaged chloroplast degradation under non-stress conditions. We identified differentially expressed genes in Ppatg3 involved in numerous biosynthetic and metabolic pathways, such as chlorophyll biosynthesis, lipid metabolism, reactive oxygen species removal and the recycling of unnecessary proteins that might have led to the premature senescence of this mutant due to defective autophagy. Our study provides new insights into the role of autophagy in preventing senescence to increase longevity in basal land plants.

Transcriptome Analysis of ppdnmt2 and Identification of Superoxide Dismutase as a Novel Interactor of DNMT2 in the Moss Physcomitrella patens

DNMT2 is a DNA/tRNA cytosine methyltransferase that is highly conserved in structure and function in eukaryotes. In plants however, limited information is available on the function of this methyltransferase. We have previously reported that in the moss Physcomitrella patens, DNMT2 plays a crucial role in stress recovery and tRNA^Asp transcription/stability under salt stress. To further investigate the role of PpDNMT2 at genome level, in this study we have performed RNA sequencing of ppdnmt2. Transcriptome analysis reveals a number of genes and pathways to function differentially and suggests a close link between PpDNMT2 function and osmotic and ionic stress tolerance. We propose PpDNMT2 to play a pivotal role in regulating salt tolerance by affecting molecular networks involved in stress perception and signal transduction that underlie maintenance of ion homeostasis in cells. We also examined interactome of PpDNMT2 using affinity purification (AP) coupled to mass spectrometry (AP-MS). Quantitative proteomic analysis reveals several chloroplast proteins involved in light reactions and carbon assimilation and proteins involved in stress response and some not implicated in stress to co-immunoprecipitate with PpDNMT2. Comparison between transcriptome and interactome datasets has revealed novel association between PpDNMT2 activity and the antioxidant enzyme Superoxide dismutase (SOD), protein turnover mediated by the Ubiquitin-proteasome system and epigenetic gene regulation. PpDNMT2 possibly exists in complex with CuZn-SODs in vivo and the two proteins also directly interact in the yeast nucleus as observed by yeast two-hybrid assay. Taken together, the work presented in this study sheds light on diverse roles of PpDNMT2 in maintaining molecular and physiological homeostasis in P. patens. This is a first report describing transcriptome and interactome of DNMT2 in any land plant.

The DEAD-box RNA helicase eIF4A regulates plant development and interacts with the hnRNP LIF2L1 in Physcomitrella patens

eIF4A is a RNA-stimulated ATPase and helicase. Besides its key role in regulating cap-dependent translation initiation in eukaryotes, it also performs specific functions in regulating cell cycle progression, plant growth and abiotic stress tolerance. Flowering plants encode three eIF4A paralogues, eIF4A1, eIF4A2 and eIF4A3 that share conserved sequence motifs but differ in functions. To date, however, no information is available on eIF4A in basal land plants. In this study we report that genome of the moss Physcomitrella patens encodes multiple eIF4A genes. The encoded proteins possess the highly conserved motifs characteristic of the DEAD box helicases. Spatial expression analysis shows these genes to be ubiquitously expressed in all tissue types with Pp3c6_1080V3.1 showing high expression in filamentous protonemata. Targeted deletion of conserved core motifs in Pp3c6_1080V3.1 slowed protonemata growth and resulted in dwarfing of leafy gametophores suggesting a role for Pp3c6_1080V3.1 in regulating cell division/elongation. Rapid and strong induction of Pp3c6_1080V3.1 under salt stress and slow recovery of knockout plants upon exposure to high salt further suggest Pp3c6_1080V3.1 to be involved in stress management in P. patens. Protein-protein interaction studies that show Pp3c6_1080V3.1 to interact with the Physcomitrella heterogenous ribonucleoprotein, LIF2L1, a transcriptional regulator of stress-responsive genes in Arabidopsis. The results presented in this study provide insight into evolutionary conserved functions of eIF4A and shed light on the novel link between eIF4A activities and stress mitigation pathways/RNA metabolic processes in P. patens.

Transgenerational Response to Heat Stress in the Form of Differential Expression of Noncoding RNA Fragments in Brassica rapa Plants

Epigenetic regulations in the form of changes in differential expression of noncoding RNAs (ncRNAs) are an essential mechanism of stress response in plants. Previously we showed that heat treatment in L. results in the differential processing and accumulation of ncRNA fragments (ncRFs) stemming from transfer RNAs (tRNAs), ribosomal RNAs (rRNAs), small nuclear RNAs (snRNAs), and small nucleolar RNAs (snoRNAs). In this work, we analyzed whether ncRFs are differentially expressed in the progeny of heat-stressed plants. We found significant changes in the size of tRF reads and a significant decrease in the percentage of tRFs mapping to tRNA-Ala, tRNA-Arg, and tRNA-Tyr and an increase in tRFs mapping to tRNA-Asp. The enrichment analysis showed significant differences in processing of tRFs from tRNA, tRNA, tRNA, tRNA, tRNA, and tRNA isoacceptors. Analysis of potential targets of tRFs showed that they regulate brassinosteroid metabolism, the proton pump ATPase activity, the antiporter activity, the mRNA decay activity as well as nucleosome positioning and the epigenetic regulation of transgenerational response. Gene ontology term analysis of potential targets demonstrated a significant enrichment in tRFs that potentially targeted a cellular component endoplasmic reticulum (ER) and in small nucleolar RNA fragments (snoRFs), the molecular function protein binding. To summarize, our work demonstrated that the progeny of heat-stressed plants exhibit changes in the expression of tRFs and snoRFs but not of small nuclear RNA fragments (snRFs) or ribosomal RNA fragments (rRFs) and these changes likely better prepare the progeny of stressed plants to future stress encounters.

Cytosine methylation by DNMT2 facilitates stability and survival of HIV-1 RNA in the host cell during infection

The enigmatic methyltransferase, DNMT2 (DNA methyltransferase 2), structurally resembles a DNA methyltransferase, but has been shown to be a tRNA methyltransferase targeting cytosine within a specific CpG in different tRNA molecules. We had previously shown that, during environmental stress conditions, DNMT2 is re-localized from the nucleus to the cytoplasmic stress granules (SGs) and is associated with RNA-processing proteins. In the present study, we show that DNMT2 binds and methylates various mRNA species in a sequence-independent manner and gets re-localized to SGs in a phosphorylation-dependent manner. Importantly, our results indicate that HIV-1 enhances its survivability in the host cell by utilizing this RNA methylation capability of DNMT2 to increase the stability of its own genome. Upon infection, DNMT2 re-localizes from the nucleus to the SGs and methylates HIV-1 RNA. This DNMT2-dependent methylation provided post-transcriptional stability to the HIV-1 RNA. Furthermore, DNMT2 overexpression increased the HIV-1 viral titre. This would suggest that HIV hijacks the RNA-processing machinery within the SGs to ensure its own survival in the host cell. Thus, our findings provide for a novel mechanism by which virus tries to modulate the host cell machinery to its own advantage.

Cross-Talk between Dnmt2-Dependent tRNA Methylation and Queuosine Modification

Enzymes of the Dnmt2 family of methyltransferases have yielded a number of unexpected discoveries. The first surprise came more than ten years ago when it was realized that, rather than being DNA methyltransferases, Dnmt2 enzymes actually are transfer RNA (tRNA) methyltransferases for cytosine-5 methylation, foremost C38 (m⁵C38) of tRNA^Asp. The second unanticipated finding was our recent discovery of a nutritional regulation of Dnmt2 in the fission yeast Schizosaccharomyces pombe. Significantly, the presence of the nucleotide queuosine in tRNA^Asp strongly stimulates Dnmt2 activity both in vivo and in vitro in S. pombe. Queuine, the respective base, is a hypermodified guanine analog that is synthesized from guanosine-5’-triphosphate (GTP) by bacteria. Interestingly, most eukaryotes have queuosine in their tRNA. However, they cannot synthesize it themselves, but rather salvage it from food or from gut microbes. The queuine obtained from these sources comes from the breakdown of tRNAs, where the queuine ultimately was synthesized by bacteria. Queuine thus has been termed a micronutrient. This review summarizes the current knowledge of Dnmt2 methylation and queuosine modification with respect to translation as well as the organismal consequences of the absence of these modifications. Models for the functional cooperation between these modifications and its wider implications are discussed.

Overexpression of GhSARP1 encoding a E3 ligase from cotton reduce the tolerance to salt in transgenic Arabidopsis

Ubiquitination plays a very important role in the response to abiotic stresses of plant. To identify key regulators of salt stress, a gene GhSARP1(Salt-Associated Ring finger Protein)encoding C3H2C3-type E3 ligase, was cloned from cotton. Transcription level of GhSARP1 was high in leaf, flower and fiber of 24,27 and 27DPA (Days Post-Anthesis), but low in root and stem. Except PEG6000 treatment, the expression of GhSARP1 was down-regulated by NaCl, cold and ABA after being treated for 1 h. GhSARP1-GFP fusion protein located on the plasma membrane, which was dependent on trans-membrane motif. In vitro ubiquitination assay showed that GhSARP1 had E3 ligase activity. Heterogeneous overexpression of GhSARP1reduced salt tolerance of transgenic Arabidopsis in germination and post-germination stage. Our results suggested that the GhSARP1 might negatively regulate the response to salt stress mediated by the ubiquitination in cotton.

GpDSR7, a Novel E3 Ubiquitin Ligase Gene in Grimmia pilifera Is Involved in Tolerance to Drought Stress in Arabidopsis

The growth and development of plants under drought stress depends mainly on the expression levels of various genes and modification of proteins. To clarify the molecular mechanism of drought-tolerance of plants, suppression subtractive hybridisation cDNA libraries were screened to identify drought-stress-responsive unigenes in Grimmia pilifera, and a novel E3 ubiquitin ligase gene, GpDSR7, was identified among the 240 responsive unigenes. GpDSR7 expression was induced by various abiotic stresses, particularly by drought. GpDSR7 displayed E3 ubiquitin ligase activity in vitro and was exclusively localised on the ER membrane in Arabidopsis mesophyll protoplasts. GpDSR7-overexpressing transgenic Arabidopsis plants showed a high water content and survival ratio under drought stress. Moreover, the expression levels of some marker genes involved in drought stress were higher in the transgenic plants than in wild-type plants. These results suggest that GpDSR7, an E3 ubiquitin ligase, is involved in tolerance to drought stress at the protein modification level.