Skip to main content
Advertisement
  • Loading metrics

Mechanical Stress Induces Biotic and Abiotic Stress Responses via a Novel cis-Element

  • Justin W Walley,

    Affiliation Section of Plant Biology, University of California Davis, Davis, California, United States of America

  • Sean Coughlan,

    ¤ Current address: DuPont, Wilmington, Delaware, United State of America

    Affiliation Agilent Technologies, Wilmington, Delaware, United States of America

  • Matthew E Hudson,

    Affiliation Department of Crop Sciences, University Of Illinois, Urbana, Illinois, United States of America

  • Michael F Covington,

    Affiliation Section of Plant Biology, University of California Davis, Davis, California, United States of America

  • Roy Kaspi,

    Affiliation Section of Plant Biology, University of California Davis, Davis, California, United States of America

  • Gopalan Banu,

    Affiliation Genomic Medicine, Cleveland Clinic, Cleveland, Ohio, United States of America

  • Stacey L Harmer,

    Affiliation Section of Plant Biology, University of California Davis, Davis, California, United States of America

  • Katayoon Dehesh

    To whom correspondence should be addressed. E-mail: kdehesh@ucdavis.edu

    Affiliation Section of Plant Biology, University of California Davis, Davis, California, United States of America

Abstract

Plants are continuously exposed to a myriad of abiotic and biotic stresses. However, the molecular mechanisms by which these stress signals are perceived and transduced are poorly understood. To begin to identify primary stress signal transduction components, we have focused on genes that respond rapidly (within 5 min) to stress signals. Because it has been hypothesized that detection of physical stress is a mechanism common to mounting a response against a broad range of environmental stresses, we have utilized mechanical wounding as the stress stimulus and performed whole genome microarray analysis of Arabidopsis thaliana leaf tissue. This led to the identification of a number of rapid wound responsive (RWR) genes. Comparison of RWR genes with published abiotic and biotic stress microarray datasets demonstrates a large overlap across a wide range of environmental stresses. Interestingly, RWR genes also exhibit a striking level and pattern of circadian regulation, with induced and repressed genes displaying antiphasic rhythms. Using bioinformatic analysis, we identified a novel motif overrepresented in the promoters of RWR genes, herein designated as the Rapid Stress Response Element (RSRE). We demonstrate in transgenic plants that multimerized RSREs are sufficient to confer a rapid response to both biotic and abiotic stresses in vivo, thereby establishing the functional involvement of this motif in primary transcriptional stress responses. Collectively, our data provide evidence for a novel cis-element that is distributed across the promoters of an array of diverse stress-responsive genes, poised to respond immediately and coordinately to stress signals. This structure suggests that plants may have a transcriptional network resembling the general stress signaling pathway in yeast and that the RSRE element may provide the key to this coordinate regulation.

Author Summary

Plants are sessile organisms constantly challenged by a wide spectrum of biotic and abiotic stresses. These stresses cause considerable losses in crop yields worldwide, while the demand for food and energy is on the rise. Understanding the molecular mechanisms driving stress responses is crucial to devising targeted strategies to engineer stress-tolerant plants. To identify primary stress-responsive genes we examined the transcriptional profile of plants after mechanical wounding, which was used as a brief, inductive stimulus. Comparison of the ensemble of rapid wound response transcripts with published transcript profiles revealed a notable overlap with biotic and abiotic stress-responsive genes. Additional quantitative analyses of selected genes over a wounding time-course enabled classification into two groups: transient and stably expressed. Bioinformatic analysis of rapid wound response gene promoter sequences enabled us to identify a novel DNA motif, designated the Rapid Stress Response Element. This motif is sufficient to confer a rapid response to both biotic and abiotic stresses in vivo, thereby confirming the functional involvement of this motif in the primary transcriptional stress response. The genes we identified may represent initial components of the general stress-response network and may be useful in engineering multi-stress tolerant plants.

Introduction

Plants are persistently challenged with numerous biotic and abiotic environmental stresses. To cope with environmental stresses plants have evolved phytohormones such as jasmonic acid, salicylic acid, ethylene, and abscisic acid, which are utilized to regulate plant responses to both abiotic and biotic stresses with considerable signaling crosstalk [1,2]. While these phytohormone pathways have been well studied, knowledge of stress perception and initial signaling events, aside from plant pathogen interactions, are less defined. It is known that application of insect oral secretions containing protein fragments of chloroplastic ATP synthase or application of purified oligouronides (OGAs) derived from the plant cell wall are capable of inducing plant defense responses, although a receptor has not yet been identified [35]. Additionally, a cellulose synthase (CESA3) mutant cev1 shows enhanced resistance to powdery mildew as a result of constitutive increase in jasmonic acid levels in these plants [6]. This has led to the hypothesis that mechanical disruption of the cell wall may result in stress signaling [3,7]. The perception of cold stress has been hypothesized to be mediated through the detection of changes in membrane fluidity and protein conformation [810]. Finally, secondary messengers such as Ca2+, reactive oxygen species (ROS), and phosphatidic acid have been implicated in initial signaling cascades in response to both abiotic and biotic stresses [1117].

One mechanism of response to stress that has been studied extensively in yeast and animals is the general stress response (GSR) (also referred to as the cellular stress response) [18]. The GSR acts in a transient manner in response to a diverse array of stresses. The GSR is initiated in response to strain imposed by environmental forces on macromolecules such as membrane lipids, proteins, and/or DNA. A critical aspect of the GSR, downstream of perception of macromolecular damage, is generation of ROS [19]. Furthermore, key molecular components of the GSR are evolutionarily conserved in all organisms [18].

To better understand plant stress responses, transcript profiling experiments have been successfully employed for many different abiotic and biotic stresses [1,20,21]. One common emerging theme from these experiments is that abiotic and biotic stresses regulate different but overlapping sets of genes [1]. For example, cDNA–amplified fragment length polymorphism analysis of the Avr9- and Cf-9-mediated defense response in tobacco cell culture revealed overlap between race-specific resistance and response to wounding [22,23]. Additionally, partial genome microarray analysis of the Arabidopsis wound response revealed that a number of wound-responsive genes encode proteins known to be involved in pathogen defense [24]. Examination of the AtGenExpress abiotic datasets demonstrates that the initial transcriptional abiotic stress response may comprise a core set of multi-stress-responsive genes. The abiotic stress response then becomes stress specific at later time points [2527]. Finally, recent analysis of the AtGenExpress abiotic and biotic datasets has uncovered ∼200 genes that are expressed in response to a broad range of stresses, which may represent the GSR of Arabidopsis [28].

Recently, a shift in stress tolerance engineering has been proposed that transfers the focus from pathway endpoints to factors governing upstream reactions. Focusing on upstream signaling components may enable the engineering of multi-stress tolerance [20,29]. Identification of cis-regulatory elements for use in synthetic promoters to confer stress tolerance has also recently been proposed [30,31]. Towards this aim, we have utilized mechanical wounding, as it uniquely confers an instantaneous and synchronous stimulus, to identify primary stress-responsive transcripts. Comparison of the 5 min rapid wound response (RWR) genes we identified with published transcript profiles demonstrated a large overlap with previously identified abiotic and biotic stress-responsive genes. Notably, RWR genes also exhibit a striking level and pattern of circadian regulation. Further investigation via real-time quantitative RT-PCR (RT-qPCR) of a wounding time course revealed genes that are expressed rapidly and transiently as well as rapidly and stably. Two rapidly and transiently expressed genes, ETHYLENE RESPONSE FACTOR #018 (ERF#018; AT1G74930) and CCR4-ASSOCIATED FACTOR 1 (CAF1-like; AT3G44260), were confirmed as wound and biotic stress inducible in vivo using stable transgenic lines expressing transcriptional luciferase fusions. Detailed analysis of the RWR promoters identified a novel cis-regulatory element we term the rapid stress response element (RSRE), which is sufficient to confer reporter gene induction in response to abiotic and biotic stress. RWR genes identified in this study may represent initial components of the GSR and be useful in engineering multi-stress tolerance.

