- State Key Laboratory of Crop Biology, College of Agronomy, Shandong Agricultural University, Taian, China
Wheat growth and nitrogen (N) uptake gradually decrease in response to high NH4+/NO3– ratio. However, the mechanisms underlying the response of wheat seedling roots to changes in NH4+/NO3– ratio remain unclear. In this study, we investigated wheat growth, transcriptome, and proteome profiles of roots in response to increasing NH4+/NO3– ratios (Na: 100/0; Nr1: 75/25, Nr2: 50/50, Nr3: 25/75, and Nn: 0/100). High NH4+/NO3– ratio significantly reduced leaf relative chlorophyll content, Fv/Fm, and ΦII values. Both total root length and specific root length decreased with increasing NH4+/NO3– ratios. Moreover, the rise in NH4+/NO3– ratio significantly promoted O2– production. Furthermore, transcriptome sequencing and tandem mass tag-based quantitative proteome analyses identified 14,376 differentially expressed genes (DEGs) and 1,819 differentially expressed proteins (DEPs). The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis indicated that glutathione metabolism and phenylpropanoid biosynthesis were the main two shared enriched pathways across ratio comparisons. Upregulated DEGs and DEPs involving glutathione S-transferases may contribute to the prevention of oxidative stress. An increment in the NH4+/NO3– ratio induced the expression of genes and proteins involved in lignin biosynthesis, which increased root lignin content. Additionally, phylogenetic tree analysis showed that both A0A3B6NPP6 and A0A3B6LM09 belong to the cinnamyl-alcohol dehydrogenase subfamily. Fifteen downregulated DEGs were identified as high-affinity nitrate transporters or nitrate transporters. Upregulated TraesCS3D02G344800 and TraesCS3A02G350800 were involved in ammonium transport. Downregulated A0A3B6Q9B3 is involved in nitrate transport, whereas A0A3B6PQS3 is a ferredoxin-nitrite reductase. This may explain why an increase in the NH4+/NO3– ratio significantly reduced root NO3–-N content but increased NH4+-N content. Overall, these results demonstrated that increasing the NH4+/NO3– ratio at the seedling stage induced the accumulation of reactive oxygen species, which in turn enhanced root glutathione metabolism and lignification, thereby resulting in increased root oxidative tolerance at the cost of reducing nitrate transport and utilization, which reduced leaf photosynthetic capacity and, ultimately, plant biomass accumulation.
Introduction
Wheat (Triticum aestivum L.) is one of the most important cereal foods worldwide. The tremendous increase in wheat crop productivity over the past several decades has made China the largest wheat consumer and producer in the world (Xu et al., 2013), mainly because of the substantial increase in the use of synthetic fertilizers, specifically, nitrogen (N) (Zhu and Chen, 2002; Yuan and Peng, 2017).
Plants have evolved sophisticated root regulatory mechanisms and N transport systems for soil N uptake (Luo et al., 2020). Although previous studies showed that urea was directly absorbed by urea transporters, such as DUR3 and MIP (Wang et al., 2012; Liu et al., 2015), it is widely accepted that plants absorb mainly ammonium or nitrate generated by microbial conversion or via their soil application as fertilizers (Witte, 2011). Two soil nitrate and ammonium uptake systems have been identified in higher plants, namely, the high-affinity transport system (HATS) and the low-affinity transport system (LATS) (Kiba and Krapp, 2016; Gojon, 2017; Plett et al., 2018). NRT1 and NRT2 are operated by two families of nitrate transporters, whereas members of the ammonium transporter (AMT) family mediate ammonium transport (Kraiser et al., 2011) and are involved in diverse aspects of plant growth and development. Although N transporters have also been reported to have specific effects on lateral root initiation (Remans et al., 2006; Motte et al., 2019), there are few reports on the effects of different NH4+/NO3– ratios on changes in N transport and uptake and metabolism in roots of young wheat seedlings.
Farmers apply excessive amounts of N fertilizers, such as urea, ammonium bicarbonate, and ammonium nitrate, simply following the conventional concept of “the more fertilizer, the higher yield.” However, studies have shown that grain yield does not increase accordingly with the excess of applied N (Ju et al., 2009; Cui et al., 2010; Liu et al., 2016). Indeed, although large amounts of N fertilizer are regularly applied as base fertilizer to increase the NH4+-N concentration and, thus, change the ratio of NH4+ and NO3– in the soil (Husain et al., 2019; Kirschke et al., 2019; Wang et al., 2020), a high-NH4+ concentration may be a disadvantage for plant growth. Wheat was shown to have a specific preference for NO3–, and in fact, showed toxicity symptoms under high NH4+ concentration (Imran et al., 2019). Consistently, plant biomass was significantly reduced when N fertilizer was supplied as NH4+ (Wang et al., 2016a), likely because NH4+ reduced plant photosynthetic capacity due to damage to the electron transport chain (Wang et al., 2019). On the other hand, NH4+-N promoted over-accumulation of glutamate, which can inhibit pyruvate kinase activity, causing a reduction in tricarboxylic acid flux (Wang et al., 2020). Furthermore, NH4+-N caused oxidative stress due to the enhanced generation of reactive oxygen species (ROS), resulting in root growth inhibition (Liu et al., 2021). Fortunately, plants possess defense mechanisms to cope with oxidative stress, whereby cells are protected against potential oxidative damage caused by ROS.
Glutathione metabolism is one of the defense mechanisms that play a vital role in plant resistance to oxidative stress (Gong et al., 2018). Plant glutathione and glutathione S-transferases (GSTs) are widely involved in ROS detoxification (Zechmann, 2014; Ding et al., 2017). As a potent antioxidant, glutathione can control cellular ROS levels (Gill et al., 2013). GSTs catalyze the nucleophilic fusion of reduced glutathione with electrophilic and hydrophobic toxic molecules generated under stress to convert them into non-toxic soluble conjugates (Tiwari et al., 2016). Lignin deposition is another defense mechanism that responds to a wide range of abiotic stress conditions (Cesarino, 2019) to which plants are often exposed during their life cycle (Moura et al., 2010). For example, salt and drought stress upregulated the expression of the phenylalanine ammonia-lyase (PAL), 4-coumarate-CoA ligase (4CL), cinnamoyl-CoA reductase (CCR), and cinnamyl alcohol dehydrogenase (CAD) genes, resulting in an increase in lignin biosynthesis, enhancing stress resistance (Liu et al., 2018; Xu et al., 2020). Although glutathione metabolism and lignin biosynthesis have been reported during plant development and under abiotic stress conditions, limited information is available regarding the involvement of these two defense strategies in wheat roots under NH4+ stress. In particular, the mechanisms underlying the responses of glutathione metabolism and lignin biosynthesis in wheat root growth to different NH4+/NO3– ratios are still far from being clearly understood. Therefore, the objectives of the present study were to (i) identify the effects of different NH4+/NO3– ratios on plant and root growth at the seedling stage, (ii) investigate the wheat root transcriptome and proteome profiles related to glutathione metabolism and lignin biosynthesis in response to different NH4+/NO3– ratios, and (iii) provide new insights into the wheat root response to N for the improvement of N fertilizer application based on a scientific N management strategy.
Materials and Methods
Plant Materials and Culture Conditions
Experiments were carried out at the experimental station of Shandong Agricultural University, Tai’an, China (36°09′N, 117°09′E, 128 m above sea level). The Jimai 22 (JM22) winter wheat cultivar was grown under controlled conditions. Variations in temperature, illumination intensity, and relative humidity are shown in Supplementary Figure 1. Seeds were surface-sterilized with 1% NaClO3 for 30 min, rinsed 5 × with sterile water, and germinated on wet filter paper in the dark at 25°C for 2 days. Germinated seeds were sowed in plastic pots (10 cm × 10 cm × 10 cm; 10 plants/pot) filled with perlite. The pots were placed on plastic trays (55 cm × 45 cm × 5 cm; 20 pots/tray).
Treatments and Experimental Design
Five NH4+/NO3– ratio treatments with a 6-mM total N concentration were tested in a controlled climate chamber experiment. Treatments included: 100/0 (Na), 75/25 (Nr1), 50/50 (Nr2), 25/75 (Nr3), and 0/100 (Nn). Ca(NO3)2, NH4Cl, or NH4NO3 were used to set the ratios. The NH4+/NO3– ratio preparation method is shown in Supplementary Table 1. The concentration of all other nutrient elements was referred to a modified Hoagland’s nutrient solution with the following chemical composition: 4 mmol L–1 Ca2+, 5 mmol L–1 KCl, 4 mmol L–1 MgSO4, 1 mmol L–1 KH2PO4, 50 μmol L–1 Fe-ethylenediaminetetraacetic acid (EDTA), 0.5 μmol L–1 H3BO3, 0.74 μmol L–1 MnSO4, 0.27 μmol L–1 ZnSO4, 0.001 μmol L–1 CuSO4, 0.001 μmol L–1 CoCl2, 0.005 μmol L–1 Na2MoO4, and 0.025 μmol L–1 KI; and pH 6.5. The nutrient solutions were supplemented once a day (1.5 L/tray). Wheat growth was evaluated from sowing (0 days) to 42 days after sowing (DAS). We observed that the first leaf emerged through the coleoptile at 7 DAS and unfolded at 14 DAS. The morphological characteristics of 42-day-old wheat plants under Nn treatment showed that four leaves were unfolded, and the fifth leaf was visible (Supplementary Figure 2). However, only three leaves were completely unfolded on 42-day-old wheat plants under Na treatment.
Measurements of the Shoot and Root Biomass Weight and Its Nitrogen Content Measurement
Ten wheat plants of each treatment were randomly sampled from 14 to 42 DAS at 7-day intervals. Sampled plants were divided into shoot and root. Samples were dried under 60°C for plant biomass weight measurement [g/plant dry weight (DW)].
The dried samples were pulverized, and their N content was measured using the semi-micro Kjeldahl method (Zou, 2000). The powder (0.1 g) was mixed with 3.2 g of the catalyzer (CuSO4/K2SO4, w/w: 1/15) and digested in 10 ml H2SO4 at 420°C for 2 h and cooled. Distillation was used as a FOSS Kjeltec 8200 analyzer (FOSS, Sweden). The boric acid is used as the receiving solution. The solution was titrated with sulfuric acid. The N content was calculated according to the following formula:
Where Vsample is the volume of used sulfuric acid for the sample (ml), Vblank is the volume of used sulfuric acid for blank (ml), N of acid represents normality of the used acid (mol L–1), and W is the sample dry weight (g).
Measurements of Leaf Relative Chlorophyll Content and Chlorophyll Fluorescence Analysis
Eight intact fully expanded leaves of each treatment at 42 DAS were selected to measure relative chlorophyll content and the chlorophyll fluorescence analysis using the MultispeQ device (PhotosynQ, United States) according to the method described by Kuhlgert et al. (2016).
Measurements of Root Growth
Root growth was evaluated by image analysis. Eight wheat plants of each treatment were sampled at 42 DAS. Roots were rinsed 3× with sterile water. Root sample images were captured with an EPSON V750 scanner (Seiko Epson Corp., Suwa, Nagano, Japan) at a resolution of 300 dpi. The images were analyzed in WinRHIZO (v. 4.0; Regent Instruments, Sainte-Foy, Quebec, Canada) which provided the total root length (TRL), root number (RN), and root average diameter (RAD). The root length and root biomass data were used to calculate the specific root length (SRL, TRL per unit of root dry weight, m g–1 DW).
Root Nitrate and Ammonium Content Measurement
Roots of 10 wheat plants under each treatment were sampled from 14 to 42 DAS at 7-day intervals. The fresh root nitrate content was determined using the salicylic acid method, as described by Cataldo et al. (1975). Three replicates of 0.1 g wheat root from each treatment were grounded to a powder in liquid N2 and homogenized in 1 ml redistilled water. Samples were boiled for 30 min, cooled, and centrifuged at 12,000 × g for 15 min. For each sample, 0.1-ml supernatant was mixed with 0.4 ml salicylic acid–sulfate acid (5 g salicylic acid in 100 ml sulfate acid). The reactions were incubated at room temperature for 20 min, and 9.5 ml of 8% NaOH solution was added. After cooling the tube to room temperature, the optical density (OD)410 value was measured.
Ammonium concentrations in roots were measured using the O-phthalaldehyde (OPA) fluorometric method, as described by Schjoerring et al. (2002) and Wang et al. (2016b). Roots (0.5 g) were homogenized in liquid N2 and added 5 ml of 10 mM formic acid to extract NH4+. The root homogenates were centrifuged at 10,000 × g (4°C) for 15 min. Supernatants were transferred to 5 ml polypropylene tubes with 0.45-μm organic ultrafiltration membranes and centrifuged at 50,000 × g (4°C) for 10 min. A total of 100 ml of OPA reagent was prepared by combining 0.2 M potassium phosphate buffer (composed of equimolar amounts of potassium dihydrogen phosphate and potassium monohydrogen phosphate), 3.75 mM OPA, and 2 mM 2-mercaptoethanol, 1 day before use. Prior to adding 2-mercaptoethanol, the solution pH was adjusted to 7.0 using 1 M NaOH and filtered through two layers of filter paper. Extract (10 μl) was added to 3 ml of OPA reagent, and color was allowed to develop for 30 min in the absence of light (25°C) before sample absorbance was measured at 410 nm.
Superoxide Anion Radical Measurement
Wheat roots under each treatment were sampled from 14 to 42 DAS at 7-day intervals. Anion radical (O2–) content and production rate were determined according to the method described by Bai et al. (2015). Three biological replicates of 0.5 g root samples for each treatment were homogenized in 3 ml of 65 mM phosphate buffer (pH 7.8), and the homogenate was centrifuged at 10,000 × g for 10 min. About 2 ml supernatant was added to 0.5 ml of 50 mM phosphate buffer (pH 7.8) and 0.1 ml of 10 mM hydroxylamine hydrochloride. After 20 min at 25°C, the mixture was added to 1 ml of 58 mM sulfanilamide and 1 ml of 7 mM α-naphthylamine, and let at 30°C for 30 min. After this period, the absorbance was measured at 530 nm.
Root Lignin Content Measurement
Wheat roots of each treatment were sampled from 14 to 42 DAS at 7-day intervals. Root lignin contents were determined according to a modified method described by Zheng et al. (2017). Three replicates of 0.1 g root for each treatment were grounded to a powder in liquid N2 and washed five times with 80% ethanol to remove soluble metabolites, followed by acetone wash and dried in a drying oven. Samples were transferred to 10-ml centrifuge tubes with 2.5 ml acetyl bromide and acetic acid solution (v/v, 1:4) and incubated at 70°C for 1 h. After cooling to room temperature, 0.9 ml of 2 M NaOH was added to terminate the reaction, then mixed and added 0.1 ml of 7.5 M hydrochloride. Finally, 4 ml acetic acid was added to each tube. The solution absorbance was measured at 280 nm using a spectrophotometer. The lignin contents were calculated according to the following linear calibration curve:
Lignin content (mg g–1 FW) = (ΔA280 − 0.0068) ÷ 0.035 × V × 10–3 ÷ W × T.
Where ΔA280 is the absorbance at 280 nm, W is sample fresh weight in g, and T is the dilution factor. Calibration curves were generated with increasing concentrations of standard lignin (Sigma-Aldrich, Inc., St. Louis, MO, United States) (4, 8, 12, 16, 20, and 24 mg L–1), which were processed by the same method used for plant samples.
Lignin Biosynthesis-Related Enzyme Activity Determination
Wheat roots of each treatment were sampled from 14 to 42 DAS at 7-day intervals. About 0.5 g root samples with three replicates from each treatment were ground to a powder in liquid N, and extracted with pH 7.0 phosphate buffer. The activities of PAL (EC:4.3.1.24), 4CL (EC:6.2.1.12), CCR (EC:1.2.1.44), CAD (EC:1.1.1.195), caffeic acid 3-O-methyltransferase (EC:2.1.1.68, COMT), and peroxidase (EC:1.11.1.7, POD) were measured using enzyme activity assay kits (Caobenyuan Biotechnology Co., Ltd., Nanjing, China).
RNA Sequencing and Data Analysis
Roots were sampled in each treatment group at 42 DAS, immediately frozen in liquid N, and stored at −80°C until the RNA extraction and proteome analysis. Total RNA was extracted using the TRIzol reagent (Invitrogen, Carlsbad, CA, United States). Total RNA quantity and purity were determined using a NanoDrop and Agilent Bioanalyzer 2100 system (Agilent Technologies). Approximately 10 μg of total RNA representing a specific adipose type was subjected to isolation of poly (A) mRNA with poly T oligo-attached magnetic beads (Invitrogen, Carlsbad, CA, United States). Following purification, mRNA was fragmented into small pieces using divalent cations at high temperatures. Cleaved RNA fragments were reverse-transcribed to create the final cDNA library in accordance with the protocol by the manufacturer for the mRNA-Seq Sample Preparation Kit (Illumina, San Diego, CA, United States); the average insert size for the paired-end libraries was 300 bp (±50 bp). Finally, we performed paired-end sequencing on an Illumina HiSeq 4000 at the LC Sciences, San Diego, CA, United States, following the protocol recommended by the vendor.
The adaptor reads were removed from raw data using the Cutadapt software1. Clean data were obtained by removing the low-quality and repeat reads. Sequence quality was verified using the FastQC software2, including Q20 (phred quality score 20, percentage of bases whose base call accuracy exceeds 99%), Q30 (phred quality score 30, percentage of bases whose base call accuracy exceeds 99.9%), and guanine-cytosine (GC) content. We aligned the sample reads to the wheat reference genome of the International Wheat Genome Sequencing Consortium (IWGSC) RefSeq v1.0 (IWGSC 2018) using HISAT2. Mapped reads from each sample were assembled using the StringTie software. The value for fragment per kilobase of transcript per million mapped reads (FPKM) was calculated to quantify its expression abundance and variations using the StringTie software (version v1.3.3)3. DESeq2 software4 was used to analyze differentially expressed genes (DEGs) selected with the parameter false discovery rate (FDR) < 0.05, and an absolute fold change of ≥2. Gene ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) databases were used to reveal the functional enrichment of DEGs.
Protein Extraction and Identification
Root samples were ground under liquid N, and 0.8 g of the powder obtained was transferred to a 5-ml centrifuge tube, homogenized with 1 ml lysis buffer, and placed on ice for 10 min. Tris–phenol (1 ml) was added to the tube, homogenized, and placed again on ice for 10 min. Samples were then centrifuged at 16,000 × g for 10 min at 4°C. The supernatant was transferred into a new tube, and 4× volumes of cold ammonium acetate/methanol were added, followed by incubation at −20°C overnight. Then, samples were centrifuged at 16,000 × g for 10 min at 4°C, the supernatant was discarded, cold methanol was added, and the mixture was incubated at − 20°C for 1 h. After centrifugation at 16,000 × g for 10 min at 4°C, the sediment was collected, resuspended in cold acetone, washed two times, and air-dried. Samples were dissolved in 0.6 ml lysis buffer [8 M urea, 50 mM pH 8.0 Tris–HCl, 1% NP40, 1% sodium deoxycholate (NaDOC), 10 mM EDTA, 5 mM dithiothreitol (DTT), 1% protease inhibitor mixture, and 1% phosphatase inhibitor mixture], and centrifuged at 20,000 × g for 10 min at 4°C. The supernatant was collected, and the protein concentration was quantified using the Bio-Rad Protein Assay Kit (Bio-Rad, Hercules, CA, United States) following the instructions of the manufacturer. Then, 30 μg of protein from each sample was used to determine the quality of the proteins using 12% sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE).
About 500 μg of protein from each sample was transferred to a new tube, adjusted to an equal volume with lysis buffer, and incubated at 30°C for 30 min. After cooling to room temperature, 30 mM iodoacetamide (IAM) was added, and the mixture was incubated at room temperature for 45 min in the dark, before adding 5× volumes of cold acetone, followed by incubation at −20°C overnight. After centrifugation at 20,000 × g for 10 min at 4°C, the sediment was collected, resuspended in cold 80% acetone, washed two times, and air-dried, followed by the addition of 300 μl of 0.1 M triethylammonium bicarbonate (TEAB) to the proteins and sonication on ice for 5 min. Trypsin (20 μg) was added to the samples and incubated at 37°C overnight. After incubation, 1% trifluoroacetic acid (TFA) was added to the mixture to terminate the reaction. The peptides were desalted using a C18 SPE column (Phenomenex, Torrance, CA, United States) and vacuum-dried prior to reconstitution in 0.5 M TEAB and processing according to the protocol of the manufacturer for tandem mass tag (TMT) kit (Thermo Fisher Scientific, Waltham, MA, United States). Tryptic peptides were fractionated by high-pH reverse-phase high-performance liquid chromatography (HPLC) using a Waters XBridge Shield C18 RP column (3.5 μm particles, 4.6 mm × 250 mm; Waters Corp., Milford, MA, United States). Peptides were separated into 60 fractions, combined into 20 fractions, and vacuum centrifugation-dried. The fractions were dissolved in 0.1% formic acid, centrifuged at 20,000 × g for 2 min at 4°C, and the supernatant was transferred to a sampler vial. Liquid chromatography with tandem mass spectrometry (LC–MS/MS) analysis was performed by Micrometer Biotech Company (Hangzhou, China).
Acquired MS/MS spectra were searched against the Triticum aestivum UniProt database (143,807 sequences). At least one unique peptide with an FDR of ≤1% was required for analysis of the protein identification and quantification data. For analysis of differentially expressed proteins (DEPs), only proteins with quantitative information from at least two biological replicates were used. The average of three biological replicates was used to indicate final protein abundance; finally, proteins showing an average protein-abundance change significantly greater than 1.3-fold (ANOVA test, p ≤ 0.05) were defined as DEPs.
Statistical Analysis and Processing
Data for leaf chlorophyll content and chlorophyll fluorescence parameters, TRL, RAD, SRL, N content, O2– content, lignin content, and lignin biosynthesis-related enzyme activity were processed in DPS v. 7.05 (Hangzhou, China). Multiple comparisons were performed after a preliminary F-test. The means were tested with a least significant difference test, and the significance was set at the probability level of 0.05. Graphs were plotted using OriginPro 2017 (OriginLab Inc., Northampton, MA, United States).
Results
Wheat Growth Responses to Increasing NH4+/NO3– Ratios
Overall, NH4+/NO3– ratio affected shoot and root biomass and morphological traits of wheat seedlings (Figure 1 and Supplementary Figure 2). Shoot and root biomass decreased with the increment of NH4+/NO3– ratio. In particular, the shoot and root biomass exhibited significant differences among the five treatments at 35 and 42 DAS. For example, compared with the Nn treatment, shoot biomass in Na, Nr1, Nr2, and Nr3 decreased by 40.56, 28.14, 19.66, and 7.79% at 42 DAS, respectively, whereas root biomass decreased by 50.66, 39.67, 23.93, and 9.84%, respectively (Figure 1).
Figure 1. The effects of different NH4+/NO3– ratios on the wheat shoot (A) and root biomass (B) from 14 to 42 days after sowing (DAS). Segments represent the SEM (n = 10). Different letters indicate significant differences among each treatment, p < 0.05.
As the NH4+/NO3– ratio increased, the leaf relative chlorophyll content at 42 DAS decreased significantly (Figure 2A). In particular, the relative chlorophyll content decreased by 57.18% in the Na treatment compared with the Nn treatment. In addition, the chlorophyll fluorescence analysis showed that both the maximal quantum efficiency of photosystem II (PSII) (Fv/Fm) and the steady-state efficiency (ΦII) decreased with increasing NH4+/NO3– ratios (Figures 2B,C).
Figure 2. The effects of different NH4+/NO3– ratios on leaf relative chlorophyll content (A), Fv/Fm (B), ΦII (C), total root length (D), specific root length (E), and root average diameter (F). Segments represent the SEM (n = 8). Different letters indicate significant differences among each treatment, p < 0.05.
TRL and SRL decreased significantly as the NH4+/NO3– ratio increased. Furthermore, the Nn treatment resulted in longer TRL and SRL but less RAD (Figures 2D–F). In contrast, the Na treatment significantly (p < 0.05) reduced TRL and SRL by 77.98% and 48.59%, respectively, compared with the Nn treatment. Meanwhile, RAD increased with increasing NH4+/NO3– ratios. Indeed, the RAD of the Nr1 treatment increased by 280.21% compared to Nn.
Root Transcriptomic and Proteomic Responses to Increasing NH4+/NO3– Ratios
Transcriptome sequencing generated from 3.80 to 5.50 million raw data. After quality control, over 97% of the reads with Q30 were saved for further analysis. On average, more than 70% of these reads were specific to wheat and uniquely aligned to the reference genome (Supplementary Table 2).
A total of 7,650 (upregulated) and 6,726 (downregulated) DEGs were identified in the four pairwise comparisons listed in Supplementary Table 3. Concomitantly, a total of 1,042 (upregulated) and 777 (downregulated) DEPs were identified (Supplementary Table 4). Compared with the Nn treatment, increasing the NH4+/NO3– ratio significantly increased the number of DEGs and DEPs (Figures 3A, 4A). A total of 2,834 DEGs (2,066 upregulated and 768 downregulated) and 357 DEPs (267 upregulated and 90 downregulated) were identified in the Nr3/Nn group. Additionally, a total of 11,922 DEGs (6,406 upregulated and 5,516 downregulated) and 1,471 DEPs (848 upregulated and 623 downregulated) were detected in the Na/Nn group. The results of the Venn diagrams showed that 1,108 upregulated DEGs and 1,372 downregulated DEGs, and 128 upregulated DEPs and 182 downregulated DEPs were unique to Na/Nn (Figures 3B, 4B). Furthermore, a total of 266 upregulated DEGs and 117 downregulated DEGs, and 37 upregulated DEPs and 13 downregulated DEPs were unique to Nr3/Nn. The four groups shared 1,486 upregulated DEGs and 494 downregulated DEGs, and 178 upregulated DEPs and 46 downregulated DEPs.
Figure 3. (A) Identification and comparison analysis of differentially expressed genes (DEGs) in the wheat root samples of different NH4+/NO3– ratio treatments. (B) Venn diagrams represent up and downregulated genes in four comparisons.
Figure 4. (A) Identification and comparison analysis of differentially expressed proteins (DEPs) in the wheat root samples of different NH4+/NO3– ratio treatments. (B) Venn diagrams represent up and downregulated proteins in four comparisons.
NH4+/NO3– Ratio Changed Kyoto Encyclopedia of Genes and Genomes Enrichment Pattern of Differentially Expressed Genes and Differentially Expressed Proteins in Wheat Seedling Roots
Thirty-seven significant KEGG pathways were enriched for DEGs in the Na/Nn transcriptome comparison group (Supplementary Table 5), whereas 33 KEGG pathways were enriched in the Nr3/Nn comparison group. Additionally, 24 pathways were shared, including glutathione metabolism (ko00480) and phenylpropanoid biosynthesis (ko00940) (Supplementary Figure 3A and Supplementary Table 5). Concomitantly, 13 significant KEGG pathways were enriched in the Na/Nn proteome comparison group, but only six pathways were enriched in the Nr3/Nn group (Supplementary Table 6). Four pathways (ko00330, arginine, and proline metabolism; ko00620, pyruvate metabolism; ko00950, isoquinoline alkaloid biosynthesis; ko00960, tropane, piperidine, and pyridine alkaloid biosynthesis) were unique to Na/Nn, and two pathways (ko00900, terpenoid backbone biosynthesis, and ko00780, biotin metabolism) were unique to Nr3/Nn (Supplementary Figure 3B and Supplementary Table 6). Additionally, only glutathione metabolism (ko00480) and phenylpropanoid biosynthesis (ko00940) were shared among the four comparison groups.
Glutathione Metabolism-Related Differentially Expressed Genes and Differentially Expressed Proteins Responses to Increasing NH4+/NO3– Ratios
The four comparison groups identified a total of 236 DEGs and 53 DEPs related to the glutathione metabolism pathway (00480) (Supplementary Tables 7, 8). Thirty-one DEGs and nine DEPs were unique to Na/Nn (Supplementary Figure 4). Additionally, 82 DEGs and 12 DEPs were shared among the four comparison groups, all of which were upregulated, in which 81 DEGs were GSTs (EC:2.5.1.18) and only one gene (TraesCS1A02G095800) was an AP domain-containing protein (EC:3.4.11.1). Although only two proteins (A0A3B6K9P2 and Q8RW00) were GST, eight uncharacterized genes (A0A3B5XZ69, A0A3B6ER15, A0A3B6H0S6, A0A3B6I5M9, A0A3B6ILA6, A0A3B6JJE1, A0A3B6LUE0, and W5CW67) had GST molecular function, accordingly with the GO enrichment analysis. Finally, A0A3B6MRF8 may play a role in zinc-ion binding, whereas A0A3B6JL78 has POD activity.
Phenylpropanoid Biosynthesis-Related Differentially Expressed Genes and Differentially Expressed Proteins Responses to Increasing NH4+/NO3– Ratio
A total of 598 DEGs and 70 DEPs related to the phenylpropanoid biosynthesis pathway (00940) were identified in the four comparison groups (Supplementary Tables 9, 10), and 99 DEGs and 12 DEPs were shared among the four comparison groups (Supplementary Figure 5). Based on the gene annotation, TraesCS2D02G581200 was defined as 4CL (EC: 6.2.1.12). Seven upregulated genes were associated with shikimate O-hydroxycinnamoyl transferase (EC: 2.3.1.133, HCT). Fourteen upregulated genes and one downregulated gene were associated with CCR (EC: 1.2.1.44). One up and two downregulated genes were associated with COMT (EC: 2.1.1.68). Four upregulated and two downregulated genes were associated with CAD (EC: 1.1.1.195). Furthermore, 19 upregulated and 28 downregulated genes were associated with a POD (EC: 1.11.1.7). Based on the protein annotation, two PKS_ER domain-containing proteins (A0A3B6NPP6 and A0A3B6LM09) were obtained, whose amino acid sequences were aligned to those of the corresponding proteins in the wheat genome. Results of the phylogenetic tree constructed showed that both A0A3B6NPP6 and A0A3B6LM09 belong to the CAD subfamily (Supplementary Figure 6). Six upregulated and three downregulated proteins were associated with POD. Related genes encoding phenylalanine ammonia-lyase (PAL), 4CL, CCR, CAD, COMT, and POD were upregulated as NH4+/NO3– ratio increase, along with the related proteins (Figure 5).
Figure 5. Heatmaps represent the candidate genes expression (A) and proteins abundance (B) involved in lignin biosynthesis in different NH4+/NO3– ratios treatments. The gene expression was scaled using the fragment per kilobase of transcript per million mapped reads (FPKM) based on a mean value of three biological replicates. Different colors represent the upregulation and downregulation difference of genes and proteins among treatments. PAL, phenylalanine ammonia-lyase; 4CL, 4-coumarate-CoA ligase; CCR, cinnamoyl-CoA reductase; CAD, cinnamyl-alcohol dehydrogenase; COMT, caffeic acid 3-O-methyltransferase; POD, peroxidase.
Nitrogen Metabolism and Transport-Related Differentially Expressed Genes and Differentially Expressed Proteins in Response to Increasing NH4+/NO3– Ratios
A total of 68 DEGs and 12 DEPs involved in the N metabolism pathway (00910) were identified in all four comparison groups (Supplementary Tables 11, 12). Seven DEGs and three DEPs were unique to Na/Nn, and one DEG and one DEP were unique to Nr1/Nn (Supplementary Figure 7). Additionally, the four comparison groups shared 27 DEGs and three DEPs. Based on the gene annotation and description, two upregulated DEGs (TraesCS3D02G344800 and TraesCS3A02G350800) were involved in ammonium transport (AMT) in the Na/Nn, Nr1/Nn, and Nr2/Nn comparison groups (Supplementary Table 13). In the four comparison groups, 15 downregulated DEGs were high-affinity nitrate transporters or nitrate transporters, including NRT2.1 and NRT2.4. Furthermore, the upregulated TraesCS7D02G449400 in the four comparison groups belongs to the NRT1/PTR Family 5.1. TraesCS1B02G224900 (NRT1/PTR Family 6.3) was downregulated in the Nr3/Nn group but upregulated in the Na/Nn, Nr1/Nn, and Nr2/Nn comparison groups (Figure 6 and Supplementary Table 13). Five DEGs were involved in nitrate reduction, and seven showed carbonic anhydrase activity. Based on protein annotation, the downregulated A0A3B6Q9B3 gene was involved in the nitrate transport in the four comparison groups and was successfully mapped to TraesCS6D02G035600. A0A3B6PQS3 was involved in the ferredoxin-nitrite reduction and mapped to TraesCS6B02G364600. A0A3B6TX28 showed to have carbonic anhydrase activity and was mapped to TraesCS7D02G443400.
Figure 6. Heatmaps represent the candidate genes expression involved in nitrogen transport in different NH4+/NO3– ratios treatments. The gene expression was scaled using the FPKM based on a mean value of three biological replicates. Different colors represent the upregulation and downregulation difference of genes and proteins among treatments. AMT2, ammonium transporter 2; NRT1/PTR6.3, nitrate transporter 1/peptide transporter family 6.3; NRT1/PTR5.1, nitrate transporter 1/peptide transporter family 5.1; NRT2.1, high-affinity nitrate transporter 2.1; NRT2.4, high-affinity nitrate transporter 2.4.
Effect of Increasing NH4+/NO3– Ratios on Nitrogen, NO3–-N, and NH4+-N Contents in Wheat Seedlings
The result showed that the shoot N content of Na and Nr1 treatments gradually increased but sharply decreased from 35 to 42 DAS (Figure 7A). Nevertheless, the shoot N content of Nn, Nr2, and Nr3 treatments increased from 14 to 42 DAS. Finally, the shoot N content decreased significantly with increasing NH4+/NO3– ratios at 42 DAS. For example, the shoot N content decreased by 6.64 and 13.52% after Nr2 and Nr1 treatments compared with Nn, respectively. Similarly, the root N content of Nn and Nr3 treatments increased from 14 to 42 DAS (Figure 7B), whereas the root N content of Na, Nr1, and Nr2 treatments gradually increased, and they later decreased from 35 to 42 DAS. The root N content significantly decreased with increasing NH4+/NO3– ratios. Relative to the Nr3 treatment, the root N content decreased by 13.83 and 17.37% at 42 DAS in Nr1 and Na treatments, respectively.
Figure 7. The effects of different NH4+/NO3– ratios on the nitrogen content of the shoot and (A) root (B) from 14 to 42 days after sowing (DAS). Segments represent the SEM (mean ± SE, n = 3). Different letters indicate significant differences among each treatment, p < 0.05.
The root NO3–-N content gradually decreased with plant growth. Conversely, the NH4+-N content gradually increased with plant growth (Figure 8). Increasing NH4+/NO3– ratios significantly reduced the NO3–-N content in wheat seedling roots (Figure 8A). Thus, the Na treatment yielded the lowest root NO3–-N content, whereas the Nn treatment showed the highest NO3–-N content. Indeed, at 42 DAS, the NO3–-N content was 71.07% lower in the Nr1 treatment than in the Nr3 treatment. Increasing the NH4+/NO3– ratios significantly increased the NH4+-N content in roots. The highest NH4+-N content was observed for the Na treatment, which was 96.35% higher than that of the Nn treatment at 42 DAS (Figure 8B).
Figure 8. The effects of different NH4+/NO3– ratios on root NO3–-N (A) and NH4+-N content (B) from 14 to 42 days after sowing (DAS). Segments represent the SEM (mean ± SE, n = 3). Different letters indicate significant differences among each treatment, p < 0.05.
Effects of the NH4+/NO3– Ratio on O2– Generation Rate and Content in Wheat Seedling Roots
Increasing the NH4+/NO3– ratio significantly promoted O2– generation and increased O2– content in wheat seedling roots (Figure 9). For example, the Na treatment yielded the highest O2– generation rate and content in the roots. Thus, compared with the Nn treatment, O2– content was enhanced by 104.20% in the Na treatment at 42 DAS; consistently, compared with the Nr3 treatment, the O2– generation rate was 50.41% higher under the Nr3 treatment at 42 DAS. Relative to the Nr3 treatment, the O2– content in the Nr1 treatment increased by 46.85% at 42 DAS.
Figure 9. The effects of different NH4+/NO3– ratios on root O2– generation rate (A) and O2– content (B) from 14 to 42 days after sowing (DAS). Segments represent the SEM (mean ± SE, n = 3). Different letters indicate significant differences among each treatment, p < 0.05.
Effects of the NH4+/NO3– Ratio on Lignin Content and the Activity of Lignin Biosynthesis-Related Enzymes in Wheat Seedling Roots
The root lignin content increased gradually with increasing growth (Figure 10); furthermore, it increased significantly with increasing NH4+/NO3– ratio. For instance, the lignin level for the Na treatment was 90.80% higher than that for the Nn treatment at 42 DAS. Similarly, compared with the Nr3 treatment, the lignin content was 45.46% higher in the Nr1 treatment at 42 DAS. The observed increased lignin contents may be due to the altered activities of lignin biosynthesis-related enzymes. Similarly, an increasing NH4+/NO3– ratio significantly increased PAL, 4CL, CCR, CAD, COMT, and POD activities (Figure 11). For instance, at 42 DAS, root PAL, 4CL, CCR, CAD, COMT, and POD activities in the Nr1 treatment increased by 36.43, 48.72, 21.55, 56.24, 11.57, and 44.79%, respectively, relative to the Nr3 treatment.
Figure 10. The effects of different NH4+/NO3– ratios on root lignin content from 14 to 42 days after sowing (DAS). Segments represent the SEM (mean ± SE, n = 3). Different letters indicate significant differences among each treatment, p < 0.05.
Figure 11. The effects of different NH4+/NO3– ratios on root PAL (A), 4CL (B), CCR (C), CAD (D), COMT (E), and POD (F) activities from 14 to 42 days after sowing (DAS). Segments represent the SEM (mean ± SE, n = 3). Different letters indicate significant differences among each treatment, p < 0.05. PAL, phenylalanine ammonia-lyase; 4CL, 4-coumarate-CoA ligase; CCR, cinnamoyl-CoA reductase; CAD, cinnamyl-alcohol dehydrogenase; COMT, caffeic acid 3-O-methyltransferase; POD, peroxidase.
Discussion
Effects of Increasing NH4+/NO3– Ratios on Wheat Seedling Growth
Nitrogen is an essential nutrient for plant growth, and consequently, the absorption of soil N is the basis for dry matter formation and accumulation. NH4+ and NO3– are the two main N forms plants can absorb and utilize (Miller and Cramer, 2005). As new roots and leaves emerge during growth, wheat seedlings become increasingly sensitive to N availability (Yang et al., 2017). In previous studies, the high NH4+-N concentration could produce apparent stress, inhibiting plant growth and reducing biomass (Wang et al., 2016b; Vega-Mas et al., 2019). However, the wheat growth response to different NH4+/NO3– ratios have different assessments. For example, Ijato et al. (2021) reported that ammonium-nitrate (NH4+/NO3– ratio: 1/1) was better for wheat growth. Conversely, the results of our study were not in line with this finding. In this study, we conducted experiments with different NH4+/NO3– ratio treatments to investigate the effects of the NH4+/NO3– ratio on wheat seedling growth. We found that the NO3–-N treatment resulted in the highest wheat growth. As the NH4+/NO3– ratio increased, significant growth inhibition occurred (Supplementary Figure 2). Although the sensitivity of different wheat cultivars to ammonium stress has been reported, wheat has been described as a species that is particularly sensitive to high NH4+ concentrations (Setién et al., 2013; Wang et al., 2016b; Liu et al., 2017). From this aspect, we used a wheat variety that is sensitive to the presence of ammonium. Indeed, wheat plants exhibited severe symptoms of ammonium toxicity. Moreover, the leaf chlorophyll content was significantly decreased with the increment in the NH4+/NO3– ratio (Figure 2A). This may be due to the activation of a series of genes by NH4+, including AMOS1/EGY1, which participated in the regulation of NH4+-stress signaling (Li et al., 2012; Li B. et al., 2014). High NH4+ damaged chloroplast ultrastructure and decreased the abundance of proteins and transcripts in chloroplasts, such as PsaA (Shi et al., 2020; Zhou et al., 2021). The chlorophyll content is a proxy of the leaf photosynthetic capacity (Croft et al., 2017), and our results showed that the lower chlorophyll content was unfavorable by the photosynthate accumulation. Furthermore, the chlorophyll fluorescence analysis indicated that leaf Fv/Fm and ΦII were significantly reduced by an increasing NH4+ supply (Figures 2B,C). Chlorophyll fluorescence traits have been used to evaluate the effects of abiotic stress factors on photosynthesis (Sharma et al., 2017). Thus, Fv/Fm reflects the maximum quantum efficiency of the photochemical reactions mediated by PSII (Sharma et al., 2015). In turn, ΦII correlates well with net photosynthetic rate and effective quantum yield of CO2 fixation (Krall and Edwards, 1992). Ammonium directly accelerates the photodamage of PSII and affects the repair of photodamaged PSII (Dai et al., 2014). Ammonium also stimulates electron transport to photorespiration to detriment of Rubisco carboxylation (Alencar et al., 2019). In this study, Fv/Fm and ΦII decrease, indicating that a high NH4+/NO3– ratio induces severe damage to the photosynthetic apparatus.
On the other hand, root growth was influenced by NH4+/NO3– ratio. Our study showed that a high NH4+/NO3– ratio changed root biomass and morphology (Figure 1B and Supplementary Figure 2). Furthermore, TRL and SRL were significantly reduced by increasing NH4+/NO3– ratios despite an increase in RAD (Figures 2D–F). Root architecture plays an essential role in N uptake (Meister et al., 2014). SRL and root diameter are considered useful root traits to evaluate nutrient absorption (Li H. et al., 2014; Sharma et al., 2021). Moreover, SRL positively correlated with the total amount of N uptake; in contrast, average root diameter negatively correlated with N uptake capacity (Hong et al., 2018). These results showed that increasing NH4+/NO3– ratios resulted in lower leaf photochemical efficiency and root uptake capacity, thereby reducing shoot and biomass and leading to stunted growth (Figure 1A and Supplementary Figure 2).
Increasing NH4+/NO3– Ratios Enhanced Root Glutathione Metabolism and Lignification, Thus Increasing Root Tolerance to Oxidative Damage but Reducing Nitrogen Transport and Utilization
In this study, we combined transcriptome with TMT-based proteome profiling to detect DEGs and DEPs and used functional predictions to explore enriched pathways that contribute to the root development under increasing NH4+/NO3– ratio conditions. A total of 14,376 DEGs and 1,819 DEPs were identified across the comparisons (Supplementary Tables 3, 4). Multiple approaches have shown that high ammonium concentrations are strongly phytotoxic and promote ROS formation, leading to oxidative stress (Bittsánszky et al., 2015; Podgórska and Szal, 2015). Increasing NH4+ supply may play a role in the oxidative stress by enhancing transcripts level and the NADPH oxidase activity, which is known as respiratory-burst oxidase homologs (RBOHs), which can oxidize cytosolic NADPH and transfer an electron across the cell membrane to generate O2– in the cytosol and cell wall (Swanson and Gilroy, 2010; Podgórska et al., 2021). In our study, we found that two upregulated DEGs were involved in the respiratory-burst oxidase expression (Supplementary Table 3). TraesCS1B02G295300 was upregulated in all four comparison groups. Compared with Nn, TraesCS1B02G295200 was upregulated in the Na, Nr1, and Nr2 treatments. This explains why increasing the NH4+/NO3– ratios significantly accelerated O2– generation and increased O2– content in the wheat roots (Figure 9). ROS act as signaling molecules regulating lateral root emergence and primary root elongation (Tsukagoshi, 2016; Liu and von Wirén, 2017). However, oxidative stress occurs when the enhanced production of ROS exceeds their degradation (Choudhary et al., 2020), thereby triggering several defense mechanisms in plants. Our KEGG pathway enrichment analysis showed that glutathione metabolism (ko00480) and phenylpropanoid biosynthesis (ko00940) were the main two shared pathways across comparisons (Supplementary Tables 5, 6). A total of 236 DEGs and 53 DEPs were involved in glutathione metabolism (Supplementary Tables 7, 8). Previous studies have found that GSTs are ubiquitous and multifunctional enzymes encoded by large gene families and have diverse roles in plant development and stress tolerance (Nianiou-Obeidat et al., 2017; Gullner et al., 2018). In our study, 81 upregulated genes were shared across comparisons and encoded GSTs involved in glutathione metabolism. According to the results of our proteome analysis, two shared upregulated proteins (A0A3B6K9P2 and Q8RW00) were GSTs. Previous studies reported that GSTs regulate root development by participating in auxin binding or transport (Jiang et al., 2010; Nongmaithem et al., 2021). GSTs also possess glutathione peroxidase activity, quench ROS with the addition of glutathione, and protect the cell from oxidative damage (Halušková et al., 2009; Kumar and Trivedi, 2018). We found that A0A3B6JL78 was described as a POD induced by a high NH4+/NO3– ratio. Thus, the increased GST gene and protein expression under a high NH4+/NO3– ratio may contribute to regulating root growth and avoiding O2– accumulation.
Recent reports have shown that root lignification is another defense mechanism that plays a critical role in adapting to oxidative stress (Cesarino, 2019; Yan et al., 2019). Higher root lignification contributes to increasing ammonium tolerance in the cell (Royo et al., 2019). In this study, we found that a high NH4+/NO3– ratio significantly increased root lignin content (Figure 10), likely because ROS directly modify the structure of the root cell wall by their involvement in lignin formation (Tsukagoshi, 2016). Indeed, ROS facilitates POD activity in the presence of NH4+ (Podgórska et al., 2021). Furthermore, we found that increasing the NH4+/NO3– ratios significantly increased the activities of PAL, 4CL, CCR, CAD, COMT, and POD (Figure 11), which are key enzymes involved in lignin biosynthesis. Moreover, the transcriptome and proteome analyses showed 598 DEGs and 70 DEPs related to phenylpropanoid biosynthesis (Supplementary Tables 9, 10). DEGs and DEPs related to PAL, 4CL, CCR, CAD, COMT, and POD were upregulated as the NH4+/NO3– ratio increased (Figure 5). In particular, we found two proteins belonging to the CAD subfamily, A0A3B6NPP6 and A0A3B6LM09, which were detected in the DEP data through phylogenetic comparison of the members of CAD among wheat genotypes (Supplementary Figure 6). Six upregulated proteins associated with POD were identified in each NH4+/NO3– ratio treatment. These results indicated that increasing the NH4+/NO3– ratios increased the expression of genes and proteins related to lignin biosynthesis, which contributed to reinforcing the root cell wall.
Excess lignin deposition in roots reduces the diameter of xylem vessels and restricts xylem extensibility, thereby limiting nutrient transport (Khandal et al., 2020). The increment in lignin deposition in the cell wall induces the formation of a root barrier that effectively reduces solute permeability and leads to reducing N transport in the roots (Ranathunge et al., 2016; Plett et al., 2020). Our results showed that increasing root NH4+/NO3– ratios increased NH4+ but reduced NO3– content in roots (Figure 8). According to the results of the transcriptome analysis, TraesCS3D02G344800 and TraesCS3A02G350800 (ammonium transporter) were upregulated under the high NH4+ conditions (Na, Nr1, and Nr2 treatments) (Figure 6). This could explain why the root had more NH4+-N content with the NH4+/NO3– ratio increment. NRT2 is a high-affinity transporter, especially under low N availability conditions (Tong et al., 2020). However, a previous study reported that TaNRT2.1 could be expressed under 5 and 10 mmol–1 N (Wang et al., 2011). A recent study also reported that Ta6A (TraesCS6A02G030800) was at or near saturation at 0.5 mM N, but it was still taking up nitrate at 10 mM nitrate (Li et al., 2021). In our study, 15 shared downregulated DEGs were identified in the N metabolism pathway, which were involved in high-affinity nitrate transport or nitrate transport, including NRT2.1 and NRT2.4 (Figure 6 and Supplementary Tables 11, 12). This explained reduced root NO3–-N content with increasing NH4+/NO3– ratios. Furthermore, a total of 12 DEPs involved in N metabolism were identified, among which A0A3B6Q9B3 was a nitrate transporter, and A0A3B6PQS3 was a ferredoxin-nitrite reductase. As a result, root N accumulation decreased significantly with increasing NH4+/NO3– ratios (Figure 7B). Correlation analysis indicated that root lignin accumulation negatively correlated with NO3– content but positively correlated with NH4+ content (Supplementary Figure 8). Based on the above, the increase in lignin deposition may decrease the N absorption ability of roots. These results demonstrate that increasing the NH4+/NO3– ratios in roots of wheat seedlings were unfavorable for root N uptake and transport due to excessive root lignification.
Conclusion
In the present study, the increment in the NH4+/NO3– ratios reduced leaf chlorophyll content, Fv/Fm, and ΦII values in the ammonium sensitive wheat, indicating that a high NH4+/NO3– ratio reduces leaf photosynthetic capacity. In addition, increasing the NH4+/NO3– ratios caused O2– generation in the roots and increased the activity of enzymes involved in lignin biosynthesis, leading to root growth inhibition. Under a high NH4+/NO3– ratio, we identified upregulated DEGs and DEPs related to glutathione metabolism and phenylpropanoid biosynthesis in roots, which may contribute to an increase of root tolerance to oxidative damage. Furthermore, increasing the NH4+/NO3– ratios downregulated DEGs and DEPs related to N metabolism and nitrate transport, leading to a decrease in root N content. These results show that high NH4+/NO3– ratios in young wheat seedlings were unfavorable for root N uptake and transport due to excess lignification of the root cell wall.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://figshare.com/articles/dataset/_/16826116.
Author Contributions
DY conceived and designed the experiments, analyzed the data, prepared figures and tables, authored and reviewed drafts of the manuscript, and approved the final draft. JZ, CB, and LL performed the experiments. ZW contributed reagents, materials, analysis tools, and authored and reviewed drafts of the manuscript. All the authors contributed to the article and approved the submitted version.
Funding
The research report in this study was supported by the National Natural Science Foundation of China (31801295), Shandong Province Natural Science Foundation (ZR2017BC106), China Postdoctoral Science Foundation funded Project (2018M632701), and China Postdoctoral Science Special Foundation (2019T120600).
Conflict of Interest
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.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.797260/full#supplementary-material
Footnotes
- ^ https://cutadapt.readthedocs.io/en/stable/
- ^ http://www.bioinformatics.babraham.ac.uk/projects/fastqc/
- ^ http://ccb.jhu.edu/software/stringtie/
- ^ https://bioc.ism.ac.jp/packages/3.0/bioc/html/DESeq2.html
References
Alencar, V. T. C. B., Lobo, A. K. M., Carvalho, F. E. L., and Silveira, J. A. G. (2019). High ammonium supply impairs photosynthetic efficiency in rice exposed to excess light. Photosynth. Res. 140, 321–335. doi: 10.1007/s11120-019-00614-z
Bai, X. Y., Dong, Y. J., Wang, Q. H., Xu, L. L., Kong, J., and Liu, S. (2015). Effects of lead and nitric oxide on photosynthesis, antioxidative ability, and mineral element content of perennial ryegrass. Biol. Plantarum. 59, 163–170. doi: 10.1007/s10535-014-0476-8
Bittsánszky, A., Pilinszky, K., Gyulai, G., and Komives, T. (2015). Overcoming ammonium toxicity. Plant Sci. 231, 184–190. doi: 10.1016/j.plantsci.2014.12.005
Cataldo, D. A., Naroon, M., Schrader, L. E., and Youngs, V. L. (1975). Rapid colorimetric determination of nitrate in plant tissue by nitration of salicylic acid. Commun. Soil. Sci. Plant. Anal. 6, 71–80. doi: 10.1080/00103627509366547
Cesarino, I. (2019). Structural features and regulation of lignin deposited upon biotic and abiotic stresses. Curr. Opin. Biotech. 56, 209–214. doi: 10.1016/j.copbio.2018.12.012
Choudhary, A., Kumar, A., and Kaur, N. (2020). ROS and oxidative burst: roots in plant development. Plant Diversity 42, 33–43. doi: 10.1016/j.pld.2019.10.002
Croft, H., Chen, J. M., Luo, X. Z., Bartlett, P., Chen, B., and Staebler, R. M. (2017). Leaf chlorophyll content as a proxy for leaf photosynthetic capacity. Global Change Biol. 23, 3513–3524. doi: 10.1111/gcb.13599
Cui, Z. L., Chen, X. P., and Zhang, F. S. (2010). Current nitrogen management status and measures to improve the intensive wheat–maize system in China. Ambio 39, 376–384. doi: 10.1007/s13280-010-0076-6
Dai, G. Z., Qiu, B. S., and Forchhammer, K. (2014). Ammonium tolerance in the cyanobacterium Synechocystis sp. strain PCC 6803 and the role of the psbA multigene family. Plant. Cell. Environ. 37, 840–851. doi: 10.1111/pce.12202
Ding, N., Wang, A. M., Zhang, X. J., Wu, Y. X., Wang, R. Y., Cui, H. H., et al. (2017). Identification and analysis of glutathione S-transferase gene family in sweet potato reveal divergent GST-mediated networks in aboveground and underground tissues in response to abiotic stresses. BMC Plant Biol. 17:225. doi: 10.1186/s12870-017-1179-z
Gill, S. S., Anjum, N. A., Hasanuzzaman, M., Trivedi, D. K., Ahmad, I., Pereira, E., et al. (2013). Glutathione and glutathione reductase: a boon in disguise for plant abiotic stress defense operations. Plant Physiol. Biochem. 70, 204–212. doi: 10.1016/j.plaphy.2013.05.032
Gojon, A. (2017). Nitrogen nutrition in plants: rapid progress and new challenges. J. Exp. Bot. 68, 2457–2462. doi: 10.1093/jxb/erx171
Gong, B., Sun, S. S., Yan, Y. Y., Jing, X., and Shi, Q. H. (2018). “Glutathione metabolism and its function in higher plants adapting to stress,” in Antioxidants and Antioxidant Enzymes in Higher Plants, eds D. K. Gupta, J. M. Palma, and F. J. Corpas (Cham: Springer International Publishing), 181–205.
Gullner, G., Komives, T., Király, L., and Schröder, P. (2018). Glutathione S-transferase enzymes in plant-pathogen interactions. Front. Plant Sci. 9:1836. doi: 10.3389/fpls.2018.01836
Halušková, L., Valentovičová, K., Huttová, J., Mistrík, J., and Tamás, L. (2009). Effect of abiotic stresses on glutathione peroxidase and glutathione S-transferase activity in barley root tips. Plant Physiol. Biochem. 47, 1069–1074. doi: 10.1016/j.plaphy.2009.08.003
Hong, J. T., Ma, X. X., Yan, Y., Zhang, X. K., and Wang, X. D. (2018). Which root traits determine nitrogen uptake by alpine plant species on the Tibetan Plateau? Plant Soil 424, 63–72. doi: 10.1007/s11104-017-3434-3
Husain, A., Muneer, M. A., Fan, W., Yin, G. F., Shen, S. Z., Wang, F., et al. (2019). Application of optimum N through different fertilizers alleviate NH4+ -N, NO3–-N and total nitrogen losses in the surface runoff and leached water and improve nitrogen use efficiency of rice crop in Erhai Lake Basin, China. Commun. Soil. Sci. Plan. 50, 716–738. doi: 10.1080/00103624.2019.1589482
Ijato, T., Porras-Murillo, R., Ganz, P., Ludewig, U., and Neuhäuser, B. (2021). Concentration-dependent physiological and transcriptional adaptations of wheat seedlings to ammonium. Physiol. Plantarum. 171, 328–342. doi: 10.1111/ppl.13113
Imran, M., Sun, X. C., Hussain, S., Ali, U., Rana, M. S., Rasul, F., et al. (2019). Molybdenum-induced effects on nitrogen metabolism enzymes and elemental profile of winter wheat (Triticum aestivum L.) under different nitrogen sources. Int. J. Mol. Sci. 20:3009. doi: 10.3390/ijms20123009
Jiang, H. W., Liu, M. J., Chen, I. C., Huang, C. H., Chao, L. Y., and Hsieh, H. L. (2010). A glutathione S-transferase regulated by light and hormones participates in the modulation of Arabidopsis seedling development. Plant Physiol. 154, 1646–1658. doi: 10.1104/pp.110.159152
Ju, X. T., Xing, G. X., Chen, X. P., Zhang, S. L., Zhang, L. J., Liu, X. J., et al. (2009). Reducing environmental risk by improving N management in intensive Chinese agricultural systems. Proc. Natl. Acad. Sci. U.S.A. 106, 3041–3046. doi: 10.1073/pnas.0813417106
Khandal, H., Singh, A. P., and Chattopadhyay, D. (2020). The MicroRNA397b -LACCASE2 module regulates root lignification under water and phosphate deficiency. Plant Physiol. 182, 1387–1403. doi: 10.1104/pp.19.00921
Kiba, T., and Krapp, A. (2016). Plant nitrogen acquisition under low availability: regulation of uptake and root architecture. Plant Cell Physiol. 57, 707–714. doi: 10.1093/pcp/pcw052
Kirschke, T., Spott, O., and Vetterlein, D. (2019). Impact of urease and nitrification inhibitor on NH4+ and NO3– dynamic in soil after urea spring application under field conditions evaluated by soil extraction and soil solution sampling. J. Plant Nutr. Soil Sci. 182, 441–450. doi: 10.1002/jpln.201800513
Kraiser, T., Gras, D. E., Gutiérrez, A. G., González, B., and Gutiérrez, R. A. (2011). A holistic view of nitrogen acquisition in plants. J. Exp. Bot. 62, 1455–1466. doi: 10.1093/jxb/erq425
Krall, J. P., and Edwards, G. E. (1992). Relationship between photosystem II activity and CO2 fixation in leaves. Physiol. Plantarum. 86, 180–187. doi: 10.1111/j.1399-3054.1992.tb01328.x
Kuhlgert, S., Austic, G., Zegarac, R., Osei-Bonsu, I., Hoh, D., Chilvers, M. I., et al. (2016). MultispeQ Beta: a tool for large-scale plant phenotyping connected to the open PhotosynQ network. Roy. Soc. Open Sci. 3:160592. doi: 10.1098/rsos.160592
Kumar, S., and Trivedi, P. (2018). Glutathione S-transferases: role in combating abiotic stresses including Aarsenic detoxification in plants. Front. Plant Sci. 9:751. doi: 10.3389/fpls.2018.00751
Li, B., Li, G. J., Kronzucker, H. J., Baluška, F., and Shi, W. M. (2014). Ammonium stress in Arabidopsis: signaling, genetic loci, and physiological targets. Trends Plant Sci. 19, 107–114. doi: 10.1016/j.tplants.2013.09.004
Li, H., Ma, Q. H., Li, H. G., Zhang, F. S., Rengel, Z., and Shen, J. B. (2014). Root morphological responses to localized nutrient supply differ among crop species with contrasting root traits. Plant Soil. 376, 151–163. doi: 10.1007/s11104-013-1965-9
Li, B. H., Li, Q., Xiong, L. M., Kronzucker, H. J., Krämer, U., and Shi, W. M. (2012). Arabidopsis plastid AMOS1/EGY1 integrates abscisic acid signaling to regulate global gene expression response to ammonium stress. Plant Physiol. 160, 2040–2051. doi: 10.1104/pp.112.206508
Li, M. J., Tian, H., and Gao, Y. J. (2021). A genome-wide analysis of NPF and NRT2 transporter gene families in bread wheat provides new insights into the distribution, function, regulation and evolution of nitrate transporters. Plant Soil 465, 47–63. doi: 10.1007/s11104-021-04927-8
Liu, G. W., Sun, A. L., Li, D. Q., Athman, A., Gilliham, M., and Liu, L. H. (2015). Molecular identification and functional analysis of a maize (Zea mays) DUR3 homolog that transports urea with high affinity. Planta 241, 861–874. doi: 10.1007/s00425-014-2219-7
Liu, H., Wang, Z. H., Yu, R., Li, F. C., Li, K. Y., Cao, H. B., et al. (2016). Optimal nitrogen input for higher efficiency and lower environmental impacts of winter wheat production in China. Agr. Ecosyst. Environ. 224, 1–11. doi: 10.1016/j.agee.2016.03.022
Liu, W., Jin, Y. Z., Li, M. M., Dong, L. J., Guo, D., Lu, C., et al. (2018). Analysis of CmCADs and three lignifying enzymes in oriental melon (‘CaiHong7’) seedlings in response to three abiotic stresses. Sci. Hortic. 237, 257–268. doi: 10.1016/j.scienta.2018.04.024
Liu, Y., Li, Y. X., Li, Y. X., Tian, Z. W., Hu, J. J., Steve, A., et al. (2021). Changes of oxidative metabolism in the roots of wheat (Triticum aestivum L.) seedlings in response to elevated ammonium concentrations. J. Integr. Agr. 20, 1216–1228. doi: 10.1016/S2095-3119(20)63216-6
Liu, Y., Sun, J. Y., Tian, Z. W., Hakeem, A., Wang, F., Jiang, D., et al. (2017). Physiological responses of wheat (Triticum aestivum L.) germination to elevated ammonium concentrations: reserve mobilization, sugar utilization, and antioxidant metabolism. Plant Growth Regul. 81, 209–220. doi: 10.1007/s10725-016-0198-3
Liu, Y., and von Wirén, N. (2017). Ammonium as a signal for physiological and morphological responses in plants. J. Exp. Bot. 68, 2581–2592. doi: 10.1093/jxb/erx086
Luo, L., Zhang, Y. L., and Xu, G. H. (2020). How does nitrogen shape plant architecture? J. Exp. Bot. 71, 4415–4427. doi: 10.1093/jxb/eraa187
Meister, R., Rajani, M. S., Ruzicka, D., and Schachtman, D. P. (2014). Challenges of modifying root traits in crops for agriculture. Trends Plant Sci. 19, 779–788. doi: 10.1016/j.tplants.2014.08.005
Miller, A. J., and Cramer, M. D. (2005). Root nitrogen acquisition and assimilation. Plant Soil 274, 1–36. doi: 10.1007/s11104-004-0965-1
Motte, H., Vanneste, S., and Beeckman, T. (2019). Molecular and environmental regulation of root development. Annu. Rev. Plant Biol. 70, 465–488. doi: 10.1146/annurev-arplant-050718-100423
Moura, J. C. M. S., Bonine, C. A. V., De Oliveira Fernandes Viana, J., Dornelas, M. C., and Mazzafera, P. (2010). Abiotic and biotic stresses and changes in the lignin content and composition in plants. J. Integr. Plant Biol. 52, 360–376. doi: 10.1111/j.1744-7909.2010.00892.x
Nianiou-Obeidat, I., Madesis, P., Kissoudis, C., Voulgari, G., Chronopoulou, E., Tsaftaris, A., et al. (2017). Plant glutathione transferase-mediated stress tolerance: functions and biotechnological applications. Plant Cell Rep. 36, 791–805. doi: 10.1007/s00299-017-2139-7
Nongmaithem, S., Ponukumatla, R., Sreelakshmi, Y., Frasse, P., Bouzayen, M., and Sharma, R. (2021). Enhanced polar auxin transport in tomato polycotyledon mutant seems to be related to glutathione levels. J. Plant Growth Regul. 40, 761–773. doi: 10.1007/s00344-020-10139-8
Plett, D. C., Holtham, L. R., Okamoto, M., and Garnett, T. P. (2018). Nitrate uptake and its regulation in relation to improving nitrogen use efficiency in cereals. Semin. Cell Dev. Biol. 74, 97–104. doi: 10.1016/j.semcdb.2017.08.027
Plett, D. C., Ranathunge, K., Melino, V. J., Kuya, N., Uga, Y., and Kronzucker, H. J. (2020). The intersection of nitrogen nutrition and water use in plants: new paths toward improved crop productivity. J. Exp. Bot. 71, 4452–4468. doi: 10.1093/jxb/eraa049
Podgórska, A., Burian, M., Dobrzyńska, K., Rasmusson, A. G., and Szal, B. (2021). Respiratory burst oxidases and apoplastic peroxidases facilitate ammonium syndrome development in Arabidopsis. Environ. Exp. Bot. 181:104279. doi: 10.1016/j.envexpbot.2020.104279
Podgórska, A., and Szal, B. (2015). “The role of reactive oxygen species under ammonium nutrition,” in Reactive Oxygen and Nitrogen Species Signaling and Communication in Plants, eds K. J. Gupta and A. U. Igamberdiev (Basel: Springer International Publishing), 133–153.
Ranathunge, K., Schreiber, L., Bi, Y. M., and Rothstein, S. J. (2016). Ammonium-induced architectural and anatomical changes with altered suberin and lignin levels significantly change water and solute permeabilities of rice (Oryza sativa L.) roots. Planta 243, 231–249. doi: 10.1007/s00425-015-2406-1
Remans, T., Nacry, P., Pervent, M., Girin, T., Tillard, P., Lepetit, M., et al. (2006). A central role for the nitrate transporter NRT2.1 in the integrated morphological and physiological responses of the root system to nitrogen limitation in Arabidopsis. Plant Physiol. 140, 909–921. doi: 10.1104/pp.105.075721
Royo, B., Esteban, R., Buezo, J., Santamaríac, E., Fernández-Irigoyenc, J., Becker, D., et al. (2019). The proteome of Medicago truncatula in response to ammonium and urea nutrition reveals the role of membrane proteins and enzymes of root lignification. Environ. Exp. Bot. 162, 168–180. doi: 10.1016/j.envexpbot.2019.02.010
Schjoerring, J. K., Husted, S., Mäck, G., and Mattsson, M. (2002). The regulation of ammonium translocation in plants. J. Exp. Bot. 53, 883–890. doi: 10.1093/jexbot/53.370.883
Setién, I., Fuertes-Mendizabal, T., González, A., Aparicio-Tejo, P. M., González-Murua, C., González-Moro, M. B., et al. (2013). High irradiance improves ammonium tolerance in wheat plants by increasing N assimilation. J. Plant. Physiol. 170, 758–771. doi: 10.1016/j.jplph.2012.12.0
Sharma, D. K., Andersen, S. B., Ottosen, C. O., and Rosenqvist, E. (2015). Wheat cultivars selected for high Fv/Fm under heat stress maintain high photosynthesis, total chlorophyll, stomatal conductance, transpiration and dry matter. Physiol. Plantarum. 153, 284–298. doi: 10.1111/ppl.12245
Sharma, D. K., Torp, A. M., Rosenqvist, E., Ottosen, C. O., and Andersen, S. B. (2017). QTLs and potential candidate genes for heat stress tolerance identified from the mapping populations specifically segregating for Fv/Fm in wheat. Front. Plant Sci. 8:1668. doi: 10.3389/fpls.2017.01668
Sharma, M., Pang, J. Y., Wen, Z. H., De Borda, A., Kim, H. S., Liu, Y. F., et al. (2021). A significant increase in rhizosheath carboxylates and greater specific root length in response to terminal drought is associated with greater relative phosphorus acquisition in chickpea. Plant Soil 460, 51–68. doi: 10.1007/s11104-020-04776-x
Shi, D. L., Zhuang, K., Chen, Y. H., Xu, F. L., Hu, Z. B., and Shen, Z. G. (2020). Effects of excess ammoniacal nitrogen (NH4+-N) on pigments, photosynthetic rates, chloroplast ultrastructure, proteomics, formation of reactive oxygen species and enzymatic activity in submerged plant Hydrilla verticillata (L.f.) Royle. Aquat. Toxicol. 226:105585. doi: 10.1016/j.aquatox.2020.105585
Swanson, S., and Gilroy, S. (2010). ROS in plant development. Physiol. Plantarum. 138, 384–392. doi: 10.1111/j.1399-3054.2009.01313.x
Tiwari, V., Patel, M. K., Chaturvedi, A. K., Mishra, A., and Jha, B. (2016). Functional characterization of the Tau class glutathione-S-transferases gene (SbGSTU) promoter of Salicornia brachiata under salinity and osmotic stress. PLoS One 11:e0148494. doi: 10.1371/journal.pone.0148494
Tong, J. F., Walk, T. C., Han, P. P., Chen, L. Y., Shen, X. J., Li, Y. S., et al. (2020). Genome-wide identification and analysis of high-affinity nitrate transporter 2 (NRT2) family genes in rapeseed (Brassica napus L.) and their responses to various stresses. BMC Plant Biol. 20:464. doi: 10.1186/s12870-020-02648-1
Tsukagoshi, H. (2016). Control of root growth and development by reactive oxygen species. Curr. Opin. Plant Biol. 29, 57–63. doi: 10.1016/j.pbi.2015.10.012
Vega-Mas, I., Cukier, C., Coleto, I., González-Murua, C., Limami, A. M., González-Moro, M. B., et al. (2019). Isotopic labelling reveals the efficient adaptation of wheat root TCA cycle flux modes to match carbon demand under ammonium nutrition. Sci. Rep. 9:8925. doi: 10.1038/s41598-019-45393-8
Wang, F., Gao, J. W., Liu, Y., Tian, Z. W., Muhammad, A., Zhang, Y. X., et al. (2016a). Higher ammonium transamination capacity can alleviate glutamate inhibition on winter wheat (Triticum aestivum L.) root growth under high ammonium stress. PLoS One 11:e0160997. doi: 10.1371/journal.pone.0160997
Wang, F., Gao, J. W., Tian, Z. W., Liu, Y., Abid, M., Jiang, D., et al. (2016b). Adaptation to rhizosphere acidification is a necessary prerequisite for wheat (Triticum aestivum L.) seedling resistance to ammonium stress. Plant. Physiol. Bioch. 108, 447–455. doi: 10.1016/j.plaphy.2016.08.011
Wang, F., Gao, J. W., Shi, S. M., He, X. H., and Dai, T. B. (2019). Impaired electron transfer accounts for the photosynthesis inhibition in wheat seedlings (Triticum aestivum L.) subjected to ammonium stress. Physiol. Plantarum. 167, 159–172. doi: 10.1111/ppl.12878
Wang, F., Gao, J. W., Yong, J. W. H., Liu, Y. F., Cao, D., and He, X. H. (2020). Glutamate over-accumulation may serve as an endogenous indicator of tricarboxylic acid (TCA) cycle suppression under NH4+ nutrition in wheat (Triticum aestivum L.) seedlings. Environ. Exp. Bot. 177:104130. doi: 10.1016/j.envexpbot.2020.104130
Wang, P., Wang, Z. L., Cai, R. G., Li, Y., Chen, X. G., and Yin, Y. P. (2011). Physiological and molecular response of wheat roots to nitrate supply in seedling stage. Agr. Sci. China 10, 695–704. doi: 10.1016/S1671-2927(11)60052-7
Wang, W. H., Köhler, B., Cao, F. Q., Liu, G. W., Gong, Y. Y., Sheng, S., et al. (2012). Rice DUR3 mediates high-affinity urea transport and plays an effective role in improvement of urea acquisition and utilization when expressed in Arabidopsis. New Phytol. 193, 432–444. doi: 10.1111/j.1469-8137.2011.03929.x
Witte, C. P. (2011). Urea metabolism in plants. Plant Sci. 180, 431–438. doi: 10.1016/j.plantsci.2010.11.010
Xu, W. Y., Tang, W. S., Wang, C. X., Ge, L. H., Sun, J. C., Qi, X., et al. (2020). SiMYB56 confers drought stress tolerance in transgenic rice by regulating lignin biosynthesis and ABA signaling pathway. Front. Plant Sci. 11:785. doi: 10.3389/fpls.2020.00785
Xu, Z. Z., Yu, Z. W., and Zhao, J. Y. (2013). Theory and application for the promotion of wheat production in China: past, present and future. J. Sci. Food Agr. 93, 2339–2350. doi: 10.1002/jsfa.6098
Yan, J., Su, P. S., Li, W., Xiao, G. L., Zhao, Y., Ma, X., et al. (2019). Genome-wide and evolutionary analysis of the class III peroxidase gene family in wheat and Aegilops tauschii reveals that some members are involved in stress responses. BMC Genomics 20:666. doi: 10.1186/s12864-019-6006-5
Yang, D. Q., Dong, W. H., Luo, Y. L., Song, W. T., Cai, T., Li, Y., et al. (2017). Effects of exogenous 6-BA on photosynthetic characteristics and endogenous hormone content in wheat leaves under two nitrogen application levels at seedling stage. Sci. Agr. Sin. 50, 3871–3884. doi: 10.3864/j.issn.0578-1752.2017.20.004
Yuan, S., and Peng, S. B. (2017). Exploring the trends in nitrogen input and nitrogen use efficiency for agricultural sustainability. Sustainability 9:1905. doi: 10.3390/su9101905
Zechmann, B. (2014). Compartment-specific importance of glutathione during abiotic and biotic stress. Front. Plant Sci. 5:566. doi: 10.3389/fpls.2014.00566
Zheng, M. J., Chen, J., Shi, Y. H., Li, Y. X., Yin, Y. P., Yang, D. Q., et al. (2017). Manipulation of lignin metabolism by plant densities and its relationship with lodging resistance in wheat. Sci. Rep. 7:41805. doi: 10.1038/srep41805
Zhou, H., Zhou, Y., Zhai, F. C., Wu, T., Xie, Y. J., Xu, G. H., et al. (2021). Rice seedlings grown under high ammonia do not show enhanced defence responses. Food Energy Secur. e331. doi: 10.1002/fes3.331
Zhu, Z. L., and Chen, D. L. (2002). Nitrogen fertilizer use in China- Contributions to food production, impacts on the environment and best management strategies. Nutr. Cycl. Agroecosys. 63, 117–127. doi: 10.1023/A:1021107026067
Keywords: root length, ROS, glutathione transferase, lignin, nitrogen transport, transcriptome, proteomic
Citation: Yang D, Zhao J, Bi C, Li L and Wang Z (2022) Transcriptome and Proteomics Analysis of Wheat Seedling Roots Reveals That Increasing NH4+/NO3– Ratio Induced Root Lignification and Reduced Nitrogen Utilization. Front. Plant Sci. 12:797260. doi: 10.3389/fpls.2021.797260
Received: 18 October 2021; Accepted: 15 December 2021;
Published: 13 January 2022.
Edited by:
Min Huang, Hunan Agricultural University, ChinaReviewed by:
Daniel Marino, University of the Basque Country, SpainIdoia Ariz, Public University of Navarre, Spain
Copyright © 2022 Yang, Zhao, Bi, Li and Wang. 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) and the copyright owner(s) 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: Dongqing Yang, chengyang2364@126.com; Zhenlin Wang, zlwang@sdau.edu.cn