How Athila retrotransposons survive in the Arabidopsis genome
© Marco and Marín; licensee BioMed Central Ltd. 2008
Received: 05 December 2007
Accepted: 14 May 2008
Published: 14 May 2008
Transposable elements are selfish genetic sequences which only occasionally provide useful functions to their host species. In addition, models of mobile element evolution assume a second type of selfishness: elements of different familes do not cooperate, but they independently fight for their survival in the host genome.
We show that recombination events among distantly related Athila retrotransposons have led to the generation of new Athila lineages. Their pattern of diversification suggests that Athila elements survive in Arabidopsis by a combination of selfish replication and of amplification of highly diverged copies with coding potential. Many Athila elements are non-autonomous but still conserve intact open reading frames which are under the effect of negative, purifying natural selection.
The evolution of these mobile elements is far more complex than hitherto assumed. Strict selfish replication does not explain all the patterns observed.
Mobile elements are selfish genomic parasites that only rarely benefit their hosts [1–4]. They belong to two main classes, with or without RNA intermediates, and most eukaryotic genomes contain several types or families of elements of each class [5–7]. A family is a set of very similar sequences that generally include some active elements plus a variable number of non-autonomous, defective copies derived from the active ones. Classical mobile element evolution models suggested that selfishness drives the evolution of each family. Altruistically amplifying either defective copies or elements of other families would decrease the likelihood of long-term survival for a family of elements [5, 6, 8]. The available data for the Saccharomyces cerevisiae and Drosophila melanogaster genomes, among others, in which the rule is to find families of recent origin, composed by almost identical and highly active elements [9, 10], agrees well with those models. However, whether elements that pervade other genomes, especially those with larger amounts of repetitive sequences, follow the same dynamics has been less extensively studied. In fact, the replication of some types of non-autonomous sequences (e. g. SINEs, MITEs, probably several types of retrotransposon-derived plant repeats) present in large numbers in some genomes depend on mobile elements (reviewed in [11–13]). It is not obvious what kind of advantage may obtain the mobile elements involved, and therefore those non-autonomous sequences are considered to replicate parasitically. However, it is possible to envisage situations in which the amplification of non-autonomous elements contributes to the survival of active elements, a possibility that remains largely unexplored. Some evidence for such type of cooperation within a family is available. For example, active Drosophila P elements may improve their likelihood of survival by replicating particular types of defective elements that negatively control the transposition rates of the active ones, thus diminishing the harmful effects on the host (reviewed in ; see also  for related examples).
Athila is one of the best-known plant long-terminal-repeat (LTR) retrotransposons [16–20]. It belongs to the Ty3/Gypsy group, evolutionary closely related to mammalian retroviruses . Actually, some Athila retrotransposons and a few related plant elements are structurally identical to simple retroviruses. They have, in addition to their gag and pol genes, a third ORF, generally absent in other LTR retrotransposons. It may encode an envelope (Env) protein, potentially able to allow the generation of viral infective particles ([17, 19, 20]; see review ). However, whether Athila behaves as an infective retrovirus is still unknown. The evolution of Athila retrotransposons has been traced back using phylogenetic analyses based on their reverse transcriptase (RT) sequences, which are part of the pol gene [17, 18, 20, 22, 23]. These analyses demonstrated that Athila elements are highly heterogenous. Particularly, our group showed that Athila RTs are more variable than those of other eight lineages of Arabidopsis Ty3/Gypsy retrotransposons and that there is no relationship between the degree of similarity among elements and the pattern of presence or absence of env sequences, suggesting that Athila evolution follows a complex pattern .
In this study, we show that the combined analyses of Athila gag, env and pol sequences provides a novel view of the evolutionary forces acting on these retrotransposons in the Arabidopsis genome. We determine that most Athila elements lack pol sequences and therefore are non-autonomous. Some of these elements have however retained intact ORFs that encode for Gag and Env proteins. These ORFs are under the effect of negative, purifying selection and therefore they must be functional. Moreover, diversification and survival of Athila elements in Arabidopsis has often involved recombination among distantly-related elements. In one particular case, recombination involving non-autonomous elements has contributed to generate an active element that moreover has acquired a typical retroviral structure. These results are not compatible with the simplistic view of selfish amplification of independent Athila families.
Arabidopsis Athila elements can be divided into ancient families, many of them exclusively composed by non-autonomous elements
As already indicated above, the evolutionary analyses performed so far on Athila elements have been focused on comparing RT sequences. However, when we deeply examined the diversity of Athila elements, we detected that the analysis of pol-derived sequences may offer at most a partial view of the patterns of evolution of these elements. We found that many Athila elements are characterized by either of two alternative structures, typical of non-autonomous retrotransposons: 1) LTRs plus a single ORF encoding Gag proteins, or, 2) LTRs plus two ORFs, encoding Gag and Env proteins. We also found that all potentially autonomous Athilas, those with pol sequences (including RTs), also have gag sequences, although they may or may not have env sequences.
These results led us to the idea of reassessing Athila evolution from the point of view of their gag sequences. We reasoned that gag sequences, common to all types of both complete and non-autonomous elements, would provide the most precise picture of the evolutionary history of Athila retrotransposons. We thus built phylogenetic trees based on Athila gag sequences. We must note here that in a previous study, based on RT sequences, Athila and the closest relative of Athila, the env-lacking retrotransposon that we named Little Athila  were confounded . However, the recent addition of many novel sequences allowed us to confirm that Athila and Little Athila elements are not only often structurally different (Athilas often contain env sequences, while Little Athilas always lack env), but also possess very different sequences and thus are better defined as two different elements. Particularly, we found that they appear as two separate lineages not only in Arabidopsis, but also in species of the Brassica genus. This result demonstrates that Athila and Little Athila split at least 15–20 millions of years ago (our results are summarized in ). This result was also found by Zhang and Wessler  in their general comparison of the elements present in Arabidopsis and Brassica. Those authors also considered Athila and Little Athila as two different elements. Thus, all the subsequent results shown here refer solely to Athila elements, as defined by Marín and Lloréns  and Zhang and Wessler .
Canonical Athila retroelements. The numbers refer to the nucleotides of each sequence that correspond to Athila ORFs or LTRs.
Locations of the ORFs
Locations of the LTRs
Insertion range‡ (Myr)
0.07 ± 0.03 – 1.60 ± 0.33
0.33 ± 0.10 – 1.63 ± 0.20
0.80 ± 0.13 – 2.07 ± 0.33
1.77 ± 0.20 – 2.07 ± 0.23
1.00 ± 0.13 – 2.40 ± 0.23
0.17 ± 0.07 – 1.07 ± 0.13
0.20 ± 0.10 – 1.63 ± 0.20
0.07 ± 0.03 – 1.73 ± 0.23
0.80 ± 0.13 – 2.33 ± 0.33
0.20 ± 0.10 – 1.83 ± 0.23
0.17 ± 0.07 – 0.93 ± 0.13
Most significantly, only four of the eleven gag-defined families (I, IIIb, IVb and VII) contained elements with pol sequences (see Table 1). These results suggest that most Athila elements, in fact complete families, are non-autonomous, and that they propagate by using the enzymatic machinery provided by elements of other families. Comparative analyses of LTRs demonstrated that non-autonomous families have been multiplying in the genome for periods of time of up to 2 millions of years (Table 1). Obviously, these results also show that all accounts of Athila element evolution published so far, based on RT sequences, offered a very incomplete view of the evolutionary dynamics of this complex ensemble of retroelements.
Activity of Athila retrotransposons
If we assume that the available sequences correctly represent the diversity of the Athila elements present in Arabidopsis thaliana, we may infer the degree of activity of the different families of elements by their number of active copies. Significantly, most Athila sequences are non functional. Out of the almost 200 sequences of Athila elements analyzed, we detected only 10 potentially active elements, which contained ORFs without any frameshifts or stop codons. These elements also contain all characteristic conserved amino acids of Athila Gag proteins and, those that contain pol sequences, also contain the typical motifs of the active centers of reverse transcriptases and integrases. The 10 elements belonged to seven different families, as follows: 1 element from family II, 1 element from family IVb, 2 elements from family IVc, 1 element from family Va, 1 element from family Vb, 3 elements from family VI and 1 element from family VII codons (see arrows in Figure 1). Interestingly, no element in four of those families (II, IVc, Va and VI) has pol sequences, they just contain gag or gag + env sequences. Thus, only three copies among all Athila elements found so far are potentially autonomous, pol-containing copies. Two of them are from family IVb – corresponding to the "Athila4" element already characterized as potentially autonomous by Marín and Lloréns  and Wright and Voytas  – and family VII, respectively. These two families contain other elements with pol sequences, albeit defective. The third one belongs, according to its gag sequence, to family Vb, but, surprisingly, all but two elements in this family lack pol sequences. These peculiar elements, named Va-rec in Figure 1, will be discussed in detail in the next section.
LTR comparative analyses showed that the youngest elements in three of the families without potentially active copies, IIIa, IIIb and IVa, retrotransposed 0.8, 1.8 and 1.0 millions of years ago respectively (Table 1). This result may imply that these families are currently extinct. However, the presence of active Athila elements of these families in other Arabidopsis genomes cannot be excluded. For the fourth family without active copies (family I), a very recent insertion (estimated to have occurred 0.07 ± 0.03 millions of years ago; Table 1) was detected, suggesting that this family is still active. On the other hand, the most recent copies of the families with potentially active elements are in general quite young (average: 0.28 ± 0.09 millions of years) suggesting that most or perhaps all of them are still currently replicating.
Results of selective regime analyses. "x", "y" and "z" refer to the three elements analyzed, with "z" being the one with coding potential, "y" a very close relative and "x" a more distant relative. In all cases except Va-rec, all elements in each analysis belong to the same family. For Va-rec, elements of the two families that give rise to the element were used. In this case, the gag sequences were not analyzed, due to the fact that they are of recombinant origin.
ω in each branch
Recombination among elements of distantly related families
2) Acquisition by a family IIIb element of env sequences originated from a family Va element. This event explains the shift in the position of family IIIb elements in the gag- and env-based trees (Figures 4A, 5B).
3) Recombination between elements of the IVc and VII families, to give rise to family VI. Family VI elements have LTRs and part of the gag sequences that are extremely similar to family VII elements, while the rest of the gag and the env sequences are very similar to those in family IVc (Figure 5C)
4) Acquisition of some family IIIa elements of an env of uncertain origin, generating an additional branch of elements in the env-based tree, that we named IIIa-rec (Figure 5D).
In summary, these results demonstrate that recombination between elements of different families has occurred frequently in the past: at least 4 of the 13 lineages observed in this study (i. e. the 11 families described in Table 1 plus the IIIa-rec and Va-rec lineages, which cannot be detected in gag-based trees), are of recombinant origin. This is probably an understimate, because ancient recombination events or those involving short sequences would remain undetected with our methods. In any case, recombination has been so frequent that none of the phylogenetic trees obtained properly reflected the diversity of Athila retrotransposons. Only tree comparisons allowed us to understand the evolution of these elements.
Selective pressures acting on Athila retroelements
The fact that 70% of the potentially active copies encode for Gag or Env proteins but not for Pol proteins raises the question of whether the pol-less elements are simply parasites of the pol-containing copies or, alternatively, they may be contributing to their own propagation or to the propagation of other Athila elements. This contribution would require the production of active Gag or Env proteins by the non-autonomous elements. Of course, to conclude that these elements may contribute functional proteins is not enough to find out that the non-autonomous copies contain potentially coding ORFs or finding ESTs derived from these elements. Even then, all they could be propagating strictly in a parasitic way, i. e. depending solely on proteins provided in trans by other elements with their own genes being non-functional or fully repressed.
Summary of ESTs derived from Athila elements
Most similar genomic DNA
EG477448, EG477439, EG477443, EG477438, EG477426
EG526171, EG526175, EG526169, EG526179, EG526177, EG526176, EG526150
EG462911, EG455888, EG462908
EG462905, EG462904, EG455886
EG463619, EG463628, EG463610, EG484314, EG484306, EG484300, EG461948, EG461965, EG461956, EG484310, EG463626, EG484303, EG461951, EG484304, EG484316, EG484314, EG484313, EG484299, EG484318, EG463627, EG463617, EG461953, EG459865, EG459891, EG461968, EG461958, EG461955, EG461949, EG463621, EG463625, EG463616, EG461967, EG484315, EG484298, EG459890, EG463611
EG452894, EG452887, EG452883, EG452874, EG452872, EG452892, EG452885, EG452879, EG452900, EG452890, EG452873, EG452871
EG491254, EG491238, EG491236, EG491239, EG491247, EG491249, EG491252, EG491244, EG491246, EG491253, EG491242
EG447146, EG446192, EG448096, EG418344
EG526117, EG526154, EG526119
EG479658, EG504857, EG479617, EG504777, EG479662, EG479648, EG479652, EG479650, EG479663, EG479611, EG479654, EG479607, EG504604, EG504604, EG479605, EG479657, EG504582, EG479655, EG479604, EG479649, EG479646, EG479644, EG479661, EG504856, EG479643, EG479651, EG479608
EG472931, EG472931, EG459693, EG459691, BP823996, BP822138, BP819107, BP826056, EG491235, EG491241, EG423786, EG491243, EG491248, EG491250, EG491245, EG491237, EG459228, EG459225, EG479658, EG504857, EG479617, EG504777, EG479662, EG479648, EG479652, EG479650, EG479663, EG479611, EG479654, EG479607, EG504604, EG479605, EG479657, EG504582, EG479655, EG479604, EG479649, EG479646, EG479644, EG479661, EG504856, EG479643, EG479651, EG479608, EG459688, EG459694, EG459692, EG459690
We may now recapitulate the observations described in the previous chapter. First, we have shown that Athila is composed by at least 11 different families, defined as monophyletic groups of closely related elements (Figure 1). Most of these families emerged in the distant past. We dated the splits between families as having occurred at least 2.7 millions of years ago. Even considering all types of elements, autonomous or not, Athilas are not present in large numbers. Our results agree very well with a previous estimation of about 200 structurally intact copies of Athila per genome . There are no predominant families, so the number of elements for each family is low, ranging from 3 to 31 in our dataset and with an average of 13.2 ± 2.8 (see also Figure 1). Finally, there are only 3 potentially autonomous, active copies in the whole dataset. All these results, together with the low number of ESTs found, suggest that Athila activity is very low. If our results can be extrapolated to other Arabidopsis genomes, we can conclude that Athila as a whole is not a particularly succesful parasite, i. e. it survives at low numbers, and that individual Athila families are at the verge of extinction, at least in individual genomes (although perhaps they are doing fine in the whole species).
Our second main result is that we have shown that families that contain very similar elements without pol sequences have been replicating in the Arabidopsis genome for more than 2 millions of years. These elements are characterized by containing gag or gag + env sequences and several copies have kept apparently intact ORFs with coding potential and which are under a purifying selection regime. Thus, the proteins derived from elements of these families may be contributing to its own replication or to the replication of other Athila elements. Finally, the third main result is that recombination between elements of distantly related families is relatively frequent.
These results are quite different from those observed for most other LTR retrotransposons. For example, in the thoroughly analyzed Saccharomyces cerevisiae, Caenorhabditis elegans or Drosophila melanogaster genomes, most LTR retrotransposons are active, and there are no descriptions of abundant non-autonomous copies with coding potential [9, 10, 26]. Recombination between elements of different families of LTR retrotransposons (Ty1 and Ty2) leading to the generation of a new lineage (Ty1/Ty2) was first observed in S. cerevisiae [27, 28]. However, differently for what we have found for recombinant Athilas, all Ty1/Ty2 recombinant copies are recent  so their long-term evolutionary potential is unclear. Similar cases, in which new lineages of active elements are produced by recombination, have been described in other species [30, 31]. Finally, a case in which a novel non-autonomous element that retains coding potential has emerged by recombination has been described in Hordeum [32, 33], but, again, this element is very young and therefore their ability to propagate for long periods of time is unknown. Significantly, recombination leading to ORFs encoding for "hybrid" proteins of mixed origin, as occurred in Athila families Va-rec and VI (Figures 5A, 5C) was not found in any of these cases.
We may now ask what are the evolutionary processes that explains the particular pattern of evolution observed for Athila elements. First, we may consider whether our results are compatible with the hypothesis of full evolutionary independence of Athila families. To consider fully independent those families for which we have found only non-autonomous copies, we ought to hypothesize that hitherto undiscovered autonomous copies exist for those lineages. This is formally possible but very unlikely. These copies should be promoting the expansion of highly similar, structurally identical defective copies for periods of millions of years while not leaving any detectable pol-containing remnant in the genome. The best argument against this happening is that such peculiar pattern is never observed for the families that do have pol sequences. That is, although many copies in families with pol-containing elements are defective – having accumulated stop codons and frameshifts –, we never observed pol-less elements within those families. We may thus reason that, if in families for which pol-containing elements are known, we never detect a set of related pol-less elements, it is highly unlikely that precisely in those families for which we have not detected pol-containing elements, they actually exist. Therefore, the simplest explanation for the observed pattern is that pol-less elements are mobilized, at least in part, in trans, by enzymes provided by elements that belong to different, pol-containing, families.
We may then ask whether this is just another case of parasitism in which non-autonomous copies use the enzymatic machinery of the active ones without providing any compensation or, alternatively, some kind of cooperation between autonomous and non-autonomous elements might exist. There are two ways in which such cooperation may arise. First, non-autonomous elements with coding potential could contribute to the replication of autonomous copies. To demonstrate this process would require direct biochemical analyses, which is beyond the scope of this work. Our data show however that two necessary conditions for the process to occur are present: 1) there are non-autonomous elements with coding potential, with proteins which are under negative selective pressures; and, 2) the products of distant Athilas are biochemically compatible, as it is demonstrated by the emergence, by recombination, of new families with genes of different origin.
The second way in which cooperation might arise is indirect: generation of coding, non-autonomous copies could be advantageous for the long-term survival of Athila elements as a whole, if the non-autonomous copies occassionally contribute to the generation of novel successful families. Our results demonstrate that this type of event has occurred. We have shown that Athila autonomous and non-autonomous families are linked by recombination events and that several successful recombinant Athila lineages, defined as lineages able to replicate and survive for long periods of time, have arisen. They are of three different types: 1) novel autonomous lineages such as the Va-rec elements; 2) non-autonomous recombinant lineages that have survived while one or perhaps both progenitor families have become extinct, as seems to be the case for the family that provided the env sequences now found in IIIa-rec elements (Figure 5D); or, 3) simply recombinant non-autonomous lineages that are able to propagate in the genome as efficiently as autonomous ones (e. g. family VI, which has been replicating for at least two millions of years). Among all these results, it is most interesting that we may have detected the birth of a new evolutionary entity: if env sequences indeed provide Athila elements with the possibility of becoming infective, Va-rec elements would be an example of how recombination between an autonomous retrotransposon (env-less, from family VII) and a non-autonomous element (pol-less, from family Va) generates a novel active retrovirus (with gag, pol and env; Figure 5A). In any case, these events demonstrate that non-autonomous copies are not strictly parasitic. They are contributing to the long-term survival of Athila elements in Arabidopsis.
In summary, our results suggest that distant Athila families may be cooperating to survive in the Arabidopsis genome. Cooperation among other type of mobile elements, bacterial IS elements, has recently received attention, with the conclusion that it may appear under precise selective regimes . Recent models also suggest situations in which mutualism may occur . In fact, we think that the accepted view that all elements behave strictly selfishly may be due to the great difficulties involved in discovering patterns of sequence evolution compatible with cooperative processes. It is possible that other mobile elements follow dynamics similar to the one we have just described. For example, evidence for a related pattern of interchange to generate novel lineages is also available for human endogenous retroviruses ([36, 37]; see also discussion in ). Therefore, this may be the first formal description of a widely-used survival strategy for eukaryotic mobile elements.
Data mining and phylogenetic analyses
We built databases of Gag, reverse transcriptase and Env proteins using BlastP and TblastN searches against the databases available at the National Center for Biotechnology Information (NCBI). We used as queries multiple representative Athila elements until the searches become saturated. After each search, we aligned the sequences obtained and removed duplicates and partial sequences. All alignments were performed with ClustalX 1.83  using default parameters. Alignments were manually corrected when necessary with GeneDoc 2.6 . We used two methods of phylogenetic inference, neighbor-joining and maximum parsimony, implemented in MEGA2  following the methods described in . For both methods, statistical support for the branches was assessed performing 1000 bootstrap replicates.
To determine the structure of Athila elements, we first used Blast2sequences searches , comparing each element with known gag, RT, integrase and env Athila sequences. ORFs were detected with ORF finder . LTR locations were determined by looking for similarity within an element, also with Blast2sequences.
Estimation of the insertion time or divergence time between elements
We estimated the insertion time for an element or the divergence time for elements of two different families following the strategy described by San Miguel et al. . The nucleotide sequences of either both LTRs of each element or single LTRs of two different elements were aligned and the Kimura two-parameter distance  was estimated using MEGA2. The distance obtained was then divided by two (because it refers to changes accumulated in both LTRs) and then again divided by the substitution rate at synonymous sites estimated for the brassicaceae Chs and Adh genes , that is, 1.5 10-8 per site per year.
Characterization of recombination events
Recombination events were deduced from incongruent phylogenetic positions of two proteins of a same element or group of elements . We searched for the recombination breakpoints by analyzing pairwise alignments of amino acidic sequences and also, at the nucleotide level, following a sliding-window approach implemented in SimPlot , which utilizes the DNAPARS and NEIGHBOR programs of the Phylip package .
Estimations fo synonymous and nonsynonymous nucleotide substitutions
The rate of synonymous and nonsynonymous substitution were estimated using the PAML3.1 package , following the strategy of branch-dependent analyses described previously in . Three models were analyzed. M0 refers to a model in which all branches are assumed to have the same rate. M1 is a model in which all branches are assumed to evolve at different rates. Finally, M2 is a model in which the branch that leads to the active element ("z" in Table 2) is assumed to evolve at a different rate than the rest.
Research supported by grant 200720I021 (Proyectos intramurales especiales, CSIC. Spain).
- Hurst GD, Werren JH: The role of selfish genetic elements in eukaryotic evolution. Nat Rev Genet. 2001, 2: 597-606. 10.1038/35084545.PubMedView ArticleGoogle Scholar
- Kidwell MG, Lisch DR: Transposable elements as sources of genomic variation. Mobile DNA II. Edited by: Craig NL, Craigie R, Gellert M, Lambowitz AM. 2002, Washington: ASM Press, 59-90.View ArticleGoogle Scholar
- Labrador M, Corcés V: Interactions between transposable elements and the host genome. Mobile DNA II. Edited by: Craig NL, Craigie R, Gellert M, Lambowitz AM. 2002, Washington: ASM Press, 1008-1023.View ArticleGoogle Scholar
- Brookfield JFY: The ecology of the genome – Mobile DNA elements and their hosts. Nat Rev Genet. 2005, 6: 128-136. 10.1038/nrg1524.PubMedView ArticleGoogle Scholar
- Charlesworth B, Langley CH: The population genetics of Drosophila transposable elements. Annu Rev Genet. 1989, 23: 251-287. 10.1146/annurev.ge.23.120189.001343.PubMedView ArticleGoogle Scholar
- Charlesworth B, Sniegowski P, Stephan W: The evolutionary dynamics of repetitive DNA in eukaryotes. Nature. 1994, 371: 215-220. 10.1038/371215a0.PubMedView ArticleGoogle Scholar
- Kazazian HHJ: Mobile elements: drivers of genome evolution. Science. 2004, 303: 1626-1632. 10.1126/science.1089670.PubMedView ArticleGoogle Scholar
- Kaplan N, Darden T, Langley CH: Evolution and extinction of transposable elements in mendelian populations. Genetics. 1985, 109: 459-480.PubMedPubMed CentralGoogle Scholar
- Lerat E, Rizzon C, Biémont C: Sequence divergence within transposable element families in the Drosophila melanogaster genome. Genome Res. 2003, 13: 1889-1896.PubMedPubMed CentralGoogle Scholar
- Lesage P, Todeschini AL: Happy together: the life and times of Ty retrotransposons and their hosts. Cytogenet Genome Res. 2005, 110: 70-90. 10.1159/000084940.PubMedView ArticleGoogle Scholar
- Batzer MA, Deininger PL: Alu repeats and human genomic diversity. Nat Rev Genet. 2002, 3: 370-379. 10.1038/nrg798.PubMedView ArticleGoogle Scholar
- Feschotte C, Jiang N, Wessler SR: Plant transposable elements: where genetics meets genomics. Nat Rev Genet. 2002, 3: 329-341. 10.1038/nrg793.PubMedView ArticleGoogle Scholar
- Sabot F, Schulman AH: Parasitism and the retrotransposon life cycle in plants: a hitchhiker's guide to the genome. Heredity. 2006, 96: 381-388. 10.1038/sj.hdy.6800903.View ArticleGoogle Scholar
- Rio DR: P transposable elements in Drosophila melanogaster. Mobile DNA II. Edited by: Craig NL, Craigie R, Gellert M, Lambowitz AM. 2002, Washington: ASM Press, 484-518.View ArticleGoogle Scholar
- Leonardo TE, Nuzhdin SV: Intracellular battlegrounds: conflict and cooperation between transposable elements. Genet Res. 2002, 80: 155-161. 10.1017/S0016672302009710.PubMedView ArticleGoogle Scholar
- Pélissier T, Tutois S, Deragon JM, Tourmente S, Genestier S, Picard G: Athila, a new retroelement from Arabidopsis thaliana. Plant Mol Biol. 1995, 29: 441-452. 10.1007/BF00020976.PubMedView ArticleGoogle Scholar
- Wright DA, Voytas DF: Potential retroviruses in plants: tat1 is related to a group of Arabidopsis thaliana Ty3/gypsy retrotransposons that encode envelope-like proteins. Genetics. 149: 703-715.
- Marín I, Lloréns C: Ty3/gypsy retrotransposons: description of new Arabidopsis thaliana elements and evolutionary perspectives derived from comparative genomic data. Mol Biol Evol. 2000, 17: 1040-1049.PubMedView ArticleGoogle Scholar
- Vicient CM, Kalendar R, Schulman AH: Envelope-class retrovirus-like elements are widespread, transcribed and spliced, and insertionally polymorphic in plants. Genome Res. 2001, 11: 2041-2049. 10.1101/gr.193301.PubMedPubMed CentralView ArticleGoogle Scholar
- Wright DA, Voytas DF: Athila4 of Arabidopsis and calypso of soybean define a lineage of endogenous plant retroviruses. Genome Res. 2002, 12: 122-131. 10.1101/gr.196001.PubMedPubMed CentralView ArticleGoogle Scholar
- Xiong Y, Eickbush TH: Origin and evolution of retroelements based upon their reverse transcriptase sequences. EMBO J. 1990, 9: 3353-3362.PubMedPubMed CentralGoogle Scholar
- Marco A, Marín I: Retrovirus-like elements in plants. Recent Res Devel Plant Sci. 2005, 3: 15-24.Google Scholar
- Zhang X, Wessler SR: Genome-wide comparative analysis of the transposable elements in the related species Arabidopsis thaliana and Brassica oleracea. Proc Natl Acad Sci USA. 2004, 101: 5589-5594. 10.1073/pnas.0401243101.PubMedPubMed CentralView ArticleGoogle Scholar
- Chan SW-L, Henderson IR, Jacobsen SE: Gardening the genome: DNA methylation in Arabidopsis thaliana. Nat Rev Genet. 2005, 6: 351-360. 10.1038/nrg1601.PubMedView ArticleGoogle Scholar
- Pereira V: Insertion bias and purifying selection of retrotransposons in the Arabidopsis thaliana genome. Genome Biol. 2004, 5: R79-10.1186/gb-2004-5-10-r79.PubMedPubMed CentralView ArticleGoogle Scholar
- Bowen NJ, McDonald JF: Genomic analysis of Caenorhabditis elegans revelas ancient families of retroviral-like elements. Genome Res. 1999, 9: 924-935. 10.1101/gr.9.10.924.PubMedView ArticleGoogle Scholar
- Jordan IK, McDonald JF: Evidence for the role of recombination in the regulatory evolution of Saccharomyces cerevisiae Ty elements. J Mol Evol. 1998, 47: 14-20. 10.1007/PL00006358.PubMedView ArticleGoogle Scholar
- Jordan IK, McDonald JF: Phylogenetic perspective reveals abundant Ty1/Ty2 hybrid elements in the Saccharomyces cerevisiae genome. Mol Biol Evol. 1999, 16: 419-422.PubMedView ArticleGoogle Scholar
- Jordan IK, McDonald JF: Tempo and mode of Ty element evolution in Saccharomyces cerevisiae. Genetics. 1999, 151: 1341-1351.PubMedPubMed CentralGoogle Scholar
- Kelly FD, Levin HL: The evolution of transposons in Schizosaccharomyces pombe. Cytogenet Genome Res. 2005, 110: 566-574. 10.1159/000084990.PubMedView ArticleGoogle Scholar
- Mugnier N, Biémont C, Vieira C: New regulatory regions of Drosophila 412 retrotransposable element generated by recombination. Mol Biol Evol. 2004, 22: 747-757. 10.1093/molbev/msi060.PubMedView ArticleGoogle Scholar
- Vicient CM, Kalendar R, Schulman AH: Variability, recombination, and mosaic evolution of the barley BARE-1 retrotransposon. J Mol Evol. 2005, 61: 275-291. 10.1007/s00239-004-0168-7.PubMedView ArticleGoogle Scholar
- Tanskanen JA, Sabot F, Vicient C, Schulman AH: Life without GAG: the BARE-2 retrotransposon as a parasite's parasite. Gene. 2007, 390: 166-174. 10.1016/j.gene.2006.09.009.PubMedView ArticleGoogle Scholar
- Wagner A: Cooperation is fleeting in the world of transposable elements. PloS Comput Biol. 2006, 2: e162-10.1371/journal.pcbi.0020162.PubMedPubMed CentralView ArticleGoogle Scholar
- Le Rouzic A, Dupas S, Capy P: Genome ecosystem and transposable element species. Gene. 2007, 390: 214-220. 10.1016/j.gene.2006.09.023.PubMedView ArticleGoogle Scholar
- Bénit L, Dessen P, Heidmann T: Identification, phylogeny, and evolution of retroviral elements based on their envelope genes. J Virol. 2001, 75: 11709-11719. 10.1128/JVI.75.23.11709-11719.2001.PubMedPubMed CentralView ArticleGoogle Scholar
- Jordan IK, McDonald JF: A biologically active family of human endogenous retroviruses (HERVs) evolved from an ancient inactive lineage. Genome Lett. 2002, 2: 105-109. 10.1166/gl.2002.011.View ArticleGoogle Scholar
- Bannert N, Kurth R: The evolutionary dynamics of human endogenous retroviral families. Annu Rev Genomics Hum Genet. 2006, 7: 149-173. 10.1146/annurev.genom.7.080505.115700.PubMedView ArticleGoogle Scholar
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25: 4876-4882. 10.1093/nar/25.24.4876.PubMedPubMed CentralView ArticleGoogle Scholar
- Nicholas KB, Nicholas HB: Genedoc: analysis and visualization of genetic variation. Distributed by the authors. 1997, [http://www.genedoc.us]Google Scholar
- Kumar S, Tamura K, Jakobsen IB, Nei M: Mega2: molecular evolutionary genetics analysis software. Bioinformatics. 2001, 17: 1244-1245. 10.1093/bioinformatics/17.12.1244.PubMedView ArticleGoogle Scholar
- Marco A, Cuesta A, Pedrola L, Palau F, Marín I: Evolutionary and structural analyses of GDAP1, involved in Charcot-Marie-Tooth disease, characterize a novel class of glutathione transferase-related genes. Mol Biol Evol. 2004, 21: 176-187. 10.1093/molbev/msh013.PubMedView ArticleGoogle Scholar
- Tatusova TA, Madden TL: Blast 2 sequences, a new tool for comparing protein and nucleotide sequences. FEMS Microbiol Lett. 1999, 174: 247-250. 10.1111/j.1574-6968.1999.tb13575.x.PubMedView ArticleGoogle Scholar
- ORF finder web page. [http://www.ncbi.nlm.nih.gov/projects/gorf]
- SanMiguel P, Gaut BS, Tikhonov A, Nakajima Y, Bennetzen JL: The paleontology of intergene retrotransposons of maize. Nat Genet. 1998, 20: 43-45. 10.1038/1695.PubMedView ArticleGoogle Scholar
- Kimura M: A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980, 16: 111-120. 10.1007/BF01731581.PubMedView ArticleGoogle Scholar
- Koch MA, Haubold B, Mitchell-Olds T: Comparative evolutionary analysis of chalcone synthase and alcohol dehydrogenase loci in Arabidopsis, Arabis, and related genera (brassicaceae). Mol Biol Evol. 2000, 17: 1483-1498.PubMedView ArticleGoogle Scholar
- Posada D, Crandall KA, Holmes EC: Recombination in evolutionary genomics. Annu Rev Genet. 2002, 36: 75-97. 10.1146/annurev.genet.36.040202.111115.PubMedView ArticleGoogle Scholar
- Lole KS, Bollinger RC, Paranjape RS, Gadkari D, Kulkarni SS, Novak NG, Ingersoll R, Sheppard HW, Ray SC: Full-length human immunodeficiency virus type 1 genomes from subtype c-infected seroconverters in india, with evidence of intersubtype recombination. J Virol. 1999, 73: 152-160.PubMedPubMed CentralGoogle Scholar
- Felsenstein J: Phylip (phylogeny inference package) version 3.5c. distributed by the author. 1993, Department of Genetics, University of Washington, SeattleGoogle Scholar
- Yang Z: Paml: a program package for phylogenetic analysis by maximum likelihood. Comput Appl Biosci. 1997, 13: 555-556.PubMedGoogle Scholar
- Fares MA, Bezemer D, Moya A, Marín I: Selection on coding regions determined Hox7 genes evolution. Mol Biol Evol. 2003, 20: 2104-2112. 10.1093/molbev/msg222.PubMedView ArticleGoogle Scholar