Tempo and Mode of Transposable Element Activity in Drosophila

Pathogen infection alters the gene expression landscape of transposable elements in Drosophila melanogaster

Transposable elements make up substantial proportions of eukaryotic genomes and many are thought to be remnants of ancient viral infections. Current research has begun to highlight the role transposable elements can play in the immune system response to infections. However, most of our knowledge about transposable element expression during infection is limited by the specific host and pathogen factors from each study, making it difficult to compare studies and develop broader patterns regarding the role of transposable elements during infection. Here, we use the tools and resources available in the model, Drosophila melanogaster, to analyze multiple gene expression datasets of flies subject to bacterial, fungal, and viral infections. We analyzed differences in pathogen species, host genotype, host tissue, and sex to understand how these factors impact transposable element expression during infection. Our results highlight both shared and unique transposable element expression patterns between pathogens and suggest a larger effect of pathogen factors over host factors for influencing transposable element expression.

Causes and Consequences of Varying Transposable Element Activity: An Evolutionary Perspective

Transposable elements (TEs) are genomic parasites found in nearly all eukaryotes, including humans. This evolutionary success of TEs is due to their replicative activity, involving insertion into new genomic locations. TE activity varies at multiple levels, from between taxa to within individuals. The rapidly accumulating evidence of the influence of TE activity on human health, as well as the rapid growth of new tools to study it, motivated an evaluation of what we know about TE activity thus far. Here, we discuss why TE activity varies, and the consequences of this variation, from an evolutionary perspective. By studying TE activity in nonhuman organisms in the context of evolutionary theories, we can shed light on the factors that affect TE activity. While the consequences of TE activity are usually deleterious, some have lasting evolutionary impacts by conferring benefits on the host or affecting other evolutionary processes.

DNA Transposons Favor De Novo Transcript Emergence Through Enrichment of Transcription Factor Binding Motifs

De novo genes emerge from noncoding regions of genomes via succession of mutations. Among others, such mutations activate transcription and create a new open reading frame (ORF). Although the mechanisms underlying ORF emergence are well documented, relatively little is known about the mechanisms enabling new transcription events. Yet, in many species a continuum between absent and very prominent transcription has been reported for essentially all regions of the genome. In this study, we searched for de novo transcripts by using newly assembled genomes and transcriptomes of seven inbred lines of Drosophila melanogaster, originating from six European and one African population. This setup allowed us to detect sample specific de novo transcripts, and compare them to their homologous nontranscribed regions in other samples, as well as genic and intergenic control sequences. We studied the association with transposable elements (TEs) and the enrichment of transcription factor motifs upstream of de novo emerged transcripts and compared them with regulatory elements. We found that de novo transcripts overlap with TEs more often than expected by chance. The emergence of new transcripts correlates with regions of high guanine-cytosine content and TE expression. Moreover, upstream regions of de novo transcripts are highly enriched with regulatory motifs. Such motifs are more enriched in new transcripts overlapping with TEs, particularly DNA TEs, and are more conserved upstream de novo transcripts than upstream their 'nontranscribed homologs'. Overall, our study demonstrates that TE insertion is important for transcript emergence, partly by introducing new regulatory motifs from DNA TE families.

Rapid evolution of piRNA clusters in the Drosophila melanogaster ovary

The piRNA pathway is a highly conserved mechanism to repress transposable element (TE) activity in the animal germline via a specialized class of small RNAs called piwi-interacting RNAs (piRNAs). piRNAs are produced from discrete genomic regions called piRNA clusters (piCs). Although the molecular processes by which piCs function are relatively well understood in Drosophila melanogaster, much less is known about the origin and evolution of piCs in this or any other species. To investigate piC origin and evolution, we use a population genomic approach to compare piC activity and sequence composition across eight geographically distant strains of D. melanogaster with high-quality long-read genome assemblies. We perform annotations of ovary piCs and genome-wide TE content in each strain. Our analysis uncovers extensive variation in piC activity across strains and signatures of rapid birth and death of piCs. Most TEs inferred to be recently active show an enrichment of insertions into old and large piCs, consistent with the previously proposed "trap" model of piC evolution. In contrast, a small subset of active LTR families is enriched for the formation of new piCs, suggesting that these TEs have higher proclivity to form piCs. Thus, our findings uncover processes leading to the origin of piCs. We propose that piC evolution begins with the emergence of piRNAs from individual insertions of a few select TE families prone to seed new piCs that subsequently expand by accretion of insertions from most other TE families during evolution to form larger "trap" clusters. Our study shows that TEs themselves are the major force driving the rapid evolution of piCs.

Genomes of historical specimens reveal multiple invasions of LTR retrotransposons in Drosophila melanogaster during the 19th century

Transposable element invasions have a profound impact on the evolution of genomes and phenotypes. It is thus an important open question how often such TE invasions occur. To address this question, we utilize the genomes of historical specimens, sampled about 200 y ago. We found that the LTR retrotransposons Blood, Opus, and 412 spread in Drosophila melanogaster in the 19th century. These invasions constitute second waves, as degraded fragments were found for all three TEs. The composition of Opus and 412, but not of Blood, shows a pronounced geographic heterogeneity, likely due to founder effects during the invasions. Finally, we identified species from the Drosophila simulans complex as the likely origin of the TEs. We show that in total, seven TE families invaded D. melanogaster during the last 200y, thereby increasing the genome size by up to 1.2Mbp. We suggest that this high rate of TE invasions was likely triggered by human activity. Based on the analysis of strains and specimens sampled at different times, we provide a detailed timeline of TE invasions, making D. melanogaster the first organism where the invasion history of TEs during the last two centuries could be inferred.

Orthoptera-TElib: a library of Orthoptera transposable elements for TE annotation

Transposable elements (TEs) are a major component of eukaryotic genomes and are present in almost all eukaryotic organisms. TEs are highly dynamic between and within species, which significantly affects the general applicability of the TE databases. Orthoptera is the only known group in the class Insecta with a significantly enlarged genome (0.93-21.48 Gb). When analyzing the large genome using the existing TE public database, the efficiency of TE annotation is not satisfactory. To address this limitation, it becomes imperative to continually update the available TE resource library and the need for an Orthoptera-specific library as more insect genomes are publicly available. Here, we used the complete genome data of 12 Orthoptera species to de novo annotate TEs, then manually re-annotate the unclassified TEs to construct a non-redundant Orthoptera-specific TE library: Orthoptera-TElib. Orthoptera-TElib contains 24,021 TE entries including the re-annotated results of 13,964 unknown TEs. The naming of TE entries in Orthoptera-TElib adopts the same naming as RepeatMasker and Dfam and is encoded as the three-level form of "level1/level2-level3". Orthoptera-TElib can be directly used as an input reference database and is compatible with mainstream repetitive sequence analysis software such as RepeatMasker and dnaPipeTE. When analyzing TEs of Orthoptera species, Orthoptera-TElib performs better TE annotation as compared to Dfam and Repbase regardless of using low-coverage sequencing or genome assembly data. The most improved TE annotation result is Angaracris rhodopa, which has increased from 7.89% of the genome to 53.28%. Finally, Orthoptera-TElib is stored in Sqlite3 for the convenience of data updates and user access.

The Complex Landscape of Structural Divergence Between the Drosophila pseudoobscura and D. persimilis Genomes

