Short and long-term evolutionary dynamics of subtelomeric piRNA clusters in Drosophila

The composition of piRNA clusters in Drosophila melanogaster deviates from expectations under the trap model

BACKGROUND

It is widely assumed that the invasion of a transposable element (TE) in mammals and invertebrates is stopped when a copy of the TE jumps into a piRNA cluster (i.e., the trap model). However, recent works, which for example showed that deletion of three major piRNA clusters has no effect on TE activity, cast doubt on the trap model.

RESULTS

Here, we test the trap model from a population genetics perspective. Our simulations show that the composition of regions that act as transposon traps (i.e., potentially piRNA clusters) ought to deviate from regions that have no effect on TE activity. We investigated TEs in five Drosophila melanogaster strains using three complementary approaches to test whether the composition of piRNA clusters matches these expectations. We found that the abundance of TE families inside and outside of piRNA clusters is highly correlated, although this is not expected under the trap model. Furthermore, the distribution of the number of TE insertions in piRNA clusters is also much broader than expected.

CONCLUSIONS

We found that the observed composition of piRNA clusters is not in agreement with expectations under the simple trap model. Dispersed piRNA producing TE insertions and temporal as well as spatial heterogeneity of piRNA clusters may account for these deviations.

Modeling early germline immunization after horizontal transfer of transposable elements reveals internal piRNA cluster heterogeneity

BACKGROUND

A fraction of all genomes is composed of transposable elements (TEs) whose mobility needs to be carefully controlled. In gonads, TE activity is repressed by PIWI-interacting RNAs (piRNAs), a class of small RNAs synthesized by heterochromatic loci enriched in TE fragments, called piRNA clusters. Maintenance of active piRNA clusters across generations is secured by maternal piRNA inheritance providing the memory for TE repression. On rare occasions, genomes encounter horizontal transfer (HT) of new TEs with no piRNA targeting them, threatening the host genome integrity. Naïve genomes can eventually start to produce new piRNAs against these genomic invaders, but the timing of their emergence remains elusive.

RESULTS

Using a set of TE-derived transgenes inserted in different germline piRNA clusters and functional assays, we have modeled a TE HT in Drosophila melanogaster. We have found that the complete co-option of these transgenes by a germline piRNA cluster can occur within four generations associated with the production of new piRNAs all along the transgenes and the germline silencing of piRNA sensors. Synthesis of new transgenic TE piRNAs is linked to piRNA cluster transcription dependent on Moonshiner and heterochromatin mark deposition that propagates more efficiently on short sequences. Moreover, we found that sequences located within piRNA clusters can have different piRNA profiles and can influence transcript accumulation of nearby sequences.

CONCLUSIONS

Our study reveals that genetic and epigenetic properties, such as transcription, piRNA profiles, heterochromatin, and conversion efficiency along piRNA clusters, could be heterogeneous depending on the sequences that compose them. These findings suggest that the capacity of transcriptional signal erasure induced by the chromatin complex specific of the piRNA cluster can be incomplete through the piRNA cluster loci. Finally, these results have revealed an unexpected level of complexity that highlights a new magnitude of piRNA cluster plasticity fundamental for the maintenance of genome integrity.

The histone demethylase Kdm3 prevents auto-immune piRNAs production in Drosophila

Genome integrity of the animal germline is protected from transposable element activity by PIWI-interacting RNAs (piRNAs). While piRNA biogenesis is intensively explored, little is known about the genetical determination of piRNA clusters, the genomic sources of piRNAs. Using a bimodal epigenetic state piRNA cluster (BX2), we identified the histone demethylase Kdm3 as being able to prevent a cryptic piRNA production. In the absence of Kdm3, dozens of coding gene-containing regions become genuine germline dual-strand piRNA clusters. Eggs laid by Kdm3 mutant females show developmental defects phenocopying loss of function of genes embedded into the additional piRNA clusters, suggesting an inheritance of functional ovarian "auto-immune" piRNAs. Antagonizing piRNA cluster determination through chromatin modifications appears crucial to prevent auto-immune genic piRNAs production.

Evolutionary dynamics of piRNA clusters in Drosophila

