Ectopic assembly of heterochromatin in Drosophila melanogaster triggered by transposable elements

Varying recombination landscapes between individuals are driven by polymorphic transposable elements

Meiotic recombination is a prominent force shaping genome evolution, and understanding the causes for varying recombination landscapes within and between species has remained a central, though challenging, question. Recombination rates are widely observed to negatively associate with the abundance of transposable elements (TEs), selfish genetic elements that move between genomic locations. While such associations are usually interpreted as recombination influencing the efficacy of selection at removing TEs, accumulating findings suggest that TEs could instead be the cause rather than the consequence. To test this prediction, we formally investigated the influence of polymorphic, putatively active TEs on recombination rates. We developed and benchmarked a novel approach that uses PacBio long-read sequencing to efficiently, accurately, and cost-effectively identify crossovers (COs), a key recombination product, among large numbers of pooled recombinant individuals. By applying this approach to Drosophila strains with distinct TE insertion profiles, we found that polymorphic TEs, especially RNA-based TEs and TEs with local enrichment of repressive marks, reduce the occurrence of COs. Such an effect leads to different CO frequencies between homologous sequences with and without TEs, contributing to varying CO maps between individuals. The suppressive effect of TEs on CO is further supported by two orthogonal approaches-analyzing the distributions of COs in panels of recombinant inbred lines in relation to TE polymorphism and applying marker-assisted estimations of CO frequencies to isogenic strains with and without transgenically inserted TEs. Our investigations reveal how the constantly changing mobilome can actively modify recombination landscapes, shaping genome evolution within and between species.

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.

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.

Species-specific chromatin landscape determines how transposable elements shape genome evolution

Transposable elements (TEs) are selfish genetic parasites that increase their copy number at the expense of host fitness. The ‘success’, or genome-wide abundance, of TEs differs widely between species. Deciphering the causes for this large variety in TE abundance has remained a central question in evolutionary genomics. We previously proposed that species-specific TE abundance could be driven by the inadvertent consequences of host-direct epigenetic silencing of TEs—the spreading of repressive epigenetic marks from silenced TEs into adjacent sequences. Here, we compared this TE-mediated local enrichment of repressive marks, or ‘the epigenetic effect of TEs’, in six species in the Drosophila melanogaster subgroup to dissect step-by-step the role of such effect in determining genomic TE abundance. We found that TE-mediated local enrichment of repressive marks is prevalent and substantially varies across and even within species. While this TE-mediated effect alters the epigenetic states of adjacent genes, we surprisingly discovered that the transcription of neighboring genes could reciprocally impact this spreading. Importantly, our multi-species analysis provides the power and appropriate phylogenetic resolution to connect species-specific host chromatin regulation, TE-mediated epigenetic effects, the strength of natural selection against TEs, and genomic TE abundance unique to individual species. Our findings point toward the importance of host chromatin landscapes in shaping genome evolution through the epigenetic effects of a selfish genetic parasite.

All the instructions required for life are encoded in the set of DNA present in a cell. It therefore seems natural to think that every bit of this genetic information should serve the organism. And yet most species carry parasitic ‘transposable’ sequences, or transposons, whose only purpose is to multiply and insert themselves at other positions in the genome.

It is possible for cells to suppress these selfish elements. Chemical marks can be deposited onto the DNA to temporarily ‘silence’ transposons and prevent them from being able to move and replicate. However, this sometimes comes at a cost: the repressive chemical modifications can spread to nearby genes that are essential for the organism and perturb their function.

Strangely, the prevalence of transposons varies widely across the tree of life. These sequences form the majority of the genome of certain species – in fact, they represent about half of the human genetic information. But their abundance is much lower in other organisms, forming a measly 6% of the genome of puffer fish for instance. Even amongst fruit fly species, the prevalence of transposable elements can range between 2% and 25%. What explains such differences?

Huang et al. set out to examine this question through the lens of transposon silencing, systematically comparing how this process impacts nearby regions in six species of fruit flies. This revealed variations in the strength of the side effects associated with transposon silencing, resulting in different levels of perturbation on neighbouring genes. A stronger impact was associated with the species having fewer transposons in its genome, suggesting that an evolutionary pressure is at work to keep the abundance of transposons at a low level in these species. Further analyses showed that the genes which determine how silencing marks are distributed may also be responsible for the variations in the impact of transposon silencing. They could therefore be the ones driving differences in the abundance of transposons between species.

Overall, this work sheds light on the complex mechanisms shaping the evolution of genomes, and it may help to better understand how transposons are linked to processes such as aging and cancer.

All Quiet on the TE Front? The Role of Chromatin in Transposable Element Silencing

Transposable elements (TEs) are mobile genetic elements that constitute a sizeable portion of many eukaryotic genomes. Through their mobility, they represent a major source of genetic variation, and their activation can cause genetic instability and has been linked to aging, cancer and neurodegenerative diseases. Accordingly, tight regulation of TE transcription is necessary for normal development. Chromatin is at the heart of TE regulation; however, we still lack a comprehensive understanding of the precise role of chromatin marks in TE silencing and how chromatin marks are established and maintained at TE loci. In this review, I discuss evidence documenting the contribution of chromatin-associated proteins and histone marks in TE regulation across different species with an emphasis on Drosophila and mammalian systems.

Meiotic, genomic and evolutionary properties of crossover distribution in Drosophila yakuba

The number and location of crossovers across genomes are highly regulated during meiosis, yet the key components controlling them are fast evolving, hindering our understanding of the mechanistic causes and evolutionary consequences of changes in crossover rates. Drosophila melanogaster has been a model species to study meiosis for more than a century, with an available high-resolution crossover map that is, nonetheless, missing for closely related species, thus preventing evolutionary context. Here, we applied a novel and highly efficient approach to generate whole-genome high-resolution crossover maps in D. yakuba to tackle multiple questions that benefit from being addressed collectively within an appropriate phylogenetic framework, in our case the D. melanogaster species subgroup. The genotyping of more than 1,600 individual meiotic events allowed us to identify several key distinct properties relative to D. melanogaster. We show that D. yakuba, in addition to higher crossover rates than D. melanogaster, has a stronger centromere effect and crossover assurance than any Drosophila species analyzed to date. We also report the presence of an active crossover-associated meiotic drive mechanism for the X chromosome that results in the preferential inclusion in oocytes of chromatids with crossovers. Our evolutionary and genomic analyses suggest that the genome-wide landscape of crossover rates in D. yakuba has been fairly stable and captures a significant signal of the ancestral crossover landscape for the whole D. melanogaster subgroup, even informative for the D. melanogaster lineage. Contemporary crossover rates in D. melanogaster, on the other hand, do not recapitulate ancestral crossovers landscapes. As a result, the temporal stability of crossover landscapes observed in D. yakuba makes this species an ideal system for applying population genetic models of selection and linkage, given that these models assume temporal constancy in linkage effects. Our studies emphasize the importance of generating multiple high-resolution crossover rate maps within a coherent phylogenetic context to broaden our understanding of crossover control during meiosis and to improve studies on the evolutionary consequences of variable crossover rates across genomes and time.

