Skip to main content

Transcriptome analysis of sweet potato responses to potassium deficiency

Abstract

Background

As one of three essential nutrients, potassium is regarded as a main limiting factor for growth and development in plant. Sweet potato (Ipomoea batatas L.) is one of seven major food crops grown worldwide, and is both a nutrient-rich food and a bioenergy crop. It is a typical ‘K-favoring’ crop, and the level of potassium ion (K+) supplementation directly influences its production. However, little is known about the transcriptional changes in sweet potato genes under low-K+ conditions. Here, we analyzed the transcriptomic profiles of sweet potato roots in response to K+ deficiency to determine the effect of low-K+ stress on this economically important crop.

Results

The roots of sweet potato seedlings with or without K+ treatment were harvested and used for transcriptome analyses. The results showed 559 differently expressed genes (DEGs) in low and high K+ groups. Among the DEGs, 336 were upregulated and 223 were downregulated. These DEGs were involved in transcriptional regulation, calcium binding, redox-signaling, biosynthesis, transport, and metabolic process. Further analysis revealed previously unknow genes involved in low-K+ stress, which could be investigated further to improve low K+ tolerance in plants. Confirmation of RNA-sequencing results using qRT-PCR displayed a high level of consistency between the two experiments. Analysis showed that many auxin-, ethylene- and jasmonic acid-related genes respond to K+ deficiency, suggesting that these hormones have important roles in K+ nutrient signaling in sweet potato.

Conclusions

According to the transcriptome data of sweet potato, various DEGs showed transcriptional changes in response to low-K+ stress. However, the expression level of some kinases, transporters, transcription factors (TFs), hormone-related genes, and plant defense-related genes changed significantly, suggesting that they have important roles during K+ deficiency. Thus, this study identifies potential genes for genetic improvement of responses to low-K+ stress and provides valuable insight into the molecular mechanisms regulating low K+ tolerance in sweet potato. Further research is required to clarify the function of these DEGs under low-K+ stress.

Peer Review reports

Background

Potassium ions (K+) are an essential macronutrients for plant growth and development. Unlike nitrate (NO3) and phosphate (H2PO4), it can assimilate into organic matter. K+ is the most abundant cation in higher plants, in which it participates in many fundamental processes, including enzyme activation, membrane potential, membrane transport, and turgor maintenance, as well as determining the yield and quality of crops [1,2,3]. As one of the most important cations, it comprises 2–10% of the plant dry weight and its concentration in fresh tissues is in the range 10–100 mM [4, 5].

To improve their K utilization efficiency (KUE), plants have evolved multiple complex uptake mechanisms that operate at low external and high environmental K+ concentrations [6,7,8]. K+ is absorbed from the soil into plant cells through K+ transport components, such as K+ transporters and channels [9, 10]. AKT1 from the plant Shaker family was the first K+ channel identified and functions at a wide range of external K+ concentrations (approximately 0.01–10 mM); it is considered the main channel mediating K+ influx into root cells [11,12,13]. HAK5 from the HAK/KUP/KT (high-affinity K+/K+ uptake/K+ transporter) family is a high-affinity K+ transporter that contributes approximately half the K+ absorption under K+-deficient conditions in Arabidopsis [14, 15]. AKT1 and HAK5 are the most important transport proteins mediating almost all K+ absorption in Arabidopsis [2, 16, 17].

K+ is involved in multiple abiotic and biotic stresses signal pathways in plants [18,19,20]. As potassium transporters, some HAK/KUP/KT family genes were reported to respond to drought and salt stresses [21, 22]. For example, the HAK/KUP/KT family K+ transport ZmHAK4 has distinct roles in promoting shoot Na+ exclusion and salt tolerance [23]. OsHAK1, a high-affinity K+ transporter, positively regulates response to drought stress in rice [24]. KUP9 is a member of the HAK/KUP/KT family, and controls primary root growth in Arabidopsis [25]. In addition, K+ has important physiological roles in cell division, as well as sugar translocation to fruits via the phloem, fruit coloration, shelf life, and the shipping quality of many horticultural crops [26,27,28]. Low K+ reduces photosynthesis, impairs the supply of photosynthates to sink organs, and promotes sucrose export to phloem from source leaves. This affects not only yield formation, but also quality parameters in, for example, wheat, potato, and grape [29]. In sweet potato (Ipomoea batatas L.), it was reported that K+ application promoted starch accumulation and storage root yield through regulating the activity and transcription of genes involved in sucrose-to-starch conversion [30]. Several studies have shown various effects of different forms of K on fruit yield; for example, the increasing yield of K-treated trees was correlated with an increase in fruit weight [28, 31]. K+ is the major osmotic solute imported into fibers and elongating fiber cells require abundant K to maintain cell turgor pressure, which they acquire via K+ transporter genes that are preferentially expressed in elongating fibers [32, 33]. K+ deficiency reduced cotton fiber strength and fiber length, and also accelerated cotton senescence and reduce production [34,35,36]. K-deficient plants are also easily infected by diseases and pests [37,38,39].

Sweet potato (Ipomoea batatas L.; Convolvulaceae) is a tuber crop and one of the most important crops in China and many other countries [40, 41]. It has 2n = 6x = 90 chromosomes and the haplotype-resolved genome sequence of sweet potato enables the investigation of gene functions [42]. However, the genome of sweet potato is highly heterozygous and generally self-incompatible, which poses numerous challenges for conventional breeding [40]. Sweet potato is a typical ‘K-favoring’ crop, and has high nutritive value. It is a cheap source of energy, carbohydrates, vitamins, K, iron, fiber, and protein [43]; however, little progress has been made in understanding changes in the transcription of sweet potato genes under low-K+ stress. Given that, of the three major nutritional elements (N, P and K), K is the nutrient in highest demand, research is required into the effect of its deficiency on sweet potato.

RNA sequencing (RNA-seq) technologies have become an important tool to identify differentially expressed genes (DEGs) at the transcriptome level and provide valuable information to accelerate studies in plants, especially for crops with complex genomes [44]. These technologies can be used to help understand all aspects of K+ management in plants, including processes related to growth, development, metabolism, and stress resistance [45]. In the present study, second sequencing technology was used to monitor the transcriptome profiles of sweet potato roots in response to K+ deficiency and to establish a useful database for transcriptome sequencing and DEGs under low-K+ stress. These results will further our understanding of the molecular mechanisms involved in, and provide novel insights into, the response of sweet potato to low-K+ stress, as well as the potential genes involved.

Results

Effect of low-K+ stress on the phenotype, K+ content, and relative expression of marker genes in sweet potato seedlings

To identify the phenotype of sweet potato seedlings under K+-deficient conditions, low-K+ stress (LK, 0 mM K+) was applied at the two-leaf seedling stage. After 14 days, the K+-deficient phenotypes were evaluated. Compared with control seedlings (HK, 1 mM K+), K+-deficient seedlings showed smaller shoots and chlorosis (Fig. 1a) and the K+ content was significantly reduced (Fig. 1b). The phenotype of sweet potato roots under low-K+ stress have been shown in Fig. S1a. The roots of sweet potato seedlings showed no obviously phenotype under normal and low-K+ stress and it may need more time to show obvious difference. These observed phenotypes were consistent with biomass and chlorophyll content measurements (Fig. S1b, c and e). In addition, these potassium deficiency phenotypes were consistent with those reported for other crops [45, 46]. In order to identify the potassium deficiency phenotype, the relative expression of the low-K+ stress marker genes was determined to further support the phenotype. CBL-interacting protein kinase 23 (IbCIPK23) and high-affinity K+ transporter 5 (IbHAK5) are two important genes involved in the response to low-K+ stress, and have crucial roles particularly under K+-deficient conditions [15, 47]. Thus genes were selected as controls for qPCR. The results showed that both genes were significantly upregulated in K+ deficient sweet potato roots (Fig. 1c and d). Phenotype and K+ content analyses showed that the sweet potato seedlings suffered low-K+ stress (Fig. 1a and b). Therefore, these materials were shown to be valid for obtaining transcriptome profiles of sweet potato under K+-deficient conditions.

