Tolerance of Transplastomic Tobacco Plants Overexpressing a Theta Class Glutathione Transferase to Abiotic and Oxidative Stresses

Chloroplasts are organelles subjected to extreme oxidative stress conditions. Biomolecules produced in the chloroplasts act as signals guiding plant metabolism toward stress tolerance and play a major role in regulating gene expression in the nucleus. Herein, we used transplastomic plants as an alternative approach to expression of transgenes in the nucleus for conferring stress tolerance to abiotic stresses and herbicides. To investigate the morphophysiological and molecular mechanisms and the role of plastid expressed GSTs in tobacco stress detoxification and stress tolerance, we used transplastomic tobacco lines overexpressing a theta class glutathione transferase (GST) in chloroplasts. The transplastomic plants were tested under drought (0, 100, and 200 mM mannitol) and salinity (0, 150, and 300 mM NaCl) in vitro, and under herbicide stress (Diquat). Our results suggest that pt^AtGSTT lines were tolerant to herbicide-induced oxidative and salinity stresses and showed enhanced response tolerance to mannitol-induced osmotic stress compared to WT plants. Overexpression of the Arabidopsis thaliana AtGSTT in the chloroplasts resulted in enhanced photo-tolerance and turgor maintenance under stress. Whole-genome transcriptome analysis revealed that genes related to stress tolerance, were upregulated in pt^AtGSTT2a line under both control and high mannitol stress conditions. Transplastomic plants overexpressing the pt^AtGSTT2a in the chloroplast showed a state of acclimation to stress, as only limited number of genes were upregulated in the pt^AtGSTT2a transplastomic line compared to WT under stress conditions while at the same time genes related to stress tolerance were upregulated in pt^AtGSTT2a plants compared to WT in stress-free conditions. In parallel, the metabolic profile indicated limited perturbations of the metabolic homeostasis in the transplastomic lines and greater accumulation of mannitol, and soluble sugars under high mannitol stress. Therefore, transplastomic lines seem to be in a state of acclimation to stress under stress-free conditions, which was maintained even under high mannitol stress. The results help to elucidate the role of GSTs in plant abiotic stress tolerance and the underlying mechanisms of the GSTs expressed in the chloroplast, toward environmental resilience of cultivated crops.

Introduction

Developing crop plants, able to yield better under abiotic stresses or plants with multiple herbicide resistance, is a prerequisite for improved crop production. The chloroplast, abundant in plant cells and eukaryotic algae, is the site of photosynthesis, providing the primary source of the world's food productivity (Verma and Daniell, 2007). As chloroplasts are the organelles responsible for photosynthesis they are also a source of reactive oxygen species (ROS) in plants (Foyer and Shigeoka, 2011). Furthermore, environmental stresses have been found to produce an excess of excitation energy in chloroplasts, resulting in the production of ROS, thus they are also considered to be implicated in the regulation of stress responses or even act as a sensor of cellular stress (Mullineaux and Karpinski, 2002). Genetic transformation of chloroplasts has been used as an alternative approach to the expression of transgenes in the nucleus (Wang et al., 2009). The transplastomic system has three main advantages: (i) prevents gene flow via pollen through transgene containment due to maternal inheritance, (ii) has highly active chloroplast transcription and translation machineries, and (iii) a lack of epigenetic interference allows stable transgene expression (Bock, 2014). Chloroplast engineering has been applied for the development of resistant crops to various abiotic and biotic stresses (Clarke and Daniell, 2011), production of biopharmaceuticals, metabolic pathway engineering and advances on RNA editing (reviewed in Wang et al., 2009) and phytoremediation (reviewed in Verma and Daniell, 2007).

Understanding the adaptation of plants to different climatic conditions, such as high temperatures, water logging, and drought is essential for addressing climate change challenges. Improving the resilience of chloroplasts through plastid engineering may provide a solution toward the improvement of crop productivity (Clarke and Daniell, 2011). To date, there are a limited number of studies regarding the development of transplastomic plants and their response to abiotic stress. Transplastomic tobacco plants expressing a choline monooxygenase (BvCMO) from Beta vulgaris demonstrated increased tolerance to salt (100 and 150 mM NaCl) and drought (300 mM mannitol) stresses (Zhang et al., 2008). Genetic engineering of carrot chloroplast genome expressing the Betaine-aldehyde dehydrogenase (badh) gene also improved tolerance to high salinity (400 mM L⁻¹ NaCl) (Kumar et al., 2004). Similarly, transplastomic Nicotiana benthamiana plants expressing multiple defense genes encoding protease inhibitors and chitinase were more tolerant to 200 mM NaCl and 3% PEG compared to the wild type plants and were able to maintain greater root growth activity due to transgene expression in the leucoplasts of roots (Chen et al., 2014). Transplastomic tobacco lines overexpressing an A. thaliana γ-tocopherol methyltransferase (Atγ-tmt) gene accumulated higher levels of α-tocopherol when grown in 400 mM NaCl, compared to wild-type plants, which accumulated higher starch and total soluble sugars, but transplastomic plants better regulated sugar transport (Jin and Daniell, 2014). Genetically engineered plastomes have provided a generation of herbicide-tolerant plants demonstrated in tobacco for tolerance to glyphosate (Ye et al., 2001; Chin et al., 2003), phosphinothricin (Iamtham and Day, 2000; Lutz et al., 2001) sulcotrione (Falk et al., 2005), isoxaflutole (IFT) (Dufourmantel et al., 2007) and paraquat (methyl-viologen) (Poage et al., 2011; Chen et al., 2014).

Plant glutathione S-transferases (GSTs) have been shown to modulate redox homeostasis by alterations in GSH content and redox state (Sappl et al., 2009), conferring tolerance to a wide range of abiotic stresses (Kumar et al., 2013; Csiszár et al., 2014; Kissoudis et al., 2015b; Kayum et al., 2018) including herbicides (Kissoudis et al., 2015a; Lo Cicero et al., 2015, 2017). Glutathione transferases (GSTs; EC 2.5.1.18) are a superfamily of multifunctional proteins that in plants, have evolved into six discreet groups classified as the zeta (Z), theta (T), phi (F), tau (U), lambda (L), and dehydroascorbate reductase (DHAR) classes, respectively, (Dixon and Edwards, 2010). Functions ascribed to date include the detoxification of herbicides (phi and tau), tyrosine degradation (zeta), the reduction of intermediates involved in redox cycling (DHAR and lambda), and acting as glutathione peroxidases toward organic hydroperoxides (theta). In the case of the theta enzymes (GSTTs), this ability to use glutathione to reduce organic hydroperoxides is conserved between plants and animals and is thought to be important in oxidative metabolism, most notably through the processing of phytotoxic oxidized lipids in the peroxisomes (Dixon et al., 2009; Dixon and Edwards, 2010).

GSTs have been used before in chloroplast transformation; The SjGST26 (EC:2.5.1.18), from Schistomosoma japonicum (Smith and Johnson, 1988) and His-tagged derivative of the maltose- binding protein (His₆-MBP) were expressed in tobacco chloroplasts to be used as affinity tags for the rapid purification of chloroplast-expressed proteins (Ahmad et al., 2012). Transplastomic tobacco lines overexpressing glutathione reductase (GR) alone or combined with GST were more tolerant under 10°C, whereas lines overexpressing dehydroascorbate reductase (DHAR) alone or in combination with GR were more sensitive compared to wild type plants (Grant et al., 2014). When these lines were chilled at 4°C and under relatively high photosynthetically active radiation (PAR), all lines were more sensitive compared to wild type plants, indicating that overexpression of the ROS-scavenging enzymes may be dependent on the interaction of light and cold stress (Grant et al., 2014). Transplastomic seedlings expressing either DHAR or an Escherichia coli GST B1-1, which has been shown to exhibit a GSH-dependent peroxidase activity against cumene hydroperoxide (Nishida et al., 1994) and proved to be important for bacterial resistance to hydrogen peroxide- induced oxidative stress (Kanai et al., 2006), or a combination of DHAR:GR and GST:GR in chloroplasts were less sensitive to salt (200 mM NaCl) and cold (4°C) compared to wild type seedlings (Le Martret et al., 2011). However, only the simultaneous expression of DHAR:GR and GST:GR conferred tolerance to methyl viologen (MV) (Le Martret et al., 2011). Transplastomic tobacco lines expressing GR in combination with either DHAR or GST (from E. coli) exhibited better tolerance to supplemental UV-B than wild type plants (Czégény et al., 2016). The expression of GSTs in compartments where they are not normally found in, can reveal new insights into their functions. For example, the expression of GSTs in cellular compartments (recombinant bacteria, plant chloroplasts) producing porphyrins has revealed their ability to bind to porphyrinogen intermediates (Dixon et al., 2008). In the case of the ZmGSTU1-ZmGSTU2, the transplastomics ability to protect plants against herbicides that inhibit porphyrin synthesis in the chloroplast shed light into the functional role of the engineered chimeric enzyme (Dixon et al., 2008).

The GSTs are predominantly not targeted for expression in the chloroplast, however, if they are expressed in this organelle, they could deliver some of their key antioxidant and detoxification functions, such as metabolizing photosystem herbicides, and reducing lipid hydroperoxides generated by ROS formed during photosynthesis. In addition to the efficiency of transplastomic expression, we were also interested in how the protective functions of GSTs could be manifested in an organelle where they are not normally targeted for expression. GSTs are important enzymes of the antioxidant pathway and when expressed in the plastome we hypothesized that the leaf physiology and performance would be enhanced under stress compared to the non-transformed wild type plants. None of the above-mentioned examples were performed with plant derived GSTs from the Theta or Tau classes. Therefore, to investigate whether the overexpression of these GSTs in the chloroplast enhances tolerance to salinity, drought, and herbicide induced oxidative stress we used T1 transplastomic tobacco lines overexpressing a theta class GST from Arabidopsis thaliana AtGSTT1 (At5g41210), an enzyme normally only expressed in the peroxisomes which is highly active as a glutathione peroxidase toward organic hydroperoxide substrates or a Zea mays tau class chimeric ZmGSTU1/ZmGSTU2 enzyme (EFD6-115A), which has been previously shown to protect the transformed plants from herbicide injury through its ability to detoxify fluorodifen (Dixon et al., 2003) and in subsequent studies it was confirmed that the chimera had the additional ability to bind porphyrinogen intermediates formed during chlorophyll biosynthesis, a trait shared with its ZmGSTU parent proteins (Dixon et al., 2008). To assess plant tolerance to abiotic stresses, we investigated the morphophysiological parameters, and the metabolic and transcriptomic reactions involved in the response of transplastomic tobacco lines. Herein, we approach plant stress tolerance from an alternative perspective via chloroplast engineering to (i) mitigate the oxidative stress imposed under various abiotic and anthropogenic stress conditions and (ii) unravel the complex networks of molecular interactions controlling plant acclimation to field conditions.

