Evaluation of protein pattern changes in roots and leaves of Zea mays plants in response to nitrate availability by two-dimensional gel electrophoresis analysis

Background Nitrogen nutrition is one of the major factors that limit growth and production of crop plants. It affects many processes, such as development, architecture, flowering, senescence and photosynthesis. Although the improvement in technologies for protein study and the widening of gene sequences have made possible the study of the plant proteomes, only limited information on proteome changes occurring in response to nitrogen amount are available up to now. In this work, two-dimensional gel electrophoresis (2-DE) has been used to investigate the protein changes induced by NO3- concentration in both roots and leaves of maize (Zea mays L.) plants. Moreover, in order to better evaluate the proteomic results, some biochemical and physiological parameters were measured. Results Through 2-DE analysis, 20 and 18 spots that significantly changed their amount at least two folds in response to nitrate addition to the growth medium of starved maize plants were found in roots and leaves, respectively. Most of these spots were identified by Liquid Chromatography Electrospray Ionization Tandem Mass Spectrometry (LC-ESI-MS/MS). In roots, many of these changes were referred to enzymes involved in nitrate assimilation and in metabolic pathways implicated in the balance of the energy and redox status of the cell, among which the pentose phosphate pathway. In leaves, most of the characterized proteins were related to regulation of photosynthesis. Moreover, the up-accumulation of lipoxygenase 10 indicated that the leaf response to a high availability of nitrate may also involve a modification in lipid metabolism. Finally, this proteomic approach suggested that the nutritional status of the plant may affect two different post-translational modifications of phosphoenolpyruvate carboxylase (PEPCase) consisting in monoubiquitination and phosphorylation in roots and leaves, respectively. Conclusion This work provides a first characterization of the proteome changes that occur in response to nitrate availability in leaves and roots of maize plants. According to previous studies, the work confirms the relationship between nitrogen and carbon metabolisms and it rises some intriguing questions, concerning the possible role of NO and lipoxygenase 10 in roots and leaves, respectively. Although further studies will be necessary, this proteomic analysis underlines the central role of post-translational events in modulating pivotal enzymes, such as PEPCase.

