Evolution, expression, and possible function of a master gene for amplification of an interspersed repeated DNA family in rodents

Subfamily-specific quantification of endogenous mouse L1 retrotransposons by droplet digital PCR

Long interspersed element type 1 (LINE-1; L1) mobilizes during early embryogenesis, neurogenesis, and germ cell development, accounting for 25% of disease-causing heritable insertions and 98% of somatic insertions in cancer. To better understand the regulation and impact of L1 mobilization in the genome, reliable methods for measuring L1 copy number variation (CNV) are needed. Here we present a comprehensive analysis of a droplet digital PCR (ddPCR) based method for quantifying endogenous mouse L1. We provide experimental evidence that ddPCR assays can be designed to target specific L1 subfamilies using diagnostic single nucleotide polymorphisms (SNPs). The target and off-target L1 subfamilies form distinct droplet clusters, which were experimentally verified using both synthetic gene fragments and endogenous L1 derived plasmid clones. We further provide a roadmap for in silico assay design and evaluation of target specificity, ddPCR testing, and optimization for L1 CNV quantification. The assay can achieve a sensitivity of 5% CNV with 8 technical replicates. With 24 technical replicates, it can detect 2% CNV because of the increased precision. The same approach will serve as a guide for the development of ddPCR based assays for quantifying human L1 copy number and any other high copy genomic target sequences.

Tracing the history of LINE and SINE extinction in sigmodontine rodents

Background

L1 retrotransposons have co-evolved with their mammalian hosts for the entire history of mammals and currently compose ~ 20% of a mammalian genome. B1 retrotransposons are dependent on L1 for retrotransposition and span the evolutionary history of rodents since their radiation. L1s were found to have lost their activity in a group of South American rodents, the Sigmodontinae, and B1 inactivation preceded the extinction of L1 in the same group. Consequently, a basal group of sigmodontines have active L1s but inactive B1s and a derived clade have both inactive L1s and B1s. It has been suggested that B1s became extinct during a long period of L1 quiescence and that L1s subsequently reemerged in the basal group.

Results

Here we investigate the evolutionary histories of L1 and B1 in the sigmodontine rodents and show that L1 activity continued until after the L1-extinct clade and the basal group diverged. After the split, L1 had a small burst of activity in the former group, followed by extinction. In the basal group, activity was initially low but was followed by a dramatic increase in L1 activity. We found the last wave of B1 retrotransposition was large and probably preceded the split between the two rodent clades.

Conclusions

Given that L1s had been steadily retrotransposing during the time corresponding to B1 extinction and that the burst of B1 activity preceding B1 extinction was large, we conclude that B1 extinction was not a result of L1 quiescence. Rather, the burst of B1 activity may have contributed to L1 extinction both by competition with L1 and by putting strong selective pressure on the host to control retrotransposition.

Electronic supplementary material

The online version of this article (10.1186/s13100-019-0164-5) contains supplementary material, which is available to authorized users.

Exaptation at the molecular genetic level

The realization that body parts of animals and plants can be recruited or coopted for novel functions dates back to, or even predates the observations of Darwin. S.J. Gould and E.S. Vrba recognized a mode of evolution of characters that differs from adaptation. The umbrella term aptation was supplemented with the concept of exaptation. Unlike adaptations, which are restricted to features built by selection for their current role, exaptations are features that currently enhance fitness, even though their present role was not a result of natural selection. Exaptations can also arise from nonaptations; these are characters which had previously been evolving neutrally. All nonaptations are potential exaptations. The concept of exaptation was expanded to the molecular genetic level which aided greatly in understanding the enormous potential of neutrally evolving repetitive DNA-including transposed elements, formerly considered junk DNA-for the evolution of genes and genomes. The distinction between adaptations and exaptations is outlined in this review and examples are given. Also elaborated on is the fact that such distinctions are sometimes more difficult to determine; this is a widespread phenomenon in biology, where continua abound and clear borders between states and definitions are rare.

Origin and evolution of SINEs in eukaryotic genomes

Short interspersed elements (SINEs) are one of the two most prolific mobile genomic elements in most of the higher eukaryotes. Although their biology is still not thoroughly understood, unusual life cycle of these simple elements amplified as genomic parasites makes their evolution unique in many ways. In contrast to most genetic elements including other transposons, SINEs emerged de novo many times in evolution from available molecules (for example, tRNA). The involvement of reverse transcription in their amplification cycle, huge number of genomic copies and modular structure allow variation mechanisms in SINEs uncommon or rare in other genetic elements (module exchange between SINE families, dimerization, and so on.). Overall, SINE evolution includes their emergence, progressive optimization and counteraction to the cell's defense against mobile genetic elements.

