Identification of Regulatory Pathways Controlling ... - Plant Physiology

Identification of Regulatory Pathways Controlling Gene Expression of Stress-Responsive Mitochondrial Proteins in Arabidopsis1[W][OA] Lois H.M. Ho2, Estelle Giraud2, Vindya Uggalla, Ryan Lister, Rachel Clifton, Angela Glen, Dave Thirkettle-Watts, Olivier Van Aken, and James Whelan* Australian Research Council Centre of Excellence in Plant Energy Biology, University of Western Australia, Crawley 6009, Western Australia, Australia

In this study we analyzed transcript abundance and promoters of genes encoding mitochondrial proteins to identify signaling pathways that regulate stress-induced gene expression. We used Arabidopsis (Arabidopsis thaliana) alternative oxidase AOX1a, external NADP H-dehydrogenase NDB2, and two additional highly stress-responsive genes, At2g21640 and BCS1. As a starting point, the promoter region of AOX1a was analyzed and functional analysis identified 10 cis-acting regulatory elements (CAREs), which played a role in response to treatment with H2O2, rotenone, or both. Six of these elements were also functional in the NDB2 promoter. The promoter region of At2g21640, previously defined as a hallmark of oxidative stress, shared two functional CAREs with AOX1a and was responsive to treatment with H2O2 but not rotenone. Microarray analysis further supported that signaling pathways induced by H2O2 and rotenone are not identical. The promoter of BCS1 was not responsive to H2O2 or rotenone, but highly responsive to salicylic acid (SA), whereas the promoters of AOX1a and NDB2 were unresponsive to SA. Analysis of transcript abundance of these genes in a variety of defense signaling mutants confirmed that BCS1 expression is regulated in a different manner compared to AOX1a, NDB2, and At2g21640. These mutants also revealed a pathway associated with programmed cell death that regulated AOX1a in a manner distinct from the other genes. Thus, at least three distinctive pathways regulate mitochondrial stress response at a transcriptional level, an SA-dependent pathway represented by BCS1, a second pathway that represents a convergence point for signals generated by H2O2 and rotenone on multiple CAREs, some of which are shared between responsive genes, and a third pathway that acts via EDS1 and PAD4 regulating only AOX1a. Furthermore, posttranscriptional regulation accounts for changes in transcript abundance by SA treatment for some genes.

The alternative oxidase (AOX) of plant mitochondria is widely used as a model to study the regulation of genes encoding mitochondrial proteins in response to stress or mitochondrial dysfunction (Vanlerberghe and McIntosh, 1997; Finnegan et al., 2004; Clifton et al., 2006; Rhoads et al., 2006). In addition to AOX several other genes are also induced by treatments that perturb mitochondrial function. A protein encoded at the locus At2g21640 has been shown to be induced by oxidative stress (Sweetlove et al., 2002) and subsequently defined as one of five hallmarks of oxidative stress (Gadjev et al., 2006). An analysis of a large number of microarray data indicated that a gene encoding a protein called BCS1, an ortholog of a 1

This work was supported by the Australian Research Council Centre of Excellence in Plant Energy Biology (CEO561495). 2 These authors contributed equally to the article. * Corresponding author; e-mail [email protected]. The author responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (www.plantphysiol.org) is: James Whelan ([email protected]). [W] The online version of this article contains Web-only data. [OA] Open Access articles can be viewed online without a subscription. www.plantphysiol.org/cgi/doi/10.1104/pp.108.121384 1858

protein involved in assembly of cytochrome bc1 in yeast (Saccharomyces cerevisiae; Nobrega et al., 1992), changes transcript abundance in response to at least as many treatments as AOX1a in Arabidopsis (Arabidopsis thaliana; Clifton et al., 2006). In addition, a gene encoding an external NADPH dehydrogenase, NDB2, follows a similar pattern of transcript change to AOX1a in Arabidopsis (Clifton et al., 2005). However, it is unknown if or how the signaling pathways that regulate the transcript abundance for these genes interact and how many signaling pathways exist. An initial microarray analysis of AOX induction suggested significant overlap in the global pattern of transcript abundance changes between chemical inhibition of mitochondrial function, and abiotic and biotic stresses. This suggests an overlap in the pathways responsible for induction of AOX and stress-responsive genes (Yu et al., 2001). A comprehensive analysis of over 250 arrays supported this initial observation, indicating that induction of AOX transcript occurs under a variety of stress treatments, with a greater response to abiotic stresses (Clifton et al., 2006). In addition, coexpression of NDB2, was observed in almost all cases of AOX1a induction, suggesting coregulation (Allocco et al., 2004; Clifton et al., 2005, 2006; Elhafez et al., 2006). Biochemical studies suggest there are at least two pathways that signal the induction of

Plant Physiology, August 2008, Vol. 147, pp. 1858–1873, www.plantphysiol.org Ó 2008 American Society of Plant Biologists

Regulation of Stress-Induced Mitochondrial Gene Expression

AOX. These pathways include a reactive oxygen species (ROS) independent pathway, based on the ability of added citrate to induce AOX transcript and protein without any observed increase in ROS in tobacco (Nicotiana tabacum) suspension cells, and an ROSdependent pathway (Vanlerberghe et al., 1998; Gray et al., 2004). In soybean (Glycine max), a similar conclusion was reached based on the ability of added citrate to induce the expression of AOX1, which was inhibited by the protein kinase inhibitor staurosporine, whereas induction of AOX1 by antimycin A was not sensitive to the addition of staurosporine (Djajanegara et al., 2002). Although much research has been carried out into the various parameters that can induce AOX (activity, protein, or transcript), there is a scarcity of information of how these varied treatments lead to the induction of AOX at a molecular level (Vanlerberghe and McIntosh, 1997; Finnegan et al., 2004; Rhoads and Subbaiah, 2007). The regulatory regions of genes, generally referred to as promoters, are the end point of any signaling pathway for induction at a transcriptional level. An analysis of the AOX1a promoter from Arabidopsis to determine which regions were responsible for stress response concluded that a 93-bp region, termed the mitochondrial retrograde response (MRR) region, responded to two different treatments, antimycin A and monofluoroacetate (Dojcinovic et al., 2005). This is of interest because antimycin A inhibits complex III and monoflouroacetate inhibits the tricarboxylic acid cycle, thus these two treatments might mimic the ROS and citrate pathways that have been proposed above. It is not known if these treatments act on the same or distinct cis-acting regulatory elements (CAREs) in this 93-bp region and if similar CAREs are present in other genes encoding mitochondrial proteins that are induced under similar circumstances. In addition to the fact that promoters are the end point of signaling pathways, the response of the promoter reveals the transcriptional response to the treatment, whereas analysis of transcript abundance, protein, or activity levels also incorporate a variety of posttranscriptional regulatory mechanisms. To gain a better understanding of the induction of genes encoding mitochondrial proteins at a transcriptional level, we analyzed the promoter region of Arabidopsis AOX1a to provide a basis to compare the promoters of other stress-responsive genes. AOX1a is the highest expressed AOX gene under untreated conditions in Arabidopsis (Thirkettle-Watts et al., 2003) and the most stress-responsive AOX gene with respect to the number of stresses and the magnitude of response (Clifton et al., 2005). The occurrence and function of CAREs defined in the AOX1a promoter was then analyzed in three other genes encoding mitochondrial proteins that have been observed to be induced with AOX1a after stress treatment (Clifton et al., 2006), NDB2, a gene at locus At2g21640, known as UPOX (up-regulated by oxidative stress), and BCS1. Also the Plant Physiol. Vol. 147, 2008

response in transcript abundance of these genes in various defense signaling mutants was determined to investigate which signaling pathways are involved. A cellular context for the induction of AOX1a was obtained by analyzing the occurrence of these elements in the promoters of all genes in Arabidopsis.

RESULTS The Transcript Abundance of Nuclear Genes Encoding Mitochondrial Proteins Change in Response to a Variety of Treatments

The alterations in transcript abundance for AOX1a and NDB2 upon treatment with H2O2, rotenone, salicylic acid (SA), and citrate have been previously documented (Clifton et al., 2005, 2006). Microarray analyses suggest that transcript abundance for genes encoding other mitochondrial proteins are also altered in response to these treatments (Clifton et al., 2006). The changes in transcript abundance upon treatment for UPOX, BCS1, PR1a (pathogen-related protein 1a, a positive control for response to SA [Shah, 2003]) and ubiquitin (UBC, a gene that is unresponsive to a wide variety of treatments [Czechowski et al., 2005]) were analyzed using quantitative reverse transcription PCR (QRT-PCR; Fig. 1). The data for AOX1a and NDB2 are replotted from Clifton et al. (2005) for comparison and the same samples were used to determine the transcript abundance for the genes outlined above. The QRT-PCR data demonstrate that transcript abundance for UPOX and BCS1 increase with these treatments, however, there are differences in the magnitude and kinetics of the response. In the case of AOX1a, the response to all treatments (except citrate; see below) peaked at 3 h, with transcript abundance increasing approximately 6-fold for treatment with H2O2, rotenone, and SA. The transcript abundance for NDB2 was most similar to that of AOX1a, except that it peaked at 12 h, with increases of 4- to 8-fold evident with H2O2, rotenone, and SA. For both AOX1a and NDB2, transcript abundance decreased at 24 h compared to the peak observed at 3 or 12 h, respectively. In the case of UPOX, transcript abundance for treatment with H2O2 and rotenone was not significantly different at 3 h, but increased afterward with a 20-fold increase observed with rotenone at 24 h. BCS1 differed in that treatment with SA caused, by far, the greatest increase in transcript abundance at all time points, with greater than a 25-fold increase observed at 12 and 24 h. Thus, even though analysis of microarray data at single time points suggest that these four genes are all coexpressed in response to a variety of treatments (Clifton et al., 2006), analysis of the kinetics and magnitude suggest significant differences in the response that may be overlooked in global analysis. This suggests that the induction of transcript abundance for these genes is not identical and that a different regulatory mechanism may be involved. 1859

Ho et al.

The transcript abundance of all the genes tested was largely unresponsive to citrate (Fig. 1), with only NDB2 showing barely a 2-fold increase in transcript abundance after 24 h. Identification of Sequence Elements That Play a Role in the Stress Response in the Promoter of AOX1a

