Quantitative genetic analysis of copia retrotransposon activity in inbred Drosophila melanogaster lines

Genetic variation in P-element dysgenic sterility is associated with double-strand break repair and alternative splicing of TE transcripts

The germline mobilization of transposable elements (TEs) by small RNA mediated silencing pathways is conserved across eukaryotes and critical for ensuring the integrity of gamete genomes. However, genomes are recurrently invaded by novel TEs through horizontal transfer. These invading TEs are not targeted by host small RNAs, and their unregulated activity can cause DNA damage in germline cells and ultimately lead to sterility. Here we use hybrid dysgenesis-a sterility syndrome of Drosophila caused by transposition of invading P-element DNA transposons-to uncover host genetic variants that modulate dysgenic sterility. Using a panel of highly recombinant inbred lines of Drosophila melanogaster, we identified two linked quantitative trait loci (QTL) that determine the severity of dysgenic sterility in young and old females, respectively. We show that ovaries of fertile genotypes exhibit increased expression of splicing factors that suppress the production of transposase encoding transcripts, which likely reduces the transposition rate and associated DNA damage. We also show that fertile alleles are associated with decreased sensitivity to double-stranded breaks and enhanced DNA repair, explaining their ability to withstand high germline transposition rates. Together, our work reveals a diversity of mechanisms whereby host genotype modulates the cost of an invading TE, and points to genetic variants that were likely beneficial during the P-element invasion.

The transposition rate has little influence on the plateauing level of the P-element

The popular trap model assumes that the invasions of transposable elements (TEs) in mammals and invertebrates are stopped by piRNAs that emerge after insertion of the TE into a piRNA cluster. It remains, however, still unclear which factors influence the dynamics of TE invasions. The activity of the TE (i.e., transposition rate) is one frequently discussed key factor. Here we take advantage of the temperature-dependent activity of the P-element, a widely studied eukaryotic TE, to test how TE activity affects the dynamics of a TE invasion. We monitored P-element invasion dynamics in experimental Drosophila simulans populations at hot and cold culture conditions. Despite marked differences in transposition rates, the P-element reached very similar copy numbers at both temperatures. The reduction of the insertion rate upon approaching the copy number plateau was accompanied by similar amounts of piRNAs against the P-element at both temperatures. Nevertheless, we also observed fewer P-element insertions in piRNA clusters than expected, which is not compatible with a simple trap model. The ping-pong cycle, which degrades TE transcripts, becomes typically active after the copy number plateaued. We generated a model, with few parameters, that largely captures the observed invasion dynamics. We conclude that the transposition rate has at the most only a minor influence on TE abundance, but other factors, such as paramutations or selection against TE insertions are shaping the TE composition.

Regulation of transposable elements: Interplay between TE-encoded regulatory sequences and host-specific trans-acting factors in Drosophila melanogaster

Transposable elements (TEs) are mobile genetic elements that can move around the genome, and their expression is one precondition for this mobility. Because the insertion of TEs in new genomic positions is largely deleterious, the molecular mechanisms for transcriptional suppression have been extensively studied. In contrast, very little is known about their primary transcriptional regulation. Here, we characterize the expression dynamics of TE families in Drosophila melanogaster across a broad temperature range (13-29°C). In 71% of the expressed TE families, the expression is modulated by temperature. We show that this temperature-dependent regulation is specific for TE families and strongly affected by the genetic background. We deduce that TEs carry family-specific regulatory sequences, which are targeted by host-specific trans-acting factors, such as transcription factors. Consistent with the widespread dominant inheritance of gene expression, we also find the prevailing dominance of TE family expression. We conclude that TE family expression across a range of temperatures is regulated by an interaction between TE family-specific regulatory elements and trans-acting factors of the host.

Tempo and Mode of Transposable Element Activity in Drosophila

