MGEScan-non-LTR: computational identification and classification of autonomous non-LTR retrotransposons in eukaryotic genomes

Computational methods for genome-wide identification of mobile genetic elements (MGEs) have become increasingly necessary for both genome annotation and evolutionary studies. Non-long terminal repeat (non-LTR) retrotransposons are a class of MGEs that have been found in most eukaryotic genomes, sometimes in extremely high numbers. In this article, we present a computational tool, MGEScan-non-LTR, for the identification of non-LTR retrotransposons in genomic sequences, following a computational approach inspired by a generalized hidden Markov model (GHMM). Three different states represent two different protein domains and inter-domain linker regions encoded in the non-LTR retrotransposons, and their scores are evaluated by using profile hidden Markov models (for protein domains) and Gaussian Bayes classifiers (for linker regions), respectively. In order to classify the non-LTR retrotransposons into one of the 12 previously characterized clades using the same model, we defined separate states for different clades. MGEScan-non-LTR was tested on the genome sequences of four eukaryotic organisms, Drosophila melanogaster, Daphnia pulex, Ciona intestinalis and Strongylocentrotus purpuratus. For the D. melanogaster genome, MGEScan-non-LTR found all known ‘full-length’ elements and simultaneously classified them into the clades CR1, I, Jockey, LOA and R1. Notably, for the D. pulex genome, in which no non-LTR retrotransposon has been annotated, MGEScan-non-LTR found a significantly larger number of elements than did RepeatMasker, using the current version of the RepBase Update library. We also identified novel elements in the other two genomes, which have only been partially studied for non-LTR retrotransposons.     

APE-type non-LTR retrotransposons: determinants involved in target site recognition

Non-long terminal repeat (Non-LTR) retrotransposons represent a diverse and widely distributed group of transposable elements and an almost ubiquitous component of eukaryotic genomes that has a major impact on evolution. Their copy number can range from a few to several million and they often make up a significant fraction of the genomes. The members of the dominating subtype of non-LTR retrotransposons code for an endonuclease with homology to apurinic/apyrimidinic endonucleases (APE), and are thus termed APE-type non-LTR retrotransposons. In the last decade both the number of identified non-LTR retrotransposons and our knowledge of biology and evolution of APE-type non-LTR retrotransposons has increased tremendously.     

Non-LTR retrotransposons in fungi

Non-long terminal repeat (non-LTR) retrotransposons have contributed to shaping the structure and function of genomes. Fungi have small genomes, usually with limited amounts of repetitive DNA. In silico approach has been used to survey the non-LTR elements in 57 fungal genomes. More than 100 novel non-LTR retrotransposons were found, which belonged to five diverse clades. The present survey identified two novel clades of fungal non-LTR retrotransposons. The copy number of non-LTR retroelements varied widely. Some of the studied species contained a single copy of non-LTR retrotransposon, whereas others possessed a great number of non-LTR retrotransposon copies per genome. Although evolutionary relationships of most elements are congruent with phylogeny of host species, a new case of possible horizontal transfer was found between Eurotiomycetes and Sordariomycetes. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Getting closer to a pre-vertebrate genome: the non-LTR retrotransposons of Branchiostoma floridae

Non-LTR retrotransposons are common in vertebrate genomes and although present in invertebrates they appear at a much lower frequency. The cephalochordate amphioxus is the closest living relative to vertebrates and has been considered a good model for comparative analyses of genome expansions during vertebrate evolution. With the aim to assess the involvement of transposable elements in these events, we have analysed the non-LTR retrotransposons of Branchiostoma floridae. In silico searches have allowed to reconstruct non-LTR elements of six different clades (CR1, I, L1, L2, NeSL and RTE) and assess their structural features. According to the estimated copy number of these elements they account for less than 1% of the haploid genome, which reminds of the low abundance also encountered in the urochordate Ciona intestinalis. Amphioxus (B. floridae) and Ciona share a pre-vertebrate-like organization for the non-LTR retrotransposons (<150 copies, < 1% of the genome) versus the complexity associated to higher vertebrates (Homo sapiens >1.3.10(6) copies, > 20% of the genome).     

LINEs, SINEs and processed pseudogenes: parasitic strategies for genome modeling

Two major classes of retrotransposons have invaded eukaryotic genomes: the LTR retrotransposons closely resembling the proviral integrated form of infectious retroviruses, and the non-LTR retrotransposons including the widespread, autonomous LINE elements. Here, we review the modeling effects of the latter class of elements, which are the most active in humans, and whose enzymatic machinery is subverted to generate a large series of "secondary" retroelements. These include the processed pseudogenes, naturally present in all eukaryotic genomes possessing non-LTR retroelements, and the very successful SINE elements such as the human Alu sequences which have evolved refined parasitic strategies to efficiently bypass the original "protectionist" cis-preference of LINEs for their own retrotransposition.     

Retroelement dynamics and a novel type of chordate retrovirus-like element in the miniature genome of the tunicate Oikopleura dioica

Retrotransposable elements have played an important role in shaping eukaryotic DNA, and their activity and turnover rate directly influence the size of genomes. With approximately 15,000 genes within 65-75 megabases, the marine tunicate Oikopleura dioica, a nonvertebrate chordate, has the smallest and most compact genome ever found in animals. Consistent with a massive elimination of retroelements, only one apparently novel clade of non-long terminal repeat (non-LTR) retrotransposons was detected within 41 megabases of nonredundant genomic sequences. In contrast, at least six clades of non-LTR elements were identified in the less compact genome of the tunicate Ciona intestinalis. Unexpectedly, Ty3/gypsy-related Tor LTR retrotransposons presented an astonishing level of diversity in O. dioica. They were generally poorly or apparently not corrupted, indicating recent activity. Both Tor3 and Tor4b families bore an envelope-like open reading frame, suggesting possible horizontal acquisition through infection. The Tor4b envelope-like gene might have been obtained from a paramyxovirus (RNA virus). Tor3 and Tor4b are phylogenetically clearly distinct from vertebrate retroviruses (Retroviridae) and are more reminiscent of certain insect and plant sequences. Tor elements potentially represent a so far unknown, ancient type of infectious retroelement in chordates. Their distribution and transmission dynamics in tunicates and other chordates deserve further study.