Recurrent and recent selective sweeps in the piRNA pathway

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.

Molecular Insights into Female Hybrid Sterility in Interspecific Crosses between Drosophila melanogaster and Drosophila simulans

Species of the genus Drosophila have served as favorite models in speciation studies; however, genetic factors of interspecific reproductive incompatibility are under-investigated. Here, we performed an analysis of hybrid female sterility by crossing Drosophila melanogaster females and Drosophila simulans males. Using transcriptomic data analysis and molecular, cellular, and genetic approaches, we analyzed differential gene expression, transposable element (TE) activity, piRNA biogenesis, and functional defects of oogenesis in hybrids. Premature germline stem cell loss was the most prominent defect of oogenesis in hybrid ovaries. Because of the differential expression of genes encoding piRNA pathway components, rhino and deadlock, the functional RDCmel complex in hybrid ovaries was not assembled. However, the activity of the RDCsim complex was maintained in hybrids independent of the genomic origin of piRNA clusters. Despite the identification of a cohort of overexpressed TEs in hybrid ovaries, we found no evidence that their activity can be considered the main cause of hybrid sterility. We revealed a complicated pattern of Vasa protein expression in the hybrid germline, including partial AT-chX piRNA targeting of the vasasim allele and a significant zygotic delay in vasamel expression. We arrived at the conclusion that the hybrid sterility phenotype was caused by intricate multi-locus differences between the species.

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.

Off-target piRNA gene silencing in Drosophila melanogaster rescued by a transposable element insertion

Transposable elements (TE) are selfish genetic elements that can cause harmful mutations. In Drosophila, it has been estimated that half of all spontaneous visible marker phenotypes are mutations caused by TE insertions. Several factors likely limit the accumulation of exponentially amplifying TEs within genomes. First, synergistic interactions between TEs that amplify their harm with increasing copy number are proposed to limit TE copy number. However, the nature of this synergy is poorly understood. Second, because of the harm posed by TEs, eukaryotes have evolved systems of small RNA-based genome defense to limit transposition. However, as in all immune systems, there is a cost of autoimmunity and small RNA-based systems that silence TEs can inadvertently silence genes flanking TE insertions. In a screen for essential meiotic genes in Drosophila melanogaster, a truncated Doc retrotransposon within a neighboring gene was found to trigger the germline silencing of ald, the Drosophila Mps1 homolog, a gene essential for proper chromosome segregation in meiosis. A subsequent screen for suppressors of this silencing identified a new insertion of a Hobo DNA transposon in the same neighboring gene. Here we describe how the original Doc insertion triggers flanking piRNA biogenesis and local gene silencing. We show that this local gene silencing occurs in cis and is dependent on deadlock, a component of the Rhino-Deadlock-Cutoff (RDC) complex, to trigger dual-strand piRNA biogenesis at TE insertions. We further show how the additional Hobo insertion leads to de-silencing by reducing flanking piRNA biogenesis triggered by the original Doc insertion. These results support a model of TE-mediated gene silencing by piRNA biogenesis in cis that depends on local determinants of transcription. This may explain complex patterns of off-target gene silencing triggered by TEs within populations and in the laboratory. It also provides a mechanism of sign epistasis among TE insertions, illuminates the complex nature of their interactions and supports a model in which off-target gene silencing shapes the evolution of the RDC complex.

Genome-wide methylation in the panmictic European eel (Anguilla anguilla)

The role of methylation in adaptive, developmental and speciation processes has attracted considerable interest, but interpretation of results is complicated by diffuse boundaries between genetic and non-genetic variation. We studied whole genome genetic and methylation variation in the European eel, distributed from subarctic to subtropical environments, but with panmixia precluding genetically based local adaptation beyond single-generation responses. Overall methylation was 70.9%, with hypomethylation predominantly found in promoters and first exons. Redundancy analyses involving juvenile glass eels showed 0.06% and 0.03% of the variance at SNPs to be explained by localities and environmental variables, respectively, with GO terms of genes associated with outliers primarily involving neural system functioning. For CpGs 2.98% and 1.36% of variance was explained by localities and environmental variables. Differentially methylated regions particularly included genes involved in developmental processes, with hox clusters featuring prominently. Life stage (adult versus glass eels) was the most important source of inter-individual variation in methylation, likely reflecting both ageing and developmental processes. Demethylation of transposable elements relative to pure European eel was observed in European X American eel hybrids, possibly representing postzygotic barriers in this system characterized by prolonged speciation and ongoing gene flow. Whereas the genetic data are consistent with a role of single-generation selective responses, the methylation results underpin the importance of epigenetics in the life cycle of eels and suggests interactions between local environments, development and phenotypic variation mediated by methylation variation. Eels are remarkable by having retained eight hox clusters, and the results suggest important roles of methylation at hox genes for adaptive processes.

Dynamics and Impacts of Transposable Element Proliferation in the Drosophila nasuta Species Group Radiation

Transposable element (TE) mobilization is a constant threat to genome integrity. Eukaryotic organisms have evolved robust defensive mechanisms to suppress their activity, yet TEs can escape suppression and proliferate, creating strong selective pressure for host defense to adapt. This genomic conflict fuels a never-ending arms race that drives the rapid evolution of TEs and recurrent positive selection of genes involved in host defense; the latter has been shown to contribute to postzygotic hybrid incompatibility. However, how TE proliferation impacts genome and regulatory divergence remains poorly understood. Here, we report the highly complete and contiguous (N50 = 33.8-38.0 Mb) genome assemblies of seven closely related Drosophila species that belong to the nasuta species group-a poorly studied group of flies that radiated in the last 2 My. We constructed a high-quality de novo TE library and gathered germline RNA-seq data, which allowed us to comprehensively annotate and compare TE insertion patterns between the species, and infer the evolutionary forces controlling their spread. We find a strong negative association between TE insertion frequency and expression of genes nearby; this likely reflects survivor bias from reduced fitness impact of TEs inserting near lowly expressed, nonessential genes, with limited TE-induced epigenetic silencing. Phylogenetic analyses of insertions of 147 TE families reveal that 53% of them show recent amplification in at least one species. The most highly amplified TE is a nonautonomous DNA element (Drosophila INterspersed Element; DINE) which has gone through multiple bouts of expansions with thousands of full-length copies littered throughout each genome. Across all TEs, we find that TEs expansions are significantly associated with high expression in the expanded species consistent with suppression escape. Thus, whereas horizontal transfer followed by the invasion of a naïve genome has been highlighted to explain the long-term survival of TEs, our analysis suggests that evasion of host suppression of resident TEs is a major strategy to persist over evolutionary times. Altogether, our results shed light on the heterogenous and context-dependent nature in which TEs affect gene regulation and the dynamics of rampant TE proliferation amidst a recently radiated species group.

