- Research article
- Open Access
Family expansion and gene rearrangements contributed to the functional specialization of PRDM genes in vertebrates
BMC Evolutionary Biology volume 7, Article number: 187 (2007)
Progressive diversification of paralogs after gene expansion is essential to increase their functional specialization. However, mode and tempo of this divergence remain mostly unclear. Here we report the comparative analysis of PRDM genes, a family of putative transcriptional regulators involved in human tumorigenesis.
Our analysis assessed that the PRDM genes originated in metazoans, expanded in vertebrates and further duplicated in primates. We experimentally showed that fast-evolving paralogs are poorly expressed, and that the most recent duplicates, such as primate-specific PRDM7, acquire tissue-specificity. PRDM7 underwent major structural rearrangements that decreased the number of encoded Zn-Fingers and modified gene splicing. Through internal duplication and activation of a non-canonical splice site (GC-AG), PRDM7 can acquire a novel intron. We also detected an alternative isoform that can retain the intron in the mature transcript and that is predominantly expressed in human melanocytes.
Our findings show that (a) molecular evolution of paralogs correlates with their expression pattern; (b) gene diversification is obtained through massive genomic rearrangements; and (c) splicing modification contributes to the functional specialization of novel genes.
The expansion of gene families in particular evolutionary groups as well as in single lineages significantly contributes to different aspects of genome innovation. However, not all classes of genes undergo equal rate of expansion. Genes controlling basic cellular functions, such as DNA modification, cytoskeletal organization, protein and RNA metabolism are remarkably conserved both in sequence and copy number across large evolutionary distances [1, 2]. Other functional categories progressively expand during evolution, with mode and tempo of this expansion largely depending on the effects on organism fitness. Genes coding for structural proteins and those controlling pathogen and stress response usually undergo lineage or organism specific expansions [3, 4]. In these cases environment-dependent regulation of gene dosage and expression is required. Genes coding for adhesion molecules, extracellular signalling or those involved in gene regulation expand preferentially during transitions towards higher complexity, defined as increase in cell types . The accretion of Zn-Fingers in vertebrates is among the most positively correlated with increase in complexity . The reason for the massive expansion of Zn-Fingers resides in the ability to bind DNA and regulate gene transcription in a tissue or cell-specific fashion. How this specificity is reached after gene expansion remains still an open question. According to the widely accepted model of evolution by gene duplication , in the first phases of their life novel genes are particularly free to change due to relaxed pressure acting on them. The most likely outcome is the acquisition of harmful modifications leading to loss of function (nonfunctionalization). However, in very rare cases, novel genes can acquire advantageous mutations differentiating their function from that of the ancestor (neofunctionalization). In a third scenario, needed to explain the abundance of gene duplications and lineage-specific expansions, mutations may occur in both genes, ancestor and duplicate, that specialize in complementary functions (subfunctionalization) [7–9].
Here we report a comparative analysis of the PRDM genes in vertebrates with the aim of connecting the expansion of this family of transcriptional regulators to their progressive functional specialization. The results of this study also contribute to unravel the general principles that drive gene family evolution. PRDM genes code for proteins controlling critical aspects of cell integrity, spanning from cell commitment and differentiation  to cell growth and apoptosis [11–13]. These genes also play a role in human cancer, where they mainly act as tumor suppressors . PRDM proteins share a characteristic domain organization with an N-terminal PR domain followed by a variable number of Zn-Finger repeats. The only exception is PRDM11, for which no Zn-Fingers are detectable (Figure 1A). The PR (P RDI-BF1 and R IZ) domain is 20–30% identical to the SET module, which is directly responsible for the catalytic activity of several histone lysine-methyltransferases . Some members of the PRDM family show intrinsic methyltransferase activity [16, 17], while others play an indirect role by recruiting chromatin remodelling enzymes [18, 19]. The number and features of PRDM genes in other species besides human are poorly known. Except for mouse, most PRDM orthologs in other organisms remain uncharacterized. Such lack of information has so far prevented from clarifying the reciprocal evolutionary relationships between different members of the PRDM family.
We started our analysis by searching for orthologs of the 17 human PRDM genes in representatives of fully sequenced eukaryotic genomes. We assessed that, after their appearance in metazoans, PRDM genes expanded in vertebrates and further duplicated in primates. By comparing the evolutionary features of PRDM genes with their expression in human tissues, we detected a tendency of newer genes to acquire tissue specificity. Slowly evolving paralogs also produce more abundant transcripts. We found that the massive structural rearrangements of primate-specific PRDM7 gene greatly contributed to modify gene function. In particular, the duplication of an internal segment conferred peculiar splicing pattern to the new gene, with the gain of a novel intron or, alternatively, its retention in the mature RNA.
Results and discussion
PRDM genes expanded in vertebrates and further duplicated in primates
By searching for PRDM orthologs in a variety of eukaryotic genomes, we did not identify any PRDM gene either in fungi or in plants, confirming that they first appeared in metazoans. In invertebrates, the orthology assignment yielded 2 genes for nematodes and 3 genes for arthropods (Figure 1B and Additional file 1). The additional sequence found in insects corresponds to Hamlet, a transcription factor involved in dendrite morphogenesis , and acting as negative regulator of S-phase entry in drosophila . Interestingly, also vertebrates PRDM2 , PRDM5 , and PRDM9  are involved in cell cycle regulation. Among vertebrates, we identified 17 putative PRDM orthologs in primates and 16 putative orthologs in rodents, birds and amphibians (Figure 1B). In fugu we only detected 15 PRDM genes, as we could not identify the PRDM17 ortholog possibly due to the presence of unresolved regions in the genome assembly.
Our results highlighted the expansion of the PRDM family in vertebrates. We tried to assess a possible timing for this expansion by checking for putative PRDM orthologs in protochordates and other invertebrates. We were able to identify only 2 PRDM genes in C. intestinalis (Figure 1B), likely because of the poor quality of the ciona genome assembly. The retrieval of 7 putative PRDM orthologs in the echinoderm S. purpuratus indicated that at least a partial expansion of the family preceded the last common ancestor of vertebrates (Figure 1B). Our data also showed that PRDM gene expansion continued in primates, as one single rodent co-ortholog corresponds to the primate in-paralog pair formed by PRDM7 and PRDM9 (see Additional file 1).
Expression pattern of PRDM genes reflects their evolutionary relationships
We rebuilt the phylogenetic relationships between PRDM genes based on the multiple alignment of the PR-domain of all PRDM proteins identified in vertebrates (Figure 2A). The trees obtained by applying both Maximum Likelihood and Bayesian methods are topologically comparable and well supported, the only exception being the internal relationships within the subfamily composed of PRDM7/9, PRDM11 and PRDM17. Overall, five subfamilies are clearly identifiable and the most divergent member is PRDM17. Despite the large variability in terms of gene structure and size among human PRDM paralogs, genes lying in sister branches of the tree maintain similar gene organization (Figure 2A). Such conservation partially reflects similar splicing patterns, as in the case of the subfamily composed of PRDM2, 3, 16. All these genes can code for both a full-length protein and an isoform without the PR-domain, which is often overexpressed in cancer [22, 23].
We observed interesting correlations between gene divergence and expression patterns. In human, all PRDM genes are generally expressed at very low levels, as assessed by qPCR (see Additional file 2). However, substantial differences are detectable among the 17 paralogs in terms of breadth and levels of expression. By ordering the PRDM proteins according to their evolutionary speed, we observed that less divergent paralogs are more expressed than faster-evolving genes (Figure 2B). This result assesses a negative correlation between rate of sequence evolution and gene expression levels (Pearson correlation, r = -0.705, P = 0.01). In addition, some paralogs, such as PRDM12, 13, 14, 7 and 9 acquire tissue specificity, since the corresponding transcripts are detectable only in very few of the analyzed samples (see Additional file 2). One of the youngest members of the family, PRDM7, shows an expression pattern mainly restricted to melanocytes, suggesting a progressive specialization and/or tighter regulation of its functions.
PRDM7 arose from PRDM9 duplication and acquired a peculiar splicing pattern
To rebuild the possible rearrangements following primate-specific duplication, we compared genetic and genomic features of the primate PRDM7 and PRDM9 genes with those of their co-ortholog PRDM7/9 in non-primates. We observed higher conservation of the genomic block around PRDM9 (Figure 3A) and higher similarity of human PRDM9 gene structure and protein domain organization with the non-primate co-ortholog (Figure 3B). These results suggest that PRDM9 likely retained the features of the ancestral locus, while PRDM7 acquired primate-specific features.
After segmental duplication, the new locus experienced massive rearrangements that significantly modified the structure of the resulting PRDM7 gene (Figure 3B). In human, the last exon coding for the Zn-Fingers underwent partial deletion and the resulting protein bears only 4 repeats instead of 14 (Figures 1A and 3B). In addition, an 89-nucleotide long segment within ancestral exon 3, which codes for the PR-domain C-terminal part, tandemly duplicated (Figure 4A). This led to the acquisition of a complex pattern of splicing variants (Figure 3C, see below). Interestingly, the duplication of the same 89-nucleotides occurred 8 times in the chimp genome (Figure 3B). However, the possible effect of this repetition on gene splicing is unknown since no chimp mRNA is available.
For human PRDM9, both a shorter (AF275816) and a longer (DQ388610) isoform have been reported. In the latter, the corresponding protein add a N-terminal KRAB domain to the usual PRDM architecture (Figure 3C). Meisetz, the mouse ortholog of the long PRDM9 isoform, is able to activate the progression into meiosis through the trimethylation of the lysine 4 on histone H3 [16, 24]. Recently, the KRAB domain of Meisetz has been proposed to be the ancestor of all modern KRAB domains, being traceable in the last common ancestor of vertebrates . In principle the genomic locus of PRDM7 could also code for a long isoform with the KRAB domain, but all mRNA sequences so far available corresponded to the short variant. One of them (AF274347) resembles the sequence of the short form of PRDM9 implying that the duplication inside exon 3 is coupled with an event of intron gain (Figures 3C and 4A). By using specific primers, we cloned the long isoform of PRDM7 and found that it can retain the 89-nucleotide long duplicon within the sequence of the mature mRNA (Figures 3C and 4A). This retention introduces a frameshift, and the resulting protein replaces the Zn-Fingers with an alternative C-terminal region (AM690991, Figure 3C). To assess the frequency of the two PRDM7 isoforms in human tissues, we amplified the corresponding region (Figure 4A) in a variety of RNA samples from normal tissues and cancer cell lines. In most cases, two bands of 289 bp and 200 bp are detectable (Figure 4B). Sequence analysis of the RT-PCR products confirmed that the higher band corresponds to the PRDM7-specific amplicon with retention of the internal duplication. The lower band could correspond to either PRDM7 or PRDM9 since the two genes have identical sequences in this region. The widespread detection of the isoform with the retained duplicon suggests that this might be a general mechanism to produce an alternative PRDM7 protein without Zn-Fingers. Interestingly, in mouse two isoforms of the Meisetz gene generated by alternative splicing also lack the Zn-Finger repeats and code for a protein with only the KRAB and the PR domains . By using several combinations of primers, we identified other long isoforms of PRDM7 with alternative patterns of exon retention (AM690992 and AM690993, Figure 3C), but we did not retrieve the long isoform coding for the Zn-Fingers. One explanation could be that the long isoform of PRDM7 with the Zn-Fingers is not produced. Another possibility is that it is extremely rare, being detectable only in specific tissues and/or cellular phases. This is the case of Meisetz, which is only expressed in fetal germ lines entering meiotic prophase .
The PRDM gene family of putative tumor suppressors constitutes an interesting example of vertebrate-specific expansion of transcriptional regulators. Except for PRDM11, all PRDM genes code for a variable number of C2H2 Zn-Finger repeats and, although the precise molecular function is mostly unknown, the ability to bind DNA seems to be a common feature. In addition, some of the PRDM genes control transcription through either direct or indirect histone methylation [16, 18, 19, 26]. Thus, rebuilding the progressive specialization of the PRDM ancestral function in different cell types and/or in diverse stages of development may yield general conclusions concerning paralog differentiation of transcriptional regulators in vertebrates.
PRDM genes show an inverse relationship between sequence divergence and expression level. Paralogs that evolve slowly, such as PRDM8, 2, and 3 tend to be more expressed than faster-evolving genes like PRDM15, 10 and 17 (Figure 2B). Although there is not a straightforward biological explanation for it [27, 28], this inverse correlation reflects an asymmetrical evolution of PRDM paralogs and is compatible with their progressive sub- and neo-functionalization [7, 8]. Both slowly (e.g. PRDM2) and fast-evolving (e.g. PRDM6) paralogs show methyltransferase activity, which implies the repartition of the putative ancestral function among several duplicates. Signs of neo-functionalization are also detectable, such as the different roles of PR domain in histone methylation and the specific involvement of PRDM9 in meiotic progression [16, 24]. The acquisition of novel functions is particularly evident in the primate-specific gene. As most newborn genes [29, 30], PRDM7 underwent remarkable structural rearrangements. Some modifications, such as the reduction of encoded Zn-Fingers, are recurrent within the PRDM family (Figure 1A), while the internal duplication of ancestral exon 3 is unusual for at least two reasons. First, the region undergoing duplication codes for the C-terminal part of the PR-domain, which is the only module always present in PRDM proteins. Usually, such well-conserved regions are under strong selective constraints that prevent harmful modifications. Second, the ancestral 89-nucleotide long segment does not contain any repetitive element to justify duplication through replication slippage and/or non-homologous recombination. Thus, the internal duplication seems to reflect adaptation of PRDM7 to novel functions. The effect of the duplication on PRDM7 splicing may constitute the first evidence for such adaptation. The duplication of a cryptic splice site (AG-GC) in the ancestral sequence results in the acquisition of a non-canonical splice site (GC-AG). This site can activate the splicing of an 89-nucleotide long intron partly deriving from the ancestral sequence and partly from the duplicon (Figure 4A). Such mechanism for intron gain was theoretically described by Rogers in 1989 , but has so far been observed very rarely [32, 33] and never in human genes. The resulting mature mRNA encodes for a protein with the usual domain architecture of PRDM genes (PR domain and Zn-Fingers). We found that the duplicated segment can be alternatively retained into the mature mRNA and eventually encodes for a frameshifted protein without Zn-Fingers (Figure 3C). Usually, the introduction of a frameshift is a dangerous event for a gene, and indeed mechanisms such as nonsense-mediated decay have evolved to prevent the generation of harmful proteins . However, a few cases of genes duplicated in vertebrates that use frameshift as a mechanism to diversify their function have been previously reported . The majority of these frameshifts originated through alternative splicing able to activate a different reading frame . In the case of PRDM7, instead, the frameshift arose through the duplication of an internal sequence. As for the other reported cases , also the mRNA transcribed by the PRDM7 gene is not a typical target for mRNA decay. It is around 2000 nucleotide long and can code for a protein with KRAB and PR domains and a C-terminal region of around 100 residues before encountering a stop codon (Figure 3C). The isoform with the retention is detectable in most of the tested samples, and it is particularly abundant in melanocytes and melanoma cell lines (Figure 4B). Although more data are needed to confirm actual functional differences, one possible consequence of the lack of Zn-Fingers could be a modified cellular localization. The PRDM9 isoforms without Zn-Fingers are indeed not able to reach cell nucleus . In addition, the alternative C-terminal region of PRDM7 shows no similarity to any sequence so far reported in sequence databases. It would be interesting to check whether this region is associated to the acquisition of any primate-specific function of the corresponding protein.
Detection of PRDM orthologs
To avoid the occurrence of aspecific hits due to the large number of Zn-Fingers, we only used the N-terminal portion of the PRDM proteins for our analysis. As a first scan of the presence of the PRDM genes in the 3 domains of life, we used each human PRDM protein as query for BLASTp  on the non-redundant database. To further confirm the presence of PRDMs only in metazoans, we searched the SMART database  for the entire repertoire of proteins with the same domain architecture. Finally, we searched the yeast genomes (sacSer1) using the 17 human PR domains as queries for tBLASTn . None of these analyses returned a significant hit for the PRDM genes in non-metazoan genomes.
We used different strategies to derive the orthologs of the human PRDM genes in other 13 metazoans (Pan troglodytes, Macaca mulatta, Mus musculus, Rattus norvegicus, Gallus gallus, Xenopus tropicalis, Takifugu rubripes, Ciona intestinalis, Strongylocentrotus purpuratus, Drosophila melanogaster, Drosophila pseudoscura, Caenorhabditis elegans, and Caenorhabditis briggsae), in order to conduct an orthology assignment as complete as possible. We started by extracting the reciprocal best hits for each PRDM proteins from the output of all all-against-all pairwise BLASTp  between human proteins and the protein set of each other species (RefSeq Release 17, April 2006; WormPep v.150; FlyBase 4.3). From the obtained results, we excluded the hits matching on the same chromosomal locus representing likely splice variants or mispredictions. To further enlarge the set of orthologs, we used the sequences of each human PR domain as queries for tBLASTn  against the different genome assemblies (human: hg18; chimp: panTro2; macaque: rheMac2; mouse: mm8; rat: rn4; chicken: galGal3; frog: xenTro2; fugu: fr1; ciona: ci2; sea urchin: Spur_0.5; fly: dm2; D. pseudoscura: dp4; worm: ce2; C. briggsae: cb1). The genomic region corresponding to each best hit was extracted and used for the complete gene prediction with GeneWise . Each predicted sequence was finally blasted back on the human protein set to confirm the orthology relationship through the best reciprocal method. In the cases where the ortholog was still missing, we defined the syntenic blocks around the PRDM genes in the UCSC human genome browser by identifying the genes at their 5' and 3'. We then extracted the corresponding regions in the target genome, and predicted the putative ortholog using GeneWise  and Blast2Gene . The complete set of the putative orthologs is reported in the Additional file 1.
Concatenated alignment and tree reconstruction
We based the evolutionary reconstruction of the PRDM family on the multiple alignment of the amino acidic sequence of the PR-domains. These are, together with the Zn-fingers, the only regions common to all members. Zn-fingers are however not suitable for the tree reconstruction given their repetitive nature. From the entire set of orthologs we only considered the sequences from seven vertebrates (Human, Chimp, Mouse, Rat, Chicken, Fugu, and Frog) and, among those, only the sequences for which a complete PR domain was retrievable. The final set of proteins used for the analysis was composed of 109 sequences, with representatives of each PRDM member in all 7 species except for PRDM3 (missing in fugu), PRDM17 (missing in fugu and chicken), and the non-primate PRDM7/9 (missing in chicken and frog).
For each PR domain of the different PRDM orthologous groups, we built a multiple sequence alignment (MSA) using MUSCLE , with the maximum number of iterations set to 100. In the resulting seventeen MSAs, the poorly aligned and divergent regions were eliminated with Gblocks . The reduced alignments were then treated as separate profiles and re-aligned using MUSCLE . The resulting global supermatrix is composed of 181 positions (see Additional file 3).
We derived Maximum Likelihood (ML) phylogenetic inferences using PHYML , and applying the JTT matrix. Our model of sequence evolution assumed that there were two classes of sites, one class being invariable and the other class being free to change. The rate variation across these sites was assumed to follow a gamma shape distribution calculated using a discrete approximation with eight categories of sites. Support for the hypotheses of relationships was assessed using 100 bootstrap replicates (see Additional file 4).
For the Bayesian inference, we used MrBayes  choosing the Jones substitution model and a gamma-shaped variation across the sites with four categories. The analysis was run for 500000 generations sampling the chain every 100 generations. To derive the final tree, we discarded the first 30% of the samples after checking for the convergence of the analysis after 150000 generations (see Additional file 5).
Screening of PRDM gene expression in human tissues
Total RNAs derived from normal human tissues (bladder, brain, breast, colon, fetal brain, fetal liver, kidney, liver, lung, ovary, pancreas, placenta, prostate, spleen, stomach, testis, thymus, thyroid and uterus) were obtained from Clontech (Mountain View, Cal., USA). In the case of hemopoietic precursor cells and human melanocytes, the RNA was extracted from CD34+ (from healthy donors) purified using standard procedures and from primary melanocyte cultures, respectively.
The expression levels of the PRDM genes were assessed using a custom TaqMan® Low Density Array (Applied Biosystems, Foster City, Cal., USA). For each PRDM gene a number of assays covering different regions of the transcript were used with 18S and GAPDH as reference genes (see Additional file 2). PRDM12, 13, 14, 7 and 9 are expressed only in very few tissues and, therefore, were not considered for further analysis. In order to reduce the technical variability due to the different efficiency of the assays, for all other genes the Ct value of each assay was plotted, and the corresponding trend in the different tissues was assessed (see Additional file 6). For each gene, the mean value only of the assays with a comparable trend was calculated and normalized using the reference genes. Finally the 2-ΔCt was calculated and normalized to the median of expression across each tissue (see Additional file 7). The resulting values were ordered according to the evolutionary divergence of the PRDM genes. We assessed the statistical significance of the negative association between gene expression level and speed of molecular evolution by plotting the logarithm of the maximum expression value for each PRDM gene against the corresponding distance to the root (Pearson correlation, r = -0.705, P = 0.01). To verify that the inverse relationship between expression data and evolutionary divergence was not due to an artifact, we treated the data from each single array in a comparable manner. In all cases comparable trends were obtained (data not shown).
RT-PCR and cloning of PRDM7-specific transcripts
Reverse transcription for subsequent PCR reactions was performed with Superscript III Reverse Transcriptase (Invitrogen) according to manufacturer's instructions. Amplification of PRDM7 long isoform was done using the high fidelity Phusion DNA polymerase (Fynnzyme) and different combinations of primers located in the 5'UTR and 3' UTR of the gene. Forward primers: PR7_9L-F1 (5'-TTCTAGACAGTCCCAGCACCATGA-3'), PR7L-F1 (5'-ACTCAGGGGCCCTTCCCACA-3'). Reverse primers: PR7L-R1 (5'-TGCAAGTGTGTGGTGATCACGT-3'), PR7L-R2 (5'-GTGTGTGGTGATCACGTTTGTCTTCT-3'), PR7L-R3 (5'-TTACTCATCCTTCCTGCAGAC-3'). A collection of normal tissues and commercial cancer cell lines was screened together with Reference RNA (Clontech). This is a mixture of total RNAs from a collection of adult human tissues, chosen to represent a broad range of expressed genes.
The PCR products were purified on Wizard SV GEL and PCR clean-up system (Promega) and, after overhangs addition with DyNAzyme II DNA polymerase (Fynnzyme), they were cloned into the pCR2.1-TOPO vector and transformed into TOP10 cells (TOPO TA cloning Kit, Invitrogen). Positive clones were screened by PCR using universal primers and sequenced. The region encompassing the duplicated sequence of PRDM7 was amplified using the following set of primers: forward 5'-TCGGCCAACTGGATGAGGTATGT-3' and reverse 5'-TCCCTGCCATGAGCTCTTTCTT-3'.
Rubin GM, Yandell MD, Wortman JR, Gabor Miklos GL, Nelson CR, Hariharan IK, Fortini ME, Li PW, Apweiler R, Fleischmann W: Comparative Genomics of the Eukaryotes. Science. 2000, 287 (5461): 2204-2215. 10.1126/science.287.5461.2204.
Chervitz SA, Aravind L, Sherlock G, Ball CA, Koonin EV, Dwight SS, Harris MA, Dolinski K, Mohr S, Smith T: Comparison of the Complete Protein Sets of Worm and Yeast: Orthology and Divergence. Science. 1998, 282 (5396): 2022-2028. 10.1126/science.282.5396.2022.
Lespinet O, Wolf YI, Koonin EV, Aravind L: The Role of Lineage-Specific Gene Family Expansion in the Evolution of Eukaryotes. Genome Res. 2002, 12 (7): 1048-1059. 10.1101/gr.174302.
Vogel C, Chothia C: Protein Family Expansions and Biological Complexity. PLoS Computational Biology. 2006, 2 (5): e48-10.1371/journal.pcbi.0020048.
Consortium IHGS: Initial sequencing and analysis of the human genome. Nature. 2001, 409 (6822): 860-921. 10.1038/35057062.
Ohno S: Evolution by gene duplication. 1970, Berlin-Heidelberg-New York: Springer-Verlang
Force A, Lynch M, Pickett FB, Amores A, Yan YL, Postlethwait J: Preservation of duplicate genes by complementary, degenerative mutations. Genetics. 1999, 151 (4): 1531-1545.
Stoltzfus A: On the possibility of constructive neutral evolution. J Mol Evol. 1999, 49 (2): 169-181. 10.1007/PL00006540.
Lynch M, Force A: The probability of duplicate gene preservation by subfunctionalization. Genetics. 2000, 154 (1): 459-473.
Turner CA, Mack DH, Davis MM: Blimp-1, a novel zinc finger-containing protein that can drive the maturation of B lymphocytes into immunoglobulin-secreting cells. Cell. 1994, 77 (2): 297-306. 10.1016/0092-8674(94)90321-2.
Yan J, Jiang J, Lim CA, Wu Q, Ng H-H, Chin K-C: BLIMP1 regulates cell growth through repression of p53 transcription. PNAS. 2007, 104 (6): 1841-1846. 10.1073/pnas.0605562104.
Deng Q, Huang S: PRDM5 is silenced in human cancers and has growth suppressive activities. Oncogene. 2004, 23 (28): 4903-4910. 10.1038/sj.onc.1207615.
He L, Yu JX, Liu L, Buyse IM, Wang MS, Yang QC, Nakagawara A, Brodeur GM, Shi YE, Huang S: RIZ1, but not the alternative RIZ2 product of the same gene, is underexpressed in breast cancer, and forced RIZ1 expression causes G2-M cell cycle arrest and/or apoptosis. Cancer Res. 1998, 58 (19): 4238-4244.
Huang S: Histone methyltransferases, diet nutrients and tumour suppressors. Nat Rev Cancer. 2002, 2 (6): 469-476. 10.1038/nrc819.
Jenuwein T: Re-SET-ting heterochromatin by histone methyltransferases. Trends Cell Biol. 2001, 11 (6): 266-273. 10.1016/S0962-8924(01)02001-3.
Hayashi K, Yoshida K, Matsui Y: A histone H3 methyltransferase controls epigenetic events required for meiotic prophase. Nature. 2005, 438 (7066): 374-378. 10.1038/nature04112.
Kim KC, Geng L, Huang S: Inactivation of a histone methyltransferase by mutations in human cancers. Cancer Res. 2003, 63 (22): 7619-7623.
Ancelin K, Lange UC, Hajkova P, Schneider R, Bannister AJ, Kouzarides T, Surani MA: Blimp1 associates with Prmt5 and directs histone arginine methylation in mouse germ cells. Nat Cell Biol. 2006, 8 (6): 623-630. 10.1038/ncb1413.
Davis CA, Haberland M, Arnold MA, Sutherland LB, McDonald OG, Richardson JA, Childs G, Harris S, Owens GK, Olson EN: PRISM/PRDM6, a transcriptional repressor that promotes the proliferative gene program in smooth muscle cells. Mol Cell Biol. 2006, 26 (7): 2626-2636. 10.1128/MCB.26.7.2626-2636.2006.
Moore AW, Jan LY, Jan YN: hamlet, a Binary Genetic Switch Between Single- and Multiple- Dendrite Neuron Morphology. Science. 2002, 297 (5585): 1355-1358. 10.1126/science.1072387.
Brumby A, Secombe J, Horsfield J, Coombe M, Amin N, Coates D, Saint R, Richardson H: A Genetic Screen for Dominant Modifiers of a cyclin E Hypomorphic Mutation Identifies Novel Regulators of S-Phase Entry in Drosophila. Genetics. 2004, 168 (1): 227-251. 10.1534/genetics.104.026617.
Steele-Perkins G, Fang W, Yang XH, Van Gele M, Carling T, Gu J, Buyse IM, Fletcher JA, Liu J, Bronson R: Tumor formation and inactivation of RIZ1, an Rb-binding member of a nuclear protein-methyltransferase superfamily. Genes Dev. 2001, 15 (17): 2250-2262. 10.1101/gad.870101.
Morishita K: Leukemogenesis of the EVI1/MEL1 gene family. Int J Hematol. 2007, 85 (4): 279-286. 10.1532/IJH97.06174.
Hayashi K, Matsui Y: Meisetz, a novel histone tri-methyltransferase, regulates meiosis-specific epigenesis. Cell Cycle. 2006, 5 (6): 615-620.
Birtle Z, Ponting CP: Meisetz and the birth of the KRAB motif. Bioinformatics. 2006, 22 (23): 2841-2845. 10.1093/bioinformatics/btl498.
Gyory I, Wu J, Fejer G, Seto E, Wright KL: PRDI-BF1 recruits the histone H3 methyltransferase G9a in transcriptional silencing. Nat Immunol. 2004, 5 (3): 299-308. 10.1038/ni1046.
Pal C, Papp B, Hurst LD: Highly expressed genes in yeast evolve slowly. Genetics. 2001, 158 (2): 927-931.
Drummond DA, Bloom JD, Adami C, Wilke CO, Arnold FH: Why highly expressed proteins evolve slowly. PNAS. 2005, 102 (40): 14338-14343. 10.1073/pnas.0504070102.
Ciccarelli FD, von Mering C, Suyama M, Harrington ED, Izaurralde E, Bork P: Complex genomic rearrangements lead to novel primate gene function. Genome Res. 2005, 15 (3): 343-351. 10.1101/gr.3266405.
Bailey JA, Eichler EE: Primate segmental duplications: crucibles of evolution, diversity and disease. Nat Rev Genet. 2006, 7 (7): 552-564. 10.1038/nrg1895.
Rogers JH: How were introns inserted into nuclear genes?. Trends Genet. 1989, 5 (7): 213-216. 10.1016/0168-9525(89)90084-X.
Venkatesh B, Ning Y, Brenner S: Late changes in spliceosomal introns define clades in vertebrate evolution. Proc Natl Acad Sci USA. 1999, 96 (18): 10267-10271. 10.1073/pnas.96.18.10267.
Knowles DG, McLysaght A: High rate of recent intron gain and loss in simultaneously duplicated Arabidopsis genes. Mol Biol Evol. 2006, 23 (8): 1548-1557. 10.1093/molbev/msl017.
Wagner E, Lykke-Andersen J: mRNA surveillance: the perfect persist. J Cell Sci. 2002, 115 (15): 3033-3038.
Raes J, Van de Peer Y: Functional divergence of proteins through frameshift mutations. Trends in Genetics. 2005, 21 (8): 428-431. 10.1016/j.tig.2005.05.013.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.
Letunic I, Copley RR, Pils B, Pinkert S, Schultz J, Bork P: SMART 5: domains in the context of genomes and networks. Nucl Acids Res. 2006, 34 (suppl_1): D257-260. 10.1093/nar/gkj079.
Birney E, Thompson JD, Gibson TJ: PairWise and SearchWise: finding the optimal alignment in a simultaneous comparison of a protein profile against all DNA translation frames. Nucleic Acids Res. 1996, 24 (14): 2730-2739. 10.1093/nar/24.14.2730.
Suyama M, Torrents D, Bork P: BLAST2GENE: a comprehensive conversion of BLAST output into independent genes and gene fragments. Bioinformatics. 2004, 20 (12): 1968-1970. 10.1093/bioinformatics/bth172.
Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucl Acids Res. 2004, 32 (5): 1792-1797. 10.1093/nar/gkh340.
Castresana J: Selection of Conserved Blocks from Multiple Alignments for Their Use in Phylogenetic Analysis. Mol Biol Evol. 2000, 17 (4): 540-552.
Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52 (5): 696-704. 10.1080/10635150390235520.
Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19 (12): 1572-1574. 10.1093/bioinformatics/btg180.
Letunic I, Bork P: Interactive Tree Of Life (iTOL): an online tool for phylogenetic tree display and annotation. Bioinformatics. 2007, 23 (1): 127-128. 10.1093/bioinformatics/btl529.
We thank Luisa Lanfrancone (IFOM-IEO Campus) for melanocyte primary cell cultures, Loris Bernard for useful discussion, Andrea Ciliberto and Anna De Grassi for critical comments on the manuscript. This work was funded by AIRC-NUSUG to F.D. Ciccarelli and by AIRC and Fondazione Cariplo to M. Alcalay.
The author(s) declares that there are no competing interests.
IF, NM, MA, FDC conceived and designed the experiments. IF and FDC performed sequence and the evolutionary analyses. NM and GS did the experimental work. MA, FDC, DR analyzed the data. FDC wrote the manuscript. All authors read and approved the final form of the manuscript.
Irene Fumasoni, Natalia Meani contributed equally to this work.
Electronic supplementary material
Additional file 1: PRDM orthologs in representatives of metazoans. Detailed description of the PRDM orthologs in 13 metazoans (Homo sapiens, Pan troglodytes, Macaca mulatta, Mus musculus, Rattus norvegicus, Gallus gallus, Xenopus tropicalis, Takifugu rubripes, Ciona intestinalis, Drosophila melanogaster, Drosophila pseudoscura, Caenorhabditis elegans, and Caenorhabditis briggsae). (PDF 93 KB)
Additional file 2: Expression data used for the analysis. The reported data represent the Ct values for each PRDM gene in 21 tissues obtained by qPCR. (XLS 50 KB)
Additional file 3: Multiple alignment used for the phylogenetic reconstruction. Multiple alignment of the vertebrate PRDM orthologs in phylip format. (PHY 26 KB)
Additional file 4: Tree data obtained with maximum likelihood. Phylogenetic tree of the PRDM genes in vertebrates using maximum likelihood in newick format. (TXT 3 KB)
Additional file 5: Tree data obtained with Bayesian inference. Phylogenetic tree of the PRDM genes in vertebrates using Bayesian inference in newick format. (TXT 3 KB)
Additional file 6: Trends of expression with different assays. The plots report the trend of expression of each human PRDM gene in 21 different tissues. (PDF 1 MB)
About this article
Cite this article
Fumasoni, I., Meani, N., Rambaldi, D. et al. Family expansion and gene rearrangements contributed to the functional specialization of PRDM genes in vertebrates. BMC Evol Biol 7, 187 (2007). https://doi.org/10.1186/1471-2148-7-187
- Intron Gain
- Orthology Assignment
- Krab Domain
- Ciona Genome