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.     

Amplification dynamics of miniature inverted-repeat transposable elements and their impact on rice trait variability

Transposable elements (TEs) are a rich source of genetic variability. Among TEs, miniature inverted-repeat TEs (MITEs) are of particular interest as they are present in high copy numbers in plant genomes and are closely associated with genes. MITEs are deletion derivatives of class II transposons, and can be mobilized by the transposases encoded by the latter through a typical cut-and-paste mechanism. However, MITEs are typically present at much higher copy numbers than class II transposons. We present here an analysis of 103 109 transposon insertion polymorphisms (TIPs) in 738 Oryza sativa genomes representing the main rice population groups. We show that an important fraction of MITE insertions has been fixed in rice concomitantly with its domestication. However, another fraction of MITE insertions is present at low frequencies. We performed MITE TIP-genome-wide association studies (TIP-GWAS) to study the impact of these elements on agronomically important traits and found that these elements uncover more trait associations than single nucleotide polymorphisms (SNPs) on important phenotypes such as grain width. Finally, using SNP-GWAS and TIP-GWAS we provide evidence of the replicative amplification of MITEs.     

Transposable elements in mosquitoes

We describe the current state of knowledge about transposable elements (TEs) in different mosquito species. DNA-based elements (class II elements), non-LTR retrotransposons (class I elements), and MITEs (Miniature Inverted Repeat Transposable Elements) are found in the three genera, Anopheles, Aedes and Culex, whereas LTR retrotransposons (class I elements) are found only in Anopheles and Aedes. Mosquitoes were the first insects in which MITEs were reported; they have several LTR retrotransposons belonging to the Pao family, which is distinct from the Gypsy-Ty3 and Copia-Ty1 families. The number of TE copies shows huge variations between classes of TEs within a given species (from 1 to 1000), in sharp contrast to Drosophila, which shows only relatively minor differences in copy number between elements (from 1 to 100). The genomes of these insects therefore display major differences in the amount of TEs and therefore in their structure and global composition. We emphasize the need for more population genetic data about the activity of TEs, their distribution over chromosomes and their frequencies in natural populations of mosquitoes, to further the current attempts to develop a transgenic mosquito unable to transmit malaria that is intended to replace the natural populations.     

How do mammalian transposons induce genetic variation? A conceptual framework

In this essay, we discuss new insights into the wide‐ranging impacts of mammalian transposable elements (TE) on gene expression and function. Nearly half of each mammalian genome is comprised of these mobile, repetitive elements. While most TEs are ancient relics, certain classes can move from one chromosomal location to another even now. Indeed, striking recent data show that extensive transposition occurs not only in the germline over evolutionary time, but also in developing somatic tissues and particular human cancers. While occasional germline TE insertions may contribute to genetic variation, many other, similar TEs appear to have little or no impact on neighboring genes. However, the effects of somatic insertions on gene expression and function remain almost completely unknown. We present a conceptual framework to understand how the ages, allele frequencies, molecular structures, and especially the genomic context of mammalian TEs each can influence their various possible functional consequences. Editor's suggested further reading in BioEssays Evolution of eukaryotic genome architecture: Insights from the study of a rapidly evolving metazoan, Oikopleura dioica Abstract     

Evidence for degeneration of the Y chromosome in the dioecious plant Silene latifolia

The human Y--probably because of its nonrecombining nature--has lost 97% of its genes since X and Y chromosomes started to diverge [1, 2]. There are clear signs of degeneration in the Drosophila miranda neoY chromosome (an autosome fused to the Y chromosome), with neoY genes showing faster protein evolution [3-6], accumulation of unpreferred codons [6], more insertions of transposable elements [5, 7], and lower levels of expression [8] than neoX genes. In the many other taxa with sex chromosomes, Y degeneration has hardly been studied. In plants, many genes are expressed in pollen [9], and strong pollen selection may oppose the degeneration of plant Y chromosomes [10]. Silene latifolia is a dioecious plant with young heteromorphic sex chromosomes [11, 12]. Here we test whether the S. latifolia Y chromosome is undergoing genetic degeneration by analyzing seven sex-linked genes. S. latifolia Y-linked genes tend to evolve faster at the protein level than their X-linked homologs, and they have lower expression levels. Several Y gene introns have increased in length, with evidence for transposable-element accumulation. We detect signs of degeneration in most of the Y-linked gene sequences analyzed, similar to those of animal Y-linked and neo-Y chromosome genes.     

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.     

MITE-Hunter: a program for discovering miniature inverted-repeat transposable elements from genomic sequences