Two primate-specific small non-protein-coding RNAs in transgenic mice: neuronal expression, subcellular localization and binding partners

In a rare occasion a single chromosomal locus was targeted twice by independent Alu-related retroposon insertions, and in both cases supported neuronal expression of the respective inserted genes encoding small non-protein coding RNAs (npcRNAs): BC200 RNA in anthropoid primates and G22 RNA in the Lorisoidea branch of prosimians. To avoid primate experimentation, we generated transgenic mice to study neuronal expression and protein binding partners for BC200 and G22 npcRNAs. The BC200 gene, with sufficient upstream flanking sequences, is expressed in transgenic mouse brain areas comparable to those in human brain, and G22 gene, with upstream flanks, has a similar expression pattern. However, when all upstream regions of the G22 gene were removed, expression was completely abolished, despite the presence of intact internal RNA polymerase III promoter elements. Transgenic BC200 RNA is transported into neuronal dendrites as it is in human brain. G22 RNA, almost twice as large as BC200 RNA, has a similar subcellular localization. Both transgenically expressed npcRNAs formed RNP complexes with poly(A) binding protein and the heterodimer SRP9/14, as does BC200 RNA in human. These observations strongly support the possibility that the independently exapted npcRNAs have similar functions, perhaps in translational regulation of dendritic protein biosynthesis in neurons of the respective primates.

Dendritic BC1 RNA in translational control mechanisms

Translational control at the synapse is thought to be a key determinant of neuronal plasticity. How is such control implemented? We report that small untranslated BC1 RNA is a specific effector of translational control both in vitro and in vivo. BC1 RNA, expressed in neurons and germ cells, inhibits a rate-limiting step in the assembly of translation initiation complexes. A translational repression element is contained within the unique 3' domain of BC1 RNA. Interactions of this domain with eukaryotic initiation factor 4A and poly(A) binding protein mediate repression, indicating that the 3' BC1 domain targets a functional interaction between these factors. In contrast, interactions of BC1 RNA with the fragile X mental retardation protein could not be documented. Thus, BC1 RNA modulates translation-dependent processes in neurons and germs cells by directly interacting with translation initiation factors.

SINE extinction preceded LINE extinction in sigmodontine rodents: implications for retrotranspositional dynamics and mechanisms

Short Interspersed Nuclear Elements, or SINEs, retrotranspose despite lacking protein-coding capability. It has been proposed that SINEs utilize enzymes produced in trans by Long Interspersed Nuclear Elements, or LINEs. Strong support for this hypothesis is found in LINE and SINE pairs that share sequence homology; however, LINEs and SINEs in primates and rodents are only linked by an insertion site motif. We have now profiled L1 LINE and B1 SINE activity in 24 rodent species including candidate taxa for the first documented L1 extinction. As expected, there was no evidence for recent activity of B1s in species that also lack L1 activity. However, B1 silencing appears to have preceded L1 extinction, since B1 activity is also lacking in the genus most closely related to those lacking active L1s despite the presence of active L1s in this genus. A second genus with active L1s but inactive B1s was also identified.

The human L1 promoter: variable transcription initiation sites and a major impact of upstream flanking sequence on promoter activity

Human L1 elements are non-LTR retrotransposons that comprise approximately 17% of the human genome. Their 5'-untranslated region (5'-UTR) serves as a promoter for L1 transcription. Now we find that transcription initiation sites are not restricted to nucleotide +1 but vary considerably in both downstream and upstream directions. Transcription initiating upstream explains additional nucleotides often seen between the 5'-target site duplication and the L1 start site. A higher frequency of G nucleotides observed upstream from the L1 can be explained by reverse transcription of the L1 RNA 5'-CAP, which is further supported by extra Gs seen for full-length HERV-W pseudogenes. We assayed 5'-UTR promoter activities for several full-length human L1 elements, and found that upstream flanking cellular sequences strongly influence the L1 5'-UTR promoter. These sequences either repress or enhance the L1 promoter activity. Therefore, the evolutionary success of a human L1 in producing progeny depends not only on the L1 itself, but also on its genomic integration site. The promoter mechanism of L1 is reminiscent of initiator (Inr) elements that are TATA-less promoters expressing several cellular genes. We suggest that the L1 5'-UTR is able to form an Inr element that reaches into upstream flanking sequence.

