Penelope retroelements from Drosophila virilis are active after transformation of Drosophila melanogaster

Experimental evolution reveals hyperparasitic interactions among transposable elements

Transposable elements (TEs) are repeated DNA sequences that can constitute a substantial part of genomes. Studying TEs' activity, interactions, and accumulation dynamics is thus of major interest to understand genome evolution. Here, we describe the transposition dynamics of cut-and-paste mariner elements during experimental (short- and longer-term) evolution in Drosophila melanogaster Flies with autonomous and nonautonomous mariner copies were introduced in populations containing no active mariner, and TE accumulation was tracked by quantitative PCR for up to 100 generations. Our results demonstrate that (i) active mariner elements are highly invasive and characterized by an elevated transposition rate, confirming their capacity to spread in populations, as predicted by the "selfish-DNA" mechanism; (ii) nonautonomous copies act as parasites of autonomous mariner elements by hijacking the transposition machinery produced by active mariner, which can be considered as a case of hyperparasitism; (iii) this behavior resulted in a failure of active copies to amplify which systematically drove the whole family to extinction in less than 100 generations. This study nicely illustrates how the presence of transposition-competitive variants can deeply impair TE dynamics and gives clues to the extraordinary diversity of TE evolutionary histories observed in genomes.

What happens when Penelope comes?: An unusual retroelement invades a host species genome exploring different strategies

Transposable elements (TEs) are ubiquitous residents in eukaryotic genomes. They can cause dramatic changes in gene expression and lead to gross rearrangements of chromosome structure, providing the basis for rapid evolution. The virilis species group of Drosophila contains certain species that can be crossed under experimental conditions and their phylogeny is thoroughly investigated. We have shown that Drosophila virilis, the most primitive karyotypically and probably the ancestral species of the group, is in the process of colonization by a very unusual retroelement Penelope which apparently repeatedly invaded the species of the group in the past. However, the molecular mechanisms and evolutionary consequences of such invasions are poorly understood. In this commentary, we discuss the implications of our recent investigation into the response of the RNA silencing system to Penelope invasion of a new host genome which can be achieved in different ways.

Whole genome sequencing in Drosophila virilis identifies Polyphemus, a recently activated Tc1-like transposon with a possible role in hybrid dysgenesis

BACKGROUND

Hybrid dysgenic syndromes in Drosophila have been critical for characterizing host mechanisms of transposable element (TE) regulation. This is because a common feature of hybrid dysgenesis is germline TE mobilization that occurs when paternally inherited TEs are not matched with a maternal pool of silencing RNAs that maintain transgenerational TE control. In the face of this imbalance TEs become activated in the germline and can cause F1 sterility. The syndrome of hybrid dysgenesis in Drosophila virilis was the first to show that the mobilization of one dominant TE, the Penelope retrotransposon, may lead to the mobilization of other unrelated elements. However, it is not known how many different elements contribute and no exhaustive search has been performed to identify additional ones. To identify additional TEs that may contribute to hybrid dysgenesis in Drosophila virilis, I analyzed repeat content in genome sequences of inducer and non-inducer lines.

RESULTS

Here I describe Polyphemus, a novel Tc1-like DNA transposon, which is abundant in the inducer strain of D. virilis but highly degraded in the non-inducer strain. Polyphemus expression is also increased in the germline of progeny of the dysgenic cross relative to reciprocal progeny. Interestingly, like the Penelope element, it has experienced recent re-activation within the D. virilis lineage.

CONCLUSIONS

Here I present the results of a comprehensive search to identify additional factors that may cause hybrid dysgenesis in D. virilis. Polyphemus, a novel Tc1-like DNA transposon, has recently become re-activated in Drosophila virilis and likely contributes to the hybrid dysgenesis syndrome. It has been previously shown that the Penelope element has also been re-activated in the inducer strain. This suggests that TE co-reactivation within species may synergistically contribute to syndromes of hybrid dysgenesis.

Endonuclease-containing Penelope retrotransposons in the bdelloid rotifer Adineta vaga exhibit unusual structural features and play a role in expansion of host gene families

Background

Penelope-like elements (PLEs) are an enigmatic group of retroelements sharing a common ancestor with telomerase reverse transcriptases. In our previous studies, we identified endonuclease-deficient PLEs that are associated with telomeres in bdelloid rotifers, small freshwater invertebrates best known for their long-term asexuality and high foreign DNA content. Completion of the high-quality draft genome sequence of the bdelloid rotifer Adineta vaga provides us with the opportunity to examine its genomic transposable element (TE) content, as well as TE impact on genome function and evolution.

Results

