Expression of genes for selected plant aminoacyl-tRNA synthetases in the abiotic stress

Introduction

As sessile organisms, plants have evolved elaborate mechanisms to adapt to adverse environmental changes and challenges. The molecular responses in plants to abiotic stresses are probably more advanced and prominent than in animals (Qin et al. 2011). Research on how plants adjust to various environmental stresses is crucial not only to plant biologists but also to agronomists, because abiotic stress has a particularly negative impact on crop production (Cramer et al. 2011). The need to understand the mechanisms involved in plant abiotic stress responses and the generation of plants tolerant to stress has accordingly gained much attention in recent years (Qin et al. 2011). In plants, as in all eukaryotes, gene expression reprogramming plays a pivotal role in the response to abiotic stress. This reprogramming is a highly complex process, including transcription, mRNA processing, mRNA transport and stability, translation or protein turnover; each participating but having a different weight in the final protein levels (Yángüez et al. 2013). One of the earliest plant metabolic responses to abiotic stresses is the inhibition of protein biosynthesis and an increase in chaperone levels controlling protein folding and processing. Repression of global protein biosynthesis is often accompanied with selective translation of mRNAs encoding proteins that are vital for cell survival and recovery from stress (Liu and Qian 2014, Merchante et al. 2017).

In the first step of protein biosynthesis, cells need to aminoacylate tRNAs with their respective amino acids in a process catalyzed by aminoacyl-tRNA synthetases (aaRSs) (Perona and Gruic-Sovulj 2014). During protein biosynthesis on the ribosome, amino acid is transferred from the tRNA onto a growing peptide according to the genetic code. Aminoacyl-tRNA synthetases are therefore housekeeping proteins that play an important role in the expression of genes to create proteins. Besides that, aaRSs may be involved in a myriad of other cellular processes exerting their noncanonical function beyond translation (Guo and Schimmel 2013, Mocibob et al. 2016). It is well documented that aaRSs participate in cellular stress responses in bacteria and metazoa. For example, vertebrate tyrosyl-tRNA synthetase has been shown to translocate to the nucleus during oxidative stress and protects against DNA damage by activating transcription factor E2F1, which promotes expression of DNA repair genes (Wei et al. 2014, Cao et al. 2017). During oxidative stress, phosphorylated methionyl-tRNA synthetase (MetRS) enhances the mischarging of methionine on non-methionyl tRNAs. Methionine carried by noncognate tRNAs is incorporated into growing polypeptides during translation and is used as scavanger of reactive oxygen species, protecting cells from oxidative damage and apoptosis (Lee et al. 2014). During oxidative stress, oxidized phenylalanyl-tRNA synthetase (PheRS) exhibits increased proofreading activity against cytotoxic tyrosine isomers, noncognate amino acids that are increased during oxidative stress. Increased proofreading of PheRS maintains high fidelity of translation even under conditions of oxidative stress (Steiner et al. 2019).

Research on plant aaRSs has revealed their cellular localization (Rokov-Plavec et al. 2008, Duchêne et al. 2009, Kekez et al. 2016), substrate specificity (Rokov et al. 1998, Rokov and Weygand-Durasevic 1999, Rokov-Plavec et al. 2002, 2004, Aldinger et al. 2012), fidelity (Rokov-Plavec et al. 2013, Lee et al. 2016, Hoffman et al. 2019) and protein interactors (Yang et al. 2018, Kekez et al. 2019). Thus far, only one crystal structure of plant aaRS has been reported (Kekez et al. 2019). As enzymes essential for protein biosynthesis in the cytosol, chloroplast and mitochondria, plant aaRSs appeared important for plant growth and development, and disruption of their function was shown to be either lethal or cause severe defects early in plant development (Berg et al. 2005, Kim et al. 2005, Kitagawa et al. 2019, Zheng et al. 2019). Participation of aaRSs in plant stress responses is poorly characterized. A few reports have indicated aaRS involvement in plant responses to stress, but the exact molecular mechanisms were not revealed. Wheat MetRS was shown to be highly expressed during biotic stress imposed by the fungal pathogen Fusarium graminearum resulting in significant resistance of the plant to this pathogen (Zuo et al. 2016). A recent genetic screen identified aspartyl-tRNA synthetase (AspRS) as a new master regulator of resistance to pathogens induced by beta-aminobutyric acid (BABA). This non-proteinaceous amino acid is a well-known priming agent that protects many plants from various biotic and abiotic stresses. BABA interferes with AspRS canonical activity resulting in the activation of cellular defense mechanisms (Luna et al. 2014). High throughput proteomic screen for proteins involved in the early responses of Arabidopsis thaliana L. to oxidation stress due to cadmium exposure identified, among other proteins, several aminoacyl-tRNA synthetases (Sarry et al. 2006).

