Telomere maintenance is crucial for preventing the linear eukaryotic chromosome ends from being mistaken for DNA double-strand breaks, thereby avoiding chromosome fusions and the loss of genetic material. Unlike most eukaryotes that use telomerase for telomere maintenance, Drosophila relies on retrotransposable elements-specifically HeT-A, TAHRE, and TART (collectively referred to as HTT)-which are regulated and precisely targeted to chromosome ends. Drosophila telomere protection is mediated by a set of fast-evolving proteins, termed terminin, which bind to chromosome termini without sequence specificity, balancing DNA damage response factors to avoid erroneous repair mechanisms. This unique telomere capping mechanism highlights an alternative evolutionary strategy to compensate for telomerase loss. The modulation of recombination and transcription at Drosophila telomeres offers insights into the diverse mechanisms of telomere maintenance. Recent studies at the population level have begun to reveal the architecture of telomere arrays, the diversity among the HTT subfamilies, and their relative frequencies, aiming to understand whether and how these elements have evolved to reach an equilibrium with the host and to resolve genetic conflicts. Further studies may shed light on the complex relationships between telomere transcription, recombination, and maintenance, underscoring the adaptive plasticity of telomeric complexes across eukaryotes.

Although repetitive DNA forms much of the human genome, its study is challenging due to limitations in assembly and alignment of repetitive short-reads. We have deployed k-Seek, software that detects tandem repeats embedded in single reads, on 2,504 human genomes from the 1,000 Genomes Project to quantify the variation and abundance of simple satellites (repeat units <20 bp). We find that the ancestral monomer of Human Satellite 3 makes up the largest portion of simple satellite content in humans (mean of ∼8 Mb). We discovered ∼50,000 rare tandem repeats that are not detected in the T2T-CHM13v2.0 assembly, including undescribed variants of telomericand pericentromeric repeats. We find broad homogeneity of the most abundant repeats across populations, except for AG-rich repeats which are more abundant in African individuals. We also find cliques of highly similar AG- and AT-rich satellites that are interspersed and form higher-order structures that covary in copy number across individuals, likely through concerted amplification via unequal exchange. Finally, we use pericentromeric polymorphisms to estimate centromeric genetic relatedness between individuals and find a strong predictive relationship between centromeric lineages and pericentromeric simple satellite abundances. In particular, ancestral monomers of Human Satellite 2 and Human Satellite 3 abundances correlate with clusters of centromeric ancestry on chromosome 16 and chromosome 9, with some clusters structured by population. These results provide new descriptions of the population dynamics that underlie the evolution of simple satellites in humans.

Meiotic drive biases the transmission of alleles in heterozygous individuals, such that Mendel's law of equal segregation is violated. Most examples of meiotic drive have been discovered over the past century based on causing sex ratio distortion or the biased transmission of easily scoreable genetic markers that were linked to drive alleles. More recently, several approaches have been developed that attempt to identify distortions of Mendelian segregation genome wide. Here, we test a candidate female meiotic drive locus in Drosophila melanogaster, identified previously as causing a ∼54:46 distortion ratio using sequencing of large pools of backcross progeny. We inserted fluorescent visible markers near the candidate locus and scored transmission in thousands of individual progeny. We observed a small but significant deviation from the Mendelian expectation; however, it was in the opposite direction to that predicted based on the original experiments. We discuss several possible causes of the discrepancy between the 2 approaches, noting that subtle viability effects are particularly challenging to disentangle from potential small-effect meiotic drive loci. We conclude that pool sequencing approaches remain a powerful method to identify candidate meiotic drive loci but that genotyping of individual progeny at early developmental stages may be required for robust confirmation.

Speciation is driven by traits that can act to prevent mating between nascent lineages, including male courtship and female preference for male traits. Mating barriers involving these traits evolve quickly because there is strong selection on males and females to maximize reproductive success, and the tight co-evolution of mating interactions can lead to rapid diversification of sexual behaviour. Populations of Drosophila melanogaster show strong asymmetrical reproductive isolation that is correlated with geographic origin. Using strains that capture natural variation in mating traits, we ask two key questions: which specific male traits are females selecting, and are these traits under divergent sexual selection? These questions have proven extremely challenging to answer, because even in closely related lineages males often differ in multiple traits related to mating behaviour. We address these questions by estimating selection gradients for male courtship and cuticular hydrocarbons for two different female genotypes. We identify specific behaviours and particular cuticular hydrocarbons that are under divergent sexual selection and could potentially contribute to premating reproductive isolation. Additionally, we report that a subset of these traits are plastic; males adjust these traits based on the identity of the female genotype they interact with. These results suggest that even when male courtship is not fixed between lineages, ongoing selection can act on traits that are important for reproductive isolation.