The low bio-availability of nitrogen in the pedosphere with respect to the request of the crops has spawned a dramatic increase in fertilization that has detrimental consequences on environment such as water eutrophication and increase in NH 3 and N 2 O in the atmosphere [6,8].
Moreover, this side-effect is severe in the case of cereals, which account for 70% of food production worldwide. Indeed, in these crops the grain yield is strictly correlated with N supply but the use efficiency is not higher than 50% [9].
Because of the economical relevance, the feasibility to combine extensive physiological, agronomic and genetic studies as well as the high metabolic efficiency of C 4 plants, maize (Zea mays L.) was proposed as the model species to study N nutrition in cereals [10].
Among nitrogen inorganic molecules, nitrate is the predominant form in agricultural soils, where it can reach concentrations three or more orders of magnitude higher than in natural soils [11,12].
The first step of nitrate assimilation, that occurs in both roots and shoots, involves its reduction to ammonia by nitrate reductase (NR) and nitrite reductase (NiR) enzymes, followed by transfer of ammonia to α-chetoglutaric acid by the action of glutamine synthetase (GS) and glutamate synthase (GOGAT) [23][24][25]. The pathway is induced in the presence of nitrate and shows many connections with other cellular traits, among which carbohydrate and amino acid metabolism, redox status and pH homeostasis [6,19,26,27]. Hence, nitrate and carbon metabolisms appear strictly linked and co-regulated, both locally and at long distance for the reciprocal root/leaf control, in response to the nutritional status of the plant and environmental stimuli [3,6,[26][27][28].
In the last years, some transcriptomic analyses have been conducted to shed light on the molecular basis of these regulatory mechanisms. Wang and co-workers studied the transcriptomic changes occurring after exposure to low and high nitrate concentrations in whole plants of Arabidopsis thaliana, by means of microarray and RNA gel blot analysis [29]. Besides the genes already known to be regulated by the presence of nitrate, the authors found new candidate genes encoding for regulatory proteins such as a MYB transcription factor, a calcium antiporter, putative protein kinases and several metabolic enzymes. Another study conducted by Scheible and co-workers [7] reports a comparative transcriptomic analysis of Arabidopsis thaliana seedlings grown in sterile liquid culture under nitrogen-limiting and nitrogen-replete conditions by using Affymetrix ATH1 arrays and (RT)-PCR. The authors observed that the response to nitrogen availability involved a deep reprogramming of primary and secondary metabolisms. These data well describe the complexity of nitrogen pathway as well as the direct and/or indirect consequences that nitrogen availability exerts on the whole metabolism of the plant.
Starting from these results it should be now desirable to deepen the knowledge about the changes at translational and post-translational levels in response to nitrogen availability. In the last decade, the improvement in technologies for protein study and the widening of gene sequences made possible the study of the plant proteomes [30][31][32][33][34].
In this context, the availability of a large EST assembly and the efforts in sequencing maize genome [35] contributed to improve the use of maize, as highlighted by a large number of studies conducted on this species, among which the proteomic characterizations of leaf [36], of chloroplasts in bundle sheath and mesophyll cells [37] and of pericycle cells of primary roots [38].
At the present time, to the best of our knowledge no studies on nitrogen nutrition in maize were conducted by this approach. The only two proteomic works regarding this issue in cereals are based on the use of 2-DE to compare the leaves [39] and the roots [40] of two wheat varieties exposed to different levels of nitrogen. These works pointed out some significant differences, correlated to N availability during the plant growth, in the protein profiles of both organs.
In order to obtain further information, in this work we investigated protein accumulation changes induced by nitrate in both roots and leaves of Zea mays plants. The attention was focused on the changes in the pattern of protein soluble fractions caused by the addition of 10 mM nitrate to the hydroponic solution, after a period in which the plants were grown in the absence of nitrogen. Firstly, the changes of some biochemical parameters were measured to describe the physiological response occurring after nitrate addition and were used to define the sampling time for proteomic analysis. These experiments led to compare the proteomes of plants previously grown for 17 days in absence of nitrogen and incubated for further 30 h without the nutrient or in the presence of 10 mM nitrate. Through 2-DE and LC-ESI-MS/MS analyses a first characterization of the proteome changes occurring in maize plants in response to an increase in nitrate availability was obtained. The results show how many of these changes were related to enzymes of the nitrate assimilation or metabolic pathways strictly linked to it (e.g. pentose phosphate pathway and photosynthesis), but also reveal new proteins that may play a role in the nitrate responses.

Experimental design and biochemical parameters
The aim of this work was to apply a proteomic approach to study the changes in protein patterns of root and leaf organs of maize plants in the first phase of exposure to high availability of nitrate, comparable to agricultural conditions, after a growth period under nitrogen starvation. This is a typical condition in which the addition of nitrate induces an increase in uptake and assimilation of this nutrient [5,28].
The need for a simultaneous analysis of the root and the leaf organs of starved plants, with completely developed but not stressed leaf apparatus, led to the definition of the experimental design showed in Figure 1. Briefly, seedlings were transferred into a hydroponic system after 3 days of germination and grown for further 14 days in a solution deprived of nitrogen. After that, at the beginning of the light period (T 0 ), some plants were maintained in the same nutritional condition (control, C) whereas others were transferred in a nutrient solution containing 10 mM NO 3 -(N). In order to define the sampling time for pro-teomic analysis, the changes of biochemical parameters in response to NO 3were firstly evaluated. Roots and leaves were collected at T 0 time and after 6, 30 and 54 h of nitrate exposure.
At these sampling times, the plants achieved the developmental stage corresponding to the complete expansion of the third leaf (pictures of harvested plants are showed in Additional file 1). The qualitative comparison between the C and N plants revealed some morphological differences. In particular, while the plants appeared very similar at the T 0 sampling time, after 30 h the expansion of the fourth leaf was slightly more evident in N plants with respect to the C ones. This trend was more pronounced at 54 h and, only in C plants, was accompanied by the comparison of faint yellow areas in the leaf blades. In the tested conditions, no significant differences were observed in root system.
In order to characterize the physiological status of the plants, the changes in nitrate content and NR activity (Figure 2) as well as the levels of proteins, amino acids, reducing sugars, sucrose and chlorophyll were evaluated ( Figure  3).
In roots and leaves of starved plants, both nitrate and NR activity were undetectable. After the addition of the nutrient to hydroponic solution the levels of nitrate progressively increased in plant tissues, reaching a level of 32.6 and 10.3 μmol of NO 3 g -1 FW after 54 h in roots and leaves respectively (Figure 2A). A parallel dramatic increase of NR activity was measured until the 30 th h of NO 3exposure, while at the longest time considered (54 h) a decreased activity was observed ( Figure 2B). This trend was more evident in the roots in which a more rapid and large availability of nitrate took place. The total protein levels did not change significantly in all the conditions tested ( Figure 3A and 3B), while a sharp increase in free amino acids was detected in both organs after nitrate addition ( Figure 3C and 3D). Moreover, the levels of amino acids were higher in the leaves than in the roots. Although many factors are involved in the overall amino acid levels, these results may suggest a contribution of translocation of nitrogen compounds between the two organs. Nitrate exposure also induced a decrease in reducing sugars in both organs ( Figure 3E and Figure 3F), while only in the roots of the plants exposed for 54 h to 10 mM NO 3 a drop of sucrose took place ( Figure 3G).
Taken together, these results well describe the induction trend of NO 3assimilation pathway, as suggested by the increase of NR activity and amino acids accompanied by the consequent decrease of reducing sugars, the main source of carbon skeletons [41]. In roots, where photosynthesis cannot satisfy this request and/or the demand of Experimental design Figure 1 Experimental design. Zea mays seeds were germinated in the dark. After 3 days, the seedlings were transferred in a hydroponic system and grown for 14 days in the absence of nitrogen (T 0 ), afterwards the plants were incubated for further 54 h in the same condition (Control, C) or in the presence of 10 mM KNO 3 (N). For details see the methods section.
carbon skeleton is high, sucrose pool was also affected.
The changes in carbohydrate availability and the increase of amino acid levels also explain the decrease in NR activity observed in roots at the 54 th h. In fact, these data are in agreement with the inhibitory effect on NR evocated by an increase of some amino acids, mainly asparagine and glutamine [5,42]. Moreover, it is know that NR activity increases after sucrose addition whilst the low sugar content, condition that we observed in the roots of N plants, affects the nitrate reduction system [5,42,43]. The results suggested that this feedback mechanism was activated in roots of the plants exposed for 54 h to 10 mM NO 3 -. Finally, only at the 54 th h, a significant decrease in chlorophyll content ( Figure 3I) was measured in the leaves of starved plants, thus suggesting that the first symptoms of stress were appearing.

2-DE analysis and protein identification
The biochemical and physiological data showed that the plants incubated for the last 30 h in the presence of 10 mM NO 3 were in a condition in which nitrogen metabolism is completely activated in both root and leaf organs and that, at the same time, no stress symptoms were detectable in the control plants. Starting from these results, the proteomic study was conducted by analyzing the soluble protein fractions extracted from roots and leaves of plants incubated for the last 30 h in the absence or in the presence of 10 mM NO 3 -.
The ratio between dry and fresh weight as well as the total protein content appeared similar both in the roots and in the leaves of C and N samples ( Table 1). The adopted protocol permitted to obtain an extraction yield of soluble proteins of about 14% and 20% for roots and leaves, respectively. Moreover, no significant differences were observed between C and N plants.  Figure 4. The electrophoretic analyses detected about 1100 and 1300 spots in roots and leaves gels, respectively. To ascertain the quantitative changes in the proteomic maps, the relative spot volumes (%Vol) were evaluated by software-assisted analysis. The Student's t-test (p < 0.05), coupled with a threshold of two-fold change in the amount, revealed that 20 spots in roots and 18 spots in leaves were affected by nitrogen availability.
The analysis of these spots by LC-ESI-MS/MS allowed to identify 15 and 14 proteins in root and leaf patterns, respectively. These proteins and the changes in their accumulation are shown in Tables 2 and 3, while further information of mass spectrometry (MS) analysis are reported in the Additional files 2 and 3.

Functional role and quantitative change of the proteins identified in roots
Many of the spots identified in roots were enzymes involved in nitrogen and carbon metabolisms ( Table 2). According to the induction of the NO 3assimilation pathway, in the roots of the plants incubated for the last 30 h in the presence of the nutrient, we observed an increase in the accumulation of nitrite reductase (spot 268, NiR) and of glutamine synthetase plastidial isoform (spot 483, GS2). Total proteins, amino acids, reducing sugars, sucrose and chlorophyll content Figure 3 Total proteins, amino acids, reducing sugars, sucrose and chlorophyll content. Time course of the changes in the content of total proteins, amino acids, reducing sugars and sucrose in roots (A, C, E and G) and leaves (B, D, F, and H) and chlorophyll content in leaves (I) of Zea mays plants, previously grown in the absence of nitrogen for 17 days (T 0 ) and incubated for further 6, 30 and 54 h in the absence (C) or presence of 10 mM NO 3 -(N). Values are the mean ± SE of three independent biological samples analyzed in triplicate (n = 9). Samples indicated with the same letters do not differ significantly according to Tukey's test (p < 0.01).

Nitrate content and nitrate reductase activity
Moreover, in response to the demand of carbon skeletons and NADPH, which is used in non-green tissues for ferredoxin reduction [44], an increase in the levels of phosphoglycerate mutase (spot 216, PGAM-1), glucose-6phosphate dehydrogenase (spot 1162, G6PD) and 6phospho-gluconate dehydrogenase (spot 392, 6PGD) took place. These results well agree with previous array data that describe the responses to nitrate exposure in Arabidopsis and tomato [7,29,45].
An increase in accumulation of the cytosolic isoform of glutamine synthetase (spot 538, GS1-1) was also detected in roots of N plants. On the basis of identified peptides by MS analysis it was possible to discriminate among the 5 GS1 isoforms known in Zea mays (SwissProt reviewed database) and to restrict the possible identification to 2 of them (GS1-1 [Swiss-Prot:P38559] and GS1-5 [Swiss-Prot:P38563] [46]). The fact that Li and co-workers [46], through a Northern blot hybridization analysis, found that the transcript of GS1-1 gene was the only one expressed in roots, conducted to the specific identification of GS1-1 protein. Moreover, Sakakibara and co-workers [47] showed that GS1-1 transcript was the only induced by NO 3 -. The proteomic approach used in the present work allows to confirm these results at the translational level, demonstrating that in maize roots a cytosolic ammonia assimilation pathway can be activated also in response to nitrate. In the table, the fresh/dry weight (FW/DW), the content of total protein (mg g -1 FW) and the % yield of the extraction of soluble proteins (% of extracted soluble proteins respect to the total content) for the roots and the leaves of the plants compared by proteomic analysis are reported. The fresh weight of the roots was 0.56 ± 0.03 and 0.60 ± 0.04 g in C and N plants, respectively. The fresh weight of the leaves was 0.79 ± 0.03 and 0.86 ± 0.04 g in C and N plants, respectively.  Other spots that were found to increase their relative volumes in response to nitrate were a non-symbiotic hemoglobin and a monodehydroascorbate reductase (spot 960, Hb2 and spot 390, MDHAR). In a previous work on Arabidopsis, it was found that NO 3induced AtHB1 and AtHB2, two genes that encode for non-symbiotic hemoglobins [7,29]. Scheible and co-workers [7] suggested that these proteins could change their abundance in relation to the redox status, whereas Wang and co-workers [29] speculated on the possibility that the induction of hemoglobin could aim at reducing oxygen concentration during NR synthesis, since molybdenium can be sensitive to oxygen. Besides, hemoglobin and MDHAR are known to be involved in the scavenging of NO that can be produced by cytosolic and/or plasmamembrane nitrate reductase when nitrite is used as substrate [48,49]. NO is a signaling molecule which is involved in many biochemical and physiological processes [50]. It has been reported that in plant roots, NO plays a role in growth, development and in some responses to environmental conditions, such as hypoxia [51]. Recently, a pos-sible involvement of NO in the mediation of nitratedependent root growth in maize has been suggested [52]. According to this work, that describes a reduction of endogenous NO at high external NO 3concentration, the observed concomitant up-accumulation of Hb2 and MDHAR in our experimental condition supports the hypothesis that they might contribute in controlling NO levels in root tissues after exposure to NO 3 - [48,49,52].

2-DE maps
The last protein found to be present in higher amount in N plants was a pyruvate decarboxylase (spot 231, PDC). This enzyme catalyzes the decarboxylation of pyruvic acid into acetaldehyde, the first step of the alcoholic fermentation. In particular, we identified the PDC isoenzyme 3 that has been previously found to be induced in hypoxia condition [53]. Although further studies are required to understand why PDC is induced by NO 3 -, we can observe that fermentation pathways are induced in response to redox status changes and that this condition could be also linked to the activation of the Hb/NO cycle (see above) [49,54]. Among the spots identified in roots, six showed a downaccumulation in N plants ( Table 2). Three of them were identified as phenylalanine ammonia-lyase (spots 171, 172 and 1160, PAL-a, PAL-b and PAL-c). The MS analysis indicated for all three spots the same protein [Gen-Bank:AAL40137] while the electrophoretic data showed some differences in M r and pI, suggesting that post-translational modification events may have occurred. It has been shown as low nitrogen availability induces transcripts encoding enzymes of phenylpropanoid and flavonoid metabolism, such as PAL, chalcone synthase and 4coumarate:coenzyme A ligase, whilst after nitrogen repletion these activities are down-regulated [7,55]. Our proteomic data appear to be in agreement with these studies.
Previously, it was found that under low nitrogen availability four proteases (e.g. serine, aspartate/metalloproteases and two cysteine proteases) increased their activity to degrade non-essential proteins in order to remobilize this nutrient [56]. In this work, we found an aspartic protease belonging to the A1 family (spot 707, AP) that was downregulated after NO 3exposure. Moreover, the experimental M r appeared lower with respect to that expected for this protein, thus suggesting that this spot is referable to the active form of the enzyme [57]. These data support a new possible role for A1 protease family [57,58].
Phosphoenolpyruvate carboxylase activity is known to increase during nitrate assimilation, having a role in cell pH homeostasis and an anaplerotic function [14,19,[59][60][61]. In addition, the monoubiquitination of this enzyme was recently well described in germinating castor oil seeds by Uhrig and co-workers [62]. It was found that this event is non-destructive and that this reversible post-translational modification of the enzyme reduces its affinity for PEP and its sensitivity to allosteric activators and inhibi-  [62]. These results suggest that in maize roots the modulation of PEP-Case activity in response to nitrogen availability could occur also through reversible monoubiquitination.
The last spot identified in roots that was down-regulated by NO 3was the 14-3-3-like protein GF14-12 (spot 774, GF14-12). Previously, it was found that this protein is localized in the nucleus where it binds the DNA at the Gbox regions in association with transcription factors and that it is involved in the regulation of gene expression [63,64]. More recently, it was described an interaction of 14-3-3 proteins with some transcription factors such as VP1, EmBP1, TBP and TFIIB [65]. Further studies are required to clarify the effective role of GF14-12, for which the functional information are still lacking.

Functional role and quantitative change of the proteins identified in leaves
Many of the spots identified in leaves by LC-ESI-MS/MS analysis were proteins linked to the NO 3assimilation as well as to the photosynthetic activity ( Table 3).
The activity of NR can be modulated also at post-translational level through a phosphorylation event followed by binding of inhibitory 14-3-3 protein [66,67]. One of the spots analyzed in the leaves was identified as TaWIN2  (Table 3, spot 1094, TaWIN2), that was previously described to be involved in the NR inactivation [67]. We found that the level of this protein decreased in leaves of N plants, where NR activity was induced (Table 3).
According to the well known relationships existing between nitrogen and carbon metabolism, the changes in accumulation of some spots after NO 3addition are consistent with an increase of photosynthesis rate. Two spots that raise after NO 3 addition were identified as CPN-60α and CPN-60β (spot 467 and 462, CPN-60α and CPN-60β, respectively), that are chaperonin proteins involved in folding of ribulose-1,5-bisphosphate carboxylase [68]. Moreover, a chloroplastic triosephosphate isomerase was up-regulated by NO 3 -(spot 1065, TIM), while a cytosolic 6-phosphogluconate dehydrogenase (spot 650, 6PGD) was down-regulated, as expected when the request of reducing power could be satisfied by the increase in photosynthetic activity [69].
Spot 500 was identified as the α subunit of the chloroplastic ATP synthase (ATPsyn α), but unexpectedly it was more abundant in leaves of C plants. Although only a speculative interpretation of this result can be made, we could hypothesize that in leaves of the N plants ATP synthase should be activated and this process requires the reconstitution of the enzymatic complex in the thylakoid membranes [70]. Hence, to clarify this point, it should be necessary to investigate if the decrease of ATPsyn α observed in the soluble fraction of N plants is effectively accompanied by an increase of this protein in the membrane fraction.
Thiamine (i.e. vitamin B 1 ) is required in many pathways, such as the Calvin cycle, the branched-chain amino acid pathway and pigment biosynthesis [71]. Along with higher request of this vitamin in leaves of N plants, where the activation of these pathways could take place, we identified, among the spots up-regulated by N, the thiazole biosynthetic enzyme (spot 999, TH1-1) that is known to be involved in thiamine biosynthesis [71].
Two spots were identified as PEPCase (spot 134 and 138, PEPCase-a and PEPCase-b respectively). In C 4 plants such as maize, this enzyme plays a central role in photosynthesis, because it catalyses the primary fixation of atmospheric CO 2 [72]. The catalytic activity and sensitivity of this enzyme are mediated by a reversible phosphorylation [73]. The experimental pIs of the spots 134 and 138 were 5.8 and 5.7, respectively. Moreover, these two PEPCase forms showed opposite changes in abundance in the leaves of plants grown in the last 30 h in the presence of NO 3 with respect to the controls. The results obtained in our work suggest that the two spots of PEPCase are referable to the phosphorylated (spot 138) and to the unphosphorylated (spot 134) form with a predicted pI of 5.7 and 5.8, respectively, that are known to correspond to the more and less active states of this enzyme [74]. Interestingly, despite the fact that data suggest an increase in the photosynthetic activity, the phosphorylated form was more abundant in the proteomic map of C plants. These results support the immunological observation by Ueno and co-workers [73] that the diurnal regulation of phosphorylation state of PEPCase appears delayed in nitrogenlimited conditions, suggesting that the circadian control of PEPcase is affected by nitrogen starvation.
Two of the spots down-regulated in leaves of N plants were a phenylalanine ammonia-lyase (spot 313, PAL) and a methionine synthase (spot 254, MetS). The decrease of PAL, observed also in root tissue (see above), is a further evidence that phenylpropanoid and flavonoid metabo-lisms are affected by nitrogen availability [7,55]. On the other hand, the change in accumulation of MetS is contrasting with a recent proteomic study performed on wheat by Bahrman and co-workers [39]. These authors found that the induction of this enzyme was positively related to nitrogen availability. This discrepancy could be associated to different genetic traits of the two species, as well as it could be linked to different experimental approaches adopted in the two studies. Nevertheless, it should be observed that in both these works a single spot referable to MetS was detected, while further information on total level and/or on activity of this enzyme is necessary to clarify this point.
The spot 219, which considerably increased in N plants, was a lipoxygenase (LOX). In particular, the analysis of the MS spectra identified the LOX codified by ZmLOX10 gene, which was found to be a plastidic type 2 linoleate 13-LOX [75]. The expression analysis of this gene revealed that its transcript was abundant in leaves and was regulated by a circadian rhythm with a trend strictly linked to the photosynthetic activity. Moreover, it has been proposed that ZmLOX10 is involved in the hydroperoxide lyase-mediated production of C 6 -aldehydes and alcohols and not in the biosynthesis of JA [75]. Although some evidences suggest a role of ZmLOX10 in the responses to (a)biotic stresses, its involvement in the diurnal lipid metabolism was also proposed [75,76].
At the same time, we identified two proteins as an oxygenevolving enhancer protein 2 (spot 1244, OEE2) and a 23 kDa polypeptide of photosystem II (spot 1612, 23pPSII), which were down-accumulated in leaves of N plants (Table 3). Both have been classified as members of PsbP family that is one of the three extrinsic protein families composing the oxygen-evolving complex (OEC) of photosystem II in higher plants [77][78][79]. In addition, it was recently demonstrated that PsbP proteins are essential for the normal function of PSII and play a crucial role in stabilizing the Mn cluster in vivo [80]. Moreover, the stability of this class of protein seems related to the lipid composition of chloroplastic membranes that is also affected by nitrogen availability [81,82].
In order to elucidate the physiological meaning of these variations and to verify if they could be related to a stress status or to an alteration in photosynthetic performance, changes of both maximum quantum yield of photosystem II (F V /F M ; dark adapted plants) and effective quantum yield of photosystem II (Φ II ; light adapted plants), dry weight and MDA levels of shoot were measured ( Figure  5). Although the F V /F M parameter, measured on overnight dark adapted plants at time points 0, 24 and 48 hours, resulted in very similar values between C and N plants (about 0.80; see also Figure 5A), the Φ II values showed a very slight decrease in C plants during the second period of illumination (C plants Φ II , 0.71 versus N plants, 0.73) and the difference became more marked between 48 and 54 hours of nitrogen starvation. Similar data could be obtained by monitoring biomass production at the different time points (Figure 5B), indicating that photosynthetic performances are highly impaired in C plants after 48-54 hours of treatment. Nevertheless, no changes in MDA were detected in all the conditions tested ( Figure 5C).
Taken together these results indicate that at the 30 th h, the time point chosen for proteomic analysis, plants start feeling the different nitrogen content in the growth media without developing major stress symptoms and the associated pleiotropic effects.
These data sustain the hypothesis that ZmLOX10 could be involved in lipid metabolism of the chloroplast that is strictly depending on photosynthetic activity [75,76]. Further analyses are needed to unravel this possible intriguing role of ZmLOX10.
Considering the PsbP proteins, the change in accumulation of OEE2 and 23pPSII could indicate that OEC stability is affected by the N availability. Through time-course experiments, it will be possible to better correlate the relationship among N nutritional status, lipid metabolism, PsbP protein levels and PSII functionality.

Conclusion
Many of the proteins found to change in accumulation in response to NO 3were directly involved in the assimilation of this mineral nutrient. Moreover, the results underline the strict relationship between nitrogen and carbon metabolisms. The experimental design chosen for this proteomic study allows to emphasize some intriguing metabolic activities in both organs. Besides a dramatic increase of NO 3assimilation pathway, the exposure to a high NO 3concentration after a starvation period seems to induce a modification in NO metabolism in roots, that could depend on the need of responding to the new nutritional status. In leaves, many proteins were found to be (in)directly involved in the photosynthesis reactivation and in the maintenance of the chloroplastic functionality.
In addition, this proteomic analysis confirms the modulation by phosphorylation of the PEPCase in the leaves, suggesting that nitrogen availability could affect the circadian rhythms, as well as it shows that the form of this enzyme operating in roots could be modulated by monoubiquitination. Although further efforts are required to elucidate these results, the present study underlines the central role of post-translational events to modulate pivotal enzymes in plant metabolic response to NO 3 -.
F V /F M , Φ II , dry weight and MDA in leaves  Figure 1). Values are the mean ± SE of three independent biological samples analyzed in triplicate (n = 9). Samples indicated with the same letters do not differ significantly according to Tukey's test (p < 0.01).

Plant material and growth conditions
Maize (Zea mays L.) seeds of T250 inbred line, kindly provided by Prof. Zeno Varanini of Udine University -Italy, were germinated in the dark at 26°C on blotting paper saturated with deionized water. After 72 h, seedlings were transferred to a hydroponic system placed in a growth chamber with a day/night regime of 16

Levels of nitrate
Nitrate was extracted from the tissues by homogenizing the samples previously boiled in 4 volumes of distilled water for 15 min. The homogenate was centrifuged at 12,000 g for 20 min to obtain a clarified supernatant. Nitrate content was measured by adding 0.8 ml of 5% (w/ v) salicylic acid in concentrated sulfuric acid solution to 0.2 ml of the supernatant. The mixture was stirred vigorously and allowed to react over 20 min, afterwards 19 ml of 2 N NaOH were slowly added and the resulting colour was read at 410 nm [83].
Nitrate reductase activity NR was extracted by using 4 volumes of ice-cold 50 mM pH 7.8 MOPS-KOH buffer containing 5 mM EDTA, 5 mM NaF, 2 mM MSH, 1 mM PMSF, 10 μM FAD, 1 μM leupeptin and 10 μM chymostatin. The homogenates were centrifuged at 13,000 g for 15 min at 4°C. NR activity was measured as described by Ferrario-Méry et al. [84] using a reaction mixture consisting of 50 mM pH 7.5 MOPS-KOH buffer, 1 mM NaF, 10 mM KNO 3 , 0.17 mM NADH, 10 mM MgCl 2 and 5 mM EDTA. The reaction was blocked after 10 or 20 min by adding an equal volume of sulphanilamide (1%, w/v in 3 M HCl) followed by n-naphtylethylethylenediamine dihydrochloride (0.02%, w/v). 30 min later, the concentration of NO 2 was determined spectrophotometrically at 540 nm. The protein concentration was determined by 2-D Quant Kit (GE Healthcare).

Determination of reducing sugars, sucrose, amino acids, total proteins and chlorophyll
Reducing sugars, sucrose and amino acids were extracted by homogenizing frozen tissues in 4 volumes of ice-cold 0.5 M perchloric acid (PCA). The homogenate was centrifuged for 10 min at 13,000 g at 4°C and the resulting pellet was washed with the same volume of PCA and then centrifuged again in the same conditions. KOH was added to the collected supernatant (to pH 7.6) to remove excess PCA. Reducing sugars were measured according to the colorimetric method by Nelson [85]. Total soluble sugars were determined by the same method boiling an aliquot of PCA extract for 1 h before neutralization. Sucrose was estimated from the difference between total soluble and reducing sugars. Total amino acids were measured by the ninhydrin method [86].
Total proteins were extracted as previously described by Martínez and co-workers [87] by homogenizing the samples, previously powdered in liquid nitrogen, in 4 volumes of a 125 mM pH 8.8 Tris-HCl buffer containing 1% (w/v) SDS, 10% (w/v) glycerol, 50 mM Na 2 S 2 O 5 . The homogenate was centrifuged at 13,000 g for 20 min to obtain a clarified supernatant. The protein content was measured by using 2-D Quant Kit (GE Healthcare).
Chlorophyll was extracted by homogenizing the leaves, previously powdered in liquid nitrogen, in 4 volumes of 80% pre-cooled acetone (v/v). The homogenate was centrifuged at 13,000 g for 20 min at 4°C to obtain a clarified supernatant. Chlorophyll concentration was measured according to Lichtenthaler [88].

Determination of malondialdehyde and chlorophyll fluorescence of the leaves
Malondialdehyde (MDA) was assayed by the method of Heath & Packer [89]. Frozen samples were homogenized with 4 volumes of ice-cold 0.1% (w/v) trichloroacetic acid (TCA) and centrifuged at 13,000 g for 20 min at 4°C. An equal volume of 20% (w/v) TCA plus 0.5% (w/v) thiobarbituric acid was added to the supernatants, which were subsequently heated at 95°C for 30 min. The extracts were then clarified by centrifugation at 13,000 g for 10 min, and the difference between the absorbance at 532 and 600 nm was measured. The MDA equivalent was calculated from the resulting difference using the extinction coefficient of 155 mM -1 cm -1 .
In order to determine the photosynthetic performance, the chlorophyll fluorescence was measured by using a portable continuous-excitation type fluorometer (Handy-PEA, Hansatech Instrument). The maximum quantum efficiency of photosystem II (F V /F M ) was calculated on over-night dark adapted plants, according to the equation (F M -F 0 )/F M , where F 0 and F M are the fluorescence levels when plastoquinone electron acceptor pool (Qa) is fully oxidized and transiently fully reduced, respectively [90]. The photosynthetic performance of light adapted plants was evaluated by monitoring the effective quantum yield of photosystem II (Φ II ) defined as (F M '-F 0 ')/F M ' [91], where F M 'and F 0 ' represent the maximal and minimal fluorescence emission of photosystem II under light conditions.

Statistical analyses of biochemical and physiological measurements
For all the biochemical and physiological measurements, the experimental design consisted in three independent biological samples each analyzed in triplicate (n = 9).
One-way analysis of variance (ANOVA) followed by the post hoc Tukey test (p < 0.01) was used to verify the significance of the variations measured among all the tested parameters. This statistical analysis was performed using the software STATISTICA 7.

Extraction of protein samples for 2-DE analysis
Three independent biological replicates were extracted for each condition. Frozen samples, each composed by leaves or roots of 6 plants, were finely powdered in liquid nitrogen using a pestle and mortar, added with PVPP [0.5% and 1% (w/w) for roots and leaves samples, respectively], homogenized in 4 volumes of extraction buffer [0.5 M Tris-HCl pH 8, 0.7 M sucrose, 10 mM EDTA, 1 mM PMSF, 1 μM leupeptin, 0.1 mg mL -1 Pefabloc (Fluka), 0.2% (v/v) MSH] and centrifuged at 13,000 g at 4°C for 20 min. The resultant supernatant was centrifuged at 100,000 g at 4°C for 38 min to obtain the soluble fraction. Proteins were then purified using the method previously described by Hurkman and Tanaka [92] by adding an equal volume of ice-cold Tris buffered phenol (pH 8) to the supernatant. Samples were shaken for 30 min at 4°C, incubated for 2 h at 4°C and finally centrifuged at 5,000 g for 20 min at 4°C to separate the phases. Proteins, grouped in the upper phenol phase, were precipitated by the addition of five volumes of -20°C pre-cooled 0.1 M ammonium acetate in methanol and the incubation at -20°C overnight. Precipitated proteins were recovered by centrifuging at 13,000 g at 4°C for 30 min and then washed again with cold methanolic ammonium acetate and three times with cold 80% (v/v) acetone. The final pellet was dried under vacuum and dissolved in IEF buffer [7 M urea, 2 M thiourea, 3% (w/v) CHAPS, 1% (v/v) octylphenoxy polyethoxy ethanol (NP-40), 50 mg mL -1 DTT and 2% (v/v) IPG Buffer pH 3-10 (GE Healthcare)] by vortexing and incubating for 1 h at room temperature. Samples were centrifuged at 10,000 g for 10 min and the supernatants stored at -80°C until further use. Protein concentration was determined by 2-D Quant Kit (GE Healthcare).
Proteins were stained using the colloidal Coomassie Brilliant Blue G-250 (cCBB) procedure, as previously described by Neuhoff and co-workers [94]. The gels were scanned in an Epson Expression 1680 Pro Scanner and analyzed with ImageMaster 2-D Platinum Software (GE Healthcare). Automatic matching was complemented by manual matching. The molecular weights of the spots were deduced on the basis of the migration of SigmaMark-ers™ wide range (MW 6.500 -205.000), while pIs were determined according to the strip manufacturer's instructions (GE Healthcare) reporting on the reference gel of the software-assisted analysis the values of pI predicted for any given length of the strip. Both M r and pI of the spots of interest were then determined by using software-automated algorithm.
Relative spot volumes (%Vol) of the six replicate gels per condition were compared and were analyzed according to the Student's t-test to verify whether the changes were statistically significant (p < 0.05). This analysis was performed by using SigmaStat software. Only spots showing at least a two-fold change in their relative volumes were considered for successive analyses.

Protein in-gel digestion and LC-ESI-MS/MS analysis
Spots excised from gels stained with cCBB were digested as described by Magni and co-workers [95] with some refine-ments. In detail, after the destaining procedure, spots were dried under vacuum on a centrifugal evaporator and incubated in 10 mM DTT, 100 mM NH 4 HCO 3 for 45 min at 56°C. The solution was replaced with 55 mM iodoacetamide, 100 mM NH 4 HCO 3 and the spots were incubated for 30 min in the dark at room temperature. After that, spots were briefly washed with 100 mM NH 4 HCO 3 and again incubated for 15 min in 50% (v/v) acetonitrile (ACN), for 3 min in 100% ACN, for 3 min in 100 mM NH 4 HCO 3 , for 15 min in 50 mM NH 4 HCO 3 in 50% (v/v) ACN and finally dried under vacuum. The following phases consisting in the protein digestion with trypsin [Sequencing grade modified Trypsin V5111, Promega, Madison] and in the recovery of peptides were carried out as described in the article above cited.
The LC-ESI-MS/MS experiments were conducted using a Surveyor (MS pump Plus) HPLC system directly connected to the ESI source of a Finnigan LCQ DECA XP MAX ion trap mass spectrometer (ThermoFisher Scientific Inc., Waltham, USA). Chromatography separations were obtained on a BioBasic C18 column (180 μm I.D × 150 mm length, 5 μm particle size), using a linear gradient from 5% to 80% solvent B [solvent A: 0.1% (v/v) formic acid; solvent B: ACN containing 0.1% (v/v) formic acid] with a flow of 2.5 μl/min. ESI was performed in positive ionization mode with spray voltage and capillary temperature set at 3 kV and at 220°C, respectively. Data were collected in the full-scan and data dependent MS/MS mode with collision energy of 35% and a dynamic exclusion window of 3 min. Spectra were searched by TurboSEQUEST ® incorporated in BioworksBrowser 3.2 software (ThermoFisher Scientific Inc., Waltham, USA) against the Zea mays protein subset, Zea mays EST subset and against the protein NCBI-nr database, all downloaded from the National Center for Biotechnology Information [96]. The searches were carried out assuming parent ion and fragment ion mass tolerance of ± 2 Da and ± 1 Da, respectively, two possible missed cleavages per peptide, fixed carboxyamidomethylation of cysteine and variable methionine oxidation. Positive hits were filtered on the basis of peptides scores [Xcorr ≥ 1.5 (+1 charge state), ≥ 2.0 (+2 charge state), ≥ 2.5 (≥ 3 charge state), ΔCn ≥ 0.1, peptide probability < 1 × 10 -3 and Sf ≥ 0.70] [97]. If needed, identified peptides were used in protein similarity search performed by alignment analyses against the NCBI-nr database using the FASTS algorithm [98]. Physical properties of the characterized proteins were predicted by in silico tools at ExPASy [99].