We performed an exhaustive search of the A. vaga genome assembly, aimed at identification of canonical PLEs combining both the reverse transcriptase (RT) and the GIY-YIG endonuclease (EN) domains. We find that the RT/EN-containing Penelope families co-exist in the A. vaga genome with the EN-deficient RT-containing Athena retroelements. Canonical PLEs are present at very low copy numbers, often as a single-copy, and there is no evidence that they might preferentially co-mobilize EN-deficient PLEs. We also find that Penelope elements can participate in expansion of A. vaga multigene families via trans-action of their enzymatic machinery, as evidenced by identification of intron-containing host genes framed by the Penelope terminal repeats and characteristic target-site duplications generated upon insertion. In addition, we find that Penelope open reading frames (ORFs) in several families have incorporated long stretches of coding sequence several hundred amino acids (aa) in length that are highly enriched in asparagine residues, a phenomenon not observed in other retrotransposons.

Conclusions

Our results show that, despite their low abundance and low transcriptional activity in the A. vaga genome, endonuclease-containing Penelope elements can participate in expansion of host multigene families. We conclude that the terminal repeats represent the cis-acting sequences required for mobilization of the intervening region in trans by the Penelope-encoded enzymatic activities. We also hypothesize that the unusual capture of long N-rich segments by the Penelope ORF occurs as a consequence of peculiarities of its replication mechanism. These findings emphasize the unconventional nature of Penelope retrotransposons, which, in contrast to all other retrotransposon types, are capable of dispersing intron-containing genes, thereby questioning the validity of traditional estimates of gene retrocopies in PLE-containing eukaryotic genomes.

Evolution and dynamics of small RNA response to a retroelement invasion in Drosophila

Although small RNAs efficiently control transposition activity of most transposons in the host genome, such an immune system is not always applicable against a new transposon's invasions. Here, we explored a possibility to introduce potentially mobile copy of the Penelope retroelement previously implicated in hybrid dysgenesis syndrome in Drosophila virilis into the genomes of two distant Drosophila species. The consequences of such introduction were monitored at different phases after experimental colonization as well as in D. virilis species, which is apparently in the process of ongoing Penelope invasion. We investigated the expression of Penelope and biogenesis of Penelope-derived small RNAs in D. virilis and D. melanogaster strains originally lacking active copies of this element after experimental Penelope invasion. These strains were transformed by constructs containing intact Penelope copies. We show that immediately after transformation, which imitates the first stage of retroelement invasion, Penelope undergoes transposition predominantly in somatic tissues, and may produce siRNAs that are apparently unable to completely silence its activity. However, at the later stages of colonization Penelope copies may jump into one of the piRNA-clusters, which results in production of homologous piRNAs that are maternally deposited and can silence euchromatic transcriptionally active copies of Penelope in trans and, hence, prevent further amplification of the invader in the host genome. Intact Penelope copies and different classes of Penelope-derived small RNAs were found in most geographical strains of D. virilis collected throughout the world. Importantly, all strains of this species containing full-length Penelope tested do not produce gonadal sterility in dysgenic crosses and, hence, exhibit neutral cytotype. To understand whether RNA interference mechanism able to target Penelope operates in related species of the virilis group, we correlated the presence of full-length and potentially active Penelope with the occurrence of piRNAs homologous to this transposable element in the ovaries of species comprising the group. It was demonstrated that Penelope-derived piRNAs are present in all virilis group species containing full-length but transcriptionally silent copies of this element that probably represent the remnants of its previous invasions taking place in the course of the virilis species divergent evolution.

Expression of Drosophila virilis retroelements and role of small RNAs in their intrastrain transposition

Transposition of two retroelements (Ulysses and Penelope) mobilized in the course of hybrid dysgenesis in Drosophila virilis has been investigated by in situ hybridization on polytene chromosomes in two D. virilis strains of different cytotypes routinely used to get dysgenic progeny. The analysis has been repeatedly performed over the last two decades, and has revealed transpositions of Penelope in one of the strains, while, in the other strain, the LTR-containing element Ulysses was found to be transpositionally active. The gypsy retroelement, which has been previously shown to be transpositionally inactive in D. virilis strains, was also included in the analysis. Whole mount is situ hybridization with the ovaries revealed different subcellular distribution of the transposable elements transcripts in the strains studied. Ulysses transpositions occur only in the strain where antisense piRNAs homologous to this TE are virtually absent and the ping-pong amplification loop apparently does not take place. On the other hand small RNAs homologous to Penelope found in the other strain, belong predominantly to the siRNA category (21nt), and consist of sense and antisense species observed in approximately equal proportion. The number of Penelope copies in the latter strain has significantly increased during the last decades, probably because Penelope-derived siRNAs are not maternally inherited, while the low level of Penelope-piRNAs, which are faithfully transmitted from mother to the embryo, is not sufficient to silence this element completely. Therefore, we speculate that intrastrain transposition of the three retroelements studied is controlled predominantly at the post-transcriptional level.