Materials and Methods

Plant Material and Experimental Design

For the experiments we used homoplastomic, transplastomic tobacco lines pt^AtGSTT and pt^EFD6−115A overexpressing the AtGSTT (lines 2a and 6-1) or a ZmGSTU1-ZmGSTU2 chimera in chloroplasts, respectively (Dixon et al., 2008). The seeds of the pt^AtGSTT and pt^EFD6−115A T1 lines were initially grown on MS selection medium supplemented with Streptomycin Sulfate (500 mg L⁻¹) and Spectinomycin Dihydrochloride (250 mg L⁻¹) (Duchefa Biochemie, The Netherlands), whereas the wild-type (WT) tobacco seeds were placed on plain MS medium. After selection, the plantlets were transferred to MS media for further growth and when they reached four true leaves were tested in vitro under drought (0, 100, and 200 mM mannitol; AppliChem-PanReac, Germany) and salinity (0, 150, and 300 mM NaCl; Centralchem, Slovakia) conditions (n = 6). The experiments lasted for 35 and 20 days, respectively.

The in vivo herbicide experiment was performed in a controlled glasshouse environment with a photoperiod of 14/10 h light/dark. The temperature was between 20 and 27°C, with a mean temperature of 23°C. Plantlets undergone acclimatization for 3 weeks and Diquat a non-selective contact herbicide, was applied as Reglone 20 SL formulation (Syngenta Hellas) at 1 and 2 L of Reglone/ha (200-low dose; Diq_L and 400 -high dose; Diq_H, g ai of diquat per hectare, respectively). Herbicide treatments were performed with a portable field plot sprayer (AZO-SPRAYERS, P.O. Box 350-6710 BJ EDE, The Netherlands) using flat-fan nozzles (Teejet Spray System Co., P.O. Box 7900, Wheaton, IL 60188) and calibrated to deliver 300 L/ha of water at 280 kPa pressure.Diquat (REGLONE^® Desiccant, Syngenta Canada Inc), Control plants were sprayed with the same volume of water only (no herbicide). All pots were placed in a randomized complete block design (transgenic lines: n = 15 and WT: n = 9). The experiment lasted for 2 days after herbicide application.

Morphophysiological Measurements

Dark-adapted chlorophyll a fluorescence measurements were performed on the youngest fully developed leaf on the adaxial leaf surface using the OS30p+ chlorophyll fluorometer (Opti-Sciences Inc., Hudson, USA) following dark adaptation of 30 min. Relative chlorophyll content was measured according to Stavridou et al. (2016) on one leaf per plant with three averaged measurements using a CCM-200 plus chlorophyll content meter (Opti-Sciences Inc., Hudson, USA). Harvested plants were separated into leaves, stems, and roots and the final morphological parameters, such as stem length, root length, number of leaves, and plant fresh matter (M_F) were measured. The plant dry matter (M_D) was obtained after drying at 60°C until constant weight.

Transcriptomic Analysis

Total RNA from whole plant tissue of pt^AtGSTT2a line and WT plants under control and high mannitol stress conditions in vitro was isolated using the Monarch Total RNA Miniprep kit (BioLabs Inc., UK) and their concentration was determined spectrophotometrically. The RNA sequencing was performed by the BGI (Denmark). The RNA results were compared as follows: pt^AtGSTT2a and WT in control conditions (groups 1 and 3) and pt^AtGSTT2a and WT in high mannitol (200 mM) stress (groups 2 and 4) (Table 1).

TABLE 1

Table 1. Plants used for RNA extraction and transcriptomics analysis.

Agilent 2100 Bioanalyzer (Agilent RNA 6000 Nano Kit) was used for the total RNA sample QC: RNA concentration, RIN value, 28S/18S and the fragment length distribution. We use NanoDrop^TM to identify the purity of the RNA samples. The first step in the workflow involves purifying the poly-A containing mRNA molecules using poly-T oligo-attached magnetic beads. Following purification, the mRNA is fragmented into small pieces using divalent cations under elevated temperature. The cleaved RNA fragments are copied into first strand cDNA using reverse transcriptase (Takara Bio Inc.) and random primers. This is followed by second strand cDNA synthesis using DNA Polymerase I and RNase H (Takara Bio Inc.). These cDNA fragments then have the addition of a single “A” base and subsequent ligation of the adapter. The products are then purified and enriched with PCR amplification. The PCR yield was quantified by Qubit and the samples were pooled together to make a single strand DNA circle (ssDNA circle), which gave the final library. DNA nanoballs (DNBs) were generated with the ssDNA circle by rolling circle replication (RCR) to enlarge the fluorescent signals at the sequencing process. The DNBs were loaded into the patterned nanoarrays and pair-end reads of 100 bp were read through on the BGISEQ-500 platform for the following data analysis study. For this step, the BGISEQ-500 platform combines the DNA nanoball-based nano arrays and stepwise sequencing using Combinational Probe-Anchor Synthesis Sequencing Method.

Bioinformatics Workflow

The reads were filtered for low-quality reads (>20% of the bases qualities are lower than 10), reads with adaptors and reads with unknown bases (N bases more than 5%) to get the clean reads using SOAPnuke software. Then we mapped the clean reads onto reference genome, followed by novel gene prediction, SNP & INDEL calling and gene splicing detection. Finally, we identified DEGs (differentially expressed genes) between samples and performed clustering analysis and functional annotations. After filtering, the remaining reads are called “Clean Reads” and stored in FASTQ format.

Regarding Genome Mapping, HISAT (Hierarchical Indexing for Spliced Alignment of Transcripts) was used to do the mapping step (Kim et al., 2015). The StringTie (Pertea et al., 2015) was used to reconstruct transcripts and Cuffcompare [Cufflinks tools- Trapnell et al. (2012)] to compare reconstructed transcripts to reference annotation. After that, the “u,” “I,” “o,” ”j” class code types were used as novel transcripts followed by a support vector machine-based classifier, named Coding Potential Calculator (CPC) (Kong et al., 2007) to predict coding potential of novel transcripts, then the coding novel transcripts were merged with reference transcripts to get a complete reference, and downstream analysis was based on this reference. The clean reads were mapped to reference using Bowtie2 (Langmead and Salzberg, 2012), and then gene expression level was calculated with RSEM (Li and Dewey, 2011), a software package for estimating gene and isoform expression levels from RNA-Seq data. After calculating Pearson's correlation between all samples using cor, hierarchical clustering was performed between all samples using hclust, and PCA analysis with all samples using princomp, and the diagrams were drawn with ggplot2 in R (R Core Team, 2016). The detection of DEGs was performed with NOIseq, which is based on noisy distribution model, as described in Tarazona et al. (2011). The Hierarchical Clustering Analysis of DEGs was performed using heat map function in R. With the GO annotation result, DEGs were classified according to official classification, and GO functional enrichment was also performed using p hyper in R. The p-value calculating formula in hypergeometric test is Equation 1:

$\begin{array}{l}P=1-{\displaystyle \sum _{i=0}^{m-1}}\frac{\left(\begin{array}{c}M\\ i\\ \end{array}\right)\left(\begin{array}{c}N-M\\ n-i\\ \end{array}\right)}{\left(\begin{array}{c}N\\ n\\ \end{array}\right)}\end{array}$

Then the false discovery rate (FDR) for each p-value was calculated and in general, the terms which FDR was not larger than 0.01 were defined as significantly enriched. With the KEGG annotation result, we classified DEGs according to official classification, and we also performed pathway functional enrichment using phyper in R with the same p-value calculating formula in Equation 1 and the FDR was calculated as described above.

To find the ORF of each DEG the getorf function was used. For plants, ORF were aligned to TF domains (from PlntfDB) using hmmsearch (Mistry et al., 2013). DIAMOND (Buchfink et al., 2014) was used to map the DEGs to the STRING database (von Mering et al., 2005) to obtain the interaction between DEG-encoded proteins using homology with known proteins. The top 100 interaction networks were selected to unfold the pathways involved and for the entire interaction result we provide an input file that can be imported directly into Cytoscape for complex network analysis and visualization.

Metabolite Extraction, Derivatization, and GC–MS Analysis

