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.     

Patterns of selection against transposons inferred from the distribution of Tc1, Tc3 and Tc5 insertions in the mut-7 line of the nematode Caenorhabditis elegans

To identify the factors (selective or mutational) that affect the distribution of transposable elements (TEs) within a genome, it is necessary to compare the pattern of newly arising element insertions to the pattern of element insertions that have been fixed in a population. To do this, we analyzed the distribution of recent mutant insertions of the Tc1, Tc3, and Tc5 elements in a mut-7 background of the nematode Caenorhabditis elegans and compared it to the distribution of element insertions (presumably fixed) within the sequenced genome. Tc1 elements preferentially insert in regions with high recombination rates, whereas Tc3 and Tc5 do not. Although Tc1 and Tc3 both insert in TA dinucleotides, there is no clear relationship between the frequency of insertions and the TA dinucleotide density. There is a strong selection against TE insertions within coding regions: the probability that a TE will be fixed is at least 31 times lower in coding regions than in noncoding regions. Contrary to the prediction of theoretical models, we found that the selective pressure against TE insertions does not increase with the recombination rate. These findings indicate that the distribution of these three transposon families in the genome of C. elegans is determined essentially by just two factors: the pattern of insertions, which is a characteristic of each family, and the selection against insertions within coding regions.     

S-element insertions are associated with the evolution of the Hsp70 genes in Drosophila melanogaster

The "selfish DNA" theory postulates that transposable elements (TEs) are intragenomic parasites, and that natural selection against deleterious effects associated with their presence is the main force preventing their genomic spread in natural populations. In agreement with this model, TEs in Drosophila melanogaster populations are usually found at low frequencies in most genomic locations. Only a few cases of fixation of TE insertions have been reported, usually in regions of low recombination, where selection is expected to be less effective. Here, we report a population genetics study on the apparent fixation of an S-element in a highly recombining region in two natural populations of D. melanogaster. Three similar fragments of an S-element are inserted into the 5' regions of three members of a heat shock gene family, Hsp70 (Hsp70Aa and Hsp70Ab in polytene chromosome band 87A, and Hsp70Bb in 87C). A PCR-based analysis suggests that the insertions are fixed or at high frequencies in the entire species. A population survey of the levels of nucleotide sequence variation at the insertion site in 87C in two natural populations of D. melanogaster provided evidence for reduced levels of variation in the region, normal levels of recombination, and selection, reflected in a significant departure from neutrality of the variant frequency spectrum. This was particularly strong for the S-element inverted repeats (IRs) and suggests that these are of functional significance for the host.     

Population‐specific dynamics and selection patterns of transposable element insertions in European natural populations

Transposable elements (TEs) are ubiquitous sequences in genomes of virtually all species. While TEs have been investigated for several decades, only recently we have the opportunity to study their genome‐wide population dynamics. Most of the studies so far have been restricted either to the analysis of the insertions annotated in the reference genome or to the analysis of a limited number of populations. Taking advantage of the European Drosophila population genomics consortium (DrosEU) sequencing data set, we have identified and measured the dynamics of TEs in a large sample of European Drosophila melanogaster natural populations. We showed that the mobilome landscape is population‐specific and highly diverse depending on the TE family. In contrast with previous studies based on SNP variants, no geographical structure was observed for TE abundance or TE divergence in European populations. We further identified de novo individual insertions using two available programs and, as expected, most of the insertions were present at low frequencies. Nevertheless, we identified a subset of TEs present at high frequencies and located in genomic regions with a high recombination rate. These TEs are candidates for being the target of positive selection, although neutral processes should be discarded before reaching any conclusion on the type of selection acting on them. Finally, parallel patterns of association between the frequency of TE insertions and several geographical and temporal variables were found between European and North American populations, suggesting that TEs can be potentially implicated in the adaptation of populations across continents.     

Transposons but not retrotransposons are located preferentially in regions of high recombination rate in Caenorhabditis elegans

We analyzed the distribution of transposable elements (TEs: transposons, LTR retrotransposons, and non-LTR retrotransposons) in the chromosomes of the nematode Caenorhabditis elegans. The density of transposons (DNA-based elements) along the chromosomes was found to be positively correlated with recombination rate, but this relationship was not observed for LTR or non-LTR retrotransposons (RNA-based elements). Gene (coding region) density is higher in regions of low recombination rate. However, the lower TE density in these regions is not due to the counterselection of TE insertions within exons since the same positive correlation between TE density and recombination rate was found in noncoding regions (both in introns and intergenic DNA). These data are not compatible with a global model of selection acting against TE insertions, for which an accumulation of elements in regions of reduced recombination is expected. We also found no evidence for a stronger selection against TE insertions on the X chromosome compared to the autosomes. The difference in distribution of the DNA and RNA-based elements along the chromosomes in relation to recombination rate can be explained by differences in the transposition processes.     

