Genome-wide and expression analysis of protein phosphatase 2C in rice and Arabidopsis
BMC Genomics volume 9, Article number: 550 (2008)
The protein phosphatase 2Cs (PP2Cs) from various organisms have been implicated to act as negative modulators of protein kinase pathways involved in diverse environmental stress responses and developmental processes. A genome-wide overview of the PP2C gene family in plants is not yet available.
A comprehensive computational analysis identified 80 and 78 PP2C genes in Arabidopsis thaliana (AtPP2Cs) and Oryza sativa (OsPP2Cs), respectively, which denotes the PP2C gene family as one of the largest families identified in plants. Phylogenic analysis divided PP2Cs in Arabidopsis and rice into 13 and 11 subfamilies, respectively, which are supported by the analyses of gene structures and protein motifs. Comparative analysis between the PP2C genes in Arabidopsis and rice identified common and lineage-specific subfamilies and potential 'gene birth-and-death' events. Gene duplication analysis reveals that whole genome and chromosomal segment duplications mainly contributed to the expansion of both OsPP2Cs and AtPP2Cs, but tandem or local duplication occurred less frequently in Arabidopsis than rice. Some protein motifs are widespread among the PP2C proteins, whereas some other motifs are specific to only one or two subfamilies. Expression pattern analysis suggests that 1) most PP2C genes play functional roles in multiple tissues in both species, 2) the induced expression of most genes in subfamily A by diverse stimuli indicates their primary role in stress tolerance, especially ABA response, and 3) the expression pattern of subfamily D members suggests that they may constitute positive regulators in ABA-mediated signaling pathways. The analyses of putative upstream regulatory elements by two approaches further support the functions of subfamily A in ABA signaling, and provide insights into the shared and different transcriptional regulation machineries in dicots and monocots.
This comparative genome-wide overview of the PP2C family in Arabidopsis and rice provides insights into the functions and regulatory mechanisms, as well as the evolution and divergence of the PP2C genes in dicots and monocots. Bioinformatics analyses suggest that plant PP2C proteins from different subfamilies participate in distinct signaling pathways. Our results have established a solid foundation for future studies on the functional divergence in different PP2C subfamilies.
The reversible phosphorylation of proteins is a fundamental mechanism by which living organisms modulate cellular processes including cell cycle events, growth factor responses, hormone and environmental stimuli responses, metabolic control, and developmental processes [1–5]. Protein phosphatases (PPs), by reversing the action of protein kinases, provide modulations of protein phosphoregulation. PPs can be grouped into two major classes based on substrate specificity: protein tyrosine phosphatases and protein serine/threonine phosphatases . Based on enzymological criteria using specific inhibitors and activators, most serine/threonine protein phosphatases have been placed into two subcategories: the protein phosphatase 1 (PP1) and protein phosphatase 2 (PP2) . PP2 proteins were further distinguished by metal ion requirements: while protein phosphatases 2A (PP2As) have no ion requirement, PP2Bs require Ca2+ and PP2Cs require Mg2+ or Mn2+. Alternatively, based on both distinct amino acid sequences and crystal structures, protein serine/threonine phosphatases can be classified into the Phosphor-protein phosphatase (PPP) and the Mg2+- or Mn2+- dependent protein phosphatase (PPM) families. The PPP family covers the PP1, PP2A and PP2B groups, whereas the PPM family includes the PP2C group and also pyruvate dehydrogenase phosphatases .
Members of the PP2C group from various organisms have been implicated in acting as negative modulators of protein kinase pathways activated by diverse environmental stresses or developmental signaling cascades [8–15]. In plants, several PP2Cs have been described as negative regulators within the abscisic acid (ABA) mediated signaling network [16–20]. Two dominant ABA insensitive mutants, abi1-1 and abi2-1, are impaired in PP2C genes (ABI1 and ABI2). Proteins encoded by these genes are important for plant tolerance to several abiotic stresses including salt, drought, and freezing [16, 21–30]. Also HAB1 and PP2CA, homologous PP2Cs from Arabidopsis, are negative regulators involved in ABA signaling [19, 20, 23, 31–33]. Anti-sense mediated down-regulation of PP2CA genes accelerates plant development and leads to freezing tolerance . Two PP2C genes (FsPP2C1 and FsPP2C2) were isolated from beech seeds (Fagus sylvatica). Ectopic expression of FsPP2C1 in Arabidopsis, which shows high sequence similarity with ABI1 and ABI2, confers ABA insensitivity during seed dormancy and germination under unfavorable conditions, indicating that FsPP2C1 also acts as a negative regulator of ABA signaling . However, ectopic expression of FsPP2C2 in Arabidopsis resulted in enhanced sensitivity to ABA and abiotic stress in seeds and vegetative tissues and caused a dwarf phenotype and delayed flowering, which suggests that FsPP2C2 is a positive regulator of ABA .
Other plant PP2C members, such as MP2C from Mesembryanthemum crystallinum, act as counterparts of the mitogen-activated protein kinase pathway during salt stress or wounding . A kinase-associated protein phosphatase 2C (KAPP) was originally identified as an interacting protein of the receptor-like kinase CLAVATA1 (CLV1). KAPP acts as a negative modulator of CLV1 and is required for the proper balance between cell proliferation and differentiation in Arabidopsis shoots and floral meristems [36, 37]. KAPP is also a component of a novel Na+ acclimation pathway . Moreover, POLTERGEIST (POL), encoding another PP2C in Arabidopsis, was also identified as a modulator of the CLAVATA pathways controlling stem cell identity in meristems .
The availability of the complete Arabidopsis and Oryza sativa genome sequences allows a timely, systematic genome-wide comparative analysis of gene families between monocots and eudicots [40–44]. Previously, 76 genes from Arabidopsis were identified as PP2C candidates ; however, important features like gene structures, protein motifs and gene duplication events were not explored, which is critical to understanding the characteristics of the Arabidopsis PP2C gene (AtPP2C) family. Moreover, the rice PP2C gene (OsPP2C) family has not been analyzed in detail, and the phylogenetic relationship of the PP2C genes between Arabidopsis and rice remains poorly understood.
In this study, we first identified all the PP2C genes from Arabidopsis and rice, and then conducted phylogenetic analyses to divide them into subfamilies. We further compared the rice and Arabidopsis PP2C genes to identify both shared and specific subfamilies, as well as potential gene birth-and-death events. Following the analysis of gene structure and protein motifs, we traced gene duplication events that likely contributed to the expansion of the PP2C family. The expression patterns of the PP2Cs from rice and Arabidopsis were analyzed and compared, in order to shed light on their functional roles, especially in ABA-signaling and stress tolerance. The upstream regulatory regions of the subfamily A members were finally analyzed in order to provide insight into the transcriptional regulation machineries that define the temporal and spatial expression and the stimuli responses.
Results and discussion
Identification of PP2C genes in Arabidopsis and rice
To uncover the entire family of genes coding for AtPP2C proteins in Arabidopsis, we performed multiple keyword searches to obtain PP2C protein sequences from several databases. Sequences from these hits were used as queries employing BLAST algorithms to search against the corresponding Arabidopsis sequence data sets (see Methods for details). A custom Perl program was used to parse the results and to exclude the redundant entries from the initial data set. Subsequently, the SMART  and Pfam  motif identification tools were used to identify predicted PP2C domains in all these candidate proteins. As a result of this filtering, 109 non-redundant putative proteins with PP2C domains encoded by 80 genes were extracted from the Arabidopsis genome annotation (Table 1).
To explore the occurrence and size of the PP2C gene family in rice, all predicted AtPP2C protein sequences were initially used as query sequences to search against the TIGR rice database  using BLASTP. Multiple BLAST searches were then performed in several rice databases using rice protein sequences with significant hits as queries (see Methods for details). After confirming the presence of PP2C catalytic domains using SMART and Pfam, 111 putative proteins encoded by 78 genes were identified as PP2C members in rice (Table 2). While scanning for protein catalytic domains, additional specific motifs were found in the regions outside of the PP2C catalytic domain in both Arabidopsis and rice PP2Cs. Transmembrane spanning regions were found in AT1G43900, Os01g37130 and Os02g05630, protein kinase catalytic domains were present in AT2G40860, Os01g36080 and Os11g37540, and a forkhead-associated domain (FHA) was identified in At5G19280 (KAPP).
Previously, 76 genes have been identified as AtPP2C candidate genes , all of which were present in our gene set except for AT1G75010, which does not contain a PP2C catalytic domain according to SMART and Pfam, and was thus excluded from our analysis. Twenty-five and twenty one alternative splicing variants were found in the Arabidopsis and rice PP2C gene families, respectively (Table 1 and 2). Because alternative splice variants share a common coding sequence, they are closely related to each other in multiple sequence alignment and phylogenetic analysis. Thus we selected only a single variant for further analysis. Despite the differences in genome sizes between Arabidopsis and rice (125 Mb and 389 Mb, respectively), the two plant species appear to have a very similar number of genes encoding PP2C proteins. Compared to other gene families analyzed in rice and Arabidopsis, the PP2C gene family is one of the largest families in the plant kingdom. Only six PP2Cs were found in the yeast (Saccharomyces cerevisiae) genome, and no more than fifteen are present in the human genome [48, 49]. The high number of PP2C genes present in the Arabidopsis and rice genome highlights the importance of the PP2C class of protein phosphatases in plants, and suggests that plant PP2Cs may have narrower substrate range than other organisms [16, 34].
Phylogenetic and comparative analysis of PP2C family in Arabidopsis and rice
Using the 80 PP2C full length amino acid sequences from Arabidopsis, a phylogenetic analysis was performed to construct a phylogenetic tree (Figure 1). The results showed that all but 7 AtPP2Cs can be divided into 13 subfamilies. Previously, 76 genes have been divided into 10 subfamilies (A-J) [4, 50]. To be consistent with this past study, we adopted the same name for the nine common subfamilies (A, B, C, D, E, G, H, I, J). Since the members of subfamily F are separated into two subfamilies in our study, we named those two subfamilies F1 and F2. High bootstrap support is observed for each subfamily clade (85% for F2 and > 90% for the remainders), suggesting that the genes in the same subfamily share the same evolutionary origin. Three of the five novel members identified in this study (AT2G30170, AT4G16580, AT5G66720) form a new subfamily, which is designated K. The fourth (AT4G03415) is placed into subfamily E, and the last (AT4G33500) forms a single branch. Subfamily A, G and E are further divided into subgroups with at least 50% bootstrap support. Moreover, subfamilies C and D as well as subfamilies E, L and H constitute sister clades in a monophyletic cluster with high bootstrap support (83% and 86%, respectively), suggesting close evolutionary relationships between the respective subfamilies. Fifty of the AtPP2Cs form 25 paralogous sister pairs, 18 of which have very strong bootstrap support (> 90%).
In order to evaluate the evolutionary relationship among the 78 rice PP2C proteins, we also performed a phylogenetic analysis based on their full length amino acid sequences (Figure 2). Eleven subfamilies are identified as monophyletic clades with at least 50% bootstrap support, and additional seven genes form single branches (Figure 2). Twenty-five sister pairs of paralogous OsPP2Cs are found, eleven of which had very strong bootstrap support (> 90%). Out results suggest a clear paralogous pattern of PP2C gene divergence by gene duplication both in monocots and eudicots.
Structural examination reveals that untranslated regions (UTRs) are present in 71 AtPP2C genes and 70 OsPP2C genes (Figure 1 and 2). The number of introns ranged from zero to twelve in AtPP2Cs and from zero to eighteen in OsPP2Cs. Most members within the same subfamily, both in Arabidopsis and rice, share a similar intron/exon structure and gene length. This conserved intron/exon structure in each subfamily supports their close evolutionary relationship and the introduced classification of subfamilies.
We further investigated the relationship between the Arabidopsis and rice PP2C protein families by generating an alignment of the full-length protein sequences of 78 OsPP2Cs and 80 AtPP2Cs followed by the construction of a phylogenetic tree (Figure 3). While bootstrap values for internal nodes are low, the outer nodes allow the clustering of the PP2C genes of Arabidopsis and rice into 13 subfamilies with higher support values (> 98%). The tree topology, as well as the subfamily organization, resembles the individual trees from rice and Arabidopsis (Figure 1 and 2). The clades defining all 13 subfamilies contain members from both rice and Arabidopsis. Subfamilies J and L each contain only a single rice member, and are thus defined as singletons instead of families in rice. The remaining 11 rice subfamilies are named according to the PP2C subfamily nomenclature in Arabidopsis (Figure 2). Consistent with the AtPP2C gene family analysis, pretty high bootstrap support was also found for a clade combining subfamilies C and D in rice (62%), confirming the close evolutionary relationship between these two subfamilies in both eudicots and monocots. The rice and Arabidopsis genes in each subfamily are thus more closely related to each other than to the PP2C genes of the same species from a different subfamily, suggesting that an ancestral set of PP2C genes defining each subfamily already existed before the monocot-eudicot divergence. Twenty-six putative Arabidopsis/rice orthologs are identified in the tree with high statistical support (indicated by the orange squares in Figure 3). Furthermore, the phylogenetic construction of the PP2C family fits well with the model of birth-and-death evolution in the flowering plant lineage [51, 52]. The branches with more than one AtPP2C or OsPP2C gene likely experienced gene birth due to gene duplication events, whereas the branches with only Arabidopsis or only rice PP2C members likely had gene death. Similarly, other large gene families also experienced birth-and-death evolution, such as the MADS-box gene family and the basic/helix-loop-helix (bHLH) transcription factor family analyzed in both rice and Arabidopsis [53, 54].
Evolution and divergence of PP2C genes in Arabidopsis and rice
Large segmental duplications of chromosomal regions during evolution, followed by gene loss, small-scale duplications and local rearrangements, have created the present complexities of plant genome . Indications of several whole genome duplication have been found in the Arabidopsis genome [56–58]. Likewise, two rounds of large scale genome duplications have occurred in the lineage leading to rice, which cover most rice chromosomes [40, 59].
To investigate the relationship between the genetic divergence within the PP2C family and gene duplication events in Arabidopsis, we located the 80 AtPP2C genes on the Arabidopsis chromosomes based on the location information provided in the TAIR database . The genes are distributed across all five chromosomes except for the upper regions of chromosome 2 (Figure 4). Only one gene cluster was identified on chromosome 3, composed of AT3G63320 and AT3G63340, which share high level of sequence similarity and have thus likely resulted from a tandem gene duplication event. Thirteen paralogous AtPP2C pairs in the AtPP2C gene family were found on segmentally duplicated chromosomal regions (Figure 4A, blue lines), which was originated from a polyploidy event occurred around 24 to 40 million years ago, probably close to the emergence of the crucifer family [61–64]. The putative protein sequences of those 13 pairs were clustered with bootstrap value higher than 70% in the phylogenetic tree (Figure 1). An additional 11 paralogous gene pairs were identified based on phylogenetic relationship and structural similar (Figure 4A, orange squares). Although the chromosomal regions hosting these 11 gene pairs were not included by Blanc et al  as recently duplicated segmental regions, they were considered to be duplicated genes. Together, these genes represent 60% of all AtPP2C s, indicating that large-scale segmental duplication events have contributed largely to the current complexity of the PP2C gene family in Arabidopsis. Only 2.5% of AtPP2Cs (one pair) are organized in clusters, and likely have evolved via local or tandem duplications.
To determine the genomic distribution of the OsPP2C genes, the DNA sequence of each OsPP2C gene was used to search the rice genome database using BLASTN. Although all of the 12 rice chromosomes do contain OsPP2C genes, the distribution appears not random (Figure 4B). Relatively high densities of OsPP2C genes were found in specific chromosomal regions, including the top of chromosome 3, and the bottom of chromosome 2. In contrast, several large chromosomal regions are devoid of OsPP2C genes, such as the top parts of chromosomes 8, 9, and 10 and the central sections of chromosomes 7, 8, and 12. To detect possible relationships between OsPP2C genes and potential genome duplication events, we mapped 17 paralogous pairs of OsPP2C genes. Twelve of these pairs are linked to each of at least 8 previously predicted chromosomal/segmental duplications (Figure 4B, green bars, linked with blue lines) [54, 65, 66]. Five additional pairs of OsPP2C genes are close paralogs (Figure 4B, orange squares), but are located in the regions that have not been proposed to be the result of rice genome duplication events. Thus, 44% of the OsPP2C family might have evolved from putative rice genome duplication events. Four OsPP2C gene clusters contain the members with high levels of sequence similarity (Figure 4B, red rectangles), representing a total of 13 OsPP2C genes. Thus 17% of OsPP2Cs are organized in clusters and likely have evolved via local or tandem duplications.
Taken together, it appears that the PP2C gene families expanded differently in rice and Arabidopsis. Whole genome and chromosomal segment duplications mainly contributed to the expansion of both OsPP2Cs and AtPP2Cs giving rise to 44% and 60% of PP2C genes, respectively. But tandem or local duplication occurred less frequently in Arabidopsis (one pair, 2.5% of all PP2Cs) than rice (four clusters, 17%). It is noteworthy that the members of divergent PP2C subfamilies in both Arabidopsis and rice are located within the same chromosomal region, whereas the members of the same subfamily are distributed in different chromosomal regions, suggesting that PP2C genes were distributed widely in the genome of the common ancestor of monocots and eudicots. Furthermore, in both species the duplicated pairs of PP2C genes appear to have been preferentially retained compared with other genes, since the density of duplicated genes retained in recently duplicated chromosomal segments was estimated to be only 28% on average . This highlights the indispensable roles of PP2Cs in both monocots and eudicots.
Protein motifs of PP2Cs in Arabidopsis and rice
To identify the motifs shared among related proteins within the PP2C family, the MEME motif search tool was employed . We identified 11 motifs and their sequences are shown in Table 3; the distribution of these motifs in PP2C proteins in Arabidopsis and rice is shown in Figure 1 and 2. Some motifs are widespread among PP2C proteins (e.g. motifs 1, 2, 3, 4, and 8). In contrast, other motifs are specific to only one or two subfamilies. For example, motifs 5 and 11 are specific to subfamilies C and D, respectively, while motif 9 is found almost exclusively in subfamily D members. Vice versa, motifs 6 and 7 are present in subfamily A, but were not found in subfamily D. Subfamilies C and D have especially distinctive characteristics in terms of motif distribution within the PP2C family.
The different motif distribution in the respective protein sequences may provide an anchor to study the divergence of gene function in different subfamilies. Furthermore, the fact that proteins in the same subfamily exhibit highly similar motif distribution supports their close evolutionary relationship.
Expression patterns of PP2C genes in Arabidopsis and rice
Gene expression patterns can provide important clues for gene function. To date, a comparative analysis of the expression patterns of the whole PP2C gene family has not been performed. Therefore, we firstly examined the expression of AtPP2C genes in inflorescence, leaf, root and silique tissues using microarray data from Genevestigator , massively parallel signature sequencing (MPSS) data  and EST abundance data from NCBI  (Figure 5A). Expression profiles for 67, 73, and 68 PP2C genes were extracted from the Genevestigator, MPSS, and EST databases, respectively. After integrating the three data sets, we found that most of the genes have a very broad expression spectrum, and only ten genes (AT1G07160, AT2G40180, AT3G09400, AT3G27140, AT3G23360, AT4G03415, AT4G08260, AT4G11040, AT5G24940, and AT5G51760) were not detected in any of those four tissues. However, transcript evidence for all genes is present in the mix tissue samples from Arabidopsis with the sole exception of AT4G11040 and AT3G27140, indicating that they might be expressed only at highly specific developmental stages. Since no evidence for the expression of AT3G27140 was found in full-length cDNA datasets at NCBI either (Table 1), it might be a pseudogene. We also summarized the expression of rice PP2C genes in leaf, flower, root, seed and stem tissues using the same approach (Figure 5B). Expression signatures for 71 OsPP2C genes were detected in the MPSS database, and 69 members had matching ESTs. No expression evidence was found in any database for three OsPP2C genes (Os02g39480, Os02g42250, Os02g42270) (Table 2), suggesting that they might be pseudogenes. Twenty OsPP2C genes had no detectable tissue specific expression, but had matching transcript sequences (Figure 5B), which might thus be expressed at specific developmental stages or under special conditions only.
A total of 49 AtPP2C genes were expressed in all tissues examined. Only three AtPP2C genes (A2G35350, AT3G05640 and AT5G02440) displayed highly preferential expression in certain tissue (Figure 5). Twenty-seven OsPP2C genes were expressed in all tissues examined, and six OsPP2C genes showed tissue-specific expression. Most AtPP2C genes (84%) and OsPP2C genes (72%) were expressed in more than 2 tissues examined, suggesting that most PP2Cs in plants play functional roles at multiple developmental processes, and only a few members act in tissue-specific biological processes. It is interesting to note that 19 of 29 paralogous gene pairs (18 duplicated gene pairs in Arabidopsis and 11 in rice) show differential expression patterns, suggesting functional diversification of most duplicated gene pairs. This result is consistent with the previous findings that the expression of paralogs frequently diverged to become part of parallel networks, and that the expression of each gene is strongly correlated with other non-homologous genes acting in the corresponding networks, suggesting that neo-functionalization is a major feature of the long-term evolution of polyploids .
The phytohormone abscisic acid (ABA) plays important regulatory roles in many stress and developmental responses throughout the plant life cycle, particularly in the ability to sense and respond to various unfavorable environmental conditions, including drought, salinity, cold shock, and wounding during vegetative growth [21, 33, 70–74]. To investigate possible involvement of PP2Cs in ABA mediated signaling, we examined the expression of all AtPP2Cs in response to ABA treatment and several environmental stimuli using the Genevestigator microarray dataset  (Figure 6). The expression of 16 PP2Cs was found to be up-regulated in response to ABA treatment, including all 9 members of subfamily A. Five AtPP2C genes belonging to subfamily A have been characterized as negative regulators in ABA-mediated processes, including AT4G26080 (ABI1), AT5G57050 (ABI2), AT1G72770 (HAB1), AT1G17550 (HAB2), and AT3G11410 (ATPP2CA) [16–20]. Unlike the PP2C genes from other subfamilies, all members from subfamily A were induced by several stimuli (Figure 6). To explore whether also OsPP2C genes from subfamily A play functional roles in ABA mediated signaling pathways and stress responses, we analyzed the expression of 10 OsPP2C genes in rice plants treated with ABA, salt, osmotic (mannitol) and cold stress by reverse transcriptase polymerase chain reaction (RT-PCR). The results of our experiment indicate elevated expression levels of all OsPP2C genes from subfamily A upon treatment with ABA and salt (Figure 7). Also mannitol and cold treatment led to induced expression of several subfamily A members, but expression and induction levels differed among these genes (Figure 7). Thus, the expression pattern of the OsPP2C subfamily A genes is in good agreement with the microarray data for Arabidopsis subfamily A members, suggesting that the members of this subfamily play foremost roles in ABA-mediated processes related to stress responses both in monocots and eudicots.
Interestingly, the expression of 7 PP2C genes in Arabidopsis was suppressed by ABA treatment, and 2 of them are the members of subfamily D (Fig. 6). In previous studies, a PP2C from Fagus sylvatica (FsPP2C2) was described as a positive regulator of the ABA signal transduction pathway, contrary to most other PP2Cs described so far . Based on phylogenetic reconstructions, FsPP2C2 belongs to subfamily D and is most closely related to AT3G51370 (data no shown). This result suggests that at least some members of subfamily D might have a positive role in ABA signaling pathways.
The expression of 30 AtPP2C genes was not influenced by different environmental stimuli (Fig. 6), suggesting that 38% of AtPP2Cs might not play functional roles in response to environmental stress signals. Among them, five genes are from subfamily C, which is composed of POLTERGEIST (POL, AT2G46920), five POL-like genes (PLL1-PLL5, gene IDs in Table 1) and AT3G16560. POL and PLL1 have been shown to regulate meristem development and pedicel length, while PLL4 and PLL5 were found regulating leaf development [39, 75]. These results indicate that at least some of these 38% AtPP2Cs (e.g. some members of subfamily C) might be involved in developmental processes.
In summary, our expression analyses provide evidences that PP2C genes in different subfamilies might play different functional roles in distinct signaling pathways.
Upstream regulatory regions of subfamily A PP2C genes in Arabidopsis and rice
Promoters in the upstream region of genes play key roles in conferring developmental and/or environmental regulation of gene expression. To further elucidate the transcriptional regulation machineries of the 19 PP2C genes in subfamily A in Arabidopsis and rice, two approaches were attempted in this study.
First, we analyzed the 1000-base-pair (bp) upstream promoter sequences of these AtPP2C genes using the PlantCARE database , and validated the results using the PLACE databases . Based upon this approach, a total of 20 types of cis-regulatory elements (CREs) were identified among the 19 subfamily A AtPP2Cs, which do not contain the common eukaryotic regulatory elements, such as TATA and CAAT boxes. The names and functions of these 20 motifs are shown in Table 4. Their amounts and positions in every subfamily A member obtained by the first method are displayed in Additional Files 1 and 2.
Second, we used these 20 motifs as a base for constructing their frequency matrix following a frequency distribution approach (see Methods and ), to compare the frequency of these motifs in the promoter regions of subfamily A with that in the whole genome, in order to explore whether these 20 motifs are overrepresented in subfamily A in Arabidopsis and rice. The frequency matrices of 9 motifs (ABRE, CE3, G-box, GCN4 motif, HSE, O2-site, TC-rich repeats, TCA-element and TGACG motif) were generated. Only a very limited number of sequences containing the remaining 11 motifs (ACE, ARE, AuxRR-core, CAT-box, CCGTCC-box, CGTCA motif, circadian, LTR, motif 1, Skn-1 motif, 5UTR Py-rich stretch) are present in the PLACE, Plantcare and other databases, therefore their frequency matrices could not be generated in this analysis. Thus, only the frequency matrices of 9 motifs were used to identify all single occurrences of these motifs in the promoter regions of the subfamily A members and all protein-coding genes using the genomic sequences of Arabidopsis (genome release 7.0, January 2007)  and rice (release 4.0, January 2006)  (Table 5). Our results suggest that in both species, 5 out of 9 motifs (GCN4, HSE, O2-site, TC-rich repeats and TCA element) are not over-represented in the upstream regions of PP2C subfamily A members. In contrast, in both species the frequency of ABREs (involved in ABA responsiveness) is much higher in the promoter regions of subfamily A than in the whole genome (Tab. 5). The CE3 (involved in ABA and VP1 responsiveness) and G-box (involved in light responsiveness) motifs were found only in subfamily A promoters from rice but absent in those from Arabidopsis. Our results support a previous finding that CE3s are abundant in rice but almost absent in Arabidopsis . In reverse, the TGACG motif (involved in MeJA responsiveness) appeared relatively abundant in subfamily A promoters from Arabidopsis while absent in those from rice. The different motif distributions in Arabidopsis and rice point to evolutionary divergent patterns of PP2C gene regulatory networks. The abundance of ABA responsive elements in subfamily A PP2C gene promoters provides additional evidences for their role in ABA signaling, and reveals shared (ABRE) as well as different (CE3) gene regulation machineries for ABA responsiveness in dicots and monocots. The occurrence of other motifs like G-box suggests that some subfamily A PP2C genes might be involved also in other biological processes, such as light-responsive pathways.
Additional information on the frequency distributions of ABRE, CE3 and G-box in the 1000 bp upstream regions and their nucleotide frequency distributions in Arabidopsis and rice are provided in Fig. 8. ABRE, CE3 and G-box are widely distributed in the upstream regions of subfamily A as well as in the whole genome of Arabidopsis and rice (data no shown).
Our results also indicate that although the mapping of known motifs to the promoter regions of genes (i.e. our first approach) can reveal interesting patterns and insights, it is necessary to compare the frequency of these motifs in the genes against their genome-wide frequency (i.e. our second approach), in order to explore whether these motifs are truly enriched in the promoter regions of these genes and to eliminate false positive motif instances.
We provide here a comparative genome-wide overview of the PP2C family of Arabidopsis and rice, one of the largest gene families in plants. Our computational studies show that despite the big differences in genome sizes of Arabidopsis and rice, the two species appear to have very similar numbers of genes encoding PP2C proteins (80 and 78, respectively), which are much higher than in yeast and human. The AtPP2Cs and Os PP2cs are divided into 13 and 11 subfamilies, respectively, each of which shares a monophyletic origin, common protein motifs and intron/exon structures. We identified both common and lineage-specific subfamilies and potential 'gene birth-and-death' events in the PP2cs of these two species. Our analyses suggest that whole genome and chromosomal segment duplications mainly contributed to the expansion of both OsPP2Cs and AtPP2Cs, but tandem or local duplication occurred less frequently in Arabidopsis than rice. Our protein motif analysis detected the motifs widespread among the PP2C proteins as well as the motifs specific to only one or two subfamilies. Our computational and experimental expression analyses suggest that 1) most PP2C genes play functional roles in multiple tissues in both species, 2) most genes in subfamily A play primary roles in stress tolerance, especially ABA response, 3) subfamily D members may constitute positive regulators in ABA-mediated signaling pathways. We further employed two approaches to analyze the putative upstream regulatory elements of PP2C subfamily members. The results obtained by both methods support the function of subfamily A in ABA responses, and provide insights into the shared and different transcriptional regulation machineries as well as evolutionary divergent patterns of PP2C gene regulatory networks in dicots and monocots. Overall, our study provides comprehensive information and novel insights into the functions and regulatory mechanisms, as well as the evolution and divergence of the PP2C genes in dicots and monocots. The results presented here lay a solid foundation for future studies on the functional divergence of different PP2c subfamilies.
Materials and methods
Identification of PP2C genes in Arabidopsis and rice
To collect all members of the PP2C gene family in Arabidopsis, an initial keyword search was performed using "protein phosphatase 2C" to screen six databases, i.e. NCBI , DATF , MAtDB , TAIR , TIGR  and PlantsP . The resulting protein hits were pooled together and redundant sequences were removed using a custom Perl program. Remaining proteins sequences were then used as a query to perform multiple database searches against proteome and genome files downloaded from TAIR. Stand-alone versions of BLASTP and TBLASTN available from NCBI were used with the e-value cutoff set to 1e-003. After eliminating redundant sequences, the obtained protein sequences were screened for the existence of a PP2C catalytic domain by the hidden Markov model provided with the SMART  and Pfam  tools. Proteins without a PP2C catalytic domain were eliminated from the dataset.
To collect all members of the PP2C gene family in rice, all predicted AtPP2C protein sequences were initially used as query sequences to search against the TIGR rice database  using BLASTP. Multiple BLAST searches using protein sequences of significant hits as the queries were then performed in several rice databases, including NCBI , TIGR , Rice Genome Database-japonica of the Rice Genome Research Program  and the International Rice Genome Sequencing Project (IRSGP) . After removal of redundant sequences from our data set, SMART and Pfam were used to confirm the presence of PP2C catalytic domains as described above.
Alignment and phylogenetic analysis of PP2C sequences
Multiple alignments of amino acid sequences were generated using ClustalX (1.83) and were manually corrected. Phylogenetic trees were constructed by neighbor-joining algorithms of Mega (4.0). Bootstrapping was performed 1000 times to obtain support values for each branch.
Chromosomal localization and duplication of PP2C genes
To determine the location of PP2C genes on the five Arabidopsis chromosomes, the Chromosome Map Tool from TAIR was used . 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 . For rice, all the sequenced contigs of japonica cultivar Nipponbare has been physically reconstructed as pseudomolecules by the IRGSP representing the 12 rice chromosomes, which are available in GeneBank. Each of the rice PP2C genes was positioned on these rice chromosome pseudomolecules using BLASTN. For the detection of large segmental duplications, we refer to the duplicated blocks map generated by Xiong et al. .
Motifs identification of PP2C proteins in Arabidopsis and rice
To identify motifs shared among related proteins within PP2C gene family, we used the MEME motif search tool  with default settings except that the maximum number of motifs to find was defined as 20 and the maximum width was set to 300.
In silico expression analysis of PP2C genes
cDNA signatures from MPSS  and EST abundance data from NCBI  were used to detect the expression patterns of Arabidopsis and rice PP2C genes in different organs. Microarray data from Genevestigiator  were used to determine the expression profile of Arabidopsis PP2C genes in different organs and under multiple environmental stimuli.
RT-PCR analysis of PP2C gene expression
Plant tissues of the rice japonica cultivar Nipponbare were ground in liquid nitrogen. For reverse transcriptase (RT) mediated PCR analysis, total RNA was isolated using the Trizol (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's instructions. RNA preparations were then treated with DNase I (Promega, Madison, WI, USA). First strand cDNA synthesis was performed by using oligo (dT) primer and M-MLV RT (Promega, Madison, WI, USA). 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 3.
A total of 32,041 1000-bp upstream sequences from Arabidopsis were downloaded from TAIR 7.0 . All upstream sequences that correspond to genes that do not code for proteins or oligopeptides were discarded (5012 sequences). Then, the reverse complementary sequences were obtained. The final database of Arabidopsis 1-kb upstream sequences consisted of 54,058 individual sequences. A database of 66,170 1000-bp upstream sequences from rice was downloaded from TIGR 4.0 . Upstream sequences were cleared of transposable elements and related gene models using the gene identifiers of transposable-element-related gene models provided by TIGR (15,424 gene identifiers). After obtaining reverse complementary sequences, the file used for the screening consisted of 82,144 sequences. The 1000-bp upstream sequences of PP2C genes in subfamily A were extracted from the final databases from Arabidopsis and rice.
Randomized sequences were generated by shuffling the extracted upstream sequences using the EMBOSS program "shuffleseq" . During the shuffling procedure, the single-nucleotide frequency of each sequence is maintained, but the nucleotide order is randomly shuffled. The reshuffling process was repeated 10 times to generate 10 different data sets of randomized 1-kb upstream sequences for Arabidopsis and rice, respectively.
Searching single occurrences of CREs
The subfamily A, genomic and randomized data sets were screened for single occurrences of CREs on the plus and minus strands, using the frequency matrices and the program "Profit" from EMBOSS . For each screened data set the program generates a list of matches, including the sequence identifier, element, score and position for every match. The score is calculated from the frequency matrix. An arbitrary 95% cut-off was used. The position of a match is expressed in terms of the distance to the translation initiation codon (ATG).
Andreeva AV, Kutuzov MA: RdgC/PP5-related phosphatases: novel components in signal transduction. Cellular signalling. 1999, 11 (8): 555-562. 10.1016/S0898-6568(99)00032-7.
Chernoff J: Protein tyrosine phosphatases as negative regulators of mitogenic signaling. Journal of cellular physiology. 1999, 180 (2): 173-181. 10.1002/(SICI)1097-4652(199908)180:2<173::AID-JCP5>3.0.CO;2-Y.
den Hertog J: Protein-tyrosine phosphatases in development. Mech Dev. 1999, 85: 1-2.
Kerk D, Bulgrien J, Smith DW, Barsam B, Veretnik S, Gribskov M: The complement of protein phosphatase catalytic subunits encoded in the genome of Arabidopsis. Plant Physiol. 2002, 129 (2): 908-925. 10.1104/pp.004002.
Schweighofer A, Kazanaviciute V, Scheikl E, Teige M, Doczi R, Hirt H, Schwanninger M, Kant M, Schuurink R, Mauch F: The PP2C-type phosphatase AP2C1, which negatively regulates MPK4 and MPK6, modulates innate immunity, jasmonic acid, and ethylene levels in Arabidopsis. The Plant cell. 2007, 19 (7): 2213-2224. 10.1105/tpc.106.049585.
Cohen P: The structure and regulation of protein phosphatases. Annual review of biochemistry. 1989, 58: 453-508. 10.1146/annurev.bi.58.070189.002321.
Cohen PT: Novel protein serine/threonine phosphatases: variety is the spice of life. Trends in biochemical sciences. 1997, 22 (7): 245-251. 10.1016/S0968-0004(97)01060-8.
Hanada M, Kobayashi T, Ohnishi M, Ikeda S, Wang H, Katsura K, Yanagawa Y, Hiraga A, Kanamaru R, Tamura S: Selective suppression of stress-activated protein kinase pathway by protein phosphatase 2C in mammalian cells. FEBS letters. 1998, 437 (3): 172-176. 10.1016/S0014-5793(98)01229-0.
Takekawa M, Maeda T, Saito H: Protein phosphatase 2Calpha inhibits the human stress-responsive p38 and JNK MAPK pathways. The EMBO journal. 1998, 17 (16): 4744-4752. 10.1093/emboj/17.16.4744.
Maeda T, Wurgler-Murphy SM, Saito H: A two-component system that regulates an osmosensing MAP kinase cascade in yeast. Nature. 1994, 369 (6477): 242-245. 10.1038/369242a0.
Gaits F, Shiozaki K, Russell P: Protein phosphatase 2C acts independently of stress-activated kinase cascade to regulate the stress response in fission yeast. The Journal of biological chemistry. 1997, 272 (28): 17873-17879. 10.1074/jbc.272.28.17873.
Sheen J: Ca2+-dependent protein kinases and stress signal transduction in plants. Science. 1996, 274 (5294): 1900-1902. 10.1126/science.274.5294.1900.
Madrid M, Nunez A, Soto T, Vicente-Soler J, Gacto M, Cansado J: Stress-activated protein kinase-mediated down-regulation of the cell integrity pathway mitogen-activated protein kinase Pmk1p by protein phosphatases. Mol Biol Cell. 2007, 18 (11): 4405-4419. 10.1091/mbc.E07-05-0484.
Ruan H, Yan Z, Sun H, Jiang L: The YCR079w gene confers a rapamycin-resistant function and encodes the sixth type 2C protein phosphatase in Saccharomyces cerevisiae. FEMS Yeast Res. 2007, 7 (2): 209-215. 10.1111/j.1567-1364.2006.00160.x.
Dai J, Zhang J, Sun Y, Wu Q, Sun L, Ji C, Gu S, Feng C, Xie Y, Mao Y: Characterization of a novel human protein phosphatase 2C family member, PP2Ckappa. Int J Mol Med. 2006, 17 (6): 1117-1123.
Merlot S, Gosti F, Guerrier D, Vavasseur A, Giraudat J: The ABI1 and ABI2 protein phosphatases 2C act in a negative feedback regulatory loop of the abscisic acid signalling pathway. Plant J. 2001, 25 (3): 295-303. 10.1046/j.1365-313x.2001.00965.x.
Tahtiharju S, Palva T: Antisense inhibition of protein phosphatase 2C accelerates cold acclimation in Arabidopsis thaliana. Plant J. 2001, 26 (4): 461-470. 10.1046/j.1365-313X.2001.01048.x.
Gonzalez-Garcia MP, Rodriguez D, Nicolas C, Rodriguez PL, Nicolas G, Lorenzo O: Negative regulation of abscisic acid signaling by the Fagus sylvatica FsPP2C1 plays a role in seed dormancy regulation and promotion of seed germination. Plant Physiol. 2003, 133 (1): 135-144. 10.1104/pp.103.025569.
Saez A, Apostolova N, Gonzalez-Guzman M, Gonzalez-Garcia MP, Nicolas C, Lorenzo O, Rodriguez PL: Gain-of-function and loss-of-function phenotypes of the protein phosphatase 2C HAB1 reveal its role as a negative regulator of abscisic acid signalling. Plant J. 2004, 37 (3): 354-369. 10.1046/j.1365-313X.2003.01966.x.
Yoshida T, Nishimura N, Kitahata N, Kuromori T, Ito T, Asami T, Shinozaki K, Hirayama T: ABA-hypersensitive germination3 encodes a protein phosphatase 2C (AtPP2CA) that strongly regulates abscisic acid signaling during germination among Arabidopsis protein phosphatase 2Cs. Plant Physiol. 2006, 140 (1): 115-126. 10.1104/pp.105.070128.
Leung J, Giraudat J: Abscisic Acid Signal Transduction. Annu Rev Plant Physiol Plant Mol Biol. 1998, 49: 199-222. 10.1146/annurev.arplant.49.1.199.
Quintero FJ, Ohta M, Shi H, Zhu JK, Pardo JM: Reconstitution in yeast of the Arabidopsis SOS signaling pathway for Na+ homeostasis. Proceedings of the National Academy of Sciences of the United States of America. 2002, 99 (13): 9061-9066. 10.1073/pnas.132092099.
Sheen J: Mutational analysis of protein phosphatase 2C involved in abscisic acid signal transduction in higher plants. Proceedings of the National Academy of Sciences of the United States of America. 1998, 95 (3): 975-980. 10.1073/pnas.95.3.975.
Leung J, Bouvier-Durand M, Morris PC, Guerrier D, Chefdor F, Giraudat J: Arabidopsis ABA response gene ABI1: features of a calcium-modulated protein phosphatase. Science. 1994, 264 (5164): 1448-1452. 10.1126/science.7910981.
Meyer K, Leube MP, Grill E: A protein phosphatase 2C involved in ABA signal transduction in Arabidopsis thaliana. Science. 1994, 264 (5164): 1452-1455. 10.1126/science.8197457.
Leung J, Merlot S, Giraudat J: The Arabidopsis ABSCISIC ACID-INSENSITIVE2 (ABI2) and ABI1 genes encode homologous protein phosphatases 2C involved in abscisic acid signal transduction. The Plant cell. 1997, 9 (5): 759-771. 10.2307/3870430.
Rodriguez PL, Benning G, Grill E: ABI2, a second protein phosphatase 2C involved in abscisic acid signal transduction in Arabidopsis. FEBS letters. 1998, 421 (3): 185-190. 10.1016/S0014-5793(97)01558-5.
Soderman E, Mattsson J, Engstrom P: The Arabidopsis homeobox gene ATHB-7 is induced by water deficit and by abscisic acid. Plant J. 1996, 10 (2): 375-381. 10.1046/j.1365-313X.1996.10020375.x.
Strizhov N, Abraham E, Okresz L, Blickling S, Zilberstein A, Schell J, Koncz C, Szabados L: Differential expression of two P5CS genes controlling proline accumulation during salt-stress requires ABA and is regulated by ABA1, ABI1 and AXR2 in Arabidopsis. Plant J. 1997, 12 (3): 557-569. 10.1046/j.1365-313X.1997.00557.x.
Murata Y, Pei ZM, Mori IC, Schroeder J: Abscisic acid activation of plasma membrane Ca(2+) channels in guard cells requires cytosolic NAD(P)H and is differentially disrupted upstream and downstream of reactive oxygen species production in abi1-1 and abi2-1 protein phosphatase 2C mutants. The Plant cell. 2001, 13 (11): 2513-2523. 10.2307/3871591.
Leonhardt N, Kwak JM, Robert N, Waner D, Leonhardt G, Schroeder JI: Microarray expression analyses of Arabidopsis guard cells and isolation of a recessive abscisic acid hypersensitive protein phosphatase 2C mutant. The Plant cell. 2004, 16 (3): 596-615. 10.1105/tpc.019000.
Saez A, Robert N, Maktabi MH, Schroeder JI, Serrano R, Rodriguez PL: Enhancement of abscisic acid sensitivity and reduction of water consumption in Arabidopsis by combined inactivation of the protein phosphatases type 2C ABI1 and HAB1. Plant Physiol. 2006, 141 (4): 1389-1399. 10.1104/pp.106.081018.
Kuhn JM, Boisson-Dernier A, Dizon MB, Maktabi MH, Schroeder JI: The protein phosphatase AtPP2CA negatively regulates abscisic acid signal transduction in Arabidopsis, and effects of abh1 on AtPP2CA mRNA. Plant Physiol. 2006, 140 (1): 127-139. 10.1104/pp.105.070318.
Reyes D, Rodriguez D, Gonzalez-Garcia MP, Lorenzo O, Nicolas G, Garcia-Martinez JL, Nicolas C: Overexpression of a protein phosphatase 2C from beech seeds in Arabidopsis shows phenotypes related to abscisic acid responses and gibberellin biosynthesis. Plant Physiol. 2006, 141 (4): 1414-1424. 10.1104/pp.106.084681.
Meskiene I, Baudouin E, Schweighofer A, Liwosz A, Jonak C, Rodriguez PL, Jelinek H, Hirt H: Stress-induced protein phosphatase 2C is a negative regulator of a mitogen-activated protein kinase. The Journal of biological chemistry. 2003, 278 (21): 18945-18952. 10.1074/jbc.M300878200.
Stone JM, Trotochaud AE, Walker JC, Clark SE: Control of meristem development by CLAVATA1 receptor kinase and kinase-associated protein phosphatase interactions. Plant Physiol. 1998, 117 (4): 1217-1225. 10.1104/pp.117.4.1217.
Trotochaud AE, Hao T, Wu G, Yang Z, Clark SE: The CLAVATA1 receptor-like kinase requires CLAVATA3 for its assembly into a signaling complex that includes KAPP and a Rho-related protein. The Plant cell. 1999, 11 (3): 393-406. 10.2307/3870868.
Manabe Y, Bressan RA, Wang T, Li F, Koiwa H, Sokolchik I, Li X, Maggio A: The Arabidopsis kinase-associated protein phosphatase regulates adaptation to Na+ stress. Plant Physiol. 2008, 146 (2): 612-622. 10.1104/pp.107.109009.
Yu LP, Miller AK, Clark SE: POLTERGEIST encodes a protein phosphatase 2C that regulates CLAVATA pathways controlling stem cell identity at Arabidopsis shoot and flower meristems. Curr Biol. 2003, 13 (3): 179-188. 10.1016/S0960-9822(03)00042-3.
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.
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.
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.
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.
Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature. 2000, 408 (6814): 796-815. 10.1038/35048692.
Schultz J, Copley RR, Doerks T, Ponting CP, Bork P: SMART: a web-based tool for the study of genetically mobile domains. Nucleic acids research. 2000, 28 (1): 231-234. 10.1093/nar/28.1.231.
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.
Stark MJ: Yeast protein serine/threonine phosphatases: multiple roles and diverse regulation. Yeast. 1996, 12 (16): 1647-1675. 10.1002/(SICI)1097-0061(199612)12:16<1647::AID-YEA71>3.0.CO;2-Q.
Cheng A, Kaldis P, Solomon MJ: Dephosphorylation of human cyclin-dependent kinases by protein phosphatase type 2C alpha and beta 2 isoforms. The Journal of biological chemistry. 2000, 275 (44): 34744-34749. 10.1074/jbc.M006210200.
Schweighofer A, Hirt H, Meskiene I: Plant PP2C phosphatases: emerging functions in stress signaling. Trends in plant science. 2004, 9 (5): 236-243. 10.1016/j.tplants.2004.03.007.
Nei M, Gu X, Sitnikova T: Evolution by the birth-and-death process in multigene families of the vertebrate immune system. Proceedings of the National Academy of Sciences of the United States of America. 1997, 94 (15): 7799-7806. 10.1073/pnas.94.15.7799.
Nei M, Rogozin IB, Piontkivska H: Purifying selection and birth-and-death evolution in the ubiquitin gene family. Proceedings of the National Academy of Sciences of the United States of America. 2000, 97 (20): 10866-10871. 10.1073/pnas.97.20.10866.
Nam J, dePamphilis CW, Ma H, Nei M: Antiquity and evolution of the MADS-box gene family controlling flower development in plants. Molecular biology and evolution. 2003, 20 (9): 1435-1447. 10.1093/molbev/msg152.
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, 14: 1167-1184. 10.1104/pp.106.080580.
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.
Lynch M, Conery JS: The evolutionary fate and consequences of duplicate genes. Science. 2000, 290 (5494): 1151-1155. 10.1126/science.290.5494.1151.
Simillion C, Vandepoele K, Van Montagu MC, Zabeau M, Peer Van de Y: The hidden duplication past of Arabidopsis thaliana. Proceedings of the National Academy of Sciences of the United States of America. 2002, 99 (21): 13627-13632. 10.1073/pnas.212522399.
Raes J, Vandepoele K, Simillion C, Saeys Y, Peer Van de Y: Investigating ancient duplication events in the Arabidopsis genome. J Struct Funct Genomics. 2003, 3: 1-4. 10.1023/A:1022666020026.
Yu J, Wang J, Lin W, Li S, Li H, Zhou J, Ni P, Dong W, Hu S, Zeng C: The Genomes of Oryza sativa: a history of duplications. PLoS Biol. 2005, 3 (2): e38-10.1371/journal.pbio.0030038.
Blanc G, Hokamp K, Wolfe KH: A recent polyploidy superimposed on older large-scale duplications in the Arabidopsis genome. Genome research. 2003, 13 (2): 137-144. 10.1101/gr.751803.
Kent WJ, Baertsch R, Hinrichs A, Miller W, Haussler D: Evolution's cauldron: duplication, deletion, and rearrangement in the mouse and human genomes. Proceedings of the National Academy of Sciences of the United States of America. 2003, 100 (20): 11484-11489. 10.1073/pnas.1932072100.
Cannon SB, Mitra A, Baumgarten A, Young ND, 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-10.1186/1471-2229-4-10.
Mehan MR, Freimer NB, Ophoff RA: A genome-wide survey of segmental duplications that mediate common human genetic variation of chromosomal architecture. Hum Genomics. 2004, 1 (5): 335-344.
Guyot R, Keller B: Ancestral genome duplication in rice. Genome. 2004, 47 (3): 610-614. 10.1139/g04-016.
Xiong Y, Liu T, Tian C, Sun S, Li J, Chen M: Transcription factors in rice: a genome-wide comparative analysis between monocots and eudicots. Plant molecular biology. 2005, 59 (1): 191-203. 10.1007/s11103-005-6503-6.
Fan LM, Zhao Z, Assmann SM: Guard cells: a dynamic signaling model. Current opinion in plant biology. 2004, 7 (5): 537-546. 10.1016/j.pbi.2004.07.009.
Hetherington AM, Woodward FI: The role of stomata in sensing and driving environmental change. Nature. 2003, 424 (6951): 901-908. 10.1038/nature01843.
Finkelstein RR, Gampala SS, Rock CD: Abscisic acid signaling in seeds and seedlings. The Plant cell. 2002, 14 (Suppl): S15-45.
Himmelbach A, Yang Y, Grill E: Relay and control of abscisic acid signaling. Current opinion in plant biology. 2003, 6 (5): 470-479. 10.1016/S1369-5266(03)00090-6.
Liu PF, Chang WC, Wang YK, Chang HY, Pan RL: Signaling pathways mediating the suppression of Arabidopsis thaliana Ku gene expression by abscisic acid. Biochim Biophys Acta. 2008, 1779 (3): 164-174.
Song SK, Clark SE: POL and related phosphatases are dosage-sensitive regulators of meristem and organ development in Arabidopsis. Developmental biology. 2005, 285 (1): 272-284. 10.1016/j.ydbio.2005.06.020.
Bailey TL, Elkan C: The value of prior knowledge in discovering motifs with MEME. Proceedings/International Conference on Intelligent Systems for Molecular Biology; ISMB. 1995, 3: 21-29.
Lescot M, Dehais P, Thijs G, Marchal K, Moreau Y, Peer Van de Y, Rouze P, Rombauts S: PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic acids research. 2002, 30 (1): 325-327. 10.1093/nar/30.1.325.
Higo K, Ugawa Y, Iwamoto M, Korenaga T: Plant cis-acting regulatory DNA elements (PLACE) database: 1999. Nucleic acids research. 1999, 27 (1): 297-300. 10.1093/nar/27.1.297.
Gomez-Porras JL, Riano-Pachon DM, Dreyer I, Mayer JE, Mueller-Roeber B: Genome-wide analysis of ABA-responsive elements ABRE and CE3 reveals divergent patterns in Arabidopsis and rice. BMC genomics. 2007, 8: 260-10.1186/1471-2164-8-260.
Sasaki T, Yano M, Kurata N, Yamamoto K: The Japanese Rice Genome Research Program. Genome research. 1996, 6 (8): 661-666. 10.1101/gr.6.8.661.
Sasaki T, Burr B: International Rice Genome Sequencing Project: the effort to completely sequence the rice genome. Current opinion in plant biology. 2000, 3 (2): 138-141. 10.1016/S1369-5266(99)00047-3.
Rice P, Longden I, Bleasby A: EMBOSS: the European Molecular Biology Open Software Suite. Trends Genet. 2000, 16 (6): 276-277. 10.1016/S0168-9525(00)02024-2.
We are grateful to Prof. Jingdong Ye from Central Florida University and Dr. Xinzheng Li from Shandong Agricultural University for their suggestions on writing this article. We are also thankful to the very helpful suggestions from the two reviewers of this manuscript. This work was supported by Shandong Agricultural University in China, the National Natural Science Foundation of China (Grant No. 30571308), the National Basic Research Program (Grant No. 2006CB1001006) of China and Millersville University of Pennsylvania in the US.
TX and DW carried out all the data analyses. TX performed the RT-PCR experiments. FN and YZ contributed the PP2C gene family background knowledge. YZ and CZ provided financial support for the study. TX and SZ jointly wrote the draft manuscript. SJ, SZ and FN double checked the presented data. XT, YZ, JE and CZ participated in the writing of the final version of the manuscript.
Tongtong Xue, Dong Wang contributed equally to this work.
Electronic supplementary material
Additional file 1:Table S1 Distribution of putative regulatory cis-elements in 5'-upstream regions of subfamily A members.(DOC 68 KB)
Additional file 2: Figure S1 Putative cis -elements distribution in upstream regulatory regions in PP2C subfamily A and D members from Arabidopsis and rice. 1000 bp genomic DNA sequences upstream of the first exon of each gene were extracted and the positions of the upstream regulatory regions are indicated. The relative positions of cis-elements are labeled with capital letters, which are annotated at the bottom. (DOC 957 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Xue, T., Wang, D., Zhang, S. et al. Genome-wide and expression analysis of protein phosphatase 2C in rice and Arabidopsis. BMC Genomics 9, 550 (2008). https://doi.org/10.1186/1471-2164-9-550