Genetic regulation by NLA and microRNA827 for maintaining nitrate-dependent phosphate homeostasis in arabidopsis.

Plants need abundant nitrogen and phosphorus for higher yield. Improving plant genetics for higher nitrogen and phosphorus use efficiency would save potentially billions of dollars annually on fertilizers and reduce global environmental pollution. This will require knowledge of molecular regulators for maintaining homeostasis of these nutrients in plants. Previously, we reported that the NITROGEN LIMITATION ADAPTATION (NLA) gene is involved in adaptive responses to low-nitrogen conditions in Arabidopsis, where nla mutant plants display abrupt early senescence. To understand the molecular mechanisms underlying NLA function, two suppressors of the nla mutation were isolated that recover the nla mutant phenotype to wild type. Map-based cloning identified these suppressors as the phosphate (Pi) transport-related genes PHF1 and PHT1.1. In addition, NLA expression is shown to be regulated by the low-Pi induced microRNA miR827. Pi analysis revealed that the early senescence in nla mutant plants was due to Pi toxicity. These plants accumulated over five times the normal Pi content in shoots specifically under low nitrate and high Pi but not under high nitrate conditions. Also the Pi overaccumulator pho2 mutant shows Pi toxicity in a nitrate-dependent manner similar to the nla mutant. Further, the nitrate and Pi levels are shown to have an antagonistic crosstalk as displayed by their differential effects on flowering time. The results demonstrate that NLA and miR827 have pivotal roles in regulating Pi homeostasis in plants in a nitrate-dependent fashion.

Introduction

High yielding crops require the application of large amounts of nitrogen (N) and phosphorus (P) fertilizers. However, most of the crop plants are able to take up less than 40% of the applied N and P fertilizers and the rest of it is lost to the environment. This leads to an increase in crop production cost and significant global environmental damage by eutrophication of marine and fresh water ecosystems and gaseous loss to the atmosphere [1]–[3]. For instance, a 1% increase in N use efficiency worldwide would save ∼$1.1 billion annually. In addition, a ∼2.5-fold increase in N- and P-driven eutrophication of water bodies is expected by the year 2050 given current trends [1]. Therefore, developing crop varieties with higher nutrient use efficiency to restrict the excessive use of N and P fertilizer is required. For this, a comprehensive knowledge of molecular mechanisms regulating N and P homeostasis in plants is a prerequisite.

P is an essential structural component of nucleic acids and phospholipids and is a key constituent of high energy phosphate compounds such as ATP and ADP [2]. Despite the integral role of P for normal plant growth, development and yield, P availability in soil is usually the lowest of the macronutrients [4]. Even though the total P content in soil is high, its availability for plant uptake is largely restricted due to its adsorption in soil, precipitation by other cations and conversion into organic forms by microbes [2], [5]. To maintain internal P homeostasis, plants have evolved a series of adaptive responses that include induction of inorganic phosphate (Pi) transporters, change of root architecture, secretion of phosphatase and symbiosis with mycorrhizal fungi [1], [5], [6]. Genetic and molecular approaches have revealed several genes involved in Pi transport, homeostasis and adaptive responses in plants [6]. Among the regulatory genes for Pi homeostasis, PHOSPHATE STARVATION RESPONSE1 (PHR1) is a MYB transcription factor which acts in the Pi starvation signaling pathway by regulating a group of Pi starvation induced genes [7]. PHR1 has sumoylation sites and is a target of the SUMO E3 ligase SIZ1 which is a controller of Pi starvation dependent responses in Arabidopsis [8]. A SEC12-related PHOSPHATE TRANSPORTER TRAFFIC FACILITATOR1 (PHF1) gene facilitates the trafficking of a key high affinity Pi transporter, PHOSPHATE TRANSPORTER1.1 (PHT1.1), which is involved in Pi acquisition [9], [10]. The pht1.1 mutant shows reduced Pi uptake in Arabidopsis [10] and its overexpression increased the Pi uptake in tobacco cells [11]. Another Pi transporter PHT1.4 might have a role in Pi uptake under high Pi conditions and a pht1.1 x pht1.4 double mutant shows a significant reduction in Pi uptake and shoot Pi content [10]. The low affinity Pi transporter PHT2.1 facilitates Pi allocation within the plant between roots and shoots and is required for Pi remobilization in old and young leaves [12]. Twenty genes with SPX domains (SYG1, Pho81 and XPR1) respond to Pi levels and SPX1 and SPX3 are proposed to regulate the expression of Pi starvation response genes [13]. Another SPX-domain containing gene PHO1 has a role in Pi loading into the xylem and possibly also in Pi signaling [14], [15]. The pho1 mutant has significantly lower levels of shoot Pi, but has normal root Pi content [14], [16]. A mutation in a ubiquitin conjugase gene PHO2 results in overaccumulation of Pi that causes Pi toxicity in Arabidopsis [17], [18]. The PHO2 gene is a target of the microRNA, miR399 [18], [19]. Like other nutrient elements, uptake of Pi also depends upon external pH and ion competition. At low pH, Pi in the H2PO4− form is present at a high proportion whereas at high pH HPO42− dominates [2]. Competition between Pi and arsenate is well known and both are taken up by the same transport system [2], [9].

N components are not only required for plant growth, but also serve as regulators of various metabolic and developmental pathways. Plants take up N mainly as nitrate (NO3−) and ammonium (NH4+), with NO3− being the predominant form in most agricultural soils [20]. NO3− and Pi are the two most important anions required for plant growth and development. However, the interaction and balance between NO3− and Pi in plants is not well studied. Here, we demonstrate the crosstalk between NO3− and Pi and the role of genes encoding an E3 ubiquitin ligase, NITROGEN LIMITATION ADAPTATION (NLA) and an E2 ubiquitin conjugase, PHO2, and a microRNA- miR827, in maintaining Pi homeostasis in Arabidopsis thaliana in a NO3− dependent manner.

