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.     

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.     

The <Emphasis Type="Italic">Juan</Emphasis> non-LTR retrotransposon in mosquitoes: genomic impact,vertical transmission and indications of recent and widespread activity

Background  

In contrast to DNA-mediated transposable elements (TEs), retrotransposons, particularly non-long terminal repeat retrotransposons (non-LTRs), are generally considered to have a much lower propensity towards horizontal transfer. Detailed studies on site-specific non-LTR families have demonstrated strict vertical transmission. More studies are needed with non-site-specific non-LTR families to determine whether strict vertical transmission is a phenomenon related to site specificity or a more general characteristic of all non-LTRs. Juan is a Jockey clade non-LTR retrotransposon first discovered in mosquitoes that is widely distributed in the mosquito family Culicidae. Being a non-site specific non-LTR, Juan offers an opportunity to further investigate the hypothesis that non-LTRs are genomic elements that are primarily vertically transmitted.     

Low diversity,activity, and density of transposable elements in five avian genomes

In this study, we conducted the activity, diversity, and density analysis of transposable elements (TEs) across five avian genomes (budgerigar, chicken, turkey, medium ground finch, and zebra finch) to explore the potential reason of small genome sizes of birds. We found that these avian genomes exhibited low density of TEs by about 10% of genome coverages and low diversity of TEs with the TE landscapes dominated by CR1 and ERV elements, and contrasting proliferation dynamics both between TE types and between species were observed across the five avian genomes. Phylogenetic analysis revealed that CR1 clade was more diverse in the family structure compared with R2 clade in birds; avian ERVs were classified into four clades (alpha, beta, gamma, and ERV-L) and belonged to three classes of ERV with an uneven distributed in these lineages. The activities of DNA and SINE TEs were very low in the evolution history of avian genomes; most LINEs and LTRs were ancient copies with a substantial decrease of activity in recent, with only LTRs and LINEs in chicken and zebra finch exhibiting weak activity in very recent, and very few TEs were intact; however, the recent activity may be underestimated due to the sequencing/assembly technologies in some species. Overall, this study demonstrates low diversity, activity, and density of TEs in the five avian species; highlights the differences of TEs in these lineages; and suggests that the current and recent activity of TEs in avian genomes is very limited, which may be one of the reasons of small genome sizes in birds.     

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).     

NeSL-1, an ancient lineage of site-specific non-LTR retrotransposons from Caenorhabditis elegans

Phylogenetic analyses of non-LTR retrotransposons suggest that all elements can be divided into 11 lineages. The 3 oldest lineages show target site specificity for unique locations in the genome and encode an endonuclease with an active site similar to certain restriction enzymes. The more "modern" non-LTR lineages possess an apurinic endonuclease-like domain and generally lack site specificity. The genome sequence of Caenorhabditis elegans reveals the presence of a non-LTR retrotransposon that resembles the older elements, in that it contains a single open reading frame with a carboxyl-terminal restriction-like endonuclease domain. Located near the N-terminal end of the ORF is a cysteine protease domain not found in any other non-LTR element. The N2 strain of C. elegans appears to contain only one full-length and several 5' truncated copies of this element. The elements specifically insert in the Spliced leader-1 genes; hence the element has been named NeSL-1 (Nematode Spliced Leader-1). Phylogenetic analysis confirms that NeSL-1 branches very early in the non-LTR lineage and that it represents a 12th lineage of non-LTR elements. The target specificity of NeSL-1 for the spliced leader exons and the similarity of its structure to that of R2 elements leads to a simple model for its expression and retrotransposition.     

Identification of rDNA-specific non-LTR retrotransposons in Cnidaria

Ribosomal RNA genes are abundant repetitive sequences in most eukaryotes. Ribosomal DNA (rDNA) contains many insertions derived from mobile elements including non-long terminal repeat (non-LTR) retrotransposons. R2 is the well-characterized 28S rDNA-specific non-LTR retrotransposon family that is distributed over at least 4 bilaterian phyla. R2 is a large family sharing the same insertion specificity and classified into 4 clades (R2-A, -B, -C, and -D) based on the N-terminal domain structure and the phylogeny. There is no observation of horizontal transfer of R2; therefore, the origin of R2 dates back to before the split between protostomes and deuterostomes. Here, we in silico identified 1 R2 element from the sea anemone Nematostella vectensis and 2 R2-like retrotransposons from the hydrozoan Hydra magnipapillata. R2 from N. vectensis was inserted into the 28S rDNA like other R2, but the R2-like elements from H. magnipapillata were inserted into the specific sequence in the highly conserved region of the 18S rDNA. We designated the Hydra R2-like elements R8. R8 is inserted at 37 bp upstream from R7, another 18S rDNA-specific retrotransposon family. There is no obvious sequence similarity between targets of R2 and R8, probably because they recognize long DNA sequences. Domain structure and phylogeny indicate that R2 from N. vectensis is the member of the R2-D clade, and R8 from H. magnipapillata belongs to the R2-A clade despite its different sequence specificity. These results suggest that R2 had been generated before the split between cnidarians and bilaterians and that R8 is a retrotransposon family that changed its target from the 28S rDNA to the 18S rDNA.     

Relationships Among Powered Flight,Metabolic Rate,Body Mass,Genome Size,and the Retrotransposon Complement of Volant Birds

