The Foldback-like element Galileo belongs to the P superfamily of DNA transposons and is widespread within the Drosophila genus

Evolutionary History of the DD41D Family of Tc1/Mariner Transposons in Two Mayetiola Species

Tc1/mariner elements are ubiquitous in eukaryotic genomes including insects. They are diverse and divided into families and sub-families. The DD34D family including mauritiana and irritans subfamilies have already been identified in two closely related species of Cecidomyiids M. destructor and M. hordei. In the current study the de novo and similarity-based methods allowed the identification for the first time of seven consensuses in M. destructor and two consensuses in M. hordei belonging to DD41D family whereas the in vitro method allowed the amplification of two and three elements in these two species respectively. Most of identified elements accumulated different mutations and long deletions spanning the N-terminal region of the transposase. Phylogenetic analyses showed that the DD41D elements were clustered in two groups belonging to rosa and Long-TIR subfamilies. The age estimation of the last transposition events of the identified Tc1/mariner elements in M. destructor showed different evolutionary histories. Indeed, irritans elements have oscillated between periods of silencing and reappearance while rosa and mauritiana elements have shown regular activity with large recent bursts. The study of insertion sites showed that they are mostly intronic and that some recently transposed elements occurred in genes linked to putative DNA-binding domains and enzymes involved in metabolic chains. Thus, this study gave evidence of the existence of DD41D family in two Mayetiola species and an insight on their evolutionary history.

Horizontal transfer and the widespread presence of Galileo transposons in Drosophilidae (Insecta: Diptera)

Galileo is a transposon notoriously involved with inversions in Drosophila buzzatii by ectopic recombination. Although widespread in Drosophila, little is known about this transposon in other lineages of Drosophilidae. Here, the abundance of the canonical Galileo and its evolutionary history in Drosophilidae genomes was estimated and reconstructed across genera within its two subfamilies. Sequences of this transposon were masked in these genomes and their transposase sequences were recovered using BLASTn. Phylogenetic analyses were employed to reconstruct their evolutionary history and compare it to that of host genomes. Galileo was found in nearly all 163 species, however, only 37 harbored nearly complete transposase sequences. In the remaining, Galileo was found highly fragmented. Copies from related species were clustered, however horizontal transfer events were detected between the melanogaster and montium groups of Drosophila, and between the latter and the Lordiphosa genus. The similarity of sequences found in the virilis and willistoni groups of Drosophila was found to be a consequence of lineage sorting. Therefore, the evolution of Galileo is primarily marked by vertical transmission and long-term inactivation, mainly through the deletion of open reading frames. The latter has the potential to lead copies of this transposon to become miniature inverted-repeat transposable elements.

Adaptation, ancestral variation and gene flow in a 'Sky Island' Drosophila species

Over time, populations of species can expand, contract, fragment and become isolated, creating subpopulations that must adapt to local conditions. Understanding how species maintain variation after divergence as well as adapt to these changes in the face of gene flow is of great interest, especially as the current climate crisis has caused range shifts and frequent migrations for many species. Here, we characterize how a mycophageous fly species, Drosophila innubila, came to inhabit and adapt to its current range which includes mountain forests in south-western USA separated by large expanses of desert. Using population genomic data from more than 300 wild-caught individuals, we examine four populations to determine their population history in these mountain forests, looking for signatures of local adaptation. In this first extensive study, establishing D. innubila as a key genomic "Sky Island" model, we find D. innubila spread northwards during the previous glaciation period (30-100 KYA) and have recently expanded even further (0.2-2 KYA). D. innubila shows little evidence of population structure, consistent with a recent establishment and genetic variation maintained since before geographic stratification. We also find some signatures of recent selective sweeps in chorion proteins and population differentiation in antifungal immune genes suggesting differences in the environments to which flies are adapting. However, we find little support for long-term recurrent selection in these genes. In contrast, we find evidence of long-term recurrent positive selection in immune pathways such as the Toll signalling system and the Toll-regulated antimicrobial peptides.

