A branching-process model for the evolution of transposable elements incorporating selection

We have formulated a very general mathematical model to analyze the evolution of transposable genetic elements in prokaryotic populations. Transposable genetic elements are DNA sequences able to replicate and insert copies of themselves at new locations in the genome. This work characterizes the equilibrium distribution of copy number under the influence of copy number-dependent selection, transposition and deletion. Our principal results concern the equilibrium distribution of copy number in response to various selective regimes. For particular transposition patterns (e.g. unregulated transposition or copy number-dependent transposition), equilibrium distributions are calculated numerically for a variety of specific selection patterns. Selection is quantified through specification of the expected number of offspring for individuals of each type, which is generally a non-increasing function of copy number, in accord with the usual evolutionary speculations.     

Evolutionary forces generating sequence homogeneity and heterogeneity within retrotransposon families

Genome projects allow us to sample copies of a retrotransposon sequence family residing in a host genome. The variation in DNA sequence between these individual copies will reflect the evolutionary process that has spread the sequences through the genome. Here I review quantitatively the expected diversity of elements belonging to a transposable genetic element family. I use a simple neutral model for replicative mobile DNAs such as retrotransposons to predict the extent of sequence variability between members of a single family of transposable elements, both within and between species. The effects of horizontal transfer are also explored. I also consider the impact on these distributions of an increase in transposition rate arising from a mutational change in copy of the sequence. In addition, I consider the question of the interaction between retrotransposons and their hosts, and the causes of the abundance of transposable elements in the genomes that they occupy.     

Population genetics models of competition between transposable element subfamilies

Transposable elements are one of the major components of genomes. Some copies are fully efficient; i.e., they are able to produce the proteins needed for their own transposition, and they can move and duplicate into the genome. Other copies are mutated. They may have lost their moving ability, their coding capacity, or both, thus becoming pseudogenes slowly eliminated from the genome through deletions and natural selection. Little is known about the dynamics of such mutant elements, particularly concerning their interactions with autonomous copies. To get a better understanding of the transposable elements' evolution after their initial invasion, we have designed a population genetics model of transposable elements dynamics including mutants or nonfunctional sequences. We have particularly focused on the case where these sequences are nonautonomous elements, known to be able to use the transposition machinery produced by the autonomous ones. The results show that such copies generally prevent the system from achieving a stable transposition-selection equilibrium and that nonautonomous elements can invade the system at the expense of autonomous ones. The resulting dynamics are mainly cyclic, which highlights the similarities existing between genomic selfish DNA sequences and host-parasite systems.     

Genome size as a mutation-selection-drift process

A novel method for estimating neutral rates and patterns of DNA evolution in Drosophila takes advantage of the propensity of non-LTR retrotransposable elements to create nonfunctional, transpositionally inactive copies as a product of transposition. For many LINE elements, most copies present in a genome at any one time are nonfunctional "dead-on-arrival" (DOA) copies. Because these are off-shoots of active, transpositionally competent "master" lineages, in a gene tree of a LINE element from multiple samples from related species, the DOA lineages are expected to map to the terminal branches and the active lineages to the internal branches, the primary exceptions being when the sample includes DOA copies that are allelic or orthologous. Analysis of nucleotide substitutions and other changes along the terminal branches therefore allows estimation of the fixation process in the DOA copies, which are unconstrained with respect to protein coding; and under selective neutrality, the fixation process estimates the underlying mutational pattern. We have studied the retroelement Helena in Drosophila. An unexpectedly high rate of DNA loss was observed, yielding a half-life of unconstrained DNA sequences approximately 60-fold faster in Drosophila than in mammals. The high rate of DNA loss suggests a straightforward explanation of the seeming paradox that Drosophila has many fewer pseudogenes than found in mammalian species. Differential rates of deletion in different taxa might also contribute to the celebrated C-value paradox of why some closely related organisms can have very different DNA contents. New data presented here rule out the possibility that the transposition process itself is highly mutagenic, hence the observed linear relation between number of deletions and number of nucleotide substitutions is most easily explained by the hypothesis that both types of changes accumulate in unconstrained sequences over time.     

Maintenance of transposable element copy number in natural populations of Drosophila melanogaster and D. simulans

To investigate the main forces controlling the containment of transposable elements (TE) in natural populations, we analyzed the copia, mdg1, and 412 elements in various populations of Drosophila melanogaster and D. simulans. A lower proportion of insertion sites on the X chromosome in comparison with the autosomes suggests that selection against the detrimental effects of TE insertions is the major force containing TE copies in populations of Drosophila. This selection effect hypothesis is strengthened by the absence of the negative correlation between recombination rate and TE copy number along the chromosomes, which was expected under the alternative ectopic exchange model (selection against the deleterious rearrangements promoted by recombination between TE insertions). A cline in 412 copy number in relation to latitude was observed among the natural populations of D. simulans, with very high numbers existing in some local populations (around 60 copies in a sample from Canberra, Australia). An apparent absence of selection effects in this Canberra sample and a value of transposition rate equal to 1–2 × 10^-3 whatever the population and its copy number agree with the idea of recent but temporarily drastic TE movements in local populations. The high values of transposition rate in D. simulans clearly disfavor the hypothesis that the low amount of transposable elements in this species could result from a low transposition rate. This revised version was published online in August 2006 with corrections to the Cover Date.     

