Email updates

Keep up to date with the latest news and content from BMC Systems Biology and BioMed Central.

Open Access Highly Accessed Research article

A system biology approach highlights a hormonal enhancer effect on regulation of genes in a nitrate responsive "biomodule"

Damion Nero1, Gabriel Krouk1, Daniel Tranchina12 and Gloria M Coruzzi1*

Author Affiliations

1 Center for Genomics and Systems Biology, Department of Biology, New York University, 100 Washington Square East, 1009 Main Building, New York, 10003, USA

2 Courant Institute of Mathematical Sciences, New York, 251 Mercer St, New York, NY, 10012, USA

For all author emails, please log on.

BMC Systems Biology 2009, 3:59  doi:10.1186/1752-0509-3-59


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1752-0509/3/59


Received:26 March 2009
Accepted:6 June 2009
Published:6 June 2009

© 2009 Nero 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

Nitrate-induced reprogramming of the transcriptome has recently been shown to be highly context dependent. Herein, a systems biology approach was developed to identify the components and role of cross-talk between nitrate and hormone signals, likely to be involved in the conditional response of NO3- signaling.

Results

Biclustering was used to identify a set of genes that are N-responsive across a range of Nitrogen (N)-treatment backgrounds (i.e. nitrogen treatments under different growth conditions) using a meta-dataset of 76 Affymetrix ATH1 chips from 5 different laboratories. Twenty-one biclusters were found to be N-responsive across subsets of this meta-dataset. N-bicluster 9 (126 genes) was selected for further analysis, as it was shown to be reproducibly responsive to NO3- as a signal, across a wide-variety of background conditions and datasets. N-bicluster 9 genes were then used as "seed" to identify putative cross-talk mechanisms between nitrate and hormone signaling. For this, the 126 nitrate-regulated genes in N-bicluster 9 were biclustered over a meta-dataset of 278 ATH1 chips spanning a variety of hormone treatments. This analysis divided the bicluster 9 genes into two classes: i) genes controlled by NO3- only vs. ii) genes controlled by both NO3- and hormones. The genes in the latter group showed a NO3- response that is significantly enhanced, compared to the former. In silico analysis identified two Cis-Regulatory Elements candidates (CRE) (E2F, HSE) potentially involved the interplay between NO3- and hormonal signals.

Conclusion

This systems analysis enabled us to derive a hypothesis in which hormone signals are proposed to enhance the nitrate response, providing a potential mechanistic explanation for the link between nitrate signaling and the control of plant development.

Background

Higher plants acquire nitrogen mainly as NO3-. The soil concentration of this mineral ion can fluctuate dramatically in the rhizosphere, often resulting in limited growth and yield [1]. Thus, nitrate signaling constitutes a key point of plant adaptation to environment. This is why nitrate signaling has so far been intensively studied by transcriptomic assays, involving more than 75 ATH1 chips in various background conditions and treatments. Taken together these transcriptomic data showed that NO3--responses are very context dependent [2,3], suggesting that evolution probably built very adaptable and robust networks involved in the integration of NO3- with other signals including light, sugar, and hormones. For instance, as sessile organisms, plants have developed a strong capacity to modulate growth according to nutrient availability. On a molecular scale, this coordination between nutrition and growth can be mediated by the co-control of metabolism and hormonal signaling. For instance, a recent work reports that molecular reprogramming induced by nutritional starvation treatments significantly involve hormone regulated genes [4]. Moreover, it has also been shown that such cross-controls exist between NO3- and: cytokinin (for review see [5]), auxin [2,6,7], and ABA [8]. To date, molecular players underlying those events are still under investigation. One striking example of such coordination at a molecular level is presented by the role of the iso-pentenyl-transferase 3 (IPT3) involved in the critical step of cytokinin biosynthesis. Transcription of the NO3- induced gene IPT3 has been shown to be involved in the production of NO3- induced cytokinins, hypothesized to coordinate shoot growth in response to NO3- provision [9-14].

Root architecture is also under the coordinated control of nutrient availability and hormone signaling [15]. For instance, NO3- controls root branching under various pathways (for review see [16,17]). Hormones have been shown to play important roles in the adaptation of root development to NO3- availability. Indeed, NO3- triggers root colonization in NO3- rich patch of the soil. Zhang et al [18] have shown that this adaptation could involve AXR4, a gene initially demonstrated to be involved in auxin signaling. Later, AXR4 was shown to be involved in targeting the auxin influx transporter AUX1 to the plasma membrane [19]. Thus AXR4 may provide a molecular link between the NO3- signal and auxin signaling through regulating auxin transport. Furthermore, the dual affinity (high and low affinity NO3- uptake) NO3- transporter NRT1.1/CHL1, hypothesized to be a part of the NO3- sensing system [20-23], was previously shown to be regulated by auxin [7]. This evidence uncovers one facet of how the NO3- sensing system is likely tuned by a hormonal/growth signal.

The complexity of the NO3- effect on root development is further complicated by the fact that high NO3- concentrations (50 mM) trigger an almost complete repression of the lateral root development (LRD). Abscisic acid (ABA) seems to be required for this effect, since the NO3- inhibitory effect on LRD is reduced by mutating either the ABI4 or ABI5 genes [8].

Despite these striking examples, very little is known concerning the transcriptional gene regulatory networks involved in NO3- /hormonal cross-talk. Here, as a step towards understanding such transcriptional co-control, we present a computational biology approach (Figure 1) designed to discover genes that are regulated in response to nitrate treatments across a range of background conditions. Further analysis identified a refined subset of gene clusters to be under the control of both NO3- and hormone signaling. We found that genes controlled both by hormones and NO3- are more responsive to NO3- and/or have a significantly higher level of baseline expression, than genes controlled by NO3- alone. This analysis lead us to generate a new hypothesis that hormone signals play a role in enhancing the effects of the NO3- signal. Moreover, we identified cis-regulatory elements (CREs) in the promoters of these genes that are candidates for enhancing the nitrate-regulation of gene expression.

thumbnailFigure 1. Schematic representation of the in silico strategy used to decipher a N/hormone crosstalk module. (A) The strategy was to discover via sequential biclustering filtering stages, co-regulated genes in response to: i) N treatments, and also ii) hormone treatments. (B) This delimitation of gene groups (whole genome, NO3--regulated, Bicluster9 NO3--regulated and NO3-/hormone co-regulated) and the comparison of their behavior allowed us to hypothesize how NO3- and hormonal signaling pathways are connected in gene networks.

Results

Bicluster analysis identifies a nitrate-responsive "biomodule"

In a previous study, meta-analysis of transcriptomes of NO3- treated plants revealed that gene responses to nitrogen were very context-dependent, and only a very small number of core genes are regulated by NO3- in a context-independent manner [2]. The underlying rules of such coordination/context dependence between signals had recently been proposed at a genome wide level concerning the interaction of carbon, nitrate, and light [3]. Moreover, in light of the context dependent nature of the N-response, mono-dimension clustering algorithms will miss genes that are co-regulated by N across a subset of treatment conditions. By contrast, an approach also known for decades [24] called biclustering can be used to identify nitrate responsive genes that are co-regulated, as a group, in response to a subset of nitrogen treatments across a matrix of meta-data (Figure 1) [25], likely susceptible to tackle the context dependence response to NO3-. Thus, detected biclusters are subsets of the studied genes exhibiting consistent patterns over a subset of N-treatment conditions. Such sets of genes would not be found using mono-dimension clustering approaches, which require that the genes in the cluster behave the same across all treatments. We used this biclustering method to analyze five microarray data sets from N-treatments of Arabidopsis generated by three different laboratory groups: the Crawford lab [26,27] (16 Affymetrix chips including controls), the Stitt lab [28] (14 Affymetrix chips including controls) and the Coruzzi lab [29,30] (46 Affymetrix chips including controls). This combined meta-data set resulting from N-treatments corresponded to a total of 76 microarray chips with controls [For details see Additional File 1].