Structural and functional liaisons between transposable elements and satellite DNAs

Transposable elements (TEs) and satellite DNAs (satDNAs) are typically identified as major repetitive DNA components in eukaryotic genomes. TEs are DNA segments able to move throughout a genome while satDNAs are tandemly repeated sequences organized in long arrays. Both classes of repetitive sequences are extremely diverse, and many TEs and satDNAs exist within a genome. Although they differ in structure, genomic organization, mechanisms of spread, and evolutionary dynamics, TEs and satDNAs can share sequence similarity and organizational patterns, thus indicating that complex mutual relationships can determine their evolution, and ultimately define roles they might have on genome architecture and function. Motivated by accumulating data about sequence elements that incorporate features of both TEs and satDNAs, here we present an overview of their structural and functional liaisons.

Drosophila muller f elements maintain a distinct set of genomic properties over 40 million years of evolution

The Muller F element (4.2 Mb, ~80 protein-coding genes) is an unusual autosome of Drosophila melanogaster; it is mostly heterochromatic with a low recombination rate. To investigate how these properties impact the evolution of repeats and genes, we manually improved the sequence and annotated the genes on the D. erecta, D. mojavensis, and D. grimshawi F elements and euchromatic domains from the Muller D element. We find that F elements have greater transposon density (25–50%) than euchromatic reference regions (3–11%). Among the F elements, D. grimshawi has the lowest transposon density (particularly DINE-1: 2% vs. 11–27%). F element genes have larger coding spans, more coding exons, larger introns, and lower codon bias. Comparison of the Effective Number of Codons with the Codon Adaptation Index shows that, in contrast to the other species, codon bias in D. grimshawi F element genes can be attributed primarily to selection instead of mutational biases, suggesting that density and types of transposons affect the degree of local heterochromatin formation. F element genes have lower estimated DNA melting temperatures than D element genes, potentially facilitating transcription through heterochromatin. Most F element genes (~90%) have remained on that element, but the F element has smaller syntenic blocks than genome averages (3.4–3.6 vs. 8.4–8.8 genes per block), indicating greater rates of inversion despite lower rates of recombination. Overall, the F element has maintained characteristics that are distinct from other autosomes in the Drosophila lineage, illuminating the constraints imposed by a heterochromatic milieu.

Structural and sequence diversity of the transposon Galileo in the Drosophila willistoni genome

Background

Galileo is one of three members of the P superfamily of DNA transposons. It was originally discovered in Drosophila buzzatii, in which three segregating chromosomal inversions were shown to have been generated by ectopic recombination between Galileo copies. Subsequently, Galileo was identified in six of 12 sequenced Drosophila genomes, indicating its widespread distribution within this genus. Galileo is strikingly abundant in Drosophila willistoni, a neotropical species that is highly polymorphic for chromosomal inversions, suggesting a role for this transposon in the evolution of its genome.

Results

We carried out a detailed characterization of all Galileo copies present in the D. willistoni genome. A total of 191 copies, including 133 with two terminal inverted repeats (TIRs), were classified according to structure in six groups. The TIRs exhibited remarkable variation in their length and structure compared to the most complete copy. Three copies showed extended TIRs due to internal tandem repeats, the insertion of other transposable elements (TEs), or the incorporation of non-TIR sequences into the TIRs. Phylogenetic analyses of the transposase (TPase)-encoding and TIR segments yielded two divergent clades, which we termed Galileo subfamilies V and W. Target-site duplications (TSDs) in D. willistoni Galileo copies were 7- or 8-bp in length, with the consensus sequence GTATTAC. Analysis of the region around the TSDs revealed a target site motif (TSM) with a 15-bp palindrome that may give rise to a stem-loop secondary structure.

Conclusions

