Next Article in Journal
The Electrical Conductivity and Drying Method Changed the Secondary Metabolite Content and Photoprotective Effects of Hydroponically Cultivated Agastache rugosa Kuntze
Next Article in Special Issue
Revealing Key Genes and Pathways in Potato Scab Disease Resistance through Transcriptome Analysis
Previous Article in Journal
Advances and Developments in Monitoring and Inversion of the Biochemical Information of Crop Nutrients Based on Hyperspectral Technology
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Agronomy
Volume 13
Issue 8
10.3390/agronomy13082164
agronomy-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptome Analysis of Nitrogen-Deficiency-Responsive Genes in Two Potato Cultivars

by
Qiaorong Wei
1,2,3,†,
Yanbin Yin
1,2,4,†,
Bin Deng
1,†,
Xuewei Song
1,2,
Zhenping Gong
1,* and
Ying Shi
1,2,*
1
College of Agriculture, Northeast Agricultural University, Harbin 150038, China
2
Key Laboratory of Germplasm Enhancement, Physiology and Ecology of Food Crops in Cold Region, Ministry of Education, Harbin 150038, China
3
National Key Laboratory of Smart Farm Technology and System, Key Laboratory of Soybean Biology in Chinese Education Ministry, Northeast Agricultural University, Harbin 150038, China
4
Rice Research Institute, Shenyang Agricultural University, Shenyang 110065, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Agronomy 2023, 13(8), 2164; https://doi.org/10.3390/agronomy13082164
Submission received: 11 July 2023 / Revised: 12 August 2023 / Accepted: 16 August 2023 / Published: 18 August 2023
(This article belongs to the Special Issue Applications of Genomics or Transcriptomics Approaches in Phytopathology)
Browse Figures
Review Reports Versions Notes

Abstract

:
The potato is the third largest food crop, and nitrogen fertilizer is important for increasing potato yields; however, the shallow root system of potatoes causes the nitrogen fertilizer utilization rate to be low, which results in waste and environmental pollution, meaning that high nitrogen efficiency breeding is highly significant for potatoes. In the high nitrogen efficiency breeding of potatoes, genes with a nitrogen-deficient response should first be identified, and RNA-seq is an efficient method for identifying nitrogen-deficiency-response genes. In this study, two potato cultivars, Dongnong 322 (DN322) and Dongnong 314 (DN314), were utilized, and two nitrogen fertilizer application rates (N0 and N1) were set for both cultivars. Through the determination of physiological indicators, we identified that DN314 is more sensitive to nitrogen fertilizer, while DN322 is relatively insensitive to nitrogen fertilizer. Samples were taken at the seedling and tuber formation stage. At the seedling stage, DN322 and DN314 had 573 and 150 differentially expressed genes (DEGs), while at the tuber formation stage, they had 59 and 1905 DEGs, respectively. A total of three genes related to a low-nitrogen response were obtained via the combined analysis of differentially expressed genes (DEGs) and weighted correlation network analysis (WGCNA), of which two genes were obtained at the tuber formation stage and one gene in the seedling stage, providing theoretical guidance for the high nitrogen efficiency breeding of potatoes.
Keywords:
potatoes; RNA-seq; DEGs; WGCNA

1. Introduction

The potato (Solanum tuberosum L.) is the third largest food crop [1], and is rich in starch, protein, vitamins, and other nutrients [2]. In China, potatoes are both vegetables and staples, so increasing the potato yield is particularly important to meet individual nutritional needs. The use of nitrogen fertilizers is important for increasing potato production [3], but due to the shallow root system of potatoes, the nitrogen fertilizer utilization rate of potatoes is low, generally being less than 50% [4,5], which causes significant waste and produces environmental pollution. The environmental pollution of nitrogen fertilizer is mainly divided into the following aspects: first, the fertilizer pollutes water resources, resulting in water eutrophication; second, it pollutes the soil, causing soil compaction; and third, it increases N2O emissions, polluting the air [6,7,8,9]. In the potato production process, reducing nitrogen input and improving nitrogen fertilizer efficiency are important research topics. It has been proven that a reasonable combination of water and fertilizer [10], or the use of slow-release nitrogen fertilizer, can improve nitrogen fertilizer utilization [11,12]. However, these two methods require producers to have a greater level of experience, and the cost of slow-release fertilizer is relatively high; therefore, high nitrogen efficiency breeding seems to be a better way to improve potatoes’ nitrogen fertilizer utilization.
In rice and maize, genes associated with nitrogen uptake, utilization, and NUE have been cloned that can increase crop yields or improve crop quality. The nitrate transporter genes NRT1.1b [13], NRT2.1 [14,15], NRT2.3a, and NRT2.3b [16,17] were reported as being able to transfer NO3− ions, while the nitrate reductase OsNR2 was reported as being able to reduce NO3− ions; these genes increase rice yield by improving NUE. Some transcriptional regulators have also been reported to increase rice yield by improving NUE; the transcriptional regulator Ghd7 directly binds to the promoter of ARE1 to improve NUE and increase rice yields [18], and OsDREB1C directly binds to the promoters of NR2, OsNRT2.4, and OsNRT1.1B [19] and activates their expression. The overexpression lines of OsDREB1C increased yields by 41.3–68.3% [19]. In maize, THP9 increases the corn protein content by improving NUE [20].
However, in potatoes, fewer genes have been reported to affect NUE; In this study, the RNA-seq method was used to identify low-nitrogen-response genes at seedling and tuber formation stages in potatoes. First, a weighted correlation network analysis (WGCNA) was used to create a co-expression module; second, the phenotype and each module were analyzed to determine the nitrogen-deficiency-response module; and third, we conducted an analysis of nitrogen-deficiency-response genes by combining the nitrogen-deficiency-response expression modules with the differentially expressed genes (DEGs). A total of three nitrogen-deficient genes were identified, including one gene at the seedling stage and two genes at the tuber formation stage, which provided theoretical knowledge for the high nitrogen efficiency breeding of potatoes. In this study, we used WGCNA and DEGs to identify three genes related to the nitrogen-deficient response; these are new genes which can be used for potato breeding to improve potato nitrogen use efficiency and enrich the potato nitrogen-response gene pathways.

2. Materials and Methods

2.1. Plant Materials and Experimental Treatment

Two potato cultivars, Dongnong322 (DN322, with a growth period of 90 days) and Dongnong314 (DN314, with a growth period of 90 days), were planted in pots with soil at the Potato Research Institute of Northeast Agricultural University in Harbin, China. The size of the pots was 33 cm × 60 cm; the volume was 0.51 m3. Detoxified seed potatoes (of about 50 g) of DN322 and DN314 were provided by the Potato Research Institute of Northeast Agricultural University. The two potato cultivars were divided into the N0 (nitrogen deficient) and N1 (normal nitrogen) groups for nitrogen treatment. Plants in the N0 group were treated with 0 g N (urea, N: 46%), 6 g P (P2O5: 44%), and 6 g K (K2O 50%); plants in the N1 group were treated with 4.5 g N (urea, N: 46%), 6 g P (P2O5: 44%), and 6 g K (K2O 50%). At the seedling and tuber formation stages, leaves were sampled and frozen in liquid nitrogen and stored at −80 °C for RNA-seq.

2.2. The Measurement of Plant Dry Matter Weight, Nitrogen Content, and Other Physiological Parameters

