Quantitatively Increased Somatic Transposition of Transposable Elements in Drosophila Strains Compromised for RNAi

In Drosophila melanogaster, small RNAs homologous to transposable elements (TEs) are of two types: piRNA (piwi-interacting RNA) with size 23-29nt and siRNA (small interfering RNA) with size 19-22nt. The siRNA pathway is suggested to silence TE activities in somatic tissues based on TE expression profiles, but direct evidence of transposition is lacking. Here we developed an efficient FISH (fluorescence in Situ hybridization) based method for polytene chromosomes from larval salivary glands to reveal new TE insertions. Analysis of the LTR-retrotransposon 297 and the non-LTR retroposon DOC shows that in the argonaut 2 (Ago2) and Dicer 2 (Dcr2) mutant strains, new transposition events are much more frequent than in heterozygous strains or wild type strains. The data demonstrate that the siRNA pathway represses TE transposition in somatic cells. Nevertheless, we found that loss of one functional copy of Ago2 or Dcr2 increases somatic transpositions of the elements at a lower level depending on the genetic background, suggesting a quantitative role for RNAi core components on mutation frequency.     

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.     

Comparative analysis of transposable elements in the melanogaster subgroup sequenced genomes

Transposable elements (TEs) are indwelling components of genomes, and their dynamics have been a driving force in genome evolution. Although we now have more information concerning their amounts and characteristics in various organisms, we still have little data from overall comparisons of their sequences in very closely-related species. While the Drosophila melanogaster genome has been extensively studied, we have only limited knowledge regarding the precise TE sequences in the genomes of the related species Drosophila simulans, Drosophila sechellia and Drosophila yakuba. In this study we analyzed the number and structure of TE copies in the sequenced genomes of these four species. Our findings show that, unexpectedly, the number of TE insertions in D. simulans is greater than that in D. melanogaster, but that most of the copies in D. simulans are degraded and in small fragments, as in D. sechellia and D. yakuba. This suggests that all three species were invaded by numerous TEs a long time ago, but have since regulated their activity, as the present TE copies are degraded, with very few full-length elements. In contrast, in D. melanogaster, a recent activation of TEs has resulted in a large number of almost-identical TE copies. We have detected variants of some TEs in D. simulans and D. sechellia, that are almost identical to the reference TE sequences in D. melanogaster, suggesting that D. melanogaster has recently been invaded by active TE variants from the other species. Our results indicate that the three species D. simulans, D. sechellia, and D. yakuba seem to be at a different stage of their TE life cycle when compared to D. melanogaster. Moreover, we show that D. melanogaster has been invaded by active TE variants for several TE families likely to come from D. simulans or the ancestor of D. simulans and D. sechellia. The numerous horizontal transfer events implied to explain these results could indicate introgression events between these species.     

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.     

Intracellular battlegrounds: conflict and cooperation between transposable elements

Transposable elements (TEs) are genomic parasites that amplify their own representation on hosts' chromosomes by inserting into new positions. It is traditionally thought that their copy number is regulated by purifying selection that eliminates hosts with higher than average TE abundance. Here, we stress that selection due to beneficial or harmful interactions between TEs introduces a whole new dimension, with implications for TE evolutionary trajectories and TE loads on hosts. This framework poses new questions requiring conceptual and experimental advances. Considering primarily Drosophila data, we make a case for within host selection on TEs by thinking expansively about the lifecycle of several TE families.     

Considering transposable element diversification in de novo annotation approaches

Transposable elements (TEs) are mobile, repetitive DNA sequences that are almost ubiquitous in prokaryotic and eukaryotic genomes. They have a large impact on genome structure, function and evolution. With the recent development of high-throughput sequencing methods, many genome sequences have become available, making possible comparative studies of TE dynamics at an unprecedented scale. Several methods have been proposed for the de novo identification of TEs in sequenced genomes. Most begin with the detection of genomic repeats, but the subsequent steps for defining TE families differ. High-quality TE annotations are available for the Drosophila melanogaster and Arabidopsis thaliana genome sequences, providing a solid basis for the benchmarking of such methods. We compared the performance of specific algorithms for the clustering of interspersed repeats and found that only a particular combination of algorithms detected TE families with good recovery of the reference sequences. We then applied a new procedure for reconciling the different clustering results and classifying TE sequences. The whole approach was implemented in a pipeline using the REPET package. Finally, we show that our combined approach highlights the dynamics of well defined TE families by making it possible to identify structural variations among their copies. This approach makes it possible to annotate TE families and to study their diversification in a single analysis, improving our understanding of TE dynamics at the whole-genome scale and for diverse species.     