Structural genomic variants are key drivers of phenotypic evolution. They can span hundreds to millions of base pairs and can thus affect large numbers of genetic elements. Although structural variation is quite common within and between species, its characterization depends upon the quality of genome assemblies and the proportion of repetitive elements. Using new high-quality genome assemblies, we report a complex and previously hidden landscape of structural divergence between the genomes of Drosophila persimilis and D. pseudoobscura, two classic species in speciation research, and study the relationships among structural variants, transposable elements, and gene expression divergence. The new assemblies confirm the already known fixed inversion differences between these species. Consistent with previous studies showing higher levels of nucleotide divergence between fixed inversions relative to collinear regions of the genome, we also find a significant overrepresentation of INDELs inside the inversions. We find that transposable elements accumulate in regions with low levels of recombination, and spatial correlation analyses reveal a strong association between transposable elements and structural variants. We also report a strong association between differentially expressed (DE) genes and structural variants and an overrepresentation of DE genes inside the fixed chromosomal inversions that separate this species pair. Interestingly, species-specific structural variants are overrepresented in DE genes involved in neural development, spermatogenesis, and oocyte-to-embryo transition. Overall, our results highlight the association of transposable elements with structural variants and their importance in driving evolutionary divergence.

The impact of paramutations on the invasion dynamics of transposable elements

According to the prevailing view, the trap model, the activity of invading transposable elements (TEs) is greatly reduced when a TE copy jumps into a piRNA cluster, which triggers the emergence of piRNAs that silence the TE. One crucial component in the host defence are paramutations. Mediated by maternally deposited piRNAs, paramutations convert TE insertions into piRNA producing loci, thereby transforming selfish TEs into agents of the host defence. Despite this significant effect, the impact of paramutations on the dynamics of TE invasions remains unknown. To address this issue, we performed extensive forward simulations of TE invasions with piRNA clusters and paramutations. We found that paramutations significantly affect TE dynamics, by accelerating the silencing of TE invasions, reducing the number of insertions accumulating during the invasions and mitigating the fitness cost of TEs. We also demonstrate that piRNA production induced by paramutations, an epigenetically inherited trait, may be positively selected. Finally, we show that paramutations may account for three important open problems with the trap model. Firstly, paramutated TE insertions may compensate for the insufficient number of insertions in piRNA clusters observed in previous studies. Secondly, paramutations may explain the discrepancy between the observed and the expected abundance of different TE families in Drosophila melanogaster. Thirdly, piRNA clusters may be crucial to trigger the host defence, but paramutations render the clusters dispensable once the defence has been established. This could account for the lack of TE activation when three major piRNA clusters were deleted in a previous study.

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.

The genomic distribution of transposable elements is driven by spatially variable purifying selection

It is widely accepted that the genomic distribution of transposable elements (TEs) mainly reflects the outcome of purifying selection and insertion bias (1). Nevertheless, the relative importance of these two evolutionary forces could not be tested thoroughly. Here, we introduce an experimental system, which allows separating purifying selection from TE insertion bias. We used experimental evolution to study the TE insertion patterns in Drosophila simulans founder populations harboring 1040 insertions of an active P-element. After 10 generations at a large population size, we detected strong selection against P-element insertions. The exception were P-element insertions in genomic regions for which a strong insertion bias has been proposed (2-4). Because recurrent P-element insertions cannot explain this pattern, we conclude that purifying selection, with variable strength along the chromosomes, is the major determinant of the genomic distribution of P-elements. Genomic regions with relaxed purifying selection against P-element insertions exhibit normal levels of purifying selection against base substitutions. This suggests that different types of purifying selection operate on base substitutions and P-element insertions. Our results highlight the power of experimental evolution to understand basic evolutionary processes, which are difficult to infer from patterns of natural variation alone.

A Quantitative, Genome-Wide Analysis in Drosophila Reveals Transposable Elements' Influence on Gene Expression Is Species-Specific

Transposable elements (TEs) are parasite DNA sequences that are able to move and multiply along the chromosomes of all genomes. They can be controlled by the host through the targeting of silencing epigenetic marks, which may affect the chromatin structure of neighboring sequences, including genes. In this study, we used transcriptomic and epigenomic high-throughput data produced from ovarian samples of several Drosophila melanogaster and Drosophila simulans wild-type strains, in order to finely quantify the influence of TE insertions on gene RNA levels and histone marks (H3K9me3 and H3K4me3). Our results reveal a stronger epigenetic effect of TEs on ortholog genes in D. simulans compared with D. melanogaster. At the same time, we uncover a larger contribution of TEs to gene H3K9me3 variance within genomes in D. melanogaster, which is evidenced by a stronger correlation of TE numbers around genes with the levels of this chromatin mark in D. melanogaster. Overall, this work contributes to the understanding of species-specific influence of TEs within genomes. It provides a new light on the considerable natural variability provided by TEs, which may be associated with contrasted adaptive and evolutionary potentials.

Population genomics reveals mechanisms and dynamics of de novo expressed open reading frame emergence in Drosophila melanogaster

Novel genes are essential for evolutionary innovations and differ substantially even between closely related species. Recently, multiple studies across many taxa showed that some novel genes arise de novo, that is, from previously noncoding DNA. To characterize the underlying mutations that allowed de novo gene emergence and their order of occurrence, homologous regions must be detected within noncoding sequences in closely related sister genomes. So far, most studies do not detect noncoding homologs of de novo genes because of incomplete assemblies and annotations, and long evolutionary distances separating genomes. Here, we overcome these issues by searching for de novo expressed open reading frames (neORFs), the not-yet fixed precursors of de novo genes that emerged within a single species. We sequenced and assembled genomes with long-read technology and the corresponding transcriptomes from inbred lines of Drosophila melanogaster, derived from seven geographically diverse populations. We found line-specific neORFs in abundance but few neORFs shared by lines, suggesting a rapid turnover. Gain and loss of transcription is more frequent than the creation of ORFs, for example, by forming new start and stop codons. Consequently, the gain of ORFs becomes rate limiting and is frequently the initial step in neORFs emergence. Furthermore, transposable elements (TEs) are major drivers for intragenomic duplications of neORFs, yet TE insertions are less important for the emergence of neORFs. However, highly mutable genomic regions around TEs provide new features that enable gene birth. In conclusion, neORFs have a high birth-death rate, are rapidly purged, but surviving neORFs spread neutrally through populations and within genomes.

Rapid evolution of piRNA clusters in the Drosophila melanogaster ovary

Animal genomes are parasitized by a horde of transposable elements (TEs) whose mutagenic activity can have catastrophic consequences. The piRNA pathway is a conserved mechanism to repress TE activity in the germline via a specialized class of small RNAs associated with effector Piwi proteins called piwi-associated RNAs (piRNAs). piRNAs are produced from discrete genomic regions called piRNA clusters (piCs). While piCs are generally enriched for TE sequences and the molecular processes by which they are transcribed and regulated are relatively well understood in Drosophila melanogaster, much less is known about the origin and evolution of piCs in this or any other species. To investigate piC evolution, we use a population genomics approach to compare piC activity and sequence composition across 8 geographically distant strains of D. melanogaster with high quality long-read genome assemblies. We perform extensive annotations of ovary piCs and TE content in each strain and test predictions of two proposed models of piC evolution. The 'de novo' model posits that individual TE insertions can spontaneously attain the status of a small piC to generate piRNAs silencing the entire TE family. The 'trap' model envisions large and evolutionary stable genomic clusters where TEs tend to accumulate and serves as a long-term "memory" of ancient TE invasions and produce a great variety of piRNAs protecting against related TEs entering the genome. It remains unclear which model best describes the evolution of piCs. Our analysis uncovers extensive variation in piC activity across strains and signatures of rapid birth and death of piCs in natural populations. Most TE families inferred to be recently or currently active show an enrichment of strain-specific insertions into large piCs, consistent with the trap model. By contrast, only a small subset of active LTR retrotransposon families is enriched for the formation of strain-specific piCs, suggesting that these families have an inherent proclivity to form de novo piCs. Thus, our findings support aspects of both 'de novo' and 'trap' models of piC evolution. We propose that these two models represent two extreme stages along an evolutionary continuum, which begins with the emergence of piCs de novo from a few specific LTR retrotransposon insertions that subsequently expand by accretion of other TE insertions during evolution to form larger 'trap' clusters. Our study shows that piCs are evolutionarily labile and that TEs themselves are the major force driving the formation and evolution of piCs.

