VisualTE: a graphical interface for transposable element analysis at the genomic scale

Background

Transposable elements are mobile DNA repeat sequences, known to have high impact on genes, genome structure and evolution. This has stimulated broad interest in the detailed biological studies of transposable elements. Hence, we have developed an easy-to-use tool for the comparative analysis of the structural organization and functional relationships of transposable elements, to help understand their functional role in genomes.

Results

We named our new software VisualTE and describe it here. VisualTE is a JAVA stand-alone graphical interface that allows users to visualize and analyze all occurrences of transposable element families in annotated genomes. VisualTE reads and extracts transposable elements and genomic information from annotation and repeat data. Result analyses are displayed in several graphical panels that include location and distribution on the chromosome, the occurrence of transposable elements in the genome, their size distribution, and neighboring genes’ features and ontologies. With these hallmarks, VisualTE provides a convenient tool for studying transposable element copies and their functional relationships with genes, at the whole-genome scale, and in diverse organisms.

Conclusions

VisualTE graphical interface makes possible comparative analyses of transposable elements in any annotated sequence as well as structural organization and functional relationships between transposable elements and other genetic object. This tool is freely available at: http://lcb.cnrs-mrs.fr/spip.php?article867.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1351-5) contains supplementary material, which is available to authorized users.     

High copy numbers of multiple transposable element families in an Australian population of Drosophila simulans

Sudden mobilization of transposable elements in Drosophila is a well-reported phenomenon but one that usually affects no more than a few elements (one to four). We report here the existence of a D. simulans natural population (Canberra) from Australia, which had high copy numbers for various transposable elements (transposons, LTR retrotransposons and non-LTR retrotransposons). The impact of transposable elements on the host genome and populations is discussed.     

TARGeT: a web-based pipeline for retrieving and characterizing gene and transposable element families from genomic sequences

Gene families compose a large proportion of eukaryotic genomes. The rapidly expanding genomic sequence database provides a good opportunity to study gene family evolution and function. However, most gene family identification programs are restricted to searching protein databases where data are often lagging behind the genomic sequence data. Here, we report a user-friendly web-based pipeline, named TARGeT (Tree Analysis of Related Genes and Transposons), which uses either a DNA or amino acid ‘seed’ query to: (i) automatically identify and retrieve gene family homologs from a genomic database, (ii) characterize gene structure and (iii) perform phylogenetic analysis. Due to its high speed, TARGeT is also able to characterize very large gene families, including transposable elements (TEs). We evaluated TARGeT using well-annotated datasets, including the ascorbate peroxidase gene family of rice, maize and sorghum and several TE families in rice. In all cases, TARGeT rapidly recapitulated the known homologs and predicted new ones. We also demonstrated that TARGeT outperforms similar pipelines and has functionality that is not offered elsewhere.     

Detection of new transposable element families in Drosophila melanogaster and Anopheles gambiae genomes

The techniques that are usually used to detect transposable elements (TEs) in nucleic acid sequences rely on sequence similarity with previously characterized elements. However, these methods are likely to miss many elements in various organisms. We tested two strategies for the detection of unknown elements. The first, which we call "TBLASTX strategy," searches for TE sequences by comparing the six-frame translations of the nucleic acid sequences of known TEs with the genomic sequence of interest. The second, "repeat-based strategy," searches genomic sequences for long repeats and clusters them in groups of similar sequences. TE copies from a given family are expected to cluster together. We tested the Drosophila melanogaster genomic sequence and the recently sequenced Anopheles gambiae genome in which most TEs remain unknown. We showed that the "TBLASTX strategy" is very efficient as it detected at least 332 new TE families in D. melanogaster and 400 in A. gambiae. This was unexpected in Drosophila as TEs of this organism have been extensively studied. The "repeat-based strategy" appeared to be very inefficient because of two problems: (i) TE copies are heavily deleted and few copies share homologous regions, and (ii) segmental duplications are frequent and it is not easy to distinguish them from TE copies.     

Genome evolution: sex and the transposable element

Mating systems are thought to play an important role in determining the fate of genomic parasites such as transposable elements. This is supported by recent studies which indicate that asexual genomes may be structured very differently to those of sexual species.     

Regulation of the transposable element mariner

The mariner/ Tc1 superfamily of transposable elements is widely distributed in animal genomes and is especially prevalent in insects. Their wide distribution results from their ability to be disseminated among hosts by horizontal transmission and also by their ability to persist in genomes through multiple speciation events. Although a great deal is known about the molecular mechanisms of transposition and excision, very little is known about the mechanisms by which transposition is controlled within genomes. The issue of mariner/Tc1 regulation is critical in view of the great interest in these elements as vectors for germline transformation of insect pests and vectors of human disease. Several potentially important regulatory mechanisms have been identified in studies of genetically engineered mariner elements. One mechanism is overproduction inhibition, in which excessive wild-type transposase reduces the rate of excision of a target element. A second mechanism is mediated by certain mutant transposase proteins, which antagonize the activity of the wild-type transposase. The latter process may help explain why the vast majority of MLEs in nature undergo ‘vertical inactivation’ by multiple mutations and, eventually, stochastic loss. Another potential mechanism of regulation may result from transposase titration by defective elements that retain their DNA binding sites and ability to transpose. There is also evidence that some mariner/Tc1 elements can be mobilized in a type of hybrid dysgenesis. This revised version was published online in August 2006 with corrections to the Cover Date.     