Fig. 1
figure 1

K+ deficiency phenotype and relation expression of low potassium maker gene. a Phenotype comparison between HK (Sufficient potassium) and LK (Low potassium) of sweet potato seedling. Sweet potato stems were cut into 5 cm fragments and growth in plot full with vermiculite to two leaf stage. Seedlings were then treated with HK (1 mM K) and LK (0 mM K) Hoagland solution for 14 days, Bar = 2 cm. b K+ content measurement. Student’s t test (**P < 0.01) was used to analyze statistical significance, n = 4. c and d Relative expression of maker gene IbHAK5 and IbCIPK23 response to low-K+ stress, n = 4, HK-S: sufficient potassium shoot, HK-R: sufficient potassium root, LK-S: low potassium shoot, LK-R: low potassium root

Transcriptome profiling in response to K+ deficiency

To assess gene expression in sweet potato roots under K+-deficient condition, a high-throughput transcriptome sequence analysis was conducted without a genome [44]. To ensure the reliability of the transcriptome data, three biological replicates were analyzed for each treatment. The roots of seedlings of ‘Taizhong6’, which is the variety cultivated in China, were treated under sufficient K+ and low-K+ stress for 14 days and RNA was extracted to obtain transcriptomic profiles. The quality of the sequencing data was sufficient to support further transcriptome analysis by paired-end sequencing using an Illumina NovaSeq6000 sequencing platform. Clean reads were obtained by removing adaptor sequences, unknown sequences, and low-quality reads from raw data. At the same time, Q20, Q30, GC-content, and the sequence duplication level of the clean data were calculated (Table S1). All analyses were based on high-quality clean reads. After quality control, ~ 22,350,759 and ~ 26,153,929 clean reads were obtained under low K+ and high K+ conditions, respectively (Table S1).

Trinity-v2.4.0 software was used to assemble clean reads and obtain unigenes. After analyzing the open-reading frame (ORF) findings, single-nucleotide polymorphisms (SNPs), and simple sequence repeats (SSRs), the genome was mapped for annotation and expression analysis. The quality of the assembly is closely related to the length and number of unigenes. In this study, the unique gene length showed wide range and those longer than 3000 bp accounted for 3.06% of genes (Fig. S2a). Unigenes were annotated in multiple databases, including NR database (the NCBI non-redundant protein sequences database); Swissprot (a manually annotated and non-redundant protein sequence database); Kyoto Encyclopedia of Genes and Genomes (KEGG); and Clusters of Orthologous Groups of proteins (KOG). In total, 91,690 unigenes were obtained (Table S2) and 38,516 unigenes were annotated (Table S3); all unigene sequences, annotation and expression data were collected in Additional files 2 and 5. The Venn diagram showing the annotated genes for each database is detailed in Fig. S2b.

To investigate differences among the transcriptome sequence data, the hierarchical cluster analysis of different expression genes during K+ sufficient and K+ deficient conditions were conducted (Fig. 2a), indicating that there were obvious difference between HK (High K+) and LK (Low K+) condition. Statistical results showed that approximately 336 genes were observably upregulated and approximately 223 were downregulated in the low-K+ treatment (Fig. 2b). The numbers and fold changes of different expression genes under low K+ were clearly displayed in the MA plot, and the red and green dots indicate differently expressed genes (Fig. 2c).

Fig. 2
figure 2

The heat map of DEGs and the number under low K+ stress. a Hierarchical clustering of all differentially expressed genes in the root of sweet potato seedlings, three biological parallels were used in each group. b The numbers of up-regulation and down-regulation genes of DEGs under K+ deficiency condition. c MA plot of the differently expressed genes in HK and LK. The x-axis shows the counts of differently expressed genes, and the y-axis shows fold change value of gene expression in LK. Blue dots indicate genes without significant differential expression. The red dots show significantly up-regulated DEGs, and the green dots show significantly down-regulated DEGs. HK: sufficient K; LK: Low-K+ treat

Transcriptome sequencing identified DEGs in sweet potato seedlings under low-K+ stress

In this study, 559 DEGs were observed under low-K+ stress, some of which were analyzed further (Table 1). These genes were selected based on their p-value and function, including transporters, TFs, cell wall-related genes, disease resistance-related genes, kinases, hormones, E3 ubiquitin ligases and uclacyanin (Table 1). Most of these genes are closely involved in low-K+ stress. For example, K+ and NO3 are reported to be cooperatively involved in K+ uptake [48]. In the current study, high-affinity nitrate transporter 2.4 was upregulated, which implys this nitrate transporter may be involved in K+ homeostasis.

Table 1 List of DEGs from sweet potato roots (LK vs HK) identified in this study

To validate the quality of the gene activity profiles, 11 high DEGs were selected according to their p-values to compare their fragments per kilobase of exon per million fragments mapped (FPKM) reads and qPCR data (Table 2). Real-time PCR analyses further confirmed that the expression of the selected genes and the trends in gene expression changes determined by the two different approaches were largely consistent (Fig. 3). Thus, the DEGs determined in this study can be considered highly accurate.

Table 2 Selected top differently expressed gene in root under low and sufficient potassium condition
Fig. 3
figure 3

Confirmation the selected genes from transcriptomic profiles by real-time PCR. The sweet potato seedlings roots conducted with HK and LK for 14 d and used for qPCR verification. The experiments were repeated three times. The error bars represent ± SE. n = 3, HK: sufficient K; LK: Low-K+ treat; R: Root

Gene ontology analysis of DEGs in sweet potato seedlings in response to K+ deficiency

To evaluate the potential functions of these DEGs in response to K+ deficiency, gene ontology (GO) analysis was performed by mapping each DEG into the records of the GO database with an adjusted p-value < 0.05 as significant enrichment. We chose the top 20 GO terms for further analysis of biological, cellular component and molecular function processes (Fig. 4). In the GO classification of biological processes, cellular catabolic processes and oxidation-reduction process were significantly enriched, indicating that the ‘Taizhong6’ variety has a wide range of catabolic and oxidation-reduction activities under low-K+ stress. Within the ‘cellular component’ classification, intracellular part, cytoplasm, intracellular organelle, and organelle part were prominently represented. Within the ‘molecular function’ classification, the main functional group of DEGs were those related to metal ion binding, protein binding, DNA binding, and nucleic acid binding of TFs (Additional file 3). Most of these processes were closely related to a response to K+ deficiency.

Fig. 4
figure 4

The numbers of DEGs in each GO term was significantly enriched. Functional categorization of genes based on the biological process of gene ontology. Different color indicates an enriched GO category, including BP (biological process); CC (cellular component); MF (molecular function). Each component shows the top 20 of GO term. The x-axis shows the counts of differently expressed genes, and the y-axis shows GO term of gene enriched in each biological process

KEGG analysis of DEGs in sweet potato seedlings in response to K+ deficiency

To identify the pathways in which the DEGs were likely to be involved, we performed the pathway analysis with the KEGG. The top 20 pathways are shown in Fig. 5. The analysis showed that the retinol metabolism, tyrosine metabolism, metabolism of xenobiotics by cytochrome P450, steroid hormone biosynthesis, and other glycan degradations were significantly enriched (Table S5). Previous studies reported that low-K+ stress affects the color of fruits and the quality of crops [26, 28, 31]. The cultivar ‘Taizhong6’ used in this study is a carotenoid-rich sweet potato variety and, usually, this orange-fleshed sweet potato is rich in vitamin A. In the transcriptome data from the current study, the retinol (vitamin A) metabolism pathway was the most significant pathway in the KEGG pathway under K+-deficient condition (Fig. 5 and Additional file 4). Carotenoids are important source for vitamin A. In this study, the total carotenoids were measured (Fig. S1d), the result showed that there were no significantly difference under normal and K+-deficient conditions. The vitamin A of root tuber will have obviously phenotype when K+-deficient time is extended. These annotations provide a valuable resource for investigating the response mechanism of sweet potato under low K+.

