Common and unique elements of the ABA-regulated transcriptome of Arabidopsis guard cells

Background In the presence of drought and other desiccating stresses, plants synthesize and redistribute the phytohormone abscisic acid (ABA). ABA promotes plant water conservation by acting on specialized cells in the leaf epidermis, guard cells, which border and regulate the apertures of stomatal pores through which transpirational water loss occurs. Following ABA exposure, solute uptake into guard cells is rapidly inhibited and solute loss is promoted, resulting in inhibition of stomatal opening and promotion of stomatal closure, with consequent plant water conservation. There is a wealth of information on the guard cell signaling mechanisms underlying these rapid ABA responses. To investigate ABA regulation of gene expression in guard cells in a systematic genome-wide manner, we analyzed data from global transcriptomes of guard cells generated with Affymetrix ATH1 microarrays, and compared these results to ABA regulation of gene expression in leaves and other tissues. Results The 1173 ABA-regulated genes of guard cells identified by our study share significant overlap with ABA-regulated genes of other tissues, and are associated with well-defined ABA-related promoter motifs such as ABREs and DREs. However, we also computationally identified a unique cis-acting motif, GTCGG, associated with ABA-induction of gene expression specifically in guard cells. In addition, approximately 300 genes showing ABA-regulation unique to this cell type were newly uncovered by our study. Within the ABA-regulated gene set of guard cells, we found that many of the genes known to encode ion transporters associated with stomatal opening are down-regulated by ABA, providing one mechanism for long-term maintenance of stomatal closure during drought. We also found examples of both negative and positive feedback in the transcriptional regulation by ABA of known ABA-signaling genes, particularly with regard to the PYR/PYL/RCAR class of soluble ABA receptors and their downstream targets, the type 2C protein phosphatases. Our data also provide evidence for cross-talk at the transcriptional level between ABA and another hormonal inhibitor of stomatal opening, methyl jasmonate. Conclusions Our results engender new insights into the basic cell biology of guard cells, reveal common and unique elements of ABA-regulation of gene expression in guard cells, and set the stage for targeted biotechnological manipulations to improve plant water use efficiency.


Background
Drought stress imposes one of the greatest limitations to crop growth and yield worldwide [1]. Limitations in fresh water availability are predicted to become an increasing problem due to industrialization and continuing global climate change [2][3][4]. One of the central mechanisms in plant drought tolerance is a reduction in stomatal apertures at the leaf surface, with consequent increase in water use efficiency [5][6][7]. Abscisic acid (ABA) is a major plant hormone which inhibits growth and promotes tolerance of abiotic stresses such as drought, salinity and cold [8][9][10]. Under drought conditions, plants synthesize and redistribute ABA which triggers cellular responses in guard cells, specialized cells that flank the stomatal pores, leading to inhibition of stomatal opening and promotion of stomatal closure and thereby reducing plant water loss. Thus, an improved understanding of guard cell responses to ABA is relevant for development of cultivars with improved productivity under drought conditions. Guard cells also have become an advanced model system for understanding plant signal transduction [11][12][13]. Much of the research on guard cell responses has focused on the signaling cascades underlying rapid responses to ABA as well as to other inputs such as light and CO 2 [14]. However, in addition to these posttranscriptional responses, guard cells also respond to ABA at the level of the transcriptome. For example, Li et al. showed in Vicia faba guard cell protoplasts that these cells respond with an increase in dehydrin transcript levels after just 15 min. of ABA exposure [15], suggesting that some of the cellular signaling processes that are observed could be influenced by transcriptional events. Leonhardt et al. published an expression analysis of the ABA-related transcriptome of guard cell protoplasts [16], utilizing the Affymetrix AG GeneChip, which represents~8200 Arabidopsis genes. In those experiments, whole plants were treated with 100 μM ABA for four hours, and transcriptional inhibitors were applied during the protoplasting process, with the idea of stabilizing transcript levels. The authors demonstrated ABA-regulation of approximately 150 genes in guard cell protoplasts. To extend the investigation of gene expression in guard cells, Yang et al. isolated guard cell promoter candidates based on 4 Arabidopsis guard cell microarrays from 4 different treatments using wholegenome Affymetrix ATH1 chips [17], which profilẽ 24,000 genes. To date, however, there has been no systematic replicated study that investigates ABA regulation of gene expression in guard cells in a genome-wide manner. A comprehensive analysis of the global ABA-related transcriptome of guard cells and its relationship to known guard-cell signaling cascades promises to engender new insights into the basic cell biology of guard cells. In addition, the wealth of microarray data regarding ABAregulated transcriptomes of other tissues, including whole plants [18,19], seedlings [20][21][22][23][24][25][26], and seeds [27], provides an opportunity to assess the extent to which the specialized guard cells employ unique vs. common promoter motifs and unique vs. common gene targets in ABA regulation of gene expression. Given growing evidence of hormonal cross-talk in physiological responses [10], another interesting question is whether or not ABA-regulated transcripts of guard cells are regulated by other hormones in other plant tissues. Finally, manipulation of transcriptional regulation of key genes is one of the most common approaches utilized in applied research aimed at improving plant stress tolerance [28][29][30], so genome-wide analysis of the Arabidopsis guard cell transcriptome may yield candidate genes for translational research. In this study, we performed a systematic microarray analysis of gene expression profile data from global transcriptomes of guard cells generated with Affymetrix ATH1 chips to explore the above topics.

Identification of ABA-responsive genes in guard cells and leaves
We and others [16] have observed that when guard cell protoplasts were used as the source of RNA, stressinduced genes could be induced by protoplasting, even without ABA treatment. Therefore, we employed a protocol [31] whereby isolated epidermal peels with guard cells as the only intact cell type ( Figure 1A and Figure  1B) were used as the source of guard cell RNA for microarray analysis.
For each tissue by treatment combination (guard cell vs. leaf; ABA vs. ethanol as solvent control), three independent biological samples were assessed using Affymetrix whole genome ATH1 chips [31]. For the present study, to reliably identify ABA-responsive genes, we adopted two different methods -a Boolean method [31] and a standard linear model method [32] -and focused on the set of ABA-regulated genes jointly confirmed by both approaches (see Methods for details). Through this overlapping analysis we identified 1173 ABA-responsive genes in Col (Arabidopsis Col-0) guard cells (696 ABAinduced genes and 477 ABA-repressed genes) and 1037 ABA-responsive genes in Col leaves (596 ABA-induced genes and 441 ABA-repressed genes). The ABA-regulated genes in guard cells and leaves are listed in Additional file 1 and Additional file 2, respectively. As shown in Figure 1C and 1D, among all the ABA-induced genes, only 177 genes are ABA-induced in both guard cells and leaves, while the rest show differential regulation with respect to tissue type. The overlap of guard cell ABArepressed genes and leaf ABA-repressed genes is even smaller with only 90 genes common to both guard cells and leaves. There are only 4 genes that are oppositely regulated by ABA in guard cells vs. leaves: AT5G20840, EXL3 (AT5G51550), and CIL (AT4G25990) are induced by ABA in guard cells but repressed by ABA in leaves, and KCS5 (AT1G25450) is induced by ABA in leaves but repressed by ABA in guard cells.