Epigenetic and chromosomal features drive transposon insertion in Drosophila melanogaster

Transposons are mobile genetic elements prevalent in the genomes of most species. The distribution of transposons within a genome reflects the actions of two opposing processes: initial insertion site selection, and selective pressure from the host. By analyzing whole-genome sequencing data from transposon-activated Drosophila melanogaster, we identified 43 316 de novo and 237 germline insertions from four long-terminal-repeat (LTR) transposons, one LINE transposon (I-element), and one DNA transposon (P-element). We found that all transposon types favored insertion into promoters de novo, but otherwise displayed distinct insertion patterns. De novo and germline P-element insertions preferred replication origins, often landing in a narrow region around transcription start sites and in regions of high chromatin accessibility. De novo LTR transposon insertions preferred regions with high H3K36me3, promoters and exons of active genes; within genes, LTR insertion frequency correlated with gene expression. De novo I-element insertion density increased with distance from the centromere. Germline I-element and LTR transposon insertions were depleted in promoters and exons, suggesting strong selective pressure to remove transposons from functional elements. Transposon movement is associated with genome evolution and disease; therefore, our results can improve our understanding of genome and disease biology.

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.

Transposable element expansion and low-level piRNA silencing in grasshoppers may cause genome gigantism

Background

Transposable elements (TEs) have been likened to parasites in the genome that reproduce and move ceaselessly in the host, continuously enlarging the host genome. However, the Piwi-interacting RNA (piRNA) pathway defends animal genomes against the harmful consequences of TE invasion by imposing small-RNA-mediated silencing. Here we compare the TE activity of two grasshopper species with different genome sizes in Acrididae (Locusta migratoria manilensis♀1C = 6.60 pg, Angaracris rhodopa♀1C = 16.36 pg) to ascertain the influence of piRNAs.

Results

We discovered that repetitive sequences accounted for 74.56% of the genome in A. rhodopa, more than 56.83% in L. migratoria, and the large-genome grasshopper contained a higher TEs proportions. The comparative analysis revealed that 41 TEs (copy number > 500) were shared in both species. The two species exhibited distinct “landscapes” of TE divergence. The TEs outbreaks in the small-genome grasshopper occurred at more ancient times, while the large-genome grasshopper maintains active transposition events in the recent past. Evolutionary history studies on TEs suggest that TEs may be subject to different dynamics and resistances in these two species. We found that TE transcript abundance was higher in the large-genome grasshopper and the TE-derived piRNAs abundance was lower than in the small-genome grasshopper. In addition, we found that the piRNA methylase HENMT, which is underexpressed in the large-genome grasshopper, impedes the piRNA silencing to a lower level.

Conclusions

Our study revealed that the abundance of piRNAs is lower in the gigantic genome grasshopper than in the small genome grasshopper. In addition, the key gene HENMT in the piRNA biogenesis pathway (Ping-Pong cycle) in the gigantic genome grasshopper is underexpressed. We hypothesize that low-level piRNA silencing unbalances the original positive correlation between TEs and piRNAs, and triggers TEs to proliferate out of control, which may be one of the reasons for the gigantism of grasshopper genomes.

A novel transposable element-mediated mechanism causes antiviral resistance in Drosophila through truncating the Veneno protein

Hosts are continually selected to evolve new defenses against an ever-changing array of pathogens. To understand this process, we examined the genetic basis of resistance to the Drosophila A virus in Drosophila melanogaster. In a natural population, we identified a polymorphic transposable element (TE) insertion that was associated with an ∼19,000-fold reduction in viral titers, allowing flies to largely escape the harmful effects of infection by this virulent pathogen. The insertion occurs in the protein-coding sequence of the gene Veneno, which encodes a Tudor domain protein. By mutating Veneno with CRISPR-Cas9 in flies and expressing it in cultured cells, we show that the ancestral allele of the gene has no effect on viral replication. Instead, the TE insertion is a gain-of-function mutation that creates a gene encoding a novel resistance factor. Viral titers remained reduced when we deleted the TE sequence from the transcript, indicating that resistance results from the TE truncating the Veneno protein. This is a novel mechanism of virus resistance and a new way by which TEs can contribute to adaptation.

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.

Epistatic selection on a selfish Segregation Distorter supergene - drive, recombination, and genetic load

Meiotic drive supergenes are complexes of alleles at linked loci that together subvert Mendelian segregation resulting in preferential transmission. In males, the most common mechanism of drive involves the disruption of sperm bearing one of a pair of alternative alleles. While at least two loci are important for male drive-the driver and the target-linked modifiers can enhance drive, creating selection pressure to suppress recombination. In this work, we investigate the evolution and genomic consequences of an autosomal, multilocus, male meiotic drive system, Segregation Distorter (SD) in the fruit fly, Drosophila melanogaster. In African populations, the predominant SD chromosome variant, SD-Mal, is characterized by two overlapping, paracentric inversions on chromosome arm 2R and nearly perfect (~100%) transmission. We study the SD-Mal system in detail, exploring its components, chromosomal structure, and evolutionary history. Our findings reveal a recent chromosome-scale selective sweep mediated by strong epistatic selection for haplotypes carrying Sd, the main driving allele, and one or more factors within the double inversion. While most SD-Mal chromosomes are homozygous lethal, SD-Mal haplotypes can recombine with other, complementing haplotypes via crossing over, and with wildtype chromosomes via gene conversion. SD-Mal chromosomes have nevertheless accumulated lethal mutations, excess non-synonymous mutations, and excess transposable element insertions. Therefore, SD-Mal haplotypes evolve as a small, semi-isolated subpopulation with a history of strong selection. These results may explain the evolutionary turnover of SD haplotypes in different populations around the world and have implications for supergene evolution broadly.

Synergistic epistasis of the deleterious effects of transposable elements

The replicative nature and generally deleterious effects of transposable elements (TEs) raise an outstanding question about how TE copy number is stably contained in host populations. Classic theoretical analyses predict that, when the decline in fitness due to each additional TE insertion is greater than linear, or when there is synergistic epistasis, selection against TEs can result in a stable equilibrium of TE copy number. While several mechanisms are predicted to yield synergistic deleterious effects of TEs, we lack empirical investigations of the presence of such epistatic interactions. Purifying selection with synergistic epistasis generates repulsion linkage between deleterious alleles. We investigated this population genetic signal in the likely ancestral Drosophila melanogaster population and found evidence supporting the presence of synergistic epistasis among TE insertions, especially TEs expected to exert large fitness impacts. Even though synergistic epistasis of TEs has been predicted to arise through ectopic recombination and TE-mediated epigenetic silencing mechanisms, we only found mixed support for the associated predictions. We observed signals of synergistic epistasis for a large number of TE families, which is consistent with the expectation that such epistatic interaction mainly happens among copies of the same family. Curiously, significant repulsion linkage was also found among TE insertions from different families, suggesting the possibility that synergism of TEs' deleterious fitness effects could arise above the family level and through mechanisms similar to those of simple mutations. Our findings set the stage for investigating the prevalence and importance of epistatic interactions in the evolutionary dynamics of TEs.

