Expression of ROS-responsive genes and transcription factors after metabolic formation of H2O2 in chloroplasts

Glycolate oxidase (GO) catalyses the oxidation of glycolate to glyoxylate, thereby consuming O2 and producing H2O2. In this work, Arabidopsis thaliana plants expressing GO in the chloroplasts (GO plants) were used to assess the expressional behavior of reactive oxygen species (ROS)-responsive genes and transcription factors (TFs) after metabolic induction of H2O2 formation in chloroplasts. In this organelle, GO uses the glycolate derived from the oxygenase activity of RubisCO. Here, to identify genes responding to an abrupt production of H2O2 in chloroplasts we used quantitative real-time PCR (qRT-PCR) to test the expression of 187 ROS-responsive genes and 1880 TFs after transferring GO and wild-type (WT) plants grown at high CO2 levels to ambient CO2 concentration. Our data revealed coordinated expression changes of genes of specific functional networks 0.5 h after metabolic induction of H2O2 production in GO plants, including the induction of indole glucosinolate and camalexin biosynthesis genes. Comparative analysis using available microarray data suggests that signals for the induction of these genes through H2O2 may originate in the chloroplast. The TF profiling indicated an up-regulation in GO plants of a group of genes involved in the regulation of proanthocyanidin and anthocyanin biosynthesis. Moreover, the upregulation of expression of TF and TF-interacting proteins affecting development (e.g., cell division, stem branching, flowering time, flower development) would impact growth and reproductive capacity, resulting in altered development under conditions that promote the formation of H2O2.