Prompted by findings that various types of stress significantly decrease global translation rates (Merchante et al. 2017, Zhao et al. 2019) we sought to determine stress related expression levels of aaRSs, the crucial enzymes in protein biosynthesis. We have analyzed promoters of genes encoding cytosolic seryl-tRNA synthetase (SerRS), cytosolic AspRS and cytosolic cysteinyl-tRNA synthetase (CysRS), and examined SerRS, AspRS and CysRS gene expression in A. thaliana seedlings exposed to different abiotic stressors for 3, 6, 8 and 24 hours. Although stress induces global translation repression, our results show that expression of plant aaRS genes is not decreased. Moreover, during exposure to stress conditions we detected increased AspRS and CysRS transcript levels, while SerRS gene expression did not change. This indicates that aaRSs are required for translation of specific stress-related mRNAs, as well as for resumption of translation after stress. Furthermore, plant aaRSs may participate directly in plant stress responses through their noncanonical activities.

Materials and methods

Analysis of transcription factor binding sites in promoter region

The occurrence of transcription factor (TF) binding site motifs in the promoter of a gene may be used to predict whether the gene is regulated by a specific TF. Here we analyzed promoter sequences of genes for cytosolic SerRS (At5g27470), cytosolic AspRS (At4g31180) and cytosolic CysRS (At5g38830). EPD, AGRIS AtcisDB and AthaMap databases were used for promoter analysis and detection of TF binding sites (Yilmaz et al. 2011, Hehl et al. 2016, Dreos et al. 2017). The databases are suitable for a comprehensive analysis of genomic sequences for potential TF binding sites and provide information about experimentally validated cis-regulatory elements, as well as predicted motifs present in a specific promoter.

Plant material and stress conditions

Wild-type seeds of Arabidopsis thaliana (ecotype Columbia) were surface-sterilized for one minute in 70% (v/v) ethanol and then 10 minutes in 1% (w/v) Izosan (Pliva, Croatia), rinsed 5 times with distilled water. Seeds were planted on half-strength Murashige and Skoog (MS5519, Sigma, USA), 0.5% (w/v) sucrose and 0.8% (w/v) agar medium. After three days of stratification, seeds were grown for eight days in a growth chamber under 16 h light / 8 h dark conditions (120 to 130 μmol m-2 s-1) at 24 °C. Half of the eight-day-old seedlings were planted on the same medium as described above for control and on the medium with added stressor (200 μM CdCl2, 200 mM mannitol or 200 mM NaCl). After 8 h and 24 h the seedlings were collected and snap-frozen. The other half of the eight-day-old seedlings were placed in liquid half-strength Murashige and Skoog medium with 0.5% (w/v) sucrose as control and in liquid half-strength Murashige and Skoog medium with 0.5% (w/v) sucrose and added stressor (200 μM CdCl2, 200 mM mannitol or 200 mM NaCl). After 3 h and 6 h the seedlings were collected and snap-frozen. To evaluate stressor effect on seedlings we applied short term exposure in liquid medium, and long term exposure on solidified medium. Treatment in liquid medium is more reliable and reproducible since seedlings are completely immersed in the stressor. However, long term incubation of A. thaliana seedlings in a liquid medium per se triggers stress response. To avoid this stress-related effect of liquid medium, for any longer exposure of seedlings to stressors we used a solidified medium. For each stressor, the entire experiments were performed three times using different batches of seedlings (three biological replicates).

RNA isolation, quality control and cDNA synthesis

To isolate the total RNA, 50 mg of plant material was frozen with liquid nitrogen and homogenized with glass beads. Isolation of total RNA using RNAzol (Sigma-Aldrich, USA) was performed according to the manufacturer’s protocol. Total RNA was resuspended in 50 μL of DEPC-treated water. The integrity of RNA was checked by measuring the absorbance using NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies) and using agarose gel-electrophoresis. DNA was removed with TURBO DNA-freeTM Kit (Ambion) according to the manufacturer’s protocol. Up to 12.5 µL of total RNA was incubated with 2 U of TURBO DNase for 30 minutes at 37 °C in 1x TURBO DNase buffer, followed by inactivation of TURBO DNase. The integrity of DNase-treated RNA was checked by measuring the absorbance using a NanoDrop ND-1000 spectrophotometer and using agarose gel-electrophoresis. cDNA was synthesised using High Capacity cDNA Reverse Transcription Kit (Ambion) according to the manufacturer’s protocol. Reverse transcription was performed with 1 µg of DNase-treated total RNA in a thermal cycler (Eppendorf) under the following conditions: 10 minutes at 25 °C, 120 minutes at 37 °C and 5 minutes at 85 °C.