Crosses between Drosophila simulans females and Drosophila melanogaster males produce viable F1 sons and poorly viable F1 daughters. Unlike most hybrid incompatibilities, this hybrid incompatibility violates Haldane's rule, the observation that incompatibilities preferentially affect the heterogametic sex. Furthermore, it has a different genetic basis than hybrid lethality in the reciprocal cross, with the causal allele in Drosophila melanogaster being a large species-specific block of complex satellite DNA on its X chromosome known as the 359-bp satellite, rather than a protein-coding locus. The causal allele(s) in Drosophila simulans are unknown but likely involve maternally expressed genes or factors since the F1 females die during early embryogenesis. The maternal haploid (mh) gene is an intriguing candidate because it is expressed maternally and its protein product localizes to the 359-bp repeat. We found that this gene has diverged extensively between Drosophila melanogaster and Drosophila simulans. This observation led to the hypothesis that Drosophila melanogaster mh may have coevolved with the 359-bp repeat and that hybrid incompatibility thus results from the absence of a coevolved mh allele in Drosophila simulans. We tested for the functional divergence of mh by creating matched transformants of Drosophila melanogaster and Drosophila simulans orthologs in both Drosophila melanogaster and Drosophila simulans strains. Surprisingly, we find that Drosophila simulans mh fully complements the female sterile phenotype of Drosophila melanogaster mh mutations. Contrary to our hypothesis, we find no evidence that adding a Drosophila melanogaster mh gene to Drosophila simulans increases hybrid viability.

Germline stem cells (GSCs) are the progenitor cells of the germline for the lifetime of an animal. In Drosophila, these cells reside in a cellular niche that is required for both their maintenance (self-renewal) and differentiation (asymmetric division resulting in a daughter cell that differs from the GSC). The stem cell-daughter cell transition is tightly regulated by a number of processes, including an array of proteins required for genome stability. The germline stem-cell maintenance factor Stonewall (Stwl) associates with heterochromatin, but its molecular function is poorly understood. We performed RNA-Seq on stwl mutant ovaries and found significant derepression of many transposon families but not heterochromatic genes. We also discovered inappropriate expression of multiple classes of genes. Most prominent are testis-enriched genes, including the male germline sex-determination switch Phf7, the differentiation factor bgcn, and a large testis-specific gene cluster on chromosome 2, all of which are upregulated or ectopically expressed in stwl mutant ovaries. Surprisingly, we also found that RNAi knockdown of stwl in somatic S2 cells results in ectopic expression of these testis genes. Using parallel ChIP-Seq and RNA-Seq experiments in S2 cells, we discovered that Stwl localizes upstream of transcription start sites and at heterochromatic sequences including repetitive sequences associated with telomeres. Stwl is also enriched at bgcn, suggesting that it directly regulates this essential differentiation factor. Finally, we identify Stwl binding motifs that are shared with known insulator binding proteins. We propose that Stwl affects gene regulation, including repression of male transcripts in the female germline, by binding insulators and establishing chromatin boundaries.

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

Drosophila telomeres have been maintained by three families of active transposable elements (TEs), HeT-A, TAHRE, and TART, collectively referred to as HTTs, for tens of millions of years, which contrasts with an unusually high degree of HTT interspecific variation. While the impacts of conflict and domestication are often invoked to explain HTT variation, the telomeres are unstable structures such that neutral mutational processes and evolutionary tradeoffs may also drive HTT evolution. We leveraged population genomic data to analyze nearly 10,000 HTT insertions in 85 Drosophila melanogaster genomes and compared their variation to other more typical TE families. We observe that occasional large-scale copy number expansions of both HTTs and other TE families occur, highlighting that the HTTs are, like their feral cousins, typically repressed but primed to take over given the opportunity. However, large expansions of HTTs are not caused by the runaway activity of any particular HTT subfamilies or even associated with telomere-specific TE activity, as might be expected if HTTs are in strong genetic conflict with their hosts. Rather than conflict, we instead suggest that distinctive aspects of HTT copy number variation and sequence diversity largely reflect telomere instability, with HTT insertions being lost at much higher rates than other TEs elsewhere in the genome. We extend previous observations that telomere deletions occur at a high rate, and surprisingly discover that more than one-third do not appear to have been healed with an HTT insertion. We also report that some HTT families may be preferentially activated by the erosion of whole telomeres, implying the existence of HTT-specific host control mechanisms. We further suggest that the persistent telomere localization of HTTs may reflect a highly successful evolutionary strategy that trades away a stable insertion site in order to have reduced impact on the host genome. We propose that HTT evolution is driven by multiple processes, with niche specialization and telomere instability being previously underappreciated and likely predominant.

