- Research article
- Open Access
Retrotransposition and mutation events yield Rap1 GTPases with differential signalling capacity
- Tomasz Zemojtel†1, 2Email author,
- Marlena Duchniewicz†4,
- Zhongchun Zhang†5,
- Taisa Paluch5,
- Hannes Luz1,
- Tobias Penzkofer2, 3,
- Jürgen S Scheele4 and
- Fried JT Zwartkruis5Email author
© Zemojtel et al; licensee BioMed Central Ltd. 2010
Received: 6 January 2010
Accepted: 19 February 2010
Published: 19 February 2010
Retrotransposition of mRNA transcripts gives occasionally rise to functional retrogenes. Through acquiring tempero-spatial expression patterns distinct from their parental genes and/or functional mutations in their coding sequences, such retrogenes may in principle reshape signalling networks.
Here we present evidence for such a scenario, involving retrogenes of Rap1 belonging to the Ras family of small GTPases. We identified two murine and one human-specific retrogene of Rap1A and Rap1B, which encode proteins that differ by only a few amino acids from their parental Rap1 proteins. Markedly, human hRap1B-retro and mouse mRap1A-retro1 acquired mutations in the 12th and 59th amino acids, respectively, corresponding to residues mutated in constitutively active oncogenic Ras proteins. Statistical and structural analyses support a functional evolution scenario, where Rap1 isoforms of retrogenic origin are functionally distinct from their parental proteins. Indeed, all retrogene-encoded GTPases have an increased GTP/GDP binding ratio in vivo, indicating that their conformations resemble that of active GTP-bound Rap1. We furthermore demonstrate that these three Rap1 isoforms exhibit distinct affinities for the Ras-binding domain of RalGDS. Finally, when tested for their capacity to induce key cellular processes like integrin-mediated cell adhesion or cell spreading, marked differences are seen.
Together, these data lend strong support for an evolution scenario, where retrotransposition and subsequent mutation events generated species-specific Rap1 isoforms with differential signaling potential. Expression of the constitutively active human Rap1B-retro in cells like those derived from Ramos Burkitt's lymphoma and bone marrow from a patient with myelodysplastic syndrome (MDS) warrants further investigation into its role in disease development.
Development of new functional genes is central to the evolution of species-specific traits . Two major mechanisms allow for creation of new genes based on existing templates: duplication of chromosomal regions containing gene loci and retrotransposition of mRNA transcripts into the genome - a process which in mammals is facilitated by active LINE-1 retrotransposons (L1s) [2, 3]. Novel functions for processed genes, termed retrogenes, may come from the acquirement of new spatio-temporal expression patterns, dictated by the genomic context of cDNA insertion [4–8]. Furthermore, functional evolution of their regulatory sequences and coding sequences, driven by natural selection, could result in new functions within signaling networks [5, 9, 10].
Evidence for the latter mode of gene evolution is scarce. First of all, expression of retrogenes has only recently begun to be explored. Although expression of processed genes seems to be more common than previously anticipated, surveys are frequently limited to a single species [9–11]. In addition, the physiological consequences of altered gene expression patterns [4, 5, 12] and/or the biochemical consequences of amino acid substitutions are rarely investigated [1, 5].
In the present study, we provide evidence for the functional evolution of retrogenes derived from the Ras-like GTPase, Rap1. The activity of Rap1 is dynamically controlled by many different receptor proteins, including tyrosine receptor kinases and G-protein coupled receptors . These receptors regulate the activity of guanine nucleotide exchange factors (GEFs), resulting in GTP-binding of Rap1, or GTPase activating proteins (GAPs) that strongly increase the hydrolyzing capacity of Rap1 and thereby inactivate Rap1 . Rap1 is involved in many different morphogenetic processes e.g. by affecting the activity of cell-matrix or cell-cell adhesion receptors. In addition, Rap1 plays important physiological roles in established structures like for example the vertebrate blood vessel system. Here, an acute rise in cAMP results in immediate enhancement of the level of GTP-bound Rap1 that is crucial for an increase in vascular permeability .
Rap1A and Rap1B are paralogs and differ only by 9 amino acids, that are mostly located at the very C-terminal part . Rap1 proteins are highly conserved in evolution and in fact identical Rap1A as well as Rap1B proteins are found in human, chimp, macaque, cow and mouse. While both proteins are ubiquitously expressed in tissues, a differential distribution of isoforms has been observed, e.g. Rap1B protein is predominant in platelets . Possibly, the high level of sequence conservation and differential gene expression reflects specific functions for each of these GTPases. This would be consistent with the outcome of studies in mice with targeted disruption of Rap1 genes. Loss of Rap1B impairs platelet aggregation due to a diminished activation of αIIbβ3 integrins. In addition, defects in endothelial cells may cause bleedings, resulting in a high percentage of embryonic lethality . Rap1A knockout mice are fully viable and fertile, but haematopoietic cells from the spleen and thymus adhere less efficiently to fibronectin or collagen as compared to wild type cells . In agreement with earlier in vitro observations , a recent work reported a reduced superoxide production in Rap1A -/- neutrophils . Double Rap1A/Rap1B knockout mice die during early embryogenesis suggesting functional overlap between the two Rap1 isoforms .
Our current study suggests that additional complexity resides within the Rap1 signaling network. We characterize here two murine-specific Rap1A retrogenes and one human-specific Rap1B retrogene that acquired limited substitutions in amino acids, including the residues G12 and A59 that were found to be mutated in oncogenic Ras isoforms. By using computational and experimental approaches we collected evidence that the species-specific events of retrotransposition and mutation created new Rap1 isoforms whose signaling capacity differs from that of their parental isoforms. We discuss potential roles for Rap1 retrogenes in cell physiology.
Identification of Rap1A and Rap1B retrogenes with limited mutations
Importantly, we previously pointed out that also in human there is evidence for expression of a Rap1B retrogene (G12R, K42E), named hRap1B-retro . More specifically, database searches revealed expression of this retrogene in Ramos Burkitt's lymphoma, neuroblastoma cell lines and placental tissue, while a study by Gyan et al. detected expression in a bone marrow sample from a patient with myelodysplastic syndrome (MDS) . As indicated above for Rap1A retrogenes, also the hRap1B-retro may encode for a Rap1 protein with elevated GTP content because of a mutation at position twelve (G12R), but no functional characterization of the encoded protein was so far undertaken.
Genomic localization and origin of Rap1A and Rap1B retrogenes
Moreover, phylogenetic analysis revealed that the group of the two Rap1A retrogenes and the mouse Rap1A gene form a monophyletic group in the tree with other mammalian Rap1A genes (with bootstrap support of 94%, see [Additional file 3: Supplemental Figure S3]). Thus, Rap1A retrogenes likely split off the mouse lineage. Together, this strongly implies that an event of insertion of a Rap1A retrogene into the 14th intron of Gtf2 h1 was followed by a duplication event and that the both events date back to times after rodent speciation. After the duplication, both Rap1A retrogenes diverged by non-synonymous mutations (i.e. mutations leading to amino acid changes) in their ORF regions (see the next subchapter).
We resolved the location of the hRap1B-retro retrogene to the region within the antisense strand of the 2nd intron of SV2C (synaptic vesicle protein 2 isoform C) on chromosome 5 (Figure 3B). This Rap1B retrogene was only identified in the human genome sequence (i.e. was absent from the chimp genome), as revealed by genomic blast searches (see Methods).
Interestingly, hRap1B-retro's mRNA sequence maps to 2 exons that are separated by an intron, which was created by L1-mediated insertion of the SVA element into a preferred target site sequence (TTTT/AA)  present in human Rap1B mRNA (at 1304 nt, Accession: 58219793). We could reconstruct the insertion site for the human Rap1B retrogene based on the comparison with the chimp genome (Figure 3). The 2nd intron of SV2C gene in the chimp genome bears the insertion site 5'-TTTT/A-3', typical of the site recognized by L1's endonuclease. Presence of target site duplications (TSDs) flanking Rap1A and Rap1B retrogenes provides further evidence for involvement of the LINE-1 retrotransposition machinery (Figure 3). Phylogenetic analysis established that hRap1-retro and the Rap1B in human and chimp form a monophyletic group (bootstrap support of 95%, see [Additional file 3: Supplemental Figure S3]). This implies that hRap1B-retro arose after divergence of macaca from the ancestor of chimp and human.
Functional evolution scenario
We tested whether it is appropriate to assume two different ω-ratios in the tree within the clades model by computing the likelihood under the simpler M0-model of codeml and performing a Likelihood Ratio Test (LRT). The simpler M0-model, which assumes only one ω-ratio for the whole tree, is a special case of the clades model. Because the M0 model is nested in the clades model, the likelihood under the clades model is by definition larger or equal than the likelihood under the M0-model for every alignment. The test statistic Δ of the LRT is two times the difference in log likelihoods, Δ = 2 ln L(clades-model) - ln L(M0). The null assumption H0 of the LRT corresponds to the assumption that sequences in the multiple alignment had evolved according to the simpler M0-model. Then, under the null assumption H0, the expected distribution of Δ is a χ2 distribution with one degree of freedom. The LRT reveals for both alignments that the clades model with two ω-ratios performs significantly better than the M0-model of codeml, computed p-values were prap1A < 10-6 and prap1B = 0.0013.
The fact that ω-ratios are larger for retrogenes suggests that the selective pressure on the replacement of amino acids was relaxed during the evolution of the retrogenes. However, ω2 estimates are at least two standard deviations smaller than one for both, for the alignment of Rap1A genes as well as for the alignment of Rap1B genes. Thus, synonymous substitutions, which are assumed to be neutral in Yang's codon substitution model, occurred significantly more often than non-synonymous substitutions and we argue that the evolution of retrogenes was subject to moderate selective constraints.
Interestingly, Lynch et al.  compared measures of ω in paralogous genes of several eukaryotic genomes. The authors demonstrated that ω systematically decreases in time and concluded that "early in their evolutionary history, duplicated genes tend to be under moderate selective constraints" where the range for ω lies between values of 0.37 and 0.46. We thus argue, that the values of ω2(Rap1A) = 0.26 and ω2(Rap1B) = 0.33 indicate that the evolution of mRap1A-retro1, -retro2 genes and hRap1B-retro is constrained by a regime of at least "moderate selection".
The tree dendrogram of Figure 4B shows the evolution of the two mouse Rap1A retrogenes after a retrotransposition event and a subsequent duplication event. If we assume the expected number of neutral synonymous substitutions under the 2-clades codon evolution model (described above) to be proportional to the elapse of time, then it can be speculated that the duplication event occurred only shortly after the retrotransposition event. Moreover, we found evidence for a scenario according to which the mutations at codon 35 resulting in synonymous and non-synonymous (T35 M) substitutions in mRap1A-retro1 and mRap1A-retro2, respectively, were generated by deamination of methylated CpG (Figure 4b).
We further employed structural analysis and functional assays to uncover evidence for the molecular basis of evolutionary pressure to conserve these genes.
Structure-related functional analysis
In order to functionally characterize the Rap1 retrogene encoded proteins, we first tested for their ability to bind guanine nucleotides. mRap1A-retro1 and hRap1B-retro contain mutations in 59th (in G3 motif) and 12th (in G1 motif) amino acid, respectively, that are involved in the control of interactions with the GTP ligand (Figure 1B). Since mutations at position 12 and 59 are transforming mutations in Ras [29–31], it is expected that human Rap1B-retro and mouse mRap1A-retro1 are more GTP-bound than wild type Rap1.
Previously, it has been reported that V12 mutations in Rap1 are still sensitive to Rap1GAP .
GAPs (GTPase activating proteins) are proteins that mediate negative regulation of Ras-like GTPases by catalyzing the hydrolysis of GTP. Co-transfection of Rap1GAP with any of the retrogenes caused a strong decrease in the GTP/GDP ratio. However, this decrease was clearly less pronounced for mRap1A-retro2 (Figure 5C) underlining the importance of T35 for GTP hydrolysis but also its involvement in control of interactions with target proteins. It was shown that mutations in T35 residue perturb a dynamic behavior of the switch I region  and this can affect interactions with effector proteins but also regulatory factors like GAP which bind to this region .
One of the best-documented functions of Rap1 is to enhance integrin-dependent adhesion to extra-cellular matrices . We assayed the ability of retrogene isoforms to increase adhesion of HB6 cells (Jurkat T cells expressing Epac1).
Together these results demonstrate that the proteins encoded by the three Rap1 retrogenes differ in their biological activities.
Reverse transcription of mRNAs followed by insertion into the genome is an event, which in germ cells can yield additional copies of genes for future generations. In most cases they will develop into transcriptionally silent pseudogenes, which harbor mutations in their coding regions, preventing the production of functional proteins. However, if retrotransposition yields processed genes that have a distinct tempero-spatial expression pattern and/or have acquired neomorphic functions as a consequence of mutation, their presence may be advantageous to the organism and expression may be maintained .
Indeed, a number of reports have recently appeared, suggesting that expression of processed genes (termed retrogenes) might be more common than previously anticipated. For example, the DPPA3 gene, which has been used as a marker for pluripotency, was found to have a number of related retrogenes, of which at least one was found to be expressed in adult tissues . Expressed retrogenes are not simple evolutionary relics, but can obtain essential functions as was shown for the Utp14B protein. Mutation of this autosomal gene results in juvenile spermatogonial depletion .
In the course of analysis of our Rap1A knockout mice, we isolated cDNAs from various tissues, corresponding to two distinct Rap1A retrogenes. The existence of Rap1-like retrogenes appears not to be unique to mice. Recently, we and others reported the existence of hRap1B-retro retrogene in human, whose expression might be linked to development of myelodysplastic syndromes and other types of leukemia [23, 24]. Genomic comparisons indicated that the Rap1A and Rap1B retrogenes are species-specific and that the Rap1A retrogenes are paralogs. Interestingly, the recently emerged Rap1 retrogenes acquired limited amino acid-changing substitutions and our statistical analysis favors a scenario whereby they evolved under selective constraints on coding sequence. In this respect, hRap1B-retro (G12R, K42E) and mRap1A-retro1 (A59V) were of special interest since they had substitutions at positions 12 and 59 (Figure 1), which resembled those found in constitutively active oncogenic Ras proteins. Both retrogene-encoded isoforms had an increased GTP/GDP binding ratio and interacted strongly with the Ras-binding domain of RalGDS in pull-down assays. Strikingly, only hRap1B-retro increased integrin-mediated cell adhesion, whereas mRap1A-retro1 was practically inactive in this assay. Cell spreading on the other hand was increased upon over-expression of both hRap1B-retro and mRap1-retro1. These results are reminiscent of those obtained with various Rap1 mutants, generated via site directed mutagenesis . In that study, Rap1A G12V appeared to be more powerful in suppressing growth of tumorigenic HT1080 cells than wild type Rap1, but its capacity to induce a flat revertant phenotype of in Ras-transformed DT cells was not strongly increased. In contrast, a Rap1A A59T mutant showed opposite characteristics with a mild effect on HT1080 proliferation but a clearly enhanced effect in the flat revertant assay. Together, these results indicate that the proteins encoded by the Rap1 retrogenes have different biological activities, and uncouple cell adhesion and cell spreading, possibly by activating distinct (sets of) Rap1 effectors. Indeed, Rap1 effectors have been described, which preferentially increase cell spreading like VAV  or both cell spreading and cell-matrix interaction like RIAM . It will be interesting to see if the Rap1 retrogene products indeed selectively interact with such effectors. The fact that Rap1 retrogene encoded proteins are in a GTP-bound state suggest that retrogene expression would allow for activation of Rap1-dependent processes in a stimulus-independent way. This phenomenon is well known in the case of Ras, where somatic mutations result in oncogenic proteins, which stimulate the ERK and PI-3K pathways in the absence of growth factors. Likewise, expression of Rap1 retrogenes may activate the above-mentioned effectors VAV or RIAM or any of the other Rap1-effectors that include diverse proteins like phospholipase Cε or the scaffold protein AF6 (reviewed in ). The precise mode by which retrogene encoded Rap1 isoforms interact with known Rap1 effectors or perhaps additional proteins will determine the physiological consequences. These may include effects on cell adhesion, important for neural pathfinding and immune surveillance or cell-cell adhesion, contributing e.g. to vascular permeability.
Collectively, the data above strongly indicate that mouse mRap1A-retro1 and human hRap1B-retro retrogenes represent recent examples of functional evolution. However, it has to be taken into consideration that not all of the described retrogenes will be selectively maintained and neo-functionalized. Indeed, we have not been able to detect a biological effect of mRap1A-retro2 over-expression. This may lend support for the scenario where mRap1A-retro2 is unlikely to escape non-functionalization.
Although it is hard to rule out that the other two retrogenes are on their way to non-functionalization, the fact that mRap1A-retro1 solely stimulates cell spreading but not adhesion, points to a unique property of this isoform. In addition, the fact that it like hRap1B-retro encodes for constitutively active protein, indicates that it can function largely independently of GEFs. Although one could argue that our current study does not reveal truly novel functions to warrant the term neo-functionalization, future studies that should try to identify downstream effectors for e.g. mRap1A-retro2 should shed further light on this issue.
It is of note that the Ras-like protein family contains at least two other members of retrogenic origin, a H-Ras retrogene named Eras, and a representative of Rap subfamily, Rap2B. Interestingly, mouse and human ERas, acquired mutations at positions 12, 59, 62 that render constitutively active proteins . ERas is absent from genomes of fish and amphibians and appears to be exclusively present in mammalian genomes. Its primary function was linked to its highly specific domain of expression in mouse ES cells , yet expression of ERas was not detected in human ES cells . On the other hand, the Rap2B retrogene, originally discovered in human platelets , is vastly conserved in vertebrate genomes. By contrast, Rap1 retrogenes are products of recent evolutionary events that created species-specific isoforms with differential signaling capacity. Although their physiological significance remains to be established, the fact that all of them are GTP-loaded, yet have distinct signaling capacities, shows how limited amino acid substitutions may turn recently emerged retrogenes into functionally novel genetic entities.
Sequences of Rap1 orthologs
The sequences of Rap1A and Rap1B orthologs were extracted using human and mouse Rap1A/B sequences (NM_145541, NM_015646) as queries in tblastx searches against the chimp, macaca, rat, dog and cow genomes.
In order to identify orthologous loci of mouse Rap1A and human Rap1B retrogenes in the rat and in the chimp genomes, respectively, we utilized nucleotide and protein blast searches. We checked for presence/absence of Rap1B and Rap1A retrogenes in the sequences of human, chimp, orangutan, rhesus, marmoset, mouse, rat, guinea pig, cat, dog, horse, cow, opossum and platypus genomes (as made available at the UCSC genome browser). The multiple alignments reported in the [Additional file 2: Supplemental Figure S2] and in the [Additional file 3: Supplemental Figure S3] were computed with ClustalX and visually inspected afterwards.
Given the alignment of all Rap1 genes, we wanted to infer a gene tree according to which the Rap1 genes putatively evolved. We used jModelTest (version 0.1.1) to identify the nucleotide substitution model that best fits the alignment  and PhyML  to infer possible gene trees. The "HKY85 + I" nucleotide substitution model (which parameterizes substitution rates with a parameter for the transition-transversion ratio, accounts for different nucleotide frequencies, and assumes that a fraction of sites is invariable) was always selected as the best-fit model according to the Akaike Information Criterion (with and without second order correction) and the Bayes Information Criterion.
We then ran PhyML with "HKY85 + I" as nucleotide substitution model on the alignment as well on 100 bootstrap replicates of the alignment. The resulting tree with annotated bootstrap values is shown in [Additional file 3: Supplemental Figure S3].
The ML tree does not reflect the accepted species phylogeny. However, the Shimodaira-Hasegawa (SH)  test to compare different assumptions about the tree revealed that the maximal likelihood of the tree computed with PhyML is not significantly different from the likelihood under the tree that reflects the species phylogeny (shown in Figure 4).
Compounding both, the SH test and weak bootstrap support values at basal branches, we concluded that the alignment is too conserved to carry enough signal to resolve basal branching patterns. We therefore assumed that orthologs in the alignment have evolved along the branches of the established species phylogeny (Figure 4) and confirmed that the computation of ω values is not sensitive to different basal branching patterns. For example we considered the topology ((rodents, (cow, dog)), primates) and compared results to the ones under ((rodents, primates), (cow, dog)) which represents the established species phylogeny. The ω values shown in Table in Figure 4A are the same under the two alternative tree topologies.
RNA isolation and reverse transcription
RNA was isolated using standard protocols and deprived of DNA using the DNAseI kit from Promega Corporation according to the manufacturer's recommendations.
cDNA was made using the SuperscriptTM II Reverse Transcripatase kit from Invitrogen according to the manufacturer's protocol. For each RNA sample a control reaction was prepared w/o RT polymerase to check whether Dnase I digestion had been successful.
RT- PCR reaction and cloning of retrogenes
For PCR reactions, the Taq PCR core kit fom Qiagene was used. 2-4 μl of RT reaction was mixed with PCR buffer, dNTPs, water, Taq polymerase and 10 pmol of 5' and 3' specific primers. The following PCR conditions were used: 96°C, 2 min. 25-35 cycles: 96°C-20 s, Tm-50°C -20 s, 72°C-1 min.: 72°C, 10 min. Using the following primer set we amplified a 1186 bp fragment of Rap1A as well as mRap1A-retro1 and mRap1A-retro2 genes: 5'-GCGGGATTGTCAATATTTAAAC-3' and 5'-GCCATAGAAATCAGTTATCCC-3'. This primer set amplifies the whole ORF and a part of the 3'UTR. To specifically amplify Rap1A retrogenes from Rap1A -/- tissues one primer was set in the targeted second exon: 5'-TCAGGAGGCGTGGGGAAG-3' and the other in the 3' UTR: 5'-GCCATAGAAATCAGTTATCCC-3'. PCR products from multiple tissues (i.e. thymus, liver, spleen, testis) were cloned into vector from TOPO TA cloning kit (Invitrogen) and then sequenced using standard primers recommended by manufacturer.
For the expression and activity studies both retrogenes were amplified using Pfu polymerase and primers compatible for cloning full length proteins into pMT2HA.
Gene Racer reaction - full length RNA ligase-mediated rapid amplification of 5' cDNA ends (RLM-RACE PCR)
The Gene Racer kit (Invitrogen) was used according to the manufacturer's recommendation. Briefly, RNA was dephosphorylated to eliminate truncated mRNA and non-mRNA from subsequent ligation with the GeneRacer RNA Oligo. Next, the cap structure from intact, full-length mRNA was removed using Tap enzyme. This treatment generates a 5'phosphate required for ligation to the "GeneRacer RNA Oligo". Subsequently, ligation of RNA Oligo to decapped mRNA was performed utilizing T4 RNA ligase. Then, cDNA was produced on this RNA using random primers. Hereafter, PCR was performed in order to obtain 5' ends. For this purpose, the first-strand cDNA was amplified using "GeneRacer 5' Primer" (homologous to the GeneRacer RNA Oligo): 5'-CGACTGGAGCACGAGGACACTGA-3' and a reverse gene-specific primer located in the targeted second exon: 5'-CTTCCCCACGCCTCCTGAACCAA-3'. Only mRNA that ligated to ''GeneRacer RNA Oligo" and is completely reverse transcribed will be amplified. This PCR was prepared the standard way using Taq polymerase from Quiagen. Subsequently, this PCR reaction was further used for nested PCR. Here, "5'GeneRacer Nest primer": 5'-GGACACTGACATGGACTGAAGGAGTA-3' and reverse nested gene specific primer located in the targeted second exon: 5'-CTTCCCCACGCCTCCTGAACCAA-3' were used. The products obtained as single dominant bands were then cloned into pCR4-TOPO vector (TOPO TA cloning kit Invitrogen) and sequenced. The sequence derived from wild type material, was a 5' sequence of mRap1A mRNA. The PCR product from knockout material was corresponding to a 5' UTR of mRap1A-retro1. The 5'UTR region of mRap1A-retro1 is longer (by ~550 nt) than the 5'UTR of mRap1A and different from the 5' UTR of mRap1A (as also evident from the genomic sequence of mRap1A-retro1).
Sequencing of Cosmid DNA harboring mRap1A-retro2
Cosmid MPMGc121C10270Q2, obtained from RZPD (German Science Centre for Genome Research) was identified during the screen for the Rap1A gene. We sequenced 1741 nt of this cosmid harboring mRap1A-retro2, using the following oligos:
5'-GCGGGATTGTCAATATTTAAAC-3'; 5'-CTTCCCCACGCCTCCTGAACCAA-3'; 5'-CACCAGTGGAAAAGAAGAAGCC-3' and 5'-GGCACAGTTACACCACTGTCTTG-3'.
Cos-7 cells were grown in Dulbecco's modified Eagles medium (DMEM) containing 10% fetal calf serum (FCS) and serum-starved in DMEM at 1.5% FCS. A431 cells expressing human p120 short interfering RNA (siRNA) was described previously . Cells were cultured in DMEM containing L-glutamine (Hyclone), 10% fetal bovine serum (FBS) (Hyclone), and 1% penicillin/streptomycin (Gibco/Invitrogen).
Cloning of hRap1B-retro into pMT2HA vector
The fragment containing the ORF of hRap1Bretro was amplified from human genomic DNA using primers: 5'-ATGTCGACCATGCGTGAGTATAAGCTAGTCGT-3' and 5'-ATGCGGCCGCAGACCTGGCTCAGAGCTACA-3'. We used Advantage polymerase from Clontech with PCR conditions: 1 cycle: 95°C 1 min./35 cycles: 95°C 15 sec, 68°C 1 min/1 cycle: 68°C 10 min. The PCR product was cloned into vector from TOPO TA cloning kit (Invitrogen) and then re-cloned into SalI and NotI of pMT2HA vector.
Activity measurements for Rap1 and related proteins
Determination of the GDP/GTP ratio using [32P]-labeling of cells or using a non-radio active pull-down assay were done as described in Zwartkruis et al. 1998 . Haemaglutamin (HA)-tagged Rap1, Rap1V12 (HA-RapV12), Rap1GAP (HA-RapGAP I) have previously been described .
Cell spreading assay
To measure cell spreading, A431 p120 RNAi cells  were plated on glass cover slips. The next day cells were transfected with either empty vector or pMT2HA vectors expressing the various Rap1 cDNAs and histon2B::GFP using Fugene. After 24 hours, cells were fixed and stained either for the HA-tag to confirm that GFP positive cells also expressed Rap1 protein or with TRITC phalloidin. Cell surface area was measured using the ImageJ program version 1.36b (NIH).
For adhesion assays, transiently transfected Jurkat cells, serum starved overnight, were harvested, washed, and suspended in TSM buffer (20 mM Tris-HCl, pH 8.0, 150 mM NaCl, 1 mM CaCl2, 2 mM MgCl2) at a concentration of 5 × 105 cells/ml. 24-Well Nunc Maxisorp plates (Corning, NY, USA) were coated with fibronectin (5 μg/ml) overnight at 4°C, washed, and blocked for 1 hour at 37°C with 1% bovine serum albumin (BSA), TSM. After washing, 200 μl of TSM was added per well with or without the indicated stimuli. Subsequently, 200 μl cell suspension was added per well. Cells were allowed to adhere for 1 hour at 37°C and non-adherent cells were removed with warmed 0.5% BSA, TSM. Adherent cells were lysed and subjected to a luciferase assay as described previously . Expression of transfected constructs was confirmed by immunoblotting of total cell lysates. Adherent cells were calculated, and the cell numbers were corrected for transfection efficiency and non-specific effects of constructs by measuring luciferase activity of total input cells ((counts in cells bound/counts in total input cells) × 100%).
We thank Johannes L. Bos and Martin Vingron for valuable discussions and a critical reading of the manuscript. We thank Dr. A. Reynolds for the A431 p120 RNAi cells. We thank Bas van Breukelen for help with statistical analysis of the cell spreading experiments. This work was supported by the Netherlands Genomics Initiative through the Cancer Genomics Center (to Z.Z.) and the Center for Biomedical Genetics (TP and FZ). JS received founding from Deutsche Krebshilfe (10-728 sche2). TZ received founding from the European Commission within its FP6 Programme "Biosapiens", under the thematic area 'Life sciences, genomics and biotechnology for health', contract number LHSG-CT-2003-503265.
- Zhang J, Zhang YP, Rosenberg HF: Adaptive evolution of a duplicated pancreatic ribonuclease gene in a leaf-eating monkey. Nat Genet. 2002, 30: 411-5. 10.1038/ng852.View ArticlePubMedGoogle Scholar
- Esnault C, Maestre J, Heidmann T: Human LINE retrotransposons generate processed pseudogenes. Nat Genet. 2000, 24: 363-7. 10.1038/74184.View ArticlePubMedGoogle Scholar
- Penzkofer T, Dandekar T, Zemojtel T: L1Base: from functional annotation to prediction of active LINE-1 elements. Nucleic Acids Res. 2005, 33: D498-500. 10.1093/nar/gki044.PubMed CentralView ArticlePubMedGoogle Scholar
- Krasnov AN, Kurshakova MM, Ramensky VE, Mardanov PV, Nabirochkina EN, Georgieva SG: A retrocopy of a gene can functionally displace the source gene in evolution. Nucleic Acids Res. 2005, 33: 6654-61. 10.1093/nar/gki969.PubMed CentralView ArticlePubMedGoogle Scholar
- Burki F, Kaessmann H: Birth and adaptive evolution of a hominoid gene that supports high neurotransmitter flux. Nat Genet. 2004, 36: 1061-3. 10.1038/ng1431.View ArticlePubMedGoogle Scholar
- Long M, Betran E, Thornton K, Wang W: The origin of new genes: glimpses from the young and old. Nat Rev Genet. 2003, 4: 865-75. 10.1038/nrg1204.View ArticlePubMedGoogle Scholar
- Brosius J: Retroposons--seeds of evolution. Science. 1991, 251: 753-10.1126/science.1990437.View ArticlePubMedGoogle Scholar
- Brosius J, Gould SJ: On "genomenclature": a comprehensive (and respectful) taxonomy for pseudogenes and other "junk DNA". Proc Natl Acad Sci USA. 1992, 89: 10706-10. 10.1073/pnas.89.22.10706.PubMed CentralView ArticlePubMedGoogle Scholar
- Marques AC, Dupanloup I, Vinckenbosch N, Reymond A, Kaessmann H: Emergence of young human genes after a burst of retroposition in primates. PLoS Biol. 2005, 3: e357-10.1371/journal.pbio.0030357.PubMed CentralView ArticlePubMedGoogle Scholar
- Vinckenbosch N, Dupanloup I, Kaessmann H: Evolutionary fate of retroposed gene copies in the human genome. Proc Natl Acad Sci USA. 2006, 103: 3220-5. 10.1073/pnas.0511307103.PubMed CentralView ArticlePubMedGoogle Scholar
- Emerson JJ, Kaessmann H, Betran E, Long M: Extensive gene traffic on the mammalian X chromosome. Science. 2004, 303: 537-40. 10.1126/science.1090042.View ArticlePubMedGoogle Scholar
- McCarrey JR, Thomas K: Human testis-specific PGK gene lacks introns and possesses characteristics of a processed gene. Nature. 1987, 326: 501-5. 10.1038/326501a0.View ArticlePubMedGoogle Scholar
- Zwartkruis FJ, Wolthuis RM, Nabben NM, Franke B, Bos JL: Extracellular signal-regulated activation of Rap1 fails to interfere in Ras effector signalling. Embo J. 1998, 17: 5905-12. 10.1093/emboj/17.20.5905.PubMed CentralView ArticlePubMedGoogle Scholar
- Bos JL, Rehmann H, Wittinghofer A: GEFs and GAPs: critical elements in the control of small G proteins. Cell. 2007, 129: 865-77. 10.1016/j.cell.2007.05.018.View ArticlePubMedGoogle Scholar
- Fukuhra S, Sakurai A, Yamagishi A, Sako K, Mochizuki N: Vascular endothelial cadherin-mediated cell-cell adhesion regulated by a small GTPase, Rap1. J Biochem Mol Biol. 2006, 39: 132-9.View ArticlePubMedGoogle Scholar
- Bokoch GM: Biology of the Rap proteins, members of the ras superfamily of GTP-binding proteins. Biochem J. 1993, 289 (Pt 1): 17-24.PubMed CentralView ArticlePubMedGoogle Scholar
- Klinz FJ, Seifert R, Schwaner I, Gausepohl H, Frank R, Schultz G: Generation of specific antibodies against the rap1A, rap1B and rap2 small GTP-binding proteins. Analysis of rap and ras proteins in membranes from mammalian cells. Eur J Biochem. 1992, 207: 207-13. 10.1111/j.1432-1033.1992.tb17039.x.View ArticlePubMedGoogle Scholar
- Chrzanowska-Wodnicka M, Smyth SS, Schoenwaelder SM, Fischer TH, White GC: Rap1b is required for normal platelet function and hemostasis in mice. J Clin Invest. 2005, 115: 680-7.PubMed CentralView ArticlePubMedGoogle Scholar
- Duchniewicz M, Zemojtel T, Kolanczyk M, Grossmann S, Scheele JS, Zwartkruis FJ: Rap1A-deficient T and B cells show impaired integrin-mediated cell adhesion. Mol Cell Biol. 2006, 26: 643-53. 10.1128/MCB.26.2.643-653.2006.PubMed CentralView ArticlePubMedGoogle Scholar
- Bokoch GM, Quilliam LA, Bohl BP, Jesaitis AJ, Quinn MT: Inhibition of Rap1A binding to cytochrome b558 of NADPH oxidase by phosphorylation of Rap1A. Science. 1991, 254: 1794-6. 10.1126/science.1763330.View ArticlePubMedGoogle Scholar
- Li Y, Yan J, De P, Chang HC, Yamauchi A, Christopherson KW, Paranavitana NC, Peng X, Kim C, Munugulavadla V, Kapur R, Chen H, Shou W, Stone JC, Kaplan MH, Dinauer MC, Durden DL, Quilliam LA: Rap1a null mice have altered myeloid cell functions suggesting distinct roles for the closely related rap1a and 1b proteins. J Immunol. 2007, 179: 8322-31.PubMed CentralView ArticlePubMedGoogle Scholar
- Chrzanowska-Wodnicka M, Kraus AE, Gale D, White GC, Vansluys J: Defective angiogenesis, endothelial migration, proliferation, and MAPK signaling in Rap1b-deficient mice. Blood. 2008, 111: 2647-56. 10.1182/blood-2007-08-109710.PubMed CentralView ArticlePubMedGoogle Scholar
- Zemojtel T, Penzkofer T, Duchniewicz M, Zwartkruis FJ: HRap1B-retro: a novel human processed rap1B gene blurs the picture?. Leukemia. 2006, 20: 145-6. 10.1038/sj.leu.2404040.View ArticlePubMedGoogle Scholar
- Gyan E, Porteu F, Fontenay M: Reply to 'hRap1B-retro: a novel human processed Rap1B gene blurs the picture?' by Zemojtel et al. Leukemia. 2006, 20: 146-147. 10.1038/sj.leu.2404042.View ArticlePubMedGoogle Scholar
- Cost GJ, Boeke JD: Targeting of human retrotransposon integration is directed by the specificity of the L1 endonuclease for regions of unusual DNA structure. Biochemistry. 1998, 37: 18081-93. 10.1021/bi981858s.View ArticlePubMedGoogle Scholar
- Yang Z: PAML: a program package for phylogenetic analysis by maximum likelihood. Comput Appl Biosci. 1997, 13: 555-6.PubMedGoogle Scholar
- Goldman N, Yang Z: A codon-based model of nucleotide substitution for protein-coding DNA sequences. Mol Biol Evol. 1994, 11: 725-36.PubMedGoogle Scholar
- Lynch M, Conery JS: The evolutionary fate and consequences of duplicate genes. Science. 2000, 290: 1151-5. 10.1126/science.290.5494.1151.View ArticlePubMedGoogle Scholar
- Bos JL: ras oncogenes in human cancer: a review. Cancer Res. 1989, 49: 4682-9.PubMedGoogle Scholar
- Corradini P, Ladetto M, Voena C, Palumbo A, Inghirami G, Knowles DM, Boccadoro M, Pileri A: Mutational activation of N- and K-ras oncogenes in plasma cell dyscrasias. Blood. 1993, 81: 2708-13.PubMedGoogle Scholar
- Barbacid M: ras genes. Annu Rev Biochem. 1987, 56: 779-827. 10.1146/annurev.bi.56.070187.004023.View ArticlePubMedGoogle Scholar
- Vossler MR, Yao H, York RD, Pan MG, Rim CS, Stork PJ: cAMP activates MAP kinase and Elk-1 through a B-Raf- and Rap1-dependent pathway. Cell. 1997, 89: 73-82. 10.1016/S0092-8674(00)80184-1.View ArticlePubMedGoogle Scholar
- Nassar N, Horn G, Herrmann C, Block C, Janknecht R, Wittinghofer A: Ras/Rap effector specificity determined by charge reversal. Nat Struct Biol. 1996, 3: 723-9. 10.1038/nsb0896-723.View ArticlePubMedGoogle Scholar
- Brinkmann T, Daumke O, Herbrand U, Kuhlmann D, Stege P, Ahmadian MR, Wittinghofer A: Rap-specific GTPase activating protein follows an alternative mechanism. J Biol Chem. 2002, 277: 12525-31. 10.1074/jbc.M109176200.View ArticlePubMedGoogle Scholar
- Spoerner M, Herrmann C, Vetter IR, Kalbitzer HR, Wittinghofer A: Dynamic properties of the Ras switch I region and its importance for binding to effectors. Proc Natl Acad Sci USA. 2001, 98: 4944-9. 10.1073/pnas.081441398.PubMed CentralView ArticlePubMedGoogle Scholar
- Polakis P, McCormick F: Structural requirements for the interaction of p21ras with GAP, exchange factors, and its biological effector target. J Biol Chem. 1993, 268: 9157-60.PubMedGoogle Scholar
- Spaargaren M, Bischoff JR: Identification of the guanine nucleotide dissociation stimulator for Ral as a putative effector molecule of R-ras, H-ras, K-ras, and Rap. Proc Natl Acad Sci USA. 1994, 91: 12609-13. 10.1073/pnas.91.26.12609.PubMed CentralView ArticlePubMedGoogle Scholar
- Bos JL: Linking Rap to cell adhesion. Curr Opin Cell Biol. 2005, 17: 123-8. 10.1016/j.ceb.2005.02.009.View ArticlePubMedGoogle Scholar
- Arthur WT, Quilliam LA, Cooper JA: Rap1 promotes cell spreading by localizing Rac guanine nucleotide exchange factors. J Cell Biol. 2004, 167: 111-22. 10.1083/jcb.200404068.PubMed CentralView ArticlePubMedGoogle Scholar
- Baertsch R, Diekhans M, Kent WJ, Haussler D, Brosius J: Retrocopy contributions to the evolution of the human genome. BMC Genomics. 2008, 9: 466-10.1186/1471-2164-9-466.PubMed CentralView ArticlePubMedGoogle Scholar
- Elliman SJ, Wu I, Kemp DM: Adult tissue-specific expression of a Dppa3-derived retrogene represents a postnatal transcript of pluripotent cell origin. J Biol Chem. 2006, 281: 16-9. 10.1074/jbc.C500415200.View ArticlePubMedGoogle Scholar
- Bradley J, Baltus A, Skaletsky H, Royce-Tolland M, Dewar K, Page DC: An X-to-autosome retrogene is required for spermatogenesis in mice. Nat Genet. 2004, 36: 872-6. 10.1038/ng1390.View ArticlePubMedGoogle Scholar
- Kitayama H, Matsuzaki T, Ikawa Y, Noda M: Genetic analysis of the Kirsten-ras-revertant 1 gene: potentiation of its tumor suppressor activity by specific point mutations. Proc Natl Acad Sci USA. 1990, 87: 4284-8. 10.1073/pnas.87.11.4284.PubMed CentralView ArticlePubMedGoogle Scholar
- Lafuente EM, van Puijenbroek AA, Krause M, Carman CV, Freeman GJ, Berezovskaya A, Constantine E, Springer TA, Gertler FB, Boussiotis VA: RIAM, an Ena/VASP and Profilin ligand, interacts with Rap1-GTP and mediates Rap1-induced adhesion. Dev Cell. 2004, 7: 585-95. 10.1016/j.devcel.2004.07.021.View ArticlePubMedGoogle Scholar
- Raaijmakers JH, Bos JL: Specificity in Ras and Rap signaling. J Biol Chem. 2009, 284: 10995-9. 10.1074/jbc.R800061200.PubMed CentralView ArticlePubMedGoogle Scholar
- Takahashi K, Mitsui K, Yamanaka S: Role of ERas in promoting tumour-like properties in mouse embryonic stem cells. Nature. 2003, 423: 541-5. 10.1038/nature01646.View ArticlePubMedGoogle Scholar
- Kameda T, Thomson JA: Human ERas gene has an upstream premature polyadenylation signal that results in a truncated, noncoding transcript. Stem Cells. 2005, 23: 1535-40. 10.1634/stemcells.2005-0054.View ArticlePubMedGoogle Scholar
- Farrell FX, Ohmstede CA, Reep BR, Lapetina EG: cDNA sequence of a new ras-related gene (rap2b) isolated from human platelets with sequence homology to rap2. Nucleic Acids Res. 1990, 18: 4281-10.1093/nar/18.14.4281.PubMed CentralView ArticlePubMedGoogle Scholar
- Posada D, Buckley TR: Model selection and model averaging in phylogenetics: advantages of akaike information criterion and bayesian approaches over likelihood ratio tests. Syst Biol. 2004, 53: 793-808. 10.1080/10635150490522304.View ArticlePubMedGoogle Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.View ArticlePubMedGoogle Scholar
- Shimodaira H, Hasegawa M: CONSEL: for assessing the confidence of phylogenetic tree selection. Bioinformatics. 2001, 17: 1246-7. 10.1093/bioinformatics/17.12.1246.View ArticlePubMedGoogle Scholar
- Mariner DJ, Davis MA, Reynolds AB: EGFR signaling to p120-catenin through phosphorylation at Y228. J Cell Sci. 2004, 117: 1339-50. 10.1242/jcs.01001.View ArticlePubMedGoogle Scholar
- Medema RH, de Laat WL, Martin GA, McCormick F, Bos JL: GTPase-activating protein SH2-SH3 domains induce gene expression in a Ras-dependent fashion. Mol Cell Biol. 1992, 12: 3425-30.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.