Complex Genetic Interactions between Piwi and HP1a in the Repression of Transposable Elements and Tissue-Specific Genes in the Ovarian Germline

Insertions of transposable elements (TEs) in eukaryotic genomes are usually associated with repressive chromatin, which spreads to neighbouring genomic sequences. In ovaries of Drosophila melanogaster, the Piwi-piRNA pathway plays a key role in the transcriptional silencing of TEs considered to be exerted mostly through the establishment of H3K9me3 histone marks recruiting Heterochromatin Protein 1a (HP1a). Here, using RNA-seq, we investigated the expression of TEs and the adjacent genomic regions upon Piwi and HP1a germline knockdowns sharing a similar genetic background. We found that the depletion of Piwi and HP1a led to the derepression of only partially overlapping TE sets. Several TEs were silenced predominantly by HP1a, whereas the upregulation of some other TEs was more pronounced upon Piwi knockdown and, surprisingly, was diminished upon a Piwi/HP1a double-knockdown. We revealed that HP1a loss influenced the expression of thousands of protein-coding genes mostly not adjacent to TE insertions and, in particular, downregulated a putative transcriptional factor required for TE activation. Nevertheless, our results indicate that Piwi and HP1a cooperatively exert repressive effects on the transcription of euchromatic loci flanking the insertions of some Piwi-regulated TEs. We suggest that this mechanism controls the silencing of a small set of TE-adjacent tissue-specific genes, preventing their inappropriate expression in ovaries.

Hinfp is a guardian of the somatic genome by repressing transposable elements

Germ cells possess the Piwi-interacting RNA pathway to repress transposable elements and maintain genome stability across generations. Transposable element mobilization in somatic cells does not affect future generations, but nonetheless can lead to pathological outcomes in host tissues. We show here that loss of function of the conserved zinc-finger transcription factor Hinfp causes dysregulation of many host genes and derepression of most transposable elements. There is also substantial DNA damage in somatic tissues of Drosophila after loss of Hinfp. Interference of transposable element mobilization by reverse-transcriptase inhibitors can suppress some of the DNA damage phenotypes. The key cell-autonomous target of Hinfp in this process is Histone1, which encodes linker histones essential for higher-order chromatin assembly. Transgenic expression of Hinfp or Histone1, but not Histone4 of core nucleosome, is sufficient to rescue the defects in repressing transposable elements and host genes. Loss of Hinfp enhances Ras-induced tissue growth and aging-related phenotypes. Therefore, Hinfp is a physiological regulator of Histone1-dependent silencing of most transposable elements, as well as many host genes, and serves as a venue for studying genome instability, cancer progression, neurodegeneration, and aging.

Population Genomics on the Fly: Recent Advances in Drosophila

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

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.

Adaptation of gene loci to heterochromatin in the course of Drosophila evolution is associated with insulator proteins

Pericentromeric heterochromatin is generally composed of repetitive DNA forming a transcriptionally repressive environment. Dozens of genes were embedded into pericentromeric heterochromatin during evolution of Drosophilidae lineage while retaining activity. However, factors that contribute to insusceptibility of gene loci to transcriptional silencing remain unknown. Here, we find that the promoter region of genes that can be embedded in both euchromatin and heterochromatin exhibits a conserved structure throughout the Drosophila phylogeny and carries motifs for binding of certain chromatin remodeling factors, including insulator proteins. Using ChIP-seq data, we demonstrate that evolutionary gene relocation between euchromatin and pericentric heterochromatin occurred with preservation of sites of insulation of BEAF-32 in evolutionarily distant species, i.e. D. melanogaster and D. virilis. Moreover, promoters of virtually all protein-coding genes located in heterochromatin in D. melanogaster are enriched with insulator proteins BEAF-32, GAF and dCTCF. Applying RNA-seq of a BEAF-32 mutant, we show that the impairment of BEAF-32 function has a complex effect on gene expression in D. melanogaster, affecting even those genes that lack BEAF-32 association in their promoters. We propose that conserved intrinsic properties of genes, such as sites of insulation near the promoter regions, may contribute to adaptation of genes to the heterochromatic environment and, hence, facilitate the evolutionary relocation of genes loci between euchromatin and heterochromatin.

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

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

Various modes of HP1a interactions with the euchromatic chromosome arms in Drosophila ovarian somatic cells

Heterochromatin protein 1a (HP1a) is a well-known component of pericentromeric and telomeric heterochromatin in Drosophila. However, its role and the mechanisms of its binding in the chromosome arms (ChAs) remain largely unclear. Here, we identified HP1a-interacting domains in the somatic cells of Drosophila ovaries using a DamID-seq approach and compared them with insertion sites of transposable elements (TEs) revealed by genome sequencing. Although HP1a domains cover only 13% of ChAs, they non-randomly associate with 42% of TE insertions. Furthermore, HP1a on average propagates at 2-kb distances from the TE insertions. These data confirm the role of TEs in formation of HP1a islands in ChAs. However, only 18% of HP1a domains have adjacent TEs, indicating the existence of other mechanisms of HP1a domain formation besides spreading from TEs. In particular, many TE-independent HP1a domains correspond to the regions attached to the nuclear pore complexes (NPCs) or contain active gene promoters. However, HP1a occupancy on the promoters does not significantly influence expression of corresponding genes. At the same time, the steady-state transcript level of many genes located outside of HP1a domains was altered upon HP1a knockdown in the somatic cells of ovaries, thus pointing to the strong indirect effect of HP1a depletion. Collectively, our results support an existence of at least three different mechanisms of HP1a domain emergence in ChAs: spreading from TE insertions, transient interactions with the chromatin located near NPCs, and targeting to the promoters of moderately expressed genes.

Taming the Turmoil Within: New Insights on the Containment of Transposable Elements

Transposable elements (TEs) are mobile genetic parasites that can exponentially increase their genomic abundance through self-propagation. Classic theoretical papers highlighted the importance of two potentially escalating forces that oppose TE spread: regulated transposition and purifying selection. Here, we review new insights into mechanisms of TE regulation and purifying selection, which reveal the remarkable foresight of these theoretical models. We further highlight emergent connections between transcriptional control enacted by small RNAs and the contribution of TE insertions to structural mutation and host-gene regulation. Finally, we call for increased comparative analysis of TE dynamics and fitness effects, as well as host control mechanisms, to reveal how interconnected forces shape the differential prevalence and distribution of TEs across the tree of life.

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

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

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

