Transcriptome Profiling of the Potato (Solanum tuberosum L.) Plant under Drought Stress and Water-Stimulus Conditions

Drought stress can seriously affect tuberization, yield and quality of potato plant. However, the precise molecular mechanisms governing potato stolon’s response to drought stress and water supply are not very well understood. In this work, a potato (Solanum tuberosum L.) variant, Ningshu 4, was subjected to severe drought stress treatment (DT) and re-watering treatment (RWT) at tuber bulking stage. Strand-specific cDNA libraries of stolon materials were constructed for paired-end transcriptome sequencing analyses and differentially expressed gene (DEG) examination. In comparison to untreated-control (CT) plants, 3189 and 1797 DEGs were identified in DT and RWT plants and 4154 solely expressed DEGs were screened out from these two comparison groups. Interestingly, 263 genes showed opposite expression patterns in DT and RWT plants. Among them, genes homologous to Protein Phosphatase 2C (PP2C), Aspartic protease in guard cell 1 (ASPG1), auxin-responsive protein, Arabidopsis pseudo response regualtor 2 (APRR2), GA stimulated transcripts in Arabidopsis 6 (GASA6), Calmodulin-like protein 19 (CML19), abscisic acid 8'-hydroxylases and calcium-transporting ATPase, et al. were related with drought-stress and water stimulus response. Sixteen DEGs involved in starch synthesis, accumulation and tuber formation exhibited significantly different expression upon re-watering. In addition, 1630, 1527 and 1596 transcription factor encoding genes were detected in CT, DT and RWT. DEGs of ERF, bHLH, MYB, NAC, WRKY, C2H2, bZIP and HD-ZIP families accounted for 50% in three comparison groups, respectively. Furthermore, characteristics of 565 gene ontology (GO) and 108 Kyoto Encyclopedia of Genes and Genomes pathways (KEGG) were analyzed with the 4154 DEGs. All these results suggest that the drought- and water-stimulus response could be implemented by the regulated expression of metabolic pathway DEGs, and these genes were involved in the endogenous hormone biosynthesis and signal transduction pathways. Our data provide more direct information for future study on the interaction between key genes involved in various metabolic pathways under drought stress in potato.


Introduction
Global climate change is reducing the reliability of rainfall, the availability of soil-water, and consequently, limiting plant production. Potato (Solanum tuberosum L.) ranks as the fourth most predominant non-grain food crop in the world. Due to their shallow root system, which limits water extraction from soil, potato plants are sensitive to drought stress. Leaf size, photosynthesis rate, tuber number, yield and quality were all severely limited when grown under a drought stress condition [1,2]. Therefore, in order to breed new potato varieties with improved drought resistance, the molecular mechanisms governing water use, tuber starch synthesis and accumulation under drought conditions must be defined.
In higher plants, drought stress can induce a wide range of physiological and biochemical responses regulated by genes encoding functional and/or regulatory proteins, to maintain a normal homeostasis and avoid wilting death [3]. Fast and efficient recovery from water stress will allow agricultural crops to adapt to the changes of meteorological conditions and maximize their growth and production during drought stress. In a diploid potato population, 47 quantitative trait loci were identified, out of them, 28 were drought-specific, and 17 were specifically expressed during recovery from drought stress [4]. Recently, the published potato genome [5] and transcriptome [6] with corresponding gene structure, location, and functional annotation provide powerful resources for the understanding of the complex responses to drought stress in potato. However, data currently available on transcriptomic changes upon rehydration in potato is insufficient, and the precise pathways involved in drought responses and recovery from water deficit still remain to be deduced.
Transcriptome sequencing technologies provide a framework dataset for researches related to construction of transcriptome map, determination of metabolic pathways, clarification of gene expression patterns, and mining of new genes. To explore the molecular regulation mechanisms of plants in response to different stresses, transcriptome analysis has been performed on a variety of crops, including sunflower (Helianthus annuus) [7], grape (Vitis vinifera) [8], maize (Zea mays) [9] and sorghum (Sorghum bicolor L. Moench) [10]. Data reported in these studies have provided information on the networks of regulatory and functional genes at different levels.
Upon the completion of whole genome and transcriptome sequencing of potato DM1-3 516R44 [5,6], a number of genes related to drought-stress response and water use efficiency have been identified by quantitative trait locus mapping [2,4], and transcriptome analyses of stress response in potato leaves has been performed with plants grown under a single condition. However, molecular mechanism of potato plant, especially the tuber or stolon responsing to water-stimulus are not very well understood with the results available now. Tubers are the main "sink" organ and it also the most important agronomic trait of potato plants. In addition, the tuber bulking stage is a critical period of time for potato tuberization, during which plants are extremely sensitive to water deficit. Therefore, we explored the molecular mechanisms underpinning drought stress responses in the stolon tissue of a potato variety. Plant cultivated under severe drought stress condition during tuber bulking stage were used for next-generation transcriptome sequencing to fully resolve gene expression profiles in response to water deficit and re-watering.