Evolutionary dynamics and conserved function of the Tudor domain-containing (TDRD) proteins in teleost fish

Tudor domain-containing (TDRD) proteins, the germline enriched protein family, play essential roles in the process of gametogenesis and genome stability through their interaction with the PIWI-interacting RNA (piRNA) pathway. Several studies have suggested the rapid evolution of the piRNA pathway in teleost lineages with striking reproductive diversity. However, there is still limited information about the function and evolution of Tdrd genes in teleost species. In this study, through genome wide screening, 13 Tdrd family genes were identified in economically important aquaculture fish, including spotted sea bass (Lateolabrax maculatus), Asian sea bass (Lates calcarifer), and tongue sole (Cynoglossus semilaevis). With copy number, structure, phylogeny, and synteny analysis, duplication of Tdrd6 and Tdrd7, as well as loss of Stk31 and Tdrd10, were characterized in teleost lineages. Codon based molecular evolution analysis indicated faster evolution of teleost Tdrd genes than that in mammals, potentially associated with the accelerated evolution of the piRNA pathway in teleost lineages. The evolutionary diversity of Tdrd genes was also detected between different teleost lineages. RNA-seq analysis showed that most teleost Tdrd genes were dominantly expressed in gonads, particularly highly expressed in testis, such as Tdrd6, Tdrd7a, Tdrd9, Ecat8, and Tdrd15. The varied expression and evolutionary pattern between the duplicated Tdrd6 and Tdrd7 in teleosts may indicate their functional diversification. All these results suggest a conserved function of teleost Tdrd family in gametogenesis and the piRNA pathway, which could lay a foundation for the evolution of Tdrd genes and be helpful for further deciphering of Tdrd functions in teleosts.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-021-00118-7.

A transposon expression burst accompanies the activation of Y-chromosome fertility genes during Drosophila spermatogenesis

Transposable elements (TEs) must replicate in germline cells to pass novel insertions to offspring. In Drosophila melanogaster ovaries, TEs can exploit specific developmental windows of opportunity to evade host silencing and increase their copy numbers. However, TE activity and host silencing in the distinct cell types of Drosophila testis are not well understood. Here, we reanalyze publicly available single-cell RNA-seq datasets to quantify TE expression in the distinct cell types of the Drosophila testis. We develop a method for identification of TE and host gene expression modules and find that a distinct population of early spermatocytes expresses a large number of TEs at much higher levels than other germline and somatic components of the testes. This burst of TE expression coincides with the activation of Y chromosome fertility factors and spermatocyte-specific transcriptional regulators, as well as downregulation of many components of the piRNA pathway. The TEs expressed by this cell population are specifically enriched on the Y chromosome and depleted on the X chromosome, relative to other active TEs. These data suggest that some TEs may achieve high insertional activity in males by exploiting a window of opportunity for mobilization created by the activation of spermatocyte-specific and Y chromosome-specific transcriptional programs.

The birth of piRNAs: how mammalian piRNAs are produced, originated, and evolved

PIWI-interacting RNAs (piRNAs), small noncoding RNAs 24–35 nucleotides long, are essential for animal fertility. They play critical roles in a range of functions, including transposable element suppression, gene expression regulation, imprinting, and viral defense. In mammals, piRNAs are the most abundant small RNAs in adult testes and the only small RNAs that direct epigenetic modification of chromatin in the nucleus. The production of piRNAs is a complex process from transcription to post-transcription, requiring unique machinery often distinct from the biogenesis of other RNAs. In mice, piRNA biogenesis occurs in specialized subcellular locations, involves dynamic developmental regulation, and displays sexual dimorphism. Furthermore, the genomic loci and sequences of piRNAs evolve much more rapidly than most of the genomic regions. Understanding piRNA biogenesis should reveal novel RNA regulations recognizing and processing piRNA precursors and the forces driving the gain and loss of piRNAs during animal evolution. Such findings may provide the basis for the development of engineered piRNAs capable of modulating epigenetic regulation, thereby offering possible single-dose RNA therapy without changing the genomic DNA. In this review, we focus on the biogenesis of piRNAs in mammalian adult testes that are derived from long non-coding RNAs. Although piRNA biogenesis is believed to be evolutionarily conserved from fruit flies to humans, recent studies argue for the existence of diverse, mammalian-specific RNA-processing pathways that convert precursor RNAs into piRNAs, perhaps associated with the unique features of mammalian piRNAs or germ cell development. We end with the discussion of major questions in the field, including substrate recognition and the birth of new piRNAs.

Large Drosophila germline piRNA clusters are evolutionarily labile and dispensable for transposon regulation

PIWI proteins and their guiding Piwi-interacting small RNAs (piRNAs) are crucial for fertility and transposon defense in the animal germline. In most species, the majority of piRNAs are produced from distinct large genomic loci, called piRNA clusters. It is assumed that germline-expressed piRNA clusters, particularly in Drosophila, act as principal regulators to control transposons dispersed across the genome. Here, using synteny analysis, we show that large clusters are evolutionarily labile, arise at loci characterized by recurrent chromosomal rearrangements, and are mostly species-specific across the Drosophila genus. By engineering chromosomal deletions in D. melanogaster, we demonstrate that the three largest germline clusters, which account for the accumulation of >40% of all transposon-targeting piRNAs in ovaries, are neither required for fertility nor for transposon regulation in trans. We provide further evidence that dispersed elements, rather than the regulatory action of large Drosophila germline clusters in trans, may be central for transposon defense.

P-elements strengthen reproductive isolation within the Drosophila simulans species complex

Determining mechanisms that underlie reproductive isolation (RI) is key to understanding how species boundaries are maintained in nature. Transposable elements (TEs) are ubiquitous across eukaryotic genomes. However, the role of TEs in modulating the strength of RI between species is poorly understood. Several species of Drosophila have been found to harbor P-elements (PEs), yet only D. simulans is known to be currently polymorphic for their presence in wild populations. PEs can cause RI between PE-containing (P) and PE-lacking (M) lineages of the same species. However, it is unclear whether they also contribute to the magnitude of RI between species. Here, we use the simulans species complex to assess whether differences in PE status between D. simulans and its sister species, which do not harbor PEs, contribute to multiple barriers to gene flow between species. We show that crosses involving a P D. simulans father and an M mother from a sister species exhibit lower F1 female fecundity than crosses involving an M D. simulans father and an M sister-species mother. We also find that another TE, I-element, might play a minor role in determining the frequency of dysgenesis between species. Our results suggest that the presence of PEs in a species can strengthen isolation from its sister species, providing evidence that TEs can play a role in RI.