The Telomere Paradox: Stable Genome Preservation with Rapidly Evolving Proteins

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

"What You Need, Baby, I Got It": Transposable Elements as Suppliers of Cis-Operating Sequences in Drosophila

Transposable elements (TEs) are constitutive components of both eukaryotic and prokaryotic genomes. The role of TEs in the evolution of genes and genomes has been widely assessed over the past years in a variety of model and non-model organisms. Drosophila is undoubtedly among the most powerful model organisms used for the purpose of studying the role of transposons and their effects on the stability and evolution of genes and genomes. Besides their most intuitive role as insertional mutagens, TEs can modify the transcriptional pattern of host genes by juxtaposing new cis-regulatory sequences. A key element of TE biology is that they carry transcriptional control elements that fine-tune the transcription of their own genes, but that can also perturb the transcriptional activity of neighboring host genes. From this perspective, the transposition-mediated modulation of gene expression is an important issue for the short-term adaptation of physiological functions to the environmental changes, and for long-term evolutionary changes. Here, we review the current literature concerning the regulatory and structural elements operating in cis provided by TEs in Drosophila. Furthermore, we highlight that, besides their influence on both TEs and host genes expression, they can affect the chromatin structure and epigenetic status as well as both the chromosome's structure and stability. It emerges that Drosophila is a good model organism to study the effect of TE-linked regulatory sequences, and it could help future studies on TE-host interactions in any complex eukaryotic genome.

The evolutionary arms race between transposable elements and piRNAs in Drosophila melanogaster

BACKGROUND

The piwi-interacting RNAs (piRNAs) are small non-coding RNAs that specifically repress transposable elements (TEs) in the germline of Drosophila. Despite our expanding understanding of TE:piRNA interaction, whether there is an evolutionary arms race between TEs and piRNAs was unclear.

RESULTS

Here, we studied the population genomics of TEs and piRNAs in the worldwide strains of D. melanogaster. By conducting a correlation analysis between TE contents and the abundance of piRNAs from ovaries of representative strains of D. melanogaster, we find positive correlations between TEs and piRNAs in six TE families. Our simulations further highlight that TE activities and the strength of purifying selection against TEs are important factors shaping the interactions between TEs and piRNAs. Our studies also suggest that the de novo generation of piRNAs is an important mechanism to repress the newly invaded TEs.

CONCLUSIONS

Our results revealed the existence of an evolutionary arms race between the copy numbers of TEs and the abundance of antisense piRNAs at the population level. Although the interactions between TEs and piRNAs are complex and many factors should be considered to impact their interaction dynamics, our results suggest the emergence, repression specificity and strength of piRNAs on TEs should be considered in studying the landscapes of TE insertions in Drosophila. These results deepen our understanding of the interactions between piRNAs and TEs, and also provide novel insights into the nature of genomic conflicts of other forms.

The SUMO Ligase Su(var)2-10 Controls Hetero- and Euchromatic Gene Expression via Establishing H3K9 Trimethylation and Negative Feedback Regulation

Сhromatin is critical for genome compaction and gene expression. On a coarse scale, the genome is divided into euchromatin, which harbors the majority of genes and is enriched in active chromatin marks, and heterochromatin, which is gene-poor but repeat-rich. The conserved molecular hallmark of heterochromatin is the H3K9me3 modification, which is associated with gene silencing. We found that in Drosophila, deposition of most of the H3K9me3 mark depends on SUMO and the SUMO ligase Su(var)2-10, which recruits the histone methyltransferase complex SetDB1/Wde. In addition to repressing repeats, H3K9me3 influences expression of both hetero- and euchromatic host genes. High H3K9me3 levels in heterochromatin are required to suppress spurious transcription and ensure proper gene expression. In euchromatin, a set of conserved genes is repressed by Su(var)2-10/SetDB1-induced H3K9 trimethylation, ensuring tissue-specific gene expression. Several components of heterochromatin are themselves repressed by this pathway, providing a negative feedback mechanism to ensure chromatin homeostasis.

Phylogenetic Shifts in Gene Body Methylation Correlate with Gene Expression and Reflect Trait Conservation

A subset of genes in plant genomes are labeled with DNA methylation specifically at CG residues. These genes, known as gene-body methylated (gbM), have a number of associated characteristics. They tend to have longer sequences, to be enriched for intermediate expression levels, and to be associated with slower rates of molecular evolution. Most importantly, gbM genes tend to maintain their level of DNA methylation between species, suggesting that this trait is under evolutionary constraint. Given the degree of conservation in gbM, we still know surprisingly little about its function in plant genomes or whether gbM is itself a target of selection. To address these questions, we surveyed DNA methylation across eight grass (Poaceae) species that span a gradient of genome sizes. We first established that genome size correlates with genome-wide DNA methylation levels, but less so for genic levels. We then leveraged genomic data to identify a set of 2,982 putative orthologs among the eight species and examined shifts of methylation status for each ortholog in a phylogenetic context. A total of 55% of orthologs exhibited a shift in gbM, but these shifts occurred predominantly on terminal branches, indicating that shifts in gbM are rarely conveyed over time. Finally, we found that the degree of conservation of gbM across species is associated with increased gene length, reduced rates of molecular evolution, and increased gene expression level, but reduced gene expression variation across species. Overall, these observations suggest a basis for evolutionary pressure to maintain gbM status over evolutionary time.

Aging in the Drosophila ovary: contrasting changes in the expression of the piRNA machinery and mitochondria but no global release of transposable elements

Background

Evolutionary theory indicates that the dynamics of aging in the soma and reproductive tissues may be distinct. This difference arises from the fact that only the germline lineage establishes future generations. In the soma, changes in the landscape of heterochromatin have been proposed to have an important role in aging. This is because redistribution of heterochromatin during aging has been linked to the derepression of transposable elements and an overall loss of somatic gene regulation. A role for changes in the chromatin landscape in the aging of reproductive tissues is less well established. Whether or not epigenetic factors, such as heterochromatin marks, are perturbed in aging reproductive tissues is of interest because, in special cases, epigenetic variation may be heritable. Using mRNA sequencing data from late-stage egg chambers in Drosophila melanogaster, we characterized the landscape of altered gene and transposable element expression in aged reproductive tissues. This allowed us to test the hypothesis that reproductive tissues may differ from somatic tissues in their response to aging.

Results

We show that age-related expression changes in late-stage egg chambers tend to occur in genes residing in heterochromatin, particularly on the largely heterochromatic 4th chromosome. However, these expression differences are seen as both decreases and increases during aging, inconsistent with a general loss of heterochromatic silencing. We also identify an increase in expression of the piRNA machinery, suggesting an age-related increased investment in the maintenance of genome stability. We further identify a strong age-related reduction in the expression of mitochondrial transcripts. However, we find no evidence for global TE derepression in reproductive tissues. Rather, the observed effects of aging on TEs are primarily strain and family specific.

