Open Access Research article

Co-expression and promoter content analyses assign a role in biotic and abiotic stress responses to plant natriuretic peptides

Stuart Meier13, René Bastian1, Lara Donaldson2, Shane Murray1, Vladimir Bajic3 and Chris Gehring1*

Author Affiliations

1 Department of Biotechnology, University of the Western Cape, Private Bag X17, Cape Town - Bellville 7535, South Africa

2 Department of Molecular and Cell Biology, University of Cape Town, Private Bag, Rondebosch 7701, South Africa

3 South African National Bioinformatics Institute, University of the Western Cape, Private Bag X17, Cape Town - Bellville 7535, South Africa

For all author emails, please log on.

BMC Plant Biology 2008, 8:24  doi:10.1186/1471-2229-8-24

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2229/8/24


Received:10 September 2007
Accepted:29 February 2008
Published:29 February 2008

© 2008 Meier et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Plant natriuretic peptides (PNPs) are a class of systemically mobile molecules distantly related to expansins. While several physiological responses to PNPs have been reported, their biological role has remained elusive. Here we use a combination of expression correlation analysis, meta-analysis of gene expression profiles in response to specific stimuli and in selected mutants, and promoter content analysis to infer the biological role of the Arabidopsis thaliana PNP, AtPNP-A.

Results

A gene ontology analysis of AtPNP-A and the 25 most expression correlated genes revealed a significant over representation of genes annotated as part of the systemic acquired resistance (SAR) pathway. Transcription of these genes is strongly induced in response to salicylic acid (SA) and its functional synthetic analogue benzothiadiazole S-methylester (BTH), a number of biotic and abiotic stresses including many SA-mediated SAR-inducing conditions, as well as in the constitutive SAR expressing mutants cpr5 and mpk4 which have elevated SA levels. Furthermore, the expression of AtPNP-A was determined to be significantly correlated with the SAR annotated transcription factor, WRKY 70, and the promoters of AtPNP-A and the correlated genes contain an enrichment in the core WRKY binding W-box cis-elements. In constitutively expressing WRKY 70 lines the expression of AtPNP-A and the correlated genes, including the SAR marker genes, PR-2 and PR-5, were determined to be strongly induced.

Conclusion

The co-expression analyses, both in wild type and mutants, provides compelling evidence that suggests AtPNP-A may function as a component of plant defence responses and SAR in particular. The presented evidence also suggests that the expression of AtPNP-A is controlled by WRKY transcription factors and WRKY 70 in particular. AtPNP-A shares many characteristics with PR proteins in that its transcription is strongly induced in response to pathogen challenges, it contains an N-terminal signalling peptide and is secreted into the extracellular space and along with PR-1, PR-2 and PR-5 proteins it has been isolated from the Arabidopsis apoplast. Based on these findings we suggest that AtPNP-A could be classified as a newly identified PR protein.

Background

Natriuretic Peptide (NP) systems have been identified in mammals, fish, amphibians, birds and reptiles. NPs and their receptors are commonly associated with organs involved in cardiac and osmoregulatory homeostasis. In amphibians, birds and fish, NPs have been shown to play a critical role in the regulation of blood fluid volume and composition [1].

The first indication that NPs function in plants came from radio-immuno assays on plant tissue extracts from Florida beauty [2] and it was shown that the rate of transpiration, solute flow and solute uptake in carnation and chrysanthemum was rapidly and significantly increased after exogenous application of synthetic human atrial NP (ANP) [3]. Subsequently it was demonstrated that rat ANP can induce stomatal opening in a concentration dependent manner [4] and this effect appears to be dependent on the intracellular second messenger cGMP (guanosine 3',5'-cyclic monophosphate) since it is inhibited by the guanylate cyclase inhibitor LY 83583, but can be induced by the synthetic cell permeant cGMP analogue 8-Br-cGMP [5-7]. Binding experiments of ANP to isolated leaf membranes provide evidence for specific receptor ligand interactions [8].

A plant NP (PNP) from A. thaliana (AtPNP-A) and several closely related sequences in different species have since been identified [9,10]. AtPNP-A, its most closely related sequence AtPNP-B, and orthologues in other higher plant species, share a family-45 glucosidase domain with the cell wall loosening expansins [11] and are related to expansins on the basis of this structural homology [9,12]. AtPNP-A (At2g18660) is a small protein of 126 amino acids in length (MW: 14016 kD; pI: 9.22) that is encoded by a gene with a single intron of 100 bp. The region most conserved between PNPs from different plant species has also been shown to be the key to its physiological activity [13]. Evidence for systemic mobility of PNPs comes from the structure and processing of the molecules [9]. The protein contains an N-terminal 24 amino acid signal peptide (MW: 2249) that directs the molecule into the extracellular space and PNPs that are recognized by anti-human atrial natriuretic polypeptide rabbit serum have been localised in situ in conductive tissue [14] and were isolated from xylem sap [15], and proteomics studies have identified the AtPNP-A protein in the apoplastic space in A. thaliana [16].

A number of physiological and biochemical studies have implicated AtPNP-A with a role in the regulation of ion and solute homeostasis. Immuno-reactive PNP (irPNP) extracts and recombinant AtPNP-A have been shown to induce swelling in leaf mesophyll protoplasts [17] and in protoplasts isolated from Arabidopsis cell suspension cultures respectively [15]. Further, irPNP rapidly and specifically induced a transient elevation of cGMP levels in the conductive stele tissue of maize roots [6] and in stomatal guard cell protoplasts [7] and recently recombinant AtPNP-A was shown to stimulate protoplast swelling in a cGMP dependent manner [18]. IrPNPs also modulate ion fluxes across plant membranes [19] and recombinant AtPNP-A induced spatially dependent H+, K+ and Na+ fluxes in A. thaliana roots [20]. Endogenous levels of irPNP are increased in response to NaCl stress in whole-plants and in Arabidopsis suspension culture cells exposed to high salt or osmoticum [15]. Collectively, these studies indicate that PNP-like molecules may function as extracellular signalling molecules that directly affect water and solute transport in response to stress. Based on biochemical and physiological data we propose mechanisms of action for AtPNP-A at the cellular level as summarised in Figure 1 (adapted from [21]).

thumbnailFigure 1. Model of AtPNP-A action at the cellular level. The model proposes that AtPNP-A can dock to receptor-like molecules that directly act as particulate guanylyl cyclases (pGCs) or indirectly activate soluble GCs (sGCs). GCs catalyse the reaction from GTP to cGMP. The latter acts as second messenger affecting cytosolic Ca2+ levels, modulating ion channels, activating phosphorylation through kinases and influences the transcriptome. Phosphodiesterases (PDEs) in turn metabolise cGMP to GMP (adapted from [21]).

Despite an increasing body of physiological and biochemical data [21], the biological role of this systemically mobile peptide has remained elusive. In order to infer a biological role for AtPNP-A, we made use of the large repositories of A. thaliana microarray data to study the expression profiles of AtPNP-A and the 25 most expression correlated genes in response to various treatments as well as in mutants. We further analysed the promoters of these genes for known regulatory motifs. The results of our study predict a function for AtPNP-A in plant abiotic and biotic stress responses, and in particular in systemic acquired resistance (SAR). Furthermore, we demonstrate how computational analyses that link regulatory potential as encoded by promoter elements and expression data can provide novel insights into the function of a specific gene as well as groups of genes.

