Next Article in Journal
Chemical Composition and Biological Activities of Tunisian Ziziphus lotus Extracts: Evaluation of Drying Effect, Solvent Extraction, and Extracted Plant Parts
Next Article in Special Issue
Influence of Nitrogen Nutrition on Fatty Acids in Oilseed Rape (Brassica napus L.)
Previous Article in Journal
Plant Response to Mechanically-Induced Stress: A Case Study on Specialized Metabolites of Leafy Vegetables
Previous Article in Special Issue
Nitrogen Assimilation Related Genes in Brassicanapus: Systematic Characterization and Expression Analysis Identified Hub Genes in Multiple Nutrient Stress Responses
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Genes Associated with Nitrogen Stress Responses in Apple Leaves

Apple Research Institute, National Institute of Horticultural and Herbal Science, Rural Development Administration, 107, Soboangye-ro, Gunwi 39000, Korea
*
Author to whom correspondence should be addressed.
These authors are co-first authors of this work.
Plants 2021, 10(12), 2649; https://doi.org/10.3390/plants10122649
Submission received: 6 November 2021 / Revised: 30 November 2021 / Accepted: 30 November 2021 / Published: 2 December 2021
(This article belongs to the Special Issue Plant Nutrition Volume II)

Abstract

:
Nitrogen (N) is an essential macronutrient that regulates diverse physiological processes for plant survival and development. In apple orchards, inappropriate N conditions can cause imbalanced growth and subsequent physiological disorders in trees. In order to investigate the molecular basis underlying the physiological signals for N stress responses, we examined the metabolic signals responsive to contrasting N stress conditions (deficient/excessive) in apple leaves using transcriptome approaches. The clustering of differentially expressed genes (DEGs) showed the expression dynamics of genes associated with each N stress group. Functional analyses of gene ontology and pathway enrichments revealed the potential candidates of metabolic signals responsible for N-deficient/excessive stress responses. The functional interactions of DEGs in each cluster were further explored by protein–protein interaction network analysis. Our results provided a comprehensive insight into molecular signals responsive to N stress conditions, and will be useful in future research to enhance the nutrition tolerance of tree crops.

1. Introduction

Nitrogen (N) is a fundamental macronutrient that plays an essential role in regulating a wide range of physiological metabolisms in plants for survival and development. Plants acquire N from the soil as inorganic forms via membrane transporters, such as ammonium transporters (AMTs) and nitrate transporters (NRTs) (or nitrate transporter/peptide transporters; NPFs) [1,2]. After uptake, N sources are utilized in plants, followed by N assimilation, which involves the reduction of nitrate to ammonium and conversion into amino acids [3,4]. The N availability and metabolic signaling of N assimilates are not only coordinated and balanced with internal carbon (C) metabolism during the developmental stages as a source–sink relationship [5,6,7], but are also primarily affected by the external conditions of N supply, including starvation/overabundance [8,9,10,11].
For many horticultural crops such as apple, the nutritional status of N affects many physiological processes. While insufficient N fertilization can cause decreased tree growth, reduction of sugar metabolism, or poor bud development [12,13,14], while the over-fertilization of N can cause a reduction in floral bud differentiation, alternate fruit-bearing, decrease in fruit quality, and elongation of vegetative growth [15,16,17,18,19].
Among the main horticultural crops, apple (Malus domestica Borkh.) occupies a large portion of the fruit industry in the world. In apple orchards, proper N fertilization is a critical component of farmers’ tree management, and during the growing season, the stress responses of apple trees induced by various N fertilization can result in unbalanced growth and physiological disorders. Global transcriptome data from previous studies [20,21] and recent studies have suggested examining the different responses to N tolerance in apple tissues [22,23]. However, little is known about the underlying molecular mechanisms of responsive metabolic signals associated with various N stress conditions in apple.
In order to investigate this further, we investigated molecular signals in response to different N stress (deficient/excessive) conditions using transcriptome analysis of apple leaves because the leaf tissue reflects the overall N nutritional status of apple trees during the growing season [24]. We identified the clusters of 2212 differentially expressed genes (DEGs), showing distinctive expression dynamics unique to each N stress. Functional enrichment and gene interaction network analyses indicated potential candidates for metabolic signals differentially associated with N stress conditions. Our results provided useful information about metabolic signals and the genes involved in N stress responses in apple.

2. Results and Discussion

2.1. Transcriptome Profiling and the Leaf N Content

In this study, we investigated the transcriptome data of ‘Fuji’/M9 apple leaves from four different N fertilization groups consisting of two N stress conditions and the control: N-deficient (0×), control (1×), and N-excessive (2×, 4×) (Figure 1a). We observed significant differences (Duncan’s multiple range test, p < 0.05) in foliar N content among the four groups from the leaves collected from one-year-old shoots (Figure 1b). The N-excessive 4× group exhibited the highest foliar N content, followed by the 2×, control, and N-deficient 0× groups in a sequential manner, confirming that N fertilization treatments correlated with the divergent leaf N content in Fuji/M9 apple trees.
To study the metabolic signals associated with N stress responses in apple, we constructed RNA-seq profiles of leaf tissues in four N fertilization groups. A total of 647 million reads were obtained for 12 libraries across all groups and used for mapping against 45,116 protein-coding genes from the GDDH13v1.1 apple reference genome [25] (Tables S1 and S2). In the principal component analysis (PCA) on the data set of 12 libraries, PC1 explained 55.6% of the total variance, separating the difference among three N fertilization groups: N-deficient, control, and N-excessive, followed by the separation between 2× and 4× samples within the N-excessive group in PC2, showing 15.4% of the variance (Figure 1c). We did not observe a clear separation between the two groups of N-excessive (2×, 4×) in the first two PCs. However, this was expected, given that the leaf samples were collected from only the same Fuji genotype. We then identified a total of 2212 DEGs by comparing samples in N stress groups (0×, 2×, 4×) against the control (1×) group with a threshold of |fold change| ≥ 1.5 and FDR ≤ 0.05 (Figure 1d). Details of the DEGs in the N stress groups are presented in Table S3. The highest number of DEGs was detected in the N-deficient condition. The expression of identified DEGs was validated by qRT-PCR on a set of eight selected genes, where the relative expression pattern of qRT-PCR in the four N fertilization groups was compared with reads per kilobase per million (RPKM) data from RNA-seq (Figure S1; Tables S4 and S5). Together these results suggest that our transcriptome data reflect the samples collected from the contrasting N fertilization groups.

2.2. Analysis of Marker Genes Known to Be Associated with N Uptake, Assimilation, and Transport Pathways in Apple