Results

Identification of the nla Mutant and nla-Suppressors, phf1 and pht1.1

The nla mutant was identified based on its altered growth response to N limitation. nla mutant plants failed to show several adaptive responses to low-N conditions, such as the inability to accumulate anthocyanin and abrupt early senescence compared to wild type (WT) plants [21]. The nla phenotype was specific to the low-N growth condition, since at optimum-N and under other abiotic stresses, both nla and WT plants were similar in phenotype [21]. Later, it was revealed that in nla mutant plants under low-N conditions, the substrate in phenylpropanoid pathway was channeled towards lignin biosynthesis instead of anthocyanin synthesis, resulting in low anthocyanin accumulation [22]. Anthocyanin protects plants from photoinhibition damage under various stress conditions [22], [23]. A lack of anthocyanin accumulation in the nla mutant under low-N conditions means an absence of the photo-protective screen, which might result in the early senescence phenotype in nla mutant plants.

This evidence suggests that the NLA gene regulates the Arabidopsis adaptive responses under low-N. The knowledge of N regulatory genes is limited and the NLA gene might be a key component for the regulation of plant adaptation to low-N. Therefore, we initiated further studies to understand in more detail the physiological and molecular role of NLA. One approach was to generate and identify suppressor mutations in the nla mutant, which would restore the nla mutant phenotype to WT. For this nla mutant seeds were chemically mutagenised and screened at low-N supply. Two suppressor plants (nla-sup1 and nla-sup2) were identified which were phenotypically similar to WT in that the nla-suppressor plants did not show early senescence under low-N conditions (Figure 1A, 1B and Figure S1A, S1B). These suppressor genes were identified separately using map-based cloning approaches. The nla-sup1 locus was mapped to the lower arm of chromosome 3, in a ∼132 kb region with 40 annotated genes (Figure 1C). Sequencing this genomic region for each gene revealed that in PHF1 (At3g52190) a transition of a single G/C to A/T occurred at the 3′ end of intron 7 (position 2218). This transition destroyed the conserved dinucleotide, ‘AG’, which is required for proper splice recognition at the 3′ end of an intron. The resulting altered splice site led to the inclusion of the adjacent exonic ‘G’ as part of the intron (Figure 1C) leading to a 1 bp deletion in the processed mRNA. The cDNA sequence also revealed that in the nla-sup1 one ‘G’ nucleotide was missing compared to WT which causes a frameshift and generation of a premature stop codon resulting in a truncated PHF1 protein sequence (Figure 1C). The second suppressor, nla-sup2, was mapped to a location on chromosome 5. Sequencing of the genomic regions confirmed a mutation in a single G/C to A/T (position 229) in the first exon of PHT1.1 (At5g43350) gene. This transition results in a switch from alanine to threonine (position 73) in the amino acid sequence of the PHT1.1 protein (Figure 1D). PHT1.1 is a high affinity Pi uptake transporter [10] and PHF1 facilitates the trafficking of PHT1.1 [9]. To confirm that these altered genes corresponded to the suppressor mutations, their respective T-DNA insertion mutants were used to make nla x phf1 and nla x pht1.1 double mutants. Like the suppressor mutations, these double mutants also recovered the nla mutant phenotype to WT.

Figure 1

Recovery of the nla mutant phenotype by nla-suppressor1 and mapping of the nla-suppressors, PHF1 and PHT1.1.

(A) WT, nla and nla-sup1 (nla-phf1) plants grown at 10 mM NO3−-10 mM Pi. (B) Plants grown at 3 mM NO3−-10 mM Pi. (C) Mapping and partial genomic sequence of WT and nla-sup1 (phf1). The transition at 3′ end nucleotide ‘g’ to ‘a’, indicated by red arrow, disrupted the conserved dinucleotide, AG (underlined), which is required for splicing. In the phf1 mutation, the next available ‘G’ from the 5′ end of exon became part of the intron to form ‘AG’, causing loss of one G from the phf1 coding sequence. This cryptic splice site resulted in a frameshift and creation of a premature stop codon. (D) The mutation in the second nla suppressor was missense in PHT1.1 gene. In the first exon, nucleotide ‘G’ at position 229 was switched to ‘A’ which changes the codon GCC (encodes alanine) to ACC (for threonine).

MicroRNA827 Suppresses the NLA Transcript Level

MicroRNAs (miRNAs) are noncoding small RNAs and suppress the expression of genes that have nearly complementary sequences by mRNA cleavage [24], [25]. The NLA gene is a putative target of a miRNA, miR827 (At3g59884) based on complementary sequences (Figure 2A and 2B). By identifying phf1 and pht1.1 as nla-suppressors, we have shown that NLA directly or indirectly targets PHF1 and PHT1.1. To know how NLA expression is regulated and whether the NLA transcript is a target of miR827, we analyzed the expression pattern of NLA and miR827 under different N and Pi regimes. Also the expression of NLA transcript was analyzed in the miR827 overexpresser (OX) and T-DNA mutant lines. The expression of miR827 is up-regulated by low-Pi conditions (Figure 2C) [26]. In contrast, the transcript level of NLA was down-regulated under low-Pi, the exact opposite response to the expression of miR827 (Figure 2D). Figure 2E shows that miR827 was overexpressed by 35-times in the OX-miR827 line and there was very low expression in the miR827-mutant plants as compared to WT. The NLA transcript level was over 3-fold down-regulated in the OX-miR827 and over 2-fold up-regulated in miR827 mutant plants compared to WT (Figure 2F). These results show that NLA is a target of miR827. Similarly, the PHO2 gene has been reported to be a target of miR399 and their expression level is also Pi dependent [18], [19].

Figure 2

Regulation of NLA expression by miR827.