Plant materials and drought treatments
The potato (Solanum tuberosum L.) strain Ningshu 4 was from our own lab. All plant materials were collected from our Guyuan Experimental Field (Guyuan, Ningxia Hui Autonomous Region, China), which does not involve any endangered or protected species. Before tuber bulking stage, 40 healthy and uniform ball-plants were transplanted into buckets and cultivated indoors at room temperature under 12-hr light/dark cycles at 4,000 lux. Experimental soil was collected from a test site in Guyuan City, Ningxia Hui Autonomous Region, with the maximum field water-holding capacity of 21%.
At flowering stage, the above plants were split in to three groups with ten plants in each group. The first group of plants was exposed to severe drought stress treatment (DT) for three days, during which relative soil moisture was controlled at 35-40% of the maximum field waterholding capacity. The second group of plants was exposed to severe drought stress treatment for 3 days, followed by re-watering to the control level and grown for 3 days (RWT). And the last group of plants was maintained under normal watering condition throughout the time period as control (CT). Then, stolon tips from three randomly selected plants of each group were pooled together, immediately grounded into fine powder in liquid nitrogen for sequencing analysis.

Total RNA extraction
Total RNA was extracted from plant materials using a TransZol Plant kit (ET121-01, Transgen Biotech, Beijing, China) according to the manufacturer's instruction. Digestion of gDNA in each sample was conducted using DNase I (2212, TaKaRa, Dalian, China). The integrity and purity of RNA was confirmed using 1.2% non-denaturing agarose gel electrophoresis and the Agilent 2100 RNA 6000 Kit (Agilent Technologies, Santa Clara, CA, USA).

Strand-specific cDNA library construction and sequencing
From each sample, 3 μg of total RNA was taken to construct strand-specific cDNA libraries using the NEBNext Ultra RNA Library Prep Kit (New England labs, Ipswich, MS, USA). The quality of cDNA libraries was tested, and quantified cDNA was subjected to purification, elution, end repair, poly(A) addition and adaptor ligation. Then, agarose gel electrophoresis was used for size selection of gene fragments, followed by PCR amplification. The established cDNA libraries were sequenced using the Illumina Hiseq 2000 platform (Illumina Inc., San Diego, CA, USA) to generate 100-bp paired-end reads.

Sequence analysis
Image data obtained from cDNA sequencing was converted to the corresponding nucleotide sequence data. Raw reads were firstly processed through in-house perl scripts. In this step, clean data (clean reads) were obtained by removing reads containing adapter, reads containing ploy-N and low quality reads from raw data. At the same time, Q20, Q30, GC content and sequence duplication level of the clean data were calculated. All the downstream analyses were based on the clean data with high quality. Then, index of the reference genome (PGSC_DM_v3_2.1.10) was built using Bowtie v2.1.0 and paired-end clean reads were aligned to the reference genome using TopHat v2.0.9 (Broad Institute, Boston, MA) [11]. Next, HTSeq v0.5.3 was used to count the reads numbers mapped to each gene. And RPKM [12] of each gene was calculated based on the length of the gene and reads count mapped to this gene. All the clean data have been deposited in the Short Read Archive (SRA) at the NCBI database with the project accession number SRP056128.
Screening, clustering and functional annotation of differentially expressed genes (DEGs) described previously. A q value <0.005 and a fold change |log2| !1 were used as the screening cutoffs for extremely significant differential gene expression between two samples.
Following pairwise comparison of samples, gene ontology (GO) classification of DEGs was performed using the GO seq software [14] with corrected p-Value <0.05 as the screening cutoff for enriched GO terms. GO enrichment was analyzed in the categories of molecular function, cellular component and biological process. Enrichment analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathways was performed using KOBAS [15].