Small RNAs produced from transposable element (TE)-rich sections of the genome, termed piRNA clusters, are a crucial component in the genomic defence against selfish DNA. In animals, it is thought the invasion of a TE is stopped when a copy of the TE inserts into a piRNA cluster, triggering the production of cognate small RNAs that silence the TE. Despite this importance for TE control, little is known about the evolutionary dynamics of piRNA clusters, mostly because these repeat-rich regions are difficult to assemble and compare. Here, we establish a framework for studying the evolution of piRNA clusters quantitatively. Previously introduced quality metrics and a newly developed software for multiple alignments of repeat annotations (Manna) allow us to estimate the level of polymorphism segregating in piRNA clusters and the divergence among homologous piRNA clusters. By studying 20 conserved piRNA clusters in multiple assemblies of four Drosophila species, we show that piRNA clusters are evolving rapidly. While 70%-80% of the clusters are conserved within species, the clusters share almost no similarity between species as closely related as D. melanogaster and D. simulans. Furthermore, abundant insertions and deletions are segregating within the Drosophila species. We show that the evolution of clusters is mainly driven by large insertions of recently active TEs and smaller deletions mostly in older TEs. The effect of these forces is so rapid that homologous clusters often do not contain insertions from the same TE families.

Mod(mdg4) variants repress telomeric retrotransposon HeT-A by blocking subtelomeric enhancers

Telomeres in Drosophila are composed of sequential non-LTR retrotransposons HeT-A, TART and TAHRE. Although they are repressed by the PIWI-piRNA pathway or heterochromatin in the germline, the regulation of these retrotransposons in somatic cells is poorly understood. In this study, we demonstrated that specific splice variants of Mod(mdg4) repress HeT-A by blocking subtelomeric enhancers in ovarian somatic cells. Among the variants, we found that the Mod(mdg4)-N variant represses HeT-A expression the most efficiently. Subtelomeric sequences bound by Mod(mdg4)-N block enhancer activity within subtelomeric TAS-R repeats. This enhancer-blocking activity is increased by the tandem association of Mod(mdg4)-N to repetitive subtelomeric sequences. In addition, the association of Mod(mdg4)-N couples with the recruitment of RNA polymerase II to the subtelomeres, which reinforces its enhancer-blocking function. Our findings provide novel insights into how telomeric retrotransposons are regulated by the specific variants of insulator proteins associated with subtelomeric sequences.

The transposition rate has little influence on the plateauing level of the P-element

The popular trap model assumes that the invasions of transposable elements (TEs) in mammals and invertebrates are stopped by piRNAs that emerge after insertion of the TE into a piRNA cluster. It remains, however, still unclear which factors influence the dynamics of TE invasions. The activity of the TE (i.e., transposition rate) is one frequently discussed key factor. Here we take advantage of the temperature-dependent activity of the P-element, a widely studied eukaryotic TE, to test how TE activity affects the dynamics of a TE invasion. We monitored P-element invasion dynamics in experimental Drosophila simulans populations at hot and cold culture conditions. Despite marked differences in transposition rates, the P-element reached very similar copy numbers at both temperatures. The reduction of the insertion rate upon approaching the copy number plateau was accompanied by similar amounts of piRNAs against the P-element at both temperatures. Nevertheless, we also observed fewer P-element insertions in piRNA clusters than expected, which is not compatible with a simple trap model. The ping-pong cycle, which degrades TE transcripts, becomes typically active after the copy number plateaued. We generated a model, with few parameters, that largely captures the observed invasion dynamics. We conclude that the transposition rate has at the most only a minor influence on TE abundance, but other factors, such as paramutations or selection against TE insertions are shaping the TE composition.

Novel quality metrics allow identifying and generating high-quality assemblies of piRNA clusters

In most animals, it is thought that the proliferation of a transposable element (TE) is stopped when the TE jumps into a piRNA cluster. Despite this central importance, little is known about the composition and the evolutionary dynamics of piRNA clusters. This is largely because piRNA clusters are notoriously difficult to assemble as they are frequently composed of highly repetitive DNA. With long reads, we may finally be able to obtain reliable assemblies of piRNA clusters. Unfortunately, it is unclear how to generate and identify the best assemblies, as many assembly strategies exist and standard quality metrics are ignorant of TEs. To address these problems, we introduce several novel quality metrics that assess: (a) the fraction of completely assembled piRNA clusters, (b) the quality of the assembled clusters and (c) whether an assembly captures the overall TE landscape of an organisms (i.e. the abundance, the number of SNPs and internal deletions of all TE families). The requirements for computing these metrics vary, ranging from annotations of piRNA clusters to consensus sequences of TEs and genomic sequencing data. Using these novel metrics, we evaluate the effect of assembly algorithm, polishing, read length, coverage, residual polymorphisms and finally identify strategies that yield reliable assemblies of piRNA clusters. Based on an optimized approach, we provide assemblies for the two Drosophila melanogaster strains Canton-S and Pi2. About 80% of known piRNA clusters were assembled in both strains. Finally, we demonstrate the generality of our approach by extending our metrics to humans and Arabidopsis thaliana.