ModuleOrganizer: detecting modules in families of transposable elements

Background  

Most known eukaryotic genomes contain mobile copied elements called transposable elements. In some species, these elements account for the majority of the genome sequence. They have been subject to many mutations and other genomic events (copies, deletions, captures) during transposition. The identification of these transformations remains a difficult issue. The study of families of transposable elements is generally founded on a multiple alignment of their sequences, a critical step that is adapted to transposons containing mostly localized nucleotide mutations. Many transposons that have lost their protein-coding capacity have undergone more complex rearrangements, needing the development of more complex methods in order to characterize the architecture of sequence variations.     

Mixed-mode loading of the cement-bone interface: a finite element study

While including the cement-bone interface of complete cemented hip reconstructions is crucial to correctly capture their response, its modelling is often overly simplified. In this study, the mechanical mixed-mode response of the cement-bone interface is investigated, taking into account the effects of the well-defined microstructure that characterises the interface. Computed tomography-based plain strain finite element analyses models of the cement-bone interface are built and loaded in multiple directions. Periodic boundaries are considered and the failure of the cement and bone fractions by cracking of the bulk components are included. The results compare favourably with experimental observations. Surprisingly, the analyses reveal that under shear loading no failure occurs and considerable normal compression is generated to prevent interface dilation. Reaction forces, crack patterns and stress fields provide more insight into the mixed-mode failure process. Moreover, the cement-bone interface analyses provide details which can serve as a basis for the development of a cohesive law.     

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.     

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     

The Tol2 transposable element of the medaka fish: an active DNA-based element naturally occurring in a vertebrate genome

Several DNA-based transposable elements are known to be present in vertebrate genomes, but few of them have been demonstrated to be active. The Tol2 element of the medaka fish is one such element and, therefore, is potentially useful for developing a gene tagging system and other molecular biological tools applicable to vertebrates. Towards this goal, analyses of the element at the molecular, cellular and population levels are in progress. Results so far obtained are described here.     

Nucleotide sequence of the maize transposable element Mul

A cloned DNA fragment from the maize allele Adhl-S3034 contains all of Mul, an insertion element involved in Robertson's Mutator activity. The element is 1367 base pairs (bp) long and is flanked by nine bp direct repeats of insertion site DNA. It has inverted terminal repeats of 215 and 213 bp showing 95% homology. Within the element are two direct repeats of 104 bp showing 96% homology. Four open reading frames (ORFs) were found, two in each DNA strand. Mul can be divided into two halves, each containing one terminal inverted repeat, an internal direct repeat, and two overlapping ORFs. The GC content of each half is high (70%), while that of a central 60 base portion of the element is low (26%). The central region contains the only sequence resembling the TAATA Goldberg and Hogness eukaryotic promoter signal. Multiple copies of DNA sequences related to Mul found in Mutator maize plants are generally similar in organization to the cloned element. A larger version containing a discrete 300 to 400 base pair insertion was found in some Mutator lines.     

DNA sequence of the transposable element IS1

Summary The nucleotide sequence of an IS1 element recently transposed into the lacI gene is reported. This sequence is nearly identical to one previously reported for another IS1 element (Ohtsubo and Ohtsubo, 1978). The implications of this similarity are discussed. The sizes of potential polypeptides encoded in the IS1 DNA have been determined and possible roles for these peptides in the illegitimate recombination events mediated by the element are considered.     

Evidence for interspecific transfer of the transposable element mariner betweenDrosophila andZaprionus

Summary The transposable element mariner occurs widely in themelanogaster species group ofDrosophila. However, in drosophilids outside of themelanogaster species group, sequences showing strong DNA hybridization with mariner are found only in the genusZaprionus. the mariner sequence obtained fromZaprionus tuberculatus is 97% identical with that fromDrosophila mauritiana, a member of themelanogaster species subgroup, whereas a mariner sequence isolated fromDrosophila tsacasi is only 92% identical with that fromD. mauritiana. BecauseD. tsacasi is much more closely related toD. mauritiana than isZaprionus, the presence of mariner inZaprionus may result from horizontal transfer. In order to confirm lack of a close phylogenetic relationship between the genusZaprionus and themelanogaster species group, we compared the alcohol dehydrogenase (Adh) sequences among these species. The results show that the coding region of Adh is only 82% identical betweenZ. tuberculatus andD. mauritiana, as compared with 90% identical betweenD. tsacasi andD. mauritiana. Furthermore, the mariner gene phylogeny obtained by maximum likelihood and maximum parsimony analyses is discordant with the species phylogeny estimated by using the Adh genes. The only inconsistency in the mariner gene phylogeny is in the placement of theZaprionus mariner sequence, which clusters with mariner fromDrosophila teissieri andDrosophila yakuba in themelanogaster species subgroup. These results strongly suggest horizontal transfer.     

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.     

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.