Complex evolution of gypsy in Drosophilid species

In an endeavor to contribute to the comprehension of the evolution of transposable elements (TEs) in the genome of host species, we investigated the phylogenetic relationships of sequences homologous to the retrotransposon gypsy of Drosophila melanogaster in 19 species of Drosophila, in Scaptodrosophila latifasciaeformis, and in Zaprionus indianus. This phylogenetic study was based on approximately 500 base pairs of the env gene. Our analyses showed considerable discrepancy between the phylogeny of gypsy elements and the relationship of their host species, and they allow us to infer a complex evolutionary pattern that could include ancestral polymorphism, vertical transmission, and several cases of horizontal transmission.     

Evolution of the Transposable Element Mariner in Drosophila Species

The distribution of the transposable element mariner was examined in the genus Drosophila. Among the eight species comprising the melanogaster species subgroup, the element is present in D. mauritiana, D. simulans, D. sechellia, D. yakuba and D. teissieri, but it is absent in D. melanogaster, D. erecta and D. orena. Multiple copies of mariner were sequenced from each species in which the element occurs. The inferred phylogeny of the elements and the pattern of divergence were examined in order to evaluate whether horizontal transfer among species or stochastic loss could better account for the discontinuous distribution of the element among the species. The data suggest that the element was present in the ancestral species before the melanogaster subgroup diverged and was lost in the lineage leading to D. melanogaster and the lineage leading to D. erecta and D. orena. This inference is consistent with the finding that mariner also occurs in members of several other species subgroups within the overall melanogaster species group. Within the melanogaster species subgroup, the average divergence of mariner copies between species was lower than the coding region of the alcohol dehydrogenase (Adh) gene. However, the divergence of mariner elements within species was as great as that observed for Adh. We conclude that the relative sequence homogeneity of mariner elements within species is more likely a result of rapid amplification of a few ancestral elements than of concerted evolution. The mariner element may also have had unequal mutation rates in different lineages.     

Transposable Element Dynamics in Two Sibling Species: Drosophila Melanogaster and Drosophila Simulans

Transposable elements (TEs) in the two sibling species, Drosophila melanogaster and D. simulans, differ considerably in amount and dynamics, with D. simulans having a smaller amount of TEs than D. melanogaster. Several hypotheses have been proposed to explain these differences, based on the evolutionary history of the two species, and claim differences either in the effective size of the population or in genome characteristics. Recent data suggest, however, that the higher amount of TEs in D. melanogaster could be associated with the worldwide invasion of D. melanogaster a long time ago while D. simulans is still under the process of such geographical spread. Stresses due to new environmental conditions and crosses between migrating populations could explain the mobilization of TEs while the flies colonize. Colonization and TE mobilization may be strong evolutionary forces that have shaped and are still shaping the eukaryote genomes.     

Horizontal transfer of hobo transposable elements within the Drosophila melanogaster species complex: evidence from DNA sequencing.

The hobo family of transposable elements, one of three transposable-element families that cause hybrid dysgenesis in Drosophila melanogaster, appears to be present in all members of the D. melanogaster species complex: D. melanogaster, D. simulans, D. mauritiana, and D. sechellia. Some hobo-hybridizing sequences are also found in the other members of the melanogaster subgroup and in many members of the related montium subgroup. Surveys of older isofemale lines of D. melanogaster suggest that complete hobo elements were absent prior to 50 years ago and that hobo has recently been introduced into the species by horizontal transfer. To test the horizontal transfer hypothesis, the 2.6-kb XhoI fragments of hobo elements from D. melanogaster, D. simulans, and D. mauritiana were cloned and sequenced. The DNA sequences reveal an extremely low level of divergence and support the conclusion that the active hobo element has been horizontally transferred into or among these species in the recent past.     

Molecular evolution of P transposable elements in the genus Drosophila. III. The melanogaster species group