The evolutionary dynamics of transposable element (TE) insertions have been of continued interest since TE activity has important implications for genome evolution and adaptation. Here, we infer the transposition dynamics of TEs by comparing their abundance in natural D. melanogaster and D. simulans populations. Sequencing pools of more than 550 South African flies to at least 320-fold coverage, we determined the genome wide TE insertion frequencies in both species. We suggest that the predominance of low frequency insertions in the two species (>80% of the insertions have a frequency <0.2) is probably due to a high activity of more than 58 families in both species. We provide evidence for 50% of the TE families having temporally heterogenous transposition rates with different TE families being affected in the two species. While in D. melanogaster retrotransposons were more active, DNA transposons showed higher activity levels in D. simulans. Moreover, we suggest that LTR insertions are mostly of recent origin in both species, while DNA and non-LTR insertions are older and more frequently vertically transmitted since the split of D. melanogaster and D. simulans. We propose that the high TE activity is of recent origin in both species and a consequence of the demographic history, with habitat expansion triggering a period of rapid evolution.

Co-evolution between transposable elements and their hosts: a major factor in genome size evolution?

Most models of genome size evolution emphasize changes in relative rates of and/or the efficacy of selection on insertions and deletions. However, transposable elements (TEs) are a major contributor to genome size evolution, and since they experience their own selective pressures for expansion, genome size changes may in part be driven by the dynamics of co-evolution between TEs and their hosts. Under this perspective, predictions about the conditions that allow for genome expansion may be altered. In this review, we outline the evidence for TE-host co-evolution, discuss the conditions under which these dynamics can change, and explore the possible contribution to the evolution of genome size. Aided partly by advances in our understanding of the mechanisms of TE silencing via small RNAs, there is growing evidence that the evolution of transposition rates can be important in driving genome expansion and contraction. Shifts in genome size and transposon abundance associated with interspecific hybridization and changes in mating system are consistent with an important role for transposition rate evolution, although other possible explanations persist. More understanding of the potential for the breakdown of host silencing mechanisms and/or the potential for TEs to evade host immune responses will improve our understanding of the importance of changes in TE activity in driving genome size evolution.

Genetic variation of copia suppression in Drosophila melanogaster

Transposable elements (TEs) are genomic parasites that propagate by exploiting its host reproductive machinery. However, some hosts have evolved the ability to silence TE activity, whereas others have not. We are investigating the population dynamics of TE host-silencing pathways, particularly copia long terminal repeat retrotransposon in Drosophila melanogaster. Here, we identify large effect genes involved in copia suppression by using a semi-quantitative analysis to assay levels of copia plasmids (believed to be an intermediate of transposition) in 98 recombinant inbred lines constructed from a line exhibiting high copia transpositions and a line exhibiting no transpositions. The results revealed that the influence of copia copy number and transcription level on copia plasmid concentrations are weak and that genomic factors, presumably encoded by the host, have stronger effects on transposition rates. We mapped a QTL affecting copia plasmid concentration within the 33A-43E cytological region of the second chromosome and applied a quantitative deficiency complementation analysis on this chromosomal region. One out of the two large effect deficiencies on copia plasmid concentrations corresponded to the vasa gene, an important component of the nuage-piwi RNA TE-silencing machinery. We hypothesize that copia suppression occurs by the joint action of several post-transcriptional mechanisms with at least one of the blocks taking place in the nuage.

Evolutionary dynamics of transposable elements in a small RNA world

Transposable elements (TEs) are selfish elements that cause harmful mutations, contribute to the structure of regulatory networks and shape the architecture of genomes. Natural selection against their harmful effects has long been considered the dominant force limiting their spread. It is now clear that a genome defense system of RNA-mediated silencing also plays a crucial role in limiting TE proliferation. A full understanding of TE evolutionary dynamics must consider how these forces jointly determine their proliferation within genomes. Here I consider these forces from two perspectives - dynamics within populations and evolutionary games within the germline. The analysis of TE dynamics from these two perspectives promises to provide new insight into their role in evolution.

Inferring the strength of selection in Drosophila under complex demographic models

