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.

Meiosis is a potentially important source of germline mutations, as sites of meiotic recombination experience recurrent double-strand breaks (DSBs). However, evidence for a local mutagenic effect of recombination from population sequence data has been equivocal, likely because mutation is only one of several forces shaping sequence variation. By sequencing large numbers of single crossover molecules obtained from human sperm for two recombination hotspots, we find direct evidence that recombination is mutagenic: Crossovers carry more de novo mutations than nonrecombinant DNA molecules analyzed for the same donors and hotspots. The observed mutations were primarily CG to TA transitions, with a higher frequency of transitions at CpG than non-CpGs sites. This enrichment of mutations at CpG sites at hotspots could predominate in methylated regions involving frequent single-stranded DNA processing as part of DSB repair. In addition, our data set provides evidence that GC alleles are preferentially transmitted during crossing over, opposing mutation, and shows that GC-biased gene conversion (gBGC) predominates over mutation in the sequence evolution of hotspots. These findings are consistent with the idea that gBGC could be an adaptation to counteract the mutational load of recombination. meiotic recombination | crossover | sequence evolution | biased gene conversion | mutation M eiotic recombination, localized in recombination hotspots, not only increases genetic diversity via the formation of new haplotypes but is also an important driver of sequence evolution. The binding sites used by the human recombination machinery involving PRDM9 (PR domain containing 9) are more eroded in humans than the same sequences in chimps, given that PRDM9 in chimps uses different binding sites (1). Moreover, regions in close vicinity to these PRDM9 binding sites also showed a significant enrichment of polymorphisms in humans (2). In addition, within-and between-species sequence diversity positively correlates with regions of high recombination activity in humans (3-7) and other eukaryotes (reviewed in refs. 8-10).One process recognized as a major evolutionary force reshaping the genomic nucleotide landscape at recombination hotspots, as shown in humans (6), chimpanzees (6, 11), mice (12), yeast (13), and metazoans (14), is GC-biased gene conversion (gBGC). In gBGC, the repair of heteroduplex tracts formed during meiotic recombination leads to the nonMendelian segregation of alleles favoring GC over AT variants. The precise molecular mechanisms leading to gBGC have yet to be unraveled, but experimental evidence has shown that in crossovers (COs) of fission yeast, GC alleles can be overtransmitted within ∼1-2 kb in length of the double-strand break (DSB) region (13), implicating mismatch repair (15).However, it is also plausible that the higher sequence variation observed at recombination hotspots is a result of a mutagenic effect of recombination: meiotic recombination is initiated by DSBs, which are associated with an increased mutati...

Abstract Genetic variation is the fuel of evolution, with standing genetic variation especially important for short-term evolution and local adaptation. To date, studies of spatiotemporal patterns of genetic variation in natural populations have been challenging, as comprehensive sampling is logistically difficult, and sequencing of entire populations costly. Here, we address these issues using a collaborative approach, sequencing 48 pooled population samples from 32 locations, and perform the first continent-wide genomic analysis of genetic variation in European Drosophila melanogaster. Our analyses uncover longitudinal population structure, provide evidence for continent-wide selective sweeps, identify candidate genes for local climate adaptation, and document clines in chromosomal inversion and transposable element frequencies. We also characterize variation among populations in the composition of the fly microbiome, and identify five new DNA viruses in our samples.

Various approaches can be applied to uncover the genetic basis of natural phenotypic variation, each with their specific strengths and limitations. Here, we use a replicated genome-wide association approach (Pool-GWAS) to fine-scale map genomic regions contributing to natural variation in female abdominal pigmentation in Drosophila melanogaster, a trait that is highly variable in natural populations and highly heritable in the laboratory. We examined abdominal pigmentation phenotypes in approximately 8000 female European D. melanogaster, isolating 1000 individuals with extreme phenotypes. We then used whole-genome Illumina sequencing to identify single nucleotide polymorphisms (SNPs) segregating in our sample, and tested these for associations with pigmentation by contrasting allele frequencies between replicate pools of light and dark individuals. We identify two small regions near the pigmentation genes tan and bric-à-brac 1, both corresponding to known cis-regulatory regions, which contain SNPs showing significant associations with pigmentation variation. While the Pool-GWAS approach suffers some limitations, its cost advantage facilitates replication and it can be applied to any non-model system with an available reference genome.

Population genetic theory shows that the efficacy of natural selection is limited by linkage-selection at one site interferes with selection at linked sites. Such interference slows adaptation in asexual genomes and may explain the evolutionary advantage of sex. Here, we test for two signatures of constraint caused by linkage in a sexual genome, by using sequence data from 255 Drosophila melanogaster and Drosophila simulans loci. We find that (i) the rate of protein adaptation is reduced in regions of low recombination, and (ii) evolution at strongly selected amino acid sites interferes with optimal codon usage at weakly selected, tightly linked synonymous sites. Together these findings suggest that linkage limits the rate and degree of adaptation even in recombining genomes. N atural selection is imperfect. To become fixed, beneficial mutations must overcome both stochastic loss and interference from selection at linked loci. In asexual genomes, where linkage is complete, two kinds of interference compromise adaptation. The first, ''ruby-in-the-rubbish'' interference, occurs because beneficial mutations often appear on genetic backgrounds loaded with segregating deleterious mutations. Since deleterious mutations are, on average, probably more strongly selected than favorable ones, adaptation is mostly limited to those lucky few beneficial mutations that arise on unloaded backgrounds (1-4). The second form of interference, ''clonal'' interference, is caused by competition among multiple segregating beneficial mutations (1, 5-7). Because only one asexual genome can be fixed at a time, adaptive substitutions are forced to be nearly sequential. Both kinds of interference limit the rate of adaptation in asexuals (8-10).The effects of both kinds of interference can be thought of as a reduction in the effective population size (N e ) caused by selection at linked loci (11-13). Recombination, by alleviating interference between linked sites, alleviates this reduction in N e . Consequently, genomic regions that differ in recombination rate also differ in effective size-indeed, this is the basis of the well-known correlation between recombination rate and levels of neutral polymorphism (4N e times the neutral mutation rate) seen in the genomes of Drosophila, humans, and others (14-18). That variation in linkage affects levels of neutral polymorphism suggests that it may also affect rates of nonneutral substitution. In particular, adaptive evolution may be limited in regions of low recombination (i.e., where N e is reduced) or in situations of extreme linkage (e.g., among sites within the same gene).Here we ask whether linkage systematically constrains adaptation in the Drosophila genome. We use divergence estimates from 255 Drosophila melanogaster and Drosophila simulans loci. These data are unique in that they include a large number of rapidly evolving genes, many of which are candidate male accessory gland proteins (Acps) and thus likely targets of sexual selection (19-21). We have reason to believe a priori not only t...