- Research article
- Open Access
Global transcription profiling reveals differential responses to chronic nitrogen stress and putative nitrogen regulatory components in Arabidopsis
BMC Genomics volume 8, Article number: 281 (2007)
A large quantity of nitrogen (N) fertilizer is used for crop production to achieve high yields at a significant economic and environmental cost. Efforts have been directed to understanding the molecular basis of plant responses to N and identifying N-responsive genes in order to manipulate their expression, thus enabling plants to use N more efficiently. No studies have yet delineated these responses at the transcriptional level when plants are grown under chronic N stress and the understanding of regulatory elements involved in N response is very limited.
To further our understanding of the response of plants to varying N levels, a growth system was developed where N was the growth-limiting factor. An Arabidopsis whole genome microarray was used to evaluate global gene expression under different N conditions. Differentially expressed genes under mild or severe chronic N stress were identified. Mild N stress triggered only a small set of genes significantly different at the transcriptional level, which are largely involved in various stress responses. Plant responses were much more pronounced under severe N stress, involving a large number of genes in many different biological processes. Differentially expressed genes were also identified in response to short- and long-term N availability increases. Putative N regulatory elements were determined along with several previously known motifs involved in the responses to N and carbon availability as well as plant stress.
Differentially expressed genes identified provide additional insights into the coordination of the complex N responses of plants and the components of the N response mechanism. Putative N regulatory elements were identified to reveal possible new components of the regulatory network for plant N responses. A better understanding of the complex regulatory network for plant N responses will help lead to strategies to improve N use efficiency.
Nitrogen (N) is the most important inorganic nutrient for plant growth . It affects many aspects of plant growth and development, such as N and carbon (C) allocation, root branching, leaf growth and flowering time [2–5]. The production of high-yielding crops is associated with the application of a large quantity of fertilizers at a substantial cost , and N pollution is becoming a threat to global ecosystems [7, 8]. Efforts have been directed to understanding the molecular basis of plant responses to N and to identifying N-responsive genes in order to manipulate their expression and enable plants to use N more efficiently .
Nitrate is the major source of N in agricultural soils . It serves as a nutrient and as a signal [11, 12]. As a nutrient, it is taken up by the low and high affinity nitrate transporter gene family members (NRT1 and NRT2), reduced to nitrite by nitrate reductase (NR), and to ammonium by nitrite reductase (NiR). Ammonium is then incorporated into amino acids, catalyzed primarily by glutamine synthetase (GS) and glutamate synthase (GOGAT) [11–14]. As a signal, it can induce the expression of a number of genes including NRT1, NRT2, NR and NiR [4, 12, 15, 16]. The expression of the ammonium assimilatory genes, GS and GOGAT, are also induced by the addition of nitrate [12, 13], as are genes involved in sugar metabolism and other metabolic pathways . In addition to these metabolic genes, expression of some regulatory genes is also affected by N levels. For example, the Arabidopsis MADS-box gene ANR1, which regulates lateral root development, responds to N, but not potassium or phosphate .
There have been several studies of plant N-responses based on microarray gene expression profiling. Wang et al  studied the response of seedlings grown on ammonium to the addition of low or high levels of nitrate. They used the Arabidopsis GEM1 microarrays, which contained 7942 cDNA clones corresponding to 5524 unique genes, and identified 25 and 49 N-responsive genes to low or high nitrate induction, respectively . Subsequently, Wang et al  used the Arabidopsis whole-genome Affymetrix ATH1 microarray containing 22,626 genes, to study the addition of the low level of nitrate to discover more N-responsive genes. Scheible et al  also used the ATH1 microarray combined with real-time RT-PCR of > 1,400 transcription factor genes to identify genes affected by N-deprivation or N-induction after 30 min or 3 hr from N-starved seedlings. Since C and N metabolism are very closely linked and tightly regulated [21, 22], Price et al  used the ATH1 microarray to identify the individual contributions of nitrogen, sugar, and nitrogen plus sugar on global gene expression. Recently, Lian et al reported expression profiles of 10,422 genes at an early stage of low N stress in rice seedling .
So far, these studies have provided valuable insights into N response and its linkage to other biological pathways. However, no studies have yet delineated the responses at the transcriptional level when plants are grown under chronic N stress. Earlier reports involved investigations of transient changes in gene expression when nitrate is added to nitrate-starved seedlings [18–20, 23]. As such, many questions still remain unanswered. For example, it is not clear how mature plants would respond to different degrees of chronic N stress. A growth system was developed where N was the growth-limiting factor, so the transcriptional changes of genes that were most affected by different degrees of N limitation could be investigated. Additionally, the global transient changes during N induction were examined by transferring mature plants from low N to high N instead of using N-starved seedlings. In addition, putative regulatory elements involved in N response were also identified as an important first step toward understanding N regulatory networks. Previous progress in this area is limited to nitrate induction from deletion analysis of NR or NiR promoters [25–29], and a study on the interaction between C and N signaling .
Results and discussion
Developing defined nitrogen growth conditions for expression profiling
To apply chronic N stress, it is important to develop a defined nitrogen growth condition. It is difficult to maintain a constant N level in a soil system due to the different size and affinity of soil particles for nutrients. A "pure" hydroponic system could control the N level well but the root system is not supported by a substrate and it is bathed directly in the nutrient solution with poor aeration. To overcome those shortcomings, a hydroponic system using rockwool as the growing substrate was adopted . The key quality of high grade rockwool is uniform wetting. It provides a buffering reservoir of nutrient solution in the root zone, while maintaining an adequate volume of air (oxygen). A detailed description of the system is provided in the Methods section. Wild-type Arabidopsis (Columbia ecotype) plants at 3 weeks of age were evaluated for shoot biomass under different N conditions, ranging from 0.1 mM to 10 mM. Under these N conditions, 3 mM nitrate was found to give maximal growth (data not shown). For subsequent experiments, 3 mM nitrate was used to produce the N-sufficient condition, 1 mM nitrate to produce the mild N-limiting condition, under which plant growth measured by shoot biomass was reduced to approximately 80% of that at 3 mM nitrate, and 0.3 mM nitrate to produce the severe N-limiting condition, under which plant growth measured by shoot biomass was further reduced to approximately 35% of that at 3 mM nitrate (Figure 1a). Under these conditions, the nitrate levels in the leaves ranged from 0.24 mg/g FW to 1.64 mg/g FW (Figure 1b). Plants have a nitrate transport system that could have a very high affinity for nitrate  and provide a certain capacity for nitrate absorption at low external nitrate concentrations. However, in our system, the amount of nitrate in one rockwool cube that was available to support the growth of one plant was apparently very limited under the low N conditions. The consequence was the obvious difference in growth under these N conditions. Apart from reduced shoot biomass due to N deficiency, a reduction in chlorophyll level was also observed (data not shown).
Identification of differentially expressed genes by expression profiling
Shoots of 3-week-old wild-type plants grown under the N-sufficient condition (3 mM nitrate), the mild N-limiting condition (1 mM nitrate) and the severe N-limiting condition (0.3 mM nitrate) were harvested for the profiling experiment to compare the baseline gene expression levels under different, but stable N conditions. Since it is known that the expression levels of nitrate assimilation genes are high early in the day and low late in the day, all samples were taken at the middle of the day to minimize diurnal changes in C and N metabolism . Each condition had three biological replicates. Also, at 2 hrs or 24 hrs before the harvest, some of the plants grown under the severe N-limiting condition (0.3 mM nitrate) were transferred to the N-sufficient condition (3 mM nitrate), and then collected at the same time as those grown under different stable N conditions to assess gene expression changes after 2 hr or 24 hr N induction. Nitrate levels in the leaves increased from 0.24 mg/g FW to 0.49 mg/g FW 2 hr after addition of higher levels of nitrate and to 1.09 mg/g FW after 24 hr (Figure 1b). Again, three biological repeats were collected for each time point. RNA was extracted and hybridized to a custom designed Arabidopsis GeneChip® whole genome microarray. This microarray contains 26,412 probe sets representing 26,367 Arabidopsis known and predicted genes. On average, each gene contains 15 perfect match probes, selected from the 3' end of 133,397 exons .
Significance analysis of microarray (SAM; academic version 2.23B)  was used to identify genes differentially expressed between the following: (1) the N-sufficient to the mild N-limiting samples (3 mM N to 1 mM N); (2) the N-sufficient to the severe N-limiting samples (3 mM N to 0.3 mM N); (3) severe N (0.3 mM N) to the 2 hr N induction samples; and (4) severe N (0.3 mM N) to the 24 hr N induction samples. Different numbers of significant genes could be obtained depending on the SAM median false discovery rate (FDR). The final number of significant genes of the following analyses was based on a median FDR of 0% and a minimum of a 2-fold expression difference (Table 1). To confirm the results of the microarray analysis, the expression of 14 significant genes (7 up-regulated and 7 down-regulated) was tested by semi-quantitative RT-PCR and found to be consistent with the microarray data [see Additional file 1].
Functional classification of significant genes under mild and severe chronic N stress
Although plant growth measured by shoot biomass was reduced by 20% under mild N stress, the baseline expression levels of most genes in the genome remained similar. Only 52 genes were identified with expression levels significantly different, of which 51 were induced and one suppressed in the mild N-limiting condition compared to the N-sufficient condition. However, when the degree of this chronic N stress increased, these numbers increased dramatically. A total of 461 genes showed significant response, of which 271 genes were up-regulated and 190 genes down-regulated under the severe N-limiting condition, when plant growth was reduced by 65% measured by shoot biomass. Severe N stress introduced a much greater transcriptional changes compared to mild N stress. This result is similar to the microarray study on plants under different degrees of phosphate deficiency . In that case, 72 genes were induced and four genes were suppressed during short-term Pi deprivation, and the numbers increased to 291 genes induced and 34 genes suppressed during medium-term Pi starvation, and further increased to 501 genes induced and 231 genes suppressed during long-term deficiency .
The entire gene lists under the mild and severe N-limiting conditions are provided in Additional file 2 and 3 respectively, together with their functional and pathway assignments. Functional assignments are defined by Gene Ontology (GO) terms , which provide broad functional classifications for genes and gene products representing their corresponding biological processes, molecular function, and cellular localization. Pathway assignments are derived from the Kyoto Encyclopedia of Genes and Genomes (KEGG) . Gene ontology on the Arabidopsis Information Resource website , and a browser-based functional classification program  were also consulted.
Under the mild N-limiting condition when plant growth was slightly reduced, none of the genes directly involved in nitrate assimilation showed significant changes in expression (Table 2). However, when plants were grown under the severe N-limiting condition, many primary metabolism genes including those involved in N assimilation were significantly down-regulated. In terms of the nitrate assimilatory genes, the baseline expressions of NR1, NR2 and NiR were clearly correlated with the three N levels – the lower the N concentration, the lower the baseline expression (Table 2). This indicates that the degree of the reduction of N assimilation gene expression is correlated to the amount of N available. As for the ammonium assimilatory genes, GS1-1 and GS1-4, the opposite trend of the nitrate assimilation genes (NR1, NR2 and NiR) was observed. GS1-1 and GS1-4 were up-regulated under severe N stress (Table 2). GS has a very high affinity for ammonia and catalyzes the ATP-dependent conversion of glutamate into glutamine by incorporating a molecule of ammonia. The individual isoenzymes of GS and GOGAT have been proposed to play roles in three major ammonium assimilation processes: (1) primary nitrogen assimilation; (2) re-assimilation of photorespiratory ammonia; and (3) re-assimilation of "recycled" nitrogen . Since primary metabolism such as nitrogen assimilation, photosynthesis and photorespiration is repressed under severe N stress, it is likely that the increased expression of cytosolic GS1 is related to the re-assimilation of ammonia released from protein degradation which usually occurs under N deficiency. With regard to the N uptake genes, different trends were observed among members of the low affinity nitrate transporter gene family NRT1. NRT1.1 baseline expression was higher, while NRT1.3 and NRT1.4 baselines expression were lower under N stress (Table 2), suggesting a lack of functional redundancy among family members. In addition to the N assimilation genes, those involved in the nitrate storage were deeply repressed by N deficiency. An anion channel protein gene, which was recently demonstrated to mediate nitrate accumulation in plant vacuoles , was down-regulated almost 7-fold under severe N shortage while no significant difference in expression was observed under mild N limitation (Table 3).
N stress led to marked changes in the expression of genes involved in carbon metabolism. Genes involved in photosynthesis, including those that code for the photosystem II family protein, the oxygen evolving enhancer 3 PsbQ family protein, the PSI type II chlorophyll a b-binding protein, and the photosystem II reaction centre W PsbW family protein were significantly down-regulated under the severe N-limiting condition, but not the mild N-limiting condition (Table 3). In addition, genes involved in the oxidative pentose-phosphate pathway that provide reducing power (NADPH) and pentose phosphates were down-regulated, including the gene coding for glucose-6-phosphate 1-dehydrogenase, which is the rate-limiting step in this process (Table 3). On the other hand, genes involved in the accumulation of starch, including ADP-glucose pyrophosphorylase which catalyzes the first, rate-limiting step in starch biosynthesis, were significantly up-regulated under severe N stress (Table 3).
N stress induces profound changes to growth and development. In order to adapt to the growth arrest due to severe N deficiency, many ribosomal genes involved in protein biosynthesis were down-regulated, including genes coding for structural constituents of ribosome and translation elongation factors [see Additional file 3]. At the same time, genes involved in protein degradation were up-regulated under severe N stress, such as the Cys peptidase (SAG12) (Table 3). It is believed that during nutrient deficiency, plants transport a series of cytosolic proteins into the vacuole, where various proteases are located, to be degraded into amino acids and exported from senescing tissues for cell reuse.
As expected, there were marked changes in the expression of genes involved in hormone synthesis and sensing since they largely control plant growth. Auxin stimulates cell devision and elongation, and many auxin-induced genes were down-regulated under severe N stress [see Additional file 3]. Cytokinin regulates cell proliferation and differentiation, and many cytokinin-responsive genes were down-regulated under severe N stress, including two transcription regulators (At3g48100, At2g18300) (Table 3).
Leaf yellowing is one of the typical responses plants have when N deficiency occurs . While chlorophyll levels were slightly reduced (~5% reduction) under the mild N- limiting condition, none of the chlorophyll biosynthesis genes changed appreciably at the transcriptional level, indicating a disconnection between transcriptional responses and phenotypic changes. Under the severe N-limiting condition, the chlorophyll levels were significantly reduced (~30% reduction). Genes involved in chlorophyll metabolism were down-regulated, including protochlorophyllide reductase, protochlorophyllide reductase and hydroxymethylbilane synthase (Table 3).
N stress induced a number of plant stress responses. Many genes involved in various stress responses were up-regulated, including some peroxidase genes and glutathioneS-transfeerase (GST) genes under both N-limiting conditions [see Additional file 2 and 3]. Toxic nitrogen compounds are usually generated under various stresses. Nitrilase 4, of which the purified enzyme has been shown to be involved in the nitrogen compound detoxification pathway, was up-regulated around 3-fold under mild N stress, but up-regulated over 11-fold under severe N stress (Table 3).
The presence of the purple flavonoid anthocyanin is an indicator of stress , and N deficiency may lead to increased anthocyanin synthesis . Plants were obviously purple in color under the severe N-limiting condition (not shown). Genes involved in anthocyanin synthesis, such as leucoanthocyanidin dioxygenase and dihydroflavonol reductase, were up-regulated just over 2-fold under mild N stress, but increased to about 14- and 18-fold under severe N stress (Table 3). Chalcone synthase, which participates in the biosynthetic pathway for all flavonoids and is required for the accumulation of anthocyanins, was up-regulated only under severe N stress (Table 3). The expression of the Myb transcription factor PAP1 was increased about 6-fold under severe N stress, and the expression of the Myb transcription factor PAP2 was increased around 3-fold under the mild N-limiting condition, but dramatically increased to over 26-fold under the severe N-limiting condition (Table 3). The two Myb genes have been shown to be able to regulate flavonoid and anthocyanin biosynthesis and over-expression of the Myb genes resulted in an enhanced accumulation of lignin and flavonoids, including various anthocyanins that produce purple color . It is clear that the accumulation of anthocyanins is closely correlated with the level of N stress, and PAP1 and PAP2 are involved in the regulation of this process. However, it is unclear how specific this regulation is for N stress response.
As N deficiency responses are known to be regulated at the transcriptional level, transcriptional factor genes were looked up from the entire gene lists [see Additional file 2 and 3]. Three transcription factors genes were induced under mild N limitation, including a WRKY transcription factor (At1g80840), an ethylene responsive element binding factor 2 (ATERF-2, At5g47220) and the aforementioned Myb transcription factor (PAP2) gene [see Additional file 2]. Plant proteins containing WRKY or ATERF domains are known regulators of abiotic and biotic stress responses [46, 47]. Detailed characterization of the WRKY gene (At1g80840) is lacking, but the ATERF-2 gene (At5g47220) is known to be regulated by other abiotic stress conditions . Thirty-nine transcription factor genes were significantly induced or repressed under severe N stress, including the Myb factors PAP1 and PAP2, the two cytokinin-responsive transcription regulators (At3g48100, At2g18300) mentioned earlier, and three genes (At1g30500, At2g34720 and At3g05690) coding for CCAAT binding factor complex (Table 3). While members in this CCAAT family are known to be involved in developing the tolerance of plants to various abiotic stress conditions [49–51], mechanistic details about these CCAAT genes and their possible specific involvement in N stress tolerance need further investigation. Some Myb transcription factors were down-regulated under severe N stress, including two Myb factors (At1g25550 and At1g68670), and some were up-regulated (PAP1 and PAP2), indicating the divergent function of individual members within this transcription factor gene family.
Overall, mild N stress triggered only a small set of genes to be expressed significantly differently. These genes are largely involved in various stress responses, including genes involved in anthocyanin biosynthesis. There is a possible disconnect between transcriptional responses and phenotypic changes, as no genes specifically involved in growth arrest or chlorophyll synthesis have significantly altered expression levels. Thus, it appears that there was a compensatory response to this limited level of N stress. Under severe N stress, however, plant responses were much more pronounced. Table 4 summarizes some of the biological processes with significant genes over-represented in the two N stress responses (P < 0.01). Among the up-regulated genes under mild N stress, 16 genes, or 40% of the significant gene list (excluding genes of unknown functions), belong to the "response to abiotic stimulus" process according to GO. Genome-wide, there are 1091 member genes in this process, or 6% of the genome. Apparently, mild N stress invoked a significant response from abiotic stimulus related pathways (P-value = 4.2E-10). Fourteen genes, or 35% of the list, are responsive to different type of stress, including 5 of those responsive to oxidative stress. Six genes, or 15% of the list, are involved in secondary metabolism, including two for flavonoid biosynthesis and two for anthocyanin biosynthesis. Under severe N stress, up-regulated genes are still over-represented by those involved in the response to abiotic stimulus and other types of stress as well as secondary metabolism. However, the relative percentage of these "stress genes" actually decreased when many more genes responded transcriptionally to severe N stress and many of them are involved in other processes, such as the response to abscisic acid stimulus, phenylpropanoid biosynthesis, senescence and starch biosynthesis (Table 4). Genes involved in primary metabolism had a preponderance of down-regulated genes under severe N stress. Fifty-eight genes, or 43% of the list, belong to that category, including 7 involved in nitrogen compound metabolism and 4 in the main pathways of carbohydrate metabolism (Table 4). Twenty-eight biosynthetic genes, or about 20% of the list, are down-regulated, with 17 of these involved in protein biosynthesis. Down-regulated genes are also involved in photosynthesis, the response to auxin and cytokinin and the generation of precursor metabolites and energy (Table 4). The number of genes with significantly altered expression in each of these processes and their degrees of over-representation is noted in Table 4.
Around 50% of the genes with significantly altered expression from the mild N-limiting condition overlap with the genes from the severe N-limiting condition. The overlapping gene list is provided in Additional file 4. Many of these genes are involved in the response to abiotic or biotic stress, such as peroxidases, GST genes and genes involved in anthocyanin synthesis, pointing to possible interactions between general stress and N stress. Two transcription factors PAP2 and ATERF2 (At5g47220) are present in both lists. Some of the overlapping genes such as PAP2 responded much more strongly at the transcript level under severe N stress, while others such as the ATERF2 gene responded at a similar level under both conditions (Table 3). Therefore, plants might have a unique scaled response system for different degrees of N stress by a gradual activation of genes from mild to severe N stress.
Functional classification of significant genes after short- and long-term N availability increase
Three-week-old Arabidopsis plants grown under severe N stress were transferred to the sufficient N condition and harvested after 2 and 24 hrs. After 2 hrs of increased N availability, a total of 77 genes were significantly up-regulated and 13 genes were significantly down-regulated. Twenty-four hrs of increased N availability invoked a greater plant response as the expression levels of 481 genes were significantly up-regulated and 356 were significantly down-regulated. The significant gene lists are provided in Additional file 5 and 6 respectively, together with their functional and pathway assignments. The short-term N availability increase triggered a rapid response of many genes directly involved in plant growth and development. The nitrate assimilatory genes, NR1, NR2 and NiR, were up-regulated significantly after a 2 hr N induction, followed by a drop after 24 hr (Table 2), a pattern consistent with previous reports. The UPM1 gene, which encodes an enzyme that makes a cofactor for NiR, had a very similar trend as NiR (Table 2). Two FNR genes, involved in supplying reductant to NiR, also had a transient increase (Table 2). With regard to the N uptake genes, the low affinity nitrate transporter gene NRT1.1 showed some increase after 2 hr N induction, while NRT1.4 showed significant induction after 24 hr (Table 2), suggesting a different role played by individual transporters in N response. Surprisingly, none of the seven NRT2 members showed significant change (Table 2). High affinity nitrate transporters (eg. NRT2.1) were induced significantly from N-starved seedlings [19, 20], but not in our experiments, possibly due to the fact that plants in this study were grown at a constant low level of nitrate, rather than being N-starved prior to the addition of a high level of nitrate. The expression of GS and GOGAT, the principle players in ammonium assimilation, didn't change significantly after 2 hr. GS1-1 and GS1-4 were repressed after the 24 hr induction (Table 2), which is similar to the previous report . Asparagine synthetase (AS) plays an important role in amino acid synthesis . Among the three AS genes, AS2 (At5g65010) showed a very similar pattern to the nitrate reduction genes (NR1, NR2 and NiR) (Table 2). Genes involved in nitrate storage were significantly up-regulated. The anion channel protein gene, mediating nitrate accumulation in plant vacuoles , was up-regulated around 3-fold after 2 hr of N addition and continued to increase to 6-fold after 24 hr (Table 3).
Many genes involved in carbon metabolism were up-regulated with increased N availability. Sugar transporter genes were up-regulated after 2 hr, together with some amino acid transporters [see Additional file 5]. The expression of the genes coding for the rate-limiting step of the oxidative pentose-phosphate pathway glucose-6-phosphate 1-dehydrogenase was induced. However, the two genes coding for the two isozymes showed a different pattern with one being induced only after 2 hrs and the other only after 24 hrs of increased N availability (Table 3). Genes involved in photosynthesis (photosystem II family protein and photosystem II reaction centre W PsbW family protein) were significantly up-regulated, and genes involved in the accumulation of starch (ADP-glucose pyrophosphorylase) were down-regulated only in the 24 hr induction samples (Table 3). Phytohormone genes coding for cell growth and expansion were induced, such as adenylate isopentenyl transferase 3 (cytokinin synthase 3, At3g63110) (Table 3). In addition, genes involved in sulphate metabolism and in iron acquisition, transport and homeostasis were all up-regulated 2 hrs after increased N availability (see Additional file 5). Seven transcription factor genes were significantly induced or repressed after short-term N induction (2 hr) [see Additional file 5], including up-regulated Myb transcriptional factor genes (At3g46130, At1g25550 and At1g68670) and down-regulated PAP1 and PAP2 genes.
Table 5 summarizes some biological processes over-represented in genes responding to short- and long-term N availability increase (P < 0.01). When N was transiently increased, genes involved in transporter activities were up-regulated. Eleven of these genes, or 20% of the significant gene list (excluding unknown function genes), belong to the "transport" process according to GO. Genome-wide, there are 1398 member genes in this process, or close to 8% of the genome, indicating a clear increase in N stimulated nutrient transport activities (P-value = 0.002). Genes involved in nitrogen assimilation, carbon metabolism, sulfate assimilation and cell homeostasis were simultaneously up-regulated after a short-term N increase (Table 5). The up-regulated genes are over-represented in the process "response to cytokinin stimulus" after both the short- and longer-term N increase (Table 5). After the 24 hr N increase, 203 genes involved in primary metabolism, or 54.7% of the significant gene list, were up-regulated, and 121 genes involved in protein biosynthesis, or 32.61% of the list, were up-regulated (Table 5). Other processes with over-represented up-regulated genes include those involved in nitrogen compound biosynthesis, chlorophyll biosynthesis, organelle organization and biogenesis, DNA packaging, nucleosome assembly, ribonucleotide biosynthesis, translation, and the responses to auxin stimulus (Table 5). On the other hand, down-regulated genes are over-represented in processes such as flavonoid biosynthesis, phenylpropanoid biosynthesis, anthocyanin biosynthesis and starch biosynthesis (Table 5), reflecting the shift from secondary to primary metabolism when more N became available. Down-regulation of anthocyanin biosynthetic genes was also observed after the shorter-term N increase (Table 5). The expression of many genes was significantly altered only after longer-term N increase, suggesting some of the changes may result from a secondary response. The expression levels of many genes fluctuated, some quite dramatically, in response to different degrees of N stress or different lengths of N induction. The most striking change was found with the Myb factor PAP2 (At1g66390), which was up-regulated over 3-fold under mild N stress, up-regulated 26-fold under severe N stress, down-regulated ~7-fold after a 2 hr N increase and continued to fall by over 16-fold after the 24 hr N increase (Table 3), indicating the close correlation between anthocyanin biosynthesis and N availability.
Comparison with published microarray data
Wang et al  published microarray data for roots and shoots of Arabidopsis seedlings grown in a "pure" hydroponic system for 10 d on ammonium and then supplied with 0.25 mM nitrate for 20 min. They found that a total of 76 genes were up-regulated and 2 genes were down-regulated using a cutoff value of 2-fold change in shoots. Scheible et al  revealed genes affected by N-induction after 30 min or 3 hr from N-starved seedlings. The induction patterns and the majority of the nitrate-induced genes of our study overlap with their findings [19, 20]. Most of the nitrate assimilation genes were induced in a similar fashion, such as NR1, NR2, NiR, UPM1 and two of the six FNR genes . However, high affinity nitrate transport gene (NRT2.1) was induced significantly from N-starved seedlings [19, 20], but not in our experiments. This is likely due to the different system used. Rather than being N-starved prior to the addition of a high level of nitrate, plants used in the current study were grown on a constant low level of nitrate, and there was some storage of free nitrate in the plants (Figure 1b). Our data presented here were from samples taken in the middle of the day to minimize diurnal effects, as nitrate levels in the plants could still fluctuate at dawn and at night and thus cause some corresponding transcriptional changes .
Wang et al  identified 6 transcription factor genes in shoots at least 2-fold induced or repressed after 20 min induction. In N-starved seedlings, Scheible et al.  identified 37 and 42 transcription factor genes, with altered expression after 30 min or 3 hr N induction respectively, with at least a 3-fold change. We identified 7 transcription factor genes significantly induced or repressed after short-term N induction (2 hr), all of which overlap with Scheible's gene list .
Additionally, Scheible et al  discovered 55 transcription factor genes at least 3-fold induced or repressed in the N-starved seedlings. We identified 3 transcription factor genes significantly induced or repressed under mild N stress and 39 under severe N stress. However, ~70% of those genes (two of the three from mild N stress, and 27 of the 39 from severe N stress) do not overlap with Scheible's gene list , indicating that plants have a very different regulatory system to cope with N starvation versus chronic N stress. These transcription factor genes not reported previously were highlighted in the Additional Files 2 and 3.
Recently, Lian et al reported expression profiles of 10,422 genes at an early stage of low N stress in rice seedling . They germinated and grew rice hydroponically with normal nutrient conditions until the emergence of the fourth leaf, then transferred the seedlings to a low N solution and harvested shoot and roots separately 20 min, 1 and 2 hr after the treatment. They found the down-regulation of genes involved in photosynthesis and energy metabolism as well as extensive cross-talk between the responses to low N stress and those to biotic and abiotic stresses. However, the genes involved in N uptake and assimilation remained unchanged, which is different from findings in our study, probably due to the fact that their low N treatment was not severe enough to cause immediate N shortage as plant vacuoles have a storage capacity for nitrogen .
Identification of putative nitrogen regulatory elements
The promoters of co-expressed genes are likely to share common regulatory motifs and are potentially regulated by a common set of transcription factors. Therefore, identification of cis-regulatory elements in the promoter regions becomes an important first step in uncovering new facets of transcriptional regulation networks for N responses. Target genes were selected from the three main stages of N assimilation, namely nitrate uptake, nitrate reduction and ammonium assimilation. Included were one of the major nitrate transporters (NRT1.1), nitrate reductase and nitrite reductase (NR1, NR2 and NiR), glutamine synthetase (GS1-1 and GS1-4) and asparagine synthetase (AS2). Data from the five N conditions (15 microarray chips) were used for clustering genes with the expression pattern most similar to those of the selected target genes. The cluster 1 (22 genes) included the target gene NRT1.1, as well as other transporter genes such as CLC-b chloride channel protein (At3g27170), the glutamine transporter (At3g56200), the sulphate transporter protein (At1g23090), the carbohydrate transporter (At1g08930 and At4g17550), and the sodium-dicarboxylate cotransporter (At5g47560), implying the coordinated regulation of different aspects of metabolism. The genes in cluster 1 had higher baseline expression levels under N stress. They were up-regulated after 2 hr N induction followed by an expression decline after 24 hrs (Figure 2). The cluster 2 (26 genes) included four target genes: NR1, NR2, NiR and AS2, as well as genes involved in glucose metabolism (At1g24280) and trehalose metabolism (At5g51460). Since C and N metabolism are very closely linked, it is interesting to note that some key enzymes in C and N metabolism fall into the same cluster. The genes in cluster 2 had a lower baseline under N stress. They were all up-regulated after 2 hr N induction and then decreased after 24 hrs (Figure 2). The cluster 3 (60 genes) included the target gene GS1-1. Genes involved in anthocyanin synthesis (At4g22880), and the rate-limiting enzyme in flavonol and anthocyanin biosynthesis chalcone synthase (At5g13930) were also in this cluster. The cluster 4 (50 genes) included the target gene GS1-4, as well as genes involved in anthocyanin synthesis (At5g42800), and the key enzyme in phenylprepanoid biosynthesis phenylalanine ammonia lyase (At2g37040), suggesting a close linkage between these metabolic pathways. The pattern of cluster 3 and 4 was similar in that they had higher baseline expression under N stress and they were repressed after the 2 hr N induction and then further decreased after 24 hrs (Figure 2). The entire list of genes in the four clusters is provided in Additional file 7, together with their functional and pathway assignments.
The upstream regions 1 kb onward from the ATG start sites of the genes present in the four clusters were analyzed for cis elements. A motif-module discovery algorithm, CisModule , was used to detect motif patterns and the 5 most highly over-represented patterns were found (Table 6). Those candidate motifs were subsequently searched against the PLACE database  to identify those containing previously published plant cis elements in both forward and reverse strand. GATA motifs were found in cluster 1, 2 and 4, Dof motifs were found in all four clusters, Myb motifs found in cluster 1 and 2, and WRKY and CCAAT motifs found in cluster 4 (Table 6).
To compare the results obtained from CisModule, another algorithm BioProspector was used to look for regulatory sequence motifs . The top 5 over-represented patterns found are listed in Table 6. Various motif finding software programs typically yield candidate cis elements not entirely consistent with each other due to the differences in their underlying statistical algorithms. CisModule is based on a Bayesian hierarchical mixture model to infer motifs by their sequence background. BioProspector, on the other hand, uses a Gibbs sampling strategy to identify motifs. Less than five significant patterns were found for cluster 1, 3 and 4 by BioProspector, based on its significance level judged by its probability (p-value < 10-5 in this case). These candidate motifs were subsequently searched against the PLACE database . Again, GATA motifs were found in both cluster 2 and 4, Dof motifs in clusters 1 and 2, and Myb in cluster 2 (Table 6).
GATA motifs have been identified in the regulatory regions of many genes involved in nitrate assimilation such as nitrate reductase, nitrite reductase and glutamine synthetase [29, 56, 57]. Previously we identified regions of the spinach nitrite reductase (NiR) promoter that are involved in N regulation [28, 29, 58]. Footprinting assay suggests that GATA factors play a role in NiR gene regulation . Therefore, the presence of GATA enriched motifs discovered by both algorithms in a NiR gene cluster lends substantial support to the validity of these putative motifs. Interestingly, the GATA motif was found in the NiR cluster as well as in the GS1-4 cluster by both algorithms. However, the expression patterns of these two clusters are very different, suggesting there are positive and negative GATA regulators. Also, the GATA motif was found by both algorithms only in the GS1-4 cluster, but not in the GS1-1 cluster, although the expression patterns of these two clusters are very similar, suggesting a complex regulatory system. Arabidopsis has 30 GATA transcription factor genes [39, 59], but only a few of them have been functionally characterized [60–63], and none has been shown to be directly involved in N regulation.
The Dof transcription factors belong to the same C2C2 zinc finger transcription factor family as the GATA factors [64, 65]. They are associated with expression of multiple genes involved in carbon metabolism in maize , but not in nitrogen metabolism. Dof factors were identified in all four clusters by CisModule and in two clusters by BioProspector. It is tempting to speculate that Dof factors could play an important role in nitrogen regulation. This speculation is supported by the fact that improved nitrogen assimilation and growth under low-nitrogen conditions could be achieved in Dof1 over-expressed plants .
In addition to the Dof factors, binding sequences for other factors such as CCAAT, Myb and WRKY were found in these clusters either by one of the algorithms or by both. Members of the CCAAT, Myb and WRKY factors have been shown to be involved in various abiotic stress regulation/tolerance [46, 49–51, 68]. Whether some of them are more specifically involved in N regulation requires further investigation. The putative regulatory elements identified in these clusters along with several other known motifs involved in N and C metabolism as well as stress response revealed a complex picture of plant N regulation.
We used an Arabidopsis whole genome array for a global evaluation of gene expression under different N conditions. The differentially expressed genes identified provide additional insights into the coordination of the complex N responses of plants and the components of the N response mechanism. Putative N regulatory elements were identified along with several previously known motifs involved in N, C, and stress responses. A better understanding of the complex regulatory network for plant N responses will ultimately lead to strategies to improve N use efficiency in crop plants. Much remains to be done in order to fully construct the regulatory networks underlying this critical aspect of plant biology.
Microarray data availability
All expression data were collected in compliance with the MIAME standards [see Additional file 8] and are available through NASCArrays database with reference number NASCARRAYS-408.
Plant growth conditions
Wild type Arabidopsis thaliana (Columbia ecotype) plants were grown under hydroponic conditions. Plant seeds were sown in rockwool cubes (25 × 25 × 40 mm; Fibrex Insulations Inc, ON, Canada), which were equilibrated with a balanced nutrient solution containing 10 mM KH2PO4, 2 mM MgSO4, 1 mM CaCl2, 0.1 mM Fe-EDTA, 50 μM H3BO4, 12 μM MnSO4, 1 μM ZnCl2, 1 μM CuSO4, 0.2 μM Na2MoO4 and varying levels of nitrate. Around 200 rockwool cubes were then placed in a tray which was connected to an 18-liter nutrient solution reservoir for sub-irrigation. Nutrient solution was pumped once a day to the tray from the reservoir to wet the rockwool completely and flowed back to the reservoir by an adjustable pump (Aquarium systems, OH, USA). Each of the above described hydroponic systems was one experimental unit used for the different nitrate levels. The electrical conductivity (EC) and pH levels of the nutrient solutions in the reservoir were monitored by TDSTestr 10 (OAKTON Instruments, IL, USA) and pH Pro Meter (Spectrum Technologies Inc. IL, USA) and adjusted to their target levels twice a week. There were three nitrate levels (treatments) in the experiment with the 3 mM nitrate as sufficient, 1 mM nitrate as the mild N-limiting condition and 0.3 mM nitrate as the severe N-limiting condition. The targeted EC was 670, 770, and 840 ppm respectively, and pH was around 6.0. Nitrate levels of the nutrient solutions in the reservoir were monitored by Cardy NO3 Nitrate Meter (Spectrum Technologies Inc. IL, USA). For the N induction, rockwool cubes with plants grown under 0.3 mM N were transferred to the tray with a connection to the 3 mM N reservoir. Nutrient solution was then pumped to the tray from each unit at 11 AM and shoots were harvested at 1 PM. The method to detect nitrate levels in the plants was described in Peng et al . Plants were grown in a growth room with 16 hr lighting per day under fluorescent lamps (with a photosynthetically active radiation of 150 μmolm-2s-1) at 23°C and 8 hr dark at 18°C for three weeks before shoots were harvested. Shoot biomass was taken from an average of 6 to 8 plants and nitrate concentration from an average of 3 samples (each sample from a pool of 3 plants). In order to harvest sufficient material for RNA extraction, the number of plants pooled to run the microarray experiment was different. It was approximately 8, 10, and 16 plants respectively for the N-sufficient, mild N-stress and severe N-stress conditions. Plants used for RNA extraction were different from those for fresh weight and nitrate concentration, but they were grown under same conditions.
Five μg of total RNA from each sample was used to synthesize double-stranded cDNAs. Labeled complementary RNA, synthesized from the cDNA was hybridized to a custom designed Arabidopsis whole genome exon GeneChip array (SYNG002) as previously described . The hybridization signals of the arrays were acquired by the GeneChip scanner 3000 and quantified by MAS 5.0 (Affymetrix). Each probe set measurement was summarized as a value of weighted average of all probes in a set, subtracting bottom 5% of average intensity of the entire array using a custom algorithm. The overall intensities of all probe sets of each array were further scaled to a target intensity of 100 to enable direct comparison.
Microarray data analysis
Out of a total of 26412 × 15 = 396180 data points in the five treatment groups, the custom algorithm flagged 20149 or 5% as 'Absent' (A), 153373 or 39% as 'Marginal' (M), and 222658 or 56% as 'Present' (P). Stepwise gene filtering was conducted in GeneSpring (Agilent, CA, USA). First, for each of the five treatment groups, each gene must have either a 'P' or 'M' flag in all 3 replicate samples. This was followed by a filtering second step requiring that at least one of the three samples had a 'P' flag. This essentially guaranteed that every gene remaining in a group would be 'PMM', 'PPM', or 'PPP' among the three replicates. For pairwise group comparisons determining differentially expressed genes, genes common to both groups were identified and data exported, log2 transformed, and analyzed in Statistical Analysis of Microarray (SAM) as two class unpaired t-test and at a permutation number of 500 .
For clustering analysis, those genes present in at least 2 out of the 5 treatment groups, totaling 17138, were exported into Gene Cluster 3.0 with the following clustering parameters selected: log2 transformation, Hierachical clustering, Correlation (centered), and Average linkage.
Arabidopsis motif analyses
Gene clusters containing several target genes were identified with correlation coefficients within a cluster greater than 0.9. The corresponding 1 kb upstream sequences for each gene were downloaded from The Arabidopsis Information Resource  as of June 14, 2006 and used as input files for motif search in BioProspector  and CisModule . For BioProspector, the entire 1 kb upstream sequence of Arabidopsis genome was used as the background distribution. The other parameters were set as follows: motif width, 6 to 16 bps; number of Monte Carlo simulations, 100; number of times to search for motifs, 100; number of top motifs to report, 10. For CisModule, all upstream sequences were masked by Repeatmasker  prior to motif analysis. The program settings were: number of motifs, 3 or 5; degenerate to independent motif sampler; total number of iterations, 1000; search both strands; and run dataset 10 times. Motif outputs from both BioProspector and CisModule were subsequently searched against PLACE database  as of August 31, 2006 to identify those containing previously published plant cis-elements.
Lea PJ, Morot-Gaudry JF: Plant nitrogen. 2001, Berlin ; London , Springer, xi, 407 p.-
Scheible WR, Lauerer M, Schulze ED, Caboche M, Stitt M: Accumulation of nitrate in the shoot acts as a signal to regulate shoot-root allocation in tobacco. Plant Journal. 1997, 11 (4): 671-691. 10.1046/j.1365-313X.1997.11040671.x.
Zhang HM, Forde BG: An Arabidopsis MADS box gene that controls nutrient-induced changes in root architecture. Science. 1998, 279 (5349): 407-409. 10.1126/science.279.5349.407.
Stitt M, Krapp A: The interaction between elevated carbon dioxide and nitrogen nutrition: the physiological and molecular background. Plant Cell and Environment. 1999, 22 (6): 583-621. 10.1046/j.1365-3040.1999.00386.x.
Marschner H: Mineral nutrition of higher plants. 1995, London , Academic, xv,889p.-2nd ed.
Good AG, Shrawat AK, Muench DG: Can less yield more? Is reducing nutrient input into the environment compatible with maintaining crop production?. Trends Plant Sci. 2004, 9 (12): 597-605. 10.1016/j.tplants.2004.10.008.
Nosengo N: Fertilized to death. Nature. 2003, 425 (6961): 894-895. 10.1038/425894a.
Giles J: Nitrogen study fertilizes fears of pollution. Nature. 2005, 433 (7028): 791-791. 10.1038/433791a.
Crawford NM, Forde BG: Molecular and Developmental Biology of Inorganic Nitrogen Nutrition. The Arabidopsis Book. Edited by: Meyerowitz EM. 2002, Rockville, MD , American Society of Plant Biologists
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.
Crawford NM: Nitrate - Nutrient and Signal for Plant-Growth. Plant Cell. 1995, 7 (7): 859-868. 10.1105/tpc.7.7.859.
Stitt M: Nitrate regulation of metabolism and growth. Current Opinion in Plant Biology. 1999, 2 (3): 178-186. 10.1016/S1369-5266(99)80033-8.
Lam HM, Coschigano KT, Oliveira IC, Melo-Oliveira R, Coruzzi GM: The Molecular-Genetics of Nitrogen Assimilation into Amino Acids in Higher Plants. Annu Rev Plant Physiol Plant Mol Biol. 1996, 47: 569-593. 10.1146/annurev.arplant.47.1.569.
Campbell WH: NITRATE REDUCTASE STRUCTURE, FUNCTION AND REGULATION: Bridging the Gap between Biochemistry and Physiology. Annu Rev Plant Physiol Plant Mol Biol. 1999, 50: 277-303. 10.1146/annurev.arplant.50.1.277.
Cheng CL, Dewdney J, Kleinhofs A, Goodman HM: Cloning and nitrate induction of nitrate reductase mRNA. Proc Natl Acad Sci U S A. 1986, 83 (18): 6825-6828. 10.1073/pnas.83.18.6825.
Crawford NM, Campbell WH, Davis RW: Nitrate reductase from squash: cDNA cloning and nitrate regulation. Proc Natl Acad Sci U S A. 1986, 83 (21): 8073-8076. 10.1073/pnas.83.21.8073.
Scheible WR, GonzalezFontes A, Lauerer M, MullerRober B, Caboche M, Stitt M: Nitrate acts as a signal to induce organic acid metabolism and repress starch metabolism in tobacco. Plant Cell. 1997, 9 (5): 783-798. 10.1105/tpc.9.5.783.
Wang R, Guegler K, LaBrie ST, Crawford NM: Genomic analysis of a nutrient response in Arabidopsis reveals diverse expression patterns and novel metabolic and potential regulatory genes induced by nitrate. Plant Cell. 2000, 12 (8): 1491-1509. 10.1105/tpc.12.8.1491.
Wang RC, Okamoto M, Xing XJ, 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 Physiology. 2003, 132 (2): 556-567. 10.1104/pp.103.021253.
Scheible WR, 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 Physiology. 2004, 136 (1): 2483-2499. 10.1104/pp.104.047019.
Coruzzi G, Bush DR: Nitrogen and carbon nutrient and metabolite signaling in plants. Plant Physiol. 2001, 125 (1): 61-64. 10.1104/pp.125.1.61.
Coruzzi GM, Zhou L: Carbon and nitrogen sensing and signaling in plants: emerging 'matrix effects'. Curr Opin Plant Biol. 2001, 4 (3): 247-253. 10.1016/S1369-5266(00)00168-0.
Price J, Laxmi A, St Martin SK, Jang JC: Global transcription profiling reveals multiple sugar signal transduction mechanisms in Arabidopsis. Plant Cell. 2004, 16 (8): 2128-2150. 10.1105/tpc.104.022616.
Lian X, Wang S, Zhang J, Feng Q, Zhang L, Fan D, Li X, Yuan D, Han B, Zhang Q: Expression profiles of 10,422 genes at early stage of low nitrogen stress in rice assayed using a cDNA microarray. Plant Mol Biol. 2006, 60 (5): 617-631. 10.1007/s11103-005-5441-7.
Lin Y, Hwang CF, Brown JB, Cheng CL: 5' proximal regions of Arabidopsis nitrate reductase genes direct nitrate-induced transcription in transgenic tobacco. Plant Physiol. 1994, 106 (2): 477-484. 10.1104/pp.106.2.477.
Neininger A, Back E, Bichler J, Schneiderbauer A, Mohr H: Deletion Analysis of a Nitrite-Reductase Promoter from Spinach in Transgenic Tobacco. Planta. 1994, 194 (2): 186-192. 10.1007/BF01101677.
Dorbe MF, Truong HN, Crete P, Daniel-Vedele F: Deletion analysis of the tobacco Nii1 promoter in Arabidopsis thaliana. Plant Science. 1998, 139 (1): 71-82. 10.1016/S0168-9452(98)00181-2.
Back E, Dunne W, Schneiderbauer A, de Framond A, Rastogi R, Rothstein SJ: Isolation of the spinach nitrite reductase gene promoter which confers nitrate inducibility on GUS gene expression in transgenic tobacco. Plant Mol Biol. 1991, 17 (1): 9-18. 10.1007/BF00036801.
Rastogi R, Bate NJ, Sivasankar S, Rothstein SJ: Footprinting of the spinach nitrite reductase gene promoter reveals the preservation of nitrate regulatory elements between fungi and higher plants. Plant Mol Biol. 1997, 34 (3): 465-476. 10.1023/A:1005842812321.
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-
Naito S, Hirai MY, Chino M, Komeda Y: Expression of a Soybean (Glycine max [L.] Merr.) Seed Storage Protein Gene in Transgenic Arabidopsis thaliana and Its Response to Nutritional Stress and to Abscisic Acid Mutations. Plant Physiol. 1994, 104 (2): 497-503.
Robinson D: What limits nitrate uptake from soil?. Plant Cell Environ. 1991, 14: 77-85. 10.1111/j.1365-3040.1991.tb01373.x.
Scheible WR, Krapp A, Stitt M: Reciprocal diurnal changes of phosphoenolpyruvate carboxylase expression and cytosolic pyruvate kinase, citrate synthase and NADP-isocitrate dehydrogenase expression regulate organic acid metabolism during nitrate assimilation in tobacco leaves. Plant Cell and Environment. 2000, 23 (11): 1155-1167. 10.1046/j.1365-3040.2000.00634.x.
Zhu T: Global analysis of gene expression using GeneChip microarrays. Current Opinion in Plant Biology. 2003, 6 (5): 418-425. 10.1016/S1369-5266(03)00083-9.
Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci U S A. 2001, 98 (9): 5116-5121. 10.1073/pnas.091062498.
Misson J, Raghothama KG, Jain A, Jouhet J, Block MA, Bligny R, Ortet P, Creff A, Somerville S, Rolland N, Doumas P, Nacry P, Herrerra-Estrella L, Nussaume L, Thibaud MC: A genome-wide transcriptional analysis using Arabidopsis thaliana Affymetrix gene chips determined plant responses to phosphate deprivation. Proc Natl Acad Sci U S A. 2005, 102 (33): 11934-11939. 10.1073/pnas.0505266102.
Gene Ontology database. [http://geneontology.org]
Kyoto Encyclopedia of Genes and Genomes . [http://www.genome.ad.jp/kegg]
Arabidopsis Information Resource website . [http://www.arabidopsis.org/tools/bulk/go/index.jsp]
Provart NJ, Gil P, Chen W, Han B, Chang HS, Wang X, Zhu T: Gene expression phenotypes of Arabidopsis associated with sensitivity to low temperatures. Plant Physiol. 2003, 132 (2): 893-906. 10.1104/pp.103.021261.
Coruzzi G: Primary N-assimilation into amino acids in Arabidopsis. The Arabidopsis book. Edited by: Meyerowitz EM. 2003, Rockville, MD , American Society of Plant Biologists
De Angeli A, Monachello D, Ephritikhine G, Frachisse JM, Thomine S, Gambale F, Barbier-Brygoo H: The nitrate/proton antiporter AtCLCa mediates nitrate accumulation in plant vacuoles. Nature. 2006, 442 (7105): 939-942. 10.1038/nature05013.
Diaz U, Saliba-Colombani V, Loudet O, Belluomo P, Moreau L, Daniel-Vedele F, Morot-Gaudry JF, Maselaux-Daubresse U: Leaf yellowing and anthocyanin accumulation are two genetically independent strategies in response to nitrogen limitation in Arabidopsis thaliana. Plant and Cell Physiology. 2006, 47 (1): 74-83. 10.1093/pcp/pci225.
Dixon RA, Paiva NL: Stress-Induced Phenylpropanoid Metabolism. Plant Cell. 1995, 7 (7): 1085-1097. 10.1105/tpc.7.7.1085.
Borevitz JO, Xia YJ, Blount J, Dixon RA, Lamb C: Activation tagging identifies a conserved MYB regulator of phenylpropanoid biosynthesis. Plant Cell. 2000, 12 (12): 2383-2393. 10.1105/tpc.12.12.2383.
Eulgem T, Rushton PJ, Robatzek S, Somssich IE: The WRKY superfamily of plant transcription factors. Trends in Plant Science. 2000, 5 (5): 199-206. 10.1016/S1360-1385(00)01600-9.
Nakano T, Suzuki K, Ohtsuki N, Tsujimoto Y, Fujimura T, Shinshi H: Identification of genes of the plant-specific transcription-factor families cooperatively regulated by ethylene and jasmonate in Arabidopsis thaliana. Journal of Plant Research. 2006, 119 (4): 407-413. 10.1007/s10265-006-0287-x.
Fujimoto SY, Ohta M, Usui A, Shinshi H, Ohme-Takagi M: Arabidopsis ethylene-responsive element binding factors act as transcriptional activators or repressors of GCC box-mediated gene expression. Plant Cell. 2000, 12 (3): 393-404. 10.1105/tpc.12.3.393.
Chinnusamy V, Schumaker K, Zhu JK: Molecular genetic perspectives on cross-talk and specificity in abiotic stress signalling in plants. Journal of Experimental Botany. 2004, 55 (395): 225-236. 10.1093/jxb/erh005.
Alonso-Blanco C, Gomez-Mena C, Llorente F, Koornneef M, Salinas J, Martinez-Zapater JM: Genetic and molecular analyses of natural variation indicate CBF2 as a candidate gene for underlying a freezing tolerance quantitative trait locus in Arabidopsis. Plant Physiology. 2005, 139 (3): 1304-1312. 10.1104/pp.105.068510.
Lee BH, Henderson DA, Zhu JK: The Arabidopsis cold-responsive transcriptome and its regulation by ICE1. Plant Cell. 2005, 17 (11): 3155-3175. 10.1105/tpc.105.035568.
Lejay L, Tillard P, Lepetit M, Olive F, Filleur S, Daniel-Vedele F, Gojon A: Molecular and functional regulation of two NO3- uptake systems by N- and C-status of Arabidopsis plants. Plant J. 1999, 18 (5): 509-519. 10.1046/j.1365-313X.1999.00480.x.
Zhou Q, Wong WH: CisModule: de novo discovery of cis-regulatory modules by hierarchical mixture modeling. Proc Natl Acad Sci U S A. 2004, 101 (33): 12114-12119. 10.1073/pnas.0402858101.
Higo K, Ugawa Y, Iwamoto M, Korenaga T: Plant cis-acting regulatory DNA elements (PLACE) database: 1999. Nucleic Acids Res. 1999, 27 (1): 297-300. 10.1093/nar/27.1.297.
Liu X, Brutlag DL, Liu JS: BioProspector: discovering conserved DNA motifs in upstream regulatory regions of co-expressed genes. Pac Symp Biocomput. 2001, 127-138.
Jarai G, Truong HN, Daniel-Vedele F, Marzluf GA: NIT2, the nitrogen regulatory protein of Neurospora crassa, binds upstream of nia, the tomato nitrate reductase gene, in vitro. Curr Genet. 1992, 21 (1): 37-41. 10.1007/BF00318652.
Oliveira IC, Coruzzi GM: Carbon and amino acids reciprocally modulate the expression of glutamine synthetase in Arabidopsis. Plant Physiol. 1999, 121 (1): 301-310. 10.1104/pp.121.1.301.
Rastogi R, Dulson J, Rothstein SJ: Cloning of tomato (Lycopersicon esculentum Mill.) arginine decarboxylase gene and its expression during fruit ripening. Plant Physiol. 1993, 103 (3): 829-834. 10.1104/pp.103.3.829.
Reyes JC, Muro-Pastor MI, Florencio FJ: The GATA family of transcription factors in Arabidopsis and rice. Plant Physiol. 2004, 134 (4): 1718-1732. 10.1104/pp.103.037788.
Nishii A, Takemura M, Fujita H, Shikata M, Yokota A, Kohchi T: Characterization of a novel gene encoding a putative single zinc-finger protein, ZIM, expressed during the reproductive phase in Arabidopsis thaliana. Biosci Biotechnol Biochem. 2000, 64 (7): 1402-1409. 10.1271/bbb.64.1402.
Zhao Y, Medrano L, Ohashi K, Fletcher JC, Yu H, Sakai H, Meyerowitz EM: HANABA TARANU is a GATA transcription factor that regulates shoot apical meristem and flower development in Arabidopsis. Plant Cell. 2004, 16 (10): 2586-2600. 10.1105/tpc.104.024869.
Bi YM, Zhang Y, Signorelli T, Zhao R, Zhu T, Rothstein S: Genetic analysis of Arabidopsis GATA transcription factor gene family reveals a nitrate-inducible member important for chlorophyll synthesis and glucose sensitivity. Plant Journal. 2005, 44 (4): 680-692. 10.1111/j.1365-313X.2005.02568.x.
Liu PP, Koizuka N, Martin RC, Nonogaki H: The BME3 (Blue Micropylar End 3) GATA zinc finger transcription factor is a positive regulator of Arabidopsis seed germination. Plant J. 2005, 44 (6): 960-971. 10.1111/j.1365-313X.2005.02588.x.
Riechmann JL, Heard J, Martin G, Reuber L, Jiang C, Keddie J, Adam L, Pineda O, Ratcliffe OJ, Samaha RR, Creelman R, Pilgrim M, Broun P, Zhang JZ, Ghandehari D, Sherman BK, Yu G: Arabidopsis transcription factors: genome-wide comparative analysis among eukaryotes. Science. 2000, 290 (5499): 2105-2110. 10.1126/science.290.5499.2105.
Yanagisawa S: The Dof family of plant transcription factors. Trends Plant Sci. 2002, 7 (12): 555-560. 10.1016/S1360-1385(02)02362-2.
Yanagisawa S: Dof1 and Dof2 transcription factors are associated with expression of multiple genes involved in carbon metabolism in maize. Plant J. 2000, 21 (3): 281-288. 10.1046/j.1365-313x.2000.00685.x.
Yanagisawa S, Akiyama A, Kisaka H, Uchimiya H, Miwa T: Metabolic engineering with Dof1 transcription factor in plants: Improved nitrogen assimilation and growth under low-nitrogen conditions. Proc Natl Acad Sci U S A. 2004, 101 (20): 7833-7838. 10.1073/pnas.0402267101.
Chen YH, Yang XY, He K, Liu MH, Li JG, Gao ZF, Lin ZQ, Zhang YF, Wang XX, Qiu XM, Shen YP, Zhang L, Deng XH, Luo JC, Deng XW, Chen ZL, Gu HY, Qu LJ: The MYB transcription factor superfamily of arabidopsis: Expression analysis and phylogenetic comparison with the rice MYB family. Plant Molecular Biology. 2006, 60 (1): 107-124. 10.1007/s11103-005-2910-y.
Peng M, Hannam C, Gu H, Bi YM, Rothstein SJ: A mutation in NLA, which encodes a RING-type ubiquitin ligase, disrupts the adaptability of Arabidopsis to nitrogen limitation. Plant J. 2007, 50 (2): 320-337. 10.1111/j.1365-313X.2007.03050.x.
Zhu T, Budworth P, Han B, Brown D, Chang HS, Zou GZ, Wang X: Toward elucidating the global gene expression patterns of developing Arabidopsis: Parallel analysis of 8 300 genes by a high-density oligonucleotide probe array. Plant Physiology and Biochemistry. 2001, 39 (3-4): 221-242. 10.1016/S0981-9428(00)01234-1.
Smit AFA, Hubley R, Green P: Repeatmasker Open-3.0. 1996, http://www.repeatmasker.org
We acknowledge M. Ono and C. Chilcott (Syngenta) for conducting microarray experiments; Dr. J. Colasanti (University of Guelph) for critical reading of the manuscript; Dr. Y. Zheng (University of Guelph) for advice on hydroponic system; J. Zhang (University of Guelph) for help with GeneSpring; T. Signorelli and R. Zhao (University of Guelph) for help growing plants. This work was supported by the Natural Sciences and Engineering Research Council of Canada and Syngenta Biotechnology Inc. to S.J.R. It has been reviewed by the US EPA and approved for publication.
YMB designed and conducted the experiment, analyzed data and drafted the manuscript. RLW performed data analysis, searched putative regulatory elements and provided critical revisions to the manuscript. TZ coordinated the microarray experiments and provided insights on the motif search. SJR designed and coordinated the experiment and finalized the manuscript. All authors read and approved the final manuscript.