The lack of recombination drives the fixation of transposable elements on the fourth chromosome of Drosophila melanogaster

In regions of suppressed recombination, where selection is expected to be less efficient in removing slightly deleterious mutations, transposable element (TE) insertions should be more likely to drift to higher frequencies, and even to reach fixation. In the absence of excision events, once a TE is fixed it cannot be eliminated from the population, and accumulation of elements thus should become an irreversible process. In the long term, this can drive the degeneration of large non-recombining fractions of the genomes. Chromosome 4 of Drosophila melanogaster has very low levels of recombination, if any, and this could be causing its degeneration. Here we report the results of a PCR-based analysis of the population frequencies of TE insertions in a sample from three African natural populations. We investigated 27 insertions from 12 TE families, located in regions of either suppressed or free recombination. Our results suggest that TE insertions tend to be fixed in the non-recombining regions, particularly on the fourth chromosome. We have also found that this involves all types of elements, and that fixed insertions are significantly shorter and more divergent from the canonical sequence than those segregating in the sample (28.1% vs 86.3% of the canonical length, and average nucleotide divergence (D(XY)) = 0.082 vs 0.008, respectively). Finally, DNA-based elements seem to show a greater tendency to reach fixation than retrotransposons. Implications of these findings for the population dynamics of TEs, and the evolutionary forces that shape the patterns of genetic variation in regions of reduced recombination, are discussed.     

Population‐specific dynamics and selection patterns of transposable element insertions in European natural populations

Transposable elements (TEs) are ubiquitous sequences in genomes of virtually all species. While TEs have been investigated for several decades, only recently we have the opportunity to study their genome‐wide population dynamics. Most of the studies so far have been restricted either to the analysis of the insertions annotated in the reference genome or to the analysis of a limited number of populations. Taking advantage of the European Drosophila population genomics consortium (DrosEU) sequencing data set, we have identified and measured the dynamics of TEs in a large sample of European Drosophila melanogaster natural populations. We showed that the mobilome landscape is population‐specific and highly diverse depending on the TE family. In contrast with previous studies based on SNP variants, no geographical structure was observed for TE abundance or TE divergence in European populations. We further identified de novo individual insertions using two available programs and, as expected, most of the insertions were present at low frequencies. Nevertheless, we identified a subset of TEs present at high frequencies and located in genomic regions with a high recombination rate. These TEs are candidates for being the target of positive selection, although neutral processes should be discarded before reaching any conclusion on the type of selection acting on them. Finally, parallel patterns of association between the frequency of TE insertions and several geographical and temporal variables were found between European and North American populations, suggesting that TEs can be potentially implicated in the adaptation of populations across continents.     

Size matters: non-LTR retrotransposable elements and ectopic recombination in Drosophila

The Drosophila melanogaster genome contains approximately 100 distinct families of transposable elements (TEs). In the euchromatic part of the genome, each family is present in a small number of copies (5-150 copies), with individual copies of TEs often present at very low frequencies in populations. This pattern is likely to reflect a balance between the inflow of TEs by transposition and the removal of TEs by natural selection. The nature of natural selection acting against TEs remains controversial. We provide evidence that selection against chromosome abnormalities caused by ectopic recombination limits the spread of some TEs. We also demonstrate for the first time that some TE families in the Drosophila euchromatin appear to be only marginally affected by purifying selection and contain many copies at high population frequencies. We argue that TEs in these families attain high population frequencies and even reach fixation as a result of low family-wide transposition rates leading to low TE copy numbers and consequently reduced strength of selection acting on individual TE copies. Fixation of TEs in these families should provide an upward pressure on the size of intergenic sequences counterbalancing rapid DNA loss through small deletions. Copy-number-dependent selection on TE families caused by ectopic recombination may also promote diversity among TEs in the Drosophila genome.     

Sequencing of pooled DNA samples (Pool-Seq) uncovers complex dynamics of transposable element insertions in Drosophila melanogaster

