Department of Biochemistry and Molecular Biology, Joef Stefan Institute, Ljubljana, Slovenia;
Department of Chemistry and Biochemistry, Faculty of Chemistry and Chemical Technology, University of Ljubljana, Ljubljana, Slovenia
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Nonlong terminal repeat (LTR) retrotransposons, one of the oldest groups of retrotransposons, have been classified into 12 clades (Malik, Burke, and Eickbush 1999
; Malik and Eickbush 2000
). Vertebrate genomes contain representatives from the L1, CR1, and RTE clades. Representatives of the remaining nine clades are distributed in invertebrates, fungi, and plants (Malik, Burke, and Eickbush 1999
).
The evolutionary dynamics of TEs in their hosts differ, but in most cases they involve vertical transmission and occasionally also horizontal transmission. The life cycles of different groups of TEs consist of vertical inactivation of active elements, stochastic loss of TEs from populations, and eventual reintroduction by horizontal transfer (Lohe et al. 1995
; Hartl, Lohe, and Lozovskaya 1997
).
The phenomenon of horizontal transfer (HT) has been observed and well documented in most TE classes (Kidwell 1993
; Capy, Anxolabehere, and Langin 1994
; Hartl, Lohe, and Lozovskaya 1997
). The most-studied and best-understood examples are those from DNA transposons represented by mariner (Robertson 1993
; Gueiros-Filho and Beverley 1997
; Hartl, Lohe, and Lozovskaya 1997
) and the P element (Houck et al. 1991
; Kidwell 1993
; Clark and Kidwell 1997
; Silva and Kidwell 2000
). In the most primitive RT-containing elements, the mobile introns, a very large number of HTs have been discovered in higher plant groups (Cho et al. 1998
). In the LTR retrotransposons, several examples of HT have been documented between closely related species (mostly among Drosophila), such as copia (Jordan and McDonald 1998
; Jordan, Matyunina, and McDonald 1999
) and gypsy (Terzian et al. 2000)
, and between more distantly related species, represented by SURL elements from sea urchins (Gonzalez and Lessios 1999
). In contrast to the general acceptance of HT of DNA transposons, mobile introns, and LTR-retrotransposons, several proposed cases of HT in non-LTR retrotransposons have met with strong criticism (Malik, Burke, and Eickbush 1999
). This issue was further investigated, and new evidence for HT of non-LTR retrotransposons is presented.
Members of the RTE clade are the shortest of the non-LTR retrotransposons and have a highly conserved structural organization (Malik and Eickbush 1998
; Malik, Burke, and Eickbush 1999
). They have been studied in nematodes and insects (Malik and Eickbush 1998
), mammals (Szemraj et al. 1995
; Okada and Hamada 1997
; Malik and Eickbush 1998
), trematodes (Drew et al. 1999
), and teleost fishes (Volff et al. 1999
). On the basis of discontinuous distribution, extreme nucleotide sequence conservation, genetic distances, and evolutionary relationships, the HT of Bov-B long interspersed nuclear elements (LINEs) 4050 MYA from the ancestral snake lineage (Boidae) to the ancestor of ruminants has been demonstrated (Kordi
and Guben
ek 1998, 1999a, 1999b
). At present, only one full-length bovine Bov-B LINE element is known (Szemraj et al. 1995
; Okada and Hamada 1997
; Malik and Eickbush 1998
), and no reptilian species have so far been investigated. Here we present the first full-length reptilian Bov-B LINE element and examine the evolutionary dynamics of Bov-B LINEs in vertebrates and the distribution and evolutionary relationships of the RTE clade in eukaryotes.
![]() |
Materials and Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
PCR Amplification of the 1.8-kb Fragment of Bov-B LINEs
All experiments were performed in parallel with negative and positive controls. Different rooms, reagents, equipment, and positive displacement pipettes were used, according to the general precautions for PCR performance. The sense (ME1: 5'-CACRGTRATYCAAGYCTAYRCYCCAAC-3') and antisense (ME2: 5'-CWGCAWATCTGAGGTTMKTKAKATTTCT-3') degenerate oligonucleotide primers were based on conserved regions in the endonuclease (TVIQVYAPTME1) and RT (RNINNLRYAME2) domains of Bos taurus Bov-B LINE and Caenorhabditis elegans RTE-1 element (Malik and Eickbush 1998
). PCR amplification of the 1.8-kb fragment of Bov-B LINEs was performed in a 100-µl volume with 1 µg of genomic DNA, each dNTP at 200 µM, 50 pmol of each primer, 2.5 U of AmpliTaq DNA polymerase, and buffer provided by the supplier of the enzyme (Perkin-Elmer). After an initial denaturation step of 5 min at 95°C, the PCR reactions were subjected to 30 cycles of amplification consisting of 2 min denaturation at 95°C, 2 min annealing at 60°C, and 2 min extension at 72°C, with a 5-min final extension at 72°C. The resulting PCR products were directly ligated into a pGEM vector using a pGEM-T-easy cloning kit (Promega) for sequence determination. The inserts were sequenced on both strands with an ABI fluorescent sequencing kit on an ABI 310 sequencer.
Isolation of a Full-Length Bov-B LINE Element from a Vipera ammodytes Genomic Library
The GEM-12 (Kordi
and Guben
ek 1996
) genomic library was screened with the 35S-labeled 1.8-kb fragment of the Bov-B LINE using the plaque hybridization method. Hybridization was carried out at 42°C for 20 h in a mixture of 6 x standard saline citrate (SSC), 5 x Denhardt's solution, 0.5% SDS, and denatured herring sperm DNA at 100 µg/ml in 50% formamide. The filters were washed successively with 6 x SSC and 2 x SSC at 35°C for 20 min each. The positive clones were rescreened by the same procedure. Phage DNA was prepared from plate lysates (Sambrook, Fritsch, and Maniatis 1989
, p. 2.64) and digested with BamHI, EcoRI, SacI, and XhoI restriction enzymes. The resulting fragments were separated by gel electrophoresis on 0.7% agarose, transferred to Hybond-N membranes (Amersham), and hybridized at 42°C with the same probe as described above. Positive genomic fragments were subcloned into pUC19.
DNA- and Protein-Based Sequence Analyses
Computer-based nucleotide and protein searches of the GenBank databases with reptilian Bov-B LINEs and several representatives from the RTE clade were performed with the different BLAST (Altschul et al. 1990
) search programs of the NCBI. Specialized EBI databases were searched with different Fasta programs (Pearson 1990
). The following sequence databases were searched at NCBI (www.ncbi.nlm.nih.gov): nonredundant (NR), dbEST (expressed sequence tags database), dbSTS (sequence tagged sites database), dbGSS (genome survey sequences database), and HTGS (unfinished high-throughput genomic sequencesphases 0, 1, and 2); while at EBI (www.ebi.ac.uk), we searched parasite genomes (www.ebi.ac.uk/parasite-genome.html) and eukaryotic genomes (www.ebi.ac.uk/fasta33/genomes.html). Evolutionary rates were estimated by standard methods (Nei and Kumar 2000
, p. 20). Poisson correction distances (d) were estimated by the equation d = -ln(1 - p), where p represents the proportion of different amino acids. The rate of amino acid substitution (r) was estimated by the standard equation r = d/2T, where T is the divergence time of the last common ancestor (LCA) of the species compared. Amino acid distances used in divergence-versus-age analysis were calculated from sequences of the complete RT domain using the MEGA2 program (Kumar et al. 2000)
.
Phylogenetic Analyses
The amino acid sequences of Bov-B LINE and other RTE representatives were aligned using CLUSTAL W (Thompson, Higgins, and Gibson 1994
) with some manual refinements. Phylogenetic trees were inferred using the neighbor-joining (NJ) method (Saitou and Nei 1987
) as implemented in TREECON and MEGA2 (Van de Peer and De Wachter 1994
; Kumar et al. 2000
). The significance of the various phylogenetic lineages was assessed by bootstrap analysis. Phylogenetic analyses were performed on separate regions of the RTE elements, such as apurinic/apyrimidinic endonuclease (AP-EN), the palm/fingers subdomain of RT, and the thumb subdomain of RT. Several potential outgroups were tested, and group II intron RT (from Neurospora, accession number S07649), Schistosoma SR2 element AP-EN, and the thumb subdomain of RT (Drew et al. 1999
) were selected, since the phylogenetic trees were in accord with the species phylogeny and taxonomy. Gaps in aligned sequences were removed for the purpose of analysis.
![]() |
Results and Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
On the basis of the nucleotide and amino acid sequences of the V. ammodytes Bov-B LINE, it is possible to reconstruct the amino acid sequence of the single previously available full-length ruminant Bov-B LINE element (Szemraj et al. 1995
; Okada and Hamada 1997
). This element was previously reconstructed on the basis of the C. elegans RTE-1 element (Malik and Eickbush 1998
) but contained some errors, such as the part of AP-EN region from C. elegans RTE-2 element. Nematode RTE-1 and bovine Bov-B LINE show only 25%30% identity at the amino acid level, in contrast to the 75% identity between the full-length Vam and B. taurus Bov-B LINE elements.
Elements from the RTE clade are the shortest among the non-LTR retrotransposons and have a highly conserved structural organization (Malik and Eickbush 1998
; Malik, Burke, and Eickbush 1999
). Several representatives from the RTE clade have been analyzed, with some of them being full-length, such as the C. elegans RTE-1 element (Malik and Eickbush 1998
), the B. taurus Bov-B LINE (BDDF) element (Szemraj et al. 1995
; Okada and Hamada 1997
; Malik and Eickbush 1998
), the Schistosoma mansoni SR2 element (Drew et al. 1999
), the Rex3 element from teleost fishes (Volff et al. 1999
), the JAM1 element from Aedes (Malik and Eickbush 1998
), and now the full-length reptilian V. ammodytes Bov-B LINE (fig. 2
). We discovered several novel nearly full length elements in GenBank databases, such as Bombyx mori Bov-B/RTE, Strongylocentrotus Bov-B/RTE, Oryzias RTE, and plant RTEs (fig. 2 ), and these will be described below in more detail.
|
We cloned 1.8-kb fragments and sequenced several (two to four) independent clones from each of 10 selected Squamata species. Intraspecies variability among several clones was very low and the same has been found for interspecies variability. Genetic distances were even smaller than those previously obtained for the extreme C-terminal part of the Bov-B LINEs (Kordi and Guben
ek 1998
). In some sequences, we found a few stop codons and short deletions; these result in frameshift mutations or in frame stop mutations.
It has been reported that marsupials, but not monotremes, also contain Bov-B LINEs (Gilbert and Labuda 1999, 2000
), but no experimental data were provided. We amplified the 1.8-kb fragment and a shorter 0.5-kb (thumb subdomain) fragment from a marsupial (Macropus) and from the ruminants tested. Since the monotremes and marsupials form a sister group, Marsupionta (Janke, Xu, and Arnason 1997
), the acquisition of Bov-B LINEs in marsupials or their loss in monotremes needs additional study with sequence and phylogenetic analyses in order to clarify their origin in mammals. Bov-B LINEs may have originally been present in an ancestor of the mammals, or, alternatively, they could have been introduced on two separate occasions by HT into marsupial and ruminant genomes.
It is important to note that all other mammalian species testedrepresentatives of different mammalian orders: Cetartiodactyla (Tragulus [PCR was negative, although the AP-EN domain is already known], Hippopotamus, Tursiops, Sus), Carnivora (Canis), Rodentia (Mus), and Primates (Homo)were negative by our PCR analyses. Additional support for the restricted and discontinuous presence of Bov-B LINEs in mammals provided an examination of mammalian sequences in GenBank databases. The currently known distribution of Bov-B LINEs in mammals is highly discontinuous; they are absent in most of the mammalian orders. Hybridization and PCR amplification data (Buntjer, Hoff, and Lenstra 1997
; Kordi
and Guben
ek 1997, 1998
; Shimamura et al. 1999
; Gilbert and Labuda 2000
) also do not support their widespread distribution in mammals, or even panvertebrate distribution, as suggested recently (Malik and Eickbush 1998
; Gilbert and Labuda 1999
). In order to determine their evolutionary origins in mammalian genomes, a much wider study in the context of the latest mammalian phylogeny (Madsen et al. 2001
; Murphy et al. 2001)
is needed. Marsupial Bov-B LINEs are not inconsistent with the hypothesis of their horizontal origin in Ruminantia.
Our preliminary sequence data of the marsupial Bov-B LINE encoding the thumb domain of RT (fig. 3
) shows up to 78% identity at the DNA level and 65% at amino acid level with the ancestral snake lineage, a level of conservation similar to that previously observed for ruminants (Kordi and Guben
ek 1998
). This may indicate an additional horizontal transfer of Bov-B LINEs into the ancestor of marsupials at least 120 MYA. However, more extensive analyses of the monotremes, marsupials, and placental mammals will be required to conclusively answer the question of the origin of Bov-B LINEs and their current distribution in mammals.
|
|
|
Wide Distribution of the RTE Clade
The distribution of the RTE clade of non-LTR retrotransposons was reanalyzed and found to be much wider than previously thought (table 2
). By searching specific databases (NR, GSS, HTGS, EST, STS, parasite genomes, and completed eukaryote genomes) and by searching within specific taxonomic groups using the full-length sequence of the Vam Bov-B LINE and other novel representatives of the RTE clade, we discovered several new representatives (table 3
).
|
|
We also detected short fragments of RT in brown algae (Laminaria); the size of the translated sequence was relatively short, encoding only 130 amino acids. Shorter fragments with significant similarity were also detected in annelids (Lumbricus), mollusks (Helix pomatia), and crustaceans (Artemia, Pennaeus). Additionally, there was a relatively high level of conservation (40%) between Laminaria and sea urchin RTEs, in contrast to the lower level between echinoderm RTEs and teleost Rex3 elements, which showed only 30% amino acid identity.
In previous analyses of the RTE clade (Malik and Eickbush 1998
; Drew et al. 1999
; Volff et al. 1999
), representatives from nematodes, arthropods, mollusks, trematodes, and vertebrates have been described. These data were extended in the present study to cover a much wider taxonomic distribution. The RTE clade is seen to originate very early in the evolution of eukaryotes, since they are present in stramenopiles, in plants, and apparently in most metazoan phyla. Most of the novel representatives recognized by TBLASTN searching were their C-terminal domains only. Available sequence data of the RTE clade show that the level of conservation on a global scale is quite low, as was found for other clades of non-LTR retrotransposons (Malik and Eickbush 1998
; Malik, Burke, and Eickbush 1999
). The availability of these new basal eukaryotic RTE representatives will enable the design of novel oligonucleotide primers and consequently much denser sampling that will clarify their long-term evolution in eukaryotes. The RTE and L1 clades are the most widespread of the known eukaryotic non-LTR retrotransposons.
Evolutionary Relationships Among Bov-B LINEs and Other Representatives of the RTE Clade
To determine the evolutionary relationships of the Bov-B LINE elements and other representatives of the RTE clade, the amino acid sequences of AP-EN, RT, and the C-terminal part of RT (thumb subdomain) were subjected to phylogenetic analysis using the NJ algorithm (Saitou and Nei 1987
). Because the level of sequence divergence within each species was insignificant compared with that between species, either a consensus or a single Bov-B LINE sequence was used to represent the elements of each species. Most Squamata Bov-B LINE elements are seen to be grouped in phylogenetic trees according to their species of origin (fig. 5 ).
|
NJ phylogenetic trees, based on amino acid sequence alignments of AP-EN, the palm/fingers subdomain of RT, and the thumb subdomain of RT domains, produce the same general pattern. In all of these trees, the Ruminantia Bov-B LINEs group together with the ancestral snake lineage (Boidae), as if they were the closest relatives (fig. 5
). The same feature has previously been recognized at the DNA level, from the thumb subdomain of RT (Kordi and Guben
ek 1998;, 1999a, 1999b
). The unusual position of mammalian (ruminant) Bov-B LINEs inside the reptilian cluster is thus confirmed and is a clear indication of HT. These evolutionary relationships show that potential HT of Bov-B LINEs or any other TE can easily be recognized from phylogenetic trees based on either DNA or amino acid sequences of the short parts of RT or AP-EN, or even from analyses of the full-length elements.
In the phylogenetic analyses of the RTE clade, we included mostly elements without a large number of stop codons or difficulties in translation. Frameshift mutations and stop codons were frequently observed during TBLASTN searching, but in most species and RTE lineages the elements are relatively well conserved and can be easily translated. As an outgroup, we used the very ancient group II introns (Neurospora) for the palm/fingers subdomain of RT, while for AP-EN and the thumb subdomain of RT this ancient sequence was too divergent and not useful for alignment or phylogenetic analyses. Schistosoma SR2 elements AP-EN and the thumb subdomain were therefore used as outgroups in such trees. Phylogenetic reconstruction of evolutionary relationships among other RTE representatives showed several well-separated groups, such as B. mori Bov-B/RTE, sea urchin (Strongylocentrotus) Bov-B/RTE, nematode RTE, Rex3 elements from teleost fishes, and plant RTEs, together with Oryzias (medaka fish) RTE. The positions of trematoda SR2 and mosquito JAM1 are not well resolved, since they show only about 30% amino acid identity with all other representatives of the RTE clade. These phylogenetic analyses of the RTE clade used larger numbers of RTE representatives, from plants to vertebrates, than previous phylogenetic analyses (Malik and Eickbush 1998
; Drew et al. 1999
; Volff et al. 1999
).
Divergence-Versus-Age Analysis of the RTE Clade
Amino acid distances between the RT domains of the RTE clade representatives are plotted against estimates of host divergence time in figure 6
. The RTE lineages show a regular pattern of increased divergence with time that approaches saturation at the time of the split of Tetrapoda and Teleostei at 400 MYA (Kumar and Hedges 1998
). Comparisons of species separated by more than 600 Myr have little resolution, as has been observed for R1 and R2 elements (Malik, Burke, and Eickbush 1999
). Comparisons of amino acid distances versus host divergence time were made within and between plant and metazoan RTE lineages. One of the comparisons between monocot and dicot plants fell above the curve and probably represented a paralogous comparison. Three examples in which points fell markedly below the RTE curve were Squamata versus Ruminantia Bov-B LINEs (at 310 Myr), B. mori Bov-B/RTE versus vertebrate Bov-B LINEs (at 670 Myr), and plant RTE versus Oryzias (medaka fish) RTE (at 1,200 Myr). It has been proposed (Malik, Burke, and Eickbush 1999
) that such points can be interpreted as an indication of HT occurrence. This finding has important implications for the existence of HT in non-LTR retrotransposons.
|
The level of amino acid divergence between vertebrate Bov-B LINEs and invertebrate B. mori Bov-B-like elements seems relatively low with regard to their taxonomic position and their LCA with vertebrates (around 670 MYA). How is it possible that invertebrate Bov-B-like elements are much more similar to vertebrate Bov-B LINEs than they are to the teleost Rex3 elements? The most parsimonious explanation, as supported by the anomalous position in figure 6 , would be separate horizontal transfers of Bov-B LINEs of invertebrate origin into the ancestor of Squamata, the ancestor of Ruminantia and the ancestor of marsupials.
Regarding the amino acid divergence observed between medaka fish (Oryzias) and the plant RTE lineage with the age of the LCA between the plant and metazoan kingdoms being >1,200 Myr (Wang, Kumar, and Hedges 1999
), their evolutionary relationships are explainable only by HT. Medaka fish, however, could not be the donor of the RTE element in the potential HT for the plants, since plant RTE elements are distributed throughout the Embryophyta.
Evolutionary Rates in the RTE Clade of Non-LTR Retrotransposons
Evolutionary rates for reptilian and mammalian Bov-B LINEs and the other representatives of the RTE clade have been estimated (table 4
), leading to several interesting observations. First, the evolutionary rates in Bov-B LINEs are very low, as in Rex3 elements from bony fishes, another representative of the RTE clade in vertebrates. Second, evolutionary rates of representatives of the RTE clade in vertebrates are slower than those in invertebrates. Third, evolutionary rates between endothermic (Ruminantia) and ectothermic (Squamata, Teleostei) vertebrates are different; slower rates of evolution were found in ectotherms. Fourth, long-term evolutionary rates in the RTE clade correlate with divergence times, with a significantly lower evolutionary rate suggesting the occurrence of HT.
|
Multiple Lineages of RTE Elements Are Present in Some Genomes
Previous studies have confirmed that multiple RTE lineages coexist in nematode genomes, such as the RTE-1 and RTE-2 lineages in C. elegans (Malik and Eickbush 1998
). In a search of the sea urchin genome database, one group of elements was found to be more closely related to the C. elegans RTE-1 elements, while another lineage was much more similar to the Bov-B LINEs (table 3 ). Another example of multiple lineages is demonstrated in the genome of medaka fish (Oryzias), where the Rex3 (Volff et al. 1999
) and Oryzias RTE elements coexist, showing only 30%35% identity, but where the Oryzias RTE element surprisingly shows the strongest similarity to the plant RTE lineage.
Evolutionary History of Bov-B LINEs in Deuterostomia
By searching the GSS database, we found that Bov-B-like elements are present in the most ancestral lineage of Deuterostomia, Echinodermata, represented by the sea urchin Strongylocentrotus. These elements exhibit 40%50% amino acid identity with the Squamata and Ruminantia Bov-B LINEs. However, they have not been retained in Urochordata (e.g., Ciona [GSS database] and Halocynthia [EST database]), Cephalochordata (Branchiostoma [HTGS database]), hagfishes (Myxine), lampreys (Petromyzontiformes), bony fishes (Teleosteithey contain another lineage, the Rex3 elements), cartilaginous fishes (Chondrichthyes), or Amphibia (PCR data and Xenopus [EST database]). They appear only in the most ancestral reptiles, in Squamata (Janke et al. 2001)
, but are absent from the turtles and Archosauria (crocodiles and birds). In mammals, they appear only in marsupials and ruminants. Why, then, are Bov-B LINEs not present in all vertebrates?
Discontinuous distribution of Bov-B LINEs in Deuterostomia can be explained by two alternative hypotheses. The first proposes that Bov-B LINEs were present in the LCA of Deuterostomia and have been transmitted vertically. A complex evolutionary scenario involving an unreasonably high number of evolutionary losses is necessary to explain the very sporadic distribution of Bov-B LINEs. Even the Bov-B LINE lineage is not retained in any of the vertebrate lineages leading to reptiles and mammals. Once lost from a lineage, it would be absent from all descendant taxa. The second hypothesis involves the possibility of HT. Three events of HT in the ancestors of Squamata, Ruminantia, and Marsupialia can quite simply explain their current distribution, providing the most parsimonious explanation for the discontinuous distribution.
Horizontal Transfer Versus Alternative Hypotheses
Evolutionary studies of TEs often yield element-generated phylogenies that are incongruent with the host species phylogenies, such as those we report, but these inconsistencies are not always indicative of HT (Capy, Anxolabehere, and Langin 1994
). Many factors can obscure phylogenetic reconstruction of multicopy TEs. Comparison of paralogous copies of elements and varying rates of sequence evolution of TE copies within and between species are factors which can yield incongruent phylogenies even under conditions of strict vertical transmission. Ancestral polymorphism, coupled with independent assortment of copies into the descendant species, inequality of substitution rates in TE sequences in different species, and the stochastic loss of TEs from a few taxa, could obscure phylogenetic reconstruction of TEs and lead to incongruence in phylogenetic reconstruction. Alternative explanations are often hard to dismiss conclusively, especially when TEs from closely related taxa are compared.
HT has been traditionally inferred when a high degree of similarity exists between TEs, coupled with a long divergence time of their respective host species. Incongruence between TE and host phylogenies, or the absence of the TE in question from taxa closely related to that into which the TE was supposedly transferred horizontally, can help corroborate that inference.
By comparing the divergence of TE amino acid sequences with those observed for host genes evolving under similar or stronger selective constraints, an HT event can be inferred whenever the divergence among TE sequences is significantly lower than that observed for the host proteins. We found that the Bov-B LINEs compared between Ruminantia and Squamata (LCA 310 MYA) show much lower divergence than the typical cellular and mitochondrial proteins from the same taxa, strongly indicating HT. No other vertebrate non-LTR retrotransposons, evolving by strict vertical transmission, show such an anomaly.
Selection and recent HT have very different consequences both for the level of congruence expected between TE and host phylogenies and for the proportionality of TE divergence relative to host gene divergence when sister groups of different ages are compared (Silva and Kidwell 2000)
. If the low TE divergence resulted from selection and the elements were transmitted vertically, TE and host phylogenies would be expected to be congruent. TE and host gene divergence should be correlated (Silva and Kidwell 2000)
. We observed this pattern in the Squamata Bov-B LINEs, as shown by very low sequence divergence and an element phylogeny that is in accord with the host phylogeny.
If recent HT caused the low TE divergence, then the TE phylogeny is not necessarily expected to be congruent with that of the host species. Also, TE and host gene divergence should not be correlated (Silva and Kidwell 2000)
. The erroneous phylogenetic position of Ruminantia and the comparison with typical cellular and mitochondrial proteins shows this pattern and strongly supports the occurrence of HT.
Evolutionary analysis of non-LTR retrotransposons (Malik, Burke, and Eickbush 1999
) suggests that vertical transmission is the most common mode of inheritance of these elements, but, according to the present study, one should not exclude the possibility of occasional HT events. Only five HT events, between Squamata/Ruminantia, Marsupialia/Ruminantia, invertebrates/vertebrates, Oryzias fishes/plants, and Laminaria/sea urchin, are needed to explain all distributional and phylogenetic discrepancies in the RTE clade. It is obvious that Bov-B LINEs and some other RTE clade representatives are evolving at a very slow rate, at least in reptilian genomes (Bov-B LINEs) and teleost fishes (Rex3). Taken together, the data from the present study suggest several explanations which are not necessarily mutually exclusive to explain the apparent discrepancies between RTE clade and species evolution.
Mechanism of HT
One of the problems in dealing with the HT of TEs is that very little is known about the actual mechanisms through which this transmission could occur. Viruses, bacteria, fungi, parasites, or symbionts could act as potential vectors. All that the vector would need is a broad host range and the ability to somehow gain access to the germ line (Kidwell 1993
). The potential HT events in the RTE clade have occurred between different classes, phyla, or kingdoms. For all of these cases, a possible ecological connection exists between them (parasitizing or feeding). The ability of any element from the RTE clade to function in such diverse hosts indicates that they are not dependent on specific host factors for their activity.
Even though we cannot directly show experimentally how the HT events may have occurred, we propose that the evidence for HT in our data indicates that the non-LTR retrotransposons (at least among several RTE lineages) do occasionally transfer horizontally. Previous conclusions to the contrary (Malik, Burke, and Eickbush 1999
) are based on the relative rarity of HT events, which makes it necessary that many species be sampled in order to detect them.
Conclusions
The results presented here show that reptiles contain full-length Bov-B LINE elements, which are highly conserved in Squamata. A very high level of sequence conservation exists between the ancestral snake lineage and Ruminantia. Bov-B LINEs are considerably more conserved than the typical cellular and mitochondrial proteins. The RTE clade, together with the L1 clade, has been shown to be one of the most widespread non-LTR retrotransposons, originating very early in the evolution of eukaryotes. Evolutionary rates in some RTE lineages are very low. Using different independent methods, including the divergence-versus-age analysis, several examples of HT in the RTE clade have been recognized, with important implications for the existence of HT in non-LTR retrotransposons.
![]() |
Supplementary Material |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Acknowledgements |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Footnotes |
---|
1 Abbreviations: AP-EN, apurinic/apyrimidinic endonuclease; HT, horizontal transfer; LCA, last common ancestor; LINE, long interspersed nuclear element; LTR, long terminal repeat; NJ, neighbor joining; ORF, open reading frame; RT, reverse transcriptase; TE, transposable element.
2 Keywords: non-LTR retrotransposon
RTE clade
Bov-B LINE
horizontal transfer
evolutionary dynamics
3 Address for correspondence and reprints: Duan Kordi
, Department of Biochemistry and Molecular Biology, Jo
ef Stefan Institute, Jamova 39, 1001 Ljubljana, Slovenia. dusan.kordis{at}ijs.si
.
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Altschul S. F., W. Gish, W. Miller, E. W. Myers, D. J. Lipman, 1990 Basic local alignment search tool J. Mol. Biol 215:403-410[ISI][Medline]
Bennetzen J. L., 2000 Transposable element contributions to plant gene and genome evolution Plant Mol. Biol 42:251-269[ISI][Medline]
Buntjer J. B., I. A. Hoff, J. A. Lenstra, 1997 Artiodactyl interspersed DNA repeats in cetacean genomes J. Mol. Evol 45:66-69[ISI][Medline]
Capy P., D. Anxolabehere, T. Langin, 1994 The strange phylogenies of transposable elements: are horizontal transfers the only explanation? Trends Genet 10:7-12[ISI][Medline]
Cho Y., Y. L. Qiu, P. Kuhlman, J. D. Palmer, 1998 Explosive invasion of plant mitochondria by a group I intron Proc. Natl. Acad. Sci. USA 95:14244-14249
Clark J. B., M. G. Kidwell, 1997 A phylogenetic perspective on P transposable element evolution in Drosophila Proc. Natl. Acad. Sci. USA 94:11428-11433
Drew A. C., D. J. Minchella, L. T. King, D. Rollinson, P. J. Brindley, 1999 SR2 elements, non-long terminal repeat retrotransposons of the RTE-1 lineage from the human blood fluke Schistosoma mansoni Mol. Biol. Evol 16:1256-1269[Abstract]
Fedoroff N., 2000 Transposons and genome evolution in plants Proc. Natl. Acad. Sci. USA 97:7002-7007
Feng D., G. Cho, R. F. Doolittle, 1997 Determining divergence times with a protein clock: update and reevaluation Proc. Natl. Acad. Sci. USA 94:13028-13033
Gilbert N., D. Labuda, 1999 CORE-SINEs: eukaryotic short interspersed retroposing elements with common sequence motifs Proc. Natl. Acad. Sci. USA 96:2869-2874
. 2000 Evolutionary inventions and continuity of CORE-SINEs in mammals J. Mol. Biol 298:365-377[ISI][Medline]
Gonzalez P., H. A. Lessios, 1999 Evolution of sea urchin retroviral-like (SURL) elements: evidence from 40 echinoid species Mol. Biol. Evol 16:938-952[Abstract]
Gueiros-Filho F. J., S. M. Beverley, 1997 Trans-kingdom transposition of the Drosophila element mariner within the protozoan Leishmania Science 276:1716-1719
Hartl D. L., A. R. Lohe, E. R. Lozovskaya, 1997 Modern thoughts on an ancyent marinere: function, evolution, regulation Annu. Rev. Genet 31:337-358[ISI][Medline]
Hedges S. B., L. L. Poling, 1999 A molecular phylogeny of reptiles Science 283:998-1001
Houck M. A., J. B. Clark, K. R. Peterson, M. G. Kidwell, 1991 Possible horizontal transfer of Drosophila genes by the mite Proctolaelaps regalis Science 253:1125-1128[ISI][Medline]
Janke A., D. Erpenbeck, M. Nilsson, U. Arnason, 2001 The mitochondrial genomes of the iguana (Iguana iguana) and the caiman (Caiman crocodylus): implications for amniote phylogeny Proc. R. Soc. Lond. B Biol. Sci 268:623-631[ISI][Medline]
Janke A., X. Xu, U. Arnason, 1997 The complete mitochondrial genome of the wallaroo (Macropus robustus) and the phylogenetic relationships among Monotremata, Marsupialia, and Eutheria Proc. Natl. Acad. Sci. USA 94:1276-1281
Jordan I. K., J. F. McDonald, 1998 Evolution of the copia retrotransposon in the Drosophila melanogaster species subgroup Mol. Biol. Evol 15:1160-1171[Abstract]
Jordan I. K., L. V. Matyunina, J. F. McDonald, 1999 Evidence for the recent horizontal transfer of long terminal repeat retrotransposon Proc. Natl. Acad. Sci. USA 96:12621-12625
Kidwell M. G., 1993 Lateral transfer in natural populations of eukaryotes Annu. Rev. Genet 27:235-256[ISI][Medline]
Kidwell M. G., D. Lisch, 1997 Transposable elements as sources of variation in animals and plants Proc. Natl. Acad. Sci. USA 94:7704-7711
. 2000 Transposable elements and host genome evolution Trends Ecol. Evol 15:95-99[ISI][Medline]
Kordi D., F. Gubenek, 1996 Ammodytoxin C gene helps to elucidate the irregular structure of Crotalinae group II phospholipase A2 genes Eur. J. Biochem 240:83-90[Abstract]
. 1997 Bov-B long interspersed repeated DNA (LINE) sequences are present in Vipera ammodytes phospholipase A2 genes and in genomes of Viperidae snakes Eur. J. Biochem 246:772-779[Abstract]
. 1998 Unusual horizontal transfer of a long interspersed nuclear element between distant vertebrate classes Proc. Natl. Acad. Sci. USA 95:10704-10709
. 1999a. Horizontal transfer of non-LTR retrotransposons in vertebrates Genetica 107:121-128[ISI][Medline]
. 1999b. Molecular evolution of Bov-B LINEs in vertebrates Gene 238:171-178[ISI][Medline]
Kumar A., J. L. Bennetzen, 1999 Plant retrotransposons Annu. Rev. Genet 33:479-532[ISI][Medline]
Kumar S., S. B. Hedges, 1998 A molecular timescale for vertebrate evolution Nature 392:917-920[ISI][Medline]
Kumar S., K. Tamura, I. Jakobsen, M. Nei, 2000 MEGA 2: molecular evolutionary genetics analysis program Version 2.0. Pennsylvania State University, University Park, Penn.; and Arizona State University, Tempe, Ariz
Lohe A. R., E. N. Moriyama, D. A. Lidholm, D. L. Hartl, 1995 Horizontal transmission, vertical inactivation, and stochastic loss of mariner-like transposable elements Mol. Biol. Evol 12:62-72[Abstract]
Madsen O., M. Scally, C. J. Douady, D. J. Kao, R. W. Debry, R. Adkins, H. M. Amrine, M. J. Stanhope, W. W. De Jong, M. S. Springer, 2001 Parallel adaptive radiations in two major clades of placental mammals Nature 409:610-614[ISI][Medline]
Malik H. S., W. D. Burke, T. H. Eickbush, 1999 The age and evolution of non-LTR retrotransposable elements Mol. Biol. Evol 16:793-805[Abstract]
Malik H. S., T. H. Eickbush, 1998 The RTE class of non-LTR retrotransposons is widely distributed in animals and is the origin of many SINEs Mol. Biol. Evol 15:1123-1134[Abstract]
. 2000 NeSL-1, an ancient lineage of site-specific non-LTR retrotransposons from Caenorhabditis elegans Genetics 154:193-203
Murphy W. J., E. Eizirik, W. E. Johnson, Y. P. Zhang, O. A. Ryder, S. J. O'Brien, 2001 Molecular phylogenetics and the origin of placental mammals Nature 409:614-618[ISI][Medline]
Nei M., S. Kumar, 2000 Molecular evolution and phylogenetics Oxford University Press, New York
Okada N., M. Hamada, 1997 The 3' ends of tRNA-derived SINEs originated from the 3' ends of LINEs: a new example from the bovine genome J. Mol. Evol 44: (Suppl. 1) S52-S56[ISI][Medline]
Pearson W. R., 1990 Rapid and sensitive sequence comparison with FASTP and FASTA Meth. Enzymol 183:63-98[ISI][Medline]
Robertson H. M., 1993 The mariner transposable element is widespread in insects Nature 362:241-245[ISI][Medline]
Saitou N., M. Nei, 1987 The neighbor-joining method: a new method for reconstructing phylogenetic trees Mol. Biol. Evol 4:406-425[Abstract]
Sambrook J., E. F. Fritsch, T. Maniatis, 1989 Molecular cloning: a laboratory manual 2nd edition. Cold Spring Harbor Laboratory, Cold Spring Harbor, N.Y
Shapiro J. A., 1999 Transposable elements as the key to a 21st century view of evolution Genetica 107:171-179[ISI][Medline]
Shimamura M., H. Abe, M. Nikaido, K. Ohshima, N. Okada, 1999 Genealogy of families of SINEs in cetaceans and artiodactyls: the presence of a huge superfamily of tRNAGlu-derived families of SINEs Mol. Biol. Evol 16:1046-1060[Abstract]
Silva J. C., M. G. Kidwell, 2000 Horizontal transfer and selection in the evolution of P elements Mol. Biol. Evol 17:1542-1557
Szemraj J., G. Plucienniczak, J. Jaworski, A. Plucienniczak, 1995 Bovine Alu-like sequences mediate transposition of a new site-specific retroelement Gene 152:261-264[ISI][Medline]
Terzian C., C. Ferraz, J. Demaille, A. Bucheton, 2000 Evolution of the gypsy endogenous retrovirus in the Drosophila melanogaster subgroup Mol. Biol. Evol 17:908-914
Thompson J. D., D. G. Higgins, T. J. Gibson, 1994 CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice Nucleic Acids Res 22:4673-4680[Abstract]
Van de Peer Y., R. De Wachter, 1994 TREECON for Windows: a software package for the construction and drawing of evolutionary trees for the Microsoft Windows environment Comput. Appl. Biosci 10:569-570[Medline]
Volff J. N., C. Korting, K. Sweeney, M. Schartl, 1999 The non-LTR retrotransposon Rex3 from the fish Xiphophorus is widespread among teleosts Mol. Biol. Evol 16:1427-1438[Abstract]
Wang D. Y., S. Kumar, S. B. Hedges, 1999 Divergence time estimates for the early history of animal phyla and the origin of plants, animals and fungi Proc. R. Soc. Lond. B Biol. Sci 266:163-171[ISI][Medline]