This Article Abstract Full Text (PDF) Supplemental Research Data All Versions of this Article: gr.6320607v1 17/8/1139    most recent Alert me when this article is cited Alert me if a correction is posted Citation Map Services Email this article to a friend Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Krull, M. Articles by Schmitz, J. Search for Related Content PubMed PubMed Citation Articles by Krull, M. Articles by Schmitz, J. Social Bookmarking                   What's this?

Letter

Functional persistence of exonized mammalian-wide interspersed repeat elements (MIRs)

¹ Institute of Experimental Pathology (ZMBE), University of Münster, Münster, Germany; ² Institute of Molecular Evolutionary Genetics and Department of Biology, Pennsylvania State University, University Park, Pennsylvania 16802, USA; ³ Institute of Bioinformatics, University of Münster, Münster, Germany

	Abstract

Top Abstract Results Discussion Conclusion Methods Acknowledgments References

 
Exonization of retroposed mobile elements, a process whereby new exons are generated following changes in non-protein-coding regions of a gene, is thought to have great potential for generating proteins with novel domains. Our previous analysis of primate-specific Alu-short interspersed elements (SINEs) showed, however, that during their 60 million years of evolution, SINE exonizations occurred in some primates, only to be lost again in some of the descendent lineages. This dynamic gain and loss makes it difficult to ascertain the contribution of exonization to genomic novelty. It was speculated that Alu-SINEs are too young to reveal persistent protein exaptation. In the present study we examined older mobile elements, mammalian-wide interspersed repeats (MIRs) that underwent active retroposition prior to the placental mammalian radiation 130 million years ago, to determine their contribution to protein-coding sequences. Of 107 potential cases of MIR exonizations in human, an analysis of splice sites substantiates a mechanism that benefits from 3' splice site selection in MIR sequences. We retraced in detail the evolution of five MIR elements that exonized at different times during mammalian evolution. Four of these are expressed as alternatively spliced transcripts; three in species throughout the mammalian phylogenetic tree and one solely in primates. The fifth is the first experimentally verified, constitutively expressed retroposed SINE element in mammals. This pattern of highly conserved, alternatively and constitutively spliced MIR sequences evinces the potential of exonized transposed elements to evolve beyond the transient state found in Alu-SINEs and persist as important parts of functional proteins.

Following the path of Alu exonizations along a phylogenetic tree of primates indicates a dynamic gain and loss of exonizations, processes that may embrace >60 million years (Myr) of evolution (Krull et al. 2005). Thus, it appears that the time frame of primate evolution might not yet be sufficiently long to ascertain whether exonized elements have been exapted as persistent modules of the proteome (Gotea and Makalowski 2006). Although transcribed, exonized sequences are destined to negotiate a grueling course along the way to becoming part of functional proteins that only a very few ever survive. Various mRNA surveillance mechanisms might degrade alternatively exonized transcripts before they become essential parts of functioning proteins (Wagner and Lykke-Andersen 2002; Hillman et al. 2004). Survivors of the degradation mechanism are further exposed to natural selection in maintaining their status as alternative splice products or in replacing the ancestral splice variant. A number of Alu exonizations in primates might not have reached that status yet. A retrospective view of much older retropositions, ones active during the Mesozoic era of mammalian evolution, promises a more definitive picture of the exonization process, especially regarding their contribution to protein plasticity.

Good candidates for illuminating these older exonization processes are the retroposed mammalian-wide interspersed repeat (MIR) elements. MIR elements amplified 130 million years ago (Mya), and they number, for example, in human, 368,000 discernible copies (Lander et al. 2001). Their age and their abundant distribution in mammals make MIRs an excellent source to elucidate the history of their exonizations. MIRs are usually truncated at either or both ends (Smit and Riggs 1995). The 5' tRNA part of the MIR is fused to a tRNA-unrelated sequence, and the 3' end features a 50-nt fragment that is similar to the 3' end of long interspersed elements (LINEs), a likely binding site of the LINE reverse transcriptase needed for retroposition (Kapitonov et al. 2004). The conserved central domain includes a 15-nt core sequence (Fig. 1). However, it is not known why the core region is highly conserved in so-called CORE-SINEs like MIRs. Possibly, CORE-SINEs have some as yet not readily understood genomic function (Gilbert and Labuda 1999; Nishihara et al. 2006). MIR elements underwent active retroposition before monotremes and placental mammals diverged, but this activity was silenced after the required LINE-specific reverse transcriptase was lost.

