Tx1: a transposable element from Xenopus laevis with some unusual properties

A transposable element annotation pipeline and expression analysis reveal potentially active elements in the microalga Tisochrysis lutea

BACKGROUND

Transposable elements (TEs) are mobile DNA sequences known as drivers of genome evolution. Their impacts have been widely studied in animals, plants and insects, but little is known about them in microalgae. In a previous study, we compared the genetic polymorphisms between strains of the haptophyte microalga Tisochrysis lutea and suggested the involvement of active autonomous TEs in their genome evolution.

RESULTS

To identify potentially autonomous TEs, we designed a pipeline named PiRATE (Pipeline to Retrieve and Annotate Transposable Elements, download: https://doi.org/10.17882/51795 ), and conducted an accurate TE annotation on a new genome assembly of T. lutea. PiRATE is composed of detection, classification and annotation steps. Its detection step combines multiple, existing analysis packages representing all major approaches for TE detection and its classification step was optimized for microalgal genomes. The efficiency of the detection and classification steps was evaluated with data on the model species Arabidopsis thaliana. PiRATE detected 81% of the TE families of A. thaliana and correctly classified 75% of them. We applied PiRATE to T. lutea genomic data and established that its genome contains 15.89% Class I and 4.95% Class II TEs. In these, 3.79 and 17.05% correspond to potentially autonomous and non-autonomous TEs, respectively. Annotation data was combined with transcriptomic and proteomic data to identify potentially active autonomous TEs. We identified 17 expressed TE families and, among these, a TIR/Mariner and a TIR/hAT family were able to synthesize their transposase. Both these TE families were among the three highest expressed genes in a previous transcriptomic study and are composed of highly similar copies throughout the genome of T. lutea. This sum of evidence reveals that both these TE families could be capable of transposing or triggering the transposition of potential related MITE elements.

CONCLUSION

This manuscript provides an example of a de novo transposable element annotation of a non-model organism characterized by a fragmented genome assembly and belonging to a poorly studied phylum at genomic level. Integration of multi-omics data enabled the discovery of potential mobile TEs and opens the way for new discoveries on the role of these repeated elements in genomic evolution of microalgae.

A transposable element annotation pipeline and expression analysis reveal potentially active elements in the microalga Tisochrysis lutea

Background

Transposable elements (TEs) are mobile DNA sequences known as drivers of genome evolution. Their impacts have been widely studied in animals, plants and insects, but little is known about them in microalgae. In a previous study, we compared the genetic polymorphisms between strains of the haptophyte microalga Tisochrysis lutea and suggested the involvement of active autonomous TEs in their genome evolution.

Results

To identify potentially autonomous TEs, we designed a pipeline named PiRATE (Pipeline to Retrieve and Annotate Transposable Elements, download: 10.17882/51795), and conducted an accurate TE annotation on a new genome assembly of T. lutea. PiRATE is composed of detection, classification and annotation steps. Its detection step combines multiple, existing analysis packages representing all major approaches for TE detection and its classification step was optimized for microalgal genomes. The efficiency of the detection and classification steps was evaluated with data on the model species Arabidopsis thaliana. PiRATE detected 81% of the TE families of A. thaliana and correctly classified 75% of them. We applied PiRATE to T. lutea genomic data and established that its genome contains 15.89% Class I and 4.95% Class II TEs. In these, 3.79 and 17.05% correspond to potentially autonomous and non-autonomous TEs, respectively. Annotation data was combined with transcriptomic and proteomic data to identify potentially active autonomous TEs. We identified 17 expressed TE families and, among these, a TIR/Mariner and a TIR/hAT family were able to synthesize their transposase. Both these TE families were among the three highest expressed genes in a previous transcriptomic study and are composed of highly similar copies throughout the genome of T. lutea. This sum of evidence reveals that both these TE families could be capable of transposing or triggering the transposition of potential related MITE elements.

Conclusion

This manuscript provides an example of a de novo transposable element annotation of a non-model organism characterized by a fragmented genome assembly and belonging to a poorly studied phylum at genomic level. Integration of multi-omics data enabled the discovery of potential mobile TEs and opens the way for new discoveries on the role of these repeated elements in genomic evolution of microalgae.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4763-1) contains supplementary material, which is available to authorized users.

piggyBac Transposon