In order to understand the dynamics of our transcriptome data associated with N stress responses in apple, we investigated the expression pattern of genes known to be involved in N uptake, transport, assimilation, and signaling pathways conserved in plant [1,2,26,27] In the ammonium uptake category, we observed that MdAMT1.1 was notably upregulated in the N-deficient (0×) group. In contrast, its expression was low in control (1×) and N-excessive (2×) (Figure 2). Similarly, among the genes involved in nitrate uptake and transport, most genes of MdNRTs and MdPTRs were upregulated in the N-deficient group compared to others, including one MdPTR5 gene that was also detected in our 2212 DEG set. On the other hand, the expression level of MdNRT1.2 was rather maintained higher in the N-excessive (2×) group than in the other three groups, and MdNitr was downregulated in both the N-deficient and N-excessive groups, while the control (1×) showed upregulation of its expression. After the uptake of N from soil as ammonium or nitrate, plants utilize N through assimilation processes [3]. In our transcriptome data, there were clear differences in the expression patterns of N assimilation. Although MdNR and MdNiR, both of which are known to be involved in the primary nitrate assimilation pathway [28,29], were upregulated in the N-deficient group and downregulated in N-excessive groups, the expression of MdGS1 and MdGOGAT genes, which are both involved in glutamate conversion from ammonia assimilation [30], showed a contrasting trend, contrary to the expression pattern of genes mediating nitrate assimilation with downregulation in N-deficient and upregulation in N-excessive groups.
After uptake into the cytoplasm, the process of nitrate signaling is mediated by Ca2+ [31]. As a primary nitrate response, Ca2+-CPK-NLP signaling is activated [32]. Our RNA-seq data showed a discrepancy among N stress groups. Although the expression pattern of MdNLP6 and MdNLP7 differed, MdNLP6 rather showed a similar upregulation in the N-deficient group as that of the main CPK-NLP complex component MdCPK10. There were many TFs enriched instead of regulating Ca2+-mediated nitrate signaling that was found to be upregulated in the N-excessive group. In particular, these included the upregulation of MdANR1, a member of MADS box TFs acting as regulators of lateral root elongation [33], and MdSPL9, a TF that was proposed as a potential hub regulator in the plant N regulatory network [34] in N-excessive 2× followed by 4×, compared with the two remaining groups. Additionally, the expression of MdNAC4, another nitrate-responsive TF that also participates in nitrate response [35], was upregulated in the N-excessive group (highest in 4× followed by 2×). As downstream genes of the nitrate response, LBD (LATERAL ORGAN BOUNDARY DOMAIN) family MdLBD37-39, TFs that are known to negatively regulate anthocyanin biosynthesis [36] were enriched. The expression of MdLBD37 and MdLBD38 was upregulated in the N-excessive group, whereas the upregulation of MdLBD39 was found in the N-deficient group. These results may imply that the response of N uptake is relatively triggered more under N starvation stress conditions, whereas the responses of N assimilation, signaling, and regulatory mechanisms are more likely to be activated under excessive N stress conditions.

2.3. Gene Expression Dynamics in N Fertilization Groups Associated with Different N Stress Responses

To uncover the expression patterns of genes responsive to N stress, we examined K-means clustering of 2212 DEG datasets, and five clusters were obtained (C1–C5, details are presented in Figure S2 and Table S6). Among the gene clusters, four clusters revealed the expression dynamics unique to each N stress condition (C1 and C4 specific to N-deficient; C2 and C5 specific to N-excessive), whereas cluster C3 was positively correlated with leaf N gradient (Figure 3a). In each cluster, we investigated the enriched GO terms (FDR < 0.05) (a detailed analysis is presented in Table S7). In cluster C1 (up in N-deficient), enrichment of the secondary metabolic process (GO: 0019748) and glycolipid metabolic process (GO: 0006664) were discovered (Figure 3b). In addition, the GO term response to stimulus (GO: 0050896) was also enriched, which included responses to osmotic stress and nutrients. Notably, phosphate starvation (GO: 0016036) was enriched in the form of clusters. This result may indicate that N-deficient nutritional conditions could trigger the defection of phosphorus (P) utilization in apple leaves stressed from N starvation, given that the uptake of N/P is known to be coordinated with each other mediated by integrative signals such as NIGT1 and SPX4 [37,38,39]. Compared to C1, cluster C2 (down in N-excessive) showed distinct GO enrichments different from those of C1. The main difference between C1 and C2 was the expression pattern of genes in control (1×). In C2, the main enrichments of GO terms were regulation of transcription (GO: 0006355) and detoxification (GO: 0098754) (Figure 3c). Furthermore, we found that the enrichment of auxin polar transport (GO: 0009926) was expected, considering that high nitrate conditions are known to repress auxin transport modulated by nitrate transporters [40]. Similarly, the enriched term floral whorl development (GO: 0048438) was discovered, suggesting that the process of floral organ development is downregulated under N-excessive stress conditions, which was partially supported by a previous study showing that NRT1.1 controlled flowering time through the interaction with FLC signaling in plants [41].
In cluster C3 (positive N correlation), various GO terms were discovered, especially those related to stress responses, including abiotic stresses: water (GO: 0009415), cold (GO: 0009409), and chemicals: sucrose (GO: 0009744) and jasmonic acid (GO: 0009753) (Figure 3d). As previously reported, N fertilization could broadly affect stress responses, such as drought and defense in plants [42,43,44]. Since 469 DEGs in C3 represented the correlation to leaf N gradient in expression levels, we further investigated their correlation patterns in expression through a Pearson’s correlation analysis of gene expression with leaf N content across each N fertilization group (Table S8). Of the 469 DEGs, 347 genes showed a high correlation in their expression patterns with leaf N content (r ≥ 0.9).
On the other hand, no significant enriched terms (biological processes) were found in C4 (down in N-deficient). In cluster C5 (up in N-excessive), the main enriched GO terms contained photosynthesis (GO: 0015979) and response to carbohydrates (GO: 0009743) (Figure 3e), indicating that the process of generation and consumption of the energy cycle was strongly activated when N was over-fertilized. In terms of the relationship between photosynthesis and N response, the expression of the nitrate reductase (NR) gene is known to decrease when enough sucrose is produced by photosynthesis [45]. Our RNA-seq data indicated that a couple of NR genes, MdNR and MdNiR, were downregulated in the N-excessive group compared to the others (Figure 2), which may also support the enrichment result in C5.

2.4. Functional Annotation of DEGs in Five Clusters Associated with N Fertilization Stresses

We then investigated the five main categories of functionally annotated DEGs in each cluster, including photosynthesis, cell wall, secondary metabolism, hormone metabolism, and transport (Figure 4, Table S9).
In the category of photosynthesis, C5 (up in N-excessive) showed the largest enrichment among the five clusters, with more than 30 relevant genes associated with photosystems (I, II), electron carrier, Calvin cycle, and photorespiration. In comparison, the remaining four clusters showed enrichment of fewer than 10 genes (Figure 4a), indicating the functional importance of photosynthesis as an upregulated response in the N-excessive group. This result was concurrent with the discovery of the photosynthesis term (GO: 0015979) found in C5 (Figure 3e).
In the cell wall category, C3 (positive N correlation) and C5 revealed abundant enrichment relative to the others, showing cell wall modification and degradation-related genes occupying more than half of the entire gene count in each cluster (Figure 4b).
As for secondary metabolism, C1 (up in N-deficient) showed the largest enrichment of genes, most of which were involved in flavonoid and lignin biosynthesis (Figure 4c), which was concurrent with the GO enrichment of secondary metabolism (GO: 0019748) found in C1 (Figure 3b). Concerning flavonoids, we demonstrated that the expression of nitrate-responsive TFs MdLBD37-38 was downregulated in N-deficient but upregulated in N-excessive groups (Figure 2), which repress anthocyanin biosynthesis [36]. Among the enriched genes related to lignin, MdPAL and Md4CL were found, the expression of which is also known to be induced under N-deficient conditions in tobacco as previously shown in [46] (Table S9).
In terms of hormone metabolism, we found a global enrichment of various hormones in most of the DEG clusters (Figure 4d). C1 showed the largest number of gene enrichments among the five clusters, while the remaining clusters showed less than 20 enriched genes. Gene ontology analysis revealed the enrichment of auxin in C2 and jasmonic acid in C3 (Figure 3c,d). The functional DEG annotation of C2 and C3 also showed a large proportion of auxin and jasmonic acid in each cluster, respectively.
In the transport group, various transporter-related genes were detected in the five clusters (Figure 4e). The DEG cluster exhibiting the largest enrichment of genes was C1, separated by 13 different target molecules in transport, but when compared to C4 (the other cluster representing N-deficient unique expression pattern contrary to C1), nitrate, potassium, and cation transporters were found to be unique to C1. On the other hand, comparing the list of enriched transporters between the pair of two contrasting N-excessive clusters (C2 and C5), ammonium, amino acids, cations, phosphate, and sulfate were C2-specific, while metabolites were unique to C5. These data suggested the possibility that global metabolic signals are differentially induced in response to nitrogen fertilization stress conditions in apple.

