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.     

Drosophila Interspecific Hybrids Phenocopy piRNA-Pathway Mutants

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.     

Patterns of puffing activity in the salivary gland chromosomes of Drosophila

The autosomal salivary gland chromosome puffing patterns of Drosophila simulans are described and compared with the puffing patterns of the sibling species D. melanogaster. During the late third larval instar and the prepupal period the patterns of puffing activity of these two species are similar — approximately 50% of the puffs common to both species showing identical activities. The remaining puffs differ in their timing of activity, or in their mean sizes, or in both of these parameters. A number of puffs (14) found in D. simulans have not been regularly observed in the Oregon stock of D. melanogaster but are active in other D. melanogaster strains. One puff (46 A) of D. melanogaster was absent from D. simulans and forms a heterozygous puff in hybrids, when the homologous chromosomes are synapsed. When the homologues are asynapsed a puff at 46 A is restricted to the melanogaster homologue. The puff at 63E on chromosome arm 3L is considerably smaller in D. simulans than in D. melanogaster and this size difference is autonomous in hybrids. Other puffs not common to both species behave non-autonomously in the species hybrid, even when the homologous chromosomes are asynapsed.     

An Age-of-Allele Test of Neutrality for Transposable Element Insertions

How natural selection acts to limit the proliferation of transposable elements (TEs) in genomes has been of interest to evolutionary biologists for many years. To describe TE dynamics in populations, previous studies have used models of transposition–selection equilibrium that assume a constant rate of transposition. However, since TE invasions are known to happen in bursts through time, this assumption may not be reasonable. Here we propose a test of neutrality for TE insertions that does not rely on the assumption of a constant transposition rate. We consider the case of TE insertions that have been ascertained from a single haploid reference genome sequence. By conditioning on the age of an individual TE insertion allele (inferred by the number of unique substitutions that have occurred within the particular TE sequence since insertion), we determine the probability distribution of the insertion allele frequency in a population sample under neutrality. Taking models of varying population size into account, we then evaluate predictions of our model against allele frequency data from 190 retrotransposon insertions sampled from North American and African populations of Drosophila melanogaster. Using this nonequilibrium neutral model, we are able to explain ∼80% of the variance in TE insertion allele frequencies based on age alone. Controlling for both nonequilibrium dynamics of transposition and host demography, we provide evidence for negative selection acting against most TEs as well as for positive selection acting on a small subset of TEs. Our work establishes a new framework for the analysis of the evolutionary forces governing large insertion mutations like TEs, gene duplications, or other copy number variants.     

Transposable elements in individual genotypes of Drosophila simulans

Transposable elements are abundant, dynamic components of the genome that affect organismal phenotypes and fitness. In Drosophila melanogaster, they have increased in abundance as the species spread out of Africa, and different populations differ in their transposable element content. However, very little is currently known about how transposable elements differ between individual genotypes, and how that relates to the population dynamics of transposable elements overall. The sister species of D. melanogaster, D. simulans, has also recently become cosmopolitan, and panels of inbred genotypes exist from cosmopolitan and African flies. Therefore, we can determine whether the differences in colonizing populations are repeated in D. simulans, what the dynamics of transposable elements are in individual genotypes, and how that compares to wild flies. After estimating copy number in cosmopolitan and African D. simulans, I find that transposable element load is higher in flies from cosmopolitan populations. In addition, transposable element load varies considerably between populations, between genotypes, but not overall between wild and inbred lines. Certain genotypes either contain active transposable elements or are more permissive of transposition and accumulate copies of particular transposable elements. Overall, it is important to quantify genotype‐specific transposable element dynamics as well as population averages to understand the dynamics of transposable element accumulation over time.     

Clonal analysis in hybrids between Drosophila melanogaster and Drosophila simulans