Transposable elements (TEs) constitute a substantial fraction of the genomes of many species, and it is thus important to understand their population dynamics. The strength of natural selection against TEs is a key parameter in understanding these dynamics. In principle, the strength of selection can be inferred from the frequencies of a sample of TEs. However, complicated demographic histories, such as found in Drosophila melanogaster, could lead to a substantial distortion of the TE frequency distribution compared with that expected for a panmictic, constant-sized population. The current methodology for the estimation of selection intensity acting against TEs does not take into account demographic history and might generate erroneous estimates especially for TE families under weak selection. Here, we develop a flexible maximum likelihood methodology that explicitly accounts both for demographic history and for the ascertainment biases of identifying TEs. We apply this method to the newly generated frequency data of the BS family of non-long terminal repeat retrotransposons in D. melanogaster in concert with two recent models of the demographic history of the species to infer the intensity of selection against this family. We find the estimate to differ substantially compared with a prior estimate that was made assuming a model of constant population size. Further, we find there to be relatively little information about selection intensity present in the derived non-African frequency data and that the ancestral African subpopulation is much more informative in this respect. These findings highlight the importance of accounting for demographic history and bear on study design for the inference of selection coefficients generally.

The dynamics of the roo transposable element in mutation-accumulation lines and segregating populations of Drosophila melanogaster

We estimated the number of copies for the long terminal repeat (LTR) retrotransposable element roo in a set of long-standing Drosophila melanogaster mutation-accumulation full-sib lines and in two large laboratory populations maintained with effective population size approximately 500, all of them derived from the same isogenic origin. Estimates were based on real-time quantitative PCR and in situ hybridization. Considering previous estimates of roo copy numbers obtained at earlier stages of the experiment, the results imply a strong acceleration of the insertion rate in the accumulation lines. The detected acceleration is consistent with a model where only one (maybe a few) of the approximately 70 roo copies in the ancestral isogenic genome was active and each active copy caused new insertions with a relatively high rate ( approximately 10(-2)), with new inserts being active copies themselves. In the two laboratory populations, however, a stabilized copy number or no accelerated insertion was found. Our estimate of the average deleterious viability effects per accumulated insert [E(s) < 0.003] is too small to account for the latter finding, and we discuss the mechanisms that could contain copy number.

Direct determination of the effects of genotype and extreme temperature on the transposition of roo in long-term mutation accumulation lines of Drosophila melanogaster

Transposable elements (TEs) are mobile repetitive DNA sequences that constitute a structurally dynamic component of genomes. In order to understand the dynamics of TEs it is necessary to have information about the control of transposition and its dependence of environmental factors. After a great deal of previous work on transposition conducted on long-term mutation accumulation (MA) lines of Drosophila melanogaster started in 1987, only roo out of 16 families was found active in this genotype. Here we test the effect of the modification of the genetic background by introducing a Cy chromosome, and the effect of extreme temperature (28 degrees C) on the transposition rate of roo. Thermal stress did not affect the transposition rate, whereas the presence of a Cy chromosome in heterozygosis lowered it. There was an excess of insertions in the X chromosome, with respect to autosomes, and in the proximal and distal regions of chromosome arms that can be interpreted as target site preference. One of the control lines became highly unstable with mean insertion and excision rates of 3.0 x 10(-3) and 8.5 x 10(-4), respectively. Instability arose spontaneously during generations of mutation accumulation, and can be attributed to "de novo" mutation. Transposition in the unstable line could be directly studied on the progeny of individual males and females, from where we deduced that transposition occurs mainly, if not exclusively, in males, with a rate of 1.125 insertions per gamete. In situ hybridization with an LTR probe showed that most excisions (12 out of 14) were precise. Our data show the prominent role of genotype in transposition control and can explain rapid turnovers in the genome without increasing the number of copies.

The fate of transposable elements in asexual populations