2.5. Identification of Candidate Metabolic Signals for N Stress Responses in Apple

We compared the KEGG pathway enrichments in each DEG cluster to further uncover the metabolic signals mediating N stress responses (Table S10). Among these, some metabolic pathways were commonly enriched in N-deficient/excessive clusters (e.g., nitrogen metabolism, carbon fixation, and α-linolenic acid metabolism), whereas others were specific to each DEG cluster (e.g., flavonoid biosynthesis, fatty acid elongation, and photosynthesis).
Among the N-deficient stress-associating clusters C1 (up in N-deficient) and C4 (down in N-deficient), flavonoid biosynthesis (ath00941) and carbon fixation (ath00710) were mainly enriched in C1, while C4 showed pathway enrichment of nitrogen metabolism (ath00910), α-linolenic acid metabolism (ath00592), and plant–pathogen interaction (ath04626). In particular, most of the enriched genes in the carbon fixation pathway were identified as the PPC phosphoenol pyruvate carboxylase (PEPC)-coding gene family (MdPPC1, MdPPC3, MdPPC5) (Table S11).
Among the enriched pathways found in N-excessive responsive DEGs C2 (down in N-excessive) and C5 (up in N-excessive), carbon fixation (ath00710) was enriched in both C2 and C5, and enrichment was also detected in C1. The functions of enriched genes within this category differed between C1 and C2. While MdPPCs were unique to C1, a couple of ribulose biphosphate carboxylase/oxygenase (Rubisco) small subunit (RBCS) family genes (MdRBCS1A and MdRBCS3B) were found to be specific to C2. Within the carbon fixation pathway, the expression of MdRBCS1A and MdRBCS3B genes is known to regulate the content of rubisco [47], the enzyme activity of which was reported to be reduced when its content was sufficient for photosynthesis under high N conditions in tree crops [48]. Considering these results, the enrichment of RBCS genes in N-excessive groups could be the result of the regulatory response triggered by the over-N condition to balance the excessive photosynthesis and N storage. On the other hand, apart from its anaplerotic function in carbon fixation, PEPC is also known to play a role in mediating the balance of carbon (C) and N metabolism by affecting both sugar accumulation and ammonium assimilation [49], suggesting that the enrichment of PPCs in the N-deficient group may reflect the physiological status of imbalanced C/N metabolism under N starvation stress condition. Furthermore, C5 showed enrichment of two photosynthesis-related pathways (ath00195, ath00196), ribosome (ath03010), N metabolism (ath00910), and hormone signal transduction (ath04075). Within the enrichment of photosynthesis-related pathways in C5, genes encoding photosystem I and II subunits (MdPnsLs, MdPSAs, MdPSBs) and antenna proteins (MdLHCAs, MdLHCBs) were upregulated in the N-excessive group, implying the maintenance of leaf photosynthesis under high N fertilization conditions.
Notably, although all clusters showed enrichment of hormone metabolism in the functional annotation of DEGs, most of which were biosynthesis-related (Figure 4e; Table S9), KEGG pathway analysis revealed that C5 exhibited a significant enrichment of hormone signal transduction (Tables S10 and S11). Auxin, cytokinin, and abscisic acid (ABA) hormones are involved in this signal transduction pathway, all of which are closely related to N signaling [50,51]. In the auxin signaling pathway, MdTIR1, a major auxin receptor that activates auxin responses [52] and a couple of Aux/IAA family genes (MdIAA7 and MdIAA16) were discovered. The enrichment of cytokinin-signaling-related genes (MdAARs, MdAHPs) in N-excessive groups could be functionally linked to upregulation of MdGS1 and MdGOGAT in our transcriptome data (Figure 2), provided that glutamate metabolism of N assimilation is known to be activated by cytokinin [53]. From the pathway enrichment of ABA signaling, we found that two genes were enriched in N-excessive groups, the main regulators of ABA signaling: a negative regulator MdPP2CA as the protein phosphatase 2C family, and a positive regulator MdOST1 as the SnRK2 family. In addition to this, in our RNA-seq data, we demonstrated the downregulation of MdCBL1 and MdCIPK23 and the upregulation of MdABI2 in N-excessive groups from the expression profile of genes involved in nitrate uptake and transport (Figure 2), suggesting that N uptake was suppressed under high N fertilization conditions. A recent study demonstrated that SnRK2-dependent ABA signaling repressed nitrate uptake [54], and the enrichment of MdOST1 in the ABA signaling pathway may be concurrent with these expression data.
In cluster C3 (positive N correlation), we focused on the highly correlated 347 DEGs to analyze the KEGG pathway enrichment, since the expression pattern of these genes is expected to reflect the physiological changes according to the gradient of N fertilization conditions (Table S8). Interestingly, among the enriched metabolic pathways of highly correlated C3 DEGs, a couple of GAD family genes MdGAD1 and MdGAD4, were found in the butanoate metabolism pathway that catalyzes the decarboxylation of glutamine to generate γ-aminobutyrate (GABA), implying that the conversion of glutamine to GABA might be functionally correlated with the external N gradient. GABA plays a role in balancing carbon and nitrogen metabolism [55,56]. Considering that GABA can also negatively control nitrate uptake under high N conditions [57], this butanoate pathway enriched in the highly correlated C3 DEG cluster might be linked to the downregulation of many nitrate-uptake-related genes, including the MdPTR family, in N-excessive groups compared to their upregulation in the N-deficient group (Figure 2).
The metabolic candidate signals in response to N-deficient/excessive stress conditions are illustrated in Figure 5, integrated with the functional category of metabolic pathways. Our results provided a comprehensive overview of molecular characteristics for a detailed understanding of the underlying N nutritional signals in apple, implying the candidate genes to potentially target for further research to enhance N nutrition in apple.

2.6. PPI Network Analysis Provides the Additional Clue for Functionally Important Metabolic Signals Responsive to N Stresses