Transposable elements (TEs) are mobile genetic elements that parasitize genomes by semi-autonomously increasing their own copy number within the host genome. While TEs are important for genome evolution, appropriate methods for performing unbiased genome-wide surveys of TE variation in natural populations have been lacking. Here, we describe a novel and cost-effective approach for estimating population frequencies of TE insertions using paired-end Illumina reads from a pooled population sample. Importantly, the method treats insertions present in and absent from the reference genome identically, allowing unbiased TE population frequency estimates. We apply this method to data from a natural Drosophila melanogaster population from Portugal. Consistent with previous reports, we show that low recombining genomic regions harbor more TE insertions and maintain insertions at higher frequencies than do high recombining regions. We conservatively estimate that there are almost twice as many "novel" TE insertion sites as sites known from the reference sequence in our population sample (6,824 novel versus 3,639 reference sites, with on average a 31-fold coverage per insertion site). Different families of transposable elements show large differences in their insertion densities and population frequencies. Our analyses suggest that the history of TE activity significantly contributes to this pattern, with recently active families segregating at lower frequencies than those active in the more distant past. Finally, using our high-resolution TE abundance measurements, we identified 13 candidate positively selected TE insertions based on their high population frequencies and on low Tajima's D values in their neighborhoods.     

Remarkable site specificity of local transposition into the Hsp70 promoter of Drosophila melanogaster

Heat-shock genes have numerous features that ought to predispose them to insertional mutagenesis via transposition. To elucidate the evolvability of heat-shock genes via transposition, we have exploited a local transposition technique and Drosophila melanogaster strains with EPgy2 insertions near the Hsp70 gene cluster at 87A7 to produce numerous novel EPgy2 insertions into these Hsp70 genes. More than 50% of 45 independent insertions were made into two adjacent nucleotides in the proximal promoter at positions -96 and -97, and no insertions were into a coding or 3'-flanking sequence. All inserted transposons were in inverse orientation to the starting transposon. The frequent insertion into nucleotides -96 and -97 is consistent with the DNase hypersensitivity, absence of nucleosomes, flanking GAGA-factor-binding sites, and nucleotide sequence of this region. These experimental insertions recapitulated many of the phenotypes of natural transposition into Hsp70: reduced mRNA expression, less Hsp70 protein, and decreased inducible thermotolerance. The results suggest that the distinctive features of heat-shock promoters, which underlie the massive and rapid expression of heat-shock genes upon heat shock, also are a source of evolutionary variation on which natural selection can act.     

Response to natural and laboratory selection at the Drosophila hsp70 genes

Abstract.— To determine whether and how laboratory and natural selection act on the hsp70 (70-Kd heat-shock protein) genes of Drosophila melanogaster , we examined hsp70 allele frequencies in two sets of populations. First, five populations reared at different temperatures for more than 20 years differentially fixed both a large insertion/deletion (indel) polymorphism at the 87A7 hsp70 cluster ("56H8"/"122") and a single nucleotide polymorphism at the 87C1 hsp70 cluster. In both cases, the 18°C and 25°C populations fixed one allele and the 28°C populations the other, consistent with previously described evolved differences among these populations in Hsp70 expression and thermo-tolerance. Second, we examined 56H8 and 122 frequencies in a set of 11 populations founded from flies collected along a latitudinal transect of eastern Australia. The 56H8 allele frequencies are positively associated with latitude, consistent with maintenance of the 56H8/122 polymorphism by natural selection. Thermal extremes and average values are negatively correlated with latitude. These results suggest that natural selection imposed by temperature and thermal variability may affect hsp70 allele frequencies.     

Genome-wide variation in the human and fruitfly: a comparison

Average levels of nucleotide diversity are ten-fold lower in humans than in the fruitfly, Drosophila melanogaster. Despite this difference, apparently as a result of a lower population size, patterns of genomic diversity are strikingly similar in being correlated with local rates of recombination, and influenced by similar interactions between positive natural selection and recombination. Both species also show lower levels of variation on average in non-African compared to African populations, reflecting a similar evolutionary history and perhaps both natural selection and founder effects in new environments.     

Transposable Element Numbers in Cosmopolitan Inversions from a Natural Population of Drosophila Melanogaster

Population studies of the distribution of transposable elements (TEs) on the chromosomes of Drosophila melanogaster have suggested that their copy number increase due to transposition is balanced by some form of natural selection. Theory suggests that, as a consequence of deleterious ectopic meiotic exchange between TEs, selection can favor genomes with lower TE copy numbers. This predicts that TEs should be less deleterious, and hence more abundant, in chromosomal regions in which recombination is reduced. To test this, we surveyed the abundance and locations of 10 families of TEs in recombination-suppressing chromosomal inversions from a natural population. The sample of 49 chromosomes included multiple independent isolates of seven different inversions and a corresponding set of standard chromosomes. For all 10 TE families pooled, copy numbers were significantly higher overall within low frequency inversions than within corresponding regions of standard chromosomes. TEs occupied chromosomal sites at significantly higher frequencies within the In(3R)M0 and In(3R)K inversions than within the corresponding regions of standard 3R chromosomes. These results are consistent with the predictions of the ectopic exchange model.     