Protein-Protein Interactions Shape Genomic Autoimmunity in the Adaptively Evolving Rhino-Deadlock-Cutoff Complex

The Piwi-interacting RNA (piRNA) pathway is a genomic defense system that controls the movement of transposable elements (TEs) through transcriptional and post-transcriptional silencing. Although TE defense is critical to ensuring germline genome integrity, it is equally critical that the piRNA pathway avoids autoimmunity in the form of silencing host genes. Ongoing cycles of selection for expanded control of invading TEs, followed by selection for increased specificity to reduce impacts on host genes, are proposed to explain the frequent signatures of adaptive evolution among piRNA pathway proteins. However, empirical tests of this model remain limited, particularly with regards to selection against genomic autoimmunity.

I examined three adaptively evolving piRNA proteins, Rhino, Deadlock, and Cutoff, for evidence of interspecific divergence in autoimmunity between Drosophila melanogaster and Drosophila simulans. I tested a key prediction of the autoimmunity hypothesis that foreign heterospecific piRNA proteins will exhibit enhanced autoimmunity, due to the absence of historical selection against off-target effects. Consistent with this prediction, full-length D. simulans Cutoff, as well as the D. simulans hinge and chromo domains of Rhino, exhibit expanded regulation of D. melanogaster genes. I further demonstrate that this autoimmunity is dependent on known incompatibilities between D. simulans proteins or domains and their interacting partners in D. melanogaster. My observations reveal that the same protein–protein interaction domains that are interfaces of adaptive evolution in Rhino and Cutoff also determine their potential for autoimmunity.

Dynamic Evolution of De Novo DNA Methyltransferases in Rodent and Primate Genomes

Transcriptional silencing of retrotransposons via DNA methylation is paramount for mammalian fertility and reproductive fitness. During germ cell development, most mammalian species utilize the de novo DNA methyltransferases DNMT3A and DNMT3B to establish DNA methylation patterns. However, many rodent species deploy a third enzyme, DNMT3C, to selectively methylate the promoters of young retrotransposon insertions in their germline. The evolutionary forces that shaped DNMT3C’s unique function are unknown. Using a phylogenomic approach, we confirm here that Dnmt3C arose through a single duplication of Dnmt3B that occurred ∼60 Ma in the last common ancestor of muroid rodents. Importantly, we reveal that DNMT3C is composed of two independently evolving segments: the latter two-thirds have undergone recurrent gene conversion with Dnmt3B, whereas the N-terminus has instead evolved under strong diversifying selection. We hypothesize that positive selection of Dnmt3C is the result of an ongoing evolutionary arms race with young retrotransposon lineages in muroid genomes. Interestingly, although primates lack DNMT3C, we find that the N-terminus of DNMT3A has also evolved under diversifying selection. Thus, the N-termini of two independent de novo methylation enzymes have evolved under diversifying selection in rodents and primates. We hypothesize that repression of young retrotransposons might be driving the recurrent innovation of a functional domain in the N-termini on germline DNMT3s in mammals.

Adaptive Evolution Targets a piRNA Precursor Transcription Network

In Drosophila, transposon-silencing piRNAs are derived from heterochromatic clusters and a subset of euchromatic transposon insertions, which are bound by the Rhino-Deadlock-Cutoff complex. The HP1 homolog Rhino binds to Deadlock, which recruits TRF2 to promote non-canonical transcription from both genomic strands. Cuff function is less well understood, but this Rai1 homolog shows hallmarks of adaptive evolution, which can remodel functional interactions within host defense systems. Supporting this hypothesis, Drosophila simulans Cutoff is a dominant-negative allele when expressed in Drosophila melanogaster, in which it traps Deadlock, TRF2, and the conserved transcriptional co-repressor CtBP in stable complexes. Cutoff functions with Rhino and Deadlock to drive non-canonical transcription. In contrast, CtBP suppresses canonical transcription of transposons and promoters flanking the major germline clusters, and canonical transcription interferes with downstream non-canonical transcription and piRNA production. Adaptive evolution thus targets interactions among Cutoff, TRF2, and CtBP that balance canonical and non-canonical piRNA precursor transcription.

Telomeric TART elements target the piRNA machinery in Drosophila

Coevolution between transposable elements (TEs) and their hosts can be antagonistic, where TEs evolve to avoid silencing and the host responds by reestablishing TE suppression, or mutualistic, where TEs are co-opted to benefit their host. The TART-A TE functions as an important component of Drosophila telomeres but has also reportedly inserted into the Drosophila melanogaster nuclear export factor gene nxf2. We find that, rather than inserting into nxf2, TART-A has actually captured a portion of nxf2 sequence. We show that TART-A produces abundant Piwi-interacting small RNAs (piRNAs), some of which are antisense to the nxf2 transcript, and that the TART-like region of nxf2 is evolving rapidly. Furthermore, in D. melanogaster, TART-A is present at higher copy numbers, and nxf2 shows reduced expression, compared to the closely related species Drosophila simulans. We propose that capturing nxf2 sequence allowed TART-A to target the nxf2 gene for piRNA-mediated repression and that these 2 elements are engaged in antagonistic coevolution despite the fact that TART-A is serving a critical role for its host genome.

Co-evolution between transposable elements (TEs) and their hosts can be antagonistic, where TEs evolve to avoid silencing and the host responds by re-establishing TE suppression, or mutualistic, where TEs are co-opted to benefit their host. This study shows that a specialized Drosophila retrotransposon that functions as a telomere has captured a portion of a host piRNA gene which may allow it to evade silencing.

Establishment and evolution of heterochromatin

The eukaryotic genome is packaged into transcriptionally active euchromatin and silent heterochromatin, with most studies focused on the former encompassing the majority of protein-coding genes. The recent development of various sequencing techniques has refined this classic dichromatic partition and has better illuminated the composition, establishment, and evolution of this genomic and epigenomic "dark matter" in the context of topologically associated domains and phase-separated droplets. Heterochromatin includes genomic regions that can be densely stained by chemical dyes, which have been shown to be enriched for repetitive elements and epigenetic marks, including H3K9me2/3 and H3K27me3. Heterochromatin is usually replicated late, concentrated at the nuclear periphery or around nucleoli, and usually lacks highly expressed genes; and now it is considered to be as neither genetically inert nor developmentally static. Heterochromatin guards genome integrity against transposon activities and exerts important regulatory functions by targeting beyond its contained genes. Both its nucleotide sequences and regulatory proteins exhibit rapid coevolution between species. In addition, there are dynamic transitions between euchromatin and heterochromatin during developmental and evolutionary processes. We summarize here the ever-changing characteristics of heterochromatin and propose models and principles for the evolutionary transitions of heterochromatin that have been mainly learned from studies of Drosophila and yeast. Finally, we highlight the role of sex chromosomes in studying heterochromatin evolution.

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.

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.