Phylogenetic relationships were determined for 76 partial P-element sequences from 14 species of the melanogaster species group within the Drosophila subgenus Sophophora. These results are examined in the context of the phylogeny of the species from which the sequences were isolated. Sequences from the P-element family fall into distinct subfamilies, or clades, which are often characteristic for particular species subgroups. When examined locally among closely related species, the evolution of P elements is characterized by vertical transmission, whereby the P-element phylogeny traces the species phylogeny. On a broader scale, however, the P-element phylogeny is not congruent with the species phylogeny. One feature of P-element evolution in the melanogaster group is the presence of more than one P-element subfamily, differing by as much as 36%, in the genomes of some species. Thus, P elements from several individual species are not monophyletic, and a likely explanation for the incongruence between P-element and species phylogenies is provided by the comparison of paralogous sequences. In certain instances, horizontal transfer seems to be a valid alternative explanation for lack of congruence between species and P-element phylogenies. The canonical P-element subfamily, which represents the active, autonomous transposable element, is restricted to D. melanogaster. Thus, its origin clearly lies outside of the melanogaster species group, consistent with the earlier conclusion of recent horizontal transfer.      

New regulatory regions of Drosophila 412 retrotransposable element generated by recombination

There are no doubts that transposable elements (TEs) have greatly influenced genomes evolution. They have, however, evolved in different ways throughout mammals, plants, and invertebrates. In mammals they have been shown to be widely present but with low transposition activity; in plants they are responsible for large increases in genome size. In Drosophila, despite their low amount, transposition seems to be higher. Therefore, to understand how these elements have evolved in different genomes and how host genomes have proposed to go around them, are major questions on genome evolution. We analyzed sequences of the retrotransposable elements 412 in natural populations of the Drosophila simulans and D. melanogaster species that greatly differ in their amount of TEs. We identified new subfamilies of this element that were the result of mutation or insertion-deletion process, but also of interfamily recombinations. These new elements were well conserved in the D. simulans natural populations. The new regulatory regions produced by recombination could give rise to new elements able to overcome host control of transposition and, thus, become potential genome invaders.     

Population genetics and phylogenetics of DNA sequence variation at multiple loci within the Drosophila melanogaster species complex

Two regions of the genome, a 1-kbp portion of the zeste locus and a 1.1- kbp portion of the yolk protein 2 locus, were sequenced in six individuals from each of four species: Drosophila melanogaster, D. simulans, D. mauritiana, and D. sechellia. The species and strains were the same as those of a previous study of a 1.9-kbp region of the period locus. No evidence was found for recent balancing or directional selection or for the accumulation of selected differences between species. Yolk protein 2 has a high level of amino acid replacement variation and a low level of synonymous variation, while zeste has the opposite pattern. This contrast is consistent with information on gene function and patterns of codon bias. Polymorphism levels are consistent with a ranking of effective population sizes, from low to high, in the following order: D. sechellia, D. melanogaster, D.mauritiana, and D. simulans. The apparent species relationships are very similar to those suggested by the period locus study. In particular, D. simulans appears to be a large population that is still segregating variation that arose before the separation of D. mauritiana and D. sechellia. It is estimated that the separation of ancestral D. melanogaster from the other species occurred 2.5-3.4 Mya. The separations of D. sechellia and D. mauritiana from ancestral D. simulans appear to have occurred 0.58- 0.86 Mya, with D. mauritiana having diverged from ancestral D. simulans 0.1 Myr more recently than D. sechellia.      

hobo transposable elements in Drosophila melanogaster and D. simulans.

Genomic patterns of occurrence of the transposable element hobo are polymorphic in the sibling species Drosophila melanogaster and D. simulans. Most tested strains of both species have apparently complete (3.0 kb) and smaller hobo elements (H lines), but in both species some strains completely lack such canonical hobo elements (E lines). The occurrence of H and E lines in D. simulans as well as in D. melanogaster implies that an hypothesis of recent introduction in the latter species is inadequate to explain the phylogenetic occurrence of hobo. Particular internally deleted elements, the approximately 1.5 kb Th1 and Th2 elements, are abundant in many lines of D. melanogaster, and an analogous 1.1 kb internally deleted element, h del sim, is abundant in most lines of D. simulans. Besides the canonical hobo sequences, both species (and their sibling species D. sechellia and D. mauritiana) have many hobo-hybridizing sequences per genome that do not appear to be closely related to the canonical hobo sequence.     