At the seedling and tuber formation stages, samples were taken from three plants from both DN322 and DN314 for every N treatment, following which the plants were divided into the root, stem, leaf, and tuber. Methods for determining plant dry matter were based on previous reports [2]; in brief, the divided parts of the potato were heated at 105 °C for 30 min in an oven, and then dried to a constant weight at 70 °C. Total nitrogen content was measured using previously defined methods [21]; in brief, the dry matter of each tissue was ground into powder and passed through an 800-mesh sieve. A 0.2 g sample was taken, and the nitrogen content was measured using a concentrated sulfuric acid–hydrogen peroxide solution, where the total nitrogen content = nitrogen content of different tissue × total dry matter of different tissue. N utilization efficiency (NUE) = (N1NC − N0NC)/TNA. N1NC represents the nitrogen accumulation in nitrogen-treated plants, N0NC represents the nitrogen accumulation in non-nitrogen-treated plants, and TNA represents the total nitrogen application. The net photosynthetic rate was measured using a method described previously [22]; in brief, on sunny days from 9:00 a.m. to 11:30 a.m., the net photosynthetic rate was determined with a Yaxin photosynthesis analyzer (XY-1101). Root vitality was also measured using previously defined methods [23]: First, the standard curve was developed using different concentrations of 1,3,5-triphenylformazan (TTF), and then 0.4% 2,3,5-triphenyl tetrazolium chloride (TTC) and a phosphoric acid buffer were configured into a mixture, where the roots were immersed in the mixture and kept warm for 1–3 h at 37 °C in the dark. Then, 1M sulfuric acid was added to react for 10 min, and finally, the roots were ground with ethyl acetate and a small amount of quartz sand. Root viability was determined according to the standard curve, where root vitality = TTC reduction amount/(root weight × kept warm time).

2.3. RNA Extraction, RNA Sequencing, and RT-qPCR Analysis

RNA extraction and RT-qPCR analysis were analyzed via methods described previously [24,25]. In brief, an RNA extraction kit (Promega, Cat.LS1040) was used to extract the RNA according to the manufacturer’s method, and 3 μg was used to synthesize the first-strand cDNA (Thermo Fisher Scientific, Cat. K1621). Controlling RNA quantity and quality was performed in the same manner as in a previous study [26]; in brief, the quantity of RNA was measured with the Nanorop2000, whereas the RNA quality was verified using a 1% agarose gel, and after 150 V electrophoresis for 5 min, two sharp bands formed on the gel, indicating that the RNA quality was good. The RNA sequencing method performed was the same as that employed in another previous study [26]; in brief, the RNA sequencing process used was library establishment–library quality control sequencing. This study used the Illumina platform for sequencing, where the number of sequencing bases was Q30 > 91%, and the sequencing quality was ideal. For RT-qPCR, TB green premix Ex Taq (Takara, Cat. RR420A) and a QuantStudio 5 real-time PCR instrument (Applied Biosystems) were used according to the manufacturers’ protocols, and EIF-3e was used as the internal control gene [27]. In RT-qPCR analysis, we used three biological replicates for different treatments.

2.4. Gene Expression Levels and Differentially Expressed Gene (DEG) Analysis

Gene expression levels were estimated using fragments per kilobase of transcripts per million fragments mapped (FPKM), as previously described [28,29]. DEGs between N0 and N1 in DN322 or DN314 were determined with the DESeq2 [30] using a model based on the negative binomial distribution. The resulting p-values were adjusted using Benjamini and Hochberg’s approach for controlling the false discovery rate (FDR). An adjusted FDR ≤ 0.01 and log2 fold changes ≥1.5 were used as the thresholds of differential expression.
Gene ontology (GO) enrichment analysis of the low-nitrogen-response module genes was implemented via the GOseq R package based on the Wallenius non-central hyper-geometric distribution [31], and KOBAS software [32] was used to analyze the low-nitrogen-response module genes for KEGG enrichment [33].

2.5. Weighted Gene Co-Expression Network Analysis (WGCNA)

A total of 20,832 and 21,061 genes with an average FPKM > 1 from the seedling and tuber formation stages were selected for the WGCNA network analysis. The WGCNA (v1.29) package in R was used, as in previous studies [34,35]. The module similarity threshold was 0.25, and the minimum number of genes for the module was 30.

3. Results

3.1. Nitrogen-Deficiency Stress Affects Potato Growth and Development

We first investigated the effect of nitrogen-deficiency stress on dry matter accumulation in potatoes. The dry matter accumulation of DN322 and DN314 at the seedling stage decreased by 60.7% and 78.6%, and that at the tuber formation stage decreased by 45.1% and 80.8%, respectively (Figure 1B,C; Supplementary Table S1). Among them, during the potato tuber formation stage, the dry matter weight of DN322 and DN314 decreased by 45.1% and 80.8%, respectively (Figure 1D), and the dry matter weight of the roots, stems, and leaves of DN322 and DN314 also decreased significantly (Supplementary Figure S1A–F). In addition, nitrogen-deficiency stress also reduced nitrogen accumulation in the roots, stems, leaves, and tubers of DN322 and DN314 (Supplementary Figure S1G–L). We also investigated the effects of low-nitrogen stress on the physiological indicators of DN322 and DN314. At the seedling stage, the root viability of DN322 and DN314 decreased by 15.4% and 41.8%, and the tuber formation stage decreased by 24.0% and 42.5%, respectively (Figure 1E,F; Supplementary Table S1). Photosynthesis is of great significance for potato yield formation, and we investigated the effect of low-nitrogen stress on the net photosynthetic rate of DN322 and DN314, which decreased by 40.1% and 39.7% at the seedling stage and by 41.0% and 40.7% at the tuber formation stage, respectively (Figure 1G,H; Supplementary Table S1). Low-nitrogen stress also significantly reduced the plant height for DN322 and DN314 (Figure 1I).We measured the N utilization efficiency (NUE) of DN322 and DN314, demonstrating that DN314 had a higher NUE (Supplementary Figure S2H). In conclusion, nitrogen-deficiency stress affected potato growth and development, and DN314 is more sensitive to nitrogen than DN322.

3.2. High-Throughput RNA-seq to Analyze the Effect of Nitrogen-Deficiency Stress on DN322

To further investigate the effect of nitrogen-deficiency stress on potato molecular mechanisms, high-throughput RNA-seq was performed to analyze the effect of nitrogen-deficiency stress on DN322 at the seedling and tuber formation stages. The clean reads ranged from 20,457,525 to 31,881,882, and the Q30 ranged from 91.25% to 94.32%; the sequencing quality was determined to be ideal (Supplementary Table S2), enabling further analysis. The leaves of DN322 had a total of 21,431 gene expressions at the seedling stage and 573 differentially expressed genes (DEGs) under nitrogen-deficiency stress, including 298 upregulated expression genes and 275 downregulated expression genes (Figure 2A, Supplementary Table S4). We also randomly selected six DEGs and performed RT-qPCR to verify the accuracy of RNA-seq; we found that the expression trend of RT-qPCR was the same as that of RNA-seq under nitrogen-deficiency stress (Supplementary Figure S2A–F). The identified DEGs were annotated with 46 GO terms: 17 biological processes, 17 cellular components, and 12 molecular functions (Figure 2B). Similarly, KEGG annotation was performed, and the identified DEGs were mainly enriched via steroid biosynthesis, diterpenoid biosynthesis, flavonoid biosynthesis, and the valine, leucine, and isoleucine degradation pathways (Figure 2C). During the tuber formation stage, there were a total of 19,905 gene expressions under nitrogen-deficiency stress and 51 differentially expressed genes, of which 15 genes were upregulated and 36 were downregulated (Figure 2D, Supplementary Table S5). The DEGs were annotated with 23 GO terms: 11 biological processes, 7 cellular components, and 5 molecular functions (Figure 2E). The DEGs were also annotated with KEGG, and mainly photosynthesis and oxidative phosphorylation were found to be enriched (Figure 2F).