The piggyBac transposon was originally isolated from the cabbage looper moth, Trichoplusia ni, in the 1980s. Despite its early discovery and dissimilarity to the other DNA transposon families, the piggyBac transposon was not recognized as a member of a large transposon superfamily for a long time. Initially, the piggyBac transposon was thought to be a rare transposon. This view, however, has now been completely revised as a number of fully sequenced genomes have revealed the presence of piggyBac-like repetitive elements. The isolation of active copies of the piggyBac-like elements from several distinct species further supported this revision. This includes the first isolation of an active mammalian DNA transposon identified in the bat genome. To date, the piggyBac transposon has been deeply characterized and it represents a number of unique characteristics. In general, all members of the piggyBac superfamily use TTAA as their integration target sites. In addition, the piggyBac transposon shows precise excision, i.e., restoring the sequence to its preintegration state, and can transpose in a variety of organisms such as yeasts, malaria parasites, insects, mammals, and even in plants. Biochemical analysis of the chemical steps of transposition revealed that piggyBac does not require DNA synthesis during the actual transposition event. The broad host range has attracted researchers from many different fields, and the piggyBac transposon is currently the most widely used transposon system for genetic manipulations.

Evolution of the Xenopus piggyBac Transposon Family TxpB: Domesticated and Untamed Strategies of Transposon Subfamilies

A new family, termed TxpB, of DNA transposons belonging to the piggyBac superfamily was found in 3 Xenopus species (Xenopus tropicalis, Xenopus laevis, and Xenopus borealis). Two TxpB subfamilies of Kobuta and Uribo1 were found in all the 3 species, and another subfamily termed Uribo2 was found in X. tropicalis. Molecular phylogenetic analyses of their open reading frames (ORFs) revealed that TxpB transposons have been maintained for over 100 Myr. Both the Uribo1 and the Uribo2 ORFs were present as multiple copies in each genome, and some of them were framed by terminal inverted repeat sequences. In contrast, all the Kobuta ORFs were present as a single copy in each genome and exhibited high evolutionary conservation, suggesting domestication of Kobuta genes by the host. Genomic insertion polymorphisms of the Uribo1 and Uribo2 transposons (nonautonomous type) were observed in a single species of X. tropicalis, indicating recent transposition events. Transfection experiments in cell culture revealed that an expression vector construct for the intact Uribo2 ORF caused precise excision of a nonautonomous Uribo2 element from the target vector construct but that for the Kobuta ORF did not. The present results support our viewpoint that some Uribo2 members are naturally active autonomous transposons, whereas Kobuta members may be domesticated by hosts.

Diversity and clustered distribution of retrotransposable elements in the compact genome of the pufferfish Tetraodon nigroviridis

We report the characterization and chromosomal distribution of retroelements in the compact genome of the pufferfish Tetraodon nigroviridis. We have reconstructed partial/complete retroelement sequences, established their phylogenetic relationship to other known eukaryotic retrotransposons, and performed double-color FISH analyses to gain new insights into their patterns of chromosomal distribution. We could identify 43 different reverse transcriptase retrotransposons belonging to the three major known subclasses (14 non-LTR retrotransposons from seven clades, 25 LTR retrotransposons representing the five major known groups, and four Penelope-like elements), and well as two SINEs (non-autonomous retroelements). Such a diversity of retrotransposable elements, which seems to be relatively common in fish but not in mammals, is astonishing in such a compact genome. The total number of retroelements was approximately 3000, roughly representing only 2.6% of the genome of T. nigroviridis. This is much less than in other vertebrate genomes, reflecting the compact nature of the genome of this pufferfish. Major differences in copy number were observed between different clades, indicating differential success in invading and persisting in the genome. Some retroelements displayed evidence of recent activity. Finally, FISH analysis showed that retrotransposable elements preferentially accumulate in specific heterochromatic regions of the genome of T. nigroviridis, revealing a degree of genomic compartmentalization not observed in the human genome.

Evolutionary dynamics in a novel L2 clade of non-LTR retrotransposons in Deuterostomia

The evolution of the novel L2 clade of non-long terminal repeat (LTR) retrotransposons and their evolutionary dynamics in Deuterostomia has been examined. The short-term evolution of long interspersed nuclear element 2s (LINE2s) has been studied in 18 reptilian species by analysis of a PCR amplified 0.7-kb fragment encoding the palm/fingers subdomain of reverse transcriptase (RT). Most of the reptilian LINE2s examined are inactive since they contain multiple stop codons, indels, or frameshift mutations that disrupt the RT. Analysis of reptilian LINE2s has shown a high degree of sequence divergence and an unexpectedly large number of deletions. The evolutionary dynamics of LINE2s in reptiles has been found to be complex. LINE2s are shown to form a novel clade of non-LTR retrotransposons that is well separated from the CR1 clade. This novel L2 clade is more widely distributed than previously thought, and new representatives have been discovered in echinoderms, insects, teleost fishes, Xenopus, Squamata, and marsupials. There is an apparent absence of LINE2s from different vertebrate classes, such as cartilaginous fishes, Archosauria (birds and crocodiles), and turtles. Whereas the LINE2s are present in echinoderms and teleost fishes in a conserved form, in most tetrapods only highly degenerated pseudogenes can be found. The predominance of inactive LINE2s in Tetrapoda indicates that, in the host genomes, only inactive copies are still present. The present data indicate that the vertical inactivation of LINE2s might have begun at the time of Tetrapoda origin, 400 MYA. The evolutionary dynamics of the L2 clade in Deuterostomia can be described as a gradual vertical inactivation in Tetrapoda, stochastic loss in Archosauria and turtles, and strict vertical transmission in echinoderms and teleost fishes.

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.