Transposable Elements: Distribution, Polymorphism, and Climate Adaptation in Populus

Transposable elements (TEs) are a class of mobile genetic elements that make effects on shaping rapid phenotypic traits of adaptive significance. TE insertions are usually related to transcription changes of nearby genes, and thus may be subjected to purifying selection. Based on the available genome resources of Populus, we found that the composition of Helitron DNA family were highly variable and could directly influence the transcription of nearby gene expression, which are involving in stress-responsive, programmed cell death, and apoptosis pathway. Next, we indicated TEs are highly enriched in Populus trichocarpa compared with three other congeneric poplar species, especially located at untranslated regions (3'UTRs and 5'UTRs) and Helitron transposons, particularly 24-nt siRNA-targeted, are significantly associated with reduced gene expression. Additionally, we scanned a representative resequenced Populus tomentosa population, and identified 9,680 polymorphic TEs loci. More importantly, we identified a Helitron transposon located at the 3'UTR, which could reduce WRKY18 expression level. Our results highlight the importance of TE insertion events, which could regulate gene expression and drive adaptive phenotypic variation in Populus.

Patterns of piRNA Regulation in Drosophila Revealed through Transposable Element Clade Inference

Transposable elements (TEs) are self-replicating "genetic parasites" ubiquitous to eukaryotic genomes. In addition to conflict between TEs and their host genomes, TEs of the same family are in competition with each other. They compete for the same genomic niches while experiencing the same regime of copy-number selection. This suggests that competition among TEs may favor the emergence of new variants that can outcompete their ancestral forms. To investigate the sequence evolution of TEs, we developed a method to infer clades: collections of TEs that share SNP variants and represent distinct TE family lineages. We applied this method to a panel of 85 Drosophila melanogaster genomes and found that the genetic variation of several TE families shows significant population structure that arises from the population-specific expansions of single clades. We used population genetic theory to classify these clades into younger versus older clades and found that younger clades are associated with a greater abundance of sense and antisense piRNAs per copy than older ones. Further, we find that the abundance of younger, but not older clades, is positively correlated with antisense piRNA production, suggesting a general pattern where hosts preferentially produce antisense piRNAs from recently active TE variants. Together these findings suggest a pattern whereby new TE variants arise by mutation and then increase in copy number, followed by the host producing antisense piRNAs that may be used to silence these emerging variants.

Transposable element variants and their potential adaptive impact in urban populations of the malaria vector Anopheles coluzzii

Anopheles coluzzii is one of the primary vectors of human malaria in sub-Saharan Africa. Recently, it has spread into the main cities of Central Africa threatening vector control programs. The adaptation of An. coluzzii to urban environments partly results from an increased tolerance to organic pollution and insecticides. Some of the molecular mechanisms for ecological adaptation are known, but the role of transposable elements (TEs) in the adaptive processes of this species has not been studied yet. As a first step toward assessing the role of TEs in rapid urban adaptation, we sequenced using long reads six An. coluzzii genomes from natural breeding sites in two major Central Africa cities. We de novo annotated TEs in these genomes and in an additional high-quality An. coluzzii genome, and we identified 64 new TE families. TEs were nonrandomly distributed throughout the genome with significant differences in the number of insertions of several superfamilies across the studied genomes. We identified seven putatively active families with insertions near genes with functions related to vectorial capacity, and several TEs that may provide promoter and transcription factor binding sites to insecticide resistance and immune-related genes. Overall, the analysis of multiple high-quality genomes allowed us to generate the most comprehensive TE annotation in this species to date and identify several TE insertions that could potentially impact both genome architecture and the regulation of functionally relevant genes. These results provide a basis for future studies of the impact of TEs on the biology of An. coluzzii.

Transposable elements expression in Rhinella marina (cane toad) specimens submitted to immune and stress challenge

Transposable elements (TEs) are important components of eukaryotic genomes and compose around 30% of the genome of Rhinella marina, an invasive toad species. Considering the possible role of TEs in the adaptation of populations, we have analyzed the expression of TEs in publicly available spleen tissue transcriptomic data generated for this species after immune and stress challenge. By analyzing the transcriptome assembly, we detected a high number of TE segments. Moreover, some distinct TE families were differentially expressed in some conditions. Our result shows that several TEs are capable of being transcribed in R. marina and they could help to generate a rapid response of specimens to the environment. Also, we can suggest that these TEs could be activated in the germinative cells as well producing variability to be selected and shaped by the evolutionary processes behind the success of this invasive species. Thus, the TEs are important targets for investigation in the context of R. marina adaptation.

Repetitive genomic regions and the inference of demographic history

Inference of demographic histories using whole-genome datasets has provided insights into diversification, adaptation, hybridization, and plant-pathogen interactions, and stimulated debate on the impact of anthropogenic interventions and past climate on species demography. However, the impact of repetitive genomic regions on these inferences has mostly been ignored by masking of repeats. We use the Populus trichocarpa genome (Pop_tri_v3) to show that masking of repeat regions leads to lower estimates of effective population size (Ne) in the distant past in contrast to an increase in Ne estimates in recent times. However, in human datasets, masking of repeats resulted in lower estimates of Ne at all time points. We demonstrate that repeats affect demographic inferences using diverse methods like PSMC, MSMC, SMC++, and the Stairway plot. Our genomic analysis revealed that the biases in Ne estimates were dependent on the repeat class type and its abundance in each atomic interval. Notably, we observed a weak, yet consistently significant negative correlation between the repeat abundance of an atomic interval and the Ne estimates for that interval, which potentially reflects the recombination rate variation within the genome. The rationale for the masking of repeats has been that variants identified within these regions are erroneous. We find that polymorphisms in some repeat classes occur in callable regions and reflect reliable coalescence histories (e.g., LTR Gypsy, LTR Copia). The current demography inference methods do not handle repeats explicitly, and hence the effect of individual repeat classes needs careful consideration in comparative analysis. Deciphering the repeat demographic histories might provide a clear understanding of the processes involved in repeat accumulation.

Diverse Eukaryotic CGG-Binding Proteins Produced by Independent Domestications of hAT Transposons

The human transcription factor (TF) CGGBP1 (CGG-binding protein) is conserved only in amniotes and is believed to derive from the zf-BED and Hermes transposase DNA-binding domains (DBDs) of a hAT DNA transposon. Here, we show that sequence-specific DNA-binding proteins with this bipartite domain structure have resulted from dozens of independent hAT domestications in different eukaryotic lineages. CGGBPs display a wide range of sequence specificity, usually including preferences for CGG or CGC trinucleotides, whereas some bind AT-rich motifs. The CGGBPs are almost entirely nonsyntenic, and their protein sequences, DNA-binding motifs, and patterns of presence or absence in genomes are uncharacteristic of ancestry via speciation. At least eight CGGBPs in the coelacanth Latimeria chalumnae bind distinct motifs, and the expression of the corresponding genes varies considerably across tissues, suggesting tissue-restricted function.