Sexual reproduction and recombination are important for maintaining a stable copy number of transposable elements (TEs). In sexual populations, elements can be contained by purifying selection against host carriers with higher element copy numbers; however, in the absence of sex and recombination, asexual populations could be driven to extinction by an unchecked proliferation of TEs. Here we provide a theoretical framework for analyzing TE dynamics under asexual reproduction. Analytic results show that, in an infinite asexual population, an equilibrium in copy number is achieved if no element excision is possible, but that all TEs are eliminated if there is some excision. In a finite population, computer simulations demonstrate that small populations are driven to extinction by a Muller's ratchet-like process of element accumulation, but that large populations can be cured of vertically transmitted TEs, even with excision rates well below transposition rates. These results may have important consequences for newly arisen asexual lineages and may account for the lack of deleterious retrotransposons in the putatively ancient asexual bdelloid rotifers.

A rapid decrease in number of the complete ninja element and concomitant increase of the defective element in a strain of Drosophila simulans

The ninja element, originally isolated from an unstable white mutant strain white-milky (w(mky)) of Drosophila simulans, is a member of the retrotransposon family with long terminal repeats (LTRs). We show that ninja is present in high copy numbers in the w(mky)-derivative sublines white-chocolate (w(cho)) and white-persimmonl (w(psm1)), in a low copy number in another derivative subline white-milky 3 (w(mky3)), and in only a few copies in a wild type strain. We have cloned the ninja elements from these sublines and examined their structures. Most of the elements cloned (38 out of 41 independent clones) from w(cho) were full length. In contrast, only 9 of 23 independent clones from w(mky3) were full length. We hypothesize that ninja elements were integrated and lost frequently in the w(mky) strain and its derivative genomes, and that a rapid decrease in numbers of the ninja element was caused not by an increased rate of loss but by a reduction of integration of full length ninja elements in w(mky3). Each defective element had a unique deletion and/or an insertion except for the three from w(mky3), which had exactly the same 81-bp deletion in each of the 5' and 3' LTRs. The 5' and 3' ends of the deletion appeared to represent sequences similar to those of Drosophila consensussplicing sites. Ectopic splicing may have produced these defective ninja elements.

Estimates of the genomic mutation rate for detrimental alleles in Drosophila melanogaster

The net rate of mutation to deleterious but nonlethal alleles and the sizes of effects of these mutations are of great significance for many evolutionary questions. Here we describe three replicate experiments in which mutations have been accumulated on chromosome 3 of Drosophila melanogaster by means of single-male backcrosses of heterozygotes for a wild-type third chromosome. Egg-to-adult viability was assayed for nonlethal homozygous chromosomes. The rates of decline in mean and increase in variance (DM and DV, respectively) were estimated. Scaled up to the diploid whole genome, the mean DM for homozygous detrimental mutations over the three experiments was between 0.8 and 1.8%. The corresponding DV estimate was approximately 0.11%. Overall, the results suggest a lower bound estimate of at least 12% for the diploid per genome mutation rate for detrimentals. The upper bound estimates for the mean selection coefficient were between 2 and 10%, depending on the method used. Mutations with selection coefficients of at least a few percent must be the major contributors to the effects detected here and are likely to be caused mostly by transposable element insertions or indels.

Characterization of active R2 retrotransposition in the rDNA locus of Drosophila simulans

The rRNA gene (rDNA) loci of all arthropod lineages contain non-LTR retrotransposable elements that have evolved to specifically insert into the 28S rRNA genes. Extensive in vitro experiments have been conducted to investigate the mechanism of R2 retrotransposition but little is known of the insertion frequency or cellular factors that might regulate R2 activity. In this article, isofemale lines obtained from a population of Drosophila simulans were surveyed for recent R2 insertions. Within most lines, all individuals showed the same collection of R2 insertions, providing no evidence for recent R2 activity. However, in a few of the isofemale lines, virtually all individuals differed in their R2 insertion profiles. The descendants of individual pairs of flies from these "active lines" rapidly accumulated new insertions. The frequent insertion of new R2 elements was associated with the elimination of old R2 elements from the rDNA locus. The existence of lines in which R2 retrotransposes frequently and lines in which the elements appear dormant suggests that cellular mechanisms that can regulate the activity of R2 exist. Retrotransposition activity was correlated with the number of full-length R2 elements but not with the size of the rDNA locus or the number of uninserted units.