Small RNA-based silencing strategies for transposons in the process of invading Drosophila species

Colonization of a host by an active transposon can increase mutation rates or cause sterility, a phenotype termed hybrid dysgenesis. As an example, intercrosses of certain Drosophila virilis strains can produce dysgenic progeny. The Penelope element is present only in a subset of laboratory strains and has been implicated as a causative agent of the dysgenic phenotype. We have also introduced Penelope into Drosophila melanogaster, which are otherwise naive to the element. We have taken advantage of these natural and experimentally induced colonization processes to probe the evolution of small RNA pathways in response to transposon challenge. In both species, Penelope was predominantly targeted by endo-small-interfering RNAs (siRNAs) rather than by piwi-interacting RNAs (piRNAs). Although we do observe correlations between Penelope transcription and dysgenesis, we could not correlate differences in maternally deposited Penelope piRNAs with the sterility of progeny. Instead, we found that strains that produced dysgenic progeny differed in their production of piRNAs from clusters in subtelomeric regions, possibly indicating that changes in the overall piRNA repertoire underlie dysgenesis. Considered together, our data reveal unexpected plasticity in small RNA pathways in germ cells, both in the character of their responses to invading transposons and in the piRNA clusters that define their ability to respond to mobile elements.

Evidence of multiple retrotransposons in two litopenaeid species

Retrotransposons encompass a specific class of mobile genetic elements that are widespread across eukaryotic genomes. The impact of the varied types of retrotransposons on these genomes is just beginning to be deciphered. In a step towards understanding their role in litopenaeid shrimp, we have herein identified nine non-LTR retrotransposons, among which several appear to exist outside the standard defined clades. Two Litopenaeus stylirostris elements were discovered through degenerate PCR amplification using previously defined non-LTR degenerate primers, and through primers designed from a RAPD-derived sequence. A third genomic L. stylirostris element was identified using specific priming from an amplification protocol. These three PCR-derived sequences showed conserved domains of the non-LTR reverse transcriptase gene. In silico searching of genome databases and subsequent contig construction yielded six non-LTR retrotransposons (both genomic and expressed) in the Litopenaeus vannamei genome that also exhibited the highly conserved domains found in our PCR-derived sequences. Phylogenetic placement among representatives from all non-LTR clades showed a possibly novel monophyletic group that included five of our nine sequences. This group, which included elements from both L. stylirostris and L. vannamei, appeared most closely related to the highly active RTE clade. Our remaining four sequences placed in the CR1 and I clades of retrotransposons, with one showing strong similarity to ancient Penelope elements. This research describes three newly discovered retrotransposons in the L. stylirostris genome. Phylogenetic analysis clusters these in a monophyletic grouping with retrotransposons previously described from two closely related species, L. vannamei and Penaeus monodon.

Transcriptome Analysis of ESTs from a Chaetognath Reveals a Deep-Branching Clade of Retrovirus-Like Retrotransposons

Chaetognaths constitute a small marine phylum exhibiting several characteristic which are highly unusual in animal genomes, including two classes of both rRNA and protein ribosomal genes. As in this phylum presence of retrovirus-like elements has never been documented, analysis of a published expressed sequence tag (EST) collection of the chaetognath Spadella cephaloptera has been made. Twelve sequences representing transcript sections of reverse transcriptase domain of active retrotransposons were isolated from~11,000 ESTs. Five of them are originated from Gypsy retrovirus-like elements, whereas the other are transcripts from a Bel-Pao LTR-retrotransposon, a Penelope-like element and LINE retrotransposons. Moreover, a part of a putative integrase has also been found. Phylogenetic analyses suggest a deep-branching clade of the retrovirus-like elements, which is in agreement with the probably Cambrian origin of the phylum. Moreover, retrotransposons have not been found in telomeric-like transcripts which are probably constituted by both vertebrate and arthropod canonical repeats.

Molecular dissection of Penelope transposable element regulatory machinery