Germline stem cells (GSCs) are the progenitor cells of the germline for the lifetime of an animal. In Drosophila , these cells reside in a cellular niche that is required for both their maintenance (self-renewal) and differentiation (asymmetric division resulting in a daughter cell that differs from the GSC). The stem cell-daughter cell transition is tightly regulated by a number of processes, including an array of proteins required for genome stability. The germline stem-cell maintenance factor Stonewall (Stwl) associates with heterochromatin, but its molecular function is poorly understood. We performed RNA-Seq on stwl mutant ovaries and found significant derepression of many transposon families but not heterochromatic genes. We also discovered that testis-enriched genes, including the differentiation factor bgcn and a large testis-specific cluster on chromosome 2, are upregulated or ectopically expressed in stwl mutant ovaries. Surprisingly, we also found that RNAi knockdown of stwl in somatic S2 cells results in ectopic expression of these genes. Using parallel ChIP-Seq and RNA-Seq experiments in S2 cells, we discovered that Stwl binds upstream of transcription start sites and localizes to heterochromatic sequences. We also find that Stwl is enriched at repetitive sequences associated with telomeres. Finally, we identify Stwl binding motifs that are shared with known insulator binding proteins. We propose that Stwl affects gene regulation by binding insulators and establishing chromatin boundaries.

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

AbstractSpeciation is driven by traits that can act to prohibit mating between nascent lineages, including male courtship and female preference for male traits. Mating barriers involving these traits evolve quickly because there is strong selection on males and females to maximize reproductive success, and the tight co-evolution of mating interactions can lead to rapid diversification of sexual behavior. Using lineages ofD. melanogasterthat show strong asymmetrical reproductive isolation, we ask two key questions: which specific male traits are females selecting, and are these traits under divergent sexual selection? These questions have proven extremely challenging to answer, because even in closely related lineages males often differ in multiple traits related to mating behavior. We address these questions by estimating selection gradients for male courtship and cuticular hydrocarbons for two different female genotypes. We identify specific behaviors and particular cuticular hydrocarbons that are under divergent sexual selection and likely contribute to reproductive isolation. Additionally, we discovered that a subset of these traits are plastic; males adjust these traits based on the identity of the female genotype they interact with. These results suggest that even when male courtship is not fixed between lineages, ongoing selection can contribute to reproductive isolation.

leads to phenotypic divergence and speciation.

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.

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

Repetitive DNAs are ubiquitous in eukaryotic genomes and, in many species, comprise the bulk of the genome. Repeats include transposable elements that can self-mobilize and disperse around the genome and tandemly-repeated satellite DNAs that increase in copy number due to replication slippage and unequal crossing over. Despite their abundance, repetitive DNAs are often ignored in genomic studies due to technical challenges in identifying, assembling, and quantifying them. New technologies and methods are now allowing unprecedented power to analyze repetitive DNAs across diverse taxa. Repetitive DNAs are of particular interest because they can represent distinct modes of genome evolution. Some repetitive DNAs form essential genome structures, such as telomeres and centromeres, that are required for proper chromosome maintenance and segregation, while others form piRNA clusters that regulate transposable elements; thus, these elements are expected to evolve under purifying selection. In contrast, other repeats evolve selfishly and cause genetic conflicts with their host species that drive adaptive evolution of host defense systems. However, the majority of repeats likely accumulate in eukaryotes in the absence of selection due to mechanisms of transposition and unequal crossing over. However, even these "neutral" repeats may indirectly influence genome evolution as they reach high abundance. In this Special Issue, the contributing authors explore these questions from a range of perspectives.

