The nitrogen responsive transcriptome in potato (Solanum tuberosum L.) reveals significant gene regulatory motifs

Agro-physiological and transcriptome profiling reveal key genes associated with potato tuberization under different nitrogen regimes in aeroponics

Nitrogen (N) is a crucial nutrient for the growth and development of potatoes. However, excessive use of nitrogen fertilizers can have detrimental effects on human health, aquatic ecosystems, and the environment. Therefore, understanding the genes involved in nitrogen metabolism is essential for developing future strategies to improve nitrogen use efficiency (NUE) in plants. This study aimed to identify genes associated with high tuber yield in two contrasting potato varieties Kufri Jyoti (N inefficient) and Kufri Pukhraj (N efficient) grown under low and high nitrogen regimes using an aeroponics system. Both varieties were grown in aeroponics with two nitrogen doses (low N: 0.5 mM N; high N: 5 mM N) using a completely randomized design (CRD) with three replications over two years. The phenotypic results confirmed that Kufri Pukhraj was more nitrogen use efficient compared to Kufri Jyoti, particularly under low nitrogen conditions. Additionally, transcriptome analysis produced high-quality data ( ≥ Q20), ranging from 4.35 to 5.46 Gb per sample. Statistically significant genes (p ≤  0.05) were identified based on the reference potato genome. Differentially expressed genes (DEGs) were categorized as either up-regulated or down-regulated in leaf and tuber tissues. Transcriptome profiling of both tuber and leaf tissues revealed genes associated with traits contributing to high tuber yield under both high and low nitrogen conditions. The DEGs were further characterized through gene ontology (GO) annotation and KEGG pathway analysis. Selected genes were validated through real-time quantitative polymerase chain reaction (RT-qPCR) analysis. In summary, several genes were identified as being involved in high tuber yield component traits in potatoes under different nitrogen conditions. These included glutaredoxin, transcription factors (BTB/POZ, AP2/ERF, and MYB), nitrate transporter, aquaporin TIP1;3, glutamine synthetase, aminotransferase, GDSL esterase/lipase, sucrose synthase, UDP-glycosyltransferases, osmotin, xyloglucan endotransglucosylase/hydrolase, and laccases. Additionally, we identified overexpressed genes including cysteine protease inhibitor 1, miraculin, sterol desaturase, and pectinesterase in Kufri Pukhraj under low N stress. Our study highlights these genes' roles in enhancing tuber yield in potatoes cultivated under both high and low nitrogen in aeroponics.

Transcriptomics combined with physiological analysis and metabolomics revealed the response of potato tuber formation to nitrogen

The absorption of the essential element nitrogen by plants affects various aspects of plant physiological activity, including gene expression, metabolite content and growth. However, the molecular mechanism underlying the potato tuberization response to nitrogen remains unclear. Potato plants were subjected to pot experiments under nitrogen deficiency, normal nitrogen levels and nitrogen sufficiency. A comprehensive analysis of the physiological responses, transcriptomic profiles, and metabolic pathways of potato stolons subjected to nitrogen stress was conducted. Transcriptomic analysis revealed 2756 differentially expressed genes (DEGs) associated with nitrogen stress. Metabolomic analysis identified a total of 600 differentially accumulated metabolites (DAMs). Further correlation analysis of the major DEGs and DAMs revealed that 9 key DEGs were associated with alpha-linolenic acid metabolism, 16 key DEGs with starch and sucrose metabolism, 7 key DEGs with nitrogen metabolism, and 16 key DEGs with ABC transporters. Nitrogen deficiency significantly increased the sucrose, GDP-glucose and L-glutamic acid levels and promoted stolon growth by increasing the expression of AMY (alpha-amylase), BE (1,4-alpha-glucan branching enzyme), SS (starch synthase), SPS (sucrose‒phosphate synthase) and AGPS (glucose‒1-phosphate adenylyltransferase). However, high nitrogen levels had the opposite effect. In addition, high nitrogen levels upregulated EG (endoglucanase), SUS (sucrose synthase) and GDH (glutamate dehydrogenase) and led to significant accumulation of 9-Hydroperoxy-10,12,15-octadecatrienoate (9(S)-HpOTrE), (13 S)-Hydroperoxyoctadeca-9,11,15-trienoate (13 (S)-HpOTrE) and L-glutamine, ultimately affecting the balance between plant growth and defense. Overall, our comprehensive study revealed the co-expressed genes and potential pathways related to potato tuber formation under different nitrogen conditions. These data provide a better understanding needed for improving potato tuber traits at the molecular and metabolic levels.

Multi-omics analysis of excessive nitrogen fertilizer application: Assessing environmental damage and solutions in potato farming