(A) NLA gene organization, exons (black boxes), introns (dotted lines), and UTR (empty boxes). Arrow indicates putative target site of miR827. (B) Partial sequences of NLA and miR827, the complementary sequences are underlined. (C) Relative expression of miR827, and (D) NLA in WT plants grown at HNHP (10 mM NO3−-10 mM Pi), LNHP (3 mM NO3−-10 mM Pi), HNLP (10 mM NO3−-0.5 mM Pi), and LNLP (3 mM NO3−-0.5 mM Pi). (E) Relative expression of miR827, and (F) NLA in WT, miR827-overexpressor, and miR827-mutant plants grown at 10 mM NO3−-10 mM Pi conditions. Expression levels were determined by real-time PCR. Shown are mean ± SD.

(A) NLA gene organization, exons (black boxes), inpan>trons (dotted linpan>es), and UTR (empty boxes). Arrow inpan>dicates putative target site of miR827. (B) Partial sequences of NLA and miR827, the complementary sequences are underlined. (C) Relative expression of miR827, and (D) NLA in WT plants grown at HNHP (10 mM NO3−-10 mM Pi), LNHP (3 mM NO3−-10 mM Pi), HNLP (10 mM NO3−-0.5 mM Pi), and LNLP (3 mM NO3−-0.5 mM Pi). (E) Relative expression of miR827, and (F) NLA in WT, miR827-overexpressor, and miR827-mutant plants grown at 10 mM NO3−-10 mM Pi conditions. Expression levels were determined by real-time PCR. Shown are mean ± SD.

The nla Mutant Overaccumulates Pi in Shoots Specifically at Low NO3− Supply

It has been shown previously that phf1 and pht1.1 mutants accumulate less Pi than WT [9], [10]. Therefore, the identification of the nla-suppressors as phf1 and pht1.1 mutations (Figure 1) and the transcript change of NLA with Pi levels (Figure 2C) indicate that Pi accumulation in the nla mutant might be impaired. Pi analysis reveals that under sufficient N-P conditions of 10 mM NO3−-10 mM Pi, where the nla mutant is phenotypically similar to WT (Figure 1A), Pi content in the nla mutant shoots was ∼1.8-fold higher than WT (Table 1). Interestingly, under the relatively low NO3− (3 mM NO3−-10 mM Pi) regime where the nla mutant shows early senescence (Figure 1B), Pi content in the nla mutant shoots increased ∼6.6-fold compared to WT (Table 1). P usually makes up ∼0.2% of plant dry matter and visual symptoms of Pi toxicity appear when P constitutes >1% of dry matter [2], [27]. Pi content in the nla mutant shoots at 3 mM NO3−-10 mM Pi regime makes up ∼2% of dry matter and its senescence phenotype likely is due to Pi toxicity, which leads to chlorosis or necrosis starting from the leaf margins (Figure 1B). Whereas, the nla mutant grown at low-Pi regimes accumulated Pi below P toxicity limits (Table 1). Also the nla suppressors, nla-sup1 and nla-sup2 as well as OX-NLA had a visual phenotype and Pi content similar to WT (Figure 1, Figure S1 and Table 1). Interestingly, the previously described Pi overaccumulator, pho2 mutant [17], [28], shows a phenotype and Pi content similar to the nla mutant (Table 1, Figure S1C and S1D). Pi toxicity in the nla and pho2 mutants was evident only when NO3− supply was relatively low with plants grown on 3 mM NO3−-10 mM Pi, while at the 10 mM NO3−-10 mM Pi regime the plants were phenotypically similar to WT (Figure S1C and S1D). This led to the obvious assumption that NO3− is also playing a role in Pi accumulation. This crosstalk is evident when Pi supply was constant and NO3− supplies were variable, where Pi content in WT shoots increased with decreasing NO3− applications (Table 1 and Table 2). In the nla and pho2 mutants, Pi accumulation accelerated and the appearance of Pi toxicity symptoms occurred earlier in accordance with decreasing NO3− supply at a given Pi level (Table 2). Pi toxicity in these mutants occurred only under low NO3− availability. This can be ascribed to the additive effects of lack of negative regulation by NLA or PHO2 and minimal suppression of NO3− on Pi accumulation.

Table 1

Pi content in Arabidopsis shoots (nmole/mg fresh weight).

Treatment	WT	nla	pho2	nla-sup1	nla-sup2	OX-NLA
NO3−-Pi (mM)
10-10	17.3±1.9gh	30.9±3.1bc	35.2±3.2b	10.6±1.1i	16.9±1.8gh	15.6±1.4h
3-10	28.7±2.9cd	191±23.1a	203±22.3a	19.9±1.7fg	36.4±4.1b	24.2±2.1de
10-0.5	3.1±0.3m	7.3±0.9k	8.4±0.8jk	2.6±0.3n	3.8±0.3l	2.7±0.2n
3-0.5	9.1±1.1ij	21.4±1.9ef	24.2±2.3de	7.2±0.7k	10.6±1.2i	8.2±0.9jk

Pi content was measured at 27 days after sowing. Shown are mean ± SD. Values with different letters indicate significant difference at P<0.05.

Table 2

Pi content (nmole/mg fresh weight) and Pi toxicity in Arabidopsis shoots.

Treatment NO3−-Pi (mM)	WT	nla	pho2	Pi toxicity in nla and pho2
	Pi content
10-10	16.8±1.4kl	30.2±3.1fg	34.9±3.5ef	No
3-10	28.1±2.2gh	189±17.4b	199±21.2b	24 DAS
1-10	32.1±3.1fg	246±21.2a	254±24.2a	19 DAS
10-3	12.1±1.1m	22.4±2.1i	24.8±2.3hi	No
3-3	20.4±1.9ij	75.4±6.8d	80.4±8.2d	26 DAS
1-3	23.3±2.0i	212±21.2ab	220±21.8ab	20 DAS
10-1	7.5±0.8n	15.7±1.6kl	17.8±1.9jk	No
3-1	15.3±1.2l	35.4±3.3ef	39.6±3.5e	No
1-1	18.7±1.5jk	112±10.4c	116±11.5c	25 DAS