Fig. 5
figure 5

The significantly enriched KEGG pathway of DEGs in sweet potato seedlings root. KEGG pathway analysis shows the top 20 significantly enriched pathway. The x-axis shows different KEGG pathways, and the y-axis shows the numbers of differently expressed genes enriched in each pathway

Discussion

K is an essential macronutrient for physiological processes in plants, and the K supply has a significant effect on crop yield through the regulation of photosynthesis rate and assimilate transport [3]. However, the K+ content in each plant is significantly different. Some plants need to take up more K+ for their normal life activities, whereas others need little K+, suggesting that the mechanisms of K+ uptake and assimilation are different. K+-enriched plants function especially well when soil K+ availability is low and fluctuating, and the plant cells are able to maintain a relative high K+ concentration [49].

Low-K+ stress affects the growth and development of sweet potato, which requires K to promote storage root formation and bulking [50,51,52,53]. The production of sweet potato is remarkably promoted when provided with sufficient K [30]. Genome-wide gene expression profiles have been performed in rice [54, 55], wheat [56], pear [57, 58], cotton [59], and tomato [60]. However, little progress has been made in understanding the transcriptional changes and the K+ uptake and utilization mechanism in sweet potato under low-K+ conditions. In this study, sweet potato seedling phenotypes under low K+ were investigated (Fig. 1a) and genome-wide gene expression profiles were determined under low-K+ stress.

Ca2+, ethylene and reactive oxygen species (ROS) signature are involved in the response to nutrient deficiency and regulate K+ absorption and plant growth. A Ser/Thr protein kinase (CIPK23), which is activated by calcineurin B-like proteins CBL1/9, phosphorylates the K+ transporters AKT1 and HAK5 for K+ acquisition under low-K+ conditions, as mechanism that is found in many crops [9, 47, 61]. In this study, the Ca2+ signal related genes are involed in low-K+ stress, including genes coding for the calcium/calmodulin-dependent protein kinase and calmodulin-binding transcription activator were both upregulated (Additional file 2). ROS have important roles in many signal transduction pathways by regulating transporter activity [62,63,64]. Many studies demonstrate that ROS and ethylene accumulated in a discrete region of roots when the plants were deprived of K+ [65, 66]. A previous report showed that ethylene acts upstream of the ROS response to potassium deprivation by changing root hair and primary root growth in Arabidopsis [67]. A recent study revealed a mechanism whereby plants sense K+ deprivation and how this translates into spatially defined ROS signals to govern specific downstream responses. K+ deprivation triggered rapid K+ and Ca2+ signals in the root K+-sensing niche (KSN), meanwhile CIF peptides were induced by Ca2+ signals to activate SGN3-LKS4/SGN1 receptor complexes. LKS4 phosphorylated and activated RBOHC/D/F for ROS signal formation, and ROS signals conveyed HAK5 induction and accelerated casparian strip formation [68]. The ethylene-responsive AP2/ERF transcription factor RAP2.11 was reported to bind to the AtHAK5 promoter and positively regulated AtHAK5 expression to enhance absorption of K+ under K+-deprived conditions [69]. In the current study, ethylene-responsive transcription factor, such as g5338.t1 (ethylene-responsive transcription factor 4) and g38192.t1 (ethylene-responsive transcription factor ERF26) were upregulated at least four times under low-K+ stress. To verified whether transcription factor IbERF can bind to IbHAK5 promoter, we conducted Y1H assay. We cloned IbHAK5 gene and analysed its promoter, and we found there are two ERF binding sites in the promoter, expectively ERE box element and ATTTTAAA-motif (Fig. S3 and S4). Yeast-one-hybrid assays showing that IbERF can bind to the fragments of IbHAK5 promoter (Fig. S3). It implies that IbERF is an important transcription factor in the conditon of potassium deficiency, and further research is under way.

Auxin and jasmonic acid (JA) are also involved in the K signal transduction pathway [45, 70,71,72]. The auxin response factor ARF2 was reported to bind to the AtHAK5 promoter to maintain the low expression level of AtHAK5 under K+-sufficient conditions. When low K signals are perceived by plants, ARF2 is phosphorylated, relieving its repression of HAK5 transcription [73]. The Arabidopsis transcription factor MYB77, which modulates auxin signal transduction is down regulated by K+ deprivation to modulate lateral root development [74]. MYB77 positively regulates the expression of HAK5 by binding to the HAK5 promoter and enhances K+ uptake by roots [75]. TRH1 is a member of the AtKT/AtKUP/AtHAK family of K carriers and has been identified as an important part of the auxin transport system, affecting root gravitropic behaviour in Arabidopsis [76, 77]. Another member of the AtKT/AtKUP/AtHAK family, KUP9, which is highly expressed in quiescent center cells in the root tips maintains Arabidopsis root meristem activity and root growth by regulating K+ and auxin homeostasis in response to low-K+ stress [25]. In the current study, the auxin-related DEGs g3636 (homologous gene of nodulin MtN21) and g30544 (homologous gene of AtPIN8) were both downregulated under K+-deficient conditions. We hypothesize that these auxin-related DEGs are down regulated to reduce plant growth and adapt to the low-K+ conditions. MYC2, an important transcription factor related to the JA-responsive signaling pathway, is involved in multiple stress resistances pathways, such as against drought, salt, and fungal pathogens [78,79,80,81]. In the current study, a transcription factor similar to MYC2 was significantly upregulated, suggesting that the JA-response gene MYC2 is related to low-K+ stress (Table 1).

In addition to hormone-related genes involved in K+ deficiency, many other genes were found to be changed at the transcriptional level in sweet potato in response to low K+. The cell wall is important in plant defense against pathogens [38], and it is hypothesized that crop demands for K+ are closely related to cell wall growth. Cell wall-related genes involved in the response to low-K+ stress hanve been reported in Arabidopsis [45]. However, little published work supports this hypothesis. In our analysis, cell wall-related transcripts, such as expansin-A2-like, wall-associated receptor kinase-like, and lignin-forming anionic peroxidase, showed significant changes under low-K+ conditions (Table 1). These genes were highly expressed under low-K+ stress, indicated that K+ is associated with cell wall-related activity. K+/NO3 are generally absorbed and transported in a coordinated manner. CIPK23 was reported to mediate K+ uptake by phosphorylating and activating AtAKT1 [47]. Moreover, CIPK23 is involved in regulating nitrate uptake; thus it might act as a connecting node for K+ and NO3 during early stages of nutrient perception [1, 82]. The nitrate transporter NPF7.3/NRT1.5 functions as a proton-coupled H+/K+ antiporter with an essential role in K+ translocation from root to shoot, and is involved in the coordination of K+/NO3 distribution in Arabidopsis [48]. In the current study, a high-affinity nitrate transporter was upregulated under low K+ stress (Table 1), suggesting that this gene is involved in K+ homeostasis in sweet potato under K+-deficiency conditions. K is essential for cell turgor and sugar accumulation, especially in fruits, and low K+ stress can promote sugar accumulation in the root [58, 83]. In previous reports, low K+ stress was closely related to sugar metabolism and signaling [59, 71, 84]. In cassava, increasing the supply of K+ led to increased sucrose in the leaves, which is loaded in the phloem and transported to the roots for storage as starch, enhancing root crop yield [85]. Sugars will be eventually effluxed transporter (SWEETs) mediates the transport of sugars across the plasma membrane or the tonoplast [86]. Here, a sugar transport (Table 1) was downregulated in sweet potato under low K+ stress and we hypothesize that it might be involved in sugar metabolism and signaling. However, further research is required to explore the function of these key DEGs under K deficiency in sweet potato and other crops.

Conclusions

In conclusion, the current study analyzed the transcriptome profiles of sweet potato roots, showing that expression of hundreds of genes changed under low K+ stress. Many physiological processes changed in response to low K+ stress in sweet potato, and some previously unknow genes were revealed. These genes were associated with the kinase, phytohormone, ROS, Ca2+, transcription factors, plant defense, transporters and cell wall (Fig. 6 and Table 1). Those TFs that were significantly altered under low K+ condition might positively or negatively regulate K+ homeostasis in sweet potato. The network of K+ uptake and assimilation in sweet potato is complex, and further research is required to clarify the differences in transcriptional regulation. Thus, this research provides a good reference for the study of the relationship between low-K+ stress and other signal pathways in sweet potato and other important crops.