TEMP: a computational method for analyzing transposable element polymorphism in populations

Insertions and excisions of transposable elements (TEs) affect both the stability and variability of the genome. Studying the dynamics of transposition at the population level can provide crucial insights into the processes and mechanisms of genome evolution. Pooling genomic materials from multiple individuals followed by high-throughput sequencing is an efficient way of characterizing genomic polymorphisms in a population. Here we describe a novel method named TEMP, specifically designed to detect TE movements present with a wide range of frequencies in a population. By combining the information provided by pair-end reads and split reads, TEMP is able to identify both the presence and absence of TE insertions in genomic DNA sequences derived from heterogeneous samples; accurately estimate the frequencies of transposition events in the population and pinpoint junctions of high frequency transposition events at nucleotide resolution. Simulation data indicate that TEMP outperforms other algorithms such as PoPoolationTE, RetroSeq, VariationHunter and GASVPro. TEMP also performs well on whole-genome human data derived from the 1000 Genomes Project. We applied TEMP to characterize the TE frequencies in a wild Drosophila melanogaster population and study the inheritance patterns of TEs during hybrid dysgenesis. We also identified sequence signatures of TE insertion and possible molecular effects of TE movements, such as altered gene expression and piRNA production. TEMP is freely available at github: https://github.com/JialiUMassWengLab/TEMP.git.     

THE TEMPO AND MODE OF EVOLUTION OF TRANSPOSABLE ELEMENTS AS REVEALED BY MOLECULAR PHYLOGENIES RECONSTRUCTED FROM MOSQUITO GENOMES

Although many mathematical models exist predicting the dynamics of transposable elements (TEs), there is a lack of available empirical data to validate these models and inherent assumptions. Genomes can provide a snapshot of several TE families in a single organism, and these could have their demographics inferred by coalescent analysis, allowing for the testing of theories on TE amplification dynamics. Using the available genomes of the mosquitoes Aedes aegypti and Anopheles gambiae, we indicate that such an approach is feasible. Our analysis follows four steps: (1) mining the two mosquito genomes currently available in search of TE families; (2) fitting, to selected families found in (1), a phylogeny tree under the general time‐reversible (GTR) nucleotide substitution model with an uncorrelated lognormal (UCLN) relaxed clock and a nonparametric demographic model; (3) fitting a nonparametric coalescent model to the tree generated in (2); and (4) fitting parametric models motivated by ecological theories to the curve generated in (3).     

What makes transposable elements move in the Drosophila genome?

Transposable elements (TEs), by their capacity of moving and inducing mutations in the genome, are considered important drivers of species evolution. The successful invasions of TEs in genomes, despite their mutational properties, are an apparent paradox. TEs' transposition is usually strongly regulated to low value, but in some cases these elements can also show high transposition rates, which has been associated sometimes to changes in environmental conditions. It is evident that factors susceptible to induce transpositions in natural populations contribute to TE perpetuation. Different factors were proposed as causative agents of TE mobilization in a wide range of organisms: biotic and abiotic stresses, inter- and intraspecific crosses and populational factors. However, there is no clear evidence of the factors capable of inducing TE mobilization in Drosophila, and data on laboratory stocks show contradictory results. The aim of this review is to have an update critical revision about mechanisms promoting transposition of TEs in Drosophila, and to provide to the readers a global vision of the dynamics of these genomic elements in the Drosophila genome.     

Abundance and chromosomal distribution of six Drosophila buzzatii transposons: BuT1, BuT2, BuT3, BuT4, BuT5, and BuT6

The abundance and chromosomal distribution of six class-II transposable elements (TEs) of Drosophila buzzatii have been analyzed by Southern blotting and in situ hybridization. These six transposons had been previously found at the breakpoints of inversions 2j and 2q ⁷ of D. buzzatii. These two polymorphic inversions were generated by an ectopic recombination event between two copies of Galileo, a Foldback element. The four breakpoints became hotspots for TE insertions after the generation of the inversion and the transposons analyzed in this work are considered to be secondary invaders of these regions. Insertions of the six transposons are present in the euchromatin but show an increased density in the pericentromeric euchromatin–heterochromatin transition region and the dot chromosome. They are also more abundant in the inverted segments of chromosome 2 rearrangements. We further observed that the accumulation of TE insertions varies between elements and is correlated between dot, proximal regions, and inverted segments. These observations fully agree with previous data in Drosophila melanogaster and support recombination rate as the chief force explaining the chromosomal distribution of TEs.Electronic Supplementary Material Supplementary material is available in the online version of this article at and is accessible for authorized users.Sequence data from this article have been deposited in the EMBL/GenBank Data Libraries under accession number DQ402469.     