Population genomics of transposable elements in Drosophila melanogaster

Transposable elements (TEs) are the primary contributors to the genome bulk in many organisms and are major players in genome evolution. A clear and thorough understanding of the population dynamics of TEs is therefore essential for full comprehension of the eukaryotic genome evolution and function. Although TEs in Drosophila melanogaster have received much attention, population dynamics of most TE families in this species remains entirely unexplored. It is not clear whether the same population processes can account for the population behaviors of all TEs in Drosophila or whether, as has been suggested previously, different orders behave according to very different rules. In this work, we analyzed population frequencies for a large number of individual TEs (755 TEs) in five North American and one sub-Saharan African D. melanogaster populations (75 strains in total). These TEs have been annotated in the reference D. melanogaster euchromatic genome and have been sampled from all three major orders (non-LTR, LTR, and TIR) and from all families with more than 20 TE copies (55 families in total). We find strong evidence that TEs in Drosophila across all orders and families are subject to purifying selection at the level of ectopic recombination. We showed that strength of this selection varies predictably with recombination rate, length of individual TEs, and copy number and length of other TEs in the same family. Importantly, these rules do not appear to vary across orders. Finally, we built a statistical model that considered only individual TE-level (such as the TE length) and family-level properties (such as the copy number) and were able to explain more than 40% of the variation in TE frequencies in D. melanogaster.     

On the abundance and distribution of transposable elements in the genome of Drosophila melanogaster

The abundance and distribution of transposable elements (TEs) in a representative part of the euchromatic genome of Drosophila melanogaster were studied by analyzing the sizes and locations of TEs of all known families in the genomic sequences of chromosomes 2R, X, and 4. TEs contribute to up to 2% of the sequenced DNA, which corresponds roughly to the euchromatin of these chromosomes. This estimate is lower than that previously available from in situ data and suggests that TEs accumulate in the heterochromatin more intensively than was previously thought. We have also found that TEs are not distributed at random in the chromosomes and that their abundance is more strongly associated with local recombination rates, rather than with gene density. The results are compatible with the ectopic exchange model, which proposes that selection against deleterious effects of chromosomal rearrangements is a major force opposing element spread in the genome of this species. Selection against insertional mutations also influences the observed patterns, such as an absence of insertions in coding regions. The results of the analyses are discussed in the light of recent findings on the distribution of TEs in other species.     

The effects of recombination rate on the distribution and abundance of transposable elements

Transposable elements (TEs) often accumulate in regions of the genome with suppressed recombination. But it is unclear whether this pattern reflects a reduction in the efficacy of selection against deleterious insertions or a relaxation of ectopic recombination. Discriminating between these two hypotheses has been difficult, because no formal model has investigated the effects of recombination under the deleterious insertion model. Here we take a simulation-based approach to analyze this scenario and determine the conditions under which element accumulation is expected in low recombination regions. We show that TEs become fixed as a result of Hill-Robertson effects in the form of Muller's ratchet, but only in regions of extremely low recombination when excision is effectively absent and synergism between elements is weak. These results have important implications for differentiating between the leading models of how selection acts on TEs and should help to interpret emerging population genetic and genomic data.     

Impact of polymorphic transposable elements on linkage disequilibrium along chromosomes

