Skip to main content

RAWUL: A new ubiquitin-like domain in PRC1 Ring finger proteins that unveils putative plant and worm PRC1 orthologs

Abstract

Background

Polycomb group (PcG) proteins are a set of chromatin-modifying proteins that play a key role in epigenetic gene regulation. The PcG proteins form large multiprotein complexes with different activities. The two best-characterized PcG complexes are the PcG repressive complex 1 (PRC1) and 2 (PRC2) that respectively possess histone 2A lysine 119 E3 ubiquitin ligase and histone 3 lysine 27 methyltransferase activities. While PRC2-like complexes are conserved throughout the eukaryotic kingdoms, PRC1-like complexes have only been described in Drosophila and vertebrates. Since both complexes are required for the gene silencing mechanism in Drosophila and vertebrates, how PRC1 function is realized in organisms that apparently lack PRC1 such as plants, is so far unknown. In vertebrates, PRC1 includes three proteins, Ring1B, Ring1A, and Bmi-1 that form an E3 ubiquitin ligase complex. These PRC1 proteins have an N-terminally located Ring finger domain associated to a poorly characterized conserved C-terminal region.

Results

We obtained statistically significant evidences of sequence similarity between the C-terminal region of the PRC1 Ring finger proteins and the ubiquitin (Ubq)-like family proteins, thus defining a new Ubq-like domain, the RAWUL domain. In addition, our analysis revealed the existence of plant and worm proteins that display the conserved combination of a Ring finger domain at the N-terminus and a RAWUL domain at the C-terminus.

Conclusion

Analysis of the conserved domain architecture among PRC1 Ring finger proteins revealed the existence of long sought PRC1 protein orthologs in these organisms, suggesting the functional conservation of PRC1 throughout higher eukaryotes.

Background

Polycomb group (PcG) proteins are epigenetic gene regulators implicated in important cellular and developmental processes. In animals, PcG proteins were primarily known for their role in maintaining cell identity during the establishment of the body plan [1], but recently they have also been implicated in other important processes like stem cell self-renewal and cancer [2, 3]. In plants, the PcG proteins are involved in the repression of flowering during vegetative development [4, 5], seed development [6], and in the vernalization response [7]. Hence, PcG proteins play an essential role in the proper development of multicellular organisms.

Biochemical and genetic characterizations of PcG proteins have revealed that they exist in distinct multi-protein complexes [8, 9], of which the two best characterized are the PcG repressive complex (PRC)1 and PRC2. The core components of PRC2 are conserved throughout higher eukaryotes. The complex mediates histone H3 methylation at lysine 27 [10–12]. Similar complexes have been found in Caenorhabditis elegans [13, 14], vertebrates [10–12, 15, 16] and plants [4, 5].

Unlike PRC2, PRC1-like complexes have only been described in Drosophila and vertebrates. Drosophila PRC1 core complex contains Posterior sex combs (Psc), Polyhomeotic (Ph), Drosophila Ring finger protein 1 or Sex combs extra (dRing1/Sce) and Polycomb (Pc) [17, 18]. The complex can inhibit in vitro both chromatin remodeling by the human SWItch/Sucrose NonFermentable (hSWI/SNF) complex [19, 20] and transcription of a chromatin template by RNA polymerase II [21]. The human PRC1 complex is composed of B lymphoma Mo-MLV insertion region 1 homolog (Bmi-1), PH2, PC3, and the Ring finger proteins (Ring1A/RING1 and Ring1B/RING2), homologs of Drosophila Psc, Ph, Pc, and dRing, respectively. Human Bmi1, Ring1A and Ring1B proteins form an E3 ubiquitin ligase complex that mono-ubiquitinates lysine 119 of nucleosomal histone H2A (H2A-K119) [22, 23]. The H2A-K119 monoubiquitination is required for the PcG-mediated in vivo gene silencing, but its role in the mechanism is not yet known [23]. Drosophila PRC1 presumably possesses ubiquitin E3 ligase activity as well, since dRing was found to co-localize with ubiquitinated H2A at the Ultrabithorax (Ubx) promoter [22].

In Drosophila and vertebrates the PcG-mediated gene silencing mechanism requires the action of both PRC2 and PRC1. The apparent lack of homologs of PRC1 components in organisms such as plants or worms and the conflicting information on histone ubiquitination in these organisms led to the speculation of alternative PcG-mediated gene silencing mechanism in which other proteins undertake PRC1 functions [4, 24–29]. Ubiquitination at lysine 119 of H2A in human was observed within the consensus peptide sequence PKKT [22]. Arabidopsis genome contains 13 histone H2A genes (HTA1-13) [30], among them just HTA10 displays conservation of the PKKT consensus sequence. Although H2A ubiquitination has not been detected in Arabidopsis using western blot analysis so far [31], a recent report showed the presence of ubiquitinated H2A in maize [32]. On the other hand, ubiquitinated H2A has been detected in C. elegans [29], indicating that some proteins may be involved in H2A ubiquitination.

The catalytic subunit of the mammalian E3 ubiquitin ligase complex is Ring1B, but Ring1A and Bmi-1 are also required for the in vivo H2A ubiquitination [31]. Bmi-1, Ring1B and Ring1A contain a Ring finger domain located at the N-terminal region. Biochemical analyses showed that the N-terminal Ring finger domain of Ring1B is sufficient for the E3 ligase activity in vitro [22, 33, 34] and that Bmi-1 displays no detectable ubiquitin ligase activity but it greatly stimulates the E3-ligase activity of Ring1B. It was also described that in vitro Ring1A can replace Ring1B [34]; however, the in vivo function must be distinct, as most of the ubiquitinated H2A is depleted upon lost of Ring1B [31, 35]. Structural analysis revealed that Ring1B interacts with Bmi-1 via their Ring fingers [33, 34].

The PRC1 Ring finger proteins also share a conserved C-terminal region [36, 37] that seems to be implicated in the interaction with Ph and Pc [38–40]. To gain mechanistic insight into the assembly and the enzymatic activity of the PRC1 proteins we performed a detailed computational sequence analysis of the C-terminal region of these proteins. Our findings result in the identification of a new ubiquitin-like domain that unveiled PRC1 Ring finger-like proteins in plants and worm.

Results and discussion