Determination of primary polar metabolites was performed as described by Lisec et al. (2006) and Michailidis et al. (2017) with slight modifications. Whole plant lyophilized material (~0.040 gr) from pt^AtGSTT2a and WT plants under in vitro high mannitol stress and control conditions (three biological replicates) were transferred in 2 mL screw cap tubes with 1400 μL of precooled (−20°C) pure methanol. Adonitol (100 μL of 0.2 mg mL⁻¹) was added as internal quantitative standard, and incubated for 10 min at 70°C. The supernatant was collected after centrifugation (11000 g, 4°C, 10 min) and 750 μL chloroform (−20°C) plus 1500 μL dH₂O (4°C) were added. Following centrifugation (2200 g, 4°C, 10 min), 150 μL of the upper polar phase were transferred into a 1.5 mL glass vial and placed under vacuum until drying. Dried residues were re-dissolved by gentle shaking in 40 μL of 20 mg mL⁻¹ methoxyamine hydrochloride for 120 min at 37°C, thereafter they were treated with 70 μL of N-methyl-N-(trimethylsilyl) trifluoroacetamide reagent (MSTFA), and incubated for 30 min at 37°C. GC-MS analysis was carried out in Thermo Trace Ultra GC equipped with ISQ MS and TriPlus RSH^TM auto-sampler (Thermo Fisher Scientific™, Switzerland). One μL was injected with a split ratio of 70:1. GC separation was held on a TR-5MS capillary column 30 m x 0.25 mm x 0.25 mm (Thermo Fisher Scientific™, Switzerland). Injector temperature was 220°C, ion source 230°C, and the interface 250°C. A constant flow of 1 mL min⁻¹ was used for carrier gas. The GC temperature program was held at 70°C for 2 min, then increased to 260°C (rate 8°C min⁻¹), where it remained for 18 min. Mass range of m/z 550 was recorded, after 5 min of solvent delay. The mass spectra were acquired in electron impact ionization mode. The peak area integration and chromatogram visualization was performed using the X-calibur processing program. Standards were used for peak identification or NIST11 database (Michailidis et al., 2017) in case of unknown peaks. The detected metabolites were assessed based on the relative response compared to adonitol and expressed as relative abundance.

Statistical Analysis

All the statistical analyses were performed using the computing environment R. The effects of stress treatments and the genotypes on the morpho-physiological parameters of the in vitro experiments and harvesting parameters of the in vivo experiments and treatments genotypes and time (days-where applicable) on the physiological parameters were assessed using two-way or three-way ANOVA, respectively, with the ez and afex packages (Lawrence, 2016; Singmann et al., 2018). All data were tested for normality (Shapiro test) and if normality failed and transformations were attempted. Data were also tested with Mauchly's test for sphericity, and if the assumption of sphericity was violated, the corresponding Greenhouse–Geisser corrections were performed. If significant differences were found among treatments, then the Tukey's HSD post hoc test was performed to determine specific treatment differences using the agricolae package (de Mendiburu, 2017). For metabolic data, two-way ANOVA was conducted using SPSS (SPSS v21.0., Chicago, USA) and statistically significant differences were based on Duncan's multiple range test (raw data) and Student's t-test for comparisons between genotypes or treatments at P < 0.05 (Table S1). The raw data are presented in Table S2 and the reported data are relative to the M_F of the pt^AtGSTT2a line and WT plants.

Results

Effect of GST Overexpression on Oxidative Stress Tolerance

Both low and high Diquat concentrations were severe enough to cause chlorotic lesions from day one (Figure S1) and senescence by day 2 on both WT plants and transplastomic lines (Figure S2). Transplastomic lines pt^AtGSTT2a and 6.1 showed chlorophyll content with increasing Diquat dose, compared to the WT control, except line pt^EFD6−115A, which showed reduced chlorophyll content following Diquat exposure (Table 2 and Table S3). Both Diquat doses negatively affected the maximum quantum efficiency of photosystem II (PSII) photochemistry in both WT and transplastomic lines, indicating that the oxidative stress was too severe possibly as a result of extensive free radical formation. The high diquat dose had a more severe effect on the M_F of transplastomic line pt^AtGSTT 6.1 (not statistically significant to the control) and WT plants (p < 0.05), and a less severe reduction was induced in pt^AtGSTT 2a (not statistically significant to the control) and pt^EFD6−115A (p < 0.05) (Table 2 and Table S3). M_D was not affected by any Diquat dose, indicating that any reduction in M_F was a result of turgor loss potentially, inhibiting the respiratory processes due to the function of Diquat as a rapid-acting translocated desiccant (Cronshey, 1961; McNaughton et al., 2015). An increase in M_D of pt^AtGSTT 2a was observed under both Diquat doses, however, the M_D of WT plants and transplastomic lines pt^AtGSTT 6.1 and pt^EFD6−115A was decreased with increasing Diquat concentration (Table 2).

TABLE 2

Table 2. Percent of change difference in growth (fresh-M_F and dry-M_D matter; g) and photophysiological parameters (relative chlorophyll content- Chl and Maximum quantum yield of PSII-Fv/Fm) of GST transplastomic lines and WT tobacco plants grown for 2 days in low (Diq_L) and high (Diq_H) Diquat dose compared to control conditions.

Tolerance of Transplastomic Lines Under in vitro NaCl Stress

Transplastomic pt^AtGSTT line 6-1 and 2a, when grown in 150 mM NaCl, exhibited increased tolerance compared to WT plants (Figure S3, Table 3 and Table S4), with shoot length and M_F not showing statistically significant differences with the stress-free plants. Transplastomic line pt^AtGSTT6-1 also showed a non-statistically significant decrease in the shoot length even under the double salt concentration, 300 mM NaCl. Transplastomic line pt^EFD6−115A showed reduced M_F in both salinity concentrations compared to the stress-free plants, however showed statistically significant increase in root length at low NaCl concentration (Table S3 and Table S4). Root length and maximum quantum efficiency of PSII were only reduced under the 300 mM NaCl concentration in all transplastomic lines and WT plants (Table 3). Wild-type plants showed the lowest chlorophyll content in both low and high NaCl concentrations compared to the transplastomic lines, although not significantly different. Relative chlorophyll content was maintained in 300 mM NaCl concentration in the transplastomic line pt^EFD6−115A, and it was reduced in all other genotypes including the WT compared to 150 mM NaCl concentration (Figure 1). Overall, the transplastomic line pt^AtGSTT2a demonstrated tolerance to both salt concentrations and especially at 150 mM NaCl as indicated by non-significant decrease in shoot length.

TABLE 3

Table 3. Morphological parameters and maximum quantum yield of PSII (Fv/Fm) of GST transplastomic lines and WT tobacco plants grown for 20 days in salinity stress (150 and 300 mM NaCl) in vitro.

FIGURE 1

Figure 1. Changes in the relative chlorophyll content in pt^AtGSTT (6-1 and 2a) and pt^EFD6−115A transplastomic lines, and WT plants growing under salinity (left panel) and drought (right panel) for 20 and 35 days, respectively. Different letters indicate significant differences between treatments for each genotype at P < 0.05 (n = 3 for control treatment and n = 6 for stress treatments).

Tolerance of Transplastomic Lines Under in vitro Mannitol Stress

Overexpression of the theta class AtGSTT in chloroplasts increased PS II functionality in both mannitol concentrations (100 and 200 mM) compared to the stress-free plants and relative to pt^EFD6−115A and WT plants, which only increased the Fv/Fm at low mannitol stress and reduced their quantum yield under high osmotic stress (Table 4 and Table S5). Additionally, the pt^AtGSTT lines showed increased chlorophyll content in 200 mM mannitol compared to pt^EFD6−115A and WT plants, which maintained their relative chlorophyll content in similar levels to the control plants (Figure 2). With respect to the effect of mannitol on growth, only pt^AtGSTT2a increased the shoot and root length in 200 mM mannitol, yet not significantly, in comparison to the stress-free plants, whilst the other transplastomic lines and the WT plants reduced the shoot length in high mannitol treatment. All transcriptomic lines reduced their M_F in the high mannitol stress, yet this reduction was less severe compared to the WT plants. Interestingly, in low mannitol stress only line pt^AtGSTT6.1 showed a non-significant reduction in the M_F and increased the shoot length and Fv/Fm, compared to the stress-free plants and the other transplastomic lines and WT plants (Table 4 and Table S5).

TABLE 4

Table 4. Morphological traits and maximum quantum yield of PSII (Fv/Fm) of transplastomic lines and WT tobacco plants grown for 35 days in osmotic stress (100 and 200 mM mannitol stress).

FIGURE 2

Figure 2. Effect of mannitol (100 and 200 mM) stress on growth of transplastomic lines and WT tobacco plants after 35 days in mannitol and stress-free (C) conditions.

Effect of GST Overexpression to the Transcriptome in Control and High Mannitol Stress

Based on the results the overexpression of the AtGSTT in the transplastomic 2a line resulted in enhanced tolerance to both salinity concentrations and osmotic stress (200 mM mannitol), along with tolerance to herbicide induced oxidative stress based on the increase in M_D and relative chlorophyll content. Taking into consideration the osmotic component of salinity, we have selected this line for further investigation of the changes occurred in transcriptome and metabolome level under osmotic stress as it looked to be the most promising one for acquired stress tolerance to investigate the whole transcriptome and metabolome response of this line in order to understand in a systemic way the response of the transplastomic line.

The transcriptome data were analyzed with RNA-Seq technology based on which, we performed the analysis of variance. The differential expression gene was selected according to the standard of P < 0.05 and the false discovery rate (FDR) was set to 0.001 to determine the threshold of the P-value for multiple tests. The absolute value of |log₂Ratio| ≥ 1 was used to determine the difference between the gene expression transcription group and the database. Gene function, annotation, and classification were researched by GO analysis (Figures 3A,B). The RNA analysis through next generation sequencing of the entire transcriptome of the pt^AtGSTT2a and WT plants under control conditions and under high mannitol (osmotic stress) was studied (each sample in duplicate). The comparison between the two samples in each group showed that the expression profile was similar, thus allowing their combination and their analysis. The analysis generated 47.105 million clean reads in total with a Q20 (%) 98.38%. The reads generated a total of 80.623 transcripts of which 51.879 are known genes and 28.744 are unknown genes. On average 93.55% reads are mapped, and the uniformity of the mapping result for each sample suggests that the samples are comparable. Analysis of differentially expressed genes (DEGs) (Table S6) between the transplastomic plants under control conditions (group 1) and WT under control conditions (group 3) showed that there are 80858 commonly expressed DEGs while there 4869 unique DEGs expressed in pt^AtGSTT2a and 3864 in WT (Figure 4A). Furthermore, when we applied the mannitol (osmotic stress simulating drought) we found 82765 commonly expressed DEGs, 4181 in pt^AtGSTT2a plants, and 3337 in WT plants, thus there is a difference of 1907 more common DEGs, 688 fewer DEGs in pt^AtGSTT2a plants, and 527 fewer DEGs in WT plants under stress showing a reduction in differentially expressed genes both in WT and pt^AtGSTT2a plants (Figure 4B).

