Drosophila Interspecific Hybrids Phenocopy piRNA-Pathway Mutants

The Piwi-interacting RNA (piRNA) pathway defends the germline of animals from the deleterious activity of selfish transposable elements (TEs) through small-RNA mediated silencing. Adaptation to novel invasive TEs is proposed to occur by incorporating their sequences into the piRNA pool that females produce and deposit into their eggs, which then propagates immunity against specific TEs to future generations. In support of this model, the F1 offspring of crosses between strains of the same Drosophila species sometimes suffer from germline derepression of paternally inherited TE families, caused by a failure of the maternal strain to produce the piRNAs necessary for their regulation. However, many protein components of the Drosophila piRNA pathway exhibit signatures of positive selection, suggesting that they also contribute to the evolution of host genome defense. Here we investigate piRNA pathway function and TE regulation in the F1 hybrids of interspecific crosses between D. melanogaster and D. simulans and compare them with intraspecific control crosses of D. melanogaster. We confirm previous reports showing that intraspecific crosses are characterized by derepression of paternally inherited TE families that are rare or absent from the maternal genome and piRNA pool, consistent with the role of maternally deposited piRNAs in shaping TE silencing. In contrast to the intraspecific cross, we discover that interspecific hybrids are characterized by widespread derepression of both maternally and paternally inherited TE families. Furthermore, the pattern of derepression of TE families in interspecific hybrids cannot be attributed to their paucity or absence from the piRNA pool of the maternal species. Rather, we demonstrate that interspecific hybrids closely resemble piRNA effector-protein mutants in both TE misregulation and aberrant piRNA production. We suggest that TE derepression in interspecific hybrids largely reflects adaptive divergence of piRNA pathway genes rather than species-specific differences in TE-derived piRNAs.     

Molecular evolution under increasing transposable element burden in <Emphasis Type="Italic">Drosophila</Emphasis>: A speed limit on the evolutionary arms race

Background

Genome architecture is profoundly influenced by transposable elements (TEs), and natural selection against their harmful effects is a critical factor limiting their spread. Genome defense by the piRNA silencing pathway also plays a crucial role in limiting TE proliferation. How these two forces jointly determine TE abundance is not well understood. To shed light on the nature of factors that predict TE success, we test three distinct hypotheses in the Drosophila genus. First, we determine whether TE abundance and relaxed genome-wide purifying selection on protein sequences are positively correlated. This serves to test the hypothesis that variation in TE abundance in the Drosophila genus can be explained by the strength of natural selection, relative to drift, acting in parallel against mildly deleterious non-synonymous mutations. Second, we test whether increasing TE abundance is correlated with an increased rate of amino-acid evolution in genes encoding the piRNA machinery, as might be predicted by an evolutionary arms race model. Third, we test whether increasing TE abundance is correlated with greater codon bias in genes of the piRNA machinery. This is predicted if increasing TE abundance selects for increased efficiency in the machinery of genome defense.

Results

Surprisingly, we find neither of the first two hypotheses to be true. Specifically, we found that genome-wide levels of purifying selection, measured by the ratio of non-synonymous to synonymous substitution rates (ω), were greater in species with greater TE abundance. In addition, species with greater TE abundance have greater levels of purifying selection in the piRNA machinery. In contrast, it appears that increasing TE abundance has primarily driven adaptation in the piRNA machinery by increasing codon bias.

Conclusions

These results indicate that within the Drosophila genus, a historically reduced strength of selection relative to drift is unlikely to explain patterns of increased TE success across species. Other factors, such as ecological exposure, are likely to contribute to variation in TE abundances within species. Furthermore, constraints on the piRNA machinery may temper the evolutionary arms race that would drive increasing rates of evolution at the amino acid level. In the face of these constraints, selection may act primarily by improving the translational efficiency of the machinery of genome defense through efficient codon usage.

     

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

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

The evolution of transposable elements in natural populations of self-fertilizing <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis> and its outcrossing relative <Emphasis Type="Italic">Arabidopsis lyrata</Emphasis>

Background  

Transposable Elements (TEs) make up the majority of plant genomes, and thus understanding TE evolutionary dynamics is key to understanding plant genome evolution. Plant reproductive systems are diverse and mating type variation is one factor among many hypothesized to influence TE evolutionary dynamics. Here, we collected a large TE-display data set in self-fertilizing Arabidopsis thaliana, and compared it to data gathered in outcrossing Arabidopsis lyrata. We analyzed seven TE families in four natural populations of each species to tease apart the effects of mating system, demography, transposition, and selection in determining patterns of TE diversity.     

piRNA clusters as a main source of small RNAs in the animal germline

PIWI subfamily Argonaute proteins and small RNAs bound to them (PIWI interacting RNA, piRNA) control mobilization of transposable elements (TE) in the animal germline. piRNAs are generated by distinct genomic regions termed piRNA clusters. piRNA clusters are often extensive loci enriched in damaged fragments of TEs. New TE integration into piRNA clusters causes production of TE-specific piRNAs and repression of cognate sequences. piRNAs are thought to be generated from long single-stranded precursors encoded by piRNA clusters. Special chromatin structures might be essential to distinguish these genomic loci as a source for piRNAs. In this review, we present recent findings on the structural organization of piRNA clusters and piRNA biogenesis in Drosophila and other organisms, which are important for understanding a key epigenetic mechanism that provides defense against TE expansion.     