Real-time PCR and data analysis

Gene-specific primers were designed using NCBI Primer-BLAST (Ye et al. 2012) and are listed in Tab. 1. Quantitative PCR was performed in an optical 96-well plate with 7500 Fast Real-Time PCR system (Applied Biosystems) and universal cycling conditions (10 min 95 °C, 40 cycles of 15 s at 95 °C and 60 s at 60 °C) followed by the generation of a dissociation curve to check for specificity of amplification. Reactions contained SYBR Green Master Mix (Applied Biosystems), 300 nM of a gene-specific forward and reverse primer and 2 μL of 5x diluted cDNA in each 20 μL reaction. No template controls (NTC) contained 2 μL RNase-free water instead. In qPCR, samples from various biological replicates were run in triplicate (technical replicates). PCR amplification efficiencies were calculated using LinRegPCR software (Ruijter et al. 2009). PCR amplification efficiencies for reference gene and genes of interest were comparable and are shown in Tab. 1. The relative quantification of gene expression was calculated using the comparative ΔΔCT method with actin 2 gene (At3g18780) as reference gene. The relative expression of specific genes was calculated by the following formulas: ΔCT=CT specific gene - CT reference gene and ΔΔCT=ΔCT treated sample - ΔCT untreated sample. The relative expression level was calculated as 2-ΔΔCT.

Gene Strand Sequence PCR efficiency
ACT2(At3g18780) Forward 5' - CTTGCACCAAGCAGCATGAA - 3' 1.886 ± 0.007
Reverse 5' - CCGATCCAGACACTGTACTTCCTT - 3'
SerRS(At5g27470) Forward 5' - AGCCCGTAGTTGCTGATACC - 3' 1.898 ± 0.008
Reverse 5' - AATTTCAAGAAAACAGAAGAGTCGT - 3'
AspRS(At4g31180) Forward 5' - TCCCAGAAGTCTTGGAGCAAC - 3' 1.875 ± 0.014
Reverse 5' - CAAATCCACCGTGTAGAGGC - 3'
CysRS(At5g38830) Forward 5' - CGCAGCTAGAGAGTTCGTCA - 3' 1.888 ± 0.004
Reverse 5' - TCCACCATCTTCAGACAATGC - 3'
Tab. 1.Primer pairs used for quantitative real-time PCR. ACT2 primer sequences were taken from Czechowski et al. 2005 and Remans et al. 2008.

Statistical analysis

Data analysis was conducted with Student’s t-test to detect the differences between control and treated groups (Yuan et al. 2006). Standard deviation (SD) was calculated with the use of the results of the 2−ΔΔCT method. The results are presented as mean ± SD from three independent experiments. Differences were considered significant at P < 0.05. Asterisk symbols (*) indicate significant differences: **P < 0.01 and *P < 0.05.

Results

We examined genomic regions upstream of the genes encoding cytoslic SerRS, AspRS and CysRS using the A. thaliana cis-regulatory element database AtcisDB, which considers a promoter to be an intergenic region upstream of the gene of interest, excluding any coding region of upstream genes (Yilmaz et al. 2011). According to the AtcisDB database, upstream intergenic regions of all three aaRS genes were identified as predicted promoters, since no experimentally validated binding of TFs to these regions is documented so far. Analysis using AtcisDB revealed that the predicted promoter of SerRS gene is only 194 base pairs long and that transcription factors of WRKY family are predicted to bind to SerRS promoter and regulate transcription of SerRS gene. The promoter sequence of AspRS gene is predicted to be 1039 bp long, with putative binding sites for transcription factor families BHLH, Homeobox, bZIPs and ABI3VP1, while CysRS promoter is the longest (2391 bp) and contains binding motifs for transcription factor families BHLH, Homeobox, WRKY, MADS, MYB, bZIP, ABI3VP1 and LFY. Furthermore, using AthaMap database (Hehl et al. 2016) we examined the region from -1,000 bp to +200 bp relative to the transcription start of a gene considering that, in general, majority of gene regulatory sequences (86%) are found inside this area (Yu et al. 2016). Analysis showed that this region in SerRS, AspRS and CysRS genes contains many potential places for binding of TFs belonging to major stress-related TF protein families AP2/ERF, NAC, MYB, bZIP and WRKY (Baillo et al. 2019) (Tab. 2). Altogether, bioinformatics analysis showed that promotors of three aaRS genes contain many potential binding sites for stress-related TFs, implying that abiotic stress may influence expression of the aaRS genes examined.

