High rate of recent transposable element-induced adaptation in Drosophila melanogaster

Evolutionary genomics of transposable elements in Saccharomyces cerevisiae

Saccharomyces cerevisiae is one of the premier model systems for studying the genomics and evolution of transposable elements. The availability of the S. cerevisiae genome led to unprecedented insights into its five known transposable element families (the LTR retrotransposons Ty1-Ty5) in the years shortly after its completion. However, subsequent advances in bioinformatics tools for analysing transposable elements and the recent availability of genome sequences for multiple strains and species of yeast motivates new investigations into Ty evolution in S. cerevisiae. Here we provide a comprehensive phylogenetic and population genetic analysis of all Ty families in S. cerevisiae based on a systematic re-annotation of Ty elements in the S288c reference genome. We show that previous annotation efforts have underestimated the total copy number of Ty elements for all known families. In addition, we identify a new family of Ty3-like elements related to the S. paradoxus Ty3p which is composed entirely of degenerate solo LTRs. Phylogenetic analyses of LTR sequences identified three families with short-branch, recently active clades nested among long branch, inactive insertions (Ty1, Ty3, Ty4), one family with essentially all recently active elements (Ty2) and two families with only inactive elements (Ty3p and Ty5). Population genomic data from 38 additional strains of S. cerevisiae show that the majority of Ty insertions in the S288c reference genome are fixed in the species, with insertions in active clades being predominantly polymorphic and insertions in inactive clades being predominantly fixed. Finally, we use comparative genomic data to provide evidence that the Ty2 and Ty3p families have arisen in the S. cerevisiae genome by horizontal transfer. Our results demonstrate that the genome of a single individual contains important information about the state of TE population dynamics within a species and suggest that horizontal transfer may play an important role in shaping the genomic diversity of transposable elements in unicellular eukaryotes.

Genetic variation within native populations of endemic silkmoth Antheraea assamensis (Helfer) from Northeast India indicates need for in situ conservation

A. assamensis is a phytophagous Lepidoptera from Northeast India reared on host trees of Lauraceae family for its characteristic cocoon silk. Source of these cocoons are domesticated farm stocks that crash frequently and/or wild insect populations that provide new cultures. The need to reduce dependence on wild populations for cocoons necessitates assessment of genetic diversity in cultivated and wild populations. Molecular markers based on PCR of Inter-simple sequence repeats (ISSR) and simple sequence repeats (SSR) were used with four populations of wild insects and eleven populations of cultivated insects. Wild populations had high genetic diversity estimates (H_i = 0.25; H_S = 0.28; H_E = 0.42) and at least one population contained private alleles. Both marker systems indicated that genetic variability within populations examined was significantly high. Among cultivated populations, insects of the Upper Assam region (H_i = 0.19; H_S = 0.18; H_E = 0) were genetically distinct (F_ST = 0.38 with both marker systems) from insects of Lower Assam (H_i = 0.24; H_S = 0.25; H_E = 0.3). Sequencing of polymorphic amplicons suggested transposition as a mechanism for maintaining genomic diversity. Implications for conservation of native populations in the wild and preserving in-farm diversity are discussed.

Population and sex differences in Drosophila melanogaster brain gene expression

Background

Changes in gene regulation are thought to be crucial for the adaptation of organisms to their environment. Transcriptome analyses can be used to identify candidate genes for ecological adaptation, but can be complicated by variation in gene expression between tissues, sexes, or individuals. Here we use high-throughput RNA sequencing of a single Drosophila melanogaster tissue to detect brain-specific differences in gene expression between the sexes and between two populations, one from the ancestral species range in sub-Saharan Africa and one from the recently colonized species range in Europe.

Results

Relatively few genes (<100) displayed sexually dimorphic expression in the brain, but there was an enrichment of sex-biased genes, especially male-biased genes, on the X chromosome. Over 340 genes differed in brain expression between flies from the African and European populations, with the inter-population divergence being highly correlated between males and females. The differentially expressed genes included those involved in stress response, olfaction, and detoxification. Expression differences were associated with transposable element insertions at two genes implicated in insecticide resistance (Cyp6g1 and CHKov1).

Conclusions

Analysis of the brain transcriptome revealed many genes differing in expression between populations that were not detected in previous studies using whole flies. There was little evidence for sex-specific regulatory adaptation in the brain, as most expression differences between populations were observed in both males and females. The enrichment of genes with sexually dimorphic expression on the X chromosome is consistent with dosage compensation mechanisms affecting sex-biased expression in somatic tissues.

Drosophila interspecific hybrids phenocopy piRNA-pathway mutants

Hybrids of two Drosophila species show transposable element derepression and piRNA pathway malfunction, revealing adaptive evolution of piRNA pathway components.

The Piwi-interacting RNA (piRNA) pathway defends the germline of animals from the deleterious activity of selfish transposable elements (TEs) through small-RNA mediated silencing. Adaptation to novel invasive TEs is proposed to occur by incorporating their sequences into the piRNA pool that females produce and deposit into their eggs, which then propagates immunity against specific TEs to future generations. In support of this model, the F1 offspring of crosses between strains of the same Drosophila species sometimes suffer from germline derepression of paternally inherited TE families, caused by a failure of the maternal strain to produce the piRNAs necessary for their regulation. However, many protein components of the Drosophila piRNA pathway exhibit signatures of positive selection, suggesting that they also contribute to the evolution of host genome defense. Here we investigate piRNA pathway function and TE regulation in the F1 hybrids of interspecific crosses between D. melanogaster and D. simulans and compare them with intraspecific control crosses of D. melanogaster. We confirm previous reports showing that intraspecific crosses are characterized by derepression of paternally inherited TE families that are rare or absent from the maternal genome and piRNA pool, consistent with the role of maternally deposited piRNAs in shaping TE silencing. In contrast to the intraspecific cross, we discover that interspecific hybrids are characterized by widespread derepression of both maternally and paternally inherited TE families. Furthermore, the pattern of derepression of TE families in interspecific hybrids cannot be attributed to their paucity or absence from the piRNA pool of the maternal species. Rather, we demonstrate that interspecific hybrids closely resemble piRNA effector-protein mutants in both TE misregulation and aberrant piRNA production. We suggest that TE derepression in interspecific hybrids largely reflects adaptive divergence of piRNA pathway genes rather than species-specific differences in TE-derived piRNAs.

Eukaryotic genomes contain large quantities of transposable elements (TEs), short self-replicating DNA sequences that can move within the genome. The selfish replication of TEs has potentially drastic consequences for the host, such as disruption of gene function, induction of sterility, and initiation or exacerbation of some cancers. Like the adaptive immune system that defends our bodies against pathogens, the Piwi-interacting RNA (piRNA) pathway defends animal genomes against the harmful effects of TEs. Fundamental to piRNA-mediated defense is the production of small noncoding RNAs that act like antibodies to target replicating TEs for destruction by piRNA-effector proteins. piRNAs are expected to diverge rapidly between species in response to genome infection by increasingly disparate TEs. Here, we tested this hypothesis by examining how differences in piRNAs between two species of fruit fly relate to TE “immunity” in their hybrid offspring. Because piRNAs are maternally deposited, we expected excessive replication of paternal TEs in hybrids. Surprisingly, we observe increased activity of both maternal and paternal TEs, together with defects in piRNA production that are reminiscent of piRNA effector-protein mutants. Our observations reveal that piRNA effector-proteins do not function properly in hybrids, and we propose that adaptive evolution among piRNA effector-proteins contributes to host genome defense and leads to the functional incompatibilities that we observe in hybrids.

Long-term and short-term evolutionary impacts of transposable elements on Drosophila

Transposable elements (TEs) are considered to be genomic parasites and their interactions with their hosts have been likened to the coevolution between host and other nongenomic, horizontally transferred pathogens. TE families, however, are vertically inherited as integral segments of the nuclear genome. This transmission strategy has been suggested to weaken the selective benefits of host alleles repressing the transposition of specific TE variants. On the other hand, the elevated rates of TE transposition and high incidences of deleterious mutations observed during the rare cases of horizontal transfers of TE families between species could create at least a transient process analogous to the influence of horizontally transmitted pathogens. Here, we formally address this analogy, using empirical and theoretical analysis to specify the mechanism of how host-TE interactions may drive the evolution of host genes. We found that host TE-interacting genes actually have more pervasive evidence of adaptive evolution than immunity genes that interact with nongenomic pathogens in Drosophila. Yet, both our theoretical modeling and empirical observations comparing Drosophila melanogaster populations before and after the horizontal transfer of P elements, which invaded D. melanogaster early last century, demonstrated that horizontally transferred TEs have only a limited influence on host TE-interacting genes. We propose that the more prevalent and constant interaction with multiple vertically transmitted TE families may instead be the main force driving the fast evolution of TE-interacting genes, which is fundamentally different from the gene-for-gene interaction of host-pathogen coevolution.

Maize genetic diversity and association mapping using transposable element insertion polymorphisms

Transposable elements are the major component of the maize genome and presumably highly polymorphic yet they have not been used in population genetics and association analyses. Using the Transposon Display method, we isolated and converted into PCR-based markers 33 Miniature Inverted Repeat Transposable Elements (MITE) polymorphic insertions. These polymorphisms were genotyped on a population-based sample of 26 American landraces for a total of 322 plants. Genetic diversity was high and partitioned within and among landraces. The genetic groups identified using Bayesian clustering were in agreement with published data based on SNPs and SSRs, indicating that MITE polymorphisms reflect maize genetic history. To explore the contribution of MITEs to phenotypic variation, we undertook an association mapping approach in a panel of 367 maize lines phenotyped for 26 traits. We found a highly significant association between the marker ZmV1-9, on chromosome 1, and male flowering time. The variance explained by this association is consistent with a flowering delay of +123 degree-days. This MITE insertion is located at only 289 nucleotides from the 3' end of a Cytochrome P450-like gene, a region that was never identified in previous association mapping or QTL surveys. Interestingly, we found (i) a non-synonymous mutation located in the exon 2 of the gene in strong linkage disequilibrium with the MITE polymorphism, and (ii) a perfect sequence homology between the MITE sequence and a maize siRNA that could therefore potentially interfere with the expression of the Cytochrome P450-like gene. Those two observations among others offer exciting perspectives to validate functionally the role of this region on phenotypic variation.

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.

Transposable elements (TE's) are parasitic genetic elements that spread by replicating themselves within a host genome. Most organisms are burdened with transposable elements; in fact, up to 80% of some genomes can consist of TE–derived DNA. Here, we use new sequencing technology to examine variation in genomic TE composition within a population at a finer scale and in a more unbiased fashion than has been possible before. We study a Portuguese population of D. melanogaster and find a large number of TE insertions, most of which occur in few individuals. Our analysis confirms that TE insertions are subject to purifying selection that counteracts their spread, and it suggests that the genome records waves of past TE invasions, with recently active elements occurring at low population frequency. We also find indications that TE insertions may sometimes have beneficial effects.

Spatiotemporal transcription of the P element and the 412 retrotransposon during embryogenesis of Drosophila melanogaster and D. willistoni

Transposable elements (TEs) are mobile nucleotide sequences which, through changing position in host genomes, partake in important evolutionary processes. The expression patterns of two TEs, P element transposon and 412 retrotransposon, were investigated during Drosophila melanogaster and D. willistoni embryogenesis, by means of embryo hybridization using riboprobes. Spatiotemporal transcription patterns for both TEs were similar to those of developmental genes. Although the two species shared the same P element transcription pattern, this was not so with 412 retrotransposon. These findings suggest that the regulatory sequences involved in the initial development of Drosophila spp are located in the transposable element sequences, and differences, such as in this case of the 412 retrotransposon, lead to losses or changes in their transcription patterns.

Evolution of genome content: population dynamics of transposable elements in flies and humans

Recent research is starting to shed light on the factors that influence the population and evolutionary dynamics of transposable elements (TEs) and TE life cycles. Genomes differ sharply in the number of TE copies, in the level of TE activity, in the diversity of TE families and types, and in the proportion of old and young TEs. In this chapter, we focus on two well-studied genomes with strikingly different architectures, humans and Drosophila, which represent two extremes in terms of TE diversity and population dynamics. We argue that some of the answers might lie in (1) the larger population size and consequently more effective selection against new TE insertions due to ectopic recombination in flies compared to humans; and (2) in the faster rate of DNA loss in flies compared to humans leading to much faster removal of fixed TE copies from the fly genome.

Active transposition in genomes

Transposons are DNA sequences capable of moving in genomes. Early evidence showed their accumulation in many species and suggested their continued activity in at least isolated organisms. In the past decade, with the development of various genomic technologies, it has become abundantly clear that ongoing activity is the rule rather than the exception. Active transposons of various classes are observed throughout plants and animals, including humans. They continue to create new insertions, have an enormous variety of structural and functional impact on genes and genomes, and play important roles in genome evolution. Transposon activities have been identified and measured by employing various strategies. Here, we summarize evidence of current transposon activity in various plant and animal genomes.

Heterozygote advantage as a natural consequence of adaptation in diploids

Molecular adaptation is typically assumed to proceed by sequential fixation of beneficial mutations. In diploids, this picture presupposes that for most adaptive mutations, the homozygotes have a higher fitness than the heterozygotes. Here, we show that contrary to this expectation, a substantial proportion of adaptive mutations should display heterozygote advantage. This feature of adaptation in diploids emerges naturally from the primary importance of the fitness of heterozygotes for the invasion of new adaptive mutations. We formalize this result in the framework of Fisher's influential geometric model of adaptation. We find that in diploids, adaptation should often proceed through a succession of short-lived balanced states that maintain substantially higher levels of phenotypic and fitness variation in the population compared with classic adaptive walks. In fast-changing environments, this variation produces a diversity advantage that allows diploids to remain better adapted compared with haploids despite the disadvantage associated with the presence of unfit homozygotes. The short-lived balanced states arising during adaptive walks should be mostly invisible to current scans for long-term balancing selection. Instead, they should leave signatures of incomplete selective sweeps, which do appear to be common in many species. Our results also raise the possibility that balancing selection, as a natural consequence of frequent adaptation, might play a more prominent role among the forces maintaining genetic variation than is commonly recognized.

Host-parasite relationships in the genome

Transposable elements are best interpreted as genomic parasites, proliferating in genomes through their over-replication relative to the rest of the genome. A new study examining correlations across Drosophila species between transposable element numbers and rates of host evolution has brought into focus one of the most complex questions in transposable element biology-what it is that determines the proportion of the genome that is transposable elements.

See research article: http://www.biomedcentral.com/1471-2148/11/258/

Coexistence of trichome variation in a natural plant population: a combined study using ecological and candidate gene approaches

The coexistence of distinct phenotypes within populations has long been investigated in evolutionary ecology. Recent studies have identified the genetic basis of distinct phenotypes, but it is poorly understood how the variation in candidate loci is maintained in natural environments. In this study, we examined fitness consequences and genetic basis of variation in trichome production in a natural population of Arabidopsis halleri subsp. gemmifera. Half of the individuals in the study population produced trichomes while the other half were glabrous, and the leaf beetle Phaedon brassicae imposed intensive damage to both phenotypes. The fitness of hairy and glabrous plants showed no significant differences in the field during two years. A similar result was obtained when sibling hairy and glabrous plants were transplanted at the same field site, whereas a fitness cost of trichome production was detected under a weak herbivory condition. Thus, equivalent fitness of hairy and glabrous plants under natural herbivory allows their coexistence in the contemporary population. The pattern of polymorphism of the candidate trichome gene GLABROUS1 (GL1) showed no evidence of long-term maintenance of trichome variation within the population. Although balancing selection under fluctuating biotic environments is often proposed to explain the maintenance of defense variation, the lack of clear evidence of balancing selection in the study population suggests that other factors such as gene flow and neutral process may have played relatively large roles in shaping trichome variation at least for the single population level.

Osvaldo and Isis retrotransposons as markers of the Drosophila buzzatii colonisation in Australia

Background

Transposable elements (TEs) constitute an important source of genetic variability owing to their jumping and regulatory properties, and are considered to drive species evolution. Several factors that are able to induce TE transposition in genomes have been documented (for example environmental stress and inter- and intra-specific crosses) but in many instances the reasons for TE mobilisation have yet to be elucidated. Colonising populations constitute an ideal model for studying TE behaviour and distribution as they are exposed to different environmental and new demographic conditions. In this study, the distribution of two TEs, Osvaldo and Isis, was examined in two colonising populations of D. buzzatii from Australia. Comparing Osvaldo copy numbers between Australian and Old World (reported in previous studies) colonisations provides a valuable tool for elucidating the colonisation process and the effect of new conditions encountered by colonisers on TEs.

Results

The chromosomal distributions of Osvaldo and Isis retrotransposons in two colonising populations of D. buzzatii from Australia revealed sites of high insertion frequency (>10%) and low frequency sites. Comparisons between Osvaldo insertion profiles in colonising populations from the Old World and Australia demonstrate a tendency towards a higher number of highly occupied sites with higher insertion frequency in the Old World than in Australian populations. Tests concerning selection against deleterious TE insertions indicate that Isis is more controlled by purifying selection than Osvaldo. The distribution of both elements on chromosomal arms follows a Poisson distribution and there are non-significant positive correlations between highly occupied sites and chromosomal inversions.

Conclusions

The occupancy profile of Osvaldo and Isis retrotransposons is characterised by the existence of high and low insertion frequency sites in the populations. These results demonstrate that Australian D. buzzatii populations were subjected to a founder effect during the colonisation process. Moreover, there are more sites with high insertion frequency in the Old World colonisation than in the Australian colonisation, indicating a probable stronger bottleneck effect in Australia. The results suggest that selection does not seem to play a major role, compared to demography, in the distribution of transposable elements in the Australian populations.

Molecular basis of adaptive shift in body size in Drosophila melanogaster: functional and sequence analyses of the Dca gene

Latitudinal body size clines in animals conforming to Bergmann's rule occur on many continents but isolating their underlying genetic basis remains a challenge. In Drosophila melanogaster, the gene Dca accounts for approximately 5-10% of the natural wing size variation (McKechnie SW, Blacket MJ, Song SV, Rako L, Carroll X, Johnson TK, Jensen LT, Lee SF, Wee CW, Hoffmann AA. 2010. A clinally varying promoter polymorphism associated with adaptive variation in wing size in Drosophila. Mol Ecol. 19:775-784). We present here functional evidence that Dca is a negative regulator of wing size. A significant negative latitudinal cline of Dca gene expression was detected in synchronized third instar larvae. In addition, we clarified the evolutionary history of the three most common Dca promoter alleles (Dca237-1, Dca237-2, and Dca247) and showed that the insertion allele (Dca247), whose frequency increases with latitude, is associated with larger wing centroid size and higher average cell number in male flies. Finally, we showed that the overall linkage disequilibrium (LD) was low in the Dca promoter and that the insertion/deletion polymorphism that defines the Dca alleles was in strong LD with two other upstream sites. Our results provide strong support that Dca is a candidate for climatic adaptation in D. melanogaster.

Modification of gene duplicability during the evolution of protein interaction network

Duplications of genes encoding highly connected and essential proteins are selected against in several species but not in human, where duplicated genes encode highly connected proteins. To understand when and how gene duplicability changed in evolution, we compare gene and network properties in four species (Escherichia coli, yeast, fly, and human) that are representative of the increase in evolutionary complexity, defined as progressive growth in the number of genes, cells, and cell types. We find that the origin and conservation of a gene significantly correlates with the properties of the encoded protein in the protein-protein interaction network. All four species preserve a core of singleton and central hubs that originated early in evolution, are highly conserved, and accomplish basic biological functions. Another group of hubs appeared in metazoans and duplicated in vertebrates, mostly through vertebrate-specific whole genome duplication. Such recent and duplicated hubs are frequently targets of microRNAs and show tissue-selective expression, suggesting that these are alternative mechanisms to control their dosage. Our study shows how networks modified during evolution and contributes to explaining the occurrence of somatic genetic diseases, such as cancer, in terms of network perturbations.

The struggle for life of the genome's selfish architects

Transposable elements (TEs) were first discovered more than 50 years ago, but were totally ignored for a long time. Over the last few decades they have gradually attracted increasing interest from research scientists. Initially they were viewed as totally marginal and anecdotic, but TEs have been revealed as potentially harmful parasitic entities, ubiquitous in genomes, and finally as unavoidable actors in the diversity, structure, and evolution of the genome. Since Darwin's theory of evolution, and the progress of molecular biology, transposable elements may be the discovery that has most influenced our vision of (genome) evolution. In this review, we provide a synopsis of what is known about the complex interactions that exist between transposable elements and the host genome. Numerous examples of these interactions are provided, first from the standpoint of the genome, and then from that of the transposable elements. We also explore the evolutionary aspects of TEs in the light of post-Darwinian theories of evolution.

Microarray estimation of genomic inter-strain variability in the genus Ectocarpus (Phaeophyceae)

BACKGROUND

Brown algae of the genus Ectocarpus exhibit high levels of genetic diversity and variability in morphological and physiological characteristics. With the establishment of E. siliculosus as a model and the availability of a complete genome sequence, it is now of interest to analyze variability among different species, ecotypes, and strains of the genus Ectocarpus both at the genome and the transcriptome level.

RESULTS

We used an E. siliculosus gene expression microarray based on EST sequences from the genome-sequenced strain (reference strain) to carry out comparative genome hybridizations for five Ectocarpus strains: four E. siliculosus isolates (the male genome strain, a female strain used for outcrosses with the genome strain, a strain isolated from freshwater, and a highly copper-tolerant strain), as well as one strain of the sister species E. fasciculatus. Our results revealed significant genomic differences between ecotypes of the same species, and enable the selection of conserved probes for future microarray experiments with these strains. In the two closely related strains (a male and a female strain used for crosses), genomic differences were also detected, but concentrated in two smaller genomic regions, one of which corresponds to a viral insertion site.

CONCLUSION

The high variability between strains supports the concept of E. siliculosus as a complex of cryptic species. Moreover, our data suggest that several parts of the Ectocarpus genome may have evolved at different rates: high variability was detected particularly in transposable elements and fucoxanthin chlorophyll a/c binding proteins.

T-lex: a program for fast and accurate assessment of transposable element presence using next-generation sequencing data

Transposable elements (TEs) are repetitive DNA sequences that are ubiquitous, extremely abundant and dynamic components of practically all genomes. Much effort has gone into annotation of TE copies in reference genomes. The sequencing cost reduction and the newly available next-generation sequencing (NGS) data from multiple strains within a species offer an unprecedented opportunity to study population genomics of TEs in a range of organisms. Here, we present a computational pipeline (T-lex) that uses NGS data to detect the presence/absence of annotated TE copies. T-lex can use data from a large number of strains and returns estimates of population frequencies of individual TE insertions in a reasonable time. We experimentally validated the accuracy of T-lex detecting presence or absence of 768 previously identified TE copies in two resequenced Drosophila melanogaster strains. Approximately 95% of the TE insertions were detected with 100% sensitivity and 97% specificity. We show that even at low levels of coverage T-lex produces accurate results for TE copies that it can identify reliably but that the rate of 'no data' calls increases as the coverage falls below 15×. T-lex is a broadly applicable and flexible tool that can be used in any genome provided the availability of the reference genome, individual TE copy annotation and NGS data.

Adaptive divergence between incipient species of Anopheles gambiae increases resistance to Plasmodium

The African malaria mosquito Anopheles gambiae is diversifying into ecotypes known as M and S forms. This process is thought to be promoted by adaptation to different larval habitats, but its genetic underpinnings remain elusive. To identify candidate targets of divergent natural selection in M and S, we performed genomewide scanning in paired population samples from Mali, followed by resequencing and genotyping from five locations in West, Central, and East Africa. Genome scans revealed a significant peak of M-S divergence on chromosome 3L, overlapping five known or suspected immune response genes. Resequencing implicated a selective target at or near the TEP1 gene, whose complement C3-like product has antiparasitic and antibacterial activity. Sequencing and allele-specific genotyping showed that an allelic variant of TEP1 has been swept to fixation in M samples from Mali and Burkina Faso and is spreading into neighboring Ghana, but is absent from M sampled in Cameroon, and from all sampled S populations. Sequence comparison demonstrates that this allele is related to, but distinct from, TEP1 alleles of known resistance phenotype. Experimental parasite infections of advanced mosquito intercrosses demonstrated a strong association between this TEP1 variant and resistance to both rodent malaria and the native human malaria parasite Plasmodium falciparum. Although malaria parasites may not be direct agents of pathogen-mediated selection at TEP1 in nature--where larvae may be the more vulnerable life stage--the process of adaptive divergence between M and S has potential consequences for malaria transmission.

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.

Recent and recurrent selective sweeps of the antiviral RNAi gene Argonaute-2 in three species of Drosophila

Antagonistic host–parasite interactions can drive rapid adaptive evolution in genes of the immune system, and such arms races may be an important force shaping polymorphism in the genome. The RNA interference pathway gene Argonaute-2 (AGO2) is a key component of antiviral defense in Drosophila, and we have previously shown that genes in this pathway experience unusually high rates of adaptive substitution. Here we study patterns of genetic variation in a 100-kbp region around AGO2 in three different species of Drosophila. Our data suggest that recent independent selective sweeps in AGO2 have reduced genetic variation across a region of more than 50 kbp in Drosophila melanogaster, D. simulans, and D. yakuba, and we estimate that selection has fixed adaptive substitutions in this gene every 30–100 thousand years. The strongest signal of recent selection is evident in D. simulans, where we estimate that the most recent selective sweep involved an allele with a selective advantage of the order of 0.5–1% and occurred roughly 13–60 Kya. To evaluate the potential consequences of the recent substitutions on the structure and function of AGO2, we used fold-recognition and homology-based modeling to derive a structural model for the Drosophila protein, and this suggests that recent substitutions in D. simulans are overrepresented at the protein surface. In summary, our results show that selection by parasites can consistently target the same genes in multiple species, resulting in areas of the genome that have markedly reduced genetic diversity.

High level of structural polymorphism driven by mobile elements in the Hox genomic region of the Chaetognath Spadella cephaloptera

Little is known about the relationships between genome polymorphism, mobile element dynamics, and population size among animal populations. The chaetognath species Spadella cephaloptera offers a unique perspective to examine this issue because they display a high level of genetic polymorphism at the population level. Here, we have investigated in detail the extent of nucleotide and structural polymorphism in a region harboring Hox1 and several coding genes and presumptive functional elements. Sequencing of several bacterial artificial chromosome inserts representative of this nuclear region uncovered a high level of structural heterogeneity, which is mainly caused by the polymorphic insertion of a diversity of genetic mobile elements. By anchoring this variation through individual genotyping, we demonstrated that sequence diversity could be attributed to the allelic pool of a single population, which was confirmed by detection of extensive recombination within the genomic region studied. The high average level of nucleotide heterozygosity provides clues of selection in both coding and noncoding domains. This pattern stresses how selective processes remarkably cope with intense sequence turnover due to substitutions, mobile element insertions, and recombination to preserve the integrity of functional landscape. These findings suggest that genome polymorphism could provide pivotal information for future functional annotation of genomes.

Transposable elements in natural populations of Drosophila melanogaster

Transposable elements (TEs) are families of small DNA sequences found in the genomes of virtually all organisms. The sequences typically encode essential components for the replicative transposition sequences of that TE family. Thus, TEs are simply genomic parasites that inflict detrimental mutations on the fitness of their hosts. Several models have been proposed for the containment of TE copy number in outbreeding host populations such as Drosophila. Surveys of the TEs in genomes from natural populations of Drosophila have played a central role in the investigation of TE dynamics. The early surveys indicated that a typical TE insertion is rare in a population, which has been interpreted as evidence that each TE is selected against. The proposed mechanisms of this natural selection are reviewed here. Subsequent and more targeted surveys identify heterogeneity among types of TEs and also highlight the large role of homologous and possibly ectopic crossing over in the dynamics of the Drosophila TEs. The recent discovery of germline-specific RNA interference via the piwi-interacting RNA pathway opens yet another interesting mechanism that may be critical in containing the copy number of TEs in natural populations of Drosophila. The expected flood of Drosophila population genomics is expected to rapidly advance understanding of the dynamics of TEs.

Evidence that adaptation in Drosophila is not limited by mutation at single sites

Adaptation in eukaryotes is generally assumed to be mutation-limited because of small effective population sizes. This view is difficult to reconcile, however, with the observation that adaptation to anthropogenic changes, such as the introduction of pesticides, can occur very rapidly. Here we investigate adaptation at a key insecticide resistance locus (Ace) in Drosophila melanogaster and show that multiple simple and complex resistance alleles evolved quickly and repeatedly within individual populations. Our results imply that the current effective population size of modern D. melanogaster populations is likely to be substantially larger (> or = 100-fold) than commonly believed. This discrepancy arises because estimates of the effective population size are generally derived from levels of standing variation and thus reveal long-term population dynamics dominated by sharp--even if infrequent--bottlenecks. The short-term effective population sizes relevant for strong adaptation, on the other hand, might be much closer to census population sizes. Adaptation in Drosophila may therefore not be limited by waiting for mutations at single sites, and complex adaptive alleles can be generated quickly without fixation of intermediate states. Adaptive events should also commonly involve the simultaneous rise in frequency of independently generated adaptive mutations. These so-called soft sweeps have very distinct effects on the linked neutral polymorphisms compared to the standard hard sweeps in mutation-limited scenarios. Methods for the mapping of adaptive mutations or association mapping of evolutionarily relevant mutations may thus need to be reconsidered.

Genome-wide patterns of adaptation to temperate environments associated with transposable elements in Drosophila

Investigating spatial patterns of loci under selection can give insight into how populations evolved in response to selective pressures and can provide monitoring tools for detecting the impact of environmental changes on populations. Drosophila is a particularly good model to study adaptation to environmental heterogeneity since it is a tropical species that originated in sub-Saharan Africa and has only recently colonized the rest of the world. There is strong evidence for the adaptive role of Transposable Elements (TEs) in the evolution of Drosophila, and TEs might play an important role specifically in adaptation to temperate climates. In this work, we analyzed the frequency of a set of putatively adaptive and putatively neutral TEs in populations with contrasting climates that were collected near the endpoints of two known latitudinal clines in Australia and North America. The contrasting results obtained for putatively adaptive and putatively neutral TEs and the consistency of the patterns between continents strongly suggest that putatively adaptive TEs are involved in adaptation to temperate climates. We integrated information on population behavior, possible environmental selective agents, and both molecular and functional information of the TEs and their nearby genes to infer the plausible phenotypic consequences of these insertions. We conclude that adaptation to temperate environments is widespread in Drosophila and that TEs play a significant role in this adaptation. It is remarkable that such a diverse set of TEs located next to a diverse set of genes are consistently adaptive to temperate climate-related factors. We argue that reverse population genomic analyses, as the one described in this work, are necessary to arrive at a comprehensive picture of adaptation.

The potential of geographic studies of genetic variation for the understanding of adaptation has been recognized for some time. In Drosophila, most of the available studies are based on a priori candidates giving a biased picture of the genes and traits under spatially varying selection. In this work, we performed a genome-wide scan of adaptations to temperate climates associated with Transposable Element (TE) insertions. We integrated the available information of the identified TEs and their nearby genes to provide plausible hypotheses about the phenotypic consequences of these insertions. Considering the diversity of these TEs and the variety of genes into which they are inserted, it is surprising that their adaptive effects are consistently related to temperate climate-related factors. The TEs identified in this work add substantially to the markers available to monitor the impact of climate change on populations.

Contrasting patterns of transposable element insertions in Drosophila heat-shock promoters

The proximal promoter regions of heat-shock genes harbor a remarkable number of P transposable element (TE) insertions relative to both positive and negative control proximal promoter regions in natural populations of Drosophila melanogaster. We have screened the sequenced genomes of 12 species of Drosophila to test whether this pattern is unique to these populations. In the 12 species' genomes, transposable element insertions are no more abundant in promoter regions of single-copy heat-shock genes than in promoters with similar or dissimilar architecture. Also, insertions appear randomly distributed across the promoter region, whereas insertions clustered near the transcription start site in promoters of single-copy heat-shock genes in D. melanogaster natural populations. Hsp70 promoters exhibit more TE insertions per promoter than all other genesets in the 12 species, similarly to in natural populations of D. melanogaster. Insertions in the Hsp70 promoter region, however, cluster away from the transcription start site in the 12 species, but near it in natural populations of D. melanogaster. These results suggest that D. melanogaster heat-shock promoters are unique in terms of their interaction with transposable elements, and confirm that Hsp70 promoters are distinctive in TE insertions across Drosophila.

Comparative Genomics on the Drosophila Phylogenetic Tree

With the sequencing of 12 complete euchromatic Drosophila genomes, the genus Drosophila is a leading model for comparative genomics. In this review, we discuss the novel insights into evolutionary processes afforded by the newly available genomic sequences when placed in the context of the phylogeny. We focus on three levels: insights into whole-genome content, such as changes in genome size and content across the phylogeny; insights into large-scale patterns of divergence and conservation, such as selective constraints on genes and chromosome-level evolution of sex chromosomes; and insights into finer-scale processes in individual lineages and genes, such as lineage-specific evolution in response to ecological context. As the field of comparative genomics is still young, we also discuss current challenges, such as the development of more sophisticated evolutionary models to capture nonequilibrium processes and the improvement of assembly and alignment algorithms to better capture uncertainty in the data.

Population dynamics of PIWI-interacting RNAs (piRNAs) and their targets in Drosophila

Transposable elements (TEs) are mobile DNA sequences that make up a large fraction of eukaryotic genomes. Recently it was discovered that PIWI-interacting RNAs (piRNAs), a class of small RNA molecules that are mainly generated from transposable elements, are crucial repressors of active TEs in the germline of fruit flies. By quantifying expression levels of 32 TE families in piRNA pathway mutants relative to wild-type fruit flies, we provide evidence that piRNAs can severely silence the activities of retrotransposons. We incorporate piRNAs into a population genetic framework for retrotransposons and perform forward simulations to model the population dynamics of piRNA loci and their targets. Using parameters optimized for Drosophila melanogaster, our simulation results indicate that (1) piRNAs can significantly reduce the fitness cost of retrotransposons; (2) retrotransposons that generate piRNAs (piRTs) are selectively more advantageous, and such retrotransposon insertions more easily attain high frequency or fixation; (3) retrotransposons that are repressed by piRNAs (targetRTs), however, also have an elevated probability of reaching high frequency or fixation in the population because their deleterious effects are attenuated. By surveying the polymorphisms of piRT and targetRT insertions across nine strains of D. melanogaster, we verified these theoretical predictions with population genomic data. Our theoretical and empirical analysis suggests that piRNAs can significantly increase the fitness of individuals that bear them; however, piRNAs may provide a shelter or Trojan horse for retrotransposons, allowing them to increase in frequency in a population by shielding the host from the deleterious consequences of retrotransposition.

RNA editing: a driving force for adaptive evolution?

Genetic variability is considered a key to the evolvability of species. The conversion of an adenosine (A) to inosine (I) in primary RNA transcripts can result in an amino acid change in the encoded protein, a change in secondary structure of the RNA, creation or destruction of a splice consensus site, or otherwise alter RNA fate. Substantial transcriptome and proteome variability is generated by A-to-I RNA editing through site-selective post-transcriptional recoding of single nucleotides. We posit that this epigenetic source of phenotypic variation is an unrecognized mechanism of adaptive evolution. The genetic variation introduced through editing occurs at low evolutionary cost since predominant production of the wild-type protein is retained. This property even allows exploration of sequence space that is inaccessible through mutation, leading to increased phenotypic plasticity and provides an evolutionary advantage for acclimatization as well as long-term adaptation. Furthermore, continuous probing for novel RNA editing sites throughout the transcriptome is an intrinsic property of the editing machinery and represents the molecular basis for increased adaptability. We propose that higher organisms have therefore evolved to systems with increasing RNA editing activity and, as a result, to more complex systems.

The hobo transposon and hobo-related elements are expressed as developmental genes in Drosophila

Transposable elements comprise a significant part of genomes and are involved in their evolvability. The hobo element is found as an active class II transposable element in Drosophila melanogaster that is able to induce gonadal dysgenesis. Some hobo-related sequences (hRSs) are thought to be relics of old "hobo" invasions, and are therefore ancient genomic constituents. However, some of these hRSs are still mobile. The present study analyzed the expression pattern of hobo and a particular type of hRSs, hobo(VAHS). Both elements were shown to be expressed as sense and antisense mRNA transcripts. Expression analysis in whole mount embryos revealed a pattern similar to that of some developmental regulatory genes. Here we suggest that cis-regulatory sequences similar to those in developmental genes exist in hobo sequences. Therefore, hobo mobilization may contribute to the development of new regulatory networks during genomic evolution.

The adaptive role of transposable elements in the Drosophila genome

Transposable elements (TEs) are short DNA sequences with the capacity to move between different sites in the genome. This ability provides them with the capacity to mutate the genome in many different ways, from subtle regulatory mutations to gross genomic rearrangements. The potential adaptive significance of TEs was recognized by those involved in their initial discovery although it was hotly debated afterwards. For more than two decades, TEs were considered to be intragenomic parasites leading to almost exclusively detrimental effects to the host genome. The sequencing of the Drosophila melanogaster genome provided an unprecedented opportunity to study TEs and led to the identification of the first TE-induced adaptations in this species. These studies were followed by a systematic genome-wide search for adaptive insertions that allowed for the first time to infer that TEs contribute substantially to adaptive evolution. This study also revealed that there are at least twice as many TE-induced adaptations that remain to be identified. To gain a better understanding of the adaptive role of TEs in the genome we clearly need to (i) identify as many adaptive TEs as possible in a range of Drosophila species as well as (ii) carry out in-depth investigations of the effects of adaptive TEs on as many phenotypes as possible.

A recent adaptive transposable element insertion near highly conserved developmental loci in Drosophila melanogaster

A recent genomewide screen identified 13 transposable elements that are likely to have been adaptive during or after the spread of Drosophila melanogaster out of Africa. One of these insertions, Bari-Juvenile hormone epoxy hydrolase (Bari-Jheh), was associated with the selective sweep of its flanking neutral variation and with reduction of expression of one of its neighboring genes: Jheh3. Here, we provide further evidence that Bari-Jheh insertion is adaptive. We delimit the extent of the selective sweep and show that Bari-Jheh is the only mutation linked to the sweep. Bari-Jheh also lowers the expression of its other flanking gene, Jheh2. Subtle consequences of Bari-Jheh insertion on life-history traits are consistent with the effects of reduced expression of the Jheh genes. Finally, we analyze molecular evolution of Jheh genes in both the long- and the short-term and conclude that Bari-Jheh appears to be a very rare adaptive event in the history of these genes. We discuss the implications of these findings for the detection and understanding of adaptation.

Identification of expressed transposable element insertions in the sequenced genome of Drosophila melanogaster

Transposable elements (TEs) are major components of most genomes, and their impact on genome evolution is now well documented. However, the way they affect the transcriptome is still not clearly established. Using the sequenced genome of Drosophila melanogaster and EST libraries, we describe here the TE insertions that are unequivocally transcribed, and we have determined their location in the sequenced genome. We show that most TE families are transcribed, and we have specifically identified 69 expressed TE insertions, half of which are located inside genes, mostly within introns and 5'UTRs.