Miniature inverted-repeat transposable elements (MITEs) are a special type of Class 2 non-autonomous transposable element (TE) that are abundant in the non-coding regions of the genes of many plant and animal species. The accurate identification of MITEs has been a challenge for existing programs because they lack coding sequences and, as such, evolve very rapidly. Because of their importance to gene and genome evolution, we developed MITE-Hunter, a program pipeline that can identify MITEs as well as other small Class 2 non-autonomous TEs from genomic DNA data sets. The output of MITE-Hunter is composed of consensus TE sequences grouped into families that can be used as a library file for homology-based TE detection programs such as RepeatMasker. MITE-Hunter was evaluated by searching the rice genomic database and comparing the output with known rice TEs. It discovered most of the previously reported rice MITEs (97.6%), and found sixteen new elements. MITE-Hunter was also compared with two other MITE discovery programs, FINDMITE and MUST. Unlike MITE-Hunter, neither of these programs can search large genomic data sets including whole genome sequences. More importantly, MITE-Hunter is significantly more accurate than either FINDMITE or MUST as the vast majority of their outputs are false-positives.     

DcSto: carrot Stowaway-like elements are abundant, diverse, and polymorphic

We investigated nine families of Stowaway-like miniature inverted-repeat transposable elements (MITEs) in the carrot genome, named DcSto1 to DcSto9. All of them were AT-rich and shared a highly conserved 6 bp-long TIR typical for Stowaways. The copy number of DcSto1 elements was estimated as ca. 5,000 per diploid genome. We observed preference for clustered insertions of DcSto and other MITEs. Distribution of DcSto1 hybridization signals revealed presence of DcSto1 clusters within euchromatic regions along all chromosomes. An arrangement of eight regions encompassing DcSto insertion sites, studied in detail, was highly variable among plants representing different populations of Daucus carota. All of these insertions were polymorphic which most likely suggests a very recent mobilization of those elements. Insertions of DcSto near carrot genes and presence of putative promoters, regulatory motifs, and polyA signals within their sequences might suggest a possible involvement of DcSto in the regulation of gene expression.     

Comparative analysis of Stowaway-like miniature inverted repeat transposable elements in wheat group 7 chromosomes: Abundance,composition, and evolution

Miniature inverted repeat transposable elements (MITEs) are the most ubiquitous transposable elements in eukaryotic genomes; they play a prominent role in sequence divergence and genome evolution. There are many well-characterized Stowaway-like MITE families in wheat, but their distribution, abundance, and composition at the chromosome level are still not well understood. In this study, we systematically investigated the Stowaway-like MITEs in wheat group 7 chromosomes based on the survey sequences of isolated wheat chromosomes, to compare them at the chromosome level and to reveal their evolutionary role on wheat polyploidization. In summary, 2026 MITEs were identified, of which 587, 714, and 725 were distributed on 7A, 7B, and 7D chromosomes, respectively. There are more MITEs present on 7D, compared to 7A and 7B, suggesting A and B subgenomes eliminated some repetitive elements during two hybridization processes. Furthermore, some chromosome/arm-specific MITEs were also identified, providing information on the function and evolution of MITEs in wheat genomes. The sequence diversity of the MITE insertions was also investigated. This study for the first time investigated the abundance and composition of MITEs at the chromosome level, which will be beneficial to improve our understanding of the distribution of wheat MITEs and their evolutionary role in polyploidization.     

Genome organization and gene expression shape the transposable element distribution in the Drosophila melanogaster euchromatin

The distribution of transposable elements (TEs) in a genome reflects a balance between insertion rate and selection against new insertions. Understanding the distribution of TEs therefore provides insights into the forces shaping the organization of genomes. Past research has shown that TEs tend to accumulate in genomic regions with low gene density and low recombination rate. However, little is known about the factors modulating insertion rates across the genome and their evolutionary significance. One candidate factor is gene expression, which has been suggested to increase local insertion rate by rendering DNA more accessible. We test this hypothesis by comparing the TE density around germline- and soma-expressed genes in the euchromatin of Drosophila melanogaster. Because only insertions that occur in the germline are transmitted to the next generation, we predicted a higher density of TEs around germline-expressed genes than soma-expressed genes. We show that the rate of TE insertions is greater near germline- than soma-expressed genes. However, this effect is partly offset by stronger selection for genome compactness (against excess noncoding DNA) on germline-expressed genes. We also demonstrate that the local genome organization in clusters of coexpressed genes plays a fundamental role in the genomic distribution of TEs. Our analysis shows that—in addition to recombination rate—the distribution of TEs is shaped by the interaction of gene expression and genome organization. The important role of selection for compactness sheds a new light on the role of TEs in genome evolution. Instead of making genomes grow passively, TEs are controlled by the forces shaping genome compactness, most likely linked to the efficiency of gene expression or its complexity and possibly their interaction with mechanisms of TE silencing.     