TF family SerRS AspRS CysRS TF family function (according to PlnTFDB)
AP2/ERF 88 30 46 plant-specific transcription factors, key regulators of floral organ identity determination and of leaf epidermal cell identity, key regulators of several abiotic stresses and respond to multiple hormones
NAC 68 51 47 plant-specific transcription factors, key regulators of stress perception and developmental programmes
MYB 56 24 43 involved in the control of the cell cycle, involved in control of plant secondary metabolism, as well as the identity and fate of plant cells
bZIP 21 1 8 regulate processes including pathogen defence, light and stress signalling, seed maturation and flower development
WRKY 3 1 0 play significant roles in responses to biotic and abiotic stresses, and in development
Tab. 2.Number of potential transcription factor (TF) binding sites in region -1000 bp to +200 bp from transcription start site of SerRS, AspRS and CysRS genes. Analysis was performed using AthaMap and plant transcription factor database (PlnTFDB, Pérez-Rodríguez et al. 2010).

In order to determine the impact of abiotic stress conditions on the expression of selected aaRS genes, seedlings of A. thaliana were exposed to osmotic stress imposed with 200 mM mannitol, salt stress imposed with 200 mM NaCl or heavy metal stress imposed with 200 µM cadmium. Expression of SerRS, AspRS and CysRS genes was measured at four different time points. The transcript levels in control and treated samples were measured using RT-qPCR and normalized by using the expression of actin2 gene as endogenous control. The effect of 200 mM mannitol is shown in Fig. 1. The SerRS gene expression did not differ significantly from SerRS gene expression in control, untreated seedlings (Fig. 1A). On the other hand, in certain conditions osmotic stress significantly enhanced AspRS and CysRS gene expression. Expression of AspRS gene was significantly elevated 8 hours after osmotic stress application (Fig. 1B), while CysRS gene expression rose immediately after mannitol application and remained significantly increased for 3, 6 and 8 hours (Fig. 1C). The effect of 200 mM NaCl is shown in Fig. 2. The SerRS gene expression did not significantly differ from that in controls (Fig. 2A), while the expression level of the AspRS gene significantly increased 8 hours after addition of NaCl and stayed induced even 24 hours following salt application (Fig. 2B). Salt stress imposed with 200 mM NaCl significantly enhanced CysRS expression 6 and 8 hours after application, while after prolonged exposure (24 hours) CysRS gene expression decreased (Fig. 2C).

Fig. 1.Expression of cytosolic seryl-tRNA synthetase (SerRS), cytosolic aspartyl-tRNA synthetase (AspRS) and cytosolic cysteinyl-tRNA synthetase (CysRS) under osmotic stress induced by 200 mM mannitol. The expression level for each gene in stress condition was calculated relative to its expression in the control, untreated sample. The gene expression was normalized with the housekeeping actin 2 transcript. Bars represent mean values ± SD from independent experiments (n = 3). Asterix indicate significant difference between treated and control samples, *P < 0.05, **P < 0.01.

Fig. 2.Expression of cytosolic seryl-tRNA synthetase (SerRS), cytosolic aspartyl-tRNA synthetase (AspRS) and cytosolic cysteinyl-tRNA synthetase (CysRS) under salt stress induced by 200 mM NaCl. The expression level for each gene in stress condition was calculated relative to its expression in the control, untreated sample. The gene expression was normalized with the housekeeping actin 2 transcript. Bars represent mean values ± SD from independent experiments (n = 3). Asterix indicate significant difference between treated and control samples, *P < 0.05.

Protein levels of SerRS, AspRS and CysRS were shown to be increased during early responses of A. thaliana due to cadmium exposure (Sarry et al. 2006). In order to examine whether accumulation of these proteins is a result of enhanced transcript accumulation we tested the gene expression of selected aaRSs after exposure to cadmium for 3, 6, 8 and 24 hours. The effect of 200 µM cadmium on aaRSs gene expression in the A. thaliana seedlings is shown in Fig. 3. Compared with control seedlings, the treated seedlings did not show a statistically significant change in SerRS transcript levels during different exposure times (Fig. 3A). The same results were obtained during analysis of AspRS and CysRS transcript levels (Fig. 3B, C).

Fig. 3.Expression of cytosolic seryl-tRNA synthetase (SerRS), cytosolic aspartyl-tRNA synthetase (AspRS) and cytosolic cysteinyl-tRNA synthetase (CysRS) under heavy metal stress induced by 200 µM cadmium. The expression level for each gene in stress condition was calculated relative to its expression in the control, untreated sample. The gene expression was normalized with the housekeeping actin 2 transcript. Bars represent mean values ± SD from independent experiments (n = 3).

Discussion

