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.     

TESD: a transposable element dynamics simulation environment

Various mathematical models have been used to explore the dynamics of transposable elements (TEs) within their host genomes. However, numerous factors can influence their dynamics, and we know only little about the dynamics of TEs when they first began to invade populations. In addition, the influence of population structuring has only recently been investigated. Transposable Element Simulator Dynamics, a population genomics simulation environment, has therefore been developed to provide a simple tool for analyzing the dynamics of TEs in a community based on (i) various TE parameters, such as the transposition and excision rates, the recombination rate and the coefficient of selection against TE insertions; and (ii) population parameters, such as population size and migration rates. The simulations can be used to illustrate the dynamic fate of TEs in structured populations, can be extended by using more specific molecular or demographic models, and can be useful for teaching population genetics and genomics. AVAILABILITY: TESD is distributed under GPL from the P?le Bioinformatique Lyonnais (PBIL) web server at http://pbil.univ-lyon1.fr/software/TESD     

Population and evolutionary dynamics of Helitron transposable elements in Arabidopsis thaliana

Helitrons, a recently discovered superfamily of DNA transposons that capture host gene fragments, constitute up to 2% of the Arabidopsis thaliana genome. In this study, we identified 565 insertions of a family of nonautonomous Helitrons, known as Basho elements. We aligned subsets of these elements, estimated their phylogenetic relationships, and used branch lengths to yield insight into the age of each Basho insertion. The age distribution suggests that 87% of Bashos inserted within 5 Myr, subsequent to the divergence between A. thaliana and its sister species Arabidopsis lyrata. We screened 278 of these insertions for their presence or absence in a sample of 47 A. thaliana accessions. With both phylogenetic and population frequency data, we investigated the effects of gene density, recombination rate, and element length on Basho persistence. Our analyses suggested that longer Basho copies are less likely to persist in the genome, consistent with selection against the deleterious effects of ectopic recombination between Basho elements. Furthermore, we determined that 39% of Basho elements contain fragments of expressed protein-coding genes, but all of these fragments were explained by only 5 gene-capture events. Overall, the picture of A. thaliana Helitron evolution is one of rapid expansion, relatively few gene-capture events, and weak selection correlated with element length.     

Rapid evolution at the Drosophila telomere: transposable element dynamics at an intrinsically unstable locus

Drosophila telomeres have been maintained by three families of active transposable elements (TEs), HeT-A, TAHRE, and TART, collectively referred to as HTTs, for tens of millions of years, which contrasts with an unusually high degree of HTT interspecific variation. While the impacts of conflict and domestication are often invoked to explain HTT variation, the telomeres are unstable structures such that neutral mutational processes and evolutionary tradeoffs may also drive HTT evolution. We leveraged population genomic data to analyze nearly 10,000 HTT insertions in 85  Drosophila melanogaster genomes and compared their variation to other more typical TE families. We observe that occasional large-scale copy number expansions of both HTTs and other TE families occur, highlighting that the HTTs are, like their feral cousins, typically repressed but primed to take over given the opportunity. However, large expansions of HTTs are not caused by the runaway activity of any particular HTT subfamilies or even associated with telomere-specific TE activity, as might be expected if HTTs are in strong genetic conflict with their hosts. Rather than conflict, we instead suggest that distinctive aspects of HTT copy number variation and sequence diversity largely reflect telomere instability, with HTT insertions being lost at much higher rates than other TEs elsewhere in the genome. We extend previous observations that telomere deletions occur at a high rate, and surprisingly discover that more than one-third do not appear to have been healed with an HTT insertion. We also report that some HTT families may be preferentially activated by the erosion of whole telomeres, implying the existence of HTT-specific host control mechanisms. We further suggest that the persistent telomere localization of HTTs may reflect a highly successful evolutionary strategy that trades away a stable insertion site in order to have reduced impact on the host genome. We propose that HTT evolution is driven by multiple processes, with niche specialization and telomere instability being previously underappreciated and likely predominant.     

A tiger mouse and relatives. Variants caused by an activated transposable element?

In a laboratory-bred population of wild Peruvian house mice, one male had an excessive rate of non-pairing of the X and Y chromosomes. After crossing him with laboratory stock mice, a mouse of very unusual phenotype appeared from a yellow (AyA) mother. He was yellow with black dorsal stripes; hence Tiger. He was mated to many females, and inbred F2 and F3 generations were raised. There were no more tiger phenotypes, but his F1 contained an excess of black-and-tans over yellows, showing him to be a gonosomic mosaic Ayat/atat; the homozygous cell line probably arose from the heterozygous one. The mitotic karyotype was normal. Some of Tiger's mates were of known allozyme types and their progeny were scored. The allozyme segregations were normal, except at the Es-3 locus (esterase-3), for which Tiger was typed as homozygous. Several unusual events among Tiger's close relatives included a mutation to an unstable pattern mutant, three probable translocations, and several cases of somatic defect. All unusual mice derived from Tiger's yellow mother, whose genome was one-quarter Peruvian. Yellow is associated with an ecotropic murine leukemia virus. The Peru genome is characterized by a high occurrence of mutation and aberrant karyotypes. It is suggested that something from the Peru genome in Tiger's mother caused instability of the DNA sequence associated with yellow, with related disturbance at different locations thereafter. The nature of this instability, and of the Peru genome, is discussed.     

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.     

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.     