There is a remarkable abundance and diversity of Galileo copies in the D. willistoni genome, although no functional copies were found. The TIRs in particular have a dynamic structure and extend in different ways, but their ends (required for transposition) are more conserved than the rest of the element. The D. willistoni genome harbors two Galileo subfamilies (V and W) that diverged ~9 million years ago and may have descended from an ancestral element in the genome. Galileo shows a significant insertion preference for a 15-bp palindromic TSM.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-792) contains supplementary material, which is available to authorized users.

Tetris is a foldback transposon that provided the building blocks for an emerging satellite DNA of Drosophila virilis

Transposable elements (TEs) and satellite DNAs (satDNAs) are abundant components of most eukaryotic genomes studied so far and their impact on evolution has been the focus of several studies. A number of studies linked TEs with satDNAs, but the nature of their evolutionary relationships remains unclear. During in silico analyses of the Drosophila virilis assembled genome, we found a novel DNA transposon we named Tetris based on its modular structure and diversity of rearranged forms. We aimed to characterize Tetris and investigate its role in generating satDNAs. Data mining and sequence analysis showed that Tetris is apparently nonautonomous, with a structure similar to foldback elements, and present in D. virilis and D. americana. Herein, we show that Tetris shares the final portions of its terminal inverted repeats (TIRs) with DAIBAM, a previously described miniature inverted transposable element implicated in the generation of chromosome inversions. Both elements are likely to be mobilized by the same autonomous TE. Tetris TIRs contain approximately 220-bp internal tandem repeats that we have named TIR-220. We also found TIR-220 repeats making up longer (kb-size) satDNA-like arrays. Using bioinformatic, phylogenetic and cytogenomic tools, we demonstrated that Tetris has contributed to shaping the genomes of D. virilis and D. americana, providing internal tandem repeats that served as building blocks for the amplification of satDNA arrays. The β-heterochromatic genomic environment seemed to have favored such amplification. Our results imply for the first time a role for foldback elements in generating satDNAs.

A genome-wide survey of genetic instability by transposition in Drosophila hybrids

Hybridization between species is a genomic instability factor involved in increasing mutation rate and new chromosomal rearrangements. Evidence of a relationship between interspecific hybridization and transposable element mobilization has been reported in different organisms, but most studies are usually performed with particular TEs and do not discuss the real effect of hybridization on the whole genome. We have therefore studied whole genome instability of Drosophila interspecific hybrids, looking for the presence of new AFLP markers in hybrids. A high percentage (27–90%) of the instability markers detected corresponds to TEs belonging to classes I and II. Moreover, three transposable elements (Osvaldo, Helena and Galileo) representative of different families, showed an overall increase of transposition rate in hybrids compared to parental species. This research confirms the hypothesis that hybridization induces genomic instability by transposition bursts and suggests that genomic stress by transposition could contribute to a relaxation of mechanisms controlling TEs in the Drosophila genome.

A divergent P element and its associated MITE, BuT5, generate chromosomal inversions and are widespread within the Drosophila repleta species group

The transposon BuT5 caused two chromosomal inversions fixed in two Drosophila species of the repleta group, D. mojavensis and D. uniseta. BuT5 copies are approximately 1-kb long, lack any coding capacity, and do not resemble any other transposable element (TE). Because of its elusive features, BuT5 has remained unclassified to date. To fully characterize BuT5, we carried out bioinformatic similarity searches in available sequenced genomes, including 21 Drosophila species. Significant hits were only recovered for D. mojavensis genome, where 48 copies were retrieved, 22 of them approximately 1-kb long. Polymerase chain reaction (PCR) and dot blot analyses on 54 Drosophila species showed that BuT5 is homogeneous in size and has a widespread distribution within the repleta group. Thus, BuT5 can be considered as a miniature inverted-repeat TE. A detailed analysis of the BuT5 hits in D. mojavensis revealed three partial copies of a transposon with ends very similar to BuT5 and a P-element-like transposase-encoding region in between. A putatively autonomous copy of this P element was isolated by PCR from D. buzzatii. This copy is 3,386-bp long and possesses a seven-exon gene coding for an 822-aa transposase. Exon–intron boundaries were confirmed by reverse transcriptase-PCR experiments. A phylogenetic tree built with insect P superfamily transposases showed that the D. buzzatii P element belongs to an early diverging lineage within the P-element family. This divergent P element is likely the master transposon mobilizing BuT5. The BuT5/P element partnership probably dates back approximately 16 Ma and is the ultimate responsible for the generation of the two chromosomal inversions in the Drosophila repleta species group.