Results/Discussion

Transcript Profiling Identifies Rapid Wound Response Genes

To identify primary stress-responsive transcripts we utilized Agilent microarrays to monitor gene expression changes 5 min after mechanical wounding of Arabidopsis rosette leaves. Because of the short duration of our stress treatment we hypothesized that expression changes would be low. In order to accurately detect these changes we utilized three biological replicates of pooled plants per treatment. In addition, two technical replicates, with dye swap of each technical replicate, were performed on each biological replicate. Using this approach, we found that the expression of 162 genes was upregulated and the expression of 44 genes was downregulated at least 2-fold and had a p-value ≤ 0.01 five min after mechanical wounding (Table S1). The expression level of selected RWR genes representing a range of high-to-low-fold change was then validated using RT-qPCR. The expression changes determined by RT-qPCR data are in good agreement with the fold change observed by microarray with a spearman rank order correlation coefficient of 0.927 (p-value = 0.000) (Figure 1A). RWR genes were then classified according to gene ontology (GO) terms in order to provide insight into their biological function [32]. The two largest defined classes of GO terms involve response to stress or abiotic/biotic stimuli (Figure 1B). It is also of interest to note that genes classified for an involvement in signal transduction were observed in upregulated but not downregulated RWR genes. These data indicate that 5 min of mechanical wounding was sufficient for induction of known stress-responsive genes as well as unknown genes that may play a role in multi-stress responses.

thumbnail
Figure 1. Verification and Functional Classification of the RWR Genes

(A) RT-qPCR expression analysis of selected genes normalized to the 60S ribosomal protein L14 (At4g27090) measured in the same samples. The resulting relative expression was then used to calculate fold change upon wounding. Data are means of n = 3.

(B) Functional classification of RWR genes using GO annotations.

https://doi.org/10.1371/journal.pgen.0030172.g001

In Vivo Validation of Rapid Wound Response Genes

We next created stable transgenic lines expressing transcriptional fusions of the ERF#018 and CAF1-like promoters to luciferase to validate in vivo RWR genes and to investigate their temporal expression pattern. For each construct, three independent T2 lines were imaged to control for positional effects of the transgene insertion site. Luciferase activity was then monitored following the wounding of a single leaf per plant to enable the observation of whether the induced activity occurred only locally or also systemically. The wound-induced expression of PERF18:LUC occurs rapidly and peaks ∼1 h 45 min after the wound stimulus (Figure 2A and 2B). Additionally, expression of PERF18:LUC was observed in the petiole and shoot apex. The expression of PCAF1-like:LUC was detected rapidly and peaked ∼1 h 25 min surrounding the wound site (Figure 2C and 2D). These data provide in vivo confirmation that RWR genes do respond rapidly to mechanical wounding.

thumbnail
Figure 2. In Vivo Monitoring of RWR Gene Induction Following Wounding

(A) Image of an individual PERF18:LUC #3 transgenic plant over time.

(B) Luciferase activity of three independent transgenic lines expressing transcriptional PERF18:LUC fusions. Luciferase activity was calculated in wounded leaves (single leaf per plant) and NW leaves (on NW plants). Data are means of n = 12 ± SEM.

(C) Image of an individual PCAF1L:LUC transgenic plant over time.

(D) Luciferase activity of three independent transgenic lines expressing transcriptional PCAF1L:LUC fusions. Luciferase activity was calculated in wounded leaves (single leaf per plant) and NW leaves (on NW plants). Data are means of n = 18 ± SEM.

https://doi.org/10.1371/journal.pgen.0030172.g002

Expression of RWR Genes over Time

To gain further insight into how the RWR genes may be acting, we performed a RT-qPCR time-course on selected genes. We classified genes as rapidly and stably expressed if 60 min post wounding they remained greater than 2-fold induced. In contrast, we classified genes that had decreased in expression to less than half of maximal expression by 60 min post wounding as rapidly and transiently expressed. Among the rapidly and stably expressed transcripts are genes with either a known or predicted role in stress signal transduction events (Figure 3). CML38 is a calmodulin-like gene that is predicted to be a sensor of Ca2+, a known secondary messenger of stress responses [33]. MPKK9, a MAPK signal transduction component, was also identified as rapidly and stably expressed. Additionally, the transcription factor WRKY40, which is known to be involved in pathogen defense, was identified [34]. Finally, BAP1, a negative regulator of defense responses whose binding of phospholipids is enhanced by calcium, was shown to respond rapidly and stably to wounding [35]. The upregulation of RWR genes in a rapid and stable manner may indicate that these genes play a more prolonged role in response to stress.

thumbnail
Figure 3. Selected RWR Genes Displaying a Rapid and Stable Expression Pattern Following Mechanical Wounding

Total RNA was extracted from 3-wk-old mechanically wounded tissue and subject to RT-qPCR analysis. AtMPKK9, WRKY40, AtPP2-A13, BAP1, AT1G32920, and CML38 transcripts were normalized to the 60S ribosomal protein L14 measured in the same samples. Data are means of n = 3 ± SEM.

https://doi.org/10.1371/journal.pgen.0030172.g003

We also uncovered rapidly and transiently expressed RWR genes with a wide range of functions (Figure 4). One such example is the chromatin remodeling ATPase SPLAYED (SYD), which peaks 15 min post wounding. Because of the large changes in gene expression following stress it has been hypothesized that chromatin remodeling may be required to allow for stress-induced transcription to occur [36]. Examples of stress-induced changes in histone acetylation state in plants have been described [37,38]. Rapid and transient upregulation of SYD suggests that ATP-dependent chromatin remodeling may also take place in order to facilitate downstream stress-induced transcriptional changes.

thumbnail
Figure 4. Selected RWR Genes Displaying a Rapid and Transient Expression Following Mechanical Wounding

Total RNA was extracted from 3-wk-old mechanically wounded tissue and subject to RT-qPCR analysis. ERF#018, CAF1-like, AP2C1, MPK3, and TIR-NBS transcripts were normalized to the 60S ribosomal protein L14 while SPLAYED transcripts were normalized to TIP41-like measured in the same samples. Data are means of n = 3 ± SEM.

https://doi.org/10.1371/journal.pgen.0030172.g004

Genes involved in signal transduction via reversible phosphorylation were also upregulated rapidly and transiently following wounding (Figure 4). One kinase identified was MPK3, a MAPK signal transduction component, which has been shown to function in innate immunity and stomatal development [39,40]. AP2C1, a PP2C-type phosphatase with a MAPK interaction motif, was also shown to exhibit a rapid and transient expression pattern resulting from wounding [41]. These results indicate that both phosphorylation and dephosphorylation of MAPK signaling components is involved in transduction of initial stress signaling events.

A third process implicated by genes identified in this study is that of mRNA turnover (Figure 4). Specifically, this process is demonstrated by the expression pattern of CAF1-like. In yeast, CAF1 has been shown to be a component of the major cytoplasmic deadenylase, which functions to remove the poly(A) tail, thereby initiating mRNA turnover [42]. Additionally, the mouse CAF1 ortholog has been demonstrated to function as a 3′-5′-RNase with a preference for poly(A) substrates [43]. In Arabidopsis, RNA processing appears to play a role in response to cold stress [44]. It is also of interest to note the difference in promoter and transcript expression patterns in response to wounding for ERF#018 and CAF1-like. When promoter activity was monitored using transcriptional luciferase fusions, activity peaked ∼1 h 30 min after wounding (Figure 2). In contrast, transcript abundance measured via RT-qPCR peaked 15–30 min after wounding (Figure 4). This discrepancy may be an artifact due to measuring promoter activity via luciferase protein activity, while RT-qPCR assayed transcript levels. An alternative explanation of these results is that mRNA levels of rapidly and transiently expressed genes such as ERF#018 and CAF1-like are controlled post-transcriptionally, possibly via mRNA turnover. In support of this hypothesis, CAF1-like mRNA has a half-life of 38 min and has been classified as a gene with an unstable transcript [45]. Furthermore, the mRNA of a second rapid and transient gene, MPK3, has a half-life of 43 min and is classified as a gene with an unstable transcript [45].