Results and Discussion

Expression Correlation and GO Analyses

In the first step of the analyses we extracted and ranked the 25 genes whose expressions are most tightly correlated with AtPNP-A (Table 1). The moderate correlation (r) values of the listed genes (maximum r = 0.73) may reflect that the expression of AtPNP-A is subject to complex combinatorial control via multiple promoter motifs with complex inputs from multiple, potentially antagonistic, signalling pathways.

Table 1. List of genes that are expression correlated with AtPNP-A (At2g18660)

In order to identify a functional role of AtPNP-A, the correlated genes were analysed in FatiGO+ [22-24] to identify any bias in GO functional annotation terms in the correlated list (list 1) compared to the remainder of the A. thaliana genome [see Additional file 1]. In the GO search category of biological process there is a significant (Family Wise Error Rate – FWER) adjusted p-value) enrichment in genes involved in biotic defence responses at a number of levels. The most notable bias being at level 8 with a significant (adjusted p-value 0.0000038) enrichment in genes involved in SAR (Table 1). SAR is an inducible plant defense response against local pathogen infection that gives rise to a systemic long lasting resistance to a broad range of virulent pathogens [25]. The SAR response is characterised by the accumulation of endogenous salicylic acid (SA) in infected tissues and later in distal uninfected tissues with a subsequent induction of a select group of pathogenesis-related genes (PR genes) [26].

Additional file 1. GO analysis of AtPNP-A expression correlated genes. List of significantly enriched GO terms associated with AtPNP-A (At2g18660) expression correlated genes in FatiGO+.

Format: PDF Size: 54KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

The enrichment in SAR annotated genes in our list is particularly striking considering that in FatiGO the entire A. thaliana genome contains only 21 annotated SAR genes and four of these are present in our list of correlated genes. The four correlated SAR genes include NIMIN1 (At1g02450; r = 0.61) that is involved in the transcriptional regulation of PR genes [26], PR-1 (At2g14610; r = 0.63); PR-2 (At3g57260; r = 0.73) and PR-5 (At1g75040; 0.61) whose expression is commonly used as diagnostic markers of the SA dependent SAR response [27]. An extended correlation analysis revealed that an additional 11 SAR annotated genes, including NPR1 (or NIM1; At1g64280; r = 0.52), which is an essential key positive regulator of signal transduction leading to the SAR response and expression of PR proteins, are significantly correlated (9 positive, 2 negative; p < 0.01; bivariate normal distribution) with the expression of AtPNP-A [see Additional file 2]. Other correlated genes in list 1 annotated to be involved in plant defence responses and response to biotic stimuli include a disease resistance family protein containing leucine rich-repeats (At2g32680), a stress responsive gene with triphosphatase activity (At2g04450) and a protease inhibitor (At5g55450). The GO analysis for the cellular component and molecular function category revealed no significant difference in biologically relevant labels between the two lists.

Additional file 2. Extended list of genes expression correlated with AtPNP-A. An extended list containing all genes that are expression correlated with AtPNP-A (At2g18660) and a list of all annotated SAR genes in Arabidopsis thaliana.

Format: XLS Size: 26KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

The results of the Swiss-Prot keyword search also identified a significant enrichment in genes annotated as PR proteins (adjusted p value = 0.002), involved in signalling (adjusted p value = 0.026) and associated with the apoplast (adjusted p value = 0.040) in list 1. It was noted that along with PR-1, PR-2 and PR-5, AtPNP-A is one of six genes in list 1 annotated as having signalling function.

Microarray Expression Profiles

The over representation of genes involved in defence responses, and specifically SAR is consistent with the observation that AtPNP-A and the correlated genes are most highly expressed in microarray experiments where defence responses are elicited. The treatments that induce up-regulation of AtPNP-A and the correlated genes more than two-fold include SA and other SAR inducing conditions as well as a number of abiotic stresses (Figure 2).

thumbnailFigure 2. Expression profile of AtPNP-A and correlated genes in response to selected treatments. The results presented illustrate the fold change (log2) in expression of AtPNP-A, WRKY 70 and WRKY 46 and the average fold change for the 25 correlated genes in response to abiotic stresses (A) and biotic and chemical treatments (B). (A) The treatments were: UV-B shoot 3 h (n = 2); O3 6 h (n = 3); Osmotic stress in the shoot after 3 h (n = 2); K+ starvation in the shoot after 7 days (n = 3); NaCl in the roots after 6 h (n = 2) and cold acclimation after 14 days (n = 3). (B) The treatments were: Erysiphe cichoracearum 3 days after inoculation (n = 4); Erysiphe orontii 3 days after inoculation (n = 2); Phytophthora infestans 1 day after inoculation (n = 3); BTH after 8 h (n = 3); SA after 3 h (n = 2) and cyclohexamide after 3 h (n = 2). Error bars represent standard errors of the mean.

The strong up-regulation of AtPNP-A and correlated genes by SA and benzothiadiazole S-methylester (BTH), a synthetic functional SA analogue [28], is a key indicator that these genes are involved in plant defence and specifically SAR since SA has been shown to be essential and sufficient to induce the SAR response in plants [29] (Figure 2B). In addition to AtPNP-A, all of the 25 correlated genes were significantly (ANOVA p-value < 0.05) up-regulated by more than two fold after 8 h and 24 h treatments with 60 μM BTH (Supplementary Table 3 in [28]) further linking these genes to the SAR defense pathway. Expression of AtPNP-A is also significantly correlated with the isochorismate synthase-1 (ICS-1) gene (At1g74710; r = 0.50) that is critical for SA biosynthesis [29].

The biotic stresses that induced the largest increase in expression of AtPNP-A include infection with the biotrophic pathogens Phytophthora infestans, Erysiphe cichoracearum and Erysiphe orontii that depend on living host tissue for survival (Figure 2B). Activation of AtPNP-A expression and the correlated genes in response to these pathogens is in accord with the literature informing that SA-dependent defenses generally act against biotrophs in contrast to jamonic acid- and ethylene-dependent responses that counteract necrotrophs [30].

The large increases in gene expression induced by ozone and UV-B (+6.86log2 and +8.09log2 respectively for AtPNP-A) is consistent with these genes being part of the SAR response since both these treatments have previously been shown to stimulate SA production and induce the expression of PR genes [31-34].

The expression of AtPNP-A and the correlated genes is also strongly modulated by a number of abiotic stresses including K+ starvation, osmotic, and NaCl stress as well as cold acclimation (Figure 2A). A common element of abiotic stresses is that they decrease water potential [35]. It is noteworthy that the induction of AtPNP-A in response to ion and osmotic stresses is tissue specific with the response to high Na+ being specific to root tissue with little change observed in shoots while both low K+ and high osmolarity induced elevated transcription in shoots only [36]. This is of interest since the AtPNP-A protein has been shown to affect water movement in shoots [37] and protoplasts [17] as well as ion fluxes in roots [20]. It is thus tempting to speculate that AtPNP-A may have a role in maintaining plant water and ion homeostasis under stress conditions.

K+ is the key inorganic ion required in high quantities by plants while Na+ on the other hand is toxic at high concentrations [38]. Na+ is able to compete with K+ ions for uptake and binding sites thus maintaining the correct Na+/K+ ratio in plants is of the utmost importance [39]. Decreases in K+ might cause the plant to take up more Na+ in order to maintain adequate osmotic pressure [40]. Therefore either the increase in cytosolic Na+ or a decrease in osmotic pressure as a consequence of K+ starvation or a combination of both may cause AtPNP-A induction.