piRNA and Transposon Dynamics in Drosophila: A Female Story

The germlines of metazoans contain transposable elements (TEs) causing genetic instability and affecting fitness. To protect the germline from TE activity, gonads of metazoans produce TE-derived PIWI-interacting RNAs (piRNAs) that silence TE expression. In Drosophila, our understanding of piRNA biogenesis is mainly based on studies of the Drosophila melanogaster female germline. However, it is not known whether piRNA functions are also important in the male germline or whether and how piRNAs are affected by the global genomic context. To address these questions, we compared genome sequences, transcriptomes, and small RNA libraries extracted from entire testes and ovaries of two sister species: D. melanogaster and Drosophila simulans. We found that most TE-derived piRNAs were produced in ovaries and that piRNA pathway genes were strongly overexpressed in ovaries compared with testes, indicating that the silencing of TEs by the piRNA pathway mainly took place in the female germline. To study the relationship between host piRNAs and TE landscape, we analyzed TE genomic features and how they correlate with piRNA production in the two species. In D. melanogaster, we found that TE-derived piRNAs target recently active TEs. In contrast, although Drosophila simulans TEs do not display any features of recent activity, the host still intensively produced silencing piRNAs targeting old TE relics. Together, our results show that the piRNA silencing response mainly takes place in Drosophila ovaries and indicate that the host piRNA response is implemented following a burst of TE activity and could persist long after the extinction of active TE families.

Monitoring insect transposable elements in large double-stranded DNA viruses reveals host-to-virus and virus-to-virus transposition

The mechanisms by which transposable elements (TEs) can be horizontally transferred between animals are unknown, but viruses are possible candidate vectors. Here, we surveyed the presence of host-derived TEs in viral genomes in 35 deep sequencing data sets produced from 11 host–virus systems, encompassing nine arthropod host species (five lepidopterans, two dipterans, and two crustaceans) and six different double-stranded (ds) DNA viruses (four baculoviruses and two iridoviruses). We found evidence of viral-borne TEs in 14 data sets, with frequencies of viral genomes carrying a TE ranging from 0.01% to 26.33% for baculoviruses and from 0.45% to 7.36% for iridoviruses. The analysis of viral populations separated by a single replication cycle revealed that viral-borne TEs originating from an initial host species can be retrieved after viral replication in another host species, sometimes at higher frequencies. Furthermore, we detected a strong increase in the number of integrations in a viral population for a TE absent from the hosts’ genomes, indicating that this TE has undergone intense transposition within the viral population. Finally, we provide evidence that many TEs found integrated in viral genomes (15/41) have been horizontally transferred in insects. Altogether, our results indicate that multiple large dsDNA viruses have the capacity to shuttle TEs in insects and they underline the potential of viruses to act as vectors of horizontal transfer of TEs. Furthermore, the finding that TEs can transpose between viral genomes of a viral species sets viruses as possible new niches in which TEs can persist and evolve.

Rapid Evolution of Autosomal Binding Sites of the Dosage Compensation Complex in Drosophila melanogaster and Its Association With Transcription Divergence

How pleiotropy influences evolution of protein sequence remains unclear. The male-specific lethal (MSL) complex in Drosophila mediates dosage compensation by 2-fold upregulation of the X chromosome in males. Nevertheless, several MSL proteins also bind autosomes and likely perform functions not related to dosage compensation. Here, we study the evolution of MOF, MSL1, and MSL2 biding sites in Drosophila melanogaster and its close relative Drosophila simulans. We found pervasive expansion of the MSL binding sites in D. melanogaster, particularly on autosomes. The majority of these newly-bound regions are unlikely to function in dosage compensation and associated with an increase in expression divergence between D. melanogaster and D. simulans. While dosage-compensation related sites show clear signatures of adaptive evolution, these signatures are even more marked among autosomal regions. Our study points to an intriguing avenue of investigation of pleiotropy as a mechanism promoting rapid protein sequence evolution.

piRNA-mediated gene regulation and adaptation to sex-specific transposon expression in D. melanogaster male germline

Small noncoding piRNAs act as sequence-specific guides to repress complementary targets in Metazoa. Prior studies in Drosophila ovaries have demonstrated the function of the piRNA pathway in transposon silencing and therefore genome defense. However, the ability of the piRNA program to respond to different transposon landscapes and the role of piRNAs in regulating host gene expression remain poorly understood. Here, we comprehensively analyzed piRNA expression and defined the repertoire of their targets in Drosophila melanogaster testes. Comparison of piRNA programs between sexes revealed sexual dimorphism in piRNA programs that parallel sex-specific transposon expression. Using a novel bioinformatic pipeline, we identified new piRNA clusters and established complex satellites as dual-strand piRNA clusters. While sharing most piRNA clusters, the two sexes employ them differentially to combat the sex-specific transposon landscape. We found two piRNA clusters that produce piRNAs antisense to four host genes in testis, including CG12717/pirate, a SUMO protease gene. piRNAs encoded on the Y chromosome silence pirate, but not its paralog, to exert sex- and paralog-specific gene regulation. Interestingly, pirate is targeted by endogenous siRNAs in a sibling species, Drosophila mauritiana, suggesting distinct but related silencing strategies invented in recent evolution to regulate a conserved protein-coding gene.

The worldwide invasion of Drosophila suzukii is accompanied by a large increase of transposable element load and a small number of putatively adaptive insertions

Transposable Elements (TEs) are ubiquitous and mobile repeated sequences. They are major determinants of host fitness. Here, we characterized the TE content of the spotted wing fly Drosophila suzukii. Using a recently improved genome assembly, we reconstructed TE sequences de novo, and found that TEs occupy 47% of the genome and are mostly located in gene poor regions. The majority of TE insertions segregate at low frequencies, indicating a recent and probably ongoing TE activity. To explore TE dynamics in the context of biological invasions, we studied variation of TE abundance in genomic data from 16 invasive and six native populations of D. suzukii. We found a large increase of the TE load in invasive populations correlated with a reduced Watterson estimate of genetic diversity θ̂w a proxy of effective population size. We did not find any correlation between TE contents and bioclimatic variables, indicating a minor effect of environmentally induced TE activity. A genome-wide association study revealed that ca. 2,000 genomic regions are associated with TE abundance. We did not find, however, any evidence in such regions of an enrichment for genes known to interact with TE activity (e.g. transcription factor encoding genes or genes of the piRNA pathway). Finally, the study of TE insertion frequencies revealed 15 putatively adaptive TE insertions, six of them being likely associated with the recent invasion history of the species.

Rapid divergence of the male reproductive proteins in the Drosophila dunni group and implications for postmating incompatibilities between species

Proteins involved in post-copulatory interactions between males and females are among the fastest evolving genes in many species, usually attributed to their involvement in reproductive conflict. As a result, these proteins are thought to often be involved in the formation of postmating-prezygotic incompatibilities between species. The Drosophila dunni subgroup consists of a dozen recently diverged species found across the Caribbean islands with varying levels of hybrid incompatibility. We performed experimental crosses between species in the dunni group and see some evidence of hybrid incompatibilities. We also find evidence of reduced survival following hybrid mating, likely due to postmating-prezygotic incompatibilities. We assessed rates of evolution between these species genomes and find evidence of rapid evolution and divergence of some reproductive proteins, specifically the seminal fluid proteins. This work suggests the rapid evolution of seminal fluid proteins may be associated with postmating-prezygotic isolation, which acts as a barrier for gene flow between even the most closely related species.