RWR Genes Are Regulated by Abiotic and Biotic Stresses

Examination of the RWR genes reveals a large number of known genes involved in abiotic and biotic stress responses (Table S1). Among the upregulated RWR genes were genes involved in ethylene signaling including ACC synthase 6 (ACS6) as well as 15 of the 122 ethylene response factors (ERFs) in Arabidopsis, which have been shown to be involved in the response to both biotic and abiotic stresses [46]. The transcriptional activators CBF1 (DREB1B), CBF2 (DREB1C), and CBF3 (DREB1A), which confer tolerance to cold and drought, were also among the RWR genes [4750]. Additional RWR genes known to confer tolerance to a range of abiotic stresses include STZ (ZAT10) and ZAT12 [51,52]. RWR genes also include genes with a known function in response to biotic stress. Examples include BAP1, a negative regulator of defense responses [35]. MPK3 and FLS2 which function in response to pathogen-associated molecular patterns in the Arabidopsis innate immune response [39,53]. The transcription factor TGA3 which regulates pathogenesis-related (PR) genes and is required for basal pathogen resistance was also among the RWR genes [54]. Finally, ERD15 regulates not only cold and drought tolerance but also resistance to the bacterial necrotroph Erwinia carotovora subsp. carotovora [55].

The abundance of RWR genes with known abiotic and biotic stress tolerance functions led us to examine the role of wounding as a general stress perception mechanism on a global level. For this analysis, we compared the overlap in gene lists between the RWR genes and published transcript profiles for a number of stress conditions. The statistical significance of the observed overlap in transcript profiles was then analyzed using empirical permutation tests [56]. We first compared RWR genes with published abiotic microarrays and found a strong overlap (unpublished data), which is in agreement with work recently published by Kilian et al. [25]. For example, 49% of upregulated RWR genes have been previously shown to be upregulated upon cold treatment [57,58]. Additionally, four of the nine genes (At1g27730, At5g51190, At5g47230, and At5g04340) found by Kilian et al. [25] to be upregulated by 30 min of cold, drought, UV-B, salt, osmotic stress treatment, and wounding we discovered to be upregulated within 5 min of wounding. We next compared the RWR transcript profile with published transcript profiles of plants challenged with different biotic stresses [5961]. For upregulated datasets there was a statistically significant overrepresentation of RWR genes in the transcript profile of all biotic stresses tested (Figure 5). Furthermore, the overrepresentation of RWR genes occurred only at early time points of P. rapae and OGA stressed plants. These data indicate that perception of mechanical stress may play a central role in the perception and initial response to a wide range of environmental stresses.

thumbnail
Figure 5. Comparison of RWR Genes with the Transcript Profile of Other Environmental Stimuli

An asterisk denotes a statistically significant overrepresentation of RWR genes in the transcript profile of the indicated stress (p < 0.0001). Arabidopsis plants were challenged with M. persicae for 48 h and 72 h; P. rapae early for 5 h; P. rapae late and P. syringae for 12 h and 24 h; B cinerea early and late for 18 h and 48 h, respectively; and OGA early and late for 1 h and 3 h, respectively [5961].

https://doi.org/10.1371/journal.pgen.0030172.g005

RWR Genes Respond to Biological Elicitors

Various biological compounds are known to elicit stress-signaling networks. Due to the overlap between RWR genes and biotic stresses we tested whether the RWR genes ERF#018 and CAF1-like respond to the biological elicitors OGA and insect regurgitant (IR) as well as cabbage looper (Trichoplusia ni) feeding. We first tested whether PERF18:LUC or PCAF1-like:LUC activity was induced by cabbage looper feeding. Indeed, cabbage looper feeding did result in enhanced luciferase activity, which verified that biological stress does induce ERF#018 and CAF1-like (Videos S1 and S2). We therefore proceeded to test induction resulting from OGA and IR treatment. When OGA, IR, or H2O were added to a nonwounded (NW) leaf, no induction of PERF18:LUC or PCAF1-like:LUC activity was observed (Figure 6A–6F). The lack of induction is likely due to the application method we used (single droplet per plant) rather than an actual lack of response to the elicitors. When transcript profiling was performed on liquid cultured 10-d-old plants incubated with 50 μg ml−1 OGAs, CAF1-like was shown to be induced 1 h after addition of OGA to the media [61].

thumbnail
Figure 6. In Vivo Monitoring of RWR Gene Induction upon Addition of Biotic Elicitors

One leaf per plant was either NW or wounded (W) and then treated with oligouronides (OGA), IR, or double-distilled H2O.

(A) Image of an individual PERF18:LUC #3 transgenic plant over time.

(B) Local expression in the treated leaves of PERF18:LUC #3 transgenic plants. Data are means of n = 9 ± SEM.

(C) Systemic expression of PERF18:LUC #3 plants monitored in the shoot apex. Data are means of n = 9 ± SEM.

(D) Image of an individual PCAF1L:LUC #2 transgenic plant over time.

(E) Local expression in the treated leaves of PCAF1L:LUC #2 transgenic plants. Data are means of n = 12 ± SEM.

(F) Systemic expression of PCAF1L:LUC #2 plants monitored in the shoot apex. Data are means of n = 12 ± SEM.

https://doi.org/10.1371/journal.pgen.0030172.g006

In PERF18:LUC-expressing plants, addition of OGA and IR to the wound site resulted in a significantly greater (p < 0.05) induction of luciferase activity than addition of H2O in both local and systemic tissue (Figure 6B and 6C, respectively). In contrast, addition of OGA or IR to the wound site did not result in a significant difference (p > 0.05) in PCAF1-like:LUC activity compared to addition of H2O to the wound site in local tissue (Figure 6E). However, a significantly greater induction, compared to H2O, in PCAF1-like:LUC activity resulted from addition of OGA or IR to the wound site in systemic tissue (Figure 6F).

There are a number of common second messengers downstream of mechanical wounding, cabbage looper feeding, and OGA treatment that may signal for the observed induction of RWR genes. One such secondary messenger is Ca2+, which increases in intracellular concentration rapidly following wounding as well as OGA treatment [3,16,62,63]. ROS are another secondary messenger that have been shown to increase in response to chewing insects, wounding, and OGA treatment [3,12]. Furthermore, while OGAs do not move systemically, ROS do accumulate systemically following wounding. This increase in ROS is likely through OGAs released by systemically induced polygalacturonase [14,16,6466]. Finally, OGA, chewing insects, and wounding may all have a common mechanism of perception resulting in similarly induced secondary messengers. Both chewing insects and wounding have a physical effect on the plasma membrane. The perception of OGA has also been hypothesized to be a result of its physical effect on the plasma membrane, rather than through an actual receptor [3,67].

Circadian Regulation of the RWR Genes

The circadian clock has been shown to regulate a number of environmentally regulated genes [68]. Additionally, cold-induced expression of RWR genes ZAT12, CBF1, CBF2, and CBF3 was recently reported to be gated by the circadian clock [69]. These findings led us to examine globally whether RWR genes are under circadian regulation. Towards this aim, we compared the RWR genes with genes recently identified as circadian regulated [56]. Surprisingly, not only were RWR genes rhythmically expressed but the upregulated and downregulated genes also showed unexpected phase distributions (Figure 7). Forty-two percent of RWR upregulated genes are expressed at subjective dusk while 81% of downregulated RWR genes peak at subjective dawn (p ≤ 0.0001). The circadian regulation of RWR genes may provide a mechanism to anticipate stresses caused by daily environmental changes.