Additional file 1. Additional Tables and Legends. Contains all Additional Table and legend cited above including: 1) Summary Table of C/N Datasets Used to Generate N-biclusters; 2) BioMaps Analysis of N-biclusters; 3) N-bicluster 9 gene list; 4) Results of N/H Biclustering; 5) Significant H-biclusters Over-represented Functional Categories; 6) Hormone Microarray Datasets.

Format: PDF Size: 145KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

To identify N-regulated genes in this meta-dataset, we first used a filtering step that consisted of selecting genes that are significantly regulated by at least one N-treatment across the various data sets (for filtering conditions see Methods, Figure 1). We found 3,782 such genes under the control of nitrogen in at least one of the considered N treatment conditions. This gene list was then used to generate 21 biclusters (containing 1,124 genes) using the SAMBA algorithm as implemented in the Click and Expander software package [31,32] (Figure 2). In order to identify whether the N-regulated biclusters had biological significance, we used the BioMaps analytical tool [30] to determine significantly over-represented functional categories using the MIPS annotation [33] in each of the 21 N-biclusters generated [Additional File 1]. Significant over-representation was determined based on a p-value that passed a 5% False Discovery Rate (FDR) cutoff (see Methods). Of the 21 N-responsive biclusters, we selected the bicluster # 9 (N-bicluster 9) for further investigation, based on the following criteria: i) genes in bicluster 9 show a reproducible response to N-treatment across several microarray datasets from several different studies (Figure 2), ii) the annotation of genes in N-bicluster 9 suggests that the genes in this bicluster have an important biological function since it comprises genes from all the known steps of N-uptake and N-assimilation (Table 1) [N-bicluster 9 gene list is provided in Additional file 1], and iii) N-bicluster 9 contained significant overrepresentation of MIPS functional categories spanning Amino Acid Metabolism, Carbohydrate Metabolism, and Transported Compounds, suggesting that the underlying regulatory mechanism serves to coordinate systems-wide responses to N (see Table 1). Here, it is noteworthy that this bicluster gathers relevant gene functions. As a matter of fact this coordination of a group of gene involved in the same functions can be designated as a "biomodule" [34-36]. Thus in this context, we used this terminology in the rest of our work.

thumbnailFigure 2. Identification of 21 N-responsive biclusters: Nitrogen-responsive genes sets identified across subsets of N-microarray meta-data. Heatmap of the responsive conditions of N-biclusters. The colors in the heatmap represent the kind of regulation (induced/depressed) and the degree of reproducibility for each data set captured in the N-biclusters. N1 to N5 represent the data sets from 3 different labs. N1:(Wang et al., 2003); N2: (Wang et al., 2004); N3: (Scheible et al., 2004); N4: (Palenchar et al., 2004); N5 (Gutierrez et al., 2007). The details of the experimental conditions of these data set is described in Additional File 1. N-bicluster 9 had the greatest degree of reproducibility over the greatest number of experiments.

Table 1. Over-represented functional categories in N-bicluster 9

To determine the connectivity between the genes in N-bicluster 9, we queried the Arabidopsis multinetwork as described in [30] using the Virtual Plant software package http://www.virtualplant.org webcite. The Arabidopsis multinetwork is a model that integrates information for gene interactions based on a variety of data including: Arabidopsis metabolic pathways, known protein:protein, protein:DNA, miRNA:RNA interactions, and predicted protein:protein and protein:DNA interactions [30]. From the network analysis of 126 genes in N-bicluster 9, a network comprised of information from metabolic and protein:protein connections included 179 nodes (52 genes and 127 metabolites) linked by 261 functional relationships (Figure 3). To test whether this level of connectivity was greater than expected by chance, we computed a p-value by comparing this network to networks generated from a collection of genes of the same size, randomly sampled from the 22,746 genes present on the full genome chip (p-value < 0.0001), computed as described in [37]. Thus, this analysis confirms that the 126 genes in N-bicluster 9 display a significant level of connectivity, providing added support that the biclustering approach has identified a biologically functional regulatory module or "biomodule" [34-36].

thumbnailFigure 3. N-bicluster 9 gene network: A highly connected network of biologically related gene functions. The metabolic layer of information about gene connectivity in the Arabidopsis Multinetwork [30] was queried with the 126 genes belonging to the N-bicluster 9. Nodes represent genes (colored squares) and metabolites (yellow circles) connected by edges (metabolic interactions are colored in grey, protein:protein interactions are colored in green). Genes belonging to: i) N-bicluster NO3- exclusive genes are colored in blue, ii) Significant N/H-biclusters are colored in red, iii) Other N/H-bicluster are colored in grey.

The genes in the N-bicluster 9 "biomodule" are induced by nitrate as a signal

Here, we asked which N-signals were involved in mediating the regulation of genes in N-the bicluster 9 biomodule. Nitrogen-responsive genes are known to be regulated by various endogenous and external nitrogen signals. For example, nitrate itself as been identified as a signal that regulates a large number of genes genome-wide [27]. In specific examples, nitrate can act as an inducer and a repressor of NRT2.1, which codes for a major component of the nitrate high affinity transport system [21,22,38]. Moreover, N-reduced metabolites such as NH4+, and organic forms of nitrogen (glutamine or glutamate) have been shown to also control gene expression, again at a genome-wide level [39,40], which was also demonstrated in specific examples [38,41]. Because NO3- is quickly assimilated into N-reduced and organic-N compounds, we tried to elucidate which form of nitrogen regulates the expression of genes in N-bicluster 9. For this, we mined Affymetrix microarray data obtained from nitrate treatments of a nitrate-reductase T-DNA double mutant knockout from Wang et al [27], which was not included in the biclustering data set. This microarray data set was derived from wild-type (WT) and nitrate reductase (NR) null mutant Arabidopsis plants treated with KNO3 for 2 hrs. The NR double mutant is a T-DNA knockout mutant that is deficient in both structural genes encoding nitrate reductase (nia1 and nia2) and is therefore unable to reduce nitrate to downstream forms of nitrogen [27]. Thus, any genes that respond to nitrate in wild type and the NR double mutant are likely to be controlled by nitrate itself. We found that genes in the N-bicluster 9 biomodule are regulated in response to nitrate, since their nitrate responsiveness is preserved in the nia1/nia2 NR double mutant at levels identical to wild-type (Figure 4).

thumbnailFigure 4. The 126 genes in N-bicluster 9 are NO3- regulated. Average raw MAS5 processed signal values for the 126 genes from N-bicluster 9 in the nitrate reductase double mutant data [27]. Values in these plots indicate the average expression for each condition taken from the nitrate reductase double mutant data set where plants were treated with KNO3 or an equal concentration of KCl (Control).

