Conserved motifs and dynamic aspects of the terminal inverted repeat organization within Bari-like transposons

Heterologous survey of 130 DNA transposons in human cells highlights their functional divergence and expands the genome engineering toolbox

Experimental studies on DNA transposable elements (TEs) have been limited in scale, leading to a lack of understanding of the factors influencing transposition activity, evolutionary dynamics, and application potential as genome engineering tools. We predicted 130 active DNA TEs from 102 metazoan genomes and evaluated their activity in human cells. We identified 40 active (integration-competent) TEs, surpassing the cumulative number (20) of TEs found previously. With this unified comparative data, we found that the Tc1/mariner superfamily exhibits elevated activity, potentially explaining their pervasive horizontal transfers. Further functional characterization of TEs revealed additional divergence in features such as insertion bias. Remarkably, in CAR-T therapy for hematological and solid tumors, Mariner2_AG (MAG), the most active DNA TE identified, largely outperformed two widely used vectors, the lentiviral vector and the TE-based vector SB100X. Overall, this study highlights the varied transposition features and evolutionary dynamics of DNA TEs and increases the TE toolbox diversity.

Hiker, a new family of DNA transposons encoding transposases with DD35E motifs, displays a distinct phylogenetic relationship with most known DNA transposon families of IS630-Tc1-mariner (ITm)

DNA transposons play a crucial role in determining the size and structure of eukaryotic genomes. In this study, a new family of IS630-Tc1-mariner (ITm) DNA transposons, named Hiker (HK), was identified. HK is characterized by a DD35E catalytic domain and is distinct from all previously known families of the ITm group. Phylogenetic analyses showed that DD35E/Hiker forms a monophyletic clade with DD34E/Gambol, indicating that they may represent a separate superfamily of ITm. A total of 178 Hiker species were identified, with 170 found mainly in Actinopterygii, one in Chondrichthyes, six in Anura and one in Mollusca. Gambol (GM), on the other hand, are found in invertebrates, with 18 in Arthropoda and one in Platyhelminthes. Hiker transposons have a total length ranging from 2.14 to 3.67 kb and contain a single open reading frame that encodes a protein of approximately 370 amino acids (range 311-413 aa). They are flanked by short terminal inverted repeats (TIRs) of 16-30 base pairs and two base pair (TA) target-site duplications. In contrast, most transposons of the Gambol family have a total length of 1.35-5.96 kb, encode a transposase protein of approximately 350 amino acids (range 306-374 aa), and are flanked by TIRs that range from 32 to 1097 bp in length. Both Hiker and Gambol transposases have several conserved motifs, including helix-turn-helix (HTH) motifs and a DDE domain. Our study observed multiple amplification waves and repeated horizontal transfer (HT) events of HK transposons in vertebrate genomes, indicating their role in diversifying and shaping the genomes of Actinopterygii, Chondrichthyes, and Anura. Conversely, GM transposons showed few Horizontal transfer events. According to cell-based transposition assays, most HK transposons are likely inactive due to the truncated DNA binding domains of their transposases. We present an updated classification of the ITm group based on these findings, which will enhance the understanding of both the evolution of ITm transposons and that of their hosts.

What Have We Learned in 30 Years of Investigations on Bari Transposons?

Transposable elements (TEs) have been historically depicted as detrimental genetic entities that selfishly aim at perpetuating themselves, invading genomes, and destroying genes. Scientists often co-opt "special" TEs to develop new and powerful genetic tools, that will hopefully aid in changing the future of the human being. However, many TEs are gentle, rarely unleash themselves to harm the genome, and bashfully contribute to generating diversity and novelty in the genomes they have colonized, yet they offer the opportunity to develop new molecular tools. In this review we summarize 30 years of research focused on the Bari transposons. Bari is a "normal" transposon family that has colonized the genomes of several Drosophila species and introduced genomic novelties in the melanogaster species. We discuss how these results have contributed to advance the field of TE research and what future studies can still add to the current knowledge.

Transcriptionally promiscuous "blurry" promoters in Tc1/mariner transposons allow transcription in distantly related genomes

Background

We have recently described a peculiar feature of the promoters in two Drosophila Tc1-like elements, Bari1 and Bari3. The AT-richness and the presence of weak core-promoter motifs make these promoters, that we have defined “blurry”, able to activate transcription of a reporter gene in cellular systems as diverse as fly, human, yeast and bacteria. In order to clarify whether the blurry promoter is a specific feature of the Bari transposon family, we have extended this study to promoters isolated from three additional DNA transposon and from two additional LTR retrotransposons.

Results