Transcriptional regulation of the human LINE-1 retrotransposon L1.2B

Although LINE-1 (L1) sequences constitute the most important family of retrotransposons in the human genome, their transcriptional regulation is poorly understood. Specifically, their unusual internal promoter is incompletely characterized. Current promoter prediction programs fail to identify the promoter in the 5'UTR of the active LINE-1 element L1.2B. Experimental investigation of this promoter using reporter gene assays in various human and murine cell types confirmed that the promoter consists of two segments, and demonstrated that the distal portion is essential for cell-type-independent activity. No differences in promoter activity were found between normal and transformed cells. The complete promoter was shown to possess approximately 20% of the activity of the strong early promoter of cytomegalovirus, and to be capable of directing the expression of levels of p53 sufficient to kill normal and transformed human cells. Thus, active LINE-1 elements contain highly active promoters capable of driving cell-type-independent expression, which are of potential use in mammalian expression constructs. In vitro methylation of the promoter at HpaII sites decreased its activity independently of cell type, but this repression was alleviated in MBD2-/- cells. Surprisingly, mutation of specific HpaII sites was also found to reduce promoter activity. Thus, efficient repression of the L1.2B promoter by DNA methylation may involve MBD2 binding, but at least one HpaII site also appears to be involved specifically in transcriptional activation. Since neither promoter activity nor the efficiency of repression by methylation differed between normal and tumor cells, the re-activation of LINE-1 sequences observed in tumor cells is probably caused by hypomethylation of the promoter.

Neuronal untranslated BC1 RNA: targeted gene elimination in mice

Despite the potentially important roles of untranslated RNAs in cellular form or function, genes encoding such RNAs have until now received surprisingly little attention. One such gene encodes BC1 RNA, a small non-mRNA that is delivered to dendritic microdomains in neurons. We have now eliminated the BC1 RNA gene in mice. Three independent founder lines were established from separate embryonic stem cells. The mutant mice appeared to be healthy and showed no anatomical or neurological abnormalities. The gross brain morphology was unaltered in such mice, as were the subcellular distributions of two prototypical dendritic mRNAs (encoding MAP2 and CaMKIIalpha). Due to the relatively recent evolutionary origin of the gene, we expected molecular and behavioral consequences to be subtle. Behavioral analyses, to be reported separately, indicate that the lack of BC1 RNA appears to reduce exploratory activity.

Non-coding ribonucleic acids--a class of their own?

Non-coding ribonucleic acids (RNAs) do not contain a peptide-encoding open reading frame and are therefore not translated into proteins. They are expressed in all phyla, and in eukaryotic cells they are found in the nucleus, cytoplasm, and mitochondria. Non-coding RNAs either can exert structural functions, as do transfer and ribosomal RNAs, or they can regulate gene expression. Non-coding RNAs with regulatory functions differ in size ranging from a few nucleotides to over 100 kb and have diverse cell- or development-specific functions. Some of the non-coding RNAs associate with human diseases. This chapter summarizes the current knowledge about regulatory non-coding RNAs.

Binding of p190RhoGEF to a destabilizing element on the light neurofilament mRNA is competed by BC1 RNA

The enhancement of RNA-mediated motor neuron degeneration in transgenic mice by mutating a major mRNA instability determinant in a light neurofilament (NF-L) transgene implicates cognate RNA binding factors in the pathogenesis of motor neuron degeneration. p190RhoGEF is a neuron-enriched guanine exchange factor (GEF) that binds to the NF-L-destabilizing element, to c-Jun N-terminal kinase-interactive protein-1 (JIP-1), and to 14-3-3 and may link neurofilament expression to pathways affecting neuronal homeostasis. This study was undertaken to identify additional RNA species that bind p190RhoGEF and could affect interactions of the exchange factor with NF-L transcripts. The C-terminal domain of p190RhoGEF, containing the RNA-binding site, was expressed as a glutathione S-transferase fusion protein and was used as an affinity probe to isolate interactive RNAs in rat brain extracts. As expected, NF-L mRNA was identified as an RNA specie eluted from the affinity column. In addition, BC1 RNA was also found enriched in the bound RNA fraction. BC1 is a 152-nucleotide RNA that is highly expressed but untranslated in differentiated neurons. We show that BC1 and NF-L mRNA bind to a similar site in the C-terminal domain of p190RhoGEF, and their bindings to p190RhoGEF are readily cross-competed. Moreover, we identify a novel binding site in BC1 to account for its interaction with p190RhoGEF. The findings suggest a novel role of BC1 in differentiated neurons involving RNA-protein interactions of p190RhoGEF.