Elevated expression of AtPNP-A and correlated genes, particularly defence genes and SAR annotated genes, by abiotic osmotic stresses as well as defence eliciting treatments may well reflect that both types of challenges lead to common homeostatic disturbances which in turn transcriptionally activate a set of common response genes. This concept is supported by several studies that recognise a role of SA in abiotic stresses such as drought, salinity and temperature [41,42] and the accumulation of PR proteins is in fact a common plant response to both abiotic and biotic stresses further highlighting the overlap in biotic and abiotic defence mechanisms [43].

The generation of reactive oxygen species and changes in ion fluxes have been identified as early responses to both abiotic and biotic stresses, including an influx of H+ and Ca2+ and an efflux of K+ and Cl- [35]. AtPNP-A has been shown to modulate H+, Na+ and K+ fluxes [20] thus further implicating AtPNP-A in plant stress responses as do studies which indicate that AtPNP-A signals via the intracellular second messenger cGMP [5,6,18] since cGMP has been shown to be an important signaling molecule in pathogen [44] and osmotic stress responses [45] in plants. It seems particularly relevant that the expression of a gene encoding a cyclic nucleotide-gated channel (CNGC20; At3g17700), which has been shown to be involved in the transport of Ca2+ and K+, and in some cases Na+, across cell membranes is also correlated with that of AtPNP-A (r = 0.60) [see Additional file 2] since ion conductance in these channels is regulated by cGMP as well as Ca2+ and calmodulin. These channels have also been implicated in regulating SA-dependent biotic defense responses [46].

AtPNP-A expression is also correlated with a number of Ca2+ sensing/binding proteins including, the above mentioned CNGC20, calreticulin 3 (At1g08450; r = 0.60), two calmodulin-binding proteins (At1g73800; r = 0.63 and At1g73805; r = 0.588) with family members involved in the induction of plant defense responses (NCBI sequence viewer, pfam07887) and a Ca2+-binding EF hand domain containing protein (At3g47480; r = 0.59). One of the earliest responses to biotic and abiotic stresses is an increase in cytosolic free Ca2+[47] that in turn plays a role in activating the oxidative burst after elicitor treatment [48,49] and is also linked to signaling SA-induced PR gene expression [50]. The expression of Ca2+ sensing molecules is rapidly induced in response to biotic and abiotic stresses [51] and functions to decode Ca2+signatures and/or relay signals to downstream targets, including kinases, which further amplify the Ca2+ signal by inducing downstream phosphorylation cascades [38]. The presence of three kinases (At4g04490; r = 0.63, At4g23150; r = 0.63 and At4g11890; r = 0.61) in amongst the correlated genes (Table 1) is entirely consistent with such a signaling cascade. Moreover, the expression correlation of three stress responsive mitogen-activated protein kinase (MAP kinase) genes MAPKK (At4g26070; r = 0.59), MPK 11 (At1g01560; r = 0.58) and MAPKK (At4g29810; r = 0.57) to AtPNP-A also ties in with the proposed cascades [see Additional file 2]. Activation of MAPKs has indeed been reported after exposure to pathogens [52] as well as a number of abiotic stresses [53].

While transcriptional responses to some stresses, including the osmotic, salt, UV-B and some of the biotic treatments, were measured over multiple time points, the data presented here are generally the earliest time point that induced the largest increase in AtPNP-A expression. The expression of AtPNP-A in some cases showed induction at earlier time points than considered in this study, however, in all cases the expression of AtPNP-A generally increased over time and thus high transcript levels were sustained for the duration of the stress, e.g. five days for E. orontii and 24 h for osmotic, salt and UV-B treatments (data not shown). The UV-B experiment can be distinguished from the other experiments in that the stress was not maintained for the duration of the experiment. Rather, plants were irradiated for 15 minutes before being transferred back to the standard phytochamber conditions until harvest. The expression of AtPNP-A in shoots was elevated at 30 minutes (1.84log2; data not shown), peaked at 3 h (8.09log2, > 250 fold) (Figure 2A) and remained elevated (6.46log2; data not shown) at 24 h after irradiation. This documents that AtPNP-A expression remains very high and sustained after the stress has been removed and thus may indicate that the initial damage inflicted, and not the actual presence of the stress itself, is the driving force for the maintained transcriptional activation.

The increase in expression of AtPNP-A (4.49log2) in response to the protein synthesis inhibitor cycloheximide (CHX) implies that transcription of AtPNP-A can occur independently of de novo protein synthesis and that concurs with the definition of immediate early response genes [54] that have been proposed to play important roles in the early regulation of defence responses [55]. It has been postulated that CHX induces gene expression via dual mechanisms; by preventing synthesis or activation of a short-lived transcriptional repressor or by removing specific labile transcript degrading enzymes [56]. There is evidence that the induced expression of genes encoding secreted proteins, such as AtPNP-A does not require de novo protein synthesis [57]. The ability to rapidly induce expression of AtPNP-A independently of de novo protein synthesis thus implies both an important and early role for this gene in response to specific elicitors.

Common Motifs in cis and Transcription Factors

The common expression profiles of APNP-A and the 25 correlated genes in response to both biotic and abiotic stresses suggests that these genes are under common regulatory control and are thus likely to share common cis-elements in their promoter regions. To reveal aspects of common transcriptional activation we analyzed promoter regions of these genes 1 kb upstream of the predicted transcription start site (TSS) for the presence of known plant cis-elements.

The analysis in POBO [58] indicated that the invariant core TTGAC W-box motif was present in 25/26 of our correlated genes a total number of 78 times at an average of 2.99 copies/promoter compared to the average of 2.24 across all A. thaliana promoters (t-test p-value >0.0001) (Figure 3 and [see Additional file 3]). The analysis in Athena [59] identified that the extended and more stringent TTGAC(A/T) W-box motif was present in 22/26 genes a total of 54 times at an average of 2.08 copies/promoter (p-value = 0.0037; data not shown). Although in Athena this p-value does not qualify the W-box motifs to be enriched in our correlated genes according to the stringent enrichment threshold of <10-4 (Bonferroni correction), it does show that a very high percentage of our genes contain multiple copies of the stringent W-box. Both these promoter analysis methods indicate that multiple copies of the W-box elements are present in a high percentage of our correlated genes with the core TTGAC motifs being significantly enriched compared to expected frequencies in the A. thaliana genome suggesting that they are important regulatory elements in these expression correlated genes.

Additional file 3. Promoter analysis. AtPNP-A (At2g18660) and the expression correlated genes were analysed in POBO and Athena for the presence of W-boxes.

Format: DOC Size: 39KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

thumbnailFigure 3. Frequency of occurrence of the W-box (TTGAC) core motif in artificial clusters generated in POBO for A. thaliana background promoters compared to the promoters of AtPNP-A and the 25 most positively expression correlated genes. The 1 kb upstream promoter sequences of the 26 expression correlated genes were analysed in POBO (see methods and [see Additional file 3]) to determine the frequency of occurrence of the TTGAC W-box core motif. The analysis determined that compared to the A. thaliana background (2.24 copies/promoter), there was a significant (t-test: p-value > 0.0001) enrichment in the frequency of the TTGAC motif in our dataset (2.99 copies/promoter).