thumbnail
Figure 7. Circadian Regulation of RWR Genes

These data are comprised of circadian-regulated RWR genes. Upregulated RWR genes (blue) peak at dusk while downregulated RWR genes (pink) peak at dawn. Plants were entrained in light/dark cycles for 7 d and then released into constant light. Samples were collected every 4 h after plants were moved to constant light.

https://doi.org/10.1371/journal.pgen.0030172.g007

Identification of a Novel Stress-Responsive cis-Regulatory Element

To begin dissecting the molecular mechanism underpinning the rapid stress response, we examined the promoters of the RWR genes for novel cis-regulatory elements. We identified the six-nucleotide repeat, CGCGTT, which we are terming the Rapid Stress Response Element (RSRE), as significantly overrepresented (58 hits in 47 of the 162 upregulated promoters) in the promoters of upregulated RWR genes. To determine whether the RSRE is sufficient alone to confer stress-responsive transcription, we used luciferase reporter constructs. Four tandem repeats of the RSRE and its consensus flanking sequence were separated by six nucleotides and cloned upstream of the minimal promoter region of the nopaline synthase (NOS) gene and modified luciferase coding region (4xRSRE:LUC). Additionally, to verify that the RSRE was the region conferring stress responsiveness, we mutated three of the six nucleotides in the RSRE (4xmtRSRE:LUC). The wound-induced expression of these constructs was then tested in 24 independent T1 plants to control for differences in expression resulting from the site of transgene insertion. All 24 4xRSRE:LUC transgenic plants exhibited wound-induced luciferase expression (Figure 8). Furthermore, luciferase activity increased immediately following wounding and peaked ∼80 min post wounding. Conversely, no 4xmtRSRE:LUC transgenic plant exhibited luciferase activity before or after wounding. These data demonstrate that the RSRE is sufficient to confer a rapid response to stress.

thumbnail
Figure 8. The RSRE Confers Wound-Induced Reporter Gene Expression

Independent T1 4xRSRE:LUC and 4xmtRSRE:LUC lines were mechanically wounded. Luciferase activity was then monitored in the wounded leaf. Data are means of n = 24 ± SEM.

https://doi.org/10.1371/journal.pgen.0030172.g008

Both abiotic and biotic stresses appear to share common signaling components with the RWR. We were therefore interested in whether the RSRE confers a rapid response to a range of stresses. To enable accurate quantification of the stress response we used a homozygous T3 4xRSRE:LUC line. Because RWR genes respond to both OGAs and IR (Figures 3 and 7), we tested the expression of 4xRSRE:LUC following treatment with these biotic elicitors. In the local leaf where the wound site was treated with OGA or IR, no induction over that of H2O treatment was observed (Figure 9A and 9B). Notably, in systemic tissues, a statistically significant (p < 0.05) synergistic enhancement of luciferase activity was detected in OGA- and IR-treated plants when compared to H2O treatment (Figure 9A and 9C). To further demonstrate that the RSRE responds to biotic stress, 4xRSRE:LUC plants were challenged with cabbage loopers and B. cinerea. When 4xRSRE:LUC plants were exposed to both forms of biotic stress, luciferase activity was induced, whereas no activity was observed in vector control lines (Videos S3 and S4). Additionally, in B. cinerea–infected plants, a low level of transient luciferase activity was first observed in the inoculated leaf. Luciferase activity was then observed at a greater level in systemic tissues (Video S4). These data clearly demonstrate that the RSRE responds to biotic stress.

thumbnail
Figure 9. The RSRE Responds to Biotic Elicitors

(A) Image of an individual 4xRSRE:LUC transgenic plant over time that was wounded and treated with OGA, IR, or H2O.

(B) Local expression in the treated leaves of 4xRSRE:LUC transgenic plants.

(C) Systemic expression of 4xRSRE:LUC in the shoot apex. Data are means of 12 ± SEM.

https://doi.org/10.1371/journal.pgen.0030172.g009

While the RSRE responds to the abiotic stress of mechanical wounding, we wished to further demonstrate the role of the RSRE in response to abiotic stress. Towards this aim, we exposed 4xRSRE:LUC expressing plants to 5 °C. Plants were then removed from cold treatment at the indicated time for imaging. Additionally, control 4xRSRE:LUC plants were also kept at 22 °C in equivalent light conditions and moved similarly to cold-treated plants to ensure that transfer to the imaging chamber did not result in induced luciferase activity. Induction of luciferase activity was observed after ∼2 h of cold treatment (Figure 10A and 10B). Furthermore, luciferase activity peaked after 5 h of cold treatment and then decreased towards basal expression levels. Notably, 4xRSRE:LUC expression was also observed in the roots of cold-treated plants. To ensure that the 4xRSRE:LUC plants were still competent to express luciferase, they were mechanically wounded after 120 h of cold treatment. Both cold- and 22 °C-treated 4xRSRE:LUC plants exhibited luciferase induction following wounding (unpublished data). These data demonstrate that the RSRE is cold responsive and that initial signaling events leading to a cold response dampen even in continuous exposure to cold.

thumbnail
Figure 10. The RSRE Confers Cold-Induced Reporter Gene Expression

(A) Image of an individual 4xRSRE:LUC transgenic plant incubated at 5 °C or 22 °C over time.

(B) Luciferase activity of 4xRSRE:LUC plants. Data are means of 16 ± SEM.

https://doi.org/10.1371/journal.pgen.0030172.g010

The rapid and transient response of the RSRE to multiple stress conditions is reminiscent of the yeast GSR promoter element STRE (stress response element; AGGGG) [18]. The STRE is responsible for rapid induction following various treatments such as heat, nitrogen starvation, low external pH, osmotic, and oxidative stress [7073]. Furthermore, even in the presence of continuous stress exposure, STRE-mediated gene induction dampens over time. An increase in unsaturated fatty acids upon stress appears to be responsible for the transient nature of STRE-mediated induction [74,75]. When plants are exposed to abiotic and biotic stresses (cold and P. syringae, respectively), there is an increase in unsaturated fatty acids [9,76]. Upon cold treatment, acyl-lipid desaturases are the enzymes that most efficiently introduce double bonds in membrane lipids, which results in the increased level of unsaturation [9]. Similar to cold treatment, two of the RWR upregulated genes are acyl-lipid desaturases (ADS1 and ADS2), which may increase the unsaturation of membranes upon wounding. It is therefore tempting to speculate that, as with the STRE in yeast, the increase in unsaturated membrane lipids resulting from both abiotic and biotic stresses may mediate the transient induction of RSRE-driven reporter gene expression.

Conclusions

We have shown that 5 min of mechanical stress is a sufficient amount of time for the plant to perceive the stress and mount a robust transcriptional response. The rapid transcriptional response to mechanical wounding shares a large overlap with both abiotic and biotic stresses and may therefore represent the initial GSR of Arabidopsis. In support of this view, the RWR upregulated genes comprised 25% of the genes identified as potential GSR genes via analysis of the AtGenExpress abiotic and biotic stress datasets [28]. Additionally, in mammalian cells, physical stress to membranes during osmotic and UV radiation stress result in the nonspecific clustering of growth factor receptor tyrosine kinases and cytokine receptors [18,77]. A similar nonspecific clustering of receptors during mechanical wounding and other environmental stresses may underlie the GSR of plants. We also show that the RWR genes are circadian regulated with consolidated phases of peak expression. Circadian regulation of RWR genes, which likely encompass initial components of the GSR, may enable plants to anticipate daily environmental changes and mount a general defense against these changes. Finally, we identified a cis-regulatory element (RSRE) overrepresented in the promoters of RWR genes. The RSRE confers a rapid and transient response similar to the yeast GSR promoter element (STRE) and is a novel GSR cis-regulatory element in plants.