Tirant Stealthily Invaded Natural Drosophila melanogaster Populations during the Last Century

It was long thought that solely three different transposable elements (TEs)-the I-element, the P-element, and hobo-invaded natural Drosophila melanogaster populations within the last century. By sequencing the "living fossils" of Drosophila research, that is, D. melanogaster strains sampled from natural populations at different time points, we show that a fourth TE, Tirant, invaded D. melanogaster populations during the past century. Tirant likely spread in D. melanogaster populations around 1938, followed by the I-element, hobo, and, lastly, the P-element. In addition to the recent insertions of the canonical Tirant, D. melanogaster strains harbor degraded Tirant sequences in the heterochromatin which are likely due to an ancient invasion, likely predating the split of D. melanogaster and D. simulans. These degraded insertions produce distinct piRNAs that were unable to prevent the novel Tirant invasion. In contrast to the I-element, P-element, and hobo, we did not find that Tirant induces any hybrid dysgenesis symptoms. This absence of apparent phenotypic effects may explain the late discovery of the Tirant invasion. Recent Tirant insertions were found in all investigated natural populations. Populations from Tasmania carry distinct Tirant sequences, likely due to a founder effect. By investigating the TE composition of natural populations and strains sampled at different time points, insertion site polymorphisms, piRNAs, and phenotypic effects, we provide a comprehensive study of a natural TE invasion.

The challenges of predicting transposable element activity in hybrids

Transposable elements (TEs) are ubiquitous mobile genetic elements that hold both disruptive and adaptive potential for species. It has long been postulated that their activity may be triggered by hybridization, a hypothesis that received mixed support from studies in various species. While host defense mechanisms against TEs are being elucidated, the increasing volume of genomic data and bioinformatic tools specialized in TE detection enable in-depth characterization of TEs at the levels of species and populations. Here, I borrow elements from the genome ecology theory to illustrate how knowledge of the diversity of TEs and host defense mechanisms may help predict the activity of TEs in the face of hybridization, and how current limitations make this task especially challenging.

Evolution of genome structure in the Drosophila simulans species complex

The rapid evolution of repetitive DNA sequences, including satellite DNA, tandem duplications, and transposable elements, underlies phenotypic evolution and contributes to hybrid incompatibilities between species. However, repetitive genomic regions are fragmented and misassembled in most contemporary genome assemblies. We generated highly contiguous de novo reference genomes for the Drosophila simulans species complex (D. simulans, D. mauritiana, and D. sechellia), which speciated ∼250,000 yr ago. Our assemblies are comparable in contiguity and accuracy to the current D. melanogaster genome, allowing us to directly compare repetitive sequences between these four species. We find that at least 15% of the D. simulans complex species genomes fail to align uniquely to D. melanogaster owing to structural divergence-twice the number of single-nucleotide substitutions. We also find rapid turnover of satellite DNA and extensive structural divergence in heterochromatic regions, whereas the euchromatic gene content is mostly conserved. Despite the overall preservation of gene synteny, euchromatin in each species has been shaped by clade- and species-specific inversions, transposable elements, expansions and contractions of satellite and tRNA tandem arrays, and gene duplications. We also find rapid divergence among Y-linked genes, including copy number variation and recent gene duplications from autosomes. Our assemblies provide a valuable resource for studying genome evolution and its consequences for phenotypic evolution in these genetic model species.

Transposable Element Landscape in Drosophila Populations Selected for Longevity

Transposable elements (TEs) inflict numerous negative effects on health and fitness as they replicate by integrating into new regions of the host genome. Even though organisms employ powerful mechanisms to demobilize TEs, transposons gradually lose repression during aging. The rising TE activity causes genomic instability and was implicated in age-dependent neurodegenerative diseases, inflammation, and the determination of lifespan. It is therefore conceivable that long-lived individuals have improved TE silencing mechanisms resulting in reduced TE expression relative to their shorter-lived counterparts and fewer genomic insertions. Here, we test this hypothesis by performing the first genome-wide analysis of TE insertions and expression in populations of Drosophila melanogaster selected for longevity through late-life reproduction for 50–170 generations from four independent studies. Contrary to our expectation, TE families were generally more abundant in long-lived populations compared with nonselected controls. Although simulations showed that this was not expected under neutrality, we found little evidence for selection driving TE abundance differences. Additional RNA-seq analysis revealed a tendency for reducing TE expression in selected populations, which might be more important for lifespan than regulating genomic insertions. We further find limited evidence of parallel selection on genes related to TE regulation and transposition. However, telomeric TEs were genomically and transcriptionally more abundant in long-lived flies, suggesting improved telomere maintenance as a promising TE-mediated mechanism for prolonging lifespan. Our results provide a novel viewpoint indicating that reproduction at old age increases the opportunity of TEs to be passed on to the next generation with little impact on longevity.

Population Genomics on the Fly: Recent Advances in Drosophila

Drosophila melanogaster, a small dipteran of African origin, represents one of the best-studied model organisms. Early work in this system has uniquely shed light on the basic principles of genetics and resulted in a versatile collection of genetic tools that allow to uncover mechanistic links between genotype and phenotype. Moreover, given its worldwide distribution in diverse habitats and its moderate genome-size, Drosophila has proven very powerful for population genetics inference and was one of the first eukaryotes whose genome was fully sequenced. In this book chapter, we provide a brief historical overview of research in Drosophila and then focus on recent advances during the genomic era. After describing different types and sources of genomic data, we discuss mechanisms of neutral evolution including the demographic history of Drosophila and the effects of recombination and biased gene conversion. Then, we review recent advances in detecting genome-wide signals of selection, such as soft and hard selective sweeps. We further provide a brief introduction to background selection, selection of noncoding DNA and codon usage and focus on the role of structural variants, such as transposable elements and chromosomal inversions, during the adaptive process. Finally, we discuss how genomic data helps to dissect neutral and adaptive evolutionary mechanisms that shape genetic and phenotypic variation in natural populations along environmental gradients. In summary, this book chapter serves as a starting point to Drosophila population genomics and provides an introduction to the system and an overview to data sources, important population genetic concepts and recent advances in the field.

Transposable Elements and the Evolution of Insects

Insects are major contributors to our understanding of the interaction between transposable elements (TEs) and their hosts, owing to seminal discoveries, as well as to the growing number of sequenced insect genomes and population genomics and functional studies. Insect TE landscapes are highly variable both within and across insect orders, although phylogenetic relatedness appears to correlate with similarity in insect TE content. This correlation is unlikely to be solely due to inheritance of TEs from shared ancestors and may partly reflect preferential horizontal transfer of TEs between closely related species. The influence of insect traits on TE landscapes, however, remains unclear. Recent findings indicate that, in addition to being involved in insect adaptations and aging, TEs are seemingly at the cornerstone of insect antiviral immunity. Thus, TEs are emerging as essential insect symbionts that may have deleterious or beneficial consequences on their hosts, depending on context.

Degradation of the Repetitive Genomic Landscape in a Close Relative of Caenorhabditis elegans