Rapid evolution and conserved function of the piRNA pathway

Transposons are major genome constituents that can mobilize and trigger mutations, DNA breaks and chromosome rearrangements. Transposon silencing is particularly important in the germline, which is dedicated to transmission of the inherited genome. Piwi-interacting RNAs (piRNAs) guide a host defence system that transcriptionally and post-transcriptionally silences transposons during germline development. While germline control of transposons by the piRNA pathway is conserved, many piRNA pathway genes are evolving rapidly under positive selection, and the piRNA biogenesis machinery shows remarkable phylogenetic diversity. Conservation of core function combined with rapid gene evolution is characteristic of a host–pathogen arms race, suggesting that transposons and the piRNA pathway are engaged in an evolutionary tug of war that is driving divergence of the biogenesis machinery. Recent studies suggest that this process may produce biochemical incompatibilities that contribute to reproductive isolation and species divergence.

Drosophila Interspecific Hybridization Causes A Deregulation of the piRNA Pathway Genes

Almost all eukaryotes have transposable elements (TEs) against which they have developed defense mechanisms. In the Drosophila germline, the main transposable element (TE) regulation pathway is mediated by specific Piwi-interacting small RNAs (piRNAs). Nonetheless, for unknown reasons, TEs sometimes escape cellular control during interspecific hybridization processes. Because the piRNA pathway genes are involved in piRNA biogenesis and TE control, we sequenced and characterized nine key genes from this pathway in Drosophila buzzatii and Drosophila koepferae species and studied their expression pattern in ovaries of both species and their F1 hybrids. We found that gene structure is, in general, maintained between both species and that two genes—armitage and aubergine—are under positive selection. Three genes—krimper, methyltransferase 2, and zucchini—displayed higher expression values in hybrids than both parental species, while others had RNA levels similar to the parental species with the highest expression. This suggests that the overexpression of some piRNA pathway genes can be a primary response to hybrid stress. Therefore, these results reinforce the hypothesis that TE deregulation may be due to the protein incompatibility caused by the rapid evolution of these genes, leading to a TE silencing failure, rather than to an underexpression of piRNA pathway genes.

Evolutionary dynamics of microRNA target sites across vertebrate evolution

MicroRNAs (miRNAs) control the abundance of the majority of the vertebrate transcriptome. The recognition sequences, or target sites, for bilaterian miRNAs are found predominantly in the 3' untranslated regions (3'UTRs) of mRNAs, and are amongst the most highly conserved motifs within 3'UTRs. However, little is known regarding the evolutionary pressures that lead to loss and gain of such target sites. Here, we quantify the selective pressures that act upon miRNA target sites. Notably, selective pressure extends beyond deeply conserved binding sites to those that have undergone recent substitutions. Our approach reveals that even amongst ancient animal miRNAs, which exert the strongest selective pressures on 3'UTR sequences, there are striking differences in patterns of target site evolution between miRNAs. Considering only ancient animal miRNAs, we find three distinct miRNA groups, each exhibiting characteristic rates of target site gain and loss during mammalian evolution. The first group both loses and gains sites rarely. The second group shows selection only against site loss, with site gains occurring at a neutral rate, whereas the third loses and gains sites at neutral or above expected rates. Furthermore, mutations that alter the strength of existing target sites are disfavored. Applying our approach to individual transcripts reveals variation in the distribution of selective pressure across the transcriptome and between miRNAs, ranging from strong selection acting on a small subset of targets of some miRNAs, to weak selection on many targets for other miRNAs. miR-20 and miR-30, and many other miRNAs, exhibit broad, deeply conserved targeting, while several other comparably ancient miRNAs show a lack of selective constraint, and a small number, including mir-146, exhibit evidence of rapidly evolving target sites. Our approach adds valuable perspective on the evolution of miRNAs and their targets, and can also be applied to characterize other 3'UTR regulatory motifs.

Evolution and biological significance of flaviviral elements in the genome of the arboviral vector Aedes albopictus

Since its genome details are publically available, the mosquito Aedes albopictus has become the central stage of attention for deciphering multiple biological and evolutionary aspects at the root of its success as an invasive species. Its genome of 1,967 Mb harbours an unusual high number of non-retroviral integrated RNA virus sequences (NIRVS). NIRVS are enriched in piRNA clusters and produce piRNAs, suggesting an antiviral effect. Here, we investigated the evolutionary history of NIRVS in geographically distant Ae. albopictus populations by comparing genetic variation as derived by neutral microsatellite loci and seven selected NIRVS. We found that the evolution of NIRVS was far to be neutral with variations both in their distribution and sequence polymorphism among Ae. albopictus populations. The Flaviviral elements AlbFlavi2 and AlbFlavi36 were more deeply investigated in their association with dissemination rates of dengue virus (DENV) and chikungunya virus (CHIKV) in Ae. albopictus at both population and individual levels. Our results show a complex association between NIRVS and DENV/CHIKV opening a new avenue for investigating the functional role of NIRVS as antiviral elements shaping vector competence of mosquitoes to arboviruses.

piRNA silencing contributes to interspecies hybrid sterility and reproductive isolation in Drosophila melanogaster

The piRNA pathway is an adaptive mechanism that maintains genome stability by repression of selfish genomic elements. In the male germline of Drosophila melanogaster repression of Stellate genes by piRNAs generated from Supressor of Stellate (Su(Ste)) locus is required for male fertility, but both Su(Ste) piRNAs and their targets are absent in other Drosophila species. We found that D. melanogaster genome contains multiple X-linked non-coding genomic repeats that have sequence similarity to the protein-coding host gene vasa. In the male germline, these vasa-related AT-chX repeats produce abundant piRNAs that are antisense to vasa; however, vasa mRNA escapes silencing due to imperfect complementarity to AT-chX piRNAs. Unexpectedly, we discovered AT-chX piRNAs target vasa of Drosophila mauritiana in the testes of interspecies hybrids. In the majority of hybrid flies, the testes were strongly reduced in size and germline content. A minority of hybrids maintained wild-type array of premeiotic germ cells in the testes, but in them harmful Stellate genes were derepressed due to the absence of Su(Ste) piRNAs, and meiotic failures were observed. Thus, the piRNA pathway contributes to reproductive isolation between D. melanogaster and closely related species, causing hybrid male sterility via misregulation of two different host protein factors.