The Role of piRNA-Mediated Epigenetic Silencing in the Population Dynamics of Transposable Elements in Drosophila melanogaster

The piwi-interacting RNAs (piRNA) are small RNAs that target selfish transposable elements (TEs) in many animal genomes. Until now, piRNAs’ role in TE population dynamics has only been discussed in the context of their suppression of TE transposition, which alone is not sufficient to account for the skewed frequency spectrum and stable containment of TEs. On the other hand, euchromatic TEs can be epigenetically silenced via piRNA-dependent heterochromatin formation and, similar to the widely known “Position-effect variegation”, heterochromatin induced by TEs can “spread” into nearby genes. We hypothesized that the piRNA-mediated spread of heterochromatin from TEs into adjacent genes has deleterious functional effects and leads to selection against individual TEs. Unlike previously identified deleterious effects of TEs due to the physical disruption of DNA, the functional effect we investigated here is mediated through the epigenetic influences of TEs. We found that the repressive chromatin mark, H3K9me, is elevated in sequences adjacent to euchromatic TEs at multiple developmental stages in Drosophila melanogaster. Furthermore, the heterochromatic states of genes depend not only on the number of and distance from adjacent TEs, but also on the likelihood that their nearest TEs are targeted by piRNAs. These variations in chromatin status probably have functional consequences, causing genes near TEs to have lower expression. Importantly, we found stronger selection against TEs that lead to higher H3K9me enrichment of adjacent genes, demonstrating the pervasive evolutionary consequences of TE-induced epigenetic silencing. Because of the intrinsic biological mechanism of piRNA amplification, spread of TE heterochromatin could result in the theoretically required synergistic deleterious effects of TE insertions for stable containment of TE copy number. The indirect deleterious impact of piRNA-mediated epigenetic silencing of TEs is a previously unexplored, yet important, element for the evolutionary dynamics of TEs.     

Epigenetics: A Historical Overview

Postmigratory mouse primordial germ cells (PGCs) undergo extensive epigenetic remodeling that includes DNA methylation (DM) reprogramming of imprinted genes and, surprisingly, of transposable elements (TEs). Given the danger posed by TEs to the integrity of the germline, even a brief derepression of TEs is counterintuitive and puzzling. In the male fetal gonocytes, a sophisticated repressive mechanism that uses DM and TE-targeting piRNAs has evolved to stably silence TEs. A recent study has further increased the complexity of this problem by revealing that TE silencing is alleviated specifically at the onset of meiosis in testes lacking MAEL, a piRNA pathway protein. These observations and prior work of others are consistent with existence of an additional reprogramming event, transient relaxation of transposon silencing (TRTS), at the onset of both male and female meiosis in mice. In this Point of View we propose that TE expression is inherent to mammalian meiosis and discuss potential functional significance of this phenomenon.     

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.     

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.     

Spatio-temporal requirements for transposable element piRNA-mediated silencing during Drosophila oogenesis

During Drosophila oogenesis, transposable element (TE) repression involves the Piwi-interacting RNA (piRNA) pathway which ensures genome integrity for the next generation. We developed a transgenic model to study repression of the Idefix retrotransposon in the germline. Using a candidate gene KD-approach, we identified differences in the spatio-temporal requirements of the piRNA pathway components for piRNA-mediated silencing. Some of them (Aub, Vasa, Spn-E) are necessary in very early stages of oogenesis within the germarium and appear to be less important for efficient TE silencing thereafter. Others (Piwi, Ago3, Mael) are required at all stages of oogenesis. Moreover, during early oogenesis, in the dividing cysts within the germarium, Idefix anti-sense transgenes escape host control, and this is associated with very low piwi expression. Silencing of P-element-based transgenes is also strongly weakened in these cysts. This region, termed the ‘Piwiless pocket’ or Pilp, may ensure that new TE insertions occur and are transmitted to the next generation, thereby contributing to genome dynamics. In contrast, piRNA-mediated silencing is strong in germline stem cells in which TE mobilization is tightly repressed ensuring the continued production of viable germline cysts.     

piRNA dynamics in divergent zebrafish strains reveal long-lasting maternal influence on zygotic piRNA profiles

Transposable elements (TEs) are mobile genetic elements that can have many deleterious effects on the fitness of their host. The germline-specific PIWI pathway guards the genome against TEs, deriving its specificity from sequence complementarity between PIWI-bound small RNAs (piRNAs) and the TEs. The piRNAs are derived from so-called piRNA clusters. Recent studies have demonstrated that the piRNA repertoire can be adjusted to accommodate recent TE invasions by capturing invading TEs in piRNA loci. Thus far, no information concerning piRNA divergence is available from vertebrates. We present piRNA analyses of two relatively divergent zebrafish strains. We find that significant differences in the piRNA populations have accumulated, most notably among active class I TEs. This divergence can be split into differences in piRNA abundance per element and differences in sense/antisense polarity ratios. In crosses between animals of the different strains, many of these differences are resolved in the progeny. However, some differences remain, often leaning to the maternally contributed piRNA population. These differences can be detected at least two generations later. Our data illustrate, for the first time, the fluidity of piRNA populations in vertebrates and how the established diversity is transmitted to future generations.