Next Article in Journal
Influence of Green, Red and Blue Light Emitting Diodes on Multiprotein Complex Proteins and Photosynthetic Activity under Different Light Intensities in Lettuce Leaves (Lactuca sativa L.)
Previous Article in Journal
Testosterone Reduces Knee Passive Range of Motion and Expression of Relaxin Receptor Isoforms via 5α-Dihydrotestosterone and Androgen Receptor Binding
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 15
Issue 3
10.3390/ijms15034635
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Analysis of the RNA Helicase Gene Family in Gossypium raimondii

by
Jie Chen
,
Yujuan Zhang
,
Jubo Liu
,
Minxuan Xia
,
Wei Wang
and
Fafu Shen
*
State Key Laboratory of Crop Biology, College of Agronomy, Shandong Agricultural University, Tai'an 271018, Shandong, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2014, 15(3), 4635-4656; https://doi.org/10.3390/ijms15034635
Submission received: 12 January 2014 / Revised: 25 February 2014 / Accepted: 27 February 2014 / Published: 17 March 2014
(This article belongs to the Section Biochemistry)
Browse Figures
Versions Notes

Abstract

:
The RNA helicases, which help to unwind stable RNA duplexes, and have important roles in RNA metabolism, belong to a class of motor proteins that play important roles in plant development and responses to stress. Although this family of genes has been the subject of systematic investigation in Arabidopsis, rice, and tomato, it has not yet been characterized in cotton. In this study, we identified 161 putative RNA helicase genes in the genome of the diploid cotton species Gossypium raimondii. We classified these genes into three subfamilies, based on the presence of either a DEAD-box (51 genes), DEAH-box (52 genes), or DExD/H-box (58 genes) in their coding regions. Chromosome location analysis showed that the genes that encode RNA helicases are distributed across all 13 chromosomes of G. raimondii. Syntenic analysis revealed that 62 of the 161 G. raimondii helicase genes (38.5%) are within the identified syntenic blocks. Sixty-six (40.99%) helicase genes from G. raimondii have one or several putative orthologs in tomato. Additionally, GrDEADs have more conserved gene structures and more simple domains than GrDEAHs and GrDExD/Hs. Transcriptome sequencing data demonstrated that many of these helicases, especially GrDEADs, are highly expressed at the fiber initiation stage and in mature leaves. To our knowledge, this is the first report of a genome-wide analysis of the RNA helicase gene family in cotton.

1. Introduction

RNA helicases are found in many organisms and have important roles in RNA metabolism. These enzymes participate in many metabolic pathways that involve the separation of double-stranded nucleic acid into single strands and the removal of nucleic-acid-associated proteins [1]. Usually, they translocate along the bound strand unidirectionally and their specific polarity depends upon the bound strand. RNA helicases potentially regulate cellular growth and differentiation, as well as responses to abiotic stress, by affecting nuclear mRNA export, translation initiation, mRNA decay, rRNA processing, cell cycle progression, and the initiation of transcription [2–4].
Based on variations within the DEAD (Asp-Glu-Ala-Asp) motif, the RNA helicase superfamily proteins can be classified into three subfamilies, which are defined by the presence of either the DEAD-box, DEAH-box, or DExD/H-box [5,6]. Multiple members of each subfamily have been found in several plants. Many of them, especially DEAD-box genes, are key regulators of developmental processes and responses to diverse abiotic stresses, such as salt stress, oxygen levels, light, or temperature [3]. In Arabidopsis, the DEAD-box RNA helicase LOS4 controls responses to low temperature and processes such as flowering and vernalization [7,8]. STRS1 and STRS2 mediate stress responses to various abiotic stresses [9]. RCF1 maintains proper splicing of pre-mRNAs and regulates responses to changes in temperature [10]. Tobacco VDL (for variegated and distorted leaf), a plastid DEAD-box RNA helicase, is essential for chloroplast differentiation and plant morphogenesis [11]. Maize DRH1 interacts with the nucleolar protein fibrillarin, MA16, and controls the metabolism of ribosomal RNA [12]. Rice ABP (ATP-binding protein) is upregulated in response to multiple abiotic stresses [13]. We reported that the DEAD-box RNA helicase AvDH1 (Apocynum venetum ATP-dependent DEAD-box helicase) participates in the regulation of salt tolerance in the halophyte Apocynum venetum [14]. The nucleolar DExD/H box RNA helicase AtMTR4 (mRNA transport protein) is required for normal rRNA biogenesis and development in Arabidopsis [15]. ISE2, a DEVH-box RNA helicase, has been shown to be involved in plasmodesmata function during embryogenesis in Arabidopsis [16]. The DEVH-box RNA helicase AtHELPS is important for responses to and tolerance of K+ deprivation in Arabidopsis [17]. The Arabidopsis DEAH-Box RNA helicase RID1 controls pre-mRNA splicing [18].
The recent availability of genome sequences enabled systematic investigation of the families of RNA helicase genes from Arabidopsis, rice, tomato, maize, and soybean. A total of 32 DEAD-box RNA helicases have been identified in Arabidopsis [19]; 113 and 115 RNA helicase genes have been found in Arabidopsis and rice, respectively [1]. Xu et al. [5] described the complete analysis and classification of 157 RNA helicase genes in tomato. Soon after, they reported a comparative genome-wide analysis of the RNA helicase gene families in Arabidopsis, rice, maize, and soybean [6]. The RNA helicase genes from Arabidopsis, rice, maize, and soybean were classified into three subfamilies, with the following respective numbers of genes for each subfamily: DEAD-box (50, 51, 57, and 87 genes), DEAH-box (40, 33, 31, and 48 genes), and DExD/H-box (71, 65, 50, and 78 genes) [6].
Cotton is an economically important crop grown worldwide as a source of fiber and edible oil [20,21]. RNA helicases are ubiquitous proteins that are believed to play important roles in cotton development and stress tolerance. However, only a few RNA helicases have been characterized in cotton [22,23]. Gossypium raimondii is a diploid species of cotton. Release of the G. raimondii genome enabled us to identify and analyze the family of RNA helicase genes in this species. Presently, there are a few reports on studies of cotton transcriptomics [22,24]. Here, we looked at expression of the different RNA helicases in transcriptome data obtained at Peking University (Beijing, China) and accessed through NCBI [25,26]. The present study identified 161 RNA helicase genes from the G. raimondii genome. These were classified into three subfamilies, which were defined by the presence of either a DEAD-box, DEAH-box, or DExD/H-box. Detailed information about the chromosomal locations, expansion, genomic structures, and phylogenetic relationships of the RNA helicase genes is provided. Transcript profiles of 161 RNA helicase genes in mature leaves and at the 0-day-post-anthesis (DPA) and 3-DPA ovule developmental stages were investigated using transcriptome-sequencing data. Our results show that most of these helicase proteins, especially GrDEADs, might function during the fiber initiation stage. This is the first report of a genome-wide analysis of the RNA helicase gene family of G. raimondii.

2. Results

2.1. Identification of RNA Helicase Family Genes in G. raimondii

RNA helicase usually contains a highly conserved adenosine triphosphate (ATP)-binding domain and a classical C-terminal domain [27–29]. To identify all members of the RNA helicase gene family in G. raimondii, all known Arabidopsis helicase gene sequences and all identified tomato RNA helicase gene sequences were used to query the protein database of G. raimondii using BLAST (basic local alignment search tool) analysis. This identified 161 genes, of which 151 genes had apparent helicase domains with the remaining 10 genes having either incomplete or no helicase domains. Given that the annotation of the G. raimondii genome indicated that these 10 genes encode helicase proteins, they were subjected to further analysis. This enabled us to classify all 161 putative helicase genes into three subfamilies based on their phylogenetic relationships and the structural features of the motif II region. The three subfamilies were defined by the presence of the DEAD-box (51 genes), DEAH-box (52 genes), and DExD/H-box (58 genes) motifs (Table 1, File S2). Based on their subfamily and the order in which they are found on chromosomes 1 through 13, the genes were renamed from GrDEAD1 through GrDEAD51, from GrDEAH1 through GrDEAH52, and from GrDExD/H1 through GrDExD/H58 (Table 1). The lengths, molecular weights, isoelectric points, and predicted subcellular localizations of each Gossypium raimondii helicase protein are summarized in Table 1. We found that GrDEAD genes were distinct from GrDEAH and GrDExD/H genes. Whereas the average GrDEAD protein contains 609 amino acids (aa), GrDEAH and GrDExD/H proteins each average approximately 1100 aa in length. The average theoretical isoelectric point (pI) of GrDEAD proteins, which is approximately 8.21, is higher than that for members of the two other families. Whereas, 74.5% of the helicases proteins analyzed were predicted to be located in the nucleus, 16.1% were predicted to reside in the cytoplasm, and others were mainly predicted to be located in the chloroplast and mitochondrion (Table 1).