The stress-responsive AOX1a promoter (Dojcinovic et al., 2005) was used as a starting point to identify CAREs that play a role in regulating expression under various treatments. Five prediction programs were used to identify putative regulatory cis-sequence elements in the AOX1a promoter: PlantCare (Rombauts et al., 1999), PLACE signal scan (Higo et al., 1999), AthaMap (Steffens et al., 2004), Athena (O’Connor et al., 2005), and AGRIS (Palaniswamy et al., 2006). These approaches predicted numerous putative sequence elements in the AOX1a promoter region (Clifton, 2006), therefore three methods were used to distill this number with the aim of increasing the likelihood of targeting elements involved in the H2O2 and rotenone response. Firstly, lists of genes that were coexpressed with AOX1a were constructed, allowing elements to be predicted from these coexpressed groups (Clifton, 2006). Secondly, we used a biclustering approach (Holt et al., 2006) to form a list of coexpressed genes from which motifs were predicted that could describe the patterns of transcript abundance observed. Finally, a phylogenetic approach compared the promoter regions of AOX1a to the AOX promoter regions of soybean (Thirkettle-Watts, 2004). An outline of the predicted elements for AOX1a is shown in Supplemental Figure S2. Different methods were used to predict the elements and with no significant experimental data relating to CAREs in the AOX1a promoter, it was not possible to rank or prioritize the elements. Element E, predicted by comparison with the soybean AOX promoter sequences and by the coexpression environments was identified as active in a previous study (Dojcinovic et al., 2005). We tested all the elements identified by phylogenetic comparison: elements A1, A2, B1, B2, E, J, K, and L; biclustering elements B1, B2, and H; and six elements predicted using hierarchical clustering coexpression (elements C, D, F, G, I1, and I2; Fig. 2). This group, therefore, consists of elements predicted by more than one approach: B1 and B2, predicted using biclustering and by comparison with soybean; elements predicted in several coexpression environments (I1 and C); and elements predicted by single approaches (F and G), identified only from coexpression environments; H, identified from biclustering; and elements J and K, identified by comparison to soybean.

Figure 1. Transcript abundance of AOX1a, NDB2, BCS1, UPOX, PR1, and UBC. QRT-PCR analysis of transcript abundance for various genes over 24 h in response to addition of citrate, H2O2, rotenone, and SA to 1860