To examine another layer of evidence, we analyzed the PPI networks of the five DEG clusters (Figure 6). The interactions of gene pairs were compared in each DEG cluster to build the PPI networks at a high confidence level (combined score ≥ 0.7); among the five clusters, C5 showed the largest interaction pairs, while the rest of the networks in the four clusters contained much fewer interacting nodes (Table S12).
The PPI network of C1 (up in N-deficient) showed the interacting modules of genes involved in secondary metabolism, the uptake of N/P, auxin signaling, and heat stress (Figure 6a), which was also considerably concurrent with the GO enrichment results of the C1 cluster containing functional terms with secondary metabolic process (GO: 0019748), response to stimulus (GO: 0050896), and P starvation (GO: 0016036) (Figure 2; Table S7). Among the interacting modules in the PPI network of C2 (down in N-excessive), gene interactions mediating auxin transport, transcription factor, detoxification, and organ development were discovered (Figure 6b), in which metabolism was linked to GO term enrichment, representing detoxification (GO: 0098754), auxin transport (GO: 0060918), regulation of transcription (GO: 0006355), and floral whorl development (GO: 0048439) (Figure 2; Table S7).
The network of C3 DEGs (positive N correlation) contained the interacting modules involved in N assimilation, α-linolenic acid metabolism, and flavonoid biosynthesis (Figure 6c). Among these, α-linolenic acid metabolism was also discovered in the list of KEGG pathway enrichments from highly correlated C3 DEGs, including the relevant genes MdHPL1 and MdLOX2,3 involved in this pathway (Tables S10 and S11).
Although there were no significant GO terms found in C4 (down in N-deficient) (Figure 2), we found some interacting nodes in the C4 PPI network, which are expected to potentially function in organs of chloroplasts and peroxisomes, starch homeostasis, and N assimilation (Figure 6d). The PPI network of C5 (up in N-excessive) indicated the two main highly interacting clusters of gene modules that function in photosynthesis and ribosomal components (Figure 6e), suggesting the functional importance of the enrichment of metabolic pathways, including photosynthesis, the photosynthesis antenna protein, and the ribosome detected in C5 (Tables S10 and S11).

3. Materials and Methods

3.1. Plant Materials and Sample Collection

Fuji apple grafted onto M9 rootstocks were planted randomly onto 55 cm × 55 cm × 60 cm pots and were grown under natural field conditions in the experimental plot of the Apple Research Institute (36.3° N, 128.5° E), Gunwi, Korea. Apple trees were then divided into three groups according to levels of N fertilization: N-deficient (0×; without fertilization), control (1×; 7 g of urea), and N-excessive (2×, 14 g of urea; 4×, 28 g of urea). To study N-stress-responsive molecular signals in apple leaves, we collected foliar tissues between the third and tenth midshoot positions of new one-year-old shoots from each group after 6 months (Figure 1a) with four biological replicates.

3.2. Leaf N Content