Detection of a new ubiquitin-like domain in the PRC1 Ring finger proteins

To analyse the conserved C-terminal region of the Ring finger PRC1 proteins, we first performed BLAST searches [41] against different sequence databases located at UniProt [42], NBCI [43], ENSEMBL [44] and JGI [45], starting from the human Bmi-1. Alignments were generated with T-Coffee [46] and checked manually. Additional profile-based sequence searches were performed against the Uniref50 and Uniref 90 protein sequence databases [47] with the defined global hidden Markov models using HMMer [48]. Using extensive profile to sequence (HMMer) and profile to profile comparisons analyses [49] we identified statistically significant E-values of sequence similarity between the C-terminal conserved regions of PRC1 Ring finger proteins and the C-terminal region of the WDR48-p80 protein family (Figure 1 and 2). The WD repeat domain 48 or p80 (WDR48-p80) family are WD40 repeat-containing endosomal proteins found in all eukaryotes including yeast and are present in only one copy per organism (Figure 1 and 3). Except for the previously described interaction between the human WDR48-p80 protein and a tyrosine kinase interacting protein from Herpesvirus saimiri virus (Tip), the exact function of the WDR48-p80 family remains unknown but has been suggested to be related to endosome/lysosome traffic [50–52].

Figure 1
figure 1

Representative multiple alignment of the RAWUL domain. The coloring scheme indicates average BLOSUM62 score (correlated to amino acid conservation) in each alignment column: red (greater than 1.8), violet (between 1.8 and 1) and light yellow (between 1 and 0.2). Residues that are part of the hydrofobic core of the Ubq fold [50] are indicated by green ovals above the alignment. The limits of the domains are indicated by the residue positions on each side. The main families and subfamilies are indicated by coloured bars to the left of the alignment; green (PRC1), violet (WD48-p80), red (Bmi1/Mel18) and yellow (Ring1A/B). Structure based sequence alignment [86] of representative members of UBL superfamily were included: pdb_1gnu: human gamma-aminobutyric acid receptor type A receptor-associated protein (GABARAP) [84]; pdb_1m94: yeast Ubl modifier protein homologous to ubiquitin 1 (Hub1) [87]; pdb_1yqb: human ubiquilin 3 [88]; pdb_1mg8: Ubl domain from mouse Parkin protein [89]; pdb_1aar: bovine Ubiquitin [90]. Independent PHD secondary structure predictions [55] for PRC1 Ring and WDR48-p80 families are shown below the family sequences. Consensus secondary structure of the Ubiquitin superfamily is shown below the alignment [54]. Alpha-helices and beta-strands are indicated by green cylinders and arrows, respectively. The sequences are named with their Uniprot identifications, and also, if necessary, with their common species name: Human, Homo sapiens; Drome, Drosophila melanogaster; Caeel, Caenorhabditis elegans; Arath, Arabidopsis thaliana; Orysj, Oryza sativa; Crypv, Cryptosporidium parvum; Schja, Schistosoma japonicum; Aspfu, Aspergillus fumigatus; Neucr, Neurospora crassa; Kluma, Kluyveromyces marxianus; Debha, Debaryomyces hansenii; Schpo, Schizosaccharomyces pombe; Yeast, Saccharomyces cerevisiae. The "est" prefix identifies sequence corrections using consensus sequences manually reconstructed by assembling highly similar expressed sequence tags from identical species (conceptual translations). The numbers in blue within the alignment represent sequence insertions that are not included.

Figure 2
figure 2

HMM-HMM profile comparison E-Values between the RAWUL and Ubl domains. The numbers correspond to HMM-HMM profile comparisons E-values from global profile search results [49] that connect independently each family with the others. The arrows indicate the profile search direction.

Figure 3
figure 3

Dendrogram and Domain Architecture for RAWUL domain containing proteins. The sequences are named as described in Figure 1. Vitis vinifera (Vitvi) sequences were included to construct the dendrogram. The main families (PRC1 Ring and WDR48-p80) and subfamilies (Bmi1/Mel18 and Ring1A/B) are indicated by coloured dotted lines, as in Figure 1. Only a representative set of sequences was labelled. The scale bar shows the average number of amino acid substitutions per site (0.1). The main bootstrap values for PRC1 Ring family classification are indicated. Schematics of domain architecture are represented with the localization of the RAWUL, Ring domains and WD40 repeats according to this analysis, Pfam [91] and REP web servers [92], respectively.

The statistical significance of profile to profile comparisons was evaluated in terms of an E-value. Significant E-values connected all these sequences and reciprocal searches produced convergent results (Figure 2). The Bmi1/Mel18 subfamily profile search finds the Ring1A/1B subfamily with a 3.2e-11 E-value and the global profile of the PRC1 Ring family (including both Bmi1/Mel18 and Ring1A/1B subfamilies) detected the WDR48-p80 and the Ubl families with E-values of 3.7e-10 and 4.7e-07, respectively. The RAWUL HMM profile of the WDR48 family have less capability to detect remote protein homologies and localize the ubiquitin family with poor statistical significance (with an E value of 0.21), although RAWUL domain is more conserved inside the WDR48 family than in PRC1 Ring family. Nevertheless, the homology of the WDR48 RAWUL domain with the UBL family is statistically supported, when using as intermediate the RAWUL HMM profile of PRC1 Ring family.

This result defines a new conserved domain present in the PRC1 Ring finger and WDR48 families. In all of those proteins the domain is always observed in a C-terminal location (Figure 3), suggesting that this is somehow important for its function.

Statistically significant sequence similarity offered by profile-profile comparisons linked this domain with ubiquitin superfamily (Figure 1 and 2). Ubiquitin (Ubq) is a highly conserved 76 amino acid residue protein found in all eukaryotic cells and whose sequence is extremely well conserved from protozoa to vertebrates [53]. The overall topology of the Ubq-like fold is β1, β2, α1, β3, β4, α2, β5 (where β indicates beta strand and α alpha helix); however, there may be deviations from this common core [54].