3.3. High-Throughput RNA-seq to Analyze the Effect of Nitrogen-Deficiency Stress on DN314

We also performed RNA-seq on leaves at the seedling and tuber formation stages of the potato cultivar DN314. The clean reads ranged from 19,127,622 to 22,712,674, and the Q30 ranged from 91.38% to 93.88% (Supplementary Table S3). The leaves of DN314 had a total of 20,832 expressed genes at the seedling stage and a total of 150 differentially expressed genes (DEGs) under nitrogen-deficiency stress, including 56 upregulated expression genes and 94 downregulated expression genes (Figure 3A, Supplementary Table S6). The identified DEGs were annotated with 34 GO terms: 13 biological processes, 12 cellular components, and 9 molecular functions (Figure 3B). Similarly, KEGG annotation was performed, and the identified DEGs were mainly enriched via protein processing in the endoplasmic, citrate cycle, and photosynthesis-antenna protein pathways (Figure 3C). During the tuber formation stage, there were a total of 21,061 gene expressions under nitrogen-deficiency stress and 1905 differentially expressed genes, of which 750 genes were upregulated and 1155 were downregulated (Figure 3D, Supplementary Table S7). The DEGs were annotated with 43 GO terms: 17 biological processes, 15 cellular components, and 11 molecular functions (Figure 3E). The DEGs were also annotated with KEGG regarding photosynthesis, photosynthesis-antenna proteins, carbon fixation in photosynthetic organisms, pentose and glucuronate interconversions, and starch and sucrose metabolism (Figure 3F).

3.4. Identification of Co-expressed Gene Clusters Using Two Potato Cultivars at the Seedling Stage

We used 12 leaf samples from the seedling stage of DN322 and DN314 for RNA sequencing. Using WGCNA, the expressed genes were divided into 14 modules (Figure 4A, Supplementary Table S8); among the modules, SM9 and SM10 were found to have the highest correlation (Supplementary Table S9), and the correlation analysis of each expression module gene’s plant dry weight, plant nitrogen content, and net photosynthetic rate was carried out. SM9 and SM10 were found to have a large correlation with each phenotype (Figure 4A). The proportion of DEGs in each module was statistically counted, and it was found that the proportion of nitrogen-fertilizer-response genes in SM9 and SM10 was the largest (Figure 4B); in summary, SM9 and SM10 had the greatest relationship with nitrogen-deficiency stress. Then, we functionally annotated the genes of SM9 and SM10 modules, and through the GO database annotation, we found that the genes in module SM9 were mainly involved in the metabolic process, cellular process, localization, developmental process, and other life processes, and that these genes mainly play a role in the membrane and organelle. They mainly possess binding and catalytic activity, as well as transporter activity associations (Figure 4C). Using the KEGG database, the gene function of module SM9 was annotated, and it was found that the genes of this module were mainly annotated to the signaling pathways of plant hormone signaling and starch and sucrose metabolism (Figure 4D). Similarly, we also annotated the module SM10 genes, and through GO database annotation, we found that the SM10 genes were mainly involved in the metabolic process, biological regulation, signaling, etc., and that these genes were mainly involved in life activities in the membrane and the extracellular region, where their main molecular functions were binding, catalytic activity, transporter activity, and enzyme regulator activity (Figure 4E). Module SM10 was analyzed using the KEGG database and was found to be mainly involved in plant hormone signal transduction and protein processing in the endoplasmic reticulum process (Figure 4F).

3.5. Identification of Co-Expressed Gene Clusters Using Two Potato Cultivars at the Tuber Formation Stage

We used 12 leaf RNAs in the tuber formation stage of DN322 and DN314 for RNA sequencing; using the method of WGCNA, the expression genes were divided into 14 modules (Figure 5A,B), and each module was associated with the plant dry weight, plant nitrogen content, and net photosynthetic rate in the tuber formation stage. It was found that module TM1 was the most associated with these phenotypic traits (Figure 5C). We then analyzed the differentially expressed genes in each model and found that module TM1 had the highest proportion of differentially expressed genes (Figure 5D). Then, we performed gene annotation on the genes of module TM1, and via GO database annotation analysis we found that the genes of module TM1 were mainly involved in the metabolic process, cellular process, and other pathways; these genes mainly play a role in the membrane and macromolecular complex and mainly carry out binding, catalytic activity, and transporter activity (Figure 5E). These genes of the TM1 module were analyzed using the KEGG database, and it was found that they were mainly located in the photosynthetic pathway and carbon metabolism pathway (Figure 5F).

3.6. Identification of Key Genes in the Nitrogen-Deficient Response

In order to further identify candidate genes with regard to the nitrogen-deficient response, this study used the method of a joint analysis of DEGs and WGCNA. In the seedling stage, we used DN322 and DN314 to obtain a total of nine co-expressed differential genes (Supplementary Figure S2), and then screened the co-expressed differential genes in the nitrogen-response modules SM9 and SM10. We found that one of the nine DEGs was present in both the SM9 and SM10 modules. This gene was Soltu.DM.05G003170, a transcription factor with the zinc-finger domain, which is involved in circadian rhythm. During the tuber formation stage, we used the same method; in module TM1, a total of two genes were found, Soltu.DM.04G009820 and Soltu.DM.04G016890. The molecular functions of Soltu.DM.04G009820 include post-translational modification, protein turnover, and chaperones, whereas Soltu.DM.04G016890 is involved in photosynthesis.

4. Discussion

4.1. Nitrogen-Deficiency Stress Affects Potato Growth

Nitrogen is a component of various cellular components and enzymes and plays an important role in the growth and development of plants. Nitrogen affects plant roots [36,37,38] and leaves [39], according to previous studies, and crop photosynthesis is also affected by nitrogen [30,31,32,33,34,35,36,37,38,39,40,41,42]. In this study, we first investigated the dry matter accumulation of two potato cultivars, DN322 and DN314, at the seedling and tuber formation stages under normal-nitrogen fertilizer (N1) and low-nitrogen conditions (N0), and we found that under low-nitrogen conditions, the dry matter weight of the roots, stems, leaves, and tubers of these two potato cultivars decreased significantly (Figure 1, Supplementary Figure S1, and Supplementary Table S1), indicating that nitrogen has an important effect on potato dry matter accumulation, which are the same findings as those of our previous study. We also investigated the nitrogen accumulation of DN322 and DN314 at the seedling and tuber formation stages, and found that nitrogen accumulation in the root, stem, leaf, and tuber was significantly lower under low-nitrogen conditions than that under normal-nitrogen conditions (Supplementary Table S1). Net photosynthesis and root vitality were significantly lower under low-nitrogen conditions (Figure 1E–H). In addition, DN314 was more sensitive to nitrogen fertilizer than DN322 in terms of the decrease in dry matter weight and root activity (Supplementary Figure S1). Under nitrogen-deficiency stress, the dry matter weight, nitrogen uptake, photosynthesis, and root activity of DN314 decreased more than that of DN322, indicating that DN314 was more sensitive to nitrogen fertilizer.