Potatoes (Solanum tuberosum L.) are the third largest food crop globally and are pivotal for global food security. Widespread N fertilizer waste in potato cultivation has caused diverse environmental issues. This study employed microbial metagenomic sequencing to analyze the causes behind the declining N use efficiency (NUE) and escalating greenhouse gas emissions resulting from excessive N fertilizer application. Addressing N fertilizer inefficiency through breeding has emerged as a viable solution for mitigating overuse in potato cultivation. In this study, transcriptome and metabolome analyses were applied to identify N fertilizer-responsive genes. Metagenomic sequencing revealed that excessive N fertilizer application triggered alterations in the population dynamics of 11 major bacterial phyla, consequently affecting soil microbial functions, particularly N metabolism pathways and bacterial secretion systems. Notably, the enzyme levels associated with NO3- increased, and those associated with NO and N2O increased. Furthermore, excessive N fertilizer application enhanced soil virulence factors and increased potato susceptibility to diseases. Transcriptome and metabolome sequencing revealed significant impacts of excessive N fertilizer use on lipid and amino acid metabolism pathways. Weighted gene co‑expression network analysis (WGCNA) was adopted to identify two genes associated with N fertilizer response: PGSC0003DMG400021157 and PGSC0003DMG400009544.

Integrated morphological, physiological and transcriptomic analyses reveal the responses of Toona sinensis seedlings to low-nitrogen stress

Nitrogen is one of the most essential elements for plant growth and development. In this study, the growth, physiology, and transcriptome of Toona sinensis (A. Juss) Roem seedlings were compared between low-nitrogen (LN) and normal-nitrogen (NN) conditions. These results indicate that LN stress adversely influences T. sinensis seedling growth. The activities of key enzymes related to nitrogen assimilation and phytohormone contents were altered by LN stress. A total of 2828 differentially expressed genes (DEGs) in roots and 1547 in leaves were identified between the LN and NN treatments. A differential enrichment analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways indicated that nitrogen and sugar metabolism, flavonoid biosynthesis, plant hormone signal transduction, and ABC transporters, were strongly affected by LN stress. In summary, this research provides information for further understanding the response of T. sinensis to LN stress.

Transcriptome Analysis of Two Tetraploid Potato Varieties under Water-Stress Conditions

Potato (Solanum tuberosum L.) is one of the most important crops worldwide, but due to its sensitivity to drought, its production can be affected by water availability. In this study, the varieties Agria and Zorba were used to determine the expression differences between control and water-stressed plants. For this purpose, they were sequenced by RNAseq, obtaining around 50 million transcripts for each variety and treatment. When comparing the significant transcripts obtained from control and drought-stressed plants of the Agria variety, we detected 931 genes that were upregulated and 2077 genes that were downregulated under stress conditions. When both treatments were compared in Zorba plants, 735 genes were found to be upregulated and 923 genes were found to be downregulated. Significantly more DEGs were found in the Agria variety, indicating a good stress response of this variety. "Abscisic acid and environmental stress-inducible protein TAS14-like" was the most overexpressed gene under drought conditions in both varieties, but expression differences were also found in numerous transcription factors and heat shock proteins. The principal GO term found was "cellular components", more specifically related to the cell membrane and the cell wall, but other metabolic pathways such as carbohydrate metabolism and osmotic adjustment were also identified. These results provide valuable information related to the molecular mechanisms of tolerance to water stress in order to establish the basis for breeding new, more tolerant varieties.

Construction of drought stress regulation networks in potato based on SMRT and RNA sequencing data

BACKGROUND

Potato (Solanum tuberosum) is the fourth most important food crop in the world and plays an important role in food security. Drought stress has a significantly negative impact on potato growth and production. There are several publications involved drought stress in potato, this research contributes to enrich the knowledge.

RESULTS

In this study, next-generation sequencing (NGS) and single-molecule real-time (SMRT) sequencing technology were used to study the transcription profiles in potato in response to 20%PEG6000 simulates drought stress. The leaves of the variety "Désirée" from in vitro plantlets after drought stress at six time points from 0 to 48 hours were used to perform NGS and SMRT sequencing. According to the sequencing data, a total of 12,798 differentially expressed genes (DEGs) were identified in six time points. The real-time (RT)-PCR results are significantly correlated with the sequencing data, confirming the accuracy of the sequencing data. Gene ontology and KEGG analysis show that these DEGs participate in response to drought stress through galactose metabolism, fatty acid metabolism, plant-pathogen interaction, glutathione metabolism and other pathways. Through the analysis of alternative splicing of 66,888 transcripts, the functional pathways of these transcripts were enriched, and 51,098 transcripts were newly discovered from alternative splicing events and 47,994 transcripts were functionally annotated. Moreover, 3445 lncRNAs were predicted and enrichment analysis of corresponding target genes was also performed. Additionally, Alternative polyadenylation was analyzed by TADIS, and 26,153 poly (A) sites from 13,010 genes were detected in the Iso-Seq data.