FB-NOF is a non-autonomous transposable element, expressed in Drosophila melanogaster and present only in the melanogaster group

Most foldback elements are defective due to the lack of coding sequences but some are associated with coding sequences and may represent the entire element. This is the case of the NOF sequences found in the FB of Drosophila melanogaster, formerly considered as an autonomous TE and currently proposed as part of the so-called FB-NOF element, the transposon that would be complete and fully functional. NOF is always associated with FB and never seen apart from the FB inverted repeats (IR). This is the reason why the FB-NOF composite element can be considered the complete element. At least one of its ORFs encodes a protein that has always been considered its transposase, but no detailed studies have been carried out to verify this. In this work we test the hypothesis that FB-NOF is an active transposon nowadays. We search for its expression product, obtaining its cDNA, and propose the ORF and the sequence of its potential protein. We found that the NOF protein is not a transposase as it lacks any of the motifs of known transposases and also shows structural homology with hydrolases, therefore FB-NOF cannot belong to the superfamily MuDR/foldback, as up to now it has been classified, and can be considered as a non-autonomous transposable element. The alignment with the published genomes of 12 Drosophila species shows that NOF presence is restricted only to the 6 Drosophila species belonging to the melanogaster group.

Identification of multiple binding sites for the THAP domain of the Galileo transposase in the long terminal inverted-repeats

Galileo is a DNA transposon responsible for the generation of several chromosomal inversions in Drosophila. In contrast to other members of the P-element superfamily, it has unusually long terminal inverted-repeats (TIRs) that resemble those of Foldback elements. To investigate the function of the long TIRs we derived consensus and ancestral sequences for the Galileo transposase in three species of Drosophilids. Following gene synthesis, we expressed and purified their constituent THAP domains and tested their binding activity towards the respective Galileo TIRs. DNase I footprinting located the most proximal DNA binding site about 70 bp from the transposon end. Using this sequence we identified further binding sites in the tandem repeats that are found within the long TIRs. This suggests that the synaptic complex between Galileo ends may be a complicated structure containing higher-order multimers of the transposase. We also attempted to reconstitute Galileo transposition in Drosophila embryos but no events were detected. Thus, although the limited numbers of Galileo copies in each genome were sufficient to provide functional consensus sequences for the THAP domains, they do not specify a fully active transposase. Since the THAP recognition sequence is short, and will occur many times in a large genome, it seems likely that the multiple binding sites within the long, internally repetitive, TIRs of Galileo and other Foldback-like elements may provide the transposase with its binding specificity.

Striking structural dynamism and nucleotide sequence variation of the transposon Galileo in the genome of Drosophila mojavensis

Background

Galileo is a transposable element responsible for the generation of three chromosomal inversions in natural populations of Drosophila buzzatii. Although the most characteristic feature of Galileo is the long internally-repetitive terminal inverted repeats (TIRs), which resemble the Drosophila Foldback element, its transposase-coding sequence has led to its classification as a member of the P-element superfamily (Class II, subclass 1, TIR order). Furthermore, Galileo has a wide distribution in the genus Drosophila, since it has been found in 6 of the 12 Drosophila sequenced genomes. Among these species, D. mojavensis, the one closest to D. buzzatii, presented the highest diversity in sequence and structure of Galileo elements.

Results

In the present work, we carried out a thorough search and annotation of all the Galileo copies present in the D. mojavensis sequenced genome. In our set of 170 Galileo copies we have detected 5 Galileo subfamilies (C, D, E, F, and X) with different structures ranging from nearly complete, to only 2 TIR or solo TIR copies. Finally, we have explored the structural and length variation of the Galileo copies that point out the relatively frequent rearrangements within and between Galileo elements. Different mechanisms responsible for these rearrangements are discussed.