DNA methylation reprogramming, TE derepression, and postzygotic isolation of nascent animal species

The genomic shock hypothesis stipulates that the stress associated with divergent genome admixture can cause transposable element (TE) derepression, which could act as a postzygotic isolation mechanism. TEs affect gene structure, expression patterns, and chromosome organization and may have deleterious consequences when released. For these reasons, they are silenced by heterochromatin formation, which includes DNA methylation. Here, we show that a significant proportion of TEs are differentially methylated between the "dwarf" (limnetic) and the "normal" (benthic) whitefish, two nascent species that diverged some 15,000 generations ago within the Coregonus clupeaformis species complex. Moreover, TEs are overrepresented among loci that were demethylated in hybrids, indicative of their transcriptional derepression. These results are consistent with earlier studies in this system that revealed TE transcriptional derepression causes abnormal embryonic development and death of hybrids. Hence, this supports a role of DNA methylation reprogramming and TE derepression in postzygotic isolation of nascent animal species.

Host-transposon interactions: conflict, cooperation, and cooption

Transposable elements (TEs) are mobile DNA sequences that colonize genomes and threaten genome integrity. As a result, several mechanisms appear to have emerged during eukaryotic evolution to suppress TE activity. However, TEs are ubiquitous and account for a prominent fraction of most eukaryotic genomes. We argue that the evolutionary success of TEs cannot be explained solely by evasion from host control mechanisms. Rather, some TEs have evolved commensal and even mutualistic strategies that mitigate the cost of their propagation. These coevolutionary processes promote the emergence of complex cellular activities, which in turn pave the way for cooption of TE sequences for organismal function.

PIWI-interacting RNAs: small RNAs with big functions

In animals, PIWI-interacting RNAs (piRNAs) of 21-35 nucleotides in length silence transposable elements, regulate gene expression and fight viral infection. piRNAs guide PIWI proteins to cleave target RNA, promote heterochromatin assembly and methylate DNA. The architecture of the piRNA pathway allows it both to provide adaptive, sequence-based immunity to rapidly evolving viruses and transposons and to regulate conserved host genes. piRNAs silence transposons in the germ line of most animals, whereas somatic piRNA functions have been lost, gained and lost again across evolution. Moreover, most piRNA pathway proteins are deeply conserved, but different animals employ remarkably divergent strategies to produce piRNA precursor transcripts. Here, we discuss how a common piRNA pathway allows animals to recognize diverse targets, ranging from selfish genetic elements to genes essential for gametogenesis.

PIWI pathway against viruses in insects

Piwi-interacting RNAs (piRNAs) are an animal-specific class of small non-coding RNAs that are generated via a biogenesis pathway distinct from small interfering RNAs (siRNAs) and microRNAs (miRNAs). There are variations in piRNA biogenesis that depend on several factors, such as the cell type (germline or soma), the organism, and the purpose for which they are being produced, such as transposon-targeting, viral-targeting, or gene-derived piRNAs. Interestingly, the genes involved in the PIWI/piRNA pathway are more rapidly evolving compared with other RNA interference (RNAi) genes. In this review, the role of the piRNA pathway in the antiviral response is reviewed based on recent findings in insect models such as Drosophila, mosquitoes, midges and the silkworm, Bombyx mori. We extensively discuss the special features that characterize host-virus piRNA responses with respect to the proteins and the genes involved, the viral piRNAs' sequence characteristics, the target strand orientation biases as well as the viral piRNA target hotspots across the viral genomes. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > RNAi: Mechanisms of Action Regulatory RNAs/RNAi/Riboswitches > Biogenesis of Effector Small RNAs.

On the Population Dynamics of Junk: A Review on the Population Genomics of Transposable Elements

Transposable elements (TEs) play an important role in shaping genomic organization and structure, and may cause dramatic changes in phenotypes. Despite the genetic load they may impose on their host and their importance in microevolutionary processes such as adaptation and speciation, the number of population genetics studies focused on TEs has been rather limited so far compared to single nucleotide polymorphisms (SNPs). Here, we review the current knowledge about the dynamics of transposable elements at recent evolutionary time scales, and discuss the mechanisms that condition their abundance and frequency. We first discuss non-adaptive mechanisms such as purifying selection and the variable rates of transposition and elimination, and then focus on positive and balancing selection, to finally conclude on the potential role of TEs in causing genomic incompatibilities and eventually speciation. We also suggest possible ways to better model TEs dynamics in a population genomics context by incorporating recent advances in TEs into the rich information provided by SNPs about the demography, selection, and intrinsic properties of genomes.

Function and Evolution of Nematode RNAi Pathways

Selfish genetic elements, like transposable elements or viruses, are a threat to genomic stability. A variety of processes, including small RNA-based RNA interference (RNAi)-like pathways, has evolved to counteract these elements. Amongst these, endogenous small interfering RNA and Piwi-interacting RNA (piRNA) pathways were implicated in silencing selfish genetic elements in a variety of organisms. Nematodes have several incredibly specialized, rapidly evolving endogenous RNAi-like pathways serving such purposes. Here, we review recent research regarding the RNAi-like pathways of Caenorhabditis elegans as well as those of other nematodes, to provide an evolutionary perspective. We argue that multiple nematode RNAi-like pathways share piRNA-like properties and together form a broad nematode toolkit that allows for silencing of foreign genetic elements.

Hybrid Genome Evolution by Transposition: An Update

Contrary to the view that hybrids are lineages devoid of evolutionary value, a number of case studies that have been lately reported show how hybrids are at the origin of many species. Some well-documented cases demonstrate that bursts of transposition often follow hybridization, generating new genetic variability. Studies in hybrid transposition strongly suggest that epigenetic changes and divergence in piRNA pathways drive deregulation in TE landscapes. Here, I have focused on mechanisms acting in Drosophila hybrids between two cactophilic species. The results reported here show that while hybrid instability by transposition is a genome-wide event, deregulation by TE overexpression in hybrid ovaries is not a general rule. When piRNA pools of ovaries are studied, results show that TEs with parental differences higher than 2-fold in their piRNA amounts are not more commonly deregulated in hybrids than TEs with similar levels, partially discrediting the generality of the maternal cytotype hypothesis. Some promising results on the piRNA pathway global failure hypothesis, which states that accumulated divergence of piRNA effector proteins is responsible for hybrid TE deregulation, have also been obtained. Altogether, these results suggest that TE deregulation might be driven by several interacting mechanisms. A natural scenario is proposed in which genome instability by transposition leads to hybrid genome reorganization. Small hybrid populations, subjected to natural selection helped by genetic drift, evolve new adaptations adapted to novel environments. The final step is either introgression or even a new hybrid species.