CONCLUSION

Our research greatly enhanced potato drought-induced gene annotations and provides transcriptome-wide insights into the molecular basis of potato drought resistance.

Transcriptome analysis reveals multiple effects of nitrogen accumulation and metabolism in the roots, shoots, and leaves of potato (Solanum tuberosum L.)

BACKGROUND

Nitrogen (N) is a major element and fundamental constituent of grain yield. N fertilizer plays an essential role in the roots, shoots, and leaves of crop plants. Here, we obtained two N-sensitive potato cultivars.

RESULTS

The plants were cultivated in the pots using N-deficient and N-sufficient conditions. Crop height, leaf chlorophyll content, dry matter, and N-accumulation significantly decreased under N-deficient conditions. Furthermore, we performed a comprehensive analysis of the phenotype and transcriptome, GO terms, and KEGG pathways. We used WGCNA of co-expressed genes, and 116 differentially expressed hub genes involved in photosynthesis, nitrogen metabolism, and secondary metabolites to generate 23 modules. Among those modules, six NRT gene families, four pigment genes, two auxin-related genes, and two energy-related genes were selected for qRT-PCR validation.

CONCLUSIONS

Overall, our study demonstrates the co-expressed genes and potential pathways associated with N transport and accumulation in potato cultivars' roots, shoots, and leaves under N-deficient conditions. Therefore, this study provides new ideas to conduct further research on improving nitrogen use efficiency in potatoes.

The chloroplast redox-responsive transcriptome of solanaceous plants reveals significant nuclear gene regulatory motifs associated to stress acclimation

Transcriptomes of solanaceous plants expressing a plastid-targeted antioxidant protein were analysed to identify chloroplast redox networks modulating the expression of nuclear genes associated with stress acclimation. Plastid functions depend on the coordinated expression of nuclear genes, many of them associated to developmental and stress response pathways. Plastid-generated signals mediate this coordination via retrograde signaling, which includes sensing of chloroplast redox state and levels of reactive oxygen species (ROS), although it remains a poorly understood process. Chloroplast redox poise and ROS build-up can be modified by recombinant expression of a plastid-targeted antioxidant protein, i.e., cyanobacterial flavodoxin, with the resulting plants displaying increased tolerance to multiple environmental challenges. Here we analysed the transcriptomes of these flavodoxin-expressing plants to study the coordinated transcriptional responses of the nucleus to the chloroplast redox status and ROS levels during normal growth and stress responses (drought or biotic stress) in tobacco and potato, members of the economically important Solanaceae family. We compared their transcriptomes against those from stressed and mutant plants accumulating ROS in different subcellular compartments and found distinct ROS-related imprints modulated by flavodoxin expression and/or stress. By introducing our datasets in a large-scale interaction network, we identified transcriptional factors related to ROS and stress responses potentially involved in flavodoxin-associated signaling. Finally, we discovered identical cis elements in the promoters of many genes that respond to flavodoxin in the same direction as in wild-type plants under stress, suggesting a priming effect of flavodoxin before stress manifestation. The results provide a genome-wide picture illustrating the relevance of chloroplast redox status on biotic and abiotic stress responses and suggest new cis and trans targets to generate stress-tolerant solanaceous crops.

Transcriptome analysis identifies CsNRT genes involved in nitrogen uptake in tea plants, with a major role of CsNRT2.4

Tea trees have a high demand for nitrogen (N) fertilizer to improve the yield and quality of tea. In this research, transcriptome analysis revealed the effect of N starvation and resupply upon N uptake in tea plants. We identified 4098 differentially expressed genes (DEGs) that were significantly enriched in amino acid and N metabolism and were extensively mapped to the tea genome. The CsNRT gene family plays vital roles in the nitrogen uptake of tea plants. The full CDS sequences of CsNRT1.1, CsNRT1.2, CsNRT1.5, CsNRT1.7, CsNRT2.4, CsNRT2.5, CsNRT3.1 and CsNRT3.2 were cloned. One-year-old cutting seedlings of Zhongcha302 (ZC302) were selected for hydroponic culture and were used for gene expression analysis. The seedlings were resupplied with 0.2 and 2 mM N after N starvation. The results of the gene expression under different N treatments and in various tissues indicated that the expression of CsNRT2.4 was highly expressed in tea roots and was greatly induced by N. Overexpressed CsNRT2.4 in transgenic Arabidopsis thaliana increased the root lengths and fresh weights and improved the NO₃^- uptake rate in the Arabidopsis roots at a low NO₃^- level. Thus, we inferred that CsNRT2.4 was a key gene for N uptake in tea plant roots. This study provides new insights into the molecular mechanisms of tea plant responses to N resupply and reveals hub genes for improving nitrogen usage efficiency (NUE) in tea plants.