Conclusions

Although Galileo is a transposable element with an ancient history in the D. mojavensis genome, our data indicate a recent transpositional activity. Furthermore, the dynamism in sequence and structure, mainly affecting the TIRs, suggests an active exchange of sequences among the copies. This exchange could lead to new subfamilies of the transposon, which could be crucial for the long-term survival of the element in the genome.

msechBari, a new MITE-like element in Drosophila sechellia related to the Bari transposon

A few occurrences of miniature inverted-repeat transposable elements (MITEs) have been reported in species of the genus Drosophila. Here, we describe msechBari, a MITE-like element in Drosophila sechellia. The element is short, approximately 90 bp in length, AT-rich and occurs in association with, or close to, genes, characteristics that are typical for MITEs. The identification was performed in silico using the sequenced genome of D. sechellia and confirmed in a laboratory strain. This short element is related to the Bari_DM transposon of Drosophila melanogaster, having terminal inverted repeats (TIRs) of a similar length and a high identity with the full-length Bari_DM element. The estimated recent origin of the element and the homogeneity observed between copies found in the genome suggests that msechBari could be active in D. sechellia.

Gene alterations at Drosophila inversion breakpoints provide prima facie evidence for natural selection as an explanation for rapid chromosomal evolution

Background

Chromosomal inversions have been pervasive during the evolution of the genus Drosophila, but there is significant variation between lineages in the rate of rearrangement fixation. D. mojavensis, an ecological specialist adapted to a cactophilic niche under extreme desert conditions, is a chromosomally derived species with ten fixed inversions, five of them not present in any other species.

Results

In order to explore the causes of the rapid chromosomal evolution in D. mojavensis, we identified and characterized all breakpoints of seven inversions fixed in chromosome 2, the most dynamic one. One of the inversions presents unequivocal evidence for its generation by ectopic recombination between transposon copies and another two harbor inverted duplications of non-repetitive DNA at the two breakpoints and were likely generated by staggered single-strand breaks and repair by non-homologous end joining. Four out of 14 breakpoints lay in the intergenic region between preexisting duplicated genes, suggesting an adaptive advantage of separating previously tightly linked duplicates. Four out of 14 breakpoints are associated with transposed genes, suggesting these breakpoints are fragile regions. Finally two inversions contain novel genes at their breakpoints and another three show alterations of genes at breakpoints with potential adaptive significance.

Conclusions

D. mojavensis chromosomal inversions were generated by multiple mechanisms, an observation that does not provide support for increased mutation rate as explanation for rapid chromosomal evolution. On the other hand, we have found a number of gene alterations at the breakpoints with putative adaptive consequences that directly point to natural selection as the cause of D. mojavensis rapid chromosomal evolution.

Tm1: a mutator/foldback transposable element family in root-knot nematodes

Three closely related parthenogenetic species of root-knot nematodes, collectively termed the Meloidogyne incognita-group, are economically significant pathogens of diverse crop species. Remarkably, these asexual root-knot nematodes are capable of acquiring heritable changes in virulence even though they lack sexual reproduction and meiotic recombination. Characterization of a near isogenic pair of M. javanica strains differing in response to tomato with the nematode resistance gene Mi-1 showed that the virulent strain carried a deletion spanning a gene called Cg-1. Herein, we present evidence that the Cg-1 gene lies within a member of a novel transposable element family (Tm1; Transposon in Meloidogyne-1). This element family is defined by composite terminal inverted repeats of variable lengths similar to those of Foldback (FB) transposable elements and by 9 bp target site duplications. In M. incognita, Tm1 elements can be classified into three general groups: 1) histone-hairpin motif elements; 2) MITE-like elements; 3) elements encoding a putative transposase. The predicted transposase shows highest similarity to gene products encoded by aphids and mosquitoes and resembles those of the Phantom subclass of the Mutator transposon superfamily. Interestingly, the meiotic, sexually-reproducing root-knot nematode species M. hapla has Tm1 elements with similar inverted repeat termini, but lacks elements with histone hairpin motifs and contains no elements encoding an intact transposase. These Tm1 elements may have impacts on root-knot nematode genomes and contribute to genetic diversity of the asexual species.