Penelope-like elements (PLEs) represent a new class of retroelements identified in more than 80 species belonging to at least 10 animal phyla. Penelope isolated from Drosophila virilis is the only known transpositionally active representative of this class. Although the size and structure of the Penelope major transcript has been previously described in both D. virilis and D. melanogaster transgenic strains, the architecture of the Penelope regulatory region remains unknown. In order to determine the localization of presumptive Penelope promoter and enhancer-like elements, segments of the putative Penelope regulatory region were linked to a CAT reporter gene and introduced into D. melanogaster by P-element-mediated transformation. The results obtained using ELISA to measure CAT expression levels and RNA studies, including RT-PCR, suggest that the active Penelope transposon contains an internal promoter similar to the TATA-less promoters of LINEs. The results also suggest that some of the Penelope regulatory sequences control the preferential expression in the ovaries of the adult flies by enhancing expression in the ovary and reducing expression in the carcass. The possible significance of the intron within Penelope for the function and evolution of PLEs, and the effect of Penelope insertions on adjacent genes, are discussed.

Effect of divergence time and recombination rate on molecular evolution of Drosophila INE-1 transposable elements and other candidates for neutrally evolving sites

Interspecies divergence of orthologous transposable element remnants is often assumed to be simply due to genetic drift of neutral mutations that occurred after the divergence of the species. However, divergence may also be affected by other factors, such as variation in the mutation rate, ancestral polymorphisms, or selection. Here we attempt to determine the impact of these forces on divergence of three classes of sites that are often assumed to be selectively unconstrained (INE-1 TE remnants, sites within short introns, and fourfold degenerate sites) in two different pairwise comparisons of Drosophila (D. melanogaster vs. D. simulans and D. simulans vs. D. sechellia). We find that divergence of these three classes of sites is strongly influenced by the recombination environment in which they are located, and this is especially true for the closer D. simulans vs. D. sechellia comparison. We suggest that this is mainly a result of the contribution of ancestral polymorphisms in different recombination regions. We also find that intergenic INE-1 elements are significantly more diverged than intronic INE-1 in both pairwise comparisons, implying the presence of either negative selection or lower mutation rates in introns. Furthermore, we show that substitution rates in INE-1 elements are not associated with the length of the noncoding sequence in which they are located, suggesting that reduced divergence in long noncoding sequences is not due to reduced mutation rates in these regions. Finally, we show that GC content for each site within INE-1 sequences has evolved toward an equilibrium value (approximately 33%) since insertion.

Comparison of dot chromosome sequences from D. melanogaster and D. virilis reveals an enrichment of DNA transposon sequences in heterochromatic domains

Sequencing and analysis of fosmid hybridization to the dot chromosomes of Drosophila virilis and D. melanogaster suggest that repetitive elements and density are important in determining higher-order chromatin packaging.

Background

Chromosome four of Drosophila melanogaster, known as the dot chromosome, is largely heterochromatic, as shown by immunofluorescent staining with antibodies to heterochromatin protein 1 (HP1) and histone H3K9me. In contrast, the absence of HP1 and H3K9me from the dot chromosome in D. virilis suggests that this region is euchromatic. D. virilis diverged from D. melanogaster 40 to 60 million years ago.

Results

Here we describe finished sequencing and analysis of 11 fosmids hybridizing to the dot chromosome of D. virilis (372,650 base-pairs) and seven fosmids from major euchromatic chromosome arms (273,110 base-pairs). Most genes from the dot chromosome of D. melanogaster remain on the dot chromosome in D. virilis, but many inversions have occurred. The dot chromosomes of both species are similar to the major chromosome arms in gene density and coding density, but the dot chromosome genes of both species have larger introns. The D. virilis dot chromosome fosmids have a high repeat density (22.8%), similar to homologous regions of D. melanogaster (26.5%). There are, however, major differences in the representation of repetitive elements. Remnants of DNA transposons make up only 6.3% of the D. virilis dot chromosome fosmids, but 18.4% of the homologous regions from D. melanogaster; DINE-1 and 1360 elements are particularly enriched in D. melanogaster. Euchromatic domains on the major chromosomes in both species have very few DNA transposons (less than 0.4 %).

Conclusion

Combining these results with recent findings about RNAi, we suggest that specific repetitive elements, as well as density, play a role in determining higher-order chromatin packaging.

Ribonuclease H evolution in retrotransposable elements

