TrEMOLO: accurate transposable element allele frequency estimation using long-read sequencing data combining assembly and mapping-based approaches.

Transposable Element MOnitoring with LOng-reads (TrEMOLO) is a new software that combines assembly- and mapping-based approaches to robustly detect genetic elements called transposable elements (TEs). Using high- or low-quality genome assemblies, TrEMOLO can detect most TE insertions and deletions and estimate their allele frequency in populations. Benchmarking with simulated data revealed that TrEMOLO outperforms other state-of-the-art computational tools. TE detection and frequency estimation by TrEMOLO were validated using simulated and experimental datasets. Therefore, TrEMOLO is a comprehensive and suitable tool to accurately study TE dynamics. TrEMOLO is available under GNU GPL3.0 at https://github.com/DrosophilaGenomeEvolution/TrEMOLO .

A Transposon Story: From TE Content to TE Dynamic Invasion of Drosophila Genomes Using the Single-Molecule Sequencing Technology from Oxford Nanopore.

Transposable elements (TEs) are the main components of genomes. However, due to their repetitive nature, they are very difficult to study using data obtained with short-read sequencing technologies. Here, we describe an efficient pipeline to accurately recover TE insertion (TEI) sites and sequences from long reads obtained by Oxford Nanopore Technology (ONT) sequencing. With this pipeline, we could precisely describe the landscapes of the most recent TEIs in wild-type strains of Drosophila melanogaster and Drosophila simulans. Their comparison suggests that this subset of TE sequences is more similar than previously thought in these two species. The chromosome assemblies obtained using this pipeline also allowed recovering piRNA cluster sequences, which was impossible using short-read sequencing. Finally, we used our pipeline to analyze ONT sequencing data from a D. melanogaster unstable line in which LTR transposition was derepressed for 73 successive generations. We could rely on single reads to identify new insertions with intact target site duplications. Moreover, the detailed analysis of TEIs in the wild-type strains and the unstable line did not support the trap model claiming that piRNA clusters are hotspots of TE insertions.

The Mi-2 nucleosome remodeler and the Rpd3 histone deacetylase are involved in piRNA-guided heterochromatin formation.

In eukaryotes, trimethylation of lysine 9 on histone H3 (H3K9) is associated with transcriptional silencing of transposable elements (TEs). In drosophila ovaries, this heterochromatic repressive mark is thought to be deposited by SetDB1 on TE genomic loci after the initial recognition of nascent transcripts by PIWI-interacting RNAs (piRNAs) loaded on the Piwi protein. Here, we show that the nucleosome remodeler Mi-2, in complex with its partner MEP-1, forms a subunit that is transiently associated, in a MEP-1 C-terminus-dependent manner, with known Piwi interactors, including a recently reported SUMO ligase, Su(var)2-10. Together with the histone deacetylase Rpd3, this module is involved in the piRNA-dependent TE silencing, correlated with H3K9 deacetylation and trimethylation. Therefore, drosophila piRNA-mediated transcriptional silencing involves three epigenetic effectors, a remodeler, Mi-2, an eraser, Rpd3 and a writer, SetDB1, in addition to the Su(var)2-10 SUMO ligase.

The somatic piRNA pathway controls germline transposition over generations.

Transposable elements (TEs) are parasitic DNA sequences that threaten genome integrity by replicative transposition in host gonads. The Piwi-interacting RNAs (piRNAs) pathway is assumed to maintain Drosophila genome homeostasis by downregulating transcriptional and post-transcriptional TE expression in the ovary. However, the bursts of transposition that are expected to follow transposome derepression after piRNA pathway impairment have not yet been reported. Here, we show, at a genome-wide level, that piRNA loss in the ovarian somatic cells boosts several families of the endogenous retroviral subclass of TEs, at various steps of their replication cycle, from somatic transcription to germinal genome invasion. For some of these TEs, the derepression caused by the loss of piRNAs is backed up by another small RNA pathway (siRNAs) operating in somatic tissues at the post transcriptional level. Derepressed transposition during 70 successive generations of piRNA loss exponentially increases the genomic copy number by up to 10-fold.

Piwi Is Required during Drosophila Embryogenesis to License Dual-Strand piRNA Clusters for Transposon Repression in Adult Ovaries.