FIGURE 3

Figure 3. Histogram representation of Gene Ontology classification. (A) control conditions (B) under high mannitol stress (200 mM).

FIGURE 4

Figure 4. Venn diagram of differentially expressed genes. Comparison among (A) pt^AtGSTT2a and WT in control conditions (groups 1 and 3) and (B) pt^AtGSTT2a and WT in high mannitol (200 mM) stress (groups 2 and 4).

Regarding the differentially expressed genes, between pt^AtGSTT2a and WT in control conditions (groups 1 and 3), we depicted 431 DEGs that were upregulated and 1500 downregulated (Figure 5A). Moreover, it is important to mention that between pt^AtGSTT2a and WT in high mannitol (200 mM) stress (groups 2 and 4), which are the samples under high mannitol stress only 264 were upregulated and 80 were downregulated (Figure 5B; Table S6).

FIGURE 5

Figure 5. Scatter plot of differentially expressed genes in (A) pt^AtGSTT2a and WT in control conditions (groups 1 and 3), and (B) pt^AtGSTT2a and WT in high mannitol (200 mM) stress (groups 2 and 4).

Analysis of pt^AtGSTT2a overexpressing line and WT plants before the application of the high mannitol (osmotic stress) showed that genes like alanine transaminase and glutamate decarboxylase both implicated in alanine metabolism and biosynthesis were upregulated in pt^AtGSTT2a under control conditions. Additionally, 2,3-bisphosphoglycerate-dependent, phosphoglycerate mutase, glycine hydroxymethyl transferase were also upregulated whilst phosphoserine phosphatase was downregulated in the pathway of Glycine, serine, and threonine metabolism. Thus, high glycine content should be expected in pt^AtGSTT2a plants under control conditions.