Mammalian retroelements

The eukaryotic genome has undergone a series of epidemics of amplification of mobile elements that have resulted in most eukaryotic genomes containing much more of this 'junk' DNA than actual coding DNA. The majority of these elements utilize an RNA intermediate and are termed retroelements. Most of these retroelements appear to amplify in evolutionary waves that insert in the genome and then gradually diverge. In humans, almost half of the genome is recognizably derived from retroelements, with the two elements that are currently actively amplifying, L1 and Alu, making up about 25% of the genome and contributing extensively to disease. The mechanisms of this amplification process are beginning to be understood, although there are still more questions than answers. Insertion of new retroelements may directly damage the genome, and the presence of multiple copies of these elements throughout the genome has longer-term influences on recombination events in the genome and more subtle influences on gene expression.

Active Alu element "A-tails": size does matter

Long and short interspersed elements (LINEs and SINEs) are retroelements that make up almost half of the human genome. L1 and Alu represent the most prolific human LINE and SINE families, respectively. Only a few Alu elements are able to retropose, and the factors determining their retroposition capacity are poorly understood. The data presented in this paper indicate that the length of Alu "A-tails" is one of the principal factors in determining the retropositional capability of an Alu element. The A stretches of the Alu subfamilies analyzed, both old (Alu S and J) and young (Ya5), had a Poisson distribution of A-tail lengths with a mean size of 21 and 26, respectively. In contrast, the A-tails of very recent Alu insertions (disease causing) were all between 40 and 97 bp in length. The L1 elements analyzed displayed a similar tendency, in which the "disease"-associated elements have much longer A-tails (mean of 77) than do the elements even from the young Ta subfamily (mean of 41). Analysis of the draft sequence of the human genome showed that only about 1000 of the over one million Alu elements have tails of 40 or more adenosine residues in length. The presence of these long A stretches shows a strong bias toward the actively amplifying subfamilies, consistent with their playing a major role in the amplification process. Evaluation of the 19 Alu elements retrieved from the draft sequence of the human genome that are identical to the Alu Ya5a2 insert in the NF1 gene showed that only five have tails with 40 or more adenosine residues. Sequence analysis of the loci with the Alu elements containing the longest A-tails (7 of the 19) from the genomes of the NF1 patient and the father revealed that there are at least two loci with A-tails long enough to serve as source elements within our model. Analysis of the A-tail lengths of 12 Ya5a2 elements in diverse human population groups showed substantial variability in both the Alu A-tail length and sequence homogeneity. On the basis of these observations, a model is presented for the role of A-tail length in determining which Alu elements are active.

Shared protein components of SINE RNPs

The heterogeneous, short RNAs produced from the high, copy, short mobile elements (SINEs) interact with proteins to form RNA-protein (RNP) complexes. In particular, the BC1 RNA, which is transcribed to high levels specifically in brain and testis from one locus of the ID SINE family, exists as a discrete RNP complex. We expressed a series of altered BC1, and other SINE-related RNAs, in several cell lines and tested for the mobility of the resulting RNP complexes in a native PAGE assay to determine which portions of these SINE RNAs contribute to protein binding. When different SINE RNAs were substituted for the BC1 ID sequence, the resulting RNPs exhibited the same mobility as BC1. This indicates that the protein(s) binding to the ID portion of BC1 is not sequence specific and may be more dependent upon the secondary structure of the RNA. It also suggests that all SINE RNAs may bind a similar set of cellular proteins. Deletion of the A-rich region of BC1 RNA has a marked effect on the mobility of the RNP. Rodent cell lines exhibit a slightly different mobility for this shifted complex when compared to human cell lines, reflecting evolutionary differences in one or more of the protein components. On the basis of mobility change observed in RNP complexes when the A-rich region is removed, we decided to examine poly(A) binding protein (PABP) as a candidate member of the RNP. An antibody against the C terminus of PABP is able to immunoprecipitate BC1 RNA, confirming PABP's presence in the BC1 RNP. Given the ubiquitous role of poly(A) regions in the retrotransposition process, these data suggest that PABP may contribute to the SINE retrotransposition process.