DNA polymorphism in recombining and non-recombining mating-type-specific loci of the smut fungus Microbotryum

The population-genetic processes leading to the genetic degeneration of non-recombining regions have mainly been studied in animal and plant sex chromosomes. Here, we report population genetic analysis of the processes in the non-recombining mating-type-specific regions of the smut fungus Microbotryum violaceum. M. violaceum has A1 and A2 mating types, determined by mating-type-specific ‘sex chromosomes'' that contain 1–2 Mb long non-recombining regions. If genetic degeneration were occurring, then one would expect reduced DNA polymorphism in the non-recombining regions of this fungus. The analysis of DNA diversity among 19 M. violaceum strains, collected across Europe from Silene latifolia flowers, revealed that (i) DNA polymorphism is relatively low in all 20 studied loci (π∼0.15%), (ii) it is not significantly different between the two mating-type-specific chromosomes nor between the non-recombining and recombining regions, (iii) there is substantial population structure in M. violaceum populations, which resembles that of its host species, S. latifolia, and (iv) there is significant linkage disequilibrium, suggesting that widespread selfing in this species results in a reduction of the effective recombination rate across the genome. We hypothesise that selfing-related reduction of recombination across the M. violaceum genome negates the difference in the level of DNA polymorphism between the recombining and non-recombining regions, and may possibly lead to similar levels of genetic degeneration in the mating-type-specific regions of the non-recombining ‘sex chromosomes'' and elsewhere in the genome.     

The effects of recombination rate on the distribution and abundance of transposable elements

Transposable elements (TEs) often accumulate in regions of the genome with suppressed recombination. But it is unclear whether this pattern reflects a reduction in the efficacy of selection against deleterious insertions or a relaxation of ectopic recombination. Discriminating between these two hypotheses has been difficult, because no formal model has investigated the effects of recombination under the deleterious insertion model. Here we take a simulation-based approach to analyze this scenario and determine the conditions under which element accumulation is expected in low recombination regions. We show that TEs become fixed as a result of Hill-Robertson effects in the form of Muller's ratchet, but only in regions of extremely low recombination when excision is effectively absent and synergism between elements is weak. These results have important implications for differentiating between the leading models of how selection acts on TEs and should help to interpret emerging population genetic and genomic data.     

Size matters: non-LTR retrotransposable elements and ectopic recombination in Drosophila

The Drosophila melanogaster genome contains approximately 100 distinct families of transposable elements (TEs). In the euchromatic part of the genome, each family is present in a small number of copies (5-150 copies), with individual copies of TEs often present at very low frequencies in populations. This pattern is likely to reflect a balance between the inflow of TEs by transposition and the removal of TEs by natural selection. The nature of natural selection acting against TEs remains controversial. We provide evidence that selection against chromosome abnormalities caused by ectopic recombination limits the spread of some TEs. We also demonstrate for the first time that some TE families in the Drosophila euchromatin appear to be only marginally affected by purifying selection and contain many copies at high population frequencies. We argue that TEs in these families attain high population frequencies and even reach fixation as a result of low family-wide transposition rates leading to low TE copy numbers and consequently reduced strength of selection acting on individual TE copies. Fixation of TEs in these families should provide an upward pressure on the size of intergenic sequences counterbalancing rapid DNA loss through small deletions. Copy-number-dependent selection on TE families caused by ectopic recombination may also promote diversity among TEs in the Drosophila genome.     

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.     

Active miniature transposons from a plant genome and its nonrecombining Y chromosome

Mechanisms involved in eroding fitness of evolving Y chromosomes have been the focus of much theoretical and empirical work. Evolving Y chromosomes are expected to accumulate transposable elements (TEs), but it is not known whether such accumulation contributes to their genetic degeneration. Among TEs, miniature inverted-repeat transposable elements are nonautonomous DNA transposons, often inserted in introns and untranslated regions of genes. Thus, if they invade Y-linked genes and selection against their insertion is ineffective, they could contribute to genetic degeneration of evolving Y chromosomes. Here, we examine the population dynamics of active MITEs in the young Y chromosomes of the plant Silene latifolia and compare their distribution with those in recombining genomic regions. To isolate active MITEs, we developed a straightforward approach on the basis of the assumption that recent transposon insertions or excisions create singleton or low-frequency size polymorphisms that can be detected in alleles from natural populations. Transposon display was then used to infer the distribution of MITE insertion frequencies. The overall frequency spectrum showed an excess of singleton and low-frequency insertions, which suggests that these elements are readily removed from recombining chromosomes. In contrast, insertions on the Y chromosomes were present at high frequencies. Their potential contribution to Y degeneration is discussed.     