Integrative Transcriptomic and Proteomic Analysis Reveals an Alternative Molecular Network of Glutamine Synthetase 2 Corresponding to Nitrogen Deficiency in Rice (Oryza sativa L.)

Nitrogen (N) is an essential nutrient for plant growth and development. The root system architecture is a highly regulated morphological system, which is sensitive to the availability of nutrients, such as N. Phenotypic characterization of roots from LY9348 (a rice variety with high nitrogen use efficiency (NUE)) treated with 0.725 mM NH₄NO₃ (1/4N) was remarkable, especially primary root (PR) elongation, which was the highest. A comprehensive analysis was performed for transcriptome and proteome profiling of LY9348 roots between 1/4N and 2.9 mM NH₄NO₃ (1N) treatments. The results indicated 3908 differential expression genes (DEGs; 2569 upregulated and 1339 downregulated) and 411 differential abundance proteins (DAPs; 192 upregulated and 219 downregulated). Among all DAPs in the proteome, glutamine synthetase (GS2), a chloroplastic ammonium assimilation protein, was the most upregulated protein identified. The unexpected concentration of GS2 from the shoot to the root in the 1/4N treatment indicated that the presence of an alternative pathway of N assimilation regulated by GS2 in LY9348 corresponded to the low N signal, which was supported by GS enzyme activity and glutamine/glutamate (Gln/Glu) contents analysis. In addition, N transporters (NRT2.1, NRT2.2, NRT2.3, NRT2.4, NAR2.1, AMT1.3, AMT1.2, and putative AMT3.3) and N assimilators (NR2, GS1;1, GS1;2, GS1;3, NADH-GOGAT2, and AS2) were significantly induced during the long-term N-deficiency response at the transcription level (14 days). Moreover, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis demonstrated that phenylpropanoid biosynthesis and glutathione metabolism were significantly modulated by N deficiency. Notably, many transcription factors and plant hormones were found to participate in root morphological adaptation. In conclusion, our study provides valuable information to further understand the response of rice roots to N-deficiency stress.

Transcriptome analysis reveals Nitrogen deficiency induced alterations in leaf and root of three cultivars of potato (Solanum tuberosum L.)

Nitrogen (N) is a key element for the production of potato. The N uptake efficiency, N use efficiency and increased N utilization efficiency can be decreased by N deficiency treatment. We performed this study to investigate the association between transcriptomic profiles and the efficiencies of N in potato. Potato cultivars "Yanshu 4" (short for Y), "Xiabodi" (cv. Shepody, short for X) and "Chunshu 4" (short for C) were treated with sufficient N fertilizer and deficient N fertilizer. Then, the growth parameters and tuber yield were recorded; the contents of soluble sugar and protein were measured; and the activities of enzymes were detected. Leaf and root transcriptomes were analyzed and differentially expressed genes (DEGs) in response to N deficiency were identified. The results showed that N deficiency decreased the nitrate reductase (NR), glutamine synthetase (GS) and root activity. Most of the DEGs between N-treated and N-deficiency participate the processes of transport, nitrate transport, nitrogen compound transport and N metabolism in C and Y, not in X, indicating the cultivar-dependent response to N deficiency. DEGs like glutamate dehydrogenase (StGDH), glutamine synthetase (StGS) and carbonic anhydrase (StCA) play key roles in these processes mentioned above. DEGs related to N metabolism showed a close relationship with the N utilization efficiency (UTE), but not with N use efficiency (NUE). The Major Facilitator Superfamily (MFS) members, like nitrate transporter 2.4 (StNRT2.4), 2.5 (StNRT2.5) and 2.7 (StNRT2.7), were mainly enriched in the processes associated with response to stresses and defense, indicating that N deficiency induced stresses in all cultivars.

Transcriptome Profiling of Potato (Solanum tuberosum L.) Responses to Root-Knot Nematode (Meloidogyne javanica) Infestation during A Compatible Interaction

Root-knot nematode (RKN) Meloidogyne javanica presents a great challenge to Solanaceae crops, including potato. In this study, we investigated transcriptional responses of potato roots during a compatible interaction with M. javanica. In this respect, differential gene expression of Solanum tuberosum cultivar (cv.) Mondial challenged with M. javanica at 0, 3 and 7 days post-inoculation (dpi) was profiled. In total, 4948 and 4484 genes were detected, respectively, as differentially expressed genes (DEGs) at 3 and 7 dpi. Functional annotation revealed that genes associated with metabolic processes were enriched, suggesting they might have an important role in M. javanica disease development. MapMan analysis revealed down-regulation of genes associated with pathogen perception and signaling suggesting interference with plant immunity system. Notably, delayed activation of pathogenesis-related genes, down-regulation of disease resistance genes, and activation of host antioxidant system contributed to a susceptible response. Nematode infestation suppressed ethylene (ET) and jasmonic acid (JA) signaling pathway hindering JA/ET responsive genes associated with defense. Genes related to cell wall modification were differentially regulated while transport-related genes were up-regulated, facilitating the formation of nematode feeding sites (NFSs). Several families of transcription factors (TFs) were differentially regulated by M. javanica infestation. Suggesting that TFs play an indispensable role in physiological adaptation for successful M. javanica disease development. This genome-wide analysis reveals the molecular regulatory networks in potato roots which are potentially manipulated by M. javanica. Being the first study analyzing transcriptome profiling of M. javanica-diseased potato, it provides unparalleled insight into the mechanism underlying disease development.