Roles of piRNAs in transposon and pseudogene regulation of germline mRNAs and lncRNAs

PIWI proteins, a subfamily of PAZ/PIWI Domain family RNA-binding proteins, are best known for their function in silencing transposons and germline development by partnering with small noncoding RNAs called PIWI-interacting RNAs (piRNAs). However, recent studies have revealed multifaceted roles of the PIWI-piRNA pathway in regulating the expression of other major classes of RNAs in germ cells. In this review, we summarize how PIWI proteins and piRNAs regulate the expression of many disparate RNAs, describing a highly complex global genomic regulatory relationship at the RNA level through which piRNAs functionally connect all major constituents of the genome in the germline.

Identification and Characterization of Breakpoints and Mutations on Drosophila melanogaster Balancer Chromosomes

Balancers are rearranged chromosomes used in Drosophila melanogaster to maintain deleterious mutations in stable populations, preserve sets of linked genetic elements and construct complex experimental stocks. Here, we assess the phenotypes associated with breakpoint-induced mutations on commonly used third chromosome balancers and show remarkably few deleterious effects. We demonstrate that a breakpoint in p53 causes loss of radiation-induced apoptosis and a breakpoint in Fucosyltransferase A causes loss of fucosylation in nervous and intestinal tissue—the latter study providing new markers for intestinal cell identity and challenging previous conclusions about the regulation of fucosylation. We also describe thousands of potentially harmful mutations shared among X or third chromosome balancers, or unique to specific balancers, including an Ankyrin 2 mutation present on most TM3 balancers, and reiterate the risks of using balancers as experimental controls. We used long-read sequencing to confirm or refine the positions of two inversions with breakpoints lying in repetitive sequences and provide evidence that one of the inversions, In(2L)Cy, arose by ectopic recombination between foldback transposon insertions and the other, In(3R)C, cleanly separates subtelomeric and telomeric sequences and moves the subtelomeric sequences to an internal chromosome position. In addition, our characterization of In(3R)C shows that balancers may be polymorphic for terminal deletions. Finally, we present evidence that extremely distal mutations on balancers can add to the stability of stocks whose purpose is to maintain homologous chromosomes carrying mutations in distal genes. Overall, these studies add to our understanding of the structure, diversity and effectiveness of balancer chromosomes.

piRNA Clusters Need a Minimum Size to Control Transposable Element Invasions

piRNA clusters are thought to repress transposable element (TE) activity in mammals and invertebrates. Here, we show that a simple population genetics model reveals a constraint on the size of piRNA clusters: The total size of the piRNA clusters of an organism must exceed 0.2% of a genome to repress TE invasions. Moreover, larger piRNA clusters accounting for up to 3% of the genome may be necessary when populations are small, transposition rates are high, and TE insertions are recessive. If piRNA clusters are too small, the load of deleterious TE insertions that accumulate during a TE invasion may drive populations extinct before an effective piRNA-based defense against the TE can be established. Our findings are solely based on three well-supported assumptions: 1) TEs multiply within genomes, 2) TEs are mostly deleterious, and 3) piRNA clusters act as transposon traps, where a single insertion in a cluster silences all TE copies in trans. Interestingly, the piRNA clusters of some species meet our observed minimum size requirements, whereas the clusters of other species do not. Species with small piRNA clusters, such as humans and mice, may experience severe fitness reductions during invasions of novel TEs, which is possibly even threatening the persistence of some populations. This work also raises the important question of how piRNA clusters evolve. We propose that the size of piRNA clusters may be at an equilibrium between evolutionary forces that act to expand and contract piRNA clusters.

Rapid evolution of piRNA-mediated silencing of an invading transposable element was driven by abundant de novo mutations