Since the RWR genes likely represent initial components of the GSR, they provide a valuable resource of candidate genes for engineering of multi-stress resistance. Similarly, the RSRE, which responds rapidly and transiently to abiotic and biotic stresses, may prove useful as a synthetic element for engineering of multi-stress tolerance. Finally, use of the RSRE in yeast one-hybrid and the 4xRSRE:LUC line for mutant screens should help elucidate the upstream mechanisms of stress perception and initial signal transduction.

Materials and Methods

Plant growth conditions and treatment.

Arabidopsis (Arabidopsis thaliana) ecotype Columbia-0 plants were grown in a 16 h light/8 h dark photoperiod at 22 °C. All experiments were conducted on 3-wk-old soil-grown plants unless otherwise noted. All rosette leaves (unless otherwise noted) were mechanically wounded one to two times with a hemostat (resulting in ∼20% of the leaf being damaged). Mechanical wounding was performed 4–6 h after dawn.

Cloning.

For cloning of PERF18:LUC, the 1.4 kb upstream of the translation start site of ERF#018 (At1g74930) was PCR amplified using primers listed in Table S2. For cloning of PCAF1-like:LUC, the 2-kb upstream region of the translation start site of CAF1-like (AT3G44260) was PCR amplified using primers listed in Table S2. The PCR amplicons were cloned into the pENTR/D-TOPO vector and subcloned into the Gateway destination vector pBGWL7 [78] by an LR reaction (Invitrogen, Carlsbad, CA). 4xRSRE:LUC (5′-cataaCGCGTTtttagatatcataaCGCGTTtttaggatccataaCGCGTTtttatctagaataaCGCGTTtttac-3′) and 4xmtRSRE:LUC (5′-cataaCATGCTtttagatatcataaCATGCTtttaggatccataaCATGCTtttatctagaataaCATGCTtttac-3′) constructs were created by cloning into the SacI/XhoI sites of pATM-Nos [79]. Transformations were performed into Columbia-0 plants by floral dip using Agrobacterium tumefaciens strain GV3101 [80].

Arabidopsis thaliana oligoarray and preparation of labeled targets for hybridization.

The Arabidopsis 2.0 oliogoarray chip containing 60-mer oligos and representing a total of 21,500 probes (TAIR ATH1 v4.0) was obtained from Agilent (G4137A; Wilmington, Delaware). Total RNA from leaf tissue was isolated using TRIzol reagent (Invitrogen) following the manufacturer's suggested protocol. Prior to hybridizations, the quality and quantity of the total RNA sample was confirmed by running 10-ng samples on an Agilent bioanalyzer (RNA chip), and by using a spectrophotometer. The oligoarray hybridization experiments utilized three biological replicates of pooled plants (∼40 plants per pool). Each biological replicate was comprised of two technical replicates with dye reversal. Total RNA (500 ng) was used as a template for cRNA production and cyanine dyes were incorporated using the Agilent low RNA input linear amp kit (Agilent). Normal yields from 500 ng of total RNA input using a 4-h in vitro transcription were 15 μg cRNA (15 pmol cyanine dye incorporated/μg cRNA).

Array hybridization and scanning.

One microgram of labeled cRNA (cy3- and cy5-labeled sample) was diluted to 175 μl and defragmented at 60 °C for 30 min following the Agilent hybridization protocols (Agilent). Defragmented samples were diluted to 500 μl (30% formamide final concentration) and hybridized for 20 h at 40 °C. Arrays were washed and dried and scanned on an Agilent G2565BA microarray scanner [81]. The raw TIFF images were analyzed using the Agilent Feature Extraction software v. 8.1 using the recommended default settings.

Microarray data analysis procedures.

The intensities of Cy3- and Cy5-labeled probes were normalized by comparing signal intensities of housekeeping genes (positive controls) for both dyes and using the determined ratio as a correction factor for differences in labeling efficiencies (Agilent Feature Extraction v. 8.1 software). The genes that had valid signal in all six replicates were exported to Rosetta Resolver software and analyzed according to the manufacturer's instructions (Seattle, Washington). The normalized values were used to calculate the ratio of channel intensities (Cy5/Cy3), which were then log10 transformed. The transformed ratios were plotted in a scatter plot (cy5/cy3) A ±1.7-fold increase or decrease in signal intensity (p-chance value, 0.01) “Signature” Features (>1.5-fold deviation from the median, p-chance value 0.01) was exported into Microsoft Excel for further analysis. In Excel, data were sorted and genes with a fold-change ≥ 2.0 were selected as differentially regulated. Additionally, differentially expressed genes were binned based on signal intensity.

RT-qPCR.

RT-qPCR analysis was conducted based on Yamagishi et al. [82] Total RNA was isolated by TRIzol extraction (Life Technologies, Grand Island, NY) and treated with DNase, MseI, and DdeI to control for DNA contamination. One microgram of RNA was reverse transcribed using Superscript III (Invitrogen, Carlsbad, California). PCRs were conducted using 12.5 μl of SYBR Green mix (40 mM Tris HCL pH 8.4, 100 mM KCl, 6 mM MgCl2, 8 % glycerol, 20 nM fluorescein, 0.4× SYBR Green [Molecular Probes, Eugene, Oregon], 100-fold dilution of BSA [New England Biolabs, Beverly, MA], and 1.6 mM dNTP), 2 μl of a 30-fold dilution of the RT reaction, 0.6 U iTaq DNA Polymerase (Bio-Rad Laboratories, Hercules, CA), 9.38 μl H2O, and 0.24 μM of each primer. Reactions were carried out on a Bio-Rad iCycler iQ multicolor real-time detection system (Bio-Rad Laboratories) using a two-step reaction condition (extension temperature was primer specific but was typically ∼60 °C), followed by a melt curve encompassing 80 steps of 0.5 °C from 60 °C to 100 °C. Gene-specific primers were designed using Beacon Designer software (Premier Biosoft Palo Alto, CA) and are listed in Table S2. Primary data analysis was performed with Bio-Rad iCycler iQ software. The Bio-Rad gene expression macro version 1.1 software (Bio-Rad Laboratories) was used to calculate relative RNA levels normalized to an internal control [8385]. The 60S ribosomal protein L14 (At4g27090) described in [82] or the TIP41-like gene (At4g34270) described in [86] was used as an internal control.

Comparison of transcript profiles.

For this analysis, we compared the overlap in gene lists between the RWR genes and published transcript profiles determining circadian-regulated genes as well as abiotic and biotic stress-responsive genes. The statistical significance of the observed overlap in transcript profiles was then analyzed using a recently described empirical permutation test [56] based on sampled randomization testing [87] with a p-value cutoff of p < 0.0001.

Luciferase imaging.

For luciferase imaging, 10–14-d-old plants grown on plates containing 1× Murashige and Skoog basal salt mixture (Sigma) were utilized. Plants were sprayed with 2.5 mM luciferin (Promega, Madison, WI) in 0.001% Triton X-100 ∼16–20 h prior to treatment. Mechanical wounding was performed on a single leaf per plant. For elicitor application, 5 μl of 30 mg/ml OGA, 1μl of cabbage looper (Trichoplusia ni) regurgitant, or 5 μl of sterile ddH2O was applied to a single leaf. For cold treatment, plants were placed in a 5 °C chamber under low light. Control plants were kept at 22 °C in equivalent light conditions and handled similarly to cold-treated plants to ensure that transfer to the imaging chamber did not result in induced luciferase activity. Plants were then removed from cold treatment for imaging and returned to either cold or 22 °C. Five microliters of Botrytis cinerea isolate KB2 in 1/2× grape juice was inoculated on a single leaf at a concentration of 500 spores/μl [88,89]. Plants were imaged using an Andor DU434-BV CCD camera (Andor Technology, South Windsor, CT). For image acquisition, PERF18:LUC or PCAF1-like:LUC plants were exposed for 20 min while 4xRSRE:LUC plants were exposed for 5 min. Luciferase activity was quantified for a defined area (leaf, shoot apex, or whole plant) as mean counts pixel−1 exposure time−1 using Andor Solis image analysis software (Andor Technology, South Windsor, CT). For statistical analysis of treatment effects, the area under the curve was calculated and compared by Kruskal-Wallis one way ANOVA on ranks, with pairwise multiple comparisons (Student-Newman-Keuls Method), using Sigma Stat v3.5 (San Jose, CA).