2.2. Phylogenetic Analysis of the RNA Helicase Family Genes in G. raimondii

The phylogenetic tree of each subfamily in our study showed that the DEAD-box (Figure 1A), DEAH-box (Figure 1B) and DExD/H-box (Figure 1C) subfamilies could be further classified into four, six, and six subgroups, respectively. However, the available phylogenetic analysis of RNA helicases from Arabidopsis, rice, maize, and soybean [6] places these subfamilies into many more subclades. That study classified the DEAD-box, DEAH-box and DExD/H-box RNA helicase proteins from tomato into three, three, or five large subgroups, respectively. The diversity of these subclades indicates the extent of the variation of RNA helicase genes in plants.

2.3. Chromosomal Position and Gene Duplications

The 161 RNA helicase genes of G. raimondii are distributed across all 13 chromosomes, with different densities of their distribution along different chromosomes (Figure 2). For example, whereas chromosome 7 contained 19 RNA helicase genes, only five RNA helicase genes were annotated on chromosome 2. Based on the definition of gene clusters, a chromosomal region containing two or more genes within 200 kb can be considered to be a gene cluster [31,32]. In G. raimondii, 18 helicase genes were identified in nine clusters (GrDEAD34–35, GrDExD/H46–47, GrDEAD40-GrDExD/H39, GrDEAH10-GrDExD/H15, GrDEAH6-GrDExD/H12, GrDEAD27-GrDEAH22, GrDEAH27-GrDExD/H31, GrDEAD48-GrDExD/H52, and GrDEAD47-GrDExD/H5) that were dispersed across several chromosomes. Two pairs of genes (GrDEAD40-GrDExD/H39 and GrDEAH10-GrDExD/H15) were arranged in tandem repeats on chromosomes 5 and 10, respectively. In addition, recent studies have shown that G. raimondii has undergone the hexaploidization event (γ-WGD) shared by the eudicots and a cotton-specific whole genome duplication [25]. To analyze the relationship between the RNA helicase genes and genome-wide duplications, we mapped the G. raimondii helicase genes to the duplicated blocks. We found that 62 of the 161 helicase genes (38.5%) had syntenic relationships (Figure 2, File S1). Of these, 40 genes involved only two chromosome regions, nine (GrDEAD4-18-48, GrDExD/H4-30-54, and GrDEAH7-16-39) spanned three chromosome regions, eight (GrDEAD3-8-11-24 and GrDExD/H13-33-38-58) traversed four chromosome regions, and five (GrDEAD19-32-35-39-41) crossed five chromosome regions (Figure 2, File S1).
We also examined the orthologous relationships between helicase genes from Gossypium raimondii and tomato, given that orthologs often retain equivalent functions during the course of evolution [33]. We found that 66 (40.99%) helicase genes from Gossypium raimondii have one or several putative orthologs in tomato (Figure 3, File S1). Of these, 27, 16, and 23 were assigned to the DEAD, DEAH, and DExD/H-box subfamilies, respectively (File S1). One member in DEAD-box family, GrDEAD7, is an ortholog of SIDEAD27 in tomato and STRS1 in Arabidopsis. Whereas, GrDEAD37 is an ortholog of SIDEAD34 in tomato and LOS4 in Arabidopsis, GrDExD/H35 is an ortholog of SIDExD/H21 in tomato and ISE2 in Arabidopsis (File S1).

2.4. Structures and Domain Analysis of the Putative Helicase Genes

Figure 4 shows the exon-intron structures and conserved domains of putative RNA helicase genes in each subfamily. The number and location of introns varied among subfamilies. In general, compared with the two other subfamilies, DEAD family genes had simpler structures and more conserved structural patterns than members of the two other subfamilies. The relative levels of conservation are exemplified by comparison of the high level of conservation evident in the comparisons GrDEAD21–30, GrDEAD25–31, GrDEAD27–50, GrDEAD4-18-48, GrDEAD23-34-42, GrDEAD1-37-45, GrDEAD5-12-29-51, GrDEAD3-8-11-24, and GrDEAD19-32-35-39-41 (Figure 4A), with the lower level of conservation amongst members of the GrDEAH-box and GrDExD/H-box helicase family genes evident in the comparisons GrDEAH4-18, GrDEAH38-3, GrDEAH35-28, GrDEAH25-8, GrDEAH33-29, and GrDEAH7-16-39 (Figure 4B), as well as GrDExD/H5–36, GrDExD/H43-1, GrDExD/H20–48, GrDExD/H54-30-4, GrDExD/H28-3-23, and GrDExD/H38-58-33-13-44 (Figure 4C). Gene structures within the same subgroup of all three subfamilies were also very diverse. In addition, we found that whereas genes duplicated in different parts of the genome (such as combinations GrDEAD35-41-19 and GrDEAD4-48-18) had the same or similar gene structures, genes that had been duplicated in tandem (GrDEAD40-GrDExD/H39 and GrDEAH10-GrDExD/H15) had different gene structures, especially in terms of the numbers of exons. Domain analysis indicated the presence of a highly conserved ATP-binding domain and a classical C-terminal domain in almost all of the predicted RNA helicases. Additionally, we found a Q motif in all members of the DEAD family, except for GrDEAD16. A WW domain was observed in three members of the DEAD family. We also found that many of the DEAH and DExD/H family genes were surrounded by defined folds, such as the zf-RING, dsRBDs, and HSA domains.

2.5. Expression of RNA Helicase Family Genes in G. raimondii

The full-length RNA helicase protein sequences in G. raimondii were aligned by MUSCLE (v3.8.31) [30] and analyzed using maximum likelihood (ML) method. The attribution of proteins in Figure 5 is according to the attribution of proteins in Figure 1. Bootstrap values were calculated. The abundances of the transcripts that encode selected RNA helicases was examined at the fiber initiation stage and in mature leaves. The expression patterns of 161 RNA helicase family genes are shown in Figure 5. Of the 161 predicted genes, only six genes (those that encode GrDEAD1, GrDEAD22, GrDEAD30, GrDEAD45, GrDEAH3, and GrDExD/H15) were not expressed, whereas 141 genes were expressed both in ovules and leaves. The GrDEAD genes (Figure 5A) showed more homogenous levels of expression and an overall higher level of expression both in ovules and leaves when compared with GrDEAH genes (Figure 5B) and GrDExD/H genes (Figure 5C). Seven GrDEAD genes (GrDEAD13, GrDEAD19, GrDEAD21, GrDEAD28, GrDEAD32, GrDEAD35, and GrDEAD41), one GrDEAH gene (GrDEAH23), and three GrDExD/H genes (GrDExD/H33, GrDExD/H54, and GrDExD/H58) were expressed at high levels in all three samples. A member of the DEAD-box family, GrDEAD7, which is an ortholog of Arabidopsis STRS1 was highly expressed in 0-DPA ovules and mature leaves. More than half of the genes were mainly expressed in one of the development stages and tissues tested, with 50 genes (8 GrDEADs, 24 GrDEAHs, and 18 GrDExD/Hs) being the most abundant in 0-DPA ovules, 25 (7 GrDEADs, 5 GrDEAHs, and 13 GrDExD/Hs) being the most abundant in 3-DPA ovules, and 14 (3 GrDEADs, 6 GrDEAHs, and 5 GrDExD/Hs) being the most abundant in mature leaves. For example, GrDEAD37 (an ortholog of Arabidopsis LOS4) was highly expressed in 0-DPA ovules, whereas, GrDExD/H35 (an ortholog of Arabidopsis ISE2) was expressed in mature leaves. Whereas, the similarity of the abundance profiles of transcripts encoded by duplicated genes, such as GrDEAD19, GrDEAD35, and GrDEAD41, suggested functional redundancy between these family members, the different expression patterns of members of the GrDEAD4-18-48 indicated that other G. raimondii helicase genes have been preserved by sub-functionalization.
The expressions of helicase genes were very diverse, even within one subgroup. For example, in the GrDExD/H group II, six genes (GrDExD/H4, GrDExD/H30, GrDExD/H33, GrDExD/H38, GrDExD/H54, and GrDExD/H58) were highly expressed in all three samples and clustered together. However, eight genes, including GrDExD/H13 and GrDExD/H45, could only be detected in one to three organs, and were expressed at a low level in most organs. Nonetheless, given that as members in this group had a Q motif, they may be GrDEAD genes, which were assigned to the wrong subfamily owing to inadequate classification.