The regulation of transposable element (TE) activity by small RNAs is a ubiquitous feature of germlines. However, despite the obvious benefits to the host in terms of ensuring the production of viable gametes and maintaining the integrity of the genomes they carry, it remains controversial whether TE regulation evolves adaptively. We examined the emergence and evolutionary dynamics of repressor alleles after P-elements invaded the Drosophila melanogaster genome in the mid-twentieth century. In many animals including Drosophila, repressor alleles are produced by transpositional insertions into piRNA clusters, genomic regions encoding the Piwi-interacting RNAs (piRNAs) that regulate TEs. We discovered that ∼94% of recently collected isofemale lines in the Drosophila melanogaster Genetic Reference Panel (DGRP) contain at least one P-element insertion in a piRNA cluster, indicating that repressor alleles are produced by de novo insertion at an exceptional rate. Furthermore, in our sample of approximately 200 genomes, we uncovered no fewer than 80 unique P-element insertion alleles in at least 15 different piRNA clusters. Finally, we observe no footprint of positive selection on P-element insertions in piRNA clusters, suggesting that the rapid evolution of piRNA-mediated repression in D. melanogaster was driven primarily by mutation. Our results reveal for the first time how the unique genetic architecture of piRNA production, in which numerous piRNA clusters can encode regulatory small RNAs upon transpositional insertion, facilitates the nonadaptive rapid evolution of repression.

Nuclear Ccr4-Not mediates the degradation of telomeric and transposon transcripts at chromatin in the Drosophila germline

Ccr4-Not is a highly conserved complex involved in cotranscriptional RNA surveillance pathways in yeast. In Drosophila, Ccr4-Not is linked to the translational repression of miRNA targets and the posttranscriptional control of maternal mRNAs during oogenesis and embryonic development. Here, we describe a new role for the Ccr4-Not complex in nuclear RNA metabolism in the Drosophila germline. Ccr4 depletion results in the accumulation of transposable and telomeric repeat transcripts in the fraction of chromatin-associated RNA; however, it does not affect small RNA levels or the heterochromatin state of the target loci. Nuclear targets of Ccr4 mainly comprise active full-length transposable elements (TEs) and telomeric and subtelomeric repeats. Moreover, Ccr4-Not foci localize at telomeres in a Piwi-dependent manner, suggesting a functional relationship between these pathways. Indeed, we detected interactions between the components of the Ccr4-Not complex and piRNA machinery, which indicates that these pathways cooperate in the nucleus to recognize and degrade TE transcripts at transcription sites. These data reveal a new layer of transposon control in the germline, which is critical for the maintenance of genome integrity.

The Telomere Paradox: Stable Genome Preservation with Rapidly Evolving Proteins

Telomeres ensure chromosome length homeostasis and protection from catastrophic end-to-end chromosome fusions. All eukaryotes require this essential, strictly conserved telomere-dependent genome preservation. However, recent evolutionary analyses of mammals, plants, and flies report pervasive rapid evolution of telomere proteins. The causes of this paradoxical observation - that unconserved machinery underlies an essential, conserved function - remain enigmatic. Indeed, these fast-evolving telomere proteins bind, extend, and protect telomeric DNA, which itself evolves slowly in most systems. We hypothesize that the universally fast-evolving subtelomere - the telomere-adjacent, repetitive sequence - is a primary driver of the 'telomere paradox'. Under this model, radical sequence changes in the subtelomere perturb subtelomere-dependent, telomere functions. Compromised telomere function then spurs adaptation of telomere proteins to maintain telomere length homeostasis and protection. We propose an experimental framework that leverages both protein divergence and subtelomeric sequence divergence to test the hypothesis that subtelomere sequence evolution shapes recurrent innovation of telomere machinery.

Environmentally-Induced Transgenerational Epigenetic Inheritance: Implication of PIWI Interacting RNAs

Environmentally-induced transgenerational epigenetic inheritance is an emerging field. The understanding of associated epigenetic mechanisms is currently in progress with open questions still remaining. In this review, we present an overview of the knowledge of environmentally-induced transgenerational inheritance and associated epigenetic mechanisms, mainly in animals. The second part focuses on the role of PIWI-interacting RNAs (piRNAs), a class of small RNAs involved in the maintenance of the germline genome, in epigenetic memory to put into perspective cases of environmentally-induced transgenerational inheritance involving piRNA production. Finally, the last part addresses how genomes are facing production of new piRNAs, and from a broader perspective, how this process might have consequences on evolution and on sporadic disease development.

Environmentally-induced epigenetic conversion of a piRNA cluster

Transposable element (TE) activity is repressed in animal gonads by PIWI-interacting RNAs (piRNAs) produced by piRNA clusters. Current models in flies propose that germinal piRNA clusters are functionally defined by the maternal inheritance of piRNAs produced during the previous generation. Taking advantage of an inactive, but ready to go, cluster of P-element derived transgene insertions in Drosophila melanogaster, we show here that raising flies at high temperature (29°C) instead of 25°C triggers the stable conversion of this locus from inactive into actively producing functional piRNAs. The increase of antisense transcripts from the cluster at 29°C combined with the requirement of transcription of euchromatic homologous sequences, suggests a role of double stranded RNA in the production of de novo piRNAs. This report describes the first case of the establishment of an active piRNA cluster by environmental changes in the absence of maternal inheritance of homologous piRNAs.