Most piRNAs in the Drosophila female germline are transcribed from heterochromatic regions called dual-strand piRNA clusters. Histone 3 lysine 9 trimethylation (H3K9me3) is required for licensing piRNA production by these clusters. However, it is unclear when and how they acquire this permissive heterochromatic state. Here, we show that transient Piwi depletion in Drosophila embryos results in H3K9me3 decrease at piRNA clusters in ovaries. This is accompanied by impaired biogenesis of ovarian piRNAs, accumulation of transposable element transcripts, and female sterility. Conversely, Piwi depletion at later developmental stages does not disturb piRNA cluster licensing. These results indicate that the identity of piRNA clusters is epigenetically acquired in a Piwi-dependent manner during embryonic development, which is reminiscent of the widespread genome reprogramming occurring during early mammalian zygotic development.

MicroRNA-Dependent Transcriptional Silencing of Transposable Elements in Drosophila Follicle Cells.

RNA interference-related silencing mechanisms concern very diverse and distinct biological processes, from gene regulation (via the microRNA pathway) to defense against molecular parasites (through the small interfering RNA and the Piwi-interacting RNA pathways). Small non-coding RNAs serve as specificity factors that guide effector proteins to ribonucleic acid targets via base-pairing interactions, to achieve transcriptional or post-transcriptional regulation. Because of the small sequence complementarity required for microRNA-dependent post-transcriptional regulation, thousands of microRNA (miRNA) putative targets have been annotated in Drosophila. In Drosophila somatic ovarian cells, genomic parasites, such as transposable elements (TEs), are transcriptionally repressed by chromatin changes induced by Piwi-interacting RNAs (piRNAs) that prevent them from invading the germinal genome. Here we show, for the first time, that a functional miRNA pathway is required for the piRNA-mediated transcriptional silencing of TEs in this tissue. Global miRNA depletion, caused by tissue- and stage-specific knock down of drosha (involved in miRNA biogenesis), AGO1 or gawky (both responsible for miRNA activity), resulted in loss of TE-derived piRNAs and chromatin-mediated transcriptional de-silencing of TEs. This specific TE de-repression was also observed upon individual titration (by expression of the complementary miRNA sponge) of two miRNAs (miR-14 and miR-34) as well as in a miR-14 loss-of-function mutant background. Interestingly, the miRNA defects differentially affected TE- and 3' UTR-derived piRNAs. To our knowledge, this is the first indication of possible differences in the biogenesis or stability of TE- and 3' UTR-derived piRNAs. This work is one of the examples of detectable phenotypes caused by loss of individual miRNAs in Drosophila and the first genetic evidence that miRNAs have a role in the maintenance of genome stability via piRNA-mediated TE repression.

A user-friendly chromatographic method to purify small regulatory RNAs.

The discovery of the small regulatory RNAs has changed our vision of cellular regulations. Indeed, when loaded on Argonaute proteins they form ribonucleoprotein complexes (RNPs) that target complementary sequences to achieve widespread silencing mechanisms conserved in most eukaryotes. The recent development of deep sequencing approaches highly contributed to their detection. Small RNA isolation from cells and/or tissues remains a crucial stage to generate robust and relevant sequencing data. In 2006, a novel strategy based on anion-exchange chromatography has been proposed as an alternative to the standard size-isolation purification procedure. Using bioinformatic comparative analysis, we here demonstrate that anion-exchange chromatographic RNP purification prior to small RNA extraction unbiasedly enriches datasets in bona fide reads (small regulatory RNA sequences) and depletes endogenous contaminants (ribosomal RNAs and degradation RNA products). The resulting increase in sequencing depth provides a major benefit to study rare populations. We then developed a fast and basic manual procedure to purify such small non-coding RNAs using anion-exchange chromatography at the bench. We validated the efficiency of this new method and used this strategy to purify small RNAs from various tissues and organisms. We moreover determined that our manual purification increases the output of the previously described anion-exchange chromatography procedure.

piRNA-mediated transgenerational inheritance of an acquired trait.