Pi content was measured at 27 days after sowing (DAS). Pi toxicity was recorded as the day when leaves of >50% plants start to show Pi toxicity symptoms. The WT plants did not show Pi toxicity. Shown are mean ± SD. Values with different letters indicate significant difference at P<0.05.

The Pi toxicity symptom inpan> the nla mutant appears mainly after bolting [21] and the above results are from plants at 27 days after sowing (DAS) when plants were already flowering. Therefore, Pi contents at the initiation of bolting (22 DAS) and before bolting (17 DAS) were also analyzed (Table S1). Pi accumulation at 22 and 17 DAS followed a similar trend to that of 27 DAS under varying NO3−-Pi regimes. However, net Pi accumulation in shoots at a specific NO3−-Pi regime was much lower at 17 DAS than 27 DAS (Table S1 and Table 1) and the nla mutant did not show Pi toxicity before bolting. This indicates that the accumulation of Pi accelerates with the transition of Arabidopsis from the vegetative to the reproductive stage. The roots are the initial sites of contact with Pi, but root growth and architecture are similar between nla and WT at varying NO3−-Pi regimes [21]. In conjunction with this, at low NO3− availability the fold differences for Pi content in roots of the nla mutant vs WT was much less than that observed in shoots (Table S2 and Table S3). This suggests that most of the Pi taken up by roots was not retained by them, but was transported towards the shoots.

A radiolabelled 33Pi uptake assay was conducted to assess whether the overaccumulation of Pi inpan> the nla and pho2 mutants was due to increased Pi uptake. Indeed, the nla and pho2 mutants showed a higher 33Pi uptake rate than WT, nla-sup1, nla-sup2 and OX-NLA plants (Figure 3). NLA is an E3 Ubiquitin Ligase and has two distinct domains, SPX and RING. The nla is a deletion mutation (Col background) with a disrupted RING domain [21]. To ascertain if a different type of nla mutation would exert similar effects, a T-DNA insertion mutant in the same gene but in the ecotype Wassilewskija (Ws) (FLAG_352A03) was obtained. The nla mutant (Ws background) showed a similar Pi toxicity phenotype and Pi content as that of the nla mutant (Col background) under low NO3−-high Pi conditions (Figure S3).

Figure 3

33Pi uptake activities in Arabidopsis plants.

Shown are mean ± SD (n = 10). (FW) Fresh weight.

For the detailed analysis of Pi homeostasis maintenance and the complementation study, several overexpresser (Table 3A) and double mutant (Table 3, Double Mutants) lines were generated. As shown above, miR827 suppresses NLA expression (Figure 2F). Indeed, Pi content in the OX-miR827 plants was ∼6.4-fold higher than WT. Pi toxicity and Pi content were similar in OX-miR827 and nla mutant plants. Whereas, the miR827 mutant had Pi content similar to the OX-NLA plants (Table 3). The PHF1 gene facilitates the trafficking of the key Pi uptake transporter, PHT1.1, and their mutant lines had lower Pi content than WT (Table 3C). Hence, their overexpression lines would be expected to accumulate higher Pi and might show Pi toxicity, as observed in the nla and pho2 mutants. However, the increase in Pi content in OX-PHF1 and OX-PHT1.1 lines was only ∼1.5-fold compared to WT, much lower than the Pi toxicity limit (Table 3A). Further, to confirm that phf1 and pht1.1 are nla-sup1 and nla-sup2 suppressor mutations, respectively, the WT cDNAs of PHF1 and PHT1.1 driven by the cauliflower mosaic virus 35S promoter (35S) were transformed into nla-sup1 and nla-sup2 mutants, respectively. The transformed nla-sup1-OX-PHF1 and nla-sup2-OX-PHT1.1 plants (generated using the genomic sequences of PHF1 and PHT1.1, respectively) grown at low NO3− had Pi content similar to the nla mutant (Table 3). In addition, the nla mutant was crossed with T-DNA mutants of phf1 and pht1.1. The double mutants, nla x phf1 and nla x pht1.1 had Pi contents similar to nla-sup1 and nla-sup2, respectively (Table 3, Double Mutants and Control Plants). The nla and pho2 plants were transformed with 35S:PHF1 and 35S:PHT1.1. However, these OX lines had no further increase of Pi contents over nla and pho2 mutants (Table 3). The nla x pho2 double mutant had no significant additive increase of Pi content compared to nla or pho2 mutants (Table 3). This indicates overlapping roles of nla and pho2 in Pi overaccumulation and their involvement in the same functional pathway. This is further confirmed in that pho2 x phf1 and pho2 x pht1.1 also restored Pi content to WT levels, like that of nla x phf1 and nla x pht1.1, respectively (Table 3, Double Mutants). PHT1.4 and PHT2.1 are also Pi transporters, but their mutants had Pi content similar to WT under high Pi conditions (Table 3, Control Plants), suggesting they may have roles in Pi reallocation within plants [6], [10], [12]. In agreement, nla x pht1.4 and nla x pht2.1 double mutants accumulated high Pi levels (Table 3, Double Mutants). The PHO1 gene is involved in Pi loading to the xylem [14] and the pho1 mutant shoots have very low Pi content (Table 3, Control Plants) [14]. The nla x pho1 double mutant gave expected results of having very low shoot Pi content compared to WT (Table 3, Double Mutants), showing that the pho1 mutation is able to restrict Pi loading for transport from roots towards shoots in the nla mutant. PHR1 regulates expression of some Pi-responsive genes [7]. The phr1 mutant had lower shoot Pi content than WT. Also the nla x phr1 mutant had a decrease in Pi content compared to the nla mutant (Table 3, Double Mutants and Control Plants), but not enough to avoid Pi toxicity. SIZ1 gene encodes a SUMO E3 ligase, mutation of which had slightly higher Pi content than WT (Table 3, Control Plants) [8]. However, nla x siz1 had no further increase in Pi levels over the nla mutant (Table 3, Control Plants).