In plants, W-box cis-elements are known to bind WRKY TFs [60] indicating that these TFs may be important in regulating the expression of the correlated genes. This is completely consistent with our expression analysis results since the WRKY family of TFs [60] have well established roles in regulating disease responses in plants [61]. In addition, they have also been documented to mediate abiotic plant responses to freezing [62], wounding [63], oxidative stress [64], drought, salinity, cold, and heat [65-67]. In our study, expression of AtPNP-A is moderately correlated with the expression of WRKY 70 (At3g56400; r = 0.60) and WRKY 46 (At2g46400; r = 0.56) [see Additional file 2]. When viewing the expression profiles of WRKY 70 and WRKY 46 genes it is apparent that the various treatments which induced large increases in the expression of the correlated genes in Table 1 also induced marked changes in the expression of the WRKY genes (Figure 2). This links the expression of AtPNP-A and the WRKY genes to common biological responses and raises the possibility that WRKY 70 and WRKY 46 may positively regulate AtPNP-A transcription.

While transcription of WRKY 70 and WRKY 46 is generally strongly induced in response to SAR eliciting treatments, only WRKY 46 is consistently co-expressed with AtPNP-A in response to the abiotic (ion and osmotic) stresses (Figure 2A and 2B). As previously described for AtPNP-A and the correlated genes, the induced expression of WRKY 46 is specific to shoots in response to K+ starvation and osmotic stress and to roots in response to NaCl stress.

The expression correlation between AtPNP-A and the discussed WRKY genes and the overrepresentation of W-boxes in the correlated genes prompted a manual analysis of the promoter of AtPNP-A which revealed the presence of four copies of the core TTGAC W-box motif and two occurrences of more stringent TTGAC(C/T) motif clustered in close proximity (starting at -738 and -775) relative to the predicated TSS. The result of the manual inspection coincided with the results returned form Athena and POBO. Similar frequencies of these motifs were observed in a study of 26 SAR regulated genes (termed PR-1 regulon genes) in which only W-boxes were present in the promoters of all 26 genes at an average of 4.3 copies of the core and 2.1 copies of the more stringent W-box elements within 1100 bp upstream of the predicated TSS [60]. These values represent a significant enrichment in W-boxes since these authors determined that the statistical expectation for a randomly distributed pentamer (TTGAC) was 2.1 copies and for the hexamer (TTGAC(C/T)) 1.1 copies per 1100 bp of promoter [60].

In summary, the presented evidence is entirely consistent with transcription of AtPNP-A and the correlated genes being positively regulated by WRKY TFs. The promoter of AtPNP-A and the correlated genes contain an enrichment of the core W-box motif and expression of AtPNP-A is correlated with two WRKY genes in response to various SAR eliciting and biotic and abiotic stresses. The correlation of WRKY 46 in response to ion and osmotic abiotc stresses was both treatment and tissue specific. In the light of these facts we suggest that the expression of AtPNP-A may be closely regulated by WRKY TFs in response to SAR-inducing and abiotic stresses.

Insights from AtPNP-A Expression in Mutants

The link between AtPNP-A, SA signalling and the WRKY TFs is also supported by expression profiles of AtPNP-A and the correlated genes in mutants including a WRKY 70 over-expresser and various SAR related mutants present in the mutant surveyor in Genevestigator [68].

WRKY 70 is a SAR annotated TF in FatiGO and has been shown to be an essential factor in plant defense responses necessary for the induction of PR gene expression in A. thaliana [28,69]. In a microarray study, a constitutive over-expressor of WRKY 70 was shown to induce constitutive expression of SA induced PR genes and five of our correlated genes, including PR-2 and PR-5, which are widely considered SAR marker genes, correspond to genes in this study that were up-regulated > 2.5 fold compared to controls [69]. While the 8K Affymetrix chip used for this study did not contain an AtPNP-A representative sequence, in an unpublished experiment using the 24K chip, AtPNP-A was up-regulated over 50 fold in over-expressing WRKY 70 lines and was amongst the top 20 genes that are up-regulated in this study (Figure 4; Personal communication: Gunter Brader, Faculty of Biosciences, University of Helsinki). Additionally, a strong induction in expression of the 25 correlated genes is also observed in this experiment providing further evidence that indicates that WRKY 70 may positively regulate the expression of AtPNP-A and the expression correlated genes.

thumbnailFigure 4. Expression profile of AtPNP-A and selected correlated genes in selected A. thaliana mutants. The expression profiles of AtPNP-A and the correlated genes were examined in a number of SA/SAR related mutants. The expression of AtPNP-A and selected genes is greatly elevated in WRKY 70 over-expresser lines and in mutants with elevated SA levels such as cpr5 and mpk4. Conversely, in the SA deficient mutant nahG, expression of the selected genes is markedly reduced. Error bars represent standard errors of the mean.

In the cpr5 (constitutive expresser of pathogenesis related genes) and mpk4 mutants that have elevated levels of SA and display constitutive expression of PR genes [70,71], the expression of AtPNP-A and the correlated genes was markedly elevated (Figure 4). It is of particular interest that the four listed mutants that displayed the largest increase in the expression of AtPNP-A in Genevestigator are all cpr5 mutants, being cpr5/scv1, cpr5/npr1, cpr5, cpr5/npr1/svi1 (range +7.55 to +6.21log2) (data not shown).

Conversely, the nahG mutant that is defective in SA production and signalling, is the only experiment presented in this study that documents a large reduction in the expression of AtPNP-A (-4.6log2) and the correlated genes (Figure 4). This experiment was performed in senescing leaves [72] to identify SA-dependent global gene expression patterns during developmental senescence since SA has previously been shown to be required for expression of some senescence-induced genes [73]. In this study, and in the ATGE developmental series of A. thaliana microarray experiments [74] transcript levels of AtPNP-A were elevated approximately 2.8 fold in senescing leaves compared to adult green leaves (data not shown) indicating that AtPNP-A is a senescence enhanced gene. Furthermore, since transcript levels of AtPNP-A were reduced beyond detection limits in senescing leaves in the nahG mutant this induction appears to be SA-dependent. This pattern is completely consistent with other results since premature senescence, including leaf yellowing and necrosis can be induced by biotic and aboitic stresses that stimulate SA production, including ozone [75] and UV-B [76] which also induce large increases in expression of AtPNP-A. Thus, there is evidence documenting induction of AtPNP-A expression in SA-mediated natural developmental and stress activated processes which both culminate in cell death indicating that AtPNP-A may be involved in these processes. The mutant analyses further enforce that the transcriptional regulation of AtPNP-A and the correlated genes is largely controlled by SA.

The role of a TGA TF in PR gene expression

Additional evidence for co-regulation of AtPNP-A with SAR annotated genes is provided by the observation that expression of another SAR annotated TF, the TGA3 bZIP TF (At1g22070), is correlated with that of AtPNP-A (r = 0.49). Upon induction of SAR, NIMIN1 and NPR1 form a ternary complex with TGA factors in the nucleus which enhances their binding to the positive regulatory as-1 (activator sequence 1) or as-1-like (TGACG) cis-elements that are present in the promoter of several plant genes activated during defense, including A. thaliana PR-1 [77-79]. A manual inspection of the AtPNP-A promoter identified two occurrences of the TGACG motif in close proximity to the TSS (start at -94 and +24). The correlation in expression between AtPNP-A, NPR1, NIMIN1 and TGA3 (Table 1, and [see Additional file 2]) along with the identification of TGA3 cis-elements in the promoter of AtPNP-A is strong evidence that these two factors contribute to the regulation of AtPNP-A expression.