Land plants are anchored in one place for most of their life cycle and therefore must constantly be able to adapt their growth and metabolism to abiotic stresses such as changes in light intensity and temperature, high salinity conditions, heavy metal exposure and the lack of water and essential minerals (Floris et al. 2009). Transcription factors are key regulators of gene expression that control cellular and developmental processes and the response to environmental stimuli. The number of transcription factor genes in A. thaliana genome is >1500 which is 3 times more than in animals with similar genome size (Drosophila melanogaster (Diptera) and Caenorhabditis elegans (nematode); Riechmann et al. 2000), which suggests that transcriptional regulation plays a more important role in plants then in animals. Hence, plants’ survival depends on their ability to rapidly regulate gene expression in order to adapt their physiology to their environment. To gain more insight into aaRSs in stress response we analyzed promoters of genes encoding cytosolic SerRS, AspRS and CysRS. The rationale was to reveal whether promoters of these genes contain binding site motifs for the transcription factors important in plant stress response. In addition, we determined stress related expression levels of selected aaRS genes.

We exposed A. thaliana seedlings to osmotic, salt and heavy metal stress and analyzed transcript expression patterns of selected aaRS genes (SerRS, AspRS and CysRS) 3, 6, 8 and 24 hours after the application of various stressors. Our results show that expression of none of the investigated aaRS genes changes during stress induced by exposure to cadmium (Fig. 3). This observation is surprising considering that previous proteomic investigation showed that protein levels of all three investigated aaRSs are increased during early responses of A. thaliana cells to oxidation stress due to cadmium exposure (Sarry et al. 2006). The discrepancy between previous proteomic data and our gene expression study may be related to the increased stability of aaRS proteins or mRNAs, as well as increased translation of aaRS mRNAs on ribosomes during stress. It is known that increased protein levels are not always correlated to transcript levels and selective mRNA translation is documented especially during exposure to stress (Sablok et al. 2017). It is interesting to note that results of the mentioned proteomic investigation indicated that SerRS undergoes posttranslational modification during cadmium stress. We have recently solved crystal structure of A. thaliana SerRS and identified the disulfide link in this cytosolic protein (Kekez et al. 2019). Considering that cadmium induces oxidative stress in cells, it is plausible that cysteines in SerRS are oxidized and form a disulfide link, thus creating a new, posttranslationally modified form of SerRS. This oxidized form of SerRS may possess higher stability due to introduction of a covalent disulfide link and may participate in plant stress response mechanisms. In this work we did not observe a statistically significant change in SerRS transcript levels during exposure to osmotic, salt and heavy metal stress, although promoter analysis revealed many potential binding sites for stress related transcription factors in the region from -1,000 bp to +200 bp relative to the transcription start site (Tab. 2). Although this region is regulatory for many A. thaliana genes (Yu et al. 2016), it appears that regulation of SerRS gene is controlled by its short promoter (196 bp) with only few putative TF binding sites (as predicted with AtcisDB database) and that expression of SerRS gene is not responsive to stress conditions used in this work. However, as explained above, SerRS expression in stress can be regulated posttranscriptionally, at the mRNA or protein level.

The expression level of the AspRS gene significantly increased during osmotic and salt stress, induced by the application of 200 mM mannitol or 200 mM NaCl, respectively. The increase in AspRS gene expression was observed 8 hours after stress application, and increased levels of AspRS mRNA were present even after 24 hours in the case of salt stress. The observed increased AspRS gene expression is in agreement with our detailed analysis of AspRS promoter which indicated 51, 30 and 24 binding site motifs for TFs belonging to NAC, AP2/ERF and MYB families, respectively (Tab. 2). These TF families are important regulators of stress response (Nuruzzaman at al. 2013, Baillo et al. 2019). It was previously shown that transcription of A. thaliana AspRS gene was enhanced during biotic stress imposed by pathogen Hyaloperonospora arabidopsidis (Peronosporaceae) infection and that the enzyme can switch from canonical AspRS activity to noncanonical defense activity upon pathogen infection (Luna et al. 2014). Here we show that the level of AspRS transcripts increases during osmotic and salt stress indicating that AspRS also participates in abiotic stress responses. Therefore, it appears that AspRS can serve as a general stress factor in plants.

The expression level of the CysRS gene significantly increased for 3, 6 and 8 hours after the application of 200 mM mannitol and with the application of 200 mM NaCl in duration for 6 and 8 hours. Our analysis showed that the promoter of CysRS gene contains 47, 46 and 43 possible binding sites for NAC, AP2/ERF and MYB transcription factors, respectively (Tab. 2). This high number of stress-related TF binding sites and large promoter size (2391 bp) is in accordance with prompt and enhanced CysRS gene expression after exposure to stress. It is interesting to note that an increase in free cysteine levels has been reported as a response to various abiotic stress factors, including cadmium, salt and drought stress (Zagorchev et al. 2013). Cysteine is mainly needed for the biosynthesis of sulfur-rich compounds with anti-stress activity, such as GSH, phytochelatins (oligomers of GSH) and stress-related proteins (Zagorchev et al. 2013). If cysteine is excessively used for synthesis of anti-stress compounds, cysteine levels are presumably not sufficient to support translation and therefore CysRS is overproduced to mitigate that effect. On the other hand, recent studies have shown that bacterial and human CysRS use cysteine to produce persulfides, such as CysSSH or longer chain sulfur compounds (polysulfides, including CysS/GS–(S)n–H) (Akaike et al. 2017). These persulfides have anti-oxidative effects. It is plausible to assume that plant CysRS participates in similar processes through its noncanonical function.