Patterns of insertion and deletion in contrasting chromatin domains

Transposable elements (TEs) play a fundamental role in the evolution of genomes. In Drosophila they are disproportionately represented in regions of low recombination, such as in heterochromatin. This pattern has been attributed to selection against repeated elements in regions of normal recombination, owing to either (1) the slightly deleterious position effects of TE insertions near or into genes, or (2) strong selection against chromosomal abnormalities arising from ectopic exchange between TE repeats. We have used defective non-long-terminal repeat (LTR) TEs that are "dead-on-arrival" (DOA) and unable to transpose in order to estimate spontaneous deletion rates in different constituents of chromatin. These elements have previously provided evidence for an extremely high rate of spontaneous deletion in Drosophila as compared with mammals, potentially explaining at least part of the differences in the genome sizes in these organisms. However, rates of deletion could be overestimated due to positive selection for a smaller likelihood of ectopic exchange. In this article, we show that rates of spontaneous deletion in DOA repeats are as high in heterochromatin and regions of euchromatin with low recombination as they are in regions of euchromatin with normal recombination. We have also examined the age distribution of five non-LTR families throughout the genome. We show that there is substantial variation in the historical pattern of transposition of these TEs. The overrepresentation of TEs in the heterochromatin is primarily due to their longer retention time in heterochromatin, as evidenced by the average time since insertion. Fragments inserted recently are much more evenly distributed in the genome. This contrast demonstrates that the accumulation of TEs in heterochromatin and in euchromatic regions of low recombination is not due to biased transposition but by greater probabilities of fixation in these regions relative to regions of normal recombination.     

Retroelements: tools for sex chromosome evolution

Many eukaryotic taxa inherit a heteromorphic sex chromosome pair. It is a generally accepted hypothesis that the sex chromosome pair is derived from a pair of homologous autosomes that has developed after the occurrence of a sex differentiator in an evolutionary process into two structurally and functionally different partners. In most of the analyzed systems the occurrence of the dominant sex differentiator is paralleled by the suppression of recombination within and close by that region. The recombinational isolation can spread in an evolutionary selection process from neighboring regions finally over the whole chromosome. Suppression of recombination strongly biases the distribution of retrotransposons in the genome. Our results and that from others indicate that the major force driving the evolution of Y chromosomes are retrotransposons, remodeling euchromatic chromosome structures into heterochromatic ones. In our model, intact or already eroded retrotransposons become trapped due to their inherent transposition mechanisms in non-recombining regions. The massive accumulation of retrotransposons interferes strongly with the activity of genes. We hypothesize that Y chromosome degeneration is a stepwise evolutionary process: (1) Massive accumulation of retrotransposons occurs in the non-recombining regions. (2) Heterochromatic nucleation centers are formed as a consequence of genomic defense against invasive parasitic elements; the established nucleation centers become epigenetically inherited. (3) Spreading of heterochromatin from the nucleation centers into flanking regions induces in an adaptive process gene silencing of neighbored genes that could either be still intact or in an already eroded condition, e.g., showing point mutations, deletions, insertions; the retroelements should be subjects to the same forces of deterioration as the genes themselves. (4) Constitutive silenced genes are not committed to the same genetic selection pressure as active genes and therefore more exposed to the decay process. (5) Gene dosage balance is reestablished by the parallel evolution of dosage compensation mechanisms. The evolving secondary sex chromosomes, neo-X and neo-Y, of Drosophila miranda are revealed to be a unique and potent model system to catch the evolutionary Y deterioration process in progress.     

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