AtPNP-A as PR Protein

Based on the above results we suggest that AtPNP-A could be classified as a PR protein since it possesses many of the criteria that define this class of proteins. The name "pathogenesis related protein" is a collective term that encompasses all proteins that are present at almost undetectable levels in healthy tissue but are induced at the protein level following pathogen infection. The classification of these proteins is based on their pathogen inducible expression rather than defined functional roles in defence [43]. This point is brought into focus when considering PR-1, which is the quintessential marker of the SAR response yet its biological role is largely unknown [80]. Although AtPNP-A is yet to be proven to be induced at the protein level in response to pathogens, elevated protein levels have been shown as a result of abiotic stresses [15]. In addition, transcription of AtPNP-A is low under control conditions but strongly induced in response to biotic and abiotic stresses and the protein has been identified and isolated from the A. thaliana apoplast together PR-1, PR-2 and PR-5 proteins [16]. AtPNP-A has other features characteristic of PR proteins including an N-terminal signal peptide [43] that directs the molecule into the extracellular space. Further, induction of AtPNP-A at the transcript level appears to occur independent of de novo protein synthesis characteristic of genes encoding secreted proteins [28]. The evolutionary history of AtPNP-A suggests that PNPs, like the related expansins, derived from ancestral family-45 endoglucanases that have lost their hydrolytic activity and have sub-functionalized into extracellular, systemically mobile signalling molecules [9].

Future directions

In order to determine the physiological role of AtPNP-A in A. thaliana a T-DNA insertion mutant, that is available from SALK, could be used. Phenotyping this mutant in response to SA-inducing abiotic and biotic stresses as well as during developmental senescence, will help characterise specific physiological processes in which AtPNP-A is involved. If such a mutant demonstrated a compromised SAR response, it would greatly strengthen the claim that AtPNP-A is indeed involved in the SAR response pathway. Additionally, it will be interesting to look at the expression of the correlated genes in an AtPNP-A mutant in response to SAR inducing conditions since this may enable us to determine a role for AtPNP-A in the context of a SAR response pathway.

Conclusion

AtPNP-A is an annotated "signal" molecule that is secreted into the apoplastic space and has been implicated with a role in the control of ion and solute movements in plants (Figure 1). The expression of AtPNP-A is significantly correlated with that of genes involved in the SAR defence response pathway in response to various biotic and abiotic stimuli and in mutant studies. The expression of AtPNP-A is correlated with the ICS-1 gene that is involved in SA biosynthesis and NPR1 and NIMIN1, which are key positive regulators of the SAR response, and two annotated SAR TFs TGA 3 and WRKY 70. Additionally, like the PR genes, the promoter of AtPNP-A contains as-1 and W-box cis-elements that correspond to binding sites for the TGA 3 and WRKY TFs. Further, over expressing WRKY 70 lines have been shown to cause a greater than 50-fold increase in the expression of AtPNP-A which is consistent with this TF being a positive regulator of AtPNP-A transcription. The induced expression of AtPNP-A by SAR elicitors and its secretion into the apoplast is similar to that of PR proteins and strongly implicates a role for AtPNP-A in plant SAR defence responses which may involve the modification of cellular ion and water homeostasis during stress responses.

Methods

Identification of Correlated Genes

We downloaded (01/07/2005) A. thaliana gene expression levels for 1877 experiments from the NASCArrays database [81], using the bulk data download option. Perl scripts were used to calculate non-parametric correlation coefficients (Spearman's rho) between the expression of AtPNP-A (At2g18660) and each of the approximately 22,000 genes represented on the Affymetrix array that was used to generate this data set. We ranked genes according to the correlation coefficient and reported genes that were most positively correlated with At2g18660. The p-values were calculated using the bivariate normal distribution, with p representing the probability of observing an equal or larger positive or negative correlation by chance.

Functional Classification and expression analysis of Correlated Genes

To characterise the correlated genes the web-based 'FatiGO+' program [22] was used to search for differential distributions of gene ontology (GO) and biological terms within this list. The search was conducted using AtPNP-A (At2g18660) and the 25 most positively correlated genes in Table 1 (list 1 = 26 genes). This list was compared to a reference gene list that contained the remaining genes in the entire A. thaliana genome (list 2 = 26147 genes). Statistical significance was determined using the Family Wise Error Rate (FWER) to calculate the adjusted p-value.

The expression profiles of AtPNP-A and the positively correlated genes (Table 1) were initially examined using Affymetrix public microarray data in the gene response viewer tool (GRV) in Genevestigator [68]). The analysis was performed using the ATH1: 22K array chip type and included all of the available 2507 chip sources. For better temporal and spatial response resolution we obtained normalised microarray data from the following sites: NASCArrays, Ozone-26 (reference ID); P. infestans-123; UV-B stress-144; Potassium starvation-105; BHT-392. TAIR (ATGenExpress): Salicylic acid-ME00364; E.orontii-ME00354; Salt stress-ME00328; Osmotic stress-ME00327; Cold acclimation-ME00369; Cyclohexamide-ME00361. GEO (NCBI): E. cichoracearum-GSE431.

In order to further reveal the relationship of AtPNP-A expression with that of key genes involved in the SAR response, the mutant surveyor in Genevestigator was used to compare gene expression in different types of defence related A. thaliana mutants. The genes investigated in this study included AtPNP-A, the 25 correlated genes, and WRKY 70 and WRKY 46. Normalised array data from the mutant experiments were obtained from: TAIR-ME00373 for cpr5/npr1 mutants; NASCarray-52 for the nahG mutant and array express (EBI) for the mpk4 mutant (E-MEXP-173). For the WRKY 70 over expresser, data was obtained through personal communication with Gunter Brader, Faculty of Biosciences, University of Helsinki.

Promoter analysis

The web-based Athena [59] and POBO [58]) applications were used to analyse the promoters (-1 kb upstream of the predicted TSS) of AtPNP- and the 25 top correlated genes.

In POBO [58], the 1 kb promoter sequences were uploaded and the analysis was run against A. thaliana background (clean) searching for the TTGAC W-box core motif using the default settings (number of sequences to pick-out = 50, number of samples to generate = 1000, sequence length = 1000 bps). A two-tailed p-value was calculated in the linked online GraphPad web-site using the generated t-value and degrees of freedom to determine the statistical differences between input sequences and background.

In Athena, the analysis was performed with the visualisation tool using the 26 correlated genes with settings of 1000 bp upstream and do not cut off at adjacent genes. The statistical significance of over-represented TF binding sites is automatically calculated using a hypergeometric probability model to calculate the p-value. A Bonferroni correction was automatically used in Athena to account for multiple hypotheses testing (up to 105 different TF binding sites) and determined that the p-value threshold for significant enrichment was < 10-4.

Abbreviations

NP – natriuretic peptide; PNP – plant natriuretic peptide; SAR – systemic acquired resistance; cGMP – guanosine 3',5'-cyclic monophosphate; TF – transcription factor; TSS – translation start site.

Authors' contributions

The project was conceived by CG and VB. StM, RB, LD and ShM have extracted and interpreted the data. The manuscript was written by StM and CG. All authors read and approved the final manuscript.