Transposable Element Misregulation Is Linked to the Divergence between Parental piRNA Pathways in Drosophila Hybrids

Interspecific hybridization is a genomic stress condition that leads to the activation of transposable elements (TEs) in both animals and plants. In hybrids between Drosophila buzzatii and Drosophila koepferae, mobilization of at least 28 TEs has been described. However, the molecular mechanisms underlying this TE release remain poorly understood. To give insight on the causes of this TE activation, we performed a TE transcriptomic analysis in ovaries (notorious for playing a major role in TE silencing) of parental species and their F1 and backcrossed (BC) hybrids. We find that 15.2% and 10.6% of the expressed TEs are deregulated in F1 and BC1 ovaries, respectively, with a bias toward overexpression in both cases. Although differences between parental piRNA (Piwi-interacting RNA) populations explain only partially these results, we demonstrate that piRNA pathway proteins have divergent sequences and are differentially expressed between parental species. Thus, a functional divergence of the piRNA pathway between parental species, together with some differences between their piRNA pools, might be at the origin of hybrid instabilities and ultimately cause TE misregulation in ovaries. These analyses were complemented with the study of F1 testes, where TEs tend to be less expressed than in D. buzzatii. This can be explained by an increase in piRNA production, which probably acts as a defence mechanism against TE instability in the male germline. Hence, we describe a differential impact of interspecific hybridization in testes and ovaries, which reveals that TE expression and regulation are sex-biased.

Adaptive Evolution Leads to Cross-Species Incompatibility in the piRNA Transposon Silencing Machinery

Reproductive isolation defines species divergence and is linked to adaptive evolution of hybrid incompatibility genes. Hybrids between Drosophila melanogaster and Drosophila simulans are sterile, and phenocopy mutations in the PIWI interacting RNA (piRNA) pathway, which silences transposons and shows pervasive adaptive evolution, and Drosophila rhino and deadlock encode rapidly evolving components of a complex that binds to piRNA clusters. We show that Rhino and Deadlock interact and co-localize in simulans and melanogaster, but simulans Rhino does not bind melanogaster Deadlock, due to substitutions in the rapidly evolving Shadow domain. Significantly, a chimera expressing the simulans Shadow domain in a melanogaster Rhino backbone fails to support piRNA production, disrupts binding to piRNA clusters, and leads to ectopic localization to bulk heterochromatin. Fusing melanogaster Deadlock to simulans Rhino, by contrast, restores localization to clusters. Deadlock binding thus directs Rhino to piRNA clusters, and Rhino-Deadlock co-evolution has produced cross-species incompatibilities, which may contribute to reproductive isolation.

Paternal Induction of Hybrid Dysgenesis in Drosophila melanogaster Is Weakly Correlated with Both P-Element and hobo Element Dosage

Transposable elements (TEs) are virtually ubiquitous components of genomes, yet they often impose significant fitness consequences on their hosts. In addition to producing specific deleterious mutations by insertional inactivation, TEs also impose general fitness costs by inducing DNA damage and participating in ectopic recombination. These latter fitness costs are often assumed to be dosage-dependent, with stronger effects occurring in the presence of higher TE copy numbers. We test this assumption in Drosophila melanogaster by considering the relationship between the copy number of two active DNA transposons, P-element and hobo element, and the incidence of hybrid dysgenesis, a sterility syndrome associated with transposon activity in the germline. By harnessing a subset of the Drosophila Genetic Reference Panel (DGRP), a group of fully-sequenced D. melanogaster strains, we describe quantitative and structural variation in P-elements and hobo elements among wild-derived genomes and associate these factors with hybrid dysgenesis. We find that the incidence of hybrid dysgenesis is associated with both P-element and hobo element copy number in a dosage-dependent manner. However, the relationship is weak for both TEs, suggesting that dosage alone explains only a small part of TE-associated fitness costs.

Natural variation of piRNA expression affects immunity to transposable elements

In the Drosophila germline, transposable elements (TEs) are silenced by PIWI-interacting RNA (piRNA) that originate from distinct genomic regions termed piRNA clusters and are processed by PIWI-subfamily Argonaute proteins. Here, we explore the variation in the ability to restrain an alien TE in different Drosophila strains. The I-element is a retrotransposon involved in the phenomenon of I-R hybrid dysgenesis in Drosophila melanogaster. Genomes of R strains do not contain active I-elements, but harbour remnants of ancestral I-related elements. The permissivity to I-element activity of R females, called reactivity, varies considerably in natural R populations, indicating the existence of a strong natural polymorphism in defense systems targeting transposons. To reveal the nature of such polymorphisms, we compared ovarian small RNAs between R strains with low and high reactivity and show that reactivity negatively correlates with the ancestral I-element-specific piRNA content. Analysis of piRNA clusters containing remnants of I-elements shows increased expression of the piRNA precursors and enrichment by the Heterochromatin Protein 1 homolog, Rhino, in weak R strains, which is in accordance with stronger piRNA expression by these regions. To explore the nature of the differences in piRNA production, we focused on two R strains, weak and strong, and showed that the efficiency of maternal inheritance of piRNAs as well as the I-element copy number are very similar in both strains. At the same time, germline and somatic uni-strand piRNA clusters generate more piRNAs in strains with low reactivity, suggesting the relationship between the efficiency of primary piRNA production and variable response to TE invasions. The strength of adaptive genome defense is likely driven by naturally occurring polymorphisms in the rapidly evolving piRNA pathway proteins. We hypothesize that hyper-efficient piRNA production is contributing to elimination of a telomeric retrotransposon HeT-A, which we have observed in one particular transposon-resistant R strain.

Transposon activity in the germline is suppressed by the PIWI-interacting RNA (piRNA) pathway. The resistance of natural Drosophila strains to transposon invasion varies considerably, but the nature of this variability is unknown. We discovered that natural variation in the efficiency of primary piRNA production in the germline causes dramatic differences in the susceptibility to expansion of a newly invaded transposon. A high level of piRNA production in the germline is achieved by increased expression of piRNA precursors. In one of the most transposon-resistant strains, increased content of primary piRNA is observed in both the germline and ovarian somatic cells. We suggest that polymorphisms in piRNA pathway factors are responsible for increased piRNA production. piRNA pathway proteins have been shown to be evolving rapidly under selective pressure. Our data are the first to describe a phenotype that might be caused by this kind of polymorphism. We also demonstrate a likely explanation as to why an overly active piRNA pathway can cause more harm than good in Drosophila: Highly efficient piRNA processing leads to elimination of domesticated telomeric retrotransposons essential for telomere elongation, an effect which has been observed in a natural strain that is extremely resistant to transposon invasion.