Conclusions

These results identify unique responses in somatic versus reproductive tissue with regards to aging. As in somatic tissues, female reproductive tissues show reduced expression of mitochondrial genes. In contrast, the piRNA machinery shows increased expression during aging. Overall, these results also indicate that global loss of TE control observed in other studies may be unique to the soma and sensitive to genetic background and TE family.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5668-3) contains supplementary material, which is available to authorized users.

On the Importance to Acknowledge Transposable Elements in Epigenomic Analyses

Eukaryotic genomes comprise a large proportion of repeated sequences, an important fraction of which are transposable elements (TEs). TEs are mobile elements that have a significant impact on genome evolution and on gene functioning. Although some TE insertions could provide adaptive advantages to species, transposition is a highly mutagenic event that has to be tightly controlled to ensure its viability. Genomes have evolved sophisticated mechanisms to control TE activity, the most important being epigenetic silencing. However, the epigenetic control of TEs can also affect genes located nearby that can become epigenetically regulated. It has been proposed that the combination of TE mobilization and the induced changes in the epigenetic landscape could allow a rapid phenotypic adaptation to global environmental changes. In this review, we argue the crucial need to take into account the repeated part of genomes when studying the global impact of epigenetic modifications on an organism. We emphasize more particularly why it is important to carefully consider TEs and what bioinformatic tools can be used to do so.

New Insights on the Evolution of Genome Content: Population Dynamics of Transposable Elements in Flies and Humans

Understanding the abundance, diversity, and distribution of TEs in genomes is crucial to understand genome structure, function, and evolution. Advances in whole-genome sequencing techniques, as well as in bioinformatics tools, have increased our ability to detect and analyze the transposable element content in genomes. In addition to reference genomes, we now have access to population datasets in which multiple individuals within a species are sequenced. In this chapter, we highlight the recent advances in the study of TE population dynamics focusing on fruit flies and humans, which represent two extremes in terms of TE abundance, diversity, and activity. We review the most recent methodological approaches applied to the study of TE dynamics as well as the new knowledge on host factors involved in the regulation of TE activity. In addition to transposition rates, we also focus on TE deletion rates and on the selective forces that affect the dynamics of TEs in genomes.

The Drosophila Dot Chromosome: Where Genes Flourish Amidst Repeats

The F element of the Drosophila karyotype (the fourth chromosome in Drosophila melanogaster) is often referred to as the "dot chromosome" because of its appearance in a metaphase chromosome spread. This chromosome is distinct from other Drosophila autosomes in possessing both a high level of repetitious sequences (in particular, remnants of transposable elements) and a gene density similar to that found in the other chromosome arms, ∼80 genes distributed throughout its 1.3-Mb "long arm." The dot chromosome is notorious for its lack of recombination and is often neglected as a consequence. This and other features suggest that the F element is packaged as heterochromatin throughout. F element genes have distinct characteristics (e.g, low codon bias, and larger size due both to larger introns and an increased number of exons), but exhibit expression levels comparable to genes found in euchromatin. Mapping experiments show the presence of appropriate chromatin modifications for the formation of DNaseI hypersensitive sites and transcript initiation at the 5' ends of active genes, but, in most cases, high levels of heterochromatin proteins are observed over the body of these genes. These various features raise many interesting questions about the relationships of chromatin structures with gene and chromosome function. The apparent evolution of the F element as an autosome from an ancestral sex chromosome also raises intriguing questions. The findings argue that the F element is a unique chromosome that occupies its own space in the nucleus. Further study of the F element should provide new insights into chromosome structure and function.

The Evolution of Small-RNA-Mediated Silencing of an Invading Transposable Element

Transposable elements (TEs) are genomic parasites that impose fitness costs on their hosts by producing deleterious mutations and disrupting gametogenesis. Host genomes avoid these costs by regulating TE activity, particularly in germline cells where new insertions are heritable and TEs are exceptionally active. However, the capacity of different TE-associated fitness costs to select for repression in the host, and the role of selection in the evolution of TE regulation more generally remain controversial. In this study, we use forward, individual-based simulations to examine the evolution of small-RNA-mediated TE regulation, a conserved mechanism for TE repression that is employed by both prokaryotes and eukaryotes. To design and parameterize a biologically realistic model, we drew on an extensive survey of empirical studies of the transposition and regulation of P-element DNA transposons in Drosophila melanogaster. We observed that even under conservative assumptions, where small-RNA-mediated regulation reduces transposition only, repression evolves rapidly and adaptively after the genome is invaded by a new TE in simulated populations. We further show that the spread of repressor alleles through simulated populations is greatly enhanced by two additional TE-imposed fitness costs: dysgenic sterility and ectopic recombination. Finally, we demonstrate that the adaptive mutation rate to repression is a critical parameter that influences both the evolutionary trajectory of host repression and the associated proliferation of TEs after invasion in simulated populations. Our findings suggest that adaptive evolution of TE regulation may be stronger and more prevalent than previously appreciated, and provide a framework for interpreting empirical data.

Stress affects the epigenetic marks added by natural transposable element insertions in Drosophila melanogaster

Transposable elements are emerging as an important source of cis-acting regulatory sequences and epigenetic marks that could influence gene expression. However, few studies have dissected the role of specific transposable element insertions on epigenetic gene regulation. Bari-Jheh is a natural transposon that mediates resistance to oxidative stress by adding cis-regulatory sequences that affect expression of nearby genes. In this work, we integrated publicly available ChIP-seq and piRNA data with chromatin immunoprecipitation experiments to get a more comprehensive picture of Bari-Jheh molecular effects. We showed that Bari-Jheh was enriched for H3K9me3 in nonstress conditions, and for H3K9me3, H3K4me3 and H3K27me3 in oxidative stress conditions, which is consistent with expression changes in adjacent genes. We further showed that under oxidative stress conditions, H3K4me3 and H3K9me3 spread to the promoter region of Jheh1 gene. Finally, another insertion of the Bari1 family was associated with increased H3K27me3 in oxidative stress conditions suggesting that Bari1 histone marks are copy-specific. We concluded that besides adding cis-regulatory sequences, Bari-Jheh influences gene expression by affecting the local chromatin state.

Coevolution between transposable elements and recombination