The control of copy number of IS6110 in Mycobacterium tuberculosis

Insertion sequence (IS) elements are bacterial genes that are able to transpose to different locations in the genome. These elements are often used in molecular epidemiology as genetic markers that track the spread of pathogens. Transposable elements have frequently been described as "selfish DNA" because they facilitate their own transposition, causing damage when they insert into coding regions, while contributing little if anything to the bacterial host. According to this hypothesis, the expansion of copy number of insertion sequences is opposed by negative selection against high copy numbers. From an alternative point of view, we might expect IS elements to intrinsically regulate transposition within cells, thereby limiting damage to their bacterial host. Here, we report evidence that the copy number of IS6110 in Mycobacterium tuberculosis is controlled by selection against the element. We first construct 12 different models of marker change resulting from a combination of possible transposition functions and selective regimes. We then compute the Akaike Information Criterion for each model to identify the models that best explain data consisting of serial isolates of M. tuberculosis genotyped with IS6110. We find that the best performing models all include selection against the accumulation of copies. Specifically, our analysis points to the interaction of separate copies of the element causing lethal effects. We discuss the implications of these findings for genome evolution and molecular epidemiology.     

Genetic Differentiation of Transposable Elements under Mutation and Unbiased Gene Conversion

A model is developed to predict the extent of genetic differentiation in a family of transposable elements under the combined effects of genetic drift, transposition, mutation and unbiased gene conversion. The model is based on simplifying assumptions that are valid when transposition is always to new sites and copy number per site is low. In the absence of gene conversion, the degree of differentiation as measured by the probability of identity of different elements is the same as at a single locus with the same mutation rate but in a population of effective size Nc/2, where N is the population size and c is the number of copies per individual. The inclusion of unbiased gene conversion does not significantly change this result. If, as seems to be the case, families of transposable elements are relatively homogeneous, then the model implies either that mutation rates for transposable elements are much lower than at comparable single-copy loci or that some other force, such as natural selection or biased gene conversion, is at work. Transposition is a very ineffective force for homogenizing a family of transposable elements.     

Continuous exchange of sequence information between dispersed Tc1 transposons in the Caenorhabditis elegans genome

In a genome-wide analysis of the active transposons in Caenorhabditis elegans we determined the localization and sequence of all copies of each of the six active transposon families. Most copies of the most active transposons, Tc1 and Tc3, are intact but individually have a unique sequence, because of unique patterns of single-nucleotide polymorphisms. The sequence of each of the 32 Tc1 elements is invariant in the C. elegans strain N2, which has no germline transposition. However, at the same 32 Tc1 loci in strains with germline transposition, Tc1 elements can acquire the sequence of Tc1 elements elsewhere in the N2 genome or a chimeric sequence derived from two dispersed Tc1 elements. We hypothesize that during double-strand-break repair after Tc1 excision, the template for repair can switch from the Tc1 element on the sister chromatid or homologous chromosome to a Tc1 copy elsewhere in the genome. Thus, the population of active transposable elements in C. elegans is highly dynamic because of a continuous exchange of sequence information between individual copies, potentially allowing a higher evolution rate than that found in endogenous genes.     

The evolution of mobile DNAs: when will transposons create phylogenies that look as if there is a master gene?

Some families of mammalian interspersed repetitive DNA, such as the Alu SINE sequence, appear to have evolved by the serial replacement of one active sequence with another, consistent with there being a single source of transposition: the "master gene." Alternative models, in which multiple source sequences are simultaneously active, have been called "transposon models." Transposon models differ in the proportion of elements that are active and in whether inactivation occurs at the moment of transposition or later. Here we examine the predictions of various types of transposon model regarding the patterns of sequence variation expected at an equilibrium between transposition, inactivation, and deletion. Under the master gene model, all bifurcations in the true tree of elements occur in a single lineage. We show that this property will also hold approximately for transposon models in which most elements are inactive and where at least some of the inactivation events occur after transposition. Such tree shapes are therefore not conclusive evidence for a single source of transposition.     

It takes two transposons to tango: transposable-element-mediated chromosomal rearrangements