Recurrent Innovation at Genes Required for Telomere Integrity in Drosophila

Telomeres are nucleoprotein complexes at the ends of linear chromosomes. These specialized structures ensure genome integrity and faithful chromosome inheritance. Recurrent addition of repetitive, telomere-specific DNA elements to chromosome ends combats end-attrition, while specialized telomere-associated proteins protect naked, double-stranded chromosome ends from promiscuous repair into end-to-end fusions. Although telomere length homeostasis and end-protection are ubiquitous across eukaryotes, there is sporadic but building evidence that the molecular machinery supporting these essential processes evolves rapidly. Nevertheless, no global analysis of the evolutionary forces that shape these fast-evolving proteins has been performed on any eukaryote. The abundant population and comparative genomic resources of Drosophila melanogaster and its close relatives offer us a unique opportunity to fill this gap. Here we leverage population genetics, molecular evolution, and phylogenomics to define the scope and evolutionary mechanisms driving fast evolution of genes required for telomere integrity. We uncover evidence of pervasive positive selection across multiple evolutionary timescales. We also document prolific expansion, turnover, and expression evolution in gene families founded by telomeric proteins. Motivated by the mutant phenotypes and molecular roles of these fast-evolving genes, we put forward four alternative, but not mutually exclusive, models of intra-genomic conflict that may play out at very termini of eukaryotic chromosomes. Our findings set the stage for investigating both the genetic causes and functional consequences of telomere protein evolution in Drosophila and beyond.

What Drives Positive Selection in the Drosophila piRNA Machinery? The Genomic Autoimmunity Hypothesis

In animals, PIWI-interacting RNAs (piRNAs) play a crucial role in genome defense. Moreover, because piRNAs can be maternally transmitted, they contribute to the epigenetic profile of inheritance. Multiple studies, especially in Drosophila, have demonstrated that the machinery of piRNA biogenesis is often the target of positive selection. Because transposable elements (TEs) are a form of genetic parasite, positive selection in the piRNA machinery is often explained by analogy to the signatures of positive selection commonly observed in genes that play a role in host-parasite dynamics. However, the precise mechanisms that drive positive selection in the piRNA machinery are not known. In this review, we outline several mechanistic models that might explain pervasive positive selection in the piRNA machinery of Drosophila species. We propose that recurrent positive selection in the piRNA machinery can be partly explained by an ongoing tension between selection for sensitivity required by genome defense and selection for specificity to avoid the off-target effects of maladaptive genic silencing by piRNA.

Recurrent Gene Duplication Diversifies Genome Defense Repertoire in Drosophila

Transposable elements (TEs) comprise large fractions of many eukaryotic genomes and imperil host genome integrity. The host genome combats these challenges by encoding proteins that silence TE activity. Both the introduction of new TEs via horizontal transfer and TE sequence evolution requires constant innovation of host-encoded TE silencing machinery to keep pace with TEs. One form of host innovation is the adaptation of existing, single-copy host genes. Indeed, host suppressors of TE replication often harbor signatures of positive selection. Such signatures are especially evident in genes encoding the piwi-interacting-RNA pathway of gene silencing, for example, the female germline-restricted TE silencer, HP1D/Rhino Host genomes can also innovate via gene duplication and divergence. However, the importance of gene family expansions, contractions, and gene turnover to host genome defense has been largely unexplored. Here, we functionally characterize Oxpecker, a young, tandem duplicate gene of HP1D/rhino We demonstrate that Oxpecker supports female fertility in Drosophila melanogaster and silences several TE families that are incompletely silenced by HP1D/Rhino in the female germline. We further show that, like Oxpecker, at least ten additional, structurally diverse, HP1D/rhino-derived daughter and "granddaughter" genes emerged during a short 15-million year period of Drosophila evolution. These young paralogs are transcribed primarily in germline tissues, where the genetic conflict between host genomes and TEs plays out. Our findings suggest that gene family expansion is an underappreciated yet potent evolutionary mechanism of genome defense diversification.

Expression of the Retrotransposon Helena Reveals a Complex Pattern of TE Deregulation in Drosophila Hybrids

Transposable elements (TEs), repeated mobile sequences, are ubiquitous in the eukaryotic kingdom. Their mobilizing capacity confers on them a high mutagenic potential, which must be strongly regulated to guarantee genome stability. In the Drosophila germline, a small RNA-mediated silencing system, the piRNA (Piwi-interacting RNA) pathway, is the main responsible TE regulating mechanism, but some stressful conditions can destabilize it. For instance, during interspecific hybridization, genomic stress caused by the shock of two different genomes can lead, in both animals and plants, to higher transposition rates. A recent study in D. buzatii-D. koepferae hybrids detected mobilization of 28 TEs, yet little is known about the molecular mechanisms explaining this transposition release. We have characterized one of the mobilized TEs, the retrotransposon Helena, and used quantitative expression to assess whether its high transposition rates in hybrids are preceded by increased expression. We have also localized Helena expression in the gonads to see if cellular expression patterns have changed in the hybrids. To give more insight into changes in TE regulation in hybrids, we analysed Helena-specific piRNA populations of hybrids and parental species. Helena expression is not globally altered in somatic tissues, but male and female gonads have different patterns of deregulation. In testes, Helena is repressed in F1, increasing then its expression up to parental values. This is linked with a mislocation of Helena transcripts along with an increase of their specific piRNA levels. Ovaries have additive levels of Helena expression, but the ping-pong cycle efficiency seems to be reduced in F1 hybrids. This could be at the origin of new Helena insertions in hybrids, which would be transmitted to F1 hybrid female progeny.

The evolution of gene expression and binding specificity of the largest transcription factor family in primates

The KRAB-containing zinc finger (KRAB-ZF) proteins represent the largest family of transcription factors (TFs) in humans, yet for the great majority, their function and specific genomic target remain unknown. However, it has been shown that a large fraction of these genes arose from segmental duplications, and that they have expanded in gene and zinc finger number throughout vertebrate evolution. To determine whether this expansion is linked to selective pressures acting on different domains, we have manually curated all KRAB-ZF genes present in the human genome together with their orthologous genes in three closely related species and assessed the evolutionary forces acting at the sequence level as well as on their expression profiles. We provide evidence that KRAB-ZFs can be separated into two categories according to the polymorphism present in their DNA-contacting residues. Those carrying a nonsynonymous single nucleotide polymorphism (SNP) in their DNA-contacting amino acids exhibit significantly reduced expression in all tissues, have emerged in a recent lineage, and seem to be less strongly constrained evolutionarily than those without such a polymorphism. This work provides evidence for a link between age of the TF, as well as polymorphism in their DNA-contacting residues and expression levels-both of which may be jointly affected by selection.