One of the most striking patterns of genome structure is the tight, typically negative, association between transposable elements (TEs) and meiotic recombination rates. While this is a highly recurring feature of eukaryotic genomes, the mechanisms driving correlations between TEs and recombination remain poorly understood, and distinguishing cause versus effect is challenging. Here, we review the evidence for a relation between TEs and recombination, and discuss the underlying evolutionary forces. Evidence to date suggests that overall TE densities correlate negatively with recombination, but the strength of this correlation varies across element types, and the pattern can be reversed. Results suggest that heterogeneity in the strength of selection against ectopic recombination and gene disruption can drive TE accumulation in regions of low recombination, but there is also strong evidence that the regulation of TEs can influence local recombination rates. We hypothesize that TE insertion polymorphism may be important in driving within-species variation in recombination rates in surrounding genomic regions. Furthermore, the interaction between TEs and recombination may create positive feedback, whereby TE accumulation in non-recombining regions contributes to the spread of recombination suppression. Further investigation of the coevolution between recombination and TEs has important implications for our understanding of the evolution of recombination rates and genome structure.This article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.

Chromatin-associated noncoding RNAs in development and inheritance

Noncoding RNAs (ncRNAs) have emerged as crucial players in chromatin regulation. Their diversity allows them to partake in the regulation of numerous cellular processes across species. During development, long and short ncRNAs act in conjunction with each other where long ncRNAs (lncRNAs) are best understood in establishing appropriate gene expression patterns, while short ncRNAs (sRNAs) are known to establish constitutive heterochromatin and suppress mobile elements. Additionally, increasing evidence demonstrates roles of sRNAs in several typically lncRNA-mediated processes such as dosage compensation, indicating a complex regulatory network of noncoding RNAs. Together, various ncRNAs establish many mitotically heritable epigenetic marks during development. Additionally, they participate in mechanisms that regulate maintenance of these epigenetic marks during the lifespan of the organism. Interestingly, some epigenetic traits are transmitted to the next generation(s) via paramutations or transgenerational inheritance mediated by sRNAs. In this review, we give an overview of the various functions and regulations of ncRNAs and the mechanisms they employ in the establishment and maintenance of epigenetic marks and multi-generational transmission of epigenetic traits. WIREs RNA 2017, 8:e1435. doi: 10.1002/wrna.1435 For further resources related to this article, please visit the WIREs website.

Pervasive epigenetic effects of Drosophila euchromatic transposable elements impact their evolution

Transposable elements (TEs) are widespread genomic parasites, and their evolution has remained a critical question in evolutionary genomics. Here, we study the relatively unexplored epigenetic impacts of TEs and provide the first genome-wide quantification of such effects in D. melanogaster and D. simulans. Surprisingly, the spread of repressive epigenetic marks (histone H3K9me2) to nearby DNA occurs at >50% of euchromatic TEs, and can extend up to 20 kb. This results in differential epigenetic states of genic alleles and, in turn, selection against TEs. Interestingly, the lower TE content in D. simulans compared to D. melanogaster correlates with stronger epigenetic effects of TEs and higher levels of host genetic factors known to promote epigenetic silencing. Our study demonstrates that the epigenetic effects of euchromatic TEs, and host genetic factors modulating such effects, play a critical role in the evolution of TEs both within and between species.

DOI:http://dx.doi.org/10.7554/eLife.25762.001

The DNA inside an organism encodes all the instructions needed for the organism to develop and work properly. Organisms carefully organize and maintain their DNA (collectively known as the genome) so that the genetic information remains intact and the cell can understand the instructions. However, there are some pieces of DNA that are capable of moving around the genome. For example, pieces known as transposable elements can make new copies of themselves and jump into new locations in the genome. Most transposons do not appear to have any important roles, and in fact they are usually harmful to organisms. Despite this, transposons are present in the genomes of almost all species. The number of transposons in a genome varies greatly between individuals and species, but it is not clear why this is the case.

Organisms have evolved ways to limit the damage caused by transposons. For example, many cells package regions of DNA containing transposons into a tightly packed structure known as heterochromatin. However, this type of DNA packaging sometimes spreads to neighboring sections of DNA. This is a problem because cells are not usually able to read the information contained within heterochromatin. This means that transposons can prevent some instructions from being produced when they should be. Lee and Karpen used fruit flies to investigate to what extent transposons harm organisms by changing the way DNA is packaged, and whether this influences how transposons evolve.

The experiments show that that more than half of the transposons in fruit flies cause neighboring sections of DNA to be packaged into heterochromatin. This can negatively impact up to 20% of genes in the genome. As a result, transposons that have harmful effects on DNA packaging are more likely to be lost from the fly population during evolution than transposons that do not have harmful effects. Fruit fly species containing transposons that tend to package more neighboring sections of DNA into heterochromatin generally have fewer transposons than genomes containing less harmful transposons.

The findings of Lee and Karpen provide new insight as to why the numbers of transposons vary among organisms. The next challenge is to find out whether transposons that alter how DNA is packaged are also common in primates and other animals.

DOI:http://dx.doi.org/10.7554/eLife.25762.002

Silencing of Transposable Elements by piRNAs in Drosophila: An Evolutionary Perspective

Transposable elements (TEs) are DNA sequences that can move within the genome. TEs have greatly shaped the genomes, transcriptomes, and proteomes of the host organisms through a variety of mechanisms. However, TEs generally disrupt genes and destabilize the host genomes, which substantially reduce fitness of the host organisms. Understanding the genomic distribution and evolutionary dynamics of TEs will greatly deepen our understanding of the TE-mediated biological processes. Most TE insertions are highly polymorphic in Drosophila melanogaster, providing us a good system to investigate the evolution of TEs at the population level. Decades of theoretical and experimental studies have well established “transposition-selection” population genetics model, which assumes that the equilibrium between TE replication and purifying selection determines the copy number of TEs in the genome. In the last decade, P-element-induced wimpy testis (PIWI)-interacting RNAs (piRNAs) were demonstrated to be master repressors of TE activities in Drosophila. The discovery of piRNAs revolutionized our understanding of TE repression, because it reveals that the host organisms have evolved an adaptive mechanism to defend against TE invasion. Tremendous progress has been made to understand the molecular mechanisms by which piRNAs repress active TEs, although many details in this process remain to be further explored. The interaction between piRNAs and TEs well explains the molecular mechanisms underlying hybrid dysgenesis for the I-R and P-M systems in Drosophila, which have puzzled evolutionary biologists for decades. The piRNA repression pathway provides us an unparalleled system to study the co-evolutionary process between parasites and host organisms.

Molecular Population Genetics

Molecular population genetics aims to explain genetic variation and molecular evolution from population genetics principles. The field was born 50 years ago with the first measures of genetic variation in allozyme loci, continued with the nucleotide sequencing era, and is currently in the era of population genomics. During this period, molecular population genetics has been revolutionized by progress in data acquisition and theoretical developments. The conceptual elegance of the neutral theory of molecular evolution or the footprint carved by natural selection on the patterns of genetic variation are two examples of the vast number of inspiring findings of population genetics research. Since the inception of the field, Drosophila has been the prominent model species: molecular variation in populations was first described in Drosophila and most of the population genetics hypotheses were tested in Drosophila species. In this review, we describe the main concepts, methods, and landmarks of molecular population genetics, using the Drosophila model as a reference. We describe the different genetic data sets made available by advances in molecular technologies, and the theoretical developments fostered by these data. Finally, we review the results and new insights provided by the population genomics approach, and conclude by enumerating challenges and new lines of inquiry posed by increasingly large population scale sequence data.