Transport of Neuronal BC1 RNA in Mauthner Axons

In neurons, localized RNAs have been identified in dendrites and axons; however, RNA transport in axons remains poorly understood. Here we analyzed axonal RNA transport in goldfish Mauthner neurons in vivo. BC1 RNA, a noncoding RNA polymerase III transcript that is targeted to dendrites in neurons of the rodent nervous system, was used as a probe for axonal RNA transport. Somata of Mauthner neurons were microinjected with various RNAs. Full-length BC1 RNA, but not control RNAs of similar length, was targeted to both axons and dendrites of Mauthner neurons. BC1 RNA was transported in the form of a rapidly advancing wave front that progressed along axons, in a microtubule-dependent manner, at a rate of 2 micrometer/sec. Whereas a BC1 5' segment of 65 nucleotides was transported to axons and dendrites in a way indistinguishable from full-length BC1 RNA, a BC1 3' segment of 60 nucleotides did not enter Mauthner cell processes to any significant extent. In the wake of the wave advancing through the axon, BC1 RNA was found localized to discrete, spatially delimited domains at the axonal surface. Such demarcated cortical concentrations of BC1 RNA could not be observed after disruption of F-actin organization in the axon. It is concluded that the specific delivery of BC1 RNA to spatially defined axonal target sites is a two-step process that requires the sequential participation of microtubules for long-range axial transport and of actin filaments for local radial transfer and focal accumulation in cortical domains.

Analysis of the SINE S1 Pol III promoter from Brassica; impact of methylation and influence of external sequences

Transcription is an important control point in the transposable element mobilization process. To better understand the regulation of the plant SINE (Short Interspersed Elements) S1, its promoter sequence was studied using an in vitro pol III transcription system derived from tobacco cells. We show that the internal S1 promoter can be functional although upstream external sequences were found to enhance this basal level of transcription. For one putative 'master' locus (na7), three CAA triplets (in positions -12, -7 and -2) and two overlapping TATA motifs (in positions -54 to -43) were important to stimulate transcription. For this locus, two transcription initiation regions were characterized, one centered on position + 1 (first nucleotide of the S1 element) and one centered on position - 19 independently of the internal motifs. The CAA triplets only influence transcription in + 1 and work in association with the internal motifs. We show that methylation can inhibit transcription at the na7 locus. We also observe that S1 RNA is cleaved in a smaller Poly (A) minus product by a process analogous to the maturation of mammalian SINEs.

Upstream flanking sequences and transcription of SINEs

SINEs, short interspersed repeated DNA elements, undergo amplification through retroposition and subsequent integration into a new location in the genome. Each new SINE insertion will be located in a new chromosomal environment, with different flanking sequences. Modulation of transcription by different flanking sequences may play an important role in determining which SINE elements are preferentially active in a genome. We evaluated the ability of upstream flanking sequences to regulate the transcription of three different SINEs (Alu, B2 and ID) by constructing chimeric constructs with known 5' flanking sequences of RNA polymerase III-transcribed genes. Upstream sequences from the 7SL RNA gene, U6 RNA gene, vault RNA gene, and BC1 gene increase transcription of Alu, B2 and BC1 in transient transfections of NIH3T3, HeLa, Neuro2a and C6 glioma cell lines. The 7SL sequence proved most efficient in increasing SINE transcription. The 7SL upstream fused to the BC1 RNA gene (an ID element) was used to create a transgenic mouse line. In contrast to the tissue-specific endogenous BC1 transcription, BC1 transgene transcripts were detected in all tissues tested. However, expression was much higher in those tissues that express the endogenous gene, demonstrating both transcriptional and post-transcriptional regulation. The BC1 RNA was detected in a similar ribonucleoprotein complex in the different tissues.

RNAs from all categories generate retrosequences that may be exapted as novel genes or regulatory elements