Recombination and selection drive the extent of linkage disequilibrium (LD) among loci and therefore affect the reshuffling of adaptive genetic variation. However, it is poorly known to what extent the enrichment of transposable elements (TEs) in recombinationally‐inert regions reflects their inefficient removal by purifying selection and whether the presence of polymorphic TEs can modify the local recombination rate. In this study, we investigate how TEs and recombination interact at fine scale along chromosomes and possibly support linked selection in natural populations. Whole‐genome sequencing data of 304 individuals from nearby alpine populations of Arabis alpina were used to show that the density of polymorphic TEs is specifically correlated with local LD along chromosomes. Consistent with TEs modifying recombination, the characterization of 28 such LD blocks of up to 5.5 Mb in length revealed strong evidence of selective sweeps at a few loci through either site frequency spectrum or haplotype structure. A majority of these blocks were enriched in genes related to ecologically relevant functions such as responses to cold, salt stress or photoperiodism. In particular, the S‐locus (i.e., supergene responsible for strict outcrossing) was identified in a LD block with high levels of polymorphic TEs and evidence of selection. Another such LD block was enriched in cold‐responding genes and presented evidence of adaptive loci related to photoperiodism and flowering being increasingly linked by polymorphic TEs. These results are consistent with the hypothesis that TEs modify recombination landscapes and thus interact with selection in driving blocks of linked adaptive loci in natural populations.     

Active miniature transposons from a plant genome and its nonrecombining Y chromosome

Mechanisms involved in eroding fitness of evolving Y chromosomes have been the focus of much theoretical and empirical work. Evolving Y chromosomes are expected to accumulate transposable elements (TEs), but it is not known whether such accumulation contributes to their genetic degeneration. Among TEs, miniature inverted-repeat transposable elements are nonautonomous DNA transposons, often inserted in introns and untranslated regions of genes. Thus, if they invade Y-linked genes and selection against their insertion is ineffective, they could contribute to genetic degeneration of evolving Y chromosomes. Here, we examine the population dynamics of active MITEs in the young Y chromosomes of the plant Silene latifolia and compare their distribution with those in recombining genomic regions. To isolate active MITEs, we developed a straightforward approach on the basis of the assumption that recent transposon insertions or excisions create singleton or low-frequency size polymorphisms that can be detected in alleles from natural populations. Transposon display was then used to infer the distribution of MITE insertion frequencies. The overall frequency spectrum showed an excess of singleton and low-frequency insertions, which suggests that these elements are readily removed from recombining chromosomes. In contrast, insertions on the Y chromosomes were present at high frequencies. Their potential contribution to Y degeneration is discussed.     

Genome organization and gene expression shape the transposable element distribution in the Drosophila melanogaster euchromatin

The distribution of transposable elements (TEs) in a genome reflects a balance between insertion rate and selection against new insertions. Understanding the distribution of TEs therefore provides insights into the forces shaping the organization of genomes. Past research has shown that TEs tend to accumulate in genomic regions with low gene density and low recombination rate. However, little is known about the factors modulating insertion rates across the genome and their evolutionary significance. One candidate factor is gene expression, which has been suggested to increase local insertion rate by rendering DNA more accessible. We test this hypothesis by comparing the TE density around germline- and soma-expressed genes in the euchromatin of Drosophila melanogaster. Because only insertions that occur in the germline are transmitted to the next generation, we predicted a higher density of TEs around germline-expressed genes than soma-expressed genes. We show that the rate of TE insertions is greater near germline- than soma-expressed genes. However, this effect is partly offset by stronger selection for genome compactness (against excess noncoding DNA) on germline-expressed genes. We also demonstrate that the local genome organization in clusters of coexpressed genes plays a fundamental role in the genomic distribution of TEs. Our analysis shows that—in addition to recombination rate—the distribution of TEs is shaped by the interaction of gene expression and genome organization. The important role of selection for compactness sheds a new light on the role of TEs in genome evolution. Instead of making genomes grow passively, TEs are controlled by the forces shaping genome compactness, most likely linked to the efficiency of gene expression or its complexity and possibly their interaction with mechanisms of TE silencing.     

Naturally occurring transposable elements disrupt hsp70 promoter function in Drosophila melanogaster

Naturally occurring transposable element (TE) insertions that disrupt Drosophila promoters are correlated with modified promoter function and are posited to play a significant role in regulatory evolution, but their phenotypes have not been established directly. To establish the functional consequences of these TE insertions, we created constructs with either TE-bearing or TE-lacking hsp70 promoters fused to a luciferase reporter gene and assayed luciferase luminescence in transiently transfected Drosophila cells. Each of the four TEs reduces luciferase signal after heat shock and heat inducibility of the hsp70 promoter. To test if the differences in hsp70 promoter activity are TE-sequence dependent, we replaced each of the TEs with multiple intergenic sequences of equal length. These replacement insertions similarly reduced luciferase signal, suggesting that the TEs affect hsp70 promoter function by altering promoter architecture. These results are consistent with differences in Hsp70 expression levels, inducible thermotolerance, and fecundity previously associated with the TEs. That two different varieties of TEs in two different hsp70 genes have common effects suggests that TE insertion represents a general mechanism through which selection manipulates hsp70 gene expression.     