Table 3

Inorganic phosphate (Pi) and total phosphorus (P) contents in Arabidopsis shoots (nmole/mg fresh weight).

Genotype	Pi	Total P
Overexpresser lines
WT OX-NLA	23.1±2.0ef	57.9±6.2ef
WT OX-miR827	175±18.2ab	419±38.2ab
WT OX-PHF1	42.7±3.5c	107±11.4c
WT OX-PHT1.1	40.2±3.9cd	101±10.1cd
nla-sup1 OX-PHF1	180±16.8a	437±46.2a
nla-sup1 OX-NLA	10.2±1.1i	33.8±3.1h
nla-sup2 OX-PHT1.1	167±15.2ab	416±39.8ab
nla-sup2 OX-NLA	18.4±1.9g	47.2±5.1fg
nla OX-PHF1	189±20.1a	450±43.7a
nla OX-PHT1.1	188±19.8a	446±42.2a
pho2 OX-PHF1	192±21.2a	453±40.4a
pho2 OX-PHT1.1	190±20.1a	450±42.5a
Double mutants
nla x pho2	200±19.1a	474±50.2a
nla x phf1	21.2±1.9fg	57.5±5.5ef
nla x pht1.1	37.5±3.2cd	95.2±8.9cd
nla x pht1.4	174±18.2ab	425±46.5ab
nla x pht2.1	171±17.4ab	417±38.4ab
nla x pho1	10.2±1.2i	31.1±2.7h
nla x phr1	146±15.9b	357±32.3b
nla x siz1	179±18.8ab	427±40.2ab
pho2 x phf1	23.4±2.1ef	62.2±6.6e
pho2 x pht1.1	35.2±2.9cd	91.5±10.2cd
Control plants
WT	27.4±2.2e	67.8±6.6e
nla	186±19.2a	441±42.2a
pho2	197±21.2a	461±48.4a
mut-miR827	23.7±1.9ef	59.1±5.5ef
nla-sup1	20.7±2.0fg	53.6±5.1ef
nla-sup2	39.7±4.1cd	99.6±10.9cd
phf1	12.8±1.2h	41.2±3.9g
pht1.1	17.2±1.9g	46.2±5.1fg
pht1.4	26.6±3.1e	65.9±6.2e
pht2.1	25.6±2.8ef	66.2±7.1e
pho1	7.5±0.9j	24.3±2.2i
phr1	20.6±1.6fg	52.6±5.6ef
siz1	33.9±2.7d	86.1±9.1d

Plants were grown at 3 mM NO3−-10 mM Pi and harvested at 27 days after sowing. Shown are mean ± SD. Values with different letters in each column indicate significant difference at P<0.05.

Plants were grown at 3 mM n class="Chemical">NO3−-10 mM Pi and harvested at 27 days after sowinpan>g. Shown are mean ± SD. Values with different letters inpan> each columnpan> inpan>dicate significant difference at P<0.05.

NO3− and Pi Have an Antagonistic Interaction

Accumulation of Pi increased with decreasing NO3− supply (Table 1). The N content was analyzed to confirm whether the changing Pi application would also exert similar suppression on N accumulation. In fact, N content increased with decreasing Pi supply (Table S4). However, the suppression by Pi on N accumulation was quite low (Table S4) compared to the suppression by NO3− on Pi accumulation (Table 1). The differential suppression by these ions could be ascribed to their contrasting mobility in soil, given that NO3− diffusion is 3–4 times faster than Pi [29]. N is available in both anionic and cationic forms, NO3− and NH4+, respectively, with NO3− being dominant in most soils and also the preferred form taken up by most plants, including Arabidopsis [20]. To analyze whether the antagonistic interaction of Pi was generalized to both NO3− and NH4+, the plants were grown in different NH4+-Pi regimes. Surprisingly, even at high NH4+ supply, Pi accumulation was much higher in nla and pho2 mutants (Table 4) and these mutant plants showed Pi toxicity (Figure S4). In contrast to this, equimolar concentration of NO3− suppressed Pi accumulation below toxicity limits in the nla and pho2 mutants (Table 1, Figure 1 and Figure S1). This suggests that Pi has an antagonistic crosstalk with NO3−, but not with NH4+.

Table 4

Pi content in Arabidopsis shoots grown under different NH4+ and Pi conditions.

Treatment NH4+-Pi (mM)	WT	nla	pho2
	nmole/mg fresh weight
10-10	28.8±3.1cd	171±18.1ab	177±18.8ab
3-10	30.2±3.4c	201±21.4a	206±23.2a
10-3	23.6±2.4d	148±15.8b	156±16.4b
3-3	24.5±2.7cd	162±16.8b	166±17.2ab

Pi content was measured at 27 days after sowing. Shown are mean ± SD. Values with different letters indicate significant difference at P<0.05.