Structural analysis of the entire proopiomelanocortin gene of Xenopus laevis

In the pars intermedia of the pituitary the prohormone proopiomelanocortin (POMC) is tissue-specifically processed to, among other peptides, alpha-melanotropin (alpha MSH). In the South African clawed toad Xenopus laevis this hormone mediates the process of background adaptation: release of alpha-MSH causes darkening of the animal, while inhibition of alpha-MSH release results in a pale toad. Elevated release of alpha-MSH coincides with a higher rate of POMC gene transcription. The present study aims to find possible transcriptional regulatory elements in the Xenopus POMC gene. For that purpose the complete nucleotide sequence of the POMC gene and its 5'- and 3'- flanking regions were determined and analyzed. The Xenopus POMC gene promoter contains several regions which may be regulatory DNA elements in view of their similarity with corresponding regions of mammalian POMC gene promoters. In the rat POMC gene promoter, many of these regions represent protein-binding sequences. Besides the promoter sequence and the protein-coding sequences, no other segments with significant identity between the Xenopus and human POMC genes were found. Intron A of the Xenopus POMC gene contains a simple sequence, (TATC)76, and a JH12 repetitive element, while the 3'-flanking region contains a repetitive-EcoRI-monomer-2 element. Comparison of the JH12 sequence of the POMC gene with JH12 sequences from other Xenopus genes revealed a 335-bp consensus sequence which is flanked by a 30-bp inverted repeat. This JH12 consensus sequence is significantly larger than the previously reported JH12 core region. Alignment of intron B of the Xenopus POMC gene with database sequences revealed a consensus sequence of a novel Xenopus repetitive element of 330 bp flanked by a nearly perfect inverted repeat, indicating that this element may be a transposon-like element.

Sequence comparison of cellular and viral copies of host cell DNA insertions found in Autographa californica nuclear polyhedrosis virus

The nucleotide sequence and structural characteristics of two nonhomologous host DNA insertions of Spodoptera frugiperda (fall armyworm) origin isolated from few polyhedra mutants of the baculovirus Autographa californica nuclear polyhedrosis virus (AcMNPV) have been determined. Neither of the host insertions contain long open reading frames suggesting that cellular genes were not introduced into the viral genome. One of the host DNA insertions, IFP1.6, terminated in short imperfect inverted repeats flanked by a duplication of the target sequence, TTAA. Analysis of three cellular copies of this host insertion isolated from a lambda genomic library of S. frugiperda DNA revealed the termini to be highly conserved and also flanked by the same TTAA sequence. IFP1.6 was shown to share homology with a putative host insertion described in the genome of a baculovirus exhibiting wild-type plaque morphology. The second of the host insertions, IFP2.2, had a structure unique among host insertions described in baculoviruses. It lacks terminal repeats but is flanked by duplications of the 8-bp target site sequence. The cellular copy of this insertion was conserved in comparison to the viral copy and was also flanked by a direct 8-bp repeat. This is the first report of the analysis of cellular copies of host DNA insertions frequently associated with baculovirus FP mutants.

Transposon mutagenesis of baculoviruses: analysis of TFP3 lepidopteran transposon insertions at the FP locus of nuclear polyhedrosis viruses