Quantitative real-time PCR (qRT-PCR) verification of DEGs
To test the validity and accuracy of transcriptome sequencing data, 10 genes were randomly selected from sequencing data for fluorescent qRT-PCR. Primers were designed using Primer Express software (v3.0, Applied Biosystems), and primer sequences were shown in S1 Table. The elongation factor 1-α (ef1α) was used as a reference gene [16]. The 20 μl PCR reaction buffer contains: 2 μl diluted cDNA, 1 μl each of forward and reverse primers (10 μmolÁL -1 ), 10 μl qPCR Master Mix (ABI, Foster City, CA, USA) and 6 μl ddH 2 O. Amplification of genes was performed on a LightCycler 480 Real-Time PCR machine (Roche, Basel, Switzerland) under the following conditions: 95°C for 3 min, followed by 40 cycles at 95°C for 15 sec and 60°C for 40 sec. Detection of each gene was repeated 3 times. The 2 -(ΔΔCt) method [17] was used to calculate and calibrate the expression level of target genes in different treatments.

Total RNA quality testing
The quality of total RNA extracted from the collected stolons of CT, DT and RWT plants was tested using agarose gel electrophoresis and an Agilent 2100 Bioanalyzer. All the three samples had 28S:18S ratios in a range of 1.8-2.0 with intact 28S, 18S and 5S RNA bands, high RNA purity, and a mean RNA integrity number (RIN) >8.0, which met the requirements for library construction and sequencing.

Sequencing data and DEG analysis
After transcriptome sequencing, a total of 78448498, 74764759, and 71081020 raw reads were obtained from CT, DT, and RWT cDNA libraries, respectively. From these reads, 21.74 G clean bases were obtained by filtering impurities, and 76.1% of them fell into the published potato genome. Sequencing analysis of CT, DT and RWT cDNA libraries produced 46302, 47166 and 46341 genes, and 75924, 76516 and 74678 transcripts, respectively (Table 1). Then, data obtained from the three samples were subjected to pairwise comparisons, and DEGs were identified using the preset cutoffs. A total of 9216 DEGs (including 5604 solely expressed DEGs) were identified. When DT was compared to CT, 1227 genes were up-regulated, and 1962 genes were down-regulated. When RWT was compared to CT, 658 genes were found to be up-regulated and 1139 gene down-regulated. Different from the above two comparison groups, 2352 genes were up-regulated, and 1878 genes showed down-regulated mode when RWT was compared to DT. Down-regulated transcripts accounts for 62%, 63% and 44% of the total DEGs in above three comparison groups, respectively (S2 Table).

qRT-PCR verification
To verify data validity, ten genes, including the target genes of interest, were selected from sequencing data, and qRT-PCR analyses were performed (S5 Table). In the comparison of DT and RWT, the correlation coefficients of gene expression trends in sequencing data and qRT-PCR results were 0.9267 and 0.5543, respectively (Fig 1). In particular, significant expression difference occurred in drought stress-and water stimulus-responsive genes before and after the treatment, including genes encoding MYB60 (PGSC0003DMG400033043), WRKY33 (PGSC0003DMG400011633) and NAC072 (PGSC0003DMG400015342), thus supporting the validity of DEG data obtained.