We have analysed the viability of cellular clones induced by mitotic recombination in Drosophila melanogaster/D. simulans hybrid females during larval growth. These clones contain a portion of either melanogaster or simulans genomes in homozygosity. Analysis has been carried out for the X and the second chromosomes, as well as for the 3L chromosome arm. Clones were not found in certain structures, and in others they appeared in a very low frequency. Only in abdominal tergites was a significant number of clones observed, although their frequency was lower than in melanogaster abdomens. The bigger the portion of the genome that is homozygous, the less viable is the recombinant melano-gaster/simulans hybrid clone. The few clones that appeared may represent cases in which mitotic recombination took place in distal chromosome intervals, so that the clones contained a small portion of either melanogaster or simulans chromosomes in homozygosity. Moreover, Lhr, a gene of D. simulans that suppresses the lethality of male and female melanogaster/simulans hybrids, does not suppress the lethality of the recombinant melanogaster/simulans clones. Thus, it appears that there is not just a single gene, but at least one per tested chromosome arm (and maybe more) that cause hybrid lethality. Therefore, the two species, D. melanogaster and D. simulans, have diverged to such a degree that the absence of part of the genome of one species cannot be substituted by the corresponding part of the genome of the other, probably due to the formation of co-adapted gene complexes in both species following their divergent evolution after speciation. The disruption of those coadapted gene complexes would cause the lethality of the recombinant hybrid clones.     

Divergent strategies in low temperature environment for the sibling species Drosophila melanogaster and D. simulans: overwintering in extension border areas of France and comparison with African populations

The two sibling species D. /melanogaster and D. simulans adopt different overwintering strategies in northern border areas situated in France. If the winter is mild, both species reappear in early spring to refound the population. If the winter is cold, with several weeks of temperatures below 0 °C, D. melanogaster leave their shelters in April whereas D. simulans, which do not use shelters, reappear in late June, probably after returning from further south. Here, we tried to identify life-history characteristics responsible for this difference. For this, we studied developmental duration, viability, fecundity, fertility and longevity, and compared the abilities of French and African populations to survive when food supplies were inadequate, at different temperatures (14, 11, 7 °C). These temperatures are lower than those commonly used in the laboratory but closer to real conditions encountered in the wild. When the temperature was mild (14 or 11 °C) and the food supply was adequate, D. simulans performed better than D. melanogaster: it had a higher fecundity, a longer life expectancy and the males remained fertile, allowing outdoor reproduction late in winter. However, D. simulans was less resistant in more extreme conditions. At 7 °C D. simulans survived shorter on normal medium and its ability to survive when food supplies were inadequate was insufficient to allow outdoor overwintering. In contrast, D. melanogaster could not reproduce during winter: its fecundity was low and males were sterile at 11 °C. Nevertheless, if only protein-deficient resources were available, temperate D. melanogaster could survive for longer than D. simulans at all the temperatures tested. This greater resistance to underfeeding allows the species to survive until spring, in shelters for several months. A comparison of French and African population performances showed differences in the evolution of the two species during the colonization of more northern areas. African D. simulans, which are efficient at mild temperatures, underwent few modifications. In contrast, the viability of D. melanogaster improved at low developmental temperatures. This species also displayed higher fecundity, longer survival and higher underfeeding resistance at low temperatures. The relationship between the long retention genotype and underfeeding resistance or survival ability observed in French D. melanogaster populations may not exist in African populations.     

Tirant Stealthily Invaded Natural Drosophila melanogaster Populations during the Last Century

It was long thought that solely three different transposable elements (TEs)—the I-element, the P-element, and hobo—invaded natural Drosophila melanogaster populations within the last century. By sequencing the “living fossils” of Drosophila research, that is, D. melanogaster strains sampled from natural populations at different time points, we show that a fourth TE, Tirant, invaded D. melanogaster populations during the past century. Tirant likely spread in D. melanogaster populations around 1938, followed by the I-element, hobo, and, lastly, the P-element. In addition to the recent insertions of the canonical Tirant, D. melanogaster strains harbor degraded Tirant sequences in the heterochromatin which are likely due to an ancient invasion, likely predating the split of D. melanogaster and D. simulans. These degraded insertions produce distinct piRNAs that were unable to prevent the novel Tirant invasion. In contrast to the I-element, P-element, and hobo, we did not find that Tirant induces any hybrid dysgenesis symptoms. This absence of apparent phenotypic effects may explain the late discovery of the Tirant invasion. Recent Tirant insertions were found in all investigated natural populations. Populations from Tasmania carry distinct Tirant sequences, likely due to a founder effect. By investigating the TE composition of natural populations and strains sampled at different time points, insertion site polymorphisms, piRNAs, and phenotypic effects, we provide a comprehensive study of a natural TE invasion.     