Eukaryotic and prokaryotic genomes encode either Type I or Type II Ribonuclease H (RNH) which is important for processing RNA primers that prime DNA replication in almost all organisms. This review highlights the important role that Type I RNH plays in the life cycle of many retroelements, and its utility in tracing early events in retroelement evolution. Many retroelements utilize host genome-encoded RNH, but several lineages of retroelements, including some non-LTR retroposons and all LTR retrotransposons, encode their own RNH domains. Examination of these RNH domains suggests that all LTR retrotransposons acquired an enzymatically weak RNH domain that is missing an important catalytic residue found in all other RNH enzymes. We propose that this reduced activity is essential to ensure correct processing of the polypurine tract (PPT), which is an important step in the life cycle of these retrotransposons. Vertebrate retroviruses appear to have reacquired their RNH domains, which are catalytically more active, but their ancestral RNH domains (found in other LTR retrotransposons) have degenerated to give rise to the tether domains unique to vertebrate retroviruses. The tether domain may serve to control the more active RNH domain of vertebrate retroviruses. Phylogenetic analysis of the RNH domains is also useful to "date" the relative ages of LTR and non-LTR retroelements. It appears that all LTR retrotransposons are as old as, or younger than, the "youngest" lineages of non-LTR retroelements, suggesting that LTR retrotransposons arose late in eukaryotes.

Penelope-like elements – a new class of retroelements: distribution, function and possible evolutionary significance

Here we describe a new class of retroelements termed PLE (Penelope-like elements). The only transpositionally active representative of this lineage found so far has been isolated from Drosophila virilis. This element, Penelope, is responsible for the hybrid dysgenesis syndrome in this species, characterized by simultaneous mobilization of several unrelated TE families in the progeny of dysgenic crosses. Several lines of evidence favor the hypothesis of recent Penelope invasion into D. virilis. Moreover, when D. virilisPenelope was introduced by P element-mediated transformation into the genome of D. melanogaster, it underwent extensive amplification in the new host and induced several traits of the dysgenesis syndrome, including gonadal atrophy and numerous mutations. The single ORF encoded by PLE consists of two principal domains: reverse transcriptase (RT) and endonuclease (EN), which is similar to GIY-YIG intron-encoded endonucleases. With the appearance of a large number of PLEs in genome databases from diverse eukaryotes, including amoebae, fungi, cnidarians, rotifers, flatworms, roundworms, fish, amphibia, and reptilia, it becomes possible to resolve their phylogenetic relationships with other RT groups with a greater degree of confidence. On the basis of their peculiar structural features, distinct phylogenetic placement, and structure of transcripts, we conclude that PLE constitute a novel class of eukaryotic retroelements, different from non-LTR and LTR retrotransposons.

Reverse transcriptase and endonuclease activities encoded by Penelope-like retroelements

Penelope-like elements are a class of retroelement that have now been identified in >50 species belonging to at least 10 animal phyla. The Penelope element isolated from Drosophila virilis is the only transpositionally active representative of this class isolated so far. The single ORF of Penelope and its relatives contains regions homologous to a reverse transcriptase of atypical structure and to the GIY-YIG, or Uri, an endonuclease (EN) domain not previously found in retroelements. We have expressed the single ORF of Penelope in a baculovirus expression system and have shown that it encodes a polyprotein with reverse transcriptase activity that requires divalent cations (Mn2+ and Mg2+). We have also expressed and purified the EN domain in Escherichia coli and have demonstrated that it has EN activity in vitro. Mutations in the conserved residues of the EN catalytic module abolish its nicking activity, whereas the DNA-binding properties of the mutant proteins remain unaffected. Only one strand of the target sequence is cleaved, and there is a certain degree of cleavage specificity. We propose that the Penelope EN cleaves the target DNA during transposition, generating a primer for reverse transcription. Our results show that an active Uri EN has been adopted by a retrotransposon.

Molecular paleontology of transposable elements in the Drosophila melanogaster genome

We report here a superfamily of "cut and paste" DNA transposons called Transib. These transposons populate the Drosophila melanogaster and Anopheles gambiae genomes, use a transposase that is not similar to any known proteins, and are characterized by 5-bp target site duplications. We found that the fly genome, which was thought to be colonized by the P element <100 years ago, harbors approximately 5 million year (Myr)-old fossils of ProtoP, an ancient ancestor of the P element. We also show that Hoppel, a previously reported transposable element (TE), is a nonautonomous derivate of ProtoP. We found that the "rolling-circle" Helitron transposons identified previously in plants and worms populate also insect genomes. Our results indicate that Helitrons were horizontally transferred into the fly or/and mosquito genomes. We have also identified a most abundant TE in the fly genome, DNAREP1_DM, which is an approximately 10-Myr-old footprint of a Penelope-like retrotransposon. We estimated that TEs are three times more abundant than reported previously, making up approximately 22% of the whole genome. The chromosomal and age distributions of TEs in D. melanogaster are very similar to those in Arabidopsis thaliana. Both genomes contain only relatively young TEs (<20 Myr old), constituting a main component of paracentromeric regions.