Distribution of transposable elements in Drosophila species

We present a global analysis of the distribution of 43 transposable elements (TEs) in 228 species of the Drosophila genus from our data and data from the literature. Data on chromosome localization come from in situ hybridization and presence/absence of the elements from southern analyses. This analysis shows great differences between TE distributions, even among closely related species. Some TEs are distributed according to the phylogeny of their host specie; others do not entirely follow the phylogeny, suggesting horizontal transfers. A higher number of insertion sites for most TEs in the genome of D. melanogaster is observed when compared with that in D. simulans. This suggests either intrinsic differences in genomic characteristics between the two species, or the influence of differing effective population sizes, although biases due to the use of TE probes coming mostly from D. melanogaster and to the way TEs are initially detected in species cannot be ruled out. Data on TEs more specific to the species under consideration are necessary for a better understanding of their distribution in organisms and populations. This revised version was published online in July 2006 with corrections to the Cover Date.     

Comparative analysis of transposable elements in the melanogaster subgroup sequenced genomes

Transposable elements (TEs) are indwelling components of genomes, and their dynamics have been a driving force in genome evolution. Although we now have more information concerning their amounts and characteristics in various organisms, we still have little data from overall comparisons of their sequences in very closely-related species. While the Drosophila melanogaster genome has been extensively studied, we have only limited knowledge regarding the precise TE sequences in the genomes of the related species Drosophila simulans, Drosophila sechellia and Drosophila yakuba. In this study we analyzed the number and structure of TE copies in the sequenced genomes of these four species. Our findings show that, unexpectedly, the number of TE insertions in D. simulans is greater than that in D. melanogaster, but that most of the copies in D. simulans are degraded and in small fragments, as in D. sechellia and D. yakuba. This suggests that all three species were invaded by numerous TEs a long time ago, but have since regulated their activity, as the present TE copies are degraded, with very few full-length elements. In contrast, in D. melanogaster, a recent activation of TEs has resulted in a large number of almost-identical TE copies. We have detected variants of some TEs in D. simulans and D. sechellia, that are almost identical to the reference TE sequences in D. melanogaster, suggesting that D. melanogaster has recently been invaded by active TE variants from the other species. Our results indicate that the three species D. simulans, D. sechellia, and D. yakuba seem to be at a different stage of their TE life cycle when compared to D. melanogaster. Moreover, we show that D. melanogaster has been invaded by active TE variants for several TE families likely to come from D. simulans or the ancestor of D. simulans and D. sechellia. The numerous horizontal transfer events implied to explain these results could indicate introgression events between these species.     

Paleo-demography of the Drosophila melanogaster subgroup: application of the maximum likelihood method

The species divergence times and demographic histories of Drosophila melanogaster and its three sibling species, D. mauritiana, D. simulans, and D. yakuba, were investigated using a maximum likelihood (ML) method. Thirty-nine orthologous loci for these four species were retrieved from DDBJ/EMBL/GenBank database. Both autosomal and X-linked loci were used in this study. A significant degree of rate heterogeneity across loci was observed for each pair of species. Most loci have the GC content greater than 50% at the third codon position. The codon usage bias in Drosophila loci is considered to result in the high GC content and the heterogenous rates across loci. The chi-square, G, and Fisher's exact tests indicated that data sets with 11, 23, and 9 pairs of DNA sequences for the comparison of D. melanogaster with D. mauritiana, D. simulans, and D. yakuba, respectively, retain homogeneous rates across loci. We applied the ML method to these data sets to estimate the DNA sequence divergences before and after speciation of each species pair along with their standard deviations. Using 1.6 x 10(-8) as the rate of nucleotide substitutions per silent site per year, our results indicate that the D. melanogaster lineage split from D. yakuba approximately 5.1 +/- 0.8 million years ago (mya), D. mauritiana 2.7 +/- 0.4 mya, and D. simulans 2.3 +/- 0.3 mya. It implies that D. melanogaster became distinct from D. mauritiana and D. simulans at approximately the same time and from D. yakuba no earlier than 10 mya. The effective ancestral population size of D. melanogaster appears to be stable over evolutionary time. Assuming 10 generations per year for Drosophila, the effective population size in the ancestral lineage immediately prior to the time of species divergence is approximately 3 x 10(6), which is close to that estimated for the extant D. melanogaster population. The D. melanogaster did not encounter any obvious bottleneck during the past 10 million years.     