Physiological and genome-wide RNA-sequencing analyses identify candidate genes in a nitrogen-use efficient potato cv. Kufri Gaurav

Nitrogen (N) is an important nutrient for plant growth. However, its excess application leads to environmental damage. Hence, improving nitrogen use efficiency (NUE) of plant is one of the plausible options to solve the problems. Aim of this study was to identify candidate genes involved in enhancing NUE in potato cv. Kufri Gaurav (N efficient). Plants were grown in aeroponic with two contrasting N regimes (low N: 0.75 mM, and high N: 7.5 mM). Higher NUE in Kufri Gaurav was observed in low N based on the parameters like NUE, NUpE (N uptake efficiency), NUtE (N utilization efficiency) and AgNUE (agronomic NUE). Further, global gene expression profiles in root, leaf and stolon tissues were analyzed by RNA-sequencing using Ion Proton™ System. Quality data (≥Q20) of 2.04-2.73 Gb per sample were mapped with the potato genome. Statistically significant (P ≤ 0.05) differentially expressed genes (DEGs) were identified such as 176 (up-regulated) and 30 (down-regulated) in leaves, 39 (up-regulated) and 105 (down-regulated) in roots, and 81 (up-regulated) and 694 (down-regulated) in stolons. The gene ontology (GO) terms like metabolic process, cellular process and catalytic activity were predominant. Our RT-qPCR analysis confirmed the gene expression profiles of RNA-seq. Overall, we identified candidate genes associated with improving NUE such as superoxide dismutase, GDSL esterase lipase, probable phosphatase 2C, high affinity nitrate transporters, sugar transporter, proline rich proteins, transcription factors (VQ motif, SPX domain, bHLH) etc. Our findings suggest that these candidate genes probably play crucial roles in enhancing NUE in potato.

Structural genome analysis in cultivated potato taxa

KEY MESSAGE

Twelve potato accessions were selected to represent two principal views on potato taxonomy. The genomes were sequenced and analyzed for structural variation (copy number variation) against three published potato genomes. The common potato (Solanum tuberosum L.) is an important staple crop with a highly heterozygous and complex tetraploid genome. The other taxa of cultivated potato contain varying ploidy levels (2X-5X), and structural variations are common in the genomes of these species, likely contributing to the diversification or agronomic traits during domestication. Increased understanding of the genomes and genomic variation will aid in the exploration of novel agronomic traits. Thus, sequencing data from twelve potato landraces, representing the four ploidy levels, were used to identify structural genomic variation compared to the two currently available reference genomes, a double monoploid potato genome and a diploid inbred clone of S. chacoense. The results of a copy number variation analysis showed that in the majority of the genomes, while the number of deletions is greater than the number of duplications, the number of duplicated genes is greater than the number of deleted ones. Specific regions in the twelve potato genomes have a high density of CNV events. Further, the auxin-induced SAUR genes (involved in abiotic stress), disease resistance genes and the 2-oxoglutarate/Fe(II)-dependent oxygenase superfamily proteins, among others, had increased copy numbers in these sequenced genomes relative to the references.

Identification of Regulatory Networks and Hub Genes Controlling Nitrogen Uptake in Tea Plants [Camellia sinensis (L.) O. Kuntze]

Nitrogen (N) uptake, as the first step of N metabolism, is a key limiting factor for plant growth. To understand the gene expression networks that control N absorption and metabolism in tea plants, we analyzed transcriptomes in the young roots of two groups of tea plants with significantly different growth rates under different N treatments (0, 0.2, and 2 mM). Using pairwise comparisons and weighted gene co-expression network analyses (WGCNA), we successfully constructed 16 co-expression modules. Among them, a specific module (turquoise) that substantially responded to the low N treatment was identified. Based on KEGG analysis, the relative genes that enriched in the "N metabolism" pathways were used to construct gene co-expression networks of N metabolism. Finally, a high-affinity ammonium (NH₄⁺) transporter designated CsAMT1.2 was identified as a hub gene in the N metabolism network in tea plant roots and the gene expression could be highly induced by N resupply. The gene functional analysis revealed that CsAMT1.2 could make functional complementation of MEP1, MEP2, and MEP3 genes in 31019b yeast cells and improve NH₄⁺ uptake rate in 31019b at low NH₄⁺ level. Thus, CsAMT1.2 was a key gene controlling N uptake in tea plants and might play a vital role in promoting NH₄⁺ uptake from the environment in tea roots. This study provided a useful foundation for improving the NUE in tea plantations.