Arabidopsis cell cultures. The amount of transcript prior to addition of compounds was set to 1 and changes expressed in a relative manner. An asterisk indicates a significant difference (P # 0.05) Data for AOX1a and NDB2 are redrawn from Clifton et al. (2005). Plant Physiol. Vol. 147, 2008

Regulation of Stress-Induced Mitochondrial Gene Expression

Figure 2. Summary of the predicted CAREs in the AOX1a and NDB2 promoter regions. Putative CAREs were identified in the AOX1a promoter region using a variety of prediction methods. Numbering is from the TSS, bases 5# are indicated with a ‘‘2’’ and bases after the TSS with a ‘‘1’’. The predicted translational start codon ATG is indicated. The 12 elements tested were designated A to L and if the element occurred more than once it was designated by both a letter and a number. A region termed the MRR defined in a previous study is indicated that contains element E analyzed in this study (Dojcinovic et al., 2005).

Altogether, 15 different elements were tested, representing 12 distinct sequences as three occurred twice in the promoter region analyzed. Ten of the 15 sequence elements were functional in the AOX1a promoter, including two found in multiple locations (A and I), resulting in the identification of 12 functional regulatory motifs (Fig. 3; Supplemental Fig. S3; Table I). A positive role in response to H2O2 treatment was defined for 10 elements (B2, C, D, E, F, G, H, I1, I2, and J). Deletion of these elements abolished or suppressed the increase in GUS activity driven by the AOX1a promoter in response to H2O2 treatment. The greatest fold change in this study was observed with the deletion of B2 resulting in a 5-fold increase of GUS activity. Deletion of this element also largely abolished induction of GUS activity by H2O2 with an increase in levels only 20% of the increase observed in the control. This defines element B2 as a strong repressor under normal conditions, which can be partially derepressed by H2O2 and to a lesser extent rotenone. Elements C, D, E, F, G, H, I1, I2, and J displayed a similar pattern in that the 50% increase in GUS activity observed with H2O2 treatment was largely abolished upon deletion. However, deleting these elements had little or no effect on basal levels of GUS activity, determined by comparing untreated samples. Thus, these elements were all defined as positive regulators of the response to H2O2. Elements A1, D, F, I1, I2, and J were defined as playing a role in the rotenone response in the AOX1a promoter (Fig. 3; Table I; Supplemental Fig. S3). Elements A1 and D act as repressors of the rotenone response as deletion of these elements was accompanied by a greater response to rotenone. Elements F, I1, and J were classified as positive response elements to rotenone on the basis that deletion resulted in a loss of induction of GUS activity in response to rotenone treatment. The elements A2, E, I1, and I2 acted as an activator and B2 as a repressor under untreated conditions. Deletion of the A2, E, I1, and I2 resulted in the loss of GUS activity whereas deletion of the B2 resulted in an increase in GUS activity. Plant Physiol. Vol. 147, 2008

Functional Elements in AOX1a Also Function in Other Genes That Are Coinduced with AOX1a under Various Treatments

The promoter regions of NDB2, UPOX, and BCS1 were searched for the presence of the CAREs that were functional in AOX1a, and promoter fragments containing these elements were cloned and tested for function and response to treatments. Six of the functional elements found in AOX1a were also found to be present and functional in the 1-kb upstream region of the transcriptional start site (TSS) of NDB2; elements B, C, F, G, H, and I (Fig. 2; Table I). Upon deletion of element B, untreated activity was reduced and the increase in GUS activity in response to H2O2 and rotenone was abolished, indicating this element was a positive element under normal conditions and played a role in the H2O2 and rotenone response. Deleting elements C, F, and H showed no effect on untreated activity, but abolished the H2O2 and rotenone stimulation. The deletion of element G increased basal activity by approximately 30%, defining it as a repressor, but had no effect on induction by H2O2 and rotenone. Deleting element I reduced basal activity to 50% or less and the induction observed with H2O2 and rotenone was abolished. Analysis of the promoter region of UPOX, containing overlapping B and I elements, revealed it is responsive to H2O2, increasing reporter activity by 2.5fold, however, it was not responsive to rotenone and only slightly responsive to SA (Fig. 4A). Removal of the overlapping B 1 I element reduced reporter activity to 50% and induction by H2O2, was abolished. The BCS1 promoter only contained one common element, H, and the approximately 700-bp region upstream of the TSS of BCS1 containing this element was cloned and tested. A 2.5-fold response to SA was observed, along with a small but significant response to rotenone but no response to H2O2 (Fig. 4B). Deletion of the H element revealed that it plays no role in driving expression. The AOX1a and NDB2 promoters showed no response to SA (Fig. 4C). To test for a repressor to SA 1861

Ho et al.

induction in the AOX1a and NDB2 promoters, extensive analysis of the effects of SA on the AOX1a and NDB2 promoters was performed. Each predicted CARE was deleted and various lengths of the upstream region were tested, but no response to this treatment was observed (data not shown). Thus, it was concluded using the promoters’ regions that were responsive to H2O2 and rotenone that no SA response regions were present. Finally, we tested the effect of citrate on the AOX1a promoter and also detected no increase in promoter activity (data not shown). H2O2 and Rotenone Have Overlapping But Distinct Effects on Transcript Abundance

From the above results, two points suggest that the response to rotenone and H2O2 is overlapping but not identical. Firstly, the promoter region of UPOX does not respond to rotenone, yet displays a strong response to H2O2, secondly there are CAREs in the AOX1a promoter that repress the response to rotenone but not to H2O2. To investigate this further we reexamined microarray data from Arabidopsis cell cultures treated with rotenone and H2O2 (Clifton et al., 2005) to determine the similarity of the response to each treatment. We also carried out microarray analysis of plant leaves treated with H2O2 and rotenone for 3 h to identify overlapping and/or distinct responses. This was carried out as there are some differences in alterations of transcript abundances between suspension cells and leaf tissue with reference to the abundance of transcripts from genes that encode proteins involved in ROS metabolism (Umbach et al., 2005). Studies with tobacco suspension cultures report alterations in some antioxidant components associated with mitochondria (Maxwell et al., 1999; Amirsadeghi et al., 2006), whereas with Arabidopsis leaf tissue no changes were observed (Umbach et al., 2005). Additionally tobacco suspension cultures grown under low inorganic phosphate have increased AOX, but leaves from tobacco plants grown under low inorganic

Figure 3. Functional analysis and stress response of predicted CAREs in the AOX1a and NDB2 promoters. The regulatory characteristics of each predicted CARE was tested by comparing the GUS activity driven by the AOX1a unmutated promoter to the promoter with the element deleted. The unmutated promoter represents the 1.85-kb and 1.0-kb regions upstream of the TSS in AOX1a and NDB2, respectively. The labeling of the element above the graph indicates the sequence element that was deleted from this fragment, as listed in Table I. The activity of the wild-type promoter region (normalized value set to 100%) and 1862

corresponding deletion was tested with H2O2 (green) and rotenone (blue) treatment. Activities for untreated (light), mock-treated (medium), and treated (dark) samples are shown for H2O2 and rotenone. A red asterisk indicates a significant difference (P # 0.05) between the GUS activities of untreated samples with unmutated versus mutated promoter fragments indicating the presence of a constitutive element (positive or negative). A black asterisk indicates a significant difference (P # 0.05) between the GUS activity of mock-treated versus treated samples indicating a stress response is associated with the promoter fragment being tested. A green asterisk indicates a significant difference (P # 0.05) between the GUS activities of treated samples with unmutated versus mutated promoter fragments indicating that when the element was absent the stress response was compromised. The latter two tests define a significant stress response when the element was present (black asterisk), which is affected when the element was deleted (green asterisk). Data for all six elements are shown for NDB2, the data for elements A1, A2, B1, and B2 are shown for AOX1a, and the data for the other elements tested are shown in Supplemental Figure S3. Plant Physiol. Vol. 147, 2008

Regulation of Stress-Induced Mitochondrial Gene Expression

Table I. Summary of the elements tested in the promoters analyzed in this study The CAREs in the AOX1a promoter are designated by a letter. If a sequence element occurs more than once this is indicated by a number. The position of the sequence element is given upstream of the TSS. A ‘‘1’’ or ‘‘2’’ a positive of negative element in untreated conditions, deleting a positive elements results in a loss of activity, whereas deleting a negative element results in an increase of activity. An ‘‘H,’’ ‘‘R,’’ or ‘‘SA’’ indicates a response to H2O2, rotenone, or SA, respectively, with a ‘‘1’’ or ‘‘2’’ indicating a positive or negative regulator. The position of the as1 element characterized to be responsive to SA is indicated in the BCS1 promoter, but this element was not tested for activity in this study. Motif Designation

Sequence

AOX1a

A1 A2 B1

TGAAGC TGAAGC CGTGAT

R2 (1642–1637) 1 (1594–1589) 1623–1618

B2 C

CGTGAT ATCCG

2; H1 (1581–1576) H1 (660–656)

D E F

CACACA CGGCTTT TCGTAAA

1;H1; R2 (316–311) H1 (223–217) H1; R1 (183–177)

G H

TCTCT GTCATC

H1 (169–165) H1 (143–138)

I1

ACGTG

1; H1; R1 (97–93)

I2

ACGTG

1; H1; R1 (73–69)

J K L as1

TTCGATCA TCATCAA ATATAAAC TGACG

H1; R1 (62–55) 57–51 33–26

NDB2

UPOX

1; H1; R1 (846–841)

1; H1; SA1 (95–90)

Previously Reported

Rogers et al. (2001)

H1; R1 (601–597)

H1; R1 (384–378) 2 (19–23) H1; R1 (876–871) 1; H1; R1 (52–48)

phosphate do not (Gonza`lez-Meler et al., 2001). This analysis indicated that there were more distinct responses for each treatment compared to common responses in both leaves and cells (Fig. 5; Supplemental Table S1). The overlap between H2O2 and rotenone treatments in leaves was not significantly different from what is expected in a random distribution (P , 0.05), but the overlap in cell cultures is significantly higher than expected (P , 0.0001). Altogether, this suggests that the signaling pathways triggered by both treatments may overlap, but are not identical. The difference in response between leaves and suspension cells may be due to the fact that cells may have a different metabolic equilibrium with reference to plant tissue (Umbach et al., 2005), particularly in terms of oxidative stress (Halliwell, 2003). Analysis of Transcript Abundance of Genes Encoding Mitochondrial Proteins Induced under Stress in Signaling Mutant Backgrounds

An alternative approach to characterize the pathways that regulate the expression of mitochondrial proteins under stress is to use mutants compromised in various defense signaling pathways. Transcript abundance of genes encoding mitochondrial proteins Plant Physiol. Vol. 147, 2008

BCS1

Dojcinovic et al. (2005)

(463–458)

Eulgem (2005); Kilian et al. (2007) Terzaghi and Cashmore (1995); Yamaguchi-Shinozaki and Shinozaki (2005); Kaplan et al. (2006)

(34–30)

Lam et al. (1989); Jupin and Chua (1996)

was assessed in a variety of these lines, representing the major phytohormone signaling pathways. Specifically, pad4 (phytoalexin-deficient mutant), which acts upstream of SA and is essential for SA-dependent defense pathways (Glazebrook and Ausubel, 1994; Jirage et al., 2001; Glazebrook et al., 2003); npr1 (nonexpressor of the PR1 gene) is required for the expression of defense genes and acts downstream of SA (Cao et al., 1994; Shah, 2003); eds4 (enhanced disease susceptibility), which acts downstream of SA (Gupta et al., 2000); etr1, which confers ethylene insensitivity in Arabidopsis (Chang et al., 1993); jar1, which is deficient in jasmonic acid signaling, particularly associated with defense against insects and necrotrophic pathogens (Staswick et al., 1992; Beckers and Spoel, 2006); abi3, which encodes a viviparous-like transcription factor that interacts with other abscisic acid (ABA)-responsive transcription factors (Giraudat et al., 1992; Nambara et al., 2002); and NahG, plants that express bacterial salicylate hydroxylase that breaks down SA (Lawton et al., 1995; Bowling et al., 1997). Three expression patterns were observed for the stress-responsive genes encoding mitochondrial proteins (Fig. 6). BCS1 was different from all the other genes, decreasing in eds1, pad4, and NahG plants, consistent with a role for SA signal induction, but 1863

Ho et al.

notably this does not depend on NPR1, defining its induction as SA dependent but NPR1 independent, distinguishing it from PR1 (Uquillas et al., 2004). It was notable that the reduction in transcript for PR1a and BCS1 was also greater in magnitude than for AOX1a, NDB2, and UPOX. In NahG plants, BCS1 transcripts were essentially absent, whereas PR1 was just detectable—both reduced by greater than 100-fold compared to the approximately 5-fold reduction for AOX1a, NDB2, and UPOX. This difference is reflected in the induction of the promoters by SA, the promoter of PR1a previously described to be SA responsive (Yang et al., 2000; Uquillas et al., 2004), and BCS1 in this study. AOX1a differed to NDB2 and UPOX in that it was up-regulated in abundance in pad4 and eds1. However, the increase in transcript abundance of AOX1a, NDB2, and UPOX was similar in some the signaling mutants npr1, jar1, abi3, and etr1. A Predicted Coregulatory Context for AOX1a

To explore the association between stress response and presence of stress-responsive elements identified in AOX1a, we compiled a list of the genes that contain six or more of the elements defined as functional in this study in their promoters (AGI, 2000), resulting in a list of 1,141 genes. As with any CARE, there is an expected random occurrence of these elements throughout promoter regions genomewide. Assuming there is an equal chance of the occurrence of the four bases in a promoter region (A, T, C, or G) under the search conditions used in this study, there would be an expected occurrence of elements A, B, C, G, and I of at least once in each promoter region genomewide,

Figure 4. Functional analysis of the stress response of promoters of genes coinduced with AOX1a. A, Analysis of the response of the promoter region of At2g21640 (UPOX) to H2O2, rotenone, and SA. The effect of removing the overlapping B and I element was also tested. B, Analysis of the response of the promoter region of At3g50930 (BCS1) to H2O2, rotenone, and SA. The affect of removing the H element was also tested. C, The response of the AOX1a and NDB2 to SA was tested. The labeling and asterisk are as for Figure 3. 1864

Figure 5. Analysis of genes whose transcript abundance changes in response to H2O2 and rotenone in Arabidopsis suspension cell cultures and leaves. The number of significant changes in transcript abundance upon treatment were determined using Arabidopsis ATH1 microarrays. The number of treatment-specific and common changes are indicated. Positive or negative refers to an increase or decrease in transcript abundance. Plant Physiol. Vol. 147, 2008

Regulation of Stress-Induced Mitochondrial Gene Expression

Figure 6. Expression of AOX1a, NDB2, UPOX, BCS1, PR1a, and UBC in defense signaling mutants. The transcript abundance was determined by QRT-PCR in the defense signaling mutants etr1, abi3, pad4, npr1, jar1, and eds1, and NahG and compared to Col-0 where expression was set to 100%.

31,764 promoter sequences (www.arabidopsis.org). The observed occurrence of these elements was significantly less (11,108, 15,674, 16,328, 8,268, and 17,074, respectively), with a P value ,0.01 according to a x2 test, indicating that these elements may be employed to regulate only a subset of genes. Elements H and D are expected to occur in 15,332 promoter regions, with actual occurrences of 9,225 and 13,014, respectively, significantly less than expected, at a P value ,0.01. The other elements E, F, and J do not show significant overor underrepresentation. However, these analyses need to be interpreted with caution as the assumption of equal occurrence for each of the four bases in promoter regions is not necessarily accurate, as promoters are AT rich, and thus, promoters are difficult to model statistically. Changes in transcript abundance of these 1,141 genes, containing six or more elements, in response to treatments designed to induce oxidative stress were investigated in a number of publicly available microarray data sets (Supplemental Table S2). These included treatment of Arabidopsis plants with cyclohexamide, N-octyl-3-nitro-2,4,6-trihydroxybenzamide, ABA, heat, osmotic stress, salt, ozone, UV light, quartz-filtered UV light, or 2,3,5-triiodobenzoic acid, along with Arabidopsis cell culture treated with rotenone, H2O2, or SA. Data files were downloaded (ftp://ftp.arabidopsis. org/home/tair/microarrays/datasets/) and normalized, and the fold changes in transcript abundance for these 1,141 genes were determined. Transcripts were considered to be responsive to oxidative stress if they showed a greater than 1.5-fold change in response to five or more of the treatments investigated. Of these 1,141 genes, 255 are not represented on the Arabidopsis ATH1 gene chip, leaving a set of 886 genes. Considering the criteria for defining a transcript as stress responsive, 601 out of these 886 genes were responsive to oxidative stress (Supplemental Table S2). These include genes encoding proteins previously reported to play roles in biotic and abiotic stress responses, such as Leu-rich repeat transmembrane protein kinase genes (Ausubel, 2005), several protein kinase genes Plant Physiol. Vol. 147, 2008

including members of the MAPKK and MEKK families (Nakagami et al., 2005; Pitzschke and Hirt, 2006), and transcription factors (Chen et al., 2002; McGrath et al., 2005). However, the list of 601 genes lacks other classes of genes associated with stress, such as enzymes involved in oxidative stress, with the exception of one glutathione-S transferase (At2g29440; van Loon et al., 2006; Supplemental Table S2). Of the 285 genes that contained six or more elements that did not show a 1.5-fold change in five treatments, only 121 genes did not respond to any treatment; the others responded to four, three, two, or one treatments (Supplemental Table S2). Analysis of these 121 genes revealed that 106 of them have expression levels too low to be detected across all microarray experiments, whereas the other 15 genes were too low to be detected in 75% of the microarray experiments. In addition, it is possible that some combinations of six of the 12 elements do not act in a combinatorial manner to regulate gene expression in planta. To investigate any percentile distribution differences in subcellular localization between the list of 601 stress-responsive genes and the whole genome, lists of genes encoding proteins targeted to the mitochondria and the chloroplast were generated (Supplemental Table S2). Lists of 1,025 mitochondrial and 1,407 chloroplastic proteins were isolated, corresponding to 3.22% and 4.48% of the whole genome, respectively. Both mitochondrial and chloroplastic proteins are significantly enriched in the list of 601 genes, at a P value ,0.0001 according to a x2 test, making up 10.98% and 8.4% of this list, respectively. To verify the role of the promoters of these genes in the transcript response we tested the ability of the promoters of a candidate set of 10 genes to drive the expression of GUS under H2O2 treatment. The only criterion used to select genes was that they encode proteins involved in a variety of functions. The promoters tested were for genes encoding transcription factors involved in stress responses (Eulgem, 2005), a ring finger E3 ligase that has been reported to play a role in the hypersensitive response (Kawasaki et al., 1865

Ho et al.

2005), stress-induced proteins predicted to be targeted to mitochondria and chloroplasts, and NIMIN-1, a protein that negatively regulates the activity of NPR1 (Weigel et al., 2005). Transcript analysis of the 10 candidate genes revealed six were up-regulated and four were down-regulated in response to both H2O2 and rotenone (Fig. 7), and analysis of the promoter regions of the 10 candidates revealed that they were all significantly responsive to H2O2 treatment in a manner consistent with changes in their transcript abundance (Fig. 7). Thus, it was concluded that the AOX1a promoter region contains stress-responsive regulatory elements that occur in the promoters of several other stress-responsive genes.

DISCUSSION

The expression/induction of AOX is a widely used model to study the MRR (Rhoads and Subbaiah, 2007), as it is induced by a wide variety of perturbations, including application of chemicals that directly affect mitochondrial function (Finnegan et al., 2004), oxidative stress (Rhoads et al., 2006), genetic lesions (Karpova et al., 2002), and the availability of various nutrients (Parsons et al., 1999; Yip and Vanlerberghe, 2001; Escobar et al., 2006). A number of genes are coexpressed with AOX1a under a variety of conditions; however, their mechanism of induction in relation to AOX1a remains largely unknown (Clifton et al., 2006). The aim of this study was to gain insight into the induction of stress-induced genes encoding mitochondrial proteins, to determine if induction overlapped, and to learn how this may relate to a wider cellular stress response. The promoter regions of four genes were analyzed to determine their ability to respond to a number of treatments and the role of the predicted sequence elements in the various promoters played in this response. The functional elements identified represent a molecular fingerprint of sequence elements that play a role in the response. The element identified may only represent a portion of the sequence element that mediates the response, or interact with other undefined sequence elements to mediate the response. Defining the promoter regions is a somewhat reiterative process and it cannot be ruled out that regulatory regions are missed. However, comparative studies carried out within a single promoter with various treatments and analysis of CAREs and between promoters with differential responses to treatments is informative on the nature of the signals that affect those regions. Under normal conditions, the AOX1a promoter is under strong repression, as evidenced by the 5-fold increase in activity upon removal of element B2. Removal of this element resulted in a promoter activity that dwarfed the effect of the other elements. It appears that application of stress relieves the repression of AOX1a and allows the positive regulatory activity of the other elements characterized to be 1866

exerted, consistent with the combinatorial nature of gene regulation. Elements B and I identified in this study contain an ABRE core, but differ in the flanking sequences that may give different binding specificities. Element B2 and E in AOX1a overlap with an ABI4 binding site (Koussevitzky et al., 2007) and I1 overlaps with a G-box binding site, defined as a light responsive cis-element (Menkens et al., 1995; Terzaghi and Cashmore, 1995). This provides a direct mechanism through which AOX1a can be regulated and/or coregulated with photosynthetic function, consistent with previous reports that AOX expression is under diurnal control (Dutilleul et al., 2003). This I1 element overlaps with an ABRE core, and is in close proximity to another ABRE element, I2, forming an ABRE-ABRE pair (Gomez-Porras et al., 2007), that is proposed to mediate ABA-dependent gene regulation. It has been demonstrated that the expression of 2-cys peroxiredoxin, a well-studied model for transcriptional regulation of a gene encoding a chloroplast antioxidant enzyme, correlates with the acceptor availability of PSI and analysis of the promoter activity in ABAbiosynthetic or ABA-insensitive mutants reveals cross talk between redox and ABA signaling downstream of ABI1 and ABI2 (Baier et al., 2004). Subsequently it was shown that this is mediated by a CE3 element in this promoter that binds Rap2.4a (Shaikhali et al., 2008). The presence of the ABRE-ABRE core in the AOX1a promoter (I1 and I2) and ABRE elements in the functional B2 element provide a means to link the expression of chloroplast antioxidant enzymes and AOX1a, providing a molecular link between expression of AOX1a and chloroplast function that has been revealed in many studies (Noguchi and Yoshida, 2008). As the elements in the 2-cys peroxiredoxin promoter differ in sequence to those in the AOX1a promoter, different ABA response factors may be involved in binding AOX1a. Several of the 10 distinct promoter elements (CAREs) found to be active in controlling AOX1a gene expression were also functional in the NDB2 and UPOX promoters. Many but not all of the elements characterized were responsive to both H2O2 and rotenone. The response of the UPOX promoter to H2O2 is consistent with it being defined as a marker to oxidative stress (Gadjev et al., 2006), yet this promoter was largely unresponsive to rotenone. The findings that the promoter of AOX1a has CAREs that suppress rotenone but not the H2O2 response and that all CAREs in AOX1a and UPOX are not equally responsive to rotenone and H2O2 suggests that they illicit at least some different responses. Furthermore, the AOX1a promoter contains CAREs not found in UPOX or NDB2 and two of these CAREs, A1 and D, repress the response to rotenone. As inhibition of complex I requires alternative routes to oxidize NADH but not necessarily a dramatic increase in AOX activity, the regulating elements are consistent with the biological functions of the encoded proteins. The analysis of transcriptome changes in response to rotenone and Plant Physiol. Vol. 147, 2008

Regulation of Stress-Induced Mitochondrial Gene Expression

Figure 7. The response of promoter regions that contain sequence elements defined as functional in AOX1a. Determination of the promoter-driven response to H2O2 of 10 genes, each containing six or more sequence elements defined as functional in AOX1a. Six genes whose transcript abundance is up-regulated after H2O2 treatment and four genes whose transcript abundance is down-regulated after H2O2 were tested. All promoters displayed a significant response in agreement with the trend observed in transcript abundance. A black asterisk indicates a significant difference (P # 0.05) when comparing the GUS activity between mock-treated and treated samples. Activities for untreated (light), mock-treated (medium), and treated (dark) samples are shown for H2O2. A schematic representation of each promoter is shown with the region in base pair or kilobase pair upstream of the TSS shown. The relative position of the occurrence of the different sequence elements in each promoter tested are indicated. The heat map displays the changes in transcript abundances for the genes whose promoters were tested for responsiveness to H2O2. A red color indicates up-regulation and green indicates down-regulation; the transcript changes for AOX1a, NDB2, and UPOX are included as a comparison. The genes are arranged numerically by AGI code. The At1g02450, At3g50870, and At1g72200 transcript levels were called absent after MAS5.0 normalization.

Plant Physiol. Vol. 147, 2008

1867

Ho et al.

H2O2 using microarrays reveals overlapping and distinct effects. Thus, these pathways may overlap at two levels; firstly, rotenone may trigger ROS production (van der Merwe and Dubery, 2006) and secondly, the signaling pathways converge to act on shared CAREs. As transcription factors are encoded by large gene families that can bind similar or identical CAREs, it cannot be concluded that these pathways act on the same transcription factors, only that the transcriptional responses share CAREs. Although the application of the phytohormone SA increases transcript abundance to the same extent as H2O2 and rotenone (Clifton et al., 2005), the AOX1a and NDB2 promoters were unresponsive to this treatment. The small but significant response of the UPOX promoter to SA may be due to the fact that addition of SA to cells results in an inhibition of electron transport via both the alternative and cytochrome pathways (Norman et al., 2004) and is therefore likely to produce a level of oxidative stress. The sensitivity of UPOX is consistent with the promoter being the most responsive to H2O2 and its definition as a marker for oxidative stress (Gadjev et al., 2006). The difference in the changes in transcript abundance in NahG plants between AOX1a, NDB2, and UPOX (which were reduced by approximately 5-fold) compared to BCS1 (which was reduced greater than 100-fold) support a different mode of action of SA on the respective transcripts. The promoter region of BCS1 was responsive to the SA application, indicating a direct effect on transcription initiation, whereas regulation seems to be posttranscriptional for the other genes tested. Similarly, in Sauromatum guttatum, application of SA can increase transcript abundance of AOX severalfold, but analysis of transcription with isolated nuclei indicated no

change in the rate of transcription, also leading to the conclusion that posttranscriptional mechanisms were involved (Rhoads and McIntosh, 1992). Examination of other regions upstream of the AOX1a promoter region used in these studies revealed no evidence of any characterized SA-responsive elements, such as the as1 element that is present in the BCS1 promoter region (Lam et al., 1989; Jupin and Chua, 1996). It is well established in mammalian systems that hormones have transcriptional and posttranscriptional effects on regulating transcript abundance (Ing, 2005) and that zinc finger transcription factors can bind both DNA and RNA in mammalian systems (Clemens et al., 2003; Lu et al., 2003). Finally, it cannot be ruled out that there are other elements that control the SA response for AOX1a, NDB2, and UPOX. A similar response of AOX1a, NDB2, and UPOX in the mutants etr1, abi3, jar1, and npr1 further support the model that they share common components in signaling pathways that regulate their expression, distinct from BCS1. The exact effects of these mutants on altering transcript abundance cannot be determined from these studies due to the complexity of interactions between various phytohormone signaling pathways (Fujita et al., 2006). However, it was evident that the transcript abundance of AOX1a only increased in the pad4 and eds1 mutants, indicating that it is regulated distinctly by these components compared to the other genes examined in this study. Although PAD4 and EDS1 act upstream of SA to regulate SA signaling (Brodersen et al., 2006), the fact that only the AOX1a transcript increases, suggests an SA-independent induction of AOX1a in these mutants. Such a pathway has been described for the regulation of various transcripts in pad4 and eds1 (Bartsch et al., 2006). Notably,

Figure 8. Summary of signaling pathways that regulate stress-induced transcript abundance for genes encoding mitochondrial proteins. Three distinctive pathways regulate mitochondrial stress response at a transcriptional level, an SA-dependent pathway represented by BCS1. A second pathway that represents a convergence point for signals generated by H2O2 and rotenone that act via the same CAREs. H2O2 and rotenone act at a transcriptional level, with some differences evident between genes, namely rotenone also has a repressive effect on specific regions of the AOX1a promoter, whereas rotenone has no direct effect on the promoter of UPOX. A third pathway that acts via EDS1 and PAD4 regulates only AOX1a. The AOX1a promoter is strongly repressed under normal conditions and is derepressed upon treatment. The regulation of transcript abundance upon treatment with SA is regulated at a transcriptional level for BCS1 and posttranscriptionally for AOX1a, NDB2, and UPOX. SA has a small but significant effect on the UPOX promoter and this may be indirect by triggering oxidative stress (see text for details). 1868

Plant Physiol. Vol. 147, 2008

Regulation of Stress-Induced Mitochondrial Gene Expression

these proteins play a central role in mediating programmed cell death conditioned by Toll-like receptors. A role for AOX in modulating programmed cell death is well described in tobacco lines that have altered AOX levels (Ordog et al., 2002; Robson and Vanlerberghe, 2002; Amirsadeghi et al., 2006). Thus, the up-regulation of AOX1a transcript abundance suggests that this is a unique signaling pathway that regulates AOX1a expression compared to NDB2 and UPOX. In plants deficient in the prohibitin AtPHB3, an inner mitochondrial membrane protein involved in mitochondrial respiration and morphology, the expression of AOX1a, UPOX, and BCS1, but not NDB2, is strongly induced, further demonstrating that multiple combinations of signaling pathways are possible (Van Aken et al., 2007).

CONCLUSION

The integrative use of promoter activity information, transcript analysis, and mutants in specific signaling cascades provides a solid base for dissecting the intricate mechanisms underlying gene expression under stress conditions. Based on our results, we propose that the mitochondrial stress response is mediated by at least three distinct pathways at the transcriptional level (Fig. 8): an SA-dependent pathway, a second pathway that converges on a number of CAREs, previously characterized to bind ABA-responsive transcription factors, and a third pathway that acts via EDS1 and PAD4 regulating AOX1a. Furthermore, we propose from the analysis of the promoters and transcripts in this study that SA acts to increase transcript abundance for AOX1a, NDB2, and UPOX via a posttranscriptional mechanism (Fig. 8), whereas BCS1 is under direct transcriptional control. H2O2 and rotenone can act at a transcriptional level, acting on an overlapping set of CAREs. Because rotenone cannot induce the promoter of UPOX, and transcript abundance for UPOX peaks at 24 h after treatment with rotenone, this suggests that different sets of transcription factors are involved in the targeting of CAREs by rotenone and H2O2. Additionally, because the transcript abundance of AOX1a is affected differently compared to NDB2 and UPOX in the pad4 and eds1 signaling mutants, it suggests an additional pathway, possibly involved in regulating programmed cell death. The identification of CAREs involved in the expression of these genes provides links to the expression of chloroplast antioxidant enzymes and to cellular stress responses because identical or similar CAREs have been characterized previously in these pathways.

MATERIALS AND METHODS Cloning of Arabidopsis thaliana Promoter Regions The promoter regions were cloned using standard protocols and subcloned into pLUS. The numbering of Arabidopsis (Arabidopsis thaliana) promoters are

Plant Physiol. Vol. 147, 2008

given from the TSS as determined from the SIGAL database (Yamada et al., 2003). Initially the 1.85-kb region upstream of the TSS of AOX1a was used based on the fact that small regions of this promoter had previously been shown to be responsive (Dojcinovic et al., 2005). The lengths of the other promoter regions used were based on the presence of overlapping CAREs with AOX1a. Thus, a 1-kb region, 442-bp region, and 689 bp region upstream of the TSS was used for NDB2, UPOX, and BCS1, respectively. The constructs were made as translational fusions with GUS, with the first ATG of the gene of interest used as the start codon for GUS. Thus, the different promoter regions also contained different lengths of 5# untranslated regions (UTRs), 98 bp for AOX1a (At3g22370), 131 bp for NDB2 (At4g05020), 49 bp for UPOX (At2g21640), and 291 bp for BCS1 (At3g509300; 49 bp comprises a UTR, an additional 478 bp comprises an intron). In the case of the other analyzed 10 promoters that contained six or more CAREs in common with AOX1a, a similar criterion was used to determine the size of the region upstream of the TSS cloned. At1g02460 contained a 5# UTR of 36 and 995 bp upstream of the TSS, At1g13990 contained a 54-bp 5# UTR and 1.4 kb upstream of the TSS, At1g70420 contained a 104-bp 5# UTR and 1.3 kb upstream of the TSS, At1g72200 contained a 172-bp 5# UTR and 1.2 kb upstream of the TSS, At2g23320 contained a 65-bp 5# UTR and 1.6 kb upstream of the TSS, At3g50870 contained a 122-bp 5# UTR and 1.1 kb upstream of the TSS, At1g75730 contained a 684-bp 5# UTR and 1.8 kb upstream of TSS, At2g01630 contained a 168-bp 5# UTR and 1.1 kb upstream of TSS, At2g17850 contained a 35-bp 5# UTR and 983 bp upstream of TSS, and At4g28240 contained a 108-bp 5# UTR and 1.6 kb upstream of the TSS. A full list of primers is shown in Supplemental Table S3.

Construction of pLUS The pLUS vector was created as follows. The LUC1 gene, flanked by the omega translational enhancer (5#) and E93# terminator (3#; Welsh et al., 2005), was cloned downstream of the 35S CaMV promoter in pGEM3Zf(1) (Promega). 35S CaMV GUS from pCAMBIA1301 (http://www.cambia.org/daisy/cambia/ 585.html) was then cloned into pGEM3Zf(1)-LUC1 downstream of the multiple cloning site, in the opposite orientation to 35S CaMV-LUC(1) (Supplemental Fig. S1).

Plant Material Growth and Treatments Suspension cell culture from Arabidopsis (ecotype Landsberg erecta) leaf tissue was obtained from (Sweetlove et al., 2002). Cells were grown under 16-h light conditions (100 mE m22 s21) and 8-h dark. Cells were subcultured (1:6 [v/v]) at 7-d intervals to maintain growth. Four-day-old cells were treated 1 h prior to collection for transient biolistic transformation, as described in Clifton et al. (2005). A final concentration of 40 mM was used for rotenone, 100 mM for SA, and 10 mM for H2O2 (Clifton et al., 2005). All Arabidopsis mutant lines were obtained from the Arabidopsis Biological Resource Center at The Arabidopsis Information Resource (TAIR). Defense signaling mutants are: pad4, CS3806; eds4, CS3799; npr1, CS3726; abi3, CS24; jar1, CS8072; and etr1, CS6374. NahG seed was kindly donated by Professor Xing Wang Deng from Yale University.

Biolistic Transformation and Assays for Luc and GUS Transformation was performed using the PDS-1000 system using the Hepta adaptor according to the manufacturer’s instructions (Bio-Rad). After transformation, Arabidopsis suspension cell culture was incubated at 22°C for 24 h under long day conditions of 16 h at approximately 100 mE m-2 s21 light conditions and 8 h of dark on paper discs on osmoticum media. Transiently transformed cells were harvested 24 h after bombardment, disrupted by grinding in a mortar and pestle under liquid nitrogen and cellular contents extracted with the lysis buffer and protocol supplied with the luciferase assay system kit (Roche). Luciferase activity assays were carried out according to the manufacturer’s instructions and activity was measured at 2-s intervals over 20 s, using the Polarstar Optima (BMG Labtech). GUS activity was determined using the fluorimetric GUS assay (Jefferson et al., 1987). Assay samples were taken over 1 h at 3-min intervals. Fluorescence was measured using the Polarstar Optima (BMG Labtech) with excitation at 355 nm and emission at 460 nm. The normalized GUS activity was determined by dividing the GUS fluorescence value by the luciferase activity value for each sample, thus eliminating variation due to transformation efficiency. A minimum of

1869

Ho et al.

nine replicate bombardments were carried out per construct. GUS activities measured for each construct were determined and expressed as a percentage and normalized to control conditions. In the case of the AOX1a promoter, the region 1.85 kb upstream of the TSS was taken as previously 1.3 kb upstream of the TSS had been shown to be responsive to some treatments (Dojcinovic et al., 2005). To test the function of various elements the 5- to 8-bp elements were deleted. In the case of the other promoters, regions that contained the elements defined as functional in AOX1a were cloned and elements deleted. For comparison of the GUS activities of the motif deletions with those of the unmutated promoter, a two-sample t test assuming unequal variances was also performed. Significance was defined as P # 0.05. The following comparisons were carried out to determine the activity of each element tested: 1. A comparison of the normalized GUS activity between the wild-type promoter and the mutated promoter; this determined if the element had any regulatory function in the absence of any stress treatment. Significance for this is indicated with a red asterisk. 2. A comparison between the mock-treated and treated GUS values. This determines if the promoter fragment was stress responsive, and if deleting the element resulted in a loss of a significant effect. Significance for this is indicated with a black asterisk. 3. A comparison between the GUS activities of treated samples with unmutated versus mutated promoter fragments indicating that when the element was absent the stress response was compromised. A green asterisk indicates a significant difference (P # 0.05).

primers used for the genes AOX1a (At3g22370) and NDB2 (At4g05020) have been previously described (Clifton et al., 2005). QRT-PCR primers used for the genes UPOX (At2g21640), BCS1 (At3g50930), PR1 (At2g14610), and UBC (At5g25760) have been described previously (Giraud et al., 2008).

Defining a Putative Coregulatory Environment for AOX1a For each functional sequence element identified in the AOX1a promoter in this study, a gene list was generated containing all Arabidopsis genes with these sequences in their upstream regions using the Patmatch function on the TAIR Web site (http://www.arabidopsis.org/cgi-bin/patmatch/nph-patmatch. pl). Each element nucleotide sequence was searched against Locus Upstream Sequences 1,000 bp (DNA) with the exception of elements A and B, which were searched against Locus Upstream Sequences 3,000 bp (DNA) because these elements occurred more than 1,000 bp upstream of the TSS in AOX1a. Elements were searched in both DNA strands (forward direction and reverse complement), with no mismatch allowed and a maximum of one hit per sequence. Once gene lists were generated, they were compiled such that the overlap, i.e. genes containing more than one element, could be identified. Lists of mitochondrial and chloroplast targeted proteins were generated using the SUBA database (http://www.plantenergy.uwa.edu.au/applications/suba2/ index.php), searching for organellar proteins confirmed via experimental methodology.

Calculating Expected Random Occurrences of Motifs Predicting Functional Elements in the Promoter Region of Arabidopsis AOX1a Three approaches were used to predict functional elements in the promoter region of AOX1a, defined as the sequences upstream of the TSS. These approaches were: (1) defining a common coexpression environment for AOX1a and NDB2 (marked in red; Supplemental Fig. S2), for AOX1a (marked in blue; Supplemental Fig. S2) and using the promoter regions of the genes within the coexpression environments to predict putative sequence elements. Briefly, six coexpression environments were created using three linkage methods (average, centroid, and complete) and two distance metrics (Euclidean and Pearson). Elements were identified in these coexpression environments using the PLACE signal scan (Higo et al., 1999) and PlantCare (Lescot et al., 2002). Additionally, the promoter regions of AOX1a and NDB2 alone were also used to predict elements (marked in green; Supplemental Fig. S2). This analysis was restricted to the 1,000 bp upstream of the TSS to restrict the number of elements predicted. (2) A biclustering approach was used with a list of genes encoding mitochondrial proteins (Holt et al., 2006); and (3) a phylogenetic approach compared the promoter region of AOX1a to the promoter regions for soybean (Glycine max) AOX genes (Thirkettle-Watts et al., 2003). A comparative genomics approach was used to compare the soybean and Arabidopsis promoter regions to identify putative regulatory sequence elements using MotifSampler (Thijs et al., 2002; http://www.esat.kuleuven. ac.be/;thijs/work/motifsampler.html) and Improbizer (http://www.cse. ucsc.edu/;kent/improbizer/improbizer.html). By limiting the comparisons to functional regions previously defined in the soybean promoter regions (Thirkettle-Watts et al., 2003), the background noise created by large amounts of nonfunctional DNA sequence is reduced, increasing the likelihood of identifying functional regulatory sequences (Zheng et al., 2003). These comparisons predicted different but some overlapping elements, which are shown in Figure 2.

QRT-PCR Analysis of Gene Expression QRT-PCR was performed on Arabidopsis leaf tissue from various signaling mutant lines, pad4, npr1, abi3-5, NahG, jar1, etr1, and eds4. Leaf tissue was excised from 3-week-old mutant and ecotype Columbia of Arabidopsis (Col-0) plants. Samples were taken in biological triplicate and snap frozen under liquid nitrogen. Total RNA isolation and cDNA synthesis was carried out as described previously (Lister et al., 2004). Transcript levels were assayed using the LightCycler 480 and the LightCycler 480 SYBR Green I Master (Roche). From each of the independent cDNA preparations, each transcript was analyzed twice. The SE was calculated for every data point. Transcript abundance for Col-0 untreated sample was normalized to one for each gene with all other values presented as relative transcript abundance. QRT-PCR

1870

Assuming each base, A, T, C, or G, has an equal chance of being incorporated into an element within the promoter region, the expected random occurrence of any motif of defined length in a given promoter region, as searched in this study would be: [length of DNA region searched/ (4length of element in bp)] 3 2, as the elements were searched in either direction in double-stranded DNA.

Microarray Analysis Analysis of the changes in transcript abundance in 4-week-old Arabidopsis Col-0 plants treated with 10 mM H2O2 or 40 mM rotenone was performed using Affymetrix GeneChip Arabidopsis ATH1 genome arrays (Affymetrix). Leaves excised from Col-0 plants were submerged in solutions of H2O2, rotenone, or a water control for 1 h before being snap frozen. All samples were collected in triplicate. For each replicate, total RNA was isolated from the leaves of two plants using the RNeasy Plant Mini Protocol (QIAGEN). The quality of the RNA was verified using a Bioanalyzer (Agilent Technologies) and spectrophotometric analysis to determine the A260 to A280 ratio. Preparation of labeled cRNA from 5 mg of total RNA, target hybridization, as well as washing, staining, and scanning of the arrays was carried out exactly as described in the Affymetrix GeneChip Expression Analysis Technical Manual, using an Affymetrix GeneChip Hybridization Oven 640, an Affymetrix Fluidics Station 450, and a GeneChip Scanner 3000 7G at the appropriate steps. Data quality was assessed using GCOS 1.4 before CEL files were exported into AVADIS Prophetic (Version 4.3; Strand Life Sciences) for further analysis. As an additional comparison, raw signal data from microarrays performed previously for Arabidopsis cell culture treated with 10 mM H2O2 or 40 mM rotenone along with an untreated control (Clifton et al., 2005) were analyzed using the program AVADIS Prophetic (Version 4.3; Strand Life Sciences). All CEL files were subjected to GC-robust multiarray average (GC-RMA) normalization, a variation of the RMA normalization algorithm. Fluorescence intensities were log transformed and correlation plots were examined between replicates using the Scatter Plot function, and in all cases r $ 0.9 (data not shown). Fold changes for untreated/mock-treated versus stress-treated cell culture and leaf tissue were calculated using the Differential Expression Analysis function and P values were calculated using an unpaired t test, P values were given as uncorrected and after P-value correction for estimation of false discovery rate in multiple comparisons. A method of false discovery rate correction that is based on the Benjamini and Hotchberg method was used (Benjamini and Hochberg, 1995; Nettleton, 2006). This method utilizes an add-on in the software package R to calculate q values based on P-value distributions (http://www.r-project.org; Storey, 2002). Changes were considered significant with a q value ,0.1. Transcripts significantly up- or down-regulated after treatment with H2O2 or rotenone for cells and plants were plotted on two-way

Plant Physiol. Vol. 147, 2008

Regulation of Stress-Induced Mitochondrial Gene Expression

Venn diagrams such that the overlap in transcripts with altered abundance could be determined. A x-square statistical test using the expected overlap in a random distribution (no. of significantly changed probe sets in set 1 3 no. of significantly changed probe sets in set 2 / no. of probe sets on ATH1 chips) or by random sampling was applied to test the significance of overlap between the two treatments. For the analysis of publicly available microarray data, data was downloaded from the TAIR Web site (ftp://ftp.arabidopsis.org/ home/tair/microarrays/datasets/) for the treatments N-octyl-3-nitro-2,4,6trihydroxybenzamide, cyclohexamide, ABA, 2,3,5-triiodobenzoic acid, osmotic, salt, heat, and UV-B, which are found in the following directories: ME00363, ME00361, ME00333, ME00358, ME00327, ME00328, ME00339, and ME00329, respectively. Access to the microarray data for the treatment of Arabidopsis cell culture with rotenone, SA, and H2O2 is outlined in Clifton et al. (2005) and the treatment of plants with quartz-filtered UV light is outlined in Ulm et al. (2004). CEL files for the treatment of plants with ozone were kindly provided by Alan D Shirras (Lancaster University). Details of the treatments and growth conditions used can be obtained from the original studies. CEL files were subjected to a MAS5.0 normalization to generate present/absent calls and a GC-RMA normalization to summarize fluorescent intensity values for further analysis. Probe sets that were determined as absent across all of the chips were removed and fold-change values were calculated as described above for each control-versus-treatment condition, for a set of 886 genes with six or more motifs. Additionally, a fold change cutoff of 1.5-fold was applied, and changes less dramatic than this were removed.

Supplemental Data The following materials are available in the online version of this article. Supplemental Figure S1. Diagrammatic representation of the biolistic transformation vector pLUS. Supplemental Figure S2. The 1,850 bp upstream of the AOX1a TSS. Supplemental Figure S3. Functional analysis and stress response of predicted CAREs in the AOX1a. Supplemental Table S1. Transcripts that are significantly changing at a q value of ,0.1 in Arabidopsis leaf tissue treated with either H2O2 or rotenone compared to mock-treated plants. Supplemental Table S2. Analysis of changes in transcript abundance for the list of 1,141 genes with six or more elements in the upstream regions in response to treatments designed to induce oxidative stress. Supplemental Table S3. Primers for cloning of promoter regions and deletion of CAREs. Received April 17, 2008; accepted June 11, 2008; published June 20, 2008.

LITERATURE CITED Allocco DJ, Kohane IS, Butte AJ (2004) Quantifying the relationship between co-expression, co-regulation and gene function. BMC Bioinformatics 5: 18 Amirsadeghi S, Robson CA, McDonald AE, Vanlerberghe GC (2006) Changes in plant mitochondrial electron transport alter cellular levels of reactive oxygen species and susceptibility to cell death signaling molecules. Plant Cell Physiol 47: 1509–1519 Arabidopsis Genome Initiative (2000) Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature 408: 796–815 Ausubel FM (2005) Are innate immune signaling pathways in plants and animals conserved? Nat Immunol 6: 973–979 Baier M, Stro¨her E, Dietz KJ (2004) The acceptor availability at photosystem I and ABA control nuclear expression of 2-Cys peroxiredoxin-A in Arabidopsis thaliana. Plant Cell Physiol 45: 997–1006 Bartsch M, Gobbato E, Bednarek P, Debey S, Schultze JL, Bautor J, Parker JE (2006) Salicylic acid-independent ENHANCED DISEASE SUSCEPTIBILITY1 signaling in Arabidopsis immunity and cell death is regulated by the monooxygenase FMO1 and the Nudix hydrolase NUDT7. Plant Cell 18: 1038–1051 Beckers GJ, Spoel SH (2006) Fine-tuning plant defence signalling: salicylate versus jasmonate. Plant Biol (Stuttg) 8: 1–10 Benjamini Y, Hochberg Y (1995) Controlling false discovery rate: a prac-

Plant Physiol. Vol. 147, 2008

tical and powerful approach to multiple testing. J R Stat Soc Ser B 57: 289–300 Bowling SA, Clarke JD, Liu Y, Klessig DF, Dong X (1997) The cpr5 mutant of Arabidopsis expresses both NPR1-dependent and NPR1-independent resistance. Plant Cell 9: 1573–1584 Brodersen P, Petersen M, Bjorn Nielsen H, Zhu S, Newman MA, Shokat KM, Rietz S, Parker J, Mundy J (2006) Arabidopsis MAP kinase 4 regulates salicylic acid- and jasmonic acid/ethylene-dependent responses via EDS1 and PAD4. Plant J 47: 532–546 Cao H, Bowling SA, Gordon AS, Dong X (1994) Characterization of an Arabidopsis mutant that is nonresponsive to inducers of systemic acquired resistance. Plant Cell 6: 1583–1592 Chang C, Kwok SF, Bleecker AB, Meyerowitz EM (1993) Arabidopsis ethylene-response gene ETR1: similarity of product to two-component regulators. Science 262: 539–544 Chen W, Provart NJ, Glazebrook J, Katagiri F, Chang HS, Eulgem T, Mauch F, Luan S, Zou G, Whitham SA, et al (2002) Expression profile matrix of Arabidopsis transcription factor genes suggests their putative functions in response to environmental stresses. Plant Cell 14: 559–574 Clemens KR, Wolf V, McBryant SJ, Zhang P, Liao X, Wright PE, Gottesfeld JM (2003) Molecular basis for specific recognition of both RNA and DNA by a zinc finger protein. Science 260: 530–533 Clifton R (2006) The Alternative Oxidase Gene Family in Arabidopsis: Insights from a Transcriptomic Study. University of Western Australia, Perth, Australia Clifton R, Lister R, Parker KL, Sappl PG, Elhafez D, Millar AH, Day DA, Whelan J (2005) Stress-induced co-expression of alternative respiratory chain components in Arabidopsis thaliana. Plant Mol Biol 58: 193–212 Clifton R, Millar AH, Whelan J (2006) Alternative oxidases in Arabidopsis: a comparative analysis of differential expression in the gene family provides new insights into function of non-phosphorylating bypasses. Biochim Biophys Acta 1757: 730–741 Czechowski T, Stitt M, Altmann T, Udvardi MK, Scheible WR (2005) Genome-wide identification and testing of superior reference genes for transcript normalization in Arabidopsis. Plant Physiol 139: 5–17 Djajanegara I, Finnegan PM, Mathieu C, McCabe T, Whelan J, Day DA (2002) Regulation of alternative oxidase gene expression in soybean. Plant Mol Biol 50: 735–742 Dojcinovic D, Krosting J, Harris AJ, Wagner DJ, Rhoads DM (2005) Identification of a region of the Arabidopsis AtAOX1a promoter necessary for mitochondrial retrograde regulation of expression. Plant Mol Biol 58: 159–175 Dutilleul C, Garmier M, Noctor G, Mathieu C, Chetrit P, Foyer CH, de Paepe R (2003) Leaf mitochondria modulate whole cell redox homeostasis, set antioxidant capacity, and determine stress resistance through altered signaling and diurnal regulation. Plant Cell 15: 1212–1226 Elhafez D, Murcha MW, Clifton R, Soole KL, Day DA, Whelan J (2006) Characterisation of mitochondrial alternative NAD(P)H dehydrogenases in Arabidopsis: intraorganelle location and expression. Plant Cell Physiol 47: 43–54 Escobar MA, Geisler DA, Rasmusson AG (2006) Reorganization of the alternative pathways of the Arabidopsis respiratory chain by nitrogen supply: opposing effects of ammonium and nitrate. Plant J 45: 775–88 Eulgem T (2005) Regulation of the Arabidopsis defense transcriptome. Trends Plant Sci 10: 71–78 Finnegan PF, Soole KL, Umbach AL (2004) Alternative mitochondrial electron transport proteins in higher plants. In DA Day, H Millar, J Whelan, eds, Plant Mitochondria: From Genome to Function, Vol 17. Kluwer Academic Publishers, Dordrecht, The Netherlands, pp 163–230 Fujita M, Fujita Y, Noutoshi Y, Takahashi F, Narusaka Y, YamaguchiShinozaki K, Shinozaki K (2006) Crosstalk between abiotic and biotic stress responses: a current view from the points of convergence in the stress signaling networks. Curr Opin Plant Biol 9: 436–442 Gadjev I, Vanderauwera S, Gechev TS, Laloi C, Minkov IN, Shulaev V, Apel K, Inze D, Mittler R, Van Breusegem F (2006) Transcriptomic footprints disclose specificity of reactive oxygen species signaling in Arabidopsis. Plant Physiol 141: 436–445 Giraud E, Ho LHM, Clifton R, Carroll A, Estavillo G, Tan Y-F, Howell K, Ivanova A, Pogson BJ, Millar AH, et al (2008) The absence of ALTERNATIVE OXIDASE1A in Arabidopsis results in acute sensitivity to combined light and drought stress. Plant Physiol 147: 595–610 Giraudat J, Hauge BM, Valon C, Smalle J, Parcy F, Goodman HM (1992) Isolation of the Arabidopsis ABI3 gene by positional cloning. Plant Cell 4: 1251–1261

1871

Ho et al.

Glazebrook J, Ausubel FM (1994) Isolation of phytoalexin-deficient mutants of Arabidopsis thaliana and characterization of their interactions with bacterial pathogens. Proc Natl Acad Sci USA 91: 8955–8959 Glazebrook J, Chen W, Estes B, Chang HS, Nawrath C, Metraux JP, Zhu T, Katagiri F (2003) Topology of the network integrating salicylate and jasmonate signal transduction derived from global expression phenotyping. Plant J 34: 217–228 Gomez-Porras JL, Riano-Pachon DM, Dreyer I, Mayer JE, MuellerRoeber B (2007) Genome-wide analysis of ABA-responsive elements ABRE and CE3 reveals divergent patterns in Arabidopsis and rice. BMC Genomics 8: 260 Gonza`lez-Meler MA, Giles L, Thomas RB, Siedow JN (2001) Metabolic regulation of leaf respiration and alternative pathway activity in response to phosphate supply. Plant Cell Environ 24: 205–215 Gray GR, Maxwell DP, Villarimo AR, McIntosh L (2004) Mitochondria/ nuclear signaling of alternative oxidase gene expression occurs through distinct pathways involving organic acids and reactive oxygen species. Plant Cell Rep 23: 497–503 Gupta V, Willits MG, Glazebrook J (2000) Arabidopsis thaliana EDS4 contributes to salicylic acid (SA)-dependent expression of defense responses: evidence for inhibition of jasmonic acid signaling by SA. Mol Plant Microbe Interact 13: 503–511 Halliwell B (2003) Oxidative stress in cell culture: an under-appreciated problem? FEBS Lett 540: 3–6 Higo K, Ugawa Y, Iwamoto M, Korenaga T (1999) Plant cis-acting regulatory DNA elements (PLACE) database: 1999. Nucleic Acids Res 27: 297–300 Holt KE, Millar AH, Whelan J (2006) ModuleFinder and CoReg: alternative tools for linking gene expression modules with promoter sequences motifs to uncover gene regulation mechanisms in plants. Plant Methods 2: 8 Ing NH (2005) Steroid hormones regulate gene expression posttranscriptionally by altering the stabilities of messenger RNAs. Biol Reprod 72: 1290–1296 Jefferson RA, Kavanagh TA, Bevan MW (1987) GUS fusions: betaglucuronidase as a sensitive and versatile gene fusion marker in higher plants. EMBO J 6: 3901–3907 Jirage D, Zhou N, Cooper B, Clarke JD, Dong X, Glazebrook J (2001) Constitutive salicylic acid-dependent signaling in cpr1 and cpr6 mutants requires PAD4. Plant J 26: 395–407 Jupin I, Chua NH (1996) Activation of the CaMV as-1 cis-element by salicylic acid: differential DNA-binding of a factor related to TGA1a. EMBO J 15: 5679–5689 Kaplan B, Davydov O, Knight H, Galon Y, Knight MR, Fluhr R, Fromm H (2006) Rapid transcriptome changes induced by cytosolic Ca21 transients reveal ABRE-related sequences as Ca21-responsive cis elements in Arabidopsis. Plant Cell 18: 2733–2748 Karpova OV, Kuzmin EV, Elthon TE, Newton KJ (2002) Differential expression of alternative oxidase genes in maize mitochondrial mutants. Plant Cell 14: 3271–3284 Kawasaki T, Nam J, Boyes DC, Holt BF III, Hubert DA, Wiig A, Dangl JL (2005) A duplicated pair of Arabidopsis RING-finger E3 ligases contribute to the RPM1- and RPS2-mediated hypersensitive response. Plant J 44: 258–270 Kilian J, Whitehead D, Horak J, Wanke D, Weinl S, Batistic O, D’Angelo C, Bornberg-Bauer E, Kudla J, Harter K (2007) The AtGenExpress global stress expression data set: protocols, evaluation and model data analysis of UV-B light, drought and cold stress responses. Plant J 50: 347–363 Koussevitzky S, Nott A, Mockler TC, Hong F, Sachetto-Martins G, Surpin M, Lim J, Mittler R, Chory J (2007) Signals from chloroplasts converge to regulate nuclear gene expression. Science 316: 715–719 Lam E, Benfey PN, Gilmartin PM, Fang RX, Chua NM (1989) Site-specific mutations alter in vitro factor binding and change promoter expression pattern in transgenic plants. Proc Natl Acad Sci USA 86: 7890–7894 Lawton K, Weymann K, Friedrich L, Vernooij B, Uknes S, Ryals J (1995) Systemic acquired resistance in Arabidopsis requires salicylic acid but not ethylene. Mol Plant Microbe Interact 8: 863–870 Lescot M, Dehais P, Thijs G, Marchal K, Moreau Y, Van de Peer Y, Rouze P, Rombauts S (2002) PlantCare, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res 30: 325–327 Lister R, Chew O, Lee MN, Heazlewood JL, Clifton R, Parker KL, Millar

1872

AH, Whelan J (2004) A transcriptomic and proteomic characterization of the Arabidopsis mitochondrial protein import apparatus and its response to mitochondrial dysfunction. Plant Physiol 134: 777–789 Lu D, Searles MA, Klug A (2003) Crystal structure of a zinc-finger-RNA complex reveals two modes of molecular recognition. Nature 426: 96–100 Maxwell DP, Wang Y, McIntosh L (1999) The alternative oxidase lowers mitochondrial reactive oxygen production in plant cells. Proc Natl Acad Sci USA 96: 8271–8276 McGrath KC, Dombrecht B, Manners JM, Schenk PM, Edgar CI, Maclean DJ, Scheible WR, Udvardi MK, Kazan K (2005) Repressor- and activatortype ethylene response factors functioning in jasmonate signaling and disease resistance identified via a genome-wide screen of Arabidopsis transcription factor gene expression. Plant Physiol 139: 949–959 Menkens AE, Schindler U, Cashmore AR (1995) The G-box: a ubiquitous regulatory DNA element in plants bound by the GBF family of bZIP proteins. Trends Biochem Sci 20: 506–510 Nakagami H, Pitzschke A, Hirt H (2005) Emerging MAP kinase pathways in plant stress signalling. Trends Plant Sci 10: 339–346 Nambara E, Suzuki M, Abrams S, McCarty DR, Kamiya Y, McCourt P (2002) A screen for genes that function in abscisic acid signaling in Arabidopsis thaliana. Genetics 161: 1247–1255 Nettleton D (2006) A discussion of statistical methods for design and analysis of microarray experiments for plant scientists. Plant Cell 18: 2112–2121 Nobrega FG, Nobrega MP, Tzagoloff A (1992) BCS1, a novel gene required for the expression of functional Rieske iron-sulfur protein in Saccharomyces cerevisiae. EMBO J 11: 3821–3829 Noguchi K, Yoshida K (2008) Interaction between photosynthesis and respiration in illuminated leaves. Mitochondrion 8: 87–99 Norman C, Howell KA, Millar AH, Whelan JM, Day DA (2004) Salicylic acid is an uncoupler and inhibitor of mitochondrial electron transport. Plant Physiol 134: 492–501 O’Connor TR, Dyreson C, Wyrick JJ (2005) Athena: a resource for rapid visualization and systematic analysis of Arabidopsis promoter sequences. Bioinformatics 21: 4411–4413 Ordog SH, Higgins VJ, Vanlerberghe GC (2002) Mitochondrial alternative oxidase is not a critical component of plant viral resistance but may play a role in the hypersensitive response. Plant Physiol 129: 1858–1865 Palaniswamy SK, James S, Sun H, Lamb RS, Davuluri RV, Grotewold E (2006) AGRIS and AtRegNet. a platform to link cis-regulatory elements and transcription factors into regulatory networks. Plant Physiol 140: 818–829 Parsons HL, Yip JY, Vanlerberghe GC (1999) Increased respiratory restriction during phosphate-limited growth in transgenic tobacco cells lacking alternative oxidase. Plant Physiol 121: 1309–1320 Pitzschke A, Hirt H (2006) Mitogen-activated protein kinases and reactive oxygen species signaling in plants. Plant Physiol 141: 351–356 Rhoads DM, McIntosh L (1992) Salicylic acid regulation of respiration in higher plants: alternative oxidase expression. Plant Cell 4: 1131–1139 Rhoads DM, Subbaiah CC (2007) Mitochondrial retrograde regulation in plants. Mitochondrion 7: 177–194 Rhoads DM, Umbach AL, Subbaiah CC, Siedow JN (2006) Mitochondrial reactive oxygen species. Contribution to oxidative stress and interorganellar signaling. Plant Physiol 141: 357–366 Robson CA, Vanlerberghe GC (2002) Transgenic plant cells lacking mitochondrial alternative oxidase have increased susceptibility to mitochondria-dependent and -independent pathways of programmed cell death. Plant Physiol 129: 1908–1920 Rogers HJ, Bate N, Combe J, Sullivan J, Sweetman J, Swan C, Lonsdale DM, Twell D (2001) Functional analysis of cis-regulatory elements within the promoter of the tobacco late pollen gene g10. Plant Mol Biol 45: 577–585 Rombauts S, Dehais P, Van Montagu M, Rouze P (1999) PlantCARE, a plant cisacting regulatory element database. Nucleic Acids Res 27: 295–296 Shah J (2003) The salicylic acid loop in plant defense. Curr Opin Plant Biol 6: 365–371 Shaikhali J, Heiber I, Seidel T, Stroeher E, Hiltscher H, Birkmann S, Dietz KJ, Baier M (2008) The redox-sensitive transcription factor Rap2.4a controls nuclear expression of 2-Cys peroxiredoxin A and other chloroplast antioxidant enzymes. BMC Plant Biol 8: 48 Staswick PE, Su W, Howell SH (1992) Methyl jasmonate inhibition of root growth and induction of a leaf protein are decreased in an Arabidopsis thaliana mutant. Proc Natl Acad Sci USA 89: 6837–6840

Plant Physiol. Vol. 147, 2008

Regulation of Stress-Induced Mitochondrial Gene Expression

Steffens NO, Galuschka C, Schindler M, Bulow L, Hehl R (2004) AthaMap: an online resource for in silico transcription factor binding sites in the Arabidopsis thaliana genome. Nucleic Acids Res 32: D368–372 Storey JD (2002) A direct approach to false discovery rates. J Roy Statist Soc Ser B Methodological 64: 479–498 Sweetlove LJ, Heazlewood JL, Herald V, Holtzapffel R, Day DA, Leaver CJ, Millar AH (2002) The impact of oxidative stress on Arabidopsis mitochondria. Plant J 32: 891–904 Terzaghi WB, Cashmore AR (1995) Photomorphogenesis. Seeing the light in plant development. Curr Biol 5: 466–468 Thijs G, Marchal K, Lescot M, Rombauts S, De Moor B, Rouze P, Moreau Y (2002) A Gibbs sampling method to detect overrepresented motifs in the upstream regions of coexpressed genes. J Comput Biol 9: 447–464 Thirkettle-Watts D (2004) Functional analysis of the soybean alternative oxidase upstream regions. PhD thesis. University of Western Australia, Perth, Australia Thirkettle-Watts D, McCabe TC, Clifton R, Moore C, Finnegan PM, Day DA, Whelan J (2003) Analysis of the alternative oxidase promoters from soybean. Plant Physiol 133: 1158–1169 Ulm R, Baumann A, Oravecz A, Ma´te´ Z, Ada´m E, Oakeley EJ, Scha¨fer E, Nagy F (2004) Genome-wide analysis of gene expression reveals function of the bZIP transcription factor HY5 in the UV-B response of Arabidopsis. Proc Natl Acad Sci USA 101: 1397–1402 Umbach AL, Fiorani F, Siedow JN (2005) Characterization of transformed Arabidopsis with altered alternative oxidase levels and analysis of effects on reactive oxygen species in tissue. Plant Physiol 139: 1806–1820 Uquillas C, Letelier I, Blanco F, Jordana X, Holuigue L (2004) NPR1independent activation of immediate early salicylic acid-responsive genes in Arabidopsis. Mol Plant Microbe Interact 17: 34–42 Van Aken O, Pecenkova´ T, van de Cotte B, De Rycke R, Eeckhout D, Fromm H, De Jaeger G, Witters E, Beemster GT, Inze´ D, et al (2007) Mitochondrial type-I prohibitins of Arabidopsis thaliana are required for supporting proficient meristem development. Plant J 52: 850–864

Plant Physiol. Vol. 147, 2008

van der Merwe JA, Dubery IA (2006) Benzothiadiazole inhibits mitochondrial NADH:ubiquinone oxidoreductase in tobacco. J Plant Physiol 163: 877–882 van Loon LC, Rep M, Pieterse CM (2006) Significance of inducible defenserelated proteins in infected plants. Annu Rev Phytopathol Vanlerberghe GC, McIntosh L (1997) ALTERNATIVE OXIDASE: from gene to function. Annu Rev Plant Physiol Plant Mol Biol 48: 703–734 Vanlerberghe GC, McIntosh L, Yip JY (1998) Molecular localization of a redox-modulated process regulating plant mitochondrial electron transport. Plant Cell 10: 1551–1560 Weigel RR, Pfitzner UM, Gatz C (2005) Interaction of NIMIN1 with NPR1 modulates PR gene expression in Arabidopsis. Plant Cell 17: 1279–1291 Welsh DK, Imaizumi T, Kay SA (2005) Real-time reporting of circadianregulated gene expression by luciferase imaging in plants and mammalian cells. Methods Enzymol 393: 269–288 Yamada K, Lim J, Dale JM, Chen H, Shinn P, Palm CJ, Southwick AM, Wu HC, Kim C, Nguyen M, et al (2003) Empirical analysis of transcriptional activity in the Arabidopsis genome. Science 302: 842–846 Yamaguchi-Shinozaki K, Shinozaki K (2005) Organization of cis-acting regulatory elements in osmotic- and cold-stress-responsive promoters. Trends Plant Sci 10: 88–94 Yang Y, Li R, Qi M (2000) In vivo analysis of plant promoters and transcription factors by agroinfiltration of tobacco leaves. Plant J 22: 543–551 Yip JY, Vanlerberghe GC (2001) Mitochondrial alternative oxidase acts to dampen the generation of active oxygen species during a period of rapid respiration induced to support a high rate of nutrient uptake. Physiol Plant 112: 327–333 Yu J, Nickels R, McIntosh L (2001) A genome approach to mitochondrialnuclear communication in Arabidopsis. Plant Physiol Biochem 39: 345–353 Zheng J, Wu J, Sun Z (2003) An approach to identify over-represented ciselements in related sequences. Nucleic Acids Res 31: 1995–2005

1873