Fig. 6
figure 6

A proposed model of transcription regulation involved in low K+ stress in roots. The model is based on the differently expression genes of transcriptome described in this study. When response to potassium deficiency condition, the sweet potato seedlings roots activate multiple methods to maintain the potassium homeostasis

Methods

Plant materials and growth condition

‘Taizhong6’, a carotenoid-rich sweet potato cultivar, the China national accession number of which is 2,013,003, was used for this research and its genomic was sequenced [42]. It was propagated by taking 5-cm-long stem cuttings, which were then grown in a plot filled with vermiculite. When the plant reached the two-leaf stage, the seedlings were subjected to either low or high K+ under constant illumination with a 13 h light/11 h dark photoperiod, 28 °C day/24 °C night temperatures, and 65% relative humidity for 14 days. Three biological replicates were used for all analyses. All seedlings were watered with modified Hoagland’s solutions every 3 days [87]. The complete nutrient solution (HK, 1 mM K+) contained 0.75 mM K2SO4, 2 mM Ca (NO3)2.4H2O, 0.65 mM MgSO4.7H2O, 0.1 mM KCl, 0.25 mM KH2PO4, 0.1 mM Fe-EDTA, 1 mM MnSO4, 1 mM ZnSO4.7H2O, 0.01 mM CuSO4.5H2O, 0.005 mM (NH4)6Mo7O24.4H2O, and 1 mM H3BO3, pH 6.0. For the low-K+ groups (0 mM K+), KCl and KH2PO4 were replaced with 0.1 mM NaCl and 0.25 mM NaH2PO4.2H2O; the remaining components were not altered.

Content measurements

To determine K+ content determination, root and shoot tissues were harvested separately from low K+ stress plants after 14 days of growth and dried at 80 °C for 72 h. The total K was determined through an atomic absorption spectrophotometer [59]. Three biological replicates were used for the phenotypic test and K+ content determination. Student’s t-test was used to determine the statistical significance (*P < 0.05, **P < 0.01).

To determine the chlorophyll content, sweet potato stems were cut into 5-cm fragments and grown in plot full with vermiculite to the two-leaf stage. Seedlings were then treated with HK (1 mM K) and LK (0 mM K) Hoagland solution for 14 days, and the shoots were then collected. Chlorophyll was extracted in 80% acetone (v/v) in the dark for 2 days. The absorbance of the extraction buffer at 645 and 663 nm was then measured using microplate reader (BioTek Power Wave XS2) [48]. Three biological replicates were used for this experiment.

The roots of sweet potato seedlings with or without K+ treatment were harvested and used for carotenoids determine. The roots were ground with liquid nitrogen and into powder, and 0.1 g power was mixed with 10 mL extraction solvent. The extract set aside in the dark until the residues became colorless. The absorbance of the extract was measured at 440 nm for carotenoid by UV-Visible spectrophotometry [88].

Genome assembly and transcriptome sequencing

mRNA was isolated from total RNA using magnetic beads with oligo (dT) and fragmentation buffer was added to cut the mRNA into short fragments. cDNA was synthesized using the short fragments as templates through random hexamer-primers. Further dNTPs, RNase H, DNA polymerase I, and buffer solution were added to synthesize the second cDNA. The cDNA fragments were purified with the QiaQuick PCR Purification kit (Qiagen) and resolved with EB buffer for end reparation and poly (A) addition. The ligation products were selected based on size by agarose gel electrophoresis, PCR amplication, and sequencing were performed using an Illumina NovaSeq6000 (the second generation sequencing platform) to conduct the high-throughput transcriptome analysis of the control and low K+ stress-treated root samples. Three biological replicates were performed for each treatment.

RNA-seq read processing, assembly and transcriptome sequence analysis