Although those E-values are significant per se, we performed three additional analyses. First, secondary structure predictions were performed independently for the conserved C-terminal regions of PRC1 Ring finger and WDR48-p80 families [55] (Figure 1). These predictions showed good agreement with the Ubq and Ubq-like structures (Figure 1). Second, to investigate whether fold recognition analysis generated consistent results, we submitted the Ubq-like domain of the human Bmi-1 protein (Uniprot-id: BMI1_HUMAN, residues 131 to 226) to an independent fold assignment software (mGenThreader [56]). mGenthreader consistently recognised only Ubq-like proteins as possible templates (8 out of the 10 proposed hits are Ubq-like folds). This consistency in the threading output additionally supports the relationship between Ubq and the conserved PRC1 Ring C-terminal domain. Third, we generated a structural model (see Additional file 1) of the human Bmi-1 putative Ubq-like domain based on alignments including ubiquitin superfamily sequences of known structure (Figure 1). The model was evaluated using statistical mean force potential and found to have a Z-score of -6.5, in the upper range of good model for this length [57], strengthening the assignment of the Ubq fold.

Despite the low sequence identity between this domain and the Ubq (below 20%), the hydrophobic core and key secondary structural elements of the Ubq fold are conserved (Figure 1) [54]. However, the PRC1 Ring finger family displays a higher degree of sequence variability in the C-terminal Ubq-like domain than the WDR48-p80 family. The domain secondary structural elements are rich in charged residues and are separated by loops of variable length (Figure 1 and Additional file 2), which are likely to be the reasons why this domain escaped detection until now. The fact that the domain shows insertions of variable length is not in contradiction with a stable fold [58]; even more, small folds, such as Ubiquitin [54], Immunoglobulin [59] or SBDS [60], share a high degree of structural conservation in a scenario of high sequence divergence and long loop insertions.

The most characteristic features of this domain is the α2 that is a bit larger in its C-terminus compared to classical Ubq-like proteins and contains a conserved pair of Thr/Leu that might have a structural role in capping the alpha helix (Figure 1). Preliminary structural characterization of the C-terminal region of Ring1B showed that this domain folds independently and is a combination of alpha-helix and beta-sheet secondary structures [61], which is in agreement with our analysis.

In summary, considering the profile-profile comparisons E-values (Figure 2), secondary structure predictions (Figure 1), fold assignment detection and the results of the homology model building and evaluation (see Additional file 1), we are confident that this domain is a new member of the Ubq superfamily. Since we found this new Ubq-like domain associated to two different domain architectures, Ring-finger domain (in PRC1 Ring family) and WD40 repeats (in WDR48-p80 family), we named it RAWUL (R ing-finger A nd W D40 associated U biquitin-L ike) domain.

The ubiquitin related proteins fall into two separate classes [62]. Ubq-like proteins (UBLs; eg, SUMO, NEDD8, etc.) that function as modifiers in a manner analogous to that of ubiquitin and exist either in a free form or attached covalently to other proteins by their C-termini, and the ubiquitin-domain proteins (UDPs) that bear domains related to ubiquitin but are otherwise unrelated in sequence to each other (i.e. parkin, Rad23, DSK2). In contrast to the UBL modifiers, the UDPs are not conjugated to other proteins.

RAWUL domain lacks the characteristic C-terminal diglycine motif required for enzymatic conjugation of the Ubq domain to an acceptor lysine. Neither are conserved the three most common acceptor lysine residues involved in polyubiquitination (Figure 1, represented in ubiquitin (PDB:1aar) as K29, K48, and K63). Hence, it appears unlikely that this domain functions as a traditional UBL modifier. PRC1 Ring finger and WDR48-p80 proteins therefore could be included in the group of UDPs. There are two classes of UDPs as well [63] differentiated by their location and sequence similarity to ubiquitin of their Ubq-like domains.

The Ubq-like domain (UBD) proteins have the Ubq-like domain located at or close to the N-terminus of the protein and defined by a stretch of 45–80 residues with significant sequence homology to ubiquitin, while the Ubiquitin regulatory X (UBX) domain proteins have the Ubq-like domain at the C-terminus consisting of ~80 amino acids that shares common secondary structure organization with ubiquitin despite the lack of significant sequence homology [64, 65]. The UBD proteins have been shown to bind to ubiquitin-interacting motifs (UIMs), like parkin protein. Defective parkin protein is responsible for a common familial form of Parkinson's disease. Parkin encodes an E3 ubiquitin ligase that contains two Ring finger domains at its C-terminus and one Ubq-like domain at its N-terminus [65, 66]. Binding of the parkin ubiquitin-like domain to the Eps15 UIM is required for parkin-mediated Eps15 ubiquitination [66].

Conversely, only very few UBX domain proteins have been studied in detail and no general function for the UBX domain has yet emerged. Most of the UBX domain proteins identified so far can be grouped into four evolutionary conserved protein families represented by the human Fas-associated factor-1 (FAF1), p47, Y33K, and Rep8 proteins [67]. Recent reports showed that the UBX domains may act as general binding modules for p97 and/or p97 homologs, possibly representing a first common role for UBX domains [68]. RAWUL domain seems to fall into the UBX domain protein family since it is located at the C-terminus, comprise about 80 residues and the conserved hydrophobic core of the Ubq fold. However, as far as we know, UBX domains have not been found together with a Ring finger domain or a WD40 repeat domain.

In the case of the poorly characterized WDR48-p80 family [50–52], the association of WD40 repeats and RAWUL domain could point out to a E3-ligase related function, since WD40 repeats containing proteins are commonly implicated in protein-protein interaction [69] and these repeats have been recently described as key elements in E3 ubiquitin ligase complexes [70].

On the other side, our analysis revealed that the combination of a Ring finger domain at the N-terminus and a RAWUL domain at the C-terminus is a key feature to define the PRC1 Ring finger protein family, since no other protein presents this domain architecture. The Ring finger domains of the Ring (Ring1A and 1B) and Bmi1 proteins are responsible for the E3 ubiquitin ligase activity [33, 34], whereas the region corresponding to the RAWUL domain has been shown to mediate protein-protein interaction with Pc and Ph [38–40]. The identification of an Ubq-like domain in this region raise the possibility that these interactions could be establish through a yet uncharacterised UIM present in Pc and Ph. However, due to the structural similarities between the RAWUL domain and ubiquitin, it is also possible that RAWUL could act as an auto-inhibitory domain, by direct competition with the substrate, regulating specificities and/or ubiquitination capabilities of the E3 ligase complex, as previously suggested for Ubq-like domains fused to ubiquitin-specific proteases [71], or as a bridge, linking the proteasome and histone post-translational modifications, as suggested for the Ubq-like domain in the Ubiquitin plant homeodomain Ring finger (UHRF) protein family [72].

