Construction of drought stress regulation networks in potato based on SMRT and RNA sequencing data
BMC Plant Biology volume 22, Article number: 381 (2022)
Potato (Solanum tuberosum) is the fourth most important food crop in the world and plays an important role in food security. Drought stress has a significantly negative impact on potato growth and production. There are several publications involved drought stress in potato, this research contributes to enrich the knowledge.
In this study, next-generation sequencing (NGS) and single-molecule real-time (SMRT) sequencing technology were used to study the transcription profiles in potato in response to 20%PEG6000 simulates drought stress. The leaves of the variety “Désirée” from in vitro plantlets after drought stress at six time points from 0 to 48 hours were used to perform NGS and SMRT sequencing. According to the sequencing data, a total of 12,798 differentially expressed genes (DEGs) were identified in six time points. The real-time (RT)-PCR results are significantly correlated with the sequencing data, confirming the accuracy of the sequencing data. Gene ontology and KEGG analysis show that these DEGs participate in response to drought stress through galactose metabolism, fatty acid metabolism, plant-pathogen interaction, glutathione metabolism and other pathways. Through the analysis of alternative splicing of 66,888 transcripts, the functional pathways of these transcripts were enriched, and 51,098 transcripts were newly discovered from alternative splicing events and 47,994 transcripts were functionally annotated. Moreover, 3445 lncRNAs were predicted and enrichment analysis of corresponding target genes was also performed. Additionally, Alternative polyadenylation was analyzed by TADIS, and 26,153 poly (A) sites from 13,010 genes were detected in the Iso-Seq data.
Our research greatly enhanced potato drought-induced gene annotations and provides transcriptome-wide insights into the molecular basis of potato drought resistance.
Drought is the most important environmental factor affecting crop yield. Studies have shown that drought can cause up to 50% of crop yield loss. As global warming intensifies the effect, frequency and range of drought, results in more serious yield loss . Therefore, development of breeding drought tolerant materials is an important goal of current breeding. Although lots of work on drought stress have been conducted in model plants, the molecular mechanism of drought response in potato, the fourth largest food crop, still needs to be further explored. The yield of potato is far from its physiological potential. Due to the drought-sensitivity of potato varieties, the yield is greatly reduced due to the lack of irrigation facilities or arid areas in the world. Drought stress leads to biochemical, molecular, physiological and morphological changes and significantly affects plant growth and yield . Materials with strong drought tolerance can reduce the impact of water shortage on plants. Therefore, screening and utilizing candidate drought stress response genes can stabilize potato yield in water-shortage area.
Plants have evolved very complex regulatory mechanisms to cope with the adverse environments [6, 51, 88]. Plants resist drought stress by reducing water loss (such as decreased leaf area, reduced stomata number and its conductance) and increasing water uptake (enhanced root formation, increased leaf thickness, and leaf rolling) in morphological aspect . Many studies have been conducted to explore the impact of drought stress on physiological, molecular, and biochemical levels in various crops . Aliche et al.  revealed the mechanism of potato carbon partition under drought stress. The development of new knowledge and techniques on omics high-throughput approaches is crucial to potato drought resistance , Obidiegwu et al.  presented an overview of drought stress and adaptive responses in potato, indicated that omics analysis was future research directions. Comparing transcript levels differences between drought-tolerant and sensitive cultivars can effectively understand the drought response mechanism of potato  .) identify key mechanisms and processes underlying single and combined abiotic stress tolerance by comparative analysis of tolerant and susceptible cultivars. Gervais et al.  has revealed that leaf protein content, instantaneous WUE, stomatal conductance, and transpiration rates can be used as the screening criteria for selecting the drought-tolerant potato genotypes. There have been many reports on the rapid identification of the molecular mechanisms of drought response at the whole-genome level using next-generation sequencing (NGS) platforms or gene microarrays in different species. Based on these data, key genes, proteins or pathways were identified. Genes encoding proteins which are involved in water transport (aquaporin) [2, 3], scavenging of free oxygen radicals (superoxide dismutase, catalase, and peroxidase), maintaining cellular membrane integrity (proline, mannitol, glycine, and betaine), and protecting macromolecules (chaperones and late embryogenesis abundant proteins) [21, 63] directly protects plants against drought stress. The other genes involved in signal perception or transduction may also play critical roles in response to drought stress. Key transcription factors (TFs) such as MYB, WRKY, bZIPs, kinases such as SnRK2 and small molecules such as abscisic acid (ABA) have been functionally proved to play key roles in plant drought tolerance (Shinozaki, Yamaguchi -[21, 64]). Some key TFs such as ABA-dependent MYC/MYB, AREB/ABF (ABA-responsive element binding/ABA-binding factor), DREB (ABA-independent dehydration-responsive element-binding proteins), CRT/DRE (C-repeat/drought-responsive element) have been proved in many crops in response to drought stress [28, 62, 82].
The S. tuberosum group Phureja clone DM1–3516 R 44 (doubled monoploid) and the S. tuberosum group Tuberosum RH89–039-16 (heterozygous diploid) have been sequenced and released via the Potato Genome Sequence Consortium (PGSC) in 2011 . There are new and improved sequences done for both diploid and tetraploid potato [55, 86]. These several genome sequences have been released and become in great genomic resources available, and provides a very good genomic reference information for transcriptome sequencing [19, 22, 46, 47]. For drought stress in potato, several studies have been conducted using RNA-Seq [10, 49, 66, 85]. In addition, Petek et al. compiled available and newly generated RNA-Seq datasets for three tetraploid potato genotypes (cultivar Désirée, cultivar Rywal, and breeding clone PW363) with diverse breeding pedigrees . Some studies on genomic assembly of tetraploid potato cultivar have also been reported recently [27, 69]. In this study, we assembled RNA-Seq datasets for Désirée in response drought stress.
Although NGS technology has the advantages of high throughput, high sensitivity, high accuracy and low cost, its main disadvantage is that the reading length is short, and it is difficult to ensure the accuracy of the reconstructed transcript in the assembly process. This presents a serious challenge to eukaryotic transcriptome sequencing with many homologous genes . In order to break the barrier, the third generation sequencing technology came into being. The third generation sequencing technology is a single molecule sequencing technology. In the process of sequencing, each RNA molecule is directly sequenced without PCR amplification. It has the characteristics of super long reading length and can cover the complete transcripts. The single molecular real time (SMRT) sequencing technology introduced by Pacific Bioscience (pacbio) was most widely used . Through this sequencing technology, information such as variable splicing, fusion gene, gene family and lncRNA can be accurately identified, which is widely used in genome research, transcriptome sequencing and epigenetic analysis . However, the main problem in SMRT is the high error rate (mainly reflected in insertion or deletion errors) . In order to solve this problem, the NGS with high quality and high coverage sequencing data was used to correct SMRT data, so “SMRT + NGS” sequencing joint analysis is more and more widely used [8, 83, 84].
In previous studies on potato drought stress transcriptome sequencing, although many differentially expressed genes were identified and screened, clustered and functionally annotated, their analysis was not comprehensive and detailed [16, 22, 85]. In our study, the leaves of the variety “Désirée” at the young plants stage after drought stress at six time points from 0 to 48 hours were used to perform NGS and SMRT sequencing. Our research will provide a valuable supplement to the potato genome and a broader regulatory mechanism in response to drought stress in potato.
Materials and methods
Plant materials and drought stress treatments
A drought tolerant tetraploid potato “Désirée ”[45, 49], provided by Southwest University, was used in this study. In 2 weeks old vitro plantlets of about 10 cm were transferred and stablished into hydroponic solution (1/4 Hoagland nutrient solution) for 1 week at 22/18 °C (16 h light / 8 h dark), then the materials were divided into two parts, one for control and the other for drought stress treatment to simulate the drought stress as reported by different authors, which were treated with 20% polyethylene glycol (PEG)-6000 [25, 48, 76]. The aboveground tissues were taken at 0, 1, 3, 6, 12, 24 and 48 hours after drought stress treatment and immediately frozen at least 30 min in liquid nitrogen and kept at − 80 °C until use. The samples were collected for three biological replicates. The collected samples were sent to Biomarker for full-field transcriptome sequencing and second-generation transcriptome sequencing, including sample detection, library construction and sequencing.
Library construction for Illumina and PacBio sequencing
Total RNA of all samples was extracted using EZ-10 DNAaway RNA Mini-Preps Kit (Sangon Biotech, Shanghai, China) and the quality of total RNAs was evaluated using an Agilent 2100 Bio-analyzer (SA Pathology, Adelaide, SA, Australia). High-quality RNAs (RIN > 8.0) were used for the following cDNA synthesis and library construction. For Illumina RNA sequencing, 1.0 μg RNA per sample was used for cDNA synthesis with the NEBNext® Ultra™ II RNA Library Prep Kit for Illumina (NEB, USA) and sent to Biomarker Technologies Corporation (Beijing, China) for Illumina sequencing. For SMRT sequencing, equal amount of the total RNA of all samples were pooled together as the template for cDNA synthesis with the SMARTer PCR cDNA Synthesis Kit (Clontech, USA) in Biomarker Technologies Corporation (Beijing, China).
Analysis of the Illumina and SMRT sequencing data
The raw Illumina data was filtered to generate clean reads using FastQC tool (http://www.bioinformatics.babraham.ac.uk/projects /fastqc/) (v0.11.5) and remove the adaptor contamination and low-quality reads using Trimmomatic (v0.35) . The high-quality clean reads were then mapped to the reference genome database (http://solanaceae.plantbiology.msu.edu/rh_potato_download.shtml) using HISAT2 version 2.1.0 and the default parameters . For SMRT sequencing data, the circular consensus sequences (CCS) were generated from the PacBio Iso-Seq raw reads with smrtanalysis_2.3.0 (http://www.pacb.com/) and then processed into error corrected reads of insert (ROIs) with a minimum full pass of > 3 and minimum accuracy of > 0.9. These ROIs were classified as full-length, non-chemiric (FLNC) transcripts using smrtanalysis_2.3.0 based on the presence and location of the 5′ and 3′ cDNA primers and poly (A) tail. All full-length reads were aligned to the S. tuberosum group Tuberosum RH89–039-16 reference genome (http://solanaceae.plantbiology.msu.edu/rh_potato_download.shtml) using GMAP software . Finally, all mapped reads were collapsed using a Python script in the ToFU package (https://github.com/PacificBiosciences/cDNA_primer/) with min-coverage = 85% and min-identity = 90%. All raw data was deposited in the National Center for Biotechnology Information (NCBI) sequence read archive with accession no. PRJNA728834.
Gene functional annotation, enrichment analysis, and weighted gene co-expression network (WGCNA) analysis
Gene functions was predicted using the following databases: NR (NCBI non-redundant protein sequences), Pfam (protein family), Swiss-Prot, KEGG (Kyoto Encyclopedia of Genes and Genomes), and GO (gene ontology). The GOseq package in R (version 3.4.3) and the KOBAS software were used to perform the GO and KEGG pathway enrichment, respectively. The gene expression levels were calculated using FPKM (fragments per kilobase of exon per million fragments mapped). Genes with FPKM values > 1 in each sample were used for WGCNA analysis using the R package WGCNA (v1.42).
Differential expression analysis and candidate gene determination
Differentially expressed genes (DEGs) were analyzed by DESeq and detected with the parameter FDR < 0.01 & |Log2Ratio| > 1. The plant transcription factors were predicted by iTAK software, and the differentially expressed transcription factors were screened by mapping to the DEGs database and the genes related to plant hormones in DEGs were screened based on Gene Ontology (GO) functional annotation. Twenty-nine DEGs were randomly selected and verified using the real-time (RT)-PCR to detect the accuracy of transcriptome data (Supplementary Table S1). Toolbox for Biologists v1.087  was used to construct the heatmap.
Identification of lncRNAs and AS events
In order to analyze the coding potential in the transcripts, four of the most widely used computational approaches including coding potential calculator (CPC) (Li, ), coding-non-coding Index (CNCI) , coding potential assessment tool (CPAT) , and pfam protein structure domain analysis were combined. Transcripts were regarded as lncRNA candidates using the parameter: lengths > 200 nt & > 2 exons.
To identify the AS (Alternative splicing) events, the ASTALAVISTA algorithm  based on the assembled GTF file from the PacBio Iso-Seq and Illumina data was conducted and five types of AS events including IR (Intron retention), ES (Exon skipping), A5 (Alternative 5′ splice site), A3 (Alternative 3′ splice site) and MEE (Mutually exclusive exon) were classified with the default settings.
Isoform identification and reverse transcription (RT)-PCR
To distinguish AS events, the reverse transcription (RT)-PCR was performed. Briefly, the RNA used for Iso-Seq was isolated using EZ-10 DNAaway RNA Mini-Preps Kit (Sangon, https://www.sangon.com) and synthesized to complementary DNA using the Hifair® II 1st Strand cDNA Synthesis SuperMix for qPCR gDNA digester plus (Yeasen). Isoform-specific primers were designed using Vector NTI software (Supplementary Table S1). PCR amplification was performed using 2 × Taq PCR Master Mix (K1034, https://www.apexbt.com/) and the PCR products were showed in agarose gel stained with GoodView™ Nucleic Acid Stain (Apr-13-2021, HGV-II, https://www.sbsbio.com).
BlazeTaq™ SYBR® Green qPCR Mix 2.0 (GeneCopoeia, USA) was used for RT-qPCR reactions on a CFX96 Real-time System (BIO-RAD, USA). The amplification program was briefly described in previous study . To ensure the precise and reproducible results, each sample was replicated three times and Ef1α in potato were used for each sample as an endogenous control (Supplementary Table S1). The formula F = 2 –ΔΔCt was used to calculate the relative expression levels of selected genes.
Alternative polyadenylation analysis
Alternative splicing and alternative polyadenylation (APA) of pre-mRNAs greatly contribute to transcriptome diversity, coding capacity of a genome and gene regulatory mechanisms in eukaryotes. TADIS pipeline (Transcriptome Analysis Pipeline for Isoform Sequencing) was used to further analyze the full-length non-chimeric sequence (FLNC) to identify APA  and MEME analysis was also performed to identify poly (A) sequence signals in our data.
Illumina and PacBio sequencing data analysis
To construct a comprehensive regulation network in response to drought stress, control and treated samples were processed through SMRT sequencing on the PacBio platform and next-generation sequencing (NGS) on the Illumina Hi-Seq platform. High-quality RNA was extracted from above tissues of “Désirée” after PEG-6000 treatments (0 h, 1 h, 3 h, 6 h, 12 h, 24 h and 48 h) for Illumina Hi-Seq and pool together for construction of Iso-Seq library. After removing low-quality reads and adapter sequences, 38,347,790-46,349,352 reads were obtained and more than 80% sequences were mapped to the reference genome (Table 1). The Pearson’s correlation coefficient showed that all correlation values between the three replicates ranged from 0.826 to 0.925, indicating a perfect positive correlation and that the sequencing results could be used for further analysis (Fig. 1). For Iso-Seq, the PacBio platform generated 33.24 Gb clean data. Among these data, 464,332 circular consensus (CCS) reads were obtained using Iso-seq pipeline with min full pass > 3 and min predicted accuracy > 0.9. Then, 375,387 full length non chimeric (FLNC) sequences were determined by searching for the polyA tail signal and the 5′ and 3′ cDNA primers and 121,434 high-quality consensus sequences were obtained by clustering the FLNC sequences. SMRT sequencing can offer better transcript integrity (up to 50 kb in length) than NGS and has greatly contributed in maize , populous  oilseed rape  and sorghum [1, 44], but no literature has been reported in potato. Significantly longer isoforms were detected by SMRT than by NGS assembly, demonstrating better isoform integrity (Fig. 2a). In addition, the overall expression levels of transcripts detected by SMRT were higher than those detected by NGS, indicating better quantification accuracy (Fig. 2b).
Analysis of gene expression and transcriptomic changes in response to drought treatment
DESeq was used for differential expression analysis between control versus treatments. False Discovery Rate (FDR) was used as statistical parameter for differential expression gene screening. Fold Change≥2 and FDR < 0.01 were used as screening criteria. In total, 31,637 DEGs were detected between control versus treatments, and the number of DEGs in D24h sample was the largest. (Fig. 3a and Supplementary Table S2). Moreover, almost all of DEGs in the libraries showed moderate expression levels with only a small percentage of the genes expressed in high levels (Fig. 3b and Supplementary Table S3). More than 88% of the DEGs had a − 3<fold change ≤3. Transcription factors (TFs) play important roles in plant abiotic stress responses (Fig. 3c) . To further analyze the function of the DEGs, we used iTAK software to predict transcription factors in DEGs. In total, 664 DEGs belonging to 28 TF families were identified. The TF family members were unevenly distributed. Seven TF families, bHLH (69), MYB (60), ERF (53), bZIP (46), NAC (43) WRKY (35) and HD-ZIP (33) accounted for more than 51% (339/664) of all TFs (Fig. 3d). In total, 919 DEGs belonging to eight classifications of plant hormones, and ABA, IAA and CK-related genes accounted for the majority, accounting for about 94% of all hormone-related genes (Fig. 3e). To combine the hormone profiling data with the expression profiles of genes involved in hormone metabolism and signaling, we analyzed 231 genes involved in plant hormone metabolism and visualized them as heat maps after functional classification (Fig. 4). Heat map results showed that DEGs were most involved in plant hormone signaling.
Gene profiling in response to drought stress in potato
To explore drought stress response mechanisms in potato, the expression patterns were calculated (Fig. 5a). The hierarchical cluster analysis was done and then the K means algorithm used for 98,624 genes. Fourteen clusters were obtained and 63 pathways were enriched, such as carbon metabolism, carbon fixation in photosynthetic organisms, porphyrin and chlorophyll metabolism, protein processing in endoplasmic reticulum and oxidative phosphorylation (Fig. 5b and c). Furthermore, WGCNA was also performed to screen genes response to drought stress by detecting molecular modules. Thirteen modules were generated and more than one module was significantly correlated with each sample. Blue, darkred, midnightblue, orange, cyan, lightgreen, and darkgreen were significantly correlated with DCK, D1h, D3h, D6h, D12h, D24h and D48h, respectively (Fig. 5d). Then, KEGG enrichment analysis was conducted to explore the enriched pathways in each module correlated to each sample. Galactose metabolism (ko00052) were enriched in more than two treated samples. Plant-pathogen interaction (ko04626) and glutathione metabolism (ko00480) were enriched in early stages after drought stress (Fig. 5e and Fig. S1). In addition, through the comparative analysis of the control group and the other six samples processed in different periods, a total of 12,798 differentially expressed genes after deduplication were found (Fig. S2 and Supplementary Table S4). The functional enrichment of these DEGs during the six drought treatment period was conducted (Fig. 6a) and 29 DEGs were selected to confirm the transcriptome sequencing data using the real-time (RT)-PCR (Fig. 6b).
Putative roles of AS for drought stress regulation in potato
AS plays critical roles in abiotic stress response by increasing genetic diversity in plants [36, 42]. However, little information about AS involved in drought stress response has been obtained in potato. In this study, 21,756 AS events were categorized into five patterns, including 12,453 retained intron (RI), 2493 skipping exon (SE), 236 mutually exclusive exons (MX), 4481 alternative 3′ splicing sites (A3), and 2093 alternative 5′ splicing sites (A5) using Isoform sequencing technology, which yields long reads without assembly (Fig. 7a). There is similar number of AS events for control and drought treatments, among these AS events, RI events predominated, accounting for 57.2%, followed by A3, SE, A5 and MX (Fig. 7b), which was consistent with findings of previous reports in other plant species . Four random selected genes with AS events were successfully validated by reverse transcription (RT)-PCR (Fig. 7c).
To identify the co-expression network of genes involved in AS events, the WGCNA package was used and three significant modules were identified (brown, blue and turquoise). The DCK sample correlated to the brown module (r = 0.79, P = 2e− 5), whereas the turquoise (r = 0.71, P = 3e− 4) and blue (r = 0.77, P = 5e− 5) modules were significantly correlated with samples after 6 and 24 h treatments, respectively (Fig. 7d). KEGG enrichment of the significant modules was also performed. Six, eleven and four pathways were significantly enriched in brown, turquoise and blue modules. The galactose metabolism (ko00052) was significantly enriched in turquoise and blue modules, which were treated with PEG-6000 after 6 and 24 h, respectively (Fig. 7e).
Putative roles of lncRNA involved in drought stress regulation in potato
LncRNA is a major component in transcriptome sequencing, and plays important roles in the regulation process after transcription [34, 39, 50, 58]. However, few information about lncRNA in potato has been obtained. In this study, 3345 lncRNAs with various lengths were identified using composite computational algorithm based on the PacBio sequencing data (Fig. 8a, b). More than 54% of lncRNAs were located in the intergenic region while less than 7% were located in the antisense region according to reference annotations (Fig. 8c). To identify the co-expression network of target genes of lncRNAs, the WGCNA package was used and yellow, brown, blue and turquoise modules were significantly enriched (Fig. 8d, Supplementary Table S5). The lncRNA target genes of KEGG analysis indicated that terpenoid backbone biosynthesis (ko00900) and phenylalanine metabolism (ko00360) pathway were most enriched in DCK and D6h; Glycerophospholipid metabolism (ko00564) and ubiquitin mediated proteolysis (ko04120) pathway were most enriched in D12h and D24h, respectively (Fig. 9).
Alternative polyadenylation analysis
In eukaryotes, polyadenylation is an important transcriptional modification in mRNAs [15, 61]. Alternative polyadenylation can easily increase the transcripts variety and regulate gene expression by producing different transcripts [18, 60, 62]. The iso-seq data allowed us to detect polyadenylation sites in potato using the program TAPIS . In this study, 26,153 poly (A) sites from 13,010 genes detected in the Iso-Seq data were identified (Supplementary Table S6). On average, 2.01 poly (A) sites per gene were found and 856 genes had at least five poly (A) sites (Fig. 10a). In addition, a clear nucleotide bias of the flanking region of all poly (A) sites was found. To confirm this phenomenon, we investigated the motifs of 50 nucleotides upstream of predominant sites using MEME analysis. Furthermore, three conserved polyadenylation signals, i.e. TTTTGT, GCCCCC and GGGGCG were overrepresented in upstream of poly (A) cleavage sites (Fig. 10b).
Potato has important economic value, while they are vulnerable to drought stress. Therefore, it is important to understand the drought resistance mechanism of potatoes. Transcriptome research is one of the indispensable tools for understanding life processes. However, RNA-Seq technology based on the second-generation high-throughput sequencing platform often cannot accurately obtain or assemble complete transcripts, and cannot identify Isoform, homologous genes. Full-length transcriptome sequencing based on PacBio SMRT does not need to interrupt RNA fragments, and directly reverse transcribes the full-length cDNA, which enables the feasible determination of the whole transcriptome at the isoform level . Previous study accelerated cloning of a potato late blight-resistance gene using RenSeq and SMRT sequencing . Recently, a chromosome-scale long-read reference assembly has been constructed . So far, a variety of potato species such as S. chacoense “M6” and S. pinnatisectum-derived somatic hybrid have been sequenced using Illumina platform [38, 71], and Kyriakidou et al. cultivated potato taxa using Illumina and long-read (PacBio) technologies . In our research, NGS and SMRT sequencing technologies were used to study the transcription profile of potato in response to drought stress. The leaves of young plants of the variety “Désirée” after drought stress were used for NGS and SMRT sequencing at six time points from 0 to 48 hours. Statistical analysis of sequencing data revealed the important roles of DEGs, alternative splicing, LncRNA and APA in potato under drought stress.
DEGs involved in drought stress response
In this study, SMRT and NGS technologies were combined to conduct a global survey of drought-treatment potato transcriptome. With the intensification of drought stress, the number of DEGs increased. The number reached the highest at 24 h after drought stress treatment (Fig. S1 and Supplementary Table S2). which was consistent with the findings previously reported in potato . Plant hormones, such abscisic acid (ABA), auxin (IAA), and cytokinin (CK), play vital roles in abiotic stresses . Transcription factors and plant hormone-related genes in differentially expressed genes were analyzed, and found that drought stress-related transcription factors such as bZIP and plant hormone-related genes such as ABA-related genes may played a dominant role. The functional annotation of differentially expressed genes confirmed their important role in drought stress, and 29 DEGs we verified using real-time (RT)-PCR (Supplementary Table S1 and Fig. 6b). Moreover, the WGCNA screening and KEGG enrichment analysis of drought stress response genes suggested that plant-pathogen interactions and glutathione metabolism pathways were enriched in the early stage of drought stress, and galactose metabolism pathways were enriched in at least two treatment samples, which indicated these pathways play an important role in the potato drought response process. Previous researches reported that galactose is involved in response to drought, high salt and low temperature tolerance , and transgenic studies demonstrated that glutathione helped Tortula ruralis adapt to rehydration following rapid desiccation . It was clear that the number of genes down regulated was greater that the upregulated which means that the plant stops regular active metabolism as main mechanism of response to the drought stress.
AS events involved in drought stress response
One of the advantages of PacBio sequencing is its ability to detect AS events. (Wang, [75, 87]). For different biological responses, AS is an effective mechanism to increase the complexity and flexibility of the transcriptome . Previous studies have confirmed that AS events play critical roles in potato growth and abiotic stress response [43, 52]. AS events identified during the potato drought stress response process were analyzed in this study. The proportion of different AS events was consistent with the findings previously reported in other plant species [80, 83]. Our experiment verified the AS events detected above, and the results showed that the size of the gel band fragments was consistent with the splice isoforms identified from the Iso-Seq data (Supplementary Table S1 and Fig. 7c) and the functional annotations of these genes also indicated that they were involved in drought stress response. In summary, we have verified that the four randomly selected genes undergo alternative splicing under drought stress, and their functional annotations are also related to plant stress tolerance. Our experiments results and sequencing data indicated that potatoes may increase their protein diversity under drought stress through alternative splicing, and thereby improve their environmental adaptability.
LncRNA plays an important role in drought stress response
As a type of non-coding RNA, lncRNA act as regulators in a series of biological processes . Previous research suggested that lncRNA has potential functions in potato defense response and development [34, 56]. In this study, a total of 3445 lncRNAs of different lengths were identified through four methods: CPC analysis, CNCI analysis, pfam protein domain analysis, and CPAT analysis (Fig. 8a and b). The WGCNA and KEGG analysis results of lncRNA target genes show the term ‘phenylalanine metabolism’ was enriched in the brown module, indicating that lncRNAs respond to drought stress by targeting genes related to amino acid metabolic pathways (Fig. 9), it also provided additional effective candidate genes for future functional characterization.
APA is essential for drought resistance of potato
Polyadenylation refers to the covalent linkage of polyadenylic acid to messenger RNA (mRNA) molecules. In the process of protein biosynthesis, this is part of the way to produce mature mRNA ready for translation. In eukaryotes, polyadenylation is a mechanism that interrupts mRNA molecules at their 3’ends. The polyadenylic acid tail (or poly A tail) protects mRNA from exonuclease attack, and is very important for the termination of transcription, export of mRNA from the nucleus, and translation. A comprehensive genome-wide APA map draft consisting of 26,153 poly (A) sites from 13,030 genes was constructed in this study (Supplementary Table S6). The nucleotide bias was evaluated and obvious nucleotide bias in the flanking regions of all poly (A) sites were found (Fig. 9b). These findings are consistent with those of previous studies in Brassica napus . Two conserved polyadenylation motifs in B. napus were identified, while three conserved polyadenylation motifs in potato were identified, indicating that the polyadenylation motifs are species-specific . Previous studies have shown that alternative polyadenylation of RNA is critical for gene function by increasing the complexity of the transcriptome and regulating gene expression .
Availability of data and materials
The datasets generated and analyzed during the current study are available in the NCBI Sequence Read Archive repository under Bioproject PRJNA728834 (https://www.ncbi.nlm.nih.gov/Traces/study/?acc=PRJNA728834).
Abdel-Ghany SE, Hamilton M, Jacobi JL, et al. A survey of the sorghum transcriptome using single-molecule long reads. Nat Commun. 2016;7:11706. https://doi.org/10.1038/ncomms11706.
Alexandersson E, Danielson JAH, Rade J, et al. Transcriptional regulation of aquaporins in accessions of Arabidopsis in response to drought stress. Plant J. 2010;61(4):650–60. https://doi.org/10.1111/j.1365-313X.2009.04087.x.
Alexandersson E, Fraysse L, Sjovall-Larsen S, et al. Whole gene family expression and drought stress regulation of aquaporins. Plant Mol Biol. 2005;59(3):469–84. https://doi.org/10.1007/s11103-005-0352-1.
Aliche EB, Theeuwen TPJM, Oortwijn M, et al. Carbon partitioning mechanisms in POTATO under drought stress. Plant Physiol Biochem. 2020;146:211–9. https://doi.org/10.1016/j.plaphy.2019.11.019.
Altschul SF, Madden TL, Schäffer AA, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25(17):3389–402. https://doi.org/10.1093/nar/25.17.3389.
Blum A. Crop responses to drought and the interpretation of adaptation. J Plant Growth Regul. 1996;20(2):135–48.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20. https://doi.org/10.1093/bioinformatics/btu170.
Chao Q, Gao ZF, Zhang D, et al. The developmental dynamics of the Populus stem transcriptome. Plant Biotechnol J. 2019;17(1):206–19. https://doi.org/10.1111/pbi.12958.
Chen C, Chen H, Zhang Y, et al. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020a;13(8):1194–202. https://doi.org/10.1016/j.molp.2020.06.009.
Chen YK, Li CH, Yi J, et al. Transcriptome response to drought, rehydration and re-dehydration in potato. Int J Mol Sci. 2020b;21(1). https://doi.org/10.3390/ijms21010159.
Dahal K, Li XQ, Tai H, et al. Improving potato stress tolerance and tuber yield under a climate change scenario – a current overview. Front Plant Sci. 2019;10:563. https://doi.org/10.3389/fpls.2019.00563.
Demirel U, Morris WL, Ducreux LJM, et al. Physiological, biochemical, and transcriptional responses to single and combined abiotic stress in stress-tolerant and stress-sensitive potato genotypes. Front Plant Sci. 2020;11:169. https://doi.org/10.3389/fpls.2020.00169.
Deng Y, Li J, Wu S, et al. Integrated NR database in protein annotation system and its localization. CE. 2006;32(5):71–4.
Dhindsa RS. Drought stress, enzymes of glutathione metabolism, oxidation injury, and protein synthesis in tortula ruralis. Plant Physiol. 1991;95(2):648–51. https://doi.org/10.1104/pp.95.2.648.
Elkon R, Ugalde AP, Agami R. Alternative cleavage and polyadenylation: extent, regulation and function. Nat Rev Genet. 2013;14(7):496-506. https://doi.org/10.1038/nrg3482.
Evers D, Lefevre I, Legay S, et al. Identification of drought-responsive compounds in potato through a combined transcriptomic and targeted metabolite approach. J Exp Bot. 2010;61(9):2327–43. https://doi.org/10.1093/jxb/erq060.
Foissac S, Sammeth M. ASTALAVISTA: dynamic and flexible analysis of alternative splicing events in custom gene datasets. Nucleic Acids Res. 2007;35:W297–9. https://doi.org/10.1093/nar/gkm311.
Fu HH, Yang DW, Su WY, et al. Genome-wide dynamics of alternative polyadenylation in rice. Genome Res. 2016;26(12):1753-60. https://doi.org/10.1101/gr.210757.116.
Galvez JH, Tai HH, Lague M, et al. The nitrogen responsive transcriptome in potato (Solanum tuberosum L.) reveals significant gene regulatory motifs. Sci Rep. 2016;6:26090. https://doi.org/10.1038/srep26090.
Gervais T, Creelman A, Li XQ, et al. Potato response to drought stress: physiological and growth basis. Front Plant Sci. 2021;12:698060. https://doi.org/10.3389/fpls.2021.698060.
Golldack D, Li C, Mohan H, et al. Tolerance to drought and salt stress in plants: unraveling the signaling networks. Front Plant Sci. 2014;5:151. https://doi.org/10.3389/fpls.2014.00151.
Gong L, Zhang H, Gan X, et al. Transcriptome profiling of the potato (Solanum tuberosum L.) plant under drought stress and water-stimulus conditions. PLoS One. 2015;10(5):e0128041. https://doi.org/10.1371/journal.pone.0128041.
Gong ZZ, Xiong LM, Shi HZ, et al. Plant abiotic stress response and nutrient use efficiency. Sci China Life Sci. 2020;63(5):635–74. https://doi.org/10.1007/s11427-020-1683-x.
Gupta N, Verma VK. Next-generation sequencing and its application: empowering in public health beyond reality. In: Arora P, editor. Microbial Technology for the Welfare of society. Microorganisms for sustainability, vol. 17. Singapore: Springer; 2019. https://doi.org/10.1007/978-981-13-8844-6_15.
Hella FA, El-Shabrawi HM, El-Hady MA, et al. Influence of PEG induced drought stress on molecular and biochemical constituents and seedling growth of Egyptian barley cultivars. J Genet Eng Biotechnol. 2018;16(1):203–12. https://doi.org/10.1016/j.jgeb.2017.10.009.
Heo JB, Lee YS, Sung S. Epigenetic regulation by long noncoding RNAs in plants. Chromosom Res. 2013;21(6–7):685–93. https://doi.org/10.1007/s10577-013-9392-6.
Hoopes G, Meng X, Hamilton JP, et al. Phased, chromosome-scale genome assemblies of tetraploid potato reveal a complex genome, transcriptome, and predicted proteome landscape underpinning genetic diversity. Mol Plant. 2022;15(3):520–36. https://doi.org/10.1016/j.molp.2022.01.003.
Hu HH, Xiong LZ. Genetic engineering and breeding of drought-resistant crops. Annu Rev Plant Biol. 2014;65:715–41. https://doi.org/10.1146/annurev-arplant-050213-040000.
Jian HJ, Ma JQ, Wei LJ, et al. Integrated mRNA, sRNA, and degradome sequencing reveal oilseed rape complex responses to Sclerotinia sclerotiorum (lib.) infection. Sci Rep. 2018;8(1):10987. https://doi.org/10.1038/s41598-018-29365-y.
Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28(11):1947–51. https://doi.org/10.1002/pro.3715.
Kanehisa M, Furumichi M, Sato Y. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49(D1):D545–51. https://doi.org/10.1093/nar/gkaa970.
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.
Kim D, Landmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357–60. https://doi.org/10.1038/nmeth.3317.
Kwenda S, Birch PR, Moleleki LN. Genome-wide identification of potato long intergenic noncoding RNAs responsive to Pectobacterium carotovorum subspecies brasiliense infection. BMC Genomics. 2016;17(1):614. https://doi.org/10.1186/s12864-016-2967-9.
Kyriakidou M, Achakkagari SR, Gálvez López JH, et al. Structural genome analysis in cultivated potato taxa. Theor Appl Genet. 2020;133(3):951–66. https://doi.org/10.1007/s00122-019-03519-6.
Laloum T, Martin G, Duque P. Alternative splicing control of abiotic stress responses. Trends Plant Sci. 2018;23(2):140–50. https://doi.org/10.1016/j.tplants.2017.09.019.
Lamaoui M, Jemo M, Datla R, et al. Heat and drought stresses in crops and approaches for their mitigation. Front Chem. 2018;6:26. https://doi.org/10.3389/fchem.2018.00026.
Leisner CP, Hamilton JP, Crisovan E, et al. Genome sequence of M6, a diploid inbred clone of the high-glycoalkaloid-producing tuber-bearing potato species Solanum chacoense, reveals residual heterozygosity. Plant J. 2018;94(3):562–70. https://doi.org/10.1111/tpj.13857.
Li MX, Xie S, Ma LP, et al. Transcriptome analysis reveals that long noncoding RNAs contribute to developmental differences between medium-sized ovarian follicles of Meishan and Duroc sows. Sci Rep. 2021;11(1):22510. https://doi.org/10.1038/s41598-021-01817-y.
Li XQ, Du D. Motif types, motif locations and base composition patterns around the RNA polyadenylation site in microorganisms, plants and animals. BMC Evol Biol. 2014;14:162. https://doi.org/10.1186/s12862-014-0162-7.
Li Y, Dai C, Hu C, Liu Z, et al. Global identification of alternative splicing via comparative analysis of SMRT- and Illumina-based RNA-seq in strawberry. Plant J. 2017;90(1):164–76. https://doi.org/10.1111/tpj.13462.
Ling Y, Alshareef S, Butt H, et al. Pre-mRNA splicing repression triggers abiotic stress signaling in plants. Plant J. 2017;89(2):291–309. https://doi.org/10.1111/tpj.13383.
Lu SW, Tian D, Borchardt-Wier HB, et al. Alternative splicing: a novel mechanism of regulation identified in the chorismate mutase gene of the potato cyst nematode Globodera rostochiensis. Mol Biochem Parasitol. 2008;162(1):1–15. https://doi.org/10.1016/j.molbiopara.2008.06.002.
Mace ES, Tai SS, Gilding EK, et al. Whole-genome sequencing reveals untapped genetic potential in Africa's indigenous cereal crop sorghum. Nat Commun. 2013;4:2320. https://doi.org/10.1038/ncomms3320.
Martínez I, Muñoz M, Acuña I, et al. Evaluating the drought tolerance of seven potato varieties on volcanic ash soils in a medium-term trial. Front Plant Sci. 2021;12:693060. https://doi.org/10.3389/fpls.2021.693060.
Massa AN, Childs KL, Buell CR. Abiotic and biotic stress responses in solanum tuberosum group phureja DM1-3 516 R44 as measured through whole transcriptome sequencing. Plant Genome. 2013;6:3. https://doi.org/10.3835/plantgenome2013.05.0014.
Massa AN, Childs KL, Lin HN, et al. The transcriptome of the reference potato genome solanum tuberosum group phureja clone DM1-3 516R44. PLoS One. 2011;6(10):e26801. https://doi.org/10.1371/journal.pone.0026801.
Meher PS, Shivakrishna P, Reddy KA, et al. Effect of PEG-6000 imposed drought stress on RNA content, relative water content (RWC), and chlorophyll content in peanut leaves and roots. Saudi J Biol Sci. 2018;25(2):285–9. https://doi.org/10.1016/j.sjbs.2017.04.008.
Moon KB, Ahn DJ, Park JS, et al. Transcriptome profiling and characterization of drought-tolerant potato plant (Solanum tuberosum L.). Mol Cells. 2018;41(11):979–92. https://doi.org/10.14348/molcells.2018.0312.
Muret K, Klopp C, Wucher V, et al. Long noncoding RNA repertoire in chicken liver and adipose tissue. Genet Sel Evol. 2017;49(1):6. https://doi.org/10.1186/s12711-016-0275-0.
Nejat N, Mantri N. Plant immune system: crosstalk between responses to biotic and abiotic stresses the missing link in understanding plant defence. Curr Issues Mol Biol. 2017;23:1–16. https://doi.org/10.21775/cimb.023.001.
Nicolas M, Rodriguez-Buey ML, Franco-Zorrilla JM, et al. A recently evolved alternative splice site in the BRANCHED1a gene controls potato plant architecture. Curr Biol. 2015;25(14):1799–809. https://doi.org/10.1016/j.cub.2015.05.053.
Obidiegwu JE, Sanetomo R, Flath K, et al. Genomic architecture of potato resistance to Synchytrium endobioticum disentangled using SSR markers and the 8.3k SolCAP SNP genotyping array. BMC Genet. 2015;16:38. https://doi.org/10.1186/s12863-015-0195-y.
Petek M, Zagorščak M, Ramšak Ž, et al. Cultivar-specific transcriptome and pan-transcriptome reconstruction of tetraploid potato. Sci Data. 2020;7(1):249. https://doi.org/10.1038/s41597-020-00581-4.
Pham GM, Hamilton JP, Wood JC, et al. Construction of a chromosome-scale long-read reference genome assembly for potato. Gigascience. 2020;9:9. https://doi.org/10.1093/gigascience/giaa100.
Ramirez Gonzales L, Shi L, Bergonzi SB, et al. Potato CYCLING DOF FACTOR 1 and its lncRNA counterpart StFLORE link tuber development and drought response. Plant J. 2021;105(4):855–69. https://doi.org/10.1111/tpj.15093.
Rhoads A, Au KF. PacBio sequencing and its applications. Genomics Proteomics Bioinformatics. 2015;13(5):278–89. https://doi.org/10.1016/j.gpb.2015.08.002.
Rinn JL, Chang HY. Genome regulation by long noncoding RNAs. Annu Rev Biochem. 2012;81:145–66. https://doi.org/10.1146/annurev-biochem-051410-092902.
Salehi-Lisar SY, Bakhshayeshan-Agdam H. Drought stress in plants: causes, consequences, and tolerance. In: Hossain M, Wani S, Bhattacharjee S, Burritt D, Tran LS, editors. Drought stress tolerance in plants, vol. 1. Cham: Springer; 2016. https://doi.org/10.1007/978-3-319-28899-4_1.
Sherstnev A, Duc C, Cole C, et al. Direct sequencing of Arabidopsis thaliana RNA reveals patterns of cleavage and polyadenylation. Nat Struct Mol Biol. 2012;19(8):845-52. https://doi.org/10.1038/nsmb.2345.
Shi YS. Alternative polyadenylation: new insights from global analyses. RNA. 2012;18(12):2105-17. https://doi.org/10.1261/rna.035899.112.
Shin D, Moon SJ, Han S, et al. Expression of StMYB1R-1, a novel potato single MYB-like domain transcription factor, increases drought tolerance. Plant Physiol. 2011;155(1):421–32. https://doi.org/10.1104/pp.110.163634.
Shinozaki K, Yamaguchi-Shinozaki K. Gene networks involved in drought stress response and tolerance. J Exp Bot. 2007;58(2):221–7. https://doi.org/10.1093/jxb/erl164.
Shinozaki K, Yamaguchi-Shinozaki K, Seki M. Regulatory network of gene expression in the drought and cold stress responses. Curr Opin Plant Biol. 2003;6(5):410–7. https://doi.org/10.1016/s1369-5266(03)00092-x.
Sicher RC, Timlin D, Bailey B. Responses of growth and primary metabolism of water-stressed barley roots to rehydration. J Plant Physiol. 2012;169(7):686–95. https://doi.org/10.1016/j.jplph.2012.01.002.
Sprenger H, Erban A, Seddig S, et al. Metabolite and transcript markers for the prediction of potato drought tolerance. Plant Biotechnol J. 2018;16(4):939–50. https://doi.org/10.1111/pbi.12840.
Sprenger H, Kurowsky C, Horn R, et al. The drought response of potato reference cultivars with contrasting tolerance. Plant Cell Environ. 2016;39(11):2370–89. https://doi.org/10.1111/pce.12780.
Steijger T, Abril JF, Engstrom PG, et al. Assessment of transcript reconstruction methods for RNA-seq. Nat Methods. 2013;10(12):1177–84. https://doi.org/10.1038/nmeth.2714.
Sun H, Jiao WB, Krause K, et al. Chromosome-scale and haplotype-resolved genome assembly of a tetraploid potato cultivar. Nat Genet. 2022;54(3):342–8. https://doi.org/10.1038/s41588-022-01015-0.
Taji T, Ohsumi C, Iuchi S, et al. Important roles of drought- and cold-inducible genes for galactinol synthase in stress tolerance in Arabidopsis thaliana. Plant J. 2002;29(4):417–26. https://doi.org/10.1046/j.0960-7412.2001.01227.x.
Tiwari JK, Rawat S, Luthra SK, et al. Genome sequence analysis provides insights on genomic variation and late blight resistance genes in potato somatic hybrid (parents and progeny). Mol Biol Rep. 2021;48(1):623–35. https://doi.org/10.1007/s11033-020-06106-x.
Verma V, Ravindran P, Kumar PP. Plant hormone-mediated regulation of stress responses. BMC Plant Biol. 2016;16:86. https://doi.org/10.1186/s12870-016-0771-y.
Wagner EM. Monitoring gene expression: quantitative real-time rt-PCR. Methods Mol Biol. 2013;1027:19–45. https://doi.org/10.1007/978-1-60327-369-5_2.
Wang B, Tseng E, Regulski M, et al. Unveiling the complexity of the maize transcriptome by single-molecule long-read sequencing. Nat Commun. 2016;7:11708. https://doi.org/10.1038/ncomms11708.
Wang M, Wang P, Liang F, et al. A global survey of alternative splicing in allopolyploid cotton: landscape, complexity and regulation. New Phytol. 2018;217(1):163–78. https://doi.org/10.1111/nph.14762.
Wei Q, Chen R, Wei X, et al. Genome-wide identification of R2R3-MYB family in wheat and functional characteristics of the abiotic stress responsive gene TaMYB344. BMC Genomics. 2020;21(1):792. https://doi.org/10.1186/s12864-020-07175-9.
Witek K, Jupe F, Witek A, et al. Accelerated cloning of a potato late blight–resistance gene using RenSeq and SMRT sequencing. Nat Biotechnol. 2016;34(6):656–60. https://doi.org/10.1038/nbt.3540.
Wu TD, Watanabe CK. GMAP: a genomic mapping and alignment program for mRNA and EST sequences. Bioinformatics. 2005;21(9):1859–75. https://doi.org/10.1093/bioinformatics/bti310.
Wu X, Liu M, Downie B, et al. Genome-wide landscape of polyadenylation in Arabidopsis provides evidence for extensive alternative polyadenylation. Proc Natl Acad Sci U S A. 2011;108(30):12533–8. https://doi.org/10.1073/pnas.1019732108.
Xie L, Teng K, Tan P, et al. PacBio single-molecule long-read sequencing shed new light on the transcripts and splice isoforms of the perennial ryegrass. Mol Gen Genomics. 2020;295(2):475–89. https://doi.org/10.1007/s00438-019-01635-y.
Xu X, Pan SK, Cheng SF, et al. Genome sequence and analysis of the tuber crop potato. Nature. 2011;475(7355):189–95. https://doi.org/10.1038/nature10158.
Yamaguchi-Shinozaki K, Shinozaki K. Transcriptional regulatory networks in cellular responses and tolerance to dehydration and cold stresses. Annu Rev Plant Biol. 2006;57:781–803. https://doi.org/10.1146/annurev.arplant.57.032905.105444.
Yan HL, Gao YW, Wu LL, et al. Potential use of the Pteris vittata arsenic hyperaccumulation-regulation network for phytoremediation. J Hazard Mater. 2019;368:386–96. https://doi.org/10.1016/j.jhazmat.2019.01.072.
Yao S, Liang F, Gill RA, et al. A global survey of the transcriptome of allopolyploid Brassica napus based on single-molecule long-read isoform sequencing and Illumina-based RNA sequencing data. Plant J. 2020;103(2):843–57. https://doi.org/10.1111/tpj.14754.
Zhang N, Liu B, Ma C, et al. Transcriptome characterization and sequencing-based identification of drought-responsive genes in potato. Mol Biol Rep. 2014;41(1):505–17. https://doi.org/10.1007/s11033-013-2886-7.
Zhou Q, Tang D, Huang W, et al. Haplotype-resolved genome analyses of a heterozygous diploid potato. Nat Genet. 2020;52(10):1018–23. https://doi.org/10.1038/s41588-020-0699-x.
Zhu CH, Li XF, Zheng JY. Transcriptome profiling using Illumina- and SMRT-based RNA-seq of hot pepper for in-depth understanding of genes involved in CMV infection. Gene. 2018;666:123–33. https://doi.org/10.1016/j.gene.2018.05.004.
Zhu JK. Abiotic stress signaling and responses in plants. Cell. 2016;167(2):313–24. https://doi.org/10.1016/j.cell.2016.08.029.
The authors thank South Western University (SWU) for providing all the facilities for the research work.
This study was supported by the National Key Research and Development Program of China (2018YFE0127900), the National Natural Science Foundation of China (32101659), the Science and Technology Partnership Program, Ministry of Science and Technology of China (KY201904016), the Talent Introduction Program of Southwest University Project (SWU019008) and the Germplasm Creation Research Program of Southwest University.
Ethics approval and consent to participate
The potato material used in this study has been approved by Southwest University, we comply with the IUCN Policy Statement on Research Involving Species at Risk of Extinction and the Convention on the Trade in Endangered Species of Wild Fauna and Flora, All methods for this study are in compliance with the provisions of Chinese law and international guidelines and regulations.
Consent for publication
The authors declared no conflict of interest in the authorship and publication of this document.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Volcano plots of expression comparisons in different drought treatment time. Supplementary Fig. S2. KEGG enrichment among the WGCNA modules most significantly correlated with drought-treated samples. The figure shows the heat map and KEGG enrichment of four drought-treated samples (D1h, D6h, D24h, D48h) [30,31,32]. Supplementary Table S1. Summary of primers used in this study. Supplementary Table S2. The number of differentially expressed genes. Supplementary Table S3. FPKM, fold changes, FDR and functional annotation of DEGs. Supplementary Table S4. Intersection data of DEGs in six samples. Supplementary Table S5. Prediction of lncRNA target genes based on co-expression. Supplementary Table S6. Statistics and localization of genes with poly (A) sites.
c. Original figures of gels.
About this article
Cite this article
Jian, H., Sun, H., Liu, R. et al. Construction of drought stress regulation networks in potato based on SMRT and RNA sequencing data. BMC Plant Biol 22, 381 (2022). https://doi.org/10.1186/s12870-022-03758-8
- Solanum tuberosum
- Drought stress
- Alternative splicing
- Alternative polyadenylation