Oligouronide preparation.

A 1% (w/v) solution of citrus pectin (Sigma-Aldrich) in 0.5 M HCl was refluxed for 3 h at a rolling boil. Following cooling, the sample was neutralized with NaOH, decolorized using activated charcoal, and dialyzed against water using 6,000 molecular weight tubing. The dialyzed sample was then lyophilized [90].

Detection of statistically significant promoter motifs.

Promoter sequences were defined as a fixed distance (2 kb for the purpose of motif detection) upstream of the annotated translation start codon, as described in Hudson and Quail [90]. An enhanced enumerative motif recognition algorithm was developed based on that described by Hudson and Quail [90]. The enhanced method does not require exact motif matching, and permits degenerate bases to occur in any number at any position through the search motif by searching for and enumerating all possible permutations of bases within the specified motif size limits including wildcard (N) bases. Significant associations between promoter gene lists and all permutations of motifs are then detected by comparison of per-promoter motif abundance between the promoters of the target coregulated gene list and the promoters of all genes in the Arabidopsis genome. The per-promoter binomial test described [91] is used to rank motifs in order of significance, with the exception that a more efficient factorial handling algorithm was used to perform the binomial test on every motif present in the list of coregulated promoters, rendering the preliminary chi-square filtering step unnecessary. Finally, motifs with related sequence are aligned and automatically clustered for output based on nucleotide-level identity. Programs were written in the Perl programming language utilizing the GMP Multiple Precision Arithmetic Library [92]. The Arabidopsis promoter motif search program is available for use via a web interface at http://stan.cropsci.uiuc.edu/tools.php. Source code is available from MEH on request.

Supporting Information

Table S1. List of Upregulated and Downregulated RWR Genes

https://doi.org/10.1371/journal.pgen.0030172.st001

(88 KB XLS)

Table S2. Primers Used for RT-qPCR and Cloning

https://doi.org/10.1371/journal.pgen.0030172.st002

(38 KB XLS)

Video S1. Movie of a PERF18:LUC-Expressing Plant Challenged with Cabbage Loopers for 18 h

https://doi.org/10.1371/journal.pgen.0030172.sv001

(595 KB MOV)

Video S2. Movie of PCAF1-like:LUC-Expressing Plants Challenged with Cabbage Loopers for 18 h

https://doi.org/10.1371/journal.pgen.0030172.sv002

(1.4 MB MOV)

Video S3. Movie of 4xRSRE:LUC-Expressing Plants Challenged with Cabbage Loopers for 16 h

https://doi.org/10.1371/journal.pgen.0030172.sv003

(659 KB MOV)

Video S4. Movie of a 4xRSRE:LUC-Expressing Plant Inoculated with B. cinere

Induction of 4xRSRE:LUC activity appears first in the inoculated leaf and subsequently spreads systemically. Images were taken for 5 min every 30 min. This movie consists of images 24–48 h post inoculation.

https://doi.org/10.1371/journal.pgen.0030172.sv004

(3.0 MB MOV)

Accession Numbers

The data discussed in this publication have been deposited in the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO, http://www.ncbi.nlm.nih.gov/geo/) and are accessible through GEO Series accession number GSE8740.

Acknowledgments

We thank Dior Kelley and Peter Quail for their critical review of this manuscript, Heather Rowe and Dan Kliebenstein for providing B. vcinerea, and John Labavitch for providing the OGAs.

Author Contributions

JWW and KD conceived and designed the experiments and wrote the paper. JWW, SC, and RK performed the experiments. JWW, SC, MEH, MFC, and GB analyzed the data. SC, MEH, MFC, SLH, and KD contributed reagents/materials/analysis tools.