Genetic variation of copia suppression in Drosophila melanogaster

Transposable elements (TEs) are genomic parasites that propagate by exploiting its host reproductive machinery. However, some hosts have evolved the ability to silence TE activity, whereas others have not. We are investigating the population dynamics of TE host-silencing pathways, particularly copia long terminal repeat retrotransposon in Drosophila melanogaster. Here, we identify large effect genes involved in copia suppression by using a semi-quantitative analysis to assay levels of copia plasmids (believed to be an intermediate of transposition) in 98 recombinant inbred lines constructed from a line exhibiting high copia transpositions and a line exhibiting no transpositions. The results revealed that the influence of copia copy number and transcription level on copia plasmid concentrations are weak and that genomic factors, presumably encoded by the host, have stronger effects on transposition rates. We mapped a QTL affecting copia plasmid concentration within the 33A-43E cytological region of the second chromosome and applied a quantitative deficiency complementation analysis on this chromosomal region. One out of the two large effect deficiencies on copia plasmid concentrations corresponded to the vasa gene, an important component of the nuage-piwi RNA TE-silencing machinery. We hypothesize that copia suppression occurs by the joint action of several post-transcriptional mechanisms with at least one of the blocks taking place in the nuage.

Phantom, a new subclass of Mutator DNA transposons found in insect viruses and widely distributed in animals

Transposons of the Mutator (Mu) superfamily have been shown to play a critical role in the evolution of plant genomes. However, the identification of Mutator transposons in other eukaryotes has been quite limited. Here we describe a previously uncharacterized group of DNA transposons designated Phantom identified in the genomes of a wide range of eukaryotic taxa, including many animals, and provide evidence for its inclusion within the Mutator superfamily. Interestingly three Phantom proteins were also identified in two insect viruses and phylogenetic analysis suggests horizontal movement from insect to virus, providing a new line of evidence for the role of viruses in the horizontal transfer of DNA transposons in animals. Many of the Phantom transposases are predicted to harbor a FLYWCH domain in the amino terminus, which displays a WRKY-GCM1 fold characteristic of the DNA binding domain (DBD) of Mutator transposases and of several transcription factors. While some Phantom elements have terminal inverted repeats similar in length and structure to Mutator elements, some display subterminal inverted repeats (sub-TIRs) and others have more complex termini reminiscent of so-called Foldback (FB) transposons. The structural plasticity of Phantom and the distant relationship of its encoded protein to known transposases may have impeded the discovery of this group of transposons and it suggests that structure in itself is not a reliable character for transposon classification.

Comparative Genomics on the Drosophila Phylogenetic Tree

With the sequencing of 12 complete euchromatic Drosophila genomes, the genus Drosophila is a leading model for comparative genomics. In this review, we discuss the novel insights into evolutionary processes afforded by the newly available genomic sequences when placed in the context of the phylogeny. We focus on three levels: insights into whole-genome content, such as changes in genome size and content across the phylogeny; insights into large-scale patterns of divergence and conservation, such as selective constraints on genes and chromosome-level evolution of sex chromosomes; and insights into finer-scale processes in individual lineages and genes, such as lineage-specific evolution in response to ecological context. As the field of comparative genomics is still young, we also discuss current challenges, such as the development of more sophisticated evolutionary models to capture nonequilibrium processes and the improvement of assembly and alignment algorithms to better capture uncertainty in the data.

The transposon Galileo generates natural chromosomal inversions in Drosophila by ectopic recombination

Background