After quality control, clean reads were mapped separately to the de novo genome assembly using trinity-v2.4.0. ORFs, SNPs, and SSRs were analyzed for each unigene. Each unigene was annotated with the NR (www.ncbi.nlm.nih.gov/), SWISSPROT (www.ebi.ac.uk/uniprot/) and KOG/COG databases (www.ncbi.nlm.nih.gov/research/cog/); GO enrichment analysis (http://geneontology.org/) and KEGG pathway analysis (https://www.kegg.jp/) were then performed for each annotated unigene [89, 90]. Finally, the expression levels of the unigenes were analyzed to reveal those that showed significant differential expression. The P-values obtained were adjusted using the Benjamini and Hochberg method for controlling the false discovery rate (FDR). FDR < 0.05 and log2|FoldChange| > 1 were set as the thresholds for significant differential expression. Thus genes with FDR-adjusted P-values < 0.05 were defined as DEGs [43]. Gene counts for each sample were imported into RStudio (www.rstudio.com/) with edgeR (www.bioconductor.org/packages/release/bioc/html/edgeR.html/) for DEGs analysis. Lowly expression genes were removed by filtering.

The gene ID used in this study was referenced by using the ‘Taizhong6’ database (https://ipomoea-genome.org/) [42]. In addition, sequencing data were submitted to the National Center for Biotechnology Information Sequence Read Archive (SRA) under accession number PRJNA760652 (www.ncbi.nlm.nih.gov/sra/PRJNA760652). All the annotation information of DEGs used in this research is provided in Additional files 2 and 5.

RNA isolation, gene cloning and RT-qPCR analyses

For RNA extraction, the shoot and root of sweet potato seedlings were collected after treatment with low K+ and complete nutrient solution for 14 days, immediately frozen in liquid nitrogen and stored at − 80 °C. All frozen samples were ground to a powder in liquid nitrogen and weighed (200 mg) to extract total RNA using an RNA extraction kit (EASYspin Plus Complex Plant RNA Kit). Three biological replicates were performed for each treatment. Quality characterization of RNA samples was determined and confirmed using a NanoDrop 2000 fluorospectrometer and formaldehyde denaturing gel electrophoresis.

The two maker genes, CIPK23 and HAK5, which were reported to respond to low K+ stress, were cloned from the ‘Taizhong6’ database. The sequences of 11 high DEGs were selected to conform the RNA-seq data. All the gene IDs used in this study were based on the ‘Taizhong6’ database [42]. The primers for these genes are shown in Table S4.

For RT-qPCR analysis, total RNA was treated with DNase I RNase Free (Takara) to eliminate genomic DNA contamination. Then, 10 ng cDNA and 50 nM of each primer were used for each quantitative PCR reaction, which was performed by using Powerup™ SYBR™ Green Master Mix (Applied Biosystems) on a QuantStudio 6 Flex PCR system machine (Thermo Fisher Scientific) following the manufacturers’ protocols. Thermal treatment was 10 min at 95 °C, followed by 40 cycles of 15 s at 95 °C and 1 min at 60 °C. Amplification was followed by a melt curve analysis. The 2–ΔΔCt method was used for relative quantification and IbACTIN was used as an internal reference for data normalization. Three biological replicates were used for each sample. Relative gene expression values were graphed using SigmaPlot v10.0 (Systat Software, https://systatsoftware.com/).

Yeast one-hybrid assay

The coding sequence of IbERF was constructed into the vector pB42AD. AD refers to the empty vector expressing the AD domain alone. LacZ was used as a reporter gene, driven by the fragments of IbHAK5 promoter in yeast. The pB42AD-IbERF and Pb42AD plasmid were co-transformed with the ProIbHAK5:lacZ plasmids into Saccharomyces cerevisiae strain EGY48 using standard transformation techniques respectively. After culturing on SD agar medium lacking Ura and leu (SD/−Ura/−Leu) at 30°Cfor 2 days. Yeast transformants were transferred onto plates SD agar medium lacking Ura and leu containing X-gal (5-bromo-4-chloro-3-indolyl-b-d-galactopyranoside) for blue colour development.

Availability of data and materials

All data generated or analyzed during this study are included in this article and additional files. The transcriptome sequence was finished in Berry Genomics Corporation (https://www.berrygenomics.com/). The datasets generated in this study are available from the NCBI SRA database under accession number PRJNA760652 (https://www.ncbi.nlm.nih.gov/sra/?term=PRJNA760652). The datasets supporting the conclusions of this article are included within the article (and its additional files).

Abbreviations

DEG:

Differently expressed gene

RT-PCR:

Real-time polymerase chain reaction

KUE:

K utilization efficiency

LK:

Low potassium

HK:

High potassium

SNP:

Singel-nucleotide polymorphisms

ORF:

Open-reading frame

TFs:

Transcription factors

ROS:

Reactive oxygen species

JA:

Jasmonic acid

FC:

Fold change

GO:

Gene ontology

KEGG:

Kyoto Encyclopedia of Genes and Genomes

FDR:

False discovery rate

FPKM:

Fragments per kilobase of exon per million fragments mapped

References

  1. Tsay YF, Ho CH, Chen HY, Lin SH. Integration of nitrogen and potassium signaling. Annu Rev Plant Biol. 2011;62:207–26. https://doi.org/10.1146/annurev-arplant-042110-103837.

    CAS  Article  PubMed  Google Scholar 

  2. Wang Y, Wu WH. Potassium transport and signaling in higher plants. Annu Rev Plant Biol. 2013;64:451–76. https://doi.org/10.1146/annurev-arplant-050312-120153.

    CAS  Article  PubMed  Google Scholar 

  3. Pettigrew WT. Potassium influences on yield and quality production for maize, wheat, soybean and cotton. Physiol Plant. 2008;133(4):670–81. https://doi.org/10.1111/j.1399-3054.2008.01073.x.

    CAS  Article  PubMed  Google Scholar 

  4. Leigh RA, Jones RGW. A hypothesis relating critical potassium concentrations for growth to the distribution and functions of this ion in the plant cell. New Phytol. 1984;97:1–13. https://doi.org/10.2307/2434189.

    CAS  Article  Google Scholar 

  5. Glass A. Plant mineral nutrition. An introduction to current concepts. Q Rev Biol. 1989;64(4):499. https://doi.org/10.1109/ICIMTR.2012.6236355.

    Article  Google Scholar 

  6. Fu HH, Luan S. AtKuP1: a dual-affinity K+ transporter from Arabidopsis. Plant Cell. 1998;10:63–73. https://doi.org/10.1105/tpc.10.1.63.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  7. Epstein E, Rains DW, Elzam OE. Resolution of dual mechanisms of potassium absorption by barley roots. Proc Natl Acad Sci U S A. 1963;49:684–92. https://doi.org/10.1073/pnas.49.5.684.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  8. Wang Y, Wu WH. Genetic approaches for improvement of the crop potassium acquisition and utilization efficiency. Curr Opin Plant Biol. 2015;25:46–52. https://doi.org/10.1016/j.pbi.2015.04.007.

    CAS  Article  PubMed  Google Scholar 

  9. Wang Y, Chen YF, Wu WH. Potassium and phosphorus transport and signaling in plants. J Integr Plant Biol. 2021;63(1):34–52. https://doi.org/10.1111/jipb.13053.

    CAS  Article  PubMed  Google Scholar 

  10. Sustr M, Soukup A, Tylova E. Potassium in root growth and development. Plants (Basel). 2019;8(10):435. https://doi.org/10.3390/plants8100435.

    CAS  Article  Google Scholar 

  11. Hirsch RE, Lewis BD, Spalding EP, Sussman MR. A role for the AKT1 potassium channel in plant nutrition. Science. 1998;280(5365):918–21. https://doi.org/10.1126/science.280.5365.918.

    CAS  Article  PubMed  Google Scholar 

  12. Lagarde D, Basset M, Lepetit M, Conejero G, Gaymard F, Astruc S, et al. Tissue-specific expression of Arabidopsis AKT1 gene is consistent with a role in K+ nutrition. Plant J. 1996;9(2):195–203. https://doi.org/10.1046/j.1365-313x.1996.09020195.x.

    CAS  Article  PubMed  Google Scholar 

  13. Spalding EP, Hirsch RE, Lewis DR, Qi Z, Sussman MR, Lewis BD. Potassium uptake supporting plant growth in the absence of AKT1 channel activity: inhibition by ammonium and stimulation by sodium. J Gen Physiol. 1999;113(6):909–18. https://doi.org/10.1085/jgp.113.6.909.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  14. Ragel P, Raddatz N, Leidi EO, Quintero FJ, Pardo JM. Regulation of K+ nutrition in plants. Front Plant Sci. 2019;10:281. https://doi.org/10.3389/fpls.2019.00281.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Gierth M, Maser P, Schroeder JI. The potassium transporter AtHAK5 functions in K(+) deprivation-induced high-affinity K(+) uptake and AKT1 K(+) channel contribution to K(+) uptake kinetics in Arabidopsis roots. Plant Physiol. 2005;137(3):1105–14. https://doi.org/10.1104/pp.104.057216.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  16. Pyo YJ, Gierth M, Schroeder JI, Cho MH. High-affinity K(+) transport in Arabidopsis: AtHAK5 and AKT1 are vital for seedling establishment and postgermination growth under low-potassium conditions. Plant Physiol. 2010;153(2):863–75. https://doi.org/10.1104/pp.110.154369.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  17. Rubio F, Nieves-Cordones M, Aleman F, Martinez V. Relative contribution of AtHAK5 and AtAKT1 to K+ uptake in the high-affinity range of concentrations. Physiol Plant. 2008;134(4):598–608. https://doi.org/10.1111/j.1399-3054.2008.01168.x.

    CAS  Article  PubMed  Google Scholar 

  18. Chérel I, Gaillard I. The complex fine-tuning of K+ fluxes in plants in relation to osmotic and ionic abiotic stresses. Int J Mol Sci. 2019;20(3):715. https://doi.org/10.3390/ijms20030715.

    CAS  Article  PubMed Central  Google Scholar 

  19. Hasanuzzaman M, Bhuyan M, Nahar K, Hossain M, Mahmud J, Hossen M, et al. Potassium: a vital regulator of plant responses and tolerance to abiotic stresses. Agron J. 2018;8(3):31. https://doi.org/10.3390/agronomy8030031.

    CAS  Article  Google Scholar 

  20. Zhang X, Wang G, Xue H, Zhang J, Wang Q, Zhang Z, et al. Metabolite profile of xylem sap in cotton seedlings is changed by K deficiency. Front Plant Sci. 2020;11:592591. https://doi.org/10.3389/fpls.2020.592591.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Liang M, Gao Y, Mao T, Zhang X, Zhang S, Zhang H, Song Z. Song. Characterization and expression of KT/HAK/KUP transporter family genes in willow under potassium deficiency, drought, and salt stresses. Biomed Res Int 2020; 2020: 2690760. https://doi.org/10.1155/2020/2690760.

  22. Rubio F, Nieves-Cordones M, Horie T, Shabala S. Doing ‘business as usual’ comes with a cost: evaluating energy cost of maintaining plant intracellular K+ homeostasis under saline conditions. New Phytol. 2020;225(3):1097–104. https://doi.org/10.1111/nph.15852.

    CAS  Article  PubMed  Google Scholar 

  23. Zhang M, Liang X, Wang L, Cao Y, Song W, Shi J, et al. A HAK family Na+ transporter confers natural variation of salt tolerance in maize. Nat Plants. 2019;5(12):1297–308. https://doi.org/10.1038/s41477-019-0565-y.

    CAS  Article  PubMed  Google Scholar 

  24. Chen G, Liu C, Gao Z, Zhang Y, Jiang H, Zhu L, et al. OsHAK1, a high-affinity potassium transporter, positively regulates responses to drought stress in rice. Front Plant Sci. 2017;8:1885. https://doi.org/10.3389/fpls.2017.01885.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Zhang ML, Huang PP, Ji Y, Wang S, Wang SS, Li Z, et al. KUP9 maintains root meristem activity by regulating K+ and auxin homeostasis in response to low K. EMBO Rep. 2020;21(6):e50164. https://doi.org/10.15252/embr.202050164.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  26. Okba SK, Mazrou Y, Elmenofy HM, Ezzat A, Salama AM. New insights of potassium sources impacts as foliar application on ‘Canino’ apricot fruit yield, fruit anatomy, quality and storability. Plants (Basel). 2021;10(6):1163. https://doi.org/10.3390/plants10061163.

    CAS  Article  Google Scholar 

  27. Kanai S, Ohkura K, Adu-Gyamfi JJ, Mohapatra PK, Nguyen NT, Saneoka H, et al. Depression of sink activity precedes the inhibition of biomass production in tomato plants subjected to potassium deficiency stress. J Exp Bot. 2007;58(11):2917–28. https://doi.org/10.1093/jxb/erm149.

    CAS  Article  PubMed  Google Scholar 

  28. Lester GE, Jifon JL, Makus DJ. Impact of potassium nutrition on postharvest fruit quality: melon (Cucumis melo L) case study. Plant Soil. 2009;335(1–2):117–31. https://doi.org/10.1007/s11104-009-0227-3.

    CAS  Article  Google Scholar 

  29. Zorb C, Senbayram M, Peiter E. Potassium in agriculture--status and perspectives. J Plant Physiol. 2014;171(9):656–69. https://doi.org/10.1016/j.jplph.2013.08.008.

    CAS  Article  PubMed  Google Scholar 

  30. Gao Y, Tang Z, Xia H, Sheng M, Liu M, Pan S, et al. Potassium fertilization stimulates sucrose-to-starch conversion and root formation in sweet potato (Ipomoea batatas (L.) lam.). Int J Mol Sci. 2021;22(9):4826. https://doi.org/10.3390/ijms22094826.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  31. Zayan MA, Mikhael GB, Okba SK. Treatments for improving tree growth, yield and fruit quality and for reducing double fruit and deep suture incidence in “desert red” peach trees. Int J Hortic Sci. 2016;22(3–4):7–19. https://doi.org/10.31421/IJHS/22/3-4./1187.

    Article  Google Scholar 

  32. Guo K, Tu L, He Y, Deng J, Wang M, Huang H, et al. Interaction between calcium and potassium modulates elongation rate in cotton fiber cells. J Exp Bot. 2017;68(18):5161–75. https://doi.org/10.1093/jxb/erx346.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  33. Ruan YL, Llewellyn DJ, Furbank RT. The control of single-celled cotton fiber elongation by developmentally reversible gating of plasmodesmata and coordinated expression of sucrose and K+ transporters and expansin. Plant Cell. 2001;13(1):47–60. https://doi.org/10.1105/tpc.13.1.47.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  34. Yang JS, Hu W, Zhao W, Meng Y, Chen B, Wang Y, et al. Soil potassium deficiency reduces cotton fiber strength by accelerating and shortening fiber development. Sci Rep. 2016;6:28856. https://doi.org/10.1038/srep28856.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  35. Yang J, Hu W, Zhao W, Chen B, Wang Y, Zhou Z, et al. Fruiting branch K(+) level affects cotton fiber elongation through osmoregulation. Front. Plant Sci. 2016;7:13. https://doi.org/10.3389/fpls.2016.00013.

    Article  Google Scholar 

  36. Hu W, Lv X, Yang J, Chen B, Zhao W, Meng Y, et al. Effects of potassium deficiency on antioxidant metabolism related to leaf senescence in cotton (Gossypium hirsutum L.). Field Crop Res. 2016;191:139–49. https://doi.org/10.1016/j.fcr.2016.02.025.

    Article  Google Scholar 

  37. Römheld V, Kirkby EA. Research on potassium in agriculture: needs and prospects. Plant Soil. 2010;335(1–2):155–80. https://doi.org/10.1007/s11104-010-0520-1.

    CAS  Article  Google Scholar 

  38. Amtmann A, Troufflard S, Armengaud P. The effect of potassium nutrition on pest and disease resistance in plants. Physiol Plant. 2008;133(4):682–91. https://doi.org/10.1111/j.1399-3054.2008.01075.x.

    CAS  Article  PubMed  Google Scholar 

  39. Zhou L, He H, Liu R, Han Q, Shou H, Liu B. Overexpression of GmAKT2 potassium channel enhances resistance to soybean mosaic virus. BMC Plant Biol. 2014;14:154. https://doi.org/10.1186/1471-2229-14-154.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  40. Liu Q. Improvement for agronomically important traits by gene engineering in sweetpotato. Breed Sci. 2017;67(1):15–26. https://doi.org/10.1270/jsbbs.16126.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  41. Nguyen HC, Chen CC, Lin KH, Chao PY, Lin HH, Huang MY. Bioactive compounds, antioxidants, and health benefits of sweet potato leaves. Molecules. 2021;26(7):1820. https://doi.org/10.3390/molecules26071820.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  42. Yang J, Moeinzadeh MH, Kuhl H, Helmuth J, Xiao P, Haas S, et al. Haplotype-resolved sweet potato genome traces back its hexaploidization history. Nat Plants. 2017;3(9):696–703. https://doi.org/10.1038/s41477-017-0002-z.

    CAS  Article  PubMed  Google Scholar 

  43. Arisha MH, Ahmad MQ, Tang W, Liu Y, Yan H, Kou M, et al. RNA-sequencing analysis revealed genes associated drought stress responses of different durations in hexaploid sweet potato. Sci Rep. 2020;10(1):12573. https://doi.org/10.1038/s41598-020-69232-3.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  44. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644–52. https://doi.org/10.1038/nbt.1883.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  45. Armengaud P, Breitling R, Amtmann A. The potassium-dependent transcriptome of Arabidopsis reveals a prominent role of jasmonic acid in nutrient signasling. Plant Physiol. 2004;136(1):2556–76. https://doi.org/10.1104/pp.104.046482.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  46. Qin YJ, Wu WH, Wang Y. ZmHAK5 and ZmHAK1 function in K+ uptake and distribution in maize under low K+ conditions. J Integr Plant Biol. 2019;61(6):691–705. https://doi.org/10.1111/jipb.12756.

    CAS  Article  PubMed  Google Scholar 

  47. Xu J, Li HD, Chen LQ, Wang Y, Liu LL, He L, et al. A protein kinase, interacting with two calcineurin B-like proteins, regulates K+ transporter AKT1 in Arabidopsis. Cell. 2006;125(7):1347–60. https://doi.org/10.1016/j.cell.2006.06.011.

    CAS  Article  PubMed  Google Scholar 

  48. Li H, Yu M, Du XQ, Wang ZF, Wu WH, Quintero FJ, et al. NRT1.5/NPF7.3 functions as a proton-coupled H+/K+ antiporter for K+ loading into the xylem in Arabidopsis. Plant Cell. 2017;29(8):2016–26. https://doi.org/10.1105/tpc.16.00972.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  49. Cherel I, Lefoulon C, Boeglin M, Sentenac H. Molecular mechanisms involved in plant adaptation to low K(+) availability. J Exp Bot. 2014;65(3):833–48. https://doi.org/10.1093/jxb/ert402.

    CAS  Article  PubMed  Google Scholar 

  50. O'Sullivan JN, Asher CJ, Blarney FPC. Nutrient disorders of sweet potato. Monographs. 1997:27–32. https://doi.org/10.22004/ag.econ.117165.

  51. Wang J, Wang H, Zhang Y, Zhou J, Chen X. Intraspecific variation in potassium uptake and utilization among sweet potato (Ipomoea batatas L.) genotypes. Field Crop Res. 2015;170:76–82. https://doi.org/10.1016/j.fcr.2014.10.007.

    Article  Google Scholar 

  52. Gajanayake B, Raja Reddy K, Shankle MW, Arancibia RA. Growth, developmental, and physiological responses of two sweetpotato (Ipomoea batatas L. [lam]) cultivars to early season soil moisture deficit. Sci Hortic. 2014;168:218–28. https://doi.org/10.1016/j.scienta.2014.01.018.

    Article  Google Scholar 

  53. Liu M, Zhang A-J, Chen X-G, Jin R, Li H-M, Tang Z-H. Effects of potassium deficiency on root morphology, ultrastructure and antioxidant enzyme system in sweet potato (Ipomoea batatas [L.] lam.) during early growth. Acta Physiol Plant. 2017;39(9):211. https://doi.org/10.1007/s11738-017-2512-8.

    CAS  Article  Google Scholar 

  54. Ma TL, Wu WH, Wang Y. Transcriptome analysis of rice root responses to potassium deficiency. BMC Plant Biol. 2012;12:161. https://doi.org/10.1186/1471-2229-12-161.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  55. Zhang X, Jiang H, Wang H, Cui J, Wang J, Hu J, et al. Transcriptome analysis of rice seedling roots in response to potassium deficiency. Sci Rep. 2017;7(1):5523. https://doi.org/10.1038/s41598-017-05887-9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  56. Zhao Y, Sun R, Liu H, Liu X, Xu K, Xiao K, et al. Multi-omics analyses reveal the molecular mechanisms underlying the adaptation of wheat (Triticum aestivum L.) to potassium deprivation. Front. Plant Sci. 2020;11:588994. https://doi.org/10.3389/fpls.2020.588994.

    Article  Google Scholar 

  57. Yang H, Li Y, Jin Y, Kan L, Shen C, Malladi A, et al. Transcriptome analysis of pyrus betulaefolia seedling root responses to short-term potassium deficiency. Int J Mol Sci. 2020;21(22):8857. https://doi.org/10.3390/ijms21228857.

    CAS  Article  PubMed Central  Google Scholar 

  58. Shen C, Wang J, Shi X, Kang Y, Xie C, Peng L, et al. Transcriptome analysis of differentially expressed genes induced by low and high potassium levels provides insight into fruit sugar metabolism of pear. Front Plant Sci. 2017;8:938. https://doi.org/10.3389/fpls.2017.00938.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Yang D, Li F, Yi F, Eneji AE, Tian X, Li Z. Transcriptome analysis unravels key factors involved in response to potassium deficiency and feedback regulation of K+ uptake in cotton roots. Int J Mol Sci. 2021;22(6):3133. https://doi.org/10.3390/ijms22063133.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  60. Zhao X, Liu Y, Liu X, Jiang J. Comparative transcriptome profiling of two tomato genotypes in response to potassium-deficiency stress. Int J Mol Sci. 2018;19(8):2402. https://doi.org/10.3390/ijms19082402.

    CAS  Article  PubMed Central  Google Scholar 

  61. Ragel P, Rodenas R, Garcia-Martin E, Andres Z, Villalta I, Nieves-Cordones M, et al. The CBL-interacting protein kinase CIPK23 regulates HAK5-mediated high-affinity K+ uptake in Arabidopsis roots. Plant Physiol. 2015;169(4):2863–73. https://doi.org/10.1104/pp.15.01401.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  62. Pei ZM, Murata Y, Benning G, Thomine S, Klusener B, Allen GJ, et al. Calcium channels activated by hydrogen peroxide mediate abscisic acid signalling in guard cells. Nature. 2000;406(6797):731–4. https://doi.org/10.1038/35021067.

    CAS  Article  PubMed  Google Scholar 

  63. Mustilli AC, Merlot S, Vavasseur A, Fenzi F, Giraudat J. Arabidopsis OST1 protein kinase mediates the regulation of stomatal aperture by abscisic acid and acts upstream of reactive oxygen species production. Plant Cell. 2002;14(12):3089–99. https://doi.org/10.1105/tpc.007906.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  64. Neill SJ. Hydrogen peroxide and nitric oxide as signalling molecules in plants. J Exp Bot. 2002;53(372):1237–47. https://doi.org/10.1093/jexbot/53.372.1237.

    CAS  Article  PubMed  Google Scholar 

  65. Shin R, Schachtman DP. Hydrogen peroxide mediates plant root cell response to nutrient deprivation. Proc Natl Acad Sci U S A. 2004;101(23):8827–32. https://doi.org/10.1073/pnas.0401707101.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  66. Shin R, Berg RH, Schachtman DP. Schachtman, reactive oxygen species and root hairs in Arabidopsis root response to nitrogen, phosphorus and potassium deficiency. Plant Cell Physiol. 2005;46(8):1350–7. https://doi.org/10.1093/pcp/pci145.

    CAS  Article  PubMed  Google Scholar 

  67. Jung JY, Shin R, Schachtman DP. Ethylene mediates response and tolerance to potassium deprivation in Arabidopsis. Plant Cell. 2009;21(2):607–21. https://doi.org/10.1105/tpc.108.063099.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  68. Wang FL, Tan YL, Wallrad L, Du XQ, Eickelkamp A, Wang ZF, et al. A potassium-sensing niche in Arabidopsis roots orchestrates signaling and adaptation responses to maintain nutrient homeostasis. Dev Cell. 2021;56(6):781–94. https://doi.org/10.1016/j.devcel.2021.02.027.

    CAS  Article  PubMed  Google Scholar 

  69. Kim MJ, Ruzicka D, Shin R, Schachtman DP. The Arabidopsis AP2/ERF transcription factor RAP2.11 modulates plant response to low-potassium conditions. Mol Plant. 2012;5(5):1042–457. https://doi.org/10.1093/mp/sss003.

    CAS  Article  PubMed  Google Scholar 

  70. Li G, Wu Y, Liu G, Xiao X, Wang P, Gao T, et al. Large-scale proteomics combined with transgenic experiments demonstrates an important role of Jasmonic acid in potassium deficiency response in wheat and rice. Mol Cell Proteomics. 2017;16(11):1889–905. https://doi.org/10.1074/mcp.RA117.000032.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  71. Tahir MM, Chen S, Ma X, Li S, Zhang X, Shao Y, et al. Transcriptome analysis reveals the promotive effect of potassium by hormones and sugar signaling pathways during adventitious roots formation in the apple rootstock. Plant Physiol Biochem. 2021;165:123–36. https://doi.org/10.1016/j.plaphy.2021.05.015.

    CAS  Article  PubMed  Google Scholar 

  72. Deepika SA. Expression dynamics indicate the role of Jasmonic acid biosynthesis pathway in regulating macronutrient (N, P and K+) deficiency tolerance in rice (Oryza sativa L.). Plant Cell Rep. 2021;40(8):1495–512. https://doi.org/10.1007/s00299-021-02721-5.

    CAS  Article  PubMed  Google Scholar 

  73. Zhao S, Zhang ML, Ma TL, Wang Y. Phosphorylation of ARF2 relieves its repression of transcription of the K+ transporter gene HAK5 in response to low potassium stress. Plant Cell. 2016;28(12):3005–19. https://doi.org/10.1105/tpc.16.00684.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  74. Shin R, Burch AY, Huppert KA, Tiwari SB, Murphy AS, Guilfoyle TJ, et al. The Arabidopsis transcription factor MYB77 modulates auxin signal transduction. Plant Cell. 2007;19(8):2440–53. https://doi.org/10.1105/tpc.107.050963.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  75. Feng CZ, Luo YX, Wang PD, Gilliham M, Long Y. MYB77 regulates high-affinity potassium uptake by promoting expression of HAK5. New Phytol. 2021;232(1):176–89. https://doi.org/10.1111/nph.17589.

    CAS  Article  PubMed  Google Scholar 

  76. Vicente-Agullo F, Rigas S, Desbrosses G, Dolan L, Hatzopoulos P, Grabov A. Potassium carrier TRH1 is required for auxin transport in Arabidopsis roots. Plant J. 2004;40(4):523–35. https://doi.org/10.1111/j.1365-313X.2004.02230.x.

    CAS  Article  PubMed  Google Scholar 

  77. Rigas S, Ditengou FA, Ljung K, Daras G, Tietz O, Palme K, et al. Root gravitropism and root hair development constitute coupled developmental responses regulated by auxin homeostasis in the Arabidopsis root apex. New Phytol. 2013;197(4):1130–41. https://doi.org/10.1111/nph.12092.

    CAS  Article  PubMed  Google Scholar 

  78. Dombrecht B, Xue GP, Sprague SJ, Kirkegaard JA, Ross JJ, Reid JB, et al. MYC2 differentially modulates diverse jasmonate-dependent functions in Arabidopsis. Plant Cell. 2007;19(7):2225–45. https://doi.org/10.1105/tpc.106.048017.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  79. Abe H, Urao T, Ito T, Seki M, Shinozaki K, Yamaguchi-Shinozaki K. Arabidopsis AtMYC2 (bHLH) and AtMYB2 (MYB) function as transcriptional activators in abscisic acid signaling. Plant Cell. 2003;15(1):63–78. https://doi.org/10.1105/tpc.006130.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  80. Verma D, Jalmi SK, Bhagat PK, Verma N, Sinha AK. A bHLH transcription factor, MYC2, imparts salt intolerance by regulating proline biosynthesis in Arabidopsis. FEBS J. 2020;287(12):2560–76. https://doi.org/10.1111/febs.15157.

    CAS  Article  PubMed  Google Scholar 

  81. Du M, Zhao J, Tzeng DTW, Liu Y, Deng L, Yang T, et al. MYC2 orchestrates a hierarchical transcriptional cascade that regulates jasmonate-mediated plant immunity in tomato. Plant Cell. 2017;29(8):1883–906. https://doi.org/10.1105/tpc.16.00953.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  82. Ho CH, Lin SH, Hu HC, Tsay YF. CHL1 functions as a nitrate sensor in plants. Cell. 2009;138(6):1184–94. https://doi.org/10.1016/j.cell.2009.07.004.

    CAS  Article  PubMed  Google Scholar 

  83. Armengaud P, Sulpice R, Miller AJ, Stitt M, Amtmann A, Gibon Y. Multilevel analysis of primary metabolism provides new insights into the role of potassium nutrition for glycolysis and nitrogen assimilation in Arabidopsis roots. Plant Physiol. 2009;150(2):772–85. https://doi.org/10.1104/pp.108.133629.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  84. Siddiqui MH, Khan MN, Mukherjee S, Alamri S, Basahi RA, Al-Amri AA, et al. Hydrogen sulfide (H2S) and potassium (K+) synergistically induce drought stress tolerance through regulation of H+-ATPase activity, sugar metabolism, and antioxidative defense in tomato seedlings. Plant Cell Rep. 2021;40(8):1543–64. https://doi.org/10.1007/s00299-021-02731-3.

    CAS  Article  PubMed  Google Scholar 

  85. Omondi JO, Lazarovitch N, Rachmilevitch S, Kukew T, Yermiyahu U, Yasuor H. Potassium and storage root development: focusing on photosynthesis, metabolites and soluble carbohydrates in cassava. Physiol Plant. 2020;169(2):169–78. https://doi.org/10.1111/ppl.13060.

    CAS  Article  PubMed  Google Scholar 

  86. Chen HY, Huh JH, Yu YC, Ho LH, Chen LQ, Tholl D, et al. The Arabidopsis vacuolar sugar transporter SWEET2 limits carbon sequestration from roots and restricts Pythium infection. Plant J. 2015;83(6):1046–58. https://doi.org/10.1111/tpj.12948.

    CAS  Article  PubMed  Google Scholar 

  87. Hoagland DR, Arnon DI. The water-culture method for growing plants without soil. California department of agriculture experimental station sircular, vol. 347; 1950. p. 1–32. https://doi.org/10.1016/S0140-6736(00)73482-9.

    Book  Google Scholar 

  88. Li X, Sun J, Chen Z, Jiang J, Jackson A. Characterization of carotenoids and phenolics during fruit ripening of Chinese raspberry (Rubus chingii Hu). RSC Adv. 2021;11(18):10804–13. https://doi.org/10.1039/d0ra10373j.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  89. Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49(D1):D545–51. https://doi.org/10.1093/nar/gkaa970.

    CAS  Article  PubMed  Google Scholar 

  90. Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30. https://doi.org/10.1093/nar/28.1.27.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgments

The transcriptome sequence was carried in Berry Genomics Corporation (https://www.berrygenomics.com/). We thank International Science Editing (http://www.internationalscienceediting.com/) for editing this manuscript.

Funding

This work was supported by Talented Introduction Project of Sichuan Academy of Agricultural Sciences (20213P10) and China Agriculture Research System of MOF and MARA (CARS-10-BO4).

Author information

Authors and Affiliations

Authors

Contributions

Fang Wang conducted the experiments, analyzed the results. Fang Wang and Wen-Fang Tan performed the figures and wrote the paper. Wei Song, Song-Tao Yang and Shuai Qiao assisted in the experimental procedures. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Fang Wang.

Ethics declarations

Ethics approval and consent to participate

The sweet potato material ‘Taizhong6’ used in this study was come from China Agriculture Research System, and in compliance with the appropriate permissions and/or licenses for the collection of seed specimens. Experimental research on plants (either cultivated or wild), including the collection of plant material, comply with relevant institutional, national, and international guidelines and legislation.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Fig. S1.

Phenotypic and quantitative data of sweet potato under normal and K+-deficient conditions. a. Phenotype of the root. b. Fresh weight of the shoot and root. c. The dry weight of the shoot and root. d. The content of carotenoids in root under low-K+ stress. e. Chlorophyll content of shoot. Sweet potato stems were cut into 5 cm fragments and growth in plot full with vermiculite to two leaf stage. Seedlings were then treated with HK (1 mM K) and LK (0 mM K) Hoagland solution for 14 days, and then the weight, carotenoids and chlorophyll content were detected. Data are shown as means ± SE (n = 4). Student’s t test (*P < 0.05) was used to analyze statistical significance, HK-S: sufficient potassium shoot, HK-R: sufficient potassium root, LK-S: low potassium shoot, LK-R: low potassium root, bar=5 cm. Fig. S2. The numbers of different unigenes length and venn diagrams of different annotation database. a. Distribution of the numbers of different length unigenes. b. Venn diagrams of transcriptions in each annotation database. Fig. S3. Transcription factor ERF binds to IbHAK5 promoter. a. Diagram of the IbHAK5 promoter. The adenine residue of the translational start codon ATG was assigned position +1, and the two ERF binding motif were showed in different color. Relative positions of the two motif were indicated by red and green lines. The scale length is 200 bp. b. Transient expression of the ProIbHAK5:lacZ fusion together with IbERF in yeast. AD together with ProIbHAK5:lacZ was taken as negative control. Observe the color and then photoes were taken. Fig. S4. Amino acid sequence alignment of AtHAK5 and IbHAK5. Table S1. Summary of RNA-seq quality information. Table S2. Unigenes statistical table. Table S3. Statistical table of unigenes annotation. Table S4. primers used in this study. Table S5. Differently expressed unigenes in KEGG.

Additional file 2: Supplement Table 6.

The raw expression data and annotation of all DEGs under low-K+ strss and high potassium condition.

Additional file 3: Supplement Table 7.

Biological process of genes based on the gene ontology. Supplement Table 8. Cellular component of genes based on the gene ontology. Supplement Table 9. Molecular function of genes based on gene ontology.

Additional file 4: Supplement Table 10.

The top 20 KEGG pathways of DEGs under low-K+ strss and high potassium condition.

Additional file 5.

The sequence of all unigenes.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Wang, F., Tan, WF., Song, W. et al. Transcriptome analysis of sweet potato responses to potassium deficiency. BMC Genomics 23, 655 (2022). https://doi.org/10.1186/s12864-022-08870-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-022-08870-5

Keywords

  • Sweet potato
  • Low-K+ stress
  • Transcriptomic
  • DEGs