Editorial note

This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

Insights Into an Unexplored Component of the Mosquito Repeatome: Distribution and Variability of Viral Sequences Integrated Into the Genome of the Arboviral Vector Aedes albopictus

The Asian tiger mosquito Aedes albopictus is an invasive mosquito and a competent vector for public-health relevant arboviruses such as Chikungunya (Alphavirus), Dengue and Zika (Flavivirus) viruses. Unexpectedly, the sequencing of the genome of this mosquito revealed an unusually high number of integrated sequences with similarities to non-retroviral RNA viruses of the Flavivirus and Rhabdovirus genera. These Non-retroviral Integrated RNA Virus Sequences (NIRVS) are enriched in piRNA clusters and coding sequences and have been proposed to constitute novel mosquito immune factors. However, given the abundance of NIRVS and their variable viral origin, their relative biological roles remain unexplored. Here we used an analytical approach that intersects computational, evolutionary and molecular methods to study the genomic landscape of mosquito NIRVS. We demonstrate that NIRVS are differentially distributed across mosquito genomes, with a core set of seemingly the oldest integrations with similarity to Rhabdoviruses. Additionally, we compare the polymorphisms of NIRVS with respect to that of fast and slow-evolving genes within the Ae. albopictus genome. Overall, NIRVS appear to be less polymorphic than slow-evolving genes, with differences depending on whether they occur in intergenic regions or in piRNA clusters. Finally, two NIRVS that map within the coding sequences of genes annotated as Rhabdovirus RNA-dependent RNA polymerase and the nucleocapsid-encoding gene, respectively, are highly polymorphic and are expressed, suggesting exaptation possibly to enhance the mosquito's antiviral responses. These results greatly advance our understanding of the complexity of the mosquito repeatome and the biology of viral integrations in mosquito genomes.

Molecular dissection of a natural transposable element invasion

The first tracking of the dynamics of a natural invasion by a transposable element (TE) provides unprecedented details on the establishment of host defense mechanisms against TEs. We captured a D. simulans population at an early stage of a P-element invasion and studied the spread of the TE in replicated experimentally evolving populations kept under hot and cold conditions. We analyzed the factors controlling the invasion by NGS, RNA-FISH, and gonadal dysgenesis assays. Under hot conditions, the P-element spread rapidly for 20 generations, but no further spread was noted later on. This plateauing of the invasion was mediated by the rapid emergence of P-element-specific piRNAs. Under cold conditions, we observed a lower expression of the P-element and a slower emergence of the piRNA defense, resulting in a three times slower invasion that continued beyond 40 generations. We conclude that the environment is a major factor determining the evolution of TEs in their host.

Targeted identification of TE insertions in a Drosophila genome through hemi-specific PCR

BACKGROUND

Transposable elements (TEs) are major components of eukaryotic genomes and drivers of genome evolution, producing intraspecific polymorphism and interspecific differences through mobilization and non-homologous recombination. TE insertion sites are often highly variable within species, creating a need for targeted genome re-sequencing (TGS) methods to identify TE insertion sites.

METHODS

We present a hemi-specific PCR approach for TGS of P-elements in Drosophila genomes on the Illumina platform. We also present a computational framework for identifying new insertions from TGS reads. Finally, we describe a new method for estimating the frequency of TE insertions from WGS data, which is based precise insertion sites provided by TGS annotations.

RESULTS

By comparing our results to TE annotations based on whole genome re-sequencing (WGS) data for the same Drosophilamelanogaster strain, we demonstrate that TGS is powerful for identifying true insertions, even in repeat-rich heterochromatic regions. We also demonstrate that TGS offers enhanced annotation of precise insertion sites, which facilitates estimation of TE insertion frequency.

CONCLUSIONS

TGS by hemi-specific PCR is a powerful approach for identifying TE insertions of particular TE families in species with a high-quality reference genome, at greatly reduced cost as compared to WGS. It may therefore be ideal for population genomic studies of particular TE families. Additionally, TGS and WGS can be used as complementary approaches, with TGS annotations identifying more annotated insertions with greater precision for a target TE family, and WGS data allowing for estimates of TE insertion frequencies, and a broader picture of the location of non-target TEs across the genome.