Validation of ABA-responsive genes from microarray analysis by Q-PCR
To evaluate the reliability of ABA-regulation derived from the microarray data and validate identified ABAresponsive genes, we performed Q-PCR experiments on 3 genes identified as not responsive to ABA and 14 genes that showed differential expression in response to ABA in either guard cells, leaves, or both tissues by our microarray analysis ( Figure 2). The three negative control genes were Actin (Actin 2: AT3G18780 and Actin 8: AT1G49240), a b-tubulin like protein TUB9 (AT4G20890), and a putative nematode resistant protein Arabidopsis ortholog of sugar beet HS1 PRO-1 2 (HSPRO2, AT2G40000). Among the ABA-responsive genes, we chose two genes that exhibited ABA-upregulation in both guard cells and leaves in our microarray analysis: ABI2 (AT5G57050), which encodes one of several PP2C-type protein phosphatases that negatively regulate many aspects of ABA signaling, including ABA-induced stomatal closure [33][34][35], and the non-specific lipid transfer protein 3 (LTP3, AT5G59320), which is a droughtresponsive gene [36]. In addition, we analyzed two genes that showed ABA-repression in both guard cell and leaf microarrays: a heat shock protein (AT4G21870) and a somatic embryogenesis receptor-like kinase 2 (SERK2, AT1G34210). For the category of ABA-induction in guard cell but not leaf microarrays, five genes were chosen for Q-PCR evaluation: LTI78/RD29A (AT5G52310), identified as a cold and salt stress-induced gene in a leaf transcription study [37,38], two genes encoding phospholipases known to be involved in ABA signaling, namely a phospholipase C gene, PI-PLC (AT3G55940) and a PLDδ homolog (AT4G35790), as well as GBF2 (AT4G01120), and a protein kinase (AT1G77280). For genes that showed ABA-repression in guard cell but not leaf microarrays, we performed Q-PCR assessment of a vacuolar anion channel gene ATCLC-A (AT5G40890) and the ion channel gene KAT1 (AT5G46240), which mediates K + uptake during stomatal opening. For genes that were ABA-responsive only in the leaf microarrays, we tested ECERIFERUM 1 (CER1, AT1G02205) and Ca 2+ -binding protein 1 (ATCP1, AT5G49480) for ABA-induced genes and WALL-ASSOCIATED KINASE 2 (WAK2, AT1G21270) for ABA-repressed genes. agreement between the microarray analysis and the Q-PCR experiments is comparable or better than that observed in other microarray studies [39]. The discrepancy most probably stems from the different techniques employed to measure gene expression. Several studies have noted that microarray analysis is less sensitive than Q-PCR and can underestimate the true extent of transcript accumulation [40,41], so observing this result for a minority of the genes (2 out of 17) is not surprising. cis-acting regulatory elements in ABA-responsive genes Several cis-regulatory elements responsible for mediating ABA-, drought-, and cold-induced gene expression have been identified by biochemical and molecular techniques [42]. We were interested to determine whether these known motifs would be enriched in our gene sets, which would support our identification of these genes as ABA-regulated. In addition, we were interested to assess whether ABA-regulated genes of guard cells and leaves would exhibit differential usage of these promoter motifs, or differential positioning of the motifs within the promoter region. The G-box-containing ABA-responsive elements (ABREs) and coupling element 3 (CE3) are classical cisacting elements in ABA-responsive gene expression [43][44][45][46]. In addition, the core motif of the dehydrationresponsive element/C-repeat (DRE/CRT) is a cis-acting element found in many drought-, high-salt-and coldresponsive genes in Arabidopsis and rice, and provides the binding site for DRE-related transcription factors such as DRE-Binding proteins (DREBs) [46,47]. MYB and MYC binding motifs are also found in the promoters of drought-inducible genes. Proteins of the MYB and MYC families function in mediating drought-and ABA-regulated gene expression [33,48] and some members of these families such as AtMYC2 and AtMYB2 are synthesized following accumulation of endogenous ABA [49]. Other related cis-acting elements include the binding site of C-repeat factors (CBFs) and the low temperature responsive element (LTRE) [46].
We performed a statistical analysis (see Methods) of the promoter sequences (defined as the 1000-bp upstream regions) of ABA-responsive genes in guard cells and leaves, to check whether the known motifs would be enriched in our gene sets. The results of this analysis are shown in Table 1. We can see that ABREs are significantly enriched in guard cell ABA-induced genes. In particular, more than 50% of ABA-induced genes in guard cells have ABREs in their promoter regions, supporting the reliability of our microarray methods and predictions. DREs and LTREs are also highly enriched in guard cell ABA-induced genes. ABREs and the cold-related cis-acting element LTRE, but not DRE/CRT elements are also highly enriched in leaf ABA-induced genes. A motif over-represented in light-induced promoters, GCCAC (SORLIP1AT in the PLACE database [46]), also is found to be highly enriched in guard cell ABA-regulated genes but not in leaves. Compared with other motifs, CE3 has only a few matches in the ABA-regulated genes of either sample type. This is consistent with the findings that CE3 is well-represented in rice but almost absent in Arabidopsis [50]. None of the cis-acting elements enriched in ABA-induced genes is significantly enriched in ABA-repressed genes of guard cells or leaves, except the MYC binding site and LTRE. The mechanisms of ABA-repression are less well characterized than those for induction [33] and motifs specific for ABA-repressed genes may exist but have not been revealed yet (see next section).
We also examined the positional preference of the cisacting regulatory elements ABRE, DRE, CBF and LTRE in a 1000-bp region upstream of identified ABA-regulated genes of guard cells and leaves ( Figure 3). There is no significant difference between ABRE locations in guard cell vs. leaf genes ( Figure 3A); the same is true for the cis-acting elements DRE, CBF and LTRE locations in guard cell genes and leaf genes (not shown). The positional distribution of these four elements upstream of guard cell ABA-induced genes is shown in Figure 3B. ABREs are more likely than the other three motifs to be located in regions close to transcription start sites (TSS), with around 28% and 20% ABREs located in the -100~-50-bp region upstream of the TSS in guard cell and leaf ABA-induced genes, respectively.