The maintenance of genome integrity is an essential trait to the successful transmission of genetic information. In animal germ cells, piRNAs guide PIWI proteins to silence transposable elements (TEs) in order to maintain genome integrity. In insects, most TE silencing in the germline is achieved by secondary piRNAs that are produced by a feed-forward loop (the ping-pong cycle), which requires the piRNA-directed cleavage of two types of RNAs: mRNAs of functional euchromatic TEs and heterochromatic transcripts that contain defective TE sequences. The first cleavage that initiates such an amplification loop remains poorly understood. Taking advantage of the existence of strains that are devoid of functional copies of the LINE-like I-element, we report here that in such Drosophila ovaries, the initiation of a ping-pong cycle is exclusively achieved by secondary I-element piRNAs that are produced in the ovary and deposited in the embryonic germline. This unusual secondary piRNA biogenesis, detected in the absence of functional I-element copies, results from the processing of sense and antisense transcripts of several different defective I-element. Once acquired, for instance after ancestor aging, this capacity to produce heterochromatic-only secondary piRNAs is partially transmitted through generations via maternal piRNAs. Furthermore, such piRNAs acting as ping-pong initiators in a chromatin-independent manner confer to the progeny a high capacity to repress the I-element mobility. Our study explains, at the molecular level, the basis for epigenetic memory of maternal immunity that protects females from hybrid dysgenesis caused by transposition of paternally inherited functional I-element.

Maternal mRNA deadenylation and decay by the piRNA pathway in the early Drosophila embryo.

Piwi-associated RNAs (piRNAs), a specific class of 24- to 30-nucleotide-long RNAs produced by the Piwi-type of Argonaute proteins, have a specific germline function in repressing transposable elements. This repression is thought to involve heterochromatin formation and transcriptional and post-transcriptional silencing. The piRNA pathway has other essential functions in germline stem cell maintenance and in maintaining germline DNA integrity. Here we uncover an unexpected function of the piRNA pathway in the decay of maternal messenger RNAs and in translational repression in the early embryo. A subset of maternal mRNAs is degraded in the embryo at the maternal-to-zygotic transition. In Drosophila, maternal mRNA degradation depends on the RNA-binding protein Smaug and the deadenylase CCR4, as well as the zygotic expression of a microRNA cluster. Using mRNA encoding the embryonic posterior morphogen Nanos (Nos) as a paradigm to study maternal mRNA decay, we found that CCR4-mediated deadenylation of nos depends on components of the piRNA pathway including piRNAs complementary to a specific region in the nos 3' untranslated region. Reduced deadenylation when piRNA-induced regulation is impaired correlates with nos mRNA stabilization and translational derepression in the embryo, resulting in head development defects. Aubergine, one of the Argonaute proteins in the piRNA pathway, is present in a complex with Smaug, CCR4, nos mRNA and piRNAs that target the nos 3' untranslated region, in the bulk of the embryo. We propose that piRNAs and their associated proteins act together with Smaug to recruit the CCR4 deadenylation complex to specific mRNAs, thus promoting their decay. Because the piRNAs involved in this regulation are produced from transposable elements, this identifies a direct developmental function for transposable elements in the regulation of gene expression.

piRNA-mediated nuclear accumulation of retrotransposon transcripts in the Drosophila female germline.

Germline silencing of transposable elements is essential for the maintenance of genome integrity. Recent results indicate that this repression is largely achieved through a RNA silencing pathway that involves Piwi-interacting RNAs (piRNAs). However the repressive mechanisms are not well understood. To address this question, we used the possibility to disrupt the repression of the Drosophila I element retrotransposon by hybrid dysgenesis. We show here that the repression of the functional I elements that are located in euchromatin requires proteins of the piRNA pathway, and that the amount of ovarian I element piRNAs correlates with the strength of the repression in the female germline. Antisense RNAs, which are likely used to produce antisense piRNAs, are transcribed by heterochromatic defective I elements, but efficient production of these antisense small RNAs requires the presence in the genome of euchromatic functional I elements. Finally, we demonstrate that the piRNA-induced silencing of the functional I elements is at least partially posttranscriptional. In a repressive background, these elements are still transcribed, but some of their sense transcripts are kept in nurse cell nuclear foci together with those of the Doc retrotransposon. In the absence of I element piRNAs, either in dysgenic females or in mutants of the piRNA silencing pathway, sense I element transcripts are transported toward the oocyte where retrotransposition occurs. Our results indicate that piRNAs are involved in a posttranscriptional gene-silencing mechanism resulting in RNA nuclear accumulation.

The flamenco locus controls the gypsy and ZAM retroviruses and is required for Drosophila oogenesis.