Another important finding was the interaction between the NO3− and Pi applications for their effect on flowerinpan>g time inpan> Arabidopsis. For this, the plants were grown under varying relative supplies of NO3− and Pi as described in the materials and methods. Plants flowered earlier when grown under low NO3− compared to high NO3− conditions. In contrast, growth with low-Pi delayed flowering compared to high Pi. Further, the interaction between these two nutrients had a cumulative effect given that the plants grown at 3 mM NO3−-10 mM Pi flowered significantly earlier than plants at 10 mM NO3−-3 mM Pi (Figure 4A, Figure S2). In addition, the nla mutant, which had higher Pi accumulation, flowered significantly earlier compared to the phf1 mutant which had low Pi accumulation (Figure 4A). Expression levels of various genes which regulate flowering time in Arabidopsis were also changed in a similar fashion by NO3−-Pi levels. FLOWERING LOCUS C (FLC) is a floral repressor [30] and its expression was lowest at 3 mM NO3−-10 mM, where plants flowered earliest. Further, the expression of FLC was lower in nla mutant leaves than in phf1 (Figure 4B). FLOWERING LOCUS T (FT), LEAFY (LFY) and APETALA1 (AP1) are positive regulators of flowering time [30], and their expression was significantly higher in leaves of plants grown on the 3 mM NO3−-10 mM Pi than at 10 mM NO3−-3 mM Pi. Also the expression of these genes was higher in the nla than in the phf1 mutant (Figure 4C–4E). The effect of varying NO3− and Pi regimes on change of expression of FLC, FT, LFY and AP1 genes indicates that the autonomous and gibberellin pathways for flowering are affected. Since, Pi and gibberellin levels are positively correlated in Arabidopsis [31] and gibberellin is known to promote flowering [32], which was reflected by early flowering at high Pi content (Figure 4A).

Figure 4

Effect of nitrate and phosphate regimes on flowering in Arabidopsis.

(A) Flowering time (days after sowing), (B, C, D, and E) The relative expression of FLC, FT, LFY, and AP1 genes, respectively, in leaves of WT, nla and phf1 plants. Total RNA was extracted from leaves of plants harvested at 18 days after sowing. Relative transcript levels were determined by real-time PCR. Shown are mean ± SD. Bars with different letters indicate significant difference at P<0.05.

Discussion

The identification of the nla suppressors as phf1 and pht1.1 and analysis of Pi levels in the nla mutant grown under different N and Pi regimes demonstrated that NLA has a role in regulating Pi homeostasis in plants. Further, the expression of NLA itself is regulated by a miR827 in a Pi dependent manner. Interestingly, comparing the phenotype of the nla and pho2 mutants has shown a similar phenotypic pattern under varying NO3− and Pi regimes and both mutants overaccumulate Pi under low NO3− and high Pi conditions. Figure 5 shows the proposed model for the maintenance of Pi homeostasis via the regulation by the NLA and PHO2 genes and the crosstalk between NO3− and Pi. The NLA and PHO2 genes, along with higher NO3− applications, suppress Pi uptake in Arabidopsis WT plants. Plants grow normally as long as both NO3− and Pi supplies are sufficient, whereas, if either of them is limiting, plants start to accumulate anthocyanin as a stress symptom (Figure 5A). In the nla or pho2 mutants, uptake of Pi is higher than in WT. However, when suppression by NO3− is also minimal when plants are grown under low NO3− and high Pi conditions, Pi accumulation is accelerated to toxic levels leading to necrotic senescence (Figure 5B). Thus, these mutants can grow normally even when Pi supply is sub-optimal, without showing Pi deficiency symptoms (Figure 5B) and would have optimal internal Pi content in shoots (Table 1).

Figure 5

Hypothetical model for the role of NLA and PHO2 genes and crosstalk between NO3− and Pi.

(A) In WT plants, NLA and PHO2 act as repressors for Pi uptake and Pi has a negative interaction with external NO3−. The condition of low NO3− and high Pi is most conducive for higher Pi uptake. (B) In nla or pho2 mutants, the NLA or PHO2 gene is not functional and these plants have higher Pi uptake compared to WT plants. Low NO3− and high Pi conditions accelerate Pi accumulation causing Pi toxicity in these mutants. Further, NO3− and Pi content has antagonistic affects on flowering time. Purple leaves show Pi deficiency leading to anthocyanin accumulation, orange leaves indicate Pi toxicity, solid red circles indicate the Pi uptake system and dotted lines indicate less suppression. The size of the letter ‘Pi’ in leaves and arrow thickness in roots is correlated with Pi content in plants.

NO3− and Pi supply had antagonistic inpan>teractions for their accumulation inpan> plants. However, the suppression by a higher supply of NO3− on Pi accumulation was larger than the suppression by higher supply of Pi on N content (Table 1 and Table S4). This could be due to two reasons, first that NO3− might have preferential uptake over Pi, since the requirement for N (>1% of dry matter) by most plants is usually higher than for P (<0.5% of dry matter). Second, the NO3− uptake system might not require negative regulators, since overaccumulation of N is not toxic to the plants and excess NO3− is stored in the vacuoles [33]. In contrast, NLA and PHO2 serve as negative regulators of the Pi uptake system and are in turn regulated by miRNAs, which is necessary for the plants to avoid Pi toxicity caused by its overaccumulation. Further, the antagonistic interaction between NO3− and Pi was reflected by their inverse effect on flowering time in Arabidopsis, given that the plants flowered early under low-N conditions and it took a longer time to flower under low-Pi conditions (Figure 4A and Figure 5). Delayed flowering by low-Pi level could be via gibberellin signalling, since low gibberellin levels are known to delay flowering [32] and Pi and gibberellin levels are positively correlated [31]. In contrast, the plants flowered early under low NO3− conditions which could be either through a general stress management strategy to complete the life cycle earlier or through direct N affects on the flowering pathways. In agreement with the antagonistic crosstalk between NO3− and Pi for their uptake and inverse effects on flowering time observed here, the contrasting effects of NO3− and Pi on root architecture in Arabidopsis were reported earlier [34]. In that case, the primary root length decreased with increasing NO3− supply and increased with increasing Pi availability.