4.2. RNA-seq Is an Effective Method for the Identification of Potato Nitrogen-Deficiency-Response Genes

For seed plants, such as rice, maize, and soybeans, genes can be mapped using the QTL method by either constructing a genetic population [43,44,45,46], or using GWAS analysis methods on the natural population [47,48]. However, as potatoes are highly heterozygous and have self-incompatibility [27,49], cloning genes using QTL and GWAS methods is relatively difficult. With the development of sequencing technology, RNA sequencing technology has become an efficient method to identify candidate potato genes. In previous studies, researchers used low-nitrogen vs. normal-nitrogen or high-nitrogen vs. normal-nitrogen treatments and took samples at a certain period for RNA sequencing to obtain the DEGs, and then annotated these DEGs to determine gene function and to analyze the nitrogen-fertilizer-response genes of potatoes [50,51]. In this study, we selected two potato cultivars and selected the seedling and tuber formation stages to observe the response mechanism of potatoes to low nitrogen before and after potato tuber formation. We found that DN322 and DN314 had a total of 150 and 573 DEGs in the seedling stage, and a total of 1905 and 51 DEGs in the tuber formation stage, respectively. These DEGs can guide the study of potato nitrogen-fertilizer-response genes.

4.3. Effects of Nitrogen-Deficiency Stress on the Molecular Functions of Different Nitrogen-Sensitive Cultivars

Previous studies have demonstrated how nitrogen deficiency affects many pathways, such as nitrogen metabolism, transcription factors, and hormone signaling [52,53,54]. In this study, DN314 was more sensitive to nitrogen fertilizers than DN322. To explore its molecular mechanisms, RNA sequencing was utilized. At the seedling stage, the two varieties of enrichment pathways observed were not identical, as DN322 is mainly enriched via several biosynthetic pathways, such as steroid biosynthesis, while DN314 is mainly enriched via protein processing in the endoplasmic reticulum and photosynthesis-antenna proteins (Figure 2C and Figure 3C, respectively). Photosynthesis [55,56,57] is of great significance for plant growth and development, and the DEGs of DN314 are mainly enriched via the photosynthesis-antenna protein pathway, which may be an important reason underlying why DN314 is more sensitive to nitrogen fertilizer. During the tuber formation stage, the pathways with the highest degree of enrichment of DN322 and DN314 DEGs were photosynthetic pathways, while the differential genes of DN314, which are sensitive to nitrogen fertilizer, were also enriched via the photosynthesis-antenna protein pathway and the starch and sucrose metabolism pathways. In summary, nitrogen-deficiency stress mainly affects the photosynthetic and synthesis pathways, while DN314 is more sensitive to nitrogen fertilizer as more DEGs were enriched via the photosynthetic and sucrose metabolism pathways.

4.4. WGCNA and DEGs Were Combined to Obtain Candidate Genes for the Potato Nitrogen-Deficiency Response

In a previous study, identifying DEGs with different nitrogen treatments provided a lot of useful information, but there were too many DEGs to make an efficient selection and these DEGs were not able to be combined with phenotypic traits [50,58]. WGCNA is an effective means for identifying co-expressed genes and for analyzing the relationship between the modules and phenotypes to obtain the key candidate genes. In this study, we performed WGCNA on genes expressing FPKM > 1 at the seedling and tuber formation stages, divided them into different modules, and correlated the modules with phenotypes. A total of 14 expression modules were obtained at both the seedling and tuber formation stages, and, interestingly, a higher proportion of DEGs were also present among the modules with a higher phenotypic correlation, indicating that these modules were associated with a low-nitrogen response. We designated these modules as low-nitrogen-response modules. Then, GO and KEGG analyses were conducted, and we found that plant hormone signal transduction and photosynthesis were enriched, as was also found in a previous study [2]. We also combined WGCNA and DEGs to obtain candidate genes with a low-nitrogen response. There were nine and four co-DEGs between DN322 and DN314 (Figure 3E) at the seedling and tuber formation stage, respectively. Then, co-DEGs and low-nitrogen-response modules were combined to obtain one and three genes associated with low-nitrogen stress; these genes will be the focus of future research.

5. Conclusions

At the seedling stage, low-nitrogen stress mainly affects potato growth by affecting the photosynthetic pathway and the steroid biosynthesis pathway. The tuber formation period mainly affects the photosynthetic pathway and the starch and sucrose metabolism pathways. In this study, four genes associated with low-nitrogen stress were identified through WGCNA combined with the study of differentially expressed genes (Figure 6), which can provide theoretical guidance for the high nitrogen efficiency breeding of potatoes.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/agronomy13082164/s1. Figure S1. Effects of nitrogen-deficiency stress on the potato. Figure S2. The RT-qPCR verification of DN322 and DN314 and DEGs analysis. Table S1. Nitrogen-deficiency stress phenotype. Table S2. RNA-seq quality of DN322. Table S3. RNA-seq quality of DN314. Table S4. The DEGs of DN322 between N0 and N1 at the seedling stage. Table S5. The DEGs of DN322 between N0 and N1 at the tuber formation stage. Table S6. The DEGs of DN314 between N0 and N1 at the seedling stage. Table S7. The DEGs of DN314 between N0 and N1 at the tuber formation stage. Table S8. Proportion of DEGs in each WGCNA module at the seedling stage. Table S9. Proportion of DEGs in each WGCNA module at the tuber formation stage.

Author Contributions

Z.G. and Q.W. designed this study; Q.W. and B.D. determined the dry matter weight and physiological indices; B.D. and X.S. took the sequencing samples; Y.Y. analyzed the data and wrote the paper; and Y.S. provided the potato materials. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the earmarked fund for the China Agriculture Research System (CARS-09).

Data Availability Statement

All data related to this manuscript can be found within this paper and its Supplementary data.

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as potential conflict of interest.