View larger version (31K): [in this window] [in a new window]   Figure 1. Distribution and orientation of exonized MIR sequences located in protein-coding sequences (CDS). The 107 human exonized MIR sequence regions are aligned against a schematic of a MIR consensus element. The tRNA-related part, including the internal promoter A and B boxes (black), is shown without shading (white). The region shown in dark gray comprises the 70 nt conserved central domain including a highly conserved 15-bp core sequence. The location of the MIR internal, antisense cryptic splice site AG is indicated on the sense strands (CT) by a vertical line. The 9-nt natural 3' MIR splice site (5'-ATTTTACAG-3') is shown as the inverse consensus sequence (Supplemental Fig. S3). The exonized MIR regions are represented as black lines; intronic or untranslated region (UTR) portions of the MIRs are not shown. The five experimentally analyzed examples are indicated (CHRNA1, Zfp384, ZNF639, LAS1L, and NTRK3).

	Results

Top Abstract Results Discussion Conclusion Methods Acknowledgments References

 
Selection of the data set and examples
To identify MIR elements in protein-coding sequences of human, mouse, and rat, we screened a compilation of mammalian mRNAs presumably harboring transposable element-cassettes assembled by Makalowski and co-workers (Genomic ScrapYard; http://warta.bio.psu.edu/SYDB/database.html), performing separate searches for the listed species. From 372 ScrapYard records with potential MIR element-cassettes identified in this initial search, 126 MIR sequences were present in protein-coding sequences (CDS), the remaining were either redundant entries or otherwise artifactual (Supplemental Table S1; Supplemental Fig. S1). Of these 126 loci, 107 were found in human (one of which was initially identified in rat) and were used to investigate the distribution and orientation of exonized MIR sequences (Fig. 1).

Five of the above 126 cases were (1) supported by ESTs or other indications of expression, (2) flanked by conserved sequence regions facilitating mammalian wide PCR amplification of fragments not exceeding 2 kb, (3) expressed in available tissue, and (4) embedded in introns, and thereby suitable for intensive phylogenetic reconstruction (Fig. 2; Supplemental Data Set S1). The five MIR cassettes were found in the following genes: (1) neurotrophic tyrosine kinase receptor type 3 gene (NTRK3; GenBank accession no. BT007291), (2) zinc finger protein 639 gene (ZNF639; NM_016331), (3) LAS1-like gene (LAS1L; NM_031206), (4) zinc finger protein 384 gene (Zfp384; AF216807), and (5) cholinergic receptor nicotinic alpha 1 gene (CHRNA1; NM_001039523). The mammalian-wide evolutionary distribution and state of expression of these exonized sequences is summarized in Figure 3.

View larger version (28K): [in this window] [in a new window]   Figure 2. Structures of five selected genes harboring internal exonized MIR elements. Numbered, thick gray boxes represent protein-coding sequences (CDS). The 5' and 3' untranslated regions (UTRs) are shown as medium thick gray bars. Introns are indicated as black lines; double slashes denote gaps in larger introns. Orientations and recognizable lengths of the MIR elements are indicated by arrows. The short white boxes represent the intronic parts of the MIR elements, and the adjacent tall white boxes the exonized regions. In three cases the exonization exceeds the MIR boundaries and includes anonymous intronic sequences (black bars adjacent to the exonized regions). (A) Human neurotrophic tyrosine kinase receptor type 3 gene (NTRK3; BT007291). (B) Human zinc finger protein 639 gene (ZNF639; NM_016331). (C) Human LAS1-like gene (LAS1L; NM_031206). (D) Rat zinc finger protein 384 gene (Zfp384; AF216807). (E) Human cholinergic receptor nicotinic alpha1 gene (CHRNA1; NM_001039523). The substructures in exon 7 of ZNF639 and in exons 6 and 7 of Zfp384 represent sequences encoding zinc finger domains.

View larger version (49K): [in this window] [in a new window]   Figure 3. Phylogenetic tree of investigated mammals showing the species distributions of the five internal MIR exonizations. In DNA, "+" indicates the presence of the intronic MIR element leaving the splice sites and reading frame intact, "(+)" its presence, without leaving the reading frame intact, and "–" its absence. In RNA, "+" indicates the presence of the exonized MIR sequence in cDNA derived by RT-PCR, and "–" the presence of a cDNA that does not include the exonized MIR sequence. Both symbols indicate that both splice forms are expressed. EST, expressed sequence tag, "+" including or "–" without the exonized MIR region, available at GenBank (http://www.ncbi.nlm.nih.gov/BLAST/Blast.cgi). Empty boxes indicate that the corresponding sequence information is not available or not PCR amplifiable; n/a, not analyzed. 1A third splice variant exists, derived from the human MegaMan Human Transcriptome Library (Stratagene). 2Five splice variants exist, three of which include the exonized MIR element. 3RT-PCR performed in the nine-banded armadillo instead of sloth. 4Not shown are "+" symbols under the DNA for chimpanzee and orangutan. 5Unrelated exonized intronic sequence in chicken and a corresponding exonization in the ostrich.

Twenty of the 64 antisense exonized MIR elements feature a MIR-contributed AG splice site that is preceded by a MIR-contributed oligopyrimidine tract (Fig. 1). This configuration is similar to the prevalent 3' splice site in the right arm of antisense-oriented Alu elements (Makalowski et al. 1994; Lev-Maor et al. 2003). Interestingly, the MIR-contributed splice sites are located just 3' adjacent to their highly conserved core regions (Fig. 1). The prevalent singular consensus splice site for the 20 antisense exonized MIR elements is 5'-ATTTTACAG-3' (Fig. 1). The splice site is similar to the major dual proximal and distal Alu consensus splice sites (5'-TTTTTTGAG-3' and 5'-TTGAGACAG-3', respectively) (Lev-Maor et al. 2003). However, we detected no other prevalent splice sites in MIR sequences. Two of the phylogenetically analyzed loci (ZNF639 and LAS1L) display the MIR-contributed splice site plus the oligopyrimidine tract.

Alternative expression of MIR sequences
Three of the five experimentally analyzed MIR elements (NTRK3, LAS1L, Zfp384) exhibit stable alternative splice forms in representative species of the major mammalian branches evidenced both in DNA (conservation of splice sites and maintenance of ORFs) and mRNA, demonstrating that the theoretical splice variants actually exist (Fig. 3). In the case of CHRNA1, stable alternative splice forms in mRNA were found in human and gorilla. However, available DNA sequence data also suggest the exonization of the intronic MIR element in chimpanzee and orangutan but not in gibbon, mandrill, or marmoset, suggesting a great ape specificity of this exonization. Demonstrable expression in human and gorilla and conservation in other great apes indicate selection pressure and probably functionality for this exonization >100 Myr after insertion (Supplemental Fig. S2). Because of the low sequence divergence between exonized MIRs of the compared species, calculated K_a/K_s values, frequently used to describe functional selection pressure for alternative splicing events (Xing and Lee 2005), were not meaningful for these examples (Supplemental Data Set S2).

Constitutive expression of a MIR sequence
Zhang and Chasin (2006) proposed that advantageous new exons are eventually expressed constitutively. For SINEs, we provide the first evidence of this prediction. Of the five genes containing internal MIR exonizations, ZNF639 is the only one in which the MIR sequence is constitutively exonized in protein-coding sequences. All analyzed mammalian representatives of Eutheria (Supraprimates: human and mouse; Laurasiatheria: mole, dog, and cow; Afrotheria: manatee; and Xenarthra: sloth), Metatheria (Didelphimorphia: opossum), and Monotremata (Ornithorhynchidae: platypus) express ZNF639 mRNAs only with, and none without, the MIR contribution. The K_a/K_s value of 0.19 indicates purifying selection. The exonized MIR component encodes 45 amino acids located within the 205 amino acid N-terminal part of the ZNF639 protein. This region (amino acids 58–102; Imoto et al. 2003) contains no recognizable protein motifs, while the 280 C-terminal amino acids contain nine zinc finger motifs (Fig. 2). The protein is exclusively expressed in cell nuclei and probably has transcription repressor activity (Bogaerts et al. 2005). ZNF639 is overexpressed in squamous cell carcinomas and over-expression correlates with a shorter survival of these patients. However, the stable MIR inclusion in the 56-kDa protein in all mammalian lineages indicates a significant function of the corresponding 45 additional amino acids.

For ZNF639, all necessary events from insertion of the MIR element to recruitment of parts of the MIR as a novel protein-coding exon occurred on the phylogenetic branch leading from the common ancestor of amniotes (mammals, reptiles, dinosaurs, and birds) to the mammalian ancestor. However, there are no living animals that diverged in this time period that would enable us to reconstruct, step by step, the successive evolution of molecular changes that were necessary to facilitate the 5' functional splice site and an intact ORF, or to show possible intermediate alternative splice variants.

	Discussion

Top Abstract Results Discussion Conclusion Methods Acknowledgments References

 
Mammalian-wide detection of MIR sequence exonizations
From a total of 372 potential exonized MIR elements, we investigated 126 elements with strong indications of exonization in human, mouse, or rat (Supplemental Fig. S1). On five of these, we performed an extensive retrospective analysis reconstructing >100 Myr of mammalian history. To establish a complete mammalian exonization pattern, we limited our analyses to those loci that could be amplified in the greatest number of species, sampling for both DNA and RNA analyses (introns 2 kb). From our extensive experience analyzing 160,000 genomic loci in the evolutionary history of mammalian species (Kriegs et al. 2006), we concluded that larger introns in some compared species render PCR amplification and comparative analyses in the rest of the species difficult. However, in addition to the five experimentally accessible loci presented in this study, we attempted to amplify additional exonizations with larger intronic regions (up to 5.8 kb; Supplemental Table S1, e.g., NM_145065, NM_015424, NM_010058). None of these cases was amplifiable in all essential key mammalian species. To gain full species sampling for those examples and to add some additional loci with lower conservation of flanking regions, we must await the upcoming data from genome sequencing projects.

Expression of MIR sequences
Relative to gene orientation, significantly more MIR exonized sequences are located in the antisense orientation (60%); this preference is even more pronounced for Alu exonizations (85%; Sorek et al. 2002). The difference might be explained by the significant predominance of antisense intronic Alu insertions compared to MIR insertions that are found preferentially in the sense orientation. Furthermore, there are two prevalent 3' splice sites in the right arm of exonized Alu sequences (proximal and distal; Lev-Maor et al. 2003), whereas we could detect only one such prevalent 3' splice site for the potential MIR exonizations (Fig. 1; Supplemental Figs. S2, S3).

An intriguing observation is that the natural, MIR-specific splice site is just 3' adjacent to the MIR conserved core region. However, there are hundreds of thousands of MIR elements that are not associated with splice sites but still bear the conserved core region. Given their origins dating back at least 130 Mya, it is still unclear why, in contrast to the rest of the MIR sequence, the core sequences remain so highly conserved. In coding sequences the 5' domain of MIR elements appears to be preferentially exonized (Fig. 1). On the one hand, this bias may be due to the naturally occurring splice site in the consensus sequence, but it also might reflect a more efficient recognition of exonized conserved 5' MIR parts compared with more diverged, but possibly as well exonized, 3' parts by the RepeatMasker. The documented persistence of MIR element exonizations appears to be different from that of Alu-SINEs. In our detailed analysis of selected Alu exonizations, we found examples in which apparently functional splice sites in ancestral species were not conserved in all related primate lineages, indicating a loss of functional splice forms (Krull et al. 2005). This implies that many Alu-derived exons still have a selection status close to neutrality; in other words, they are not yet indispensable. Due to the relatively short time that has transpired since their exonization, the alternative exons might not yet be capable of making a significant evolutionary contribution. Gotea and Makalowski (2006) propose that such "young" exonized sequences might have functions other than protein coding, such as regulatory functions, including nonsense-mediated decay (NMD) to control gene expression. Another explanation may be the loss of temporarily exapted (historically constrained) alternative splice forms in response to divergent functional needs in different taxa, comparable to the loss of historically constrained duplicated genes (Greenberg et al. 2006).

Time point of MIR exonization
MIR elements were actively mobile prior to the mammalian radiation 130 Mya (Smit and Riggs 1995). In four of the five phylogenetically analyzed cases, exonization of MIR elements also took place very early in the evolution of mammals (Fig. 3). However, in the case of CHRNA1 we found a more recent, great ape-specific exonization pattern. Although Saini et al. (2005) suggested, based on RT-PCR fragment size not verified by sequencing, that the CHRNA1 isoform (containing the MIR element) is also present in mouse, we can clearly exclude its expression in mouse. By sequencing both DNA and RNA, we clearly observed that the human orthologous exonized MIR region is missing in mouse (Supplemental Fig. S2; Supplemental Data Set S3). This underscores once more that exaptation in general and exonization in particular of retroposed sequences, can occur at any age and consequently at any stage of decay (Brosius 2005). In the other four cases, five of eight splice sites were part of the MIR consensus sequence, a fact that may have facilitated their exonization immediately after integration. Depending on the evolutionary time point of exonization or the degree of initial positive selection, the exonized part of the retroposed SINE more or less reflects the ancestral sequence of the MIR element.

Nearly 30 years ago, Walter Gilbert recognized in alternative splicing a process that allows evolution to try out new solutions without destroying the old (Gilbert 1978), a variation of a "strategy" that had been recognized in gene and genome duplication (Bridges 1936; Ohno 1970). Particularly in mammalian genomes, exonized transposed elements contribute significantly to alternative splicing.

Evolutionary time seems to be a critical factor in establishing essential key mutations required for exonization. Three examples endorse this assumption (Fig. 4): (1) Alu exonizations that are present in specific primate lineages but not in others (RPE2–1, C-rel-2, MTO1–3, survivin) (Krull et al. 2005; Mola et al. 2007). A relatively short evolutionary period (<60 Myr) led to exonizations, the gain and loss of which fluctuated greatly in various species of given taxonomic groups. During early stages of exonization, the newly acquired exons (perhaps present only in a small percentage of the alternative transcripts) probably ranged from slightly detrimental to slightly advantageous. In any event, selection pressure was perhaps either not present or too modest for persistence of the exaptation. (2) MIR exonizations (NTRK3, LAS1L, and Zfp384) are present in all mammals. Selection over longer time periods favored exonizations that are now stably established and expressed in all representative taxa. In many cases, both the ancestral and the new alternative splice forms were advantageous and therefore expression of both forms persisted. (3) The MIR exonization in ZNF639 is constitutively expressed in all analyzed mammalian species. Expression shifted completely, possibly via co-existence of two alternative splice forms for a time, from 100% of the ancient form (prior to MIR exonization) to 100% of the form harboring the MIR exon, such that the latter is now constitutive. We cannot discern whether the constitutive inclusion of the exonized sequence in ZNF639 was achieved via a transient alternative splicing or directly after acquisition of the necessary splice sites (Fig. 4).

View larger version (22K): [in this window] [in a new window]   Figure 4. Evolutionary scenario after exonization of sequences, from the ancestral constitutive to alternative and, occasionally, to novel constitutive splicing. Intensities of selection are indicated at the left and right bars. For evolution of alternative splice sites, the intensities of selection vary over time. Initially, inclusion of the novel exon might be slightly deleterious, neutral, or slightly advantageous. Subsequently, selection changes, perhaps during a period of positive selection to increasing levels of negative selection (bar at the left with a gradient from white to red). The numeration refers to intronic insertion of a transposed element (red box) (1), followed by acquisition of necessary components for alternative splicing (2). Under relaxed selective pressure, the exonized condition might revert or become fixed in different lineages (vertical double arrow). In an extreme case, (3), as shown for ZNF639, the novel alternative splice form might, over time, completely replace the ancestral constitutive splice form, thus representing a different constitutive form (bar at the right side, with gradients from black via white to red). However, the constitutive inclusion of the exonized sequence in ZNF639 could also have been acquired directly after acquisition of the necessary splice sites. Examples corresponding to different stages of exonization taken from previous works (Singer et al. 2004; Krull et al. 2005; Mola et al. 2007) or from the present work (bold letters) are presented in ovals.

Although K_a/K_s values indicate that the exonized part of the MIR element itself is under moderate selection pressure, the 291-nt adjacent protein coding flanks (parts of exons 5 and 7) show even lower K_a/K_s values (0.17 vs. 0.19 for the exonized MIR). This, at most, suggests a possibly lower selection pressure on the MIR exonized sequence than on the flanks. This difference is even greater when compared to the nine functional zinc finger domains of ZNF639 in exon 7 (582 nt; K_a/K_s = 0.03; data not shown). However, more information about functional domains of the N-terminal region of the protein is necessary to present more conclusive information about a potential spacer function of the exonized MIR sequence. There is also another report of two independent exonizations, although of different lengths and origin, in the same intron of the ADARB1 gene in different taxonomic groups (human and mouse; Slavov and Gardiner 2002).

In theory, the first transition from random insertion to alternative splicing is reversible and is either not at all under purifying selection or under relaxed negative selection (Xing and Lee 2006). Novel exons might be generated and shaped for "testing" under such relaxed selection (this period might include a phase of positive selection as well). Through a gradient of purifying selection, the second transition leads to stable alternative splice forms and in rare cases to constitutive splicing that includes the exonized form. ZNF639 is the first gene experimentally verified to contain such a constitutively expressed MIR exonization and exaptation. Of course, at this time we cannot exclude the possibility that the ZNF639 exonization was constitutive from the start and did not pass through this transitional state.

	Conclusion

Top Abstract Results Discussion Conclusion Methods Acknowledgments References

 
The contribution of transposed elements to gene structures is more or less coincidental. Their persistence is usually transient. If not deleted, they fade beyond recognition over longer evolutionary periods. However, a notable fraction of transposed elements escapes transience, for example, by integrating into protein-coding parts of genes, facilitated by internal components providing splice sites and oligopyrimidine tracts and "reprogramming" the splicing system of a targeted gene. Once proven worthy in the struggle of survival, they endure recognizably over hundreds of millions of years and contribute to significant tasks. We have identified and analyzed some of these candidates, thus shedding light on their >100 million-year-old evolutionary histories that show they have clearly stood the test of evolutionary time and persisted in mammalian lineages. We showed that 3' splice site selection in exonized transposed elements is not restricted to Alu elements but seems to be an older and significant mechanism for MIR exonization as well. Alternative splicing of exonized MIRs is exemplarily shown to be a stable process retained >100 Myr in all major groups of mammals. Functional persistence shown by constitutive splicing of an exonized MIR sequence was evidenced for the first time and demonstrates that this evolutionary pathway is not necessary correlated with genetic disorder, as has been suggested for Alu exonizations by Lev-Maor et al. (2003). The present data are ample evidence of the value of exonization, given enough evolutionary time, in generating genomic novelties.

	Methods

Top Abstract Results Discussion Conclusion Methods Acknowledgments References

 
Data selection
To identify MIR elements in protein-coding sequences, we screened a compilation of mammalian mRNAs with transposable element cassettes (Genomic ScrapYard; http://warta.bio.psu.edu/SYDB/database.html), performing independent searches for human, mouse, and rat. Out of 1091, 472, and 201 matching records we selected 314, 38, and 20 cases that included MIR elements in human, mouse, and rat, respectively. These 372 cases were then scrutinized to filter out duplications and other artifacts (246 cases; Supplemental Fig. S1). To search for cases that were suitable for experimental evaluation of their phylogeny, we screened the remaining 126 potential exonizations by applying the following criteria: (1) The MIR-derived sequence should not be located in the first or last protein-coding exon as they usually lack highly conserved flanking regions and hence are difficult to amplify by polymerase chain reaction (PCR). (2) To enable manageable PCRs of up to 2 kb in highly diverged mammalian species, the MIR-derived sequences should be flanked by conserved sequences. Conservation was determined by comparison of available genomic sequences (e.g., of human, mouse, and dog). (3) Indication for exonization should be supplied by such available information as EST data or other published information. (4) Relevant tissues should be available for specific alternative transcripts. Five of the 126 potential exonizations fulfilled all four criteria and were subjected to further phylogenetic examination using PCR, RT-PCR, and sequence analyses. Note that the ScrapYard database consists of GenBank entries from January 15, 2002. It is expected that an updated database would facilitate the recovery of additional cases of MIR exonizations.

Available sequence information was obtained from the UCSC Genome Browser (http://genome.ucsc.edu/cgi-bin/hgBlat) and the NCBI trace archive (http://www.ncbi.nlm.nih.gov/blast/tracemb.shtml).

Experimental procedures concerning DNA and RNA extraction, PCR amplification, and reverse transcription are given in the Supplemental Protocol S1. PCR primers are listed in Supplemental Table S2 and are illustrated in Supplemental Fig. S5.

Sequence analyses
For detection and classification of the inserted and partially exonized MIR sequences, we used the RepeatMasker server (A.F.A. Smit, R. Hubley, and P. Green, unpubl.; RepeatMasker at http://www.repeatmasker.org/cgi-bin/WEBRepeatMasker). The available mammalian sequences of the five relevant loci were derived from the NCBI database, and all additionally amplified sequences were manually aligned. GenBank entries are compiled in Supplemental Table S3.

K_a/K_s values
We used the Yang–Nielsen maximum likelihood method (Yang and Nielsen 2000), implemented in the YN00 program of the PAML package (Yang 1997), to calculate the K_a/K_s values for aligned exonized parts of the MIR elements. K_a/K_s values were averaged over all investigated species pairs. It should be noted that in most cases the significance of the derived values was limited by the short length and low variation of the exonized regions.

	Acknowledgments

Top Abstract Results Discussion Conclusion Methods Acknowledgments References

 
We thank Frank Grützner, Rodney L. Honeycutt, Uwe Joite, Jan Ole Kriegs, Jörg Molten, Bernhard Neurohr, Christian Roos, Gertrud Scheele, Heike Weber, and Anja Zemann for providing us with tissue samples and Marsha Bundman for editorial assistance. We thank Valer Gotea for his help in selecting data from the Genomic ScrapYard database and Michael Haberl for his comments. We thank Django Sussman for introducing us to methods for analyzing structural features of the ZNF639 protein. J.S. thanks Matthias Schmitz for all his personal support. This work was supported by the Nationales Genomforschungsnetz (NGFN) (0313358A to J.B. and J.S.), the European Union (EU) (LSHG-CT-2003-503022 to J.B.), and the Deutsche Forschungsgemeinschaft (DFG) (SCHM1469 to J.S. and J.B.).

	Footnotes

 
⁴ Corresponding authors.

E-mail jueschm{at}uni-muenster.de; fax 49-251-8352134.

E-mail RNA.world{at}uni-muenster.de; fax 49-251-8358512.

[Supplemental material is available online at www.genome.org. The sequence data from this study have been submitted to NCBI under accession nos: DQ323592–DQ323661, DQ507223–DQ507235, DQ855908–DQ855913, EF418572, EF422277, EF520683, and EF520684.]

Article published online before print. Article and publication date are at http://www.genome.org/cgi/doi/10.1101/gr.6320607

	References

Top Abstract Results Discussion Conclusion Methods Acknowledgments References

 

Bejerano, G., Lowe, C.B., Ahituv, N., King, B., Siepel, A., Salama, S.R., Rubin, E.M., Kent, W.J., and Haussler, D. 2006. A distal enhancer and an ultraconserved exon are derived from a novel retroposon. Nature 441: 87–90.[CrossRef][Medline]

Bogaerts, S., Vanlandschoot, A., van Hengel, J., and van Roy, F. 2005. Nuclear translocation of N-catenin by the novel zinc finger transcriptional repressor ZASC1. Exp. Cell Res. 311: 1–13.[CrossRef][Medline]

Bridges, C.B. 1936. Genes and chromosomes. Teaching Biol. Nov: 17–23.

Brosius, J. 2005. Echoes from the past—Are we still in an RNP world? Cytogenet. Genome Res. 110: 8–24.[CrossRef][Medline]

Brosius, J. and Gould, S.J. 1992. On "genomenclature": a comprehensive (and respectful) taxonomy for pseudogenes and other "junk DNA". Proc. Natl. Acad. Sci. 89: 10706–10710.[Abstract/Free Full Text]

Gilbert, W. 1978. Why genes in pieces? Nature 271: 501.[CrossRef][Medline]

Gilbert, N. and Labuda, D. 1999. CORE-SINEs: Eukaryotic short interspersed retroposing elements with common sequence motifs. Proc. Natl. Acad. Sci. 96: 2869–2874.[Abstract/Free Full Text]

Gotea, V. and Makalowski, W. 2006. Do transposable elements really contribute to proteomes? Trends Genet. 22: 260–267.[CrossRef][Medline]

Greenberg, A.J., Moran, J.R., Fang, S., and Wu, C.-I. 2006. Adaptive loss of an old duplicated gene during incipient speciation. Mol. Biol. Evol. 23: 401–410.[Abstract/Free Full Text]

Hillman, R.T., Green, R.E., and Brenner, S.E. 2004. An unappreciated role for RNA surveillance. Genome Biol. 5: R8. doi: 10.1186/gb-2004-5-2-r8.[CrossRef][Medline]

Imoto, I., Yuki, Y., Sonoda, I., Ito, T., Shimada, Y., Imamura, M., and Inazawa, J. 2003. Identification of ZASC1 encoding a Krüppel-like zinc finger protein as a novel target for 3q26 amplification in esophageal squamous cell carcinomas. Cancer Res. 63: 5691–5696.[Abstract/Free Full Text]

Kapitonov, V.V., Pavlicek, A., and Jurka, J. 2004. Anthology of human repetitive DNA. In Encyclopedia of molecular cell biology and molecular medicine (ed. R.A. Meyers), pp. 251–305. Wiley-VCH, Weinheim.

Kriegs, J.O., Churakov, G., Kiefmann, M., Jordan, U., Brosius, J., and Schmitz, J. 2006. Retroposed elements as archives for the evolutionary history of placental mammals. PLoS Biol. 4: e91. doi: 10.1371/journal.pbio.0040091.[CrossRef][Medline]

Krull, M., Brosius, J., and Schmitz, J. 2005. Alu-SINE exonization: En route to protein-coding function. Mol. Biol. Evol. 22: 1702–1711.[Abstract/Free Full Text]

Lander, E.S., Linton, L.M., Birren, B., Nusbaum, C., Zody, M.C., Baldwin, J., Devon, K., Dewar, K., Doyle, M., FitzHugh, W., et al. 2001. Initial sequencing and analysis of the human genome. Nature 409: 860–921.[CrossRef][Medline]

Lev-Maor, G., Sorek, R., Shomron, N., and Ast, G. 2003. The birth of an alternatively spliced exon: 3' Splice-site selection in Alu exons. Science 300: 1288–1291.[Abstract/Free Full Text]

Makalowski, W., Mitchell, G.A., and Labuda, D. 1994. Alu sequences in the coding regions of mRNA: A source of protein variability. Trends Genet. 10: 188–193.[CrossRef][Medline]

Mola, G., Vela, E., Fernández-Figueras, M.T., Isamat, M., and Munoz-Mármol, A.M. 2007. Exonization of Alu-generated splice variants in the Survivin gene of human and non-human primates. J. Mol. Biol. 366: 1055–1063.[CrossRef][Medline]

Nishihara, H., Smit, A.F.A., and Okada, N. 2006. Functional noncoding sequences derived from SINEs in the mammalian genome. Genome Res. 16: 864–874.[Abstract/Free Full Text]

Ohno, S. 1970. Evolution by gene duplication. Springer, New York.

Saini, S.S., Tüzün, E., and Christadoss, P. 2005. The cDNA of mouse skeletal muscle transcribe for both isoforms 1 and 2 of acetylcholine receptor subunit. J. Neuroimmunol. 169: 177–179.[CrossRef][Medline]

Singer, S.S., Maennel, D.N., Hehlgans, T., Brosius, J., and Schmitz, J. 2004. From "junk" to gene: Curriculum vitae of a primate receptor isoform gene. J. Mol. Biol. 341: 883–886.[CrossRef][Medline]

Slavov, D. and Gardiner, K. 2002. Phylogenetic comparison of the pre-mRNA adenosine deaminase ADAR2 genes and transcripts: Conservation and diversity in editing site sequence and alternative splicing patterns. Gene 299: 83–94.[CrossRef][Medline]

Smit, A.F. and Riggs, A.D. 1995. MIRs are classic, tRNA-derived SINEs that amplified before the mammalian radiation. Nucleic Acids Res. 23: 98–102.[Abstract/Free Full Text]

Sorek, R., Ast, G., and Graur, D. 2002. Alu-containing exons are alternatively spliced. Genome Res. 12: 1060–1067.[Abstract/Free Full Text]

Wagner, E. and Lykke-Andersen, J. 2002. mRNA surveillance: The perfect persist. J. Cell Sci. 115: 3033–3038.[Abstract/Free Full Text]

Xing, Y. and Lee, C. 2005. Evidence of functional selection pressure for alternative splicing events that accelerate evolution of protein subsequences. Proc. Natl. Acad. Sci. 102: 13526–13531.[Abstract/Free Full Text]

Xing, Y. and Lee, C. 2006. Alternative splicing and RNA selection pressure—Evolutionary consequences for eukaryotic genomes. Nat. Rev. Genet. 7: 499–509.[CrossRef][Medline]

Yang, Z. 1997. PAML: A program package for phylogenetic analysis by maximum likelihood. Comput. Appl. Biosci. 13: 555–556.[Free Full Text]

Yang, Z. and Nielsen, R. 2000. Estimating synonymous and nonsynonymous substitution rates under realistic evolutionary models. Mol. Biol. Evol. 17: 32–43.[Abstract/Free Full Text]

Zhang, X.H.-F. and Chasin, L.A. 2006. Comparison of multiple vertebrate genomes reveals the birth and evolution of human exons. Proc. Natl. Acad. Sci. 103: 13427–13432.[Abstract/Free Full Text]

Received January 24, 2007; accepted in revised format May 8, 2007.

CiteULike

Connotea

Del.icio.us

Digg

Reddit

Technorati

This article has been cited by other articles:

				J. Schmitz, A. Zemann, G. Churakov, H. Kuhl, F. Grutzner, R. Reinhardt, and J. Brosius Retroposed SNOfall--A mammalian-wide comparison of platypus snoRNAs Genome Res., June 1, 2008; 18(6): 1005 - 1010. [Abstract] [Full Text] [PDF]

				R. Sorek The birth of new exons: Mechanisms and evolutionary consequences RNA, October 1, 2007; 13(10): 1603 - 1608. [Abstract] [Full Text] [PDF]

This Article Abstract Full Text (PDF) Supplemental Research Data All Versions of this Article: gr.6320607v1 17/8/1139    most recent Alert me when this article is cited Alert me if a correction is posted Citation Map Services Email this article to a friend Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Krull, M. Articles by Schmitz, J. Search for Related Content PubMed PubMed Citation Articles by Krull, M. Articles by Schmitz, J. Social Bookmarking                   What's this?