Transcriptome analysis of potato shoots, roots and stolons under nitrogen stress

Potato crop requires high dose of nitrogen (N) to produce high tuber yield. Excessive application of N causes environmental pollution and increases cost of production. Hence, knowledge about genes and regulatory elements is essential to strengthen research on N metabolism in this crop. In this study, we analysed transcriptomes (RNA-seq) in potato tissues (shoot, root and stolon) collected from plants grown in aeroponic culture under controlled conditions with varied N supplies i.e. low N (0.2 milli molar N) and high N (4 milli molar N). High quality data ranging between 3.25 to 4.93 Gb per sample were generated using Illumina NextSeq500 that resulted in 83.60-86.50% mapping of the reads to the reference potato genome. Differentially expressed genes (DEGs) were observed in the tissues based on statistically significance (p ≤ 0.05) and up-regulation with ≥ 2 log2 fold change (FC) and down-regulation with ≤ -2 log2 FC values. In shoots, of total 19730 DEGs, 761 up-regulated and 280 down-regulated significant DEGs were identified. Of total 20736 DEGs in roots, 572 (up-regulated) and 292 (down-regulated) were significant DEGs. In stolons, of total 21494 DEG, 688 and 230 DEGs were significantly up-regulated and down-regulated, respectively. Venn diagram analysis showed tissue specific and common genes. The DEGs were functionally assigned with the GO terms, in which molecular function domain was predominant in all the tissues. Further, DEGs were classified into 24 KEGG pathways, in which 5385, 5572 and 5594 DEGs were annotated in shoots, roots and stolons, respectively. The RT-qPCR analysis validated gene expression of RNA-seq data for selected genes. We identified a few potential DEGs responsive to N deficiency in potato such as glutaredoxin, Myb-like DNA-binding protein, WRKY transcription factor 16 and FLOWERING LOCUS T in shoots; high-affinity nitrate transporter, protein phosphatase-2c, glutaredoxin family protein, malate synthase, CLE7, 2-oxoglutarate-dependent dioxygenase and transcription factor in roots; and glucose-6-phosphate/phosphate translocator 2, BTB/POZ domain-containing protein, F-box family protein and aquaporin TIP1;3 in stolons, and many genes of unknown function. Our study highlights that these potential genes play very crucial roles in N stress tolerance, which could be useful in augmenting research on N metabolism in potato.

SWATH-MS-facilitated proteomic profiling of fruit skin between Fuji apple and a red skin bud sport mutant

BACKGROUND

Apple is one of the most popular fruit crops world-wide and its skin color is an important quality consideration essential for commercial value. However, the strategy on genetic breeding for red skin apple and the genetic basis of skin color differentiation is very limited and still largely unknown.

RESULTS

Here, we reported a bud sport mutant of Fuji apple with red skin color and enhanced anthocyanins accumulation. Quantitative SWATH-MS (sequential window acquisition of all theoretical spectra-mass spectrometry) proteomics investigations revealed proteome changes in the apple red skin bud mutation and a total of 451 differentially expressed proteins were identified in apple skin. The mutant showed significantly increased expression levels of photosynthesis-related proteins, stress-related proteins as well as anthocyanins biosynthesis pathway. On the other hand, substantial downregulation of mitogen-activated protein kinase 4 (MAPK4) and mevalonate kinase (MVK) were detected, indicating a promising role for the red skin color development in the mutant. Furthermore, we also hypothesize that a post-transcriptional regulation of the skin color formation occurs in the mutant through the advanced SWATH-MS analysis.

CONCLUSION

Our work provides important information on the application of proteomic methods for analysing proteomes changes in Fuji apple and highlights a clade of regulatory proteins potentially contributing for the molecular breeding of fruit skin color.

Map-Based Functional Analysis of the GhNLP Genes Reveals Their Roles in Enhancing Tolerance to N-Deficiency in Cotton