3. Discussion

RNA helicases play important roles in plant development and responses to stress. However, only a few RNA helicases have been identified in plants. Genome-wide analysis is the first step to elucidating the biological roles of members of the RNA helicase family members in certain plant species. The recent availability of genome sequences has enabled systematic investigation of this family of genes in Arabidopsis, rice, tomato, maize, and soybean. However, no RNA helicases have been characterized in cotton. This study involved a complete analysis of the RNA helicase gene family in the G. raimondii genome, including gene classification and the analysis of chromosomal locations, gene expansion, phylogenetic relationships, and structures of the genes, as well as their expression profiles at the fiber initiation stage and in mature leaves under normal growth conditions.

3.1. RNA Helicases in G. raimondii

We identified 161 RNA helicase genes in the diploid genome of the cotton species G. raimondii (Table 1). This is close to the estimated number of RNA helicase predicted from the analysis of the tomato (157), Arabidopsis (161), rice (149), maize (136), and soybean (213) genomes. The relatively large number of cotton RNA helicase genes identified likely reflects the detailed method we used to identify the RNA helicase genes. The presence of a large RNA helicase gene family in all of these species underscores that RNA helicases likely play important regulatory roles in various processes during plant growth and development. We classified these helicases into three subfamilies, which include the DEAD-box (51 genes), DEAH-box (52 genes), and DExD/H-box (58 genes) gene families (Table 1). Whereas, each of the three subfamilies of cotton RNA helicase genes consists of a similar number of genes, Xu et al. [6] have shown that the DEAD-box and DExD/H-box subfamily were larger than DEAH-box subfamily in Arabidopsis, rice, maize, and soybean. We also analyzed the predicted lengths, molecular weights, isoelectric points, and subcellular localizations of each putative helicase protein identified in the G. raimondii genome. We found that DEAD-box RNA helicases were distinct from DEAH-box and DExD/H-box RNA helicases. The apparently higher isoelectric points and smaller sizes of DEAD helicase proteins might be related to their relatively simple gene structures compared with other classes of RNA helicase. In addition, most of DEAD-box and DExD/H-box RNA helicase proteins were predicted to be located in the nucleus and cytoplasm while most DEAH-box RNA helicase proteins were predicted to reside in the nucleus (Table 1). Thus, we suggested DEAH-box RNA helicase proteins may mainly function in nuclear RNA processing. Linder et al. [34] have reported two DEAH-box RNA helicases, ESP3 and MUT6 which were located in nucleus perform different roles-RNA splicing and decay-in nuclear RNA processing. Several DEAD-box RNA helicase proteins and DExD/H-box RNA helicase proteins were predicted to be located in the chloroplast and mitochondria. Recently, there have been a few reports about RNA helicase proteins in chloroplast or mitochondria in plant. Asakura et al. [35] demonstrated that chloroplast RH3 DEAD-box RNA helicases in maize and Arabidopsis function in splicing of specific group II introns and affect chloroplast ribosome biogenesis. He et al. [36] demonstrated mitochondria ABO6 DExH-box RNA helicase in Arabidopsis involves in regulating the splicing of several genes of complex I.

3.2. Expansion of the RNA Helicase Gene Family in G. raimondii

Recent studies have shown that the G. raimondii genome has undergone at least two rounds of genome-wide duplication [25]. To detect possible relationships between RNA helicase genes and genome duplication events, we mapped 36 paralogous gene pairs (38.5%) of RNA helicase genes in G. raimondii (Figure 2, File S1). A similar percentage of paralogous pairs of RNA helicase genes was observed in Arabidopsis, rice, maize and soybean (35, 27, 25, and 62 pairs, respectively) [6]. These results suggest that the expansion of RNA helicase gene family is associated with whole-genome duplication events. However, the 38.5% value is much lower than that of the family of genes that encode NAC (NAM/ATAF/CUC) transcription factors (NAC-TFs) in G. raimondii. Of the 127 G. raimondii NAC-TF genes, 76.37% were within 307 identified syntenic blocks [37]. Given that segmental duplication events occur more often in the more slowly evolving gene families [38], the RNA helicase gene family may be evolving more rapidly than most other gene families in the cotton genome. In addition to the whole-genome duplication event, gene families can also arise through tandem amplification. For instance, in Chinese plum, tandem duplications played a key role in the expansion of the AP2/ERF family [39]. In G. raimondii, Shang et al. [37] also detected 20 tandem duplications associated with NAC-TF genes. However, only two gene pairs in our study showed evidence of having participated in tandem duplication. The mechanisms that supported the expansion of the RNA helicase gene family may be more complicated than is suggested by our classification; the specific mechanisms involved require further investigation. Marchat et al. [27] reported that the evolution of RNA helicase proteins involved gene fusion. Moreover, we speculate that given that the majority of RNA helicase gene family members play vital and diverse role in plants, there would be strong selection against variation in their copy numbers. The DEAD genes showed the greatest degree of duplication. Given the more simple and conserved gene structures and domains, as well as higher expression of members of this subfamily than those of others, we propose that DEAD helicase genes may evolve more slowly and may play more basic roles in plant growth and development than the members of other subfamilies of RNA helicases.

3.3. Phylogenetic Analysis and Gene Structural Organization

The full-length RNA helicase protein sequences in tomato and G. raimondii were aligned by MUSCLE and analyzed using the more accurate maximum likelihood (ML) method (File S2). Members of the DEAD-box, DEAH-box, and DExD/H-box subfamilies were further classified into four, six, and six subgroups, respectively (Figure 1), although Xu et al. [6] placed them into many more subclades following phylogenetic analysis of the RNA helicases of Arabidopsis, rice, maize, and soybean. Those workers classified the DEAD-box, DEAH-box, and DExD/H-box RNA helicase proteins from tomato into three, three, or five large subgroups, respectively. The diversity in the number and compositions of the subclades from different species indicate variation in the compositions of different RNA helicase gene families from different plant species. In addition, analysis of the exon-intron structures and sequences of the conserved domains can provide insights into the evolution of gene families. Our results showed that the numbers and locations of introns varied among subfamilies. Members of the DEAD-box subfamily had comparatively simple and conserved structural patterns. Genes from the GrDEAH-box and GrDExD/H-box subfamilies were less conserved and more diverse than those from the DEAD-box subfamily. It is noteworthy that genome-wide duplicated genes had the similar or same gene structures, whereas, tandem duplicated genes did not; this requires further analysis. Domain analysis indicated that an ATP-binding domain and a C-terminal domain were highly conserved in these putative RNA helicase proteins, whereas, more diversity was evident amongst the other domains present in each subfamily. The most characteristic feature of the DEAD-box family is the conserved Q motif [40]. We found Q motifs in all members of DEAD genes, except for GrDEAD16. In addition, three of the family members were found to have a WW domain. The WW domain has been implicated in mediating protein-protein interactions and linking cell signaling to the membrane cytoskeleton [41]. Whereas, DEAH genes lack a Q motif, most members of group II of DExD/H family have a Q motif. This might be attributed to the not-very-strict classification criteria used to distinguish between members of the DEAD and DExD/H subfamily. Many of the DEAH and DExD/H family genes were surrounded by defined folds, such as the zf-RING, dsRBD, and HAS domains, which may extend the length of the helicase. These regions influence or even define the function of a helicase [42]. The great diversity in these helicases may allow them to regulate many specific pathways in plants.

3.4. Expression Analysis Based on Transcriptome Sequencing Data

RNA helicases rearrange RNA secondary structure, potentially playing roles in any cellular process that involves RNA metabolism [3]. Of the 161 predicted genes in our study, 141 (87.6%) RNA helicase genes were expressed both in ovules and leaves. Xu et al. [6] reported that more than 80% RNA helicase genes in Arabidopsis, rice, and maize were expressed in at least one of the development stages and tissues tested. The high expression level of this RNA helicase gene family in all of these species further indicates that the RNA helicases may play important roles in various processes. The GrDEAD genes were more homogenous in terms of their level of expression and were expressed at higher levels both in ovules and leaves when compared with GrDEAH genes and GrDExD/H genes (Figure 5). The DEAD-box family member GrDEAD7, which is an ortholog of Arabidopsis STRS1, was highly expressed in 0-DPA ovules and mature leaves. The RNA helicases STRS1 and STRS2 have been shown to be involved in responses to various abiotic stresses [9]. GrDEAD genes may be more important in diverse cellular processes than GrDEAH genes and GrDExD/H genes. However, Xu et al. [6] have shown the DEAH-box RNA helicase genes higher proportion of the development stages and tissues in Arabidopsis, rice, and maize than DEAD-box RNA helicase genes and DExD/H-box RNA helicase genes. This phenomenon might be attributed to diversity between the different crops investigated or the diversity of the development stages and tissues examined. Moreover, tissue specificity of the expression of genes is commonly observed in plants. For example, cyclin dependent kinases-like proteins (CKL) of Arabidopsis show strong tissue specificity of expression, with CKL12 being root specific, CKL3 induced in stem tumors and callus, and CKL6 showing strong expression only in leaf tissue [43]. Though Arabidopsis ISE2 has been shown to be involved in posttranscriptional gene silencing, and the absence of the ISE2 affects a critical factor required for correct plasmodesmata (PD) formation and function [16], the duration of PD closure was positively correlated with the final fiber length attained [44]. Thus, GrDExD/H35 is not expressed at the initiation fiber stage.

4. Experimental Section

4.1. Identification of RNA Helicase Genes

The latest version (V2.0) of the Gossypium raimondii genome and protein sequences was downloaded from CottonGen [45]. To identify the members of the RNA helicase gene family in Gossypium raimondii, all known Arabidopsis helicase gene sequences and the identified tomato RNA helicase gene sequences [5] were used as queries to perform multiple database searches using BLASTP [46]. They were downloaded from The Arabidopsis Information Resource (TAIR) [47] and the plantGDB database [48], respectively. After selection of G. raimondii proteins with at least 50% identity with the query sequence, the candidate helicases proteins were aligned with each other to ensure that no gene was represented multiple times. All of the remaining protein sequences were examined using the domain analysis program PROSITE [49] with the default cutoff parameters. The three fields (length, molecular weight, and isoelectric point) of each Gossypium raimondii helicase protein were calculated using the online ExPasy program [50]. Subcellular localization was analyzed using the CELLO v2.5 server [51].

4.2. Phylogenetic Analysis

The full-length protein sequences of the helicase genes were aligned using the MUSCLE (v3.8.31) [30] program with the default settings. Phylogenetic trees were constructed by employing the maximum-likelihood (ML) method of the phyML (20120412) program [52] with the WAG (Whelan and Goldman) substitution model. Bootstrap values were calculated using the aLRT (approximate likelihood ratio test) model with the default cutoff parameters.

4.3. Chromosome Localization and Gene Duplications

Synteny analysis was conducted locally using a method similar to that developed for the Plant Genome Duplication Database [53]. Mcscan [54] was employed to identify homologous regions, and syntenic blocks were evaluated using Circos-0.64 [55]. Default parameters were used in all steps. Tandem duplication was characterized as multiple genes of one family located within the same or neighboring intergenic region [39].

4.4. Gene Structure and Domain Analysis

All of the putative protein sequences were analyzed using the domain analysis program PROSITE [56] at ExPASy [57]. Exon-intron structure information of these helicase genes was parsed from the GFF (Generic Feature Format) file downloaded along with the genomic data. Gene structures of the helicase genes were generated using the GSDS (Gene Structure Display Server) algorithm [58].

4.5. Gene Expression Analyses

The expression pattern of the helicase genes was analyzed using transcriptome sequencing data from mature leaves, 0-DPA ovules, and 3-DPA ovules of G. raimondii. These data were obtained from the NCBI Sequence Read Archive (SRA) [25]. The accession numbers were: SRX111367, SRX111365, and SRX111366, respectively. The search was performed using nucleotide signatures at least 20 nucleotide long. Reads mapping were performed by BWA (0.7.5a-r405) [59] with the default parameters except that the seed length was set to be 31. Sequenced reads that were mapped on these helicase sequences were converted to RPKM in order to estimate gene expression levels [60,61]. The formula used was:
R P K M = 10 6 C / ( N L / 10 3 )
where C is the number of reads that were uniquely aligned to the transcript, N is the total number of reads that were uniquely aligned to all the transcripts in a specific sample and L is number of bases in the transcript.

5. Conclusions

Our study has reported a genome-wide analysis of the important RNA helicase gene family in G. raimondii. Based on their expression analysis, we hypothesize that GrDEAD37 and GrDExD/H35 might function during the fiber initiation stage. Only 38.5% of the putative RNA helicase genes were mapped to the previously identified syntenic blocks. The specific mechanisms used during the expansion of the RNA helicase family might be more complicated than suggested by our analysis, and require further investigation. Of the subfamilies of RNA helicases from G. raimondii, the GrDEADs have undergone the greatest degree of duplication and have the most conserved structural patterns and highest levels of expression, when measured at the level of transcript abundance. This suggests that the GrDEAD gene subfamily may evolve more slowly than the two other subfamilies, and that GrDEAD genes play a more important and basic role in cotton than DEAH or DExD/H genes. This study should provide a solid foundation for future functional studies and for guiding future experimental work on helicase genes in plants.

Acknowledgments

