1. Introduction
Nitrogen (N) is a crucial component for crop growth, biomass accumulation and yield formation [
1], because it is a constituent of amino acids, which are the elemental units of protein and enzymes. However, nitrogen fertilizer application increases crop yield, often at the rates far greater than the rates at which crops can consume, resulting in a leftover of N in soils [
2]. Improper application of nitrogen fertilizers not only reduces crop nitrogen utilization efficiency, but also has the negative impact on the environment [
3]. Therefore, precise management of crop nitrogen status are important to improve crop yield and reduce environmental pollution [
4].
Nitrogen is mainly available to plants in the form of inorganic nitrogen such as nitrate (NO
3−) and ammonium (NH
4+), but organic sources such urea, free amino acids, and short peptides are also acquired by plants. These nitrogen sources acquisitions by plant depends mainly on pant habitat and time, soil heterogeneity, microbial community, agronomic practices and other environmental conditions [
5]. In aerated soils, NO
3− is the most dominant nitrogen source with concentrations range between 1 to 5 mM, while NH
4+ concentration usually ranges between 20 to 200 μM [
6]. On the contrary free amino acids and urea are usually present in small concentrations, ranging from 1 to 150 µM and <70 µM, respectively [
6,
7,
8]. Urea is the major nitrogen (N) form supplied as fertilizer in agriculture, instead, it is also important N metabolite in plants [
9]. The low manufacturing cost and the high nitrogen content of urea (46% N) make this organic compound a common N fertilizer used in arable soils [
10]. Currently in terrestrial N cycle, organic N uptake by plants is considered to be a critical component [
11]. Before, the uptake of organic N was considered as secondary source of N in low inorganic N soils, such as arctic or boreal ecosystems [
12,
13,
14] but conversely, recent studies on agricultural [
15,
16,
17] and sub-tropical plants [
18] confirmed that organic nitrogen uptake accounts for a significant portion of the N budget in plants.
Different N forms can affect the physiological and metabolic processes of plants, such as nutrient uptake, enzyme activity, photosynthesis rate, respiration rate, water balance, and signaling pathways, thus eventually influencing plant growth and crop yield [
19,
20,
21,
22]. NH
4+ as a sole N source provokes negative effects on NH
4+−sensitive plants. The effects of NH
4+ and NO
3− nutrition on plant growth have been intensively studied; however, the results are not consistent and depend mainly on plant species. Wheat and maize [
23], sugar beet [
24], beans [
25]
,, tobacco [
26] and canola [
27], grow preferentially on NO
3− nutrition; whereas, rice [
28], pine and larch [
29], grow preferentially on NH
4+ nutrition. Urea fertilizer containing 46% of increased N utilization as compared to other fertilizer sources as well as better than other reduced N sources and appeared to be the best source which gave the highest dry matter yield compared to other nitrogen sources [
30,
31]. Several studies on Arabidopsis [
32], pakchoi [
33] and tomato [
34] revealed that utilization of glycine could result in high biomass production compared with plants grown under deficiency of N supply. Wang et al. [
35] once found that glycine as a nitrogen source influences a series of physiological and biochemical reactions, including plant defense systems, redox homeostasis, amino acids metabolism and protein synthesis.
To understand the adaptive mechanisms of crops to nutrient deficiency, a large number of nutrient-responsive genes or proteins have been identified and characterized through high-throughput omics techniques, such as genomics, transcriptomics, and proteomics [
36,
37]. Due to the changes in gene transcripts, protein levels, and enzyme activities, metabolites, are regarded as the direct readouts of plant growth or developmental status [
38]. Metabolome analysis has been performed to categorize nitrogen deficiency responsive metabolites and their respective pathways, in many crops, such as rice, maize, wheat, barley, soybean, tomato, and rapeseed [
39,
40,
41,
42,
43,
44,
45,
46]. Previously, a metabolomics study performed with
Arabidopsis leaf tissues, found alteration in plant phenotype, and biomass, along with regulation of enzymes and major metabolites under moderate N starvation [
47]. Wang et al. [
35] found varying physiological responses of the two cultivars to glycine as a nitrogen source and relate it with the up-regulation or down-regulation of glycine specific responsive proteins.
Wild barley (
Hordeum spontaneum) is the progenitor of cultivated barley (
Hordeum vulgare) and provides a rich source of genetic variations for barley improvement [
48]. In our previous findings, wild genotype XZ16 had shown much better growth at urea, and was less affected under glycine-nitrogen treatment as compared to Hua30 in comparison with the plants exposed to NO
3− fertilizer. The objective of this study was to find out the shoot metabolites responding to four N-form nutritions, to reveal the specific metabolites differing among N-form nutritions and the two barley genotypes.
2. Materials and Methods
2.1. Germination and Growth Conditions
Seeds of two barley genotypes, wild (XZ16, low-N tolerant) and cultivated (Hua30, low-N sensitive) were surface sterilized with 2% H2O2 for 20 min and then thoroughly washed with distilled water five times. Seeds were soaked for five hours and put in sterilized sand for germination. After germination, light was supplied with fluorescent lamps at 250 μmol m−2s−1 until two-leaf stage. Uniform seedlings were then transferred to 1 L black pots containing modified Hoagland solutions with KNO3 as a sole nitrogen source, and aerated with air pumps. The solution was renewed every day and the pH of the solution was adjusted to 5.8±1. Seedlings were grown in a controlled growth room at temperature 22℃ of 14 h day/18℃ of 10 h night until three-leaf stage, and then treatments were performed.
2.2. Treatment and Sampling
Barley seedlings were treated with four nitrogen forms (NO
3-, NH
4+, urea and glycine) at N level of 2 mM based on Hoagland solution. NO
3- and NH
4+ were supplied as KNO
3 and (NH
4)
2SO
4, respectively, and K concentration in other three N treatments without addition of KNO
3 was adjusted with KCl. Basic nutrient solution for barley was modified to make the nitrogen level same for all the treatments [
49]. Shoots were sampled in five replicates after 7 days treatments. Samples were frozen in liquid nitrogen immediately, and stored in a refrigerator at −80℃ for use in metabolomics analysis..
2.3. Elemental Analysis
For element analysis, shoots were collected after 7 days of treatments. Samples were washed with dH2O and blotted on filter paper. To attain a constant weight, samples were dried in an oven at 80°C for 72 hours. Approximately 0.1 g dry samples were used for element analysis. Tissue digestion was carried out in 5 ml HNO3 at 120°C using a microwave digester (MCA 3000, Anton Paar, Beijing, China) and final volume was made 15 ml by adding dH2O to the tube. Finally inductively coupled plasma optical emission spectrometer (ICP-OES, Optima 8000DV, PerkinElmer, USA) was used for mineral concentration determination.
2.4. Metabolites Extraction
Frozen shoots, weighing 10±1 mg were placed into a 2 mL test tube. Then 500 μL of pre-cold extraction mixture (methanol/dH2O (v:v) = 3:1) along with 10 μL internal standard (adonitol, 0.5 mg/mL stock) were added. The samples were vortexed for 30 s and homogenized with ball mill for 4 min at 45 Hz, followed by ultra-sonication for 5 min in ice water. After centrifugation at 4 ℃ for 15 min at 12000 rpm, 200 μL supernatant was added to a fresh test tube. To prepare the QC (Quality control) samples, 60 μL of each sample was taken out and combined together. After evaporation in a vacuum concentrator, 80 μL of methoxyamination hydrochloride (20 mg/mL in pyridine) was added and incubated at 80 ℃ for 30 min, and then derivatized by 100 μL of BSTFA regent (1% TMCS, v/v) at 70 ℃ for 1.5 h. Samples were gradually cooled to room temperature and 5 μL of FAMEs (in chloroform) was added. All samples were then analyzed by a gas chromatograph system (Agilent, Santa Clara, CA) coupled with a time-of-flight mass spectrometer (GC-TOF-MS) (Leco, St. Joseph, MI).
2.5. GC-TOF-MS Analysis
GC-TOF-MS analysis was performed using an Agilent 7890 gas chromatograph coupled with a time-of-flight mass spectrometer. The system utilized a DB-5MS capillary column. Aliquot of 1 μL extraction was injected in a splitless mode. Helium was used as the carrier gas. Front inlet purge flow was 3 mL min−1, and the gas flow rate through the column was 1 mL min−1. The initial temperature was kept at 50 °C for 1 min, and then raised to 310 °C at a rate of 10 °C min−1, finally kept for 8 min at 310 °C. The injection, transfer line, and ion source temperatures were 280, 280 and 250 °C, respectively. The energy was -70 eV in electron impact mode. The mass spectrometry data were acquired in full-scan mode with the m/z range of 50-500 at a rate of 12.5 spectra per second after a solvent delay of 6.30 min.
2.6. Data Preprocessing and Annotation
Chroma TOF (V 4.3x, LECO) software was used for raw data analysis, including peak extraction, baseline adjustment, deconvolution, alignment and integration. LECO-Fiehn Rtx5 database was used for metabolite identification by matching the mass spectrum and retention index. Lastly, the peaks detected in less than half of QC samples or RSD > 30% in QC samples was removed [
50].
2.7. Data Analysis
Differentially expressed metabolites were presented by Venn diagrams drawn by online jvenn software (
http://genoweb.toulouse.inra.fr:8091/app/index.html). Enriched KEGG pathway mapping was performed by the online database (
http://www.genome.jp/kegg/). PCA, PLSDA, and heatmap analysis were done by Metaboanalyst 6.0 (
https://www.metaboanalyst.ca). The metabolites with VIP > 1 and p < 0.05 (student t test) were considered as significantly changed metabolites. For element analysis, statistical analysis was performed with statistical software package (Statistix 8.1v) using two-way analysis of variance (ANOVA). The significant difference between treatments was observed using Fisher’s least significant difference (LSD) test at 95 % probability.
4. Discussion
Leaf nutrient concentration and resorption play an important role in determining plant nutrient-use strategies [
51]. High leaf nutrient concentration and low leaf nutrient resorption generally reflect a “resource spending” nutrient use strategy for plants to quickly grow and develop [
52]. However, low leaf nutrient concentration and high leaf nutrient resorption usually reflect a “conservative consumption” nutrient use strategy for plant survival and reproduction [
53]. Leaf nutrient concentration and stoichiometry also reflect plant nutrient status and the strength of nutrient limitation in plants [
54,
55]. In our study, we found dramatic reduction in shoot concentrations of macro and micro nutrients in glycine and ammonium treatments relative to other two N-form treatments. Still there was significant difference in the two genotypes in nutrient concentrations in each N treatment. XZ16 had higher concentrations of nutrient ions in glycine treatment as compared to Hua30, but Hua30 accumulated more nutrients in ammonium treatment (
Table 1). These results are consistent with our previous results [
49] where some parameters like leaf surface area and tiller per plant increased in ammonium treatment for the cultivated barley as compared to wild genotypes. Similarly higher biomass of wild genotypes was found in urea treatment in our previous study and here also we found higher nutrient concentrations in both the genotypes in urea treatment in comparison with other two treatments, with XZ16 being higher than Hua30. Dubey et al. [
56] found that ammonium nitrogen not only induced disorder in nutrient balance but also reduced calcium (Ca) uptake. In an another study more ammonium uptake resulted in reduced plants shoot Ca and Mg concentrations and increased Mn and Al uptake [
57,
58].
Appropriate nitrogen fertilizer input can significantly increase the activity of enzymes related to nitrogen metabolism in leaves [
59], and promote photosynthesis and biosynthesis of organic acids, amino acids, proteins and other nitrogen-containing secondary metabolites in plants [
60]. In the present study, we found different pattern of carbohydrates, amino acids and organic acids among four nitrogen forms. In root metabolomics we had found down regulation of most carbohydrates, amino acids and organic acids in glycine treatment, but in shoot metabolomics, we found the opposite trend. Carbohydrate metabolism in shoots was greatly increased in ammonium and urea treatments in both the genotypes, with urea treatment being higher than ammonium treatment. But we can observe the difference between two genotypes (
Table S1). In ammonium treatment some commonly regulated sugar metabolites such as turanose, allose, D-(glycerol 1-phosphate), 3,6-anhydro-d-galactose had greater fold change in Hua30 (11.13, 3.23, 1.84, 8.12) as compared to XZ16 (2.57, 1.68, 1.94, 1.62), respectively (
Table S1). As osmolytes and signaling molecules, soluble sugars participate in the response and adaptation of plants to environmental stresses [
61]. The beta-mannosylglycerate is a potential protein thermostabilizer [
62] and significantly up-regulation of this metabolite in Hua30 may increase leaves protein stability in ammonium treatment [
63]. More accumulation of soluble sugar may confer higher resistance to abiotic stresses in ammonium treatment to Hua30 as compared to XZ16. From figure 4, it is clear that glycolysis intermediate (fructose-6-P) was more up-regulated and pyruvate was less down-regulated in XZ16 as compared to Hua30, thus adding the property of more glycine utilization in XZ16. Similarly, TCA cycle intermediates (aconitic acid and isocitrate) were more up-regulated and citrate, alpha ketoglutarate, and fumarate were less down-regulated in XZ16 as compared to Hua30, which may also be the contributive metabolites in glycine for more uptakes (
Figure 4). Overall, glycine reduced TCA intermediate in both the genotypes and these results are consistent with former studies conducted by Wang et al. [
64] and Liu et al. [
65], which partly explains the lower N use efficiency in response to glycine [
64].
The amino acid pool of plants is characterized by high flexibility and varies among species, cell types, and physiological situations [
39,
66,
67,
68]. In the environments of nutritional, low temperature and drought stresses, it might stimulate the synthesis of some amino acids [
60,
69]. Amino acids up-regulation may be also due to degradation of proteins to single amino acids, which can provide carbon skeleton for TCA cycle for enhancement of energy under stress [
70]. In the current study, we found more accumulation of amino acids in ammonium and glycine treatments in both the genotypes (
Table S1) with Hua30 being more affected than XZ16 under glycine treatment. Perchlik and Tegeder [
71] also found the increased amino acid delivery to
aap2 leaves, resulting in an enhanced sink C supply for higher NUE under N limitation. Some studies support that amino acid transport processes in the shoot and/or associated changes in the shoot N status exert regulatory control over N uptake in roots [
72,
73,
74,
75]. Stress responsive amino acids such as proline, asparagine, homoserine were greatly induced in ammonium and glycine treatments but less responsive in urea treatment. Accumulation of asparagine is regarded as a self-regulation mechanism for plant adaptation to organic N sources [
32,
76], this effect was also observed in our study (
Figure 4). Intermediates of arginine and proline metabolism (L-ornithine and putrescine) and 4-aminobutanoic acid are associated with plant stress responses [
77,
78,
79]. Glycine treatment led to higher concentrations of these components (
Figure 4), suggesting that supply of glycine as an N source may increase stress resistance in barley as reported in lettuce earlier [
80]. Considering these changes in plant metabolism, it is expected that glycine fertilizers could enhance the levels of antioxidants which might improve tolerance to abiotic stresses [
81].
Certain metabolites; including carbon compounds from photosynthesis, starch and sucrose metabolism plus oxidative pentose phosphate pathway, tricarboxylic acid (TCA) cycle and glycolysis metabolites and N-containing metabolites such as glutamine, have been seen to be at low levels in plants undergoing high growth. This suggests that these metabolites are providing the major building blocks for bio-macromolecules such as proteins and that growth drives metabolism [
82]. In urea treatment, amino acid and organic acids were greatly down-regulated in XZ16 as compared to Hua30. Urea is a good N fertilizer form to barley only after nitrate [
83]. Thus for the barley seedlings exposed to urea treatment, it is not necessary for them to accumulate more amino acids or organic acids as stress responsive. Also the amino acids produced may be directly used for protein synthesis that may be the reason for less accumulation of amino acids. Furthermore more down regulation of these metabolites in XZ16 may be related to the more biomass production in urea treatment as reported earlier [
83].
Glyoxylate and dicarboxylate metabolism is an important abiotic stress-related pathway, providing a metabolic balance to improve tolerance [
84]. In glycine treatment this was greatly affected (
Table S6). Oxalic acid, isocitrate, and glutamin were significantly up-regulated in XZ16 as compared to Hua30 in glycine of this pathway which may be the reason for more utilization of glycine in XZ16. Similarly arginine and proline pathway is related to nitrogen metabolism in plants, essential for production of nucleic acids and proteins. Arginine is a precursor of polyamines and has a role in proline biosynthesis when glutamate is not available. In ammonium treatment this pathway is greatly affected (
Table S4). The influence of ammonium treatment is highly expected in this pathway. Neto et al. [
85] also found these pathways as stress responsive.