Here we show that the blurry promoter is also a feature of two vertebrate transposable elements, Sleeping Beauty and Hsmar1, belonging to the Tc1/mariner superfamily. In contrast, this feature is not shared by the promoter of the hobo transposon, which belongs to the hAT superfamily, nor by LTR retrotransposon-derived promoters, which, in general, do not activate transcription when introduced into non-related genomes.

Conclusions

Our results suggest that the blurry promoter could be a shared feature of the members of the Tc1/mariner superfamily with possible evolutionary and biotechnological implications.

Electronic supplementary material

The online version of this article (10.1186/s13100-019-0155-6) contains supplementary material, which is available to authorized users.

Does the Promoter Constitute a Barrier in the Horizontal Transposon Transfer Process? Insight from Bari Transposons

The contribution of the transposons’ promoter in the horizontal transfer process is quite overlooked in the scientific literature. To shed light on this aspect we have mimicked the horizontal transfer process in laboratory and assayed in a wide range of hosts (fly, human, yeast and bacteria) the promoter activity of the 5′ terminal sequences in Bari1 and Bari3, two Drosophila transposons belonging to the Tc1-mariner superfamily. These sequences are able to drive the transcription of a reporter gene even in distantly related organisms at least at the episomal level. By combining bioinformatics and experimental approaches, we define two distinct promoter sequences for each terminal sequence analyzed, which allow transcriptional activity in prokaryotes and eukaryotes, respectively. We propose that the Bari family of transposons, and possibly other members of the Tc1-mariner superfamily, might have evolved “blurry promoters,” which have facilitated their diffusion in many living organisms through horizontal transfer.

Identification of Bari Transposons in 23 Sequenced Drosophila Genomes Reveals Novel Structural Variants, MITEs and Horizontal Transfer

Bari elements are members of the Tc1-mariner superfamily of DNA transposons, originally discovered in Drosophila melanogaster, and subsequently identified in silico in 11 sequenced Drosophila genomes and as experimentally isolated in four non-sequenced Drosophila species. Bari-like elements have been also studied for their mobility both in vivo and in vitro. We analyzed 23 Drosophila genomes and carried out a detailed characterization of the Bari elements identified, including those from the heterochromatic Bari1 cluster in D. melanogaster. We have annotated 401 copies of Bari elements classified either as putatively autonomous or inactive according to the structure of the terminal sequences and the presence of a complete transposase-coding region. Analyses of the integration sites revealed that Bari transposase prefers AT-rich sequences in which the TA target is cleaved and duplicated. Furthermore evaluation of transposon’s co-occurrence near the integration sites of Bari elements showed a non-random distribution of other transposable elements. We also unveil the existence of a putatively autonomous Bari1 variant characterized by two identical long Terminal Inverted Repeats, in D. rhopaloa. In addition, we detected MITEs related to Bari transposons in 9 species. Phylogenetic analyses based on transposase gene and the terminal sequences confirmed that Bari-like elements are distributed into three subfamilies. A few inconsistencies in Bari phylogenetic tree with respect to the Drosophila species tree could be explained by the occurrence of horizontal transfer events as also suggested by the results of dS analyses. This study further clarifies the Bari transposon’s evolutionary dynamics and increases our understanding on the Tc1-mariner elements’ biology.

Structural and sequence diversity of the transposon Galileo in the Drosophila willistoni genome

Background

Galileo is one of three members of the P superfamily of DNA transposons. It was originally discovered in Drosophila buzzatii, in which three segregating chromosomal inversions were shown to have been generated by ectopic recombination between Galileo copies. Subsequently, Galileo was identified in six of 12 sequenced Drosophila genomes, indicating its widespread distribution within this genus. Galileo is strikingly abundant in Drosophila willistoni, a neotropical species that is highly polymorphic for chromosomal inversions, suggesting a role for this transposon in the evolution of its genome.

Results

We carried out a detailed characterization of all Galileo copies present in the D. willistoni genome. A total of 191 copies, including 133 with two terminal inverted repeats (TIRs), were classified according to structure in six groups. The TIRs exhibited remarkable variation in their length and structure compared to the most complete copy. Three copies showed extended TIRs due to internal tandem repeats, the insertion of other transposable elements (TEs), or the incorporation of non-TIR sequences into the TIRs. Phylogenetic analyses of the transposase (TPase)-encoding and TIR segments yielded two divergent clades, which we termed Galileo subfamilies V and W. Target-site duplications (TSDs) in D. willistoni Galileo copies were 7- or 8-bp in length, with the consensus sequence GTATTAC. Analysis of the region around the TSDs revealed a target site motif (TSM) with a 15-bp palindrome that may give rise to a stem-loop secondary structure.