Dominance of mutations affecting viability in Drosophila melanogaster

There have been several attempts to estimate the average dominance (ratio of heterozygous to homozygous effects) of spontaneous deleterious mutations in Drosophila melanogaster, but these have given inconsistent results. We investigated whether transposable element (TE) insertions have higher average dominance for egg-to-adult viability than do point mutations, a possibility suggested by the types of fitness-depressing effects that TEs are believed to have. If so, then variation in dominance estimates among strains and crosses would be expected as a consequence of variation in TE activity. As a first test, we estimated the average dominance of all mutations and of copia insertions in a set of lines that had accumulated spontaneous mutations for 33 generations. A traditional regression method gave a dominance estimate for all mutations of 0.17, whereas average dominance of copia insertions was 0.51; the difference between these two estimates approached significance (P = 0.08). As a second test, we reanalyzed Ohnishi 1974 data on dominance of spontaneous and EMS-induced mutations. Because a considerable fraction of spontaneous mutations are caused by TE insertions, whereas EMS induces mainly point mutations, we predicted that average dominance would decline with increasing EMS concentration. This pattern was observed, but again fell short of formal significance (P = 0.07). Taken together, however, the two results give modest support for the hypothesis that TE insertions have greater average dominance in their viability effects than do point mutations, possibly as a result of deleterious effects of expression of TE-encoded genes.

Female age and sperm competition: last-male precedence declines as female age increases

Until very recently, most studies of sperm competition have focused on variation in male competitive ability. However, we now know that a number of reproductive traits, including oviposition rate, use of stored sperm and receptivity to mating, vary with female condition. Because females can play an active part in the movement of sperm within their reproductive tract, sperm competition may be influenced by female condition. Existing studies of sperm competition in fruitflies ignore the effects of female condition, using females that are 3-4 days old and in their reproductive prime. But condition will decline as a female senesces. Here, we examine the effect of female age on the outcome of sperm competition in three strains of the fruitfly, Drosophila melanogaster. Previous studies have shown that female age influences preference for mates and male ejaculation strategies. In this study, we find that when males are mated to females that are older than 17 days, last-male sperm precedence decreases significantly. These results could lead to a greater understanding of the physiological mechanisms that regulate the outcome of sperm competition.

Tissue-specificity of 412 retrotransposon expression in Drosophila simulans and D. melanogaster

We analyse the expression of the retrotransposon 412 in the soma, testes, and ovaries in populations of Drosophila simulans and D. melanogaster, using RT-PCR and in situ hybridization. We find that expression of 412 is highly variable in the soma, confirming previous findings based on Northern blots. No 412RNA is detected in the ovaries by either in situ hybridization or RT-PCR, in any population of either species. Transcripts are, however, detected in the male germline, which show a very characteristic spatial pattern of 412 expression in primary spermatocytes. There is no relationship between expression of the 412 element in the soma and in the testes in the populations. These findings show that the expression of 412 is independently regulated in the soma and the testes, and this raises the question of the real influence of the somatic transcripts on the organism and on the transposition rate.

Genomics in pure and applied entomology

Genomics is the study of the structure and function of the genome: the set of genetic information encoded in the DNA of the nucleus and organelles of an organism. It is a dynamic field that combines traditional paths of inquiry with new approaches that would have been impossible without recent technological developments. Much of the recent focus has been on obtaining the sequence of entire genomes, determining the order and organization of the genes, and developing libraries that provide immediate physical access to any desired DNA fragment. This has enabled functional studies on a genome-wide level, including analysis of the genetic basis of complex traits, quantification of global patterns of gene expression, and systematic gene disruption projects. The successful contribution of genomics to problems in applied entomology requires the cooperation of the private and public sectors to build upon the knowledge derived from the Drosophila genome and effectively develop models for other insect Orders.