References

  1. Tiwari, J.K.; Buckseth, T.; Devi, S.; Varshney, S.; Sahu, S.; Patil, V.U.; Zinta, R.; Ali, N.; Moudgil, V.; Singh, R.K.; et al. Physiological and genome-wide RNA-sequencing analyses identify candidate genes in a nitrogen-use efficient potato cv. Kufri Gaurav. Plant Physiol. Biochem. 2020, 154, 171–183. [Google Scholar] [CrossRef] [PubMed]
  2. Guo, H.; Pu, X.; Jia, H.; Zhou, Y.; Ye, G.; Yang, Y.; Na, T.; Wang, J. Transcriptome analysis reveals multiple effects of nitrogen accumulation and metabolism in the roots, shoots, and leaves of potato (Solanum tuberosum L.). BMC Plant Biol. 2022, 22, 282. [Google Scholar] [CrossRef] [PubMed]
  3. Garnett, T.; Plett, D.; Heuer, S.; Okamoto, M. Genetic approaches to enhancing nitrogen-use efficiency (NUE) in cereals: Challenges and future directions. Funct. Plant Biol. FPB 2015, 42, 921–941. [Google Scholar] [CrossRef] [PubMed]
  4. Bundy, L.G.; Andraski, T.W. Recovery of Fertilizer Nitrogen in Crop Residues and Cover Crops on an Irrigated Sandy Soil. Soil. Sci. Soc. Am. J. 2005, 69, 640–648. [Google Scholar] [CrossRef]
  5. Gao, X.; Li, C.; Zhang, M.; Wang, R.; Chen, B. Controlled release urea improved the nitrogen use efficiency, yield and quality of potato (Solanum tuberosum L.) on silt loamy soil. Field Crop. Res. 2015, 181, 60–68. [Google Scholar] [CrossRef]
  6. Shoji, S.; Delgado, J.; Mosier, A.R.; Miura, Y. Use of controlled release fertilizers and nitrification inhibitors to increase nitrogen use efficiency and to conserve air andwater quality. Commun. Soil Sci. Plant Anal. 2001, 32, 1051–1070. [Google Scholar] [CrossRef]
  7. Galvez, J.H.; Tai, H.H.; Lague, M.; Zebarth, B.J.; Stromvik, M.V. The nitrogen responsive transcriptome in potato (Solanum tuberosum L.) reveals significant gene regulatory motifs. Sci. Rep. 2016, 6, 26090. [Google Scholar] [CrossRef]
  8. Dungait, J.A.; Cardenas, L.M.; Blackwell, M.S.; Wu, L.; Withers, P.J.; Chadwick, D.R.; Bol, R.; Murray, P.J.; Macdonald, A.J.; Whitmore, A.P.; et al. Advances in the understanding of nutrient dynamics and management in UK agriculture. Sci. Total Environ. 2012, 434, 39–50. [Google Scholar] [CrossRef]
  9. Naqqash, T.; Malik, K.A.; Imran, A.; Hameed, S.; Shahid, M.; Hanif, M.K.; Majeed, A.; Iqbal, M.J.; Qaisrani, M.M.; van Elsas, J.D. Inoculation with Azospirillum spp. Acts as the Liming Source for Improving Growth and Nitrogen Use Efficiency of Potato. Front. Plant Sci. 2022, 13, 929114. [Google Scholar] [CrossRef]
  10. Li, W.; Xiong, B.; Wang, S.; Deng, X.; Yin, L.; Li, H. Regulation Effects of Water and Nitrogen on the Source-Sink Relationship in Potato during the Tuber Bulking Stage. PLoS ONE 2016, 11, e0146877. [Google Scholar] [CrossRef]
  11. Francis Zvomuya, C.J.R. Nitrate Leaching and Nitrogen Recovery Following Application of Polyolefin-Coated Urea to Potato. Ground Water Qual. 2003, 32, 480–489. [Google Scholar]
  12. Zebarth, B.J.; Rosen, C.J. Research Perspective on Nitrogen BMP Development for Potato. Am. J. Potato Res. 2007, 84, 3–18. [Google Scholar] [CrossRef]
  13. Hu, B.; Wang, W.; Ou, S.; Tang, J.; Li, H.; Che, R.; Zhang, Z.; Chai, X.; Wang, H.; Wang, Y.; et al. Variation in NRT1.1B contributes to nitrate-use divergence between rice subspecies. Nat. Genet. 2015, 47, 834–838. [Google Scholar] [CrossRef] [PubMed]
  14. Chen, J.; Zhang, Y.; Tan, Y.; Zhang, M.; Zhu, L.; Xu, G.; Fan, X. Agronomic nitrogen-use efficiency of rice can be increased by driving OsNRT2.1 expression with the OsNAR2.1 promoter. Plant Biotechnol. J. 2016, 14, 1705–1715. [Google Scholar] [CrossRef] [PubMed]
  15. Chen, J.; Fan, X.; Qian, K.; Zhang, Y.; Song, M.; Liu, Y.; Xu, G.; Fan, X. pOsNAR2.1:OsNAR2.1 expression enhances nitrogen uptake efficiency and grain yield in transgenic rice plants. Plant Biotechnol. J. 2017, 15, 1273–1283. [Google Scholar] [CrossRef]
  16. Chen, J.; Liu, X.; Liu, S.; Fan, X.; Zhao, L.; Song, M.; Fan, X.; Xu, G. Co-Overexpression of OsNAR2.1 and OsNRT2.3a Increased Agronomic Nitrogen Use Efficiency in Transgenic Rice Plants. Front. Plant Sci. 2020, 11, 1245. [Google Scholar] [CrossRef]
  17. Fan, X.; Tang, Z.; Tan, Y.; Zhang, Y.; Luo, B.; Yang, M.; Lian, X.; Shen, Q.; Miller, A.J.; Xu, G. Overexpression of a pH-sensitive nitrate transporter in rice increases crop yields. Proc. Natl. Acad. Sci. USA 2016, 113, 7118–7123. [Google Scholar] [CrossRef]
  18. Wang, Q.; Su, Q.; Nian, J.; Zhang, J.; Guo, M.; Dong, G.; Hu, J.; Wang, R.; Wei, C.; Li, G.; et al. The Ghd7 transcription factor represses ARE1 expression to enhance nitrogen utilization and grain yield in rice. Mol. Plant 2021, 14, 1012–1023. [Google Scholar] [CrossRef]
  19. Wei, S.; Li, X.; Lu, Z.; Zhang, H.; Ye, X.; Zhou, Y.; Li, J.; Yan, Y.; Pei, H.; Duan, F.; et al. A transcriptional regulator that boosts grain yields and shortens the growth duration of rice. Science 2022, 377, eabi8455. [Google Scholar] [CrossRef]
  20. Huang, Y.; Wang, H.; Zhu, Y.; Huang, X.; Li, S.; Wu, X.; Zhao, Y.; Bao, Z.; Qin, L.; Jin, Y.; et al. THP9 enhances seed protein content and nitrogen-use efficiency in maize. Nature 2022, 612, 292–300. [Google Scholar] [CrossRef]
  21. Li, B.; Xin, W.; Sun, S.; Shen, Q.; Xu, G. Physiological and Molecular Responses of Nitrogen-starved Rice Plants to Re-supply of Different Nitrogen Sources. Plant Soil. 2006, 287, 145–159. [Google Scholar] [CrossRef]
  22. Basu, P.S.; Sharma, A.; Sukumaran, N.P. Changes in net photosynthetic rate and chlorophyll fluorescence in potato leaves induced by water stress. Photosynthetica 1998, 35, 13–19. [Google Scholar] [CrossRef]
  23. Zhang, J.; Wang, Y.; Zhao, Y.; Zhang, Y.; Zhang, J.; Ma, H.; Han, Y. Transcriptome analysis reveals Nitrogen deficiency induced alterations in leaf and root of three cultivars of potato (Solanum tuberosum L.). PLoS ONE 2020, 15, e0240662. [Google Scholar] [CrossRef] [PubMed]
  24. Wang, X.; Wu, F.; Xie, Q.; Wang, H.; Wang, Y.; Yue, Y.; Gahura, O.; Ma, S.; Liu, L.; Cao, Y.; et al. SKIP is a component of the spliceosome linking alternative splicing and the circadian clock in Arabidopsis. Plant Cell 2012, 24, 3278–3295. [Google Scholar] [CrossRef]
  25. Cui, Z.; Tong, A.; Huo, Y.; Yan, Z.; Yang, W.; Yang, X.; Wang, X.X. SKIP controls flowering time via the alternative splicing of SEF pre-mRNA in Arabidopsis. BMC Biol. 2017, 15, 80. [Google Scholar] [CrossRef]
  26. Pandey, A.; Khan, M.K.; Hamurcu, M.; Brestic, M.; Topal, A.; Gezgin, S. Insight into the Root Transcriptome of a Boron-Tolerant Triticum zhukovskyi Genotype Grown under Boron Toxicity. Agronomy 2022, 12, 2421. [Google Scholar] [CrossRef]
  27. Ye, M.; Peng, Z.; Tang, D.; Yang, Z.; Li, D.; Xu, Y.; Zhang, C.; Huang, S. Generation of self-compatible diploid potato by knockout of S-RNase. Nat. Plants 2018, 4, 651–654. [Google Scholar] [CrossRef]
  28. Putri, G.H.; Anders, S.; Pyl, P.T.; Pimanda, J.E.; Zanini, F. Analysing high-throughput sequencing data in Python with HTSeq 2.0. Bioinformatics 2022, 38, 2943–2945. [Google Scholar] [CrossRef]
  29. Li, B.; Dewey, C.N. RSEM: Accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 2011, 12, 323. [Google Scholar] [CrossRef]
  30. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef]
  31. Young, M.D.; Wakefield, M.J.; Smyth, G.K.; Oshlack, A. Gene ontology analysis for RNA-seq: Accounting for selection bias. Genome Biol. 2010, 11, R14. [Google Scholar] [CrossRef]
  32. Mao, X.; Cai, T.; Olyarchuk, J.G.; Wei, L. Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics 2005, 21, 3787–3793. [Google Scholar] [CrossRef] [PubMed]
  33. Kanehisa, M.; Araki, M.; Goto, S.; Hattori, M.; Hirakawa, M.; Itoh, M.; Katayama, T.; Kawashima, S.; Okuda, S.; Tokimatsu, T.; et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008, 36, D480–D484. [Google Scholar] [CrossRef] [PubMed]
  34. Yang, L.; Ma, T.J.; Zhang, Y.B.; Wang, H.; An, R.H. Construction and Analysis of lncRNA-miRNA-mRNA ceRNA Network Identify an Eight-Gene Signature as a Potential Prognostic Factor in Kidney Renal Papillary Cell Carcinoma (KIRP). Altern. Ther. Health Med. 2022, 28, 42–51. [Google Scholar] [PubMed]
  35. Zhang, X.; Huang, N.; Mo, L.; Lv, M.; Gao, Y.; Wang, J.; Liu, C.; Yin, S.; Zhou, J.; Xiao, N.; et al. Global Transcriptome and Co-Expression Network Analysis Reveal Contrasting Response of Japonica and Indica Rice Cultivar to γ Radiation. Int. J. Mol. Sci. 2019, 20, 4358. [Google Scholar] [CrossRef] [PubMed]
  36. Ruiz Herrera, L.F.; Shane, M.W.; López-Bucio, J. Nutritional regulation of root development. Dev. Biol. 2015, 4, 431–443. [Google Scholar] [CrossRef]
  37. Xin, W.; Zhang, L.; Zhang, W.; Gao, J.; Yi, J.; Zhen, X.; Du, M.; Zhao, Y.; Chen, L. An Integrated Analysis of the Rice Transcriptome and Metabolome Reveals Root Growth Regulation Mechanisms in Response to Nitrogen Availability. Int. J. Mol. Sci. 2019, 20, 2349. [Google Scholar] [CrossRef]
  38. Walch-Liu, P.; Ivanov, I.I.; Filleur, S.; Gan, Y.; Remans, T.; Forde, B.G. Nitrogen regulation of root branching. Ann. Bot. 2006, 97, 875–881. [Google Scholar] [CrossRef]
  39. Mu, X.; Chen, Q.; Chen, F.; Yuan, L.; Mi, G. A RNA-Seq Analysis of the Response of Photosynthetic System to Low Nitrogen Supply in Maize Leaf. Int. J. Mol. Sci. 2017, 18, 2624. [Google Scholar] [CrossRef]
  40. Maxwell, K.; Johnson, G.N. Chlorophyll fluorescence—A practical guide. J. Exp. Bot. 2000, 51, 659–668. [Google Scholar] [CrossRef]
  41. Ghannoum, O.; Evans, J.R.; Chow, W.S.; Andrews, T.J.; Conroy, J.P.; von Caemmerer, S. Faster Rubisco is the key to superior nitrogen-use efficiency in NADP-malic enzyme relative to NAD-malic enzyme C4 grasses. Plant Physiol. 2005, 137, 638–650. [Google Scholar] [CrossRef]
  42. Ding, L.; Wang, K.J.; Jiang, G.M.; Biswas, D.K.; Xu, H.; Li, L.F.; Li, Y.H. Effects of nitrogen deficiency on photosynthetic traits of maize hybrids released in different years. Ann. Bot. 2005, 96, 925–930. [Google Scholar] [CrossRef] [PubMed]
  43. Du, A.; Tian, W.; Wei, M.; Yan, W.; He, H.; Zhou, D.; Huang, X.; Li, S.; Ouyang, X. The DTH8-Hd1 Module Mediates Day-Length-Dependent Regulation of Rice Flowering. Mol. Plant 2017, 10, 948–961. [Google Scholar] [CrossRef]
  44. Xue, W.; Xing, Y.; Weng, X.; Zhao, Y.; Tang, W.; Wang, L.; Zhou, H.; Yu, S.; Xu, C.; Li, X.; et al. Natural variation in Ghd7 is an important regulator of heading date and yield potential in rice. Nat. Genet. 2008, 40, 761–767. [Google Scholar] [CrossRef] [PubMed]
  45. Han, X.; Qin, Y.; Sandrine, A.M.N.; Qiu, F. Fine mapping of qKRN8, a QTL for maize kernel row number, and prediction of the candidate gene. Theor. Appl. Genet. 2020, 133, 3139–3150. [Google Scholar] [CrossRef]
  46. Lu, S.; Zhao, X.; Hu, Y.; Liu, S.; Nan, H.; Li, X.; Fang, C.; Cao, D.; Shi, X.; Kong, L.; et al. Natural variation at the soybean J locus improves adaptation to the tropics and enhances yield. Nat. Genet. 2017, 49, 773–779. [Google Scholar] [CrossRef] [PubMed]
  47. Sun, J.; Zhang, G.; Cui, Z.; Kong, X.; Yu, X.; Gui, R.; Han, Y.; Li, Z.; Lang, H.; Hua, Y.; et al. Regain flood adaptation in rice through a 14-3-3 protein OsGF14h. Nat. Commun. 2022, 13, 5664. [Google Scholar] [CrossRef]
  48. Duan, Z.; Zhang, M.; Zhang, Z.; Liang, S.; Fan, L.; Yang, X.; Yuan, Y.; Pan, Y.; Zhou, G.; Liu, S.; et al. Natural allelic variation of GmST05 controlling seed size and quality in soybean. Plant Biotechnol. J. 2022, 20, 1807–1818. [Google Scholar] [CrossRef]
  49. Kao, T.H.; McCubbin, A.G. How flowering plants discriminate between self and non-self pollen to prevent inbreeding. Proc. Natl. Acad. Sci. USA 1996, 93, 12059–12065. [Google Scholar] [CrossRef]
  50. Tiwari, J.K.; Buckseth, T.; Zinta, R.; Saraswati, A.; Singh, R.K.; Rawat, S.; Dua, V.K.; Chakrabarti, S.K. Transcriptome analysis of potato shoots, roots and stolons under nitrogen stress. Sci. Rep. 2020, 10, 1152. [Google Scholar] [CrossRef] [PubMed]
  51. Zhang, X.; Huang, G.; Bian, X.; Zhao, Q. Effects of root interaction and nitrogen fertilization on the chlorophyll content, root activity, photosynthetic characteristics of intercropped soybean and microbial quantity in the rhizosphere. Plant Soil Environ. 2013, 59, 80–88. [Google Scholar] [CrossRef]
  52. Ma, P.; Zhang, X.; Luo, B.; Chen, Z.; Gao, S. Transcriptomic and genome-wide association study reveal long noncoding RNAs responding to nitrogen deficiency in maize. BMC Plant Biol. 2021, 21, 93. [Google Scholar] [CrossRef]
  53. Chen, Z.; Jiang, Q.; Jiang, P.; Zhang, W.; Lu, R. Novel low-nitrogen stressresponsive long non-coding RNAs (lncRNA) in barley landrace B968 (Liuzhutouzidamai) at seedling stage. BMC Plant Biol. 2020, 20, 142. [Google Scholar] [CrossRef] [PubMed]
  54. Wang, C.; Li, Y.; Li, M.; Zhang, K.; Ma, W.; Zheng, L.; Xu, H.; Cui, B.; Liu, R.; Yang, Y.; et al. Functional assembly of root-associated microbial consortia improves nutrient efficiency and yield in soybean. J. Integr. Plant Biol. 2021, 63, 1021–1035. [Google Scholar] [CrossRef] [PubMed]
  55. Esper, B.; Badura, A.; Rögner, M. Photosynthesis as a power supply for (bio-)hydrogen production. Trends Plant Sci. 2006, 11, 543–549. [Google Scholar] [CrossRef] [PubMed]
  56. Goltsev, V.; Zaharieva, I.; Chernev, P.; Strasser, R.J. Delayed fluorescence in photosynthesis. Photosynth. Res. 2009, 101, 217–232. [Google Scholar] [CrossRef]
  57. Nabity, P.D.; Zavala, J.A.; DeLucia, E.H. Indirect suppression of photosynthesis on individual leaves by arthropod herbivory. Ann. Bot. 2009, 103, 655–663. [Google Scholar] [CrossRef] [PubMed]
  58. Sinha, S.K.; Sevanthi, V.A.; Chaudhary, S.; Tyagi, P.; Venkadesan, S.; Rani, M.; Mandal, P.K. Transcriptome Analysis of Two Rice Varieties Contrasting for Nitrogen Use Efficiency under Chronic N Starvation Reveals Differences in Chloroplast and Starch Metabolism-Related Genes. Genes 2018, 9, 206. [Google Scholar] [CrossRef]