The abundance, diversity, and genomic distribution of repetitive elements is highly variable among species. These patterns are thought to be driven in part by reproductive mode and the interaction of selection and recombination, and recombination rates typically vary by chromosomal position. In the nematode Caenorhabditis elegans, repetitive elements are enriched at chromosome arms and depleted on centers, and this mirrors the chromosomal distributions of other genomic features such as recombination rate. How conserved is this genomic landscape of repeats, and what evolutionary forces maintain it? To address this, we compared the genomic organization of repetitive elements across five Caenorhabditis species with chromosome-level assemblies. As previously reported, repeat content is enriched on chromosome arms in most Caenorhabditis species, and no obvious patterns of repeat content associated with reproductive mode were observed. However, the fig-associated C. inopinata has experienced repetitive element expansion and reveals no association of global repeat density with chromosome position. Patterns of repeat superfamily specific distributions reveal this global pattern is driven largely by a few repeat superfamilies that in C. inopinata have expanded in number and have weak associations with chromosome position. Additionally, 15% of predicted protein-coding genes in C. inopinata align to transposon-related proteins. When these are excluded, C. inopinata has no enrichment of genes in chromosome centers, in contrast to its close relatives who all have such clusters. Forward evolutionary simulations reveal that chromosomal heterogeneity in recombination rate alone can generate structured repetitive genomic landscapes when insertions are weakly deleterious, whereas chromosomal heterogeneity in the fitness effects of transposon insertion can promote such landscapes across a variety of evolutionary scenarios. Thus, patterns of gene density along chromosomes likely contribute to global repetitive landscapes in this group, although other historical or genomic factors are needed to explain the idiosyncrasy of genomic organization of various transposable element taxa within C. inopinata. Taken together, these results highlight the power of comparative genomics and evolutionary simulations in testing hypotheses regarding the causes of genome organization.

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.

Transposable elements in Drosophila

Drosophila has been studied as a biological model for many years and many discoveries in biology rely on this species. Research on transposable elements (TEs) is not an exception. Drosophila has contributed significantly to our knowledge on the mechanisms of transposition and their regulation, but above all, it was one of the first organisms on which genetic and genomic studies of populations were done. In this review article, in a very broad way, we will approach the TEs of Drosophila with a historical hindsight as well as recent discoveries in the field.

Double-edged sword: The evolutionary consequences of the epigenetic silencing of transposable elements

Transposable elements (TEs) are genomic parasites that selfishly replicate at the expense of host fitness. Fifty years of evolutionary studies of TEs have concentrated on the deleterious genetic effects of TEs, such as their effects on disrupting genes and regulatory sequences. However, a flurry of recent work suggests that there is another important source of TEs' harmful effects-epigenetic silencing. Host genomes typically silence TEs by the deposition of repressive epigenetic marks. While this silencing reduces the selfish replication of TEs and should benefit hosts, a picture is emerging that the epigenetic silencing of TEs triggers inadvertent spreading of repressive marks to otherwise expressed neighboring genes, ultimately jeopardizing host fitness. In this Review, we provide a long-overdue overview of the recent genome-wide evidence for the presence and prevalence of TEs' epigenetic effects, highlighting both the similarities and differences across mammals, insects, and plants. We lay out the current understanding of the functional and fitness consequences of TEs' epigenetic effects, and propose possible influences of such effects on the evolution of both hosts and TEs themselves. These unique evolutionary consequences indicate that TEs' epigenetic effect is not only a crucial component of TE biology but could also be a significant contributor to genome function and evolution.

Adaptive evolution among cytoplasmic piRNA proteins leads to decreased genomic auto-immunity

In metazoan germlines, the piRNA pathway acts as a genomic immune system, employing small RNA-mediated silencing to defend host DNA from the harmful effects of transposable elements (TEs). Expression of genomic TEs is proposed to initiate self regulation by increasing the production of repressive piRNAs, thereby “adapting” piRNA-mediated control to the most active TE families. Surprisingly, however, piRNA pathway proteins, which execute piRNA biogenesis and enforce silencing of targeted sequences, evolve rapidly and adaptively in animals. If TE silencing is ensured through piRNA biogenesis, what necessitates changes in piRNA pathway proteins? Here we used interspecific complementation to test for functional differences between Drosophila melanogaster and D. simulans alleles of three adaptively evolving piRNA pathway proteins: Armitage, Aubergine and Spindle-E. In contrast to piRNA-mediated transcriptional regulators examined in previous studies, these three proteins have cytoplasmic functions in piRNA maturation and post-transcriptional silencing. Across all three proteins we observed interspecific divergence in the regulation of only a handful of TE families, which were more robustly silenced by the heterospecific piRNA pathway protein. This unexpected result suggests that unlike transcriptional regulators, positive selection has not acted on cytoplasmic piRNA effector proteins to enhance their function in TE repression. Rather, TEs may evolve to “escape” silencing by host proteins. We further discovered that D. simulans alleles of aub and armi exhibit enhanced off-target effects on host transcripts in a D. melanogaster background, as well as modest reductions in the efficiency of piRNA biogenesis, suggesting that promiscuous binding of D. simulans Aub and Armi proteins to host transcripts reduces their participation in piRNA production. Avoidance of genomic auto-immunity may therefore be a critical target of selection. Our observations suggest that piRNA effector proteins are subject to an evolutionary trade-off between defending the host genome from the harmful effect of TEs while also minimizing collateral damage to host genes.

Transposable elements are mobile fragments of selfish DNA that burden host genomes with deleterious mutations and incite genome instability. Host cells employ a specialized small-RNA mediated silencing pathway, the piRNA pathway, to act as a genomic immune system suppressing the mobilization of TEs. Changes in genomic TE content are met with rapid changes in the piRNA pool, thereby maintaining host control over transposition. However, piRNA pathway proteins—which enact piRNA biogenesis and silence target TEs—also evolve adaptively. To isolate forces that underlie this adaptive evolution, we examined functional divergence between two Drosophila species for three adaptively evolving piRNA pathway proteins. To our surprise, we found very few differences in TE regulation, suggesting that evolution has not generally acted to enhance control of TE parasites. Rather, we discovered interspecific differences in the regulation of host mRNAs for two proteins, which suggested that proteins evolve to avoid off-target silencing of host transcripts. We propose that the avoidance of such “genomic autoimmunity” is an important and underappreciated force driving the adaptive evolution of piRNA proteins.

Paternally Inherited P-Element Copy Number Affects the Magnitude of Hybrid Dysgenesis in Drosophila simulans and D. melanogaster

Transposable elements (TEs) are repetitive regions of DNA that are able to self-replicate and reinsert themselves throughout host genomes. Since the discovery of TEs, a prevalent question has been whether increasing TE copy number has an effect on the fitness of their hosts. P-elements (PEs) in Drosophila are a well-studied TE that has strong phenotypic effects. When a female without PEs (M) is crossed to a male with them (P), the resulting females are often sterile, a phenomenon called hybrid dysgenesis (HD). Here, we used short- and long-read sequencing to infer the number of PEs in the genomes of dozens of isofemale lines from two Drosophila species and measured whether the magnitude of HD was correlated with the number of PEs in the paternal genome. Consistent with previous reports, we find evidence for a positive correlation between the paternal PE copy number and the magnitude of HD in progeny from ♀M × ♂ P crosses for both species. Other crosses are not affected by the number of PE copies. We also find that the correlation between the strength of HD and PE copy number differs between species, which suggests that there are genetic differences that might make some genomes more resilient to the potentially deleterious effects of TEs. Our results suggest that PE copy number interacts with other factors in the genome and the environment to cause HD and that the importance of these interactions is species specific.