GO classification and KEGG pathway annotation analysis of DEGs
In order to investigate their functions and highlight metabolic pathways potentially related to drought tolerance, GO classification and KEGG pathway annotation analyses were performed on the identified DEGs. A total of 639 unique GO functional annotation terms were assigned to the identified DEGs (Fig 2). Of them, 979 were within the biological process category, accounting for 73.4% of the total GO functional annotation categories, and the rest were within the cellular component (217) and molecular function categories (137). Intra-group analysis ( Table 6)   significantly enriched GO functional annotation categories were translation (GO:0006412), ribosome (GO:000 5840) and structural constituent of ribosome (GO:0003735). These included 274 DEGs, of which, 256 were up-regulated genes. Clearly, DEGs were mainly down-regulated in DT and up-regulated in RWT. Similar to the RWT vs. CT group, the top 3 significantly enriched GO functional annotation categories in the comparison of DT and RWT were translation (GO:0006412), ribosome (GO:0005840) and structural constituent of ribosome (GO:0003735) either. However, each of the GO category include much more DEGs. Totally, this comparison group includes 713 DEGs, of which, 675 were up-regulated genes. A directed acyclic graph (DAG) was used to illustrate the affiliation and functional scope of the 10 most significantly enriched GO categories. In the comparison of CT and DT (Fig 3A), DEGs were intensively enriched in the following five GO categories: metabolic process (GO:0008152), stress response (GO:0006950), organic substance metabolic process (GO:0071704), primary metabolic process (GO:0044238) and cellular macromolecular complex assembly (GO:0034622), which included 4076 DEGs. In the comparison of CT and RWT (Fig 3B), the top 3 significantly enriched GO categories were metabolic process (GO:0008152), cellular metabolic process (GO:0044237) and translation (GO:0006412). To further identify the key regulatory genes related to drought-and water-stimulus responses, as well as genes involved in starch and sucrose metabolism, we focused on DEGs of the same GO terms but showed opposite expression regulation patterns between the two comparison treatments, that is, transcripts that were regulated in an opposite direction in DT and RWT compared to CT. A total number of 5 GO terms were related to starch and sucrose metabolism, and 3 GO terms were related to drought-stress and water-stimulus responses (S6 Table). Among these 8 GO terms, 60 DEGs showed opposite regulation patterns between DT and RWT. Further analysis revealed that 41 genes were down-regulated in DT and up-regulated in RWT, of which 37 were related to starch and sucrose metabolism, including homologous  genes encoding beta-galactosidase 1 (BGAL1), UDP-glucose 4-epimerase and 6-phosphogluconate dehydrogenase. Additionally, 19 genes were up-regulated in DT and down-regulated in RWT, of which 13 were related to drought-stress and water-stimulus responses.
In the pathway annotation analysis, data was compared with the published KOBAS database. A total number of 4549 DEGs were enriched in 314 (114 unique) pathways. Intra-group analysis of each comparison was then performed to identify the 10 most significantly enriched pathways (Fig 4). Revealing DEGs were intensively enriched in 14 pathways. The metabolic pathways with the most significant DEGs enrichment were metabolic pathways (sly01100), biosynthesis of secondary metabolites (sly01110), ribosome (sly03010) and carbon cycle metabolism (sly01200), which included 868, 446, 297, and 143 DEGs, respectively.
Furthermore, we performed a detailed classification of KEGG pathways for the 263 genes which showed opposite expression regulation patterns between the two comparison treatments, revealing 66 DEGs annotated to KEGG pathways (S7 Table) including 16 metabolic pathways (Table 7). Of them, the most frequently associated pathways were metabolic pathways (21), followed by ribosome (18). Only 3 enriched DEGs were associated with carotenoid biosynthesis and plant hormone signal transduction pathways (S1A and S1B Fig). These included PGSC0003DMG402018475 (encoding abscisic acid 8'-hydroxylase 1-like, CYP707A2), PGSC0003DMG400001960 (encoding abscisic acid 8'-hydroxylase 1-like) and PGSC0003DMG400028897 (encoding beta-carotene hydroxylase) involved in carotenoid biosynthesis pathway, and PGSC0003DMG400008011 (encoding AREB-like protein, AREB), PGSC0003DMG40 1025624 (encoding two-component response regulator ARR2-like) and PGSC0003DM G400006093 (encoding auxin-responsive protein, AUX) associated with plant hormone signal transduction pathway. The above results suggest that DEGs with opposite expression regulation patterns were mainly involved in drought-and water-stimulus response and regulation via the carotenoid biosynthesis pathway, which affect the synthesis and signal transduction pathways of other endogenous hormones (e.g., abscisic acid and indole acetic acid).

