Transcriptome profiling analysis for two Tibetan wild barley genotypes in responses to low nitrogen
BMC Plant Biology volume 16, Article number: 30 (2016)
Nitrogen (N) is the most common limiting factor for crop productivity worldwide. An effective approach to solve N deficiency is to develop low N (LN) tolerant crop cultivars. Tibetan annual wild barley is well-known for its wide genetic diversity and high tolerance to poor soil fertility. Up to date, no study has been done to illustrate the mechanism of LN tolerance underlying the wild barley at transcriptional level.
In this study, we employed Illumina RNA-Sequencing to determine the genotypic difference in transcriptome profile using two Tibetan wild barley genotypes differing in LN tolerance (XZ149, tolerant and XZ56, sensitive). A total of 1469 differentially expressed genes (DEGs) were identified in the two genotypes at 6 h and 48 h after LN treatment. Genetic difference existed in DEGs between XZ149 and XZ56, including transporters, transcription factors (TFs), kinases, antioxidant stress and hormone signaling related genes. Meanwhile, 695 LN tolerance-associated DEGs were mainly mapped to amino acid metabolism, starch and sucrose metabolism and secondary metabolism, and involved in transporter activity, antioxidant activities, and other gene ontology (GO). XZ149 had a higher capability of N absorption and use efficiency under LN stress than XZ56. The higher expression of nitrate transporters and energy-saving assimilation pattern could be attributed to its more N uptake and higher LN tolerance. In addition, auxin (IAA) and ethylene (ETH) response pathways may be also related to the genotypic difference in LN tolerance.
The responses of XZ149 and XZ56 to LN stress differed dramatically at transcriptional level. The identified candidate genes related to LN tolerance may provide new insights into comprehensive understanding of the genotypic difference in N utilization and LN tolerance.
Nitrogen (N) is an essential mineral nutrient element for plant growth and development , which serves as a constituent of many important macro-molecules, including proteins, enzymes, and several plant hormones [2–4]. On the whole, it is a major limiting factor for crop production in the world. Commonly, plants are often exposed to N deficient conditions with the situation becoming more severe due to increasingly-declined soil fertility and widely planted high-yield crop cultivars. Hence, a great amount of N fertilizer is applied to meet N requirement by crops , which not only increases the cost for farmers but also brings the environmental problems. In addition it is well documented that only less than half of the applied N fertilizer is used by crops , with the surplus contributing to severe environmental pollution. Therefore, it is extremely imperative and important to develop crop cultivars with high LN tolerance or N use efficiency (NUE), as it is a basic and also the most efficient approach for coping with low N availability in the soil and insufficient N fertilizer supply.
Actually, plants have evolved many adaptive responses for coping with LN condition. Moreover, it was found that such N limitation adaptability in crops is closely associated with their yield performance [6, 7]. Thus, genetic improvement of LN tolerance in crops would be of significance for developing sustainable agriculture. On the other hand, it has been well documented that NUE is a genetically controlled trait, differing dramatically among genotypes, such as in Arabidopsis , as well as in crops including wheat, rice, maize and barley [9–12]. However, narrower genetic diversity in cultivated barley has become a bottleneck for genetic improvement . On the other hand, the Tibetan annual wild barley, growing in the Tibetan Plateau, which has been proved to be one of the centers of cultivated barley , contains accessions with high tolerance to abiotic stresses, such as drought, salinity and potassium deficiency [15–17]. Meanwhile, the wild barley shows generally better adaption to poor soil fertility, and we also have identified some wild barley genotypes with high LN tolerance in a previous study . Therefore, it may be assumed that the wild barley contains some unique mechanisms of LN stress tolerance.
RNA-Seq, one of next-generation high-throughput sequencing technologies, has been widely used recently, due to low background noise, high sensitivity and reproducibility, great dynamic range of expression and base-pair resolution for transcription profiling . Using this technique, transcriptomic profiles of many plants have been dissected under any given conditions, including biotic and abiotic stresses, such as heat stress  and drought stress . Recently, this method has been also performed on some crops, viz rice, sorghum and cucumber [22–24], to investigate the mechanism of LN stress tolerance.
In our previous research, we found a large genetic variation in LN tolerance among the wild barley accessions . However, up to date, no study has been done to analyze the genetic difference of transcriptomic profiles in response to LN stress. It is imperative for us to reveal the underlying mechanism or to explore the relevant genes of LN tolerance in wild barley. In this study, a comprehensive transcriptome analysis was conducted on two wild barley accessions (XZ149, LN-tolerant and XZ56, LN-sensitive), according to a previous study . The main objectives of the current study were to understand the mechanisms of LN tolerance existed in wild barley; and to determine the signaling pathways and regulatory networks of LN tolerance.
Plant materials and treatments
The experiment was carried out in a greenhouse at Zijingang Campus, Zhejiang University, China. Healthy seeds of two wild barley genotypes (XZ149, LN-tolerant and XZ56, LN-sensitive) were disinfected with 2 % H2O2 for 30 min, rinsed thoroughly with distilled water 5 times, and then soaked for 6 h at room temperature. Then the seeds were moved onto moist filter papers in germination boxes, which were placed in a growth chamber (22/18 °C, day/night) at dark for 3 days, and then incubated for another 7 days with light. At the second leaf stage (10-day old), the similar seedlings were transplanted into black plastic pots (5 L). Endosperm was removed away from the seedlings when they were transplanted, so as to eliminate nutrient supply from seeds. The used full-strength nutrient solution was the same as that reported by Yang et al. . The pH of the solution was adjusted to 5.8 ± 0.1 with NaOH or HCl as required, and was continuously aerated with pumps and renewed every five days. Low N treatment was initiated on the three-leaf-stage seedlings, with 0.2 mM N as LN treatment and 2 mM N as normal level (control).
For biomass and N content determination, the seedlings were harvested and separated into shoots and roots, at 14 d after LN treatment. All the plant samples were heated at 105 °C for 30 min, dried at 80 °C until their weight remained constant, and then dry weight was recorded. N content was determined using Foss Kjeltec 8400.
In order to know the time course of gene HvNRT2.1 expression under LN stress, the roots of XZ149 were sampled with three biological replicates at 3 h, 6 h, 12 h, 24 h, 48 h, 4 d, and 8 d after LN treatment, frozen in liquid nitrogen immediately, and stored at −80 °C for use in RNA extraction.
RNA-Seq sampling, RNA extraction and quality control
For RNA-Seq sampling, the seeds of XZ149 and XZ56 were incubated at the same condition as mentioned above. The samples were taken at 6 h and 48 h after exposure to LN stress (0.2 mM) and control (2 mM), respectively. Roots of four seedlings for each treatment were collected and mixed together at each time point. All samples [totally 16, 2 genotypes (XZ149, LN tolerant and XZ56, LN sensitive) × 2 treatments (LN stress and control) × 2 time points (6 h and 48 h) × 2 biological replications] were prepared for further RNA-Seq analysis.
Total RNA was extracted according to the instructions of miRNeasy mini kit (QIAGEN, Germany). RNA degradation and contamination were monitored on 1 % agarose gels. RNA purity was examined using the NanoPhotometer® spectrophotometer (IMPLEN, CA, USA). RNA concentration was determined using the Qubit® RNA Assay Kit in a Qubit® 2.0 Fluorometer (Life Technologies, CA, USA). RNA integrity was analyzed using the RNA 6000 Pico Assay Kit of the Agilent 2100 Bioanalyzer (Agilent Tecnologies, Santa Clara, CA, USA).
RNA-Seq library construction, sequencing and reads mapping
Sequencing libraries were generated using the Illumina TruSeq™ RNA Sample Preparation Kit (Illumina, San Diego, CA, USA) following the manufacturer`s instructions. Initially, mRNA was purified from the total RNA using poly-T oligo-attached magnetic beads. Then the purified mRNA was fragmented using divalent cations under elevated temperature in Illumina proprietary fragmentation buffer, and reversely transcribed into cDNA using SuperScript II (Invitrogen, Carlsbad, CA, USA). Illumina paired-end sequencing adapters were ligated for preparation of hybridization after making adenylation of the 3’ ends of DNA fragments. In order to preferentially select cDNA fragments of 150 bp, the library fragments were purified using an AMPure XP system (Beckman Coulter, Beverly, USA). With ligated adapters on both ends, DNA fragments were selectively amplified and enriched. Then PCR products were purified again using an AMPure XP system and quantified by Agilent Bioanalyzer 2100 system. Finally, each final cDNA library was applied on one lane of the Illumina paired-end flow cell for the cluster generation processing, and subsequently sequenced on an Illumina NextSeq 500 platform, thus generating 2 × 75 bp pair-ends reads.
The raw reads were generated through the Illumina data processing pipeline (version 1.8). For further analysis, the clean data were obtained by removing low quality bases, empty reads and adaptor sequences at the 3’ end from the raw reads. Meanwhile, the Q20, Q30, GC contents, and sequence duplication level of the clean data were calculated. We downloaded the barley genome sequence and annotation data, and used the ultra high-throughput short read aligner to align RNA-Seq reads to the barley reference genomes on TopHat (http://tophat.cbcb.umd.edu/), and then identified splice junctions between exons by analyzing the mapping results.
Identification of the DEGs and validation of RNA-Seq by quantitative RT-PCR
For gene expression analysis, FPKM (fragments per kilo-base of exon per million fragments mapped reads) was calculated at expression level . The difference in expression between control and treatment (two biological replicates per time point) was analyzed using the DESeq R package (1.10.1) . A FDR (false discovery rate) was set as 0.05 for the threshold of DEGs .
To confirm the validation of the RNA-Seq results, 1 μg total RNA of each sample for RNA-Seq was used for real-time quantitative PCR assays. After eliminating the genomic DNA contamination, first strand cDNA was synthesized with oligo dT primer and Random 6 mers in a 20 μl reaction (Takara, Japan). All quantitative RT-PCR (qRT-PCR) analyses were done in two biological replicates and three technical replicates using a CFX96 system (Bio-Rad, USA). The PCR profiles were as follows: 30 s at 95 °C for pre-denaturation, 40 cycles of 5 s at 95 °C for denaturation and 30 s at 60 °C for annealing, followed by Melt-Curve analysis (60 °C - 95 °C, 0.5 °C increment for 5 s per step) to test the amplicon specificity. For relative quantification, the comparative CT method was used . The amplification of HvGAPDH sequence was used as an endogenous reference to normalize all the data. The gene-specific primers were designed using primer-blast (http:/www.ncbi.nlm.nih.gov/tools/primer-blast/), and all the primers were presented in Additional file 1: Table S1.
Analysis of GO enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG)
GO annotation and KEGG analysis for the DEGs were performed using the Blast2GO program  and the similar steps as reported by Zeng et al. . The GOs distribution associated with DEGs were then obtained from three levels: biological process (BP), molecular function (MF) and cellular component (CC) . The KEGG maps, which contained the EC numbers and enzymatic functions in many metabolic pathways were available in a variety of formats .
Significant differences of physiological traits and gene expression among treatments and genotypes were tested using the Duncan’s Multiple Range Test (DMRT) on data processing system (DPS) statistical software, and the difference at P < 0.05 and P < 0.01 was considered as significant and highly significant, respectively.
Effect of N level on growth performances of two wild barley accessions
XZ149 and XZ56 were considered as LN tolerant and sensitive, respectively, according to a previous study . Although LN treatment caused a significant reduction of shoots dry weight for the two accessions, XZ149 was much less affected than XZ56, with XZ149 and XZ56 showing 7.0 % and 28 % reduction, respectively (Table 1). Meanwhile, LN enhanced root growth increase by 35 % and 28 % for XZ149 and XZ56, respectively (Table 1). The less reduction of shoot dry weight and the much increase in root dry weight both contributed to the higher relative total plant dry weight in XZ149. On the other hand, although there was little difference in shoot N concentration between the two accessions under normal N, XZ149 had significantly higher shoot N concentration than XZ56 under LN (Table 1). As a result, N accumulation of XZ149 was 1.58 times larger than that of XZ56 under LN (Table 1). Obviously the current results proved that XZ149 is more LN tolerant than XZ56.
Identification of DEGs
In order to investigate an appropriate time of sampling for RNA-Seq analysis, we observed the response kinetic of the N starvation responsive gene, HvNRT2.1. Relative expression of HvNRT2.1 was calculated using pair-wise comparison between LN and control at 3 h, 6 h, 12 h, 24 h, 48 h, 4 d, 6 d and 8 d after LN treatment, respectively (Additional file 2: Figure S1). The transcript level of HvNRT2.1 was slightly increased at 3 h and 6 h under LN stress, after then peaked abruptly at 12 h followed by another little peak at 48 h, and thereafter decreased slowly but still remained on a little higher level at 4 d and 8 d (Additional file 2: Figure S1). The results demonstrated that the roots were capable of sensing LN signal and activating relevant signal transduction as early as 3 h after treatment, resulting in differential expression of the relevant genes, which in turn showed highly significant differences at 12 h and 48 h between the N levels. Accordingly, in view of sampling convenience in time arrangement, we then took the samples at 6 h and 48 h for RNA-Seq analysis.
To obtain an overall view of the early LN responsive transcriptome in the two accessions, RNA samples were prepared from the roots of both accessions at 6 h and 48 h after LN treatment. Gene expression profiles of the wild barley roots under both control and LN conditions were analyzed. For each sample, two biological replicates were performed in sequencing. In total, 265597518 clean reads were obtained in the tested samples. For the most samples, 70 % of the sequenced reads could be uniquely mapped (Additional file 3: Table S2).
The transcriptional levels were normalized using the FPKM method. Meanwhile, FDR < 0.05 was used as screening thresholds to test the significance of difference in transcript abundance. Consequently, 1469 DEGs under LN stress were identified using pair-wise comparison of each accession between normal and low N conditions at each time point (Additional file 4: Table S3). These included both up-regulated (782) and down-regulated (728) genes (Fig. 1). Interestingly, DEGs in the tolerant genotype XZ149 (1203) were nearly twice as much as those in the sensitive genotype XZ56 (524) (Additional file 5: Table S4, Additional file 6: Table S5). However, there were 258 DEGs commonly found in both XZ149 and XZ56. The two accessions displayed dissimilar expression patterns, in which the amount of up-regulated DEGs in XZ56 at 48 h decreased to one third of that at 6 h, while XZ149 maintained little change between 6 h and 48 h (Fig. 1).
To confirm the validation of the RNA-Seq data, 15 responsive genes were randomly selected for quantitative RT-PCR analysis. The results from both qRT-PCR and RNA-Seq analysis showed that expressions of these genes were highly consistent, thus validating the RNA-Seq data (Fig. 2).
DEGs involved in nitrogen metabolism
Many genes involved in nitrogen absorption and assimilation were differentially expressed under LN stress relative to the normal condition. In the current study, 12 DEGs encoding nitrate transporters were detected (Table 2). The abundance of these gene transcripts were increased under LN stress. Whereas, the expression patterns of these DEGs in XZ149 differed from those in XZ56, in which transcript levels of most DEGs were increased both in XZ149 and XZ56 at 6 h, and remained being enhanced in XZ149 but declined in XZ56 at 48 h. In addition, three DEGs (MLOC_60308, MLOC_56891, MLOC_51737) were up-regulated only in XZ149, and one DEG (MLOC_65110) was down-regulated only in XZ56. Two LN responsive genes encoding ammonium transporters (MLOC_33834, MLOC_35002) were identified (Table 2), in which, MLOC_35002 was responsive to LN stress only in XZ149. Moreover, seven, one and six DEGs encoding amino acid, lysine histidine and oligopeptide transporters, respectively, were found (Table 2). There were also some DEGs encoding the key enzymes in nitrate assimilation, including two nitrate reductases (NR), one nitrite reductase (NiR) and one glutamate synthase 1 (GOGAT 1) (Table 2). Interestingly, all of these DEGs were down-regulated in XZ149 but remained little changed in XZ56.
Meanwhile, expression of many genes associated with absorption or translocation of other nutrients changed under LN stress, such as phosphate (5), potassium (4), sulfate (2), Zinc (1), iron (9) and molybdenum (1), indicating that uptake of these nutrients in barley plants is affected by N metabolism under cross-talking regulation.
Transcription factors and protein kinases
In total 89 DEGs encoding transcription factors (TFs) were identified in this study, and they belonged to different families, such as Zinc finger (30), bHLH (17), MYB (10), bZIP (10), ERF (7), NAC(7), WRKY (5) and HSF (3) (Fig. 3). XZ149 (75) had more than two fold TFs DEGs than XZ56 (34). Moreover, we found that the proteins with zinc finger domains were the most enriched among the TFs, accounting for 34 % of all these DEGs (Fig. 3).
Kinases play important roles in the development of eukaryotic cells , and adaptation to abiotic stresses. Because some of their targets are transcription factors, they also have the functions of regulating transcription . In the current study, 85 DEGs encoding kinases of different groups were identified, which include serine threonine-protein kinase (STK), leucine-rich repeat receptor-like kinase (LRR), lectin protein kinase (LPK), cysteine-rich receptor-like kinase (CRK), CBL-interacting protein kinase (CIPK), wall-associated receptor kinase and calmodulin-binding receptor-like cytoplasmic kinase (CRCK) (Additional file 7: Table S6).
DEGs related to hormone signaling
In addition to the basic roles in growth and development, phytohormones are also involved in various environmental responses, such as light, salt and drought [34, 35]. It has been proposed that some hormones coordinate demand and acquisition of nitrogen . A total of 47 hormone signaling-related DEGs were found in this study, including gibberellin (GA, 6), cytokinin (CTK, 7), auxin (IAA, 18), ethylene (ETH, 7), abscisic acid (ABA, 8), jasmonic acid (JA, 7), and brassinosteroid (BR, 5) (Additional file 8: Figure S2). Heat map clustering analysis was performed to detect these DEGs involved in hormone signaling. It was found that cytokinin dehydrogenase DEG (MLOC_58639) was up-regulated and remained unchanged in XZ149 and XZ56 under LN stress, respectively (Additional file 8: Figure S2). In addition, three aminocyclopropane-1-carboxylate (ACC) oxidase (ACO) DEGs were up-regulated in XZ149 but remained little changed in XZ56 (Additional file 5: Table S4, Additional file 6: Table S5).
GO analysis and pathway enrichment analysis of LN tolerance related DEGs
Of the total 1469 DEGs, 695 DEGs, showing significant up-regulation in XZ149, but down-regulation or unchanged in XZ56, or little change in XZ149 but down-regulation in XZ56, were selected for further investigation. Based on hierarchical clustering analysis, these DEGs could be mainly grouped into four classes. GO functional enrichment analysis were done to classify these DEGs into their corresponding biological process (BP), molecular function (MF) and cellular component (CC) (Fig. 4). The DEGs with known annotation could be categorized into 37 functional groups. In the biological process ontology, GO terms associated with ‘metabolic process’, ‘single-organism process’, and ‘cellular process’, occupied the majority (Fig. 4). The genes associated with catalytic activity and binding were the most enriched, accounting for 74.3 % of molecular function ontology (Fig. 4). Meanwhile, DEGs related to LN tolerance also acted as diverse cellular components (Fig. 4).
In addition to GO analysis, 695 DEGs were mapped to terms in KEGG pathway enrichment, and the encoded enzymes were assigned to 72 KEGG pathways (Fig. 5), including amino acid, nucleotide, lipid, carbohydrate, energy and other metabolisms. Among these pathways, DEGs involved in phenylpropanoid biosynthesis (34), phenylalanine metabolism (33), and starch and sucrose metabolism (13) were the most abundant. DEGs (MLOC_4686, MLOC_14829) (Additional file 4: Table S3) encoding the two key enzymes phenylalanine ammonia-lyase (PAL) and cinnamate 4-hydroxylase (C4H), respectively, involved in phenylalanine metabolism, as well as DEG (MLOC_6177) encoding flavonoid 3'5'-hydroxylase (F3'5'H) homolog, were up-regulated only in XZ149 but not changed in XZ56 under LN stress. In addition, the DEG (MLOC_68184) encoding chalcone isomerse (CHI) was down-regulated only in XZ56 and had little change in XZ149. Ten enzymes encoded by 13 LN tolerance-related DEGs were found to be associated with starch and sucrose metabolism.
Nitrogen is an essential mineral nutrient required for plant growth and development. However, N deficiency in soil is becoming a big issue in crop production worldwide. It is imperative to develop the crop cultivars with high LN tolerance or NUE for coping with the issue. There is a large difference among species or genotypes within a species in LN tolerance. Relative biomass or dry weight is often used as an indicator of plant tolerance to low nutrition stress [37, 38]. In the present study, the differences in growth performance between the two wild barley genotypes proved our previous result that XZ149 is more tolerant to LN stress than XZ56, and the results also indicated that XZ149 had the higher capability of N absorption and translocation than XZ56. Meanwhile, we used the RNA-Seq to reveal the differences of the transcriptome profiling between the two wild barley genotypes under low N stress. Clearly, there is a distinct difference in the transcriptional level between the two genotypes in their responses to LN.
Nitrogen metabolism genes responsive to LN stress
Nitrate transporters are responsible for nitrate absorption from soils. On the whole, N deficiency enhances the expression of high affinity transport systems for nitrate . It was reported that high affinity nitrate transporters AtNRT2.1, AtNRT2.4 and AtNRT2.5 were induced in N-starved Arabidopsis roots [40, 41]. Similarly, the current study showed that the nitrate transporter DEGs were up-regulated under LN stress (Table 2). Moreover, three nitrate transporter DEGs were up-regulated only in XZ149, and in general, the relative increase in abundance of most nitrate transporter DEG transcripts last longer in XZ149 than in XZ56. This unique and higher expression of nitrate transporters in XZ149 may contribute to higher nitrate uptake efficiency, producing more N-containing metabolites required for its survival under LN stress. Thus, it could be assumed that high LN stress tolerance observed in XZ149 is described to its more N uptake and accumulation in plants.
Nitrogen assimilation is another fundamental biological process in plants, which is very energy consuming . The energy cost is particularly larger when nitrate is used as a major N source. When inorganic N is assimilated in roots, energy and C skeletons are provided through respiration, in which sucrose must be supplied from source leaves . Thus, the location of nitrate assimilation may affect energy budget in plants. Against the backdrop, it is reasonable that N assimilation related DEGs (Table 2) in the roots of the tolerant genotype XZ149 were depressed under LN condition, so as to keep energy supply more efficient than transporting sucrose from source leaves to roots for generating ATP and NAD(P)H used in nitrate assimilation. Alternatively, its nitrate and ammonium may be assimilated in the shoots after their transportation out of the roots via the xylem. This less energy-consuming counter measure may also contribute to XZ149’s tolerance to LN stress.
Carbon metabolism and reducing equivalents
In order to survive under LN stress, some genes related to alleviation of the detrimental effect are abundantly expressed, thus resulting in enhanced stress tolerance. Currently, the transcription of soluble acid invertase (MLOC_60412) was increased in XZ149 but not in XZ56 (Additional file 5: Table S4), indicating that sucrose degradation was enhanced in the tolerant genotype. The similar findings were observed in drought and LN stresses in other studies [23, 44]. Therefore, we may assume that the enhanced invertase expression in the roots of the tolerant genotype XZ149 may stimulate the cycling of sucrose and carbon partitioning in favor of sucrose accumulation for counteracting the LN stress .
6-Phosphogluconate dehydrogenase (6PGDH) is a key enzyme involved in the pentose phosphate pathway, where it catalyzes the conversion of 6-phosphogluconate to ribulose-5-phosphate, with the generation of NADPH. In the current study, the DEGs (MLOC_80338) encoding this enzyme were up-regulated in XZ149 (Additional file 5: Table S4), but unchanged in XZ56. As one of the primary end products of the pentose phosphate pathway, NADPH is necessary for fatty acid synthesis, and is needed in response to oxidative stress. Apart from that, it also helps to maintain the reduced state of glutathione (GSH) by serving as co-substrate for glutathione reductases that reduce oxidized glutathione.
Several TFs have been described in plants exposed to limited N. For example, Peng et al. found that NAC29 showed elevated expression at N deficiency . An R2R3-type MYB TF, CmMYB1, enhances the expression of CmNRT, CmNAR, CmNIR, CmAMT, and CmGS under N-starvation . It was reported that in fungi some GATA factors were involved in regulation of N metabolism and may activate expression of N catabolic enzymes during N-deficiency . However, little has been known about their functions in plants. Here, we found one GATA type zinc finger transcription factor family protein DEG (MLOC_53547) to be up-regulated only in XZ149 (Fig. 3). In addition, three MYB and three NAC transcription factors were found only up-regulated in XZ149 or only down-regulated in XZ56 (Fig. 3). Obviously, it will be quite imperative to determine the possible roles of these TFs in LN stress tolerance in the future.
As a hormone, CTK (cytokinin) acts as a critical signalling molecule in communicating N availability in plants [48, 49]. It has been well documented that CKT is involved in the repression of high-affinity NO3 − transporter genes and root-shoot-root signaling in Arabidopsis [50, 51]. In this study, cytokinin dehydrogenase DEGs (MLOC_58639) involved in CTK degradation pathway was up-regulated and not changed in XZ149 and XZ56 under LN stress, respectively (Additional file 8: Figure S2). Obviously, the regulation of CTK homeostasis in both accessions is consistent with the expression profile of high-affinity NO3 − transporter genes.
Interestingly, we found that three ACO homologs were increased in XZ149 but remained little changed in XZ56 (Additional file 5: Table S4, Additional file 6: Table S5). It is well known that ACO, a rate-limiting enzyme in ethylene synthesis, can oxidize ACC (the direct precursor of ethylene synthesis) to ethylene,and ACO4 and another ACO homolog had responses to N deficiency [45, 52]. Recent investigations on Arabidopsis detected a negative feedback loop between NRT2.1 expression and ethylene biosynthesis under low nitrate level . Thus we deduced that it may be the finely tuning between up-regulation of NRT2.1 expression and ethylene biosynthesis in XZ149 should be beneficial for its LN tolerance.
Recent study visualized reactive oxygen species (ROS) accumulation when Arabidopsis roots were subjected to N deprivation , due to the reduction in the frequency of electron carriers of electron transport systems . Excessive accumulation of ROS in plant cells can result in oxidative stress, a major cause of cellular damage and cell death . To cope with this, plants have developed numerous strategies for the detoxification of ROS. Glutathione S-transferases (GSTs) can directly scavenge peroxides with the help of glutathione as electron acceptor . In this study, three DEGs encoding putative GSTs were up-regulated only in XZ149 (Additional file 9: Table S7). In addition, Cytochrome P450s (CYPs), which are involved in biosynthesis and detoxification of a wide variety of molecules , had higher expression in rice  and sorghum  under LN stress, indicating their specific effect of N-induced CYPs on the antioxidant level. Currently, we found nine CYPs DEGs were up-regulated in XZ149, but unchanged in XZ56, and another three unchanged in XZ149, but down-regulated in XZ56 under LN stress (Additional file 9: Table S7). The expression of one gene (MLOC_75745) for aldehyde dehydrogenase (ALDH) was up-regulated only in XZ149 under LN stress (Additional file 9: Table S7). Moreover, a large amount of peroxidases, up to 29, were up-regulated only in XZ149, apparently contributing to the higher capacity of antioxidant defense in XZ149 (Additional file 9: Table S7). Clearly, XZ149 may develop the higher ability of scavenging excessive ROS through forming a stronger antioxidant system under LN stress, facilitating higher tolerance to LN stress.
Phenylpropanoid metabolism pathway regulated by PAL under LN stress
Phenylalanine is a key amino acid at the interphase of primary and secondary metabolism, and PAL is an initial rate limiting enzyme in phenylpropanoid synthesis . In addition to its important role in plant development, PAL also acts as a key enzyme in response to stress. Its biosynthesis would be stimulated in plants exposed to pathogenic attack, low temperature, salt stress, and low nitrogen, phosphate, or iron . Therefore, PAL has generally been considered as one of the main markers of environmental stress . Possibly, the response to N deficiency could be partly altered through regulation of PAL. In our study, DEG (MLOC_4684) encoding PAL was found only in XZ149 under LN, suggesting that the phenylpropanoid metabolism pathway mediated by PAL regulation may confer to genotypic difference in LN tolerance. Furthermore, lots of DEGs, up to 34, related to LN tolerance, were assigned to phenylpropanoid metabolism pathway according to KEGG metabolic pathway enrichment (Fig. 5). It may be assumed that the enhanced phenylpropanoid metabolism observed in XZ149 under LN stress may contribute to its high tolerance.
Regulatory network of the flavoniod synthesis pathway under LN stress
LN stressed plants shows significant reduction of photosynthetic capacity and become more susceptible to oxidative damage caused by excessive light . As an adaptive strategy, synthesis of photo-protective pigments, such as anthocyanins and flavonols could be enhanced in plants exposed to LN [63, 64], because anthocyanins may act as the ‘sunscreen’ and scavenger of ROS , and flavonols could filter off damaging wavelengths of radiation .
Flavonols and anthocyanins synthesis shared the first two phases. Structural genes in anthocyanins and flavonols synthesis specially displayed higher transcriptional level under LN stress in the tolerant genotype XZ149 (Fig. 6). Starting from the first shared phase phenylpropanoid pathway, the DEGs encoding PAL, C4H and F3’5’H were up-regulated only in XZ149 and remained little changed in XZ56 under LN stress (Fig. 6). Moreover, the DEGs (MLOC_68184) encoding CHI were down-regulated only in XZ56 and had little change in XZ149 (Fig. 6).
In addition to the altered structure genes, the regulatory genes in anthocyanins and flavonols synthesis were also affected by LN stress. A MYB12 orthologous gene (MLOC_69570) was up-regulated in XZ149 under LN stress, while remain unchanged in XZ56. The MYB12 transcription factor, known to be a specific activator of flavonoid metabolism, activated the expression of chalcone synthase (CHS), CHI, flavanone 3-hydroxylase (F3H), and flavonol synthase (FLS) in Arabidopsis , and was also intersection of signaling pathways for auxin- and ethylene-mediated flavonol increases , indicating that the ethylene and auxin signalling may play an important role in low-N-induced flavonoid synthesis in barley roots. Moreover, three ACO homologs were up-regulated in XZ149 but unchanged in XZ56. Besides, two ethylene-responsive transcription factors (ERFs) (MLOC_18164, MLOC_75601) were unchanged/increased in XZ149 while decreased/unchanged in XZ56 (Fig. 3).
In short, the current results strongly suggest that higher flavonoid accumulation in the tolerant genotype XZ149, especially the accumulation of anthocyanins and flavonols, may be attributed to its higher LN tolerance.
Identification of DEG transcripts in plants under LN stress would reveal the genetic mechanism of LN tolerance. Here, the results of RNA-Seq analysis (all the sequences of the DEGs were available in Additional file 10: Table S8) demonstrated that there was the dramatic difference at transcriptional level between the two Tibetan wild barley genotypes in response to LN stress. Accordingly, a hypothetical model was developed for the mechanism of LN tolerance in XZ149 (Fig. 7). We deduced that more N absorption, less energy-consuming N assimilation pattern, more energy-producing model, which is contributed to high LN tolerance, may explain its well growth performance under LN stress. In addition, the current study identified some candidate genes related to LN tolerance, and cast a light on comprehensive understanding of the genotypic difference in N utilization and LN tolerance.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and its additional files.
p-Coumaroyl CoA synthase
Cinnamic acid 4-hydroxylase
Chalcone flavanone isomerase
CBL-interacting protein kinase
Calmodulin-binding receptor-like cytoplasmic kinase
Cysteine-rich receptor-like kinase
Differentially expressed genes
False discovery rate
- GOGAT 1:
Glutamate synthase 1
Kyoto Encyclopedia of Genes and Genomes
Lectin protein kinase
Leucine-rich repeat receptor-like kinase
Reactive oxygen species
Serine threonine-protein kinase
Horchani F, R’bia O, Hajri R, Aschi-Smiti A. Nitrogen nutrition and ammonium toxicity in higher plants. International. 2011;7(1):1–16.
Marschner H, Rimmington G. Mineral nutrition of higher plants. Plant Cell Environ. 1988;11:147–8.
Frink CR, Waggoner PE, Ausubel JH. Nitrogen fertilizer: retrospect and prospect. Proc Natl Acad Sci USA. 1999;96(4):1175–80.
Crawford NM, Forde BG. Molecular and developmental biology of inorganic nitrogen nutrition. Arabidopsis Book. 2002;1:e0011.
Socolow RH. Nitrogen management and the future of food: lessons from the management of energy and carbon. Proc Natl Acad Sci. 1999;96(11):6001–8.
McCullough DE, Girardin PH, Mihajlovic M, Aguilera A, Tollenaar M. Influence of N supply on development and dry matter accumulation of an old and a new maize hybrid. Can J Plant Sci. 1994;74(3):471–7.
Ding L, Wang KJ, Jiang GM, Biswas DK, Xu H, Li LF, et al. Effects of nitrogen deficiency on photosynthetic traits of maize hybrids released in different years. Ann Bot-London. 2005;96(5):925–30.
Chardon F, Noel V, Masclaux-Daubresse C. Exploring NUE in crops and in Arabidopsis ideotypes to improve yield and seed quality. J Exp Bot. 2012;63:3401–12.
Le Gouis J, Beghin D, Heumez E, Pluchard P. Genetic differences for nitrogen uptake and nitrogen utilisation efficiencies in winter wheat. Eur J Agron. 2000;12(3):163–73.
Anbessa Y, Juskiw P, Good A, Nyachiro J, Helm J. Genetic variability in nitrogen use efficiency of Spring Barley. Crop Sci. 2009;49(4):1259–69.
Namai S, Toriyama K, Fukuta Y. Genetic variations in dry matter production and physiological nitrogen use efficiency in rice (Oryza sativa L.) varieties. Breeding Sci. 2009;59(3):269–76.
Presterl T, Seitz G, Landbeck M, Thiemt EM, Schmidt W, Geiger HH. Improving nitrogen-use efficiency in European maize: estimation of quantitative genetic parameters. Crop Sci. 2003;43(4):1259–65.
Ellis RP, Forster BP, Robinson D, Handley LL, Gordon DC, Russell JR, et al. Wild barley: a source of genes for crop improvement in the 21st century? J Exp Bot. 2000;51(342):9–17.
Dai F, Nevo E, Wu DZ, Comadran J, Zhou MX, Qiu L, et al. Tibet is one of the centers of domestication of cultivated barley. Proc Natl Acad Sci USA. 2012;109(42):16969–73.
Wu DZ, Shen QF, Cai SG, Chen ZH, Dai F, Zhang GP. Ionomic responses and correlations between elements and metabolites under salt stress in wild and cultivated barley. Plant Cell Physiol. 2013;54(12):1976–88.
Zhang M, Jin ZQ, Zhao J, Zhang GP, Wu FB. Physiological and biochemical responses to drought stress in cultivated and Tibetan wild barley. Plant Growth Regul. 2015;75(2):567–74.
Zeng JB, He XY, Quan XY, Cai SG, Han Y, Nadira UA, et al. Identification of the proteins associated with low potassium tolerance in cultivated and Tibetan wild barley. J Proteomics. 2015;126:1–11.
Yang LN, Hu HL, Zhu B, Jin XL, Wu FB, Zhang GP. Genotypic variations of N use efficiency in Tibetan wild and cultivated barleys. J Zhejiang Univ (Agric & Life Sci). 2014;40(2):155–64.
Marioni JC, Mason CE, Mane SM, Stephens M, Gilad Y. RNA-Seq: An assessment of technical reproducibility and comparison with gene expression arrays. Genome Res. 2008;18(9):1509–17.
Gleason LU, Burton RS. RNA-seq reveals regional differences in transcriptome response to heat stress in the marine snail Chlorostoma funebralis. Mol Ecol. 2015;24(3):610–27.
Hübner S, Korol AB, Schmid KJ. RNA-Seq analysis identifies genes associated with differential reproductive success under drought-stress in accessions of wild barley Hordeum spontaneum. BMC Plant Biol. 2015;15(1):134.
Yang SY, Hao DL, Song ZZ, Yang GZ, Wang L, Su YH. RNA-Seq analysis of differentially expressed genes in rice under varied nitrogen supplies. Gene. 2015;555(2):305–17.
Gelli M, Duo Y, Konda AR, Zhang C, Holding D, Dweikat I. Identification of differentially expressed genes between sorghum genotypes with contrasting nitrogen stress tolerance by genome-wide transcriptional profiling. BMC Genomics. 2014;15(1):179.
Zhao WC, Yang XY, Yu HJ, Jiang WJ, Sun N, Liu XR, et al. RNA-Seq-based transcriptome profiling of early nitrogen deficiency response in cucumber seedlings provides new insight into the putative nitrogen regulatory network. Plant and Cell Physiol. 2015;56(3):455–67.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.
Robinson MD, Oshlack A. A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. 2010;11(3):R25.
Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, Van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat biotechnol. 2010;28(5):511–5.
Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative CT method. Nat Protoc. 2008;3(6):1101–8.
Conesa A, Go¨tz S, Garcı´a-Go´mez JM, Terol J, Talo´n M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6.
Zeng JB, He XY, Wu DZ, Zhu B, Cai SG, Nadira UA, et al. Comparative transcriptome profiling of two Tibetan wild barley genotypes in responses to low potassium. PLoS ONE. 2014;9(6):e100567.
Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Hanks SK, Quinn AM, Hunter T. The protein kinase family: conserved features and deduced phylogeny of the catalytic domains. Science. 1988;241(4861):42–52.
Hunter T, Karin M. The regulation of transcription by phosphorylation. Cell. 1992;70(3):375–87.
Zhu JK. Salt and drought stress signal transduction in plants. Annu Rev Plant Biol. 2002;53:247–73.
Halliday KJ, Martinez-Garcia JF, Josse EM. Integration of light and auxin signaling. CSH Perspect Biol. 2009;1(6):a001586.
Argueso CT, Ferreira FJ, Kieber JJ. Environmental perception avenues: the interaction of cytokinin and environmental response pathways. Plant Cell Environ. 2009;32(9):1147–60.
Broadley MR, Bowen HC, Cotterill HL, Hammond JP, Meacham MC, Mead A, et al. Phylogenetic variation in the shoot mineral concentration of angiosperms. J Exp Bot. 2004;55(396):321–36.
Hermans C, Hammond JP, White PJ, Verbruggen N. How do plants respond to nutrient shortage by biomass allocation? Trends Plant Sci. 2006;11(12):610–7.
Crawford N, Glass A. Molecular and physiological aspects of nitrate uptake in plants. Trends Plant Sci. 1998;3(10):389–95.
Kiba T, Feria-Bourrellier AB, Lafouge F, Brehaut V, Miller A, Daniel-Vedele F, et al. The Arabidopsis nitrate transporter NRT2.4 plays a double role in roots and shoots of nitrogen-starved plants. Plant Cell. 2012;24(1):245–58.
Lezhneva L, Kiba T, Feria-Bourrellier AB, Lafouge F, Boutet-Mercey S, Zoufan P, et al. The Arabidopsis nitrate transporter NRT2.5 plays a role in nitrate acquisition and remobilization in nitrogen-starved plants. Plant J. 2014;80(2):230–41.
Werf A, Kooijman A, Welschen R, Lambers H. Respiratory energy costs for the maintenance of biomass, for growth and for iron uptake in roots of Carex diandra and Carex acutiformis. Physiol Plant. 1988;72(3):483–91.
Nunes-Nesi A, Fernie AR, Stitt M. Metabolic and signaling aspects underpinning the regulation of plant carbon nitrogen interactions. Mol Plant. 2010;3(6):973–96.
Saeedipour S, Moradi F. Comparison of the drought stress responses of tolerant and sensitive wheat cultivars during grain filling: impact of invertase activity on carbon metabolism during kernel development. J Agric Sci. 2011;3(2):32–44.
Peng MS, Bi YM, Zhu T, Rothstein SJ. Genome-wide analysis of Arabidopsis responsive transcriptome to nitrogen limitation and its regulation by the ubiquitin ligase gene NLA. Plant Mol Biol. 2007;65(6):775–97.
Imamura S, Kanesaki Y, Ohnuma M, Inouye T, Sekine Y, Fujiwara T, et al. R2R3-type MYB transcription factor, CmMYB1, is a central nitrogen assimilation regulator in Cyanidioschyzon merolae. Proc Natl Acad Sci. 2009;106(30):12548–53.
Lowry JA, Atchley WR. Molecular evolution of the GATA family of transcription factors: conservation within the DNA-binding domain. J Mol Evol. 2000;50(2):103–15.
Sakakibara H. Nitrate-specific and cytokinin-mediated nitrogen signaling pathways in plants. J Plant Res. 2003;116(3):253–7.
Sakakibara H, Takei K, Hirose N. Interactions between nitrogen and cytokinin in the regulation of metabolism and development. Trends Plant Sci. 2006;11(9):440–8.
Ruffel S, Krouk G, Ristova D, Shasha D, Birnbaum KD, Coruzzi GM. Nitrogen economics of root foraging: Transitive closure of the nitrate-cytokinin relay and distinct systemic signaling for N supply vs. demand. Proc Natl Acad Sci USA. 2011;108(45):18524–9.
Nacry P, Bouguyon E, Gojon A. Nitrogen acquisition by roots: Physiological and developmental mechanisms ensuring plant adaptation to a fluctuating resource. Plant Soil. 2013;370(1–2):1–29.
Bi YM, Wang RL, Zhu T, Rothstein SJ. Global transcription profiling reveals differential responses to chronic nitrogen stress and putative nitrogen regulatory components in Arabidopsis. BMC Genomics. 2007;8(1):281.
Zheng D, Han X, An Y, Guo H, Xia X, Yin W. The nitrate transporter NRT2.1 functions in the ethylene response to nitrate deficiency in Arabidopsis. Plant Cell Environ. 2013;36(7):1328–37.
Shin R, Berg RH, Schachtman DP. Reactive oxygen species and root hairs in Arabidopsis root response to nitrogen, phosphorus and potassium deficiency. Plant Cell Physiol. 2005;46(8):1350–7.
Grossman A, Takahashi H. Macronutrient utilisation by photosynthetic eukaryotes and the fabric of interactions. Annu Rev Plant Biol. 2001;52(1):163–210.
Mittler R. Oxidative stress, antioxidants and stress tolerance. Trends Plant Sci. 2002;7(9):405–10.
Becana M, Matamoros MA, Udvardi M, Dalton DA. Recent insights into antioxidant defenses of legume root nodules. New Phytol. 2010;188(4):960–76.
Mizutani M, Ohta D. Diversification of P450 genes during land plant evolution. Annu Rev Plant Biol. 2010;61:291–315.
Cai H, Lu Y, Xie W, Zhu T, Lian X. Transcriptome response to nitrogen starvation in rice. J Biosci. 2012;37(4):731–47.
Dixon RA, Paiva NL. Stress-induced phenylpropanoid metabolism. Plant Cell. 1995;7(7):1085.
Sanchez-Ballesta MT, Lafuente MT, Zacarias L, Zacarias L, Granell A. Involvement of phenylalanine ammonia-lyase in the response of Fortune mandarin fruits to cold temperature. Physiol Plant. 2000;108(4):382–9.
Diaz C, Saliba-Colombani V, Loudet O, Belluomo P, Moreau L, Daniel-Vedele F, et al. Leaf yellowing and anthocyanin accumulation are two genetically independent strategies in response to nitrogen limitation in Arabidopsis thaliana. Plant Cell Physiol. 2006;47(1):74–83.
Peng M, Hudson D, Schofield A, Tsao R, Yang R, Gu H, et al. Adaptation of Arabidopsis to nitrogen limitation involves induction of anthocyanin synthesis which is controlled by the NLA gene. J Exp Bot. 2008;59(11):2933–44.
Stewart AJ, Chapman W, Jenkins GI, Graham I, Martin T, Crozier A. The effect of nitrogen and phosphorus deficiency on flavonol accumulation in plant tissues. Plant Cell Environ. 2001;24(11):1189–97.
Stewart AJ, Bozonnet S, Mullen W, Jenkins GI, Lean MEJ, Crozier A. Occurrence of flavonols in tomatoes and tomato-based products. J AGR Food CHEM. 2000;48(7):2663–9.
Mehrtens F, Kranz H, Bednarek P, Weisshaar B. The Arabidopsis transcription factor MYB12 is a flavonol-specific regulator of phenylpropanoid biosynthesis. Plant Physiol. 2005;138(2):1083–96.
Lewis DR, Ramirez MV, Miller ND, Vallabhaneni P, Ray WK, Helm RF, et al. Auxin and ethylene induce flavonol accumulation through distinct transcriptional networks. Plant Physiol. 2011;156(1):144–64.
We thank Prof. Dongfa Sun (Huazhong Agricultural University, China) for providing Tibetan wild barley accessions. This work was supported by Natural Science Foundation of China (31330055 and 31271656), China Agriculture Research System (CARS-05) and Jiangsu Collaborative Innovation Center for Modern Crop Production (JCIC-MCP).
The authors declare that they have no competing interests.
XQ, JZ and GZ designed research. XQ, JZ, GC, ZH and JS performed research. XQ, JZ, LY and GZ analyzed data. XQ and GZ wrote the paper. All authors have read, edited and approved the current version of the manuscript.
The primers used in quantitative real-time PCR. (DOC 51 kb)
Real-time PCR analysis of the HvHRT2.1 gene in XZ149 under low N stress. (DOCX 16 kb)
Summary of mapping reads of the RNA-Seq. (XLSX 16 kb)
The FPKM value of 1469 DEGs in XZ149 and XZ56. (XLSX 464 kb)
DEGs at 6 h and 48 h after low N treatment in XZ149. (XLSX 160 kb)
DEGs at 6 h and 48 h after low N treatment in XZ56. (XLSX 70 kb)
DEGs encoding kinases under low N stress in XZ149 and XZ56. (XLSX 14 kb)
Heat Map analysis of DEGs involved in hormone signaling in XZ149 and XZ56. (DOCX 76 kb)
DEGs related antioxidant stress under low N stress. (DOCX 16 kb)
Gene accession numbers and sequences of 1469 DEGs. (XLS 2369 kb)
About this article
Cite this article
Quan, X., Zeng, J., Ye, L. et al. Transcriptome profiling analysis for two Tibetan wild barley genotypes in responses to low nitrogen. BMC Plant Biol 16, 30 (2016). https://doi.org/10.1186/s12870-016-0721-8
- Low N tolerance
- Differentially expressed genes