Population dynamics of miniature inverted-repeat transposable elements (MITEs) in Medicago truncatula

Miniature inverted-repeat transposable elements (MITEs) are small and high copy number transposons, related to and mobilized by some class II autonomous elements. New MITE families can be identified by computer-based mining of sequenced genomes. We describe four MITE families related to MtPH transposons mined de novo in the genome of Medicago truncatula, together with one previously described family MITRAV. Different levels of their intra-family sequence diversity and insertion polymorphism indicate that they were active at different evolutionary periods. MetMIT1 and MITRAV families were uniform in sequence and produced highly polymorphic insertion sites in 26 ecotypes representing a M. truncatula core collection. A subset of insertions was present only in the reference genome of A17 ‘Jemalong’, suggesting that the two families might have been active in the course of domestication. In contrast, all investigated insertions of the MetMIT2 family were fixed, showing that it was not active after M. truncatula speciation. MetMIT1 elements were divided into three clusters, i.e. (I) relatively heterogenous copies fixed in the genome of M. truncatula, (II) uniform but also mostly fixed, and (III) uniform and polymorphic among the investigated accessions. It might reflect the evolutionary history of the MetMIT1 family, showing multiple bursts of activity. A number of MetMIT1 and MITRAV insertions were present within 1 kb upstream or downstream the ORF. A high proportion of insertions proximal to coding regions was unique to A17 ‘Jemalong’.     

Population dynamics of nitrifying bacteria in an aquatic ecosystem.

In a space environment such as Space Shuttle or Space Station, animal experiments with aquatic species in a closed system pose a crucial problem in maintaining their water quality for a long term. In nature, ammonia as an animal wastes is converted by nitrifying bacteria to nitrite or nitrate compounds, which usually become nitrogen sources for plants. Thus an application of the biological reactor with such bacteria attached on some filters has been suggested and experimentally studied for efficient waste managements of ammonia. Although some successful results were reported (Kozu et al. 1995, Nagaoka et al. 1998, Nakamura et al. 1997, 1998) in the space applications, purely empirical approaches have so far been taken to develop a biological filter having a stable nitrifying activity. In this study, we constructed a mathematical model to deal with the dynamics of the ammonia nitrifying processes in a biological reactor. The model describes population dynamics of the ammonia-oxidizing bacteria and the nitrite-oxidizing bacteria cultivated on the same filter. We estimated parameters involved in the model using the experimental data. The result shows that these estimated parameters could be applied to general cases and that the two bacteria are in a symbiotic relationship; they can better perform when both coexist, as has been empirically recognized. Based on the model analysis, we discuss how to prepare a high performance biological filter.     

Trans-activation of an artificial dTam3 transposable element in transgenic tobacco plants

In Antirrhinum majus only autonomous Tam3 transposons have been characterized. We investigated whether an artificial dTam3 element, with a deletion in the presumptive transposase coding region, can be trans-activated in tobacco by an activator Tam3 element, which was immobilized by the deletion of one inverted repeat. A phenotypic assay based on restored hygromycin resistance demonstrates that a dTam3 element harbouring a bacterial plasmid can be trans-activated with a low frequency. Molecular analysis confirms that the dTam3 element has been excised from the HPTII marker gene. Reintegration of the dTam3 element into the tobacco genome is detected only in one out of six hygromycin-resistant plants analysed. PCR analysis of empty donor sites shows that excision of the dTam3 element in tobacco results in rearrangements (deletions and additions), that have been shown to be characteristic of Tam3 excision in the original host Antirrhinum majus. This trans-activation assay allowed us to establish that, in contrast to what has been detected in Antirrhinum majus, a periodical temperature shift down to 15°C does not enhance dTam3 transposition in regenerating tobacco calli.     

The pogo transposable element family of Drosophila melanogaster.

Summary A 190 by insertion is associated with the white-eosin mutation in Drosophila melanogaster. This insertion is a member of a family of transposable elements, pogo elements, which is of the same class as the P and hobo elements of D. melanogaster. Strains typically have many copies of a 190 by element, 10–15 elements 1.1–1.5 kb in size and several copies of a 2.1 kb element. The smaller elements all appear to be derived from the largest by single internal deletions so that all elements share terminal sequences. They either always insert at the dinucleotide TA and have perfect 21 bp terminal inverse repeats, or have 22 by inverse repeats and produce no duplication upon insertion. Analysis by DNA blotting of their distribution and occupancy of insertion sites in different strains suggests that they may be less mobile than P or hobo. The DNA sequence of the largest element has two long open reading frames on one strand which are joined by splicing as indicated by cDNA analysis. RNAs of this strand are made, whose sizes are similar to the major size classes of elements. A protein predicted by the DNA sequence has significant homology with a human centrosomal-associated protein, CENP-B. Homologous sequences were not detected in other Drosophila species, suggesting that this transposable element family may be restricted to D. melanogaster.     