In Drosophila, the as yet uncloned heterochromatic locus flamenco (flam) controls mobilization of the endogenous retrovirus gypsy through the repeat-associated small interfering (rasi) RNA silencing pathway. Restrictive alleles (flamR) downregulate accumulation of gypsy transcripts in the somatic follicular epithelium of the ovary. In contrast, permissive alleles (flamP) are unable to repress gypsy. DIP1, the closest transcription unit to a flam-insertional mutation, was considered as a good candidate to be a gypsy regulator, since it encodes a dsRNA-binding protein. To further characterize the locus we analyzed P-induced flam mutants and generated new mutations by transposon mobilization. We show that flam is required somatically for morphogenesis of the follicular epithelium, the tissue where gypsy is repressed. This developmental activity is necessary to control gypsy and another retroelement, ZAM. We also show that flam is not DIP1, as none of the new permissive mutants affect the DIP1 coding sequence. In addition, two deletions removing DIP1 coding sequences do not affect any of the flamenco functions. Our results suggest that flamenco extends proximally to DIP1, spanning >130 kb of transposon-rich heterochromatin. We propose a model explaining the multiple functions of this large heterochromatic locus.

A novel repeat-associated small interfering RNA-mediated silencing pathway downregulates complementary sense gypsy transcripts in somatic cells of the Drosophila ovary.

Replication of the gypsy endogenous retrovirus involves contamination of the female germ line by adjacent somatic tissues. This is prevented by flam, an as-yet-uncloned heterochromatic pericentromeric locus, at the level of transcript accumulation in these somatic ovarian tissues. We tested the effect of a presumptive RNA silencing mechanism on the accumulation of RNAs produced by constructs containing various gypsy sequences and report that the efficiency of silencing is indeed correlated with the amount of complementary RNAs, 25 to 30 nucleotides in length, in the ovary. For instance, while these RNAs were found to display a three- to fivefold excess of the antisense strands, only the transcripts that contain the complementary sense gypsy sequences could be repressed, indicating that they are targeted at the RNA, not DNA, level. Their size and asymmetry in strand polarity are typical of the novel repeat-associated small interfering RNA (rasiRNA)-mediated pathway, recently suspected to prevent the deleterious expression of selfish DNA specifically in the germ line. Unlike microRNAs (but like rasiRNAs and, surprisingly, siRNAs as well), gypsy rasiRNAs are modified at the 3' end. The rasiRNA-associated protein Piwi (but not Aub) is required for gypsy silencing, whereas Dicer-2 (which makes siRNAs) is not. In contrast, piwi, aub, and flam do not appear to affect somatic siRNA-mediated silencing. The amount of gypsy rasiRNAs is genetically determined by the flam locus in a provirus copy number-independent manner and is triggered in the somatic tissues by some pericentromeric provirus(es), which are thereby able to protect the germ line from retroviral invasion.

Restrictive flamenco alleles are maintained in Drosophila melanogaster population cages, despite the absence of their endogenous gypsy retroviral targets.

The flamenco (flam) locus, located at 20A1-3 in the centromeric heterochromatin of the Drosophila melanogaster X chromosome, is a major regulator of the gypsy/mdg4 endogenous retrovirus. In restrictive strains, functional flam alleles maintain gypsy proviruses in a repressed state. By contrast, in permissive strains, proviral amplification results from infection of the female germ line and subsequent insertions into the chromosomes of the progeny. A restrictive/permissive polymorphism prevails in natural and laboratory populations. This polymorphism was assumed to be maintained by the interplay of opposite selective forces; on one hand, the increase of genetic load caused by proviral insertions would favor restrictive flam alleles because they make flies resistant to these gypsy replicative transpositions and, on the other, a hypothetical resistance cost would select against such alleles in the absence of the retrovirus. However, the population cage data presented in this paper do not fit with this simple resistance cost hypothesis because restrictive alleles were not eliminated in the absence of functional gypsy proviruses; on the contrary, using 2 independent flam allelic pairs, the restrictive frequency rose to about 90% in every experimental population, whatever the pair of alleles and the allelic proportions in the initial inoculum. These data suggest that the flam polymorphism is maintained by some strong balancing selection, which would act either on flam itself, independently of the deleterious effect of gypsy, or on a hypothetical flanking gene, in linkage disequilibrium with flam. Alternatively, restrictive flam alleles might also be resistant to some other retroelements that would be still present in the cage populations, causing a positive selection for these alleles. Whatever selective forces that maintain high levels of restrictive alleles independently of gypsy, this unknown mechanism can set up an interesting kind of antiviral innate immunity, at the population level.

Characterization of a nucleocapsid-like region and of two distinct primer tRNALys,2 binding sites in the endogenous retrovirus Gypsy.