Conclusions

There is a remarkable abundance and diversity of Galileo copies in the D. willistoni genome, although no functional copies were found. The TIRs in particular have a dynamic structure and extend in different ways, but their ends (required for transposition) are more conserved than the rest of the element. The D. willistoni genome harbors two Galileo subfamilies (V and W) that diverged ~9 million years ago and may have descended from an ancestral element in the genome. Galileo shows a significant insertion preference for a 15-bp palindromic TSM.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-792) contains supplementary material, which is available to authorized users.

The Drosophila mojavensis Bari3 transposon: distribution and functional characterization

BACKGROUND

Bari-like transposons belong to the Tc1-mariner superfamily, and they have been identified in several genomes of the Drosophila genus. This transposon's family has been used as paradigm to investigate the complex dynamics underlying the persistence and structural evolution of transposable elements (TEs) within a genome. Three structural Bari variants have been identified so far and can be distinguished based on the organization of their terminal inverted repeats. Bari3 is the last discovered member of this family identified in Drosophila mojavensis, a recently emerged species of the Repleta group of the genus Drosophila.

RESULTS

We studied the insertion pattern of Bari3 in different D. mojavensis populations and found evidence of recent transposition activity. Analysis of the transposase domains unveiled the presence of a functional nuclear localization signal, as well as a functional binding domain. Using luciferase-based assays, we investigated the promoter activity of Bari3 as well as the interaction of its transposase with its left terminus. The results suggest that Bari3 is transposition-competent. Finally we demonstrated transposase transcript processing when the transposase gene is overexpressed in vivo and in vitro.

CONCLUSIONS

Bari3 displays very similar structural and functional features with its close relative, Bari1. Our results strongly suggest that Bari3 is an independent element that has generated genomic diversity in D. mojavensis. It can autonomously transcribe its transposase gene, which in turn can localize in the nucleus and bind the terminal inverted repeats of the transposon. Nevertheless, the identification of an unpredicted spliced form of the Bari3 transposase transcript allows us to hypothesize a control mechanism of its mobility based on mRNA processing. These results will aid the studies on the Bari family of transposons, which is intriguing for its widespread diffusion in Drosophilids coupled with a structural diversity generated during the evolution of Bari-like elements in their host genomes.

Functional characterization of the Bari1 transposition system

The transposons of the Bari family are mobile genetic elements widespread in the Drosophila genus. However, despite a broad diffusion, virtually no information is available on the mechanisms underlying their mobility. In this paper we report the functional characterization of the Bari elements transposition system. Using the Bari1 element as a model, we investigated the subcellular localization of the transposase, its physical interaction with the transposon, and its catalytic activity. The Bari1 transposase localized in the nucleus and interacted with the terminal sequences of the transposon both in vitro and in vivo, however, no transposition activity was detected in transposition assays. Profiling of mRNAs expressed by the transposase gene revealed the expression of abnormal, internally processed transposase transcripts encoding truncated, catalytically inactive transposase polypeptides. We hypothesize that a post-transcriptional control mechanism produces transposase-derived polypeptides that effectively repress transposition. Our findings suggest further clues towards understanding the mechanisms that control transposition of an important class of mobile elements, which are both an endogenous source of genomic variability and widely used as transformation vectors/biotechnological tools.

Identification of multiple binding sites for the THAP domain of the Galileo transposase in the long terminal inverted-repeats

Galileo is a DNA transposon responsible for the generation of several chromosomal inversions in Drosophila. In contrast to other members of the P-element superfamily, it has unusually long terminal inverted-repeats (TIRs) that resemble those of Foldback elements. To investigate the function of the long TIRs we derived consensus and ancestral sequences for the Galileo transposase in three species of Drosophilids. Following gene synthesis, we expressed and purified their constituent THAP domains and tested their binding activity towards the respective Galileo TIRs. DNase I footprinting located the most proximal DNA binding site about 70 bp from the transposon end. Using this sequence we identified further binding sites in the tandem repeats that are found within the long TIRs. This suggests that the synaptic complex between Galileo ends may be a complicated structure containing higher-order multimers of the transposase. We also attempted to reconstitute Galileo transposition in Drosophila embryos but no events were detected. Thus, although the limited numbers of Galileo copies in each genome were sufficient to provide functional consensus sequences for the THAP domains, they do not specify a fully active transposase. Since the THAP recognition sequence is short, and will occur many times in a large genome, it seems likely that the multiple binding sites within the long, internally repetitive, TIRs of Galileo and other Foldback-like elements may provide the transposase with its binding specificity.

Striking structural dynamism and nucleotide sequence variation of the transposon Galileo in the genome of Drosophila mojavensis