Avian genomes are of interest because the rapid metabolic rate associated with powered flight requires small cells which constrain genome size. Consequently, flying birds tend to have small genomes relative to other vertebrates such as mammals. It thus stands to reason that flying birds should have smaller genomes than ground-dwelling birds with lower metabolic rates. Small genomes could be condensed but uncompromised in a number of ways, including smaller intergenic intervals, shorter introns, and/or a reduced transposable element (TE) complement. We evaluated genome size in light of the orthologous TE complement among 41 flying (FY) and seven ground-dwelling (GD) bird species to determine if a preponderance of deletions in orthologous TEs might explain the compact genomes of flying birds with high metabolic rates. We measured, across multiple loci in all 48 species, the lengths of 50 contemporary orthologous chicken repeat 1 (CR1, a non-LTR retrotransposon) copies relative to inferred ancestral CR1 sequences. We found genome sizes in GD birds were not different than those in FY birds, but the mean lengths of orthologous CR1 loci were significantly shorter in FY birds than in GD birds. Moreover, we observed a negative correlation between basal metabolic rate and length of orthologous CR1 loci. Finally, we observed positive correlations between body mass and both genome sizes as well as length of orthologous CR1 loci, which we expected given that body mass correlates negatively with metabolic rates. Our results support the contention that metabolism helps shape the avian TE complement and thus indirectly contributes to the compact genomes of birds.     

Multiple tandem gene duplications in a neutral lipase gene cluster in Drosophila

We have examined a highly dynamic section of the Drosophila melanogaster genome which contains neutral lipase family genes that have undergone multiple tandem duplication events. We have identified the orthologous clusters, encoding between five and eight apparently functional lipases, in other Drosophila genomes: yakuba, ananassae, pseudoobscura, virilis, mojavensis, persimilis, grimshawi and willistoni. We examined their gene structure, duplication and pseudogene formation, and the presence of transposable elements. Based on phylogenetic comparisons, the lipase genes contained in each of the clusters fall into four distinct clades. Clades I and II have distinct evolutionary constraints to clades III and IV. Multiple gene duplications have occurred in different lineages of clades I and II while clades III and IV contain a single lipase gene from each species. Compared with lipases from other clades, clade IV genes contain an additional 3' domain of tandemly repeated sequence of varying length and composition, and a substitution in the residue adjacent to the key catalytic serine in the encoded proteins. A comparison of non-synonymous to synonymous nucleotide substitution (dN/dS) rates within each clade showed the highest rate of divergence was between paralogous lipase gene pairs suggesting selection pressure on duplicated genes. Analysis of the encoded lipase protein sequences within each species using PAML identified positively selected sites; structure homology modeling based on human pancreatic lipase indicated many of these residues formed part of the active site of the enzyme. As some of the cluster lipase genes are known to be expressed in the insect midgut and respond to changes in dietary components, we propose that the lipase cluster has undergone dynamic evolutionary changes to maximize absorption of lipid nutrients from the diet.     

The Non-LTR Retrotransposon R2 in Termites (Insecta,Isoptera): Characterization and Dynamics

The full-length element of the non-LTR retrotransposon R2 is here characterized in three European isopteran species: the more primitive Kalotermes flavicollis (Kalotermitidae), including two highly divergent mitochondrial lineages, and the more derived Reticulitermes lucifugus and R. urbis (Rhinotermitidae). Partial 3′ sequences for R. grassei and R. balkanensis were also analyzed. The essential structural features of R2 elements are conserved in termites. Phylogenetic analysis revealed that termite elements belong to the same clade and that their phylogeny is fully compatible with the phylogeny of their host species. The study of the number and the frequency of R2 insertion variants in four R. urbis colonies suggests a greatly reduced, or completely absent, recent element activity.     

Retrotransposon populations of <Emphasis Type="Italic">Vicia</Emphasis> species with varying genome size

The (non-LTR) LINE and Ty3-gypsy-type LTR retrotransposon populations of three Vicia species that differ in genome size (Vicia faba, Vicia melanops and Vicia sativa) have been characterised. In each species the LINE retrotransposons comprise a complex, very heterogeneous set of sequences, while the Ty3-gypsy elements are much more homogeneous. Copy numbers of all three retrotransposon groups (Ty1-copia, Ty3-gypsy and LINE) in these species have been estimated by random genomic sequencing and Southern hybridisation analysis. The Ty3-gypsy elements are extremely numerous in all species, accounting for 18–35% of their genomes. The Ty1-copia group elements are somewhat less abundant and LINE elements are present in still lower amounts. Collectively, 20–45% of the genomes of these three Vicia species are comprised of retrotransposons. These data show that the three retrotransposon groups have proliferated to different extents in members of the Vicia genus and high proliferation has been associated with homogenisation of the retrotransposon population.Electronic Supplementary Material Supplementary material is available for this article at .     

Teleost Fish Genomes Contain a Diverse Array of L1 RetrotransposonLineages That Exhibit a Low Copy Number and High Rate of Turnover

Retrotransposable elements exhibit a wide range of variation in population dynamics, abundance, and lineage diversity among host genomes across taxa. This range of diversity is illustrated by a single well-defined constituent monophyletic clade of L1 non-LTR retrotransposons that is shared between mammalian and teleost fish genomes. Despite the clear phylogenetic relationships that exist between mammalian and teleost L1 sequences, these elements exhibit markedly different dynamics within their respective taxa. While mammalian genomes typically contain a single, abundant lineage of L1 elements that traces millions of years of evolution, the zebraflsh genome was recently shown to exhibit a high diversity of ancient lineages coexisting at a very low copy number and apparently exhibiting a high rate of turnover. In the present study, a combination of degenerate PCR, lineage-specific PCR, and genomic Southern blot analysis is utilized to demonstrate high L1 lineage diversity, low copy number, and a high proportion of polymorphic inserts in the genomes of the killifish species, Fundulus heteroclitus. Additional species surveyed by degenerate PCR include Cyprinodon variegatus, Rivulus marmoratus, and Menidia beryllina. These results further support the generality of the differences that exist in host–element dynamics between teleost fish and mammalian genomes with regard to L1 retrotransposons.Reviewing Editor: Dr. Axel Meyer