Comparative analysis of genome composition in Triticeae reveals strong variation in transposable element dynamics and nucleotide diversity

A 454 sequencing snapshot was utilised to investigate the genome composition and nucleotide diversity of transposable elements (TEs) for several Triticeae taxa, including Triticum aestivum, Hordeum vulgare, Hordeum spontaneum and Secale cereale together with relatives of the A, B and D genome donors of wheat, Triticum urartu (A), Aegilops speltoides (S) and Aegilops tauschii (D). Additional taxa containing the A genome, Triticum monococcum and its wild relative Triticum boeoticum, were also included. The main focus of the analysis was on the genomic composition of TEs as these make up at least 80% of the overall genome content. Although more than 200 TE families were identified in each species, approximately 50% of the overall genome comprised 12–15 TE families. The BARE1 element was the largest contributor to all genomes, contributing more than 10% to the overall genome. We also found that several TE families differ strongly in their abundance between species, indicating that TE families can thrive extremely successfully in one species while going virtually extinct in another. Additionally, the nucleotide diversity of BARE1 populations within individual genomes was measured. Interestingly, the nucleotide diversity in the domesticated barley H. vulgare cv. Barke was found to be twice as high as in its wild progenitor H. spontaneum, suggesting that the domesticated barley gained nucleotide diversity from the addition of different genotypes during the domestication and breeding process. In the rye/wheat lineage, sequence diversity of BARE1 elements was generally higher, suggesting that factors such as geographical distribution and mating systems might play a role in intragenomic TE diversity.     

Transposition of the Tol2 element, an Ac-like element from the Japanese medaka fish Oryzias latipes, in mouse embryonic stem cells

The Tol2 transposable element of the Japanese medaka fish belongs to the hAT family of transposons including hobo of Drosophila, Ac of maize, and Tam3 of snapdragon. To date, Tol2 is the only natural transposon in vertebrates that has ever been shown to encode a fully functional transposase. It has not been known, however, whether Tol2 can transpose in vertebrates other than fish. We report here transposition of Tol2 in mouse embryonic stem (ES) cells. We constructed a transposon donor plasmid containing a nonautonomous Tol2 element with the neomycin resistance gene and a helper plasmid capable of expressing the transposase and introduced the donor plasmid with various amounts of the helper plasmid by electroporation into mouse ES cells. The number of G418-resistant ES colonies increased as the amount of helper plasmid was increased, in a dose-dependent manner, indicating that the transposase activity elevated the integration efficiency. These G418-resistant ES colonies were cloned and the structure of the junction of the integrated Tol2 element and the genomic DNA was analyzed by inverse PCR. In those clones, Tol2 was surrounded by mouse genomic sequences and an 8-bp direct repeat was created adjacent to both ends of Tol2, indicating that Tol2 was integrated in the genome through transposition. The Tol2 transposon system is thus active in mouse as well as in fish. We propose that it should be used as a genetic tool to develop novel gene transfer, transgenesis, and mutagenesis methods in mammals.     

Tag1 is an autonomous transposable element that shows somatic excision in both Arabidopsis and tobacco.

Tag1 is a transposable element first identified as an insertion in the CHL1 gene of Arabidopsis. The chl1::Tag1 mutant originated from a plant (ecotype Landsberg erecta) that had been transformed with the maize transposon Activator (Ac), which is distantly related to Tag1. Genomic analysis of untransformed Landsberg erecta plants demonstrated that two identical Tag1 elements are present in the Landsberg erecta genome. To determine what provides transposase function for Tag1 transposition, we examined Tag1 excision in different genetic backgrounds. First, the chl1::Tag1 mutant was backcrossed to untransformed wild-type Arabidopsis plants to remove the Ac element(s) from the genome. F2 progeny that had no Ac elements but still retained Tag1 in the CHL1 gene were identified. Tag1 still excised in these Ac-minus progeny producing CHL1 revertants; therefore, Ac is not required for Tag1 excision. Next, Tag1 was inserted between a cauliflower mosaic virus 35S promoter and a beta-glucuronidase (GUS) marker gene and transformed into tobacco. Transformants showed blue-staining sectors indicative of Tag1 excision. Transgenic tobacco containing a defective Tag1 element, which was constructed in vitro by deleting an internal 1.4-kb EcoRI fragment, did not show blue-staining sectors. We conclude that Tag1 is an autonomous element capable of independent excision. The 35S-GUS::Tag1 construct was then introduced into Arabidopsis. Blue-staining sectors were found in cotyledons, leaves, and roots, showing that Tag1 undergoes somatic excision during vegetative development in its native host.     

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.