Mobile LTR-retroelements comprising retroviruses and LTR-retrotransposons form a large part of eukaryotic genomes. Their mode of replication and abundance favour the notion that they are major actors in eukaryote evolution. The Gypsy retroelement can spread in the germ line of the fruit fly Drosophila melanogaster via both env-independent and env-dependent processes. Thus, Gypsy is both an active retrotransposon and an infectious retrovirus resembling the gammaretrovirus MuLV. However, unlike gammaretroviruses, the Gypsy Gag structural precursor is not processed into Matrix, Capsid and Nucleocapsid (NC) proteins. In contrast, it has features in common with Gag of the ancient yeast TY1 retroelement. These characteristics of Gypsy make it a very interesting model to study replication of a retroelement at the frontier between ancient retrotransposons and retroviruses. We investigated Gypsy replication using an in vitro model system and transfection of insect cells. Results show that an unstructured domain of Gypsy Gag has all the properties of a retroviral NC. This NC-like peptide forms ribonucleoparticle-like complexes upon binding Gypsy RNA and directs the annealing of primer tRNA(Lys,2) to two distinct primer binding sites (PBS) at the genome 5' and 3' ends. Only the 5' PBS is indispensable for cDNA synthesis in vitro and in Drosophila cells.

Evidence for a piwi-dependent RNA silencing of the gypsy endogenous retrovirus by the Drosophila melanogaster flamenco gene.

In Drosophila melanogaster, the endogenous retrovirus gypsy is repressed by the functional alleles (restrictive) of an as-yet-uncloned heterochromatic gene called flamenco. Using gypsy-lacZ transcriptional fusions, we show here that this repression takes place not only in the follicle cells of restrictive ovaries, as was previously observed, but also in restrictive larval female gonads. Analyses of the role of gypsy cis-regulatory sequences in the control of gypsy expression are also presented. They rule out the hypothesis that gypsy would contain a single binding region for a putative Flamenco repressor. Indeed, the ovarian expression of a chimeric yp3-lacZ construct was shown to become sensitive to the Flamenco regulation when any of three different 5'-UTR gypsy sequences (ranging from 59 to 647 nucleotides) was incorporated into the heterologous yp3-lacZ transcript. The piwi mutation, which is known to affect RNA-mediated homology-dependent transgene silencing, was also shown to impede the repression of gypsy in restrictive female gonads. Finally, a RNA-silencing model is also supported by the finding in ovaries of short RNAs (25-27 nucleotides long) homologous to sequences from within the gypsy 5'-UTR.

A novel double-stranded RNA-binding protein, disco interacting protein 1 (DIP1), contributes to cell fate decisions during Drosophila development.

We report the identification of the Disco Interacting Protein 1 (DIP1) gene isolated in a yeast interaction trap screen using the zinc finger protein disconnected (disco) as a bait. DIP1 encodes a protein containing two double-stranded RNA binding domains (dsRBD). Consistent with the presence of dsRBD, DIP1 binds dsRNA or structured RNAs in Northwestern assays. DIP1 is found in nuclear subdomains resembling speckles known to accumulate transcription and splicing factors. In early embryos, nuclear localization of DIP1 protein coincides with the onset of zygotic gene expression. Later in development DIP1 expression is decreased in dividing cells in different tissues. Overexpression of DIP1 in the eye-antennal imaginal disc, early in embryonic and larval development, causes the formation of supernumerary structures in the head capsule. A role for DIP1 in epigenetic mechanisms that lead to the establishment and/or maintenance of cell fate specification is discussed.

Comparative and functional studies of Drosophila species invasion by the gypsy endogenous retrovirus.

Gypsy is an endogenous retrovirus of Drosophila melanogaster. Phylogenetic studies suggest that occasional horizontal transfer events of gypsy occur between Drosophila species. gypsy possesses infective properties associated with the products of the envelope gene that might be at the origin of these interspecies transfers. We report here the existence of DNA sequences putatively encoding full-length Env proteins in the genomes of Drosophila species other than D. melanogaster, suggesting that potentially infective gypsy copies able to spread between sexually isolated species can occur. The ability of gypsy to invade the genome of a new species is conditioned by its capacity to be expressed in the naive genome. The genetic basis for the regulation of gypsy activity in D. melanogaster is now well known, and it has been assigned to an X-linked gene called flamenco. We established an experimental simulation of the invasion of the D. melanogaster genome by gypsy elements derived from other Drosophila species, which demonstrates that these non- D. melanogaster gypsy elements escape the repression exerted by the D. melanogaster flamenco gene.