CellLine, a stochastic cell lineage simulator

Summary: We present CellLine, a simulator of the dynamics ofgene regulatory networks (GRN) in the cells of a lineage. Fromuser-defined reactions and initial substance quantities, itgenerates cell lineages, i.e. genealogic pedigrees of cellsrelated through mitotic division. Each cell's dynamics is drivenby a delayed stochastic simulation algorithm (delayed SSA),allowing multiple time delayed reactions. The cells of the lineage can be individually subject to ‘perturbations’,such as gene deletion, duplication and mutation. External interventions,such as adding or removing a substance at a given moment, canbe specified. Cell differentiation lineages, where differentiationis stochastically driven or externally induced, can be modeledas well. Finally, CellLine can generate and simulate the dynamicsof multiple copies of any given cell of the lineage. As examples of CellLine use, we simulate the following systems:cell lineages containing a model of the P53-Mdm2 feedback loop,a differentiation lineage where each cell contains a 4 generepressilator (a bistable circuit), a model of the differentiationof the cells of the retinal mosaic required for color visionin Drosophila melanogaster, where the differentiation pathwaydepends on one substance's concentration that is controlledby a stochastic process, and a 9 gene GRN to illustrate theadvantage of using CellLine rather than simulating multipleindependent cells, in cases where the cells of the lineage aredynamically correlated. Availability: The CellLine program, instructions and examplesare available at http://www.cs.tut.fi/~sanchesr/CellLine/CellLine.html Contact: andre.sanchesribeiro{at}tut.fi Associate Editor: Limsoon Wong     

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.     

ROBIN: a tool for genome rearrangement of block-interchanges

Summary: ROBIN is a web server for analyzing genome rearrangementof block-interchanges between two chromosomal genomes. It takestwo or more linear/circular chromosomes as its input, and computesthe number of minimum block-interchange rearrangements betweenany two input chromosomes for transforming one chromosome intoanother and also determines an optimal scenario taking thisnumber of rearrangements. The input can be either bacterial-sizesequence data or landmark-order data. If the input is sequencedata, ROBIN will automatically search for the identical landmarksthat are the homologous/conserved regions shared by all theinput sequences. Availability: ROBIN is freely accessed at http://genome.life.nctu.edu.tw/ROBIN Contact: cllu{at}mail.nctu.edu.tw     

Gross chromosome rearrangements mediated by transposable elements in Drosophila melanogaster

A combination of cytogenetic and molecular analyses has shown that several different transposable elements are involved in the restructuring of Drosophila chromosomes. Two kinds of elements, P and hobo, are especially prone to induce chromosome rearrangements. The mechanistic details of this process are unclear, but, at least some of the time, it seems to involve ectopic recombination between elements inserted at different chromosomal sites; the available data suggest that these ectopic recombination events are much more likely to occure between elements in the same chromosome than between elements in different chromosomes. Other Drosophila transposons also appear to mediate chromosome restructuring by ectopic recombination; these include the retrotransposons BEL, roo, Docand I and the foldback element FB. In addition, two retrotransposons, HeT-A and TART, have been found to be associated specifically with the ends of Drosophila chromosomes. Very limited data indicate that transposon-mediated chromosome restructuring is occurring in natural populations of Drosophila. This suggests that transposable elements may help to shape the structure of the Drosophila genome and implies that they may have a similar role in other organisms.     

A precise and scalable method for querying genes in chromosomal banding regions based on cytogenetic annotations