High rate of horizontal transfer of transposable elements in Drosophila

We have conducted molecular population genetics analyses to understand the relationships among the transposable elements (TEs) in Drosophila melanogaster, in combination with sequence comparisons of TEs from two related species, D. simulans and D. yakuba. We observed much lower than expected genetic differences among elements, clear evidence for departure from expectations for equilibrium copy numbers and little divergence between species. This suggests that a large proportion of TEs in D. melanogaster had a recent origin as a result of interspecies movement.     

Evolutionary stability of the R1 retrotransposable element in the genus Drosophila

R1 is a non-long terminal repeat (non-LTR) retrotransposable element that inserts into a specific sequence of insect 28S ribosomal RNA genes. We have previously shown that this element has been maintained through vertical transmission in the melanogaster species subgroup of Drosophila. To address whether R1 elements have been vertically transmitted for longer periods of evolutionary time, the analysis has been extended to 11 other species from four species groups of the genus Drosophila (melanogaster, obscura, testecea, and repleta). All sequenced elements appeared functional on the basis of the preservation of their open-reading frames and consistently higher rate of substitution at synonymous sites relative to replacement sites. The phylogenetic relationships of the R1 elements from all species analyzed were congruent with the species phylogenies, suggesting that the R1 elements have been vertically transmitted since the inception of the Drosophila genus, an estimated 50-70 Mya. The stable maintenance of R1 through the germ line appears to be the major mechanism for the widespread distribution of these elements in Drosophila. In two species, D. neotestecea of the testecea group and D. takahashii of the melanogaster group, a second family of R1 elements was also present that differed in sequence by 46% and 31%, respectively, from the family that was congruent with the species phylogeny. These second families may represent occasional horizontal transfers or, alternatively, they could reflect the ability of R1 elements to diverge into new families within a species and evolve independently.      

Distribution of the retrotransposable element 412 in Drosophila species

Copy numbers of sequences homologous to the Drosophila melanogaster retrotransposable element 412, their distribution between the chromosome arms and the chromocenter, and whether they contain full- size copies were analyzed for 55 species of the Drosophila genus. Element 412 insertion sites were detected on the chromosome arms of D. melanogaster, Drosophila simulans, and a few species of the obscura group, but the chromocenter was labeled in almost all species. The presence of element 412 sequences in the majority of species shows that this element has a long evolutionary history in Drosophilidae, although it may have recently invaded the chromosomes in some species, such as D. simulans. Differences in copy number between species may be due to population size or specific endogenous or environmental factors and may follow the worldwide invasion of the species. Putative full-length copies were detected in the chromocenters of some species with no copies on the chromosome arms, suggesting that the chromocenter may be a shelter for such copies and not only for deleted ones.      

Inferences on the evolutionary history of the S-element family of Drosophila melanogaster

The S-element family of transposable elements has been characterized in D. melanogaster. Attempts to find it in other Drosophila-related species have failed, suggesting that this element family may have recently invaded the D. melanogaster genome by horizontal transfer. In order to investigate its evolutionary history, we analyzed the patterns of DNA polymorphism among the S-element copies present in a sample genome (Drosophila Genome Project). The observed levels of nucleotide diversity are significantly lower than theoretical expectations based on the neutral model. This is consistent with evidence for ongoing gene conversion among copies and for purifying selection on the elements' sequences, particularly on the terminal inverted repeats. A phylogenetic analysis revealed that the members of the S-element family can be grouped into at least two genetically differentiated clusters. The level of divergence between these clusters suggests that the S elements invaded the genome of the ancestor of D. melanogaster before the speciation of the D. melanogaster complex. However, other relevant scenarios are also discussed.     

Molecular identification of the active ninja retrotransposon and the inactive aurora element in Drosophila simulans and D. melanogaster.