This research was mainly supported by the China Major Projects for Transgenic Breeding (Grant Nos. 2011ZX08005-004 and 2011ZX08005-002) and the China Key Development Project for Basic Research (973) (Grant No. 2010CB12606).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. 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, 449–465. [Google Scholar]
  2. Lohman, T.M.; Bjornson, K.P. Mechanisms of helicase-catalyzed DNA unwinding. Annu. Rev. Biochem. 1996, 65, 169–214. [Google Scholar]
  3. Owttrim, G.W. RNA helicases and abiotic stress. Nucleic Acids Res. 2006, 34, 3220–3230. [Google Scholar]
  4. Schmid, S.; Linder, P. D-E-A-D protein family of putative RNA helicases. Mol. Microbiol. 1992, 6, 283–292. [Google Scholar]
  5. 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, 128–140. [Google Scholar]
  6. 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, e78982. [Google Scholar]
  7. Gong, Z.; Dong, C.H.; Lee, H.; Zhu, J.; Xiong, L.; Gong, D.; Stevenson, B.; Zhu, J.K. A DEAD box RNA helicase is essential for mrna export and important for development and stress responses in Arabidopsis. Plant Cell Online 2005, 17, 256–267. [Google Scholar]
  8. Gong, Z.; Lee, H.; Xiong, L.; Jagendorf, A.; Stevenson, B.; Zhu, J.K. RNA helicase-like protein as an early regulator of transcription factors for plant chilling and freezing tolerance. Proc. Natl. Acad. Sci. USA 2002, 99, 11507–11512. [Google Scholar]
  9. 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, 814–830. [Google Scholar]
  10. Guan, Q.; Wu, J.; Zhang, Y.; Jiang, C.; Liu, R.; Chai, C.; Zhu, J. A DEAD vox RNA helicase is critical for pre-mrna splicing cold-responsive gene regulation and cold tolerance in Arabidopsis. Plant Cell Online 2013, 25, 342–356. [Google Scholar]
  11. Wang, Y.; Duby, G.; Purnelle, B.; Boutry, M. Tobacco VDL gene encodes a plastid DEAD box RNA Helicase and is involved in chloroplast differentiation and plant morphogenesis. Plant Cell Online 2000, 12, 2129–2142. [Google Scholar]
  12. Gendra, E.; Moreno, A.; Albà, 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, 875–886. [Google Scholar]
  13. Macovei, A.; Vaid, N.; Tula, S.; Tuteja, N. A new DEAD-box helicase ATP-binding protein (OsABP) from rice is responsive to abiotic stress. Plant Signal. Behav. 2012, 7, 1138–1143. [Google Scholar]
  14. Liu, H.; Liu, J.; Fan, S.; Song, M.; Han, X.; Liu, F.; Shen, F. Molecular cloning and characterization of a salinity stress-induced gene encoding DEAD-box helicase from the halophyte Apocynum venetum. J. Exp. Bot. 2008, 59, 633–644. [Google Scholar]
  15. Lange, H.; Sement, F.M.; Gagliardi, D. MTR4 a putative RNA helicase and exosome co-factor is required for proper rRNA biogenesis and development in Arabidopsis thaliana. Plant J. 2011, 68, 51–63. [Google Scholar]
  16. Kobayashi, K.; Otegui, M.S.; Krishnakumar, S.; Mindrinos, M.; Zambryski, P. INCREASED SIZE EXCLUSION LIMIT2 encodes a putative DEVH Box RNA helicase involved in plasmodesmata function during Arabidopsis embryogenesis. Plant Cell Online 2007, 19, 1885–1897. [Google Scholar]
  17. Xu, R.R.; Qi, S.D.; Lu, L.T.; Chen, C.T.; Wu, C.A.; Zheng, C.C. A DExD/H box RNA helicase is important for K+ deprivation responses and tolerance in Arabidopsis thaliana. FEBS J. 2011, 278, 2296–2306. [Google Scholar]
  18. Ohtani, M.; Demura, T.; Sugiyama, M. Arabidopsis ROOT INITIATION DEFECTIVE1 a DEAH-Box RNA helicase involved in pre-mrna splicing is essential for plant development. Plant Cell Online 2013. [Google Scholar] [CrossRef]
  19. Aubourg, S.; Kreis, M.; Lecharny, A. The DEAD box RNA helicase family in Arabidopsis thaliana. Nucleic Acids Res. 1999, 27, 628–636. [Google Scholar]
  20. Rathore, K.S. Cotton. In Genetic Modification of Plants; Springer: Berlin/Heidelberg, Germany, 2010; pp. 269–285. [Google Scholar]
  21. Wendel, J.F.; Brubaker, C.; Alvarez, I.; Cronn, R.; Stewart, J.M. Evolution and natural history of the cotton genus. In Genetics and Genomics of Cotton; Springer: New York, NY, USA, 2009; pp. 3–22. [Google Scholar]
  22. Hovav, R.; Udall, J.A.; Hovav, E.; Rapp, R.; Flagel, L.; Wendel, J.F. A majority of cotton genes are expressed in single-celled fiber. Planta 2008, 227, 319–329. [Google Scholar]
  23. Hovav, R.; Udall, J.A.; Chaudhary, B.; Hovav, E.; Flagel, L.; Hu, G.; Wendel, J.F. The evolution of spinnable cotton fiber entailed prolonged development and a novel metabolism. PLoS Genet. 2008, 4, e25. [Google Scholar]
  24. Yao, D.; Zhang, X.; Zhao, X.; Liu, C.; Wang, C.; Zhang, Z.; Zhang, C.; Wei, Q.; Wang, Q.; Yan, H. Transcriptome analysis reveals salt-stress-regulated biological processes and key pathways in roots of cotton (Gossypium hirsutum L). Genomics 2011, 98, 47–55. [Google Scholar]
  25. Wang, K.; Wang, Z.; Li, F.; Ye, W.; Wang, J.; Song, G.; Yue, Z.; Cong, L.; Shang, H.; Zhu, S.; et al. The draft genome of a diploid cotton Gossypium raimondii. Nat. Genet. 2012, 44, 1098–1103. [Google Scholar]
  26. NCBI. Available online: http://www.ncbi.nlm.nih.gov/ (accessed on 1 October 2013).
  27. Marchat, L.A.; Orozco, E.; Guillen, N.; Weber, C.; Lopez-Camarillo, C. Putative DEAD and DExH-box RNA helicases families in Entamoeba histolytica. Gene 2008, 424, 1–10. [Google Scholar]
  28. Rocak, S.; Linder, P. DEAD-box proteins: the driving forces behind RNA metabolism. Nat. Rev. Mol. Cell Biol. 2004, 5, 232–241. [Google Scholar]
  29. Zou, J.; Chang, M.; Nie, P.; Secombes, C. Origin and evolution of the RIG-I like RNA helicase gene family. BMC Evol. Biol. 2009, 9, 85. [Google Scholar]
  30. Edgar, R.C. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar]
  31. Holub, E.B. The arms race is ancient history in Arabidopsis the wildflower. Nat. Rev. Genet. 2001, 2, 516–527. [Google Scholar]
  32. Zhang, C.; Zhang, H.; Zhao, Y.; Jiang, H.; Zhu, S.; Cheng, B.; Xiang, Y. Genome-wide analysis of the CCCH zinc finger gene family in Medicago truncatula. Plant Cell Rep. 2013, 1–13. [Google Scholar]
  33. Altenhoff, A.M.; Dessimoz, C. Phylogenetic and functional assessment of orthologs inference projects and methods. PLoS Comput. Biol. 2009, 5, e1000262. [Google Scholar]
  34. Linder, P.; Owttrim, G.W. Plant RNA helicases: Linking aberrant and silencing RNA. Trends Plant Sci. 2009, 14, 344–352. [Google Scholar]
  35. Asakura, Y.; Galarneau, E.; Watkins, K.P.; Barkan, A.; van Wijk, K.J. Chloroplast RH3 DEAD box RNA helicases in maize and Arabidopsis function in splicing of specific group II introns and affect chloroplast ribosome biogenesis. Plant Physiol. 2012, 159, 961–974. [Google Scholar]
  36. He, J.; Duan, Y.; Hua, D.; Fan, G.; Wang, L.; Liu, Y.; Chen, Z.; Han, L.; Qu, L.J.; Gong, Z. DEXH Box RNA helicase-mediated mitochondrial reactive oxygen species production in Arabidopsis mediates crosstalk between abscisic acid and auxin signaling. Plant Cell Online 2012, 24, 1815–1833. [Google Scholar]
  37. Shang, H.; Li, W.; Zou, C.; Yuan, Y. Analyses of the NAC transcription factor gene family in Gossypium raimondii Ulbr: Chromosomal location structure phylogeny and expression patterns. J. Integr. Plant Biol. 2013, 55, 663–676. [Google Scholar]
  38. Cannon, S.B.; Mitra, A.; Baumgarten, A.; Young, N.D.; May, G. The roles of segmental and tandem gene duplication in the evolution of large gene families in Arabidopsis thaliana. BMC Plant Biol. 2004, 4, 10. [Google Scholar]
  39. Du, D.; Hao, R.; Cheng, T.; Pan, H.; Yang, W.; Wang, J.; Zhang, Q. Genome-wide analysis of the AP2/ERF gene family in Prunus mume. Plant Mol. Biol. Rep. 2013, 31, 741–750. [Google Scholar]
  40. Tanner, N.K.; Cordin, O.; Banroques, J.; Doère, M.; Linder, P. The Q motif: A newly identified motif in DEAD Box helicases may regulate ATP binding and hydrolysis. Mol. Cell 2003, 11, 127–138. [Google Scholar]
  41. Ilsleya, J.L.; Sudolb, M.; Windera, S.J. The WW domain: Linking cell signalling to the membrane cytoskeleton. Cell. Signal. 2002, 14, 183–189. [Google Scholar]
  42. Klostermeier, D.; Rudolph, M.G. A novel dimerization motif in the C-terminal domain of the Thermus thermophilus DEAD box helicase Hera confers substantial flexibility. Nucleic Acids Res. 2009, 37, 421–430. [Google Scholar]
  43. Menges, M.; de Jager, S.M.; Gruissem, W.; Murray, J.A. Global analysis of the core cell cycle regulators of Arabidopsis identifies novel genes reveals multiple and highly specific profiles of expression and provides a coherent model for plant cell cycle control. Plant J. 2005, 41, 546–566. [Google Scholar]
  44. Lucas, W.J.; Ham, B.K.; Kim, J.Y. Plasmodesmata—Bridging the gap between neighboring plant cells. Trends Cell Biol. 2009, 19, 495–503. [Google Scholar]
  45. CottonGen. Available online: http://www.cottongen.org/ (accessed on 1 October 2013).
  46. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar]
  47. The Arabidopsis Information Resource. Available online: http://www.arabidopsis.org/ (accessed on 6 October 2013).
  48. PlantGDB. Available online: http://www.plantgdb.org/ (accessed on 6 October 2013).
  49. PROSITE Program. Available online: http://www.expasy.ch/prosite/ (accessed on 25 October 2013).
  50. Compute pI/Mw Tool. ExPASy Server. Available online: http://web.expasy.org/compute_pi/ (accessed on 25 October 2013).
  51. CELLO v2.5. Available online: http://cello.life.nctu.edu.tw/ (accessed on 25 October 2013).
  52. 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, W557–W559. [Google Scholar]
  53. Tang, H.; Bowers, J.E.; Wang, X.; Ming, R.; Alam, M.; Paterson, A.H. Synteny and collinearity in plant genomes. Science 2008, 320, 486–488. [Google Scholar]
  54. Wang, Y.; Tang, H.; DeBarry, J.D.; Tan, X.; Li, J.; Wang, X.; Lee, T.-H.; Jin, H.; Marler, B.; Guo, H. MCScanX: A toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012, 40, e49–e49. [Google Scholar]
  55. Krzywinski, M.; Schein, J.; Birol, İ; Connors, J.; Gascoyne, R.; Horsman, D.; Jones, S.J.; Marra, M.A. Circos: An information aesthetic for comparative genomics. Genome Res. 2009, 19, 1639–1645. [Google Scholar]
  56. Bairoch, A. PROSITE: A dictionary of sites and patterns in proteins. Nucleic Acids Res. 1991, 19, 2241. [Google Scholar]
  57. Gasteiger, E.; Gattiker, A.; Hoogland, C.; Ivanyi, I.; Appel, R.D.; Bairoch, A. ExPASy: The proteomics server for in-depth protein knowledge and analysis. Nucleic Acids Res. 2003, 31, 3784–3788. [Google Scholar]
  58. GSDS v2.0. Available online: http://gsds.cbi.pku.edu.cn/ (accessed on 25 October 2013).
  59. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar]
  60. Mortazavi, A.; Williams, B.A.; McCue, K.; Schaeffer, L.; Wold, B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat. Methods 2008, 5, 621–628. [Google Scholar]
  61. Yin, Z.; Wang, J.; Wang, D.; Fan, W.; Wang, S.; Ye, W. The MAPKKK gene family in Gossypium raimondii: Genome-wide identification classification and expression analysis. Int. J. Mol. Sci. 2013, 14, 18740–18757. [Google Scholar]
