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.     

Amplification of the 1731 LTR retrotransposon in Drosophila melanogaster cultured cells: origin of neocopies and impact on the genome

Transposable elements (TEs), represent a large fraction of the eukaryotic genome. In Drosophila melanogaster, about 20% of the genome corresponds to such middle repetitive DNA dispersed sequences. A fraction of TEs is composed of elements showing a retrovirus-like structure, the LTR-retrotransposons, the first TEs to be described in the Drosophila genome. Interestingly, in D. melanogaster embryonic immortal cell culture genomes the copy number of these LTR-retrotransposons was revealed to be higher than the copy number in the Drosophila genome, presumably as the result of transposition of some copies to new genomic locations [Potter, S.S., Brorein Jr., W.J., Dunsmuir, P., Rubin, G.M., 1979. Transposition of elements of the 412, copia and 297 dispersed repeated gene families in Drosophila. Cell 17, 415-427; Junakovic, N., Di Franco, C., Best-Belpomme, M., Echalier, G., 1988. On the transposition of copia-like nomadic elements in cultured Drosophila cells. Chromosoma 97, 212-218]. This suggests that so many transpositions modified the genome organisation and consequently the expression of targeted genes. To understand what has directed the transposition of TEs in Drosophila cell culture genomes, a search to identify the newly transposed copies was undertaken using 1731, a LTR-retrotransposon. A comparison between 1731 full-length elements found in the fly sequenced genome (y(1); cn(1)bw(1), sp(1) stock) and 1731 full-length elements amplified by PCR in the two cell line was done. The resulting data provide evidence that all 1731 neocopies were derived from a single copy slightly active in the Drosophila genome and subsequently strongly activated in cultured cells; and that this active copy is related to a newly evolved genomic variant (Kalmykova, A.I., et al., 2004. Selective expansion of the newly evolved genomic variants of retrotransposon 1731 in the Drosophila genomes. Mol. Biol. Evol. 21, 2281-2289). Moreover, neocopies are shown to be inserted in different sets of genes in the two cell lines suggesting they might be involved in the biological and physiological differences observed between Kc and S2 cell lines.     

Identification of Two Subfamilies of micropia Transposable Element in Species of the repleta Group of Drosophila

The occurrence, number of insertion sites and antisense RNA expression of micropia transposable element were studied in 26 species that belong to three subgroups (mercatorum, mulleri and hydei) of repleta group of Drosophila. Under high specific PCR, micropia sequences were detected in 11 species, but under less stringent condition, this retrotransposon was detected in all species. The widespread distribution of micropia suggests that this element was already present at the common ancestor of the repleta group of Drosophila. Southern blot analysis showed a variation from 0 to 17 different insertion sites and the occurrence of male-specific sequences. We found that the expression of the 1.0 kb micropia antisense RNA is variable among the species and tissues (soma and testis), which suggests that more than one mechanism regulates transposition in these species. Variation of amplification by PCR and of antisense RNA expression, as well as divergence of nucleotide sequences among the species allow us to suggest that at least two subfamilies of micropia transposable element are harbored by the genome of this species group.     

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.     

Evolutionary patterns of the gypsy and bilbo retrotransposon families in the Drosophila species of the obscura group

We analyse in this paper the evolutionary patterns of two types of Drosophila retrotransposons, gypsy (a virus-like element), and bilbo (a LINE-like element), in host species from the Drosophila and Scaptomyza genus. Phylogenetic analysis of the retrotransposon sequences amplified by PCR, revealed concordance with the phylogeny of the Drosophila host species from the obscura group, which is consistent with vertical transmission during differentiation of the species. However, in the species outside of the obscura group, horizontal transmission can be considered. The amplified sequences that presented intact open reading frames were used in an analysis of the evolutionary constraints on the amino acid sequences. The analysed sequences seem to be functional, and the selective constraints are evidenced, especially when sequences from distant species are compared. Comparison of the evolutionary rates of both retrotransposons in the same species, suggests that bilbo seems to evolve more rapidly than gypsy.     

Differential detection of transposable elements between Saccharum species

Cultivars of sugarcane (Saccharum) are hybrids between species S. officinarum (x = 10, 2n = 8x = 80) and S. spontaneum (x = 8, 2n = 5 – 16x = 40 – 128). These accessions have 100 to 130 chromosomes, 80–85% of which are derived from S. officinarum, 10–15% from S. spontaneum, and 5–10% are possible recombinants between the two genomes. The aim of this study was to analyze the repetition of DNA sequences in S. officinarum and S. spontaneum. For this purpose, genomic DNA from S. officinarum was digested with restriction enzymes and the fragments cloned. Sixty-eight fragments, approximately 500 bp, were cloned, sequenced and had their identity analyzed in NCBI, and in the rice, maize, and sorghum genome databases using BLAST. Twelve clones containing partial transposable elements, one single-copy control, one DNA repetitive clone control and two genome controls were analyzed by DNA hybridization on membrane, using genomic probes from S. officinarum and S. spontaneum. The hybridization experiment revealed that six TEs had a similar repetitive DNA pattern in the genomes of S. officinarum and S. spontaneum, while six TEs were more abundant in the genome of S. officinarum. We concluded that the species S. officinarum and S. spontaneum have differential accumulation LTR retrotransposon families, suggesting distinct insertion or modification patterns.     

An Age-of-Allele Test of Neutrality for Transposable Element Insertions