Background

Galileo is a transposable element responsible for the generation of three chromosomal inversions in natural populations of Drosophila buzzatii. Although the most characteristic feature of Galileo is the long internally-repetitive terminal inverted repeats (TIRs), which resemble the Drosophila Foldback element, its transposase-coding sequence has led to its classification as a member of the P-element superfamily (Class II, subclass 1, TIR order). Furthermore, Galileo has a wide distribution in the genus Drosophila, since it has been found in 6 of the 12 Drosophila sequenced genomes. Among these species, D. mojavensis, the one closest to D. buzzatii, presented the highest diversity in sequence and structure of Galileo elements.

Results

In the present work, we carried out a thorough search and annotation of all the Galileo copies present in the D. mojavensis sequenced genome. In our set of 170 Galileo copies we have detected 5 Galileo subfamilies (C, D, E, F, and X) with different structures ranging from nearly complete, to only 2 TIR or solo TIR copies. Finally, we have explored the structural and length variation of the Galileo copies that point out the relatively frequent rearrangements within and between Galileo elements. Different mechanisms responsible for these rearrangements are discussed.

Conclusions

Although Galileo is a transposable element with an ancient history in the D. mojavensis genome, our data indicate a recent transpositional activity. Furthermore, the dynamism in sequence and structure, mainly affecting the TIRs, suggests an active exchange of sequences among the copies. This exchange could lead to new subfamilies of the transposon, which could be crucial for the long-term survival of the element in the genome.

Ancestral polymorphism and recent invasion of transposable elements in Drosophila species

BACKGROUND

During the evolution of transposable elements, some processes, such as ancestral polymorphisms and horizontal transfer of sequences between species, can produce incongruences in phylogenies. We investigated the evolutionary history of the transposable elements Bari and 412 in the sequenced genomes of the Drosophila melanogaster group and in the sibling species D. melanogaster and D. simulans using traditional phylogenetic and network approaches.

RESULTS

Maximum likelihood (ML) phylogenetic analyses revealed incongruences and unresolved relationships for both the Bari and 412 elements. The DNA transposon Bari within the D. ananassae genome is more closely related to the element of the melanogaster complex than to the sequence in D. erecta, which is inconsistent with the species phylogeny. Divergence analysis and the comparison of the rate of synonymous substitutions per synonymous site of the Bari and host gene sequences explain the incongruence as an ancestral polymorphism that was inherited stochastically by the derived species. Unresolved relationships were observed in the ML phylogeny of both elements involving D. melanogaster, D. simulans and D. sechellia. A network approach was used to attempt to resolve these relationships. The resulting tree suggests recent transfers of both elements between D. melanogaster and D. simulans. The divergence values of the elements between these species support this conclusion.

CONCLUSIONS

We showed that ancestral polymorphism and recent invasion of genomes due to introgression or horizontal transfer between species occurred during the evolutionary history of the Bari and 412 elements in the melanogaster group. These invasions likely occurred in Africa during the Pleistocene, before the worldwide expansion of D. melanogaster and D. simulans.

The role of vertical and horizontal transfer in the evolution of Paris-like elements in drosophilid species

The transposable element (TE) Paris was described in a Drosophila virilis strain (virilis species group) as causing a hybrid dysgenesis with other mobile genetic elements. Since then, the element Paris has only been found in D. buzzatii, a species from the repleta group. In this study, we performed a search for Paris-like elements in 56 species of drosophilids to improve the knowledge about the distribution and evolution of this element. Paris-like elements were found in 30 species from the Drosophila genus, 15 species from the Drosophila subgenus and 15 species from the Sophophora subgenus. Analysis of the complete sequences obtained from the complete available Drosophila genomes has shown that there are putative active elements in five species (D. elegans, D. kikkawai, D. ananassae, D. pseudoobscura and D. mojavensis). The Paris-like elements showed an approximately 242-bp-long terminal inverted repeats in the 5' and 3' boundaries (called LIR: long inverted repeat), with two 28-bp-long direct repeats in each LIR. All potentially active elements presented degeneration in the internal region of terminal inverted repeat. Despite the degeneration of the LIR, the distance of 185 bp between the direct repeats was always maintained. This conservation suggests that the spacing between direct repeats is important for transposase binding. The distribution analysis showed that these elements are widely distributed in other Drosophila groups beyond the virilis and repleta groups. The evolutionary analysis of Paris-like elements suggests that they were present as two subfamilies with the common ancestor of the Drosophila genus. Since then, these TEs have been primarily maintained by vertical transmission with some events of stochastic loss and horizontal transfer.