Original Research ARTICLE
RNA-SEQ Reveals Transcriptional Level Changes of Poplar Roots in Different Forms of Nitrogen Treatments
- 1State Key Laboratory of Tree Genetics and Breeding, School of Forestry, Northeast Forestry University, Harbin, China
- 2College of Life Science, Northeast Forestry University, Harbin, China
- 3School of Forestry, Northeast Forestry University, Harbin, China
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 -N and -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 (), S2 (NH4NO3), and S3 () 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.
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 (-N) and/or nitrate (-N) from soil. For most plants, a mixed nutrition of and is superior over sole -N or -N source (Marschner, 2011). The proportion of to 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 and , 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 as sole N source, both total dry weight and fruit dry weight were increased by 11 and 30% when adding a low concentration of (N:A ratio = 75:25) to the culture solution (Wang et al., 2003). as the sole N source resulted in lower dry weight of tobacco roots as compared with the other N forms ( and NH4NO3; 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 and , 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 in a global regulation of plant transcriptome has been extensively explored. The previous studies show that as compared with N-free samples, supplying 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 and/or NH4NO3 nutrition are less compared with the studies on that of , though several recent publications reported regulation of on gene expression in various plant systems. Fizames et al. (2004) identified 270 genes differentially expressed in Arabidopsis roots when supplied with or NH4NO3 as N source. Zhu et al. (2006) demonstrated that 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 ( or , 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−2 s−1 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 , 0.5 mM NH4NO3, or 1 mM , respectively. The different N forms (S1 [-N]; S2 [NH4NO3]; S3 [-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. 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 [-N]; S2 [NH4NO3]; S3 [-N]) for another 21 days. The total treatment period is 34 days.
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 log2 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 log2 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).
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 and NH4NO3 treated seedlings were higher than that of treated seedlings for 21 days.
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 (), S2 (NH4NO3), and S3 (), 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. 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. 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. Change in gene expression profiles among different N forms. S1, treatment; S2, NH4NO3; S3, . 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. 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 , while the genes induced or repressed by were grouped in Cluster 2 and 3 (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 treatment. In Cluster 2, the genes encoding transporters and transporters (such as AMT2;1, NRT 2.5, NRT3.1) were changed significantly under and NH4NO3 forms. Under 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. 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.
In this study, we aimed to investigate transcriptional-level changes of poplar roots under different N forms (, , and NH4NO3) for a long time (21 days). We found root length and dry weight in or NH4NO3 significantly higher than that in 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 treatment and GO term cell wall under treatment. Podgorska et al. (2013) analyzed physiological responses (up to 8 weeks) of Arabidopsis leaves to or treatment and found that 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 -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 , , or NH4NO3 for 48 h, only three genes were found to be specifically /-induced/repressed (Lopes and Araus, 2008). Hoffmann et al. (2007) studied the genes expression of Arabidopsis seedlings continuously on the medium containing or 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 -induced cluster (Cluster 1), an -repressed cluster (Cluster 2), an -induced cluster (Cluster 3), and a -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 and 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 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 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 uptake (Okamoto et al., 2006). These genes are not transporters, but have been shown necessary for 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 treatment, whereas under or NH4NO3 treatment, NRT3.1 expression level was increased significantly as compared with that with treatment. So we considered that might be an essential for NRT3.1 expression, which would promote N absorption. AMT2.1 genes that are associated with absorption were expressed in low abundance of 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 . 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 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.
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.
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).
The Supplementary Material for this article can be found online at: https://www.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.
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
Lopes, M. S., and Araus, J. L. (2008). Comparative genomic and physiological analysis of nutrient response to: , NH4NO3 and in barley seedlings. Physiol. Plantarum 134, 134–150. doi: 10.1111/j.1399-3054.2008.01114.x
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
Keywords: poplar, nitrogen metabolism, RNA-SEQ, nitrogen forms, long-time nitrogen treatment
Citation: Qu C-P, Xu Z-R, Hu Y-B, Lu Y, Yang C-J, Sun G-Y and Liu G-J (2016) RNA-SEQ Reveals Transcriptional Level Changes of Poplar Roots in Different Forms of Nitrogen Treatments. Front. Plant Sci. 7:51. doi: 10.3389/fpls.2016.00051
Received: 24 September 2015; Accepted: 13 January 2016;
Published: 02 February 2016.
Edited by:Raul Antonio Sperotto, Centro Universitário Univates, Brazil
Reviewed by:Agostino Sorgona, Università degli Studi Mediterranea, Italy
Anna Maria Mastrangelo, CRA-Centro di Ricerca per la Cerealicoltura, Italy
Meng Xu, Nanjing Forestry University, China
Zhibin Luo, Northwest A&F University of China, China
Agnieszka Bagniewska-Zadworna, Adam Mickiewicz University in Poznań, Poland
Copyright © 2016 Qu, Xu, Hu, Lu, Yang, Sun and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Guan-Jun Liu, firstname.lastname@example.org
†Present Address: Guan-Jun Liu, School of Forestry, Northeast Forestry University, Harbin, China