Identification of nitrate/hormone responsive biclusters

As stated above, the N-regulation of the bicluster 9 "biomodule" appears to have biological relevance as it contains genes covering all steps in the N-uptake/assimilation pathway, and genes integrating the regulation of genes in N-metabolism with other related metabolic processes including C-metabolism and Energy. As many of the genes in these pathways have been shown to be under control of other hormones [7,42,43], we hypothesized that hormonal control might play a role in the regulation of the genes in the N-bicluster 9 biomodule. To test this hypothesis, and to uncover the possible mechanistic basis for nitrate/hormome signaling interactions, we performed an in silico analysis of N-bicluster 9 gene regulation in response to hormones. For this analysis, we constructed a meta-data matrix of microarray experiments from the NASC repository (Nottingham Arabidospsis Stock Center, [44]) covering all available hormone and hormone inhibitor treatment experiments (including auxin, abscisic acid, ethylene, cytokinin, brassinosteroids, giberellic acid and jasmonic acid). We included all data sets having replicates and a relative control (see Methods). The data were converted to a log base 2 ratio of treatment/control, and biclustering was used to identify sets of genes co-regulated under various hormone treatments. Biclustering of the 126 genes in the N-bicluster 9 biomodule over the hormone meta-data set generated 22 hormone biclusters referred to as N/H-biclusters (Nitrogen/Hormone responsive). These N/H biclusters encompass 77 genes contained within N-bicluster 9 [see Additional File 1 for a complete list of N/H-biclusters]. Each of these N/H biclusters has an average number of 8 genes with extensive overlap between several N/H biclusters. These results show that 77 of the genes (out of 126) in N-bicluster 9 are reproducibly regulated by both nitrate and hormone treatments. Further the fact that these genes overlapped between N/H biclusters suggested that they may be responsive to several different types of hormone treatments. The 22 N/H-biclusters were classified according to filtering criteria described in Methods section (Table 2). Out of the 22 N/H-biclusters, five biclusters which showed the most extensive and reproducible responses to hormones were selected for further analysis (Figure 5). The metabolic and functionally interacting genes contained in these N/Hormone biclusters are highlighted in the 77 N/H responsive genes in the network view of the N-bicluster 9, as depicted in Figure 3.

Table 2. NO3-/Hormone Response Interaction.

thumbnailFigure 5. Hormone responsiveness of the five selected N/H-biclusters from N-bicluster 9. Centroid plots of the expression patterns for genes across the 5 significant (> 50% reproducibility; > 1.5 fold change) N/H-biclusters. Black bars below the plots indicate replicate experiments and grey bars indicate singleton experiments. N/H-biclusters 1 and 19 show a strong response to ABA treatment during seed imbibition (induction and repression respectively), N/H-biclusters 16 and 6 show an induction to cytokinin in roots and shoots of mutant and wildtype plants respectively. N/H-bicluster 20 shows a strong depression due to cytokinin in root tissues of mutant and wild type plants. Treatments are taken from the NASC data set: 1- Brassinolide 10 nM 3 hours (Seedlings), 2- ABA 3 μM 24 H (during Seed Imbibition), 3- ABA 30 μM 24 H (during Seed Imbibition), 4- Zeatin 20 μM 1 H (Seedlings), 5- Zeatin 20 μM 1 H (Shoots), 6- Zeatin 20 μM 1 H (Roots), 7- arr10/12 Zeatin 20 μM 1 H (Roots), 8- arr10/12 Zeatin 20 μM 1 H (Shoots), 9- AtIPT8/pga22 (Seedlings), 10- ga1–5 GA3 1 μM 0.5 H (Seedlings), 11- IAA 1 μM 0.5 H (Seedlings), 12- IAA 1 μM 1 H (Seedlings).

The analysis of the N/Hormone biclusters (biclusters 1, 6, 16, 19, 20) revealed that cytokinin and ABA are the main hormone treatments under which the NO3- regulated genes from N-bicluster 9 are co-regulated (Figure 3 and 5). N/H-biclusters 1 and 19 are both mainly driven by ABA treatments, although their respective regulation is in opposing directions (induced vs. repressed). N/H-biclusters 6, 16, and 20 are comprised of genes almost exclusively regulated in response to cytokinin treatment (Figure 5). Together, our results suggest that the coordinated regulation of these genes to nitrate as well as cytokinin or ABA may be part of a regulatory network that mediates the responsiveness of these genes.

Functional interactions within the nitrate/hormone biclusters

As all N/H-biclusters were derived from the 126 genes contained in N-bicluster 9, we used a modified version of BioMaps analysis to determine which if any of the five selected N/H-biclusters were enriched for specific MIPS functional categories (see Methods). This analysis demonstrated that N/H-biclusters 1, 16 and 19 had at least one over-represented MIPS category, when compared to N-bicluster 9 [see Additional File 1]. The most significantly over-represented categories from these N/H-biclusters are genes involved in metabolic pathways, suggesting that this NO3-/Hormone "crosstalk" may be directed towards the coordinate regulation of genes in interconnected metabolic pathways (see Network View of N-bicluster 9). Further, genes from N/H-bicluster 1 have several additional categories over-represented including Energy, Pentose phosphate pathway and Photosynthesis.

Elucidation of a hormone "enhancement" of gene nitrate responsiveness

With the aim of elucidating the mechanisms which mediate differences in regulation between genes that are controlled by NO3- only vs. genes that are co-regulated by NO3- and hormones, we analyzed the expression of the genes N-bicluster 9 and N/H bicluster in the NR double mutant data set [27]. This analysis uncovered a strong and unexpected difference between: a) the NO3- responsiveness of genes in all of the five N/H-biclusters vs. b) genes from N-bicluster 9 that did not appear in any hormone bicluster (termed "N-bicluster 9 exclusive") (Figure 6) (see also informatic analysis scheme Figure 1B). Based on the analysis described below, genes belonging to N/H-biclusters 1, 16, 19 and 20 are significantly more NO3- responsive as measured by the amplitude of expression in genes from N-bicluster 9 that do not fall into any N/H-bicluster (N-bicluster 9 exclusive genes) (see results of ANOVA, Table 2). Further, we validated that the enhanced nitrate responsiveness for the genes belonging to the N/H biclusters, is exhibited not only in the NR double mutant data set, but also in the entire N-treatment meta-data set used to build the N-biclusters [N1 to N5, see Additional File 1] (data not shown).

thumbnailFigure 6. Changes in NO3- responsiveness for genes in significant N/H biclusters vs. N-bicluster 9 exclusive genes. Histogram plots of average expression over NO3- treatment conditions in both WT and NR double mutant microarray data showing the nitrate responsiveness of genes that were present in significant hormone biclusters and genes that were not present in any N/H-bicluster (N-bicluster 9 exclusive (E)), data indicates root and shoot responses in WT and NR double mutant plants). Comparison between treatment and control across the different groups indicates a clear difference in both the baseline expression level and average response to nitrate for N/H-biclusters 1, 16, 19 and 20 and was confirmed with ANOVA (see text). (+) KNO3 treatment, (-) KCl treatment.