A Comprehensive Study of Genic Variation in Natural Populations of Drosophila melanogaster. III. Variations in Genetic Structure and Their Causes between Drosophila melanogaster and Its Sibling Species Drosophila simulans

The natural populations of Drosophila melanogaster and Drosophila simulans were compared for their genetic structure. A total of 114 gene-protein loci were studied in four mainland (from Europe and Africa) and an island (Seychelle) populations of D. simulans and the results were compared with those obtained on the same set of homologous loci in fifteen worldwide populations of D. melanogaster. The main results are as follows: (1) D. melanogaster shows a significantly higher proportion of loci polymorphic than D. simulans (52% vs. 39%, P<0.05), (2) both species have similar mean heterozygosity and mean number of alleles per locus, (3) the two species share some highly polymorphic loci but they do not share loci that show high geographic differentiation, and (4) D. simulans shows significantly less geographic differentiation than D. melanogaster. The differences in genetic differentiation between the two species are limited to loci located on the X and second chromosomes only; loci on the third chromosome show similar level of geographic differentiation in both species. These two species have previously been shown to differ in their pattern of variation for chromosomal polymorphisms, quantitative and physiological characters, two-dimensional electrophoretic (2DE) proteins, middle repetitive DNA and mitochondrial DNA. Variation in niche-widths and/or genetic "strategies" of adaptation appear to be the main causes of differences in the genetic structure of these two species.     

TEMP: a computational method for analyzing transposable element polymorphism in populations

Insertions and excisions of transposable elements (TEs) affect both the stability and variability of the genome. Studying the dynamics of transposition at the population level can provide crucial insights into the processes and mechanisms of genome evolution. Pooling genomic materials from multiple individuals followed by high-throughput sequencing is an efficient way of characterizing genomic polymorphisms in a population. Here we describe a novel method named TEMP, specifically designed to detect TE movements present with a wide range of frequencies in a population. By combining the information provided by pair-end reads and split reads, TEMP is able to identify both the presence and absence of TE insertions in genomic DNA sequences derived from heterogeneous samples; accurately estimate the frequencies of transposition events in the population and pinpoint junctions of high frequency transposition events at nucleotide resolution. Simulation data indicate that TEMP outperforms other algorithms such as PoPoolationTE, RetroSeq, VariationHunter and GASVPro. TEMP also performs well on whole-genome human data derived from the 1000 Genomes Project. We applied TEMP to characterize the TE frequencies in a wild Drosophila melanogaster population and study the inheritance patterns of TEs during hybrid dysgenesis. We also identified sequence signatures of TE insertion and possible molecular effects of TE movements, such as altered gene expression and piRNA production. TEMP is freely available at github: https://github.com/JialiUMassWengLab/TEMP.git.     

ADAPTATION AND SPECIALIZATION IN A TWO-RESOURCE ENVIRONMENT IN DROSOPHILA SPECIES

We assayed two components of performance (development time and survivorship), on food medium with and without ethanol, in laboratory populations of Drosophila simulans and D. melanogaster for which ethanol-medium was a novel food resource. These assays were done before and after 12 generations of rearing in either one-(regular medium only) or two-resource (regular medium and ethanol medium) environments. Initially, D. simulans was highly susceptible to ethanol, whereas D. melanogaster was relatively unaffected. After 12 generations in the two-resource environment, D. simulans showed significantly improved mean performance on ethanol medium; mean performance of D. melanogaster did not significantly change. Variation among families for both traits was higher on ethanol medium in D. simulans. Variation in D. melanogaster was not significantly affected by ethanol level, suggesting that resource quality was more important than novelty per se. In both species, the least variation was seen in populations after 12 generations in the two-resource environment. For development time in D. simulans, the decrease in variation was largely due to reduced variation within families, suggesting the evolution of canalization. Development time on the two media was not negatively correlated. In D. simulans, correlations measured before and after the experiment were not heterogeneous, suggesting that trade-offs in performance did not block diet expansion. In D. melanogaster, correlations became significantly less positive after 12 generations in the two-resource environment, supporting the view that correlations between performance on different resources may become less positive over time through selection.     

Rapid evolution at the Drosophila telomere: transposable element dynamics at an intrinsically unstable locus

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

Population genomics: whole-genome analysis of polymorphism and divergence in Drosophila simulans