Leaf samples were dried for 3 days and then ground. Contents from each ground leaf were extracted with a sulfuric acid solution with catalyst accelerators (K2SO4 + selenium (1000 Kjeltabs S/3.5, Foss, Hägersten, Sweden), heated at 420 °C for 1 h, and then cooled down for Kjeldahl digestion [58,59]. Leaf N content was measured using a Kjeldahl 8400 nitrogen analyzer (Foss, Hägersten, Sweden).

3.3. RNA Extraction and Sequencing

Total RNA was extracted using a modified CTAB method [59]. RNA quality was measured using an Agilent 2100 bioanalyzer (Palo Alto, CA, USA), and samples with a RIN of 7 or above were sent to CnK Genomics (Korea) for sequencing. There were three biological replicates in each N group used for sequencing. RNA-seq libraries were prepared using the MGIEasy RNA directional library preparation kit (MGI, Shenzhen, China), and a total of 12 libraries were constructed using the MGISEQ-2000 platform (MGI, Shenzhen, China) with 2 × 100 bp paired-end (PE) reads.

3.4. Read Processing, Mapping, and Expression Data Analysis

Raw reads were trimmed to remove low-quality and short sequences. Only PE reads with a quality score of <0.001 (equivalent to Phred Q30), ≤2 ambiguous nucleotides, and a length of ≥75 bp were retained for further analyses. Trimmed reads were aligned to the apple reference genome GDDH13v1.1 [25] in CLC Genomics Workbench 20 software [60] using the parameters of both the minimum length fraction and minimum similarity to 0.8. This genome version contained 51,000 genes, including 45,116 protein-coding genes and noncoding RNAs. Raw read counts were normalized to RPKM for each gene. Differentially expressed genes (DEGs) were determined by comparing three groups of N stress conditions (deficient: 0×; excessive: 2× and 4×) against control (1×) with a |fold change| ≥1.5, and a false discovery rate (FDR) value of ≤0.05. DEG clusters were analyzed using the K-means method in MeV software ver4.9.0 [61]. GO term enrichment was investigated using the agriGO enrichment tool [62] and clustered using the REVIGO visualization tool [63]. Pathway enrichment of DEGs was analyzed in each cluster using the DAVID bioinformatics resources database [64], based on the pairwise sequence comparison of GDDH13v1.1 transcript with NCBI Arabidopsis (TAIR10) homologs with a threshold of E-value < 1 × 10−6 provided from the GDR genome database [65].

3.5. Protein–Protein Interaction (PPI) Network Analysis

In order to build the PPI network of DEGs in each cluster, Arabidopsis homolog gene IDs were used to predict the PPI network using the STRING database [66]. The interaction score of each gene pair was obtained by combining the main sources (coexpression, lab experiments, text mining, and annotation in databases), and the PPI network was constructed with a high confidence level (combined interaction score ≥ 0.7).

3.6. Quantitative Reverse Transcription PCR (qRT-PCR)

First-strand complementary DNA was synthesized using 1.0 μg of total RNA, oligo dT primer, and Transcriptor Reverse Transcriptase (Roche, Penzberg, Germany). The qRT-PCR analysis was performed using LightCycler 480 SYBR Green I Master mix (Roche, Germany) on a Roche 480 LightCycler® (Basel, Switzerland). Complementary DNA (1:20 dilution) was used as a template (5 μL) in a 20 μL reaction volume. For each sample type, there were 4–6 technical replicates. PCR cycles were as follows: initial denaturation at 95 °C for 5 min, followed by 45 cycles of 95 °C for 10 s, 65 °C for 15 s, and 72 °C for 12 s, and a final melt curve analysis to determine the amplification of a single product. Primers were designed using Primer-Blast (http://www.ncbi.nlm.nih.gov/tools/Primer-Blast, accessed on 10 May 2021) to span an intron, if possible, with a 100–150 bp product size. MDP0000336547 (SGF29 Tudor-like domain) was selected as the reference gene [67]. Primer efficiencies and relative expression levels of targets were calculated using the Roche 480 Light Cycler software E-Method [68].

3.7. Quantification and Statistical Analyses

The heatmap of genes showing hierarchical clustering was analyzed using the Pheatmap R package (https://cran.r-project.org/web/packages/pheatmap, accessed on 20 July 2021). Functional annotations were generated using the online Mercator v3.6 tool [69]. Other statistical analyses were performed using Microsoft Excel 2016 and the R program (http://www.R-project.org, accessed on 9 August 2021).

4. Conclusions

Recently, studies have targeted apple tissues to explore global N tolerance using RNA-seq analysis [22,23]; however, the detailed metabolic signals associated with N stress responses have not been investigated well. Our transcriptome data of apple leaves collected under contrasting N stress conditions provided valuable information about the expression dynamics of genes associated with each N-deficient/excessive stress response. The functional annotations of DEG clusters and pathway enrichment analyses demonstrated the potential list of regulatory genes that may contribute to the differential responses specific to contrasting N stress conditions, following further validation with the gene interaction network approach. In order to better understand the physiological responses to N stress, we should apply it to an orchard of fruit crops with improved fertilization or tolerance to external nutrient fluctuations. We suggest a collective profile of metabolic signals responsive to contrasting N stresses, including those related to photosynthesis (PSAF, LHCB, RBCS1A), hormone signaling (TIR1, ARR, OST) unique to N-excessive signals, and those related to carbon fixation (PPC) and secondary metabolism (C4H, PAL, 4CL) unique to N-deficient signals. Our results provided a novel insight into signals associated with N stress responses, and will be useful in future research to enhance N tolerance and optimize fertilization in apple.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/plants10122649/s1, Figure S1. The qPCR results for RNA-seq validation. RPKM data were compared with relative expression values in qPCR, and the correlation coefficient r was obtained using a set of 8 selected genes. Figure S2. Figure of merit (FOM) analysis of the 2212 DEG datasets. Table S1. Summary statistics of RNA-seq data. Table S2. Expression level (RPKM) of 45,116 apple genes in our transcriptome data. Table S3. DEGs from each N fertilization group and the expression level (RPKM) of 2212 pooled DEGs across all groups. Table S4. Primers used for qPCR validation. Table S5. The qRT-PCR results compared against RNA-seq data. Table S6. The five expression clusters of the 2212 DEG datasets (C1–C5). Table S7. GO enrichment analysis of five expression clusters (C1–C5). Table S8. Pearson’s correlation analysis of cluster 3 DEGs that positively correlated with leaf N content. Table S9. Functional annotation of DEGs in five expression clusters. Table S10. KEGG pathway enrichment of five expression clusters (C1–C5). Table S11. Annotation of genes found in the enriched metabolic pathways from five expression clusters (C1–C5). Table S12. Scoring of protein–protein interaction (PPI) of DEGs in five expression clusters (C1–C5).

Author Contributions

Y.L. and V.G.D. contributed equally as first authors; Y.L. conceived the idea and performed the tissue preparation and RNA-seq library construction; Y.L. and S.K. contributed to the DEG analysis; V.G.D. contributed to qRT-PCR validation; S.K. and H.K. contributed to data curation; Y.L. and V.G.D. contributed to managing experiments and analyses; Y.L., V.G.D., S.K. and H.K. contributed to manuscript writing. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the Rural Development Administration research program (PJ01424101), Apple Research Institute, National Institute of Horticultural and Herbal Sciences, Rural Development Administration, Korea.

Data Availability Statement

All sequencing data were deposited in the National Center for Biotechnology Information Sequence Read Archive database, bearing the BioProject ID PRJNA773514.

Acknowledgments

This work was supported by the Rural Development Administration research program, Apple Research Institute, National Institute of Horticultural and Herbal Sciences, Rural Development Administration, Korea.

Conflicts of Interest

The authors declare that they have no competing interest.

References

  1. Lanquar, V.; Loqué, D.; Hörmann, F.; Yuan, L.; Bohner, A.; Engelsberger, W.R.; Lalonde, S.; Schulze, W.X.; von Wirén, N.; Frommer, W.B. Feedback inhibition of ammonium uptake by a phospho-Dependent allosteric mechanism in Arabidopsis. Plant Cell 2009, 21, 3610–3622. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Leran, S.; Varala, K.; Boyer, J.C.; Chiurazzi, M.; Crawford, N.; Daniel-Vedele, F.; David, L.; Dickstein, R.; Fernandez, E.; Forde, B.; et al. A unified nomenclature of NITRATE TRANSPORTER 1/PEPTIDE TRANSPORTER family members in plants. Trends Plant. Sci 2014, 19, 5–9. [Google Scholar] [CrossRef] [PubMed]
  3. Kwang, W.J. International Review of Cell and Molecular Biology; Academic Press: Cambridge, MA, USA, 2014; Volume 310, pp. 1–37. [Google Scholar] [CrossRef]
  4. Xu, G.; Fan, X.; Miller, A.J. Plant nitrogen assimilation and use efficiency. Annu. Rev. Plant Biol. 2012, 63, 153–182. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Nunes-Nesi, A.; Fernie, A.R.; Stitt, M. Metabolic and signaling aspects underpinning the regulation of plant carbon nitrogen interactions. Mol. Plant 2010, 3, 973–996. [Google Scholar] [CrossRef]
  6. Huarancca Reyes, T.; Scartazza, A.; Pompeiano, A.; Ciurli, A.; Lu, Y.; Guglielminetti, L.; Yamaguchi, J. Nitrate reductase modulation in response to changes in C/N balance and nitrogen source in Arabidopsis. Plant. Cell Physiol. 2018, 59, 1248–1254. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Shaw, R.; Cheung, C. A dynamic multi-Tissue flux balance model captures carbon and nitrogen metabolism and optimal resource partitioning during Arabidopsis growth. Front. Plant Sci. 2018, 9, 884. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Bi, Y.-M.; Wang, R.-L.; Zhu, T.; Rothstein, S.J. Global transcription profiling reveals differential responses to chronic nitrogen stress and putative nitrogen regulatory components in Arabidopsis. BMC Genom. 2007, 8, 1–17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Krapp, A.; Berthomé, R.; Orsel, M.; Mercey-Boutet, S.; Yu, A.; Castaings, L.; Elftieh, S.; Major, H.; Renou, J.-P.; Daniel-Vedele, F. Arabidopsis roots and shoots show distinct temporal adaptation patterns toward nitrogen starvation. Plant Physiol. 2011, 157, 1255–1282. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Chardon, F.; Barthélémy, J.; Daniel-Vedele, F.; Masclaux-Daubresse, C. Natural variation of nitrate uptake and nitrogen use efficiency in Arabidopsis thaliana cultivated with limiting and ample nitrogen supply. J. Exp. Bot. 2010, 61, 2293–2302. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Menz, J.; Range, T.; Trini, J.; Ludewig, U.; Neuhäuser, B. Molecular basis of differential nitrogen use efficiencies and nitrogen source preferences in contrasting Arabidopsis accessions. Sci. Rep. 2018, 8, 1–11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Dos Santos, T.B.; Soares, J.D.; Lima, J.E.; Silva, J.C.; Ivamoto, S.T.; Baba, V.Y.; Souza, S.G.; Lorenzetti, A.P.; Paschoal, A.R.; Meda, A.R. An integrated analysis of mRNA and sRNA transcriptional profiles in Coffea arabica L. roots: Insights on nitrogen starvation responses. Funct. Integr. Genom. 2019, 19, 151–169. [Google Scholar] [CrossRef] [Green Version]
  13. Keller, M.; Arnink, K.J.; Hrazdina, G. Interaction of nitrogen availability during bloom and light intensity during veraison. I. Effects on grapevine growth, fruit development, and ripening. Am. J. Enol. Vitic. 1998, 49, 333–340. [Google Scholar]
  14. Weber, K.; Burow, M. Nitrogen–essential macronutrient and signal controlling flowering time. Physiol. Plant. 2018, 162, 251–260. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Ferrante, A.; Spinardi, A.; Maggiore, T.; Testoni, A.; Gallina, P.M. Effect of nitrogen fertilisation levels on melon fruit quality at the harvest time and during storage. J. Sci. Food Agric. 2008, 88, 707–713. [Google Scholar] [CrossRef]
  16. Carranca, C.; Brunetto, G.; Tagliavini, M. Nitrogen nutrition of fruit trees to reconcile productivity and environmental concerns. Plants 2018, 7, 4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Albornoz, F. Crop responses to nitrogen overfertilization: A review. Sci. Hortic. 2016, 205, 79–83. [Google Scholar] [CrossRef]
  18. Sete, P.B.; Comin, J.J.; Ciotta, M.N.; Salume, J.A.; Thewes, F.; Brackmann, A.; Toselli, M.; Nava, G.; Rozane, D.E.; Loss, A. Nitrogen fertilization affects yield and fruit quality in pear. Sci. Hortic. 2019, 258, 108782. [Google Scholar] [CrossRef]
  19. Neilsen, G.; Neilsen, D.; Peryea, F.; Fallahi, E.; Fallahi, B. Effects of mineral nutrition on fruit quality and nutritional disorders in apples. In Proceedings of the VI International Symposium on Mineral Nutrition of Fruit Crops 868, Faro, Portugal, 19–21 May 2008; pp. 49–60. [Google Scholar]
  20. Zhao, H.; Sun, S.; Zhang, L.; Yang, J.; Wang, Z.; Ma, F.; Li, M. Carbohydrate metabolism and transport in apple roots under nitrogen deficiency. Plant. Physiol. Biochem. 2020, 155, 455–463. [Google Scholar] [CrossRef] [PubMed]
  21. Han, Y.; Vimolmangkang, S.; Soria-Guerra, R.E.; Rosales-Mendoza, S.; Zheng, D.; Lygin, A.V.; Korban, S.S. Ectopic expression of apple F3′ H genes contributes to anthocyanin accumulation in the Arabidopsis tt7 mutant grown under nitrogen stress. Plant Physiol. 2010, 153, 806–820. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Sun, T.; Zhang, J.; Zhang, Q.; Li, X.; Li, M.; Yang, Y.; Zhou, J.; Wei, Q.; Zhou, B. Integrative physiological, transcriptome, and metabolome analysis reveals the effects of nitrogen sufficiency and deficiency conditions in apple leaves and roots. Environ. Exp. Bot. 2021, 192, 104633. [Google Scholar] [CrossRef]
  23. Zhang, X.; Tahir, M.M.; Li, S.; Mao, J.; Nawaz, M.A.; Liu, Y.; Li, K.; Xing, L.; Niu, J.; Zhang, D. Transcriptome analysis reveals the inhibitory nature of high nitrate during adventitious roots formation in the apple rootstock. Physiol. Plant. 2021, 173, 867–882. [Google Scholar] [CrossRef]
  24. Lee, Y.; Kweon, H.J.; Park, M.-Y.; Lee, D. Field Assessment of Macronutrients and Nitrogen in Apple Leaves Using a Chlorophyll Meter. Hort. Technol. 2019, 29, 300–307. [Google Scholar] [CrossRef] [Green Version]
  25. Daccord, N.; Celton, J.-M.; Linsmith, G.; Becker, C.; Choisne, N.; Schijlen, E.; Van de Geest, H.; Bianco, L.; Micheletti, D.; Velasco, R. High-Quality de novo assembly of the apple genome and methylome dynamics of early fruit development. Nat. Genet. 2017, 49, 1099–1106. [Google Scholar] [CrossRef] [PubMed]
  26. Armijo, G.; Gutierrez, R.A. Emerging Players in the Nitrate Signaling Pathway. Mol. Plant 2017, 10, 1019–1022. [Google Scholar] [CrossRef]
  27. Undurraga, S.F.; Ibarra-Henríquez, C.; Fredes, I.; Álvarez, J.M.; Gutiérrez, R.A. Nitrate signaling and early responses in Arabidopsis roots. J. Exp. Bot. 2017, 68, 2541–2551. [Google Scholar] [CrossRef]
  28. Wang, R.; Tischner, R.; Gutiérrez, R.A.; Hoffman, M.; Xing, X.; Chen, M.; Coruzzi, G.; Crawford, N.M. Genomic analysis of the nitrate response using a nitrate reductase-Null mutant of Arabidopsis. Plant Physiol. 2004, 136, 2512–2522. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Takahashi, M.; Sasaki, Y.; Ida, S.; Morikawa, H. Nitrite reductase gene enrichment improves assimilation of NO2 in Arabidopsis. Plant Physiol. 2001, 126, 731–741. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Konishi, N.; Ishiyama, K.; Matsuoka, K.; Maru, I.; Hayakawa, T.; Yamaya, T.; Kojima, S. NADH-dependent glutamate synthase plays a crucial role in assimilating ammonium in the Arabidopsis root. Physiol. Plant. 2014, 152, 138–151. [Google Scholar] [CrossRef] [PubMed]
  31. Asim, M.; Ullah, Z.; Xu, F.; An, L.; Aluko, O.O.; Wang, Q.; Liu, H. Nitrate Signaling, Functions, and Regulation of Root System Architecture: Insights from Arabidopsis Thaliana. Genes 2020, 11, 633. [Google Scholar] [CrossRef]
  32. Liu, K.-H.; Niu, Y.; Konishi, M.; Wu, Y.; Du, H.; Chung, H.S.; Li, L.; Boudsocq, M.; McCormack, M.; Maekawa, S. Discovery of nitrate–CPK–NLP signalling in central nutrient–growth networks. Nature 2017, 545, 311–316. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Zhang, H.; Forde, B.G. An Arabidopsis MADS box gene that controls nutrient-Induced changes in root architecture. Science 1998, 279, 407–409. [Google Scholar] [CrossRef]
  34. Krouk, G.; Mirowski, P.; LeCun, Y.; Shasha, D.E.; Coruzzi, G.M. Predictive network modeling of the high-Resolution dynamic plant transcriptome in response to nitrate. Genome Biol. 2010, 11, 1–19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Vidal, E.A.; Moyano, T.C.; Riveras, E.; Contreras-López, O.; Gutiérrez, R.A. Systems approaches map regulatory networks downstream of the auxin receptor AFB3 in the nitrate response of Arabidopsis thaliana roots. Proc. Natl. Acad. Sci. USA 2013, 110, 12840–12845. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Rubin, G.; Tohge, T.; Matsuda, F.; Saito, K.; Scheible, W.-R.D. Members of the LBD family of transcription factors repress anthocyanin synthesis and affect additional nitrogen responses in Arabidopsis. Plant Cell 2009, 21, 3567–3584. [Google Scholar] [CrossRef] [Green Version]
  37. Maeda, Y.; Konishi, M.; Kiba, T.; Sakuraba, Y.; Sawaki, N.; Kurai, T.; Ueda, Y.; Sakakibara, H.; Yanagisawa, S. A NIGT1-Centred transcriptional cascade regulates nitrate signalling and incorporates phosphorus starvation signals in Arabidopsis. Nat. Commun. 2018, 9, 1376. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Hu, B.; Jiang, Z.; Wang, W.; Qiu, Y.; Zhang, Z.; Liu, Y.; Li, A.; Gao, X.; Liu, L.; Qian, Y. Nitrate–NRT1. 1B–SPX4 cascade integrates nitrogen and phosphorus signalling networks in plants. Nat. Plants 2019, 5, 401–413. [Google Scholar] [CrossRef] [PubMed]
  39. Hu, B.; Chu, C. Nitrogen–Phosphorus interplay: Old story with molecular tale. New Phytol. 2020, 225, 1455–1460. [Google Scholar] [CrossRef] [Green Version]
  40. Krouk, G.; Lacombe, B.; Bielach, A.; Perrine-Walker, F.; Malinska, K.; Mounier, E.; Hoyerova, K.; Tillard, P.; Leon, S.; Ljung, K. Nitrate-Regulated auxin transport by NRT1. 1 defines a mechanism for nutrient sensing in plants. Dev. Cell 2010, 18, 927–937. [Google Scholar] [CrossRef]
  41. Teng, Y.; Liang, Y.; Wang, M.; Mai, H.; Ke, L. Nitrate Transporter 1.1 is involved in regulating flowering time via transcriptional regulation of FLOWERING LOCUS C in Arabidopsis thaliana. Plant Sci. 2019, 284, 30–36. [Google Scholar] [CrossRef] [PubMed]
  42. Li, G.; Tillard, P.; Gojon, A.; Maurel, C. Dual regulation of root hydraulic conductivity and plasma membrane aquaporins by plant nitrate accumulation and high-Affinity nitrate transporter NRT2. 1. Plant Cell Physiol. 2016, 57, 733–742. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Araus, V.; Swift, J.; Alvarez, J.M.; Henry, A.; Coruzzi, G.M. A balancing act: How plants integrate nitrogen and water signals. J. Exp. Bot. 2020, 71, 4442–4451. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Verly, C.; Djoman, A.C.R.; Rigault, M.; Giraud, F.; Rajjou, L.; Saint-Macary, M.-E.; Dellagi, A. Plant Defense Stimulator Mediated Defense Activation Is Affected by Nitrate Fertilization and Developmental Stage in Arabidopsis thaliana. Front. Plant Sci. 2020, 11, 583. [Google Scholar] [CrossRef] [PubMed]
  45. Cheng, C.-l.; Acedo, G.N.; Cristinsin, M.; Conkling, M.A. Sucrose mimics the light induction of Arabidopsis nitrate reductase gene transcription. Proc. Natl. Acad. Sci. USA 1992, 89, 1861–1864. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Fritz, C.; Palacios-Rojas, N.; Feil, R.; Stitt, M. Regulation of secondary metabolism by the carbon–Nitrogen status in tobacco: Nitrate inhibits large sectors of phenylpropanoid metabolism. Plant J. 2006, 46, 533–548. [Google Scholar] [CrossRef] [PubMed]
  47. Izumi, M.; Tsunoda, H.; Suzuki, Y.; Makino, A.; Ishida, H. RBCS1A and RBCS3B, two major members within the Arabidopsis RBCS multigene family, function to yield sufficient Rubisco content for leaf photosynthetic capacity. J. Exp. Bot. 2012, 63, 2159–2170. [Google Scholar] [CrossRef]
  48. Warren, C.R.; Dreyer, E.; Adams, M.A. Photosynthesis-Rubisco relationships in foliage of Pinus sylvestris in response to nitrogen supply and the proposed role of Rubisco and amino acids as nitrogen stores. Trees 2003, 17, 359–366. [Google Scholar] [CrossRef]
  49. Shi, J.; Yi, K.; Liu, Y.; Xie, L.; Zhou, Z.; Chen, Y.; Hu, Z.; Zheng, T.; Liu, R.; Chen, Y. Phospho enol pyruvate carboxylase in Arabidopsis leaves plays a crucial role in carbon and nitrogen metabolism. Plant Physiol. 2015, 167, 671–681. [Google Scholar] [CrossRef] [Green Version]
  50. Kiba, T.; Kudo, T.; Kojima, M.; Sakakibara, H. Hormonal control of nitrogen acquisition: Roles of auxin, abscisic acid, and cytokinin. J. Exp. Bot. 2011, 62, 1399–1409. [Google Scholar] [CrossRef]
  51. Vega, A.; O’Brien, J.A.; Gutiérrez, R.A. Nitrate and hormonal signaling crosstalk for plant growth and development. Curr. Opin. Plant Biol. 2019, 52, 155–163. [Google Scholar] [CrossRef]
  52. Gray, W.M.; Kepinski, S.; Rouse, D.; Leyser, O.; Estelle, M. Auxin regulates SCF TIR1-Dependent degradation of AUX/IAA proteins. Nature 2001, 414, 271–276. [Google Scholar] [CrossRef]
  53. Poitout, A.; Crabos, A.; Petřík, I.; Novák, O.; Krouk, G.; Lacombe, B.; Ruffel, S. Responses to systemic nitrogen signaling in Arabidopsis roots involve trans-Zeatin in shoots. Plant Cell 2018, 30, 1243–1257. [Google Scholar] [CrossRef] [Green Version]
  54. Su, H.; Wang, T.; Ju, C.; Deng, J.; Zhang, T.; Li, M.; Tian, H.; Wang, C. Abscisic acid signaling negatively regulates nitrate uptake via phosphorylation of NRT1. 1 by SnRK2s in Arabidopsis. J. Integr. Plant Biol. 2021, 63, 597–610. [Google Scholar] [CrossRef]
  55. Fait, A.; Nesi, A.N.; Angelovici, R.; Lehmann, M.; Pham, P.A.; Song, L.; Haslam, R.P.; Napier, J.A.; Galili, G.; Fernie, A.R. Targeted enhancement of glutamate-To-γ-Aminobutyrate conversion in Arabidopsis seeds affects carbon-Nitrogen balance and storage reserves in a development-Dependent manner. Plant Physiol. 2011, 157, 1026–1042. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Bouche, N.; Fromm, H. GABA in plants: Just a metabolite? Trends Plant Sci. 2004, 9, 110–115. [Google Scholar] [CrossRef]
  57. Barbosa, J.M.; Singh, N.K.; Cherry, J.H.; Locy, R.D. Nitrate uptake and utilization is modulated by exogenous γ-Aminobutyric acid in Arabidopsis thaliana seedlings. Plant Physiol. Biochem. 2010, 48, 443–450. [Google Scholar] [CrossRef]
  58. Morikawa, H.; Takahashi, M.; Sakamoto, A.; Matsubara, T.; Arimura, G.-I.; Kawamura, Y.; Fukunaga, K.; Fujita, K.; Sakurai, N.; Hirata, T. Formation of unidentified nitrogen in plants: An implication for a novel nitrogen metabolism. Planta 2004, 219, 14–22. [Google Scholar] [CrossRef]
  59. Chang, S.; Puryear, J.; Cairney, J. A simple and efficient method for isolating RNA from pine trees. Plant Mol. Biol. Report. 1993, 11, 113–116. [Google Scholar] [CrossRef]
  60. Liu, C.-H.; Di, Y.P. Analysis of RNA sequencing data using CLC Genomics Workbench. In Molecular Toxicology Protocols; Springer: Berlin/Heidelberg, Germany, 2020; pp. 61–113. [Google Scholar]
  61. Howe, E.A.; Sinha, R.; Schlauch, D.; Quackenbush, J. RNA-Seq analysis in MeV. Bioinformatics 2011, 27, 3209–3210. [Google Scholar] [CrossRef] [Green Version]
  62. Tian, T.; Liu, Y.; Yan, H.; You, Q.; Yi, X.; Du, Z.; Xu, W.; Su, Z. agriGO v2.0: A GO analysis toolkit for the agricultural community, 2017 update. Nucleic Acids Res. 2017, 45, W122–W129. [Google Scholar] [CrossRef] [PubMed]
  63. Supek, F.; Bošnjak, M.; Škunca, N.; Šmuc, T. REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS ONE 2011, 6, e21800. [Google Scholar] [CrossRef] [Green Version]
  64. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 2009, 4, 44–57. [Google Scholar] [CrossRef] [PubMed]
  65. Jung, S.; Lee, T.; Cheng, C.-H.; Buble, K.; Zheng, P.; Yu, J.; Humann, J.; Ficklin, S.P.; Gasic, K.; Scott, K. 15 years of GDR: New data and functionality in the Genome Database for Rosaceae. Nucleic Acids Res. 2019, 47, D1137–D1145. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Szklarczyk, D.; Gable, A.L.; Nastou, K.C.; Lyon, D.; Kirsch, R.; Pyysalo, S.; Doncheva, N.T.; Legeay, M.; Fang, T.; Bork, P. The STRING database in 2021: Customizable protein–Protein networks, and functional characterization of user-Uploaded gene/measurement sets. Nucleic Acids Res. 2021, 49, D605–D612. [Google Scholar] [CrossRef] [PubMed]
  67. Bowen, J.; Ireland, H.S.; Crowhurst, R.; Luo, Z.; Watson, A.E.; Foster, T.; Gapper, N.; Giovanonni, J.J.; Mattheis, J.P.; Watkins, C. Selection of low-Variance expressed Malus x domestica (apple) genes for use as quantitative PCR reference genes (housekeepers). Tree Genet. Genomes 2014, 10, 751–759. [Google Scholar] [CrossRef]
  68. Tellmann, G. The E-Method: A highly accurate technique for gene-Expression analysis. Nat. Methods 2006, 3, i–ii. [Google Scholar] [CrossRef]
  69. Lohse, M.; Nagel, A.; Herter, T.; May, P.; Schroda, M.; Zrenner, R.; Tohge, T.; Fernie, A.R.; Stitt, M.; Usadel, B. Mercator: A Fast and Simple Web Server for Genome Scale Functional Annotation of Plant Sequence Data. Plant Cell Environ. 2014, 5, 1250–1258. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. RNA-seq experiment design for the study of apple N stress response signals and DEG analysis. (a) A scheme for apple N stress response study. Leaf tissues of ‘Fuji’/M9 apple trees treated with different N fertilization gradients showing a different level of foliar N content were sampled for transcriptome analysis. (b) Leaf N content among four N fertilization groups. Lines on the bars show SE of the mean and different letters indicate significant difference as determined by DMRT (p < 0.05). (c) Principal component analysis of RNA-seq expression data. The top 5000 genes listed according to the variance rank were chosen as informative genes and used for PCA. (d) Venn diagram of DEGs with criteria of |fold change| ≥ 1.5 and FDR ≤ 0.05.
Figure 1. RNA-seq experiment design for the study of apple N stress response signals and DEG analysis. (a) A scheme for apple N stress response study. Leaf tissues of ‘Fuji’/M9 apple trees treated with different N fertilization gradients showing a different level of foliar N content were sampled for transcriptome analysis. (b) Leaf N content among four N fertilization groups. Lines on the bars show SE of the mean and different letters indicate significant difference as determined by DMRT (p < 0.05). (c) Principal component analysis of RNA-seq expression data. The top 5000 genes listed according to the variance rank were chosen as informative genes and used for PCA. (d) Venn diagram of DEGs with criteria of |fold change| ≥ 1.5 and FDR ≤ 0.05.
Plants 10 02649 g001
Figure 2. Expression profile of apple genes involved in a conserved N uptake, assimilation, and transport pathway.
Figure 2. Expression profile of apple genes involved in a conserved N uptake, assimilation, and transport pathway.
Plants 10 02649 g002
Figure 3. Investigation of DEGs functionally associated with nitrogen stress responses. (a) k-means clustering of rootstock 2212 DEGs with criteria of fold change above 1.5 and FDR below 0.05 (k = 5). (be) Gene ontology enrichment of DEGs in each cluster (biological process). All clusters except C4 (up in N-deficient) showed significant GO terms with a criteria of FDR 0.05. (b) C1: up in N-deficient; (c) C2: down in N-excessive; (d) C3: positive N correlation; (e) C5: up in N-excessive.
Figure 3. Investigation of DEGs functionally associated with nitrogen stress responses. (a) k-means clustering of rootstock 2212 DEGs with criteria of fold change above 1.5 and FDR below 0.05 (k = 5). (be) Gene ontology enrichment of DEGs in each cluster (biological process). All clusters except C4 (up in N-deficient) showed significant GO terms with a criteria of FDR 0.05. (b) C1: up in N-deficient; (c) C2: down in N-excessive; (d) C3: positive N correlation; (e) C5: up in N-excessive.
Plants 10 02649 g003
Figure 4. Functional annotation of DEGs: (a) photosynthesis; (b) cell wall; (c) secondary metabolism; (d) hormone metabolism; (e) transport.
Figure 4. Functional annotation of DEGs: (a) photosynthesis; (b) cell wall; (c) secondary metabolism; (d) hormone metabolism; (e) transport.
Plants 10 02649 g004
Figure 5. A proposed scheme for candidate metabolic signals associated with N-deficient/excessive stress responses in apple.
Figure 5. A proposed scheme for candidate metabolic signals associated with N-deficient/excessive stress responses in apple.
Plants 10 02649 g005
Figure 6. Network of protein–protein interaction (PPI) of five DEG clusters (C1–C5). (ae) The PPI network in each DEG cluster was constructed with high confidence (0.7) using the STRING database: (a) C1, up in N-deficient; (b) C2, down in N-excessive; (c) C3, positive N correlation; (d) C4, down in N-deficient; (e) C5, up in N-excessive.
Figure 6. Network of protein–protein interaction (PPI) of five DEG clusters (C1–C5). (ae) The PPI network in each DEG cluster was constructed with high confidence (0.7) using the STRING database: (a) C1, up in N-deficient; (b) C2, down in N-excessive; (c) C3, positive N correlation; (d) C4, down in N-deficient; (e) C5, up in N-excessive.
Plants 10 02649 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lee, Y.; Do, V.G.; Kim, S.; Kweon, H. Identification of Genes Associated with Nitrogen Stress Responses in Apple Leaves. Plants 2021, 10, 2649. https://doi.org/10.3390/plants10122649

AMA Style

Lee Y, Do VG, Kim S, Kweon H. Identification of Genes Associated with Nitrogen Stress Responses in Apple Leaves. Plants. 2021; 10(12):2649. https://doi.org/10.3390/plants10122649

Chicago/Turabian Style

Lee, Youngsuk, Van Giap Do, Seonae Kim, and Hunjoong Kweon. 2021. "Identification of Genes Associated with Nitrogen Stress Responses in Apple Leaves" Plants 10, no. 12: 2649. https://doi.org/10.3390/plants10122649

APA Style

Lee, Y., Do, V. G., Kim, S., & Kweon, H. (2021). Identification of Genes Associated with Nitrogen Stress Responses in Apple Leaves. Plants, 10(12), 2649. https://doi.org/10.3390/plants10122649

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop