Stress response, behavior, and development are shaped by transposable element-induced mutations in Drosophila

Mapping genotype to phenotype is challenging because of the difficulties in identifying both the traits under selection and the specific genetic variants underlying these traits. Most of the current knowledge of the genetic basis of adaptive evolution is based on the analysis of single nucleotide polymorphisms (SNPs). Despite increasing evidence for their causal role, the contribution of structural variants to adaptive evolution remains largely unexplored. In this work, we analyzed the population frequencies of 1,615 Transposable Element (TE) insertions in 91 samples from 60 worldwide natural populations of Drosophila melanogaster . We identified a set of 300 TEs that are present at high population frequencies, and located in genomic regions with high recombination rate, where the efficiency of natural selection is high. The age and the length of these 300 TEs are consistent with relatively young and long insertions reaching high frequencies due to the action of positive selection. Indeed, we, and others, found evidence of selective sweeps and/or population differentiation for 65 of them. The analysis of the genes located nearby these 65 candidate adaptive insertions suggested that the functional response to selection is related with the GO categories of response to stimulus, behavior, and development. We further showed that a subset of the candidate adaptive TEs affect expression of nearby genes, and five of them have already been linked to an ecologically relevant phenotypic effect. Our results provide a more complete understanding of the genetic variation and the fitness-related traits relevant for adaptive evolution. Similar studies should help uncover the importance of TE-induced adaptive mutations in other species as well.