In metazoan germlines, the piRNA pathway acts as a genomic immune system, employing small-RNA-mediated silencing to defend host DNA from the harmful effects of transposable elements (TEs). In response to dynamic changes in TE content, host genomes are proposed to alter the piRNAs that they produce in order to silence the most active TE families. Surprisingly, however, piRNA pathway proteins, which execute piRNA biogenesis and enforce silencing of targeted sequences, also evolve rapidly and adaptively in animals. If TE silencing evolves through changes in piRNAs, what necessitates changes in piRNA pathway proteins? Here we used interspecific complementation to test for functional differences between Drosophila melanogaster and D. simulans alleles of three adaptively evolving piRNA pathway proteins: Armitage, Aubergine and Spindle-E. Surprisingly, we found very few differences in TE transcript levels, but global effects on piRNA biogenesis, particularly for Armitage. Adaptive divergence in piRNA biogenesis in the absence of clear downstream effects on TE targets reveals a surprising degree of robustness that suggests that piRNA-mediated silencing is an “all-or-nothing” threshold trait. These results lead us to propose that the feature of importance for adaptation is piRNA production rather than piRNA abundance, and we further suggest that piRNA biogenesis evolves rapidly to establish control over invader and escaper TEs or in response to TE antagonism. We also discovered that D. simulans alleles exhibit enhanced off-target effects on host transcripts in a D. melanogaster background, suggesting the avoidance of genomic auto-immunity as an additional target of selection. Taken together, we propose that piRNA effector proteins are subject to an evolutionary trade-off between defending the host genome from the harmful effect of TEs while also minimizing friendly-fire against host genes.

A substantial portion of the genomes of most multicellular eukaryotes consists of large arrays of tandemly repeated sequence, collectively called satellite DNA. The processes generating and maintaining different satellite DNA abundances across lineages are important to understand as satellites have been linked to chromosome mis-segregation, disease phenotypes, and reproductive isolation between species. While much theory has been developed to describe satellite evolution, empirical tests of these models have fallen short because of the challenges in assessing satellite repeat regions of the genome. Advances in computational tools and sequencing technologies now enable identification and quantification of satellite sequences genome-wide. Here, we describe some of these tools and how their applications are furthering our knowledge of satellite evolution and function.

Simple satellites are tandemly repeating short DNA motifs that can span megabases in eukaryotic genomes. Because they can cause genomic instability through nonallelic homologous exchange, they are primarily found in the repressive heterochromatin near centromeres and telomeres where recombination is minimal, and on the Y chromosome, where they accumulate as the chromosome degenerates. Interestingly, the types and abundances of simple satellites often vary dramatically between closely related species, suggesting that they turn over rapidly. However, limited sampling has prevented detailed understanding of their evolutionary dynamics. Here, we characterize simple satellites from whole-genome sequences generated from males and females of nine Drosophila species, spanning 40 Ma of evolution. We show that PCR-free library preparation and postsequencing GC-correction better capture satellite quantities than conventional methods. We find that over half of the 207 simple satellites identified are species-specific, consistent with previous descriptions of their rapid evolution. Based on a maximum parsimony framework, we determined that most interspecific differences are due to lineage-specific gains. Simple satellites gained within a species are typically a single mutation away from abundant existing satellites, suggesting that they likely emerge from existing satellites, especially in the genomes of satellite-rich species. Interestingly, unlike most of the other lineages which experience various degrees of gains, the lineage leading up to the satellite-poor D. pseudoobscura and D. persimilis appears to be recalcitrant to gains, providing a counterpoint to the notion that simple satellites are universally rapidly evolving.

Genomic divergence can cause reproductive isolation between species. The molecular mechanisms underlying reproductive isolation can thus reveal which genomic features evolve rapidly and become unstable or incompatible in hybrids. In a recent paper in Nature, Gibeaux et al. (2018) report paternal genome instability and metabolic imbalance in hybrids between frog species.

Eukaryotic genomes are replete with repeated sequences in the form of transposable elements (TEs) dispersed across the genome or as satellite arrays, large stretches of tandemly repeated sequences. Many satellites clearly originated as TEs, but it is unclear how mobile genetic parasites can transform into megabase-sized tandem arrays. Comprehensive population genomic sampling is needed to determine the frequency and generative mechanisms of tandem TEs, at all stages from their initial formation to their subsequent expansion and maintenance as satellites. The best available population resources, short-read DNA sequences, are often considered to be of limited utility for analyzing repetitive DNA due to the challenge of mapping individual repeats to unique genomic locations. Here we develop a new pipeline called ConTExt that demonstrates that paired-end Illumina data can be successfully leveraged to identify a wide range of structural variation within repetitive sequence, including tandem elements. By analyzing 85 genomes from five populations of Drosophila melanogaster, we discover that TEs commonly form tandem dimers. Our results further suggest that insertion site preference is the major mechanism by which dimers arise and that, consequently, dimers form rapidly during periods of active transposition. This abundance of TE dimers has the potential to provide source material for future expansion into satellite arrays, and we discover one such copy number expansion of the DNA transposon hobo to approximately 16 tandem copies in a single line. The very process that defines TEs-transposition-thus regularly generates sequences from which new satellites can arise.