Histone-lysine N-methyltransferase ASH1L, was down regulated in pt^AtGSTT2a plants whereas in histidine metabolism genes responsible for phosphoribosyl-ATP pyrophosphohydrolase, phosphoribosyl-AMP cyclohydrolase, histidinol dehydrogenase, and histidine decarboxylase were upregulated, suggesting that histidine should be accumulated in the pt^AtGSTT2a plants. Glutamate decarboxylase implicated in the Alanine, aspartate and glutamate metabolism and taurine and hypotaurine metabolism was upregulated in the transplastomic lines. Furthermore, glutathione S-transferases were found to be overexpressed (BGI_novel_G036737 K00799: BGI_novel_G036737 (-2.6), BGI_novel_G008526 (-2.3), BGI_novel_G023189 (-2.3).

An important metabolic pathway related to stress tolerance is starch and sucrose metabolism where trehalose 6-phosphate synthase was downregulated in pt^AtGSTT2a line under control conditions suggesting that plants were in a state of stress-priming. Transplastomic line pt^AtGSTT2a upregulated inositol polyphosphate 5-phosphatase and inositol-pentakis phosphate 2-kinase in the Inositol phosphate metabolism and Phosphatidylinositol signaling system. Alanine transaminase involved in Alanine, aspartate, and glutamate metabolism was upregulated, whilst glycine hydroxyl methyltransferase and phosphoserine phosphatase were found to be downregulated in the biosynthesis of amino acids pathway.

In the comparison of pt^AtGSTT2a vs. WT plants under the high mannitol stress the number of transcripts for pt^AtGSTT2a and WT plants was similar to those in control conditions; however, there were only 264 DEGs upregulated and 80 downregulated compared with 431 and 1500, respectively, in control conditions (Table S6). Important genes found with altered expression are glycerate dehydrogenase and hydroxypyruvate reductase upregulated in Glycine, serine and threonine metabolism as well as in DNA repair pathway, which is expected as stress produces ROS, to also affect nucleic acids. Additionally, the gene responsible for spermidine synthase implicated in glutathione metabolism, cysteine, and methionine metabolism and in arginine and proline metabolism was downregulated along with pectinesterase, an important gene implicated in Pentose and glucuronate interconversions as well as in cell wall degradation, in pt^AtGSTT2a compared to WT plants under stress. In the Phenylalanine, tyrosine and tryptophan biosynthesis pathway, the genes encoding bifunctional anthranilate synthase/indole-3-glycerol-phosphate synthase (G005943) related to tryptophane biosynthesis were upregulated as was the 5-methyltetrahydrofolate-homocysteine methyltransferase (G005943), leading to methionine.

Effect of GST Overexpression in the Chloroplast to the Metabolome in High Mannitol Stress

The response of transplastomic line pt^AtGSTT2a and the WT plants was investigated further through the induced metabolic alterations. A total of 51 polar metabolites were identified (Figure 6; Table S1), of which 11 were soluble sugars, 5 soluble alcohols, 9 organic acids, 21 amino acids, and 5 other compounds (Figure 6; Tables S1, S2). The differences between the transplastomic line and WT plants under stress-free and high mannitol stress, revealed that ~78, 66.6, and 90% of the metabolic changes occurred due to treatment, genotypic, and treatment x genotype interaction effects, respectively.

FIGURE 6

Figure 6. Heat map of primary metabolites of pt^AtGSTT2a and WT plants under high mannitol (200 mM) stress compared to WT control plants. Increase is indicated as red and decrease as blue (see color scale). Mean values of 3 independent determinations for each treatment were expressed as relative abundance compared to internal standard adonitol and are reported relative to the respective M_F. Actual data are provided in Table S1.

In stress-free conditions, the overexpression of the AtGSTT had a significant effect to the metabolic profile of transplastomic plants as indicated by the 27 out of 51 metabolites being significantly altered from the WT plants, 26 of which were down-regulated and only the benzoic acid was up-regulated (27.6-fold) (Table 5). The downregulated metabolites were mostly amino acids (13), such as proline, oxoproline, and valine, organic acids (5), such as citric, quinic, and threonic acids, soluble sugars (5), such as sucrose, fructose, and glucose and the soluble alcohols, erythritol, myo-inositol and glycerol. In contrast, the effect of high mannitol stress on plants overexpressing the GST chimera was moderate as only 16 metabolites were significantly changed. More specifically, plants overexpressing the GST chimera up-regulated only the soluble sugars threose (9.5-fold) and arabinose (0.82-fold), whilst, 14 metabolites were downregulated compared to the WT plants (Table 5). These results indicate that pt^AtGSTT2a transplastomic line was osmotolerant and able to maintain cellular homeostasis in comparison to the WT plants that required more energy to tolerate high mannitol stress.

TABLE 5

Table 5. Metabolites that were significantly (p < 0.05) altered in pt^AtGSTT2a and WT plants under high mannitol (200 mM) and control conditions.

The pt^AtGSTT2a line under high mannitol stress significantly altered more metabolites (38) compared to the 34 metabolites of the WT plants (Figures 7, 8; Table S7). The pt^AtGSTT2a line upregulated six metabolites of which four were common. In the increased metabolites two were soluble sugars, such as arabinose, which was unique for the transplastomic line, two were soluble alcohols, such as mannitol, and quinic acid (Figure 7). Erythrose and sorbitol were accumulated in greater concentrations in the WT plants than in the transplastomic line under high mannitol compared to control conditions (Table S7). Interestingly, the compatible solute mannitol was accumulated in greater concentration in the pt^AtGSTT2a by 6.84-fold compared to the 4.73-fold increase in the WT plants.

FIGURE 7

Figure 7. Venn diagram representation of metabolites commonly or differentially increased in the leaves of WT and pt^AtGSTT2a tobacco plants under high mannitol (200 mM) compared to non-stressed plants.

FIGURE 8

Figure 8. Venn diagram representation of metabolites commonly or differentially decreased in the leaves of WT and pt^AtGSTT2a tobacco plants under high mannitol (200 mM) compared to non-stressed plants.

The WT plants down-regulated 31 metabolites compared to the 28 of the pt^AtGSTT2a plants. Among the metabolites that were decreased, 24 were common (6 soluble sugars, 12 amino acids, and 2 organic acids), such as the TCA cycle intermediate citric acid and the precursor of various amino acids aspartic acid (Figure 8; Table S7). Additionally, the WT plants had more differentially decreased metabolites compared to the pt^AtGSTT2a plants (Figure 8).

Discussion

The adaptation response mechanisms of plants to adverse abiotic stresses result in the up-regulation of the reactive oxygen species (ROS) detoxification network, to mitigate the negative effects of oxidative stress, commonly induced under such conditions (Gill and Tuteja, 2010; Nianiou-Obeidat et al., 2017). Enhancing the ROS scavenging capacity in plants by direct gene expression in the chloroplast, an active cell compartment could in theory increase the photosynthetic rate and thus increase in yield, yet this is only a speculation that needs thorough investigation, which is beyond the scope of this research. The functional role of tobacco lines overexpressing the AtGSTT in chloroplasts has been previously characterized (Dixon et al., 2008), yet, the roles of this AtGSTT in plant homeostasis and response mechanisms both under abiotic and herbicide-induced oxidative stresses, and non-stress conditions are still required to be unraveled. Targeting the chloroplasts, we have assessed the osmotic, ionic, and oxidative potential of the pt^AtGSTT lines and the ZmGSTU1-ZmGSTU2 chimera overexpressing line in comparison to WT plants. Our work shows that pt^AtGSTT lines were tolerant to herbicide-induced oxidative and salinity stresses and showed enhanced response tolerance to mannitol-induced osmotic stress compared to WT plants.

The mode of bipyridiniums action is within the chloroplast by diverting electrons from photosystem I (PSI) of photosynthesis to form the Diquat radical, which in turn generates a highly destructive superoxide radical (Devine et al., 1992; Hawkes, 2014). Despite the observed reduction in maximum quantum yield of PSII, the transplastomic lines pt^AtGSTT showed differential response mechanism in the accumulation of relative chlorophyll content in both half- and recommended field dose of Diquat, possibly to alleviate the negative effect of oxidative damage on PSII, which was apparent under all levels of Diquat in WT plants. It has been observed that transplastomic tobacco plants expressing simultaneously DHAR:GR and GST:GR showed enhanced tolerance to paraquat induced oxidative stresses while expression of either single transgene did not (Le Martret et al., 2011). Transplastomic overexpression of glutathione peroxidase (GP) in tobacco plants has shown to confer moderate tolerance to paraquat (Yoshimura et al., 2004), whereas, transplastomic tobacco lines overexpressing an Escherichia coli glutathione reductase (gor) gene have not enhanced protection from paraquat induced photooxidative stress (Poage et al., 2011). In the present study pt^AtGSTT lines also had enhanced turgor maintenance in contrast to WT plants, which showed extensive dehydration, since Diquat is a rapid desiccant (Cronshey, 1961; Hawkes, 2014).

At 150 mM of NaCl concentration, both pt^AtGSTT transplastomic lines showed enhanced growth (shoot length and M_F) and pt^AtGSTT2a moderate tolerance to high salinity stress (300 mM NaCl) by maintaining the shoot length compared to WT plants. Additionally, under both NaCl concentrations, all transplastomic lines demonstrated a higher relative chlorophyll content compared to the WT plants. Similar results were observed in transplastomic tobacco plants overexpressing a choline monooxygenase (BvCMO) from beetroot which increased photosynthetic rate and apparent quantum yield of photosynthesis in the presence of 150 mM NaCl when compared to WT, and the maximal efficiency of PSII photochemistry in both wild type and transplastomic plants was not affected (Zhang et al., 2008). This is also consistent with our results indicating that the overexpression of the AtGSTT and the ZmGSTU1-ZmGSTU2 chimera can protect PSII reaction centers from damage. Transplastomic carrot plants expressing the badh gene demonstrated enhanced tolerance up to up to 400 mM NaCl compared to untransformed plants exhibiting severe growth inhibition at 200 mM NaCl (Kumar et al., 2004). Herein, the pt^AtGSTT transplastomic lines pt^AtGSTT2a and especially pt^AtGSTT6-1 demonstrated enhanced photo-tolerance when exposed to 200 mM mannitol stress demonstrating increased relative chlorophyll content and maximum yield of PSII compared to WT plants. Increasing or maintaining the chlorophyll content in transgenic chloroplasts suggests the integrity of thylakoid membranes, even in the presence of high concentrations of NaCl and mannitol, demonstrating the advantage of overexpressing the AtGSTT in the chloroplasts. Similar results were observed by Lee et al. (2003).

The transcriptomics analysis of pt^AtGSTT2a line and WT plants under control and high mannitol stress suggests that plants of the pt^AtGSTT2a overexpressing line, before the application of the high mannitol (osmotic stress) upregulated genes related to stress tolerance such as genes encoding for alanine transaminase and glutamate decarboxylase both implicated in alanine metabolism and biosynthesis. Alanine was found to be the main amino acid accumulated in Medicago truncatula seedlings under hypoxic stress (Limami et al., 2008). Furthermore, genes encoding for 2,3-bisphosphoglycerate-dependent phosphoglycerate mutase and glycine hydroxymethyl transferase were also upregulated, while phosphoserine phosphatase was downregulated in glycine, serine and threonine metabolism, which may lead to increased glycine. Therefore, high glycine content should be expected in pt^AtGSTT2a plants under control conditions but considering that high glycine content is correlated with stress tolerance and especially drought resistance (Thankur and Rai, 1982), the pt^AtGSTT2a line is probably in a stress primed state before the application of the stress. . Also, glutamate decarboxylase is upregulated which is implicated in Alanine, aspartate and glutamate metabolism and in taurine and hypotaurine metabolism and was also confirmed in metabolomics analysis. Interestingly, in the resurrection plant Sporobolus stapfianus Martinelli et al. (2007) reported that the accumulation of asparagine and glutamate might have led to its conversion to arginine and asparagine, as all of the above are considered to play important role in plant protection against drought stress (Martinelli et al., 2007). Moreover, we found that the histone-lysine N-methyltransferase ASH1L, was down regulated in pt^AtGSTT2a plants in histidine metabolism whereas, genes responsible for histidine metabolism like phosphoribosyl-ATP pyrophosphohydrolase phosphoribosyl-AMP cyclohydrolase, histidinol dehydrogenase and histidine decarboxylase were upregulated, suggesting that histidine should be accumulated in the pt^AtGSTT2a plants prior to the application of stress, which reinforces the notion that the transplastomic plants are in a primed condition as before (Tran et al., 2007; Witt et al., 2012). Similarly, trehalose 6-phosphate synthase was found to be downregulated in pt^AtGSTT2a line under control conditions highlighting the stress primed condition of the transplastomic line (Lee et al., 2003; Ilhan et al., 2015).

When the transplastomic line pt^AtGSTT2a and WT plants were exposed to osmotic stress (high mannitol), the number of transcripts for pt^AtGSTT2a and WT plants did not change compared to those in control conditions. However, important genes found with an altered expression such as those encoding for glycerate dehydrogenase and hydroxypyruvate reductase which were found to be upregulated in glycine, serine and threonine metabolism as well as in DNA repair pathway, which, as stress produces ROS is expected to affect nucleic acids. In contrast, the gene responsible for spermidine synthase implicated in glutathione metabolism, cysteine and methionine metabolism and in arginine, and proline metabolism was down regulated. Proline is important in stress tolerance and has been found to increase during different environmental stresses like salinity, drought, UV, and extreme temperatures (Ashraf and Foolad, 2007). In addition, polyamines like spermidine have been reported to play a role in inducing stress response under various stresses that produce ROS as they might serve as ROS scavengers, and as positive regulators for expression of stress response genes. Thus, polyamines like spermidine could perform as primal stress molecules in plants (Rhee et al., 2007). Additionally, pectinesterase was found to be downregulated in pt^AtGSTT2a compared to WT plants under stress. This is an important gene implicated in pentose and glucuronate interconversions as well as in cell wall degradation as it was found to be upregulated in plants exposed to permissive high temperature conditions (37°C). This parallels to acclimation in order to acquire thermotolerance as a result of the cell wall modification (Yang et al., 2006). However, the downregulation of such enzymes in the pt^AtGSTT2a line under mannitol stress suggests that these plants might be in a state of acclimation prior to the application of the stress. In the Phenylalanine, tyrosine and tryptophan biosynthesis pathway the genes encoding bifunctional anthranilate synthase/indole-3-glycerol-phosphate synthase (G005943) related to tryptophane biosynthesis were upregulated as was the 5-methyltetrahydrofolate-homocysteine methyltransferase (G005943), leading to methionine. In a rat model, actin oxidative damage by ROS was found to occur through the oxidation of cysteine, tryptophan and methionine (Fedorova et al., 2010). If this is also the case in plants, then increased amounts of these amino acids might be needed and thus, leading to the upregulation of the genes responsible for their production, as it was found herein; however, this hypothesis needs further investigation.

The metabolomics analysis was performed on transplastomic and WT plants grown under high mannitol stress and controlled conditions in vitro for 35 days. The overexpression of the AtGSTT had a significant effect to the metabolic profile of transplastomic plants, since many metabolites were downregulated under both control and drought conditions indicating limited perturbation of metabolic homeostasis in the transplastomic lines. Especially under high mannitol stress the pt^AtGSTT2a line had higher concentrations of the soluble sugars, threose and arabinose, which demonstrates the protective role against osmotic stress (Keunen et al., 2013). The soluble alcohol mannitol was accumulated in greater concentration in pt^AtGSTT2a line despite the common increase in WT plants under high mannitol stress. In contrast to our results mannitol accumulation was decreased in transgenic tobacco plants overexpressing a Gmgstu4 gene under salinity stress (Kissoudis et al., 2015b). Mannitol accumulation plays an important role in osmotic adjustment and signaling molecule enhance tolerance to water stress in various plant species (Slama et al., 2015). Additionally, the greater M_F and shoot length of pt^AtGSTT2a compared to the WT plants under mannitol stress indicates a possible relation between increase in mannitol and improved growth. Similar results were observed in peanut (Bhauso et al., 2014) and Zea mays (Nguyen et al., 2013) plants overexpressing mtlD genes, which conferred water-deficit stress tolerance by inducing the accumulation of mannitol and increase in biomass and relative water content under drought conditions.

The results above suggest that overexpression of the AtGSTT in the chloroplasts resulted in enhanced photo-tolerance and turgor maintenance under herbicide-induced oxidative (increased M_D and Relative chlorophyll content) and salinity stresses (higher chlorophyll, non-significant decrease in shoot length and M_F compared to the control plants) and enhanced response tolerance to high mannitol-induced osmotic stress (increased shoot and root length). Whole-genome transcriptome analysis revealed that genes related to stress tolerance, such as GSTs, were upregulated in pt^AtGSTT2a line under both control and high mannitol stress conditions indicating an acclimation state to stress. In parallel, the metabolic profile indicated limited perturbations of the metabolic homeostasis in the transplastomic lines and greater accumulation of mannitol and soluble sugars under high mannitol stress. We have therefore established that the transplastomic plants overexpressing the pt^AtGSTT2a in the chloroplast are probably in a state of acclimation to stress, thus, when the actual stress is applied there is limited need for overexpression of the whole array of stress tolerance mechanisms, which is imprinted in the levels of relative gene expression. As mentioned before, we found only limited genes to be upregulated in the pt^AtGSTT2a transplastomic line compared to WT under stress conditions while at the same time we have found genes related to stress tolerance upregulated in pt^AtGSTT2a plants compared to WT in stress-free conditions, strengthening the hypothesis that the AtGSTT overexpressed in plastids might have conferred plant stress tolerance.

Challenges caused by climate change will demand for quick action of the scientific community in order to develop stress tolerant varieties to secure enough food for the increasing world population. GSTs, for have proven to be enzymes involved in stress tolerance (Dixon et al., 1998, 2008, 2011; Axarli et al., 2009, 2017; Chronopoulou and Labrou, 2009; Benekos et al., 2010; Chronopoulou et al., 2011, 2012, 2014; Madesis et al., 2013; Kissoudis et al., 2015a,b; Labrou et al., 2015; Lo Cicero et al., 2015, 2017; Nianiou-Obeidat et al., 2017) might help toward the development of plant acclimation to environmental stresses. In some cases, the overexpression of a single antioxidant enzyme might not provide protection against oxidative stress whilst, simultaneous expression of multiple antioxidant enzymes is more effective than a single expression for enhancing tolerance to environmental stresses (Le Martret et al., 2011). Herein, the ZmGSTU1-ZmGSTU2 chimera was able to induce photoprotection of the photosystem II under sever salinity stress, yet it was not as tolerant as the single AtGSTTs overexpressed in the chloroplasts. Potentially, the expression of multiple defense genes encoding enzymes belonging to different classes could generate plants with enhanced stress tolerance (Zhao and Zhang, 2006) able to withstand multiple stresses, which needs to be further investigated. This study provides evidence that overexpression of both the theta class AtGSTT and the unique chimera GSTU1-GSTU2 from Zea mays in the chloroplast resulted in enhanced tolerance of the transplastomic plants to abiotic stresses. Furthermore, transcriptomics and metabolomics analysis showed that the GST overexpressing plants were in a stress tolerance priming state even before the application of the severe osmotic stress (high mannitol concentration) thus, enhancing the plant's ability to tolerate abiotic stresses.

Author Contributions

ES designed and performed part of the experiments, and the statistical analysis, and wrote the manuscript. SG and AI performed part of the experiments. SK performed part of the experiments. MM performed part of the metabolomics and wrote part of the metabolomics section. EC performed part of the experiments and read the manuscript. RE wrote and edited part of the manuscript. AD wrote and edited the manuscript. NL, IN-O and PM designed the research, wrote and edited the manuscript.

Conflict of Interest Statement

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.

Acknowledgments

We acknowledge the General Secretariat for Research and Development that financed the work described here under the project Intelligent και Automated systems for enabling the design, simulation and development of integrated processes and products –ODYSSEAS, (MIS) 5002462.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.01861/full#supplementary-material

References

Ahmad, N., Michoux, F., McCarthy, J., and Nixon, P. J. (2012). Expression of the affinity tags, glutathione-S-transferase and maltose-binding protein, in tobacco chloroplasts. Planta 235, 863–871. doi: 10.1007/s00425-011-1584-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Ashraf, M., and Foolad, M. R. (2007). Roles of glycine betaine and proline in improving plant abiotic stress resistance. Environ. Exp. Bot. 59, 206–216. doi: 10.1016/j.envexpbot.2005.12.006

CrossRef Full Text | Google Scholar

Axarli, I., Dhavala, P., Papageorgiou, A. C., and Labrou, N. E. (2009). Crystallographic and functional characterization of the fluorodifen-inducible glutathione transferase from Glycine max reveals an active site topography suited for diphenylether herbicides and a novel L-site. J. Mol. Biol. 385, 984–1002. doi: 10.1016/j.jmb.2008.10.084

PubMed Abstract | CrossRef Full Text | Google Scholar

Axarli, I., Muleta, A. W., Chronopoulou, E. G., Papageorgiou, A. C., and Labrou, N. E. (2017). Directed evolution of glutathione transferases towards a selective glutathione-binding site and improved oxidative stability. Biochim. Biophys. Acta Gen. Subj. 1861, 3416–3428. doi: 10.1016/j.bbagen.2016.09.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Benekos, K., Kissoudis, C., Nianiou-Obeidat, I., Labrou, N., Madesis, P., Kalamaki, M., et al. (2010). Overexpression of a specific soybean GmGSTU4 isoenzyme improves diphenyl ether and chloroacetanilide herbicide tolerance of transgenic tobacco plants. J. Biotechnol. 150, 195–201. doi: 10.1016/j.jbiotec.2010.07.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Bhauso, T. D., Radhakrishnan, T., Kumar, A., Mishra, G. P., Dobaria, J. R., Patel, K., et al. (2014). Overexpression of bacterial mtlD gene in peanut improves drought tolerance through accumulation of mannitol. Sci. World J. 2014:125967. doi: 10.1155/2014/125967

PubMed Abstract | CrossRef Full Text | Google Scholar

Bock, R. (2014). Genetic engineering of the chloroplast: novel tools and new applications. Curr. Opin. Biotechnol. 26, 7–13. doi: 10.1016/j.copbio.2013.06.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Buchfink, B., Xie, C., and Huson, D. H. (2014). Fast and sensitive protein alignment using diamond. Nat. Methods 12, 59–60. doi: 10.1038/nmeth.3176

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, P. J., Senthilkumar, R., Jane, W. N., He, Y., Tian, Z., and Yeh, K. W. (2014). Transplastomic Nicotiana benthamiana plants expressing multiple defence genes encoding protease inhibitors and chitinase display broad-spectrum resistance against insects, pathogens and abiotic stresses. Plant Biotechnol. J. 12, 503–515. doi: 10.1111/pbi.12157

PubMed Abstract | CrossRef Full Text | Google Scholar

Chin, H. G., Kim, G.-D., Marin, I., Mersha, F., Evans, T. C., Chen, L., et al. (2003). Protein trans-splicing in transgenic plant chloroplast: reconstruction of herbicide resistance from split genes. Proc. Natl. Acad. Sci. U.S.A. 100, 4510–4515. doi: 10.1073/pnas.0736538100

PubMed Abstract | CrossRef Full Text | Google Scholar

Chronopoulou, E., Axarli, I., Nianiou-obeidat, I., Madesis, P., and Labrou, N. E. (2011). Structure and antioxidant catalytic function of plant glutathione transferases. Curr. Chem. Biol. 5, 64–74. doi: 10.2174/187231311793564306

CrossRef Full Text | Google Scholar

Chronopoulou, E., Madesis, P., Asimakopoulou, B., Platis, D., Tsaftaris, A., and Labrou, N. E. (2012). Catalytic and structural diversity of the fluazifop-inducible glutathione transferases from Phaseolus vulgaris. Planta 235, 1253–1269. doi: 10.1007/s00425-011-1572-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Chronopoulou, E., Madesis, P., Tsaftaris, A., and Labrou, N. E. (2014). Cloning and characterization of a biotic-stress-inducible glutathione transferase from Phaseolus vulgaris. Appl. Biochem. Biotechnol. 172, 595–609. doi: 10.1007/s12010-013-0509-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Chronopoulou, E. G., and Labrou, N. E. (2009). Glutathione transferases : emerging multidisciplinary tools in red and green biotechnology. Recent Pat. Biotechnol. 3, 211–223. doi: 10.2174/187220809789389135

PubMed Abstract | CrossRef Full Text | Google Scholar

Clarke, J. L., and Daniell, H. (2011). Plastid biotechnology for crop production: present status and future perspectives. Plant Mol. Biol. 76, 211–220. doi: 10.1007/s11103-011-9767-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Cronshey, J. F. H. (1961). A review of experimental work with diquat and related compounds. Weed Res. 1, 68–77. doi: 10.1111/j.1365-3180.1961.tb00007.x

CrossRef Full Text | Google Scholar

Csiszár, J., Horváth, E., Váry, Z., Gallé, Á., Bela, K., Brunner, S., et al. (2014). Glutathione transferase supergene family in tomato: salt stress-regulated expression of representative genes from distinct GST classes in plants primed with salicylic acid. Plant Physiol. Biochem. 78, 15–26. doi: 10.1016/j.plaphy.2014.02.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Czégény, G., Le Martret, B., Pávkovics, D., Dix, P. J., and Hideg, É. (2016). Elevated ROS-scavenging enzymes contribute to acclimation to UV-B exposure in transplastomic tobacco plants, reducing the role of plastid peroxidases. J. Plant Physiol. 201, 95–100. doi: 10.1016/j.jplph.2016.05.026

PubMed Abstract | CrossRef Full Text | Google Scholar

de Mendiburu, F. (2017). agricolae: Statistical Procedures for Agricultural Research. R package version 1.2-8. Available online at: https://CRAN.R-project.org/package=agricola

Devine, M., Duke, S. O., and Fedtke, C. (1992). Physiology of herbicide action. Exp. Agric. 29:441. doi: 10.1017/S0014479700021281

CrossRef Full Text | Google Scholar

Dixon, D. P., Cummins, L., Cole, D. J., and Edwards, R. (1998). Glutathione-mediated detoxification systems in plants. Curr. Opin. Cell Biol. 1, 258–266.

PubMed Abstract | Google Scholar

Dixon, D. P., and Edwards, R. (2010). “Gutathione transferases,” in The Arabidopsis Book American Society of Plant Physiologist. doi: 10.1199/tab.0131

CrossRef Full Text

Dixon, D. P., Hawkins, T., Hussey, P. J., and Edwards, R. (2009). Enzyme activities and subcellular localization of members of the arabidopsis glutathione transferase superfamily. J. Exp. Bot. 60, 1207–1218. doi: 10.1093/jxb/ern365

PubMed Abstract | CrossRef Full Text | Google Scholar

Dixon, D. P., Lapthorn, A., Madesis, P., Mudd, E. A., Day, A., and Edwards, R. (2008). Binding and glutathione conjugation of porphyrinogens by plant glutathione transferases. J. Biol. Chem. 283, 20268–20276. doi: 10.1074/jbc.M802026200

PubMed Abstract | CrossRef Full Text | Google Scholar

Dixon, D. P., McEwen, A. G., Lapthorn, A. J., and Edwards, R. (2003). Forced evolution of a herbicide detoxifying glutathione transferase. J. Biol. Chem. 278, 23930–23935. doi: 10.1074/jbc.M303620200

PubMed Abstract | CrossRef Full Text | Google Scholar

Dixon, D. P., Sellars, J. D., and Edwards, R. (2011). The Arabidopsis phi class glutathione transferase AtGSTF2: binding and regulation by biologically active heterocyclic ligands. Biochem. J. 438, 63–70. doi: 10.1042/BJ20101884

PubMed Abstract | CrossRef Full Text | Google Scholar

Dufourmantel, N., Dubald, M., Matringe, M., Canard, H., Garcon, F., Job, C., et al. (2007). Generation and characterization of soybean and marker-free tobacco plastid transformants over-expressing a bacterial 4-hydroxyphenylpyruvate dioxygenase which provides strong herbicide tolerance. Plant Biotechnol. J. 5, 118–133. doi: 10.1111/j.1467-7652.2006.00226.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Falk, J., Brosch, M., Schäfer, A., Braun, S., and Krupinska, K. (2005). Characterization of transplastomic tobacco plants with a plastid localized barley 4-hydroxyphenyl-pyruvate dioxygenase. J. Plant Physiol. 162, 738–742. doi: 10.1016/j.jplph.2005.04.005

CrossRef Full Text | Google Scholar

Fedorova, M., Kuleva, N., and Hoffmann, R. (2010). Identification of cysteine, methionine and tryptophan residues of actin oxidized in vivo during oxidative stress. J. Proteome Res. 9, 1598–1609. doi: 10.1021/pr901099e

PubMed Abstract | CrossRef Full Text | Google Scholar

Foyer, C. H., and Shigeoka, S. (2011). Understanding oxidative stress and antioxidant functions to enhance photosynthesis. Plant Physiol. 155, 93–100. doi: 10.1104/pp.110.166181

PubMed Abstract | CrossRef Full Text | Google Scholar

Gill, S. S., and Tuteja, N. (2010). Reactive oxygen species and antioxidant machinery in abiotic stress tolerance in crop plants. Plant Physiol. Biochem. 48, 909–930. doi: 10.1016/j.plaphy.2010.08.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Grant, O. M., Brennan, D. P., Mellisho Salas, C. D., and Dix, P. J. (2014). Impact of enhanced capacity to scavenge reactive oxygen species on cold tolerance of tobacco. Int. J. Plant Sci. 175, 544–554. doi: 10.1086/675976

CrossRef Full Text | Google Scholar

Hawkes, T. R. (2014). Mechanisms of resistance to paraquat in plants. Pest Manag. Sci. 70, 1316–1323. doi: 10.1002/ps.3699

PubMed Abstract | CrossRef Full Text | Google Scholar

Iamtham, S., and Day, A. (2000). Removal of antibiotic resistance genes from transgenic tobacco plastids. Nat. Biotechnol. 18, 1172–1176. doi: 10.1038/81161

PubMed Abstract | CrossRef Full Text | Google Scholar

Ilhan, S., Ozdemir, F., and Bor, M. (2015). Contribution of trehalose biosynthetic pathway to drought stress tolerance of Capparis ovata Desf. Plant Biol. 17, 402–407. doi: 10.1111/plb.12227

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, S., and Daniell, H. (2014). Expression of γ-tocopherol methyltransferase in chloroplasts results in massive proliferation of the inner envelope membrane and decreases susceptibility to salt and metal-induced oxidative stresses by reducing reactive oxygen species. Plant Biotechnol. J. 12, 1274–1285. doi: 10.1111/pbi.12224

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanai, T., Takahashi, K., and Inoue, H. (2006). Three distinct-type glutathione S-transferases from Escherichia coli important for defense against oxidative stress. J. Biochem. 140, 703–711. doi: 10.1093/jb/mvj199

PubMed Abstract | CrossRef Full Text | Google Scholar

Kayum, M., Nath, U. K., Park, J. I., Biswas, M. K., Choi, E. K., Song, J. Y., et al. (2018). Genome-wide identification, characterization, and expression profiling of glutathione S-transferase (GST) family in pumpkin reveals likely role in cold-stress tolerance. Genes 9, 1–21. doi: 10.3390/genes9020084

PubMed Abstract | CrossRef Full Text | Google Scholar

Keunen, E., Peshev, D., Vangronsveld, J., Van Den Ende, W., and Cuypers, A. (2013). Plant sugars are crucial players in the oxidative challenge during abiotic stress: extending the traditional concept. Plant Cell Environ. 36, 1242–1255. doi: 10.1111/pce.12061

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, D., Langmead, B., and Salzberg, S. L. (2015). HISAT: a fast spliced aligner with low memory requirements. Nat. Methods 12, 357–360. doi: 10.1038/nmeth.3317

PubMed Abstract | CrossRef Full Text | Google Scholar

Kissoudis, C., Kalloniati, C., Flemetakis, E., Madesis, P., Labrou, N. E., Tsaftaris, A., et al. (2015a). Maintenance of metabolic homeostasis and induction of cytoprotectants and secondary metabolites in alachlor-treated GmGSTU4-overexpressing tobacco plants, as resolved by metabolomics. Plant Biotechnol. Rep. 9, 287–296. doi: 10.1007/s11816-015-0364-5

CrossRef Full Text | Google Scholar

Kissoudis, C., Kalloniati, C., Flemetakis, E., Madesis, P., Labrou, N. E., Tsaftaris, A., et al. (2015b). Stress-inducible GmGSTU4 shapes transgenic tobacco plants metabolome towards increased salinity tolerance. Acta Physiol. Plant. 37, 1–11. doi: 10.1007/s11738-015-1852-5

CrossRef Full Text | Google Scholar

Kong, L., Zhang, Y., Ye, Z. Q., Liu, X. Q., Zhao, S. Q., Wei, L., et al. (2007). CPC: assess the protein-coding potential of transcripts using sequence features and support vector machine. Nucleic Acids Res. 35, 345–349. doi: 10.1093/nar/gkm391

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, S., Asif, M. H., Chakrabarty, D., Tripathi, R. D., Dubey, R. S., and Trivedi, P. K. (2013). Expression of a rice Lambda class of glutathione S-transferase, OsGSTL2, in Arabidopsis provides tolerance to heavy metal and other abiotic stresses. J. Hazard. Mater. 248–249, 228–237. doi: 10.1016/j.jhazmat.2013.01.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, S., Dhingra, A., and Daniell, H. (2004). Plastid-expressed betaine aldehyde dehydrogenase gene in carrot cultured cells, roots, and leaves confers enhanced salt tolerance. Plant Physiol. 136, 2843–2854. doi: 10.1104/pp.104.045187.and

PubMed Abstract | CrossRef Full Text | Google Scholar

Labrou, N. E., Papageorgiou, A. C., Pavli, O., and Flemetakis, E. (2015). Plant GSTome: structure and functional role in xenome network and plant stress response. Curr. Opin. Biotechnol. 32, 186–194. doi: 10.1016/j.copbio.2014.12.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Langmead, B., and Salzberg, S. L. (2012). Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359. doi: 10.1038/nmeth.1923

PubMed Abstract | CrossRef Full Text | Google Scholar

Lawrence, M. A. (2016). ez: Easy Analysis and Visualization of Factorial Experiments. R package version 4.4-0. Available online at: https://CRAN.R-project.org/package=ez

Le Martret, B., Poage, M., Shiel, K., Nugent, G. D., and Dix, P. J. (2011). Tobacco chloroplast transformants expressing genes encoding dehydroascorbate reductase, glutathione reductase, and glutathione-S-transferase exhibit altered anti-oxidant metabolism and improved abiotic stress tolerance. Plant Biotechnol. J. 9, 661–673. doi: 10.1111/j.1467-7652.2011.00611.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, S. B., Kwon, H.-B., Kwon, S.-J., Park, S.-C., Jeong, M.-J., Han, S.-E., et al. (2003). Accumulation of trehalose within transgenic chloroplasts confers drought tolerance. Mol. Breed. 11, 1–13. doi: 10.1023/A:1022100404542

CrossRef Full Text | Google Scholar

Li, B., and Dewey, C. N. (2011). RSEM : accurate transcript quantification from RNA - Seq data with or without a reference genome. BMC Bioinformatics 12, 1–16. doi: 10.1186/1471-2105-12-323

CrossRef Full Text | Google Scholar

Limami, A. M., Glévarec, G., Ricoult, C., Cliquet, J. B., and Planchet, E. (2008). Concerted modulation of alanine and glutamate metabolism in young Medicago truncatula seedlings under hypoxic stress. J. Exp. Bot. 59, 2325–2335. doi: 10.1093/jxb/ern102

PubMed Abstract | CrossRef Full Text | Google Scholar

Lisec, J., Schauer, N., Kopka, J., Willmitzer, L., and Fernie, A. R. (2006). Gas chromatography mass spectrometry-based metabolite profiling in plants. Nat. Protoc. 1, 387–396. doi: 10.1038/nprot.2006.59

PubMed Abstract | CrossRef Full Text | Google Scholar

Lo Cicero, L., Catara, V., Strano, C. P., Bella, P., Madesis, P., and Lo Piero, A. R. (2017). Over-expression of CsGSTU promotes tolerance to the herbicide alachlor and resistance to Pseudomonas syringae pv. tabaci in transgenic tobacco. Biol. Plant. 61, 169–177. doi: 10.1007/s10535-016-0659-6

CrossRef Full Text | Google Scholar

Lo Cicero, L., Madesis, P., Tsaftaris, A., and Lo Piero, A. R. (2015). Tobacco plants over-expressing the sweet orange tau glutathione transferases (CsGSTUs) acquire tolerance to the diphenyl ether herbicide fluorodifen and to salt and drought stresses. Phytochemistry 116, 69–77. doi: 10.1016/j.phytochem.2015.03.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Lutz, K. A., Knapp, J. E., and Maliga, P. (2001). Expression of bar in the plastid genome confers herbicide resistance. Plant Physiol. 125, 1585–1590. doi: 10.1104/pp.125.4.1585

PubMed Abstract | CrossRef Full Text | Google Scholar

Madesis, P., Nianiou-Obeidat, I., Labrou, N. E., and Tsaftaris, A. (2013). “The Use of highly specific GSTs towards the development of stress tolerant transgenic plants,” in Glutathione: Biochemistry, Mechanisms of Action and Biotechnological Implications, eds N. E. Labrou, and E. Flemetakis (New York, NY: Nova Science Pub Inc.) 263–274. doi: 10.1146/annurev.soc.30.012703.110603

CrossRef Full Text | Google Scholar

Martinelli, T., Whittaker, A., Bochicchio, A., Vazzana, C., Suzuki, A., and Masclaux-Daubresse, C. (2007). Amino acid pattern and glutamate metabolism during dehydration stress in the “resurrection” plant Sporobolus stapfianus: a comparison between desiccation-sensitive and desiccation-tolerant leaves. J. Exp. Bot. 58, 3037–3046. doi: 10.1093/jxb/erm161

PubMed Abstract | CrossRef Full Text | Google Scholar

McNaughton, K. E., Blackshaw, R. E., Waddell, K. A., Gulden, R. H., Sikkema, P. H., and Gillard, C. L. (2015). Effect of five desiccants applied alone and in combination with glyphosate in dry edible bean (Phaseolus vulgaris L.). Can. J. Plant Sci. 95, 1235–1242. doi: 10.4141/cjps-2015-098

CrossRef Full Text | Google Scholar

Michailidis, M., Karagiannis, E., Tanou, G., Karamanoli, K., Lazaridou, A., Matsi, T., et al. (2017). Metabolomic and physico-chemical approach unravel dynamic regulation of calcium in sweet cherry fruit physiology. Plant Physiol. Biochem. 116, 68–79. doi: 10.1016/j.plaphy.2017.05.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Mistry, J., Finn, R. D., Eddy, S. R., Bateman, A., and Punta, M. (2013). Challenges in homology search: HMMER3 and convergent evolution of coiled-coil regions. Nucleic Acids Res. 41:e121. doi: 10.1093/nar/gkt263

PubMed Abstract | CrossRef Full Text | Google Scholar

Mullineaux, P., and Karpinski, S. (2002). Signal transduction in response to excess light: getting out of the chloroplast. Curr. Opin. Plant Biol. 5, 43–48. doi: 10.1016/S1369-5266(01)00226-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Nguyen, T. X., Nguyen, T., Alameldin, H., Goheen, B., Loescher, W., and Sticklen, M. (2013). Transgene pyramiding of the HVA1 and mtlD in T3 Maize (Zea mays L.) plants confers drought and salt tolerance, along with an increase in crop biomass. Int. J. Agron. 2013, 1–10. doi: 10.1155/2013/598163

CrossRef Full Text | Google Scholar

Nianiou-Obeidat, I., Madesis, P., Kissoudis, C., Voulgari, G., Chronopoulou, E., Tsaftaris, A., et al. (2017). Plant glutathione transferase-mediated stress tolerance: functions and biotechnological applications. Plant Cell Rep. 36, 791–805. doi: 10.1007/s00299-017-2139-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Nishida, M., Kong, K.-H., Inoue, H., and Takahashi, K. (1994). Molecular cloning and site-directed mutagenesis of glutathione s-transferase from Escherichia coli. J. Biol. Chem. 269, 32536–32541.

PubMed Abstract | Google Scholar

Pertea, M., Pertea, G. M., Antonescu, C. M., Chang, T. C., Mendell, J. T., and Salzberg, S. L. (2015). Stringtie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 33, 290–295. doi: 10.1038/nbt.3122

PubMed Abstract | CrossRef Full Text | Google Scholar

Poage, M., Le Martret, B., Jansen, M. A., Nugent, G. D., and Dix, P. J. (2011). Modification of reactive oxygen species scavenging capacity of chloroplasts through plastid transformation. Plant Mol. Biol. 76, 371–384. doi: 10.1007/s11103-011-9784-y

PubMed Abstract | CrossRef Full Text | Google Scholar

R Core Team (2016). R: A Language and Environment for Statistical Computing. Available online at: https://cran.r-project.org/doc/manuals/r-release/fullrefman.pdf

Rhee, H. J., Kim, E. J., and Lee, J. K. (2007). Physiological polyamines: simple primordial stress molecules. J. Cell. Mol. Med. 11, 685–703. doi: 10.1111/j.1582-4934.2007.00077.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Sappl, P. G., Carroll, A. J., Clifton, R., Lister, R., Whelan, J., Harvey Millar, A., et al. (2009). The Arabidopsis glutathione transferase gene family displays complex stress regulation and co-silencing multiple genes results in altered metabolic sensitivity to oxidative stress. Plant J. 58, 53–68. doi: 10.1111/j.1365-313X.2008.03761.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Singmann, H., Bolker, B., Westfall, J., and Aust, F. (2018). afex: Analysis of Factorial Experiments. R package version 0.19-1. Available online at: https://CRAN.R-project.org/package=afex

Slama, I., Abdelly, C., Bouchereau, A., Flowers, T., and Savouré, A. (2015). Diversity, distribution and roles of osmoprotective compounds accumulated in halophytes under abiotic stress. Ann. Bot. 115, 433–447. doi: 10.1093/aob/mcu239

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, D. B., and Johnson, K. S. (1988). Single-step purification of polypeptides expressed in Escherichia coli as fusions with glutathione S-transferase. Gene 67, 31–40. doi: 10.1016/0378-1119(88)90005-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Stavridou, E., Hastings, A., Webster, R. J., and Robson, P. R. H. (2016). The impact of soil salinity on the yield, composition and physiology of the bioenergy grass Miscanthus x giganteus. GCB Bioenergy 9, 92–104. doi: 10.1111/gcbb.12351

CrossRef Full Text | Google Scholar

Tarazona, S., García-Alcalde, F., Dopazo, J., Ferrer, A., and Conesa, A. (2011). Differential expression in RNA-seq: a matter of depth. Genome Res. 21, 2213–2223. doi: 10.1101/gr.124321.111

PubMed Abstract | CrossRef Full Text | Google Scholar

Thankur, P. S., and Rai, V. K. (1982). Dynamics of amino acid accumulation of two differentially drought resistant Zea Mays cultivars in response to osmotic stress. Environ. Exp. Bot. 22, 221–225. doi: 10.1016/0098-8472(82)90042-9

CrossRef Full Text | Google Scholar

Tran, L. P., Urao, T., Qin, F., Maruyama, K., Kakimoto, T., Shinozaki, K., et al. (2007). Functional analysis of AHK1/ATHK1 and cytokinin receptor histidine kinases in response to abscisic acid, drought, and salt stress in Arabidopsis. PNAS 104, 20623–20628. doi: 10.1073/pnas.0706547105

CrossRef Full Text | Google Scholar

Trapnell, C., Roberts, A., Goff, L., Pertea, G., Kim, D., Kelley, D. R., et al. (2012). Differential gene and transcript expression analysis of RNA-seq experiments with tophat and cufflinks. Nat. Protoc. 7, 562–578. doi: 10.1038/nprot.2012.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Verma, D., and Daniell, H. (2007). Chloroplast vector systems for biotechnology applications. Plant Physiol. 145, 1129–1143. doi: 10.1104/pp.107.106690

PubMed Abstract | CrossRef Full Text | Google Scholar

von Mering, C., Jensen, L. J., Snel, B., Hooper, S. D., Krupp, M., Foglierini, M., et al. (2005). STRING: known and predicted protein-protein associations, integrated and transferred across organisms. Nucleic Acids Res. 33, 433–437. doi: 10.1093/nar/gki005

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, H. H., Yin, W. B., and Hu, Z. M. (2009). Advances in chloroplast engineering. J. Genet. Genomics 36, 387–398. doi: 10.1016/S1673-8527(08)60128-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Witt, S., Galicia, L., Lisec, J., Cairns, J., Tiessen, A., Araus, J. L., et al. (2012). Metabolic and phenotypic responses of greenhouse-grown maize hybrids to experimentally controlled drought stress. Mol. Plant 5, 401–417. doi: 10.1093/mp/ssr102

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, K. A., Lim, C. J., Hong, J. K., Park, C. Y., Cheong, Y. H., Chung, W. S., et al. (2006). Identification of cell wall genes modified by a permissive high temperature in Chinese cabbage. Plant Sci. 171, 175–182. doi: 10.1016/j.plantsci.2006.03.013

CrossRef Full Text | Google Scholar

Ye, G. N., Hajdukiewicz, P. T., Broyles, D., Rodriguez, D., Xu, C. W., Nehra, N., et al. (2001). Plastid-expressed 5-enolpyruvylshikimate-3-phosphate synthase genes provide high level glyphosate tolerance in tobacco. Plant J. 25, 261–270. doi: 10.1046/j.1365-313X.2001.00958.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoshimura, K., Miyao, K., Gaber, A., Takeda, T., Kanaboshi, H., Miyasaka, H., et al. (2004). Enhancement of stress tolerance in transgenic tobacco plants overexpressing Chlamydomonas glutathione peroxidase in chloroplasts or cytosol. Plant J. 37, 21–33. doi: 10.1046/j.1365-313X.2003.01930.x

CrossRef Full Text | Google Scholar

Zhang, J., Tan, W., Yang, X. H., and Zhang, H. X. (2008). Plastid-expressed choline monooxygenase gene improves salt and drought tolerance through accumulation of glycine betaine in tobacco. Plant Cell Rep. 27, 1113–1124. doi: 10.1007/s00299-008-0549-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, F., and Zhang, H. (2006). Salt and paraquat stress tolerance results from co-expression of the Suaeda salsa glutathione S-transferase and catalase in transgenic rice. Plant Cell. Tissue Organ Cult. 86, 349–358. doi: 10.1007/s11240-006-9133-z

CrossRef Full Text | Google Scholar