Although transposable elements (TEs) are known to be potent sources of mutation, their contribution to the generation of recent adaptive changes has never been systematically assessed. In this work, we conduct a genome-wide screen for adaptive TE insertions in Drosophila melanogaster that have taken place during or after the spread of this species out of Africa. We determine population frequencies of 902 of the 1,572 TEs in Release 3 of the D. melanogaster genome and identify a set of 13 putatively adaptive TEs. These 13 TEs increased in population frequency sharply after the spread out of Africa. We argue that many of these TEs are in fact adaptive by demonstrating that the regions flanking five of these TEs display signatures of partial selective sweeps. Furthermore, we show that eight out of the 13 putatively adaptive elements show population frequency heterogeneity consistent with these elements playing a role in adaptation to temperate climates. We conclude that TEs have contributed considerably to recent adaptive evolution (one TE-induced adaptation every 200-1,250 y). The majority of these adaptive insertions are likely to be involved in regulatory changes. Our results also suggest that TE-induced adaptations arise more often from standing variants than from new mutations. Such a high rate of TE-induced adaptation is inconsistent with the number of fixed TEs in the D. melanogaster genome, and we discuss possible explanations for this discrepancy.     

SHORT INDELS ARE SUBJECT TO INSERTION‐BIASED GENE CONVERSION

Recombination between homologous loci is accompanied by formation of heteroduplexes. Repairing mismatches in heteroduplexes often leads to single nucleotide substitutions in a process known as gene conversion. Gene conversion was shown to be GC‐biased in different organisms; that is, a W(A or T)→S(G or C) substitution is more likely in this process than a S→W substitution. Here, we show that the insertion/deletion ratio for short noncoding indels that reach fixation between species is positively correlated with the recombination rate in Drosophila melanogaster, Homo sapiens, and Saccharomyces cerevisiae. This correlation is both due to an increase of the fixation rate of insertions and decrease of the fixation rate of deletions in regions of high recombination. Whole‐genome data on indel polymorphism and divergence in D. melanogaster rule out mutation biases and selection as the cause of this trend, pointing to insertion‐biased gene conversion as the most likely explanation. The bias toward insertions is the strongest for single‐nucleotide indels, and decreases with indel length. In regions of high recombination rate this bias leads to an up to ～5‐fold excess of fixed short insertions over deletions, and substantially affects the evolution of DNA segments.     

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.     

High Rate of Recent Transposable Element–Induced Adaptation in Drosophila melanogaster

Although transposable elements (TEs) are known to be potent sources of mutation, their contribution to the generation of recent adaptive changes has never been systematically assessed. In this work, we conduct a genome-wide screen for adaptive TE insertions in Drosophila melanogaster that have taken place during or after the spread of this species out of Africa. We determine population frequencies of 902 of the 1,572 TEs in Release 3 of the D. melanogaster genome and identify a set of 13 putatively adaptive TEs. These 13 TEs increased in population frequency sharply after the spread out of Africa. We argue that many of these TEs are in fact adaptive by demonstrating that the regions flanking five of these TEs display signatures of partial selective sweeps. Furthermore, we show that eight out of the 13 putatively adaptive elements show population frequency heterogeneity consistent with these elements playing a role in adaptation to temperate climates. We conclude that TEs have contributed considerably to recent adaptive evolution (one TE-induced adaptation every 200–1,250 y). The majority of these adaptive insertions are likely to be involved in regulatory changes. Our results also suggest that TE-induced adaptations arise more often from standing variants than from new mutations. Such a high rate of TE-induced adaptation is inconsistent with the number of fixed TEs in the D. melanogaster genome, and we discuss possible explanations for this discrepancy.     

Species and Population Level Molecular Profiling Reveals Cryptic Recombination and Emergent Asymmetry in the Dimorphic Mating Locus of C. reinhardtii

Heteromorphic sex-determining regions or mating-type loci can contain large regions of non-recombining sequence where selection operates under different constraints than in freely recombining autosomal regions. Detailed studies of these non-recombining regions can provide insights into how genes are gained and lost, and how genetic isolation is maintained between mating haplotypes or sex chromosomes. The Chlamydomonas reinhardtii mating-type locus (MT) is a complex polygenic region characterized by sequence rearrangements and suppressed recombination between its two haplotypes, MT+ and MT−. We used new sequence information to redefine the genetic contents of MT and found repeated translocations from autosomes as well as sexually controlled expression patterns for several newly identified genes. We examined sequence diversity of MT genes from wild isolates of C. reinhardtii to investigate the impacts of recombination suppression. Our population data revealed two previously unreported types of genetic exchange in Chlamydomonas MT—gene conversion in the rearranged domains, and crossover exchanges in flanking domains—both of which contribute to maintenance of genetic homogeneity between haplotypes. To investigate the cause of blocked recombination in MT we assessed recombination rates in crosses where the parents were homozygous at MT. While normal recombination was restored in MT+×MT+ crosses, it was still suppressed in MT−×MT− crosses. These data revealed an underlying asymmetry in the two MT haplotypes and suggest that sequence rearrangements are insufficient to fully account for recombination suppression. Together our findings reveal new evolutionary dynamics for mating loci and have implications for the evolution of heteromorphic sex chromosomes and other non-recombining genomic regions.     