Although stress triggers a rapid downregulation of general protein biosynthesis (Merchante et al. 2017), our results show that expression of plant aaRSs that participate in translation is not decreased. Moreover, during exposure to stress conditions we detected increased AspRS and CysRS transcript levels, while others reported increased SerRS, AspRS and CysRS protein levels (Sarry et al. 2006). Considering that during stress, translation of stress related proteins is enhanced we may assume that the levels of aaRSs are kept at steady-state to pursue translation of specific mRNAs encoding proteins that are essential for stress resistance and survival. Also, when stress conditions cease and global translation starts again, it may be advantageous to have sufficient levels of aaRSs to sustain translation. Finally, increased protein levels of plant aaRSs during stress may serve as a pool of aaRS proteins that can participate directly in stress responses through their noncanonical activities. Further research on plant aaRS role in stress is required to reveal the underlying molecular mechanisms that link aaRSs to plant stress response.

References

  1. Akaike T, Ida T, Wei FY, Nishida M, Kumagai Y, Alam MM. Cysteinyl-tRNA synthetase governs cysteine polysulfidation and mitochondrial bioenergetics.. Nat Commun. 2017; 8:1177-92. DOI | PubMed
  2. Aldinger CA, Leisinger AK, Igloi GL. The influence of identity elements on the aminoacylation of tRNA(Arg) by plant and Escherichia coli arginyl-tRNA synthetases.. FEBS J. 2012; 279:3622-38. DOI | PubMed
  3. Baillo EH, Kimotho RN, Zhang Z, Xu P. Transcription factors associated with abiotic and biotic stress tolerance and their potential for crops improvement.. Genes (Basel). 2019; 10:771. DOI | PubMed
  4. Berg M, Rogers R, Muralla R, Meinke D. Requirement of aminoacyl-tRNA synthetases for gametogenesis and embryo development in Arabidopsis.. Plant J. 2005; 44:866-78. DOI | PubMed
  5. Cao X, Lia C, Xiaoa S, Tanga Y, Huanga J, Zhaoa S. Acetylation promotes TyrRS nuclear translocation to prevent oxidative damage.. Proc Natl Acad Sci USA. 2017; 114:687-92. DOI | PubMed
  6. Cramer GR, Urano K, Delrot S, Pezzotti M, Shinozaki K. Effects of abiotic stress on plants: a systems biology perspective.. BMC Plant Biol. 2011; 11:163. DOI | PubMed
  7. Czechowski T, Stitt M, Altmann T, Udvardi MK, Scheible W-R. Genome-wide identification and testing of superior reference genes for transcript normalization in Arabidopsis.. Plant Physiol. 2005; 139:5-17. DOI | PubMed
  8. Dreos R, Ambrosini G, Groux R, Périer RC, Bucher P. The eukaryotic promoter database in its 30th year: focus on non-vertebrate organisms.. Nucleic Acids Res. 2017; 45:D51-5. DOI | PubMed
  9. Duchêne AM, Pujol C, Maréchal-Drouard L. Import of tRNAs and aminoacyl-tRNA synthetases into mitochondria.. Curr Genet. 2009; 55:1-18. DOI | PubMed
  10. Floris M, Mahgoub H, Lanet E, Robaglia C, Menand B. Post-transcriptional regulation of gene expression in plants during abiotic stress.. Int J Mol Sci. 2009; 10:3168-85. DOI | PubMed
  11. Guo M, Schimmel P. Essential nontranslational functions of tRNA synthetases.. Nat Chem Biol. 2013; 9:145-53. DOI | PubMed
  12. Hehl R, Norval L, Romanov A, Bülow L. Boosting AthaMap database content with data from protein binding microarrays.. Plant Cell Physiol. 2016; 57DOI | PubMed
  13. Hoffman KS, Vargas-Rodriguez O, Bak DW, Mukai T, Woodward LK, Weerapana E. A cysteinyl-tRNA synthetase variant confers resistance against selenite toxicity and decreases selenocysteine misincorporation.. J Biol Chem. 2019; 294:12855-65. DOI | PubMed
  14. Kekez M, Bauer N, Saric E, Rokov-Plavec J. Exclusive cytosolic localization and broad tRNASer specificity of Arabidopsis thaliana seryl-tRNA synthetase.. J Plant Biol. 2016; 59:44-54. DOI
  15. Kekez M, Zanki V, Kekez I, Baranasic J, Hodnik V, Duchêne AM. Arabidopsis seryl‐tRNA synthetase: the first crystal structure and novel protein interactor of plant aminoacyl‐tRNA synthetase.. FEBS J. 2019; 286:536-54. DOI | PubMed
  16. Kim Y-K, Lee J-Y, Cho HS, Lee SS, Ha HJ, Kim S. Inactivation of organellar glutamyl- and seryl-tRNA synthetases leads to developmental arrest of chloroplasts and mitochondria in higher plants.. J Biol Chem. 2005; 280:37098-106. DOI | PubMed
  17. Kitagawa M, Balkunde R, Bui H, Jackson D. An aminoacyl tRNA synthetase, OKI1, is required for proper shoot meristem size in Arabidopsis.. Plant Cell Physiol. 2019; 60:2597-608. DOI | PubMed
  18. Lee J, Joshi N, Pasini R, Dobson RC, Allison J, Leustek T. Inhibition of Arabidopsis growth by the allelopathic compound azetidine-2-carboxylate is due to the low amino acid specificity of cytosolic prolyl-tRNA synthetase.. Plant J. 2016; 88:236-46. DOI | PubMed
  19. Lee JY, Kim DG, Kim BG, Yang WS, Hong J, Kang T. Promiscuous methionyl-tRNA synthetase mediates adaptive mistranslation to protect cells against oxidative stress.. J Cell Sci. 2014; 127:4234-45. DOI | PubMed
  20. Liu B, Qian S-B. Translational reprogramming in stress response.. WIREs RNA. 2014; 5:301-5. DOI | PubMed
  21. Luna E, van Hulten M, Zhang Y, Berkowitz O, López A, Pétriacq P. Plant perception of b-aminobutyric acid is mediated by an aspartyl-tRNA synthetase.. Nat Chem Biol. 2014; 10:450-6. DOI | PubMed
  22. Merchante C, Stepanova AN, Alonso JM. Translation regulation in plants: an interesting past, an exciting present and a promising future.. Plant J. 2017; 90:628-53. DOI | PubMed
  23. Mocibob M, Rokov-Plavec J, Godinic-Mikulcic V, Gruic-Sovulj I. Seryl-tRNA synthetases in translation and beyond.. Croat Chem Acta. 2016; 89:261-76. DOI
  24. Nuruzzaman M, Sharoni AM, Kikuchi S. Roles of NAC transcription factors in the regulation of biotic and abiotic stress responses in plants.. Front Microbiol. 2013; 4:248. DOI | PubMed
  25. Pérez-Rodríguez P, Riano-Pachon DM, Correa LGG, Rensing SA, Kersten B, Mueller-Roeber B. PlnTFDB: updated content and new features of the plant transcription factor database.. Nucleic Acids Res. 2010; 38:D822-7. DOI | PubMed
  26. Perona JJ, Gruic-Sovulj I. Synthetic and editing mechanisms of aminoacyl-tRNA synthetases.. Top Curr Chem. 2014; 344:1-41. PubMed
  27. Qin F, Shinozaki K, Yamaguchi-Shinozaki K. Achievements and challenges in understanding plant abiotic stress responses and tolerance.. Plant Cell Physiol. 2011; 52:1569-82. DOI | PubMed
  28. Remans T, Smeets K, Opdenakker K, Mathijsen D, Vangronsveld J, Cuypers A. Normalisation of real-time RT-PCR gene expression measurements in Arabidopsis thaliana exposed to increased metal concentrations.. Planta. 2008; 227:1343-9. DOI | PubMed
  29. Riechmann JL, Heard J, Martin G, Reuber L, Jiang C, Keddie J. Arabidopsis transcription factors: genome-wide comparative analysis among eukaryotes.. Science. 2000; 290:2105-10. DOI | PubMed
  30. Rokov J, Söll D, Weygand-Durasevic I. Maize mitochondrial seryl-tRNA synthetase recognizes Escherichia coli tRNA(Ser) in vivo and in vitro.. Plant Mol Biol. 1998; 38:497-502. DOI | PubMed
  31. Rokov J, Weygand-Durasevic I. Seryl-tRNA synthesis in maize.. Period Biol. 1999; 1999:137-42.
  32. Rokov-Plavec J, Bilokapic S, Gruic-Sovulj I, Mocibob M, Glavan F, Brgles M. Unilateral flexibility in tRNASer recognition by heterologuous seryl-tRNA synthetase.. Period Biol. 2004; 106:147-54.
  33. Rokov-Plavec J, Dulic M, Duchêne AM, Weygand-Durasevic I. Dual targeting of organellar seryl-tRNA synthetase to maize mitochondria and chloroplasts.. Plant Cell Rep. 2008; 27:1157-68. DOI | PubMed
  34. Rokov-Plavec J, Lesjak S, Gruic-Sovulj I, Mocibob M, Dulic M, Weygand-Durasevic I. Substrate recognition and fidelity of maize seryl-tRNA synthetases.. Arch Biochem Biophys. 2013; 529:122-30. DOI | PubMed
  35. Rokov-Plavec J, Lesjak S, Landeka I, Mijakovic I, Weygand-Durasevic I. Maize seryl-tRNA synthetase: specificity of substrate recognition by the organellar enzyme.. Arch Biochem Biophys. 2002; 397:40-50. DOI | PubMed
  36. Ruijter JM, Ramakers C, Hoogaars WM, Karlen Y, Bakker O, van den Hoff MJ. Amplification efficiency: linking baseline and bias in the analysis of quantitative PCR data.. Nucleic Acids Res. 2009; 37DOI | PubMed
  37. Sablok G, Powell JJ, Kazan K. Emerging roles and landscape of translating mRNAs in plants.. Front Plant Sci. 2017; 8:1443. DOI | PubMed
  38. Sarry JE, Kuhn L, Ducruix C, Lafaye A, Junot C, Hugouvieux V. The early responses of Arabidopsis thaliana cells to cadmium exposure explored by protein and metabolite profiling analyses.. Proteomics. 2006; 6:2180-98. DOI | PubMed
  39. Steiner RE, Kylec AM, Ibbaa M. Oxidation of phenylalanyl-tRNA synthetase positively regulates translational quality control.. Proc Natl Acad Sci USA. 2019; 116:10058-63. DOI | PubMed
  40. Wei N, Shi Y, Truong LN, Fisch KM, Xu T, Gardiner E. Oxidative stress diverts tRNA synthetase to nucleus for protection against DNA damage.. Mol Cell. 2014; 56:323-32. DOI | PubMed
  41. Yang X, Li G, Tian Y, Song Y, Liang W, Zhang D. A rice glutamyl-tRNA synthetase modulates early anther cell division and patterning.. Plant Physiol. 2018; 177:728-44. DOI | PubMed
  42. Yángüez E, Castro-Sanz AB, Fernández-Bautista N, Oliveros JC, Castellano MM. Analysis of genome-wide changes in the translatome of Arabidopsis seedlings subjected to heat stress.. PLoS One. 2013; 8DOI | PubMed
  43. Ye J, Coulouris G, Zaretskaya I, Cutcutache I, Rozen S, Madden T. Primer-BLAST: A tool to design target-specific primers for polymerase chain reaction.. BMC Bioinformatics. 2012; 13:134. DOI | PubMed
  44. Yilmaz A, Mejia-Guerra MK, Kurz K, Liang X, Welch L, Grotewold E. AGRIS: Arabidopsis gene regulatory information server, an update.. Nucleic Acids Res. 2011; 39(Database issue):D1118-22. DOI | PubMed
  45. Yu CP, Lin JJ, Li WH. Positional distribution of transcription factor binding sites in Arabidopsis thaliana. Sci Rep. 2016; 6:25164. DOI | PubMed
  46. Yuan JS, Reed A, Chen F, Stewart CN. Statistical analysis of real-time PCR dana.. BMC Bioinformatics. 2006; 7:85. DOI | PubMed
  47. Zagorchev L, Seal CE, Kranner I, Odjakova M. A central role for thiols in plant tolerance to abiotic stress.. Int J Mol Sci. 2013; 14:7405-32. DOI | PubMed
  48. Zhao J, Qin B, Nikolay R, Spahn CMT, Zhang G. Translatomics: The global view of translation.. Int J Mol Sci. 2019; 20:1-23. PubMed
  49. Zheng H, Wang Z, Tian Y, Liu L, Lv F, Kong W. Rice albino 1, encoding a glycyl-tRNA synthetase, is involved in chloroplast development and establishment of the plastidic ribosome system in rice.. Plant Physiol Biochem. 2019; 139:495-503. DOI | PubMed
  50. Zuo D-Y, Yi S-Y, Liu R-J, Qu B, Huang T, He W-J. A deoxynivalenol-activated methionyl-tRNA synthetase gene from wheat encodes a nuclear localized protein and protects plants against Fusarium pathogens and mycotoxins.. Phytopathology. 2016; 106:614-23. DOI | PubMed