How transposable elements evolve is a key facet in understanding of spontaneous mutation and genomic rearrangements in various organisms. One of the best ways to approach this question is to study a newly evolved transposable element whose presence is restricted to a specific population or strain. The retrotransposons ninja and aurora may provide insights into the process of their evolution, because of their contrasting characteristics, even though they show high sequence identity. The ninja retrotransposon was found in a Drosophila simulans strain in high copy number and is potent in transposition. On the other hand, aurora elements are distributed widely among the species belonging to the Drosophila melanogaster species complex, but are immobile at least in D. melanogaster. In order to distinguish the two closely resembled retrotransposons by molecular means, we determined and compared DNA sequence of the elements, and identified characteristic internal deletions and nucleotide substitutions in 5'-long terminal repeats (LTR). Analyses of the structure of ninja homologs and LTR sequences amplified from both genomic and cloned DNA revealed that the actively transposable ninja elements were present only in D. simulans strains, but inactive aurora elements exist in both D. melanogaster and D. simulans.     

Evidence for horizontal transfer of the LTR retrotransposon mdg3, which lacks an env gene

Horizontal (interspecific) transfer is regarded as a possible strategy for the propagation of transposable elements through evolutionary time. To date, however, conclusive evidence that transposable elements are capable of horizontal transfer from one species to another has been limited to class II or DNA-type elements. We tested the possibility of such transfer for several Drosophila melanogaster LTR retrotransposons of the gypsy group in an experiment in which D. melanogaster and D. virilis somatic cell lines were used as donor and recipient cells, respectively. This approach was chosen in light of the high levels of LTR retrotransposon amplification and expression observed in cultured D. melanogaster cells. In the course of the experiment, parallel analysis for mdg1, mdg3, 17.6, 297, 412 and B104/roo retrotransposons was performed to detect their presence in the genome of recipient cells. Only the mdg3 retrotransposon, which lacks an env gene, was found to be transmitted into recipient cells. This model, based on the use of cultured cells, is a promising system for further investigating the mechanisms of LTR retrotransposon transfer.     

Horizontal transfer of non-LTR retrotransposons in vertebrates

Since their discovery in family Bovidae (bovids), Bov-B LINEs, believed to be order-specific SINEs, have been found in all ruminants and recently also in Viperidae snakes. The distribution and the evolutionary relationships of Bov-B LINEs provide an indication of their origin and evolutionary dynamics in different species. The evolutionary origin of Bov-B LINE elements has been shown unequivocally to be in Squamata (squamates). The horizontal transfer of Bov-B LINE elements in vertebrates has been confirmed by their discontinuous phylogenetic distribution in Squamata (Serpentes and two lizard infra-orders) as well as in Ruminantia, by the high level of nucleotide identity, and by their phylogenetic relationships. The direction of horizontal transfer from Squamata to the ancestor of Ruminantia is evident from the genetic distances and discontinuous phylogenetic distribution of Bov-B LINE elements. The ancestral snake lineage (Boidae) has been recognized as a possible donor of Bov-B LINE elements to Ruminantia. The timing of horizontal transfer has been estimated from the distribution of Bov-B LINE elements in Ruminantia and the fossil data of Ruminantia to be 40–50mya. The phylogenetic relationships of Bov-B LINE elements from the various Squamata species agrees with that of the species phylogeny, suggesting that Bov-B LINE elements have been stably maintained by vertical transmission since the origin of Squamata in the Mesozoic era. This revised version was published online in July 2006 with corrections to the Cover Date.     

The Evolutionary History of the Transposable Element Penelope in the Drosophila virilis Group of Species

We have used phylogenetic techniques to study the evolutionary history of the Penelope transposable element in the Drosophila virilis species group. Two divergent types of Penelope have been detected, one previously described, clade I, and a new one which we have termed clade III. The phylogeny of some copies of the Penelope clade I element was partially consistent with the species phylogeny of the D. montana subphylad, suggesting cospeciation and allowing the estimation of the evolutionary rate of Penelope. Divergence times of elements found in different species are younger than the age of the species, suggesting horizontal transfer events. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Dmitri Petrov]