Identification of new promoter motifs
In guard cells, 39 of the ABA-induced genes have none of the known motifs of Table 1 in their promoter regions, and this is also true for 79 of the genes induced 'All genes' means all genes represented in the Affymetrix ATH1 chips. "GC_up" and "GC_down" represent ABA-induced and ABA-repressed genes in guard cells, respectively. "LF_up" and "LF_down" mean ABA-induced and ABA-repressed genes in leaves, respectively. The numbers in brackets in the header denote the total number of ABA-regulated genes in the corresponding category. P-values smaller than 0.01 are marked in bold.
by ABA in leaves. In addition, not all genes containing a known ABA-related promoter motif show ABA-regulation in both tissue types. Accordingly, we sought to identify possible new motifs that may be involved in ABA-regulated gene expression, particularly of a cellspecific nature. We were also interested in the possibility of identifying repression-specific motifs. We identified all 5-10mer motifs highly enriched in the 1000-bp sequences upstream of ABA-regulated genes in guard cells and leaves (see Methods), and also compared their enrichment in previously published ABA-regulated gene sets. The statistically significant 5-mer motifs are summarized in Table 2. These motifs include not only the core of G-box-containing ABREs (ACGTG), subsequences of ABREs, the LTRE (CCGAC), and the lightinduced motif (GCCAC), as expected from the above analysis, but also some new, previously unidentified, motifs. The results for 5-10mer motifs also show that there are many more significant motifs in ABA-induced genes than in ABA-repressed genes. In addition, a greater number of distinct motifs were found in guard cell ABA-induced genes than in leaf ABA-induced genes.
To check whether newly identified 5-mer motifs significantly enriched in our gene sets are guard cell-or leaf-specific, we merged previously published ABA-regulated gene sets in Arabidopsis from [18][19][20][21][22][23][24][25][26][27] and compared the enrichment of these motifs. Note that none of these studies uses guard cells or rosette leaves as plant materials. We calculated the enrichment P-values of these motifs in the merged ABA-regulated gene set using all genes in the whole genome as the background, and also calculated the enrichment of these motifs in our gene sets using the merged ABA-regulated gene set as the background. The results are shown in Table S1 and Table S2 in Additional file 3, respectively. Guard cell-specific or leaf-specific motifs are designated as those that are not significantly enriched in the above merged ABA-regulated gene set, but are still significant in our guard cell or leaf ABA-regulated gene set when using the merged ABA-regulated gene set as the background. We find that GTCGG, which is not a subsequence of any known motifs, is significantly enriched in guard cell ABA-induced genes using either all genes in the whole genome or the merged ABA-regulated gene set as the background. By contrast, the enrichment of GTCGG in the merged ABA-regulated gene set is not significant (larger than 1.0E-03). We accordingly identify GTCGG as a potential guard cell-specific motif involved in ABA-induced gene expression. Similarly, TGCAA is a potential leaf-specific motif involved in ABA-repressed gene expression.
Although only a few motifs emerge as significant from the promoters of ABA-repressed genes, in our analysis we find a novel strong candidate for a repression-specific motif, CAAGTTG, which is enriched in both guard cell ABA-repressed genes and leaf ABA-repressed genes (P = 6.0E-09 and 1.6E-08, respectively). This motif is very similar to the motif CAACTTG identified in [51], but in our guard cell and leaf gene sets, the enrichment of CAAGTTG is stronger than that of CAACTTG (P = 6.6E-05 and 1.7E-04, respectively). Known ABArepressed genes KAT2 and MYB60 [14,16] contain the CAAGTTG motif. To examine whether this motif is also enriched in ABA-repressed genes in other tissues/ cell types, we collected eight ABA-repressed gene sets from [19][20][21][22][23][24][25][26], and compared this motif's enrichment in these ABA-repressed gene sets (studies that report a very small number of ABA-repressed genes are not included here). As shown in Figure 4, CAAGTTG is highly enriched in the upstream sequences of most of the eight ABA-repressed gene sets, and the combined motif CAA(G/C)TTG is consistently enriched in all ABA-repressed gene sets, which further supports its involvement in ABA-mediated gene repression.

Comparison with other ABA transcriptome studies
Many transcriptome studies have been conducted on ABA regulation of gene expression in different Arabidopsis tissues and cell types, including whole plants [18,19], seedlings [20][21][22][23][24][25][26], seeds [27], and one previous study reporting ABA responsive genes in guard cells and mesophyll cells using a partial-genome (~8K) Affymetrix chip [16]. 5789 ABA-induced genes and 5635 ABArepressed genes are obtained by merging the ABAresponsive genes in these studies. We examined how many of the ABA-regulated genes that we identified from our guard cell transcriptomes can be found in this merged ABA-regulated gene set. The overlap of our ABA-regulated genes with the previously published ABA-regulated genes is shown in Figure 5. We find that 69.4% of the ABA-induced genes and 54.5% ABArepressed genes in guard cells identified in this work are found in previous ABA transcriptome studies. In leaves, 86.6% of ABA-induced genes and 58.0% ABA-repressed genes are found in previous ABA transcriptome studies. Such significant overlap supports the reliability of our microarray analysis. In addition to the overall comparison, we also conducted a pairwise comparison of the 11 ABA transcriptome studies mentioned above and our study, which describe a total of 14 different ABA-regulated transcriptomes (experiments) (Leonhardt et al. [16] contribute two experiments, and we contribute two as well). These transcriptome experiments have been done under various experimental conditions, including different ABA doses and treatment durations, tissue/cell types, microarray platforms, plant developmental stages, and methods for identifying ABA-regulated genes. A summary of experimental conditions for these ABA transcriptome studies is given in Table S3 (Additional file  3). The pairwise overlaps are measured by percentage (Additional file 3) and representation factor (see Methods), and their significance is examined by the hypergeometric distribution. The results are shown in Figure  S1 (Additional file 3). Given the diversity of experimental parameters, it is not surprising that these studies do not have very high pairwise overlap percentages for ABA-regulated genes. However, when compared with random overlaps, all of them are significant for ABAinduced genes. For ABA-repressed gene sets, all have significant pairwise overlaps as compared with the expected overlap of two random sets except for the ABA-repressed gene sets in [16,18,27] which have nonsignificant overlaps with several other studies. Figure 6 shows the cumulative number of ABA-regulated genes that are common to different numbers of experiments. As expected from the pairwise overlap analysis, there are not many ABA-induced genes found in common as the number of experiments increases, and there are no ABA-repressed genes that are common to >11 experiments (transcriptomes). However, despite the diversity in tissue types, ABA treatments, and transcriptome platforms utilized in these experiments, we were able to find a "core set" of~50 ABA-induced genes reported in any 9 out of 14 experiments ( Figure 6A), and a similar number of "core set" ABA-repressed genes reported in any 6 out of 14 experiments ( Figure 6B). Table 3 lists the core set of ABA-induced genes, and Table 4 lists the core set of ABA-repressed genes.