INTRODUCTION
Photosynthetic organisms are confronted with reactive oxygen species (ROS), such as singlet oxygen ( 1 O 2 ), the superoxide anion radical (O − 2 ), the hydroxyl radical (OH·), and hydrogen peroxide (H 2 O 2 ), which may cause oxidative stress and damage to important biological molecules (Apel and Hirt, 2004;Møller et al., 2007). Plants in their natural environments are often exposed to sudden increases in light intensity, which results in the absorption of excitation energy in excess of that required for metabolism. In chloroplasts, when absorbed energy is in excess at photosystem II (PSII), O − 2 is produced during the Mehler reaction by Fd-NADPH oxidase at PSI and is dismutated by superoxide dismutase (SOD) to H 2 O 2 (Ort and Baker, 2002;Asada, 2006). The photorespiratory pathway consumes photosynthetic reducing energy and produces H 2 O 2 in the peroxisomes through the action of glycolate oxidase (GO) (Maurino and Peterhansel, 2010). H 2 O 2 is also produced during a variety of different reactions under stress conditions, often through the detoxification of 1 O 2 and O − 2 . The generated H 2 O 2 is scavenged by different antioxidant/enzyme reactions: the ascorbate and glutathione cycles, ascorbate peroxidase (APX), catalase, and peroxiredoxin (PRX) (Tripathi et al., 2009). ROS generated in the chloroplast have been implicated as triggers of signaling pathways that influence expression of nuclearencoded genes, which may initiate responses such as cell death or acclimation depending on the degree of the stress (Karpinski et al., 1999;Fryer et al., 2003;Op den Camp et al., 2003;Danon et al., 2005). H 2 O 2 can take part in signaling acting as messenger either directly (e.g., by reversibly modifying critical thiol groups in target proteins; Neill et al., 2002) or by using an oxidized product as a secondary messenger (Møller et al., 2007). The H 2 O 2 -scavenging enzymes APX and dehydroascorbate reductase (DHAR) may act as highly efficient initiators of oxidative signaling by generating transient bursts of reduced glutathione. This in consequence triggers glutaredoxin-mediated protein oxidation (Neill et al., 2002). Crosstalk between redox pools of different cellular compartments, possibly transmitted by a redox shift in cellular components, has also been suggested to be important for control of the expression of nuclear genes (Baier and Dietz, 2005;Leister, 2005). A generalized model of H 2 O 2 signal transduction pathways suggests that H 2 O 2 may also directly oxidize transcription factors (TFs) in either the cytosol or the nucleus. Alternatively, H 2 O 2 -mediated activation of a signaling protein such as a protein kinase may activate TFs (Mittler et al., 2004;Miao et al., 2007). TFs would interact with cognate H 2 O 2 -response elements in target gene promoters thereby modulating gene expression (Foyer and Noctor, 2005). Recently, Møller and Sweetlove (2010) put forward the hypothesis that H 2 O 2 itself is unlikely to be the signaling molecule that selectively regulates nuclear-encoded chloroplastic genes but rather that oxidized peptides deriving from proteolysis of oxidized proteins would act as second messengers during retrograde ROS signaling. On the other hand, using spin trapping EPR spectroscopy in addition to chemical assays (employing Amplex Red reagent), Mubarakshina et al. (2010) showed that 5% of the H 2 O 2 produced inside chloroplasts at high light intensities can actually be detected outside the organelles. This process may involve the pass of H 2 O 2 through aquaporins (Bienert et al., 2007) and might be sufficient to trigger signaling processes outside the chloroplasts. Desikan et al. (2001) showed that approximately 1% of the transcriptome was altered in H 2 O 2 -treated Arabidopsis thaliana (A. thaliana) cell cultures. Although H 2 O 2 -responsive promoters have been identified (Desikan et al., 2001), specific H 2 O 2regulatory DNA sequences and their cognate TFs have not been isolated and characterized. In more recent studies genes involved in H 2 O 2 signal transduction have been identified or proposed, including mitogen-activated protein kinases (MAPKs), various TFs of e.g., the NAC, ZAT, and WRKY families, miRNAs and others (Van Breusegem et al., 2008;Li et al., 2011;Petrov and Van Breusegem, 2012). Moreover, using genome-wide analysis of catalase deficient A. thaliana, H 2 O 2 was inferred to regulate the expression of genes encoding specific small heat shock proteins, several TFs and candidate regulatory proteins (Vandenabeele et al., 2004;Vanderauwera et al., 2005).
To date, it is not known to which extent the chemical specificity of the ROS and the cellular compartment of their release may contribute to the multiplicity of responses that occur in plants. A major challenge is to dissect the genetic networks that control ROS signaling and to assess specific and common responses toward different types of ROS signals. To this end, the molecular, biochemical and physiological responses of A. thaliana to elevated in planta levels of H 2 O 2 were and are being investigated in various types of model systems including mutants altered in the ROS scavenging machinery (Maurino and Flügge, 2008). However, the analysis of dynamic physiological processes using (knock-out) mutants may not always be straightforward, especially when compensatory cellular mechanisms are induced. With respect to ROS-related mutants, changing the balance of scavenging enzymes may induce compensatory mechanisms such that signaling and oxidative damage effects may not be easily separated. Moreover, invasive experimental setups like the application of oxidative stress-causing agents may induce a nonspecific oxidative stress that acts throughout the cell and triggers additional responses that may complicate the analysis of ROS signal transduction pathways (Maurino and Flügge, 2008). We have recently developed a tool to functionally dissect the action of plastid-generated H 2 O 2 , using plants overexpressing GO in plastids (GO plants;Fahnenstich et al., 2008). During photosynthesis, the oxygenase activity of ribulose 1,5-bisphosphate carboxylase/oxygenase (RubisCO) produces glycolate 2-phosphate within the chloroplasts, which is then dephosphorylated to glycolate by phosphoglycolate phosphatase (Maurino and Peterhansel, 2010). In GO plants, glycolate is oxidized to glyoxylate by the plastidic GO, with the parallel production of H 2 O 2 . When growing under moderate photon fluxes and ambient CO 2 concentration (photorespiratory conditions) the GO plants remain smaller than the wild type, presenting a reduced rosette diameter and yellowish leaves due to H 2 O 2 accumulation (Fahnenstich et al., 2008). In contrast, in non-photorespiratory conditions (e.g., at high CO 2 concentration) the oxygenase activity of RubisCO is abolished and thus, the metabolic flux through GO is suppressed, allowing GO plants to grow like wild type (Fahnenstich et al., 2008). Transferring GO plants from high to ambient CO 2 concentration specifically induces H 2 O 2 formation in the chloroplasts (Fahnenstich et al., 2008). These properties permit the modulation of plastidic produced H 2 O 2 levels by changing light intensity and/or CO 2 levels (Maurino and Flügge, 2008). Moreover, H 2 O 2 is specifically generated without a concomitant accumulation of superoxide or singlet oxygen, which are common precursors of H 2 O 2 during ROS generation in chloroplasts. A similar experimental set-up was employed in previous studies using catalase null mutants in which the production of peroxisomal H 2 O 2 is induced by changing the conditions of plant growth from non-photorespiratory to photorespiratory conditions (e.g., high light intensity) (Dat et al., 2000;Vandenabeele et al., 2004;Vanderauwera et al., 2005). The metabolic production of H 2 O 2 may avoid the pleiotropic effects discussed above but it cannot be ruled out that ROS-unrelated pleiotropic reactions may occur in both approaches due to abrupt changes in CO 2 level or light intensity.
In this work we attempted to identify genes strongly responding to an abrupt production of H 2 O 2 in chloroplasts of A. thaliana. To this end we tested the expressional changes of 187 nuclear-encoded ROS-responsive genes and 1880 TFs, using quantitative real-time (qRT)-PCR (Czechowski et al., 2004;Balazadeh et al., 2008;Wu et al., 2012) upon transfer of high CO 2 -grown GO and wild-type (WT) plants to ambient CO 2 concentration. Our data revealed a rapid and coordinated expression response of ROS-affected genes of specific functional networks in GO including an early induction of indole glucosinolate and camalexin biosynthesis genes and an up-regulation of a group of genes involved in the regulation of proanthocyanidin and anthocyanin biosynthesis. Moreover, the upregulation of expression of TF and TF-interacting proteins affecting development (e.g., cell division, stem branching, flowering time, flower development) would impact growth and reproductive capacity, resulting in altered development under conditions that promote the formation of H 2 O 2 .

