RNA-SEQ Reveals Transcriptional Level Changes of Poplar Roots in Different Forms of Nitrogen Treatments

Poplar has emerged as a model plant for better understanding cellular and molecular changes accompanying tree growth, development, and response to environment. Long-term application of different forms of nitrogen (such as ${\text{NO}}_3^{-}$-N and ${\text{NH}}_4^{+}$-N) may cause morphological changes of poplar roots; however, the molecular level changes are still not well-known. In this study, we analyzed the expression profiling of poplar roots treated by three forms of nitrogen: S1 (${\text{NH}}_4^{+}$), S2 (NH₄NO₃), and S3 (${\text{NO}}_3^{-}$) by using RNA-SEQ technique. We found 463 genes significantly differentially expressed in roots by different N treatments, of which a total of 112 genes were found to differentially express between S1 and S2, 171 genes between S2 and S3, and 319 genes between S1 and S3. A cluster analysis shows significant difference in many transcription factor families and functional genes family under different N forms. Through an analysis of Mapman metabolic pathway, we found that the significantly differentially expressed genes are associated with fermentation, glycolysis, and tricarboxylic acid cycle (TCA), secondary metabolism, hormone metabolism, and transport processing. Interestingly, we did not find significantly differentially expressed genes in N metabolism pathway, mitochondrial electron transport/ATP synthesis and mineral nutrition. We also found abundant candidate genes (20 transcription factors and 30 functional genes) regulating morphology changes of poplar roots under the three N forms. The results obtained are beneficial to a better understanding of the potential molecular and cellular mechanisms regulating root morphology changes under different N treatments.

Introduction

Nitrogen (N) element is one of macronutrients essential for plant growth, which accounts for 1.5–2% of plant dry matter and ~16% of total plant protein (Frink et al., 1999). Plant roots mainly take up inorganic nitrogen in the form of ammonium (${\text{NH}}_4^{+}$-N) and/or nitrate (${\text{NO}}_3^{-}$-N) from soil. For most plants, a mixed nutrition of ${\text{NO}}_3^{-}$ and ${\text{NH}}_4^{+}$ is superior over sole ${\text{NH}}_4^{+}$-N or ${\text{NO}}_3^{-}$-N source (Marschner, 2011). The proportion of ${\text{NO}}_3^{-}$ to ${\text{NH}}_4^{+}$ for optimal plant growth depends on plant species, developmental stage, environmental conditions, and the total concentrations of supplied N (Jackson and Caldwell, 1993; Luo et al., 2013a; Zhang et al., 2014). Morphological characters of plant roots and shoots are significantly different when it is supplied by a moderate concentration of ${\text{NO}}_3^{-}$ and ${\text{NH}}_4^{+}$, respectively, as sole N source (Schortemeyer et al., 1997; Claussen, 2002; Wang et al., 2003). For example, total dry weight of tomato plants was decreased by 32–86% when it was cultured by ${\text{NH}}_4^{+}$ as sole N source, both total dry weight and fruit dry weight were increased by 11 and 30% when adding a low concentration of ${\text{NO}}_3^{-}$ (N:A ratio = 75:25) to the culture solution (Wang et al., 2003). ${\text{NH}}_4^{+}$ as the sole N source resulted in lower dry weight of tobacco roots as compared with the other N forms (${\text{NO}}_3^{-}$ and NH₄NO₃; Zou et al., 2005). Moreover, there were significant differences in physiological characters, including activity of glutamate dehydrogenase (Acer pseudoplatanus), total amino acid concentration (soybeans), photosynthetic rates (wheat and maize), phosphoenolpyruvate carboxylase (Alfalfa), glutamine synthetase (Pea), type II NAD(P)H dehydrogenase, AOX genes and proline oxidase (Arabidopsis) when the plants were supplied with ${\text{NH}}_4^{+}$ and ${\text{NO}}_3^{-}$, respectively (Goodchild and Givan, 1990; Chaillou et al., 1991; Cramer and Lewis, 1993; Pasqualini et al., 2001; Frechilla et al., 2002; Escobar et al., 2006; Patterson et al., 2010).

Morphological and physiological changes induced by different N forms are closely linked to transcription-level changes. In recent years, the role of ${\text{NO}}_3^{-}$ in a global regulation of plant transcriptome has been extensively explored. The previous studies show that as compared with N-free samples, supplying ${\text{NO}}_3^{-}$ to Arabidopsis seedlings make transcriptional-level changes of the biological processes including transcription and RNA processing, biosynthesis of amino acids and nucleic acids, trehalose metabolism, hormone biosynthesis, and N assimilation as well as PtaNAC1, a transcription factor, which is thought to be associated with root architecture under low nitrogen (LN) conditions (Wang et al., 2000, 2003; Scheible et al., 2004; Bi et al., 2007; Gifford et al., 2008; Wei et al., 2013). Researches on molecular effects of ${\text{NH}}_4^{+}$ and/or NH₄NO₃ nutrition are less compared with the studies on that of ${\text{NO}}_3^{-}$, though several recent publications reported regulation of ${\text{NH}}_4^{+}$ on gene expression in various plant systems. Fizames et al. (2004) identified 270 genes differentially expressed in Arabidopsis roots when supplied with ${\text{NO}}_3^{-}$ or NH₄NO₃ as N source. Zhu et al. (2006) demonstrated that ${\text{NH}}_4^{+}$ as N source stimulated sulfur assimilation in rice leaves. In alfalfa, Ruffel et al. (2008) revealed that over 3000 genes expression was regulated by the status of plant N supply. Poplar has emerged as a model system for understanding molecular mechanisms of woody plants growth, development, and response to environment (Brunner et al., 2004). Some progresses have been achieved in morphological, physiological characteristics of some fast-growing poplar trees (such as P. simonii × P. nigra) and selection of stress-tolerant genes (Wang et al., 2011; Chen et al., 2012; Li et al., 2012; Luo et al., 2013a, 2015; Gan et al., 2015). However, to date, few studies focus on linkages of morphological changes of roots and transcriptional-level characteristics when poplar plants were treated by different N-forms nutrition (${\text{NH}}_4^{+}$ or ${\text{NO}}_3^{-}$, or both). In this study, we hypothesized there were potential coupling changes of N metabolism-related genes and root morphology when the poplar plants were treated by different N forms. So we examined transcriptome profiling of the P. simonii × P. nigra roots using high throughput sequencing technique and analyzed potential effects of long-term different N forms on N metabolism and root morphology-related genes of hydroponic-cultured P. simonii × P. nigra seedlings by a large-scale comparative transcriptomes analysis.