Patterns of insertion and deletion in contrasting chromatin domains

Transposable elements (TEs) play a fundamental role in the evolution of genomes. In Drosophila they are disproportionately represented in regions of low recombination, such as in heterochromatin. This pattern has been attributed to selection against repeated elements in regions of normal recombination, owing to either (1) the slightly deleterious position effects of TE insertions near or into genes, or (2) strong selection against chromosomal abnormalities arising from ectopic exchange between TE repeats. We have used defective non-long-terminal repeat (LTR) TEs that are "dead-on-arrival" (DOA) and unable to transpose in order to estimate spontaneous deletion rates in different constituents of chromatin. These elements have previously provided evidence for an extremely high rate of spontaneous deletion in Drosophila as compared with mammals, potentially explaining at least part of the differences in the genome sizes in these organisms. However, rates of deletion could be overestimated due to positive selection for a smaller likelihood of ectopic exchange. In this article, we show that rates of spontaneous deletion in DOA repeats are as high in heterochromatin and regions of euchromatin with low recombination as they are in regions of euchromatin with normal recombination. We have also examined the age distribution of five non-LTR families throughout the genome. We show that there is substantial variation in the historical pattern of transposition of these TEs. The overrepresentation of TEs in the heterochromatin is primarily due to their longer retention time in heterochromatin, as evidenced by the average time since insertion. Fragments inserted recently are much more evenly distributed in the genome. This contrast demonstrates that the accumulation of TEs in heterochromatin and in euchromatic regions of low recombination is not due to biased transposition but by greater probabilities of fixation in these regions relative to regions of normal recombination.     

Restriction-map variation with the yellow-achaete-scute region in five populations of Drosophila melanogaster

It has been proposed that the degree of recombination for a genomic region will affect the level of both nucleotide heterozygosity and the density of transposable elements. Both features of genomic diversity have been examined in a number of recent reports for regions undergoing relatively normal levels of recombination in Drosophila melanogaster. In this study the genomic variation associated with yellow-achaete- scute loci located at the tip of the X chromosome is examined by six- cutter restriction mapping. In this region, as usual for regions adjacent to telomeres, crossing-over is dramatically reduced, and published studies of visible mutants indicate extremely little restriction-map variation. Eight six-cutter restriction endonucleases were used to locate sequence variation in 14- and 16.5-kb regions in 109 lines sampled from North America, Africa, and Europe. The overall level of heterozygosity is estimated as 0.29%. Nine large insertions, all presumed to be transposable elements, were observed. Base-pair heterozygosity appears to be reduced compared with regions having normal levels of recombination. The estimated heterozygosity is much higher than reported in earlier studies of restriction-map variation among visible mutations in the complex. The incidence of large insertions is not elevated compared with that in other regions of the genome. This suggests that asymmetric synapsis and exchange is not an important mechanism for the elimination of transposable elements.      

Genomic analysis of adaptive differentiation in Drosophila melanogaster

Drosophila melanogaster shows clinal variation along latitudinal transects on multiple continents for several phenotypes, allozyme variants, sequence variants, and chromosome inversions. Previous investigation suggests that many such clines are due to spatially varying selection rather than demographic history, but the genomic extent of such selection is unknown. To map differentiation throughout the genome, we hybridized DNA from temperate and subtropical populations to Affymetrix tiling arrays. The dense genomic sampling of variants and low level of linkage disequilibrium in D. melanogaster enabled identification of many small, differentiated regions. Many regions are differentiated in parallel in the United States and Australia, strongly supporting the idea that they are influenced by spatially varying selection. Genomic differentiation is distributed nonrandomly with respect to gene function, even in regions differentiated on only one continent, providing further evidence for the role of selection. These data provide candidate genes for phenotypes known to vary clinally and implicate interesting new processes in genotype-by-environment interactions, including chorion proteins, proteins regulating meiotic recombination and segregation, gustatory and olfactory receptors, and proteins affecting synaptic function and behavior. This portrait of differentiation provides a genomic perspective on adaptation and the maintenance of variation through spatially varying selection.