Drosophila melanogaster retrotransposon and inverted repeat-derived endogenous siRNAs are differentially processed in distinct cellular locations

BACKGROUND

Endogenous small interfering (esi)RNAs repress mRNA levels and retrotransposon mobility in Drosophila somatic cells by poorly understood mechanisms. 21 nucleotide esiRNAs are primarily generated from retrotransposons and two inverted repeat (hairpin) loci in Drosophila culture cells in a Dicer2 dependent manner. Additionally, proteins involved in 3' end processing, such as Symplekin, CPSF73 and CPSR100, have been recently implicated in the esiRNA pathway.

RESULTS

Here we present evidence of overlap between two essential RNA metabolic pathways: esiRNA biogenesis and mRNA 3' end processing. We have identified a nucleus-specific interaction between the essential esiRNA cleavage enzyme Dicer2 (Dcr2) and Symplekin, a component of the core cleavage complex (CCC) required for 3' end processing of all eukaryotic mRNAs. This interaction is mediated by the N-terminal 271 amino acids of Symplekin; CCC factors CPSF73 and CPSF100 do not contact Dcr2. While Dcr2 binds the CCC, Dcr2 knockdown does not affect mRNA 3' end formation. RNAi-depletion of CCC components Symplekin and CPSF73 causes perturbations in esiRNA abundance that correlate with fluctuations in retrotransposon and hairpin esiRNA precursor levels. We also discovered that esiRNAs generated from retrotransposons and hairpins have distinct physical characteristics including a higher predominance of 22 nucleotide hairpin-derived esiRNAs and differences in 3' and 5' base preference. Additionally, retrotransposon precursors and derived esiRNAs are highly enriched in the nucleus while hairpins and hairpin derived esiRNAs are predominantly cytoplasmic similar to canonical mRNAs. RNAi-depletion of either CPSF73 or Symplekin results in nuclear retention of both hairpin and retrotransposon precursors suggesting that polyadenylation indirectly affects cellular localization of Dcr2 substrates.

CONCLUSIONS

Together, these observations support a novel mechanism in which differences in localization of esiRNA precursors impacts esiRNA biogenesis. Hairpin-derived esiRNAs are generated in the cytoplasm independent of Dcr2-Symplekin interactions, while retrotransposons are processed in the nucleus.

Helitrons shaping the genomic architecture of Drosophila: enrichment of DINE-TR1 in α- and β-heterochromatin, satellite DNA emergence, and piRNA expression

Drosophila INterspersed Elements (DINEs) constitute an abundant but poorly understood group of Helitrons present in several Drosophila species. The general structure of DINEs includes two conserved blocks that may or not contain a region with tandem repeats in between. These central tandem repeats (CTRs) are similar within species but highly divergent between species. It has been assumed that CTRs have independent origins. Herein, we identify a subset of DINEs, termed DINE-TR1, which contain homologous CTRs of approximately 150 bp. We found DINE-TR1 in the sequenced genomes of several Drosophila species and in Bactrocera tryoni (Acalyptratae, Diptera). However, interspecific high sequence identity (∼ 88 %) is limited to the first ∼ 30 bp of each tandem repeat, implying that evolutionary constraints operate differently over the monomer length. DINE-TR1 is unevenly distributed across the Drosophila phylogeny. Nevertheless, sequence analysis suggests vertical transmission. We found that CTRs within DINE-TR1 have independently expanded into satellite DNA-like arrays at least twice within Drosophila. By analyzing the genome of Drosophila virilis and Drosophila americana, we show that DINE-TR1 is highly abundant in pericentromeric heterochromatin boundaries, some telomeric regions and in the Y chromosome. It is also present in the centromeric region of one autosome from D. virilis and dispersed throughout several euchromatic sites in both species. We further found that DINE-TR1 is abundant at piRNA clusters, and small DINE-TR1-derived RNA transcripts (∼25 nt) are predominantly expressed in the testes and the ovaries, suggesting active targeting by the piRNA machinery. These features suggest potential piRNA-mediated regulatory roles for DINEs at local and genome-wide scales in Drosophila.

Horizontal Transfer Can Drive a Greater Transposable Element Load in Large Populations

Genomes are comprised of contrasting domains of euchromatin and heterochromatin, and transposable elements (TEs) play an important role in defining these genomic regions. Therefore, understanding the forces that control TE abundance can help us understand the chromatin landscape of the genome. What determines the burden of TEs in populations? Some have proposed that drift plays a determining role. In small populations, mildly deleterious TE insertion alleles are allowed to fix, leading to increased copy number. However, it is not clear how the rate of exposure to new TE families, via horizontal transfer (HT), can contribute to broader patterns of genomic TE abundance. Here, using simulation and analytical approaches, we show that when the effects of drift are weak, exposure rate to new TE families via HT can be an important determinant of genomic copy number. If population exposure rate is proportional to population size, larger populations are expected to have a higher rate of exposure to rare HT events. This leads to the counterintuitive prediction that larger populations may carry a higher TE load. We also find that increased rates of recombination can lead to greater probabilities of TE establishment. This work has implications for our understanding of the evolution of chromatin landscapes, genome defense by RNA silencing, and recombination rates.

Multiple Independent Retroelement Insertions in the Promoter of a Stress Response Gene Have Variable Molecular and Functional Effects in Drosophila

Promoters are structurally and functionally diverse gene regulatory regions. The presence or absence of sequence motifs and the spacing between the motifs defines the properties of promoters. Recent alternative promoter usage analyses in Drosophila melanogaster revealed that transposable elements significantly contribute to promote diversity. In this work, we analyzed in detail one of the transposable element insertions, named FBti0019985, that has been co-opted to drive expression of CG18446, a candidate stress response gene. We analyzed strains from different natural populations and we found that besides FBti0019985, there are another eight independent transposable elements inserted in the proximal promoter region of CG18446. All nine insertions are solo-LTRs that belong to the roo family. We analyzed the sequence of the nine roo insertions and we investigated whether the different insertions were functionally equivalent by performing 5'-RACE, gene expression, and cold-stress survival experiments. We found that different insertions have different molecular and functional consequences. The exact position where the transposable elements are inserted matters, as they all showed highly conserved sequences but only two of the analyzed insertions provided alternative transcription start sites, and only the FBti0019985 insertion consistently affects CG18446 expression. The phenotypic consequences of the different insertions also vary: only FBti0019985 was associated with cold-stress tolerance. Interestingly, the only previous report of transposable elements inserting repeatedly and independently in a promoter region in D. melanogaster, were also located upstream of a stress response gene. Our results suggest that functional validation of individual structural variants is needed to resolve the complexity of insertion clusters.