References

  1. 1. Fujita M, Fujita Y, Noutoshi Y, Takahashi F, Narusaka Y, et al. (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.
  2. 2. Bostock RM (2005) Signal crosstalk and induced resistance: straddling the line between cost and benefit. Annu Rev Phytopathol 43: 545–580.
  3. 3. Vorwerk S, Somerville S, Somerville C (2004) The role of plant cell wall polysaccharide composition in disease resistance. Trends Plant Sci 9: 203–209.
  4. 4. Schmelz EA, Carroll MJ, LeClere S, Phipps SM, Meredith J, et al. (2006) Fragments of ATP synthase mediate plant perception of insect attack. Proc Natl Acad Sci USA 103: 8894–8899.
  5. 5. Schmelz EA, LeClere S, Carroll MJ, Alborn HT, Teal PE (2007) Cowpea (Vigna unguiculata) chloroplastic ATP synthase is the source of multiple plant defense elicitors during insect herbivory. Plant Physiol 144: 793–805.
  6. 6. Ellis C, Turner JG (2001) The Arabidopsis mutant cev1 has constitutively active jasmonate and ethylene signal pathways and enhanced resistance to pathogens. Plant Cell 13: 1025–1033.
  7. 7. Ellis C, Karafyllidis I, Wasternack C, Turner JG (2002) The Arabidopsis mutant cev1 links cell wall signaling to jasmonate and ethylene responses. Plant Cell 14: 1557–1566.
  8. 8. Chinnusamy V, Zhu J, Zhu J-K (2006) Gene regulation during cold acclimation in plants. Physiol Plant 126: 52–61.
  9. 9. Murata N, Los DA (1997) Membrane fluidity and temperature perception. Plant Physiol 115: 875–879.
  10. 10. Xiong L, Schumaker KS, Zhu J-K (2002) Cell signaling during cold, drought, and salt stress. Plant Cell 14: S165–183.
  11. 11. Testerink C, Munnik T (2005) Phosphatidic acid: a multifunctional stress signaling lipid in plants. Trends Plant Sci 10: 368–375.
  12. 12. Apel K, Hirt H (2004) Reactive oxygen species: metabolism, oxidative stress, and signal transduction. Annu Rev Plant Biol 55: 373–399.
  13. 13. Mishra G, Zhang W, Deng F, Zhao J, Wang X (2006) A bifurcating pathway directs abscisic acid effects on stomatal closure and opening in Arabidopsis. Science 312: 264–266.
  14. 14. Orozco-Cardenas ML, Narvaez-Vasquez J, Ryan CA (2001) Hydrogen peroxide acts as a second messenger for the induction of defense genes in tomato plants in response to wounding, systemin, and methyl jasmonate. Plant Cell 13: 179–191.
  15. 15. de Jong CF, Laxalt AM, Bargmann BOR, de Wit PJGM, Joosten MHAJ, et al. (2004) Phosphatidic acid accumulation is an early response in the Cf-4/Avr4 interaction. Plant J 39: 1–12.
  16. 16. Leon J, Rojo E, Sanchez-Serrano JJ (2001) Wound signalling in plants. J Exp Bot 52: 1–9.
  17. 17. Torres MA, Dangl JL (2005) Functions of the respiratory burst oxidase in biotic interactions, abiotic stress and development. Curr Opin Plant Biol 8: 397–403.
  18. 18. Kultz D (2005) Molecular and evolutionary basis of the cellular stress response. Annu Rev Physiol 67: 225–257.
  19. 19. Mikkelsen R, Wardman P (2003) Biological chemistry of reactive oxygen and nitrogen and radiation-induced signal transduction mechanisms. Oncogene 22: 5734–5754.
  20. 20. Bohnert HJ, Gong Q, Li P, Ma S (2006) Unraveling abiotic stress tolerance mechanisms - getting genomics going. Curr Opin Plant Biol 9: 180–188.
  21. 21. Eulgem T (2005) Regulation of the Arabidopsis defense transcriptome. Trends Plant Sci 10: 71–78.
  22. 22. Durrant WE, Rowland O, Piedras P, Hammond-Kosack KE, Jones JDG (2000) cDNA-AFLP reveals a striking overlap in race-specific resistance and wound response gene expression profiles. Plant Cell 12: 963–977.
  23. 23. Farmer EE (2000) Adding injury to insult: pathogen detection and responses. Genome Biol 1: REVIEWS1012.
  24. 24. Cheong YH, Chang H-S, Gupta R, Wang X, Zhu T, et al. (2002) Transcriptional profiling reveals novel interactions between wounding, pathogen, abiotic stress, and hormonal responses in Arabidopsis. Plant Physiol 129: 661–677.
  25. 25. Kilian J, Whitehead D, Horak J, Wanke D, Weinl S, et al. (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.
  26. 26. Swindell WR (2006) The association among gene expression responses to nine abiotic stress treatments in Arabidopsis thaliana. Genetics 174: 1811–1824.
  27. 27. Walther D, Brunnemann R, Selbig J (2007) The regulatory code for transcriptional response diversity and its relation to genome structural properties in A. thaliana. PLoS Genet 3: e11.
  28. 28. Ma S, Bohnert H (2007) Integration of Arabidopsis thaliana stress-related transcript profiles, promoter structures, and cell-specific expression. Genome Biol 8: R49.
  29. 29. Vinocur B, Altman A (2005) Recent advances in engineering plant tolerance to abiotic stress: achievements and limitations. Curr Opin Biotechnol 16: 123–132.
  30. 30. Rushton PJ, Reinstadler A, Lipka V, Lippok B, Somssich IE (2002) Synthetic plant promoters containing defined regulatory elements provide novel insights into pathogen- and wound-induced signaling. Plant Cell 14: 749–762.
  31. 31. Wood NT (2002) Synthetic promoters illuminate roles of cis-acting elements in plant defence. Trends Plant Sci 7: 288–288.
  32. 32. Berardini TZ, Mundodi S, Reiser L, Huala E, Garcia-Hernandez M, et al. (2004) Functional annotation of the Arabidopsis genome using controlled vocabularies. Plant Physiol 135: 745–755.
  33. 33. McCormack E, Tsai Y-C, Braam J (2005) Handling calcium signaling: Arabidopsis CaMs and CMLs. Trends Plant Sci 10: 383–389.
  34. 34. Xu X, Chen C, Fan B, Chen Z (2006) Physical and functional interactions between pathogen-induced Arabidopsis WRKY18, WRKY40, and WRKY60 transcription factors. Plant Cell 18: 1310–1326.
  35. 35. Yang H, Li Y, Hua J (2006) The C2 domain protein BAP1 negatively regulates defense responses in Arabidopsis. Plant J 48: 238–248.
  36. 36. Arnholdt-Schmitt B (2004) Stress-induced cell reprogramming. A role for global genome regulation? Plant Physiol 136: 2579–2586.
  37. 37. Zhou C, Zhang L, Duan J, Miki B, Wu K (2005) HISTONE DEACETYLASE19 is involved in jasmonic acid and ethylene signaling of pathogen response in Arabidopsis. Plant Cell 17: 1196–1204.
  38. 38. Vlachonasios KE, Thomashow MF, Triezenberg SJ (2003) Disruption mutations of ADA2b and GCN5 transcriptional adaptor genes dramatically affect Arabidopsis growth, development, and gene expression. Plant Cell 15: 626–638.
  39. 39. Asai T, Tena G, Plotnikova J, Willmann MR, Chiu W-L, et al. (2002) MAP kinase signalling cascade in Arabidopsis innate immunity. Nature 415: 977–983.
  40. 40. Wang H, Ngwenyama N, Liu Y, Walker JC, Zhang S (2007) Stomatal development and patterning are regulated by environmentally responsive mitogen-activated protein kinases in Arabidopsis. Plant Cell 19: 63–73.
  41. 41. Schweighofer A, Hirt H, Meskiene I (2004) Plant PP2C phosphatases: emerging functions in stress signaling. Trends Plant Sci 9: 236–243.
  42. 42. Tucker M, Valencia-Sanchez MA, Staples RR, Chen J, Denis CL, et al. (2001) The transcription factor associated Ccr4 and Caf1 proteins are components of the major cytoplasmic mRNA deadenylase in Saccharomyces cerevisiae. Cell 104: 377–386.
  43. 43. Viswanathan P, Ohn T, Chiang Y-C, Chen J, Denis CL (2004) Mouse CAF1 can function as a processive deadenylase/3′-5′-exonuclease in vitro but in yeast the deadenylase function of CAF1 is not required for mRNA Poly(A) removal. J Biol Chem 279: 23988–23995.
  44. 44. Zhu J, Dong C-H, Zhu J-K (2007) Interplay between cold-responsive gene regulation, metabolism and RNA processing during plant cold acclimation. Curr Opin Plant Biol 10: 290–295.
  45. 45. Gutierrez RA, Ewing RM, Cherry JM, Green PJ (2002) Identification of unstable transcripts in Arabidopsis by cDNA microarray analysis: rapid decay is associated with a group of touch- and specific clock-controlled genes. Proc Natl Acad Sci USA 99: 11513–11518.
  46. 46. Nakano T, Suzuki K, Fujimura T, Shinshi H (2006) Genome-wide analysis of the ERF gene family in Arabidopsis and rice. Plant Physiol 140: 411–432.
  47. 47. Stockinger EJ, Gilmour SJ, Thomashow MF (1997) Arabidopsis thaliana CBF1 encodes an AP2 domain-containing transcriptional activator that binds to the C-repeat/DRE, a cis-acting DNA regulatory element that stimulates transcription in response to low temperature and†water†deficit. Proc Natl Acad Sci USA 94: 1035–1040.
  48. 48. Liu Q, Kasuga M, Sakuma Y, Abe H, Miura S, et al. (1998) Two transcription factors, DREB1 and DREB2, with an EREBP/AP2 DNA binding domain separate two cellular signal transduction pathways in drought- and low-temperature-responsive gene expression, respectively, in Arabidopsis. Plant Cell 10: 1391–1406.
  49. 49. Gilmour SJ, Zarka DG, Stockinger EJ, Salazar MP, Houghton JM, et al. (1998) Low temperature regulation of the Arabidopsis CBF family of AP2 transcriptional activators as an early step in cold-induced COR gene expression. Plant J 16: 433–442.
  50. 50. Medina J, Bargues M, Terol J, Perez-Alonso M, Salinas J (1999) The Arabidopsis CBF gene family is composed of three genes encoding AP2 domain-containing proteins whose expression is regulated by low temperature but not by abscisic acid or dehydration. Plant Physiol 119: 463–470.
  51. 51. Lippuner V, Cyert MS, Gasser CS (1996) Two classes of plant cDNA clones differentially complement yeast calcineurin mutants and increase salt tolerance of wild-type yeast. J Biol Chem 271: 12859–12866.
  52. 52. Davletova S, Schlauch K, Coutu J, Mittler R (2005) The zinc-finger protein Zat12 plays a central role in reactive oxygen and abiotic stress signaling in Arabidopsis. Plant Physiol 139: 847–856.
  53. 53. Gomez-Gomez L, Boller T (2000) FLS2: An LRR receptor-like kinase involved in the perception of the bacterial elicitor flagellin in Arabidopsis. Mol Cell 5: 1003–1011.
  54. 54. Kesarwani M, Yoo J, Dong X (2007) Genetic interactions of TGA transcription factors in the regulation of pathogenesis-related genes and disease resistance in Arabidopsis. Plant Physiol 144: 336–346.
  55. 55. Kariola T, Brader G, Helenius E, Li J, Heino P, et al. (2006) EARLY RESPONSIVE TO DEHYDRATION 15, a negative regulator of abscisic acid responses in Arabidopsis. Plant Physiol 142: 1559–1573.
  56. 56. Covington M, Harmer SL (2007) the circadian clock regulates auxin signaling and responses in Arabidopsis. PLoS Biol 5: e222.
  57. 57. Vogel JT, Zarka DG, Van Buskirk HA, Fowler SG, Thomashow MF (2005) Roles of the CBF2 and ZAT12 transcription factors in configuring the low temperature transcriptome of Arabidopsis. Plant J 41: 195–211.
  58. 58. Lee B-h, Henderson DA, Zhu J-K (2005) The Arabidopsis cold-responsive transcriptome and its regulation by ICE1. Plant Cell 17: 3155–3175.
  59. 59. De Vos M, Van Oosten VR, Van Poecke RMP, Van Pelt JA, Pozo MJ, et al. (2005) Signal signature and transcriptome changes of Arabidopsis during pathogen and insect attack. Mol Plant Microbe Interact 18: 923–937.
  60. 60. Little D, Gouhier-Darimont C, Bruessow F, Reymond P (2007) oviposition by pierid butterflies triggers defense responses in Arabidopsis. Plant Physiol 143: 784–800.
  61. 61. Ferrari S, Galletti R, Denoux C, De Lorenzo G, Ausubel FM, et al. (2007) Resistance to Botrytis cinerea induced in Arabidopsis by elicitors is independent of salicylic acid, ethylene, or jasmonate signaling but requires PHYTOALEXIN DEFICIENT3. Plant Physiol 144: 367–379.
  62. 62. Knight MR, Read ND, Campbell AK, Trewavas AJ (1993) Imaging calcium dynamics in living plants using semi-synthetic recombinant aequorins. J Cell Biol 121: 83–90.
  63. 63. Navazio L, Moscatiello R, Bellincampi D, Baldan B, Meggio F, et al. (2002) The role of calcium in oligogalacturonide-activated signalling in soybean cells. Planta 215: 596–605.
  64. 64. Baydoun EAH, Fry SC (1985) The immobility of pectic substances in injured tomato leaves and its bearing on the identity of the wound hormone. Planta 165: 269–276.
  65. 65. Bergey DR, Orozco-Cardenas M, de Moura DS, Ryan CA (1999) A wound- and systemin-inducible polygalacturonase in tomato leaves. Proc Natl Acad Sci USA 96: 1756–1760.
  66. 66. Orozco-Cardenas M, Ryan CA (1999) Hydrogen peroxide is generated systemically in plant leaves by wounding and systemin via the octadecanoid pathway. Proc Natl Acad Sci USA 96: 6553–6557.
  67. 67. Aldington S, McDougall GJ, Fry SC (1991) Structure-activity relationships of biologically active oligosaccharides. Plant Cell Environ 14: 625–636.
  68. 68. Harmer SL, Hogenesch JB, Straume M, Chang H-S, Han B, et al. (2000) Orchestrated transcription of key pathways in Arabidopsis by the circadian clock. Science 290: 2110–2113.
  69. 69. Fowler SG, Cook D, Thomashow MF (2005) Low temperature induction of Arabidopsis CBF1, 2, and 3 is gated by the circadian clock. Plant Physiol 137: 961–968.
  70. 70. Martínez-Pastor M, Marchler G, Schüller C, Marchler-Bauer A, Ruis H, et al. (1996) The Saccharomyces cerevisiae zinc finger proteins Msn2p and Msn4p are required for transcriptional induction through the stress response element (STRE). EMBO J 15: 2227–2235.
  71. 71. Kobayashi N, McEntee K (1993) Identification of cis and trans components of a novel heat shock stress regulatory pathway in Saccharomyces cerevisiae. Mol Cell Biol 13: 248–256.
  72. 72. Marchler G, Schüller C, Adam G, Ruis H (1993) A Saccharomyces cerevisiae UAS element controlled by protein kinase A activates transcription in response to a variety of stress conditions. EMBO J 12: 1997–2003.
  73. 73. Schüller C, Brewster J, Alexander M, Gustin M, Ruis H (1994) The HOG pathway controls osmotic regulation of transcription via the stress response element (STRE) of the Saccharomyces cerevisiae CTT1 gene. EMBO J 13: 4382–4389.
  74. 74. Chatterjee MT, Khalawan SA, Curran BP (1997) Alterations in cellular lipids may be responsible for the transient nature of the yeast heat shock response. Microbiology 143: 3063–3068.
  75. 75. Chatterjee MT, Khalawan SA, Curran BPG (2000) Cellular lipid composition influences stress activation of the yeast general stress response element (STRE). Microbiology 146: 877–884.
  76. 76. Shah J (2005) Lipids, lipases, and lipid-modifying enzymes in plant disease resistance. Annu Rev Phytopathol 43: 229–260.
  77. 77. Rosette C, Karin M (1996) Ultraviolet light and osmotic stress: activation of the JNK cascade through multiple growth factor and cytokine receptors. Science 274: 1194–1197.
  78. 78. Karimi M, De Meyer B, Hilson P (2005) Modular cloning in plant cells. Trends Plant Sci 10: 103–105.
  79. 79. Harmer SL, Kay SA (2005) Positive and negative factors confer phase-specific circadian regulation of transcription in Arabidopsis. Plant Cell 17: 1926–1940.
  80. 80. Clough SJ, Bent AF (1998) Floral dip: a simplified method for Agrobacterium-mediated transformation of Arabidopsis thaliana. Plant J 16: 735–743.
  81. 81. Coughlan S, Agrawal V, Meyers BC (2004) A comparison of global gene expression measurement technologies in Arabidopsis thaliana. Comp Funct Genomics 5: 245–252.
  82. 82. Yamagishi K, Nagata N, Yee KM, Braybrook SA, Pelletier J, et al. (2005) TANMEI/EMB2757 encodes a WD repeat protein required for embryo development in Arabidopsis. Plant Physiol 139: 163–173.
  83. 83. Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2-[Delta][Delta]CT method. Methods 25: 402–408.
  84. 84. Pfaffl MW (2001) A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 29.
  85. 85. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, et al. (2002) Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol 3: RESEARCH0034.
  86. 86. 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.
  87. 87. Sokal R, Rohlf F (1995) Biometry: the principles and practice of statistics in biological research. New York: W. H. Freeman and Co. pp. 808–810.
  88. 88. Kliebenstein DJ, Rowe HC, Denby KJ (2005) Secondary metabolites influence Arabidopsis/Botrytis interactions: variation in host production and pathogen sensitivity. Plant J 44: 25–36.
  89. 89. Rowe H, Kliebenstein DJ (2007) Elevated genetic variation within virulence-associated Botrytis cinerea polygalacturonase loci. Mol Plant Microbe Interact. In Press.
  90. 90. Campbell AD, Labavitch JM (1991) Induction and regulation of ethylene biosynthesis by pectic oligomers in cultured pear cells. Plant Physiol 97: 699–705.
  91. 91. Hudson ME, Quail PH (2003) Identification of promoter motifs involved in the network of phytochrome a-regulated gene expression by combined analysis of genomic sequence and microarray data. Plant Physiol 133: 1605–1616.
  92. 92. Granlund T (2004) GMP, the GNU multiple precision arithmetic library. Available: http://gmplib.org/. Accessed 14 September 2007.