PLANT MATERIAL
Arabidopsis thaliana (L.) Heynh. ecotype Columbia-0 (Col-0, wild-type) constitutively expressing glycolate oxidase (GO, At3g14420) in the plastids (GO plants) under the cauliflower mosaic virus 35S promoter were generated in our previous work (Fahnenstich et al., 2008). In these plants to direct the expression of GO to the chloroplats the stromal targeting presequence from Arabidopsis thaliana phosphoglucomutase (At5g51820) was used (Fahnenstich et al., 2008). WT and GO transgenic plants were grown in pots containing 3 parts of soil (Gebr. Patzer KG, Sinntal-Jossa, Germany) and one part of vermiculite (Basalt Feuerfest, Linz, Austria) under a 16 h-light/8 h-dark regime at photosynthetically active photon flux densities (PPFD) of 75 µmol quanta m −2 s −1 at 22 • C day/18 • C night temperatures and a CO 2 concentration of 3000 ppm. After 3 weeks of growth plants were transferred to ambient CO 2 concentration (380 ppm) and the same PPFD. Whole rosettes were harvested at different time points after transfer, immediately frozen in liquid nitrogen and stored at −80 • C until use for RNA isolation and H 2 O 2 measurements.

ISOLATION OF RNA AND REAL-TIME PCR ANALYSIS
For the large-scale qRT-PCR analysis, total RNA was extracted from 100 mg leaves (fresh weight) using RNeasy Plant Mini kit (Qiagen, Valencia, USA) according to the manufacturer's protocol. DNAse I digestion was performed on 20-30 µg of total RNA using TURBO DNase Kit (Ambion, Cambridgeshire, UK) according to manufacturer's instructions. RNA integrity was checked on 1% (w/v) agarose gels and concentration measured with a Nanodrop ND-1000 spectrophotometer before and after DNAse treatment. Absence of genomic DNA was confirmed subsequently by quantitative PCR using primers that amplify an intron sequence of the gene At5g65080 (forward 5 -TTTTTTGCCCCCTTCGAATC-3 and reverse 5 -ATCTTCCGCCACCACATTGTAC-3 ). First-strand cDNA was synthesized from 8 µg to 10 µg of total RNA using RevertAid™ First Strand cDNA Synthesis Kit (Fermentas, St. Leon-Rot, Germany) following the manufacturer's protocol. The efficiency of cDNA synthesis was estimated by qRT-PCR using two different primer sets annealing to the 5 and 3 ends, respectively, of a control gene (At3g26650, GAPDH, glyceraldehyde-3phosphate dehydrogenase). Primer sequences were as follows: for GAPDH3 , forward 5 -TTGGTGACAACAGGTCAAGCA-3 and reverse 5 -AAACTTGTCGCTCAATGCAATC-3 ; for GAPDH5 , forward 5 -TCTCGATCTCAATTTCGCAAAA-3 and reverse 5 -CGAAACCGTTGATTCCGATTC-3 . Transcript levels of each gene were normalized to ACTIN2 (At3g18780) transcript abundance (forward 5 -TCCCTCAGCACATTCCAGCAGAT-3 and reverse 5 -AACGATTCCTGGACCTGCCTCATC-3 ). A total of 187 ROS-responsive genes (Wu et al., 2012) and 1880TFs (Czechowski et al., 2004Balazadeh et al., 2008) were analyzed by qRT-PCR as previously described (Caldana et al., 2007;Balazadeh et al., 2008). PCR reactions were run on an ABI PRISM 7900HT sequence detection system (Applied Biosystems, Darmstadt, Germany), and amplification products were visualized using SYBR Green (Applied Biosystems).
Levels of H 2 O 2 were determined using the Amplex® Red Technology (Life Technologies, Darmstadt, Germany) following the manufacturer's instructions. Amplex Red (N-acetyl-3,7-dihydroxyphenoxazine) reacts with H 2 O 2 in the presence of horseradish peroxidase and forms the fluorescent product resorufin. For the determinations, 100 mg leaves (fresh weight) were ground in liquid nitrogen into a fine powder and resuspended with 0.15 mL extraction buffer prepared as indicated by the manufacturer. This suspension was centrifuged at 4 • C at 13,000 rpm for 15 min. Five µL of the supernatant, 45 µL distilled water and 50 µL of Amplex® Red solution were added to a microtitre plate. After 30 min incubation in the dark fluorescence was measured by excitation at 560 nm and emission reads at 590 nm. A calibration curve was established with known H 2 O 2 concentrations.