Comparative Genomic Analyses Provide New Insights into the Evolutionary Dynamics of Heterochromatin in Drosophila

The term heterochromatin has been long considered synonymous with gene silencing, but it is now clear that the presence of transcribed genes embedded in pericentromeric heterochromatin is a conserved feature in the evolution of eukaryotic genomes. Several studies have addressed the epigenetic changes that enable the expression of genes in pericentric heterochromatin, yet little is known about the evolutionary processes through which this has occurred. By combining genome annotation analysis and high-resolution cytology, we have identified and mapped 53 orthologs of D. melanogaster heterochromatic genes in the genomes of two evolutionarily distant species, D. pseudoobscura and D. virilis. Our results show that the orthologs of the D. melanogaster heterochromatic genes are clustered at three main genomic regions in D. virilis and D. pseudoobscura. In D. virilis, the clusters lie in the middle of euchromatin, while those in D. pseudoobscura are located in the proximal portion of the chromosome arms. Some orthologs map to the corresponding Muller C element in D. pseudoobscura and D. virilis, while others localize on the Muller B element, suggesting that chromosomal rearrangements that have been instrumental in the fusion of two separate elements involved the progenitors of genes currently located in D. melanogaster heterochromatin. These results demonstrate an evolutionary repositioning of gene clusters from ancestral locations in euchromatin to the pericentromeric heterochromatin of descendent D. melanogaster chromosomes. Remarkably, in both D. virilis and D. pseudoobscura the gene clusters show a conserved association with the HP1a protein, one of the most highly evolutionarily conserved epigenetic marks. In light of these results, we suggest a new scenario whereby ancestral HP1-like proteins (and possibly other epigenetic marks) may have contributed to the evolutionary repositioning of gene clusters into heterochromatin.

Repetitive sequences and epigenetic modification: inseparable partners play important roles in the evolution of plant sex chromosomes

The present review discusses the roles of repetitive sequences played in plant sex chromosome evolution, and highlights epigenetic modification as potential mechanism of repetitive sequences involved in sex chromosome evolution. Sex determination in plants is mostly based on sex chromosomes. Classic theory proposes that sex chromosomes evolve from a specific pair of autosomes with emergence of a sex-determining gene(s). Subsequently, the newly formed sex chromosomes stop recombination in a small region around the sex-determining locus, and over time, the non-recombining region expands to almost all parts of the sex chromosomes. Accumulation of repetitive sequences, mostly transposable elements and tandem repeats, is a conspicuous feature of the non-recombining region of the Y chromosome, even in primitive one. Repetitive sequences may play multiple roles in sex chromosome evolution, such as triggering heterochromatization and causing recombination suppression, leading to structural and morphological differentiation of sex chromosomes, and promoting Y chromosome degeneration and X chromosome dosage compensation. In this article, we review the current status of this field, and based on preliminary evidence, we posit that repetitive sequences are involved in sex chromosome evolution probably via epigenetic modification, such as DNA and histone methylation, with small interfering RNAs as the mediator.

Bioinformatic analyses of sense and antisense expression from terminal inverted repeat transposons in Drosophila somatic cells

Understanding regulation of transposon movement in somatic cells is important as mobile elements can cause detrimental genomic rearrangements. Generally, transposons move via one of 2 mechanisms; retrotransposons utilize an RNA intermediate, therefore copying themselves and amplifying throughout the genome, while terminal inverted repeat transposons (TIR Tns) excise DNA sequences from the genome and integrate into a new location. Our recently published work indicates that retrotransposons in Drosophila tissue culture cells are actively transcribed in the antisense direction. Our data support a model in which convergent transcription of retrotransposons from intra element transcription start sites results in complementary RNAs that hybridize to form substrates for Dicer-2, the endogenous small interfering (esi)RNA generating enzyme. Here, we extend our previous analysis to TIR Tns. In contrast to retrotransposons, our data show that antisense TIR Tn RNAs result from transcription of intronic TIR Tns oriented antisense to their host genes. Also, disproportionately less esiRNAs are generated from TIR transcripts than from retrotransposons and transcription of very few individual TIR Tns could be confirmed. Collectively, these data support a model in which TIR Tns are regulated at the level of Transposase production while retrotransposons are regulated with esiRNA post-transcriptional mechanisms in Drosophila somatic cells.

Heterochromatin-Associated Proteins HP1a and Piwi Collaborate to Maintain the Association of Achiasmate Homologs in Drosophila Oocytes

Accurate segregation of homologous chromosomes during meiosis depends on their ability to remain physically connected throughout prophase I. For homologs that achieve a crossover, sister chromatid cohesion distal to the chiasma keeps them attached until anaphase I. However, in Drosophila melanogaster wild-type oocytes, chromosome 4 never recombines, and the X chromosome fails to cross over in 6-10% of oocytes. Proper segregation of these achiasmate homologs relies on their pericentric heterochromatin-mediated association, but the mechanism(s) underlying this attachment remains poorly understood. Using an inducible RNA interference (RNAi) strategy combined with fluorescence in situ hybridization (FISH) to monitor centromere proximal association of the achiasmate FM7a/X homolog pair, we asked whether specific heterochromatin-associated proteins are required for the association and proper segregation of achiasmate homologs in Drosophila oocytes. When we knock down HP1a, H3K9 methytransferases, or the HP1a binding partner Piwi during mid-prophase, we observe significant disruption of pericentric heterochromatin-mediated association of FM7a/X homologs. Furthermore, for both HP1a and Piwi knockdown oocytes, transgenic coexpression of the corresponding wild-type protein is able to rescue RNAi-induced defects, but expression of a mutant protein with a single amino acid change that disrupts the HP1a-Piwi interaction is unable to do so. We show that Piwi is stably bound to numerous sites along the meiotic chromosomes, including centromere proximal regions. In addition, reduction of HP1a or Piwi during meiotic prophase induces a significant increase in FM7a/X segregation errors. We present a speculative model outlining how HP1a and Piwi could collaborate to keep achiasmate chromosomes associated in a homology-dependent manner.

Targeting of P-Element Reporters to Heterochromatic Domains by Transposable Element 1360 in Drosophila melanogaster