Acknowledgements

This project was supported by the South African National Research Foundation and an Ernest Oppenheimer Memorial Trust fellowship to C.G.

References

  1. Takei Y: Does the natriuretic peptide system exist throughout the animal and plant kingdom?

    Comp Biochem Physiol B: Biochem Mol Biol 2001, 129:559-573. PubMed Abstract | Publisher Full Text OpenURL

  2. Vesely DL, Giordano AT: Atrial natriuretic peptide hormonal system in plants.

    Biochem Biophys Res Commun 1991, 179:695-700. PubMed Abstract | Publisher Full Text OpenURL

  3. Vesely DL, Gower WR, Giordano AT: Atrial natriuretic peptides are present throughout the plant kingdom and enhance solute flow in plants.

    Am J Physiol 1993, 265:E465-E477. PubMed Abstract | Publisher Full Text OpenURL

  4. Gehring CA, Khalid KM, Toop T, Donald JA: Rat natriuretic peptide binds specifically to plant membranes and induces stomatal opening.

    Biochem Biophys Res Commun 1996, 228:739-744. PubMed Abstract | Publisher Full Text OpenURL

  5. Pharmawati M, Billington T, Gehring CA: Stomatal guard cell responses to kinetin and natriuretic peptides are cGMP-dependent.

    Cell Mol Life Sci 1998, 54:272-276. PubMed Abstract | Publisher Full Text OpenURL

  6. Pharmawati M, Gehring CA, Irving HR: An immunoaffinity purified plant natriuretic peptide analogue modulates cGMP level in the Zea mays root stele.

    Plant Sci 1998, 137:107-115. Publisher Full Text OpenURL

  7. Pharmawati M, Maryani MM, Nikolakopoulos T, Gehring CA, Irving HR: Cyclic GMP modulates stomatal opening induced by natriuretic peptides and immunoreactive analogues.

    Plant Physiol Biochem 2001, 39:385-394. Publisher Full Text OpenURL

  8. Suwastika IN, Toop T, Irving HR, Gehring CA: In situ and in vitro binding of natriuretic peptide hormones in Tradescantia multiflora.

    Plant Biol 2000, 2:1-3. Publisher Full Text OpenURL

  9. Ludidi NN, Heazlewood JL, Seoighe C, Irving HR, Gehring CA: Expansin-like molecules: novel functions derived from common domains.

    J Mol Evol 2002, 54:587-594. PubMed Abstract | Publisher Full Text OpenURL

  10. Nembaware V, Seoighe C, Sayed M, Gehring C: A plant natriuretic peptide-like gene in the bacterial pathogen Xanthomonas axonopodis may induce hyper-hydration in the plant host: a hypothesis of molecular mimicry.

    Bmc Evol Biol 2004, 4:10. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  11. McQueen-Mason S, Cosgrove DJ: Expansin mode of action on cell walls: analysis of wall hydrolysis, stress relaxation, and binding.

    Plant Physiol 1995, 107:87-100. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Kende H, Bradford K, Brummell D, Cho HT, Cosgrove D, Fleming A, Gehring C, Lee Y, McQueen-Mason S, Rose J, Voesenek LA: Nomenclature for members of the expansin superfamily of genes and proteins.

    Plant Mol Biol 2004, 55:311-314. PubMed Abstract | Publisher Full Text OpenURL

  13. Morse M, Pironcheva G, Gehring C: AtPNP-A is a systemically mobile natriuretic peptide immunoanalogue with a role in Arabidopsis thaliana cell volume regulation.

    Febs Lett 2004, 556:99-103. PubMed Abstract | Publisher Full Text OpenURL

  14. Maryani MM, Morse MV, Bradley G, Irving HR, Cahill DM, Gehring CA: In situ localization associates biologically active plant natriuretic peptide immuno-analogues with conductive tissue and stomata.

    J Exp Bot 2003, 54:1553-1564. PubMed Abstract | Publisher Full Text OpenURL

  15. Rafudeen S, Gxaba G, Makgoke G, Bradley G, Pironcheva G, Raitt L, Irving H, Gehring C: A role for plant natriuretic peptide immuno-analogues in NaCl- and drought-stress responses.

    Physiol Plantarum 2003, 119:554-562. Publisher Full Text OpenURL

  16. Boudart G, Jamet E, Rossignol M, Lafitte C, Borderies G, Jauneau A, Esquerre-Tugaye MT, Pont-Lezica R: Cell wall proteins in apoplastic fluids of Arabidopsis thaliana rosettes: Identification by mass spectrometry and bioinformatics.

    Proteomics 2004, 5:212-221. Publisher Full Text OpenURL

  17. Maryani MM, Bradley G, Cahill DM, Gehring CA: Natriuretic peptides and immunoreactants modify osmoticum-dependent volume changes in Solanum tuberosum L. mesophyll cell protoplasts.

    Plant Sci 2001, 161:443-452. Publisher Full Text OpenURL

  18. Wang Y, Gehring C, Cahill D, Irving H: Plant natriuretic peptide active site determination and effects on cGMP and cell volume regulation.

    Funct Plant Biol 2007, 34:645-653. Publisher Full Text OpenURL

  19. Pharmawati M, Shabala SN, Newman IA, Gehring CA: Natriuretic peptides and cGMP modulate K+, Na+, and H+ fluxes in Zea mays roots.

    Mol Cell Biol Res Commun 1999, 2:53-57. PubMed Abstract | Publisher Full Text OpenURL

  20. Ludidi N, Morse M, Sayed M, Wherrett T, Shabala S, Gehring C: A recombinant plant natriuretic peptide causes rapid and spatially differentiated K+, Na+ and H+ flux changes in Arabidopsis thaliana roots.

    Plant Cell Physiol 2004, 45:1093-1098. PubMed Abstract | Publisher Full Text OpenURL

  21. Gehring CA, Irving HR: Natriuretic peptides - a class of heterologous molecules in plants.

    Int J Biochem Cell Biol 2003, 35:1318-1322. PubMed Abstract | Publisher Full Text OpenURL

  22. Al-Shahrour F, Diaz-Uriarte R, Dopazo J: FatiGO: a web tool for finding significant associations of gene ontology terms with groups of genes.

    Bioinf 2004, 20:578-580. Publisher Full Text OpenURL

  23. Al-Shahrour F, Minguez P, Vaquerizas JM, Conde L, Dopazo J: BABELOMICS: a suite of web tools for functional annotation and analysis of groups of genes in high-throughput experiments.

    Nucleic Acids Res 2005, 33(Web Server issue):W460-464. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Al-Shahrour F, Minguez P, Tarraga J, Medina I, Alloza E, Montaner D, Dopazo J: FatiGO +: a functional profiling tool for genomic data. Integration of functional annotation, regulatory motifs and interaction data with microarray experiments.

    Nucleic Acids Res 2007, 35(Web Server issue):W91-96. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Durrant WE, Dong X: Systemic acquired resistance.

    Annu Rev Phytopathol 2004, 42:185-209. PubMed Abstract | Publisher Full Text OpenURL

  26. Weigel RR, Pfitzner UM, Gatz C: Interaction of NIMIN1 with NPR1 modulates PR gene expression in Arabidopsis.

    Plant Cell 2005, 17:1279-1291. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Delaney TP, Friedrich L, Ryals JA: Arabidopsis signal transduction mutant defective in chemically and biologically induced disease resistance.

    Proc Natl Acad Sci U S A 1995, 92:6602-6606. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Wang D, Amornsiripanitch N, Dong X: A genomic approach to identify regulatory nodes in the transcriptional network of systemic acquired resistance in plants.

    PLoS Pathog 2006, 2:e123. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Wildermuth MC, Dewdney J, Wu G, Ausubel FM: Isochorismate synthase is required to synthesize salicylic acid for plant defence.

    Nature 2001, 414:562-565. PubMed Abstract | Publisher Full Text OpenURL

  30. McDowell JM, Dangl JL: Signal transduction in the plant immune response.

    Trends Biochem Sci 2000, 25:79-82. PubMed Abstract | Publisher Full Text OpenURL

  31. Yalpani N, Enyedi AJ, Leon J, Raskin I: Ultraviolet light and ozone stimulate accumulation of salicylic acid, pathogenesis-related proteins and virus resistance in tobacco.

    Planta 1994, 193:372-376. Publisher Full Text OpenURL

  32. Yang PZ, Chen CH, Wang ZP, Fan BF, Chen ZX: A pathogen- and salicylic acid-induced WRKY DNA-binding activity recognizes the elicitor response element of the tobacco class I chitinase gene promoter.

    Plant J 1999, 18:141-149. Publisher Full Text OpenURL

  33. Mackerness SAH, Surplus SL, Blake P, John CF, Buchanan-Wollaston V, Jordan BR, Thomas B: Ultraviolet-B-induced stress and changes in gene expression in Arabidopsis thaliana: role of signalling pathways controlled by jasmonic acid, ethylene and reactive oxygen species.

    Plant Cell Environ 1999, 22:1413-1423. Publisher Full Text OpenURL

  34. Mackerness SAH, Jordan BR, Thomas B: Reactive oxygen species in the regulation of photosynthetic genes by ultraviolet-B radiation (UV-B: 280-320 nm) in green and etiolated buds of pea (Pisum sativum L.).

    J Photoch Photobio B 1999, 48:180-188. Publisher Full Text OpenURL

  35. Verslues PE, Bray EA: Role of abscisic acid (ABA) and Arabidopsis thaliana ABA-insensitive loci in low water potential-induced ABA and proline accumulation.

    J Exp Bot 2006, 57:201-212. PubMed Abstract | Publisher Full Text OpenURL

  36. Denby K, Gehring C: Engineering drought and salinity tolerance in plants: lessons from genome-wide expression profiling in Arabidopsis.

    Trends Biotech 2005, 23:547-552. Publisher Full Text OpenURL

  37. Suwastika IN, Gehring CA: Natriuretic peptide hormones promote radial water movements from the xylem of Tradescantia shoots.

    Cell Mol Life Sci 1998, 54:1161-1167. Publisher Full Text OpenURL

  38. Mahajan S, Tuteja N: Cold, salinity and drought stresses: an overview.

    Arch Biochem Biophys 2005, 444:139-158. PubMed Abstract | Publisher Full Text OpenURL

  39. Rus A, Lee BH, Munoz-Mayor A, Sharkhuu A, Miura K, Zhu JK, Bressan RA, Hasegawa PM: AtHKT1 facilitates Na+ homeostasis and K+ nutrition in planta.

    Plant Physiol 2004, 136:2500-2511. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  40. Serrano R, Rodriguez-Navarro A: Ion homeostasis during salt stress in plants.

    Curr Opin Cell Biol 2001, 13:399-404. PubMed Abstract | Publisher Full Text OpenURL

  41. Dat JF, Foyer CH, Scott IM: Changes in salicylic acid and antioxidants during induced thermotolerance in mustard seedlings.

    Plant Physiol 1998, 118:1455-1461. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  42. Kang HM, Saltveit ME: Chilling tolerance of maize, cucumber and rice seedling leaves and roots are differentially affected by salicylic acid.

    Physiol Plantarum 2002, 115:571-576. Publisher Full Text OpenURL

  43. van Loon LC, Rep M, Pieterse CMJ: Significance of inducible defense-related proteins in infected plants.

    Annu Rev Phytopathol 2006, 44:135–162. Publisher Full Text OpenURL

  44. Durner J, Wendehenne D, Klessig DF: Defense gene induction in tobacco by nitric oxide, cyclic GMP, and cyclic ADP-ribose.

    Proc Natl Acad Sci U S A 1998, 95:10328-10333. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  45. Donaldson L, Ludidi N, Knight MR, Gehring C, Denby K: Salt and osmotic stress cause rapid increases in Arabidopsis thaliana cGMP levels.

    Febs Lett 2004, 569:317-320. PubMed Abstract | Publisher Full Text OpenURL

  46. Yoshioka K, Moeder W, Kang HG, Kachroo P, Masmoudi K, Berkowitz G, Klessig DF: The chimeric Arabidopsis CYCLIC NUCLEOTIDE-GATED ION CHANNEL11/12 activates multiple pathogen resistance responses.

    Plant Cell 2006, 18:747-763. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  47. Bouche N, Yellin A, Snedden WA, Fromm H: Plant-specific calmodulin-binding proteins.

    Annu Rev Plant Biol 2005, 56:435-466. PubMed Abstract | Publisher Full Text OpenURL

  48. Chandra S, Stennis M, Low PS: Measurement of Ca2+ fluxes during elicitation of the oxidative burst in aequorin-transformed tobacco cells.

    J Biol Chem 1997, 272:28274-28280. PubMed Abstract | Publisher Full Text OpenURL

  49. Sasabe M, Takeuchi K, Kamoun S, Ichinose Y, Govers F, Toyoda K, Shiraishi T, Yamada T: Independent pathways leading to apoptotic cell death, oxidative burst and defense gene expression in response to elicitin in tobacco cell suspension culture.

    Eur J Biochem 2000, 267:5005-5013. PubMed Abstract | Publisher Full Text OpenURL

  50. Raz V, Fluhr R: Calcium requirement for ethylene-dependent responses.

    Plant Cell 1992, 4:1123-1130. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  51. Snedden WA, Fromm H: The calcium and calmodulin signaling network in plants.

    New Phytol 2001, 151:35-66. Publisher Full Text OpenURL

  52. Suzuki K, Suzuki N, Ohme-Takagi M, Shinshi H: Immediate early induction of mRNAs for ethylene-responsive transcription factors in tobacco leaf strips after cutting.

    Plant J 1998, 15:657-665. Publisher Full Text OpenURL

  53. Knetsch M, Wang M, Snaar-Jagalska BE, Heimovaara-Dijkstra S: Abscisic acid induces mitogen-activated protein kinase activation in barley aleurone protoplasts.

    Plant Cell 1996, 8:1061-1067. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  54. Abel S, Theologis A: Early genes and auxin action.

    Plant Physiol 1996, 111:9-17. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  55. Nishiuchi T, Suzuki K, Kitajima S, Sato F, Shinshi H: Wounding activates immediate early transcription of genes for ERFs in tobacco plants.

    Plant Mol Biol 2002, 49:473-482. PubMed Abstract | Publisher Full Text OpenURL

  56. Koshiba T, Ballas N, Wong LM, Theologis A: Transcriptional regulation of PS-IAA4/5 and PS-IAA6 early gene expression by indoleacetic acid and protein synthesis inhibitors in pea (Pisum sativum).

    J Mol Biol 1995, 253:396-413. PubMed Abstract | Publisher Full Text OpenURL

  57. Wang D, Weaver ND, Kesarwani M, Dong X: Induction of protein secretory pathway is required for systemic acquired resistance.

    Science 2005, 308:1036-1040. PubMed Abstract | Publisher Full Text OpenURL

  58. Kankainen M, Holm L: POBO, transcription factor binding site verification with bootstrapping.

    Nucleic Acids Res 2004, 32(Web Server issue):W222-229. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  59. O'Connor TR, Dyreson C, Wyrick JJ: Athena: a resource for rapid visualization and systematic analysis of Arabidopsis promoter sequences.

    Bioinf 2005, 21:4411-4413. Publisher Full Text OpenURL

  60. Maleck K, Levine A, Eulgem T, Morgan A, Schmidl J, Lawton KA, Dangl JL, Dietrich RA: An Arabidopsis promoter element shared among genes co-regulated during systemic acquired disease resistance.

    Nat Genet 2000, 26:403-410. PubMed Abstract | Publisher Full Text OpenURL

  61. Eulgem T, Rushton PJ, Robatzek S, Somssich IE: The WRKY superfamily of plant transcription factors.

    Trends Plant Sci 2000, 5:199-206. PubMed Abstract | Publisher Full Text OpenURL

  62. Huang T, Nicodemus J, Zarka DG, Thomashow MF, Wisniewski M, Duman JG: Expression of an insect (Dendroides canadensis) antifreeze protein in Arabidopsis thaliana results in a decrease in plant freezing temperature.

    Plant Mol Biol 2002, 50:333-344. PubMed Abstract | Publisher Full Text OpenURL

  63. Hara K, Yagi M, Kusano T, Sano H: Rapid systemic accumulation of transcripts encoding a tobacco WRKY transcription factor upon wounding.

    Mol Gen Genet 2000, 263:30-37. PubMed Abstract | Publisher Full Text OpenURL

  64. Rizhsky L, Davletova S, Liang H, Mittler R: The zinc finger protein Zat12 is required for cytosolic ascorbate peroxidase 1 expression during oxidative stress in Arabidopsis.

    J Biol Chem 2004, 279:11736-11743. PubMed Abstract | Publisher Full Text OpenURL

  65. Pnueli L, Hallak-Herr E, Rozenberg M, Cohen M, Goloubinoff P, Kaplan A, Mittler R: Molecular and biochemical mechanisms associated with dormancy and drought tolerance in the desert legume Retama raetam.

    Plant J 2002, 31:319-330. PubMed Abstract | Publisher Full Text OpenURL

  66. Seki M, Ishida J, Narusaka M, Fujita M, Nanjo T, Umezawa T, Kamiya A, Nakajima M, Enju A, Sakurai T, Satou M, Akiyama K, Yamaguchi-Shinozaki K, Carninci P, Kawai J, Hayashizaki Y, Shinozaki K: Monitoring the expression pattern of around 7,000 Arabidopsis genes under ABA treatments using a full-length cDNA microarray.

    Funct Integr Genomics 2002, 2:282-291. PubMed Abstract | Publisher Full Text OpenURL

  67. Rizhsky L, Liang H, Shuman J, Shulaev V, Davletova S, Mittler R: When defense pathways collide. The response of Arabidopsis to a combination of drought and heat stress.

    Plant Physiol 2004, 134:1683-1696. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  68. Zimmermann P, Hirsch-Hoffmann M, Hennig L, Gruissem W: GENEVESTIGATOR. Arabidopsis microarray database and analysis toolbox.

    Plant Physiol 2004, 136:2621-2632. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  69. Li J, Brader G, Palva ET: The WRKY70 transcription factor: a node of convergence for jasmonate-mediated and salicylate-mediated signals in plant defense.

    Plant Cell 2004, 16:319-331. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  70. Bowling SA, Clarke JD, Liu Y, Klessig DF, Dong X: The cpr5 mutant of Arabidopsis expresses both NPR1-dependent and NPR1-independent resistance.

    Plant Cell 1997/10/27 edition. 1997, 9:1573-1584. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  71. Petersen M, Brodersen P, Naested H, Andreasson E, Lindhart U, Johansen B, Nielsen HB, Lacy M, Austin MJ, Parker JE, Sharma SB, Klessig DF, Martienssen R, Mattsson O, Jensen AB, Mundy J: Arabidopsis map kinase 4 negatively regulates systemic acquired resistance.

    Cell 2000, 103:1111-1120. PubMed Abstract | Publisher Full Text OpenURL

  72. Buchanan-Wollaston V, Page T, Harrison E, Breeze E, Lim PO, Nam HG, Lin JF, Wu SH, Swidzinski J, Ishizaki K, Leaver CJ: Comparative transcriptome analysis reveals significant differences in gene expression and signalling pathways between developmental and dark/starvation-induced senescence in Arabidopsis.

    Plant J 2005, 42:567-585. PubMed Abstract | Publisher Full Text OpenURL

  73. Morris K, MacKerness SA, Page T, John CF, Murphy AM, Carr JP, Buchanan-Wollaston V: Salicylic acid has a role in regulating gene expression during leaf senescence.

    Plant J 2000, 23:677-685. PubMed Abstract | Publisher Full Text OpenURL

  74. Schmid M, Davison T, Henz S, Pape U, Demar M, Vingron M, Schölkopf B, Weigel D, Lohmann J: A gene expression map of Arabidopsis thaliana development.

    Nat Genet 2005, 37:501–506.. Publisher Full Text OpenURL

  75. Torsethaugen G, Pitcher LH, Zilinskas BA, Pell EJ: Overproduction of ascorbate peroxidase in the tobacco chloroplast does not provide protection against ozone.

    Plant Physiol 1997, 114:529-537. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  76. Jansen MAK, Gaba V, Greenberg BM: Higher plants and UV-B radiation: balancing damage, repair and acclimation.

    Trends Plant Sci 1998, 3:131–135. OpenURL

  77. Zhang Y, Fan W, Kinkema M, Li X, Dong X: Interaction of NPR1 with basic leucine zipper protein transcription factors that bind sequences required for salicylic acid induction of the PR-1 gene.

    Proc Natl Acad Sci U S A 1999, 96:6523-6528. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  78. Despres C, Chubak C, Rochon A, Clark R, Bethune T, Desveaux D, Fobert PR: The Arabidopsis NPR1 disease resistance protein is a novel cofactor that confers redox regulation of DNA binding activity to the basic domain/leucine zipper transcription factor TGA1.

    Plant Cell 2003, 15:2181-2191. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  79. Johnson C, Boden E, Arias J: Salicylic acid and NPR1 induce the recruitment of trans-activating TGA factors to a defense gene promoter in Arabidopsis.

    Plant Cell 2003, 15:1846-1858. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  80. Van Loon LC, Van Strien EA: The families of pathogenesis-related proteins, their activities, and comparative analysis of PR-1 type proteins.

    Physiol Mol Plant Pathol 1999, 55:85–97. Publisher Full Text OpenURL

  81. Craigon DJ, James N, Okyere J, Higgins J, Jotham J, May S: NASCArrays: a repository for microarray data generated by NASC's transcriptomics service.

    Nucleic Acids Res 2004, 32(Database issue):D575-577. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL