Skip to main content
  • Research article
  • Open access
  • Published:

Genome-wide analysis of CCCH zinc finger family in Arabidopsis and rice

Abstract

Background

Genes in the CCCH family encode zinc finger proteins containing the motif with three cysteines and one histidine residues. They have been known to play important roles in RNA processing as RNA-binding proteins in animals. To date, few plant CCCH proteins have been studied functionally.

Results

In this study, a comprehensive computational analysis identified 68 and 67 CCCH family genes in Arabidopsis and rice, respectively. A complete overview of this gene family in Arabidopsis was presented, including the gene structures, phylogeny, protein motifs, and chromosome locations. In addition, a comparative analysis between these genes in Arabidopsis and rice was performed. These results revealed that the CCCH families in Arabidopsis and rice were divided into 11 and 8 subfamilies, respectively. The gene duplication contributed to the expansion of the CCCH gene family in Arabidopsis genome. Expression studies indicated that CCCH proteins exhibit a variety of expression patterns, suggesting diverse functions. Finally, evolutionary analysis showed that one subfamily is higher plant specific. The expression profile indicated that most members of this subfamily are regulated by abiotic or biotic stresses, suggesting that they could have an effective role in stress tolerance.

Conclusion

Our comparative genomics analysis of CCCH genes and encoded proteins in two model plant species provides the first step towards the functional dissection of this emerging family of potential RNA-binding proteins.

Background

Transcription factors are important regulators of cellular processes, and the complexity of living organisms necessitates a large number of transcription factors. The zinc finger motifs, which are classified based on the arrangement of the zinc-binding amino acids, are present in many transcription factors and play critical roles in interactions with other molecules [1, 2]. A large amount of zinc-finger transcription factors are implicated in important biological processes and many of them often share common characteristic to form a family. So far, several zinc finger families have been found in plants, such as RING-finger, ERF, WRKY, DOF and LIM family [3–8]. However, most of them are identified as DNA-binding or protein-binding proteins, fewer function as RNA-binding proteins.

The CCCH zinc finger motif has been found in proteins from organisms ranging from man to yeast [9–18]. The CCCH proteins are a large family of zinc finger containing C3H-type motifs and many evidences proved that they may be RNA-binding proteins functioning in RNA processing [19–21]. In mouse, tristetraprolin, a protein containing two CCCH zinc fingers, binds directly to AU-rich elements within the 3'-untranslated region of target transcripts to facilitate mRNA degradation [22–25]. Zfp36l2, like its better-known relative TTP, is a mRNA-binding and destabilizing protein, functions in the physiological control of female fertility at the level of early embryonic development [26]. The PIE-1 is an essential regulator of Caenorhabditis elegans germ cell fate that segregates with the germ lineage by inhibition of transcription or activation of protein expression from maternal RNAs [27]. Compared to animals, only a few of plant CCCH proteins have been characterized functionally. In Arabidopsis thaliana, HUA1 has been proved to be a RNA-binding protein and likely participates in a new regulatory mechanism governing flower development [19]. AtCPSF30 encodes a small polypeptide which shares the probable ortholog of the 30-kD subunit of the mammalian cleavage and polyadenylation specificity factor, and it was shown to be nucleus-localized RNA-binding protein that binds calmodulin [20]. PEI1 is an embryo-specific CCCH zinc finger protein that plays an important role during Arabidopsis embryogenesis, functioning primarily in the apical domain of the embryo [28]. FES1 interacts genetically with FRI and FRL1 to promote the winter-annual habit of Arabidopsis and might be involved in the processing of mRNA [29]. In rice (Oryza sativa), OsDOS is a nuclear protein that delays leaf senescence by integrating developmental cues to the jasmonate (JA) pathway [30]. It was also proposed to play a role in posttranscriptional level by interacting with target RNA.

Given the potential for CCCH proteins to function diverse roles by associating with RNA and most of members containing CCCH motif remaining poorly understood, it was of considerable interest to us to characterize the plant CCCH gene family. The availability of the Arabidopsis and Oryza sativa genome sequences allows the genome-wide comparative analysis of gene families between monocot and eudicot plants [31–35].

Perl is one of the most widely used programming languages for managing and manipulating life-science information [36]. Perl has been extremely successful for connecting software applications together into sequence analysis pipelines, converting file formats, and extracting information from the output of analysis programs and other text files [37, 38].

In this study, we wrote Perl program to search against the entire genome of two species and performed a genome-wide analysis of the CCCH gene family in Oryza sativa and Arabidopsis. A total of 68 CCCH genes from Arabidopsis and 67 from Oryza sativa were identified and most of them have not been reported previously. Phylogenetic analyses revealed that the CCCH families in Arabidopsis and Oryza sativa were divided into 11 and 8 subfamilies, respectively. Expression studies indicated that CCCH proteins exhibit a variety of expression patterns, suggesting diverse functions. Detailed analysis of the higher plant-specific subfamily IX of Arabidopsis revealed its involvement in response to various stresses. Our genomics analysis provides the framework for future studies to dissect the function of this emerging family of potential RNA-binding proteins.

Results and discussion

Identification of genes coding for CCCH zinc finger protein in Arabidopsis

The previous research revealed that CCCH zinc finger protein contained 1–6 copies of C3H-type zinc-finger motifs characterized by three Cys and one His (Figure 1). Berg et al. defined that the CCCH family is a group of zinc-finger protein consisting of canonical C-X-C-X-C-X-H motif (C-X6–14-C-X4–5-C-X3-H) [21]. Based on the different amino acid spacing numbers between Cys and His in zinc finger motif, we divided the family into 18 groups (Table 1). The protein containing the motif within these groups was considered as a candidate member of CCCH family.

Figure 1
figure 1

Alignment of the CCCH zinc finger motifs from selected CCCH proteins. Black and light grey shading indicate identical and conserved amino acid residues, respectively. The three cysteine and one histidine residues putatively responsible for the zinc-finger structure are indicated.

Table 1 The potential CCCH motifs.

To uncover the entire family of genes coding for CCCH zinc finger protein in Arabidopsis, we analyzed the proteome data with a program written by Perl (Figure 2, see Additional file 1, 2, 3), in which a regular expression (C\w{6,14}C\w{4,5}C\w{3}H) was used to search for the CCCH motifs. Sequences that matched the criteria within 18 groups could be easily detected by the program and were considered as putative CCCH proteins.

Figure 2
figure 2

The program for detecting the CCCH motifs by regular expression. The code was used to detect the proteins containing CCCH motifs from local MySQL database. The whole program is available in additional file (see Additional file 5 and 7).

In addition, to verify the reliability of our result, we performed multiple BLAST algorithms to search against the corresponding data set using the known CCCH proteins got from several database as our query (see Methods). Another program written by Perl was used to parse the results and exclude the redundant entries from the initial data set.

The Perl program could identify all the amino acid sequences included the ones which do not belong to CCCH family. For instance, the sequences "CIRLEAYEDCIEYFVCLNGH" in AT5G15440 could be detected as a C-X8-C-X5-C-X3-H motif. Our comparison analysis showed that the entries obtained from Perl program contained all the entries obtained from BLAST searches. Thus, the data got from Perl program were used for further analysis.

Subsequently, we performed SMART and Pfam analysis to find the putative CCCH protein sequences [39, 40]. Seventy-nine proteins with typical CCCH motif encoded by 68 genes were detected from the original data, and five of them, AT1G19860, AT2G16470, AT2G33835, AT3G51120 and AT5G12440, contained a predicted CCCH motif with low confidence values. The proteins which were not detected by SMART and Pfam were eliminated from our study. For AT3G51180, which could not be detected by SMART and Pfam but has a sister gene in genome duplication region and contains a C-X8-C-X5-C-X3-H motif, was also considered as a member of CCCH zinc-finger family. In addition, we wrote another program to analyse the type of CCCH motifs based on the residue number between the four conserved amino acids (see Additional file 4 and 5). If a new motif was not referred in the 18 groups, we modified our Perl program with the new motif to search against the entire proteome again. For instance, when a detected C-X8-C-X6-C-X3-H motif was not mentioned by Berg et al., we rewrote new regular expression for matching the sequences of this motif. After verifying the reliability of the new motif in the protein by SMART and Pfam, we added it to the CCCH family. The motifs which were not mentioned by Berg et al. were also listed in Table 1, making a total of 23 groups.

Community microarray data, Massively Parallel Signature Sequencing searches, and ESTs database provided confirmative transcript information for all 68 genes (Table 2). From the above evidences, no pseudogene was found among the 68 genes. The CCCH genes were listed in Table 2 and their entry numbers were assigned according to the order of CCCH genes appearing on the Arabidopsis chromosomes, from the short arm to the long arm, and from chromosomes 1 to 5 (see Additional file 6).

Table 2 The CCCH gene family of Arabidopsis.

Phylogenetic analysis of the CCCH proteins in Arabidopsis

To evaluate the evolutionary relationships within the CCCH gene family in Arabidopsis, we performed a phylogenetic analysis of the 68 Arabidopsis protein sequences to construct a phylogenetic tree (Figure 3). Although the signature CCCH motif of the CCCH proteins is well conserved, the motif number of each protein and spacing amino acids between adjacent CCCH zinc-finger motifs are diverse, so it is difficult to make alignment using CCCH motifs. Accordingly, an alignment of full-length CCCH proteins was constructed using ClustalX and refined manually [41]. A phylogenetic tree was generated with the neighbor-joining method [42]. For statistical reliability, we conducted bootstrap analysis with 1000 replicates. The tree topologies are similar with different tree-building methods except in deep nodes. From the values obtained in the bootstrap analysis, it is apparent that the phylogenetic relationship is unclear and the bootstrap values are low in deep nodes. Nevertheless, in outer clades, the CCCH protein has better resolution, permitting subfamilies of proteins to be delimited. This phenomenon was also observed in the analysis of the bHLH [43], Dof [8], WRKY [44], and IQD [45] transcription factor families. Alternatively, it may be due to the divergence of the CCCH motif, especially the diverse motif numbers and spacing in the amino acid sequences between Cys and His in each protein. We could not infer evolutionary relationships between the different subfamilies of CCCH proteins because the internal nodes did not show high support. By contrast, within each subfamily, the strong amino acid sequence conservation is evident from the short branch lengths at the tips of the tress, suggestive of strong evolutionary relationships among subfamily members.

Figure 3
figure 3

An analytical view of the Arabidopsis CCCH gene family. The following parts are shown from left to right. Protein neighbor-joining tree: The unrooted tree, constructed using ClustalX (1.83), summarizes the evolutionary relationship among the 68 members of CCCH families. The neighbor-joining tree was constructed using aligned full-length amino acid sequences. The proteins are named according to their gene name (see Table 2) with the CCCH zinc finger number of each protein. The tree shows the 11 major phylogenetic subfamilies (left column, numbered I to XI and marked with different alternating tones of a gray background to make subfamily identification easier) with high predictive value. The numbers beside the branches represent bootstrap values (≥500) based on 1000 replications that were used to class the major 11 subfamilies. Gene structure: The gene structure is presented by black exon(s) and spaces between the black boxes correspond to introns. The sizes of exons and introns can be estimated using the horizontal lines. Protein structure: Each black box represents the motif in the protein, as indicated in the table on the left side. The conserved motifs outside CCCH motif are highlighted with an arranged number, and the same number referred to the same motif. The length of the motif can be estimated using the scale at top. aa, amino acids.

Based on the statistical support of each branch, we divided the Arabidopsis CCCH family into 11 subfamilies, which was supported by more than 90% bootstrap value, designated I to XI (Figure 3). The genes with low bootstrap value were not divided into subfamilies, and also were not taken into consideration for further analysis. The genes clustered pairwisely and supported by bootstrap value 1000 are mostly paralogous genes in Arabidopsis. Totally, 18 pairs of putative paralogous genes were found in Arabidopsis.

The analysis of protein motif and gene structure

To discover motifs shared among related proteins within the family, we used MEME [46], which performs motif searches in groups of related DNA or protein sequences. The search was performed separately for each of the subfamilies of proteins. Moreover, SMART and Pfam were used to annotate the motifs found by MEME. As expected, the proteins within same subfamily exhibit the common CCCH motifs, suggesting the major functional role of CCCH motif in these proteins, although few motifs outside the CCCH domain were detected by MEME as well. The schemes of the protein motifs of the individual members of the CCCH family clearly demonstrate structural similarities among the proteins within one subfamily (Figure 3). For example, the members of subfamily IV contain tandem WD40 domains. All eleven members of subfamily IX consist of three highly conserved tandem motifs (C-X5-H-X4-C-X3-H, C-X7–8-C-X5-C-X3-H and C-X5-C-X4-C-X3-H) and nine members of subfamily I contain five C-X8-C-X5-C-X3-H zinc finger motifs. Furthermore, twelve proteins contain well-defined RNA-binding domain RRM or KH motifs suggesting their potential roles involved in RNA binding [47].

The gene structures in terms of intron number and gene length were also consistent with the phylogenetic subfamilies defined in Figure 3. Most members in same subfamilies had similar intron/exon structure. The fact that they not only contain common motifs but also share very similar intron/exon structure supports their close evolutionary relationship and membership in the same subfamily. Taken together, our results validate the classification of the CCCH family genes of Arabidopsis.

Evolution and divergence of the CCCH family genes in Arabidopsis

The Arabidopsis genome has undergone several rounds of genome-wide duplication events, including polyploidy [48, 49], which has great impact on the amplification of members of a gene family in the genome. To further investigate the relationship between the genetic divergence within the CCCH family and gene duplication in Arabidopsis, the chromosomal location of each CCCH gene was determined from the genomic sequences of Arabidopsis [50]. Based on the chromosomal location information provided by the NCBI [51] and TAIR [52], we localized 68 CCCH genes in Arabidopsis chromosomes and determined that the genes are distributed across all five chromosomes (Figure 4). Relatively low densities of CCCH genes were observed in some chromosomal regions, including the top of chromosome 4, and the bottom of chromosome 1. Twenty-two genes were found in previously identified duplicated segmental regions on chromosomes that are the result of a polyploidy that occurred around 24 to 40 million years ago, probably close to the emergence of the crucifer family [53]. Another seven gene pairs (AtC3H1 and AtC3H51, AtC3H8 and AtC3H60, AtC3H12 and AtC3H28, AtC3H14 and AtC3H15, AtC3H30 and AtC3H56, AtC3H46 and AtC3H55, AtC3H59 and AtC3H62) were identified to share common gene structure and zinc finger motif. Although Blanc et al. did not include the seven gene pairs as duplicated pair genes in recently duplicated segmental chromosomes, the phylogenetic relationship and sequence similarity suggest that they are closely related to each other. Furthermore, AtC3H12 and AtC3H28 were listed as duplicated genes in the segmental duplications dataset maintained from TIGR. Therefore, they were considered to be putative duplicated genes (Figure 4). Consequently, about 53% of CCCH genes, which lie within recently duplicated segmental chromosomes, have a clear relative in these regions. Since the density of the duplicated genes in recently duplicated segmental chromosomes was reported to be 28.0% ± 7.8% [53], the duplicated pairs of CCCH genes have been preferentially retained compared with other genes. This finding is consistent with a previous study demonstrating that duplicated genes involved in signal transduction and transcription are preferentially retained [54, 55]. In summary, large-scale segmental duplication events appear to have exclusively contributed to the current complexes of the CCCH gene family.

Figure 4
figure 4

Chromosomal distribution and segmental duplication events for Arabidopsis CCCH genes. The chromosome number is indicated at the top of each chromosome. The green boxes indicate the duplicated segmental regions resulting from the most recent polyploidy. Only the duplicated regions containing CCCH genes are shown. Blue lines connect corresponding sister gene pairs in duplicated blocks. AtC3H1 and AtC3H51, AtC3H8 and AtC3H60, AtC3H12 and AtC3H28, AtC3H14 and AtC3H15, AtC3H30 and AtC3H56, AtC3H46 and AtC3H55, AtC3H59 and AtC3H62 are potential duplicated gene pairs which are marked with the same color rectangle, as described in the text.

In order to test whether CCCH family is ubiquitous in other organisms, CCCH genes in moss were also identified. Although ESTs are incomplete, a preliminary examination using the BLAST program suggested that almost all of the CCCH genes belonging to subfamily I, II, III, IV, V, VI, VII, X and XI in Arabidopsis have corresponding homologous members in moss. For subfamily VIII, since only 3 members are present in Arabidopsis, it is difficult to detect their corresponding homologous genes in the incomplete EST database of moss. However, the subfamily IX, which contains 11 genes in Arabidopsis, is not extracted in moss, suggesting that it is a particular subfamily only belonging to advanced plants or arises from the evolutionary process to adapt the circumstance.

To explore the evolutionary history of the CCCH gene family in greater detail and to testify our assumption about subfamily IX, we searched publicly available genomic and EST databases for all angiosperm species represented in the TIGR Plant Gene Indices as well as for the gymnosperm Pinus ssp. [56]. As expected, almost all of the members in Arabidopsis have their counterparters in these higher plants. These results suggest that the basis of the phylogenetic topology of the CCCH family had already been established before the divergence of vascular plants and the subfamily IX may be a special subfamily belonging only to higher plant.

Identification of CCCH genes in rice

Rice is one of the most important food crops in the world and it has been used as a major model species in plant (especially monocot) functional genomics research because of its relatively small genome size and extensively sequenced genome [32]. To explore the occurrence and size of the CCCH gene family in rice, we wrote a Perl program to analyze the genome file of rice (ssp. japonica) got from The TIGR Rice Genome Annotation Database and Resource (see Additional file 7). Multiple BLAST searches were also performed in several rice databases using the typical CCCH protein sequence as the queries to validate the primary results. After detecting the CCCH motifs by SMART and Pfam, 67 CCCH genes were identified in rice (Table 3, see Additional file 8). To determine their phylogenetic relationships, a phylogenetic tree was constructed using the full length protein sequences. Like the CCCH protein family in Arabidopsis, the phylogenetic analysis reveals 8 subfamilies supported by high bootstrap values. The members of same subfamily share the similar gene structure and zinc finger motifs (see Additional file 9A). We also detected two large subfamilies in the tree, in which the members share the same characteristic with the genes of corresponding subfamilies in Arabidopsis. Unlike the Arabidopsis CCCH family, the distribution of CCCH genes in the rice genome is clearly biased towards the chromosomes 1 to 7 (see Additional file 10). Four CCCH genes are each found on chromosomes 8, 9 and 12, and only one gene, however, is present on chromosome 10 and 11, respectively. The topology of the phylogenetic tree of rice CCCH genes indicates 13 pairs of putative duplicated genes that are close paralogs (OsC3H2 and OsC3H35, OsC3H5 and OsC3H6, OsC3H9 and OsC3H39, OsC3H10 and OsC3H37, OsC3H14 and OsC3H44, OsC3H19 and OsC3H32, OsC3H22 and OsC3H23, OsC3H24 and OsC3H67, OsC3H38 and OsC3H60, OsC3H40 and OsC3H49, OsC3H42 and OsC3H43, OsC3H54 and OsC3H51, OsC3H63 and OsC3H66). These genes represent 39% of the rice CCCH gene family members and might have evolved from putative rice genome duplication events.

Table 3 The CCCH gene family of rice.

Comparative phylogenetic analyses of the CCCH genes in Arabidopsis and rice

In order to evaluate the evolutionary relationship among the CCCH proteins, we performed a phylogenetic analysis based on the full-length amino acid sequences of Arabidopsis and rice. Because of the large number of taxa and relatively small number of characters in full-length sequences, the bootstrap values of internal nodes were low, whereas the outer nodes had more credible bootstrap values, allowing for clustering of the CCCH genes of Arabidopsis and rice into 20 subfamilies (see Additional file 9B). The tree topology, as well as the subfamily organization, resembled those from the rice and Arabidopsis individual trees. Eight subfamilies from Arabidopsis and rice clustered together into the same branch in the combined tree and 30 putative orthologs were identified in the tree (see Additional file 9B), suggesting that an ancestral set of CCCH genes already existed before the monocot-eudicot divergence. Moreover, the paralogous genes of Arabidopsis and rice were already displayed as paralogs in the respective trees.

Despite the differences in genome sizes between Arabidopsis and rice (125 Mb and 389 Mb respectively) and encoded number of genes, the two plant species appear to have a very similar number of genes encoding CCCH zinc-finger proteins (68 and 67 putative numbers, respectively). Compared with other gene families in Arabidopsis and rice, the CCCH gene family is one of the largest families in plant and includes diverse members with distinct specificities. As shown in Figure 5, the proteins containing C-X7–8-C-X5-C-X3-H motif constitute the largest groups in the CCCH zinc-finger family. In Arabidopsis, 42 members contain C-X7-C-X5-C-X3-H motif and 44 members contain C-X8-C-X5-C-X3-H motif, while the number is 35 and 36 in rice, respectively, suggesting that C-X7–8-C-X5-C-X3-H motif may be an ancestor of other CCCH motifs. In addition, five novel CCCH motifs beyond the 18 groups were identified, including C-X4-C-X5-C-X3-H, C-X5-C-X4-C-X3-H, C-X7-C-X6-C-X3-H, C-X8-C-X6-C-X3-H and C-X15-C-X5-C-X3-H. Further experiments are required to determine their binding activity to Zn2+ and biological functions.

Figure 5
figure 5

Schematic structures of Arabidopsis CCCH proteins. The figure is schematic protein structures of all CCCH-type zinc-finger proteins identified in Arabidopsis. The CCCH zinc fingers are shown by black boxes. Types of zinc fingers are indicated by A to W, which are cross-indexed in the table on the right side. Numbers of amino acids in the spacers are indicated on each spacer region for several-fingered proteins. The lengths of each protein are shown on the right of the schematic structures. Numbers of different types of CCCH zinc finger proteins are presented in the brackets.

To determine sequence features of these CCCH motifs, we performed sequence alignments of 302 CCCH motifs (152 from Arabidopsis and 150 from rice) using the Clustal X program [41]. The weblogo indicated that four (three cysteines and one histidine) amino acids are completely conserved among all the CCCH motifs and more than 75% of the CCCH motifs contain glycine and phenylalanine (Figure 6) [57]. In addition, the distributions of conserved amino acids among the CCCH motifs of both Arabidopsis and rice are similar. Given these characteristic of verified protein sequences, we developed a new criteria to objectively define those sequences to be considered as CCCH motif. The CCCH proteins are characterized by one to six C-X4–15-C-X4–6-C-X3-H motifs which are glycine-rich and phenylalanine-rich sequences.

Figure 6
figure 6

Sequence logos for the CCCH zinc finger motifs of Arabidopsis and rice. Numbers on the x-axis represent the sequence positions in zinc finger motifs. The y-axis represents the information content measured in bits. The sequence logos were derived using WebLogo [4].

The expression pattern of CCCH genes in Arabidopsis and rice

Since gene expression patterns can provide important clues for gene function, we examined the expression of Arabidopsis and rice CCCH genes in root, leaf, flower and seed tissues using massively parallel signature sequencing (MPSS) data and EST data of NCBI (Figure 7) [58, 59]. In Arabidopsis, expression profiles of 65 and 63 CCCH genes were extracted from MPSS and EST databases, respectively. After integrating two data together, we found that most of the genes have a very broad expression spectrum, and only six genes (AtC3H9, AtC3H10, AtC3H28, AtC3H51, AtC3H53 and AtC3H62) were not detected from any specific tissue according to EST and MPSS data. Except AtC3H53, five of these six genes were detected in mix tissue of Arabidopsis from EST database. Furthermore, we also summarized the expression of rice CCCH genes using the same methods. Expression of 57 rice CCCH genes was detected from MPSS database, and 58 had matching ESTs. Thirteen CCCH genes in rice had no tissue expression pattern, while eight of them contain corresponding ESTs and five (OsC3H21, OsC3H48, OsC3H57, OsC3H58 and OsC3H64) without expression information might be pseudogenes, or expressed at specific developmental stages or under special conditions.

Figure 7
figure 7

The expression patterns for Arabidopsis and rice CCCH genes from MPSS and EST data. The expression patterns for Arabidopsis (the part of left side) and rice (the part of right side) CCCH genes are shown. The letter R above the column of expression data refers to root, I refers to inflorescence, L refers to leaf, and S refers to seed (silique). A positive signal is indicated by a colored box for the following tissues: grey for roots (R), red for inflorescences (I), green for rosette leaves (L), and yellow for siliques (S). The white box indicates that no expression could be detected. The number on the left indicates the subfamily and the black points on the right show the origin of expression data for each gene. CCCH proteins are aligned in the same order as they appear in the phylogenetic trees. Subfamilies of CCCH proteins are highlighted by vertical bars next to the gene identifiers. The expression profile of AtC3H37 (HUA1) was got from RNA filter hybridization [19].

By combining the EST and MPSS results together, we found that the accumulation of CCCH gene transcripts not only is associated with different tissues, but also the expression pattern of each CCCH gene member differed. According to expression profiles, CCCH genes can be classified into three groups. The largest group is the genes that expressed in all tissues, including 33 genes in Arabidopsis and 36 genes in rice, and the second group contains 18 Arabidopsis CCCH genes and 12 rice CCCH genes whose expression were detected in the majority, but not all, of the tissues. Interestingly, the third group includes 17 genes with very specific expression. AtC3H50 and AtC3H68 are exclusively expressed in roots, while AtC3H18, AtC3H35, AtC3H43, OsC3H36 and OsC3H39 in inflorescences, AtC3H45, OsC3H28 and OsC3H30 in leaves, AtC3H2, AtC3H5, AtC3H13, AtC3H15, AtC3H54, OsC3H1 and OsC3H52 in seeds. When compared the expression pattern of the AtC3H45, which could not be detected by SMART and Pfam but contains a C-X8-C-X5-C-X3-H motif, with its sister gene AtC3H68, we found that they share different expression profiles with relatively low level (AtC3H45 could be only detected in leaves, while AtC3H68 is exclusively expressed in roots). In addition, we examined the expression of other CCCH duplicated gene pairs of both Arabidopsis and rice, and only 20 of 31 pairs (18 duplicated gene pairs in Arabidopsis and 13 in rice) share the same expression pattern. These results are consistent with the previous research by Blanc and Wolfe, that the expression profiles of the two paralogs have diverged in concert, forming two parallel networks, and the expression of each gene is strongly correlated with the other nonhomologous genes in its network but poorly correlated with its paralog in the other network, suggesting functional diversification of the surviving duplicated genes is a major feature of the long-term evolution of polyploids [55].

Expression profiling within different tissues is only the first step to understanding the function of CCCH genes. Overall, the CCCH family members show diverse expression patterns. A majority of the CCCH genes are expressed in all tissues tested, with various expression levels. It might be a common character of large transcription factor families, such as MYB family [60]. In addition, the genes with specific expression patterns can be the focus of functional studies for their possible roles in specific tissues. Because of the limitation of current information from MPSS and EST databases, expression of some genes within the tissues might not be detected. Further investigation of spatial and developmental expression patterns will need to be performed to obtain more detailed expression information.

Characteristics of subfamily IX genes of Arabidopsis and subfamily I genes of rice

Subfamily IX containing 11 members in Arabidopsis and subfamily I containing 9 members in rice are the largest CCCH subfamily in each species, respectively. All the proteins encoded by these 20 genes include two CCCH motifs. With the exception of OsC3H10, OsC3H37 and OsC3H52 in rice, the products of 17 genes commonly contain two tandem consensus motifs, C-X7–8-C-X5-C-X3-H and C-X5-C-X4-C-X3-H. In OsC3H10 and OsC3H37 proteins, the first zinc finger is replaced by a C-X10-C-X5-C-X3-H motif, while in OsC3H52 the first zinc finger is C-X15-C-X5-C-X3-H. Amino acid sequence analysis revealed that all twenty members contain highly conserved C-X5-C-X4-C-X3-H motif, which is a characteristic of these proteins. Additionally, putative consensus CHCH (C-X5-H-X4-C-X3-H) motif was also observed in these 20 proteins, implying that it may be a novel zinc finger motif and execute some biological functions (Figure 8A). Moreover, the number of amino acids between the three motifs is invariable and all the 20 genes have no introns (Figure 3, see Additional file 9A).

As shown in Figure 9A, phylogenetic analysis indicated that the subfamily IX of Arabidopsis consists of two major subgroups. Interestingly, with exception of two CCCH motifs, ankyrin (ANK) repeat motif was also identified by MEME search using amino acid sequences of all members of subgroup 2 (Figure 8B). The ANK repeat motif is one of the most common protein-protein interaction motifs in nature, and it has been found in proteins of diverse function such as transcriptional initiators, ion transporters and signal transducers [61–64]. The database search further revealed that only six genes coding proteins with ANK repeats and zinc finger domains in whole Arabidopsis genome, five of them belong to CCCH subfamily IX. In addition, to evaluate the evolutionary relationships within the 20 genes, we performed a combined phylogenetic analysis of the 20 Arabidopsis and rice amino acid sequences to obtain a joint tree. Three subgroups can be observed from the combined tree (Figure 9B). In subgroup 1 of the combined tree, all nine genes of both Arabidopsis and rice contain the ANK domain. Surprisingly, within each subgroup, the rice and Arabidopsis genes appear closely related to genes of the same species, and only one putative ortholog pairs (AtC3H66 and OsC3H33) were identified in the combined tree. One possible explanation of this phenomenon is that these genes may go through a complex and extensive evolution in each species after monocot-eudicot divergence.

Figure 8
figure 8

Amino acid sequence analysis of the members in subfamily IX of Arabidopsis. A, Multiple sequence alignment of the three zinc finger motifs of the members in the subfamily IX. The black and dashed bars represent CCCH motifs and putative CHCH zinc finger, respectively. Black and gray shading indicate identical and conserved amino acid residues present in more than 50% of the aligned sequences, respectively. B, Amino acid sequence alignment of ankyrin repeats in CCCH family. C, Amino acid sequence alignment of putative NES sequences in subfamily IX.

Figure 9
figure 9

Phylogenetic trees of genes in subfamily IX of Arabidopsis and subfamily I of rice. A, Phylogenetic tree of Arabidopsis CCCH genes in subfamily IX. The unrooted tree was inferred by MEGA 3.1 and the neighbor-joining method after the alignment of the full-length amino acid sequences of the 11 Arabidopsis genes in subfamily IX. The numbers beside the branches represent bootstrap values based on 1000 replications. Subgroups of CCCH proteins are highlighted by vertical bars next to the gene identifiers. The scale bar corresponds to 0.1 estimated amino acid substitutions per site. B, Joined phylogenetic tree of rice subfamily I and Arabidopsis subfamily IX CCCH genes. The unrooted tree was constructed using MEGA 3.1 and the neighbor-joining method after the alignment of the full-length amino acid sequences of 20 Arabidopsis and rice genes. Numbers on branches indicate the percentage of 1000 bootstrap replicates that support the adjacent node. Black braces and numerals at right indicate the three subgroups. The scale bar corresponds to 0.05 estimated amino acid substitutions per site.

Previous studies demonstrated that OsDOS, AtCPSF30 and HUA1 are nuclear-localized proteins [19, 20, 30]. The database research with full-length protein sequences of subfamily IX of Arabidopsis and subfamily I of rice suggested that all members are localized in the cell nucleus (see Methods). In animal, TTP was shown to export from the nucleus mediated by a Leucine-rich Nuclear Export Signal (NES). To our knowledge, no Leucine-rich NES in plant proteins has been identified. When we performed our program using the developed widely accepted NES consensus [LV]-x(2,3)-[LIVFM]-x(2,3)-L-x-[LIMTKD] to detect 68 Arabidopsis CCCH proteins [65], 54 proteins containing putative NES sequences were identified including all of 11 members of subfamily IX (Figure 8C, see Additional file 11 and 12). The result suggests that all subfamily IX proteins of Arabidopsis may be nucleocytoplasmic shuttle proteins involved in signal transduction events [66].

Recently, a few genes within these 20 genes have been shown to play crucial roles in abiotic and/or biotic stress-responsive gene expression. OsDOS (LOC_Os01g09620) was proved to be involved in JA pathway [30]. ZFAR1 (At2g40140) transcript was induced by Botrytis in inoculated Arabidopsis lines and zfar1 mutant was hypersensitive to ABA [67]. In order to better understand the function of these genes, we firstly examined the expression of all the genes of subfamily IX of Arabidopsis in response to multiple environmental stimuli by means of microarray data available at Genevestigator site [68]. The results revealed that all genes of subfamily IX were activated or suppressed by various stresses, including salt, cold, mannitol, ABA, hypoxia and osmotic stress (Figure 10A). Secondly, we analyzed the expression of these 11 genes in Arabidopsis plants treated with drought (mannitol), salt, cold, and ABA by RT-PCR. As shown in Figure 10B, the RT-PCR results are in better agreement with the microarray data, suggesting their involvement in abiotic and/or biotic stress.

Figure 10
figure 10

The expression patterns of subfamily IX genes in response to abiotic stress and ABA treatments. A, Microarray data for the CCCH genes of subfamily IX of Arabidopsis were extracted using Genevestigator. Up- and down-expression patterns of eleven genes under six different treatments are indicated by black and grey arrowhead, respectively. B, Expression pattern for genes in subfamily IX of Arabidopsis following stress and ABA treatments detected by RT-PCR. Plants were grown on Murashige and Skoog (1962) (MS) agar medium for 3 weeks and were treated with NaCl (300 mM), mannitol (300 mM), cold (4°C), ABA (100 μM) or water (as a control). The EF1-α gene was used as an internal control.

RNA-binding roles for CCCH proteins

Recently, the structure of the CCCH domain from the TIS11D in complex with RNA monomer of AU-rich element (ARE) have been determined [69]. Sequences alignment show that among 135 CCCH proteins in Arabidopsis and rice, AtC3H14, AtC3H15, OsC3H9 and OsC3H39 share high amino acid sequence identities with TIS11D (Figure 11A). Like the TIS11D, each of these four proteins contains two tandem (C-X8-C-X5-C-X3-H) CCCH zinc fingers and a linker of 18 residues between the two zinc finger motifs. In addition, many other residues are strictly conserved. To obtain insights of the structure of plant CCCH zinc finger, we have taken the coordinates of the TIS11D CCCH domain complex and modelled the structure of the interaction of the AtC3H14 (AT1G66810) CCCH domain peptide with the RNA nonamer (5'-UUAUUUAUU-3'), using the Swiss-Model programs [70]. It is apparent that the RNA-binding domain of the AtC3H14 is likely to be identical in structure to that of TIS11D. Compared with the structure of TIS11D, the characteristic KTEL(V) motif at the N-terminus of each zinc finger domain provides a critical part of the RNA-binding surface; each motif forms two walls of a deep pocket that accommodates the bases of U6 and U2 in fingers 1 and 2, respectively (Figure 11B). The model illustrates the aromatic stack formed from U2-Phe293-A3, U4-Tyr287-U5, U6-Phe255-A7 and U8-Tyr249-U9, which are essential for high-affinity binding. Mutation within these critical hydrophobic amino acids might make the domain abrogate RNA binding (Figure 11C). Further mutagenesis experiments would be required to establish the contributions of these residues to binding affinity.

Figure 11
figure 11

Sequence alignments and the structure of AtC3H14 tandem domains in complex with the RNA. A, Alignment of AtC3H14, AtC3H15, OsC3H9, OsC3H39, hTTP and hTIS11d. Conserved motifs are highlighted in light blue, gold and red (cysteines and histidines involved in zinc binding, respectively), and green (aromatic residues involved in base-stacking interactions with RNA). h, Homo sapiens. B, Structure of the RNA complex of AtC3H14. This proposed structure was modelled on the original nuclear magnetic resonance structure describe by Hudson et al. [69], using their pdb coordinates and the Swiss-Model program. The RNA oligonucleotide is shown in magenta (A3 and A7 are shown in yellow) with the 5' and 3' ends indicated. The protein backbone are shown to represent the structure, with alpha helices are coloured royal blue, beta sheets are coloured yellow, the zinc-coordinating ligands are colored green and the zinc residues are coloured red. C, Molecular contact surface of AtC3H14 structure showing surface topology. Blue denotes convex surfaces. The location of the motifs that form the U6 and U2 binding pockets and the aromatic residues which are essential for high-affinity binding are indicated.

In mammals, some proteins containing tandem CCCH domains, such as TTP and TIS11D, have been proved to bind to the class II ARE in the 3'-untranslated region (3'-UTR) of target mRNAs and promote their deadenylation and degradation [71]. AREs are sequence elements of 50–150 nt that are rich in adenosine and uridine bases. They are located in the 3'-UTRs of many but not all mRNAs that have a short half-life and have been identified by their capacity to provoke degradation of the host mRNA by a mechanism dependent on deadenylation. It has been estimated that 5–8% of human genes code for ARE-containing mRNAs, and the corresponding proteins perform a variety of functions implicated in numerous transient biological processes [72]. Henceforth, only two plant CCCH (HUA1 and AtCPSF30) proteins have been proved to have the RNA-binding ability by experimental evidence. Extraction of ARE-containing mRNAs in 3'-UTR of Arabidopsis genome was performed using Patmatch tools of TAIR. The 3'-UTRs were searched for the 13-bp pattern WWWU(AUUUA)UUUW. Totally, we identified 200 genes containing 258 ARE in their 3'-UTR mRNA (see Additional file 13). The potential relationship between the plant CCCH proteins and the AU-rich mRNA need been confirmed by experiment approaches in the future.

Conclusion

In this study, we developed a new method employed regular expression of Perl to identify the gene family. By this method, 68 and 67 CCCH genes were identified in Arabidopsis and rice, respectively, and they formed 20 subfamilies that were supported by phylogeny, protein motifs and intron/exon structures. Compared with other gene families in rice and Arabidopsis, the CCCH gene family is one of the largest families in plants. Gene duplication analysis revealed that during the expansion of the CCCH gene family, many subfamilies and subgroups have evolved, resulting in a high level of functional divergence. Most of the subfamilies/subgroups are present both in Arabidopsis and rice, suggesting that the appearance of many of the genes in these species predates monocot/eudicot divergence. Members within a given subfamily/subgroup may have recent common evolutionary origins and may possess common zinc finger motifs that have related molecular functions. Available data supported the hypothesis that CCCH genes in plants perform a variety of functions in different tissues at multiple developmental stages. The subfamily IX of Arabidopsis with many interesting characters was investigated in details, in which members are shown to be involved in plant stress responses. We then proposed that the genes of subfamily IX possibly play roles as nucleocytoplasmic shuttle proteins involved in signal transduction events. Although the function of most plant CCCH genes is unknown, the phylogenetic and expression analyses provide a solid foundation for future functional studies in both Arabidopsis and rice. Overall, we have identified a novel zinc finger family and results presented here can serve as useful information for guiding future experimental work and understanding the structure-function relationship of the members of the CCCH gene family.

Methods

Identification of CCCH protein in Arabidopsis and rice

A preliminary search for CCCH proteins was performed using a program written by Perl. This program could search against the CCCH proteins with the regular expression in entire Arabidopsis and rice genome file which was downloaded from the FTP of TAIR [52] and TIGR [73]. The regular expression (C\w{6,14}C\w{4,5}C\w{3}H) was designed according to the previous studies. All of the proteins which matched the expression would be considered as candidate CCCH proteins. The obtained protein sequences were then examined for the CCCH motif using the hidden Markov model of SMART/Pfam tool. The proteins without CCCH motif were eliminated from the datasets. The new type motif which was not referred by the previous studies was used to modify the regular expression. The modified regular expression was used to retrieve for another round of searching, and new results were added to the original dataset.

Another approach of retrieving CCCH genes from database was employed. Multiple database searches were performed using the Basic Local Alignment Search Tool algorithms BLASTP and TBLASTN with various published CCCH proteins as query sequences and with the E-value cutoff set as 1e-005. For Arabidopsis, the databases searched included NCBI [51], DATF [74], MAtDB [75], TAIR [52], and TIGR [73] and for rice, the databases searched included NCBI [51], TIGR [76], Rice Genome Database-japonica of the Rice Genome Research Program [77] and The International Rice Genome Sequencing Project (IRSGP) [78]. The redundant sequences with different identification numbers and the same chromosome locus were removed from our data set.

The hits obtained from all the above methods were pooled together and another program written by Perl was performed to eliminate the redundant dataset. Some programs used in this research are listed in additional file (see Additional file 1, 2, 3, 4, 5 and 7).

Sequence properties of Arabidopsis and rice CCCH genes

The amino acid sequences of the CCCH proteins were analyzed for physicochemical parameters (ProtParam) and predicted subcellular localization (SubLoc v1.0) on DBSubLoc [79]. MEME (Multiple Expectation Maximization for Motif Elicitation) was used to identify conserved motif structures of CCCH protein sequences.

Alignment and phylogenetic analysis of CCCH sequences

Multiple alignments of amino acid sequences were performed using ClustalX and were manually corrected. For generating the phylogenetic tree, we used ClustalX (1.83) and the neighbor-joining algorithm. Bootstrap analysis with 1,000 replicates was used to evaluate the significance of the nodes. Representations of the calculated trees were constructed using TreeView. The phylogenetic trees of Figure 9 were constructed by neighbor-joining algorithms of MEGA3.1. Bootstrapping was performed 1000 times to obtain support values for each branch.

The location of CCCH genes on chromosomes

To determine the location of CCCH genes on Arabidopsis chromosomes, Chromosome Map Tool at TAIR was used [50]. Gene duplications and their presence on duplicated chromosomal segments were investigated using "Paralogous in Arabidopsis" with the default parameters set to a minimum threshold for paired proteins per block above 7 [80]. For rice, all the sequenced contigs of japonica cv Nipponbare have been physically constructed as pseudomolecules by the IRGSP, representing the 12 rice chromosomes, and available in GenBank. Each of the rice CCCH genes was positioned on these rice chromosome pseudomolecules by the BLASTN search.

Other plant species

To identify members of the moss CCCH protein family, multiple database searches were performed using the stand-alone BLAST tools available on NCBI. The EST database of moss was obtained on the Physcomitrella EST Project Web site [81], and then the file of results was parsed by a program written by Perl. For CCCH members of pine, the same method was used to search against the database got from NSF Genomics of Loblolly Pine Embryogenesis Project [56].

Expression analysis of CCCH genes