How natural selection acts to limit the proliferation of transposable elements (TEs) in genomes has been of interest to evolutionary biologists for many years. To describe TE dynamics in populations, previous studies have used models of transposition–selection equilibrium that assume a constant rate of transposition. However, since TE invasions are known to happen in bursts through time, this assumption may not be reasonable. Here we propose a test of neutrality for TE insertions that does not rely on the assumption of a constant transposition rate. We consider the case of TE insertions that have been ascertained from a single haploid reference genome sequence. By conditioning on the age of an individual TE insertion allele (inferred by the number of unique substitutions that have occurred within the particular TE sequence since insertion), we determine the probability distribution of the insertion allele frequency in a population sample under neutrality. Taking models of varying population size into account, we then evaluate predictions of our model against allele frequency data from 190 retrotransposon insertions sampled from North American and African populations of Drosophila melanogaster. Using this nonequilibrium neutral model, we are able to explain ∼80% of the variance in TE insertion allele frequencies based on age alone. Controlling for both nonequilibrium dynamics of transposition and host demography, we provide evidence for negative selection acting against most TEs as well as for positive selection acting on a small subset of TEs. Our work establishes a new framework for the analysis of the evolutionary forces governing large insertion mutations like TEs, gene duplications, or other copy number variants.     

What makes transposable elements move in the Drosophila genome?

Transposable elements (TEs), by their capacity of moving and inducing mutations in the genome, are considered important drivers of species evolution. The successful invasions of TEs in genomes, despite their mutational properties, are an apparent paradox. TEs' transposition is usually strongly regulated to low value, but in some cases these elements can also show high transposition rates, which has been associated sometimes to changes in environmental conditions. It is evident that factors susceptible to induce transpositions in natural populations contribute to TE perpetuation. Different factors were proposed as causative agents of TE mobilization in a wide range of organisms: biotic and abiotic stresses, inter- and intraspecific crosses and populational factors. However, there is no clear evidence of the factors capable of inducing TE mobilization in Drosophila, and data on laboratory stocks show contradictory results. The aim of this review is to have an update critical revision about mechanisms promoting transposition of TEs in Drosophila, and to provide to the readers a global vision of the dynamics of these genomic elements in the Drosophila genome.     

Single-Nucleotide-Specific Targeting of the Tf1 Retrotransposon Promoted by the DNA-Binding Protein Sap1 of Schizosaccharomyces pombe

Transposable elements (TEs) constitute a substantial fraction of the eukaryotic genome and, as a result, have a complex relationship with their host that is both adversarial and dependent. To minimize damage to cellular genes, TEs possess mechanisms that target integration to sequences of low importance. However, the retrotransposon Tf1 of Schizosaccharomyces pombe integrates with a surprising bias for promoter sequences of stress-response genes. The clustering of integration in specific promoters suggests that Tf1 possesses a targeting mechanism that is important for evolutionary adaptation to changes in environment. We report here that Sap1, an essential DNA-binding protein, plays an important role in Tf1 integration. A mutation in Sap1 resulted in a 10-fold drop in Tf1 transposition, and measures of transposon intermediates support the argument that the defect occurred in the process of integration. Published ChIP-Seq data on Sap1 binding combined with high-density maps of Tf1 integration that measure independent insertions at single-nucleotide positions show that 73.4% of all integration occurs at genomic sequences bound by Sap1. This represents high selectivity because Sap1 binds just 6.8% of the genome. A genome-wide analysis of promoter sequences revealed that Sap1 binding and amounts of integration correlate strongly. More important, an alignment of the DNA-binding motif of Sap1 revealed integration clustered on both sides of the motif and showed high levels specifically at positions +19 and −9. These data indicate that Sap1 contributes to the efficiency and position of Tf1 integration.     

A non-LTR retroelement extinction in Spermophilus tridecemlineatus

The typical mammalian genome is dominated by two types of transposable elements (TEs), the autonomous and non-autonomous non-LTR retrotransposons, i.e. LINEs and SINEs, and with few exceptions there is a sole active LINE family (L1). During an ongoing investigation of TEs in rodents we determined that overall transposon activity has been steadily declining in Spermophilus tridecemlineatus. More specifically, the typically ubiquitous L1 activity of mammals has decreased drastically within the last 26MY. Indeed, only three L1 insertions with intact ORF1 sequences were readily identifiable and no intact ORF2 sequences were identified. The last L1 and SINE insertions date to ~5.3MYA and 4MYA, respectively. Based on our inability to computationally identify recently inserted L1 elements we suggest that S. tridecemlineatus is experiencing a quiescence or extinction of non-LTR retrotransposon activity. Such a finding represents only the fourth instance of a loss of non-LTR retrotransposon activity identified in mammals and, as such, represents an important additional data point to guide our understanding of LINE dynamics in eutherians.     

Mobilizing the genome of Lepidoptera through novel sequence gains and end creation by non-autonomous Lep1 Helitrons

Transposable elements (TEs) can affect the structure of genomes through their acquisition and transposition of novel DNA sequences. The 134-bp repetitive elements, Lep1, are conserved non-autonomous Helitrons in lepidopteran genomes that have characteristic 5'-CT and 3'-CTAY nucleotide termini, a 3'-terminal hairpin structure, a 5'- and 3'-subterminal inverted repeat (SIR), and integrations that occur between AT or TT nucleotides. Lep1 Helitrons have acquired and propagated sequences downstream of their 3'-CTAY termini that are 57-344-bp in length and have termini composed of a 3'-CTRR preceded by a 3'-hairpin structure and a region complementary to the 5'-SIR (3'-SIRb). Features of both the Lep1 Helitron and multiple acquired sequences indicate that secondary structures at the 3'-terminus may have a role in rolling circle replication or genome integration mechanisms, and are a prerequisite for novel end creation by Helitron-like TEs. The preferential integration of Lep1 Helitrons in proximity to gene-coding regions results in the creation of genetic novelty that is shown to impact gene structure and function through the introduction of novel exon sequence (exon shuffling). These findings are important in understanding the structural requirements of genomic DNA sequences that are acquired and transposed by Helitron-like TEs.