NLA is an E3 ligase and PHO2 is an E2 conjugase protein. The Arabidopsis genome contains 37 E2 and ∼1300 E3 genes [35] and their proteins have potential roles for ubiquitination of target proteins destined to be degraded by the ubiquitin-26S proteasome. This pathway is employed by plants for the degradation and removal of proteins to maintain optimal growth and development. It would be interesting to study how NLA and PHO2 proteins interact to coordinate signals that specify the target proteins involved in Pi acquisition. PHF1 and PHT1.1 are likely direct or indirect targets of NLA and PHO2. This is supported by our results that mutations of phf1 or pht1.1 in the nla and pho2 mutant background successfully reduced Pi overaccumulation and restored their Pi contents similar to WT. The nla mutant had high Pi accumulation while the OX-NLA plants do not have a significantly lower Pi than does WT (Table 1). The simplest explanation is that once a threshold level of NLA is reached, then increasing NLA does not change phenotype with regards Pi accumulation. In addition, NLA is an E3 ubiquitin ligase which is part of the E1–E2–E3 ubiquitination pathway to target proteins that are degraded by the 26S proteasome pathway. Hence, the co-ordination of both pathways might be necessary for the complete degradation of PHF1 and PHT1.1 to a level which can significantly reduce the Pi content in OX-NLA plants, compared to WT plants. The role of NLA and PHO2 in Pi homeostasis is further supported in that they are regulated by the low-Pi induced miRNAs, miR827 (Figure 2) and miR399 [18], [19], respectively. The miR827 and miR399 are expressed at low levels under high-Pi conditions with the concomitant higher expression of NLA and PHO2. This is necessary to prevent the overloading of Pi into plants. Under low-Pi supply, increased expressions of miRNAs repress NLA and PHO2 transcript levels and the repression of Pi transporters is alleviated (Figure 5). Given that Pi availability in soil is usually limited, this insures that sufficient Pi accumulation occurs in plants.

Over the next 50 years, it will be essential to increase crop yields by almost double to meet the needs of the growing world population. Increasing fertilizer use for higher crop yields would lead to large economic costs and vastly add to the environmental load of crop production. Therefore, the development of crops with improved genetics for nutrient use efficiency is absolutely crucial for the sustainability of crop production. We have demonstrated that NLA and PHO2 (their regulation by miR827 and miR399) have pivotal roles in NO3− regulated control of Pi homeostasis. In fertilized agricultural soils, Pi availability is typically low when compared to NO3−, although more fixed N is required for maximal plant growth. It should be possible by modulating the expression of NLA or PHO2 to develop crop cultivars that better utilize the available P in balance with N supply.

Materials and Methods

Plant Material and Growth Conditions

