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, many previous studies have used models of transposition-selection equilibrium that rely on the assumption of a constant rate of transposition. However, since TE invasions are known to happen in bursts through time, this assumption may not be reasonable in natural populations. 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 and have subsequently had their allele frequency estimated in a population sample. By conditioning on the age of an individual TE insertion (using information contained in the number of substitutions that have occurred within the TE sequence since insertion), we determine the probability distribution for 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 non-equilibrium model, we are able to explain about 80% of the variance in TE insertion allele frequencies based on age alone. Controlling both for 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.

💡 Research Summary

The paper tackles a long‑standing problem in evolutionary genomics: how to infer the selective forces acting on transposable element (TE) insertions when the underlying transposition process is highly non‑equilibrium. Classical approaches assume a constant transposition rate and a transposition‑selection balance, which is unrealistic because TE invasions often occur in short, intense bursts and host populations frequently experience demographic fluctuations. To overcome these limitations, the authors develop a “age‑of‑allele” test of neutrality that conditions on the estimated age of each TE insertion rather than on a steady‑state transposition rate.

The key methodological innovation is to infer the insertion age from the number of neutral substitutions that have accumulated within the TE sequence since insertion. Because a newly inserted TE is identical to its source copy, any divergence reflects mutations that have arisen after the insertion event. By counting these substitutions, the authors obtain a proxy for the number of generations elapsed since insertion. With this age estimate in hand, they derive the probability distribution of the insertion’s allele frequency in a sample under a neutral Wright‑Fisher model, explicitly incorporating time‑varying effective population size (Ne(t)). The resulting conditional frequency distribution takes the form of a beta‑binomial (or beta‑negative‑binomial) that reflects random genetic drift from an initial frequency of 1/Ne at the moment of insertion, but does not require a constant transposition rate.

To test the theory, the authors assembled a dataset of 190 retrotransposon insertions from two geographically distinct Drosophila melanogaster populations (North America and Africa). For each insertion they (i) aligned the TE sequence to its consensus, (ii) counted internal nucleotide differences to estimate age, and (iii) measured the insertion’s allele frequency in a population sample of several hundred individuals. Using maximum‑likelihood methods, they fitted demographic parameters (e.g., timing and severity of bottlenecks) and evaluated the fit of the age‑conditioned neutral model to the observed frequency spectrum. Remarkably, the model explained roughly 80 % of the variance in observed frequencies solely on the basis of age, demonstrating that insertion age is the dominant predictor of frequency when demographic history is accounted for.

Residual analysis revealed systematic deviations from the neutral expectation. The majority of insertions displayed lower frequencies than predicted, consistent with weak purifying (negative) selection against TE insertions—a pattern that aligns with the view that most TEs impose a fitness cost on the host. Conversely, a small subset of insertions showed higher frequencies than expected, suggesting episodes of positive selection. These positively selected insertions tended to be located near genes involved in stress response or other adaptive pathways, hinting that some TE insertions may be co‑opted as regulatory elements.

The authors discuss several methodological caveats. First, age estimation assumes that the counted substitutions are neutral; any selective constraint within the TE could bias age estimates. Second, the demographic model collapses complex histories into a limited number of parameters (e.g., a single bottleneck‑recovery scenario), which may oversimplify real population dynamics such as migration or admixture. Third, the analysis treats each TE independently, ignoring possible interactions among insertions (e.g., epistatic effects or competition for host suppression mechanisms). Future work could integrate more sophisticated coalescent simulations, incorporate selection on the TE sequence itself, and extend the framework to other classes of large insertion mutations such as gene duplications and copy‑number variants.

In summary, this study provides a robust, non‑equilibrium framework for testing neutrality of large insertion mutations by leveraging the molecular clock embedded within the inserted element. By conditioning on insertion age, the method sidesteps the unrealistic assumption of a constant transposition rate and directly incorporates demographic fluctuations. Applied to Drosophila retrotransposons, the approach quantifies the relative contributions of drift, purifying selection, and occasional adaptive selection, offering a powerful tool for dissecting the evolutionary forces shaping genome architecture across diverse taxa.