GENE EXPRESSION NETWORK ANALYSIS
The two genes that were most strongly induced under photorespiratory conditions in GO plants at the 0.5 and 6 h time points (At3g02840 and At1g17180, respectively) were used as baits to identify globally coexpressed genes using the ATTED-II database (http://atted.jp), which allows evaluating genes that are coexpressed under five experimental conditions (tissue, abiotic stress, biotic stress, hormones, and light conditions) (Obayashi et al., 2009).

INDUCTION OF H 2 O 2 FORMATION IN GO PLANTS
The production of H 2 O 2 in leaves of plants overexpressing GO in the plastids (Fahnenstich et al., 2008) was analyzed after activation of photorespiration by transferring high CO 2 -grown plants to ambient-CO 2 conditions. As shown in Table 1, higher levels of H 2 O 2 were determined in GO than in WT plants at 0.5 and 4 h after transfer while GO and WT plants maintained under non-photorespiratory conditions (3000 ppm CO 2 ) showed similar H 2 O 2 levels at both time points (Table 1). Note, that as the measurements were performed using whole-leaf extracts the expected differences in chloroplastic H 2 O 2 levels between GO and WT plants under photorespiratory condition may be higher than determined here.

EXPRESSION PROFILING OF ROS MARKER GENES IN GO AND WILD-TYPE PLANTS AFTER THE INDUCTION OF H 2 O 2 FORMATION IN CHLOROPLASTS
To study the impact of an abrupt production of H 2 O 2 in chloroplasts on nuclear gene expression, we analyzed transcript level changes of 187 ROS-responsive genes using a previously established qRT-PCR platform (detailed in Wu et al., 2012). The genes Gene expression was analyzed in whole rosettes of 3-week-old WT and GO plants at 0.5, 4, 6, and 12 h after shifting high-CO 2grown plants (non-photorespiratory condition) to ambient CO 2 concentration (photorespiratory condition). Expression profiling was performed in two biological replicates and log-fold change (log2 FC) ratios of expression changes were calculated for GO and WT plants by comparing gene expression levels before and after the CO 2 concentration shift. A total of 131 genes were expressed in all examined samples ( Table A1 in Appendix). The remaining 56 genes did not yield detectable PCR amplicons, indicating no or marginal expression under our experimental conditions.
Considering a 3-fold expression difference cut-off, 120 genes displayed differential expression in GO and/or WT plants upon transfer from high to ambient CO 2 concentration; the vast majority of the affected genes (116 in total) were up-regulated, and only four genes were down-regulated (Figure 1, Table A1 in Appendix). Most noticeably, expression of 58 genes was induced in GO plants already within 0.5 h after the transfer to ambient CO 2 condition, whilst only a single gene was induced in the wild type at the same time point (Figure 1). Importantly, however, many genes showed also high expression in the wild type at later time points after the CO 2 concentration shift, but the expressional changes were in most cases more pronounced in GO than WT plants (Figure 1, and section "Early Induction of Indole Glucosinolate and Camalexin Biosynthesis Genes in GO Plants"). Thus, our data indicate that similar sets of ROS-responsive genes responded to the CO 2 shift in GO and WT plants; however, the dynamics of the transcriptional responses were clearly different in the two types of plants, being faster and more prominent in the GO plants.

EARLY INDUCTION OF INDOLE GLUCOSINOLATE AND CAMALEXIN BIOSYNTHESIS GENES IN GO PLANTS
To identify transcripts responsive to metabolically produced H 2 O 2 we focused our analysis on the 0.5-and 6-h time points. Genes were considered differentially expressed when the fold change was more than 3-fold (log 2 ≥ 1.56).
At 0.5 h after shifting plants to ambient CO 2 concentration, 58 of the 131 expressed genes were induced in GO plants by more than 3-fold, whilst in the wild type the expression change was less than 3-fold, suggesting that these genes participate in early signaling steps triggered by the production of H 2 O 2 under photorespiratory conditions ( Table 2). After 6 h, seven of these genes showed WT levels of expression (below 3-fold), while 29 were further overexpressed only in GO ( Table 2). Although at 6 h after transfer to ambient CO 2 the expression fold-change (FC) of the remaining 22 genes was higher than 3 in both, GO and WT plants, the expression change between GO and WT (FC GO /FC WT ) was higher than 2 for 16 of these genes ( Later responding genes, which were affected only after 6 h under photorespiratory conditions, were also identified. From the 23 genes that showed an expression change of above 3-fold in GO, 13 were only induced in GO, while 10 genes were induced in both, GO and WT. The FC ratio in GO and WT (FC GO /FC WT ) was above 2 for the 10 genes ( Table 3), indicating that their expression in GO plants is controlled by the higher levels of H 2 O 2 , similar to the early-responsive genes.   The most highly up-regulated gene in GO plants at 0.5 h after induction of H 2 O 2 production was At3g02840 (encoding a putative U-box-type E3 ubiquitin ligase, known to respond immediately-early to fungal elicitation) ( Table 2). We used the ATTED-II database (http://atted.jp; Obayashi et al., 2009) to discover genes coexpressed with At3g02840 and observed that 45 of the 58 genes induced at 0.5 h after induction of H 2 O 2 production cluster together (Table 2), indicating that metabolically produced H 2 O 2 in GO plants induces the coordinate expression of functionally related genes. A similar analysis using the most highly expressed gene at 6 h after induction of H 2 O 2 production (At1g17180, encoding glutathione S-transferase Tau 25) indicated that another group of eight genes are coordinately expressed in GO plants at this later time point (Table 3). Recently, Inzé et al. (2012) listed the 85 most strongly H 2 O 2 -responsive genes in catalase loss-of-function mutants shifted from low-to high-light conditions, where H 2 O 2 is produced in peroxisomes by the action of photorespiratory GOs. Interestingly, 23 of the 81 genes, which changed their expression in the GO plants were also differentially expressed in catalase loss-of-function mutants (Tables 2 and 3), indicating that they respond to enhanced levels of H 2 O 2 independent of the site of its generation; the remaining genes may then represent candidates preferentially responsive to H 2 O 2 produced in chloroplasts. Many of the genes up-regulated in GO plants encode proteins or TFs of currently unknown specific functions. Interestingly, however, several of the early-responsive genes are involved in tryptophan-derived biosynthesis of the phytoanticipins camalexin and indole glucosinolates, i.e., secondary metabolites that have antifungal and insect-deterring functions (Kliebenstein et al., 2001;Bednarek et al., 2009). These genes encode (1) the transcription factor WRKY33 (At2g38470), which is involved in controlling camalexin biosynthesis (Birkenbihl et al., 2012); (2) the Myb-type transcription factor HIG1/MYB51 (At1g18570) involved in the positive regulation of indole glucosinolate biosynthesis by activating several target genes (Gigolashvili et al., 2007); (3) the O-methyltransferases IGMT1 (At1g21100) and IGMT2 (At1g21120), which catalyze the transfer of a methyl group to the hydroxy indole glucosinolate hydroxyindol-3-ylmethylglucosinolate (4 and 1OH-I3M, respectively) to form methoxyindol-3-ylmethylglucosinolate (4 and 1MO-I3M, respectively) (Pfalz et al., 2011); and (4) cytochrome P450 monooxygenase CYP81F2 (At5g57220), that is essential for the pathogen-induced accumulation of 4-methoxyindol-3ylmethylglucosinolate (4MI3G) (Bednarek et al., 2009). Our data thus show the early induction of indole glucosinolate and camalexin biosynthesis genes in GO plants after metabolic formation of H 2 O 2 through the activation of genes encoding enzymes involved in intermediate metabolite conversions and of TFs that act on several target genes of these biosynthetic pathways.

TRANSCRIPTION FACTOR PROFILING
To understand the potential effects of overexpression of GO in chloroplasts on the nuclear transcriptional program, we next broadened our analysis by testing the expression of 1880 TFs using a highly sensitive quantitative real-time PCR (qRT-PCR) platform (Czechowski et al., 2004;Balazadeh et al., 2008). Considering the data obtained from the profiling of the ROSresponsive genes, we analyzed the expression at 0.5 h after induction of H 2 O 2 production to capture the early-responsive TFs. Expression profiling was performed in two biological replicates and log-fold change (log2 FC) ratios of expression changes were calculated for GO and WT plants by comparing gene expression levels before and after the transfer of plants grown at high CO 2 to ambient CO 2 . TFs most strongly responding to H 2 O 2 were identified by comparing their expression FC in GO and WT plants. A TF was considered differentially expressed when the FC in GO was more than 3-fold (log 2 ≥ 1.56) and less than 2-fold in the wild type (log 2 ≥ 1.0) ( Table 4). Analysis of transcript profiles revealed that the expression of 1449 genes, representing 77% of all TF genes tested, could be detected ( Table A2 in Appendix). The remaining 23% (431 of the 1880 TFs) did not yield detectable PCR amplicons, indicating no or very weak expression in the tested material.
At 0.5 h after shifting plants to ambient CO 2 concentration, 78 of the 1449 genes were induced by more than 3-fold in GO plants, whereas in WT plants the expression changes of the same genes were less than 2-fold ( Table 4). Using published data the involvement/participation of the TFs in specific biological processes (Table 4) Table 4 and Figure 2). JAZ proteins bind directly to the key transcription factor MYC2 and thereby prevent JA-dependent gene transcription (Chini et al., 2007;Pauwels et al., 2010). At the same time JAZ genes are rapidly induced by JA and some are MYC2-regulated. This feedback loop regulation would provide a rapid on and off switch of the pathway involving JA. Transcriptional activation of JAZ genes was found to occur in response to several biotic www.frontiersin.org November 2012 | Volume 3 | Article 234 | 7   Figure 2). The early activation of camalexin and indole glucosinolate biosynthesis was also observed in the analysis performed with the ROS-responsive gene platform (Table 1). Finally, FG5 includes TFs with currently unknown functions (Table 4 and Figure 2).

FIGURE 2 | Pie chart representation of the five functional groups (FG) of early H 2 O 2 -responsive TFs in GO plants.
FG5 includes genes for which a distinct biological function has not been reported yet. The analysis of the transcript profiles at 0.5 h after induction of H 2 O 2 production in GO plants ( Table A2 in Appendix) also revealed a group of 13 genes that are down-regulated in GO relative to WT plants ( Table 5). The function of five of these genes is currently unknown, but interestingly, the remaining eight genes positively control developmental processes. The down-regulation of expression of these TFs in GO plants together with the upregulation of expression of TFs negatively affecting development (see FG2, Table 3) would act in concert to arrest growth and especially to delay the transition from the vegetative to the reproductive phase. Consistently, our previous results showed that GO plants growing under photorespiratory conditions are smaller than WT plants, presenting a reduced rosette diameter and a delay in flowering time (Fahnenstich et al., 2008).

CONCLUDING REMARKS
The metabolic induction of H 2 O 2 formation in chloroplasts of GO plants under photorespiratory conditions triggered a faster and more prominent transcriptional response of ROS-responsive genes in these plants than in wild type. The changes of the transcriptional activities observed in GO plants early after induction of H 2 O 2 production in chloroplasts suggest the establishment of responses that resemble those occurring at early times after wounding or herbivore attack, where H 2 O 2 is also produced (Orozco-Cardenas and Ryan, 1999). These responses include (1) the retardation of development, which in part would be linked to JA signaling, and (2) the production of the phytoanticipins indole glucosinolates and camalexin. As in the case of herbivore attack, the retardation of development such as reductions in growth and reproduction observed in GO plants could be regarded as a strategy to allow more resource allocation to support defense and tolerance responses (Zavala and Baldwin, 2006). The data also suggest that signals for the early induction of indole glucosinolate and camalexin biosynthesis genes in GO plants through H 2 O 2 may originate in chloroplasts as these genes showed no modified expression in catalase loss-of-function mutants (Inzé et al., 2012).