To quantify and statistically validate the regulation and nitrate-responsiveness of the genes in the N/H biclusters in the NR double mutant data set (Figure 6), we modeled the expression of genes from these groups using the lm, summary.lm and ANOVA functions in [R] [45]. In this analysis the gene-expression response variable was modeled as a function of 4 explanatory-variable factors: i) Treatment, with 2 levels (nitrate treatment and control); ii) Tissue, with 2 levels (roots and shoots); iii) Genotype, with 2 levels (mutant (NR double mutant) and wild-type); iv) N/H-bicluster, with 6 levels (N/H bicluster 1, 6, 16, 19, 20 and N-bicluster 9 exclusive). To avoid any ambiguity between factor levels, overlapping genes from H-biclusters were removed from the analysis. The response variable (signal values) were taken from the normalized MAS5 data from the Wang et al, dataset [27] (These data were used to build Figure 6). In our ANOVA analysis, we started with an initial model that included main effects for each of the factors and an interaction term for the Treatment and N/H-bicluster factors. We simplified the model systematically in a step-wise procedure as outlined in Crawley [46] and fully described in Additional File 2. Briefly, our results from ANOVA analysis showed that the main effects of Tissue and Genotype factors were not significant (p-values of 0.13 and 0.94, respectively). Further it revealed that the factor levels of N/H-biclusters were not all significantly different from each other. Specifically NH-bicluster 16 19, and 20, in one hand, and 6 and N-bicluster 9 exclusive, in the other hand, were not significantly different from each other so that these levels could be combined into a single compound level. The final result of our simplification procedure was a model with main effects of Treatment (with 2 levels), N/H bicluster (with 3 compound levels), and an interaction term between N/H bicluster and Treatment. The R code and output for this model simplification are fully available in Additional File 2.

Additional file 2. ANOVA [R] script and results. Contains the entire commented procedure of the ANOVA analysis, comparing N/H controlled genes vs. N-exclusive genes.

Format: PDF Size: 49KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Using this final model, we were able to show that the N/H bicluster 1 level and the compound level for N/H biclusters 16, 19 and 20 are both significantly different from the compound level of N/H bicluster 6 and N-bicluster 9 exclusive in both having a stronger baseline response and a stronger response to nitrate (see Figure 6). Taken together our combined bicluster analysis data from N-treatment and hormone-treatment meta-datasets, leads us to propose a new hypothesis that hormone signaling (specifically ABA, and/or cytokinins which are represented in N/H biclusters 1, 16, 19 and 20) may act as an enhancer of NO3- signaling/induction.

Identification of candidate cis-binding elements involved in the hormonal enhancing effect of the nitrate response and in silico validation

In order to determine candidate cis-regulatory elements (CREs) involved in the proposed 'enhancing' effect of hormones on nitrate signaling, we first scanned the ~3,000 bp parsed upstream promoter sequence of the genes in N-bicluster 9 for known transcription factor binding sites using a DNA pattern search tool [47]. We next determined the over-representation of these CREs in the 126 N-bicluster 9 genes using the Fisher Exact Test (see Methods). Based on this analysis, 23 CRE elements were found to be significantly over-represented in N-bicluster 9, compared to the genome-wide frequency of these elements. These CREs could potentially be involved in nitrate responses, or in a yet to be defined signal (or complex of signals) controlling this gene cluster as a whole (Table 3). Using a similar approach we tested for a significant difference in the CRE frequency for these CREs between two groups: one containing the N/H bicluster 1, 16, 19 and 20 genes and the other containing the N-bicluster 9 exclusive (background). N/H bicluster 6 was removed from this and all subsequent analysis as it was shown to not be significant in the previous ANOVA and contained very few genes (2 genes) when overlapping genes were removed. This analysis revealed that 2 CREs, HSE element [48] and E2F element [49]. These two CRE elements are significantly over-represented in the N/H biclusters (See Table 4) analyzed as a group compared to their frequency within the subset of N-bicluster 9 genes not present in any of the N/H-biclusters (termed NO3- exclusive). The specific over-representation of these CREs in the N/H biclusters implicates them as potential candidates to be involved in the cross-talk between hormone and NO3- response. The known physiological roles of the E2F and HSE CREs are discussed below.

Table 3. Over-represented known CREs in N-bicluster 9.

Table 4. Over-represented CREs for 4 significant N/H-biclusters.

Discussion

Biclustering identifies a "biomodule" of biologically related nitrate-regulated genes involved in metabolism and signal transduction

In a previous meta-analysis of nitrate-regulated genes, we demonstrated that a very small number of genes are nitrate regulated across a variety of background conditions, while the vast number of nitrate-regulated genes are regulated in a context-dependant manner [2]. This observation suggests that the NO3- signaling pathway is also under the influence of other (as yet) unidentified controls. Taking this observation as a starting point, we decided to use biclustering technique: an approach that clusters both genes and treatments, as a tool to discover genes that are co-regulated by nitrate across a wide variety of background conditions corresponding to a subset of the meta-data analysis. This biclustering approach allowed us to uncover a "biomodule" of 126 NO3- regulated genes that are related in expression pattern and in biological function (N-bicluster 9, Figures 2 and 3, Additional file 1, table s2). Indeed, as a group, the genes in N-bicluster 9 comprise a set of 52 metabolic genes including, for example, all steps in the pathway of nitrate uptake & reduction (NRT1.1, NRT2.1, NRT3.1, NIA1, NIA2, NIR), as well as genes involved in N-assimilation into organic form (GDH1, ASN2 and GLT1). In addition, the N-bicluster 9 also contains significant overrepresentation of genes involved in Energy, Nitrogen and Carbon metabolism (Table 2). This strong functional coherence of the genes in N-bicluster 9 is illustrated by the interactions between 52 genes in the metabolic/protein interactions shown in the subnetwork (Figure 3).

It is noteworthy that the concept of the 126 genes in bicluster 9 constituting a "biomodule" in our study is comparable to ideas that have been already developed by others in the field of systems biology. For instance, i) Baliga et al. [50] state that "a biomodule is a group of proteins that execute a particular function", and ii) Bonneau et al. [51] also used a biclustering approach (cMonkey) to define "biologically meaningful biclusters". The conjunction of both above definitions match our concept/definition of a "biomodule".

As an insight into potential TFs that regulate the genes in this network, it is noteworthy that N-regulated bicluster 9 contains 17 transcription factors (based on AGRIS transcription factor annotation) whose regulation is by definition correlated with targets in N-bicluster 9, as well as with genes from other functional and unknown categories. N-bicluster 9 also contains other regulatory genes potentially involved in signal transduction such as kinases or phosphatases (6% of the genes from this N-bicluster fall into this category) (Table 1).

Hormones enhance the NO3- responsiveness of genes within the bicluster 9 biomodule: what are the physiological consequences?

