{'fraction_non_alphanumeric': 0.03874532413603761, 'fraction_numerical': 0.028003706372902296, 'mean_word_length': 5.176345909283595, 'pattern_counts': {'":': 0, '<': 7, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Comparative analysis of tissue-specific transcriptomic responses to nitrogen stress in spinach (Spinacia oleracea)

Nitrogen (N) is critical to the growth and productivity of crops. To understand the molecular mechanisms influenced by N stress, we used RNA-Sequencing (RNA-Seq) to analyze differentially expressed genes (DEGs) in root and leaf tissues of spinach. N stress negatively influenced photosynthesis, biomass accumulation, amino acid profiles, and partitioning of N across tissues. RNA-seq analysis revealed that N stress caused most transcriptomic changes in roots, identifying 1,346 DEGs. High-affinity nitrate transporters (NRT2.1, NRT2.5) and glutamine amidotransferase (GAT1) genes were strongly induced in roots in response to N deplete and replete conditions, respectively. GO and KEGG analyses revealed that the functions associated with metabolic pathways and nutrient reservoir activity were enriched due to N stress. Whereas KEGG pathway enrichment analysis indicated the upregulation of DEGs associated with DNA replication, pyrimidine, and purine metabolism in the presence of high N in leaf tissue. A subset of transcription factors comprising bHLH, MYB, WRKY, and AP2/ERF family members was over-represented in both tissues in response to N perturbation. Interesting DEGs associated with N uptake, amino acid metabolism, hormonal pathway, carbon metabolism, along with transcription factors, were highlighted. The results provide valuable information about the underlying molecular processes in response to N stress in spinach and; could serve as a resource for functional analysis of candidate genes/pathways and enhancement of nitrogen use efficiency.


Introduction
Nitrogen (N), being a constituent of amino acids, nucleic acids, cofactors, and secondary metabolites, is central to most biological activities in plants. Effects of N on plant morphology rely on its concentration, which is a determining factor in several metabolic processes, resource allocation, growth, and development of plants [1]. N moves along a complex branching and merging pathway which interacts at numerous sites with the carbon (C) flow, ion, and assimilates at the cell and whole plant level [2]. N is typically transported in the form of nitrate , which leads to groundwater contamination, health, and environmental hazards [7,8]. In commercial spinach production, it is estimated that about 60% of N is lost through leaching [9], most likely due to its shallow root system [10] and short production cycle. Although spinach needs significant amounts of N to sustain rapid growth, it is relatively weak in NO 3 − reduction [11][12][13]. One way to tackle this crisis without compromising productivity is to improve nitrogen use efficiency (NUE) [14]. N perturbation induces series of molecular changes associated with intrinsic processes such as chlorophyll synthesis [15], root architecture [16], and metabolism of sugars and sugar phosphates [17]. In particular, N starvation results in diminished levels of N containing metabolites such as amino acids glutamate (Glu) and glutamine (Gln) [18]. As plants take up NO 3 − and ammonium (NH 4 + ) from the soil and assimilate into amino acids in roots/shoots, improving the amino acid partitioning from source to sink would be an effective strategy to optimize N utilization [19]. To date, quite a few studies have been conducted to investigate genome-wide expression analysis in response to N stress in various plants such as Arabidopsis [20][21][22], rice [23][24][25], maize [26,27], wheat [28,29], and Brassica species [30]. Although this information should help in understanding general molecular processes, each species has its unique and non-overlapping mechanisms to withstand N stress. Besides few studies focusing on mechanisms of N uptake [31, 32] and genotypic variation [6,33,34], not much is known about the transcriptomic changes in the vegetative or non-vegetative tissues of spinach in response to N perturbations.
Here we performed RNA-sequencing (RNA-Seq) of leaf and root tissues of spinach combined with physiological and metabolic analyses using two N treatments: high N (200 ppm) and low N (50 ppm). A total of 1346 and 1136 differentially expressed genes (DEGs) were identified in response to N in root and leaf tissues, respectively. The changes in the expression of the identified genes and the impacted pathways in response to N would help in advancing our understanding of tissue-specific molecular and physiological processes in spinach. The novel genes/pathways identified in this study have provided new targets for genetic manipulation or introgression breeding to improve NUE in spinach. elements were kept constant using N-free Hoagland's nutrient solution ( PO 4 , and 0.1 mg�LPPPPP −1 PPPPP ZnCl 2 R. Samples collected from 6-week old plants were frozen in liquid nitrogen and stored at −80˚C until subsequent analyses. Three independent plants were used for metabolic analysis, total RNA extractions, physiological and biochemical analysis.

Determination of physiological traits, N, and free amino acids
Leaf photosynthetic rate (Pn), transpiration rate (Tr), stomatal conductance (Cs), and intercellular CO 2 concentration were measured with LICOR-6400-XT Photosynthetic system (Lincoln, USA) from the fully expanded leaves of 6-week old plants. The chlorophyll content was measured using a portable chlorophyll meter (SPAD-502, Konica Minolta, Tokyo, Japan). The plant samples were analyzed for total N (TKN), NO 3 − , NH 4 + Rusing Easy Chem Plus analyzer (Chinchilla Scientific, Oak Brook, IL). The amino acid analysis was performed using WatersPPPPP R PPPPP Acquity H-class UPLC system coupled to WaterPPPPP R PPPPP's Xevo TQs MS-MS mass spectrometer with electrospray ionization (ESI) probe following pre-established protocol [35]. Data integration and quantitation were performed using the Water's Tar-getLynx™ software. Statistical analysis was performed using JMP Pro 14 (SAS Institute, Cary, NC, USA).
was built using Bowtie v2.2.3, and paired-end clean reads were aligned to the reference genome using TopHat v2.0.12. Gene expression quantification and DEG analysis: HTSeq v0.6.1 was used to count the reads mapped to each gene. FPKM [36] of each gene was calculated based on the length of the gene and reads count mapped to this gene. Differential expression analysis of high N and low N conditions (three biological replicates per tissue per treatment) was performed using the DESeq R package (1.18.0) [37]. Genes with Pvalue < 0.05 found by DESeq were assigned as differentially expressed.
GO and KEGG enrichment analysis. Gene Ontology (GO) [38] enrichment analysis of differentially expressed genes was implemented by the GOseq R package, in which gene length bias was corrected. GO terms with P-value less than 0.05 were considered significantly enriched by DEGs. KOBAS software was used to test the statistical enrichment of differential expression genes in the Kyoto encyclopedia of genes and genome pathways database (KEGG; http://www.genome.jp/kegg/) [39].
Validation by quantitative real-time PCR. To validate the RNA-seq data, quantitative real-time PCR (RT-qPCR) was conducted to examine the expression pattern of twenty DEGs. Primer Premier 3.0 software was used to design gene-specific primers on the basis of the selected unigene sequences (S1 Table). Total RNA was extracted with the Quick-RNA™ Miniprep Kit (Zymo Research Corporation, Irvine, CA) treated with DNase1 (Zymo Research Corporation, Irvine, CA), and subjected to reverse transcription using iScript RT Supermix (Bio-Rad Laboratories, Inc, Hercules, USA). The quality and quantity of the RNA were analyzed by a Denovix DS-11+ spectrophotometer (Wilmington, Delaware, USA). Gene expression analysis via reverse transcription-qPCR was performed using the BioRad CFX96 qPCR instrument using SsoAdv Univer SYBR GRN Master Kit (Bio-Rad Laboratories, Inc, Hercules, USA). The expression levels of selected DEGs were normalized by comparing with an internal reference gene, 18SrRNA [40]. The relative expression levels (Cq values) for each gene were normalized to that of reference genes by taking an average of three biological replicates. The relative expression levels were calculated using the 2 −ΔΔCt method [41]. All RT-qPCR were repeated in three biological and three technical replications.

Results and discussions
Validation of N stress for expression profiling in spinach N availability affects its acquisition and metabolism in plants. Hence a detailed phenotypic characterization preceded sampling of material for RNA-Seq analysis. A generalized response to N limitation was confirmed by evaluating responses of gas exchange parameters in the leaves of 6-week old spinach plants. N stress significantly reduced the net photosynthetic rate, stomatal conductance, and transpiration rate, implying compromised photosynthetic performance (Table 1). Further, elemental N and NO 3 − was significantly reduced in leaf, petiole, and  (Fig 1). NH 4 + content was also lower in leaf tissues due to N stress. As most N in plants is converted to organic form as free amino acids, we compared the amounts of six highly abundant amino acids; Asparagine (Asn), Aspartate (Asp), Gln, Glu, Serine (Ser) and Alanine (Ala) in four tissue types of spinach. The contents of all the amino acids showed a significant reduction in the roots, while foremost N carrying amino acids (Asp, Asn, Glu, Gln) were significantly decreased in the leaf tissue (Fig 2). Unlike total N and content of NO 3 − , no significant changes were observed in amino acids in petiole due to N stress.

Transcriptome sequencing and assembly of sequencing data
The transcriptomic changes induced by N stress in a spinach leaf and root tissues were analyzed by RNA-Seq. Total twelve independent libraries; six from leaves (HNL1, HNL2, HNL3 for high N and LNL1, LNL2, LNL3 for low N) and six from roots (HNR1, HNR2, HNR3 for high N and LNR1, LNR2, LNR3 for low N) were sequenced using the Illumina HiSeq platform. On average, 45.53 to 43.47 million raw reads were generated for leaf and root tissues in both N treatments ( Table 2). Across all reads, the Q20 and Q30 percentage was more than 97 and 93%, respectively (sequencing error rate was less than 0.03%), and GC content for the libraries was more than 43%. Additional significant characteristics of these libraries are summarized in S2 Table. Among all the libraries, the ratio of total mapped reads and multiple mapped reads in both tissue types was above 85% and 2.5%, respectively. Average 85% reads were uniquely

Analysis of differentially expressed genes (DEGs)
The cluster analysis confirmed that a large number of genes were differentially expressed in leaf and root tissues in both treatments (low N and high N) (S3 Fig). The DEG analysis of roots revealed that 1346 transcripts were significantly (p < 0.05) altered when plants were grown in high N. These included 726 upregulated and 620 downregulated transcripts (Fig 3,  S3 Table, S4 Table). A total of 1136 genes were differentially expressed in leaf tissue, of which 550 genes were upregulated (Fig 3, S5 Table), while 586 were downregulated in high N (Fig 3, S6 Table). Within treatments, the number of DEGs was higher in high N treatment while within tissue types, it was more abundant in root tissue in both N treatments (Fig 4). Between treatments, the DEG analysis revealed that expression of 222 genes was significantly (p < 0.05) affected by N stress, of which 88 genes were upregulated and 134 genes downregulated in the presence of high N availability (S4 Fig). Together, all the DEGs in root and leaf tissues under both treatments are presented as a single Venn diagram ( Fig 5). Total 200 and 181 DEGs were uniquely expressed in root and leaf tissue, respectively.

Functional annotation and GO enrichment analysis of DEGs
Among the DEGs that were significantly upregulated in leaf under Low N, functional annotation showed that 15 biological processes, 9 cell component metabolic pathways, and 6  PLOS ONE molecular functions were significantly over-represented (p-value < 0.05) (Fig 6). Similarly, DEGs showing down-regulation in leaf affected 21 molecular functions, followed by 8 biological processes under low N availability (Fig 6). In leaves, the low N availability mainly upregulated DEGs associated with biological processes, metabolic and organo-nitrogen metabolic processes, while the membrane and its integral components, and protein modification processes were over-represented in DEGs that were down-regulated. In the case of high N treatment, among the up-regulated DEGs, functional annotation showed 15 biological processes, 12 cellular components, and 3 molecular functions were significantly affected (Fig 7). While among the downregulated DEGs, 7 biological processes and cellular components, and 16 molecular functions were over-represented. Similar to the upregulated DEGs associated with biosynthetic processes in leaves under low N, high N treatment clearly represented processes related to the biosynthesis and metabolism of nitrogenous compounds and peptides. Among the downregulated DEGs, functions associated with membrane modifications, response to a stimulus, protein modifications, and transferase activity were significantly over-represented in leaf tissue in high N treatment.
Overall, in root tissue, the GO terms such as membrane and cellular components, biological processes associated with a single organism, response to stimulus and transferase activity were over-represented in the up-regulated DEGs while nutrient reservoir activity related function was represented in the down-regulated DEGs due to N stress (S5

KEGG pathways associated with differentially expressed genes
The enrichment of KEGG pathways was similar in leaf and root tissues in both N treatments (Fig 8), showing over-representation of DEGs associated with metabolic pathways, especially in amino acid metabolism and biosynthesis of secondary metabolites. DEGs related to amino acid metabolic pathways such as Ala, Asp, and Glu synthesis, glutathione metabolism, phenylalanine metabolism were enriched in both N treatments. KEGG pathways of DEGs identified in both N treatments in leaf Vs. root tishsue is shown in S8 Fig. Like tissue types, the analysis showed that the number of DEGs in both N treatments were also associated with metabolic and biosynthesis of secondary metabolites pathways.

Quantitative real-time-PCR validation of DEGs from RNA-Seq
To validate RNA-seq results, we selected twenty DEGs to confirm their expression patterns of the by quantitative real-time PCR. The fold-changes of the selected genes using RT-qPCR were consistent (S9 Fig) with the results obtained with RNA-Seq analysis (R 2 = 0.085, r 2 = 0.93 for roots and R 2 = 0.080, r 2 = 0.91 for leaf), indicating reproducibility and credibility of the RNA-seq data to evaluate the expression of nitrogen-responsive genes in spinach. Histograms were produced by comparing the Log2 fold-changes by transcriptome analysis and RT-qPCR (S10 Fig).

Key genes responding to N perturbations
A complex network of developmental, physiological, biochemical, and molecular responses regulates uptake and partitioning of N. We selected a subset of DEGs that were directly associated with N metabolisms, such as N uptake, assimilation, and transport to understand their role in leaf or root tissues based on N availability.

Nitrate transporters
Among the four families of transporters, Nitrate Transporter 1 (NRT1) [42] and Nitrate Transporter 2 (NRT2) [43] participate in NO 3 − uptake, distribution, or storage. As per the external NO 3 − concentration, plants use two specialized transport systems [44], a constitutive, lowaffinity transport system (LATS) and an inducible, high-affinity transport system (HATS) to maximize N uptake efficiency [45][46][47]. The LATS allows transport during high external NO 3 − while low concentrations activate HATS [44], although it also contributes to fulfilling NO 3 − demand at higher concentrations [48]. We first looked at differentially expressed NO 3 − transporters under HATS induced in roots in response to N stress. We found at least three NRT2.1 (Spo03988, Spo09968, Spo03990) members from HATS were induced in roots due to N stress. In particular, the expression of Spo03988 (NRT2.5) and Spo09968 (NRT2.1) were upregulated 4.3 and 2.6-fold higher in N stressed roots. The affinity of NRT2 transporters for now NO 3 − influx in the root has already been demonstrated for NRT2.1 and NRT2.5 [43,49,50]. NRT2.5 is a part of HATS that enables roots to absorb NO 3 − under limited availability, contributes to phloem loading in shoots, and activates nitrate-inducible genes [49,51]. Both; NRT2.1 and NRT2. 5 in Arabidopsis have been characterized as NO 3 − responsive genes [52]. Additionally, expression of Spo15883 (NRG2, Nitrate regulatory gene) was also induced in roots due to limited N availability. These genes are required for NO 3 − signaling and regulate expression of the nitrate-responsive genes, including NRT2.1 [53]. Under limited N, we also observed the induction of three ammonium transporters (Spo24676, Spo16971, Spo24677), implying alternative strategies deployed by roots to accomplish N requirement. Under high N (S3 Table)

Amino acid metabolism
In plants, inorganic N (NO 3 − and NH 4 + ) is incorporated into Glu and Gln in roots before converting to other amino acids and nitrogenous compounds. Amino acids are later transported to the sink organs by the long-distance transport systems of the plant. We found that the expression of a putative glutamine amidotransferase GAT1 (Spo00216) was almost 5 -fold (Log2) up in roots as well as 4 -fold up in leaf under high N. Root specific GAT1 gene has been shown to be highly responsive to N status in Arabidopsis [58]. It is plausible to assume that under high N conditions, GAT1 facilitates the conversion of Gln to Glu for channeling C to the TCA cycle. We found two Glutamate receptor (GLR) genes (Spo15921 and Spo01837) upregulated in response to high N in roots. GLR genes are implicated in C: N signaling, hypocotyl elongation, root growth, stress responses, and general ion transport [59][60][61][62]. Asn and Gln serve as the major N transport and storage compounds in most non-leguminous plants [63]. The expression of Asparagine synthetases (AS) (Spo05295 and Spo08792) were significantly upregulated in the roots, and leaf under high N. AS activity is regulated by N status [60,64] as it catalyzes the conversion of Asp and Gln to Asn and Glu in an ATP-dependent reaction. Aminotransferases have been suggested to play essential roles in redirecting N resources to different pathways [65]. The expression of the branched-chain amino acid (BCAA) aminotransferase genes BCAT-2 (Spo14673 and Spo13997) regulating conversion of leucine/ valine/isoleucine to Glu were significantly upregulated in roots. Among the most enriched pathway terms, several DEGs associated with degradation of branched-chain amino acids were significantly upregulated in the root (33 out of 48 reference genes annotated to the BCAA degradation pathway) as well as leaf (23 out of 48) tissues. Although BCAA catabolism likely has many functions in plants [66][67][68], its involvement in N redistribution needs further functional characterization of the genes involved in BCAA catabolism. Additionally, we found expression of serine decarboxylase (SDC; Spo13254) significantly upregulated in response to high N in roots. Most recently, SDC like gene showing remarkably high expression in young roots under sufficient N has been characterized in tea plants [69]. We performed a BLASTP homology search of the NCBI database. We found that the deduced protein sequence of the spinach SDC gene shared high similarity with characterized Tea Camellia sinensis (74.6% identity) SDC as well as other plant species such as Arabidopsis (77.26% identity), and rice (74.7% identity) (S11 Fig). Decarboxylation of Ser is the major source of ethanolamine production required for plant growth [70]. Moreover, it has been proposed that Ser could serve as a source of N in non-photosynthetic tissues like roots [71]. Nevertheless, additional functional characterization of this gene will be required to understand its role in N metabolism. Under high N, expression of putative amino acid transporters such as WAT-1 related (UMAMIT34, Usually Multiple Acids Move In and out Transporters; Spo14333), vacuolar amino acid transporters (Spo21518 and Spo21518) and ABC transporter C family member (Spo05865) were significantly upregulated (�2 -fold (Log2) in leaf implying increased N assimilation.

Nucleic acid metabolism
The development and proliferation of living cells require DNA replication, which is responsible for genome duplication in plants [72]. DNA replication is initiated and facilitated by the formation of the replication fork comprising the Minichromosome Maintenance MCM(2-7) protein helicase complex [73]. The MCMs are found to be highly expressed in dividing tissues such as shoot apex and root tips of several plants [74], and their expression is coordinated during plant development, possibly at the level of transcription [75]. Consistent with these reports, nearly half of the DEGs (24 out of 50; P < 0.00029) associated with DNA replication fork assembly were up-regulated in leaf tissue under high N. Since DNA replication is linked to cell cycle progression and DNA repair processes, upregulation of genes associated with MCM complex, RFAs, genes regulating DNA polymerase complexes (α primase, δ, and ε), DNA ligases (Lig 1, Fen1) and helicase (Dna2) were not unexpected (S12 Fig). Besides replication fork, upregulation of several DEGs associated with pyrimidine (34 out 116 reference genes annotated for the pyrimidine pathway; P <0.004; S13 Fig) and purine metabolism (41 out of 158 reference genes annotated for the purine pathway; P <0.006; S14 Fig) also supported increased nucleic acid synthesis and plant growth in response to N availability.

Hormonal pathways
Several studies have demonstrated that NO 3 − can regulate biosynthesis, degradation, transport, and signaling of phytohormones [76,77]. The auxin transport and signaling components are critical for altering root architecture in response to N availability in plants [78,79]. A relationship between the availability of NO 3 − and auxin metabolism has been validated in Arabidopsis, maize, and rice [80][81][82]. Among the DEGs associated with hormonal pathways, 11 genes coding Auxinbinding proteins (ABP) with putative auxin receptor function were up-regulated in roots under high N treatment. Although these proteins share homology with many ABP19 and ABP20 proteins from plants, a detailed functional characterization would be required to understand their significance during excess N availability. The Auxin Binding Protein 1 (ABP1) has been extensively studied as a candidate auxin receptor and has been shown to be a key regulator for auxin-mediated responses [83,84]. Due to N stress, two GA2oxs; Gibberellin 2-beta-dioxygenase (Spo11903 and Spo03697) and a Gibberellin-regulated (Spo07806) genes were up-regulated in roots. GA2oxs regulate plant growth by inactivating endogenous bioactive gibberellins. The gibberellin's impact on leaf-growth and enhanced apical dominance is well known; however, little information is available about root-specific changes in gibberellin in response to N stress. Nonetheless, reduction in endogenous GA concentrations in roots due to the down-regulation of GA biosynthesis at the transcriptional level has been shown to regulate root growth and NO 3 − uptake in cucumber plants [85].

Carbon metabolism
In high N treatment, several genes associated with C metabolism were differentially expressed in leaf tissue. In particular, the expression of Phosphoenolpyruvate carboxykinase (PEPCK; Spo20411) was up-regulated by 4-fold (Log2 scale). PEPCK is involved in multiple C metabolism-related functions such as CO 2 assimilation, relocating the consumed TCA intermediates in the biosynthesis, and assimilation of N [86]. Similarly, expression of Isocitrate lyase (ICL; Spo13898) and Malate synthase (MS; Spo16696) involved in the glyoxylate pathway were upregulated (>2-fold Log2). It has been proposed that cellular lipids can be metabolized via the glyoxylate cycle. Sucrose can be synthesized from the four C products of the glyoxylate cycle, which converts TCA cycle components to sucrose via the process of gluconeogenesis. Although detailed enzyme kinetics and tracer experiments will be required to understand the role of these genes in spinach, unusually coordinated expression of PEPCK with MS and ICL poses a possibility of the presence of gluconeogenic pathway in leaf tissue. In C4 plants, like maize, gluconeogenesis has been shown to regulate the import/mobilization of nitrogenous compounds. Upregulation of Expansins (Spo03956 and Spo11774) and Sugar transporter proteins (MST4; Spo11108 and Spo11108) and Beta-amylase 2 (Spo10422) also supports cell expansion and enhanced influx of photo-assimilates in leaf tissue under high N treatment, respectively.

Differentially expressed transcription factors (TFs) responding to N stress
We used iTAK [87] to perform the transcription factor analysis. A total of 786 transcription factor related genes were identified from spinach and the 56 TF families that expressed genes differentially between leaf Vs. root tissues under low N or high N are shown in  Table). It will be valuable to functionally characterize the identified TFs for enhancing plant performance under limited N availability and thereby providing novel targets for genetic manipulation. In roots, low N up-regulated expression (� 2-fold log2) of three ERFs-TINY (Spo16233), ERF17 (Spo20470), and ERF12 (Spo10826); and two NAC TFs (NAC43; Spo25926 and NAC86; Spo15445). The interaction between ethylene and availability of N affects numerous physiological processes, including roots architecture, leaf, reproductive organ development, as well as the synthesis of amino acids, proteins, and enzymes [88]. N starvation increases the number or affinity of root receptors, allowing roots to respond to lower concentrations of ethylene than those found in unstressed roots. The AP2/ERF-TF family has been involved in signaling processes and the responses to environmental stresses [89], such as micronutrient deficiency [90,91]. Similar to our observations, under N starvation, seven ERF genes in cucumber seedlings were also differentially regulated [92]. Under high N, the expression of Ethylene-responsive

PLOS ONE
transcription factor 13 (ERF13; Spo11940) and two of the NAC29 TFs was induced (� 2-fold Log2). We found that the expression of 1-aminocyclopropane-1-carboxylate oxidase1 (ACO; Spo03782 and Spo03777) genes involved in ethylene synthesis were also upregulated in roots. Induction of expression and enzymatic activities of ethylene biosynthetic genes, ACO, and ACC synthase (ACS) in response to high N has been demonstrated in other plants [93][94][95].

Conclusions
In conclusion, comparative analysis of the leaf and root transcriptomes indicated a coordinated regulation of multiple genes and pathways and provided a catalog of transcriptomic variation between the spinach tissues in response to N. Although transcriptomic responses to N perturbation have been studied in Arabidopsis and few other crops, to the best of our knowledge, this is the first study that provides information about the partitioning of leaf and root transcriptomic responses in spinach. The genomic resources generated from this RNA-Seq experiment can be used to characterize genes associated with NUE in spinach breeding materials and to develop informative markers specifically to select parental lines. RNA-seq analysis revealed that N stress caused most transcriptomic changes in roots, identifying 1,346 DEGs compared to 1136 DEGs in leaf. Significant upregulation of NRT2.1 (Spo03988, Spo09968, Spo03990) in roots validated its role in N acquisition under limited availability. The root-specific genetic manipulation to alter the expression of NRT2.1 could be a way forward to enhance NUE in spinach. Concomitant activation of putative glutamine amidotransferase GAT1 (Spo00216) and Asparagine synthetases (AS) (Spo05295 and Spo08792) along with Glutamate receptor (GLR) (Spo15921 and Spo01837) in response to excess N highlights the functional significance of Gln reduction in N assimilation. Consistent with the effects of higher N availability on plant growth, expression of genes associated with nucleic acid synthesis and DNA replication were significantly upregulated. The RNA-Seq data generated in this study will be a valuable resource for the identification of key genes related to N metabolism and identifying potential targets for genetic manipulation to enhance N uptake and utilization. The functional characterization of the differentially expressed genes during N stress would provide new insights into developing N use efficient varieties.