Heterochromatin is a common DNA packaging form employed by eukaryotes to constitutively silence transposable elements. Determining which sequences to package as heterochromatin is vital for an organism. Here, we use Drosophila melanogaster to study heterochromatin formation, exploiting position-effect variegation, a process whereby a transgene is silenced stochastically if inserted in proximity to heterochromatin, leading to a variegating phenotype. Previous studies identified the transposable element 1360 as a target for heterochromatin formation. We use transgene reporters with either one or four copies of 1360 to determine if increasing local repeat density can alter the fraction of the genome supporting heterochromatin formation. We find that including 1360 in the reporter increases the frequency with which variegating phenotypes are observed. This increase is due to a greater recovery of insertions at the telomere-associated sequences (∼50% of variegating inserts). In contrast to variegating insertions elsewhere, the phenotype of telomere-associated sequence insertions is largely independent of the presence of 1360 in the reporter. We find that variegating and fully expressed transgenes are located in different types of chromatin and that variegating reporters in the telomere-associated sequences differ from those in pericentric heterochromatin. Indeed, chromatin marks at the transgene insertion site can be used to predict the eye phenotype. Our analysis reveals that increasing the local repeat density (via the transgene reporter) does not enlarge the fraction of the genome supporting heterochromatin formation. Rather, additional copies of 1360 appear to target the reporter to the telomere-associated sequences with greater efficiency, thus leading to an increased recovery of variegating insertions.

Antisense Transcription of Retrotransposons in Drosophila: An Origin of Endogenous Small Interfering RNA Precursors

Movement of transposons causes insertions, deletions, and chromosomal rearrangements potentially leading to premature lethality in Drosophila melanogaster. To repress these elements and combat genomic instability, eukaryotes have evolved several small RNA-mediated defense mechanisms. Specifically, in Drosophila somatic cells, endogenous small interfering (esi)RNAs suppress retrotransposon mobility. EsiRNAs are produced by Dicer-2 processing of double-stranded RNA precursors, yet the origins of these precursors are unknown. We show that most transposon families are transcribed in both the sense (S) and antisense (AS) direction in Dmel-2 cells. LTR retrotransposons Dm297, mdg1, and blood, and non-LTR retrotransposons juan and jockey transcripts, are generated from intraelement transcription start sites with canonical RNA polymerase II promoters. We also determined that retrotransposon antisense transcripts are less polyadenylated than sense. RNA-seq and small RNA-seq revealed that Dicer-2 RNA interference (RNAi) depletion causes a decrease in the number of esiRNAs mapping to retrotransposons and an increase in expression of both S and AS retrotransposon transcripts. These data support a model in which double-stranded RNA precursors are derived from convergent transcription and processed by Dicer-2 into esiRNAs that silence both sense and antisense retrotransposon transcripts. Reduction of sense retrotransposon transcripts potentially lowers element-specific protein levels to prevent transposition. This mechanism preserves genomic integrity and is especially important for Drosophila fitness because mobile genetic elements are highly active.

An Interaction between RRP6 and SU(VAR)3-9 Targets RRP6 to Heterochromatin and Contributes to Heterochromatin Maintenance in Drosophila melanogaster

RNA surveillance factors are involved in heterochromatin regulation in yeast and plants, but less is known about the possible roles of ribonucleases in the heterochromatin of animal cells. Here we show that RRP6, one of the catalytic subunits of the exosome, is necessary for silencing heterochromatic repeats in the genome of Drosophila melanogaster. We show that a fraction of RRP6 is associated with heterochromatin, and the analysis of the RRP6 interaction network revealed physical links between RRP6 and the heterochromatin factors HP1a, SU(VAR)3-9 and RPD3. Moreover, genome-wide studies of RRP6 occupancy in cells depleted of SU(VAR)3-9 demonstrated that SU(VAR)3-9 contributes to the tethering of RRP6 to a subset of heterochromatic loci. Depletion of the exosome ribonucleases RRP6 and DIS3 stabilizes heterochromatic transcripts derived from transposons and repetitive sequences, and renders the heterochromatin less compact, as shown by micrococcal nuclease and proximity-ligation assays. Such depletion also increases the amount of HP1a bound to heterochromatic transcripts. Taken together, our results suggest that SU(VAR)3-9 targets RRP6 to a subset of heterochromatic loci where RRP6 degrades chromatin-associated non-coding RNAs in a process that is necessary to maintain the packaging of the heterochromatin.

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.

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.

The piwi-interacting RNAs (piRNAs) are small RNAs that can suppress the expression of selfish transposable elements (TEs) in many animal genomes. One mechanism by which piRNAs silence TEs is through the formation of heterochromatin, which is condensed chromatin and generally associated with repressed gene expression. Several functional studies have demonstrated that piRNA-mediated heterochromatin of TEs can spread to adjacent genes. We hypothesized that this spread of TE-induced heterochromatin influences the function of adjacent genes, ultimately resulting in selection against individual TEs. Consistent with our hypothesis, we found that sequences and genes adjacent to TEs are enriched in heterochromatic marks. We determine that this TE-induced variation in epigenetic status is probably piRNA-dependent and that this change in chromatin state influences the expression levels of adjacent genes. Importantly, TEs that lead to higher heterochromatin enrichment of adjacent genes are more strongly selected against, demonstrating the evolutionary consequences of TE-induced epigenetic silencing. In contrast to previously studied deleterious impacts of TEs, which depend on TEs’ physical disruptions of DNAs, our proposed functional effect of TEs is mediated through their epigenetic influence. Our study suggests that the piRNA-dependent epigenetic impact of TEs may play an important role in the evolutionary dynamics of TEs.

Exploring the Phenotypic Space and the Evolutionary History of a Natural Mutation in Drosophila melanogaster

A major challenge of modern Biology is elucidating the functional consequences of natural mutations. Although we have a good understanding of the effects of laboratory-induced mutations on the molecular- and organismal-level phenotypes, the study of natural mutations has lagged behind. In this work, we explore the phenotypic space and the evolutionary history of a previously identified adaptive transposable element insertion. We first combined several tests that capture different signatures of selection to show that there is evidence of positive selection in the regions flanking FBti0019386 insertion. We then explored several phenotypes related to known phenotypic effects of nearby genes, and having plausible connections to fitness variation in nature. We found that flies with FBti0019386 insertion had a shorter developmental time and were more sensitive to stress, which are likely to be the adaptive effect and the cost of selection of this mutation, respectively. Interestingly, these phenotypic effects are not consistent with a role of FBti0019386 in temperate adaptation as has been previously suggested. Indeed, a global analysis of the population frequency of FBti0019386 showed that climatic variables explain well the FBti0019386 frequency patterns only in Australia. Finally, although FBti0019386 insertion could be inducing the formation of heterochromatin by recruiting HP1a (Heterochromatin Protein 1a) protein, the insertion is associated with upregulation of sra in adult females. Overall, our integrative approach allowed us to shed light on the evolutionary history, the relevant fitness effects, and the likely molecular mechanisms of an adaptive mutation and highlights the complexity of natural genetic variants.