Transposable elements in individual genotypes of Drosophila simulans

Transposable elements are abundant, dynamic components of the genome that affect organismal phenotypes and fitness. In Drosophila melanogaster, they have increased in abundance as the species spread out of Africa, and different populations differ in their transposable element content. However, very little is currently known about how transposable elements differ between individual genotypes, and how that relates to the population dynamics of transposable elements overall. The sister species of D. melanogaster, D. simulans, has also recently become cosmopolitan, and panels of inbred genotypes exist from cosmopolitan and African flies. Therefore, we can determine whether the differences in colonizing populations are repeated in D. simulans, what the dynamics of transposable elements are in individual genotypes, and how that compares to wild flies. After estimating copy number in cosmopolitan and African D. simulans, I find that transposable element load is higher in flies from cosmopolitan populations. In addition, transposable element load varies considerably between populations, between genotypes, but not overall between wild and inbred lines. Certain genotypes either contain active transposable elements or are more permissive of transposition and accumulate copies of particular transposable elements. Overall, it is important to quantify genotype-specific transposable element dynamics as well as population averages to understand the dynamics of transposable element accumulation over time.

Pericentromeric heterochromatin is hierarchically organized and spatially contacts H3K9me2 islands in euchromatin

Membraneless pericentromeric heterochromatin (PCH) domains play vital roles in chromosome dynamics and genome stability. However, our current understanding of 3D genome organization does not include PCH domains because of technical challenges associated with repetitive sequences enriched in PCH genomic regions. We investigated the 3D architecture of Drosophila melanogaster PCH domains and their spatial associations with the euchromatic genome by developing a novel analysis method that incorporates genome-wide Hi-C reads originating from PCH DNA. Combined with cytogenetic analysis, we reveal a hierarchical organization of the PCH domains into distinct “territories.” Strikingly, H3K9me2-enriched regions embedded in the euchromatic genome show prevalent 3D interactions with the PCH domain. These spatial contacts require H3K9me2 enrichment, are likely mediated by liquid-liquid phase separation, and may influence organismal fitness. Our findings have important implications for how PCH architecture influences the function and evolution of both repetitive heterochromatin and the gene-rich euchromatin.

The three dimensional (3D) organization of genomes in cell nuclei can influence a wide variety of genome functions. However, most of our understanding of this critical architecture has been limited to the gene-rich euchromatin, and largely ignores the gene-poor and repeat-rich pericentromeric heterochromatin, or PCH. PCH comprises a large part of most eukaryotic genomes, forms 3D membraneless PCH domains in nuclei, and plays a vital role in chromosome dynamics and genome stability. In this study, we developed a new method that overcomes the technical challenges imposed by the highly repetitive PCH DNA, and generated a comprehensive picture of its 3D organization. Combined with image analyses, we reveal a hierarchical organization of the PCH domains. Surprisingly, we showed that distant euchromatic regions enriched for repressive epigenetic marks also dynamically interact with the main PCH domains. These 3D interactions are likely mediated by liquid-liquid phase separation (similar to how oil and vinegar separate in salad dressing) and the resulting liquid-like fusion events, and can influence the fitness of individuals. Our discoveries have strong implications for how seemingly “junk” DNA could impact functions in the gene-rich euchromatin.

Nanopore sequencing and Hi-C scaffolding provide insight into the evolutionary dynamics of transposable elements and piRNA production in wild strains of Drosophila melanogaster

Illumina sequencing has allowed for population-level surveys of transposable element (TE) polymorphism via split alignment approaches, which has provided important insight into the population dynamics of TEs. However, such approaches are not able to identify insertions of uncharacterized TEs, nor can they assemble the full sequence of inserted elements. Here, we use nanopore sequencing and Hi-C scaffolding to produce de novo genome assemblies for two wild strains of Drosophila melanogaster from the Drosophila Genetic Reference Panel (DGRP). Ovarian piRNA populations and Illumina split-read TE insertion profiles have been previously produced for both strains. We find that nanopore sequencing with Hi-C scaffolding produces highly contiguous, chromosome-length scaffolds, and we identify hundreds of TE insertions that were missed by Illumina-based methods, including a novel micropia-like element that has recently invaded the DGRP population. We also find hundreds of piRNA-producing loci that are specific to each strain. Some of these loci are created by strain-specific TE insertions, while others appear to be epigenetically controlled. Our results suggest that Illumina approaches reveal only a portion of the repetitive sequence landscape of eukaryotic genomes and that population-level resequencing using long reads is likely to provide novel insight into the evolutionary dynamics of repetitive elements.

Dynamics of Transposable Element Invasions with piRNA Clusters

In mammals and invertebrates, the proliferation of an invading transposable element (TE) is thought to be stopped by an insertion into a piRNA cluster. Here, we explore the dynamics of TE invasions under this trap model using computer simulations. We found that piRNA clusters confer a substantial benefit, effectively preventing extinction of host populations from a proliferation of deleterious TEs. TE invasions consist of three distinct phases: first, the TE amplifies within the population, next TE proliferation is stopped by segregating cluster insertions, and finally the TE is inactivated by fixation of a cluster insertion. Suppression by segregating cluster insertions is unstable and bursts of TE activity may yet occur. The transposition rate and the population size mostly influence the length of the phases but not the amount of TEs accumulating during an invasion. Solely, the size of piRNA clusters was identified as a major factor influencing TE abundance. We found that a single nonrecombining cluster is more efficient in stopping invasions than clusters distributed over several chromosomes. Recombination among cluster sites makes it necessary that each diploid carries, on the average, four cluster insertions to stop an invasion. Surprisingly, negative selection in a model with piRNA clusters can lead to a novel equilibrium state, where TE copy numbers remain stable despite only some individuals in a population carrying a cluster insertion. In Drosophila melanogaster, the trap model accounts for the abundance of TEs produced in the germline but fails to predict the abundance of TEs produced in the soma.

Genome-wide variation in nucleotides and retrotransposons in alpine populations of Arabis alpina (Brassicaceae)

Advances in high-throughput sequencing have promoted the collection of reference genomes and genome-wide diversity. However, the assessment of genomic variation among populations has hitherto mainly been surveyed through single-nucleotide polymorphisms (SNPs) and largely ignored the often major fraction of genomes represented by transposable elements (TEs). Despite accumulating evidence supporting the evolutionary significance of TEs, comprehensive surveys remain scarce. Here, we sequenced the full genomes of 304 individuals of Arabis alpina sampled from four nearby natural populations to genotype SNPs as well as polymorphic long terminal repeat retrotransposons (polymorphic TEs; i.e., presence/absence of TE insertions at specific loci). We identified 291,396 SNPs and 20,548 polymorphic TEs, comparing their contributions to genomic diversity and divergence across populations. Few SNPs were shared among populations and overall showed high population-specific variation, whereas most polymorphic TEs segregated among populations. The genomic context of these two classes of variants further highlighted candidate adaptive loci having a putative impact on functional genes. In particular, 4.96% of the SNPs were identified as nonsynonymous or affecting start/stop codons. In contrast, 43% of the polymorphic TEs were present next to Arabis genes enriched in functional categories related to the regulation of reproduction and responses to biotic as well as abiotic stresses. This unprecedented data set, mapping variation gained from SNPs and complementary polymorphic TEs within and among populations, will serve as a rich resource for addressing microevolutionary processes shaping genome variation.