We report the complete sequences of two representatives of the TFP3 transposable element family of the lepidopteran, Trichoplusia ni. These elements were isolated as insertions mobilized from the Lepidopteran host genome into two closely related nuclear polyhedrosis viruses (NPV) during infection. Both elements inserted within the 500-bp FP locus of the respective viral genomes (map units 36.0 to 37.0), causing a distinctive plaque morphology phenotype and the loss of a 25-kDa viral-specific protein. Both insertions occurred at the identical TTAA target site in the respective genomes, in the same relative orientation, and are flanked by 15-bp imperfect inverted repeats. The inserted elements interrupt the 25K open reading frame (ORF). One of these FP mutants undergoes spontaneous reversion. Sequence analysis at the excision site of a spontaneous revertant demonstrates that the TFP3 elements are capable of precise excision, restoring the expression of the 25-kDa protein. We also compare the sequences of the 25K genes of the Autographa californica and Galleria mellonella viruses (AcMNPV and GmMNPV, respectively). The 25K gene sequences diverge in five areas, resulting in an additional EcoRV and TaqI site within the GmMNPV 25K gene, and extension of the ORF for an additional 2 amino acids at the C-terminus of the predicted GmMNPV 25 kDa protein. The phenomenon of transposon mutagenesis of Baculovirus genomes provides a unique opportunity for analysis of transposon mobility.

Transposon mutagenesis of baculoviruses: analysis of Trichoplusia ni transposon IFP2 insertions within the FP-locus of nuclear polyhedrosis viruses

The transposable IFP2 element of Trichoplusia ni was originally isolated as a host DNA insertion in spontaneous FP mutants of Galleria mellonella or Autographa californica nuclear polyhedrosis viruses (NPVs). The termini of IFP2 insertions from five independently isolated FP mutants were sequenced. In all cases IFP2 is flanked by 13-bp terminal inverted repeats and has additional inverted repeats of 19 bp in length located asymmetrically with respect to the ends of the element. Insertion of IFP2 into the viral genome always generated a duplication of the tetranucleotide target site, TTAA. There was an apparent preference for insertion within a 12-bp A + T-rich imperfect palindromic sequence surrounding the target site. Sequence analysis of three independent IFP2 elements revealed an internal domain of 2.475 kb containing an RNA polymerase II promoter region and two large open reading frames. Primer extension analysis of IFP2-specific mRNA positioned the 5' terminus of the transcript. The element is present in DNA isolated from T. ni cell lines TN-368 and TN-5B1, but is not apparent in DNAs isolated from the TN-R2 cell line or our laboratory colony of T. ni larvae, suggesting IFP2 was recently introduced into the T. ni genome.

Composite transposable elements in the Xenopus laevis genome

Members of two related families of transposable elements, Tx1 and Tx2, were isolated from the genome of Xenopus laevis and characterized. In both families, two versions of the elements were found. The smaller version in each family (Tx1d and Tx2d) consisted largely of two types of 400-base-pair tandem internal repeats. These elements had discrete ends and short inverted terminal repeats characteristic of mobile DNAs that are presumed to move via DNA intermediates, e.g., Drosophila P and maize Ac elements. The longer versions (Tx1c and Tx2c) differed from Tx1d and Tx2d by the presence of a 6.9-kilobase-pair internal segment that included two long open reading frames (ORFs). ORF1 had one cysteine-plus-histidine-rich sequence of the type found in retroviral gag proteins. ORF2 showed more substantial homology to retroviral pol genes and particularly to the analogs of pol found in a subclass of mobile DNAs that are supposed retrotransposons, such as mammalian long interspersed repetitive sequences, Drosophila I factors, silkworm R1 elements, and trypanosome Ingi elements. Thus, the Tx1 elements present a paradox by exhibiting features of two classes of mobile DNAs that are thought to have very different modes of transposition. Two possible resolutions are considered: (i) the composite versions are actually made up of two independent elements, one of the retrotransposon class, which has a high degree of specificity for insertion into a target within the other, P-like element; and (ii) the composite elements are intact, autonomous mobile DNAs, in which the pol-like gene product collaborates with the terminal inverted repeats to cause transposition of the entire unit.

The Tc3 family of transposable genetic elements in Caenorhabditis elegans

We describe genetic and molecular properties of Tc3, a family of transposable elements in Caenorhabditis elegans. About 15 Tc3 elements are present in the genomes of several different wild-type varieties of C. elegans, but Tc3 transposition and excision are not detected in these strains. Tc3 transposition and excision occur at high frequencies, however, in strain TR679, a mutant identified because of its highly active Tc1 elements. In TR679, Tc3 is responsible for several spontaneous mutations affecting the unc-22 gene. Tc3-induced mutations are unstable, and revertants result from precise or nearly precise excision of Tc3. Although Tc3 is very active in TR679, it is not detectably active in several other mutator mutants, all of which exhibit high levels of Tc1 activity. Tc3 is 2.5 kilobases long, and except for sequences near its inverted repeat termini, it is unrelated to Tc1. The termini of Tc3 are inverted repeats of at least 70 base pairs; the terminal 8 nucleotides of Tc3 are identical to 8 of the terminal 9 nucleotides of Tc1.