Agronomy 13 02164 g001
Figure 1. Nitrogen-deficiency stress affects potato growth and development. (A) Comparison of DN322 and DN314 under low- and normal-nitrogen fertilizer conditions. (B,C) Effects of nitrogen-deficiency stress on the dry matter of DN322 and DN314 at the seedling stage and tuber formation stage, respectively. (D) Effects of nitrogen-deficiency stress on the tuber dry matter of DN322 and DN314 at the tuber formation stage. (E,F) Effects of nitrogen-deficiency stress on the root vitality of DN322 and DN314 at the seedling stage and tuber formation stage, respectively. (G,H) Effects of nitrogen-deficiency stress on the net photosynthetic rates of DN322 and DN314 at the seedling stage and tuber formation stage, respectively. (I) Effects of nitrogen-deficiency stress on the plant height of DN322 and DN314. SS represents the seedling stage, TFS represents the tuber formation stage, and DMW represents the dry matter weight. Data in Figure 1 are means ± SDs. Each black dot in the figure represents the data measured using a single potato plant. A two-tailed unpaired t-test with Welch’s correction was used for statistical analysis. Statistically significant differences (p < 0.01) are indicated with two stars.
Figure 1. Nitrogen-deficiency stress affects potato growth and development. (A) Comparison of DN322 and DN314 under low- and normal-nitrogen fertilizer conditions. (B,C) Effects of nitrogen-deficiency stress on the dry matter of DN322 and DN314 at the seedling stage and tuber formation stage, respectively. (D) Effects of nitrogen-deficiency stress on the tuber dry matter of DN322 and DN314 at the tuber formation stage. (E,F) Effects of nitrogen-deficiency stress on the root vitality of DN322 and DN314 at the seedling stage and tuber formation stage, respectively. (G,H) Effects of nitrogen-deficiency stress on the net photosynthetic rates of DN322 and DN314 at the seedling stage and tuber formation stage, respectively. (I) Effects of nitrogen-deficiency stress on the plant height of DN322 and DN314. SS represents the seedling stage, TFS represents the tuber formation stage, and DMW represents the dry matter weight. Data in Figure 1 are means ± SDs. Each black dot in the figure represents the data measured using a single potato plant. A two-tailed unpaired t-test with Welch’s correction was used for statistical analysis. Statistically significant differences (p < 0.01) are indicated with two stars.
Agronomy 13 02164 g001
Agronomy 13 02164 g002
Figure 2. Transcriptome analysis of DN322 at the seedling and tuber formation stages under nitrogen-deficiency stress. (A) Volcano map of the DEGs upregulated and downregulated between N0 and N1 of DN322 at the seedling stage. (B) GO annotation of the total DEGs of DN322 at the seedling stage. (C) KEGG enrichment analysis of the total DEGs of DN322 at the seedling stage. (D) Volcano map of the DEGs upregulated and downregulated between N0 and N1 of DN322 at the tuber formation stage. (E) GO annotation of the total DEGs of DN322 at the tuber formation stage. (F) KEGG enrichment analysis of the total DEGs of DN322 at the tuber formation stage. Rich factor indicates the ratio of the proportion of DEGs in a pathway to the proportion of genes in all genes annotated to that pathway. The green dots represent the downregulated expression genes, and the red dots represent the upregulated expression genes in A; the larger the enrichment factor, the darker the color, and the more confident the pathway results.
Figure 2. Transcriptome analysis of DN322 at the seedling and tuber formation stages under nitrogen-deficiency stress. (A) Volcano map of the DEGs upregulated and downregulated between N0 and N1 of DN322 at the seedling stage. (B) GO annotation of the total DEGs of DN322 at the seedling stage. (C) KEGG enrichment analysis of the total DEGs of DN322 at the seedling stage. (D) Volcano map of the DEGs upregulated and downregulated between N0 and N1 of DN322 at the tuber formation stage. (E) GO annotation of the total DEGs of DN322 at the tuber formation stage. (F) KEGG enrichment analysis of the total DEGs of DN322 at the tuber formation stage. Rich factor indicates the ratio of the proportion of DEGs in a pathway to the proportion of genes in all genes annotated to that pathway. The green dots represent the downregulated expression genes, and the red dots represent the upregulated expression genes in A; the larger the enrichment factor, the darker the color, and the more confident the pathway results.
Agronomy 13 02164 g002
Agronomy 13 02164 g003
Figure 3. Transcriptome analysis of DN314 at the seedling and tuber formation stages under nitrogen-deficiency stress. (A) Volcano map of DEGs upregulated and downregulated between N0 and N1 of DN314 at the seedling stage. (B) GO annotation of the total DEGs of DN314 at the seedling stage. (C) KEGG enrichment analysis of the total DEGs of DN314 at the seedling stage. (D) Volcano map of the DEGs upregulated and downregulated between N0 and N1 of DN314 at the tuber formation stage. (E) GO annotation of the total DEGs of DN314 at the tuber formation stage. (F) KEGG enrichment analysis of the total DEGs of DN314 at the tuber formation stage.
Figure 3. Transcriptome analysis of DN314 at the seedling and tuber formation stages under nitrogen-deficiency stress. (A) Volcano map of DEGs upregulated and downregulated between N0 and N1 of DN314 at the seedling stage. (B) GO annotation of the total DEGs of DN314 at the seedling stage. (C) KEGG enrichment analysis of the total DEGs of DN314 at the seedling stage. (D) Volcano map of the DEGs upregulated and downregulated between N0 and N1 of DN314 at the tuber formation stage. (E) GO annotation of the total DEGs of DN314 at the tuber formation stage. (F) KEGG enrichment analysis of the total DEGs of DN314 at the tuber formation stage.
Agronomy 13 02164 g003
Agronomy 13 02164 g004
Figure 4. Co-expression network of two potato cultivars at the seedling stage constructed using WGCNA. (A) Correlation analysis of each module at the seedling stage. (A) Correlation analysis of each module to each phenotype at the seedling stage. PDW represents plant dry weight, PNC represents plant nitrogen content, and Pho represents the net photosynthetic rate. Red represents a positive correlation, and blue represents a negative correlation. (B) Proportion of DEGs in each module at the seedling stage. (C) GO annotation of the genes in module SM9. (D) KEGG annotation of the genes in module SM9. (E) GO annotation of the genes in module SM10. (F) KEGG annotation of the genes in module SM10.
Figure 4. Co-expression network of two potato cultivars at the seedling stage constructed using WGCNA. (A) Correlation analysis of each module at the seedling stage. (A) Correlation analysis of each module to each phenotype at the seedling stage. PDW represents plant dry weight, PNC represents plant nitrogen content, and Pho represents the net photosynthetic rate. Red represents a positive correlation, and blue represents a negative correlation. (B) Proportion of DEGs in each module at the seedling stage. (C) GO annotation of the genes in module SM9. (D) KEGG annotation of the genes in module SM9. (E) GO annotation of the genes in module SM10. (F) KEGG annotation of the genes in module SM10.
Agronomy 13 02164 g004
Agronomy 13 02164 g005
Figure 5. Co-expression network of two potato cultivars under the tuber formation stage constructed using WGCNA. (A) Overview of the co-expression network of two potato cultivars in the tuber formation stage. (B) Correlation analysis of each module in the tuber formation stage. The redder the color between modules, the higher the correlation between modules. (C) Correlation analysis of each module to each phenotype in the tuber formation stage. PDW represents plant dry weight, PNC represents plant nitrogen content, and Pho represents the net photosynthetic rate. Red represents a positive correlation, and blue represents a negative correlation. (D) Proportion of differentially expressed genes in each module in the tuber formation stage. (E) GO annotation of the genes in module TM1. (F) KEGG annotation of the genes in module TM1.
Figure 5. Co-expression network of two potato cultivars under the tuber formation stage constructed using WGCNA. (A) Overview of the co-expression network of two potato cultivars in the tuber formation stage. (B) Correlation analysis of each module in the tuber formation stage. The redder the color between modules, the higher the correlation between modules. (C) Correlation analysis of each module to each phenotype in the tuber formation stage. PDW represents plant dry weight, PNC represents plant nitrogen content, and Pho represents the net photosynthetic rate. Red represents a positive correlation, and blue represents a negative correlation. (D) Proportion of differentially expressed genes in each module in the tuber formation stage. (E) GO annotation of the genes in module TM1. (F) KEGG annotation of the genes in module TM1.
Agronomy 13 02164 g005
Agronomy 13 02164 g006
Figure 6. The model employed in this study. On the left is the seedling stage: When nitrogen-deficiency stress occurs, the DEGs of DN322 are mainly enriched via the steroid biosynthesis pathway and diterpenoid biosynthesis, whereas DN314 DEGs are mainly enriched via the photosynthesis-antenna proteins and protein processing in the endoplasmic reticulum pathway. The gene Soltu.DM.05G003170 was identified in response to nitrogen-deficiency stress through combined analysis of WGCNA and DEGs. On the right is the tuber formation stage: When nitrogen-deficiency stress occurs, the DEGs of DN322 and DN314 are mainly enriched via the photosynthesis pathway. Two genes, Soltu.DM.04G009820 and Soltu.DM.04G016890, were identified in response to nitrogen-deficiency stress through combined analysis of WGCNA and DEGs.
Figure 6. The model employed in this study. On the left is the seedling stage: When nitrogen-deficiency stress occurs, the DEGs of DN322 are mainly enriched via the steroid biosynthesis pathway and diterpenoid biosynthesis, whereas DN314 DEGs are mainly enriched via the photosynthesis-antenna proteins and protein processing in the endoplasmic reticulum pathway. The gene Soltu.DM.05G003170 was identified in response to nitrogen-deficiency stress through combined analysis of WGCNA and DEGs. On the right is the tuber formation stage: When nitrogen-deficiency stress occurs, the DEGs of DN322 and DN314 are mainly enriched via the photosynthesis pathway. Two genes, Soltu.DM.04G009820 and Soltu.DM.04G016890, were identified in response to nitrogen-deficiency stress through combined analysis of WGCNA and DEGs.
Agronomy 13 02164 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wei, Q.; Yin, Y.; Deng, B.; Song, X.; Gong, Z.; Shi, Y. Transcriptome Analysis of Nitrogen-Deficiency-Responsive Genes in Two Potato Cultivars. Agronomy 2023, 13, 2164. https://doi.org/10.3390/agronomy13082164

AMA Style

Wei Q, Yin Y, Deng B, Song X, Gong Z, Shi Y. Transcriptome Analysis of Nitrogen-Deficiency-Responsive Genes in Two Potato Cultivars. Agronomy. 2023; 13(8):2164. https://doi.org/10.3390/agronomy13082164

Chicago/Turabian Style

Wei, Qiaorong, Yanbin Yin, Bin Deng, Xuewei Song, Zhenping Gong, and Ying Shi. 2023. "Transcriptome Analysis of Nitrogen-Deficiency-Responsive Genes in Two Potato Cultivars" Agronomy 13, no. 8: 2164. https://doi.org/10.3390/agronomy13082164

APA Style

Wei, Q., Yin, Y., Deng, B., Song, X., Gong, Z., & Shi, Y. (2023). Transcriptome Analysis of Nitrogen-Deficiency-Responsive Genes in Two Potato Cultivars. Agronomy, 13(8), 2164. https://doi.org/10.3390/agronomy13082164

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