We used MPSS and EST data to detect the expression patterns of CCCH genes. The locus name of CCCH genes were used to query the MPSS database containing the signature information of the CCCH genes. EST data came from UniGene of NCBI, TIGR and TAIR. We also searched the expression data of Arabidopsis CCCH genes in the database of Genevestigator [68].

RT-PCR analysis

Plant tissues of Arabidopsis were harvested and ground in liquid nitrogen. For reverse transcriptase-mediated PCR analysis, total RNA was isolated with the RNeasy mini kit (Qiagen, USA) according to the manufacturer's instructions. The RNA preparation was then treated with DNase I. First strand synthesis of cDNA was performed by using oligo (dT) primer and M-MLV RT (Promega). PCR products were fractionated on 1% agarose gels containing ethidium bromide and photographed under UV light. These experiments were independently replicated at least three times under identical conditions. Details of primers are listed in the Additional file 14.

Structure model

The model of complex of plant CCCH protein with the RNA nonamer was constructed using the Swiss-Model programs [70]. The final merged files was presented and annotated with Swiss Pdb Viewer 3.7 and Rastop 2.2.

References

  1. Takatsuji H: Zinc-finger transcription factors in plants. Cell Mol Life Sci. 1998, 54 (6): 582-596. 10.1007/s000180050186.

    Article  PubMed  CAS  Google Scholar 

  2. Moore M, Ullman C: Recent developments in the engineering of zinc finger proteins. Brief Funct Genomic Proteomic. 2003, 1 (4): 342-355. 10.1093/bfgp/1.4.342.

    Article  PubMed  CAS  Google Scholar 

  3. Freemont PS: The RING finger. A novel protein sequence motif related to the zinc finger. Ann N Y Acad Sci. 1993, 684: 174-192. 10.1111/j.1749-6632.1993.tb32280.x.

    Article  PubMed  CAS  Google Scholar 

  4. Arnaud D, Dejardin A, Leple JC, Lesage-Descauses MC, Pilate G: Genome-wide analysis of LIM gene family in Populus trichocarpa, Arabidopsis thaliana, and Oryza sativa. DNA Res. 2007, 14 (3): 103-116. 10.1093/dnares/dsm013.

    Article  PubMed  CAS  Google Scholar 

  5. Kosarev P, Mayer KF, Hardtke CS: Evaluation and classification of RING-finger domains encoded by the Arabidopsis genome. Genome Biol. 2002, 3 (4): RESEARCH0016-10.1186/gb-2002-3-4-research0016.

    Article  PubMed  Google Scholar 

  6. Nakano T, Suzuki K, Fujimura T, Shinshi H: Genome-wide analysis of the ERF gene family in Arabidopsis and rice. Plant Physiol. 2006, 140 (2): 411-432. 10.1104/pp.105.073783.

    Article  PubMed  CAS  Google Scholar 

  7. Zhang Y, Wang L: The WRKY transcription factor superfamily: its origin in eukaryotes and expansion in plants. BMC Evol Biol. 2005, 5 (1): 1-10.1186/1471-2148-5-1.

    Article  PubMed  CAS  Google Scholar 

  8. Lijavetzky D, Carbonero P, Vicente-Carbajosa J: Genome-wide comparative phylogenetic analysis of the rice and Arabidopsis Dof gene families. BMC Evol Biol. 2003, 3: 17-10.1186/1471-2148-3-17.

    Article  PubMed  Google Scholar 

  9. De J, Lai WS, Thorn JM, Goldsworthy SM, Liu X, Blackwell TK, Blackshear PJ: Identification of four CCCH zinc finger proteins in Xenopus, including a novel vertebrate protein with four zinc fingers and severely restricted expression. Gene. 1999, 228 (1–2): 133-145. 10.1016/S0378-1119(98)00617-9.

    Article  PubMed  CAS  Google Scholar 

  10. DuBois RN, McLane MW, Ryder K, Lau LF, Nathans D: A growth factor-inducible nuclear protein with a novel cysteine/histidine repetitive sequence. J Biol Chem. 1990, 265 (31): 19185-19191.

    PubMed  CAS  Google Scholar 

  11. Gomperts M, Pascall JC, Brown KD: The nucleotide sequence of a cDNA encoding an EGF-inducible gene indicates the existence of a new family of mitogen-induced genes. Oncogene. 1990, 5 (7): 1081-1083.

    PubMed  CAS  Google Scholar 

  12. Carrick DM, Lai WS, Blackshear PJ: The tandem CCCH zinc finger protein tristetraprolin and its relevance to cytokine mRNA turnover and arthritis. Arthritis Res Ther. 2004, 6 (6): 248-264. 10.1186/ar1441.

    Article  PubMed  CAS  Google Scholar 

  13. Ma Q, Wadleigh D, Chi T, Herschman H: The Drosophila TIS11 homologue encodes a developmentally controlled gene. Oncogene. 1994, 9 (11): 3329-3334.

    PubMed  CAS  Google Scholar 

  14. Mello CC, Schubert C, Draper B, Zhang W, Lobel R, Priess JR: The PIE-1 protein and germline specification in C. elegans embryos. Nature. 1996, 382 (6593): 710-712. 10.1038/382710a0.

    Article  PubMed  CAS  Google Scholar 

  15. Nie XF, Maclean KN, Kumar V, McKay IA, Bustin SA: ERF-2, the human homologue of the murine Tis11d early response gene. Gene. 1995, 152 (2): 285-286. 10.1016/0378-1119(94)00696-P.

    Article  PubMed  CAS  Google Scholar 

  16. Seydoux G, Mello CC, Pettitt J, Wood WB, Priess JR, Fire A: Repression of gene expression in the embryonic germ lineage of C. elegans. Nature. 1996, 382 (6593): 713-716. 10.1038/382713a0.

    Article  PubMed  CAS  Google Scholar 

  17. Taylor GA, Lai WS, Oakey RJ, Seldin MF, Shows TB, Eddy RL, Blackshear PJ: The human TTP protein: sequence, alignment with related proteins, and chromosomal localization of the mouse and human genes. Nucleic Acids Res. 1991, 19 (12): 3454-10.1093/nar/19.12.3454.

    Article  PubMed  CAS  Google Scholar 

  18. Thompson MJ, Lai WS, Taylor GA, Blackshear PJ: Cloning and characterization of two yeast genes encoding members of the CCCH class of zinc finger proteins: zinc finger-mediated impairment of cell growth. Gene. 1996, 174 (2): 225-233. 10.1016/0378-1119(96)00084-4.

    Article  PubMed  CAS  Google Scholar 

  19. Li J, Jia D, Chen X: HUA1, a regulator of stamen and carpel identities in Arabidopsis, codes for a nuclear RNA binding protein. Plant Cell. 2001, 13 (10): 2269-2281. 10.1105/tpc.13.10.2269.

    Article  PubMed  CAS  Google Scholar 

  20. Delaney KJ, Xu R, Zhang J, Li QQ, Yun KY, Falcone DL, Hunt AG: Calmodulin interacts with and regulates the RNA-binding activity of an Arabidopsis polyadenylation factor subunit. Plant Physiol. 2006, 140 (4): 1507-1521. 10.1104/pp.105.070672.

    Article  PubMed  CAS  Google Scholar 

  21. Berg JM, Shi Y: The galvanization of biology: a growing appreciation for the roles of zinc. Science. 1996, 271 (5252): 1081-1085. 10.1126/science.271.5252.1081.

    Article  PubMed  CAS  Google Scholar 

  22. Lai WS, Carballo E, Strum JR, Kennington EA, Phillips RS, Blackshear PJ: Evidence that tristetraprolin binds to AU-rich elements and promotes the deadenylation and destabilization of tumor necrosis factor alpha mRNA. Mol Cell Biol. 1999, 19 (6): 4311-4323.

    Article  PubMed  CAS  Google Scholar 

  23. Lai WS, Carballo E, Thorn JM, Kennington EA, Blackshear PJ: Interactions of CCCH zinc finger proteins with mRNA. Binding of tristetraprolin-related zinc finger proteins to Au-rich elements and destabilization of mRNA. J Biol Chem. 2000, 275 (23): 17827-17837. 10.1074/jbc.M001696200.

    Article  PubMed  CAS  Google Scholar 

  24. Lai WS, Kennington EA, Blackshear PJ: Tristetraprolin and its family members can promote the cell-free deadenylation of AU-rich element-containing mRNAs by poly(A) ribonuclease. Mol Cell Biol. 2003, 23 (11): 3798-3812. 10.1128/MCB.23.11.3798-3812.2003.

    Article  PubMed  CAS  Google Scholar 

  25. Lai WS, Stumpo DJ, Blackshear PJ: Rapid insulin-stimulated accumulation of an mRNA encoding a proline-rich protein. J Biol Chem. 1990, 265 (27): 16556-16563.

    PubMed  CAS  Google Scholar 

  26. Ramos SB, Stumpo DJ, Kennington EA, Phillips RS, Bock CB, Ribeiro-Neto F, Blackshear PJ: The CCCH tandem zinc-finger protein Zfp36l2 is crucial for female fertility and early embryonic development. Development. 2004, 131 (19): 4883-4893. 10.1242/dev.01336.

    Article  PubMed  CAS  Google Scholar 

  27. Tenenhaus C, Subramaniam K, Dunn MA, Seydoux G: PIE-1 is a bifunctional protein that regulates maternal and zygotic gene expression in the embryonic germ line of Caenorhabditis elegans. Genes Dev. 2001, 15 (8): 1031-1040. 10.1101/gad.876201.

    Article  PubMed  CAS  Google Scholar 

  28. Li Z, Thomas TL: PEI1, an embryo-specific zinc finger protein gene required for heart-stage embryo formation in Arabidopsis. Plant Cell. 1998, 10 (3): 383-398. 10.1105/tpc.10.3.383.

    PubMed  CAS  Google Scholar 

  29. Schmitz RJ, Hong L, Michaels S, Amasino RM: FRIGIDA-ESSENTIAL 1 interacts genetically with FRIGIDA and FRIGIDA-LIKE 1 to promote the winter-annual habit of Arabidopsis thaliana. Development. 2005, 132 (24): 5471-5478. 10.1242/dev.02170.

    Article  PubMed  CAS  Google Scholar 

  30. Kong Z, Li M, Yang W, Xu W, Xue Y: A novel nuclear-localized CCCH-type zinc finger protein, OsDOS, is involved in delaying leaf senescence in rice. Plant Physiol. 2006, 141 (4): 1376-1388. 10.1104/pp.106.082941.

    Article  PubMed  CAS  Google Scholar 

  31. Yu J, Hu S, Wang J, Wong GK, Li S, Liu B, Deng Y, Dai L, Zhou Y, Zhang X: A draft sequence of the rice genome (Oryza sativa L. ssp. indica). Science. 2002, 296 (5565): 79-92. 10.1126/science.1068037.

    Article  PubMed  CAS  Google Scholar 

  32. Goff SA, Ricke D, Lan TH, Presting G, Wang R, Dunn M, Glazebrook J, Sessions A, Oeller P, Varma H: A draft sequence of the rice genome (Oryza sativa L. ssp. japonica). Science. 2002, 296 (5565): 92-100. 10.1126/science.1068275.

    Article  PubMed  CAS  Google Scholar 

  33. Sasaki T, Matsumoto T, Yamamoto K, Sakata K, Baba T, Katayose Y, Wu J, Niimura Y, Cheng Z, Nagamura Y: The genome sequence and structure of rice chromosome 1. Nature. 2002, 420 (6913): 312-316. 10.1038/nature01184.

    Article  PubMed  CAS  Google Scholar 

  34. Feng Q, Zhang Y, Hao P, Wang S, Fu G, Huang Y, Li Y, Zhu J, Liu Y, Hu X: Sequence and analysis of rice chromosome 4. Nature. 2002, 420 (6913): 316-320. 10.1038/nature01183.

    Article  PubMed  CAS  Google Scholar 

  35. The Arabidopsis Genome Initiative: Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature. 2000, 408 (6814): 796-815. 10.1038/35048692.

    Article  Google Scholar 

  36. Jamison DC: Perl programming for biologists. 2003, Hoboken, N.J.: Wiley-Liss

    Chapter  Google Scholar 

  37. Patwardhan N, Siever E, Spainhour S: Perl in a nutshell. 2002, Sebastopol, CA: O'Reilly, 2

    Google Scholar 

  38. Brown MC: Perl programmer's reference. 1999, Berkeley, CA: Osborne/McGraw-Hill

    Google Scholar 

  39. Schultz J, Copley RR, Doerks T, Ponting CP, Bork P: SMART: a web-based tool for the study of genetically mobile domains. Nucleic Acids Res. 2000, 28 (1): 231-234. 10.1093/nar/28.1.231.

    Article  PubMed  CAS  Google Scholar 

  40. Sonnhammer EL, Eddy SR, Durbin R: Pfam: a comprehensive database of protein domain families based on seed alignments. Proteins. 1997, 28 (3): 405-420. 10.1002/(SICI)1097-0134(199707)28:3<405::AID-PROT10>3.0.CO;2-L.

    Article  PubMed  CAS  Google Scholar 

  41. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25 (24): 4876-4882. 10.1093/nar/25.24.4876.

    Article  PubMed  CAS  Google Scholar 

  42. Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4 (4): 406-425.

    PubMed  CAS  Google Scholar 

  43. Li X, Duan X, Jiang H, Sun Y, Tang Y, Yuan Z, Guo J, Liang W, Chen L, Yin J: Genome-wide analysis of basic/helix-loop-helix transcription factor family in rice and Arabidopsis. Plant Physiol. 2006, 141 (4): 1167-1184. 10.1104/pp.106.080580.

    Article  PubMed  CAS  Google Scholar 

  44. Wu KL, Guo ZJ, Wang HH, Li J: The WRKY family of transcription factors in rice and Arabidopsis and their origins. DNA Res. 2005, 12 (1): 9-26. 10.1093/dnares/12.1.9.

    Article  PubMed  CAS  Google Scholar 

  45. Abel S, Savchenko T, Levy M: Genome-wide comparative analysis of the IQD gene families in Arabidopsis thaliana and Oryza sativa. BMC Evol Biol. 2005, 5: 72-10.1186/1471-2148-5-72.

    Article  PubMed  Google Scholar 

  46. Bailey TL, Elkan C: The value of prior knowledge in discovering motifs with MEME. Proc Int Conf Intell Syst Mol Biol. 1995, 3: 21-29.

    PubMed  CAS  Google Scholar 

  47. Lorkovic ZJ, Barta A: Genome analysis: RNA recognition motif (RRM) and K homology (KH) domain RNA-binding proteins from the flowering plant Arabidopsis thaliana. Nucleic Acids Res. 2002, 30 (3): 623-635. 10.1093/nar/30.3.623.

    Article  PubMed  CAS  Google Scholar 

  48. Vision TJ, Brown DG, Tanksley SD: The origins of genomic duplications in Arabidopsis. Science. 2000, 290 (5499): 2114-2117. 10.1126/science.290.5499.2114.

    Article  PubMed  CAS  Google Scholar 

  49. Bowers JE, Chapman BA, Rong J, Paterson AH: Unravelling angiosperm genome evolution by phylogenetic analysis of chromosomal duplication events. Nature. 2003, 422 (6930): 433-438. 10.1038/nature01521.

    Article  PubMed  CAS  Google Scholar 

  50. Chromosome Map Tool. [http://www.arabidopsis.org/jsp/ChromosomeMap/tool.jsp]

  51. NCBI. [http://www.ncbi.nlm.nih.gov/]

  52. TAIR. [http://www.Arabidopsis.org/]

  53. Blanc G, Hokamp K, Wolfe KH: A recent polyploidy superimposed on older large-scale duplications in the Arabidopsis genome. Genome Res. 2003, 13 (2): 137-144. 10.1101/gr.751803.

    Article  PubMed  CAS  Google Scholar 

  54. Blanc G, Wolfe KH: Widespread paleopolyploidy in model plant species inferred from age distributions of duplicate genes. Plant Cell. 2004, 16 (7): 1667-1678. 10.1105/tpc.021345.

    Article  PubMed  CAS  Google Scholar 

  55. Blanc G, Wolfe KH: Functional divergence of duplicated genes formed by polyploidy during Arabidopsis evolution. Plant Cell. 2004, 16 (7): 1679-1691. 10.1105/tpc.021410.

    Article  PubMed  CAS  Google Scholar 

  56. NSF Genomics of Loblolly Pine Embryogenesis Project. [http://www.tigr.org/tdb/e2k1/pine/index.shtml]

  57. Crooks GE, Hon G, Chandonia JM, Brenner SE: WebLogo: a sequence logo generator. Genome Res. 2004, 14 (6): 1188-1190. 10.1101/gr.849004.

    Article  PubMed  CAS  Google Scholar 

  58. Meyers BC, Vu TH, Tej SS, Ghazal H, Matvienko M, Agrawal V, Ning J, Haudenschild CD: Analysis of the transcriptional complexity of Arabidopsis thaliana by massively parallel signature sequencing. Nat Biotechnol. 2004, 22 (8): 1006-1011. 10.1038/nbt992.

    Article  PubMed  CAS  Google Scholar 

  59. UniGene of NCBI. [http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db5unigene]

  60. Martin C, Paz-Ares J: MYB transcription factors in plants. Trends Genet. 1997, 13 (2): 67-73. 10.1016/S0168-9525(96)10049-4.

    Article  PubMed  CAS  Google Scholar 

  61. Bork P: Hundreds of ankyrin-like repeats in functionally diverse proteins: mobile modules that cross phyla horizontally?. Proteins. 1993, 17 (4): 363-374. 10.1002/prot.340170405.

    Article  PubMed  CAS  Google Scholar 

  62. Lin JH, Makris A, McMahon C, Bear SE, Patriotis C, Prasad VR, Brent R, Golemis EA, Tsichlis PN: The ankyrin repeat-containing adaptor protein Tvl-1 is a novel substrate and regulator of Raf-1. J Biol Chem. 1999, 274 (21): 14706-14715. 10.1074/jbc.274.21.14706.

    Article  PubMed  CAS  Google Scholar 

  63. Mosavi LK, Minor DL, Peng ZY: Consensus-derived structural determinants of the ankyrin repeat motif. Proc Natl Acad Sci USA. 2002, 99 (25): 16029-16034. 10.1073/pnas.252537899.

    Article  PubMed  CAS  Google Scholar 

  64. Becerra C, Jahrmann T, Puigdomenech P, Vicient CM: Ankyrin repeat-containing proteins in Arabidopsis: characterization of a novel and abundant group of genes coding ankyrin-transmembrane proteins. Gene. 2004, 340 (1): 111-121. 10.1016/j.gene.2004.06.006.

    Article  PubMed  CAS  Google Scholar 

  65. la Cour T, Gupta R, Rapacki K, Skriver K, Poulsen FM, Brunak S: NESbase version 1.0: a database of nuclear export signals. Nucleic Acids Res. 2003, 31 (1): 393-396. 10.1093/nar/gkg101.

    Article  PubMed  CAS  Google Scholar 

  66. Gama-Carvalho M, Carmo-Fonseca M: The rules and roles of nucleocytoplasmic shuttling proteins. FEBS Lett. 2001, 498 (2–3): 157-163. 10.1016/S0014-5793(01)02487-5.

    Article  PubMed  CAS  Google Scholar 

  67. AbuQamar S, Chen X, Dhawan R, Bluhm B, Salmeron J, Lam S, Dietrich RA, Mengiste T: Expression profiling and mutant analysis reveals complex regulatory networks involved in Arabidopsis response to Botrytis infection. Plant J. 2006, 48 (1): 28-44. 10.1111/j.1365-313X.2006.02849.x.

    Article  PubMed  CAS  Google Scholar 

  68. Genevestigator. [http://www.genevestigator.ethz.ch/at/]

  69. Hudson BP, Martinez-Yamout MA, Dyson HJ, Wright PE: Recognition of the mRNA AU-rich element by the zinc finger domain of TIS11d. Nat Struct Mol Biol. 2004, 11 (3): 257-264. 10.1038/nsmb738.

    Article  PubMed  CAS  Google Scholar 

  70. Guex N, Peitsch MC: SWISS-MODEL and the Swiss-PdbViewer: an environment for comparative protein modeling. Electrophoresis. 1997, 18 (15): 2714-2723. 10.1002/elps.1150181505.

    Article  PubMed  CAS  Google Scholar 

  71. Barreau C, Paillard L, Osborne HB: AU-rich elements and associated factors: are there unifying principles?. Nucleic Acids Res. 2005, 33 (22): 7138-7150. 10.1093/nar/gki1012.

    Article  PubMed  CAS  Google Scholar 

  72. Bakheet T, Williams BR, Khabar KS: ARED 3.0: the large and diverse AU-rich transcriptome. Nucleic Acids Res. 2006, D111-114. 10.1093/nar/gkj052. 34 Database

  73. TIGR. [http://www.tigr.org]

  74. DATF. [http://datf.cbi.pku.edu.cn/index.php]

  75. MAtDB. [http://mips.gsf.de/proj/plant/jsf/athal/index.jsp]

  76. TIGR Rice Genome Annotation Database and Resource. [http://www.tigr.org/tdb/e2k1/osa1/]

  77. Rice Genome Database-japonica of the Rice Genome Research Program. [http://rgp.dna.affrc.go.jp/]

  78. The International Rice Genome Sequencing Project. [http://rgp.dna.affrc.go.jp]

  79. DBSubLoc. [http://www.bioinfo.tsinghua.edu.cn/dbsubloc.html]

  80. Paralogons in Arabidopsis. [http://wolfe.gen.tcd.ie/athal/dup]

  81. Physcomitrella EST Project. [http://moss.nibb.ac.jp/]

Download references

Acknowledgements

We are grateful to Dr. Perry J. Blackshear and Brandon Cuthbertson for their suggestion to construct the structure model using the Swiss PDB Viewer. This work was supported by the National Basic Research Program (Grant No. 2006CB1001006) and the National Natural Science Foundation (Grant No. 30570144 and 30500042) in China.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Chengchao Zheng.

Additional information

Authors' contributions

D.W. carried out all the analyses and interpreted the results. D.W. and Y.Y.L. jointly wrote the manuscript. Y.H.G. and C.C.Z contributed with the CCCH gene family background knowledge and edited the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12864_2007_1238_MOESM1_ESM.pdf

Additional file 1: Figure S1. The program imports the BLAST result including the query sites and subject sites to the excel file. (PDF 37 KB)

12864_2007_1238_MOESM2_ESM.pdf

Additional file 2: Figure S2. The program reads dataset from "TAIR6_pep_20060907" (the file of Arabidopsis proteome) or "rice.pep" (the file of rice proteome) and input it into MySQL local database. (PDF 25 KB)

Additional file 3: Figure S3. The program detects the putative CCCH proteins from Arabidopsis proteome. (PDF 30 KB)

Additional file 4: Figure S4. The program sorts the CCCH motifs by their orders in CCCH protein sequences. (PDF 30 KB)

Additional file 5: Figure S5. The program parses the database to search motifs that contain each other. (PDF 28 KB)

Additional file 6: Table S1. Detailed information of CCCH gene family in Arabidopsis. (XLS 62 KB)

Additional file 7: Figure S6. The program detects the putative CCCH proteins from rice proteome. (PDF 30 KB)

Additional file 8: Table S2. Detailed information of CCCH gene family in rice. (XLS 64 KB)

12864_2007_1238_MOESM9_ESM.doc

Additional file 9: Figure S7. Phylogenetic trees of CCCH gene family. A, Phylogenetic tree of rice CCCH gene family. The unrooted tree, constructed using ClustalX (1.83), summarizes the evolutionary relationship among the members of CCCH families in rice. The neighbor-joining tree was constructed using aligned full-length amino acid sequences. The tree shows the 8 major phylogenetic subfamilies (numbered I to VIII and marked with different alternating tones of a gray background to make subfamily identification easier) with high predictive value. Numbers on branches indicate the bootstrap values (1000 replicates) that support the adjacent node. The gene structure of rice are shown on the right side (boxes represent exons and lines represent introns). B, Joined phylogenetic tree of the rice and Arabidopsis CCCH gene families. The unrooted tree was inferred by the neighbor-joining method after the alignment of the full-length amino acid sequences of the 135 Arabidopsis and rice genes listed in Table 2 and Table 3 respectively. The tree shows the 20 phylogenetic subfamilies (numbered I to XX and marked with different alternating tones of a gray background to make subfamily identification easier) with high predictive value. Numbers on branches indicate the bootstrap values (1000 replicates) that support the adjacent node. (DOC 2 MB)

Additional file 10: Figure S8. Chromosomal distribution for rice CCCH genes. (DOC 152 KB)

Additional file 11: Figure S9. The program detects the putative NES sequences form CCCH proteins. (PDF 30 KB)

12864_2007_1238_MOESM12_ESM.doc

Additional file 12: Table S3: The putative NES sequences in CCCH proteins detecting by regular expression consensus [LV]-x(2,3)-[LIVFM]-x(2,3)-L-x-[LIMTKD] (DOC 95 KB)

12864_2007_1238_MOESM13_ESM.xls

Additional file 13: Table S4. The gene locus containing WWWT(ATTTA)TTTW pattern in 3'-UTR of Arabidopsis genome. (XLS 36 KB)

Additional file 14: Table S5: Primers for RT-PCRs. (DOC 34 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Wang, D., Guo, Y., Wu, C. et al. Genome-wide analysis of CCCH zinc finger family in Arabidopsis and rice. BMC Genomics 9, 44 (2008). https://doi.org/10.1186/1471-2164-9-44

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-9-44

Keywords