Arabidopsis thaliana ecotype Columbia was used inpan> the experiments unless otherwise stated. The seeds were stratified at 4°C for 3 d and sown inpan> nutrient-free LB2 soil (SunGro Horticulture Canada Ltd, http://www.sungro.com/). The plants were grown inpan> controlled growth chambers at 16 h light/8 h dark, 23°C day/18°C night, light inpan>tensity 200 µE m−1 s−1 and 65% relative humidity. The nutrient solution was applied once a week for 4 weeks and containpan>ed 2 mM MgSO4, 1 mM CaCl2, 100 µM Fe-EDTA, 50 µM H3BO3, 12 µM MnSO4, 1 µM ZnCl2, 1 µM CuSO4 and 0.2 µM Na2MoO4. N and P were supplied in the form of KNO3 and KH2PO4 (pH 6.0), respectively, with varying concentration as mentioned in the results. The nutrient levels used in the experiments are relative levels, but not absolute levels. To avoid the discrepancies of maintaining the absolute nutrient levels in soil between different experiments, we have previously developed the defined nutrient growth conditions [21]. The application of nutrient solution once every week for 4 weeks to Arabidopsis with varying N supplies showed that 10 mM NO3− is sufficient-N for optimum plant growth and yield, 3 mM NO3− is moderately low-N where plants start to develop adaptive responses required to deal with low-N conditions and 1 mM NO3− is a severe-N limiting condition [21] (Figure S5A). The varying P levels showed that 10 mM Pi is sufficient-P for optimum plant growth and yield, 3 mM Pi is moderately low to sufficient-P with a slight appearance of red color in leaves indicating the start of anthocyanin accumulation, 1 mM Pi is a moderately low-P condition with higher anthocyanin accumulation and significantly lower plant growth and yield and 0.5 mM Pi is a severe-P limiting condition (Figure S5B). Shoots and roots were harvested separately at the indicated days, frozen in liquid nitrogen and stored at −80°C until use. The T-DNA insertion mutant lines for different genes obtained from the Arabidopsis Biological Resource Center and used in the experiments are: phf1 (SALK_144943), pht1.1 (SALK_088586), pht1.4 (SALK_103881), pht2.1 (SALK_094069), siz1 (SALK_065397), phr1 (SALK_067629), and miR827 (SALK_020837).

Isolation of nla Mutant Suppressors and Positional Cloning

The homozygous nla mutant (inpan> Col background) seeds were soaked overnpan>ight at 4°C inpan> 100 mM phosphate buffer pH 7.5 and then treated with 0.4% ethyl methanesulfonate (EMS) for 8 h at room temperature. The seeds were washed thoroughly and sown in LB2 soil supplied with 3 mM KNO3− and 10 mM KH2PO4 and other nutrients as above. The M1 plants were grouped into 120 sets, with each set having a pool of 36–40 plants. The seeds from each set were harvested separately. About 200 M2 plants within each set were grown on 3 mM KNO3− and 10 mM KH2PO4 to screen for putative nla-suppressor plants which would recover the nla mutant phenotype similar to WT. Two putative nla-suppressors (nla-sup1 and nla-sup2) plants were identified and were backcrossed twice with nla-Col mutant plants. Simultaneously, to achieve the nla mutation in Ler genetic background, the homozygous nla-Col plant was crossed with Ler WT plant and the nla-Ler plant obtained from this cross was further backcrossed nine times with Ler WT plant to produce genetically clean nla-Ler. The nla-suppressor-Col plant was crossed with nla-Ler plant. Among the segregating F2 progeny, which was grown at 3 mM KNO3− and 10 mM KH2PO4, ∼500 plants showing the WT phenotype were selected for mapping. The first round of mapping was performed as described by Luckowitz et al. [36]. The second round of fine mapping was accomplished using the simple sequence length polymorphism (SSLP) markers prepared from the Arabidopsis genome sequence (www.arabidopsis.org). The mapping schemes are shown in Figure 1. We identified two independent putative nla-suppressors, with the mapping for each done separately. One suppressor, nla-sup1, was linked to the PHF1 (At3g52190) gene and the other, nla-sup2, was linked to the PHT1.1 (At5g43350) gene.

Plant Transformation

For the complementation and overexpression analysis, the coding and/or genomic sequences of NLA (At1g02860), PHF1 and PHT1.1 genes and a 300 bp region including miR827 sequences were amplified by PCR and cloned into gateway-compatible vectors pEarleyGate 100 and/or pEarleyGate 101 [37]. The plasmids were transformed into Agrobacterium tumefaciens strain EHA105 and then transformed into WT Col, nla-sup1, nla-sup2, nla and/or pho2 plants as described by Clough and Bent [38]. The transgenic plants were selected with the herbicide selection BASTA (active ingredient glufosinate ammonium). Five independent OX lines carrying a single transgene copy for all the constructs were generated. At least two independent lines for each construct were used in the experiments and representative data of one OX line is given in the results.

Biochemical Assays

Frozen shoot and root tissue was used for the following biochemical assays. Pi and total P contents were measured as described by Chiou et al. [18], which is originally adapted from the protocol mentioned by Ames [39]. The percentage of total N inpan> the dried tissues was measured by the Micro-Dumas combustion analysis method usinpan>g a Carlo Erba NA1500 C/N analyzer, (Carlo Erba Strumentazione, Milan, Italy).

Pi Uptake Assay

The plants were grown for 12 day on agar plates (1% sucrose and 0.8% agar) with 1 mM KNO3 and 0.5 mM KH2PO4 with the other nutrient concentrations the same as in the soil experiments. Ten plants were pooled and roots were immersed in 7 ml nutrient solution, except with the replacement of KH2PO4 by [33P]orthophosphate (10.2×10−2 MBq/7 ml; Perkin-Elmer). The plants were incubated for 2 and 4 h and then rinsed with nutrient solution without 33P. The plants were weighed and lysed in 500 µl of 30% H2O2 and 200 µl of perchloric acid for 1 h at 70°C and 5 ml scintillation liquid was added to each sample, incubated over night and the activity measured by scintillation counter.

Quantitative Real-Time PCR

Total RNA was isolated usinpan>g TRIZOL reagent (Invitrogen, Carlsbad, CA, USA). To eliminate any residual genomic DNA, total RNA was treated with RQ1 RNase-free DNase (Promega, Madison, WI, USA). The cDNA was synthesized from total RNA by using the Reverse Transcription System kit (Promega). Primer Express 2.0 software (Applied Biosystems, Forster City, CA, USA) was used to design the primers. Primer sequences for each gene are given in Table S5. Real-time PCR was performed as described previously [40]. Relative quantification (RQ) values for each target gene were calculated by the 2−ΔΔCΤ method [41] using UBIQUITIN10 (UBQ10) and/or GLYCERALDEHYDE 3-PHOSPHATE DEHYDROGENASE A SUBUNIT (GAPA) as an internal reference gene. To ensure the validity of the 2−ΔΔCΤ method, twofold serial dilutions of cDNA were used to create standard curves, and the amplification efficiencies of the target and reference genes shown to be approximately equal.

Statistical Analysis

Statistical analysis was done by Fisher's protected LSD test using SAS statistical software (SAS Institute, Inc., n class="Chemical">NC). The results shown are representative of three inpan>dependent experiments and withinpan> each experiment treatments were replicated three times, unpan>less otherwise stated.

Effect of nitrate and phosphate regimes on growth of Arabidopsis WT, nla-sup2, OX-NLA, nla and pho2 plants. Plants grown at (A, C) 10 mM NO3−-10 mM Pi, and (B, D) 3 mM NO3−-10 mM Pi conditions.

(TIF)

Click here for additional data file.

Effect of nitrate and phosphate regimes on growth of Arabidopsis WT and nla plants. Plants grown at (A) 10 mM NO3−-10 mM Pi, (B) 3 mM NO3−-10 mM Pi, and (C) 10 mM NO3−-3 mM Pi conditions.

(TIF)

Click here for additional data file.

Effect of nitrate and phosphate regimes on growth of Arabidopsis WT and nla (ecotype Wassilewskija, Ws) plants. Plants grown at (A) 10 mM NO3−-10 mM Pi, and (B) 3 mM NO3−-10 mM Pi conditions.

(TIF)

Click here for additional data file.

Effect of ammonium and phosphate regimes on growth of Arabidopsis WT and nla plants. Plants at (A) 20 days after sowing, and (B) 28 days after sowing.

(TIF)

Click here for additional data file.

Effect of varying N (A) and P (B) regimes on growth of Arabidopsis WT plants.

(TIF)

Click here for additional data file.

Pi content in n class="Species">Arabidopsis shoots.

(DOC)

Click here for additional data file.

Pi content in n class="Species">Arabidopsis roots.

(DOC)

Click here for additional data file.

Pi content in n class="Species">Arabidopsis roots and inpan>florescence.

(DOC)

Click here for additional data file.

Total nitrogen content inpan> Arabidopsis WT plants.

(DOC)

Click here for additional data file.

Primers used for real-time PCR analysis of gene expression.

(DOC)

Click here for additional data file.