Transposable elements (TEs) are responsible for the generation of chromosomal inversions in several groups of organisms. However, in Drosophila and other Dipterans, where inversions are abundant both as intraspecific polymorphisms and interspecific fixed differences, the evidence for a role of TEs is scarce. Previous work revealed that the transposon Galileo was involved in the generation of two polymorphic inversions of Drosophila buzzatii.

Methodology/Principal Findings

To assess the impact of TEs in Drosophila chromosomal evolution and shed light on the mechanism involved, we isolated and sequenced the two breakpoints of another widespread polymorphic inversion from D. buzzatii, 2z³. In the non inverted chromosome, the 2z³ distal breakpoint was located between genes CG2046 and CG10326 whereas the proximal breakpoint lies between two novel genes that we have named Dlh and Mdp. In the inverted chromosome, the analysis of the breakpoint sequences revealed relatively large insertions (2,870-bp and 4,786-bp long) including two copies of the transposon Galileo (subfamily Newton), one at each breakpoint, plus several other TEs. The two Galileo copies: (i) are inserted in opposite orientation; (ii) present exchanged target site duplications; and (iii) are both chimeric.

Conclusions/Significance

Our observations provide the best evidence gathered so far for the role of TEs in the generation of Drosophila inversions. In addition, they show unequivocally that ectopic recombination is the causative mechanism. The fact that the three polymorphic D. buzzatii inversions investigated so far were generated by the same transposon family is remarkable and is conceivably due to Galileo's unusual structure and current (or recent) transpositional activity.

Evolution of Tom, 297, 17.6 and rover retrotransposons in Drosophilidae species

LTR retrotransposons are the most abundant transposable elements in Drosophila and are believed to have contributed significantly to genome evolution. Different reports have shown that many LTR retrotransposon families in Drosophila melanogaster emerged from recent evolutionary episodes of transpositional activity. To contribute to the knowledge of the evolutionary history of Drosophila LTR retrotransposons and the mechanisms that control their abundance, distribution and diversity, we conducted analyses of four related families of LTR retrotransposons, 297, 17.6, rover and Tom. Our results show that these elements seem to be restricted to species from the D. melanogaster group, except for 17.6, which is also present in D. virilis and D. mojavensis. Genetic divergences and phylogenetic analyses of a 1-kb fragment region of the pol gene illustrate that the evolutionary dynamics of Tom, 297, 17.6 and rover retrotransposons are similar in several aspects, such as low codon bias, the action of purifying selection and phylogenies that are incongruent with those of the host species. We found an extremely complex association among the retrotransposon sequences, indicating that different processes shaped the evolutionary history of these elements, and we detected a very high number of possible horizontal transfer events, corroborating the importance of lateral transmission in the evolution and maintenance of LTR retrotransposons.

Cloning and sequencing of the breakpoint regions of inversion 5g fixed in Drosophila buzzatii

Chromosomal inversions are ubiquitous in Drosophila both as intraspecific polymorphisms and interspecific differences. Many gaps still remain in our understanding of the mechanisms that generate them. Previous work has shown that in Drosophila buzzatii, three polymorphic inversions were generated by ectopic recombination between copies of the transposon Galileo. In this study, we have characterized the breakpoint regions of inversion 5g, fixed in D. buzzatii and absent in Drosophila koepferae and other closely related species. A novel approach comprising four experimental steps was used. First, D. buzzatii BAC clones encompassing the breakpoints were identified and their ends sequenced. Then, breakpoint regions were mapped at high resolution in the Drosophila mojavensis genome sequence. Finally, breakpoint regions were isolated by polymerase chain reaction in D. buzzatii and D. koepferae and sequenced. Our aim was to shed light on the mechanism that generated inversion 5g and specifically to test for an implication of the transposon Galileo. No evidence implicates Galileo or other transposable elements in the origin of inversion 5g that was generated most likely by two independent breaks and non-homologous end-joining repair. Our results show that different inversion-generating mechanisms may coexist within the same lineage and suggest a hypothesis for the evolutionary time and mode of their operation.