Materials and Methods

Plant Material and Treatments

Poplar seedlings (P. simonii × P. nigra) were germinated on LA media (Hewitt, 1966) with 2% (w/v) sucrose and grown in a growth chamber (light intensity of 200 μmol photons m⁻² s⁻¹ for 16 h per day, day/night temperature of 24/22°C, relative humidity of 50–55%) for 25 days. Then the seedlings were transferred to a 4.0-L LA nutrient medium for 10 days, containing complete nutrient solution within hydroponic boxes under a 16/8 h light/dark regime at 24/22°C and constant (60–65%) relative humidity. The seedlings were transferred to N-free medium for 3 days growth, which the time of N-free treatment was determined according to the result of Balazadeh et al. (2014) and our preliminary test (data not shown), then subsets of seedlings were transferred back to LA complete mediums containing 1 mM ${\text{NH}}_4^{+}$, 0.5 mM NH₄NO₃, or 1 mM ${\text{NO}}_3^{-}$, respectively. The different N forms (S1 [${\text{NH}}_4^{+}$-N]; S2 [NH₄NO₃]; S3 [${\text{NO}}_3^{-}$-N]) were resupplied for another 21 days. To minimize the effect of the altered N content on osmotic potential, the nutrient solution was augmented by a certain amount of sodium chloride to maintain the same cation concentration in the nutrient medium. Samples were taken from the roots with different N⁺ processing for 21 days (Urbanczyk and Fernie, 2005; Wei et al., 2013). Then samples were stored at −80°C for RNA extraction (Figure 1).

FIGURE 1

Figure 1. Schematic representation of the experimental set-up: poplar plants were grown hydroponically in complete LA liquid medium (N⁺ medium) for 10 days, then transferred to nitrogen-free medium (N⁻ medium) for 3 days; subsets of plants were transferred back to complete medium (N resupply, S1 [${\text{NH}}_4^{+}$-N]; S2 [NH₄NO₃]; S3 [${\text{NO}}_3^{-}$-N]) for another 21 days. The total treatment period is 34 days.

Root Measurements

Root length (length of taproots from stem end to root tip) was measured using a ruler (Figure S1). Subsequently, the samples were ground into fine powder in liquid N with a mortar and pestle. Frozen powder of each sample was dried at 60°C for 24 h, then the dry weight was calculated.

Illumina Sequencing and Data Processing

Total RNA was isolated by using a modified CTAB method (Chen et al., 2009; Qu et al., 2010), then sent to Beijing Genomics Institute (Shenzhen) where the libraries were produced and sequenced using the Illumina Genome Analyzer (Solexa). All the samples for Digital Gene Expression were run in two biological replicates, which each replicate is a mixture of root samples from four individual seedlings. Sequence tag preparation was done with Illumina's Digital Gene Expression Tag Profiling Kit according to the manufacturer's protocol (version 2.1B). Firstly, the total RNA samples were treated with DNase I to degrade any possible DNA contamination. Then the mRNA was enriched by using the oligo (dT) magnetic beads (for eukaryotes). After mixed with the fragmentation buffer, the mRNA was fragmented into short fragments (about 200 bp). Then the double strands of cDNA were synthesized by a series of primers, buffer, RNase H and DNA polymerase I. The double strand cDNA was purified with magnetic beads. End reparation and 3′-end single nucleotide A (adenine) addition was then performed. Finally, the fragments were ligated with the adaptors and enriched by PCR amplification, each fragment will generate millions of raw reads. Raw sequences were transformed into clean reads after certain steps of data processing, including removal of the 3′ adaptor sequence, empty reads, and low-quality reads.

All clean reads were mapped to the poplar × trichocarpa Torr. Gray contigs assembly using SOAP2 and only no more than a 2-nucleotide mismatch was allowed (Li et al., 2009). Clean reads mapped to the reference contigs assembly from multiple genes were filtered. The remaining clean reads were designed as unambiguous clean reads. The number of unambiguous clean reads for each gene was calculated and then normalized to RPKM (Reads Per Kb per Million reads), which associated the read number with gene expression levels (Morrissy et al., 2009). Differential gene expression between different nitrogen forms samples was determined by taking the log₂ ratio of RPKM.

Identification of Differentially Expressed Genes and Gene Ontology

The NOIseq was used to identify differentially expressed genes for the samples treated by different N forms. Probability ≥ 0.8 and the absolute value of log₂ Ratio > 1 were used as the threshold to judge the significance of gene expression difference (Tarazona et al., 2011). Cluster analysis of gene expression patterns was performed by Genesis based on the K-means method (Soukas et al., 2000; de Hoon et al., 2004). Gene ontology (GO) analysis was applied to predict gene function and calculate the functional category distribution frequency (Du et al., 2010). Pathway analysis was mainly based on the Mapman (Thimm et al., 2004).

Data Validation by qRT-PCR