Transposable Element Numbers in Cosmopolitan Inversions from a Natural Population of Drosophila Melanogaster

Population studies of the distribution of transposable elements (TEs) on the chromosomes of Drosophila melanogaster have suggested that their copy number increase due to transposition is balanced by some form of natural selection. Theory suggests that, as a consequence of deleterious ectopic meiotic exchange between TEs, selection can favor genomes with lower TE copy numbers. This predicts that TEs should be less deleterious, and hence more abundant, in chromosomal regions in which recombination is reduced. To test this, we surveyed the abundance and locations of 10 families of TEs in recombination-suppressing chromosomal inversions from a natural population. The sample of 49 chromosomes included multiple independent isolates of seven different inversions and a corresponding set of standard chromosomes. For all 10 TE families pooled, copy numbers were significantly higher overall within low frequency inversions than within corresponding regions of standard chromosomes. TEs occupied chromosomal sites at significantly higher frequencies within the In(3R)M0 and In(3R)K inversions than within the corresponding regions of standard 3R chromosomes. These results are consistent with the predictions of the ectopic exchange model.     

Fitness costs of Doc expression are insufficient to stabilize its copy number in Drosophila melanogaster

The stable coexistence of transposable elements (TEs) with their host genome over long periods of time suggests TEs have to impose some deleterious effect upon their host fitness. Three mechanisms have been proposed to account for the deleterious effect caused by TEs: host gene interruptions by TE insertions, chromosomal rearrangements by TE-induced ectopic recombination, and costly TE expression. However, the relative importance of these mechanisms remains controversial. Here, we test specifically if TE expression accounts for the host fitness cost imposed by TE insertions. In the retrotransposon Doc, expression requires binding of the host RNA polymerase to the internal promoter. If expression of Doc elements is deleterious to their host, Doc copies with promoters would be more strongly selected against and would persist in the population for shorter periods of time compared with Docs lacking promoters. We tested this prediction using sequence-specific amplified polymorphism (SSAP) analyses. We compared the populations of these two types of Doc elements in two sets of lines of Drosophila melanogaster: selection-free isogenic lines accumulating new Doc insertions and isogenized isofemale lines sampled from a natural population. We found that (1) there is no difference in the proportion of promoter-bearing and promoter-lacking copies between sets of lines, and (2) the site occupancy distribution of promoter-bearing copies does not skew toward lower frequency compared with that of promoter-lacking copies. Thus, selection against promoter-bearing copies does not appear to be stronger than that of promoter-lacking copies. Our results show that expression is not playing a major role in stabilizing Doc copy numbers.     

Paleogenomic record of the extinction of human endogenous retrovirus ERV9

An outstanding question of genome evolution is what stops the invasion of a host genome by transposable elements (TEs). The human genome, harboring the remnants of many extinct TE families, offers an extraordinary opportunity to investigate this problem. ERV9 is an endogenous retrovirus repeatedly mobilized during primate evolution, 15 to 6 million years ago (MYA), which left a trace of over a hundred provirus-like copies and at least 4,000 solitary long terminal repeats (LTRs) in the human genome. Then, its proliferation ceased for unknown reasons, and the family went extinct. We have made a detailed reconstruction of its last active subfamily, ERV9_XII, by examining 115 solitary LTRs from it. These insertions were grouped into 11 sets according to shared nucleotide variants, which could be placed in a sequential order of 10 to 6 MYA. At least 75% of the subfamily was produced 8 to 6 MYA, during a stage of intense proliferation. With new analytical tools, we show that the youngest and most prolific sets may have been produced by effectively instantaneous expansions of corresponding single-sequence variants. The extinction of this family apparently was not a consequence of its slow gradual degeneration, but the outcome of the fixation of specific restrictive alleles in the human-chimpanzee ancestral population. Three species-specific insertions (two in humans and one in chimpanzees) were identified, further supporting that extinction took place when these two species were beginning to diverge. These are the only fixed differences of this kind so far observed between humans and chimpanzees, apart from those belonging to the human endogenous retrovirus K family.