Transposable elements (TEs) promote various chromosomal rearrangements more efficiently, and often more specifically, than other cellular processes(1-3). One explanation of such events is homologous recombination between multiple copies of a TE present in a genome. Although this does occur, strong evidence from a number of TE systems in bacteria, plants and animals suggests that another mechanism - alternative transposition - induces a large proportion of TE-associated chromosomal rearrangements. This paper reviews evidence for alternative transposition from a number of unrelated but structurally similar TEs. The similarities between alternative transposition and V(D)J recombination are also discussed, as is the use of alternative transposition as a genetic tool.     

A cloned I-factor is fully functional in Drosophila melanogaster

Summary I-R hybrid dysgenesis in Drosophila melanogaster occurs in female progeny of crosses between reactive strain females and inducer strain males, and is controlled by transposable elements called I-factors. These are 5.3 kb elements that are structurally similar to mammalian LINE elements and other retroposons. We have tested the activity of an I-factor directly, by introducing it into the genome of a reactive strain, using P-element mediated transformation. It confers the complete inducer phenotype on the reactive strain, and can stimulate dysgenesis when transformed males are mated with reactive females. It has transposed in the transformed lines, and we have cloned one of the transposed copies. This is the first time that it has been possible to demonstrate that a particular retroposon is transposition proficient, and to compare donor and transposed elements. We propose a mechanism for I-factor transposition based on these results, and the coding capacity of these elements. We have been unable to detect either autonomous transposition of a complete I-factor from a plasmid injected into reactive strain embryos, or transposition of a marked I-factor when co-injected with a complete element.     

Copy Number Control of Tn5 Transposition

Transposition of Tn5 in Escherichia coli strains containing one or multiple copies of the transposable element was investigated. It was found that the overall frequency of transposition within a cell remained constant regardless of the number of copies of Tn5 present in that cell. Experiments measuring the transposition frequency of differentially marked Tn5s confirmed that the frequency of transposition of an individual Tn5 decreased proportionally with the total number of copies of the element present in a cell. The IS50R -encoded function, protein 2, which has previously been shown to be an inhibitor of transposition, is sufficient to mediate this inhibitory effect. The concentration of protein 2 in a cell appears to modulate the transposition of individual Tn5 elements in such a way that the overall transposition of Tn5 in a cell remains constant.     

A branching process model for the evolution of transposable elements

A discrete-time multitype branching process model is presented for the evolution of transposable elements in haploid populations. An individual is classified as type i if it possesses i copies of the TE, i0. The general model incorporates copy-dependent selection and transposition, and recursion relations are derived for the distribution of the number of individuals of the various types. The asymptotic relative proportions of individuals of the different types is studied in the neutral case. The behavior of this equilibrium distribution is examined for various patterns of regulated transposition and deletion.     

Mechanisms regulating the copy numbers of six LTR retrotransposons in the genome of Drosophila melanogaster

There has been debate over the mechanisms that control the copy number of transposable elements in the genome of Drosophila melanogaster. Target sites in D. melanogaster populations are occupied at low frequencies, suggesting that there is some form of selection acting against transposable elements. Three main theories have been proposed to explain how selection acts against transposable elements: insertions of a copy of a transposable element are selected against; chromosomal rearrangements caused by ectopic exchange between element copies are selected against; or the process of transposition itself is selected against. The three theories give different predictions for the pattern of transposable element insertions in the chromosomes of D. melanogaster. We analysed the abundance of six LTR (long terminal repeat) retrotransposons on the X and fourth chromosomes of multiple strains of D. melanogaster, which we compare with the predictions of each theory. The data suggest that no one theory can account for the insertion patterns of all six retrotransposons. Comparing our results with earlier work using these transposable element families, we find a significant correlation between studies in the particular model of copy number regulation supported by the proportion of elements on the X for the different transposable element families. This suggests that different retrotransposon families are regulated by different mechanisms.     

The interaction between mobile DNAs and their hosts in a fluctuating environment

The interaction between mobile DNA sequences and their hosts raises important questions in the context of hosts which reproduce clonally with only rare horizontal transmission between clones. The activity of some mobile DNAs as reversible mutators of genes raises the possibility that, in a fluctuating environment, cells may gain an advantage if they have mobile DNAs which mutate genes whose inactivation is favoured in one of the environments that the population encounters. Here we analyse a model of this process and ask what would be the optimal rate of transposition in a population whose elements are maintained by this mechanism. We also examine the impact of horizontal transfer on such a population. With movement of elements between cells, we can imagine elements with differing rates of transposition and host cells with differing rates of transposition. We find that evolution in the population of elements favours a rapid rate of transposition, and evolution of the host cells favours cells in which this rapid rate of element-dependent transposition results in an optimal rate of transposition per cell. However, when horizontal transfer rates are high, some unexpected features of the model are observed. In particular, a polymorphism between cell types (some with an optimal rate of transposition and some with no transposition at all from endogenous elements) can be stably maintained. We consider the possible biological predictions of this analysis.