The primers used for qRT-PCR validation are listed in Table S1. They were designed on the basis of poplar refseq mRNA sequences using the Primer-BLAST web resource at NCBI (National Center for Biotechnology Information; http://www.ncbi.nlm.nih.gov/BLAST). Quantitative RT-PCR (qRT-PCR) was performed using the ABI7500 Real Time System (Applied Biosystems). Gene expression was analyzed quantitatively using the SYBR Green detection system with melting curve. Amplification conditions were 95°C for 3 min, followed by 40 cycles of: denaturation, 95°C for 15 s; annealing (55–60°C) for 20 s; extension at 72°C for 34 s. Samples for qRT-PCR were run in three biological replicates and two technical replicates. The results were normalized using the Pfaffl method to report relative expression (Pfaffl, 2001). For normalization of gene expression, CYC063 and UBQ7 were used as internal standard (Figure S2).

Statistical Analysis of Root Morphological Parameters

All data were analyzed using SPSS 19.0 software (SPSS, Inc., Chicago, IL, USA). The root length and dry weight of poplar seedlings with different N forms were compared by one-way ANOVA on the basis of Duncan's test at the significance level of 0.05 (P < 0.05).

Results

Morphological Characters of Poplar Roots under Different N Forms

Significant difference in root length and dry weight was found in roots treated by different N forms for 21 days (Figure 2). Root length and dry weight of ${\text{NO}}_3^{-}$ and NH₄NO₃ treated seedlings were higher than that of ${\text{NH}}_4^{+}$ treated seedlings for 21 days.

FIGURE 2

Figure 2. Morphological parameters of poplar roots with different N forms for 21 days. Upper: The picture of the root under different N forms; Lower: Length of roots and the dry weighs of roots. Values are the mean of four replicates ± SE. a, b, and c indicate significant difference based on ANOVA analysis and Duncan's test (P < 0.05).

Illumina Sequencing and Aligning to the Reference Genome

We sequenced three groups of cDNA libraries, S1 (${\text{NH}}_4^{+}$), S2 (NH₄NO₃), and S3 (${\text{NO}}_3^{-}$), and generated 18,232,002 sequence reads, encompassing 1.71 Gb of sequence data (Table 1). Each treatment was represented by ~6 million reads that were sufficient for the quantitative analysis of gene expression. The sequence reads were aligned to the poplar reference genome database using soap2 software (set to allow two base mismatches). Of the total reads, 71.39% were matched either to a unique (62.50%) or multiple (8.88%) genomic locations; the remaining 28.61% were unmatched (Table 1). Only the reads aligning entirely inside exonic regions were matched, the reads from exon-exon junction regions were not matched.

TABLE 1

Table 1. Summary of read numbers based on the RNA-SEQ data from poplar roots under different N forms.

Global Analysis of Gene Expression

A total of 22,414, 25,691, and 26,170 genes, ranging from 100 to ≥2000 bp, were detected in the samples of S1, S2, and S3, respectively. As shown in Table 2, the proportion of sequences with matches to poplar databases was higher among the longer assembled sequences. Specifically, a match efficiency of 32.23% was observed for sequences longer than 2000 bp, whereas the match efficiency decreased to 16.89% for those ranging from 500 to 1000 bp, and to 5.62% for sequences between 100 and 500 bp (Table 2).

TABLE 2

Table 2. Distribution of the gene sequences detected in poplar roots treated by different forms of nitrogen by RNA-SEQ.

Gene Expression Profiles under Different Nitrogen Treatments

To obtain statistical differences in gene expression among different N treated libraries, we compared the RPKM-derived read count using a likelihood ratio test. To minimize false positives and negatives, we assumed that a statistical analysis was reliable when applied to genes with an RPKM value ≥ 2 in both of the two replication libraries. It should be noted that the statistical significances are based on expected sampling distributions. To determine the differentially expressed genes among different N treated libraries, the threshold we used is a two-fold or greater change in expression and Probability ≥ 0.8, and we obtained a set of 602 DEGs (Tables S2–S4). A total of 112 significantly changed genes were detected between the S1 and S2 libraries, with 33 up-regulated genes and 79 down-regulated genes (Figure 3; Table S3). Between the S2 and S3 libraries, a total of 171 DEGs were detected, with 76 up-regulated genes and 95 down-regulated genes (Figure 3; Table S4). This suggests that the differentially expressed genes between S1/S2 is smaller than that between S3/S2. After eliminating duplicate genes, we found 463 genes significantly differentially expressed between the N treatments. All the predicted poplar genes were assigned to different functional categories using Blast2GO (version 2.2.5; http://www.blast2go.org/; Conesa et al., 2005).

FIGURE 3

Figure 3. Change in gene expression profiles among different N forms. S1, ${\text{NH}}_4^{+}$ treatment; S2, NH₄NO₃; S3, ${\text{NO}}_3^{-}$. The number of up-regulated and down-regulated genes between S1 and S2, and between S2 and S3 are summarized. Between the S1 and S2 libraries, there are 33 up-regulated genes and 79 down-regulated genes, while 76 up-regulated genes and 95 down-regulated genes between the S3 and S2 libraries, and 130 up-regulated genes and 189 down-regulated genes between the S1 and S3 libraries.

Functional Analysis of DEGs Based on RNA-SEQ Data

Based on sequence homology, all the significantly changed genes were categorized into 26 functional groups in the three main categories (cellular component, molecular function, and biological process) of the GO classification. Among the groups, nine functional groups were significantly enriched (P ≤ 0.05), including extra cellular region (GO:0005576), oxidoreductase activity (GO:0016491), dioxygenase activity (GO:0051213), oxidoreductase activity (single donors; GO:0016701), oxidoreductase activity (paired donors; GO:0016705), ion binding (GO:0043167), cation binding (GO:0043169), plant-type cell wall modification (GO:0009827), and cell wall modification (GO:0042545; Table 3).

TABLE 3

Table 3. Summary of significantly over-represented functional groups in ammonium-regulated and nitrate-regulated gene sets.

Clustering of DEGs in the Three N-treated Conditions

Based on a similarity of gene expression profiles of two-dimensional hierarchical clustering, we classified 463 differential expression profiles into four expression cluster groups (Clusters 1, 2, 3, and 4; Figure 4). Visual inspection of these expression groups suggested diverse and complex patterns of gene regulation. Clusters 1 and 4 contained the genes induced or repressed by ${\text{NO}}_3^{-}$, while the genes induced or repressed by ${\text{NH}}_4^{+}$ were grouped in Cluster 2 and 3 (Figure 4).

FIGURE 4

Figure 4. Expression patterns clustering and gene annotation. Two-dimensional hierarchical clustering classifies 463 differential gene expression profiles into four expression cluster groups according to the similarity of their expression profiles.

In order to identify differentially expressed genes and transcription factors that are associated with different N forms, we identified 20 putative TF genes through the Mapman Transcription factors Database (pathway). Meanwhile, 30 functional genes associated with N forms are shown in Figure 4; the entire data were listed in Table S5. Figure 4 showed that the expression level of extension of cell wall related genes such as EXPA11 and EXPB3 (Cluster 1) were significantly increased under ${\text{NO}}_3^{-}$ treatment. In Cluster 2, the genes encoding ${\text{NH}}_4^{+}$ transporters and ${\text{NO}}_3^{-}$ transporters (such as AMT2;1, NRT 2.5, NRT3.1) were changed significantly under ${\text{NO}}_3^{-}$ and NH₄NO₃ forms. Under ${\text{NH}}_4^{+}$ treatment, the expression levels of the genes encoding N storage proteins are significantly increased (PLP2). Interestingly, there was no significant difference in genes expression in N metabolism pathway, mitochondrial electron transport/ATP synthesis and mineral nutrition (Figure 5; Figure S3).

FIGURE 5

Figure 5. Schematic representation of primary metabolism of poplar roots at different forms of nitrogen. Gray charts represent significantly differentially expressed genes in the metabolic pathways. White charts represent the genes without significantly different expression in the metabolic pathways. The left and right columns represent the S1/S2 and S3/S2 at the transcriptional level, respectively. Up- and down-facing triangles represent an increase and decrease in transcripts. The digital in columns represent the number of change genes. It can be found that there was no significantly difference gene in nitrogen metabolic pathway.

Discussion

In this study, we aimed to investigate transcriptional-level changes of poplar roots under different N forms (${\text{NH}}_4^{+}$, ${\text{NO}}_3^{-}$, and NH₄NO₃) for a long time (21 days). We found root length and dry weight in ${\text{NO}}_3^{-}$ or NH₄NO₃ significantly higher than that in ${\text{NH}}_4^{+}$ condition. Our results (Figure 2) are in agreement with that found in the rice (Oryza Sativa L.; Schortemeyer et al., 1997; Li et al., 2006). To better understand the changes of root morphology and growth under different N forms, we examined differentially expressed genes of poplar roots treated by different N forms (Tables S2–S4). Based on the result of GO analysis, we only found enrichment of GO term oxidative stress under ${\text{NH}}_4^{+}$ treatment and GO term cell wall under ${\text{NO}}_3^{-}$ treatment. Podgorska et al. (2013) analyzed physiological responses (up to 8 weeks) of Arabidopsis leaves to ${\text{NO}}_3^{-}$ or ${\text{NH}}_4^{+}$ treatment and found that ${\text{NH}}_4^{+}$ nutrition led to increase of leaf NAD(P)H/NAD(P)⁺ ratio, reactive oxygen species content and accumulation of biomolecules oxidized by free radicals of Arabidopsis thaliana. Patterson et al. (2010) believed that the ${\text{NH}}_4^{+}$-induced responses were primarily associated with biotic stress and cellular redox (Patterson et al., 2010). In this paper, our results are consistent with the above results.

We examined genes expression of N-form associated metabolic pathways, but did not find significantly differentially expressed genes in N metabolism pathway. In a previous study on transcript levels of barley plants supplied with ${\text{NO}}_3^{-}$, ${\text{NH}}_4^{+}$, or NH₄NO₃ for 48 h, only three genes were found to be specifically ${\text{NO}}_3^{-}$/${\text{NH}}_4^{+}$-induced/repressed (Lopes and Araus, 2008). Hoffmann et al. (2007) studied the genes expression of Arabidopsis seedlings continuously on the medium containing ${\text{NH}}_4^{+}$ or ${\text{NO}}_3^{-}$ for 15 days, and only found two genes differentially expressed under different N forms. In the study of Beatty et al. (2009), the alanine aminotransferase (AlaAT) gene was transferred into rice plants and over-expressed through a tissue-specific promoter; the authors found the transgenic plants had a strong N use efficiency but less change in the transgenic transcriptome as compared with the controls, with only 0.11 and 0.07% differentially regulated genes in roots and shoots, respectively. According to significant difference of N transport-related genes (Figure 5), we assumed that N transport related genes play an important role in the regulation of long-time N uptake, and N metabolism related genes may be in a steady state in the poplar roots under N forms treatment for 21 days.

In this study, we analyzed genes expression of N-form treated poplar roots by Mapman software (Thimm et al., 2004); and identified four gene expression clusters, a ${\text{NO}}_3^{-}$-induced cluster (Cluster 1), an ${\text{NH}}_4^{+}$-repressed cluster (Cluster 2), an ${\text{NH}}_4^{+}$-induced cluster (Cluster 3), and a ${\text{NO}}_3^{-}$-repressed cluster (Cluster 4), respectively. The significantly differentially expressed genes (Table S5) were classified to 20 metabolism pathways (Figure 5), some of which (30 functional genes) have been clearly annotated in Figure 4. We found gene families members related to the synthesis of cell wall in different clusters (Figure 4; Table S5). For example, exp3 and expA11 belong to cluster 1, and expA17 belongs to cluster 2, while paralog expB2 belongs to cluster 4 (Figure 4). EXPA and EXPB were known to have cell-wall loosening activity and to be involved in cell expansion and other developmental events during which cell-wall modification occurs (Cosgrove, 2000). EXPA and several EXPB are implicated as catalysts of “acid growth,” and regulate the expansion activity rapidly by modulating pH of cell wall (McQueen-Mason et al., 1992; Li et al., 1993, 2003; Cho and Kende, 1997; Sampedro and Cosgrove, 2005). In fine roots of Populus popularis, net fluxes of ${\text{NH}}_4^{+}$ and ${\text{NO}}_3^{-}$ are in association with H⁺ fluxes and change the pH around the root (Luo et al., 2013b). In the present study, the transcriptional level changes of exp genes may be related only to N form but not to change of extracellular pH; because we replaced the culture solution per 2 days for eliminating the effect of medium pH. It is inferred that the exp family genes may play a key role in morphogenesis of the poplar roots when they are treated by different N forms for a long time.

Nitrogen absorption related genes, for example, NRT2.5 and NRT3.1, and the genes associated with ${\text{NH}}_4^{+}$ absorption, AMT2.1, occurred in the cluster 2. The NRT3 family in Arabidopsis contained two members, AtNRT3.1 and AtNRT3.2. The NRT3 family genes in Arabidopsis play a role in ${\text{NO}}_3^{-}$ transport (Okamoto et al., 2006; Orsel et al., 2006). The two NRT3 genes appear to be closely correlated with each other, but NRT3.1 (NAR2.1) appears to play a more significant role in high-affinity ${\text{NO}}_3^{-}$ uptake (Okamoto et al., 2006). These genes are not ${\text{NO}}_3^{-}$ transporters, but have been shown necessary for ${\text{NO}}_3^{-}$ transport through interaction with the other NRT2 transporters (Plett et al., 2010). From Figure 4, we found that the expression of NRT3.1 gene was inhibited only by ${\text{NH}}_4^{+}$ treatment, whereas under ${\text{NO}}_3^{-}$ or NH₄NO₃ treatment, NRT3.1 expression level was increased significantly as compared with that with ${\text{NH}}_4^{+}$ treatment. So we considered that ${\text{NO}}_3^{-}$ might be an essential for NRT3.1 expression, which would promote N absorption. AMT2.1 genes that are associated with ${\text{NH}}_4^{+}$ absorption were expressed in low abundance of ${\text{NH}}_4^{+}$ as a sole nitrogen source.

Poplar BSP genes, which belong to the nucleoside phosphorylases gene family, are expressed before the metabolic nucleotide salvaging, and play an important role in ecophysiological adaptation for inter- and intra-seasonal N storage and cycling (Pettengill et al., 2013). In this study, we obtained a transcription read of bsp-like gene (Comments as PLP; Figure 4), which might be a kind of storage protein of a poplar, specifically induced by ${\text{NH}}_4^{+}$. The magnitude of gene expression in stems is significantly higher than that of leaves and roots, and has a high homology with the NP-like subfamily of the BSP gene family (data not shown). We speculated that the BSP protein synthesis may have a close correlation with the exogenous ${\text{NH}}_4^{+}$ concentration.

Recent publications show that myb48 (Plavcova et al., 2013) and wrky75 (Devaiah et al., 2007) are significantly differentially expressed under different N treatments. In the present study, we identified 20 transcription factors (TF) mRNAs that have close correlations with N forms. However, potential roles of the identified TF mRNAs in root morphogenesis are unknown. So functional verification of the transcription factors is necessary in future work.

In summary, there are significant differences in morphological characteristics and N transport-related genes of the poplar roots between the three N forms for 21 days, but with no significant difference in the transcription level in N-metabolism pathway (Figure 5). Further study should focus on transcriptional level changes of poplar roots treated by different N forms for a short-term period using RNA-SEQ and reveal potential molecular mechanisms. Further work is proposed to select functional genes related to root morphology and functional assignment of the transcription factors under different nitrogen forms.

Author Contributions

Conceived and designed the experiments: CQ, GL. Performed the experiments: CQ, ZX. Analyzed the data: CQ, ZX, CY, YL. Revised the paper: YH, GL, GS. Wrote the paper: CQ, GL. All authors have read and approve of the final manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This project was financially supported by the National High Technology Research and Development Program of China (863 Program, 2013AA102702), National Natural Science Foundation of China (No. 31570648), and by Fundamental Research Funds for the Central Universities (DL13EA03-01).

Supplementary Material

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2016.00051

Table S1. The genes primers used for qRT-PCR validation are listed.

Table S2–S4. To determine which of the 40,361 genes were differentially expressed among the different N-treatments libraries, we required a two-fold or greater change in expression and Probability ≥ 0.8, which resulted in a set of 602 DEGs.

Table S5. The annotation results of expression pattern clustering analysis.

Figure S1. Morphological parameters of poplar roots under different N forms for 21 days. Upper: a whole seedling under different N forms, respectively; Lower: a detailed root structure. Table- dry root weights. Values are the mean of four replicates ± SE. Different letters mean significant difference at a significance of 0.05 through ANOVA analysis and Duncan's test.

Figure S2. Validation of RNA-SEQ results by real-time PCR; r is the correlation coefficient value between the two platforms.

Figure S3. Alidation of nitrogen metabolism-related genes expression by real-time PCR.

References

Balazadeh, S., Schildhauer, J., Ara, W. L., Munn9-Bosch, S., Fernie, A. R., Proost, S., et al. (2014). Reversal of senescence by N resupply to N-starved Arabidopsis thaliana: transcriptomic and metabolomic consequences. J. Exp. Bot. 65, 3975–3992. doi: 10.1093/jxb/eru119

PubMed Abstract | CrossRef Full Text | Google Scholar

Beatty, P. H., Shrawat, A. K., Carroll, R. T., Zhu, T., and Good, A. G. (2009). Transcriptome analysis of nitrogen−efficient rice over−expressing alanine aminotransferase. Plant Biotechnol. J. 7, 562–576. doi: 10.1111/j.1467-7652.2009.00424.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Bi, Y. M., Wang, R. L., Zhu, T., and Rothstein, S. J. (2007). Global transcription profiling reveals differential responses to chronic nitrogen stress and putative nitrogen regulatory components in Arabidopsis. BMC Genomics 8:281. doi: 10.1186/1471-2164-8-281

PubMed Abstract | CrossRef Full Text | Google Scholar

Brunner, A. M., Busov, V. B., and Strauss, S. H. (2004). Poplar genome sequence: functional genomics in an ecologically dominant plant species. Trends Plant Sci. 9, 49–56. doi: 10.1016/j.tplants.2003.11.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Chaillou, S., Vessey, J. K., Morot-Gaudry, J. F., Raper, C. D., Henry, L. T., and Boutin, J. P. (1991). Expression of characteristics of ammonium nutrition as affected by pH of the root medium. J. Exp. Bot. 42, 189–196. doi: 10.1093/jxb/42.2.189

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, C. B., Zhu, S. H., and Zhou, J. (2009). Improved CTAB method for total RNA extraction of mature feicheng Peach fruit. Shandong Agric. Sci. 5, 102–104.

Google Scholar

Chen, S., Jiang, J., Li, H., and Liu, G. (2012). The salt-responsive transcriptome of Populus simonii × Populus nigra via DGE. Gene 504, 203–212. doi: 10.1016/j.gene.2012.05.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Cho, H. T., and Kende, H. (1997). Expansins in deepwater rice internodes. Plant Physiol. 113, 1137–1143. doi: 10.1104/pp.113.4.1137

PubMed Abstract | CrossRef Full Text | Google Scholar

Claussen, W. (2002). Growth, water use efficiency, and proline content of hydroponically grown tomato plants as affected by nitrogen source and nutrient concentration. Plant Soil 247, 199–209. doi: 10.1023/A:1021453432329

CrossRef Full Text | Google Scholar

Conesa, A., Gotz, S., Garcia-Gomez, J. M., Terol, J., Talon, M., and Robles, M. (2005). Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21, 3674–3676. doi: 10.1093/bioinformatics/bti610

PubMed Abstract | CrossRef Full Text | Google Scholar

Cosgrove, D. J. (2000). Loosening of plant cell walls by expansins. Nature 407, 321–326. doi: 10.1038/35030000

PubMed Abstract | CrossRef Full Text | Google Scholar

Cramer, M. D., and Lewis, O. A. M. (1993). The influence of nitrate and ammonium nutrition on the growth of wheat (Triticum aestivum) and maize (Zea mays) plants. Ann. Bot. 72, 359–365. doi: 10.1006/anbo.1993.1119

CrossRef Full Text | Google Scholar

de Hoon, M. J., Imoto, S., Nolan, J., and Miyano, S. (2004). Open source clustering software. Bioinformatics 20, 1453–1454. doi: 10.1093/bioinformatics/bth078

PubMed Abstract | CrossRef Full Text | Google Scholar

Devaiah, B. N., Karthikeyan, A. S., and Raghothama, K. G. (2007). WRKY75 transcription factor is a modulator of phosphate acquisition and root development in Arabidopsis. Plant Physiol. 143, 1789–1801. doi: 10.1104/pp.106.093971

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, Z., Zhou, X., Ling, Y., Zhang, Z., and Su, Z. (2010). agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 38(suppl. 2), W64–W70. doi: 10.1093/nar/gkq310

PubMed Abstract | CrossRef Full Text | Google Scholar

Escobar, M. A., Geisler, D. A., and Rasmusson, A. G. (2006). Reorganization of the alternative pathways of the Arabidopsis respiratory chain by nitrogen supply: opposing effects of ammonium and nitrate. Plant J. 45, 775–788. doi: 10.1111/j.1365-313X.2005.02640.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Fizames, C., Muños, S., Cazettes, C., Nacry, P., Boucherez, J., Gaymard, F., et al. (2004). The Arabidopsis root transcriptome by serial analysis of gene expression. Gene identification using the genome sequence. Plant Physiol. 134, 67–80. doi: 10.1104/pp.103.030536

PubMed Abstract | CrossRef Full Text | Google Scholar

Frechilla, S., Lasa, B., Aleu, M., Juanarena, N., Lamsfus, C., and Aparicio-Tejo, P. M. (2002). Short-term ammonium supply stimulates glutamate dehydrogenase activity and alternative pathway respiration in roots of pea plants. J. Plant Physiol. 159, 811–818. doi: 10.1078/0176-1617-00675

CrossRef Full Text | Google Scholar

Frink, C. R., Waggoner, P. E., and Ausubel, J. H. (1999). Nitrogen fertilizer: retrospect and prospect. Proc. Natl. Acad. Sci. U.S.A. 96, 1175–1180. doi: 10.1073/pnas.96.4.1175

PubMed Abstract | CrossRef Full Text | Google Scholar

Gan, H., Jiao, Y., Jia, J., Wang, X., Li, H., Shi, W., et al. (2015). Phosphorus and nitrogen physiology of two contrasting poplar genotypes when exposed to phosphorus and/or nitrogen starvation. Tree physiol. 36, 22–38. doi: 10.1093/treephys/tpv093

PubMed Abstract | CrossRef Full Text | Google Scholar

Gifford, M. L., Dean, A., Gutierrez, R. A., Coruzzi, G. M., and Birnbaum, K. D. (2008). Cell-specific nitrogen responses mediate developmental plasticity. Proc. Natl. Acad. Sci. U.S.A. 105, 803–808. doi: 10.1073/pnas.0709559105

PubMed Abstract | CrossRef Full Text | Google Scholar

Goodchild, J. A., and Givan, C. V. (1990). Influence of ammonium and extracellular pH on the amino and organic acid contents of suspension culture cells of Acer pseudoplatanus. Physiol. Plant. 78, 29–37. doi: 10.1111/j.1399-3054.1990.tb08710.x

CrossRef Full Text | Google Scholar

Hewitt, E. J. (1966). Sand and Water Culture Methods Used in the Study of Plant Nutrition, Vol. 1. Farnham Royal: Commonwealth Agricultural Bureaux.

Google Scholar

Hoffmann, A., Milde, S., Desel, C., Hümpel, A., Kaiser, H., Hammes, E., et al. (2007). N form–dependent growth retardation of Arabidopsis thaliana seedlings as revealed from physiological and microarray studies. J. Plant Nutr. Soil Sci. 170, 87–97. doi: 10.1002/jpln.200625032

CrossRef Full Text | Google Scholar

Jackson, R. B., and Caldwell, M. M. (1993). Geostatistical patterns of soil heterogeneity around individual perennial plants. J. Ecol. 81, 683–692. doi: 10.2307/2261666

CrossRef Full Text | Google Scholar

Li, B., Xin, W., Sun, S., Shen, Q., and Xu, G. (2006). Physiological and molecular responses of nitrogen-starved rice plants to re-supply of different nitrogen sources. Plant Soil 287, 145–159. doi: 10.1007/s11104-006-9051-1

CrossRef Full Text | Google Scholar

Li, H., Li, M. C., Luo, J., Cao, X., Qu, L., Gai, Y., et al. (2012). N-fertilization has different effects on the growth, carbon and nitrogen physiology, and wood properties of slow- and fast-growing Populus species. J. Exp. Bot. 63, 6173–6185. doi: 10.1093/jxb/ers271

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, L. C., Bedinger, P. A., Volk, C., Jones, A. D., and Cosgrove, D. J. (2003). Purification and characterization of four β-expansins (Zea m 1 isoforms) from maize pollen. Plant Physiol. 132, 2073–2085. doi: 10.1104/pp.103.020024

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, R., Yu, C., Li, Y., Lam, T. W., Yiu, S. M., Kristiansen, K., et al. (2009). SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics 25, 1966–1967. doi: 10.1093/bioinformatics/btp336

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Z. C., Durachko, D. M., and Cosgrove, D. J. (1993). An oat coleoptile wall protein that induces wall extension in vitro and that is antigenically related to a similar protein from cucumber hypocotyls. Plantarum 191, 349–356. doi: 10.1007/BF00195692

CrossRef Full Text | Google Scholar

Lopes, M. S., and Araus, J. L. (2008). Comparative genomic and physiological analysis of nutrient response to: ${\text{NH}}_4^{+}$, NH₄NO₃ and ${\text{NO}}_3^{-}$ in barley seedlings. Physiol. Plantarum 134, 134–150. doi: 10.1111/j.1399-3054.2008.01114.x

CrossRef Full Text

Luo, J., Li, H., Liu, T. X., Polle, A., Peng, C. H., and Luo, Z. B. (2013a). Nitrogen metabolism of two contrasting poplar species during acclimation to limiting nitrogen availability. J. Exp. Bot. 64, 4207–4224. doi: 10.1093/jxb/ert234

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, J., Qin, J., He, F., Li, H., Liu, T., Polle, A., et al. (2013b). Net fluxes of ammonium and nitrate in association with H+ fluxes in fine roots of Populus popularis. Plantarum 237, 919–931. doi: 10.1007/s00425-012-1807-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, J., Zhou, J., Li, H., Shi, W., Polle, A., Lu, M., et al. (2015). Global poplar root and leaf transcriptomes reveal links between growth and stress responses under nitrogen starvation and excess. Tree Physiology. 35, 1283–1302. doi: 10.1093/treephys/tpv091

PubMed Abstract | CrossRef Full Text | Google Scholar

Marschner, P. (ed). (2011). Marschner's Mineral Nutrition of Higher Plants. 3rd Edn. (Amsterdam: Elsevier/Academic Press).

Google Scholar

McQueen-Mason, S., Durachko, D. M., and Cosgrove, D. J. (1992). Two endogenous proteins that induce cell wall extension in plants. Plant Cell 4, 1425–1433. doi: 10.1105/tpc.4.11.1425

PubMed Abstract | CrossRef Full Text | Google Scholar

Morrissy, A. S., Morin, R. D., Delaney, A., Zeng, T., McDonald, H., Jones, S., et al. (2009). Next-generation tag sequencing for cancer gene expression profiling. Genome Res. 19, 1825–1835. doi: 10.1101/gr.094482.109

PubMed Abstract | CrossRef Full Text | Google Scholar

Okamoto, M., Kumar, A., Li, W., Wang, Y., Siddiqi, M. Y., Crawford, N. M., et al. (2006). High-affinity nitrate transport in roots of Arabidopsis depends on expression of the NAR2-like gene AtNRT3. 1. Plant Physiol. 140, 1036–1046. doi: 10.1104/pp.105.074385

PubMed Abstract | CrossRef Full Text | Google Scholar

Orsel, M., Chopin, F., Leleu, O., Smith, S. J., Krapp, A., Daniel-Vedele, F., et al. (2006). Characterization of a two-component high-affinity nitrate uptake system in Arabidopsis. Physiology and protein-protein interaction. Plant Physiol. 142, 1304–1317. doi: 10.1104/pp.106.085209

PubMed Abstract | CrossRef Full Text | Google Scholar

Pasqualini, S., Ederli, L., Piccioni, C., Batini, P., Bellucci, M., Arcioni, S., et al. (2001). Metabolic regulation and gene expression of root phosphoenolpyruvate carboxylase by different nitrogen sources. Plant Cell Environ. 24, 439–447. doi: 10.1046/j.1365-3040.2001.00692.x

CrossRef Full Text | Google Scholar

Patterson, K., Cakmak, T., Cooper, A., Lager, I., Rasmusson, A. G., and Escobar, M. A. (2010). Distinct signalling pathways and transcriptome response signatures differentiate ammonium−and nitrate−supplied plants. Plant Cell Environ 33, 1486–1501. doi: 10.1111/j.1365-3040.2010.02158.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Pettengill, E. A., Pettengill, J. B., and Coleman, G. D. (2013). Elucidating the evolutionary history and expression patterns of nucleoside phosphorylase paralogs (vegetative storage proteins) in Populus and the plant kingdom. BMC Plant Biol. 13:118. doi: 10.1186/1471-2229-13-118

PubMed Abstract | CrossRef Full Text | Google Scholar

Pfaffl, M. W. (2001). A new mathematical model for relative quantification in real-time RT–PCR. Nucleic Acids Res. 29, e45–e45. doi: 10.1093/nar/29.9.e45

PubMed Abstract | CrossRef Full Text | Google Scholar

Plavcova, L., Hacke, U. G., Almeida, M., Li, E., and Douglas, C. J. (2013). Gene expression patterns underlying changes in xylem structure and function in response to increased nitrogen availability in hybrid poplar. Plant Cell Environ. 36, 186–199. doi: 10.1111/j.1365-3040.2012.02566.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Plett, D., Toubia, J., Garnett, T., Tester, M., Kaiser, B. N., and Baumann, U. (2010). ichotomy in the NRT gene families of dicots and grass species. PLoS ONE 5:e15289. doi: 10.1371/journal.pone.0015289

PubMed Abstract | CrossRef Full Text | Google Scholar

Podgorska, A., Gieczewska, K., Łukawska, K., Rasmusson, A. G., Gardeström, P., and Szal, B. (2013). Long−term ammonium nutrition of Arabidopsis increases the extrachloroplastic NAD (P) H/NAD (P)+ ratio and mitochondrial reactive oxygen species level in leaves but does not impair photosynthetic capacity. Plant Cell Environ. 36, 2034–2045. doi: 10.1111/pce.12113

PubMed Abstract | CrossRef Full Text | Google Scholar

Qu, C. P., Xu, Z. R., Liu, G. J., Liu, C., Li, Y., Wei, Z. G., et al. (2010). Differential expression of copper-zinc superoxide dismutase gene of Polygonum sibiricum leaves, stems and underground stems, subjected to high-salt stress. Int. J. Mol. Sci. 11, 5234–5245. doi: 10.3390/ijms11125234

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruffel, S., Freixes, S., Balzergue, S., Tillard, P., Jeudy, C., Martin-Magniette, M. L., et al. (2008). Systemic signaling of the plant nitrogen status triggers specific transcriptome responses depending on the nitrogen source in Medicago truncatula. Plant Physiol. 146, 2020–2035. doi: 10.1104/pp.107.115667

PubMed Abstract | CrossRef Full Text | Google Scholar

Sampedro, J., and Cosgrove, D. J. (2005). The expansin superfamily. Genome Biol. 6:242. doi: 10.1186/gb-2005-6-12-242

PubMed Abstract | CrossRef Full Text | Google Scholar

Scheible, W. R., Morcuende, R., Czechowski, T., Fritz, C., Osuna, D., Palacios-Rojas, N., et al. (2004). Genome-wide reprogramming of primary and secondary metabolism, protein synthesis, cellular growth processes, and the regulatory infrastructure of Arabidopsis in response to nitrogen. Plant Physiol. 136, 2483–2499. doi: 10.1104/pp.104.047019

PubMed Abstract | CrossRef Full Text | Google Scholar

Schortemeyer, M., Stamp, P., and Feil, B. (1997). Ammonium tolerance and carbohydrate status in maize cultivars. Ann. Bot. 79, 25–30. doi: 10.1006/anbo.1996.0298

CrossRef Full Text | Google Scholar

Soukas, A., Cohen, P., Socci, N. D., and Friedman, J. M. (2000). Leptin-specific patterns of gene expression in white adipose tissue. Gene Dev. 14, 963–980. doi: 10.1101/gad.14.8.963

PubMed Abstract | CrossRef Full Text | Google Scholar

Tarazona, S., García-Alcalde, F., Dopazo, J., Ferrer, A., and Conesa, A. (2011). Differential expression in RNA-seq: a matter of depth. Genome Res. 21, 2213–2223. doi: 10.1101/gr.124321.111

PubMed Abstract | CrossRef Full Text | Google Scholar

Thimm, O., Bläsing, O., Gibon, Y., Nagel, A., Meyer, S., Krüger, P., et al. (2004). Mapman: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 37, 914–939. doi: 10.1111/j.1365-313X.2004.02016.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Urbanczyk, W. E., and Fernie, A. R. (2005). Metabolic profiling reveals altered nitrogen nutrient regimes have diverse effects on the metabolism of hydroponically-grown tomato (Solanum lycopersicum) plants. J. Exp. Bot. 56, 309–321. doi: 10.1093/jxb/eri059

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, G., Li, C., and Zhang, F. (2003). Effects of different nitrogen forms and combination with foliar spraying with 6-benzylaminopurine on growth, transpiration, and water and potassium uptake and flow in tobacco. Plant Soil 256, 169–178. doi: 10.1023/A:1026268722353

CrossRef Full Text | Google Scholar

Wang, L., Zhou, B., Wu, L., Guo, B., and Jiang, T. (2011). Differentially expressed genes in Populus simonii× Populus nigra in response to NaCl stress using cDNA-AFLP. Plant Sci. 180, 796–801. doi: 10.1016/j.plantsci.2011.02.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, R., Guegler, K., LaBrie, S. T., and Crawford, N. M. (2000). Genomic analysis of a nutrient response in Arabidopsis reveals diverse expression patterns and novel metabolic and potential regulatory genes induced by nitrate. Plant Cell 12, 1491–1509. doi: 10.1105/tpc.12.8.1491

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, H., Yordanov, Y. S., Georgieva, T., Li, X., and Busov, V. (2013). Nitrogen deprivation promotes Populus root growth through global transcriptome reprogramming and activation of hierarchical genetic networks. New Phytol. 200, 483–497. doi: 10.1111/nph.12375

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, S., Jiang, H., Zhao, H., Korpelainen, H., and Li, C. (2014). Sexually different physiological responses of Populus cathayana to nitrogen and phosphorus deficiencies. Tree Physiol. 34, 343–354. doi: 10.1093/treephys/tpu025

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, G. H., Zhuang, C. X., Wang, Y. Q., Jiang, L. R., and Peng, X. X. (2006). Differential expression of rice genes under different nitrogen forms and their relationship with sulfur metabolism. J. Integr. Plant Biol. 48, 1177–1184. doi: 10.1111/j.1744-7909.2006.00332.x

CrossRef Full Text | Google Scholar

Zou, C., Wang, X., Wang, Z., and Zhang, F. (2005). Potassium and nitrogen distribution pattern and growth of flue-cured tobacco seedlings influenced by nitrogen form and calcium carbonate in hydroponic culture. J. Plant Nutr. 28, 2145–2157. doi: 10.1080/01904160500320624

CrossRef Full Text | Google Scholar