To identify potential regulatory mechanisms for the NO3- responses of these genes to by other stimuli, we examined the regulation of the 126 genes in N-bicluster 9 across a metadata set of hormone microarrays. This was done in order to try to understand whether these genes, or subsets of these genes, are coregulated by hormones as well. Hormone treatments have been previously shown to have strong interactions with nitrogen signaling [5]. This analysis identified a subset of 77 genes in N-bicluster 9 that also cluster together across a subset of hormone treatments. The position of these 77 genes present in the N/Hormone biclusters are shown in the context of the metabolic/protein interaction network presented in Figure 3. This view demonstrates a strong potential effect of diverse hormonal controls on the level of response of NO3--controlled metabolic processes [See Figure 3 color coding for nodes: Red squares = significant N/H-bicluster genes (genes from N/H-biclusters 1, 16, 19 and 20), Blue squares = genes controlled only by NO3- and not co-regulated by hormone treatments, based on results of hormone biclustering, Grey squares = non-significant N/H-bicluster genes controlled by hormones (i.e. no reproducible hormone response in N/H biclusters, see Methods)]. This hormonal control of nitrate-regulated genes represents a potential mechanism to fine tune and co-ordinate response levels of genes in a biomodule so that metabolic processes (here N-assimilation, carbon metabolism, and signaling components) can be regulated according to the growth rate of the plant. This is consistent with the observations made for phosphate [52,53], sugar [54], sulfate [55], and iron metabolism [56]. In all of these previous studies, when the hormone receptor is mutated, the response of genes to the nutrient under investigation is maintained, but the hormone response of the same genes is abrogated. This implies that hormonal control of nutrition pathways has a broad effect and controls metabolism as a whole, and is distinct from nutrient signaling. Our current work supports this view and also goes a step further. Indeed, our systems approach has enabled us to derive the hypothesis that hormone signals can interact with NO3- signals to enhance the responsiveness of genes, and we have performed and in silico test of this hypothesis. This hypothesis is based on the finding that genes controlled by NO3- only, were shown to be less responsive to NO3- than genes under the control of NO3- and hormones (Figure 6, Table 2). This kind of interaction has to our knowledge never been reported, and is a particularly novel aspect of for the effect of hormones as they relate to NO3- induction. Although the effect of external hormone supply on genes belonging to NO3- assimilation pathway or sensing system has already been documented, our results propose a new dimension of interaction at the transcriptional level between hormonal and NO3- signaling. The existence of specific links between different nutrient and hormonal signals reported herein is also of particular interest and deserves further investigation.

Putative roles of the E2F and HSE Cis Regulatory Elements (CRES) in mediating cross talk between nitrate and hormone signaling

Our study has identified two putative regulatory elements that are over-represented in the four significant N/H-biclusters identified by ANOVA (Table 2). To identify the potential role of such elements in mediating the hormone enhancing effect on nitrate responsiveness, we performed an in silico analysis aiming at deciphering the potential effect of each candidate binding site. By removing all genes from N-bicluster 9 exclusive gene list that contained these CREs (E2F and HSE), we were able to "virtually" examine their respective role in the enhancement of the baseline and NO3- response by comparing these genes to genes from significant N/H-biclusters 1, 6, 16 and 20. The analysis demonstrated that E2F and HSE CREs are potentially involved in the hormonal enhancing effect of expression of these NO3- responsive genes (Table 4, 5, and 6). To date, the heat shock elements (HSE) were not shown to be involved in the control of N-regulated genes though their role in Arabidopsis in the transcriptional control of responses to heat stress has been extensively studied [57]. However, a heat shock transcription factor HsfA9 has been shown to be under hormonal control in seeds (ABA through ABI3) [58]. This observation leads to the tentative hypothesis that heat shock elements could potentially be involved in conveying a hormonal signal. Moreover, to further have insight into the HSF/hormonal connection we ran a Sungear [59] analysis to decipher if these factors are under any other hormonal controls. To do so, we queried gene annotation for HSF term. We found 21 HSF, and looked to see if they were found regulated by any hormone as reported by Nemhauser et al. [60]. Out of the 21 HSF detected we found that 6 (28%) are regulated by ABA (2 of which are also regulated by methyl-jasmonate), and 1 gene is regulated by cytokinins. This kind of co-regulation might further support the potential connection between HSF and hormonal signals.

E2F binding elements and the role of their associated transcription factors are still poorly understood in plants. However, what is known in plants as well as in other organisms is that these factors (considered in animals as oncogenes) are involved in the control of the cell cycle [61,62]. Remarkably the role of E2F in the control of gene expression related to N-assimilation has already been shown in Arabidopsis, providing an independent validation of our results. Vlieghe et al. [63] demonstrated that the over-expression of the E2Fa-DPa transcription factor leads to the induction of nitrate reductase (NIA2), glutamine synthetase (GS), glutamate synthase (GOGAT), and nitrite reductase (NIR) gene. It is noteworthy that all of these genes respond to both nitrate and hormonal signals in our analysis (Figure 3). Furthermore, several genes in the N/H biclusters that are involved in C-metabolism are also mis-regulated in plants over-expressing E2Fa-DPa. Interestingly, E2F CREs were also identified in the nitrate reductase promoter of the green algae, Chlorella vulgaris. The protein binding activity at this site was validated but was not dependent on nitrate in the media [64]. This confirms the idea that E2F CREs are involved in the interaction of the NO3- response with other signals such as hormones and may be mediating crosstalk between these signals. Finally, the cell cycle is known to be an important target of hormonal signaling. For instance the Arabidopsis E2FC-DPB transcription factor was demonstrated to be involved in the control of the cell cycle. Also, cell division (monitored by CYCB1-GUS) in plants over-expressing E2FC-DPB was found to be less sensitive to auxin than cell division in wild types plants. This supports the hypothesis that E2F transcription factors are involved in mediating hormonal control of cell division [65].

Conclusion

In conclusion, our results suggest and highlight a significant level of control of NO3- signaling by hormones. This control may allow plants to modulate biomodules of genes spanning N and C metabolism according to growth-dependant hormone signals. The systems biology approach presented herein demonstrates the inference of relationships between signals a postriori using extensive microarray data sets (76 chips for Nitrogen + 278 chips for hormones) to uncover new hypotheses for mechanisms underlying the much studied but poorly understood interactions between nutrient and hormone signaling. This in silico approach opens the door toward unraveling new biological concepts by systems analysis of existing microarray and other genome scale data sets within the public domain.

Methods

Nitrogen Microarray Meta-Data set used for bicluster analysis

