Protein kinase GCN2 mediates responses to glyphosate in Arabidopsis
© Faus et al.; licensee BioMed Central. 2015
Received: 20 June 2014
Accepted: 10 December 2014
Published: 21 January 2015
The increased selection pressure of the herbicide glyphosate has played a role in the evolution of glyphosate-resistance in weedy species, an issue that is becoming a threat to global agriculture. The molecular components involved in the cellular toxicity response to this herbicide at the expression level are still unidentified.
In this study, we identify the protein kinase GCN2 as a cellular component that fosters the action of glyphosate in the model plant Arabidopsis thaliana. Comparative studies using wild-type and gcn2 knock-out mutant seedlings show that the molecular programme that the plant deploys after the treatment with the herbicide, is compromised in gcn2. Moreover, gcn2 adult plants show a lower inhibition of photosynthesis, and both seedlings and adult gcn2 plants accumulate less shikimic acid than wild-type after treatment with glyphosate.
These results points to an unknown GCN2-dependent factor involved in the cascade of events triggered by glyphosate in plants. Data suggest either that the herbicide does not equally reach the target-enzyme in a gcn2 background, or that a decreased flux in the shikimate pathway in a gcn2 plants minimize the impact of enzyme inhibition.
KeywordsGlyphosate Gcn2 Transcriptomic Shikimate Translation Herbicide
Since its introduction in 1974, glyphosate has become the world most widely used herbicide, especially after the emergence of transgenic resistant crops in 1996. In 2007, more than 80% of the transgenic crops worldwide were engineered to be glyphosate-resistant . An increase in the application frequency of the herbicide has, however, played a role in the evolution of glyphosate-resistance in weedy species, an issue that is becoming a threat to global agriculture. At least three of the ten most conspicuous weeds have evolved resistant to glyphosate after one decade of transgenic crops . In this scenario, understanding the molecular responses to glyphosate can be useful for future biotechnological approaches.
Resistance to herbicides can be achieved by alterations in the gene encoding the target protein, causing a reduction in the efficacy of the herbicide (target-site resistance) or by any other mechanism independent of the target enzyme (non-target-site resistance, NTSR) . NTRS, that includes decreased herbicide penetration into the plant, decreased rate of herbicide translocation and increased rate of herbicide sequestration/metabolism, has been reported to be the most widespread type of resistance to glyphosate . Although several of these mechanisms have been proposed [5,6], the molecular components involved are still unidentified.
Glyphosate affects plants systemically after application to the leaf surface. The phytotoxic symptoms develop slowly, plant death requiring days or weeks depending on the dose applied . Inhibition of its target enzyme 5-enolpyruvylshikimate-3-phosphate-synthase (EPSPs; EC 220.127.116.11), inhibits the shikimate pathway, leading to a shortage in aromatic amino acids, quinones and cofactor biosynthesis. This is considered by some authors as the main cause of glyphosate toxicity, consistent with the slow development of symptoms . In contrast, others consider that shikimate accumulation, due to a decrease in feedback inhibition through the pathway, leads to an energy drain imposed by a utilization of one phosphoenolpyruvate (PEP) molecule and one ATP molecule for every molecule of shikimate-3-phosphate accumulated and 3 ATP’s for every NADPH [8,9]. Herbicides inhibiting amino acid biosynthesis also induce non-target indirect effects, as proteolysis and an increase in free amino acids . The last omics studies [11,12] reveal that the full picture of molecular disturbances after EPSP inhibition is complex and far from being totally understood. One of the biochemically best characterized cellular events after glyphosate treatment is the rapid shutdown of photosynthesis . Cessation of carbon fixation, decrease in chlorophyll content and electron transport has been reported to occur soon after herbicide application. This impairment is not only a metabolic perturbation but involves gene expression and protein contents. Turfgrass and soybean plants exposed to glyphosate repressed most genes related to photosynthesis [14,15] and repression of photosynthetic proteins was observed in rice after glyphosate treatment . Besides photosynthesis, other cellular processes such as cell cycle  cell motility (in bacteria) , cell death and redox homeostasis  are directly or indirectly affected by glyphosate.
In yeast, amino acid starvation is followed by activation of the protein kinase GCN2 by uncharged tRNAs. This enzyme phosphorylates the α subunit of eukaryotic translation initiation factor 2 (eIF2α), inhibiting the conversion of eIF2γ-GDP to eIF2γ -GTP, preventing further cycles of translation initiation and suppressing protein synthesis . Phosphorylation of eIF2α not only causes a general reduction of protein synthesis, but initiates the selective translation of GCN4, an mRNA containing short open reading frames (ORF) upstream of the long, protein-coding ORF. Its translation upon amino acid starvation produces a transcription factor that activates amino acid biosynthesis genes, helping the cell to recover from the stress. This regulatory response is known as general amino acids control (GAAC) .
Plants contain a GCN2 homologue kinase that complements the gcn2 yeast mutant strain, indicating that the GAAC response could be also operating in plants . In Arabidopsis plants with the aromatic amino acids biosynthesis impaired as imposed by glyphosate treatment, eIF2α is phosphorylated, and this phosphorylation is GCN2-dependent, as it was abolished in an insertion line in the GCN2 gene . Phosphorylation of eIF2α by GCN2 has been observed under several abiotic stresses [20,21], and a GCN2-dependent translational arrest has been observed after treatment with chlorsulfuron, an inhibitor of branched amino acid biosynthesis . Despite these data, the importance of GCN2 pathway as a regulatory mechanism in plants is still under debate [22-24]. The potential homologous gene of the yeast GCN4 has not been found in Arabidopsis, and little evidence was found for the involvement of AtGCN2 in the regulation of expression of amino acid biosynthesis genes, suggesting that this kinase could be playing another role than just regulating amino acid biosynthesis . The role of GCN2 in the plant response to glyphosate remains elusive.
The transcriptional response to glyphosate has been determined in the bacteria Escherichia coli  in the grass Festuca  and in soybean  and in all cases a large number of genes change in expression in response to the herbicide. In the model plant A. thaliana, however, only a few genes (sixteen) were detected as modulated by glyphosate in leaves . This discrepancy, and the unexpected behavior of aAtGCN2 in the regulation of amino acid biosynthetic genes, prompted us to reinvestigate the transcriptional response to glyphosate in Arabidopsis and to give some insights into the role of GCN2 in the triggering of this response.
In this study we show that glyphosate treatment triggers a complex transcriptional response in Arabidopsis plants. Surprisingly, many of these responses are not triggered or are altered in a gcn2 mutant line. We also show that shikimate accumulation in gcn2 plants after herbicide treatment is lower than in wild-type plants. All these results indicate that GCN2 is an important factor in the response of plants to glyphosate and that this protein kinase fosters the action of the herbicide by some unknown mechanism.
Glyphosate treatment causes a dramatic shift of the Arabidopsis transcriptome
Selected categories enriched in Arabidopsis seedlings after glyphosate treatment (for a complete list, see Additional file 2 )
Selected GO categories enriched in glyphosate-induced genes
Selected GO categories enriched in glyphosate-repressed genes
GO biological process
GO biological process
Response to wounding
Photosynthetic electron transport chain
Response to bacterium
Porphyrin biosynthetic process
Aromatic amino acid metabolic process
Response to osmotic stress
Electron transport chain
Response to oxidative stress
Fixation of carbon dioxide
Regulation of cell size
Response to drug
Most remarkable, the herbicide-treated plants show a dramatic down-regulation of the photosynthesis (7.06×10−29 adj. p-value), including chlorophyll biosynthesis (3.89×10−07), electron transport chain (1.26×10−08), and, to a lesser extent, CO2 fixation (4.68×10−03) (Table 1 and Additional files 3 and 4), a fact that was already observed at a later time-point (5 days after treatment) by Cebeci and Budak , in turtgrass. As examples, protochlorophyllide oxidoreductase (AT1G03630), involved in chlorophyll biosynthesis, is repressed more than 9 times; and PsbQ (AT1G14150) and PSB29 (AT2G20890), part of the PSII, are repressed more than 4 times. Interestingly, categories related to translation (including genes for the translation elongation factor EF1B (AT1G64510) and some ribosomal proteins (AT2G24090, AT1G07320 or AT1G48350, among others, repressed more than 4 times), growth and cell division (including genes encoding cyclin A1;1 (AT1G44110) or cyclin D1;1(AT1G70210), repressed more than 4 times) are also down-regulated [see Additional file 4], probably indicating that plant metabolism is being reprogrammed to cope with the stress situation .
The presence of the protein kinase Gcn2 is necessary for the deployment of early cellular responses after glyphosate treatment
GCN2 is a protein kinase that phosphorylates the α subunit of the eIF2 translation initiation factor, a key regulatory mechanism that arrest general protein synthesis and allows the re-stablishment of homeostasis in eukaryotes after several stress conditions . In Saccharomyces cerevisiae, GCN2 is activated under amino acid starvation and triggers a general translational arrest but also promotes the selective translation of the transcription factor GCN4. This factor activates transcription of several hundred genes, including those involved in amino acid biosynthesis . As glyphosate is known to block AAA biosynthesis, we wanted to know whether a similar mechanism was operating in plants, and also to identify other cellular responses to glyphosate eventually regulated by GCN2. We compared the transcriptome of 16-day-old Arabidopsis wild-type seedlings with that of gcn2 mutant seedlings of the same age , both at 6 h after glyphosate treatment. As described in Zhang et al.  and shown in Additional file 5: Figure S2, gcn2 plants are phenotypically indistinguishable from wild-type plans at this and all stages of growth. Previously, we had compared the transcriptome of the same plants under normal conditions and determined that only 24 genes were changing their expression using the selected criteria. [see Additional files 1 and 3]. As shown in Additional file 6: Figure S1, and as reported by Zhang et al. , eIF2α phosphorylation indicates that, after 6 h treatment, GCN2 kinase is activated in wild-type plants.
Selected categories enriched in wild-type vs. gcn2 Arabidopsis seedlings after glyphosate treatment (for a complete list, see Additional file 2 )
Selected GO categories enriched in wild-type as compared with gcn2 in glyphosate-induced genes
Selected GO categories enriched in gcn2 as compared with wild-type in glyphosate –repressed genes
GO biological process
GO biological process
Response to wounding
Response to bacterium
Chlorophyl biosynthetic process
Response to salt stress
Photosynthesis electron transport
Response to osmotic stress
Response to fungus
Regulation of cell size
Electron transport chain
Response to oxidative stress
M phase of cell cycle
gcn2 plants show less glyphosate-derived effects
Effect of glyphosate application on the photosynthetic rate (A N ), stomatal conductance (g s ), substomatal CO 2 concentration (C i ), transpiration rate (E), quantum efficiency of photosystem II (PhiPS2), maximum quantum yield efficiency (Fv/Fm) and SPAD index
A N (μmol m −2 s −1 )
g s (mol m −2 s −1 )
C i (μmol mol −1 )
E (mmol m −2 s −1 )
By contrast, the application of the herbicide has negligible effects on photosynthetic measurements in the gcn2 mutant plants. No significant differences were observed between mock-treated and glyphosate-treated plants on CO2 assimilation, stomatal conductance, quantum efficiency of photosystem II or maximal quantum yield (Table 3). These results further confirm the dependence on GCN2 of the glyphosate-induced repression of photosynthesis.
These results agree with the relative expression of several photosynthetic genes in both genotypes. We analyzed by RT-PCR the expression of the three genes most repressed by glyphosate in the seedlings experiment, namely two light-harvesting complex genes (LHC 2.2 and LHC 4.2) and a gene involved in the assembly of cytochrome b6f (At5g36120). These three genes were not repressed in gcn2 seedlings [see Additional file 3]. As shown in Figure 1E, expression of these genes was repressed in wild-type plants after glyphosate treatment, but the effect was not observed in the gcn2-treated plants, confirming the results obtained in seedlings.
DAB staining is used as an efficient method to detect hydroxen peroxide accumulation in plant tissues  and H2O2 increases have been described in plants treated with glyphosate . Young leaves of wild-type plants after glyphosate treatment were stained using DAB. As shown in Figure 1D, a dramatic accumulation of H2O2 is observed, preceeding the cell death observed in the following days. As expected, gcn2 plants of the same age did not show H2O2 accumulation, further confirming the dependence on GCN2 of the oxidative-stress burst observed after glyphosate treatment.
In summary, gene expression data, as well as glyphosate effects such as photosynthesis decay and oxidative stress, suggest that gcn2 plants are less affected by glyphosate than wild-type plants.
Shikimic acid accumulation is compromised in gcn2 plants after glyphosate treatment
Glyphosate inhibition of the enzyme 5-enolpyruvyl-shikimate-3-phosphate synthase leads to reduced feedback inhibition of the pathway, resulting in carbon flow to shikimate-3-phosphate, which is converted into high levels of shikimate . Given that gcn2 plants were less prone to glyphosate effects, as observed in the experiments described above, we wanted to determine if the target enzyme of glyphosate was inhibited in the same way in both lines. Arabidopsis wild-type and gcn2 seedlings were treated with glyphosate as explained above, and shikimic acid accumulation was measured three days after treatment. The same experiment was performed in four-weeks old plants, in the way explained above and in Methods. As observed in Figure 2, basal levels of shikimic acid are similar in both lines, either in seedlings or in adult plants. However, 3-days after glyphosate application, shikimic acid in wild-type seedlings has increased almost 20 times over the basal levels (46 times in adult plants treatments), whereas gcn2 only accumulates half the amount of shikimic acid that wild-type does (only 10 times in seedlings, 20 times in adult plants). This trend indicates that the target enzyme of glyphosate is not being inactivated at the same extend in gcn2 plants, or that the metabolic flux to the pathway has been reduced in gcn2 plants.
Glyphosate uptake is not affected in gcn2 plants
In the plant, EPSPS is encoded in the nucleus and translocated to the chloroplast, where the aromatic amino acids are synthesized . The highest shikimate accumulation occurs in the youngest tissues, suggesting that glyphosate enters the plant and then translocates to the active meristems to reach the target site, where its action takes place preferentially . Several reasons could then account for the less glyphosate action observed in gcn2 plants. In one hand, glyphosate uptake and/or translocation could be compromised. Less concentration of active glyphosate will be present in the plant chloroplasts, explaining the apparent resistance of gcn2 plants to the herbicide. On the other hand, detoxification mechanisms could be exacerbated in the mutant, inactivating the herbicide.
The mechanisms of glyphosate uptake into plant cells are not well understood. At least at low concentrations of herbicide (in the micromolar range, as it is used in this study) it seems to involve a phosphate transporter, as glyphosate uptake is inhibited by sodium phosphate and phosphonoformic acid, a competitor inhibitor of phosphate transport in plants [33,34]. In order to ascertain if gcn2 plants were compromised in phosphate transport, that could confer to these plants an advantage in the presence of the herbicide, we performed germination assays of wild-type and gcn2 seeds in media with phosphate deficiency. No differences were observed between wild-type and gcn2 seedling when grown in media containing 0-500 μM phosphate (data not shown), indicating that gcn2 plants are taking up phosphate in a similar way as wild-type does. The same results were observed in root-growth assays (data not shown).
Little is known about the molecular events that contribute to non-target-site based resistance (NTSR), the combination of mechanisms that limit to a non-lethal dose the amount of herbicide reaching the target-site. In this study, we identify GCN2 as a cellular component that fosters the action of glyphosate in the model plant A. thaliana. GCN2 is a conserved protein kinase responsible for the phosphorylation of the initiation factor eIF2-α after a number of stress situations. In Arabidopsis, herbicide treatments, wounding, cold treatments, UV or purine starvation has been described to activate AtGCN2  the only kinase able to phosphorylate eIF2-α in this species . Phosphorylation of eIF2-α prevents further cycles of protein translation, and it is assumed that this helps the cell to conserve metabolic resources until it has overcome the immediate biological impact of the stress . In this model, the activation of GCN2 should then be beneficial for the plant to cope with the stress. However, we have shown that the presence of GCN2 is somehow facilitating the action of the herbicide, and cellular responses to glyphosate are not triggered or attenuated in a mutant line that is not able to phosphorylate eIF2-α. This is not the first report where the lack of GCN2 is conferring an advantage when a particular stress is applied. In yeast, ScGCN2 acts also as a negative factor, conferring toxic effects on growth under NaCl stress, and being a gcn2 knock-out strain able to grow normally under 400 mM NaCl . In human tumor cells, HsGcn2 was shown to have an unexpected proapoptotic effect under glucose deficiency stress, and gcn2 knock-out cells are able to survive more than wild-type under these stress conditions . As such, these stresses and the canonical amino acid starvation may utilize distinct pathways that converge on eIF2α phosphorylation with opposing biological outcomes . It will be interesting to know whether a gcn2 background confers resistance to other abiotic stresses in Arabidopsis, and to investigate if this resistance is converging in a single factor for all stresses or in different factors for every one of them.
Known responses to glyphosate include the rapid repression of photosynthesis. Inhibition of CO2 assimilation and depletion of intermediates of the carbon reduction cycle had been documented years ago  and these effects had been attributed to an upregulated flux into the shikimate pathway due to depletion downstream of EPSPs. Recent transcriptomic assays revealed that this repression could be also genetically regulated. In a comparative study of Festuca species, Cebeci and Budak reported that, 5 days after treatment with glyphosate, a marked repression of photosynthetic genes, including chlorophyll biosynthesis, photosystems and Calvin cycle enzymes was occurring . Consistent results were obtained in proteomic assays performed in rice . We have also observed repression in gene expression in adult plants 15 days after treatment, consistent with a clear inhibition of photosynthetic rate. The dramatic repression of photosynthetic genes only 6 h after treatment to Arabidopsis seedlings, when no visible symptoms of leaf chlorosis were observed yet (Table 1), suggests that besides a likely photosynthetic decrease due to metabolic toxicity, an early genetically programmed inhibition of photosynthesis is also taking place. Consistent with our data, repression in some photosynthetic genes was observed in soybean sensitive plants only 4 h after treatment with glyphosate . This decay of photosynthesis after glyphosate treatment is not observed in gcn2 plants, nor is the oxidative stress that characterizes the herbicide effect (Table 2). Moreover, enzymes known to be involved in the detoxification of herbicides such as ABC transporters, glutathione-transferases or glycosyltransferases [38,39] are dramatically activated after glyphosate treatment in wild-type plants, but unaltered or weakly activated in gcn2 plants [see Additional file 1]. Finally, shikimate accumulation in gcn2 plants compared to wild-type clearly demonstrates that the lack of GCN2 becomes an advantage when plants are treated with this herbicide.
At present, these data do not allow to discern at what point the herbicide action benefits from the presence of GCN2 in the cell. Glyphosate is said to starve cells of aromatic amino acids, due to inhibition of the shikimate pathway and the production of chorismate, the precursor of phenylalanine, tryptophan and threonine. However, total amino acid pools after glyphosate treatment rather increase than decrease [40,41] after the first days after treatment, as a consequence of the induction of some proteolytic activities  or are not much different in glyphosate resistant soybean as compared to sensitive cultivars . The known mechanism of activation of GCN2 through uncharged t-RNAs  also works in Arabidopsis  and although not described in detail for glyphosate, treatment with chlorsulfuron, that blocks valine, leucine and isoleucine biosynthesis  yield a peak of eIF2-α phosphorylation 6 h after treatment but return to basal levels after 24 h. Starvation of amino acids and prolonged protein translational arrest via GCN2 are then unlikely to be the major cause of the slow effects of glyphosate treatment. The involvement of GCN2 in the glyphosate mode of action should then fall on the first hours after the treatment, this conditioning the final effect in the plant. How could GCN2 foster the action of the herbicide? As stated above, the early activation of the kinase in wild-type plants after glyphosate application, likely due to the initial decay in aromatic amino acids, do not avoid the expression of cellular factors involved in detoxification of xenobiotics. In some species, vacuolar sequestration is contributing to the resistance mechanism in resistant variants . If the same mechanism is working in Arabidopsis, activation of GCN2 could impair vacuolar membrane trafficking through inhibition of some important protein. Alternatively, the selective translation of certain mRNAs with uORFs in the leader sequence, in the same way of ScGCN4, in yeast, and HsATF4, in humans , could facilitate translocation of glyphosate to the young tissues, where the target enzyme is mainly expressed and the action of glyphosate is more dramatic . Reduced translocation to meristematic sinks is a major mechanism of resistance in horseweed , Italian ryegrass  or hairy fleabane . However, no homologous sequence of ScGCN4 has been found in Arabidopsis so far, and a GCN2-dependent selective translation of mRNAs is unknown. Finally, the absence of GCN2 activity in the mutant gcn2 line could provide a constitutive advantage in the mutant background that diminished herbicide effects, independently of the post-treatment GCN2 activation. One possibility was a higher uptake of glyphosate in wild-type plants. Although the mechanism of glyphosate uptake into plant cells is not well understood, the involvement of a phosphate transporter has been proposed [32,47]. gcn2 plants are not more sensitive than wild-type to growth in phosphate deficiency media, and microarray experiments over non-treated seedlings [see Additional file 1] did not reveal differences in gene expression that make suspect of phosphate transport missregulation, indicating that phosphate transport is not compromised in gcn2 plants. If phosphate transporters are involved in glyphosate uptake in Arabidopsis, then gcn2 plants should be taking up the herbicide at the same rate as wild-type plants.
Several attempts have been made to find genes involved in glyphosate resistance using mutant collections of A. thaliana [48,49]. In Brotherton et al., the same concentration of glyphosate used in this study was used in a germination assay to find EMS-mutagenized mutant lines of Arabidopsis resistant to the herbicide, but no resistant mutant was recovered . If glyphosate resistance single mutations were common, they should have been found in these saturation mutagenesis studies . The lack of GCN2 does not confer resistance to glyphosate in germination assays (data not shown). In a seedling assay performed by Zhang et al., the gcn2 mutant line showed sensitivity to glyphosate treatment . However, we were not able to find this sensitivity in a similar experiment using plants in the same developmental stage (data not shown). May be subtle differences in the experimental conditions are the cause of this apparent discrepancy. A gcn2 knock-out mutant conferring resistance or sensitivity to the same stress has been already reported in animal cell cultures, GCN2 acting as a molecular switch that shifts cells from a proapoptotic to a cytoprotective state in response to glucose deficiency , depending on the duration of the stress. In our experimental setup, we have shown that shikimate is not accumulating and gene activation is not being triggered at the same rate than in wild-type plants, indicating that this protein kinase could be an important clue to discover components involved in the resistance to this herbicide. “Glyphosate is as important to world agriculture as penicillin is to human health” stated Stephen Powles, director of the Australian Herbicide Resistance Initiative (http://www.ahri.uwa.edu.au). Given the spread of glyphosate resistance weeds around the planet and the economic importance for agriculture, understanding the mechanisms of such resistance could help to design new biotechnological approaches for a more efficient use of this important herbicide.
Several mechanisms have been proposed for non-target based resistance to glyphosate, but any study has so far identified any gene that could be directly involved or influencing the final effect of this herbicide in the plant. We have demonstrated that the translational regulator GCN2 is fostering the action of the herbicide by an unknown mechanism. The loss-of-function gcn2 mutant in the model plant A. thaliana emerge as an important tool to decipher the way glyphosate enters the plant and reach its target site. This information will help to design new strategies to preserve the use of glyphosate in the emerging glyphosate-resistance-weeds era.
Plant growth and treatments
A. thaliana accession Landsberg erecta (Ler-0) was used in this study. Genetrap GT8359, containing a Ds transposable element interrupting the GCN2 gene  was obtained from Cold Spring Harbor Laboratory, New York (http://genetrap.cshl.edu/).
Seeds were pretreated in 70% ethanol for 20 min, surface-sterilized in 2.5% bleach for 10 min, and washed with distillated water at least five times. After stratification at 4°C in the dark for 5 days, seeds were sown on 1% agar-containing MS Salts, 1% sucrose, pH 5.5, and grown at 23°C with a 16-h-light/8-h-dark cycle.
For experiments at seedling stage, 16-day-old plantlets (see Additional file 5: Figure S2 for pictures) were submerged in 200 μM glyphosate (SIGMA) or distilled water (mock) for 1 minute, and incubated for further growth on liquid MS medium, 1% sucrose, pH 5.5, under the same conditions. Samples for microarray experiments were collected 6 h after treatment. Samples for shikimate assays were collected 72 h after treatment. Visual inspection was observed during the next two weeks after treatment.
For experiments in adult plants, 10 to 15-day-old plantlets were transferred to soil and grown at 23°C under short-day conditions (8-h-light/16-h-dark). Glyphosate treatments were done three-weeks later (see Additional file 5: Figure S2 for pictures). Plants were sprayed with 200 μM glyphosate once (or distilled water for mock treatments), using a standard sprayer by applying three pulses to every plant at a distance of 10–15 cm from the rosette (aprox 250 μL per plant), and incubated for further growth under the same conditions. Expression analysis by RT-PCR, DAB staining and photosynthetic measurements were done two-weeks after treatment. Visual inspection was followed during the next four weeks after treatment.
Total RNA of glyphosate- and mock-treated 16-day-old seedlings was extracted using RNeasy kit (Qiagen). 1.5 μg of total RNA was labeled using MessageAmp II amplification kit from Life Technologies, following manufacturer instructions. Labeling was done using Cy3 and Cy5 dyes from GE (RPN5661). Before hybridization, slides were pre-hybridized at 42°C for 45 m in 5xSSC, 0.1%SDS and 0.1 mg/mL BSA. Microarray hybridizations and washings were done in manual chambers at 42°C, according to Forment et al. . Scanning and Image Analysis was performed using GenePix Pro 6.0 software (Molecular Devices).
For wild type vs. wild type + Glyphosate, two biological replicates were done. Expression ratios in both microarrays were averaged and considered for further analysis if equal trend (induction or repression) was observed in both replicates. Gene set Enrichment Analysis was done using Fatiscan  taking as significant those categories with and adjusted p-value lower than 0.05.
For gcn2 mutant + Glyphosate vs. wild type + Glyphosate, three biological replicated were done swapping the dyes in one of them. A gene was considered differentially expressed if average fold-change was higher than 2 and had a FDR < 5% after a SAM test . Functional analysis was done using FatiGO  taking as significant those categories with and adjusted p-value lower than 0.05. The same analysis was performed for gcn2 mutant vs. wild type without glyphosate.
These microarrays data have been included in the GEO Omnibus database with the reference numbers GSE 56146 and GSE 56147.
Shikimate determination was done 3 days after treatment in seedlings and 7 days after treatment in adult plants. For the seedling assay, fresh weight was annotated before freezing. From each adult leaf, three discs (4 mm diameter using a micropunch) were placed in a 2-ml tube. Seedlings and discs were frozen with liquid nitrogen and kept at −80 until use. Extraction of shikimate was performed as described in Koger et al. . Vials were removed from the freezer and 0.25 M HCl was added to each vial (1 ml per 100 mg FW or 100 μl per leaf disc). Vials were mixed by vortexing and incubated at room temperature for 1.5 h. Afterwards the solution was frozen and kept at −20°C until analysis. Shikimate was analyzed by HPLC as described before .
Gas exchange and chlorophyll fluorescence measurements were performed as described by Flexas et al. . Instantaneous determinations of net CO2 assimilation rate (AN), stomatal conductance (gs), transpiration rate (E) and substomatal CO2 concentration (Ci) were carried out at steady-state conditions under saturating light (1000 μmol m−2 s−1), a vapour pressure difference (vpd) between 1 and 2 kPa and 400 ppm CO2 with a LI-6400 (LICOR, Nebraska, USA). The actual photochemical efficiency of photosystem II (PhiPS2) was determined by measuring steady-state fluorescence (Fs) and maximum fluorescence (Fm’) during a light-saturating pulse (8000 μmol m−2 s−1) . Maximal photochemical efficiency (Fv/Fm) on dark adapted leaves was measured with a MINI PAM fluorometer (Walz, Effeltrich, Germany). SPAD values were measured with a chlorophyll meter SPAD-502 (Konica Minolta, Osaka, Japan). One measurement per plant was taken on the 8th to 10th leaf from the apex, and for each genotype and treatment, 8 different plants were measured.
In situ detection of hydrogen peroxide was performed by staining with diaminobenzidine (DAB) staining, according to Daudi et al.  with modifications. Briefly, rosette leaves were incubated in staining buffer (1 mg/mL DAB containing Tween 20 (0.05% v/v) and 50 mM sodium phosphate buffer (pH 3.8) and vacuum infiltrated applying 3 pulses of 1.5 m, and stained for 24 h at room temperature. Leaves were fixed in ethanol:glycerol:acetic acid 3:1:1 (bleaching solution) placed in a water bath at 95°C for 15 m. Bleaching solution was replaced and plants were visualized under white light and photographed.
For RT-PCR experiments, total RNA was extracted using RNeasy kit (Qiagen) and treated with DNase I to remove genomic DNA. cDNA was obtained using the Maxima First Strand cDNA Synthesis Kit (Fermentas). Quantitative real-time PCR was performed in a 7500 Fast Real-Time PCR System, from Applied Biosystems, using EvaGreen as a fluorescent reporter and Taq polymerase (Biotools). Primers were designed using PRIMER3 software. Actin 8 (At1g49240) was used as an internal control (Fw 5′-AGTGGTCGTACAACCGGTATTGT; Rv 5′- GAGGATAGCATGTGGAAGTGAGAA). Primers for LHCB 4.2 (Fw 5′- CCACTCTTGGCGCTATCAC; Rv 5′- GCCGATCACTAACACTTCGAT). Primers for LHCB 2.4 (Fw 5′- AGCGACCTCATCCAAAAGG; Rv 5′- TCCGAGAATGGTCCCAAGTA). Primers for B6F (Fw 5′- AGTGACCACCAGCTTCGTCT; Rv 5′- AAGAGACGTGGATCGATTGC). The reaction commenced at 95°C for 5 m, followed by 40 cycles of 15 s at 95°C, 30 s at 55°C, and 30 s at 72°C. Data were analyzed using 7500 Applied Biosystem proprietary software v.2.0.4.
Availability of supporting data
The data sets supporting the results of this article are available in the GEO repository (GSE56146 and GSE56147).
This work was mainly supported by the Universidad Politecnica de Valencia (PAID2011-16) and the Ministerio Español de Ciencia y Tecnología (BFU2011-22526). The work was partially supported through a grant from the Ministerio Español de Ciencia y Tecnología (AGL-2010-18621).
- James C. Global status of Commercialized Biotech/GM Crops. 2010. www.isaaa.org
- Basu C, Halfhill MD, Mueller TC, Stewart Jr CN. Weed genomics: new tools to understand weed biology. Trends Plant Sci. 2004;9(8):391–8.PubMedView ArticleGoogle Scholar
- Délye C. Unravelling the genetic bases of non-target-site-based resistance (NTSR) to herbicides: a major challenge for weed science in the forthcoming decade. Pest Manag Sci. 2013;69(2):176–87.PubMedView ArticleGoogle Scholar
- Powles SB, Yu Q. Evolution in action: plants resistant to herbicides. Annu Rev Plant Biol. 2010;61:317–47.PubMedView ArticleGoogle Scholar
- Ge X, d’Avignon DA, Ackerman JJ, Collavo A, Sattin M, Ostrander EL, et al. Vacuolar glyphosate-sequestration correlates with glyphosate resistance in ryegrass (Lolium spp.) from Australia, South America, and Europe: a 31P NMR investigation. J Agric Food Chem. 2012;60(5):1243–50.PubMedView ArticleGoogle Scholar
- Ge X, d’Avignon DA, Ackerman JJ, Sammons RD. Rapid vacuolar sequestration: the horseweed glyphosate resistance mechanism. Pest Manag Sci. 2010;66(4):345–8.PubMed CentralPubMedGoogle Scholar
- Duke SO, Powles SB. Glyphosate: a once-in-a-century herbicide. Pest Manag Sci. 2008;64(4):319–25.PubMedView ArticleGoogle Scholar
- Fisher RS, Berry A, Gaines CG, Jensen RA. Comparative action of glyphosate as a trigger of energy drain in Eubacteria. J Bacteriol. 1986;168:1147–54.Google Scholar
- De María N, Becerril JM, García-Plazaola JI, Hernández A, De Felipe MR, Fernández-Pascual M. New insights on glyphosate mode of action in nodular metabolism: Role of shikimate accumulation. J Agric Food Chem. 2006;54:2621–8.PubMedView ArticleGoogle Scholar
- Zulet A, Gil-Monreal M, Villamor JG, Zabalza A, van der Hoorn RAL, Royuela M. Proteolytic Pathways Induced by Herbicides That Inhibit Amino Acid Biosynthesis. PLoS One. 2013;8(9):e73847.PubMed CentralPubMedView ArticleGoogle Scholar
- Ahsan N, Lee DG, Lee KW, Alam I, Lee SH, Bahk JD, et al. Glyphosate-induced oxidative stress in rice leaves revealed by proteomic approach. Plant Physiol Biochem. 2008;46(12):1062–70.PubMedView ArticleGoogle Scholar
- Lu W, Li L, Chen M, Zhou Z, Zhang W, Ping S, et al. Genome-wide transcriptional responses of Escherichia coli to glyphosate, a potent inhibitor of the shikimate pathway enzyme 5-enolpyruvylshikimate-3-phosphate synthase. Mol Biosyst. 2013;9(3):522–30.PubMedView ArticleGoogle Scholar
- Servaites JC, Tucci MA, Geiger DR. Glyphosate effects on carbon assimilation, ribulose bisphosphate carboxylase activity, and metabolite levels in sugar beet leaves. Plant Physiol. 1987;85(2):370–4.PubMed CentralPubMedView ArticleGoogle Scholar
- Cebeci O, Budak H. Global Expression Patterns of Three Festuca Species Exposed to Different Doses of Glyphosate Using the Affymetrix GeneChip Wheat Genome Array. Comp Funct Genomics. 2009;Article ID 505701:10. doi:10.1155/2009/505701.Google Scholar
- Zhu J, Patzoldt WL, Shealy RT, Vodkin LO, Clough SJ, Tranel PJ. Transcriptome response to glyphosate in sensitive and resistant soybean. J Agric Food Chem. 2008;13;56(15):6355–63.View ArticleGoogle Scholar
- Marc J, Mulner-Lorillon O, Bellé R. Glyphosate-based pesticides affect cell cycle regulation. Biol Cell. 2004;96(3):245–9.PubMedView ArticleGoogle Scholar
- Wek RC, Jiang HY, Anthony TG. Coping with stress: eIF2 kinases and translational control. Biochem Soc Trans. 2006;34:7–11.PubMedView ArticleGoogle Scholar
- Hinnebusch AG. Translational regulation of GCN4 and the general amino acid control of yeast. Annu Rev Microbiol. 2005;59:407–50.PubMedView ArticleGoogle Scholar
- Zhang Y, Dickinson JR, Paul MJ, Halford NG. Molecular cloning of an Arabidopsis homologue of GCN2, a protein kinase involved in co-ordinated response to amino acid starvation. Planta. 2003;217(4):668–75.PubMedView ArticleGoogle Scholar
- Zhang Y, Wang Y, Kanyuka K, Parry MA, Powers SJ, Halford NG. GCN2-dependent phosphorylation of eukaryotic translation initiation factor-2alpha in Arabidopsis. J Exp Bot. 2008;59(11):3131–41.PubMed CentralPubMedView ArticleGoogle Scholar
- Lageix S, Lanet E, Pouch-Pelissier MN, Espagnol MC, Robaglia C, Deragon JM, et al. Arabidopsis eIF2a kinase GCN2 is essential for growth in stress conditions and is activated by wounding. BMC Plant Biol. 2008;8:134–42.PubMed CentralPubMedView ArticleGoogle Scholar
- Shaikhin SM, Smailov SK, Lee AV, Kozhanov EV, Iskakov BK. Interaction of wheat germ translation initiation factor 2 with GDP and GTP. Biochimie. 1992;74(5):447–54.PubMedView ArticleGoogle Scholar
- Krishna VM, Janaki N, Ramaiah KV. Wheat germ initiation factor 2 (WG × eIF2) decreases the inhibition in protein synthesis and eIF2B activity of reticulocyte lysates mediated by eIF2alpha phosphorylation. Arch Biochem Biophys. 1997;346(1):28–36.PubMedView ArticleGoogle Scholar
- Immanuel TM, Greenwood DR, MacDiarmid RM. A critical review of translation initiation factor eIF2α kinases in plants-regulating protein synthesis during stress. Funct Plant Biol. 2012;39:717–35.View ArticleGoogle Scholar
- Byrne EH, Prosser I, Muttucumaru N, Curtis TY, Wingler A, Powers S, et al. Overexpression of GCN2-type protein kinase in wheat has profound effects on free amino acid concentration and gene expression. Plant Biotechnol J. 2012;10(3):328–40.PubMedView ArticleGoogle Scholar
- Das M, Reichman JR, Haberer G, Welzl G, Aceituno F, Mader MT, et al. A composite transcriptional signature differentiates responses towards closely related herbicides in Arabidopsis thaliana and Brassica napus. Plant Mol Biol. 2010;72:545–56.PubMed CentralPubMedView ArticleGoogle Scholar
- Yuan JS, Tranel PJ, Stewart Jr CN. Non-target-site herbicide resistance: a family business. Trends Plant Sci. 2007;12(1):6–13.PubMedView ArticleGoogle Scholar
- Peng YH, Abercrombie LLG, Yuan JS, Riggins CW, Sammons RD, Tranel PJ, et al. Characterization of the horseweed (Conyza canadensis) transcriptome using GS-FLX 454 pyrosequencing and its application for expression analysis of candidate non-target herbicide resistance genes. Pest Manag Sci. 2010;66:1053–62.PubMedView ArticleGoogle Scholar
- Baena-González E. Energy signaling in the regulation of gene expression during stress. Mol Plant. 2010;3(2):300–13.PubMedView ArticleGoogle Scholar
- Vivancos PD, Driscoll SP, Bulman CA, Ying L, Emami K, Treumann A, et al. Perturbations of amino acid metabolism associated with glyphosate-dependent inhibition of shikimic acid metabolism affect cellular redox homeostasis and alter the abundance of proteins involved in photosynthesis and photorespiration. Plant Physiol. 2011;157(1):256–68.PubMed CentralPubMedView ArticleGoogle Scholar
- Daudi A, Cheng Z, O’Brien JA, Mammarella N, Khan S, Ausubel FM, et al. The apoplastic oxidative burst peroxidase in Arabidopsis is a major component of pattern-triggered immunity. Plant Cell. 2012;24(1):275–87.PubMed CentralPubMedView ArticleGoogle Scholar
- Schmid J, Amrhein N. The shikimate pathway. In: Singh B, editor. Plant Amino Acids, Biochemistry and Biotechnology. New York: Marcel Dekker; 1999. p. 147–70.Google Scholar
- Denis M-H, Delrot S. Carrier-mediated uptake of glyphosate in broad bean (Vicia faba) via a phosphate transporter. Physiol Plant. 1993;87:569–75.View ArticleGoogle Scholar
- Hetherington PR, Marshall G, Kirkwood RC, Warner JM. Absorption and efflux of glyphosate by cell suspensions. J Exp Bot. 1998;49:527–33.View ArticleGoogle Scholar
- Goossens A, Dever TE, Pascual-Ahuir A, Serrano R. The protein kinase Gcn2p mediates sodium toxicity in yeast. J Biol Chem. 2001;276(33):30753–60.PubMedView ArticleGoogle Scholar
- Muaddi H, Majumder M, Peidis P, Papadakis AI, Holcik M, Scheuner D, et al. Phosphorylation of eIF2α at serine 51 is an important determinant of cell survival and adaptation to glucose deficiency. Mol Biol Cell. 2010;21(18):3220–31.PubMed CentralPubMedView ArticleGoogle Scholar
- Geiger DR, Kapitan SW, Tucci MA. Glyphosate inhibits photosynthesis and allocation of carbon to starch in sugar beet leaves. Plant Physiol. 1986;82:468–72.PubMed CentralPubMedView ArticleGoogle Scholar
- Edwards R, Del Buono D, Fordham M, Skipsey M, Brazier M, Dixon DP, et al. Differential induction of glutathione transferases and glucosyltransferases in wheat, maize and Arabidopsis thaliana by herbicide safeners. Z Naturforsch. 2005;60(3–4):307–16.Google Scholar
- Cummins I, Wortley DJ, Sabbadin F, He Z, Coxon CR, Straker HE, et al. Key role for a glutathione transferase in multiple-herbicide resistance in grass weeds. Proc Natl Acad Sci U S A. 2013;110(15):5812–7.PubMed CentralPubMedView ArticleGoogle Scholar
- Orcaray L, Igal M, Marino D, Zabalza A, Royuela M. The possible role of quinate in the mode of action of glyphosate and acetolactate synthase inhibitors. Pest Manag Sci. 2010;66:262–9.PubMedView ArticleGoogle Scholar
- Orcaray L, Zulet A, Zabalza A, Royuela M. Impairment of carbon metabolism induced by the herbicide glyphosate. J Plant Physiol. 2012;169:27–33.PubMedView ArticleGoogle Scholar
- Li MW, AuYeung WK, Lam HM. The GCN2 homologue in Arabidopsis thaliana interacts with uncharged tRNA and uses Arabidopsis eIF2α molecules as direct substrates. Plant Biol (Stuttg). 2013;15(1):13–8.View ArticleGoogle Scholar
- Ge X, Davignon DA, Ackerman JJ, Sammons D. In Vivo 31P-NMR Studies of Glyphosate Uptake, Vacuolar Sequestration, and Tonoplast Pump Activity in Glyphosate Resistant Horseweed (Conyza canadensis). Plant Physiol. 2014;166:1255–68.PubMedView ArticleGoogle Scholar
- Feng PCC, Chiu T, Sammons RD. Glyphosate efficacy is contributed by its tissue concentration and sensitivity in velvetleaf (Abutilon theophrasti). Pestic Biochem Physiol. 2003;77:83–91.View ArticleGoogle Scholar
- Koger CH, Reddy KN. Role of absorption and translocation in the mechanism of glyphosate resistance in horseweed (Conyza canadensis). Weed Sci. 2005;53:84–9.View ArticleGoogle Scholar
- Perez-Jones AKW, Park J, Colquhoun C, Mallory-Smith J, Shaner D. Identification of glyphosate-resistant Italian ryegrass (Lolium multiflorum) in Oregon. Weed Sci. 2005;53:775–9.View ArticleGoogle Scholar
- Morin F, Vera V, Nurit F, Tissut M, Marigo G. Glyphosate uptake in Catharanthus roseus cells: Role of a phosphate transporter. Pestic Biochem Physiol. 1997;58:13–22.View ArticleGoogle Scholar
- Jander G, Baerson SR, Hudak JA, Gonzalez KA, Gruys KJ, Last RL. Ethylmethanesulfonate saturation mutagenesis in Arabidopsis to determine frequency of herbicide resistance. Plant Physiol. 2003;131:139–46.PubMed CentralPubMedView ArticleGoogle Scholar
- Brotherton JE, Jeschke MR, Tranel PJ, Widholm JM. Identification of Arabidopsis thaliana variants with differential glyphosate responses. J Plant Physiol. 2007;164:1337–45.PubMedView ArticleGoogle Scholar
- Forment J, Gadea J, Huerta L, Abizanda L, Agusti J, Alamar S, et al. Development of a citrus genome-wide EST collection and cDNA microarray as resources for genomic studies. Plant Mol Biol. 2005;57(3):375–91.PubMedView ArticleGoogle Scholar
- Medina I, Carbonell J, Pulido L, Madeira SC, Goetz S, Conesa A, et al. Babelomics: an integrative platform for the analysis of transcriptomics, proteomics and genomic data with advanced functional profiling. Nucleic Acids Res. 2010;38:W210–3.PubMed CentralPubMedView ArticleGoogle Scholar
- Tusher VG, Tibshirani R, Chu G. Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci U S A. 2001;98(9):5116–21.PubMed CentralPubMedView ArticleGoogle Scholar
- Koger CH, Shaner DL, Krutz LJ, Walker TW, Buehring N, Henry WB, et al. Rice (Oryza sativa) response to drift rates of glyphosate. Pest Manag Sci. 2005;61(12):1161–7.PubMedView ArticleGoogle Scholar
- Flexas J, Ortuño MF, Ribas-Carbo M, Diaz-Espejo A, Flórez-Sarasa ID, Medrano H. Mesophyll conductance to CO2 in Arabidopsis thaliana. New Phytol. 2007;175:501–11.PubMedView ArticleGoogle Scholar
- Genty B, Briantais JM, Baker NR. The relationship between the quantum yield of photosynthetic electron transport and quenching of chlorophyll fluorescence. Biochim Biophys Acta. 1989;990:87–92.View ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.