Motivation: Staining the human metaphase chromosomes revealscharacteristic banding patterns known as cytogenetic bands orcytobands. Using technologies based on metaphase chromosomes,researchers have accumulated much knowledge about the correlationsbetween human diseases and specific cytoband aberrations, indicatingthe presence of disease-associated genes in those bands. Withthe progress of human genome project and techniques such asfluorescent in situ hybridization, many genes have been assignedto the cytobands and annotated in public databases, making itpossible to find all genes in the disease-related cytobandsthrough database queries. However, finding genes in cytobandsremains an imprecise process, partly due to the insufficiencyof current methods for cytoband queries, especially for thosebased on cytogenetic annotations. Results: By transforming the cytoband annotations into numericalsegments, a new query method is developed that is able to accuratelydefine any cytogenetic ranges in human chromosomes. A querysystem (designated cytoband query sys CQS) is implemented usingcytogenetic annotations in the public domain. Judged by a performancetest, CQS executed as accurately as expected using cytogeneticannotations from NCBI Map Viewer. The new method is scalableand can be applied to genomes from other species. Availability: The CQS is freely accessible over the Internetat http://moris.csie.ncku.edu.tw/cqs/ Contact: clh9{at}mail.ncku.edu.tw Supplementary information: http://moris.csie.ncku.edu.tw/cqs/     

Updating of transposable element annotations from large wheat genomic sequences reveals diverse activities and gene associations

Triticeae species (including wheat, barley and rye) have huge and complex genomes due to polyploidization and a high content of transposable elements (TEs). TEs are known to play a major role in the structure and evolutionary dynamics of Triticeae genomes. During the last 5 years, substantial stretches of contiguous genomic sequence from various species of Triticeae have been generated, making it necessary to update and standardize TE annotations and nomenclature. In this study we propose standard procedures for these tasks, based on structure, nucleic acid and protein sequence homologies. We report statistical analyses of TE composition and distribution in large blocks of genomic sequences from wheat and barley. Altogether, 3.8 Mb of wheat sequence available in the databases was analyzed or re-analyzed, and compared with 1.3 Mb of re-annotated genomic sequences from barley. The wheat sequences were relatively gene-rich (one gene per 23.9 kb), although wheat gene-derived sequences represented only 7.8% (159 elements) of the total, while the remainder mainly comprised coding sequences found in TEs (54.7%, 751 elements). Class I elements [mainly long terminal repeat (LTR) retrotransposons] accounted for the major proportion of TEs, in terms of sequence length as well as element number (83.6% and 498, respectively). In addition, we show that the gene-rich sequences of wheat genome A seem to have a higher TE content than those of genomes B and D, or of barley gene-rich sequences. Moreover, among the various TE groups, MITEs were most often associated with genes: 43.1% of MITEs fell into this category. Finally, the TRIM and copia elements were shown to be the most active TEs in the wheat genome. The implications of these results for the evolution of diploid and polyploid wheat species are discussed. Electronic Supplementary Material Supplementary material is available for this article at     

An Analysis of Male-Recombination Elements in a Natural Population of DROSOPHILA MELANOGASTER in South Texas

Genomes from a group of Drosophila melanogaster collected from a natural population at San Benito, South Texas, in March of 1975 were analyzed for the presence of male-recombination elements. All three autosomes and both sex chromosomes were examined, with emphasis placed on the two major autosomes, the second and third chromosomes. In samples of 16 second and 16 third chromosomes, at least half, but not all, of each were found to carry male-recombination elements. It is suggested, although the data are not conclusive, that some of the fourth, X, and Y chromosomes might also be associated with male-recombination elements.—When a male-recombination element, or elements, was located in the second chromosome, relatively more male recombination was induced in the second than in the third chromosome. This situation was reversed when the element(s) was located in the third chromosome.—Distortion of transmission frequency, one of the characteristics of previously studied second chromosome lines associated with male recombination, was confirmed for these second chromosomes that carried male-recombination elements. Similar, but less pronounced, distortion was observed for the third chromosome lines that carried male-recombination elements.     

The evolution of retrotransposon regulatory regions and its consequences on the Drosophila melanogaster and Homo sapiens host genomes

It has now been established that transposable elements (TEs) make up a variable, but significant proportion of the genomes of all organisms, from Bacteria to Vertebrates. However, in addition to their quantitative importance, there is increasing evidence that TEs also play a functional role within the genome. In particular, TE regulatory regions can be viewed as a large pool of potential promoter sequences for host genes. Studying the evolution of regulatory region of TEs in different genomic contexts is therefore a fundamental aspect of understanding how a genome works. In this paper, we first briefly describe what is currently known about the regulation of TE copy number and activity in genomes, and then focus on TE regulatory regions and their evolution. We restrict ourselves to retrotransposons, which are the most abundant class of eukaryotic TEs, and analyze their evolution and the subsequent consequences for host genomes. Particular attention is paid to much-studied representatives of the Vertebrates and Invertebrates, Homo sapiens and Drosophila melanogaster, respectively, for which high quality sequenced genomes are available.