MAK,a computational tool kit for automated MITE analysis

Miniature inverted repeat transposable elements (MITEs) are ubiquitous and numerous in higher eukaryotic genomes. Analysis of MITE families is laborious and time consuming, especially when multiple MITE families are involved in the study. Based on the structural characteristics of MITEs and genetic principles for transposable elements (TEs), we have developed a computational tool kit named MITE analysis kit (MAK) to automate the processes (http://perl.idmb.tamu.edu/mak.htm). In addition to its ability to routinely retrieve family member sequences and to report the positions of these elements relative to the closest neighboring genes, MAK is a powerful tool for revealing anchor elements that link MITE families to known transposable element families. Implementation of the MAK is described, as are genetic principles and algorithms used in its derivation. Test runs of the programs for several MITE families yielded anchor sequences that retain TIRs and coding regions reminiscent of transposases. These anchor sequences are consistent with previously reported putative autonomous elements for these MITE families. Furthermore, analysis of two MITE families with no known links to any transposon family revealed two novel transposon families, namely Math and Kid, belonging to the IS5/Harbinger/PIF superfamily.     

Dynamics of <Emphasis Type="Italic">Vulmar/VulMITE</Emphasis> group of transposable elements in <Emphasis Type="Italic">Chenopodiaceae</Emphasis> subfamily <Emphasis Type="Italic">Betoideae</Emphasis>

Transposable elements are important factors driving plant genome evolution. Upon their mobilization, novel insertion polymorphisms are being created. We investigated differences in copy number and insertion polymorphism of a group of Mariner-like transposable elements Vulmar and related VulMITE miniature inverted-repeat transposable elements (MITEs) in species representing subfamily Betoideae. Insertion sites of these elements were identified using a modified transposon display protocol, allowing amplification of longer fragments representing regions flanking insertion sites. Subsequently, a subset of TD fragments was converted into insertion site-based polymorphism (ISBP) markers. The investigated group of transposable elements was the most abundant in accessions representing the section Beta, showing intraspecific insertion polymorphisms likely resulting from their recent activity. In contrast, no unique insertions were observed for species of the genus Beta section Corollinae, while a set of section-specific insertions was observed in the genus Patellifolia, however, only two of them were polymorphic between P. procumbens and P. webbiana. We hypothesize that Vulmar and VulMITE elements were inactivated in the section Corollinae, while they remained active in the section Beta and the genus Patellifolia. The ISBP markers generally confirmed the insertion patterns observed with TD markers, including presence of distinct subsets of TE insertions specific to Beta and Patellifolia.     

The lack of recombination drives the fixation of transposable elements on the fourth chromosome of Drosophila melanogaster

In regions of suppressed recombination, where selection is expected to be less efficient in removing slightly deleterious mutations, transposable element (TE) insertions should be more likely to drift to higher frequencies, and even to reach fixation. In the absence of excision events, once a TE is fixed it cannot be eliminated from the population, and accumulation of elements thus should become an irreversible process. In the long term, this can drive the degeneration of large non-recombining fractions of the genomes. Chromosome 4 of Drosophila melanogaster has very low levels of recombination, if any, and this could be causing its degeneration. Here we report the results of a PCR-based analysis of the population frequencies of TE insertions in a sample from three African natural populations. We investigated 27 insertions from 12 TE families, located in regions of either suppressed or free recombination. Our results suggest that TE insertions tend to be fixed in the non-recombining regions, particularly on the fourth chromosome. We have also found that this involves all types of elements, and that fixed insertions are significantly shorter and more divergent from the canonical sequence than those segregating in the sample (28.1% vs 86.3% of the canonical length, and average nucleotide divergence (D(XY)) = 0.082 vs 0.008, respectively). Finally, DNA-based elements seem to show a greater tendency to reach fixation than retrotransposons. Implications of these findings for the population dynamics of TEs, and the evolutionary forces that shape the patterns of genetic variation in regions of reduced recombination, are discussed.     

Insertion Preference of Maize and Rice Miniature Inverted Repeat Transposable Elements as Revealed by the Analysis of Nested Elements

A 128-bp insertion into the maize waxy-B2 allele led to the discovery of Tourist, a family of miniature inverted repeat transposable elements (MITEs). As a special category of nonautonomous elements, MITEs are distinguished by their high copy number, small size, and close association with plant genes. In maize, some Tourist elements (named Tourist-Zm) are present as adjacent or nested insertions. To determine whether the formation of multimers is a common feature of MITEs, we performed a more thorough survey, including an estimation of the proportion of multimers, with 30.2 Mb of publicly available rice genome sequence. Among the 6600 MITEs identified, >10% were present as multimers. The proportion of multimers differs for different MITE families. For some MITE families, a high frequency of self-insertions was found. The fact that all 340 multimers are unique indicates that the multimers are not capable of further amplification.     

Large-scale discovery of insertion hotspots and preferential integration sites of human transposed elements

Throughout evolution, eukaryotic genomes have been invaded by transposable elements (TEs). Little is known about the factors leading to genomic proliferation of TEs, their preferred integration sites and the molecular mechanisms underlying their insertion. We analyzed hundreds of thousands nested TEs in the human genome, i.e. insertions of TEs into existing ones. We first discovered that most TEs insert within specific ‘hotspots’ along the targeted TE. In particular, retrotransposed Alu elements contain a non-canonical single nucleotide hotspot for insertion of other Alu sequences. We next devised a method for identification of integration sequence motifs of inserted TEs that are conserved within the targeted TEs. This method revealed novel sequences motifs characterizing insertions of various important TE families: Alu, hAT, ERV1 and MaLR. Finally, we performed a global assessment to determine the extent to which young TEs tend to nest within older transposed elements and identified a 4-fold higher tendency of TEs to insert into existing TEs than to insert within non-TE intergenic regions. Our analysis demonstrates that TEs are highly biased to insert within certain TEs, in specific orientations and within specific targeted TE positions. TE nesting events also reveal new characteristics of the molecular mechanisms underlying transposition.     

Population genetics of transposable DNA elements

This paper is an attempt to bring together the various, dispersed data published in the literature on insertion polymorphism of transposable elements from various kinds of populations (natural populations, laboratory strains, isofemale and inbred lines). Although the results deal mainly with Drosophila, data on other organisms have been incorporated when necessary to illustrate the discussion. The data pertinent to the regions of insertion, the rates of transposition and excision, the copy number regulation, and the degree of heterozygosity were analysed in order to be confronted with the speculations made with various theoretical models of population biology of transposable elements. The parameters of these models are very sensitive to the values of the transposable element characteristics estimated on populations, and according to the difficulties of these estimations (population not at equilibrium, particular mutations used to estimate the transposition and excision rates, trouble with the in situ technique used to localize the insertions, undesired mobilization of TEs in crosses, spontaneous genome resetting, environmental effects, etc.) it cannot be decided accurately which model better accounts for the population dynamics of these TEs. Tendencies, however, emerge in Drosophila: the copia element shows evidence for deficiency of insertions on the X chromosomes, a result consistent with selection against mutational effects of copia insertions; the P element repartition does not significantly deviate from the neutral assumption, in spite of a systematic copy number of insertions higher on the X than on the autosomes. Data on other elements support either the neutral model of TE containment, neither of the two models, or both. Prudence in conclusion should then be de rigueur when dealing with such kind of data. Finally the potential roles of TEs in population adaptation and evalution are discussed.     

Transposable elements are found in a large number of human protein-coding genes.

To study the genome-wide impact of transposable elements (TEs) on the evolution of protein-coding regions, we examined 13 799 human genes and found 533 (approximately 4%) cases of TEs within protein-coding regions. The majority of these TEs (approximately 89.5%) reside within 'introns' and were recruited into coding regions as novel exons. We found that TE integration often has an effect on gene function. In particular, there were two mouse genes whose coding regions consist largely of TEs, suggesting that TE insertion might create new genes. Thus, there is increasing evidence for an important role of TEs in gene evolution. Because many TEs are taxon-specific, their integration into coding regions could accelerate species divergence.     

Marker utility of miniature inverted-repeat transposable elements for wheat biodiversity and evolution

Transposable elements (TEs) account for up to 80% of the wheat genome and are considered one of the main drivers of wheat genome evolution. However, the contribution of TEs to the divergence and evolution of wheat genomes is not fully understood. In this study, we have developed 55 miniature inverted-repeat transposable element (MITE) markers that are based on the presence/absence of an element, with over 60% of these 55 MITE insertions associated with wheat genes. We then applied these markers to assess genetic diversity among Triticum and Aegilops species, including diploid (AA, BB and DD genomes), tetraploid (BBAA genome) and hexaploid (BBAADD genome) species. While 18.2% of the MITE markers showed similar insertions in all species indicating that those are fossil insertions, 81.8% of the markers showed polymorphic insertions among species, subspecies, and accessions. Furthermore, a phylogenetic analysis based on MITE markers revealed that species were clustered based on genus, genome composition, and ploidy level, while 47.13% genetic divergence was observed between the two main clusters, diploids versus polyploids. In addition, we provide evidence for MITE dynamics in wild emmer populations. The use of MITEs as evolutionary markers might shed more light on the origin of the B-genome of polyploid wheat.