While the significance of middle repetitive elements had been neglected for a long time, there are again tendencies to ascribe most members of a given middle repetitive sequence family a functional role--as if the discussion of SINE (short interspersed repetitive elements) function only can occupy extreme positions. In this article, I argue that differences between the various classes of retrosequences concern mainly their copy numbers. Consequently, the function of SINEs should be viewed as pragmatic such as, for example, mRNA-derived retrosequences, without underestimating the impact of retroposition for generation of novel protein coding genes or parts thereof (exon shuffling by retroposition) and in particular of SINEs (and retroelements) in modulating genes and their expression. Rapid genomic change by accumulating retrosequences may even facilitate speciation [McDonald, J.F., 1995. Transposable elements: possible catalysts of organismic evolution. Trends Ecol. Evol. 10, 123-126.] In addition to providing mobile regulatory elements, small RNA-derived retrosequences including SINEs can, in analogy to mRNA-derived retrosequences, also give rise to novel small RNA genes. Perhaps not representative for all SINE/master gene relationships, we gained significant knowledge by studying the small neuronal non-messenger RNAs, namely BC1 RNA in rodents and BC200 RNA in primates. BC1 is the first identified master gene generating a subclass of ID repetitive elements, and BC200 is the only known Alu element (monomeric) that was exapted as a novel small RNA encoding gene.

A novel type of non-coding RNA expressed in the rat brain

We have characterized a novel type of non-coding RNA which consists of tandem repeats of similar sequences, approximately 0.9 kb in size. This RNA, termed Bsr (brain specific repetitive) RNA, is encoded at a single locus (6 q31-->q32) in the rat genome, where 100 to 150 copies of the 0.9 kb sequences are repeated in tandem. Bsr RNA is preferentially expressed in the rat central nervous system (CNS), especially in phylogenetically old structures, such as the pareo- and archicortex, amygdala, thalamus and hypothalamus. In the developing brains, Bsr RNA is expressed in the subsets of differentiating cells but not in proliferating cells. Despite the finding that Bsr RNA appears to be conserved only among the Rattus species, the specific expression pattern of Bsr RNA suggests that it might have some role in the rat CNS.

Control of genes by mammalian retroposons

Available data on possible genetic impacts of mammalian retroposons are reviewed. Most important is the growing number of established examples showing the involvement of retroposons in modulation of expression of protein-coding genes transcribed by RNA polymerase II (Pol II). Retroposons contain conserved blocks of nucleotide sequence for binding of some important Pol II transcription factors as well as sequences involved in regulation of stability of mRNA. Moreover, these mobile genes provide short regions of sequence homology for illegitimate recombinations, leading to diverse genome rearrangements during evolution. Therefore, mammalian retroposons representing a significant fraction of noncoding DNA cannot be considered at present as junk DNA but as important genetic symbionts driving the evolution of regulatory networks controlling gene expression.

Reverse transcriptase: mediator of genomic plasticity

Reverse transcription has been an important mediator of genomic change. This influence dates back more than three billion years, when the RNA genome was converted into the DNA genome. While the current cellular role(s) of reverse transcriptase are not yet completely understood, it has become clear over the last few years that this enzyme is still responsible for generating significant genomic change and that its activities are one of the driving forces of evolution. Reverse transcriptase generates, for example, extra gene copies (retrogenes), using as a template mature messenger RNAs. Such retrogenes do not always end up as nonfunctional pseudogenes but form, after reinsertion into the genome, new unions with resident promoter elements that may alter the gene's temporal and/or spatial expression levels. More frequently, reverse transcriptase produces copies of nonmessenger RNAs, such as small nuclear or cytoplasmic RNAs. Extremely high copy numbers can be generated by this process. The resulting reinserted DNA copies are therefore referred to as short interspersed repetitive elements (SINEs). SINEs have long been considered selfish DNA, littering the genome via exponential propagation but not contributing to the host's fitness. Many SINEs, however, can give rise to novel genes encoding small RNAs, and are the migrant carriers of numerous control elements and sequence motifs that can equip resident genes with novel regulatory elements [Brosius J. and Gould S.J., Proc Natl Acad Sci USA 89, 10706-10710, 1992]. Retrosequences, such as SINEs and portions of retroelements (e.g., long terminal repeats, LTRs), are capable of donating sequence motifs for nucleosome positioning, DNA methylation, transcriptional enhancers and silencers, poly(A) addition sequences, determinants of RNA stability or transport, splice sites, and even amino acid codons for incorporation into open reading frames as novel protein domains. Retroposition can therefore be considered as a major pacemaker for evolution (including speciation). Retroposons, with their unique properties and actions, form the molecular basis of important evolutionary concepts, such as exaptation [Gould S.J. and Vrba E., Paleobiology 8, 4-15, 1982] and punctuated equilibrium [Elredge N. and Gould S.J. in Schopf T.J.M. (ed). Models in Paleobiology. Freeman, Cooper, San Francisco, 1972, pp. 82-115].