piRNAs Are Associated with Diverse Transgenerational Effects on Gene and Transposon Expression in a Hybrid Dysgenic Syndrome of D. virilis

Sexual reproduction allows transposable elements (TEs) to proliferate, leading to rapid divergence between populations and species. A significant outcome of divergence in the TE landscape is evident in hybrid dysgenic syndromes, a strong form of genomic incompatibility that can arise when (TE) family abundance differs between two parents. When TEs inherited from the father are absent in the mother's genome, TEs can become activated in the progeny, causing germline damage and sterility. Studies in Drosophila indicate that dysgenesis can occur when TEs inherited paternally are not matched with a pool of corresponding TE silencing PIWI-interacting RNAs (piRNAs) provisioned by the female germline. Using the D. virilis syndrome of hybrid dysgenesis as a model, we characterize the effects that divergence in TE profile between parents has on offspring. Overall, we show that divergence in the TE landscape is associated with persisting differences in germline TE expression when comparing genetically identical females of reciprocal crosses and these differences are transmitted to the next generation. Moreover, chronic and persisting TE expression coincides with increased levels of genic piRNAs associated with reduced gene expression. Combined with these effects, we further demonstrate that gene expression is idiosyncratically influenced by differences in the genic piRNA profile of the parents that arise though polymorphic TE insertions. Overall, these results support a model in which early germline events in dysgenesis establish a chronic, stable state of both TE and gene expression in the germline that is maintained through adulthood and transmitted to the next generation. This work demonstrates that divergence in the TE profile is associated with diverse piRNA-mediated transgenerational effects on gene expression within populations.

Transposable elements (TEs) are selfish elements that copy themselves. More than half of the human genome is comprised of such elements. Studies in the fruit flies Drosophila melanogaster and D. virilis have been important in demonstrating a role for RNA silencing by PIWI-interacting RNAs (piRNAs) in protecting the genome against these harmful elements. These small RNAs are capable of recognizing TE mRNAs and mediating their destruction. They are also transmitted by the female germline to offspring in order to maintain a stable genome across generations. When males carrying a particular TE family are crossed with females lacking the element, the mother is unable to provide genome defense via complementary piRNAs that target the element. This leads to excess TE activation in the germline and sterility, a phenomenon known as hybrid dysgenesis. In this article we characterize the genomic landscape of TE destabilization that occurs in dysgenic crosses of D. virilis. We demonstrate that this mobilization is associated with an increased level of germline TE expression that persists through adulthood. In addition, we find that TE activation is associated with diverse effects on normal gene expression that are also mediated by piRNAs.

Genome-wide analysis of signatures of selection in populations of African honey bees (Apis mellifera) using new web-based tools

Background

With the development of inexpensive, high-throughput sequencing technologies, it has become feasible to examine questions related to population genetics and molecular evolution of non-model species in their ecological contexts on a genome-wide scale. Here, we employed a newly developed suite of integrated, web-based programs to examine population dynamics and signatures of selection across the genome using several well-established tests, including F_ST, pN/pS, and McDonald-Kreitman. We applied these techniques to study populations of honey bees (Apis mellifera) in East Africa. In Kenya, there are several described A. mellifera subspecies, which are thought to be localized to distinct ecological regions.

Results

We performed whole genome sequencing of 11 worker honey bees from apiaries distributed throughout Kenya and identified 3.6 million putative single-nucleotide polymorphisms. The dense coverage allowed us to apply several computational procedures to study population structure and the evolutionary relationships among the populations, and to detect signs of adaptive evolution across the genome. While there is considerable gene flow among the sampled populations, there are clear distinctions between populations from the northern desert region and those from the temperate, savannah region. We identified several genes showing population genetic patterns consistent with positive selection within African bee populations, and between these populations and European A. mellifera or Asian Apis florea.

Conclusions

These results lay the groundwork for future studies of adaptive ecological evolution in honey bees, and demonstrate the use of new, freely available web-based tools and workflows (http://usegalaxy.org/r/kenyanbee) that can be applied to any model system with genomic information.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1712-0) contains supplementary material, which is available to authorized users.

Variation in piRNA and transposable element content in strains of Drosophila melanogaster

Transposable elements (TEs) are one of the most important features of genome architecture, so their evolution and relationship with host defense mechanisms have been topics of intense study, especially in model systems such as Drosophila melanogaster. Recently, a novel small RNA-based defense mechanism in animals called the Piwi-interacting RNA (piRNA) pathway was discovered to form an adaptive defense mechanism against TEs. To investigate the relationship between piRNA and TE content between strains of a species, we sequenced piRNAs from 16 inbred lines of D. melanogaster from the Drosophila Genetic Reference Panel. Instead of a global correlation of piRNA expression and TE content, we found evidence for a host response through de novo piRNA production from novel TE insertions. Although approximately 20% of novel TE insertions induced de novo piRNA production, the abundance of de novo piRNAs was low and did not markedly affect the global pool of ovarian piRNAs. Our results provide new insights into the evolution of TEs and the piRNA system in an important model organism.

Rapid evolution of piRNA pathway in the teleost fish: implication for an adaptation to transposon diversity

The Piwi-interacting RNA (piRNA) pathway is responsible for germline specification, gametogenesis, transposon silencing, and genome integrity. Transposable elements can disrupt genome and its functions. However, piRNA pathway evolution and its adaptation to transposon diversity in the teleost fish remain unknown. This article unveils evolutionary scene of piRNA pathway and its association with diverse transposons by systematically comparative analysis on diverse teleost fish genomes. Selective pressure analysis on piRNA pathway and miRNA/siRNA (microRNA/small interfering RNA) pathway genes between teleosts and mammals showed an accelerated evolution of piRNA pathway genes in the teleost lineages, and positive selection on functional PAZ (Piwi/Ago/Zwille) and Tudor domains involved in the Piwi-piRNA/Tudor interaction, suggesting that the amino acid substitutions are adaptive to their functions in piRNA pathway in the teleost fish species. Notably five piRNA pathway genes evolved faster in the swamp eel, a kind of protogynous hermaphrodite fish, than the other teleosts, indicating a differential evolution of piRNA pathway between the swamp eel and other gonochoristic fishes. In addition, genome-wide analysis showed higher diversity of transposons in the teleost fish species compared with mammals. Our results suggest that rapidly evolved piRNA pathway in the teleost fish is likely to be involved in the adaption to transposon diversity.