Drought stress and water stimulus responsive TFs
TF-encoding genes are key factors in the process of stress signal perception and transduction. Based on the sequencing analyses, a total number of 1630, 1527 and 1596 TF-encoding genes were detected in CT, DT and RWT, respectively (Table 8). In the comparison to the control treatment, both drought and re-watering treatments reduced the number of TF-encoding genes detected. In the comparison between DT and RWT with CT, a total number of 1019 and 866 TF-encoding DEGs were identified, respectively.
Further analyses revealed that in the comparison between DT and CT, the identified 1019 DEGs fell into 50 TF families, and most of DEGs (268) belong to bHLH, ERF and MYB families (S8 Table). Similarly, in the comparison between RWT and CT, the identified DEGs were involved in a total of 52 TF families, and most of them (230) also belong to ERF, bHLH, and MYB families. When it comes to RWT vs. DT, 1018 DEGs were involved in a total of 51 TF families. bHLH (97), ERF (90), C2H2 (57) and MYB (57) were the top 4 families with the most differently expressed genes. Overall (S8 Table,  DEGs belonging to the above eight transcription factor families accounted for almost 50% of the total differentially expressed TFs detected in the DT and RWT treatments, respectively. We also analyzed the quantity relationship of up-and down-regulated TFs in 3 comparison groups (S8 Table). In the comparison between DT and CT, the number of down-regulated DEGs was approximately 2-fold that of the up-regulated ones among the 1019 differentially expressed TFs. The largest number of down-regulated DEGs was found in the bHLH family, while that of up-regulated DEGs was mostly found in the ERF family. Similar with comparison between DT and CT, the number of down-regulation DEGs in RWT vs. DT exhibited almost twice as much as the number of up-regulation TFs. The bHLH family contained 69 decreased expression of TFs and the ERF family included 31 increased expression of DEGs. In the comparison between RWT and CT, the number of up-or down-regulated DEGs was about the same among the 866 differentially expressed TFs. The largest number of down-regulated DEGs was found in the ERF family, whereas the largest number of up-regulated ones was found in the bHLH family.

Discussion
Rapid growth of global population and decrease of arable land have presented a serious challenge for food production. Potato is a non grain crop and one of the most important sources for food and industry uses worldwide. One of the most essential approaches is to breed new plant species which can grow on uncultivated lands. Understanding how plants including potato respond to water-deficit stress will facilitate the breeding of drought tolerant plants. In the present study, we employed high-throughput sequencing technology to profile the transcriptome changes in the stolon tissues of potato undergoing water-deficit stress and re-watering treatment. A total number of 3189, 1797 and 4230 DEGs (S2 Table), including 1630, 1527 and 1596 transcriptional factor-encoding DEGs (Table 8), were identified in the comparisons within CT, DT and RWT, respectively. These DEGs were classified into 639 unique gene ontology functional annotation terms and were involved in 114 Kyoto Encyclopedia of Genes and Genomes pathways.
In terms of the current reports, cDNA microarrays and target metabolite analysis were performed on leaves of two different drought-tolerance potato clones (Solanum tuberosum subsp. andigenum) [2]. Gene expression analysis revealed that photosynthesis, photorespiration and carbohydrate-related gene, coupled with SPS and SS, showed more drastically down regulated in drought-tolerance clone than the expression in drought-sensitive clone. Zhang et al. [18] reported transcriptome dynamics of potato leaves under drought stress with high throughput sequencing platform and according to their report the DEGs were mainly annotated as genes coding for metabolism alteration, osmolite adjustment, cell rescue and gene regulation. Peripheral-type benzodiazepine receptor and polygalacturonase noncatalytic subunit AroGP3 were the one of the ten most up-regulated and the ten most down-regulated genes, respectively. Unlike previous studies, we took stolon tissues, one of the major agronomic traits, as sequencing material to examine transcriptome changes under two treatments and identified more significantly altered DEGs associated with drought stress (Tables 2, 3 and 4), starch synthesis and transfer (S4 Table), and water stimulus response (S6 Table). The difference of DEGs identified in these studies could be due to the different material properties, treatments, sampling strategies and genetic backgrounds of plants.
Many functional and regulatory genes are responsive to water deficit stress [19,20]. The main role of these functional genes is to maintain intracellular water and ion homeostasis, membrane structural stability, and reconstruction of primary and secondary metabolism [21], whereas regulatory genes (e.g., calcium ions, receptor protein kinases, and TFs) are involved in stress response by regulating signal transduction and metabolic pathways [22,23]. Heat shock responsive genes was one of the functional gene that facilitate protein refolding and stabilize polypeptides and membranes under stress condition by preventing protein from aggregation, protecting non-native enzymes from degradation and assisting protein refolding [24,25]. Expression induction of these genes under drought stress has been reported in barley [26] and rice [27]. Our data showed different expression of heat shock protein-encoding genes (Tables 2 and 3), indicating that they were also directly involved in regulating drought-stress responses in potato.
Aquaporins, which facilitate transmembrane transport of small molecules such as water, have been showed to be involved in response to drought and salt stress, and to ABA treatment [28,29]. In our data, aquaporin genes (Tables 2 and 3, S4 Table) were significantly expressed in Ningshu 4 stolons under drought stress condition. Similar expression pattern was also observed in Arabidopsis [28], chickpea (Cicer arietinum L.) [30], foxtail millet (Setaria italica L.) [31] and maize [32]. Additionally, ScPIP2a, an aquaporin, was reported to act in photoperiod perception, generating flower-inducing signals in the leaves [33]. The existing researches indicated these proteins could play an important role in water stimulus response and the maintenance of intracellular osmotic potential stability and organelle integrity. Therefore, we hypothesize that aquaporin(s) may play a similar regulatory role in potato tuber bulking and growth under drought stress.
Potato tuber formation is a complex biological process governed by both environmental factors and genes [34]. Cell division and expansion, together with the associated starch and protein accumulation, are the major events of tuber bulking in potato [35]. Photoperiod perception in leaves initiates tuberization in the subapical region of the underground stolons [36]. Dioxygenases, which catalyze the oxygenation of polyunsaturated fatty acids (LOXs) such as linoleic and linolenic acids, contribute to 9(S)-hydroperoxy linolenic acid (9(S)-HPOT) and 13(S)-hydroperoxy linolenic acid (13(S)-HPOT) accumulation [37]. Jasmonic acid is synthesized from 13(S)-HPOT and then metabolized to tuberonic acid and tuberonic acid glucoside, both of which induce tuber formation [38]. In this study, expression of PHYB, GA and LOXencoding gene (S4 Table) was significantly different in drought and re-watering plants, suggesting that these genes could function in regulating starch accumulation and tuber formation of potato plants.
As an important plant hormone, ABA plays important roles in plant response to drought stress by inducing the expression of TF-, heat shock protein-, transporter-, and osmotic regulator-encoding genes downstream of stress signaling pathways in both ABA-dependent and ABA-independent manners [19]. We observed that expression of ABA anabolism-or signal transduction-related genes such as abscisic acid 8'-hydroxylase 1-like, Beta-carotene hydroxylase, AREB, and PP2C5 was all significantly altered during drought stress and re-watering treatment (Table 5, S3 Table), indicating that ABA and anabolism related genes were also important regulators for drought-stress and re-watering stimulus-response in potato. This hypothesis is also supported by previous findings in corn (Zea mays) and cotton (Gossypium hirsutum) [24,32]. In addition, expression of some gene involved in auxin, cytokinin and gibberellin anabolism was significantly affected (S3 Table, S1 Fig), which, together with ABA, jointly regulates the response of potato stolons to drought stress and re-watering. Similar results were also found in Chinese cabbage (Brassica rapa ssp. pekinensis) [39] and barley (Hordeum vulgare) [40]. Collectively, our findings demonstrate the existence of a synergistic and antagonistic cross talk at endogenous hormone level in potato during the response to drought stress and re-watering treatment.
Many TFs, such as ERF [41], MYB [42], WRKY [43] and bHLH, act as key regulators in signal transduction pathways involved in plant response to drought stress [19,20]. Consistently, we observed four TF-encoding genes, WRKY, bHLH, ERF and MYB, were differentially expressed under drought stress and re-watering conditions (Tables 2 and 3, S3 Table). In particular, genes encoding WRKY, bHLH, and ERF showed opposite expression alterations under either condition, indicating that TF-encoding genes may have important functions in regulating the response to drought stress and re-watering stimulus in potato. Meanwhile, expression of genes mediated by them, such as those encoding transporter, osmotic adaptor, and those involved in fatty acid metabolism pathway (Tables 2 and 3, S3 Table), was also up-or down-regulated. Therefore, we postulate that drought stress response in potato stolons is similar to that in cotton as proposed by Padmalatha et al [24]. That is, potato stolons mediate the synthesis and signal transduction of ABA and the other endogenous hormones, and thereby influence the expression of downstream genes related to heat shock proteins, fatty acid metabolism and starch anabolism. Consequently, this mechanism further realizes drought-and water-stimulus response and regulation.
Based on DEGs detection, enriched biological processes and metabolism pathways analysis, our results provide a probable insight into the molecular mechanism in potato tuber response to drought stress and water-stimulus conditions. Our data shown in this work will provide more direct evidence and information for future study on the interaction between key genes involved in various metabolic pathways under drought stress condition. More detailed studies on the components involved in these pathways will help to understand the precise regulation process and mechanism related to drought stress response in potato, as well as for the molecular breeding of drought tolerant plants.