Ijms 15 04635f1 1024
Figure 1. Phylogenetic analysis of RNA helicases in tomato and G. raimondii. (A) The DEAD-box RNA helicase proteins; (B) The DEAH-box RNA helicase proteins; (C) The DExD/H-box RNA helicase proteins. The amino acid sequences of the RNA helicase proteins were aligned with MUSCLE (v3.8.31) [30], and the phylogenetic tree was constructed using the maximum-likelihood (ML) method.
Figure 1. Phylogenetic analysis of RNA helicases in tomato and G. raimondii. (A) The DEAD-box RNA helicase proteins; (B) The DEAH-box RNA helicase proteins; (C) The DExD/H-box RNA helicase proteins. The amino acid sequences of the RNA helicase proteins were aligned with MUSCLE (v3.8.31) [30], and the phylogenetic tree was constructed using the maximum-likelihood (ML) method.
Ijms 15 04635f1
Ijms 15 04635f2 1024
Figure 2. Syntenic relationships among RNA helicases in G. raimondii. The G. raimondii RNA helicases genes are shown in the outer circle, with DEAD-box, DEAH-box and DExD/H-box RNA helicase genes indicated in red, green, and blue, respectively. Genome-wide duplicated genes in the DEAD-box, DEAH-box, and DExD/H-box subfamilies are connected by red, green, and blue lines, respectively.
Figure 2. Syntenic relationships among RNA helicases in G. raimondii. The G. raimondii RNA helicases genes are shown in the outer circle, with DEAD-box, DEAH-box and DExD/H-box RNA helicase genes indicated in red, green, and blue, respectively. Genome-wide duplicated genes in the DEAD-box, DEAH-box, and DExD/H-box subfamilies are connected by red, green, and blue lines, respectively.
Ijms 15 04635f2
Ijms 15 04635f3 1024
Figure 3. Syntenic relationships between RNA helicases of tomato and G. raimondii. ch, chromosome of tomato; Chr, chromosome of G. raimondii. DEAD-box, DEAH-box, and DExD/H-box genes are indicated in red, green, and blue, respectively. The putative orthologous RNA helicase genes belonging to the DEAD-box, DEAH-box, and DExD/H-box subfamilies are connected by red, green, and blue lines, respectively.
Figure 3. Syntenic relationships between RNA helicases of tomato and G. raimondii. ch, chromosome of tomato; Chr, chromosome of G. raimondii. DEAD-box, DEAH-box, and DExD/H-box genes are indicated in red, green, and blue, respectively. The putative orthologous RNA helicase genes belonging to the DEAD-box, DEAH-box, and DExD/H-box subfamilies are connected by red, green, and blue lines, respectively.
Ijms 15 04635f3
Ijms 15 04635f4a 1024Ijms 15 04635f4b 1024
Figure 4. Gene structure and domain analysis of RNA helicases in G. raimondii. (A) The DEAD-box RNA helicase genes; (B) The DEAH-box RNA helicase genes; (C) The DExD/H-box RNA helicase genes. Exons and introns of each subgroup are shown in the yellow boxes on the left. Conserved domains annotated by PROSITE are shown in colored boxes on the right.
Figure 4. Gene structure and domain analysis of RNA helicases in G. raimondii. (A) The DEAD-box RNA helicase genes; (B) The DEAH-box RNA helicase genes; (C) The DExD/H-box RNA helicase genes. Exons and introns of each subgroup are shown in the yellow boxes on the left. Conserved domains annotated by PROSITE are shown in colored boxes on the right.
Ijms 15 04635f4aIjms 15 04635f4b
Ijms 15 04635f5 1024
Figure 5. Expression patterns of RNA helicases in Gossypium raimondii. (A) The DEAD-box RNA helicase proteins; (B) The DEAH-box RNA helicase protein; (C) The DExD/H-box RNA helicase proteins. Heatmap showing the clustering of RNA helicase genes across three tissues (leaf, ovules at 0 DPA, and ovules at 3 DPA; mentioned at the top of each lane). RNA seq data (series accession number SRA048621) was obtained from the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) database [26]. The color scale at the top of the dendrogram indicates the relative expression levels. Color scale represents reads per kilobase per million normalized log2 transformed counts where light green indicates low level and red indicates high level.
Figure 5. Expression patterns of RNA helicases in Gossypium raimondii. (A) The DEAD-box RNA helicase proteins; (B) The DEAH-box RNA helicase protein; (C) The DExD/H-box RNA helicase proteins. Heatmap showing the clustering of RNA helicase genes across three tissues (leaf, ovules at 0 DPA, and ovules at 3 DPA; mentioned at the top of each lane). RNA seq data (series accession number SRA048621) was obtained from the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) database [26]. The color scale at the top of the dendrogram indicates the relative expression levels. Color scale represents reads per kilobase per million normalized log2 transformed counts where light green indicates low level and red indicates high level.
Ijms 15 04635f5
Table
Table 1. Characteristics of RNA helicases from G. raimondii. AA: Amino acid; pI: The theoretical isoelectric point of proteins; Mw: The theoretical molecular weight of proteins.
Table 1. Characteristics of RNA helicases from G. raimondii. AA: Amino acid; pI: The theoretical isoelectric point of proteins; Mw: The theoretical molecular weight of proteins.
Gene nameGene identifierGenomic positionSize (AA)MwpISubcellular localization
GrDEAD1Gorai.001G008400.1Chr01: 784950-78726247252.956.12Cytoplasmic
GrDEAD2Gorai.001G064600.1Chr01: 6445636-645139360468.988.93Nuclear
GrDEAD3Gorai.001G169800.1Chr01: 24289833-2429527549957.138.66Nuclear
GrDEAD4Gorai.001G263500.1Chr01: 53921207-5392639965570.989.19Nuclear
GrDEAD5Gorai.002G232200.1Chr02: 59155342-5915929062067.948.4Nuclear
GrDEAD6Gorai.003G106600.1Chr03: 32786814-3279382862271.268.96Nuclear
GrDEAD7Gorai.003G186400.1Chr03: 45681420-4568631756462.898.97Nuclear
GrDEAD8Gorai.004G042200.1Chr04: 3636522-364265849056.28.84Nuclear
GrDEAD9Gorai.004G050000.1Chr04: 4582045-458620677685.095.72Nuclear
GrDEAD10Gorai.004G089100.1Chr04: 11968779-1197424045049.776.53Cytoplasmic
GrDEAD11Gorai.004G157400.1Chr04: 44565986-4457100549656.998.72Nuclear
GrDEAD12Gorai.004G277500.1Chr04: 61046063-6105015562468.336.93Chloroplast
GrDEAD13Gorai.005G021200.1Chr05: 1744957-174787268378.28.94Nuclear
GrDEAD14Gorai.005G070600.1Chr05: 7665522-76732951154126.0910.07Nuclear
GrDEAD15Gorai.005G084500.1Chr05: 10403322-1040915478889.029.82Nuclear
GrDEAD16Gorai.005G118500.1Chr05: 24460627-2446477330033.49.26Cytoplasmic
GrDEAD17Gorai.005G186000.1Chr05: 54468800-5447465258063.667.68Chloroplast
GrDEAD18Gorai.006G031200.1Chr06: 8129011-813400669276.229.66Nuclear
GrDEAD19Gorai.006G095400.1Chr06: 33433098-3343636641346.845.28Cytoplasmic
GrDEAD20Gorai.006G254000.1Chr06: 49690002-4969578481491.899.16Nuclear
GrDEAD21Gorai.007G009400.1Chr07: 733772-73831369775.738.99Nuclear
GrDEAD22Gorai.007G029000.1Chr07: 1983857-198937234939.228.92Cytoplasmic
GrDEAD23Gorai.007G047000.1Chr07: 3268271-327269362368.339.83Nuclear
GrDEAD24Gorai.007G097600.1Chr07: 7210939-721719349456.518.4Nuclear
GrDEAD25Gorai.007G107800.1Chr07: 8134802-813969379288.919.46Nuclear
GrDEAD26Gorai.007G226800.1Chr07: 27332174-2733756461568.619.2Nuclear
GrDEAD27Gorai.007G304300.1Chr07: 51812182-5181839544449.548.97Cytoplasmic
GrDEAD28Gorai.007G309800.1Chr07: 52334743-5234031782090.348.33Chloroplast
GrDEAD29Gorai.007G362700.1Chr07: 59429333-5943403261767.378.16Nuclear
GrDEAD30Gorai.008G282800.1Chr08: 55902284-5590692370577.129.54Nuclear
GrDEAD31Gorai.008G245100.1Chr08: 52991359-5299614283594.629.51Nuclear
GrDEAD32Gorai.009G043400.1Chr09: 3162403-316605341246.765.37Cytoplasmic
GrDEAD33Gorai.009G126600.1Chr09: 9548083-955119458765.96.7Cytoplasmic
GrDEAD34Gorai.009G140700.1Chr09: 10637675-1064293450555.888.97Cytoplasmic
GrDEAD35Gorai.009G142800.1Chr09: 10823077-1082673641347.015.29Cytoplasmic
GrDEAD36Gorai.009G193900.1Chr09: 14896221-1490372783394.117.12Nuclear
GrDEAD37Gorai.009G410800.1Chr09: 61968396-6197255550155.795.65Cytoplasmic
GrDEAD38Gorai.009G437700.1Chr09: 68810677-6881531875284.978.96Nuclear
GrDEAD39Gorai.010G098500.1Chr10: 17146619-1714961541346.775.37Cytoplasmic
GrDEAD40Gorai.010G134100.1Chr10: 29846928-2985285562068.479.87Mitochondrial
GrDEAD41Gorai.011G063400.1Chr11: 5308916-531226541346.925.29Cytoplasmic
GrDEAD42Gorai.011G066000.1Chr11: 5584414-559053850555.729.04Mitochondrial
GrDEAD43Gorai.011G148000.1Chr11: 24134100-2413611666176.368.62Nuclear
GrDEAD44Gorai.011G261800.1Chr11: 59292567-592986421104125.795.78Nuclear
GrDEAD45Gorai.011G288300.1Chr11: 61987801-6199144949355.026.3Cytoplasmic
GrDEAD46Gorai.012G015300.1Chr12: 1700644-170544359868.179.48Nuclear
GrDEAD47Gorai.012G039900.1Chr12: 4943976-495034064369.169.51Nuclear
GrDEAD48Gorai.012G063000.1Chr12: 8891463-889664368875.169.8Nuclear
GrDEAD49Gorai.012G079600.1Chr12: 12903543-1290793259366.239.11Mitochondrial
GrDEAD50Gorai.013G142600.1Chr13: 38926934-3893227344449.598.97Cytoplasmic
GrDEAD51Gorai.013G247600.1Chr13: 56620576-5662597261666.798.66Nuclear
GrDEAH1Gorai.001G141100.1Chr01: 19085241-190935001328148.921.97Nuclear
GrDEAH2Gorai.002G190800.1Chr02: 51591530-51599433979109.029.01Nuclear
GrDEAH3Gorai.003G041600.1Chr03: 5046943-50542541750196.936.75Nuclear
GrDEAH4Gorai.004G081200.1Chr04: 9735695-974142875885.986.61Nuclear
GrDEAH5Gorai.004G164900.1Chr04: 45891531-459013801615182.547.34Nuclear
GrDEAH6Gorai.004G265100.1Chr04: 60082397-6009203271180.986.59Nuclear
GrDEAH7Gorai.005G058400.1Chr05: 6046005-60526921063123.425.87Nuclear
GrDEAH8Gorai.005G157800.1Chr05: 45129943-4513709873483.315.09Nuclear
GrDEAH9Gorai.005G161600.1Chr05: 46720561-467348001028114.745.99Nuclear
GrDEAH10Gorai.005G170100.1Chr05: 49756283-4976426170078.717.56Nuclear
GrDEAH11Gorai.005G213400.1Chr05: 59588117-595965961067119.046.37Nuclear
GrDEAH12Gorai.006G000100.1Chr06: 50440-5822982593.928.85Cytoplasmic
GrDEAH13Gorai.006G125500.1Chr06: 37710924-377275821747199.385.8Nuclear
GrDEAH14Gorai.006G135500.1Chr06: 39140704-391516721536174.798.3Nuclear
GrDEAH15Gorai.006G248800.1Chr06: 49317977-493233651112126.268.08Nuclear
GrDEAH16Gorai.006G264900.1Chr06: 50379562-503860541065123.475.93Nuclear
GrDEAH17Gorai.007G040200.1Chr07: 2796320-28061891142128.147.07PlasmaMembrane
GrDEAH18Gorai.007G136500.1Chr07: 11144557-1115042775886.087.26Nuclear
GrDEAH19Gorai.007G190100.1Chr07: 18534829-1854084573281.939.15Nuclear
GrDEAH20Gorai.007G273900.1Chr07: 46703371-467195561484168.815.48Nuclear
GrDEAH21Gorai.007G299200.1Chr07: 51082328-5108809988299.768.49Nuclear
GrDEAH22Gorai.007G305400.1Chr07: 51890852-5189567661168.078.18Nuclear
GrDEAH23Gorai.008G137800.1Chr08: 38782533-3878618272582.667.22Nuclear
GrDEAH24Gorai.008G154400.1Chr08: 41323487-41330255885100.527.91Nuclear
GrDEAH25Gorai.008G168800.1Chr08: 44099126-4410786973182.915.38Nuclear
GrDEAH26Gorai.008G176500.1Chr08: 45286639-4529451771880.58.88Nuclear
GrDEAH27Gorai.008G197200.1Chr08: 48161550-481678781067121.285.64Nuclear
GrDEAH28Gorai.009G117300.1Chr09: 8628265-86330631189135.536.56Nuclear
GrDEAH29Gorai.009G171000.1Chr09: 13221402-132300431391153.375.43Nuclear
GrDEAH30Gorai.009G186700.1Chr09: 14354401-143652071234137.838.29Nuclear
GrDEAH31Gorai.009G267400.1Chr09: 22217004-222207221006113.775.31Nuclear
GrDEAH32Gorai.010G093100.1Chr10: 15093632-151049331675187.556.05Nuclear
GrDEAH33Gorai.010G151100.1Chr10: 40671899-406804921325146.145.3Nuclear
GrDEAH34Gorai.010G228200.1Chr10: 59997688-6001389310211158.53Nuclear
GrDEAH35Gorai.011G007600.1Chr11: 569884-5744341184135.346.68Nuclear
GrDEAH36Gorai.011G032400.1Chr11: 2409712-241702269777.958.37PlasmaMembrane
GrDEAH37Gorai.011G075000.1Chr11: 7156546-7161392913103.077.79Nuclear
GrDEAH38Gorai.011G111000.1Chr11: 13326515-133342291760197.967.62Nuclear
GrDEAH39Gorai.011G187000.1Chr11: 44528408-445351611060123.25.95Nuclear
GrDEAH40Gorai.011G206500.1Chr11: 49795790-498144741767199.657.03Nuclear
GrDEAH41Gorai.012G017400.1Chr12: 1996523-20037471345151.148.44Nuclear
GrDEAH42Gorai.012G136000.1Chr12: 30913260-309206331038114.766.6Nuclear
GrDEAH43Gorai.012G176400.1Chr12: 34564794-3456732965874.126.1Nuclear
GrDEAH44Gorai.013G001200.1Chr13: 68544-7318086396.218.9Chloroplast
GrDEAH45Gorai.013G038600.1Chr13: 3074495-30837141203135.397.84Nuclear
GrDEAH46Gorai.013G043100.1Chr13: 3742541-3750413971107.926.2Nuclear
GrDEAH47Gorai.013G133100.1Chr13: 35035921-35039609990112.455.64Nuclear
GrDEAH48Gorai.013G139300.1Chr13: 37542094-375550831371152.685.82Nuclear
GrDEAH49Gorai.013G141900.1Chr13: 38687077-386953771232139.186.13Nuclear
GrDEAH50Gorai.013G218600.1Chr13: 53839220-538543302161238.256.56Nuclear
GrDEAH51Gorai.013G252300.1Chr13: 56979166-569922541773202.645.77Nuclear
GrDEAH52Gorai.013G258600.1Chr13: 57397754-574062401037115.477.45Nuclear
GrDExD/H1Gorai.001G036900.1Chr01: 3403877-34202371455165.335.25Nuclear
GrDExD/H2Gorai.001G150400.1Chr01: 20836533-208461781950218.725.96Nuclear
GrDExD/H3Gorai.001G168700.1Chr01: 24030814-240419492238251.249Nuclear
GrDExD/H4Gorai.002G016100.1Chr02: 1050235-105335440545.95.99Cytoplasmic
GrDExD/H5Gorai.002G143000.1Chr02: 25828154-258390651393158.076.64PlasmaMembrane
GrDExD/H6Gorai.002G216000.1Chr02: 56556123-565737451470163.866.4Nuclear
GrDExD/H7Gorai.003G100200.1Chr03: 30861094-308750471287144.526.5Nuclear
GrDExD/H8Gorai.003G125100.1Chr03: 37290357-373027571340149.475.67Cytoplasmic
GrDExD/H9Gorai.003G156900.1Chr03: 42640922-4264703979690.088.15Nuclear
GrDExD/H10Gorai.004G209000.1Chr04: 54157901-5416065354859.429.63Mitochondrial
GrDExD/H11Gorai.004G221400.1Chr04: 55605032-556109801225137.546.4Nuclear
GrDExD/H12Gorai.004G264200.1Chr04: 59938105-599568632042231.375.47Nuclear
GrDExD/H13Gorai.004G292500.1Chr04: 62043821-6204842142848.465.4Cytoplasmic
GrDExD/H14Gorai.005G137100.1Chr05: 35497500-35505536926104.628.77Nuclear
GrDExD/H15Gorai.005G170000.1Chr05: 49691064-4969321540145.329.41Mitochondrial
GrDExD/H16Gorai.005G194600.1Chr05: 56476340-564949841217137.047.12Nuclear
GrDExD/H17Gorai.006G005000.1Chr06: 967752-9812011063120.145.98Nuclear
GrDExD/H18Gorai.006G017900.1Chr06: 4179363-418425085197.236.43Nuclear
GrDExD/H19Gorai.006G241400.1Chr06: 48762774-4876753756363.28.72Mitochondrial
GrDExD/H20Gorai.006G257000.1Chr06: 49880967-498966982586283.316.46Nuclear
GrDExD/H21Gorai.007G006300.1Chr07: 523949-533327119614.448.96Nuclear
GrDExD/H22Gorai.007G023800.1Chr07: 1694751-169963577086.868.08Nuclear
GrDExD/H23Gorai.007G100800.1Chr07: 7468905-74801352258253.198.82Nuclear
GrDExD/H24Gorai.007G374300.1Chr07: 60633068-6063709153259.729.39Nuclear
GrDExD/H25Gorai.008G031000.1Chr08: 3714524-3722015930105.316.29Nuclear
GrDExD/H26Gorai.008G049300.1Chr08: 6883759-68985772377263.017.2Nuclear
GrDExD/H27Gorai.008G068400.1Chr08: 11250164-1125478474885.765.76Nuclear
GrDExD/H28Gorai.008G121100.1Chr08: 35933836-359447752250252.238.95Nuclear
GrDExD/H29Gorai.008G145400.1Chr08: 39788075-397950651017115.156.02PlasmaMembrane
GrDExD/H30Gorai.008G149500.1Chr08: 40485456-4048871840846.165.98Nuclear
GrDExD/H31Gorai.008G198400.1Chr08: 48292290-483211422090236.596.08Nuclear
GrDExD/H32Gorai.009G053300.1Chr09: 3872981-38865151654186.186.7Nuclear
GrDExD/H33Gorai.009G057600.1Chr09: 4145298-414968542848.495.41Cytoplasmic
GrDExD/H34Gorai.009G153000.1Chr09: 11690907-1169488051457.399.54Nuclear
GrDExD/H35Gorai.009G345500.1Chr09: 41630420-416371651179133.415.3Nuclear
GrDExD/H36Gorai.009G374900.1Chr09: 50970390-509794911395158.666.56PlasmaMembrane
GrDExD/H37Gorai.009G393800.1Chr09: 54260282-542697641035115.378.81Nuclear
GrDExD/H38Gorai.010G114300.1Chr10: 22016499-2202200442848.445.41Cytoplasmic
GrDExD/H39Gorai.010G134200.1Chr10: 29886669-2988952040945.739.2Nuclear
GrDExD/H40Gorai.010G135100.1Chr10: 30499026-30514000688778.43Cytoplasmic
GrDExD/H41Gorai.010G175700.1Chr10: 51517626-51529268991112.566.05Nuclear
GrDExD/H42Gorai.010G194700.1Chr10: 55411136-554180531268146.358.06Nuclear
GrDExD/H43Gorai.010G244100.1Chr10: 61318853-613326141462166.455.32Nuclear
GrDExD/H44Gorai.011G082500.1Chr11: 8316574-832187342748.385.4Cytoplasmic
GrDExD/H45Gorai.011G088700.1Chr11: 9302151-930591053558.48.55Nuclear
GrDExD/H46Gorai.011G165900.1Chr11: 32740783-3274500548154.199Nuclear
GrDExD/H47Gorai.011G166100.1Chr11: 32787373-327986641114128.036.23Nuclear
GrDExD/H48Gorai.011G185100.1Chr11: 43975253-439928513361365.55.33Nuclear
GrDExD/H49Gorai.011G217300.1Chr11: 52330942-523388042177247.235.53Cytoplasmic
GrDExD/H50Gorai.012G023100.1Chr12: 2830126-28460191138128.768.88Nuclear
GrDExD/H51Gorai.012G039300.1Chr12: 4882407-48921161221134.916.11Nuclear
GrDExD/H52Gorai.012G063400.1Chr12: 8945198-89529501210135.58.5Nuclear
GrDExD/H53Gorai.013G034700.1Chr13: 2677455-2684701982111.627.1Nuclear
GrDExD/H54Gorai.013G070700.1Chr13: 8189231-819286740645.985.98Nuclear
GrDExD/H55Gorai.013G146500.1Chr13: 40255533-402788591256139.067.25Nuclear
GrDExD/H56Gorai.013G197700.1Chr13: 50405230-504252372054227.785.83Nuclear
GrDExD/H57Gorai.013G211300.1Chr13: 52317728-52328526990111.416.26Cytoplasmic
GrDExD/H58Gorai.013G215200.1Chr13: 53550658-5355565742848.445.4Cytoplasmic

Share and Cite

MDPI and ACS Style

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, 4635-4656. https://doi.org/10.3390/ijms15034635

AMA Style

Chen J, Zhang Y, Liu J, Xia M, Wang W, Shen F. Genome-Wide Analysis of the RNA Helicase Gene Family in Gossypium raimondii . International Journal of Molecular Sciences. 2014; 15(3):4635-4656. https://doi.org/10.3390/ijms15034635

Chicago/Turabian Style

Chen, Jie, Yujuan Zhang, Jubo Liu, Minxuan Xia, Wei Wang, and Fafu Shen. 2014. "Genome-Wide Analysis of the RNA Helicase Gene Family in Gossypium raimondii " International Journal of Molecular Sciences 15, no. 3: 4635-4656. https://doi.org/10.3390/ijms15034635

APA Style

Chen, J., Zhang, Y., Liu, J., Xia, M., Wang, W., & Shen, F. (2014). Genome-Wide Analysis of the RNA Helicase Gene Family in Gossypium raimondii . International Journal of Molecular Sciences, 15(3), 4635-4656. https://doi.org/10.3390/ijms15034635

Article Metrics

Back to TopTop