Nitrogen is a key macronutrient needed by plants to boost their production, but the development of cotton genotypes through conventional approaches has hit a bottleneck due to the narrow genetic base of the elite cotton cultivars, due to intensive selection and inbreeding. Based on our previous research, in which the BC₂F₂ generations developed from two upland cotton genotypes, an abiotic stress-tolerant genotype, G. tomentosum (donor parent) and a highly-susceptible, and a highly-susceptible, but very productive, G. hirsutum (recurrent parent), were profiled under drought stress conditions. The phenotypic and the genotypic data generated through genotyping by sequencing (GBS) were integrated to map drought-tolerant quantitative trait loci (QTLs). Within the stable QTLs region for the various drought tolerance traits, a nodule-inception-like protein (NLP) gene was identified. We performed a phylogenetic analysis of the NLP proteins, mapped their chromosomal positions, intron-exon structures and conducted ds/dn analysis, which showed that most NLP genes underwent negative or purifying selection. Moreover, the functions of one of the highly upregulated genes, Gh_A05G3286 (Gh NLP5), were evaluated using the virus gene silencing (VIGS) mechanism. A total of 226 proteins encoded by the NLP genes were identified, with 105, 61, and 60 in Gossypium hirsutum, G. raimondii, and G. arboreum, respectively. Comprehensive Insilico analysis revealed that the proteins encoded by the NLP genes had varying molecular weights, protein lengths, isoelectric points (pI), and grand hydropathy values (GRAVY). The GRAVY values ranged from a negative one to zero, showing that proteins were hydrophilic. Moreover, various cis-regulatory elements that are the binding sites for stress-associated transcription factors were found in the promoters of various NLP genes. In addition, many miRNAs were predicted to target NLP genes, notably miR167a, miR167b, miR160, and miR167 that were previously shown to target five NAC genes, including NAC1 and CUC1, under N-limited conditions. The real-time quantitative polymerase chain reaction (RT-qPCR) analysis, revealed that five genes, Gh_D02G2018, Gh_A12G0439, Gh_A03G0493, Gh_A03G1178, and Gh_A05G3286 were significantly upregulated and perhaps could be the key NLP genes regulating plant response under N-limited conditions. Furthermore, the knockdown of the Gh_A05G3286 (GhNLP5) gene by virus-induced silencing (VIGS) significantly reduced the ability of these plants to the knockdown of the Gh_A05G3286 (GhNLP5) gene by virus-induced gene silencing (VIGS) significantly reduced the ability of the VIGS-plants to tolerate N-limited conditions compared to the wild types (WT). The VIGS-plants registered lower chlorophyll content, fresh shoot biomass, and fresh root biomass, addition to higher levels of malondialdehyde (MDA) and significantly reduced levels of proline, and superoxide dismutase (SOD) compared to the WT under N-limited conditions. Subsequently, the expression levels of the Nitrogen-stress responsive genes, GhTap46, GhRPL18A, and GhKLU were shown to be significantly downregulated in VIGS-plants compared to their WT under N-limited conditions. The downregulation of the nitrogen-stress responsive genes provided evidence that the silenced gene had an integral role in enhancing cotton plant tolerance to N-limited conditions.

mRNA transcription profile of potato (Solanum tuberosum L.) exposed to ultrasound during different stages of in vitro plantlet development

KEY MESSAGE

In response to an ultrasound pulse, several hundred DEGs, including in response to stress, were up- or down-regulated in in vitro potato plantlets. Despite this abiotic stress, plantlets survived. Ultrasound (US) can influence plant growth and development. To better understand the genetic mechanism underlying the physiological response of potato to US, single-node segments of four-week-old in vitro plantlets were subjected to US at 35 kHz for 20 min. Following mRNA purification, 10 cDNA libraries were assessed by RNA-seq. Significantly differentially expressed genes (DEGs) were categorized by gene ontology or Kyoto Encyclopedia of Genes and Genomes identifiers. The expression intensity of 40,430 genes was studied. Several hundred DEGs associated with biosynthesis, carbohydrate metabolism and catabolism, cellular protein modification, and response to stress, and which were expressed mainly in the extracellular region, nucleus, and plasma membrane, were either up- or down-regulated in response to US. RT-qPCR was used to validate RNA-seq data of 10 highly up- or down-regulated DEGs, and both Spearman and Pearson correlations between SeqMonk LFC and RT-qPCR LFC were highly positive (0.97). This study examines how some processes evolved over time (0 h, 24 h, 48 h, 1 week and 4 weeks) after an abiotic stress (US) was imposed on in vitro potato explants, and provides clues to the temporal dynamics in DEG-based enzyme functions in response to this stress. Despite this abiotic stress, plantlets survived.

Analysis of Key Genes Involved in Potato Anthocyanin Biosynthesis Based on Genomics and Transcriptomics Data