Expression values for all genes within the Arabidopsis genome present on the Affymetrix chip were taken from published data on nitrogen treatments vs. controls for all the available experiments from the data sets published in:[26-28,30]. All microarray data used in this analysis was processed and normalized using Affymetrix Suite 5.0 or MAS5 Software (as implemented in the R statistical package [45] the two normalization Methods gave equivalent results. For biclustering analysis (see below), signal values were converted to log base 2 ratios with the treatment condition compared to its relative baseline condition (control). Genes with raw signal values less 100 in their treatment or control conditions in either replicate had their signal log ratio values replaced with a non-numerical NA value which is ignored by the biclustering algorithm. Finally, Log 2 ratio data from the microarray data was analyzed to determine which genes in the genome were greater then 1.5-fold responsive in any pair of replicate experiments in this meta data set. The resulting list of 3,752 N-responsive genes was used for biclustering, as described below.

Hormone Microarray Meta Data Analysis

All hormone microarray data was taken from the MAS5-processed NASC Microarray database (Nottingham Arabidopsis Stock Center [44]) data. Data was chosen based on annotation and experimental conditions that referred to a hormone or hormone inhibitor treatment vs. a relative control, with only replicated data used for biclustering analysis. A total of 19 data sets comprising 278 microarray experiments were compiled based on these criteria. The full list of data sets and the contributing number of microarrays from each data set is provided in Additional File 1. All hormone data with signal values of < 100 had their signal values replaced with a non-numerical NA value. Further, all data was converted to log base 2 ratios prior to biclustering analysis. The 126 genes from N-bicluster 9 were biclustered over all hormone data as described (see below).

Biclustering of Microarray Data

Biclustering was performed using the SAMBA algorithm as described by [32] and as implemented in the CLICK and EXPANDER program [31]. Biclustering was performed using default parameters except as follows: amount of overlap allowed (50%), gene coverage (set to cover all genes) and the number of genes expected (set to maximum number of genes in the data set). Biclusters that were used for further analysis were chosen based on genes being ≥ 1.5 fold regulated across reproducible experiments and the presence of replicate experiments for ≥ 50% of the experiments.

BioMaps Analysis of N-biclusters

BioMaps analysis of N-biclusters was performed as described in [30] as accessed via http://www.virtualplant.org webcite. The program was run using the MIPS [33] annotation option for functional definitions. A 5% FDR (False Discovery Rate) cut-off was computed using the R statistical package [45] to determine significant p-values.

Modified BioMaps of N/Hormone Biclusters

As N/H-biclusters were derived from N-bicluster 9, over-representation of MIPS [33] functional terms were determined using the Fisher's Exact Test as implemented in the R statistical package [45] to compare the proportions of genes from N/H-biclusters containing a MIPS term vs. the proportion of genes from N-bicluster 9 not in N/H biclusters containing that same term.

Multinetwork Analysis of Bicluster 9

To understand the relationships among the 126 genes from N-bicluster 9, the Arabidopsis multinetwork analysis tool was used [30] as accessed by http://www.virtualplant.org webcite. This network contains many validated connections for gene interactions in the Arabidopsis genome. Network interactions were visualized using Cytoscape [66].

Computing the False Discovery Rate (FDR)

An FDR control method was used to determine a significance cutoff for p-values. This value was computed using a script written for the R statistical package [45]. This script was derived from the Storey and Tibshirani method [67] which determines a cut-off based on the expected proportion of false positives incurred when calling a feature significant.

Cis Regulatory Element (CRE) Detection Using Known CREs

In order to detect known CREs that may be over-represented in a group of genes (e.g. N-bicluster 9, N/H-biclusters), sequence analysis of the promoter regions of these genes was performed. We used Cis Regulatory Element (CRE) annotation from the AGRIS Database (Arabidopsis Gene Regulatory Information Server [68]) as well as our own literature search to identify biologically active CREs that have been validated by in vivo experimentation. CRE detection was performed using the DNA pattern search tool available from RSA Tools [47] upon 3,000 bp of parsed upstream promoter region (taken from the AGRIS database).

The test for over-representation of CREs was performed using Fisher's Exact Test. This test compared the proportion of promoters in which a particular CRE of interest appeared in one group with the proportion of these same CREs detected in another group. A p-value cutoff was computed using a 5% FDR cut-off for significance.

Authors' contributions

GC, GK, DT and DN designed the study. DT and DN performed statistical analysis. GK and DN wrote the paper. All authors read and approved the final version of the manuscript.

Acknowledgements

We thank Sandrine Ruffel, and Miriam Gifford for helpful discussion and critical reading of the manuscript. We thank Dennis Shasha for his help in determining the statistical significance of the connectivity of multinetworks. This work was funded by NIH NIGMS GRANT GM032877 to G.C., an NIH Pre-doctoral minority fellowship GM032877-S1 to D.N., and an NSF Arabidopsis 2010 Genome Grant (IOB 0519985) to G.C.

References

  1. Crawford NM, Glass ADM: Molecular and physiological aspects of nitrate uptake in plants.

    Trends in Plant Science 1998, 3(10):389-395. Publisher Full Text OpenURL

  2. Gutierrez RA, Gifford ML, Poultney C, Wang R, Shasha DE, Coruzzi GM, Crawford NM: Insights into the genomic nitrate response using genetics and the Sungear Software System.

    J Exp Bot 2007, 58(9):2359-2367. PubMed Abstract | Publisher Full Text OpenURL

  3. Krouk G, Tranchina D, Lejay L, Cruikshank AA, Shasha D, Coruzzi GM, Gutierrez RA: A Systems Approach Uncovers Restrictions for Signal Interactions Regulating Genome-wide Responses to Nutritional Cues in Arabidopsis.

    PLoS Comput Biol 2009, 5(3):e1000326. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Rubio V, Bustos R, Irigoyen ML, Cardona-Lopez X, Rojas-Triana M, Paz-Ares J: Plant hormones and nutrient signaling.

    Plant Mol Biol 2009, 69(4):361-373. PubMed Abstract | Publisher Full Text OpenURL

  5. Sakakibara H, Takei K, Hirose N: Interactions between nitrogen and cytokinin in the regulation of metabolism and development.

    Trends Plant Sci 2006, 11(9):440-448. PubMed Abstract | Publisher Full Text OpenURL

  6. Avery GS, Pottorf L: Auxin and nitrogen relationshinps in green plants.

    American Journal of Botany 1945, 32:666-669. Publisher Full Text OpenURL

  7. Guo FQ, Wang R, Crawford NM: The Arabidopsis dual-affinity nitrate transporter gene AtNRT1.1 (CHL1) is regulated by auxin in both shoots and roots.

    J Exp Bot 2002, 53(370):835-844. PubMed Abstract | Publisher Full Text OpenURL

  8. Signora L, De Smet I, Foyer CH, Zhang H: ABA plays a central role in mediating the regulatory effects of nitrate on root branching in Arabidopsis.

    Plant J 2001, 28(6):655-662. PubMed Abstract | Publisher Full Text OpenURL

  9. Takei K, Ueda N, Aoki K, Kuromori T, Hirayama T, Shinozaki K, Yamaya T, Sakakibara H: AtIPT3 is a key determinant of nitrate-dependent cytokinin biosynthesis in Arabidopsis.

    Plant Cell Physiol 2004, 45(8):1053-1062. PubMed Abstract | Publisher Full Text OpenURL

  10. Miyawaki K, Matsumoto-Kitano M, Kakimoto T: Expression of cytokinin biosynthetic isopentenyltransferase genes in Arabidopsis: tissue specificity and regulation by auxin, cytokinin, and nitrate.

    Plant J 2004, 37(1):128-138. PubMed Abstract | Publisher Full Text OpenURL

  11. Sakakibara H: Nitrate-specific and cytokinin-mediated nitrogen signaling pathways in plants.

    J Plant Res 2003, 116(3):253-257. PubMed Abstract | Publisher Full Text OpenURL

  12. Takei K, Takahashi T, Sugiyama T, Yamaya T, Sakakibara H: Multiple routes communicating nitrogen availability from roots to shoots: a signal transduction pathway mediated by cytokinin.

    J Exp Bot 2002, 53(370):971-977. PubMed Abstract | Publisher Full Text OpenURL

  13. Takei K, Sakakibara H, Taniguchi M, Sugiyama T: Nitrogen-dependent accumulation of cytokinins in root and the translocation to leaf: implication of cytokinin species that induces gene expression of maize response regulator.

    Plant Cell Physiol 2001, 42(1):85-93. PubMed Abstract | Publisher Full Text OpenURL

  14. Rahayu YS, Walch-Liu P, Neumann G, Romheld V, von Wiren N, Bangerth F: Root-derived cytokinins as long-distance signals for NO3--induced stimulation of leaf growth.

    J Exp Bot 2005, 56(414):1143-1152. PubMed Abstract | Publisher Full Text OpenURL

  15. Drew MC: Comparison of the effects of a localized supply of phosphate, nitrate, ammonium and potassium on the growth of the seminal root system, and the shoot, in barley.

    New Phyt 1975, 75:479-490. Publisher Full Text OpenURL

  16. Walch-Liu P, Filleur S, Gan Y, Forde BG: Signaling mechanisms integrating root and shoot responses to changes in the nitrogen supply.

    Photosynthesis research 2005, 83(2):239-250. PubMed Abstract | Publisher Full Text OpenURL

  17. Walch-Liu P, Ivanov II, Filleur S, Gan Y, Remans T, Forde BG: Nitrogen regulation of root branching.

    Ann Bot (Lond) 2006, 97(5):875-881. PubMed Abstract | Publisher Full Text OpenURL

  18. Zhang H, Jennings A, Barlow PW, Forde BG: Dual pathways for regulation of root branching by nitrate.

    Proc Natl Acad Sci USA 1999, 96(11):6529-6534. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Dharmasiri S, Swarup R, Mockaitis K, Dharmasiri N, Singh SK, Kowalchyk M, Marchant A, Mills S, Sandberg G, Bennett MJ, et al.: AXR4 is required for localization of the auxin influx facilitator AUX1.

    Science 2006, 312(5777):1218-1220. PubMed Abstract | Publisher Full Text OpenURL

  20. Remans T, Nacry P, Pervent M, Filleur S, Diatloff E, Mounier E, Tillard P, Forde BG, Gojon A: The Arabidopsis NRT1.1 transporter participates in the signaling pathway triggering root colonization of nitrate-rich patches.

    Proc Natl Acad Sci USA 2006, 103:19206-19211. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  21. Krouk G, Tillard P, Gojon A: Regulation of the high-affinity NO3- uptake system by NRT1.1-mediated NO3- demand signaling in Arabidopsis.

    Plant Physiol 2006, 142(3):1075-1086. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Muños S, Cazettes C, Fizames C, Gaymard F, Tillard P, Lepetit M, Lejay L, Gojon A: Transcript profiling in the chl1–5 mutant of Arabidopsis reveals a role of the nitrate transporter NRT1.1 in the regulation of another nitrate transporter, NRT2.1.

    Plant Cell 2004, 16(9):2433-2447. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  23. Walch-Liu P, Forde BG: Nitrate signalling mediated by the NRT1.1 nitrate transporter antagonises L-glutamate-induced changes in root architecture.

    Plant J 2008, 54(5):820-828. PubMed Abstract | Publisher Full Text OpenURL

  24. Hartigan J: Clustering algorithms. New York, NY, USA: John Wiley & Sons; 1975.

  25. Tanay ASR, Shamir R: Discovering statistically significant Biclusters that are associated with replication origins.

    Curr Biology 1999, 9:R661-R667. Publisher Full Text OpenURL

  26. Wang R, Okamoto M, Xing X, Crawford NM: Microarray analysis of the nitrate response in Arabidopsis roots and shoots reveals over 1,000 rapidly responding genes and new linkages to glucose, trehalose-6-phosphate, iron, and sulfate metabolism.

    Plant Physiol 2003, 132(2):556-567. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Wang R, Tischner R, Gutierrez RA, Hoffman M, Xing X, Chen M, Coruzzi G, Crawford NM: Genomic analysis of the nitrate response using a nitrate reductase-null mutant of Arabidopsis.

    Plant Physiol 2004, 136(1):2512-2522. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Scheible W-R, Morcuende R, Czechowski T, Fritz C, Osuna D, Palacios-Rojas N, Schindelasch D, Thimm O, Udvardi MK, Stitt M: Genome-wide reprogramming of primary and secondary metabolism, protein synthesis, cellular growth processes, and the regulatory infrastructure of Arabidopsis in response to nitrogen.

    Plant Physiol 2004, 136(1):2483-2499. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Palenchar PM, Kouranov A, Lejay LV, Coruzzi GM: Genome-wide patterns of carbon and nitrogen regulation of gene expression validate the combined carbon and nitrogen (CN)-signaling hypothesis in plants.

    Genome Biol 2004, 5(11):R91. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Gutierrez RA, Lejay LV, Dean A, Chiaromonte F, Shasha DE, Coruzzi GM: Qualitative network models and genome-wide expression data define carbon/nitrogen-responsive molecular machines in Arabidopsis.

    Genome Biol 2007, 8(1):R7. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  31. Sharan R: M-KAaSR: CLICK and EXPANDER: A System for Clustering and Visualizing Gene Expression Data.

    Bioinformatics 2003, 19(14):1787-1799. PubMed Abstract | Publisher Full Text OpenURL

  32. Tanay A, Sharan R, Kupiec M, Shamir R: Revealing modularity and organization in the yeast molecular network by integrated analysis of highly heterogeneous genomewide data.

    Proc Natl Acad Sci USA 2004, 101(9):2981-2986. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  33. Schoof H, Zaccaria P, Gundlach H, Lemcke K, Rudd S, Kolesov G, Arnold R, Mewes HW, Mayer KF: MIPS Arabidopsis thaliana Database (MAtDB): an integrated biological knowledge resource based on the first complete plant genome.

    Nucleic Acids Res 2002, 30(1):91-93. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Madar A, Bonneau R: Learning global models of transcriptional regulatory networks from data.

    Methods in molecular biology (Clifton, NJ) 2009, 541:181. OpenURL

  35. Gan RR, Yi EC, Chiu Y, Lee H, Kao YC, Wu TH, Aebersold R, Goodlett DR, Ng WV: Proteome analysis of Halobacterium sp. NRC-1 facilitated by the biomodule analysis tool BMSorter.

    Mol Cell Proteomics 2006, 5(6):987-997. PubMed Abstract | Publisher Full Text OpenURL

  36. Bonneau R: Learning biological networks: from modules to dynamics.

    Nat Chem Biol 2008, 4(11):658-664. PubMed Abstract | Publisher Full Text OpenURL

  37. Thum KE, Shin MJ, Gutierrez RA, Mukherjee I, Katari MS, Nero D, Shasha D, Coruzzi GM: An integrated genetic, genomic and systems approach defines gene networks regulated by the interaction of light and carbon signaling pathways in Arabidopsis.

    BMC systems biology 2008, 2:31. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  38. Cerezo M, Tillard P, Filleur S, Munos S, Daniel-Vedele F, Gojon A: Major alterations of the regulation of root NO3- uptake are associated with the mutation of Nrt2.1 and Nrt2.2 genes in Arabidopsis.

    Plant Physiol 2001, 127(1):262-271. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  39. Gutierrez RA, Stokes TL, Thum K, Xu X, Obertello M, Katari MS, Tanurdzic M, Dean A, Nero DC, McClung CR, et al.: Systems approach identifies an organic nitrogen-responsive gene network that is regulated by the master clock control gene CCA1.

    Proc Natl Acad Sci USA 2008, 105(12):4939-4944. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  40. Ruffel S, Freixes S, Balzergue S, Tillard P, Jeudy C, Martin-Magniette ML, Merwe MJ, Kakar K, Gouzy J, Fernie AR, et al.: Systemic signaling of the plant nitrogen status triggers specific transcriptome responses depending on the nitrogen source in Medicago truncatula.

    Plant Physiol 2008, 146(4):2020-2035. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  41. Zhuo D, Okamoto M, Vidmar JJ, Glass AD: Regulation of a putative high-affinity nitrate transporter (Nrt2;1At) in roots of Arabidopsis thaliana.

    Plant J 1999, 17(5):563-568. PubMed Abstract | Publisher Full Text OpenURL

  42. Brenner WG, Romanov GA, Kollmer I, Burkle L, Schmulling T: Immediate-early and delayed cytokinin response genes of Arabidopsis thaliana identified by genome-wide expression profiling reveal novel cytokinin-sensitive processes and suggest cytokinin action through transcriptional cascades.

    Plant J 2005, 44(2):314-333. PubMed Abstract | Publisher Full Text OpenURL

  43. Gan Y, Filleur S, Rahman A, Gotensparre S, Forde BG: Nutritional regulation of ANR1 and other root-expressed MADS-box genes in Arabidopsis thaliana.

    Planta 2005, 222(4):730-742. PubMed Abstract | Publisher Full Text OpenURL

  44. 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):D575-577. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  45. Team RDC: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria; 2007.

  46. Crawley MJ: Statistics: an introduction using R. Volume 1. 1st edition. New York: John Wiley & Sons; 2005.

  47. van Helden JAB, Collado-Vides J: A web site for the computational analysis of yeast regulatory sequences.

    Yeast 2000, 16(2):177-187. PubMed Abstract | Publisher Full Text OpenURL

  48. Littlefield O, Nelson HC: A new use for the 'wing' of the 'winged' helix-turn-helix motif in the HSF-DNA cocrystal.

    Nat Struct Biol 1999, 6(5):464-470. PubMed Abstract | Publisher Full Text OpenURL

  49. Vandepoele K, Vlieghe K, Florquin K, Hennig L, Beemster GT, Gruissem W, Peer Y, Inze D, De Veylder L: Genome-wide identification of potential plant E2F target genes.

    Plant Physiol 2005, 139(1):316-328. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  50. Baliga NS, Pan M, Goo YA, Yi EC, Goodlett DR, Dimitrov K, Shannon P, Aebersold R, Ng WV, Hood L: Coordinate regulation of energy transduction modules in Halobacterium sp. analyzed by a global systems approach.

    Proc Natl Acad Sci USA 2002, 99(23):14913-14918. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  51. Bonneau R, Facciotti MT, Reiss DJ, Schmid AK, Pan M, Kaur A, Thorsson V, Shannon P, Johnson MH, Bare JC, et al.: A predictive model for transcriptional control of physiology in a free living cell.

    Cell 2007, 131(7):1354-1365. PubMed Abstract | Publisher Full Text OpenURL

  52. Franco-Zorrilla JM, Martin AC, Leyva A, Paz-Ares J: Interaction between phosphate-starvation, sugar, and cytokinin signaling in Arabidopsis and the roles of cytokinin receptors CRE1/AHK4 and AHK3.

    Plant Physiol 2005, 138(2):847-857. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  53. Franco-Zorrilla JM, Martin AC, Solano R, Rubio V, Leyva A, Paz-Ares J: Mutations at CRE1 impair cytokinin-induced repression of phosphate starvation responses in Arabidopsis.

    Plant J 2002, 32(3):353-360. PubMed Abstract | Publisher Full Text OpenURL

  54. Moore B, Zhou L, Rolland F, Hall Q, Cheng WH, Liu YX, Hwang I, Jones T, Sheen J: Role of the Arabidopsis glucose sensor HXK1 in nutrient, light, and hormonal signaling.

    Science 2003, 300(5617):332-336. PubMed Abstract | Publisher Full Text OpenURL

  55. Maruyama-Nakashita A, Nakamura Y, Yamaya T, Takahashi H: A novel regulatory pathway of sulfate uptake in Arabidopsis roots: implication of CRE1/WOL/AHK4-mediated cytokinin-dependent regulation.

    Plant J 2004, 38(5):779-789. PubMed Abstract | Publisher Full Text OpenURL

  56. Seguela M, Briat JF, Vert G, Curie C: Cytokinins negatively regulate the root iron uptake machinery in Arabidopsis through a growth-dependent pathway.

    Plant J 2008, 55(2):289-300. PubMed Abstract | Publisher Full Text OpenURL

  57. Nover L, Bharti K, Doring P, Mishra SK, Ganguli A, Scharf KD: Arabidopsis and the heat stress transcription factor world: how many heat stress transcription factors do we need?

    Cell stress & chaperones 2001, 6(3):177-189. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  58. Kotak S, Larkindale J, Lee U, von Koskull-Doring P, Vierling E, Scharf KD: Complexity of the heat stress response in plants.

    Curr Opin Plant Biol 2007, 10(3):310-316. PubMed Abstract | Publisher Full Text OpenURL

  59. Poultney CS, Gutierrez RA, Katari MS, Gifford ML, Paley WB, Coruzzi GM, Shasha DE: Sungear: interactive visualization and functional analysis of genomic datasets.

    Bioinformatics 2007, 23(2):259-261. PubMed Abstract | Publisher Full Text OpenURL

  60. Nemhauser JL, Hong F, Chory J: Different plant hormones regulate similar processes through largely nonoverlapping transcriptional responses.

    Cell 2006, 126(3):467-475. PubMed Abstract | Publisher Full Text OpenURL

  61. Inze D: Green light for the cell cycle.

    Embo J 2005, 24(4):657-662. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  62. Blais A, Dynlacht BD: E2F-associated chromatin modifiers and cell cycle control.

    Current opinion in cell biology 2007, 19(6):658-662. PubMed Abstract | Publisher Full Text OpenURL

  63. Vlieghe K, Vuylsteke M, Florquin K, Rombauts S, Maes S, Ormenese S, Van Hummelen P, Peer Y, Inze D, De Veylder L: Microarray analysis of E2Fa-DPa-overexpressing plants uncovers a cross-talking genetic network between DNA replication and nitrogen assimilation.

    J Cell Sci 2003, 116(20):4249-4259. PubMed Abstract | Publisher Full Text OpenURL

  64. Cannons AC, Shiflett SD: Transcriptional regulation of the nitrate reductase gene in Chlorella vulgaris: identification of regulatory elements controlling expression.

    Current genetics 2001, 40(2):128-135. PubMed Abstract | Publisher Full Text OpenURL

  65. del Pozo JC, Diaz-Trivino S, Cisneros N, Gutierrez C: The balance between cell division and endoreplication depends on E2FC-DPB, transcription factors regulated by the ubiquitin-SCFSKP2A pathway in Arabidopsis.

    Plant Cell 2006, 18(9):2224-2235. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  66. Shannon PMA, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: A software environment for integrated models of biomolecular interaction networks.

    Genome Research 2003, 13:2498-2404. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  67. Storey JD, Tibshirani R: Statistical significance for genomewide studies.

    Proc Natl Acad Sci USA 2003, 100(16):9440-9445. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  68. Davuluri RV, Sun H, Palaniswamy SK, Matthews N, Molina C, Kurtz M, Grotewold E: AGRIS: Arabidopsis gene regulatory information server, an information resource of Arabidopsis cis-regulatory elements and transcription factors.

    BMC bioinformatics 2003, 4:25. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL