Skip to main content

Genome-wide systematic survey and analysis of the RNA helicase gene family and their response to abiotic stress in sweetpotato

Abstract

Sweetpotato (Ipomoea batatas (L.) Lam.) holds a crucial position as one of the staple foods globally, however, its yields are frequently impacted by environmental stresses. In the realm of plant evolution and the response to abiotic stress, the RNA helicase family assumes a significant role. Despite this importance, a comprehensive understanding of the RNA helicase gene family in sweetpotato has been lacking. Therefore, we conducted a comprehensive genome-wide analysis of the sweetpotato RNA helicase family, encompassing aspects such as chromosome distribution, promoter elements, and motif compositions. This study aims to shed light on the intricate mechanisms underlying the stress responses and evolutionary adaptations in sweetpotato, thereby facilitating the development of strategies for enhancing its resilience and productivity. 300 RNA helicase genes were identified in sweetpotato and categorized into three subfamilies, namely IbDEAD, IbDEAH and IbDExDH. The collinearity relationship between the sweetpotato RNA helicase gene and 8 related homologous genes from other species was explored, providing a reliable foundation for further study of the sweetpotato RNA helicase gene family's evolution. Furthermore, through RNA-Seq analysis and qRT-PCR verification, it was observed that the expression of eight RNA helicase genes exhibited significant responsiveness to four abiotic stresses (cold, drought, heat, and salt) across various tissues of ten different sweetpotato varieties. Sweetpotato transgenic lines overexpressing the RNA helicase gene IbDExDH96 were generated using A.rhizogenes-mediated technology. This approach allowed for the preliminary investigation of the role of sweetpotato RNA helicase genes in the response to cold stress. Notably, the promoters of RNA helicase genes contained numerous cis-acting elements associated with temperature, hormone, and light response, highlighting their crucial role in sweetpotato abiotic stress response.

Peer Review reports

Introduction

RNA helicases are found in all organisms, including prokaryotes and eukaryotes, ranging from viruses to humans [1,2,3]. They play a crucial role in unwinding double-stranded RNA by utilizing the energy derived from NTP molecules [4]. During RNA metabolism, RNA molecules undergo various modifications [5]. However, their inherent instability can lead to disruptions in RNA metabolism, ultimately impacting the development and ability of plants to withstand abiotic stress [6,7,8].

Although structurally similar, RNA helicases exhibit diverse functions. These remarkable proteins play crucial roles in multiple biological processes [9], such as RNA splicing, RNA metabolism, ribosome formation, and nuclear-cytoplasmic transport. By actively participating in these intricate molecular mechanisms, RNA helicases contribute significantly to the proper functioning and regulation of cellular activities [10,11,12]. These genes have many functions in RNA metabolism, among which RNA helicase is associated with resistance to stress. The helicases are divided into 6 superfamilies, namely SF-1, -2, -3, -4, -5, and -6. Among them, the most representative and largest family is SF2. According to the change of DEAD (Asp-Glu-Ala-Asp) motif, the RNA helicase superfamily is DEAD, DEAH and DExD / H, respectively [13, 14]. Almost all the helicase proteins contain nine conserved motifs [15]. Each of these nine motifs plays distinct roles, which are essential for helicase enzymatic activities [16, 17].

RNA helicases are found in numerous plant species, including Arabidopsis thaliana [18], Oryza sativa [19], Ipomoea trifida [20], Glycine max [21], Zea mays [21], Gossypium spp [22], soybean [23], Gossypium raimondii [22] and Solanum lycopersicum [14]. A total of 32 DEAD were initially identified in Arabidopsis thaliana [24]. Then, 113 and 115 RNA helicase genes were identified in Arabidopsis and Oryza sativa [19]. Studies have shown that RNA helicase genes in Arabidopsis, Oryza sativa, Gossypium spp, Gossypium raimondii and Zea mays are divided into three subfamilies, the numeral of genes in apiece subfamily is as follows: DEAD-box (50, 51, 87, 51, and 57 genes), DEAH-box (40, 33, 48, 52, and 31 genes), and DExDH-box (71, 65, 78, 58, and 50 genes).

Sweetpotato is a valuable food source and industrial raw material, with high economic importance [25]. Being a hexaploid plant with 90 chromosomes, high heterozygosity, and abundant repetitive sequences, the identification and functional studies of genes in sweetpotato have been challenging [26, 27]. Given that sweetpotato is susceptible to abiotic stress, which significantly affects potato chip production [28], a genome-wide identification of sweetpotato RNA helicase genes is highly warranted. Genomic-assisted breeding technology can potentially improve sweetpotato yield by developing new or improved varieties. To explore the biological basis of cold resistance in sweetpotato, it is essential to identify differentially expressed genes in response to low temperature stress and apply them in production. Proteins and enzymes play critical roles in controlling numerous metabolic pathways by regulating RNA occurrence and metabolism. RNA helicases are involved in various molecular functions, including stress tolerance and development regulation. Therefore, their identification in sweetpotato and efforts to improve sweetpotato varieties hold significant practical value and bear great significance for enhancing the productivity and resilience of this important crop.

Therefore, this study seeks to elucidate the role and involvement pathway of RNA helicase genes in sweetpotato. The main objectives of this study are to conduct a thorough and comprehensive genome-wide validation of RNA helicase genes in sweetpotato. Additionally, we aim to analyze the intricate molecular mechanisms that underlie their active participation in various biological processes specific to sweetpotato. By successfully accomplishing these goals, our study aims to significantly enhance our understanding of the functional significance of RNA helicases in sweetpotato. This newfound knowledge will ultimately contribute to advancements in sweetpotato research and provide valuable insights for crop improvement efforts.

Materials and methods

Identification of the RNA helicase genes in sweetpotato genomes

To identify the members of the Arabidopsis RNA helicase gene family, we employed BLASTP to search for all known RNA helicase gene sequences of Arabidopsis (113) and rice (115) across multiple databases. Subsequently, utilizing this as a reference, duplicate genes were eliminated, and the DEAD-box RNA helicase family genes of sweetpotato were retrieved [29]. And all the information about the RNA helicase genes in Arabidopsis (https://www.arabidopsis.org/) and rice (http://rice.plantbiology.msu.edu/) was downloaded [30]. The whole sweetpotato genome sequence was derivative by Ipomoea genome Hub (https://ipomoea-genome.org) [31]. Subsequently, all protein were covered and each member of the RNA helicase gene was verified using the Pfam database (http://pfam.xfam.org/), the CD-search (https://www.ncbi.nlm.nih.gov/cdd/Structure/cdd/wrpsb.cgi), and the PROSITE (https://prosite.expasy.org/), and members lacking typical conserved RNA helicase domains were deleted. The sequence information of all sweetpotato RNA helicase proteins can be found in the Supplementary File 1.

Phylogenetic relationships of RNA helicase proteins in sweetpotato

The Clustal X program was utilized to compare the sequence of RNA helicase, while the MUSCLE program was employed for multiple sequence alignment to support Clustal X (http://www.clustal.org/) [32]. The maximum likelihood (ML) method was utilized to construct the phylogenetic tree. The optimal evolutionary model LG + G + F was applied to DEAD-box, while the DEAH-box and DExDH-box utilized the optimal evolutionary model WAG + G + F for phylogenetic tree construction. These methodological choices ensure robust and accurate representation of the evolutionary relationships among RNA helicase genes in sweetpotato. Additionally, the bootstrap value was set to 1000 for increased accuracy and reliability of the results [33].

Protein property and conserved domain of helicase genes in sweetpotato

The physicochemical properties of sweetpotato RNA helicase protein were predicted utilizing the ExPASy online database (http://expasy.org/). To predict the subcellular site and phosphorylation site of RNA helicase genes in sweetpotato, default parameters were utilized for both the Plant-mPLoc (http://www.csbio.sjtu.edu.cn/bioinf/plant-multi/) and NetPhos (http://www.cbs.dtu.dk/services/NetPhos/) databases, respectively. Furthermore, the gene structure of sweetpotato RNA helicase was obtained by comparing its gene sequence with the genome sequence, using Tbtools [34]. Finally, MEME (https://meme-suite.org/meme/tools/meme) was employed to determine the conserved domains [35], and the maximum number of bases was set to 15 in the program settings.

Chromosomal location and collinearity analysis of the RNA helicase genes in sweetpotato

The isoelectric points and molecular weights of these proteins were determined using ExPASy (http://expasy.org/) [36]. Additionally, the structural characteristics of these sweetpotato RNA helicase genes were analyzed in conjunction with genomic data. To investigate collinearity among RNA helicase genes and other plant genes, the genome sequence information of sweetpotato, Arabidopsis thaliana, and rice was downloaded and examined. MCScanX was employed to establish gene duplication and collinearity relationships, utilizing default parameters [37]. The resulting data was visualized using TBtools [34, 38]. Default parameters were utilized throughout all steps of the analysis.

Plant materials and sample collection

The 11 sweetpotato varieties were selected as plant materials, which were collected and conserved in the Xuzhou Sweetpotato Research Center in China, and no permission was required for plant collection. Stem segments of XuShu18, exhibiting stable growth and measuring about 15 cm in length, were chosen for abiotic stress and hormone treatment, subsequently, the treated leaves were collected to detect the expression of related DEAD-box RNA helicase genes. Young leaves, mature leaves (leaves), stems, and roots were collected from ten sweetpotato varieties (PuShu 32, AoZhouZiBai, LongShu 9, XinXiang, XuZiShu 8, ShangShu 19, JiShu 26, YanShu 25, Marsali, and Hami) during the same period to analyze their tissue-specific expression.

Expression analyses of the sweetpotato RNA helicase genes

The RNA helicase gene of sweetpotato was analyzed in conjunction with the published transcriptome results under cold treatment on XuShu18 (Xie et. al (2019), the supplementary data can be found online at https://doi.org/https://doi.org/10.1016/j.ygeno.2018.05.014). XuShu18 was exposed to four abiotic stress treatments and four hormonal stress treatments at time points of 0 h, 1 h, 3 h, 6 h, 12 h, 24 h, 48 h, and 72 h. Four plant tissue parts (YL, L, S, and R) were collected from ten sweetpotato varieties, the objective was to examine the tissue-specific expression of the DEAD-box RNA helicase gene. RNA was extracted from these samples for use in subsequent experiments.

The total RNA was excavated using RNA extraction kit (TianGen, Beijing, China), and reverse transcriptions using TransScript® gDNA removal (TransGen, Beijing, China). All the sweetpotato RNA helicase promoter regions were examined by plantCARE (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/).

A. rhizogenes-mediated sweetpotato transformation

ZiShu 8 was chosen as the experimental material, and the following methods were employed for mediating penetration. The constructed IbDExDH96 overexpression vector (OE-IbDExDH96) and the control vector (CK) were introduced into A.rhizogenes K599 and cultured in 30 mL LB liquid medium containing 20 mg/L rifampicin and 50 mg/L kanamycin at 28 °C for 12 -16 h. After the OD600 value reached between 0.4—0.6, the solution was centrifuged at 4000 rpm for 10 min, and then resuspended in the infiltration buffer. Subsequently, 1 mL of the rhizobium suspension was injected into the stems of sweetpotato cuttings using a needle and syringe. Following this, the cuttings were transplanted into the soil and grown in a high-humidity environment [39]. After 3 weeks of cultivation in the greenhouse, transgenic positive plants were collected for further analysis.

Statistical analysis

Statistical analysis was conducted using Microsoft Excel 2019, Graphpad Prism 5.0, and SPSS statistical software (IBM Corp, Armonk, NY, USA). In order to ensure the biological significance of the results, a two-fold cut-off value was applied for differential gene expression analysis. The data were subjected to ANOVA (one-way analysis of variance) using the SPSS software, and differences in means were considered significant if they had a p-value < 0.05, as determined by Dunnett's test.

Results

Identification of RNA helicase family genes in sweetpotato

To identify the components of the sweetpotato RNA helicase gene, we utilized bioinformatics methods to gather information on numerous sweetpotato RNA helicases. A total of 300 sweetpotato RNA helicase proteins were identified. Based on conserved motifs within the RNA helicases, the 300 sweetpotato RNA helicase genes were categorized into the DEAD-box subfamily (53 genes), the DEAH-box subfamily (54 genes), and the DExDH-box subfamily (193 genes) (Supplementary File 1 and Supplementary Fig. 1). The amount of RNA helicase proteins in Arabidopsis and rice was 113 and 115, respectively. It is predicted that there are 300 RNA helicase proteins in sweetpotato, which is much more than the number of Arabidopsis and rice [19]. Subsequently, the 300 RNA helicase genes located on the 15 sweetpotato chromosomes were systematically named from top to bottom as IbDEAD1 ~ IbDEAD53, IbDEAH1 ~ IbDEAH54, and IbDExDH1 ~ IbDExDH193, following a consistent nomenclature convention (Supplementary Fig. 2). Subsequently, the physicochemical properties of 300 sweetpotato RNA helicase proteins were examined. The length and relative molecular mass of the RNA helicase vary greatly. The length of DEAD ranges from 323 aa (IbDEAD34) to 1301 aa (IbDEAD2), the relative molecular mass ranges from 3584.09 to 145,261.6 Da, and the isoelectric point ranges from 5.02 (IbDEAD40) to 9.87 (IbDEAD37). The length of DEAH is between 206 aa (IbDEAH5) and 2904 aa (IbDEAH44), the relative molecular mass is between 23,377.03 and 323,370.3 Da, and the isoelectric point is between 5.16 (IbDEAH22) and 9.25 (IbDEAH36). The length of DExDH ranged from 128aa (IbDExDH7) to 2801aa (IbDExDH30), the relative molecular mass ranged from 13,904.77 to 306,450.01 Da, and the isoelectric point graded from 4.81 (IbDExDH7) to 9.67 (IbDExDH45). The subcellular location showed that most of the RNA helicase proteins were positioned in the nucleus. Furthermore, the potential phosphorylation sites showed that IbDEAD contains 29 (IbDEAD28) to 141 (IbDEAD2) phosphorylation sites, IbDEAH contains 19 (IbDEAH5) to 280 (IbDEAH44) phosphorylation sites, and IbDExDH contains 18 (IbDExDH141) to 362 (IbDExDH30) phosphorylation sites, all of these sweetpotato RNA helicase proteins restrain more Ser sites than the Tyr and Thr sites (Supplementary Table 1).

Phylogenetic analysis of the RNA helicase family proteins in sweetpotato

To investigate the evolutionary patterns and classification of RNA helicases in sweetpotato, we constructed a rootless phylogenetic tree, comparing the identified RNA helicases in sweetpotato with the known RNA helicases in Arabidopsis thaliana (Fig. 1). The RNA helicases in Arabidopsis were classified into three subfamilies: DEAD, DEAH, and DExDH. Among the 53 IbDEAD proteins, except for IbDEAD5, they were further divided into nine subgroups. Similarly, the 54 IbDEAH proteins, excluding IbDEAH36, were separated into nine subgroups. Furthermore, the 193 IbDExDH proteins were categorized into 13 subgroups, with the exception of IbDExDH19 and IbDExDH23. Notably, the V subgroup in IbDEAH and the XII subgroup in IbDExDH were unique to sweetpotato. The distribution of RNA helicase proteins varied across the different subfamilies. Subgroups VI, VIII, and XI of IbDExDH contained the fewest sweetpotato RNA helicase genes. Intriguingly, IbDEAD5, IbDEAH36, IbDExDH19, and IbDExDH23 did not belong to any of the aforementioned three subfamilies, suggesting that they might possess alternative functions.

Fig. 1
figure 1

Phylogenetic tree analysis of RNA helicase in Arabidopsis and sweetpotato. A The DEAD-box RNA helicase proteins in two species. B The DEAH-box RNA helicase proteins in two species. C The DExDH-box RNA helicase proteins in two species. These three subfamilies were used to construct the maximum-likelihood (ML) phylogenetic tree by MEGA X with 1,000 bootstrap replicates. Different subgroups are named based on the reports in Arabidopsis and are distinguished with different colors. The gene names are marked at the end of the branch, the red triangle and green circle represent the sweetpotato RNA helicase and Arabidopsis RNA helicase, respectively

Gene structure and motif composition analyses in sweetpotato

To gain further insights into the structural characteristics of the sweetpotato RNA helicase family, we conducted a comparative analysis of conserved motifs and intron/exon composition based on the phylogenetic tree of sweetpotato RNA helicases (Figs. 2 and 4). The structural features of IbDEAD, IbDEAH, and IbDExDH subfamily genes exhibit notable variations across different subgroups. However, it is worth mentioning that genes within the same side branch are homologous, sharing similar genetic structures. This observation suggests a strong correlation between the exon–intron structure and phylogenetic relationships within the sweetpotato RNA helicase family (Figs. 2, 3 and 4). The genetic architecture of the RNA helicase family members in sweetpotato is remarkably intricate, with the IbDEAD and IbDEAH gene members exhibiting a higher number of exons. In contrast, among the IbDExDH gene members, only a small fraction of genes (2.6%) were found to comprise a single exon, while the rest were composed of multiple exons. To gain further insight into the structural conservation of these RNA helicases, we performed motif analysis for each member using the MEME website, revealing the presence of 15 highly conserved motifs across all sweetpotato RNA helicases. This finding underscores the potential functional significance of these motifs in regulating RNA helicase activity in sweetpotato (Figs. 2, 3 and 4). Sweetpotato RNA helicases exhibit 15 conserved motifs, with motifs 1, 2, 3, 4, 5, 7, 8, 9, and 15 being widely present in most proteins. Motif 5 is particularly noteworthy because it contains highly conserved sequences, including DEAD, DEAH, and DExDH, which are essential for RNA helicase function. Although these conserved motifs are shared among sweetpotato RNA helicases, significant variations exist in their domains and amino acid sequences. These differences suggest that despite having common structural features, sweetpotato RNA helicases may have diverse functions and regulatory mechanisms (Figs. 2, 3 and 4). In general, phylogenetic tree examination of data shows that the system development and features and the divergence of genetic structure and sequence distribution are closely related.

Fig. 2
figure 2

Phylogenetic tree, gene structure and motif distribution analysis of 53 IbDEAD-box members. A Phylogenetic tree of DEAD-box members in sweetpotato. A total of 53 IbDEAD-box proteins were used to construct the maximum-likelihood (ML) phylogenetic tree by MEGA X with 1,000 bootstrap replicates. The different groups are marked with different colors. B The exon–intron structure analysis of 53 IbDEAD-box genes. Exon and introns lengths are displayed proportionally, yellow boxes represent exons, black lines represent introns, and green boxes represent 5’ and 3’ UTR. C Motif distribution of 53 IbDEAD-box proteins. The different colored boxes represent 15 motifs, and the box lengths represent motif lengths

Fig. 3
figure 3

Phylogenetic tree, gene structure and motif distribution analysis of 54 IbDEAH-box members. A Phylogenetic tree of DEAH-box members in sweetpotato and the consistent parameters shown in Fig. 2(A). The different groups are marked with different colors. B The exon–intron structure analysis of 54 IbDEAH-box genes and the consistent parameters shown in Fig. 2(B). C Motif distribution of 54 IbDEAH-box proteins. The different colored boxes represent 15 motifs, and the box lengths represent motif lengths

Fig. 4
figure 4

Phylogenetic tree, gene structure and motif distribution analysis of 193 IbDExDH-box members. A Phylogenetic tree of DExDH-box members in sweetpotato and the consistent parameters shown in Fig. 2(A). The different groups are marked with different colors. B The exon–intron structure analysis of 193 IbDExDH-box genes and the consistent parameters shown in Fig. 2(B). C Motif distribution of 103 IbDExDH-box proteins. The different colored boxes represent 15 motifs, and the box lengths represent motif lengths

Chromosome localization of the RNA helicase family in sweetpotato

Using GFF3 genome annotation, physical location detection analysis revealed that 53 IbDEAD genes, 54 IbDEAH genes, and 193 IbDExDH genes were distributed across all 15 sweetpotato chromosomes. Notably, chromosome 6 harbored the highest number of IbDEAD genes, with 7 members, while chromosome 11 contained the largest number of IbDEAH genes, with 9 members. Strikingly, chromosome 9 had no IbDEAD or IbDEAH genes. Furthermore, most chromosomes exhibited a high abundance of IbDExDH genes, with the exception of chromosome 10, which contained only four IbDExDH genes. These findings shed light on the genomic distribution of RNA helicases in sweetpotato, suggesting potential functional diversification among different gene family members, and providing a foundation for further investigations into their biological roles and regulatory mechanisms (Figs. 5 and Supplementary Fig. 3). The distribution of the three subfamilies of RNA helicases across the 15 chromosomes exhibits notable divergence. Collinearity analysis revealed a cluster of tandem duplicated IbDEAH genes and two clusters of tandem duplicated IbDExDH genes, specifically IbDEAH2 / 24 and IbDExDH114 / 168 as well as IbDExDH162 / 185. However, no tandem duplicated genes were identified in the IbDEAD subfamily (Supplementary Table 7). These findings provide insights into the evolutionary dynamics and genomic organization of RNA helicase genes in sweetpotato, suggesting potential mechanisms underlying their expansion and functional diversification.

Fig. 5
figure 5

The interchromosomal relationship of RNA helicase genes in sweetpotato chromosomes. A IbDEAH, B IbDExDH. Circular visualizations of RNA helicase genes mapped to the LG1- LG15 chromosomes are indicated by colored rectangles. The colored curves represent duplicated RNA helicase gene pairs. The corresponding RNA helicase genes located in segmental duplications are marked with colors

Cis element analysis of the RNA helicase family gene promoters

To investigate the potential regulatory mechanisms of sweetpotato RNA helicases in response to abiotic stress and hormones, we conducted a comprehensive analysis by scanning the cis-acting elements present in the 2 kb upstream promoter region of these genes. This analysis was performed using the PlantCare database, which provides valuable insights into the putative regulatory elements involved in the modulation of gene expression (Figs. 6, 7 and 8, Supplementary Table 2).

Fig. 6
figure 6

Phylogenetic clustering and predicted stress- and hormone-related cis-elements in the promoters of IbDEAD genes. A The phylogenetic tree of 53 IbDEADs was constructed by MEGA X based on the consistent parameters used in Fig. 2(A). B Predicted cis-elements in the IbDEAD promoters. 2000 bp promoter regions of each IbDEAD gene were detected by PlantCARE database. Different colored rectangles represent different cis-elements that are potentially involved in stress or hormone regulation

Fig. 7
figure 7

Phylogenetic clustering and predicted stress- and hormone-related cis-elements in the promoters of IbDEAH genes. A The phylogenetic tree of 54 species of IbDEAH and the consistent parameters shown in Fig. 2(A). B Predicted cis-elements in the IbDEAH promoters. Different colored rectangles represent different cis-elements that are potentially involved in stress or hormone regulation

Fig. 8
figure 8

Phylogenetic clustering and predicted stress- and hormone-related cis-elements in the promoters of IbDExDH genes. A The phylogenetic tree of 193 species of IbDExDH and the consistent parameters shown in Fig. 2(A). B Predicted cis-elements in the IbDExDH promoters. Different colored rectangles represent different cis-elements that are potentially involved in stress or hormone regulation

The conclusions showed that each promoter region of RNA helicase had multiple cis-acting elements connected with stress- and / or hormones. Among others, nearly 90% of RNA helicase promoters contain multiple stress cis-elements, for instance stress response elements (TC-rich repeats), low temperature response elements (LTR), MAJA response elements (TGACG-motif), drought response elements (MBS), etc. These cis-acting elements can be associated with expression. For example, the expression of several RNA helicase genes providing IbDEAH-32 / -42 and IbDExDH-36 / -47 / -96 was increased under different stresses. Correspondingly, stress-related repeats of MBS, TC, or LTR cis elements are more numerous in their promoter regions. However, TC-rich repeats, MAJA response elements and LTR elements were found on the promoters of the IbDExDH-25 / -48 and IbDEAH53 genes, their expression was not particularly significant under salt, drought, or cold stress, but there was a certain degree of response under high-temperature stress (Fig. 9 and Supplementary Table 5). Furthermore, all sweetpotato RNA helicase promoters restrain many hormone elements, such as abscisic acid response element (ABRE), Me-JA response element (CGTCA motif and TGACG motif) or auxin response element (TGA-box) (Figs. 6, 7 and 8, Supplementary Table 2). Nevertheless, the expression of the RNA helicase genes in dissimilar tissues of 10 different sweetpotato varieties was also different (Fig. 14). Most of the RNA helicase genes are communicated in high amounts in sweetpotato stems and root tissues, in particular, the expression levels of IbDExDH36 and IbDExDH48 in the roots were lower than those in other tissues. It shows that RNA helicase is associated to plant development. Among them, the transcription levels of most RNA helicase genes do not respond significantly to hormonal treatment (Fig. 10). These data indicate that the cis-acting elements of sweetpotato RNA helicase can be concerned in both hormonal and abiotic stresses.

Fig. 9
figure 9

Relative expression levels of 8 RNA helicase genes in response to abiotic stresses detected by qRT-PCR. The abiotic stress treatments include cold (4 °C), heat (42 °C), salt (200 mM NaCl) and drought (200 mM PEG4000). The expression levels at 0 h were normalized to 1, and the Y-axis delineates the fold changes of relative expression comparing with 0 h. Bars represent the mean of three biological replicates ± SE. The two-fold threshold is presented by a dotted line. Values are presented relative to untreated plants (0 h). Asterisks indicate a significant difference (P < 0.05) between untreated and treated plants

Fig. 10
figure 10

Relative expression levels detected by qRT-PCR under different hormone treatments including GA, IAA, ABA and ZT. The expression levels at 0 h were normalized to 1, and the Y-axis delineates the fold changes of relative expression comparing with 0 h. Bars represent the mean of three biological replicates ± SE. The two-fold threshold is presented by a dotted line. Values are presented relative to untreated plants (0 h). Asterisks indicate a significant difference (P < 0.05) between untreated and treated plants

Identification of RNA helicase family cold-response genes in transcriptomes and their expression profiles under multiple abiotic stress and hormonal treatments

Numerous studies have demonstrated the significant involvement of RNA helicases in diverse abiotic stresses, including cold, salt, drought, and heat, as well as developmental processes. In order to elucidate the potential biological role of the sweetpotato RNA helicase gene under adverse stress conditions, we conducted a comprehensive analysis using our previously generated RNA-seq data. Specifically, we focused on investigating the expression pattern of the XuShu 18 RNA helicase gene under cold stress [40]. The consequence indicated that eight RNA helicase genes were screened out. Subsequently, we further studied the expression patterns of eight genes (IbDEAH-32 / -42 / -53, IbDExDH-25 / -36 / -47 / -48 / -96) screened by qRT-PCR under salt, drought, heat and cold stresses, and explore a two-fold cut off value [41]. The results showed that all selected genes were up-regulated to varying degrees after salt, drought, cold, and heat treatment. Among them, four stresses could up-regulate IbDEAH32 and IbDExDH-36 / -47 transcription, two stresses could up-regulate IbDEAH42 and IbDExDH-48 / -96 expression, and one stress could up-regulate IbDExDH25 and IbDEAH53 transcription. All RNA helicases could be up-regulated under high temperature stress and salt stress. In particular, IbDEAH32 and IbDExDH47 had the highest level of induction after cold treatment, about six times, while the level of induction was relatively low in the transcription of other RNA helicases. Except for IbDExDH48, the expression of other genes was significantly enhanced under high temperature conditions, which was 2.3 and 19.3-fold that of the control. Under salt stress conditions, the expression of IbDEAH42 and IbDExDH-47 / -48 was enhanced, and the expression of other genes was weakly induced by salt stress. In particular, the response of all genes to drought stress was not very strong, which may be affected by other factors (Fig. 9). Overall, these data suggest that multiple sweetpotato RNA helicase members can take the lead in reaction to abiotic stress.

Moreover, qRT-PCR was acclimated to further detect the transcription profiles of eight RNA helicase genes under distinct hormone treatments, encircling ABA, IAA, GA and ZT. Unexpectedly, most of the RNA helicase genes were down-regulated when we used double as the cut-off value for differential expression. Just the stress hormone ZT could prompt the expression of IbDExDH-36 / -47 / -48 / -98 and IbDEAH42 (Fig. 10). It is reported that the response of the RNA helicase to hormone treatment is not obvious, mainly related to development and abiotic stress [42, 43]. Overall, these data propose that multiple members of the sweetpotato RNA helicase can be important players in answer to hormone and / or abiotic stresses.

Collinearity analysis of the RNA helicase genes between sweetpotato and other plants

To furthermore explore the origin and evolutionary mechanism of sweetpotato RNA helicase genes, we compared the homology of 300 RNA helicase genes with 8 representative species-related genes. These species include wild diploid relatives of sweetpotato (Ipomoea trioba and Ipomoea trifida), two model plants (Arabidopsis thaliana and Oryza sativa), two cruciferous plants (Brassica rapa and Brassica oleracea) and two Solanaceae plants (Solanum lycopersicum and Capsicum annuum). Among them, 42 (79.2%) and 43 (81.1%) IbDEAD genes were homologous to genes in Ipomoea trioba and Ipomoea trifida, respectively, accompanied by Solanum lycopersicum (8), Capsicum annuum (6), Arabidopsis thaliana (6), Brassica rapa (2), and Brassica oleracea (1), but no homologous genes were found both sweetpotato and rice. 39 (72.2%) and 40 (74.1%) IbDEAH genes were also homologous to the genes in Ipomoea trioba and Ipomoea trifida, respectively, accompanied by Solanum lycopersicum (8), Capsicum annuum (3) and Arabidopsis thaliana (3), nevertheless, no similar homologous genes were noted between sweetpotato and Brassica rapa, cabbage and rice. Similarly, 89 (46.1%) and 87 (45.1%) IbDExDH genes were homologous to the genes in Ipomoea trioba and Ipomoea trifida, followed by tomato (24), pepper (10), Arabidopsis (7), Brassica rapa (3), and cabbage (2), however, no homologous genes were found with rice (Figs. 11, 12 and 13). It should be mentioned that the collinearity of sweetpotato RNA helicase genes between the Ipomoea trioba and Ipomoea trifida genes more than the extra six varieties, which can be associated to the wild diploid relationship of sweetpotato.

Fig. 11
figure 11

Synteny analyses of IbDEAD genes between sweetpotato and eight representative plant species from Ipomoea triloba and Ipomoea trifida A Solanum lycopersicum and Capsicum annuum B Brassica rapa and Brassica oleracea C and Arabidopsis thaliana and Oryza sativa D The chromosomes of different plants are distinguished with differential colors. The blue lines connecting two different chromosomes highlight the syntenic IbDEAD gene pairs within sweetpotato and other plant genomes, respectively

Fig. 12
figure 12

Synteny analyses of IbDEAH genes between sweetpotato and eight representative plant species from Ipomoea triloba and Ipomoea trifida A Solanum lycopersicum and Capsicum annuum B Brassica rapa and Brassica oleracea C and Arabidopsis thaliana and Oryza sativa D The chromosomes of different plants are distinguished with differential colors. The blue lines connecting two different chromosomes highlight the syntenic IbDEAH gene pairs within sweetpotato and other plant genomes, respectively

Fig. 13
figure 13

Synteny analyses of IbDExDH genes between sweetpotato and eight representative plant species from Ipomoea triloba and Ipomoea trifida A Solanum lycopersicum and Capsicum annuum B Brassica rapa and Brassica oleracea C and Arabidopsis thaliana and Oryza sativa D The chromosomes of different plants are distinguished with differential colors. The blue lines connecting two different chromosomes highlight the syntenic IbDExDH gene pairs within sweetpotato and other plant genomes, respectively

Moreover, the mass of the RNA helicase genes in our sweetpotato are one-to-one gene pairs with sweetpotato diploids, and a few are one sweetpotato gene corresponding to two or three diploid gene pairs. For example, the sweetpotato genes IbDEAD (10), IbDEAH (2) and IbDExDH (19) have a one-to-many relationship in Ipomoea trioba, and IbDEAD (10), IbDEAH (2) and IbDExDH (17) sweetpotato genes have a one-to-many relationship in Ipomoea trifida (Supplementary Table 3). Interestingly, we found that IbDEAH did not find some collinear gene pairs between sweetpotato and Oryza sativa / Brassica rapa / Brassica oleracea, and there was no collinear gene pair between sweetpotato RNA helicase and rice.

Tissue expression patterns of cold-response genes of the RNA helicase family of different sweetpotato cultivars

Objective to preliminarily understand the role of the RNA helicase gene in the development stage of sweetpotato, the expression contours of the RNA helicase genes in diverse tissues of 10 sweetpotato varieties were analyzed by qRT-PCR, such as young leaves, leaves, stems and roots. The expression contours of eight RNA helicase genes were grouped together on their respective heat maps (Fig. 14). The expression of the RNA helicase in various tissues is different, and the expression in different varieties of sweetpotato is also quite different. This also shows that the RNA helicase is related to growth and development [44, 45]. Some the RNA helicase genes were highly communicated in stems and roots of sweetpotato, such as IbDExDH25, IbDEAH42, IbDExDH47, and IbDExDH96. Similarly, some RNA helicase genes were also highly communicated in sweetpotato leaves, including IbDExDH36, IbDExDH48, and IbDEAH53. This shows that the expression of the RNA helicase in unequal growth stages of sweetpotato is diverse and have an influence in the development stage of plants. This is steady with the results of other plants in RNA helicase. HS3 is located in the DEAD-box RNA helicase 22 in Arabidopsis plastids, which is necessary for proper accumulation of plastid gene mRNA throughout seed germination and plant growth [46, 47]. The RNA helicase can also concern the development of plants under chilling stress. AtRH7, one of the RNA helicases in Arabidopsis, is an RNA chaperone involved in cold adaptation [48]. The mutation of rh7 affects the abnormal development of flowers in Arabidopsis thaliana, and makes the leaves of Arabidopsis thaliana smaller in chilling stress [49].

Fig. 14
figure 14

Heat map representation for the organ-specific RNA helicase gene expression profiles in different sweetpotato varieties. YL, Young leaves, L, leaves, S, stems and R, roots. The expression levels of genes are presented using fold-change values transformed to Log2 format. The Log2 (fold-change values) and the color scale are shown at the bottom of heat map. For all the columns labeled with “*,” “**,” and “***” indicate significant differences at P < 0.05, P < 0.01, and P < 0.001, respectively, compared with the YL

Analysis of cold stress related genes in IbDExDH96 overexpressing sweetpotato

From the transcriptome of sweetpotato subjected to cold stress, we screened eight RNA helicase genes. Among them, four genes, namely IbDEAH32, IbDExDH47, IbDExDH48, and IbDExDH96, exhibited responsive behavior to cold stress and were selected for further amplification. Successfully constructing the IbDExDH96 overexpression vector, we mediated its integration into sweetpotato hair roots using rhizobium. Through qRT-PCR analysis, it was observed that the overexpression of OE-IbDExDH96 reached a peak level approximately 15.5 times higher than CK (Fig. 15). Upon analyzing cold-related genes in the transgenic sweetpotato hair roots lines, significant down-regulation was observed in MYB15, ICE1, and CBF2 compared to the CK, while ICE2, CBF1, CBF3, and COR genes showed significant up-regulation in the transgenic sweetpotato hair roots lines (Fig. 16 and Supplementary Table 5). These results indicate that the overexpression of IbDExDH96 in sweetpotato significantly enhances the expression of key genes under cold stress conditions, providing a solid foundation for further exploration of cold stress tolerance in transgenic sweetpotato.

Fig. 15
figure 15

A. rhizogenes-mediated transgenic sweetpotato. A pUBI.U4::IbDExDH96-CaMV35S::DsRed expression diagram, TRs: transgenic roots; ARs: adventitious roots. B Identification of OE-IbDExDH96 transgenic lines, the columns labeled with “***” indicate significant differences at P < 0.001

Fig. 16
figure 16

Expression of cold stress-related genes in OE-IbDExDH96 transgenic lines, the columns labeled with “***” indicate significant differences at P < 0.001

Discussion

The RNA helicases are ubiquitous in miscellaneous organisms [50]. It is involved in almost all processes of RNA metabolism, such as transcription, mRNA splicing and output, mRNA translation, etc. It involves almost all aspects of organisms [4, 51, 52]. In the past, the study of the RNA helicase family in plants mainly focuses on dicotyledonous model plant Arabidopsis thaliana and monocotyledonous rice [53,54,55], these play an essential role in plant growth and stress response [56, 57]. Sweetpotato is a significant food crop, which is broadly used in food, feed and industrial raw materials. It is a hexaploid plant, the seventh largest food crop in the world, with the strong ability to adapt to the environment, varieties of high yield, strong stress resistance [58, 59]. However, only 17 RNA helicase genes have been nominated in Ipomoea trifida [20]. There is no comprehensive study on RNA helicase in sweetpotato. This study systematically identified RNA helicase genes in sweetpotato, which laid a solid basis for further study on the engagement of the RNA helicase in plant abiotic stress and development.

In the present study, we conducted a comprehensive reasoning of the sweetpotato RNA helicase gene family, containing phylogenetic tree, expression profiles of different sweetpotato varieties under common growth circumstances, and expression profiles under diverse abiotic stresses and hormone stresses. 300 RNA helicase genes were nominated in sweetpotato genome, which is a very large gene family. A larger family of the RNA helicase gene exists in sweetpotato, suggesting that RNA helicase plays a leading role in regulating environmental responses. According to the family classification of Arabidopsis and rice RNA helicases [21], they were divided into three subfamilies, including IbDEAD (53), IbDEAH (54) and IbDExDH (193) (Supplementary File 1). Studies have indicated that the DExDH subfamily in Arabidopsis, rice, and maize is larger than the DEAD and DEAH subfamilies. Additionally, being a hexaploid plant, sweetpotato possesses a greater number of genes in each subfamily compared to diploid and tetraploid plants.

The sweetpotato RNA helicase gene contains different numbers of exons and different lengths of introns (Figs. 2, 3 and 4). Intriguingly, five genes in the IbDExDH subfamily contain only one exon, while the IbDEAD and IbDEAH subfamilies do not contain this case and contain multiple exons. In fact, the length of RNA helicase family proteins varies greatly. IbDEAD is 323 to 1301 amino acids, IbDEAH is 206 to 2904 amino acids, and IbDExDH is 128 to 2801 amino acids. The highly different amino acid sequences of these sweetpotato RNA helicases induce diverse protein structures and functions in dissimilar developmental and stress-resistant environments [60].

At the same time, the sweetpotato RNA helicase motif was analyzed, the most conserved motif in different species is Asp-Glu-Ala-Asp, which is divided into three subfamilies according to its difference [16, 61, 62]. According to the structural characteristics and phylogenetic analysis of the Motif V region, the determined helicases can be moreover divided into these subfamilies, comprising IbDEAD, IbDEAH, and IbDExDH. Phylogenetic analysis showed that IbDEAD, IbDEAH and IbDExDH RNA helicase proteins can be moreover divided into nine, or thirteen great subgroups (Fig. 1). The gene structure results showed that the main RNA helicase genes in Arabidopsis were uniform to the AtRH family genes [24], but the position and length of introns were not fixed. However, this is the first genome-wide examination of the RNA helicase genes family in sweetpotato. The distinct subfamilies and gene structures of sweetpotato RNA helicase genes indirectly indicate the different functions in RNA metabolism, stress resistance and growth and development [49, 63].

In fact, we establish that most RNA helicase gene promoters contain some cis-regulatory elements, such as plant development [64], abiotic stress [65], plant hormones, and light response elements [50, 53, 66]. It is worth noting that IbDEAD and IbDEAH have more cis-regulatory elements associated to low temperature and light response, and IbDExDH has more cis-acting elements associated with abscisic acid and plant hormones (Figs. 6, 7 and 8). According to the study that temperature, abscisic acid and jasmonic acid are involved in abiotic stress processes in plants [67,68,69]. Under hormone and abiotic stress treatments, qRT-PCR data showed that RNA helicase mainly responded to abiotic stress (Figs. 9 and 10).

It was observed that under drought stress, the DEAD-box RNA helicase CmRH56 inhibited rhizome growth by disrupting gibberellic acid (GA) biosynthesis [70]. By regulating salt-tolerant genes, the DEAD-box RNA helicase BnRH6 plays a negative role in promoting the adaptation of Brassica plants to salt stress [71]. By modulating the extent of membrane lipid peroxidation, Triticum aestivum DEAD-box RNA helicase TaDEAD-box57 has the potential to enhance the tolerance of transgenic plants to drought, salt, and low-temperature stresses [72]. In Arabidopsis, DEAD-box LOS4 is able to participate in the process of cryogenic stress processes [73, 74]. Among them, in the cold stress response pathway, LOS RNA helicase plays a key role in target gene output, maturation and reaction to temperature stress. The transcription of STRS1 and STRS2 was inhibited under salt stress. The salt tolerance of mutants strs1 and strs2 was enhanced, and the expression of RD29A, DREB1A and DREB2A was enhanced [75]. The development of the germ and leaf of the Arabidopsis rh7 mutant was seriously delayed under low temperature stress [49]. The growth of rh3 mutants was severely inhibited under salt or cold stress [76]. Our previous findings showed that the tolerance of tomato SlDEAD31 was enhanced in salt and drought stress [41]. The rice SUV3 protein has DNA and RNA helicase and ATPase activities, and SUV3 expression can be induced by salt stress [77, 78]. Low temperature and high salt stress can induce the expression of GmRH in soybeans, and GmRH plays a significant in RNA processing [79]. Tobacco P68 can enhance tolerance to salt stress [80]. AvDH1 increased salt tolerance [81]. The Arabidopsis RCF1 gene plays an integral role in maintaining normal splicing of mRNA precursors, and some cold stress-induced genes were error spliced in the rcf-1 mutant [82]. Maize DRH1 can interact with the nucleoprotein fiber MA16, which is involved in ribosomal RNA metabolism [83]. DEVH-box RNA helicase AtHELPS play a key role in K+ deprivation in Arabidopsis thaliana [84]. The presence of numerous cis-regulatory elements in the promoter regions of sweetpotato RNA helicase genes highlights their crucial role in conferring stress resistance to the plant. This discovery has led us to identify novel candidate genes involved in regulating sweetpotato ability to withstand abiotic stress through RNA metabolism. Given that sweetpotato ranks as the seventh largest food crop globally, its economic significance and environmental adaptability pose significant challenges. The unveiling of stress-related genes has laid a strong groundwork for advancing our understanding of the molecular mechanisms underlying sweetpotato resistance and expediting the development of resilient sweetpotato varieties through breeding efforts.

RNA helicases play a crucial role in regulating plant development and responding to environmental cues. To gain insights into their functionality, we conducted an analysis of the expression profiles of RNA helicase genes in various tissues of different sweetpotato varieties. Additionally, their response patterns were examined under four abiotic stresses and four hormone stresses (Fig. 14). Our findings revealed significant variations in the expression of sweetpotato RNA helicase genes across different tissues, indicating their involvement in plant development. Furthermore, these genes exhibited prominent responsiveness to abiotic stresses, such as cold, heat, and salt, underscoring their significance in enhancing plant stress resistance.This is consistent with the functional studies of other plant RNA helicases, for instance Arabidopsis [48], rice [85], tomato [86], rapeseed [71], chrysanthemum [70], Zea mays [87], barley [72] and soybean [88]. Furthermore, we proceeded to develop sweetpotato lines through A.rhizogenes-mediated overexpression of IbDExDH96, an RNA helicase gene. This allowed us to conduct initial validation of its involvement in cold stress response in sweetpotato. Based on our findings, we propose that RNA helicase proteins not only serve as functional genes influencing the growth and development of sweetpotato but also act as regulatory factors under diverse environmental conditions. Hence, our study presents novel experimental insights and provides valuable leads to support the aforementioned speculation.

In summary, delving into the reception of the RNA helicase gene in sweetpotato holds the potential to enhance transgenic research aimed at improving the crop's yield and stress resistance. The exploration of the identification, classification, and phylogenetic tree construction of sweetpotato RNA helicase genes through bioinformatics has yielded valuable insights for a comprehensive investigation into the biological functions of these genes. Furthermore, these studies significantly contribute to our understanding of the molecular basis underlying several crucial agricultural traits in sweetpotato cultivation. However, the precise regulatory mechanism governing the development and stress response of sweetpotato RNA helicase genes remains elusive and warrants further investigation.

Conclusion

In this research, we conducted an exhaustive genome-wide reasoning of the sweetpotato RNA helicase family, containing chromosome distribution, promoter elements, and protein motif analysis. All of 300 RNA helicase genes were detected, containing IbDEAD, IbDEAH and IbDExDH subfamilies. The expression patterns of eight RNA helicase genes in different sweetpotato varieties and their responses to abiotic stress and hormonal stress were analyzed by qRT-PCR. The expression of the RNA helicase genes was significantly distinct in individual tissues of 10 sweetpotato varieties and notably raised under divergent abiotic stresses. The overexpression of the RNA helicase gene IbDExDH96 in transgenic sweetpotato hair roots, mediated by rhizobium, led to alterations in the expression levels of cold-related genes. This enabled us to investigate the role of sweetpotato RNA helicase in alleviating cold stress. We obtained preliminary evidence supporting its involvement in endogenous cold stress tolerance. However, further exploration is still required to elucidate the specific mechanism of the sweetpotato RNA helicase gene. The results showed that RNA helicase was complicated in the direction of extension and the resistance to stress of sweetpotato. This study supplies new inspirations into the development and exploration of the RNA helicase gene families.

Availability of data and materials

The RNA-seq data in this study was acquired from NCBI with the accession number GSE255226 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE255226). The datasets supporting the conclusions of this article are included in the article and its supplementary material.

References

  1. Heaton SM, Gorry PR, Borg NA. DExD/H-box helicases in HIV-1 replication and their inhibition. Trends Microbiol. 2023;31(4):393–404.

    Article  CAS  PubMed  Google Scholar 

  2. Tapescu I, Taschuk F, Pokharel SM, Zginnyk O, Ferretti M, Bailer PF, Whig K, Madden EA, Heise MT, Schultz DC, et al. The RNA helicase DDX39A binds a conserved structure in Chikungunya virus RNA to control infection. Mol Cell. 2023;83(22):4174-4189e4177.

    Article  CAS  PubMed  Google Scholar 

  3. You L, Shi C, Wang D, Fu ZQ. Helicases clear hurdles during plant defense protein translation. Trends Biochem Sci. 2023;S0968-0004(23):00276–1.

  4. Vashisht AA, Tuteja N. Stress responsive DEAD-box helicases: a new pathway to engineer plant stress tolerance. J Photochem Photobiol B. 2006;84(2):150–60.

    Article  CAS  PubMed  Google Scholar 

  5. Wang X, Kong W, Wang Y, Wang J, Zhong L, Lao K, Dong X, Zhang D, Huang H, Mo B, et al. Uridylation and the SKI complex orchestrate the calvin cycle of photosynthesis through RNA surveillance of TKL1 in Arabidopsis. Proc Natl Acad Sci U S A. 2022;119(38):e2205842119.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Gc K, Gyawali P, Balci H, Abeysirigunawardena S. Ribosomal RNA methyltransferase RsmC moonlights as an RNA chaperone. ChemBioChem. 2020;21(13):1885–92.

    Article  PubMed  Google Scholar 

  7. Nawaz G, Kang H. Chloroplast- or Mitochondria-targeted DEAD-Box RNA helicases play essential roles in Organellar RNA metabolism and abiotic stress responses. Front Plant Sci. 2017;8:871.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Tao GY, Xie YH, Li WF, Li KP, Sun C, Wang HM, Sun XM. LkARF7 and LkARF19 overexpression promote adventitious root formation in a heterologous poplar model by positively regulating LkBBM1. Commun Biol. 2023;6(1):372.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Xiang Y, Huang W, Tan L, Chen T, He Y, Irving PS, Weeks KM, Zhang QC, Dong X. Pervasive downstream RNA hairpins dynamically dictate start-codon selection. Nature. 2023;621(7978):423–30.

    Article  CAS  PubMed  PubMed Central  ADS  Google Scholar 

  10. de la Cruz J, Kressler D, Linder P. Unwinding RNA in Saccharomyces cerevisiae: DEAD-box proteins and related families. Trends Biochem Sci. 1999;24(5):192–8.

    Article  PubMed  Google Scholar 

  11. Lorsch JR. RNA chaperones exist and DEAD box proteins get a life. Cell. 2002;109(7):797–800.

    Article  CAS  PubMed  Google Scholar 

  12. Sahoo RK, Chandan RK, Swain DM, Tuteja N, Jha G. Heterologous overexpression of PDH45 gene of pea provides tolerance against sheath blight disease and drought stress in rice. Plant Physiol Biochem. 2022;186:242–51.

    Article  CAS  PubMed  Google Scholar 

  13. Rocak S, Linder P. DEAD-box proteins: the driving forces behind RNA metabolism. Nat Rev Mol Cell Biol. 2004;5(3):232–41.

    Article  CAS  PubMed  Google Scholar 

  14. Xu R, Zhang S, Lu L, Cao H, Zheng C. A genome-wide analysis of the RNA helicase gene family in Solanum lycopersicum. Gene. 2013;513(1):128–40.

    Article  CAS  PubMed  Google Scholar 

  15. Tuteja N, Tuteja R. Unraveling DNA helicases. Motif, structure, mechanism and function. Eur J Biochem. 2004;271(10):1849–63.

    Article  CAS  PubMed  Google Scholar 

  16. Tanner NK, Linder P. DExD/H box RNA helicases: from generic motors to specific dissociation functions. Mol Cell. 2001;8(2):251–62.

    Article  CAS  PubMed  Google Scholar 

  17. Jiang Z, Zhang K, Li Z, Li Z, Yang M, Jin X, Cao Q, Wang X, Yue N, Li D, et al. The barley stripe mosaic virus γb protein promotes viral cell-to-cell movement by enhancing ATPase-mediated assembly of ribonucleoprotein movement complexes. PLoS Pathog. 2020;16(7):e1008709.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Boudet N, Aubourg S, Toffano-Nioche C, Kreis M, Lecharny A. Evolution of intron/exon structure of DEAD helicase family genes in arabidopsis, caenorhabditis, and drosophila. Genome Res. 2001;11(12):2101–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Umate P, Tuteja R, Tuteja N. Genome-wide analysis of helicase gene family from rice and Arabidopsis: a comparison with yeast and human. Plant Mol Biol. 2010;73(4–5):449–65.

    Article  CAS  PubMed  Google Scholar 

  20. Wan R, Liu J, Yang Z, Zhu P, Cao Q, Xu T. Genome-wide identification, characterisation and expression profile analysis of DEAD-box family genes in sweet potato wild ancestor Ipomoea trifida under abiotic stresses. Genes Genomics. 2020;42(3):325–35.

    Article  CAS  PubMed  Google Scholar 

  21. Xu R, Zhang S, Huang J, Zheng C. Genome-wide comparative in silico analysis of the RNA helicase gene family in Zea mays and Glycine max: a comparison with Arabidopsis and Oryza sativa. PLoS One. 2013;8(11):e78982.

    Article  CAS  PubMed  PubMed Central  ADS  Google Scholar 

  22. Chen J, Zhang Y, Liu J, Xia M, Wang W, Shen F. Genome-wide analysis of the RNA helicase gene family in Gossypium Raimondii. Int J Mol Sci. 2014;15(3):4635–56.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Karthik S, Tuteja N, Ganapathi A, Manickavasagam M. Pea p68, a DEAD-box helicase, enhances salt tolerance in marker-free transgenic plants of soybean [Glycine max (L.) Merrill]. 3 Biotech. 2019;9(1):10.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Aubourg S, Kreis M, Lecharny A. The DEAD box RNA helicase family in Arabidopsis thaliana. Nucleic Acids Res. 1999;27(2):628–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Katayama K, Kobayashi A, Sakai T, Kuranouchi T, Kai Y. Recent progress in sweetpotato breeding and cultivars for diverse applications in Japan. Breed Sci. 2017;67(1):3–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Isobe S, Shirasawa K, Hirakawa H. Challenges to genome sequence dissection in sweetpotato. Breed Sci. 2017;67(1):35–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Yan M, Nie H, Wang Y, Wang X, Jarret R, Zhao J, Wang H, Yang J. Exploring and exploiting genetics and genomics for sweetpotato improvement: Status and perspectives. Plant Commun. 2022;3(5):100332.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Ramamoorthy P, Bheemanahalli R, Meyers SL, Shankle MW, Reddy KR. Drought, Low Nitrogen Stress, and Ultraviolet-B Radiation effects on Growth, Development, and physiology of Sweetpotato cultivars during early season. Genes (Basel). 2022;13(1):156.

  29. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.

    Article  CAS  PubMed  Google Scholar 

  30. Xu R, Zhang S, Huang J, Zheng C. Correction: genome-wide comparative in silico analysis of the RNA helicase Gene Family in Zea mays and Glycine max: a comparison with Arabidopsis and Oryza sativa. PLoS One. 2020;15(10):e0240759.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Yang J, Moeinzadeh MH, Kuhl H, Helmuth J, Xiao P, Haas S, Liu G, Zheng J, Sun Z, Fan W, et al. Haplotype-resolved sweet potato genome traces back its hexaploidization history. Nat Plants. 2017;3(9):696–703.

    Article  CAS  PubMed  Google Scholar 

  32. Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Guindon S, Lethiec F, Duroux P, Gascuel O. PHYML Online–a web server for fast maximum likelihood-based phylogenetic inference. Nucleic Acids Res. 2005;33(Web Server issue):W557-559.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Chen C, Chen H, Zhang Y, Thomas HR, Frank MH, He Y, Xia R. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020;13(8):1194–202.

    Article  CAS  PubMed  Google Scholar 

  35. Bailey TL, Johnson J, Grant CE, Noble WS. The MEME suite. Nucleic Acids Res. 2015;43(W1):W39–49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Gasteiger E, Gattiker A, Hoogland C, Ivanyi I, Appel RD, Bairoch A. ExPASy: the proteomics server for in-depth protein knowledge and analysis. Nucleic Acids Res. 2003;31(13):3784–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Wang Y, Tang H, Debarry JD, Tan X, Li J, Wang X, Lee TH, Jin H, Marler B, Guo H, et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012;40(7):e49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA. Circos: an information aesthetic for comparative genomics. Genome Res. 2009;19(9):1639–45.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Yu Y, Xuan Y, Bian X, Zhang L, Pan Z, Kou M, Cao Q, Tang Z, Li Q, Ma D, et al. Overexpression of phosphatidylserine synthase IbPSS1 affords cellular na(+) homeostasis and salt tolerance by activating plasma membrane na(+)/H(+) antiport activity in sweet potato roots. Hortic Res. 2020;7:131.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Xie Z, Zhou Z, Li H, Yu J, Jiang J, Tang Z, Ma D, Zhang B, Han Y, Li Z. High throughput sequencing identifies chilling responsive genes in sweetpotato (Ipomoea batatas Lam.) During storage. Genomics. 2019;111(5):1006–17.

    Article  CAS  PubMed  Google Scholar 

  41. Zhu M, Chen G, Dong T, Wang L, Zhang J, Zhao Z, Hu Z. SlDEAD31, a putative DEAD-Box RNA helicase gene, regulates salt and drought tolerance and stress-related genes in tomato. PLoS ONE. 2015;10(8):e0133849.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Camborde L, Kiselev A, Pel MJC, Le Ru A, Jauneau A, Pouzet C, Dumas B, Gaulin E. An oomycete effector targets a plant RNA helicase involved in root development and defense. New Phytol. 2022;233(5):2232–48.

    Article  CAS  PubMed  Google Scholar 

  43. Li Y, Du Y, Huai J, Jing Y, Lin R. The RNA helicase UAP56 and the E3 ubiquitin ligase COP1 coordinately regulate alternative splicing to repress photomorphogenesis in Arabidopsis. Plant Cell. 2022;34(11):4191–212.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Tyagi V, Parihar V, Malik G, Kalra V, Kapoor S, Kapoor M. The DEAD-box RNA helicase eIF4A regulates plant development and interacts with the hnRNP LIF2L1 in Physcomitrella patens. Mol Genet Genomics. 2020;295(2):373–89.

    Article  CAS  PubMed  Google Scholar 

  45. Liu L, Zhou X, Cheng S, Ge Y, Chen B, Shi J, Li H, Li S, Li Y, Yuan J, et al. RNA-binding protein DHX9 promotes glioma growth and tumor-associated macrophages infiltration via TCF12. CNS Neurosci Ther. 2023;29(4):988–99.

    Article  CAS  PubMed  Google Scholar 

  46. Kanai M, Hayashi M, Kondo M, Nishimura M. The plastidic DEAD-box RNA helicase 22, HS3, is essential for plastid functions both in seed development and in seedling growth. Plant Cell Physiol. 2013;54(9):1431–40.

    Article  CAS  PubMed  Google Scholar 

  47. Iglesias-Fernández R, Pastor-Mora E, Vicente-Carbajosa J, Carbonero P. A possible role of the Aleurone expressed gene HvMAN1 in the hydrolysis of the cell wall mannans of the starchy endosperm in germinating hordeum vulgare L. seeds. Front Plant Sci. 2019;10:1706.

    Article  PubMed  Google Scholar 

  48. Huang CK, Shen YL, Huang LF, Wu SJ, Yeh CH, Lu CA. The DEAD-Box RNA helicase AtRH7/PRH75 participates in Pre-rRNA processing, plant development and cold tolerance in arabidopsis. Plant Cell Physiol. 2016;57(1):174–91.

    Article  CAS  PubMed  Google Scholar 

  49. Liu Y, Tabata D, Imai R. A Cold-Inducible DEAD-Box RNA helicase from arabidopsis thaliana regulates plant growth and development under low temperature. PLoS ONE. 2016;11(4):e0154040.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Li X, Li C, Zhu J, Zhong S, Zhu H, Zhang X. Functions and mechanisms of RNA helicases in plants. J Exp Bot. 2023;74(7):2295–310.

    Article  CAS  PubMed  Google Scholar 

  51. Donsbach P, Klostermeier D. Regulation of RNA helicase activity: principles and examples. Biol Chem. 2021;402(5):529–59.

    Article  CAS  PubMed  Google Scholar 

  52. Sloan KE, Bohnsack MT. Unravelling the mechanisms of RNA helicase regulation. Trends Biochem Sci. 2018;43(4):237–50.

    Article  CAS  PubMed  Google Scholar 

  53. Nawaz G, Kang H. Rice OsRH58, a chloroplast DEAD-box RNA helicase, improves salt or drought stress tolerance in Arabidopsis by affecting chloroplast translation. BMC Plant Biol. 2019;19(1):17.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Takagi M, Iwamoto N, Kubo Y, Morimoto T, Takagi H, Takahashi F, Nishiuchi T, Tanaka K, Taji T, Kaminaka H, et al. Arabidopsis SMN2/HEN2, Encoding DEAD-Box RNA helicase, governs proper expression of the Resistance Gene SMN1/RPS6 and is involved in dwarf, autoimmune phenotypes of mekk1 and mpk4 mutants. Plant Cell Physiol. 2020;61(8):1507–16.

    Article  CAS  PubMed  Google Scholar 

  55. Banu MSA, Huda KMK, Harun-Ur-Rashid M, Parveen S, Tuteja N. A DEAD box helicase Psp68 positively regulates salt stress responses in marker-free transgenic rice plants. Transgenic Res. 2023;32(4):293–304.

    Article  CAS  PubMed  Google Scholar 

  56. Kim JS, Kim KA, Oh TR, Park CM, Kang H. Functional characterization of DEAD-box RNA helicases in Arabidopsis thaliana under abiotic stress conditions. Plant Cell Physiol. 2008;49(10):1563–71.

    Article  CAS  PubMed  Google Scholar 

  57. Li X, Li C, Zhu J, Zhong S, Zhu H, Zhang X. Functions and mechanisms of RNA helicases in plants. J Exp Bot. 2022.

  58. Yu J, Su D, Yang D, Dong T, Tang Z, Li H, Han Y, Li Z, Zhang B. Chilling and heat stress-Induced physiological changes and MicroRNA-Related mechanism in Sweetpotato (Ipomoea batatas L). Front Plant Sci. 2020;11:687.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Banu MSA, Huda KMK, Harun-Ur-Rashid M, Parveen S, Shahinul Islam SM, Tuteja N. Phenotypic and microarray analysis reveals salinity stress-induced oxidative tolerance in transgenic rice expressing a DEAD-box RNA helicase, OsDB10. Plant Mol Biol. 2023;113(1–3):19–32.

    Article  CAS  PubMed  Google Scholar 

  60. Pandey S, Prasad A, Sharma N, Prasad M. Linking the plant stress responses with RNA helicases. Plant Sci. 2020;299:110607.

    Article  CAS  PubMed  Google Scholar 

  61. Linder P, Jankowsky E. From unwinding to clamping - the DEAD box RNA helicase family. Nat Rev Mol Cell Biol. 2011;12(8):505–16.

    Article  CAS  PubMed  Google Scholar 

  62. Xu C, Zhang Z, He J, Bai Y, Cui J, Liu L, Tang J, Tang G, Chen X, Mo B. The DEAD-box helicase RCF1 plays roles in miRNA biogenesis and RNA splicing in Arabidopsis. Plant J. 2023;116(1):144–60.

    Article  CAS  PubMed  Google Scholar 

  63. Zhang XD, Han Y, Yang ZM, Sun D. DEAD-box RNA helicase 6 regulates drought and abscisic acid stress responses in rapeseed (Brassica napus). Gene. 2023;886:147717.

    Article  CAS  PubMed  Google Scholar 

  64. Luo P, Shi C, Zhou Y, Zhou J, Zhang X, Wang Y, Yang Y, Peng X, Xie T, Tang X. The nuclear-localized RNA helicase 13 is essential for chloroplast development in Arabidopsis thaliana. J Exp Bot. 2023;74(17):5057–71.

    Article  CAS  PubMed  Google Scholar 

  65. Zheng M, Song Y, Wang L, Yang D, Yan J, Sun Y, Hsu YF. CaRH57, a RNA helicase, contributes pepper tolerance to heat stress. Plant Physiol Biochem. 2023;205:108202.

    Article  CAS  PubMed  Google Scholar 

  66. Huang C-K, Sie Y-S, Chen Y-F, Huang T-S, Lu C-A. Two highly similar DEAD box proteins, OsRH2 and OsRH34, homologous to eukaryotic initiation factor 4AIII, play roles of the exon junction complex in regulating growth and development in rice. BMC Plant Biol. 2016;16(1):84.

    Article  PubMed  PubMed Central  Google Scholar 

  67. Bari R, Jones JD. Role of plant hormones in plant defence responses. Plant Mol Biol. 2009;69(4):473–88.

    Article  CAS  PubMed  Google Scholar 

  68. Bharath P, Gahir S, Raghavendra AS. Abscisic acid-induced stomatal closure: an important component of plant defense against abiotic and biotic stress. Front Plant Sci. 2021;12:615114.

    Article  PubMed  PubMed Central  Google Scholar 

  69. Nidumukkala S, Tayi L, Chittela RK, Vudem DR, Khareedu VR. DEAD box helicases as promising molecular tools for engineering abiotic stress tolerance in plants. Crit Rev Biotechnol. 2019;39(3):395–407.

    Article  CAS  PubMed  Google Scholar 

  70. Zhang L, Xu Y, Liu X, Qin M, Li S, Jiang T, Yang Y, Jiang CZ, Gao J, Hong B, et al. The chrysanthemum DEAD-box RNA helicase CmRH56 regulates rhizome outgrowth in response to drought stress. J Exp Bot. 2022;73(16):5671–81.

    Article  PubMed  Google Scholar 

  71. Zhang X, Song J, Wang L, Yang ZM, Sun D. Identification of a DEAD-box RNA helicase BnRH6 reveals its involvement in salt stress response in rapeseed (Brassica napus). Int J Mol Sci. 2022;24(1):2.

  72. Ru JN, Hou ZH, Zheng L, Zhao Q, Wang FZ, Chen J, Zhou YB, Chen M, Ma YZ, Xi YJ, et al. Genome-wide analysis of DEAD-box RNA helicase family in wheat (Triticum aestivum) and functional identification of TaDEAD-box57 in Abiotic stress responses. Front Plant Sci. 2021;12:797276.

    Article  PubMed  PubMed Central  Google Scholar 

  73. Gong Z, Dong CH, Lee H, Zhu J, Xiong L, Gong D, Stevenson B, Zhu JK. A DEAD box RNA helicase is essential for mRNA export and important for development and stress responses in Arabidopsis. Plant Cell. 2005;17(1):256–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Gong Z, Lee H, Xiong L, Jagendorf A, Stevenson B, Zhu JK. RNA helicase-like protein as an early regulator of transcription factors for plant chilling and freezing tolerance. Proc Natl Acad Sci U S A. 2002;99(17):11507–12.

    Article  CAS  PubMed  PubMed Central  ADS  Google Scholar 

  75. Kant P, Kant S, Gordon M, Shaked R, Barak S. STRESS RESPONSE SUPPRESSOR1 and STRESS RESPONSE SUPPRESSOR2, two DEAD-box RNA helicases that attenuate Arabidopsis responses to multiple abiotic stresses. Plant Physiol. 2007;145(3):814–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  76. Gu L, Xu T, Lee K, Lee KH, Kang H. A chloroplast-localized DEAD-box RNA helicaseAtRH3 is essential for intron splicing and plays an important role in the growth and stress response in Arabidopsis thaliana. Plant Physiol Biochem. 2014;82:309–18.

    Article  CAS  PubMed  Google Scholar 

  77. Tuteja N, Sahoo RK, Garg B, Tuteja R. OsSUV3 dual helicase functions in salinity stress tolerance by maintaining photosynthesis and antioxidant machinery in rice (Oryza sativa L. Cv. IR64). Plant J. 2013;76(1):115–27.

    Article  CAS  PubMed  Google Scholar 

  78. Sahoo RK, Ansari MW, Tuteja R, Tuteja N. Salt tolerant SUV3 overexpressing transgenic rice plants conserve physicochemical properties and microbial communities of rhizosphere. Chemosphere. 2015;119:1040–7.

    Article  CAS  PubMed  ADS  Google Scholar 

  79. Chung E, Cho CW, Yun BH, Choi HK, So HA, Lee SW, Lee JH. Molecular cloning and characterization of the soybean DEAD-box RNA helicase gene induced by low temperature and high salinity stress. Gene. 2009;443(1–2):91–9.

    Article  CAS  PubMed  Google Scholar 

  80. Tuteja N, Banu MS, Huda KM, Gill SS, Jain P, Pham XH, Tuteja R. Pea p68, a DEAD-box helicase, provides salinity stress tolerance in transgenic tobacco by reducing oxidative stress and improving photosynthesis machinery. PLoS One. 2014;9(5):e98287.

    Article  PubMed  PubMed Central  ADS  Google Scholar 

  81. Chen J, Wan S, Liu H, Fan S, Zhang Y, Wang W, Xia M, Yuan R, Deng F, Shen F. Overexpression of an Apocynum Venetum DEAD-Box helicase gene (AvDH1) in cotton confers salinity tolerance and increases yield in a saline field. Front Plant Sci. 2015;6:1227.

    PubMed  Google Scholar 

  82. Guan Q, Wu J, Zhang Y, Jiang C, Liu R, Chai C, Zhu J. A DEAD box RNA helicase is critical for pre-mRNA splicing, cold-responsive gene regulation, and cold tolerance in Arabidopsis. Plant Cell. 2013;25(1):342–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  83. Gendra E, Moreno A, Albà MM, Pages M. Interaction of the plant glycine-rich RNA-binding protein MA16 with a novel nucleolar DEAD box RNA helicase protein from Zea mays. Plant J. 2004;38(6):875–86.

    Article  CAS  PubMed  Google Scholar 

  84. Xu RR, Qi SD, Lu LT, Chen CT, Wu CA, Zheng CC. A DExD/H box RNA helicase is important for K + deprivation responses and tolerance in Arabidopsis thaliana. Febs j. 2011;278(13):2296–306.

    Article  CAS  PubMed  Google Scholar 

  85. Xiaomei W, Rongrong K, Ting Z, Yuanyuan G, Jianlong X, Zhongze P, Gangseob L, Dongzhi L, Yanjun D. A DEAD-box RNA helicase TCD33 that confers chloroplast development in rice at seedling stage under cold stress. J Plant Physiol. 2020;248:153138.

    Article  PubMed  Google Scholar 

  86. Capel C, Albaladejo I, Egea I, Massaretto IL, Yuste-Lisbona FJ, Pineda B, García-Sogo B, Angosto T, Flores FB, Moreno V, et al. The res (restored cell structure by salinity) tomato mutant reveals the role of the DEAD-box RNA helicase SlDEAD39 in plant development and salt response. Plant Cell Environ. 2020;43(7):1722–39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  87. Yang YZ, Ding S, Liu XY, Xu C, Sun F, Tan BC. The DEAD-box RNA helicase ZmRH48 is required for the splicing of multiple mitochondrial introns, mitochondrial complex biosynthesis, and seed development in maize. J Integr Plant Biol. 2023;65(11):2456–68.

    Article  CAS  PubMed  Google Scholar 

  88. Wang Y, Liao J, Wu J, Huang H, Yuan Z, Yang W, Wu X, Li X. Genome-wide identification and characterization of the soybean DEAD-Box Gene Family and expression response to Rhizobia. Int J Mol Sci. 2022;23(3):1120.

Download references

Funding

This work was supported by the earmarked fund for CARS-10-Sweetpotato, Postgraduate Research & Practice Innovation Program of Jiangsu Province (KYCX22_2818).

Author information

Authors and Affiliations

Authors

Contributions

ZL, MZ and TD conceived and designed the research. FM, HZ and QZ performed the research and analyzed the data. FM wrote the manuscript. MZ, LK and ZL helped to revise the manuscript. All authors read and approved the manuscript.

Corresponding author

Correspondence to Zongyun Li.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interest

The authors declare 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

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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Mu, F., Zheng, H., Zhao, Q. et al. Genome-wide systematic survey and analysis of the RNA helicase gene family and their response to abiotic stress in sweetpotato. BMC Plant Biol 24, 193 (2024). https://doi.org/10.1186/s12870-024-04824-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-024-04824-z

Keywords