The population genetic perspective is that the processes shaping genomic variation can be revealed only through simultaneous investigation of sequence polymorphism and divergence within and between closely related species. Here we present a population genetic analysis of Drosophila simulans based on whole-genome shotgun sequencing of multiple inbred lines and comparison of the resulting data to genome assemblies of the closely related species, D. melanogaster and D. yakuba. We discovered previously unknown, large-scale fluctuations of polymorphism and divergence along chromosome arms, and significantly less polymorphism and faster divergence on the X chromosome. We generated a comprehensive list of functional elements in the D. simulans genome influenced by adaptive evolution. Finally, we characterized genomic patterns of base composition for coding and noncoding sequence. These results suggest several new hypotheses regarding the genetic and biological mechanisms controlling polymorphism and divergence across the Drosophila genome, and provide a rich resource for the investigation of adaptive evolution and functional variation in D. simulans.     

Germ line abnormalities InDrosophila simulans transfected with the transposable P element

A strain ofDrosophila simulans was studied 40 generations after the transposable P element had been introduced into the genome by means of transformation. The genome also contained arosy transposon consisting of the wildtype allele of therosy gene flanked by P element DNA. During the 40 generations of evolution the number of P elements had increased to the level of 8–15 and the number ofrosy transposons to the level of 4–12. Continued transpositional activity in the germ line of the strain was evidenced by deletions occurring in therosy transposon and, in two independent sublines, by the transposition of therosy transposon from the X chromosome to the autosomes. Although at 25°C gonadal development and fertility appeared normal in both sexes, at 29°C both sexes were sterile. The sterile females had morphologically normal ovaries, but the sterile males often had shrunken, dysmorphic testes containing few or no immature sperm bundles. However, the sterility found in the transfected strain may not result directly from transpositional activity of the P element. The characteristics of theD. simulans strain infected with the P element are discussed in the context of factors that influence hybrid dysgenesis inD. melanogaster.     

Divergence of water balance mechanisms in two sibling species (Drosophila simulans and D. melanogaster): effects of growth temperatures

Drosophila simulans is more abundant under colder and drier montane habitats in the western Himalayas as compared to its sibling D. melanogaster but the mechanistic bases of such climatic adaptations are largely unknown. Previous studies have described D. simulans as a desiccation sensitive species which is inconsistent with its occurrence in temperate regions. We tested the hypothesis whether developmental plasticity of cuticular traits confers adaptive changes in water balance-related traits in the sibling species D. simulans and D. melanogaster. Our results are interesting in several respects. First, D. simulans grown at 15 °C possesses a high level of desiccation resistance in larvae (~39 h) and in adults (~86 h) whereas the corresponding values are quite low at 25 °C (larvae ~7 h; adults ~13 h). Interestingly, cuticular lipid mass was threefold higher in D. simulans grown at 15 °C as compared with 25 °C while there was no change in cuticular lipid mass in D. melanogaster. Second, developmental plasticity of body melanisation was evident in both species. Drosophila simulans showed higher melanisation at 15 °C as compared with D. melanogaster while the reverse trend was observed at 25 °C. Third, changes in water balance-related traits (bulk water, hemolymph and dehydration tolerance) showed superiority of D. simulans at 15 °C but of D. melanogaster at 25 °C growth temperature. Rate of carbohydrate utilization under desiccation stress did not differ at 15 °C in both the species. Fourth, effects of developmental plasticity on cuticular traits correspond with changes in the cuticular water loss i.e. water loss rates were higher at 25 °C as compared with 15 °C. Thus, D. simulans grown under cooler temperature was more desiccation tolerant than D. melanogaster. Finally, desiccation acclimation capacity of larvae and adults is higher for D. simulans reared at 15 °C but quite low at 25 °C. Thus, D. simulans and D. melanogaster have evolved different strategies of water conservation consistent with their adaptations to dry and wet habitats in the western Himalayas. Our results suggest that D. simulans from lowland localities seems vulnerable due to limited acclimation potential in the context of global climatic change in the western Himalayas. Finally, this is the first report on higher desiccation resistance of D. simulans due to developmental plasticity of both the cuticular traits (body melanisation and epicuticular lipid mass) when grown at 15 °C, which is consistent with its abundance in temperate regions.