The accumulation of secondary metabolites, such as anthocyanins, in cells plays an important role in colored plants. The synthesis and accumulation of anthocyanins are regulated by multiple genes, of which the R2R3-MYB transcription factor gene family plays an important role. Based on the genomic data in the Potato Genome Sequencing Consortium database (PGSC) and the transcriptome data in the SRA, this study used potato as a model plant to comprehensively analyze the plant anthocyanin accumulation process. The results indicated that the most critical step in the synthesis of potato anthocyanins was the formation of p-coumaroyl-CoA to enter the flavonoid biosynthetic pathway. The up-regulated expression of the CHS gene and the down-regulated expression of HCT significantly promoted this process. At the same time, the anthocyanins in the potato were gradually synthesized during the process from leaf transport to tubers. New transcripts of stAN1 and PAL were cloned and named stAN1-like and PAL-like, respectively, but the functions of these two new transcripts still need further study. In addition, the sequence characteristics of amino acids in the R2-MYB and R3-MYB domains of potato were preliminarily identified. The aims of this study are to identify the crucial major genes that affect anthocyanin biosynthesis through multi-omics joint analysis and to transform quantitative traits into quality traits, which provides a basis and reference for the regulation of plant anthocyanin biosynthesis. Simultaneously, this study provides the basis for improving the anthocyanin content in potato tubers and the cultivation of new potato varieties with high anthocyanin content.

Transcriptome Profiling and Characterization of Drought-Tolerant Potato Plant (Solanum tuberosum L.)

Potato (Solanum tuberosum L.) is the third most important food crop, and breeding drought-tolerant varieties is vital research goal. However, detailed molecular mechanisms in response to drought stress in potatoes are not well known. In this study, we developed EMS-mutagenized potatoes that showed significant tolerance to drought stress compared to the wild-type (WT) 'Desiree' cultivar. In addition, changes to transcripts as a result of drought stress in WT and drought-tolerant (DR) plants were investigated by de novo assembly using the Illumina platform. One-week-old WT and DR plants were treated with -1.8 Mpa polyethylene glycol-8000, and total RNA was prepared from plants harvested at 0, 6, 12, 24, and 48 h for subsequent RNA sequencing. In total, 61,100 transcripts and 5,118 differentially expressed genes (DEGs) displaying up- or down-regulation were identified in pairwise comparisons of WT and DR plants following drought conditions. Transcriptome profiling showed the number of DEGs with up-regulation and down-regulation at 909, 977, 1181, 1225 and 826 between WT and DR plants at 0, 6, 12, 24, and 48 h, respectively. Results of KEGG enrichment showed that the drought tolerance mechanism of the DR plant can mainly be explained by two aspects, the 'photosynthetic-antenna protein' and 'protein processing of the endoplasmic reticulum'. We also divided eight expression patterns in four pairwise comparisons of DR plants (DR0 vs DR6, DR12, DR24, DR48) under PEG treatment. Our comprehensive transcriptome data will further enhance our understanding of the mechanisms regulating drought tolerance in tetraploid potato cultivars.

Nitrogen Deficiency Induced Alterations in the Root Proteome of a Pair of Potato (Solanum tuberosum L.) Varieties Contrasting for their Response to Low N

Improving crop nitrogen use efficiency is important both from the economic and the environmental viewpoint. Here, the aim is to highlight differences between the proteomic response of the roots of two potato cultivars contrasting in their response to nitrogen (N) deficiency, in an effort to understand which proteins and metabolic pathways contribute to the tolerance of N deprivation. The two cultivars ''Topas'' (tolerant) and ''Lambada'' (sensitive) are grown under both an N sufficient and an N deficient regime, using an in vitro-based cultivation system. Responsive proteins are identified and quantified using label-free quantitative shotgun proteomics. The contrasting cultivars differed with respect to components of the glutamine synthetase/glutamine oxoglutarate aminotransferase pathway, tricarboxylic acid cycle, the glycolysis/gluconeogenesis pathway as well as protein and amino acid synthesis machinery. Additional differences are associated with protein catabolism and defense mechanisms.

Gene expression profiles predictive of cold-induced sweetening in potato

Cold storage (2-4 °C) used in potato production to suppress diseases and sprouting during storage can result in cold-induced sweetening (CIS), where reducing sugars accumulate in tuber tissue leading to undesirable browning, production of bitter flavors, and increased levels of acrylamide with frying. Potato exhibits genetic and environmental variation in resistance to CIS. The current study profiles gene expression in post-harvest tubers before cold storage using transcriptome sequencing and identifies genes whose expression is predictive for CIS. A distance matrix for potato clones based on glucose levels after cold storage was constructed and compared to distance matrices constructed using RNA-seq gene expression data. Congruence between glucose and gene expression distance matrices was tested for each gene. Correlation between glucose and gene expression was also tested. Seventy-three genes were found that had significant p values in the congruence and correlation tests. Twelve genes from the list of 73 genes also had a high correlation between glucose and gene expression as measured by Nanostring nCounter. The gene annotations indicated functions in protein degradation, nematode resistance, auxin transport, and gibberellin response. These 12 genes were used to build models for prediction of CIS using multiple linear regression. Nine linear models were constructed that used different combinations of the 12 genes. An F-box protein, cellulose synthase, and a putative Lax auxin transporter gene were most frequently used. The findings of this study demonstrate the utility of gene expression profiles in predictive diagnostics for severity of CIS.