Identification of plant and worm PRC1 Ring finger protein orthologs

In the course of characterizing the PRC1 Ring finger proteins, we identified a set of proteins in Arabidopsis, Oryza sativa, Vitis vinifera and worms that share the PRC1 ring finger domain architecture, an N-terminal Ring finger domain and a RAWUL domain at their C-terminal region (Figure 1, 3 and 4). The preservation of domain architecture among these proteins indicates that these plant and worm proteins are potential orthologs of animal PRC1 Ring finger proteins.

Figure 4
figure 4

Representative full sequence alignment of the PRC1 Ring family. The coloring scheme indicates average BLOSUM62 score (correlated to amino acid conservation) in each alignment column: red (greater than 4), violet (between 4 and 1.5) and light yellow (between 1.5 and 0.5). The sequences are named as described in Figure 1. Lowercase characters represent parts of the sequence that are not evolutionarily conserved. In the PRC1 Ring family the N-terminal Ring domain (blue rectangle) is more conserved than the C-terminal RAWUL domain (green rectangle).

The Ring finger domain is highly abundant in plant proteins, with as many as 469 predicted Ring finger-containing proteins in the Arabidopsis genome [73]. The large number of potential Ring finger proteins in Arabidopsis suggests that target specific ubiquitination plays an important role in protein regulation in plants. However, the combination of the Ring finger domain and the RAWUL domain, which is a specific feature of the PRC1 Ring finger proteins, is found in only a few plant proteins. In Arabidopsis, there are four PRC1 Ring finger genes, At1g03770, At5g44280, At2g30580 and At3g23060 that codify the predicted proteins Q0WX00, estQ9FKW0, Q94AY3 and Q9LS86, respectively. There are full length cDNAs for these four genes with one open reading frame predicted for each (the Arabidopsis information resource (TAIR) database [74]. Online microarray expression data proof that these genes are expressed [74].

Outside of these two domains, the plant proteins have diverged significantly from their putative animal counterparts (Figure 4), but usually PcG protein orthologs do not display sequence homology outside the conserved domains, not only in evolutionary distant organisms but also among paralogs in a given organism [75]. Taking into account the fact that convergent evolution of domain architectures is rare [76], the common domain architecture presented by the plants and animal PRC1 Ring finger proteins (Figure 3) suggest a common evolutionary origin, raising the possibility that these plant proteins could be the PRC1 Ring finger orthologs, in turn implying that the PRC1-mediated histone H2A ubiquitination might also take place in plant.

Plants were thought to have no homologs of PRC1 components. However, two recent reports (25,26) shows that Arabidopsis Like Heterochromatin Protein1 (LHP1, also called Terminal Flower2, TFL2) is functionally similar to Pc, a subunit of PRC1 that recognizes and binds the H3K27 marks created by PRC2. Similarly, PRC1 Ring finger proteins were thought to be missing in plants; however, our analysis revealed the presence of putative orthologs to these proteins. All together these data suggest that PRC1 is also conserved in plants, in agreement with our analysis.

Based on sequence similarity and biochemical data, two groups of PRC1 Ring finger proteins have been defined [18, 36], the Bmi1 subfamily that includes vertebrate Bmi1 and Mel18 and Drosophila Psc and Suppressor 2 of zeste (Su(z)2), and the Ring subfamily that includes vertebrate Ring1A, Ring1B and Drosophila dRing. In order to identify potential orthologous relationships between these plant and animal homologous genes we conducted a phylogenetic analysis of this family. Due to the low degree of sequence conservation in the C-terminal domain RAWUL (see Figure 1), we add the N-terminal Ring domain, which is highly conserved between plants and animals (Figure 4), to the phylogenetic analysis of the PRC1 Ring family. The high confidence bootstrap value (997/1000) at the connecting node (Figure 3) indicates that both subfamilies (Bmi1/Mel18 and Ring1A/1B) are clearly present in plants (Figure 3). For instance, Arabidopsis Q0WX00 and estQ9FKW0 proteins are part of the Ring1A/1B subfamily, while Q94AY3 and Q9LS86 proteins are the orthologs of the Bmi1/Mel18 subfamily. The identification of orthologous relationships between the two groups including both plant and animal proteins suggests an early gene duplication event, giving rise to Bmi1/Mel18 and Ring1A/1B subfamilies prior to the plant-animal divergence. The plant and animal proteins therefore could share a common biochemical function.

Similarly, lower metazoans, such as nematodes, are also thought to lack PRC1 homologs. In C. elegans, PcG-mediated gene repression has been proposed to require a PRC2-like complex [77] and two novel proteins, Suppressor of Pal-1 (SOP-2) and Sop-2-related protein 1 (SOR-1) [27, 28]. Despite the lack of obvious sequence similarity, several conserved properties between PRC1 and the putative SOP-2/SOR-1 complex suggested a conservation of the mechanism [27, 28]. However, these proteins are not found in other organisms, including its closest related C. briggsae [28], weakening the case for an alternative PRC1-like complex. In addition, ubiquitinated H2A has been detected by western blot analysis in C. elegans [31], indicating that some proteins may be involved in H2A ubiquitination. Interestingly, we found one PRC1 Ring finger homolog in C. elegans and in C. briggsae (Figure 1 and 3). The function of this PRC1 Ring finger homolog is so far unknown. However, the conservation of the protein architecture, an N-terminal Ring finger and a C-terminal RAWUL domain, suggest that the function of this protein could be related to histone ubiquitination.

Conclusion

We detected statistically significant sequence similarity between the ubiquitin superfamily and the conserved C-terminal regions of PRC1 Ring and WDR48-p80 families, defining a new Ubq-like domain, the RAWUL domain. The identification of an Ubq-like domain in the PRC1 Ring family offers new experimental approaches aimed at elucidating their roles in important cellular processes, such us, stem cell self-renewal and cancer. In addition, characterizing the RAWUL domain allowed us to identify putative PRC1 Ring finger proteins that were thought to be missing in the plants and worms. The possibility that a PRC1-like complex is also involved in PcG-mediated gene silencing mechanism in these organisms is intriguing and opens new avenues in PcG investigation.

Methods

Sequence analysis

We first performed BLAST sequence similarity searches [41] against different protein sequence databases: UniProt [42], NBCI [43], ENSEMBL [44] and JGI [45]. For the sequence analysis we related distant protein families via profile searches using global hidden Markov models: HMMer [48] over the UniRef non redundant sequence database [47] and HHpred over pdb70 profile database [49]. The alignments were produced with the T-Coffee and HMMer [46, 48] using default parameters, slightly refined manually and viewed with the Belvu program [78]. Dendrogram were calculated with the neighbor-joining method [79] using ClustalW [80] and was edited with Treetool [81]. The stability of different branches with respect to different choices of subsets of residue positions was checked by bootstrapping experiments (1000 replicates) [82]. For profile-profile comparison we used HHpred software (version 1.2.0) downloaded from the authors [49]. Global HMM profiles were done for each family (PRC1 Ring and WDR48-p80) and subfamily (Bmi1/Mel18 and Ring1A/B) using non-redundant alignments of the RAWUL domain (90% of sequence identity) and removing gappy columns (>70%). We include these profiles into the pdb70 HMM-database [49] before profile-profile searches were done. All the profile-profile comparisons were done using default parameters and without known (at Ubl family) or predicted secondary structure information (hhsearch option – ssm 0).

Protein-structure predictions and modeling

Secondary-structure predictions were performed with Predict Protein server [55]. Fold-recognition analyses were performed with mGenThreader [56]. We generated structural models using Modeler [83] based on pdb structure 1gnu[84]. Models were evaluated using statistical mean force potential [53]. The figure of the model was generated with Pymol program [85].

References

  1. Ringrose L, Paro R: Epigenetic regulation of cellular memory by the Polycomb and Trithorax group proteins. Annu Rev Genet. 2004, 38: 413-43.

    Article  PubMed  CAS  Google Scholar 

  2. Martinez A, Cavalli G: The role of polycomb group proteins in cell cycle regulation during development. Cell Cycle. 2006, 5: 1189-97.

    Article  PubMed  CAS  Google Scholar 

  3. Valk-Lingbeek M, Bruggeman S, van Lohuizen M: Stem cells and cancer; the polycomb connection. Cell. 2004, 118: 409-18.

    Article  PubMed  CAS  Google Scholar 

  4. Pien S, Grossniklaus U: Polycomb group and trithorax group proteins in Arabidopsis. Biochim Biophys Acta. 2007, 1769: 375-82.

    Article  PubMed  CAS  Google Scholar 

  5. Calonje M, Sung ZR: Complexity beneath the silence. Curr Opin Plant Biol. 2006, 9: 530-537.

    Article  PubMed  CAS  Google Scholar 

  6. Kohler C, Makarevich G: Epigenetic mechanisms governing seed development in plants. EMBO Rep. 2006, 7: 1223-7.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Schmitz R, Amasino R: Vernalization: a model for investigating epigenetics and eukaryotic gene regulation in plants. Biochim Biophys Acta. 2007, 1769: 269-75.

    Article  PubMed  CAS  Google Scholar 

  8. Muller J, Kassis JA: Polycomb response elements and targeting of Polycomb group proteins in Drosophila. Curr Opin Genet Dev. 2006, 16: 476-484.

    Article  PubMed  Google Scholar 

  9. Breiling A, Sessa L, Orlando V: Biology of polycomb and trithorax group proteins. Int Rev Cytol. 2007, 258: 83-136.

    Article  PubMed  CAS  Google Scholar 

  10. Cao R, Wang L, Wang H, Xia L, Erdjument-Bromage H, Tempst P, Jones RS, Zhang Y: Role of histone H3 lysine 27 methylation in Polycomb-group silencing. Science. 2002, 298: 1039-43.

    Article  PubMed  CAS  Google Scholar 

  11. Czermin B, Melfi R, McCabe D, Seitz V, Imhof A, Pirrotta V: Drosophila enhancer of Zeste/ESC complexes have a histone H3 methyltransferase activity that marks chromosomal. Polycomb sites Cell. 2002, 111: 185-96.

    PubMed  CAS  Google Scholar 

  12. Muller J, Hart CM, Francis NJ, Vargas ML, Sengupta A, Wild B, Miller EL, O'Connor MB, Kingston RE, Simon JA: Histone methyltransferase activity of a Drosophila Polycomb group repressor complex. Cell. 2002, 111: 197-208.

    Article  PubMed  CAS  Google Scholar 

  13. Ketel CS, Andersen EF, Vargas ML, Suh J, Strome S, Simon JA: Subunit contributions to histone methyltransferase activities of fly and worm polycomb group complexes. Mol Cell Biol. 2005, 25: 6857-68.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  14. Bender LB, Cao R, Zhang Y, Strome S: The MES-2/MES-3/MES-6 complex and regulation of histone H3 methylation in C. elegans. Curr Biol. 2004, 14: 1639-43.

    Article  PubMed  CAS  Google Scholar 

  15. Beisel C, Imhof A, Greene J, Kremmer E, Sauer F: Histone methylation by the Drosophila epigenetic transcriptional regulator Ash1. Nature. 2002, 419: 857-62.

    Article  PubMed  CAS  Google Scholar 

  16. Kuzmichev A, Nishioka K, Erdjument-Bromage H, Tempst P, Reinberg D: Histone methyltransferase activity associated with a human multiprotein complex containing the Enhancer of Zeste protein. Genes Dev. 2002, 16: 2893-905.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  17. Fritsch C, Beuchle D, Muller J: Molecular and genetic analysis of the Polycomb group gene Sex combs extra/Ring in Drosophila. Mech Dev. 2003, 120: 949-54.

    Article  PubMed  CAS  Google Scholar 

  18. Gorfinkiel N, Fanti L, Melgar T, García E, Pimpinelli S, Guerrero I, Vidal M: The Drosophila Polycomb group gene Sex combs extra encodes the ortholog of mammalian Ring1 proteins. Mech Dev. 2004, 121: 449-62.

    Article  PubMed  CAS  Google Scholar 

  19. Francis NJ, Saurin AJ, Shao Z, Kingston RE: Reconstitution of a functional core polycomb repressive complex. Mol Cell. 2001, 8: 545-56.

    Article  PubMed  CAS  Google Scholar 

  20. Shao Z, Raible F, Mollaaghababa R, Guyon JR, Wu CT, Bender W, Kingston RE: Stabilization of chromatin structure by PRC1, a Polycomb complex. Cell. 1999, 98: 37-46.

    Article  PubMed  CAS  Google Scholar 

  21. King IF, Francis NJ, Kingston RE: Native and recombinant polycomb group complexes establish a selective block to template accessibility to repress transcription in vitro. Mol Cell Biol. 2002, 22: 7919-28.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  22. Wang H, Wang L, Erdjument-Bromage H, Vidal M, Tempst P, Jones RS, Zhang Y: Role of histone H2A ubiquitination in Polycomb silencing. Nature. 2004, 431: 873-8.

    Article  PubMed  CAS  Google Scholar 

  23. Cao R, Tsukada Y, Zhang Y: Role of Bmi-1 and Ring1A in H2A ubiquitylation and Hox gene silencing. Mol Cell. 2005, 20: 845-54.

    Article  PubMed  CAS  Google Scholar 

  24. Mylne JS, Barrett L, Tessadori F, Mesnage S, Johnson L, Bernatavichute YV, Jacobsen SE, Fransz P, Dean C: LHP1, the Arabidopsis homologue of HETEROCHROMATIN PROTEIN1, is required for epigenetic silencing of FLC. Proc Natl Acad Sci USA. 2006, 103: 5012-7.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  25. Zhang X, Germann S, Blus BJ, Khorasanizadeh S, Gaudin V, Jacobsen SE: The Arabidopsis LHP1 protein colocalizes with histone H3 Lys27 trimethylation. Nat Struct Mol Biol. 2007, 14: 869-871.

    Article  PubMed  CAS  Google Scholar 

  26. Turck F, Roudier F, Farrona S, Martin-Magniette ML, Guillaume E, Buisine N, Gagnot S, Martienssen RA, Coupland G, Colot V: Arabidopsis TFL2/LHP1 specifically associates with genes marked by trimethylation of histone H3 lysine 27. PLoS Genet. 2007, 3: e86-

    Article  PubMed  PubMed Central  Google Scholar 

  27. Zhang H, Christoforou A, Aravind L, Emmons SW, Heuvel van den S, Haber DA: The C. elegans Polycomb gene SOP-2 encodes an RNA binding protein. Mol Cell. 2004, 14: 841-7.

    Article  PubMed  CAS  Google Scholar 

  28. Zhang T, Sun Y, Tian E, Deng H, Zhang Y, Luo X, Cai Q, Wang H, Chai J, Zhang H: RNA-binding proteins SOP-2 and SOR-1 form a novel PcG-like complex in C. elegans. Development. 2006, 133: 1023-33.

    Article  PubMed  CAS  Google Scholar 

  29. Calonje M, Sanchez R, Chen LJ, Sung ZR: EMBRYONIC FLOWER1 participates in the Polycomb group-mediated AG gene silencing in Arabidopsis. Plant Cell. 2008, 20: 277-91.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  30. Yi H, Sardesai N, Fujinuma T, Chan CW, Veena , Gelvin SB: Constitutive expression exposes functional redundancy between the Arabidopsis histone H2A gene HTA1 and other H2A gene family members. Plant Cell. 2006, 18: 1575-89.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  31. de Napoles M, Mermoud JE, Wakao R, Tang YA, Endoh M, Appanah R, Nesterova TB, Silva J, Otte AP, Vidal M, Koseki H, Brockdorff N: Polycomb group proteins Ring1A/B link ubiquitylation of histone H2A to heritable gene silencing and X inactivation. Dev Cell. 2004, 7: 663-76.

    Article  PubMed  CAS  Google Scholar 

  32. Alatzas A, Foundouli A: Distribution of ubiquitinated histone H2A during plant cell differentiation in maize root and dedifferentiation in callus culture. Plant Science. 2006, 171: 481-487.

    Article  PubMed  CAS  Google Scholar 

  33. Li Z, Cao R, Wang M, Myers MP, Zhang Y, Xu RM: Structure of a Bmi-1-Ring1B polycomb group ubiquitin ligase complex. J Biol Chem. 2006, 281: 20643-9.

    Article  PubMed  CAS  Google Scholar 

  34. Buchwald G, Stoop van der P, Weichenrieder O, Perrakis A, van Lohuizen M, Sixma TK: Structure and E3-ligase activity of the Ring-Ring complex of polycomb proteins Bmi1 and Ring1b. Embo J. 2006, 25: 2465-74.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  35. Fang J, Chen T, Chadwick B, Li E, Zhang Y: Ring1b-mediated H2A ubiquitination associates with inactive X chromosomes and is involved in initiation of X inactivation. J Biol Chem. 2004, 279: 52812-5.

    Article  PubMed  CAS  Google Scholar 

  36. Brunk BP, Martin EC, Adler PN: Drosophila genes Posterior Sex Combs and Suppressor two of zeste encode proteins with homology to the murine bmi-1 oncogene. Nature. 1991, 353: 351-3.

    Article  PubMed  CAS  Google Scholar 

  37. van Lohuizen M, Frasch M, Wientjens E, Berns A: Sequence similarity between the mammalian bmi-1 proto-oncogene and the Drosophila regulatory genes Psc and Su(z)2. Nature. 1991, 353: 353-5.

    Article  PubMed  CAS  Google Scholar 

  38. Kyba M, Brock HW: The Drosophila polycomb group protein Psc contacts ph and Pc through specific conserved domains. Mol Cell Biol. 1998, 18: 2712-20.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  39. Alkema MJ, Bronk M, Verhoeven E, Otte A, van 't Veer LJ, Berns A, van Lohuizen M: Identification of Bmi1-interacting proteins as constituents of a multimeric mammalian polycomb complex. Genes Dev. 1997, 11: 226-40.

    Article  PubMed  CAS  Google Scholar 

  40. Satijn DP, Otte AP: Polycomb group protein complexes: do different complexes regulate distinct target genes?. Biochim Biophys Acta. 1999, 1447: 1-16.

    Article  PubMed  CAS  Google Scholar 

  41. Altschul S, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-402.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  42. UniProt. [http://www.ebi.uniprot.org/]

  43. NBCI. [http://www.ncbi.nlm.nih.gov]

  44. ENSEMBL. [http://www.ensembl.org]

  45. JGI. [http://www.jgi.doe.gov/]

  46. Notredame C, Higgins D, Heringa J: T-Coffee: A novel method for fast and accurate multiple sequence alignment. J Mol Biol. 2000, 302: 205-17.

    Article  PubMed  CAS  Google Scholar 

  47. Wu CH, Apweiler R, Bairoch A, Natale DA, Barker WC, Boeckmann B, Ferro S, Gasteiger E, Huang H, Lopez R, Magrane M, Martin MJ, Mazumder R, O'Donovan C, Redaschi N, Suzek B: The Universal Protein Resource (UniProt): an expanding universe of protein information. Nucleic Acids Res. 2006, 34: D187-91.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  48. Eddy SR: Hidden Markov models. Curr Opin Struct Biol. 1996, 6: 361-365.

    Article  PubMed  CAS  Google Scholar 

  49. Soding J: Protein homology detection by HMM-HMM comparison. Bioinformatics. 2005, 21: 951-60.

    Article  PubMed  Google Scholar 

  50. Cho N, Kingston D, Chang H, Kwon EK, Kim JM, Lee JH, Chu H, Choi MS, Kim IS, Jung JU: Association of herpesvirus saimiri tip with lipid raft is essential for downregulation of T-cell receptor and CD4 coreceptor. J Virol. 2006, 80: 108-18.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  51. Park J, Lee BS, Choi JK, Means RE, Choe J, Jung JU: Herpesviral protein targets a cellular WD repeat endosomal protein to downregulate T lymphocyte receptor expression. Immunity. 2002, 17: 221-33.

    Article  PubMed  CAS  Google Scholar 

  52. Park J, Cho NH, Choi JK, Feng P, Choe J, Jung JU: Distinct roles of cellular Lck and p80 proteins in herpesvirus saimiri Tip function on lipid rafts. J Virol. 2003, 77: 9041-51.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  53. Vijay-Kumar S, Bugg CE, Wilkinson KD, Vierstra RD, Hatfield PM, Cook WJ: Comparison of the three-dimensional structures of human, yeast and oat ubiquitin. J Biol Chem. 1987, 262: 6396-9.

    PubMed  CAS  Google Scholar 

  54. Kiel C, Serrano L: The ubiquitin domain superfold: structure-based sequence alignments and characterization of binding epitopes. J Mol Biol. 2006, 355: 821-44.

    Article  PubMed  CAS  Google Scholar 

  55. Rost B: PHD: predicting one-dimensional protein structure by profile-based neural networks. Methods Enzymol. 1996, 266: 525-539.

    Article  PubMed  CAS  Google Scholar 

  56. McGuffin L, Jones D: Benchmarking secondary structure prediction for fold recognition. Proteins. 2003, 52: 166-75.

    Article  PubMed  CAS  Google Scholar 

  57. Melo F, Sanchez R, Sali A: Statistical potentials for fold assessment. Protein Sci. 2002, 11: 430-48.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  58. Bork P, Holm L, Sander C: The immunoglobulin fold. Structural classification, sequence patterns and common core. J Mol Biol. 1994, 242: 309-20.

    PubMed  CAS  Google Scholar 

  59. Scalley-Kim M, Minard P, Baker D: Low free energy cost of very long loop insertions in proteins. Protein Sci. 2003, 12: 197-206.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  60. Czypionka A, de los Paños OR, Mateu MG, Barrera FN, Hurtado-Gómez E, Gómez J, Vidal M, Neira JL: The isolated C-Terminal domain of Ring1B is a dimer made od stable, well-structured monomers. Biochemistry. 2007, 46: 12764-12776.

    Article  PubMed  CAS  Google Scholar 

  61. Jentsch S, Pyrowolakis G: Ubiquitin and its kin: how close are the family ties?. Trends Cell Biol. 2000, 10: 335-42.

    Article  PubMed  CAS  Google Scholar 

  62. Buchberger A: From UBA to UBX: new words in the ubiquitin vocabulary. Trends Cell Biol. 2002, 12: 216-21.

    Article  PubMed  CAS  Google Scholar 

  63. Buchberger A, Howard M, Proctor M, Bycroft M: The UBX domain: a widespread ubiquitin-like module. J Mol Biol. 2001, 307: 17-24.

    Article  PubMed  CAS  Google Scholar 

  64. Yuan X, Shaw A, Zhang X, Kondo H, Lally J, Freemont PS, Matthews S: Solution structure and interaction surface of the C-terminal domain from p47: a major p97-cofactor involved in SNARE disassembly. J Mol Biol. 2001, 311: 255-263.

    Article  PubMed  CAS  Google Scholar 

  65. Safadi SS, Shaw GS: A Disease State Mutation Unfolds the Parkin Ubiquitin-like Domain. Biochemistry. 2007, 46: 14162-14169.

    Article  PubMed  CAS  Google Scholar 

  66. Fallon L, Bélanger CM, Corera AT, Kontogiannea M, Regan-Klapisz E, Moreau F, Voortman J, Haber M, Rouleau G, Thorarinsdottir T, Brice A, En Henegouwen van Bergen PM, Fon EA: A regulated interaction with the UIM protein Eps15 implicates parkin in EGF receptor trafficking and PI(3)K-Akt signalling. Nat Cell Biol. 2006, 8: 834-42.

    Article  PubMed  CAS  Google Scholar 

  67. Dreveny I, Kondo H, Uchiyama K, Shaw A, Zhang X, Freemont PS: Structural basis of the interaction between the AAA ATPase p97/VCP and its adaptor protein p47. EMBO J. 2004, 23: 1030-1039.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  68. Park S, Rancour DM, Bednarek SY: Protein Domain-Domain Interactions and Requirements for the Negative Regulation of Arabidopsis CDC48/p97 by the Plant Ubiquitin Regulatory X (UBX) Domain-containing Protein, PUX1. J Biol Chem. 2007, 282: 5217-5224.

    Article  PubMed  CAS  Google Scholar 

  69. Andrade M, Perez-Iratxeta C, Ponting C: Protein repeats: structures, functions, and evolution. J Struct Biol. 2001, 134: 117-131.

    Article  PubMed  CAS  Google Scholar 

  70. Higa L, Zhang H: Stealing the spotlight: CUL4-DDB1 ubiquitin ligase docks WD40-repeat proteins to destroy. Cell Div. 2007, 6: 2-5.

    Google Scholar 

  71. Zhu X, Menard R, Sulea T: High incidence of ubiquitin-like domains in human ubiquitin-specific proteases. Proteins. 2007, 69: 1-7.

    Article  PubMed  CAS  Google Scholar 

  72. Bronner C, Achour M, Arima Y, Chataigneau T, Saya H, Schini-Kerth VB: The UHRF family: Oncogenes that are drugable targets for cancer therapy in the near future?. Pharmacol Ther. 2007, 115: 419-434.

    Article  PubMed  CAS  Google Scholar 

  73. Stone SL, Hauksdóttir H, Troy A, Herschleb J, Kraft E, Callis J: Functional analysis of the RING-type ubiquitin ligase family of Arabidopsis. Plant Physiol. 2005, 137: 13-30.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  74. TAIR. [http://www.arabidopsis.org/]

  75. Whitcomb S, Basu A, Allis C, Bernstein E: Polycomb Group proteins: an evolutionary perspective. Trends Genet. 2007, 23: 494-502.

    Article  PubMed  CAS  Google Scholar 

  76. Gough J: Convergent evolution of domain architectures (is rare). Bioinformatics. 2005, 21: 1464-71.

    Article  PubMed  CAS  Google Scholar 

  77. Fong Y, Bender L, Wang W, Strome S: Regulation of the different chromatin states of autosomes and X chromosomes in the germ line of C. elegans. Science. 2002, 296: 2235-2238.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  78. Sonnhammer EL, Hollich V: Scoredist: a simple and robust protein sequence distance estimator. BMC Bioinformatics. 2005, 6: 108-

    Article  PubMed  PubMed Central  Google Scholar 

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

    PubMed  CAS  Google Scholar 

  80. Higgins DG, Thompson JD, Gibson TJ: Using CLUSTAL for multiple sequence alignments. Methods Enzymol. 1996, 266: 383-402.

    Article  PubMed  CAS  Google Scholar 

  81. Maciukenas M: Treetool 2.0.2. 1994, [http://iubio.bio.indiana.edu/soft/molbio/unix/treetool]

    Google Scholar 

  82. Felsenstein J: Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985, 39: 783-791.

    Article  Google Scholar 

  83. Sali A, Blundell TL: Comparative protein modelling by satisfaction of spatial restraints. J Mol Biol. 1993, 234: 779-815.

    Article  PubMed  CAS  Google Scholar 

  84. Knight D, Harris R, McAlister MS, Phelan JP, Geddes S, Moss SJ, Driscoll PC, Keep NH: The X-ray crystal structure and putative ligand-derived peptide binding properties of gamma-aminobutyric acid receptor type A receptor-associated protein. J Biol Chem. 2002, 277: 5556-61.

    Article  PubMed  CAS  Google Scholar 

  85. DeLano WL: Unraveling hot spots in binding interfaces: progress and challenges. Curr Opin Struct Biol. 2002, 12: 14-20.

    Article  PubMed  CAS  Google Scholar 

  86. Holm L, Sander C: Dali: a network tool for protein structure comparison. Trends Biochem Sci. 1995, 20: 478-80.

    Article  PubMed  CAS  Google Scholar 

  87. Ramelot T, Cort JR, Yee AA, Semesi A, Edwards AM, Arrowsmith CH, Kennedy MA: Solution structure of the yeast ubiquitin-like modifier protein Hub1. J Struct Funct Genomics. 2003, 4: 25-30.

    Article  PubMed  CAS  Google Scholar 

  88. Mah A, Perry G, Smith M, Monteiro M: Identification of ubiquilin, a novel presenilin interactor that increases presenilin protein accumulation. J Cell Biol. 2000, 151: 847-62.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  89. Tashiro M, Okubo S, Shimotakahara S, Hatanaka H, Yasuda H, Kainosho M, Yokoyama S, Shindo H: NMR structure of ubiquitin-like domain in PARKIN: gene product of familial Parkinson's disease. J Biomol NMR. 2003, 25: 153-6.

    Article  PubMed  CAS  Google Scholar 

  90. Cook W, Jeffrey L, Carson M, Chen Z, Pickart C: Structure of a diubiquitin conjugate and a model for interaction with ubiquitin conjugating enzyme (E2). J Biol Chem. 1992, 267: 16467-71.

    PubMed  CAS  Google Scholar 

  91. Bateman A, Coin L, Durbin R, Finn RD, Hollich V, Griffiths-Jones S, Khanna A, Marshall M, Moxon S, Sonnhammer EL, Studholme DJ, Yeats C, Eddy SR: The Pfam protein families database. Nucleic Acids Res. 2004, 32: D138-41.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  92. Andrade MA, Ponting CP, Gibson TJ, Bork P: Protein repeats: structures, functions, and evolution. J Mol Biol. 2000, 298: 521-37.

    Article  PubMed  CAS  Google Scholar 

Download references

Acknowledgements

This work is supported by a Deutsche Forschungsgemeinschaft (DFG) grant (KO2302/8-1). ZRS is supported by a National Science Fundation (NSF) grant (IBN 0236399) and a United States Department of Agriculture (USDA) grant (03-35301-13244). DD is supported by the EU grant 3D-Repertoire (LSHG-CT-2005-5120828). We thank Kimmen Sjolander (Departments of bioengineering and Plant and Microbial Biology, University of California, Berkeley, CA, USA) for helpful suggestions of the manuscript.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Myriam Calonje.

Additional information

Authors' contributions

LSP carried out the database search, sequence analysis and drafted the manuscript. DD carried out the protein structure prediction and modelling, and helped to draft the manuscript. ZRS participated in the coordination and helped to draft the manuscript. MC conceived of the study, and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

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

Sanchez-Pulido, L., Devos, D., Sung, Z.R. et al. RAWUL: A new ubiquitin-like domain in PRC1 Ring finger proteins that unveils putative plant and worm PRC1 orthologs. BMC Genomics 9, 308 (2008). https://doi.org/10.1186/1471-2164-9-308

Download citation

  • Received:

  • Accepted:

  • Published:

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

Keywords