Identification of new guard cell ABA-responsive genes
In addition to identifying core sets of ABA-induced and ABA-repressed genes, we also discovered a number of novel ABA-responsive genes in both guard cells and leaves. In particular, our results provide a valuable comprehensive source of ABA-regulated genes in guard cells since the smaller Affymetrix~8K array was used in [16] and therefore could not provide information on~2/3 of the transcriptome. For example, among 1173 genes significantly ABA-regulated in our guard cell transcriptome, 1122 genes (658 ABA-induced and 464 ABArepressed) were not previously reported to be ABAregulated by Leonhardt et al. [16]. Of these 1173 genes, 148 ABA-induced genes and 149 ABA-repressed genes were also not reported as ABA-regulated in any of the 11 previous transcriptome studies we analyzed (see  identified ABA-induced genes of guard cells, consistent with the hypothesis that this motif participates in ABAregulation of gene expression in guard cells. We also identified a smaller set of ABA-regulated transcripts in leaves that have not been reported in the previous ABA microarray studies: 62 ABA-induced genes and 113 ABA-repressed genes. In Table 5 and Table 6 we list by P-value the top 50 new ABA-induced genes and ABArepressed genes in guard cells not previously reported by any of the 11 previous ABA microarray studies that we evaluated. Additional file 4 and Additional file 5 list all the new ABA-responsive genes in guard cells and leaves, respectively. Functional enrichment in ABA-responsive genes ABA is involved in a myriad of biological processes, including plant growth and development, and plant responses to environmental stresses such as cold, drought, and pathogens. To examine whether ABAresponsive genes identified in our study tend to be involved in such aspects of plant physiology, we used the GO functional analysis tool BiNGO 2.42 embedded in the Cytoscape project [52,53] to functionally categorize our ABA-responsive genes in guard cells and leaves. TAIR9 GO data (GOSlim_Plants) [54] for the whole Arabidopsis gene annotation were used as the reference set. The significance of overrepresentation and underrepresentation of a functional category was calculated by the hypergeometric cumulative distribution. Bonferroni familywise error rate (FWER) correction was applied for multiple testing. The functional categories with P-values less than 1.0E-03 that also have FWER -corrected Pvalues smaller than 0.05 are given in Figure 7. For the ABA-induced genes, the three most significant functional categories in both guard cells and leaves are response to abiotic stimulus, response to endogenous stimulus, and response to stress ( Figure 7A), all consistent with the fact that ABA plays important roles in the regulation of gene expression during stress responses. Lipid metabolic process, catalytic activity and secondary metabolic process are also enriched in both guard cell and leaf ABA-induced genes but the enrichment of vacuole and transcription factor activity is unique to guard cells. In addition, several functional categories are underrepresented within ABA-induced genes, including translation, RNA binding, protein metabolism, and nucleobase metabolism (see Discussion).
ABA-repressed genes show enrichment of different functional categories than ABA-induced genes ( Figure  7B). The functional categories enriched in both guard cell and leaf ABA-repressed genes are plastid, transferase activity, plasma membrane, and nucleotide binding. Response to abiotic stimulus, response to stress and several other functional categories are highly enriched in guard cell ABA-repressed genes but not in leaf ABA-repressed genes. Two functional categories, nucleobase metabolic process and protein modification process, are enriched in leaf ABA-repressed genes but not in guard cell ABArepressed genes. The functional categories protein metabolism, nucleobase metabolism and nucleotide binding which are underrepresented in leaf ABA-induced genes are overrepresented in leaf ABA-repressed genes.  AT5G66770 scarecrow transcription factor family protein 6 The first column lists the AGI locus identifiers of the genes. *Genes that appear in our guard cell ABA-repressed gene set.^Genes that appear in Leonhardt et al.'s guard cell ABA-repressed gene set [16]. The genes that appear in our ABA-repressed gene set are marked with underlines if they are not represented on the Affymetrix~8K chips used in [16]. The second column gives brief TAIR descriptions for the genes. The third column indicates the number of experiments that the genes are found common to. The fold changes are calculated by averaging three biological replicates and comparing the average expression values with ABA treatment to those without ABA treatment. The genes with underlines also appear in our leaf ABA-induced genes. The fold changes are calculated by averaging three biological replicates and comparing the average expression values with ABA treatment to those without ABA treatment. The genes with underlines also appear in our leaf ABA-repressed genes.

Guard cell signaling and development genes regulated at the transcript level by ABA
A number of genes and proteins in guard cells have been revealed to function in stomatal movements in response to stimuli, including ABA, CO 2 , light and pathogens. In [55], a total of 67 guard cell signaling genes was compiled from published literature. Li and Assmann [56] reviewed selected guard cell signaling and development genes. In addition, Kwak et al. collected 69 guard cell signaling genes [13]. We combined these sets as well as 22 additional genes with functional guard cell roles as described by Kim et al. [14] and 24 genes involved in guard cell patterning and development as described by Dong and Bergmann [57], to obtain a list of 149 known guard cell signaling or development genes/proteins (Additional file 6). Of these guard cell components, 29 are significantly regulated by ABA at the transcript level ( Figure 8).

Cross-regulation by ABA and heterotrimeric G proteins
Heterotrimeric G-proteins, composed of α, β and γ subunits, participate in a wide range of crucial signaling pathways in eukaryotes [58,59]. ABA signaling is known to interact with heterotrimeric G-protein signaling in both developmental and stress responses [60][61][62][63]. In a previous study, we found that ABA signaling also crosstalks with G protein signaling at the level of the transcriptome [31]. In the present study, we used the genome-wide guard cell microarray data from wild-type, Gα (gpa1), Gβ (agb1) and Gγ double (agb1 gpa1) mutant plants, with and without ABA treatment to systematically investigate ABA regulation of gene expression in guard cells and compared it with the ABAregulated transcriptome of leaves from the same genotypes. We identified a number of G-protein-independent ABA-regulated genes and G-protein-dependent ABAregulated genes, which will facilitate the screening and identification of novel guard cell signaling genes. Figure  9 demonstrates the heat map of expression patterns of these genes (drawn by Matrix2png [64]) and shows that in both guard cells and leaves the number of G-proteinindependent ABA-regulated genes is dominant, consistent with the fact that ABA mediates many stress and developmental signaling pathways without any known participation from the heterotrimeric G protein. The number of G-protein-dependent ABA-regulated genes in leaves is much larger than in guard cells, suggesting that the G-protein may be more involved in ABA signaling in mesophyll cells than in guard cells. Six (out of 9) guard cell G-protein-dependent ABA-induced genes and 14 (out of 58) leaf G-protein-dependent ABA-induced genes have ABREs in their promoter regions, which is statistically significant.

Cross-regulation by ABA and other hormones
The phenomenon of cross-talk between plant hormones is well-established at the physiological level [10,65,66]. Nemhauser et al. [22] performed an analysis of seven hormone-treated transcriptomes generated by Goda et al. [67] from 7-day-old Arabidopsis seedlings and identified target genes regulated by the hormones abscisic acid (ABA), gibberellic acid 3 (GA), indole-3-acetic acid (IAA, auxin), 1-amino-cyclopropane-1-carboxylic acid (ACC, ethylene precursor), zeatin (CK, cytokinin), brassinolide (BL, brassinosteroid), and methyl jasmonate (MJ, jasmonate). Six of these seven hormones both positively and negatively regulate significant numbers of genes at the transcript level, identified by stringent analysis using two different methods [22], whereas the seventh hormone, GA, regulates only a few genes identified by a low-stringency linear model [22]. We took the intersection of the hormone-regulated gene sets identified by the rank product and linear model methods used Lipid metabolic process (36,26) Response to stress (106, 107) Response to endogenous simulus (67,55) Response to abiotic stimulus (92,79) GC LF in [22] and compared our ABA-regulated genes in guard cells and leaves with these hormonally regulated genes. The overlap of our ABA-regulated genes with hormoneregulated gene sets was calculated both in the same direction, that is, genes up-(down-) regulated by a hormone are compared with genes up-(down-) regulated by another hormone, and in the antagonistic direction. Representation factor and P-values calculated by the hypergeometric distribution (see Methods) were used to evaluate the overlaps. Figure 10 shows the significance of the overlaps of our ABA-regulated genes in guard cells and leaves with the six hormone-regulated gene sets of Nemhauser et al. [22]. Our ABA-regulated genes in both guard cells and leaves have highly significant overlap with ABA-regulated and MJ-regulated genes in [22]. ABA-induced genes have a larger overlap with MJ-induced genes in leaves than in guard cells. Although the numbers are smaller, guard cell ABA-induced genes also have a significant overlap with IAA-induced genes, and leaf ABAinduced genes have a significant overlap with ACCinduced genes. Guard cell ABA-repressed genes have a significant overlap with BL-repressed genes. Overlap of our ABA-regulated genes with hormone-regulated gene sets in an antagonistic way, although not significant, is observed in each hormone-regulated gene set (Table S4 in Additional file 3). In some instances, the extent of antagonistic overlap is similar to or even larger than that corresponding to the comparison in the same direction. For example, ABA-repressed genes in guard cells and leaves have 10 and 8 genes in common with IAAinduced genes, and ABA-induced genes in guard cells and leaves have 14 and 13 common genes with MJrepressed genes. Nemhauser et al. defined marker genes as those genes specifically regulated by one hormone in a high stringency analysis and not by any other hormone even in a lower stringency analysis [22]. However, some genes designated as marker genes specific for hormones other than ABA are actually found to be ABA-regulated in our guard cell or leaf transcriptomes (Table  S5 and Table S6 in Additional file 3).

Discussion
Our identified ABA-regulated guard cell and leaf transcriptomes are supported not only by Q-PCR verification but also by their significant overlap with the ABA- The numbers over the bars represent the numbers of our ABA-regulated genes that overlap with the hormone-regulated genes in [22]. regulated transcriptomes of other tissues with respect to the promoter motifs that are utilized, the genes that are regulated, and the functional roles of the encoded gene products in stress tolerance as ascertained from GO analysis. However, as discussed in detail below, a deeper study of these new transcriptomes reveals many interesting features.

Promoters
As expected, known promoter motifs conferring ABAand stress-related gene regulation, such as ABREs and LTREs, are over-represented in the promoters of both our guard cell and leaf ABA-induced gene sets. However, promoters of ABA-regulated genes of guard cells but not leaves are enriched in MYB and DRE/CRT elements. These differences suggest differential utilization of specific transcription factors in guard cells vs. leaves. Functional support for this hypothesis is offered by published analyses of MYB family members and stomatal regulation. Two members of the Arabidopsis R2R3-MYB family, MYB15 and AtMYB44, are expressed in guard cells and their overexpression lines show enhancement of ABA-induced stomatal closure and drought tolerance [68,69]. Comparable functional data detailing guard cell roles for the DRE/CRT binding factors are presently lacking; our analysis of promoter motifs suggests that DRE binding proteins will also be found to participate in guard cell ABA responses. Analyses in other tissues have shown that genes with DRE elements in their promoters are regulated by dehydration but not by ABA, although ABA-sensitivity may be conferred by the concerted action of the DREB and ABF proteins [70]. When our guard cell ABA-regulated gene set is assessed for its complement of known transcription factors, we find that DREB1A, ABF3 and ABF4 show concerted ABA-induction in guard cells. We hypothesized that, given the specialized nature of guard cell responses to ABA, these cells might also utilize specialized promoter motifs. One motif, GTCGG, which is not a subsequence of any known motifs, was thereby provisionally identified, awaiting functional verification. In addition, a greater number of different enriched motifs are found for guard cell ABA-induced genes than for leaf ABA-induced genes, suggesting that transcription factors in guard cells may target a greater diversity of cis-acting regulatory elements in ABA responses as compared to those in leaves.
Only a few motifs, including CAAGTTG, emerged as significant from the promoters of ABA-repressed genes. For both guard cells and leaves, there are more significant 5-10mer motifs in ABA-induced genes than in ABA-repressed genes. Since the numbers of ABAinduced and ABA-repressed genes are of similar magnitude in our datasets, this result suggests that evolution may have favored a greater diversity in mechanisms for regulation of gene repression than for gene induction, thus is it more difficult for any one ABA-repression motif to achieve statistical significance.

Comparison with other ABA transcriptome studies
The fact that many classic ABA-and stress-related promoter motifs are enriched upstream of genes in our ABA-regulated gene sets led us to predict that many of the genes identified by our analyses would also be present in comparable transcriptomes derived from other tissues. Indeed, we find over half of all ABA-regulated genes in our guard cell and leaf datasets were also identified in other studies, with the highest percentage overlap (86.6%) found for ABA-induced genes in our leaf gene set, probably reflecting the fact that leaf tissue is included as a component tissue (e.g. as a component of seedlings and whole plants) in several previous studies. On the other hand, it is interesting that there is only a single ABA-regulated gene, AT1G01470: LEA14 (LATE EMBRYOGENESIS ABUNDANT 14) found in common to all ABA transcriptome studies that we analyzed. This result hints that there might be greater flexibility in how tissues achieve a common end result of stress tolerance than in how they perform other functions such as primary metabolism.
Despite the lack of genes found in common to all studies, we nevertheless identified core ABA-regulated gene sets. Many of these core ABA-induced genes have already been assigned functional roles in response to ABA or another stress, or are members of a family (e.g. PP2C phosphatases [33][34][35]) already known to play such a role. In contrast, fewer ABA-repressed genes (Table 4) have known relationships to stress signaling or tolerance. Taken together with the observation that it has been more difficult to identify promoter motifs for ABA-repression than for ABA-induction, and with the fact that we had to reduce the number of experiments to 6 to identify a shared set of~50 ABA-repressed genes, these results suggest a greater between-tissue diversity in the roles of ABA-repressed genes than those of ABA-induced genes. For example, our guard cell dataset includes 56 of the 67 core ABA-induced genes but only 24 of the 47 core ABA-repressed genes.
During drought stress, it is optimal for guard cells to lose water to promote stomatal closure while it is optimal for other cell types to retain water to maintain cellular hydration. We therefore expected to find some ABA-regulated genes unique to guard cells, and this expectation was satisfied. We identified~150 ABAinduced and~150 ABA-repressed genes in guard cells that were not identified as ABA-regulated genes in any of 11 previous studies of ABA-regulated transcriptomes (see Table S3 in Additional file 3). These genes represent potential guard cell-specific ABA-regulated genes. Consistent with this hypothesis, our predicted guard cell-specific motif GTCGG is significantly enriched in the upstream regions of these genes (P < 1.2E-04). Supporting the idea that these genes may have novel, as yet undescribed ABA-related functions in guard cells, many fewer genes in the lists of new and most significantly ABA-regulated guard cell genes (Tables 5, 6) are annotated as ABA-or stress-related than in the lists of core ABA-regulated genes (Tables 3,  4).
Additional file 4 provides the full list of the newly identified ABA-responsive genes of guard cells. Some of these genes already have been reported to be involved in stomatal movements or guard cell signaling. JAR1 (AT2G46370) is Jasmonate Resistant 1 and has been detected in guard cells [16]. This gene product has been reported to enhance ROS production and ABA sensitivity of guard cells [71] and its expression is ABA-induced in our guard cell data. Other hormone-regulated genes, such as the IAA inducible gene IAA16 (AT3G04730) and cytokinin response factor 3 CRF3 (AT5G53290), are also regulated by ABA in guard cells according to our microarray analysis. In addition, several Ca 2+ -related genes are identified as new ABA-induced genes in guard cells, such as AtCAMBP25 (AT2G41010) and CBL2 (AT5G55990). CRY1/BLU1 (AT4G08920) encodes the blue light photoreceptor cryptochrome 1 and has been reported to participate in blue light-regulation of stomatal movements [72], while NPH3 functions downstream of the PHOT1/PHOT2 blue light photoreceptors which mediate blue light-specific stomatal opening [73]. Our results indicate that these genes are ABA-repressed in guard cells, which might contribute to the phenomenon of ABA-repression of light-induced stomatal opening, particularly over a longer time scale. Sucrose is a significant osmoticum maintaining stomatal opening in the afternoon [74], and the sugar transporter STP1 has already been described as a guard cell-specific transporter that shows a peak in expression in the afternoon [75], in concert with the daily timing of sucrose accumulation into guard cells [74]. The observed repression of STP1 and the putative sugar transporter AT3G05160 by ABA might also contribute to long-term inhibition of stomatal opening.

GO functional analysis
The GO category lipid metabolic processes is over-represented in ABA-induced genes of guard cells and leaves. Little is known concerning the relationship between ABA and lipid signaling in mesophyll cells; however, the guard cell result is consistent with functional evidence for roles of inositol phosphates (produced by the action of phospholipase C) and phosphatidic acid (produced by the action of phospholipase D) in ABA inhibition of stomatal opening [76][77][78][79]. In fact, among ABA-induced genes in this category we find AT5G58700, a phosphoinositide-specific phospholipase C family protein, and AT2G22240, an inositol-3-phosphate synthase, as well as two PLD isoforms: PLDδ (AT4G35790; also confirmed in our Q-PCR analysis) and PLDζ (AT3G05630), suggesting that activation of these genes may contribute toward long-term suppression of stomatal opening under drought conditions.
Conversely, several functional categories are underrepresented within the ABA-induced gene sets, including translation, RNA binding, protein metabolism, and nucleobase metabolism, consistent with the supposition that metabolism is generally down-regulated under stress conditions [80]. However, such down-regulation may be greater for leaves than for guard cells, given that GO categories associated with several types of metabolic processes are over-represented in leaf but not in guard cell ABA-repressed gene sets. This difference may reflect a requirement for guard cells to maintain metabolic activity and active stomatal regulation under drought conditions. One interesting observation is that the GO category response to abiotic stimulus is over-represented in both ABA-induced and ABA-repressed guard cell genes. Many of these genes are genes of known function in guard cells, and are discussed in the next section.
Genes of known guard cell function, G protein regulation, and hormone cross-talk Our analysis reveals that many genes which encode proteins with known functions in guard cell physiology are also ABA-regulated at the transcript level. This phenomenon was reported less robustly in a previous analysis [16] which queried only~1/3 of the Arabidopsis genes assayed here. Changes in rates of solute transport resulting from non-transcription-related regulation of ion channels and transporters are central to rapid, osmotically-driven control of stomatal apertures [81,82]. However, our data show that many of the relevant transport proteins are regulated at the transcript level as well. Transcripts encoding the H + ATPase OST2, the inward K + channels KAT1 and KAT2, the sugar transporter STP1, the endosomal Na + /H + antiporter CHX20, and the nitrate importer CHL1 are all down-regulated in guard cells following ABA treatment, suggesting that their downregulation may be an important aspect of longer-duration inhibition of stomatal opening by ABA. Two known genes related to guard cell development, FAMA (FMA) which regulates proliferation of stomatal precursors as well as differentiation of guard mother cells [83], and ATMPK3, an environmentally responsive mitogen-activated protein kinase mediating stomatal development and patterning [84], are downregulated, as are several signaling proteins. Several other genes, functioning in the guard cell response to CO 2 (HT2, HIC) and pathogens (e.g. FLS2), are slightly repressed by ABA (not shown), although failing to meet the threshold of significance used to categorize ABA-responsive genes in our microarray analyses.
ABA up-regulated genes include zeaxanthin epoxidase (ABA1), required for ABA synthesis. In addition, up-regulation of ABA-related transcription factors ABF3, ABF4, and AtMYB44 could provide positive feedback to the ABA response by promoting transcription of stressrelated genes. A number of genes encoding known guard cell signaling proteins are up-regulated, including four protein phosphatase 2C's: PP2Ca, HAB1, ABI2, and ABI1. Since PP2C phosphatases are negative regulators of ABA signaling, their ABA-induced expression suggests instigation of negative feedback on the ABA response. This supposition is supported by the converse observation that five out of six genes encoding PYR/ PYL/RCAR type soluble ABA receptors with documented function in guard cells [85][86][87] are down-regulated by ABA (out of the 14 receptors in PYR/PYL/RCAR family, only ten are represented on the ATH1 chip), as are two other genes encoding positive transducers of the ABA signal, the MAPK MPK3, and an enzyme involved in nitric oxide production, NIA2.
Analysis of T-DNA insertional mutants has implicated heterotrimeric G proteins in ABA responses of diverse tissues [58,59,88]. Mutant plants with G-protein subunit knockouts are ABA hyposensitive in aspects of guard cell ion channel and stomatal regulation [60,63,89,90] but, unexpectedly, ABA hypersensitive in inhibition of seed germination and root growth [61,91], indicating system-specificity of G-protein effectors. Interestingly, the leaf transcriptome actually has more genes co-regulated by ABA and G-proteins than the guard cell transcriptome. This result suggests important but as yet unidentified functions of ABA/G-protein co-regulation in mesophyll cell physiology. Conversely, in guard cells, this result implies that non-ABA related functions of Gproteins should be investigated. Indeed, based on phenotypes of gpa1 mutants, G-proteins also participate in CO 2 and pathogen responses of guard cells [92,93].
Our results indicate that signaling cross-talk also occurs at the transcript level between ABA and other plant hormones. ABA-regulated and MJ-regulated genes identified by Nemhauser et al. [22] each have highly significant overlap with our ABA-regulated genes in both guard cells and leaves, and this is consistent with ABA-MJ crosstalk at the physiological level, in both guard cells [71,[94][95][96] and other tissues [97]. In addition, a few genes designated as marker genes for hormones other than ABA in [22] are found to be significantly ABAregulated in our transcriptomes. It is possible that the previous analysis [22] was confounded by the use of whole Arabidopsis seedlings as the material for microarray analysis, which may have masked hormonal co-regulation occurring in specific tissue or cell types. Although our evaluation of overlap is limited to one hormone (ABA), these results nevertheless indicate that the conclusion of [22], that plant hormones do not regulate common gene sets, merits re-assessment. Our work provides a starting point to re-assess core hormone regulatory modules and hormone marker genes.

Conclusions
This microarray analysis investigates ABA regulation of gene expression in guard cells in a systematic genomewide manner. A number of ABA-regulated genes of guard cells identified here overlap with ABA-regulated genes of other tissues while a subset of them show ABA-regulation unique to this cell type. A unique cisacting motif, GTCGG, associated with ABA-induction of gene expression specifically in guard cells, was identified at the bioinformatic level. Many of the genes known to encode ion transporters associated with stomatal opening are down-regulated by ABA, providing one mechanism for long-term maintenance of stomatal closure during drought. We also found examples of both negative and positive feedback in the transcriptional regulation by ABA of known ABA-signaling genes, particularly with regard to the PYR/PYL/RCAR class of soluble ABA receptors and their downstream targets, the type 2C protein phosphatases. In conclusion, the results of this study engender new insights into the basic cell biology of guard cells, reveal common and unique elements of ABA-regulation of gene expression in guard cells, and set the stage for targeted biotechnological manipulations to improve plant water use efficiency.

Microarray hybridizations and Q-PCR
Our microarray data utilized the Affymetrix ATH1 chip and are from hybridizations described in detail in [31]. All of these data are available in the Gene Expression Omnibus (GEO) database http://www.ncbi.nlm.nih.gov/ geo with accession no. GSE19520. To reprise, epidermal peels with guard cells as the only intact cell type were used as the source of RNA for our guard cell microarrays [31] ( Figure 1A and Figure 1B), and leaves from the first and second rosette layer of five-week old plants were the RNA source for our leaf microarrays [31]. ABA treatments were 50 μM × 3 hr. and EtOH was used as the solvent control. If an epidermal peel preparation did not pass quality control parameters for guard cell purity it was not used for RNA isolation, and if an RNA preparation was not of high quality as indicated by Bioanalyzer profiling, it was discarded. Four genotypes were evaluated: Col (wild-type), gpa1-4 mutant, agb1-2 mutant, and agb1-2 gpa1-4 double mutant, where GPA1 (AT2G26300) and AGB1 (AT4G34460) are the Arabidopsis G protein α subunit and β subunit respectively. For each type of sample (guard cells or leaves), three independent biological replicates were obtained, resulting in a total of 48 microarray hybridizations (2 sample types × 4 genotypes × two treatments × 3 replicates).
For Q-PCR analysis, cDNA was aliquoted and kept at 4°C throughout each Q-PCR experiment to avoid discrepancy in the data because of freeze-thaw cycles. Real time PCR was performed using pre-mix containing SYBR-Green intercalating dye (BioRad). Actin was used as an internal control (Actin 2: AT3G18780 and Actin 8: AT1G49240) [98]. The positions of the oligonucleotide primers used for real-time PCR was chosen so that the size of all PCR products was between 100 and 150-bp. The suitability of the oligonucleotide sequences in terms of efficiency of annealing was evaluated in advance using the Primer 3 program. Q-PCR experiments were repeated thrice independently, and the data were averaged. The data obtained were analyzed with IQ5 software (Bio-Rad).

Identification of differentially expressed genes
Our ABA-regulated gene sets consisted of the intersecting set of genes identified by two methods: a Boolean method and a linear model method ( Figure 11).
In [31], we adopted a Boolean framework F(ABA, GPA1,AGB1) = B(ABA, A(GPA1, AGB1)) + C ABA to describe the regulatory modes of ABA and the G protein and identify ABA-and/or G protein-regulated genes. In this framework, GPA1 and AGB1 are Boolean variables which can have two states: 1 denoting 'on' (not knocked out) and 0 denoting 'off' (knocked out). ABA is also a Boolean variable: ABA = 1 indicates the presence of ABA (i.e. ABA treatment) and ABA = 0 represents the absence of ABA (i.e. solvent control). According to this framework, B 4 (ABA, A) = ABA and B 13 (ABA,A) = not ABA determine two Boolean rules describing ABAregulation of gene expression independent of the heterotrimeric G protein. The rest of the Boolean rules B (ABA, A) describe G-protein-only regulation or the coregulation of ABA and the G protein on gene expression. A correlation measure is used to assign differentially expressed genes to these Boolean rules [31]. In the present study, we used 1.5 as the correlation score threshold for ABA-only regulated genes to designate them as G-protein-independent ABA-regulated genes. This threshold controls the false discovery rate (FDR) within 0.001. In addition, those ABA-G-protein co-regulated genes that also exhibit greater than 2-fold expression change (with FDR less than 0.05) in response to ABA in the wild-type are designated as G-proteindependent ABA regulated genes in the present study. In other words, all ABA-G-protein co-regulated genes described here satisfy the condition of showing ABAregulation in the wild-type Col genotype.
To make the identification of differentially expressed genes with respect to ABA more reliable, we also used linear models combined with empirical Bayes methods to determine ABA-responsive genes [32]. This method is implemented by the limma package embedded in the Bioconductor project [99]. We followed the 2 × 4 factorial design process, using +ABA/-ABA as one factor and genotypes agb1 gpa1 mutant, gpa1 mutant, agb1 mutant, and wild-type as another factor. Contrasts of interest (ABA treated wild-type versus control wild-type, control mutants versus control wild-type, ABA treated mutants versus ABA treated wild-type) were extracted and the P-values for moderate t-tests were adjusted by Benjamini and Hochberg's method to control for the false discovery rate (FDR). We chose FDR <0.001 as a cut-off to select significantly differentially expressed genes. The P-value threshold for the contrast of ABA treated wild-type versus control wild-type was set at 0.0001, and the P-value threshold for the contrast of genotypes versus wild-type was set to 0.01. Those genes that are only significantly differentially expressed with respect to ABA are designated G-protein-independent ABA-regulated genes. Those genes that are significantly differentially expressed with respect to both ABA and genotypes are called G-protein-dependent ABA-regulated genes. As with the Boolean method, all genes designated as ABA-G-protein co-regulated from the linear model approach satisfy the condition of showing ABA-regulation in the wild-type Col genotype. The Boolean method identified 1323 G-protein-independent ABA-regulated genes and 32 G-protein-dependent ABA-regulated genes in guard cells, with 1024 and 408 genes, respectively, identified in leaves. The linear model method identified 1379 G-protein-independent ABA regulated genes and 204 G-protein-dependent ABA-regulated genes in guard cells, with 1040 and 998 genes, determined in leaves. The G-protein-independent ABA-regulated genes determined by the two methods (intersection of green and red circles in Figure 11) and the ABA-regulated genes identified by both methods but with inconsistency in the existence of G-protein regulation (intersection of green and yellow circles plus intersection of red and blue circles in Figure 11) are designated G-protein-independent ABA-regulated genes in our final gene list. G-protein-dependent ABA-regulated genes that are confirmed by both methods (intersection of yellow and blue circles in Figure 11) are reported as G-protein-dependent ABA-regulated genes in our final gene list. The genes in each subsector of Figure 11A and Figure 11B are given in Additional file 7 and Additional file 8, respectively. The genes that are present in guard cells or leaves but not regulated by ABA at the transcript level (i.e. not present in any of the 8 subsectors of Figure 11A and Figure 11B) are given in Additional file 9 and Additional file 10, respectively.

Promoter motif analysis
We performed a statistical analysis of the promoter sequences of ABA-responsive genes in guard cells and leaves. The 1000-bp upstream regions of all Arabidopsis gene sequences were obtained from the FTP site of the Arabidopsis Information Resource (TAIR) (TAIR9_up-stream_1000_20090619.txt) [54]. Motif search was performed by writing custom Python scripts http://www. python.org which search a given motif pattern with an exact match without insert or mismatch. The significance of motif enrichment was computed by the hypergeometric cumulative distribution function. Specifically, when we examined the enrichment of a known cis-regulatory element, we counted the number of ABA-regulated genes in guard cells or leaves and the number of all genes in the chips whose promoter regions contain the cis-regulatory element. A P-value was obtained according to occurrences of this cis-regulatory element in ABA-regulated genes and in all genes. Positional distribution of a cis-regulatory element in ABA-regulated genes was obtained by shifting a window of the same length to scan the promoter region of each gene, and recording the occurrence positions. To identify possible new motifs that may be involved in ABA-regulated gene expression, we generated sets of 5-10mer sequences and calculated the enrichment significance of each sequence in our ABA-regulated gene sets. As there are many more significant motifs in ABA-induced genes than in ABA-repressed genes, we set 1.0E-10 and 1.0E-04 as the significance thresholds for new motifs of ABA-induction and ABA-repression respectively.

Comparison of hormone-regulated gene sets
In the comparison of our ABA-regulated genes with previous ABA/hormone transcriptome studies, we used a representation factor (RF) and its associated P-value to evaluate whether the overlap of two gene sets from the same background is significant or not. The representation factor is defined by the ratio of the size of the real overlap to the expected number of common genes between the two gene sets: where |A| is the number of genes in the gene set A, | B| is the number of genes in the gene set B, |A∩B| represents the number of genes common to gene set A and gene set B, and N is the total number of genes in the Arabidopsis genome (in this study, it is the number of genes covered by the Affymetrix ATH1 gene chip).