A system biology approach highlights a hormonal enhancer effect on regulation of genes in a nitrate responsive "biomodule"
© Nero et al; licensee BioMed Central Ltd. 2009
Received: 26 March 2009
Accepted: 06 June 2009
Published: 06 June 2009
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.
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.
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.
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 . 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 . Moreover, it has also been shown that such cross-controls exist between NO3- and: cytokinin (for review see ), auxin [2, 6, 7], and ABA . 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 . 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  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 . 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 . 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 .
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 . 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 . 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  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) , 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  (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].
Over-represented functional categories in N-bicluster 9
Observed Frequency %
Expected Frequency %
metabolism of energy reserves (e.g. glycogen, trehalose)
glycolysis and gluconeogenesis
C-compound and carbohydrate metabolism
amino acid metabolism
assimilation of ammonia, metabolism of the glutamate group
nitrogen and sulfur metabolism
TRANSPORTED COMPOUNDS (SUBSTRATES)
The genes in the N-bicluster 9 "biomodule" are induced by nitrate as a signal
Identification of nitrate/hormone responsive biclusters
NO3-/Hormone Response Interaction.
Description of Comparison
N/H-bicluster 1 vs. Exclusive +N/H-bicluster 6
N/H-Bicluster 16, 19, 20 vs. Exclusive + N/H-bicluster 6
Treated vs. Control
N/H-Bicluster 16, 19, 20: Treated
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
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] . 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  (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  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.
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
Over-represented known CREs in N-bicluster 9.
Cis Regulatory Element
% Present in Whole Genome
% Present in N-bicluster 9
Bellringer/replumless/pennywise BS1 IN AG
ATB2/AtbZIP53/AtbZIP44/GBF5BS in ProDH
AtMYB2 BS in RD22
Over-represented CREs for 4 significant N/H-biclusters.
P-value (N/H-Bicluster vs. N-bicluster 9 Exclusive)
% Present in N/H-bicluster Genes
% Present in N-bicluster 9 Exclusive Genes
% Present in N-bicluster 9 Genes
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 . 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.  state that "a biomodule is a group of proteins that execute a particular function", and ii) Bonneau et al.  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 . 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 , sulfate , and iron metabolism . 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 . However, a heat shock transcription factor HsfA9 has been shown to be under hormonal control in seeds (ABA through ABI3) . 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  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. . 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.  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 . 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 .
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.
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  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 ) 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  and as implemented in the CLICK and EXPANDER program . 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  as accessed via http://www.virtualplant.org. The program was run using the MIPS  annotation option for functional definitions. A 5% FDR (False Discovery Rate) cut-off was computed using the R statistical package  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  functional terms were determined using the Fisher's Exact Test as implemented in the R statistical package  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  as accessed by http://www.virtualplant.org. This network contains many validated connections for gene interactions in the Arabidopsis genome. Network interactions were visualized using Cytoscape .
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 . This script was derived from the Storey and Tibshirani method  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 ) 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  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.
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.
- Crawford NM, Glass ADM: Molecular and physiological aspects of nitrate uptake in plants. Trends in Plant Science. 1998, 3 (10): 389-395. 10.1016/S1360-1385(98)01311-9.View ArticleGoogle Scholar
- 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. 10.1093/jxb/erm079View ArticlePubMedGoogle Scholar
- 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- 10.1371/journal.pcbi.1000326PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1007/s11103-008-9380-yView ArticlePubMedGoogle Scholar
- 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. 10.1016/j.tplants.2006.07.004View ArticlePubMedGoogle Scholar
- Avery GS, Pottorf L: Auxin and nitrogen relationshinps in green plants. American Journal of Botany. 1945, 32: 666-669. 10.2307/2437622.View ArticleGoogle Scholar
- 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. 10.1093/jexbot/53.370.835View ArticlePubMedGoogle Scholar
- 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. 10.1046/j.1365-313x.2001.01185.xView ArticlePubMedGoogle Scholar
- 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. 10.1093/pcp/pch119View ArticlePubMedGoogle Scholar
- 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. 10.1046/j.1365-313X.2003.01945.xView ArticlePubMedGoogle Scholar
- Sakakibara H: Nitrate-specific and cytokinin-mediated nitrogen signaling pathways in plants. J Plant Res. 2003, 116 (3): 253-257. 10.1007/s10265-003-0097-3View ArticlePubMedGoogle Scholar
- 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. 10.1093/jexbot/53.370.971View ArticlePubMedGoogle Scholar
- 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. 10.1093/pcp/pce009View ArticlePubMedGoogle Scholar
- 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. 10.1093/jxb/eri107View ArticlePubMedGoogle Scholar
- 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. 10.1111/j.1469-8137.1975.tb01409.x.View ArticleGoogle Scholar
- 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. 10.1007/s11120-004-2080-9View ArticlePubMedGoogle Scholar
- 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. 10.1093/aob/mcj601View ArticleGoogle Scholar
- 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. 10.1073/pnas.96.11.6529PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1126/science.1122847View ArticlePubMedGoogle Scholar
- 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. 10.1073/pnas.0605275103PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1104/pp.106.087510PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1105/tpc.104.024380PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1111/j.1365-313X.2008.03443.xView ArticlePubMedGoogle Scholar
- Hartigan J: Clustering algorithms. 1975, New York, NY, USA: John Wiley & SonsGoogle Scholar
- Tanay ASR, Shamir R: Discovering statistically significant Biclusters that are associated with replication origins. Curr Biology. 1999, 9: R661-R667. 10.1016/S0960-9822(99)80418-7.View ArticleGoogle Scholar
- 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. 10.1104/pp.103.021253PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1104/pp.104.044610PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1104/pp.104.047019PubMed CentralView ArticlePubMedGoogle Scholar
- 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 CentralView ArticlePubMedGoogle Scholar
- 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- 10.1186/gb-2007-8-1-r7PubMed CentralView ArticlePubMedGoogle Scholar
- Sharan R: M-KAaSR: CLICK and EXPANDER: A System for Clustering and Visualizing Gene Expression Data. Bioinformatics. 2003, 19 (14): 1787-1799. 10.1093/bioinformatics/btg232View ArticlePubMedGoogle Scholar
- 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. 10.1073/pnas.0308661100PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1093/nar/30.1.91PubMed CentralView ArticlePubMedGoogle Scholar
- Madar A, Bonneau R: Learning global models of transcriptional regulatory networks from data. Methods in molecular biology (Clifton, NJ). 2009, 541: 181-View ArticleGoogle Scholar
- 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. 10.1074/mcp.M500367-MCP200View ArticlePubMedGoogle Scholar
- Bonneau R: Learning biological networks: from modules to dynamics. Nat Chem Biol. 2008, 4 (11): 658-664. 10.1038/nchembio.122View ArticlePubMedGoogle Scholar
- 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- 10.1186/1752-0509-2-31PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1104/pp.127.1.262PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1073/pnas.0800211105PubMed CentralView ArticlePubMedGoogle Scholar
- Ruffel S, Freixes S, Balzergue S, Tillard P, Jeudy C, Martin-Magniette ML, Merwe van der 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. 10.1104/pp.107.115667PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1046/j.1365-313X.1999.00396.xView ArticlePubMedGoogle Scholar
- 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. 10.1111/j.1365-313X.2005.02530.xView ArticlePubMedGoogle Scholar
- 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. 10.1007/s00425-005-0020-3View ArticlePubMedGoogle Scholar
- 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, D575-577. 32 Database
- Team RDC: A language and environment for statistical computing. 2007, R Foundation for Statistical Computing, Vienna, AustriaGoogle Scholar
- Crawley MJ: Statistics: an introduction using R. 2005, 1: New York: John Wiley & Sons, 1View ArticleGoogle Scholar
- van Helden JAB, Collado-Vides J: A web site for the computational analysis of yeast regulatory sequences. Yeast. 2000, 16 (2): 177-187. 10.1002/(SICI)1097-0061(20000130)16:2<177::AID-YEA516>3.0.CO;2-9View ArticlePubMedGoogle Scholar
- 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. 10.1038/8269View ArticlePubMedGoogle Scholar
- Vandepoele K, Vlieghe K, Florquin K, Hennig L, Beemster GT, Gruissem W, Peer Van de Y, Inze D, De Veylder L: Genome-wide identification of potential plant E2F target genes. Plant Physiol. 2005, 139 (1): 316-328. 10.1104/pp.105.066290PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1073/pnas.192558999PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1016/j.cell.2007.10.053View ArticlePubMedGoogle Scholar
- 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. 10.1104/pp.105.060517PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1046/j.1365-313X.2002.01431.xView ArticlePubMedGoogle Scholar
- 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. 10.1126/science.1080585View ArticlePubMedGoogle Scholar
- 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. 10.1111/j.1365-313X.2004.02079.xView ArticlePubMedGoogle Scholar
- 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. 10.1111/j.1365-313X.2008.03502.xView ArticlePubMedGoogle Scholar
- 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. 10.1379/1466-1268(2001)006<0177:AATHST>2.0.CO;2View ArticleGoogle Scholar
- 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. 10.1016/j.pbi.2007.04.011View ArticlePubMedGoogle Scholar
- 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. 10.1093/bioinformatics/btl496View ArticlePubMedGoogle Scholar
- Nemhauser JL, Hong F, Chory J: Different plant hormones regulate similar processes through largely nonoverlapping transcriptional responses. Cell. 2006, 126 (3): 467-475. 10.1016/j.cell.2006.05.050View ArticlePubMedGoogle Scholar
- Inze D: Green light for the cell cycle. Embo J. 2005, 24 (4): 657-662. 10.1038/sj.emboj.7600561PubMed CentralView ArticlePubMedGoogle Scholar
- Blais A, Dynlacht BD: E2F-associated chromatin modifiers and cell cycle control. Current opinion in cell biology. 2007, 19 (6): 658-662. 10.1016/j.ceb.2007.10.003PubMed CentralView ArticlePubMedGoogle Scholar
- Vlieghe K, Vuylsteke M, Florquin K, Rombauts S, Maes S, Ormenese S, Van Hummelen P, Peer Van de 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. 10.1242/jcs.00715View ArticlePubMedGoogle Scholar
- 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. 10.1007/s002940100232View ArticlePubMedGoogle Scholar
- 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. 10.1105/tpc.105.039651PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1101/gr.1239303PubMed CentralView ArticlePubMedGoogle Scholar
- Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proc Natl Acad Sci USA. 2003, 100 (16): 9440-9445. 10.1073/pnas.1530509100PubMed CentralView ArticlePubMedGoogle Scholar
- 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- 10.1186/1471-2105-4-25PubMed CentralView ArticlePubMedGoogle Scholar