Simple and fast classification of non-LTR retrotransposons based on phylogeny of their RT domain protein sequences

The Families of Non-LTR Transposable Elements within Neritimorpha and Other Gastropoda

Repeated sequences, especially transposable elements (TEs), are known to be abundant in some members of the important invertebrate class Gastropoda. TEs that do not have long terminal repeated sequences (non-LTR TEs) are frequently the most abundant type but have not been well characterised in any gastropod. Despite this, sequences in draft gastropod genomes are often described as non-LTR TEs, but without identification to family type. This study was conducted to characterise non-LTR TEs in neritimorph snails, using genomic skimming surveys of three species and the recently published draft genome of Theodoxus fluviatilis. Multiple families of non-LTR TEs from the I, Jockey, L1, R2 and RTE superfamilies were found, although there were notably few representatives of the first of these, which is nevertheless abundant in other Gastropoda. Phylogenetic analyses of amino acid sequences of the reverse transcriptase domain from the elements ORF2 regions found considerable interspersion of representatives of the four neritimorph taxa within non-LTR families and sub-families. In contrast, phylogenetic analyses of sequences from the elements' ORF1 region resolved the representatives from individual species as monophyletic. However, using either region, members of the two species of the Neritidae were closely related, suggesting their potential for investigation of phyletic evolution at the family level.

Transposable elements in Rosaceae: insights into genome evolution, expression dynamics, and syntenic gene regulation

Transposable elements (TEs) exert significant influence on plant genomic structure and gene expression. Here, we explored TE-related aspects across 14 Rosaceae genomes, investigating genomic distribution, transposition activity, expression patterns, and nearby differentially expressed genes (DEGs). Analyses unveiled distinct long terminal repeat retrotransposon (LTR-RT) evolutionary patterns, reflecting varied genome size changes among nine species over the past million years. In the past 2.5 million years, Rubus idaeus showed a transposition rate twice as fast as Fragaria vesca, while Pyrus bretschneideri displayed significantly faster transposition compared with Crataegus pinnatifida. Genes adjacent to recent TE insertions were linked to adversity resistance, while those near previous insertions were functionally enriched in morphogenesis, enzyme activity, and metabolic processes. Expression analysis revealed diverse responses of LTR-RTs to internal or external conditions. Furthermore, we identified 3695 pairs of syntenic DEGs proximal to TEs in Malus domestica cv. 'Gala' and M. domestica (GDDH13), suggesting TE insertions may contribute to varietal trait differences in these apple varieties. Our study across representative Rosaceae species underscores the pivotal role of TEs in plant genome evolution within this diverse family. It elucidates how these elements regulate syntenic DEGs on a genome-wide scale, offering insights into Rosaceae-specific genomic evolution.

Conserved and divergent DNA recognition specificities and functions of R2 retrotransposon N-terminal domains

R2 non-long terminal repeat (non-LTR) retrotransposons are among the most extensively distributed mobile genetic elements in multicellular eukaryotes and show promise for applications in transgene supplementation of the human genome. They insert new gene copies into a conserved site in 28S ribosomal DNA with exquisite specificity. R2 clades are defined by the number of zinc fingers (ZFs) at the N terminus of the retrotransposon-encoded protein, postulated to additively confer DNA site specificity. Here, we illuminate general principles of DNA recognition by R2 N-terminal domains across and between clades, with extensive, specific recognition requiring only one or two compact domains. DNA-binding and protection assays demonstrate broadly shared as well as clade-specific DNA interactions. Gene insertion assays in cells identify the N-terminal domains sufficient for target-site insertion and reveal roles in second-strand cleavage or synthesis for clade-specific ZFs. Our results have implications for understanding evolutionary diversification of non-LTR retrotransposon insertion mechanisms and the design of retrotransposon-based gene therapies.

The good, the bad and the ugly of transposable elements annotation tools

Transposable elements are repetitive and mobile DNA segments that can be found in virtually all organisms investigated to date. Their complex structure and variable nature are particularly challenging from the genomic annotation point of view. Many softwares have been developed to automate and facilitate TEs annotation at the genomic level, but they are highly heterogeneous regarding documentation, usability and methods. In this review, we revisited the existing software for TE genomic annotation, concentrating on the most often used ones, the methodologies they apply, and usability. Building on the state of the art of TE annotation software we propose best practices and highlight the strengths and weaknesses from the available solutions.

Human Endogenous Retrovirus-K (HML-2)-Related Genetic Variation: Human Genome Diversity and Disease

Human endogenous retroviruses (HERVs) comprise a significant portion of the human genome, making up roughly 8%, a notable comparison to the 2-3% represented by coding sequences. Numerous studies have underscored the critical role and importance of HERVs, highlighting their diverse and extensive influence on the evolution of the human genome and establishing their complex correlation with various diseases. Among HERVs, the HERV-K (HML-2) subfamily has recently attracted significant attention, integrating into the human genome after the divergence between humans and chimpanzees. Its insertion in the human genome has received considerable attention due to its structural and functional characteristics and the time of insertion. Originating from ancient exogenous retroviruses, these elements succeeded in infecting germ cells, enabling vertical transmission and existing as proviruses within the genome. Remarkably, these sequences have retained the capacity to form complete viral sequences, exhibiting activity in transcription and translation. The HERV-K (HML-2) subfamily is the subject of active debate about its potential positive or negative effects on human genome evolution and various pathologies. This review summarizes the variation, regulation, and diseases in human genome evolution arising from the influence of HERV-K (HML-2).

Multiple and diversified transposon lineages contribute to early and recent bivalve genome evolution

BACKGROUND

Transposable elements (TEs) can represent one of the major sources of genomic variation across eukaryotes, providing novel raw materials for species diversification and innovation. While considerable effort has been made to study their evolutionary dynamics across multiple animal clades, molluscs represent a substantially understudied phylum. Here, we take advantage of the recent increase in mollusc genomic resources and adopt an automated TE annotation pipeline combined with a phylogenetic tree-based classification, as well as extensive manual curation efforts, to characterize TE repertories across 27 bivalve genomes with a particular emphasis on DDE/D class II elements, long interspersed nuclear elements (LINEs), and their evolutionary dynamics.

RESULTS

We found class I elements as highly dominant in bivalve genomes, with LINE elements, despite less represented in terms of copy number per genome, being the most common retroposon group covering up to 10% of their genome. We mined 86,488 reverse transcriptases (RVT) containing LINE coming from 12 clades distributed across all known superfamilies and 14,275 class II DDE/D-containing transposons coming from 16 distinct superfamilies. We uncovered a previously underestimated rich and diverse bivalve ancestral transposon complement that could be traced back to their most recent common ancestor that lived ~ 500 Mya. Moreover, we identified multiple instances of lineage-specific emergence and loss of different LINEs and DDE/D lineages with the interesting cases of CR1- Zenon, Proto2, RTE-X, and Academ elements that underwent a bivalve-specific amplification likely associated with their diversification. Finally, we found that this LINE diversity is maintained in extant species by an equally diverse set of long-living and potentially active elements, as suggested by their evolutionary history and transcription profiles in both male and female gonads.

CONCLUSIONS

We found that bivalves host an exceptional diversity of transposons compared to other molluscs. Their LINE complement could mainly follow a "stealth drivers" model of evolution where multiple and diversified families are able to survive and co-exist for a long period of time in the host genome, potentially shaping both recent and early phases of bivalve genome evolution and diversification. Overall, we provide not only the first comparative study of TE evolutionary dynamics in a large but understudied phylum such as Mollusca, but also a reference library for ORF-containing class II DDE/D and LINE elements, which represents an important genomic resource for their identification and characterization in novel genomes.

How to start a LINE: 5' switching rejuvenates LINE retrotransposons in tobacco and related Nicotiana species

By contrast to their conserved mammalian counterparts, plant long interspersed nuclear elements (LINEs) are highly variable, splitting into many low-copy families. Curiously, LINE families from the retrotransposable element (RTE) clade retain a stronger sequence conservation and hence reach higher copy numbers. The cause of this RTE-typical property is not yet understood, but would help clarify why some transposable elements are removed quickly, whereas others persist in plant genomes. Here, we bring forward a detailed study of RTE LINE structure, diversity and evolution in plants. For this, we argue that the nightshade family is the ideal taxon to follow the evolutionary trajectories of RTE LINEs, given their high abundance, recent activity and partnership to non-autonomous elements. Using bioinformatic, cytogenetic and molecular approaches, we detect 4029 full-length RTE LINEs across the Solanaceae. We finely characterize and manually curate a core group of 458 full-length LINEs in allotetraploid tobacco, show an integration event after polyploidization and trace hybridization by RTE LINE composition of parental genomes. Finally, we reveal the role of the untranslated regions (UTRs) as causes for the unique RTE LINE amplification and evolution pattern in plants. On the one hand, we detected a highly conserved motif at the 3' UTR, suggesting strong selective constraints acting on the RTE terminus. On the other hand, we observed successive rounds of 5' UTR cycling, constantly rejuvenating the promoter sequences. This interplay between exchangeable promoters and conserved LINE bodies and 3' UTR likely allows RTE LINEs to persist and thrive in plant genomes.

Paleozoic Protein Fossils Illuminate the Evolution of Vertebrate Genomes and Transposable Elements

Genomes hold a treasure trove of protein fossils: fragments of formerly protein-coding DNA, which mainly come from transposable elements (TEs) or host genes. These fossils reveal ancient evolution of TEs and genomes, and many fossils have been exapted to perform diverse functions important for the host's fitness. However, old and highly-degraded fossils are hard to identify, standard methods (e.g. BLAST) are not optimized for this task, and few Paleozoic protein fossils have been found. Here, a recently optimized method is used to find protein fossils in vertebrate genomes. It finds Paleozoic fossils predating the amphibian/amniote divergence from most major TE categories, including virus-related Polinton and Gypsy elements. It finds 10 fossils in the human genome (8 from TEs and 2 from host genes) that predate the last common ancestor of all jawed vertebrates, probably from the Ordovician period. It also finds types of transposon and retrotransposon not found in human before. These fossils have extreme sequence conservation, indicating exaptation: some have evidence of gene-regulatory function, and they tend to lienearest to developmental genes. Some ancient fossils suggest "genome tectonics", where two fragments of one TE have drifted apart by up to megabases, possibly explaining gene deserts and large introns. This paints a picture of great TE diversity in our aquatic ancestors, with patchy TE inheritance by later vertebrates, producing new genes and regulatory elements on the way. Host-gene fossils too have contributed anciently-conserved DNA segments. This paves the way to further studies of ancient protein fossils.

Migrators within migrators: exploring transposable element dynamics in the monarch butterfly, Danaus plexippus

Background

Lepidoptera (butterflies and moths) are an important model system in ecology and evolution. A high-quality chromosomal genome assembly is available for the monarch butterfly (Danaus plexippus), but it lacks an in-depth transposable element (TE) annotation, presenting an opportunity to explore monarch TE dynamics and the impact of TEs on shaping the monarch genome.

Results

We find 6.21% of the monarch genome is comprised of TEs, a reduction of 6.85% compared to the original TE annotation performed on the draft genome assembly. Monarch TE content is low compared to two closely related species with available genomes, Danaus chrysippus (33.97% TE) and Danaus melanippus (11.87% TE). The biggest TE contributions to genome size in the monarch are LINEs and Penelope-like elements, and three newly identified families, r2-hero_dPle (LINE), penelope-1_dPle (Penelope-like), and hase2-1_dPle (SINE), collectively contribute 34.92% of total TE content. We find evidence of recent TE activity, with two novel Tc1 families rapidly expanding over recent timescales (tc1-1_dPle, tc1-2_dPle). LINE fragments show signatures of genomic deletions indicating a high rate of TE turnover. We investigate associations between TEs and wing colouration and immune genes and identify a three-fold increase in TE content around immune genes compared to other host genes.

Conclusions

We provide a detailed TE annotation and analysis for the monarch genome, revealing a considerably smaller TE contribution to genome content compared to two closely related Danaus species with available genome assemblies. We identify highly successful novel DNA TE families rapidly expanding over recent timescales, and ongoing signatures of both TE expansion and removal highlight the dynamic nature of repeat content in the monarch genome. Our findings also suggest that insect immune genes are promising candidates for future interrogation of TE-mediated host adaptation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-022-00263-5.

Horizontal Transposon Transfer and Its Implications for the Ancestral Ecology of Hydrophiine Snakes

Transposable elements (TEs), also known as jumping genes, are sequences able to move or copy themselves within a genome. As TEs move throughout genomes they often act as a source of genetic novelty, hence understanding TE evolution within lineages may help in understanding environmental adaptation. Studies into the TE content of lineages of mammals such as bats have uncovered horizontal transposon transfer (HTT) into these lineages, with squamates often also containing the same TEs. Despite the repeated finding of HTT into squamates, little comparative research has examined the evolution of TEs within squamates. Here we examine a diverse family of Australo-Melanesian snakes (Hydrophiinae) to examine if the previously identified, order-wide pattern of variable TE content and activity holds true on a smaller scale. Hydrophiinae diverged from Asian elapids ~30 Mya and have since rapidly diversified into six amphibious, ~60 marine and ~100 terrestrial species that fill a broad range of ecological niches. We find TE diversity and expansion differs between hydrophiines and their Asian relatives and identify multiple HTTs into Hydrophiinae, including three likely transferred into the ancestral hydrophiine from fish. These HTT events provide the first tangible evidence that Hydrophiinae reached Australia from Asia via a marine route.

A genomic survey of LINE elements in Pipidae aquatic frogs shed light on Rex-elements evolution in these genomes

The transposable elements (TE) represent a large portion of anuran genomes that act as components of genetic diversification. The LINE order of retrotransposons is among the most representative and diverse TEs and is poorly investigated in anurans. Here we explored the LINE diversity with an emphasis on the elements generically called Rex in Pipidae species, more specifically, in the genomes ofXenopus tropicalis, used as a model genome in the study of anurans,the allotetraploid sister species Xenopus laevis and theAmerican species Pipa carvalhoi. We were able to identify a great diversity of LINEs from five clades, Rex1, L2, CR1, L1 and Tx1, in these three species, and the RTE clade was lost in X. tropicalis. It is clear that elements classified as Rex are distributed in distinct clades. The evolutionary pattern of Rex1 elements denote a complex evolution with independent losses of families and some horizontal transfer events between fishes and amphibians which were supported not only by the phylogenetic inconsistencies but also by the very low Ks values found for the TE sequences. The data obtained here update the knowledge of the LINEs diversity in X. laevis and represent the first study of TEs in P. carvalhoi.

Comparative genomic and transcriptomic analyses of transposable elements in polychaetous annelids highlight LTR retrotransposon diversity and evolution

Background

With the expansion of high throughput sequencing, we now have access to a larger number of genome-wide studies analyzing the Transposable elements (TEs) composition in a wide variety of organisms. However, genomic analyses often remain too limited in number and diversity of species investigated to study in depth the dynamics and evolutionary success of the different types of TEs among metazoans. Therefore, we chose to investigate the use of transcriptomes to describe the diversity of TEs in phylogenetically related species by conducting the first comparative analysis of TEs in two groups of polychaetes and evaluate the diversity of TEs that might impact genomic evolution as a result of their mobility.

Results

We present a detailed analysis of TEs distribution in transcriptomes extracted from 15 polychaetes depending on the number of reads used during assembly, and also compare these results with additional TE scans on associated low-coverage genomes. We then characterized the clades defined by 1021 LTR-retrotransposon families identified in 26 species. Clade richness was highly dependent on the considered superfamily. Copia elements appear rare and are equally distributed in only three clades, GalEa, Hydra and CoMol. Among the eight BEL/Pao clades identified in annelids, two small clades within the Sailor lineage are new for science. We characterized 17 Gypsy clades of which only 4 are new; the C-clade largely dominates with a quarter of the families. Finally, all species also expressed for the majority two distinct transcripts encoding PIWI proteins, known to be involved in control of TEs mobilities.

Conclusions

This study shows that the use of transcriptomes assembled from 40 million reads was sufficient to access to the diversity and proportion of the transposable elements compared to those obtained by low coverage sequencing. Among LTR-retrotransposons Gypsy elements were unequivocally dominant but results suggest that the number of Gypsy clades, although high, may be more limited than previously thought in metazoans. For BEL/Pao elements, the organization of clades within the Sailor lineage appears more difficult to establish clearly. The Copia elements remain rare and result from the evolutionary consistent success of the same three clades.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-021-00252-0.

Differential expression of transposable elements in the medaka melanoma model

Malignant melanoma incidence is rising worldwide. Its treatment in an advanced state is difficult, and the prognosis of this severe disease is still very poor. One major source of these difficulties is the high rate of metastasis and increased genomic instability leading to a high mutation rate and the development of resistance against therapeutic approaches. Here we investigate as one source of genomic instability the contribution of activation of transposable elements (TEs) within the tumor. We used the well-established medaka melanoma model and RNA-sequencing to investigate the differential expression of TEs in wildtype and transgenic fish carrying melanoma. We constructed a medaka-specific TE sequence library and identified TE sequences that were specifically upregulated in tumors. Validation by qRT- PCR confirmed a specific upregulation of a LINE and an LTR element in malignant melanomas of transgenic fish.

The function of twister ribozyme variants in non-LTR retrotransposition in Schistosoma mansoni

The twister ribozyme is widely distributed over numerous organisms and is especially abundant in Schistosoma mansoni, but has no confirmed biological function. Of the 17 non-LTR retrotransposons known in S. mansoni, none have thus far been associated with ribozymes. Here we report the identification of novel twister variant (T-variant) ribozymes and their function in S. mansoni non-LTR retrotransposition. We show that T-variant ribozymes are located at the 5′ end of Perere-3 non-LTR retrotransposons in the S. mansoni genome. T-variant ribozymes were demonstrated to be catalytically active in vitro. In reporter constructs, T-variants were shown to cleave in vivo, and cleavage of T-variants was sufficient for the translation of downstream reporter genes. Our analysis shows that the T-variants and Perere-3 are transcribed together. Target site duplications (TSDs); markers of target-primed reverse transcription (TPRT) and footmarks of retrotransposition, are located adjacent to the T-variant cleavage site and suggest that T-variant cleavage has taken place inS. mansoni. Sequence heterogeneity in the TSDs indicates that Perere-3 retrotransposition is not site-specific. The TSD sequences contribute to the 5′ end of the terminal ribozyme helix (P1 stem). Based on these results we conclude that T-variants have a functional role in Perere-3 retrotransposition.

Retrotransposable Elements: DNA Fingerprinting and the Assessment of Genetic Diversity

Retrotransposable elements (RTEs) are highly common mobile genetic elements that are composed of several classes and make up the majority of eukaryotic genomes. The "copy-out and paste-in" life cycle of replicative transposition in these dispersive and ubiquitous RTEs leads to new genome insertions without excision of the original element. RTEs are important drivers of species diversity; they exhibit great variety in structure, size, and mechanisms of transposition, making them important putative components in genome evolution. Accordingly, various applications have been developed to explore the polymorphisms in RTE insertion patterns. These applications include conventional or anchored polymerase chain reaction (PCR) and quantitative or digital PCR with primers designed for the 5' or 3' junction. Marker systems exploiting these PCR methods can be easily developed and are inexpensively used in the absence of extensive genome sequence data. The main inter-repeat amplification polymorphism techniques include inter-retrotransposon amplified polymorphism (IRAP), retrotransposon microsatellite amplified polymorphism (REMAP), and Inter-Primer Binding Site (iPBS) for PCR amplification with a single or two primers.

Existence of Bov-B LINE Retrotransposons in Snake Lineages Reveals Recent Multiple Horizontal Gene Transfers with Copy Number Variation

Transposable elements (TEs) are dynamic elements present in all eukaryotic genomes. They can “jump” and amplify within the genome and promote segmental genome rearrangements on both autosomes and sex chromosomes by disruption of gene structures. The Bovine-B long interspersed nuclear element (Bov-B LINE) is among the most abundant TE-retrotransposon families in vertebrates due to horizontal transfer (HT) among vertebrate lineages. Recent studies have shown multiple HTs or the presence of diverse Bov-B LINE groups in the snake lineage. It is hypothesized that Bov-B LINEs are highly dynamic and that the diversity reflects multiple HTs in snake lineages. Partial sequences of Bov-B LINE from 23 snake species were characterized. Phylogenetic analysis resolved at least two Bov-B LINE groups that might correspond to henophidian and caenophidian snakes; however, the tree topology differed from that based on functional nuclear and mitochondrial gene sequences. Several Bov-B LINEs of snakes showed greater than 80% similarity to sequences obtained from insects, whereas the two Bov-B LINE groups as well as sequences from the same snake species classified in different Bov-B LINE groups showed sequence similarities of less than 80%. Calculation of estimated divergence time and pairwise divergence between all individual Bov-B LINE copies suggest invasion times ranging from 79.19 to 98.8 million years ago in snakes. Accumulation of elements in a lineage-specific fashion ranged from 9 × 10⁻⁶% to 5.63 × 10⁻²% per genome. The genomic proportion of Bov-B LINEs varied among snake species but was not directly associated with genome size or invasion time. No differentiation in Bov-B LINE copy number between males and females was observed in any of the snake species examined. Incongruence in tree topology between Bov-B LINEs and other snake phylogenies may reflect past HT events. Sequence divergence of Bov-B LINEs between copies suggests that recent multiple HTs occurred within the same evolutionary timeframe in the snake lineage. The proportion of Bov-B LINEs varies among species, reflecting species specificity in TE invasion. The rapid speciation of snakes, coinciding with Bov-B LINE invasion in snake genomes, leads us to better understand the effect of Bov-B LINEs on snake genome evolution.

R2 and Non-Site-Specific R2-Like Retrotransposons of the German Cockroach, Blattella germanica

The structural and functional organization of the ribosomal RNA gene cluster and the full-length R2 non-LTR retrotransposon (integrated into a specific site of 28S ribosomal RNA genes) of the German cockroach, Blattella germanica, is described. A partial sequence of the R2 retrotransposon of the cockroach Rhyparobia maderae is also analyzed. The analysis of previously published next-generation sequencing data from the B. germanica genome reveals a new type of retrotransposon closely related to R2 retrotransposons but with a random distribution in the genome. Phylogenetic analysis reveals that these newly described retrotransposons form a separate clade. It is shown that proteins corresponding to the open reading frames of newly described retrotransposons exhibit unequal structural domains. Within these retrotransposons, a recombination event is described. New mechanism of transposition activity is discussed. The essential structural features of R2 retrotransposons are conserved in cockroaches and are typical of previously described R2 retrotransposons. However, the investigation of the number and frequency of 5′-truncated R2 retrotransposon insertion variants in eight B. germanica populations suggests recent mobile element activity. It is shown that the pattern of 5′-truncated R2 retrotransposon copies can be an informative molecular genetic marker for revealing genetic distances between insect populations.

Comparative genomic analyses illuminate the distinct evolution of megabats within Chiroptera

The revision of the sub-order Microchiroptera is one of the most intriguing outcomes in recent mammalian molecular phylogeny. The unexpected sister–taxon relationship between rhinolophoid microbats and megabats, with the exclusion of other microbats, suggests that megabats arose in a relatively short period of time from a microbat-like ancestor. In order to understand the genetic mechanism underlying adaptive evolution in megabats, we determined the whole-genome sequences of two rousette megabats, Leschenault’s rousette (Rousettus leschenaultia) and the Egyptian fruit bat (R. aegyptiacus). The sequences were compared with those of 22 other mammals, including nine bats, available in the database. We identified that megabat genomes are distinct in that they have extremely low activity of SINE retrotranspositions, expansion of two chemosensory gene families, including the trace amine receptor (TAAR) and olfactory receptor (OR), and elevation of the dN/dS ratio in genes for immunity and protein catabolism. The adaptive signatures discovered in the genomes of megabats may provide crucial insight into their distinct evolution, including key processes such as virus resistance, loss of echolocation, and frugivorous feeding.

Variation in base composition underlies functional and evolutionary divergence in non-LTR retrotransposons

Background

Non-LTR retrotransposons often exhibit base composition that is markedly different from the nucleotide content of their host’s gene. For instance, the mammalian L1 element is AT-rich with a strong A bias on the positive strand, which results in a reduced transcription. It is plausible that the A-richness of mammalian L1 is a self-regulatory mechanism reflecting a trade-off between transposition efficiency and the deleterious effect of L1 on its host. We examined if the A-richness of L1 is a general feature of non-LTR retrotransposons or if different clades of elements have evolved different nucleotide content. We also investigated if elements belonging to the same clade evolved towards different base composition in different genomes or if elements from different clades evolved towards similar base composition in the same genome.

Results

We found that non-LTR retrotransposons differ in base composition among clades within the same host but also that elements belonging to the same clade differ in base composition among hosts. We showed that nucleotide content remains constant within the same host over extended period of evolutionary time, despite mutational patterns that should drive nucleotide content away from the observed base composition.

Conclusions

Our results suggest that base composition is evolving under selection and may be reflective of the long-term co-evolution between non-LTR retrotransposons and their host. Finally, the coexistence of elements with drastically different base composition suggests that these elements may be using different strategies to persist and multiply in the genome of their host.

Conversion of DNA Sequences: From a Transposable Element to a Tandem Repeat or to a Gene

Eukaryotic genomes are rich in repetitive DNA sequences grouped in two classes regarding their genomic organization: tandem repeats and dispersed repeats. In tandem repeats, copies of a short DNA sequence are positioned one after another within the genome, while in dispersed repeats, these copies are randomly distributed. In this review we provide evidence that both tandem and dispersed repeats can have a similar organization, which leads us to suggest an update to their classification based on the sequence features, concretely regarding the presence or absence of retrotransposons/transposon specific domains. In addition, we analyze several studies that show that a repetitive element can be remodeled into repetitive non-coding or coding sequences, suggesting (1) an evolutionary relationship among DNA sequences, and (2) that the evolution of the genomes involved frequent repetitive sequence reshuffling, a process that we have designated as a "DNA remodeling mechanism". The alternative classification of the repetitive DNA sequences here proposed will provide a novel theoretical framework that recognizes the importance of DNA remodeling for the evolution and plasticity of eukaryotic genomes.

MGERT: a pipeline to retrieve coding sequences of mobile genetic elements from genome assemblies

Background

Genomes of eukaryotes are inhabited by myriads of mobile genetic elements (MGEs) – transposons and retrotransposons - which play a great role in genome plasticity and evolution. A lot of computational tools were developed to annotate them either in genomic assemblies or raw reads using de novo or homology-based approaches. But there has been no pipeline enabling users to get coding and flanking sequences of MGEs suitable for a downstream analysis from genome assemblies.

Results

We developed a new pipeline, MGERT (Mobile Genetic Elements Retrieving Tool), that automates all the steps necessary to obtain protein-coding sequences of mobile genetic elements from genomic assemblies even if no previous knowledge on MGE content of a particular genome is available.

Conclusions

Using MGERT, researchers can easily find MGEs, their coding and flanking sequences in the genome of interest. Thus, this pipeline helps researchers to focus on the biological analysis of MGEs rather than excessive scripting and pipelining.

Electronic supplementary material

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

LINEs Contribute to the Origins of Middle Bodies of SINEs besides 3' Tails

Short interspersed elements (SINEs), which are nonautonomous transposable elements, require the transposition machinery of long interspersed elements (LINEs) to mobilize. SINEs are composed of two or more independently originating parts. The 5′ region is called the “head” and is derived mainly from small RNAs, and the 3′ region (“tail”) originates from the 3′ region of LINEs and is responsible for being recognized by counterpart LINE proteins. The origin of the middle “body” of SINEs is enigmatic, although significant sequence similarities among SINEs from very diverse species have been observed. Here, a systematic analysis of the similarities among SINEs and LINEs deposited on Repbase, a comprehensive database of eukaryotic repeat sequences was performed. Three primary findings are described: 1) The 5′ regions of only two clades of LINEs, RTE and Vingi, were revealed to have contributed to the middle parts of SINEs; 2) The linkage of the 5′ and 3′ parts of LINEs can be lost due to occasional tail exchange of SINEs; and 3) The previously proposed Ceph-domain was revealed to be a fusion of a CORE-domain and a 5′ part of RTE clade of LINE. Based on these findings, a hypothesis that the 5′ parts of bipartite nonautonomous LINEs, which possess only the 5′ and 3′ regions of the original LINEs, can contribute to the undefined middle part of SINEs is proposed.

Differential Effect of Selection against LINE Retrotransposons among Vertebrates Inferred from Whole-Genome Data and Demographic Modeling

Variation in LINE composition is one of the major determinants for the substantial size and structural differences among vertebrate genomes. In particular, the larger genomes of mammals are characterized by hundreds of thousands of copies from a single LINE clade, L1, whereas nonmammalian vertebrates possess a much greater diversity of LINEs, yet with orders of magnitude less in copy number. It has been proposed that such variation in copy number among vertebrates is due to differential effect of LINE insertions on host fitness. To investigate LINE selection, we deployed a framework of demographic modeling, coalescent simulations, and probabilistic inference against population-level whole-genome data sets for four model species: one population each of threespine stickleback, green anole, and house mouse, as well as three human populations. Specifically, we inferred a null demographic background utilizing SNP data, which was then exploited to simulate a putative null distribution of summary statistics that was compared with LINE data. Subsequently, we applied the inferred null demographic model with an additional exponential size change parameter, coupled with model selection, to test for neutrality as well as estimate the strength of either negative or positive selection. We found a robust signal for purifying selection in anole and mouse, but a lack of clear evidence for selection in stickleback and human. Overall, we demonstrated LINE insertion dynamics that are not in accordance to a mammalian versus nonmammalian dichotomy, and instead the degree of existing LINE activity together with host-specific demographic history may be the main determinants of LINE abundance.

Using bioinformatic and phylogenetic approaches to classify transposable elements and understand their complex evolutionary histories

In recent years, much attention has been paid to comparative genomic studies of transposable elements (TEs) and the ensuing problems of their identification, classification, and annotation. Different approaches and diverse automated pipelines are being used to catalogue and categorize mobile genetic elements in the ever-increasing number of prokaryotic and eukaryotic genomes, with little or no connectivity between different domains of life. Here, an overview of the current picture of TE classification and evolutionary relationships is presented, updating the diversity of TE types uncovered in sequenced genomes. A tripartite TE classification scheme is proposed to account for their replicative, integrative, and structural components, and the need to expand in vitro and in vivo studies of their structural and biological properties is emphasized. Bioinformatic studies have now become front and center of novel TE discovery, and experimental pursuits of these discoveries hold great promise for both basic and applied science.

Identification and characterization of mobile genetic elements LINEs from Brassica genome

Among transposable elements (TEs), the LTR retrotransposons are abundant followed by non-LTR retrotransposons in plant genomes, the lateral being represented by LINEs and SINEs. Computational and molecular approaches were used for the characterization of Brassica LINEs, their diversity and phylogenetic relationships. Four autonomous and four non-autonomous LINE families were identified and characterized from Brassica. Most of the autonomous LINEs displayed two open reading frames, ORF1 and ORF2, where ORF1 is a gag protein domain, while ORF2 encodes endonuclease (EN) and a reverse transcriptase (RT). Three of four families encoded an additional RNase H (RH) domain in pol gene common to 'R' and 'I' type of LINEs. The PCR analyses based on LINEs RT fragments indicate their high diversity and widespread occurrence in tested 40 Brassica cultivars. Database searches revealed the homology in LINE sequences in closely related genera Arabidopsis indicating their origin from common ancestors predating their separation. The alignment of 58 LINEs RT sequences from Brassica, Arabidopsis and other plants depicted 4 conserved domains (domain II-V) showing similarity to previously detected domains. Based on RT alignment of Brassica and 3 known LINEs from monocots, Brassicaceae LINEs clustered in separate clade, further resolving 4 Brassica-Arabidopsis specific families in 2 sub-clades. High similarities were observed in RT sequences in the members of same family, while low homology was detected in members across the families. The investigation led to the characterization of Brassica specific LINE families and their diversity across Brassica species and their cultivars.

Subgenome Dominance in an Interspecific Hybrid, Synthetic Allopolyploid, and a 140-Year-Old Naturally Established Neo-Allopolyploid Monkeyflower

Recent studies have shown that one of the parental subgenomes in ancient polyploids is generally more dominant, having retained more genes and being more highly expressed, a phenomenon termed subgenome dominance. The genomic features that determine how quickly and which subgenome dominates within a newly formed polyploid remain poorly understood. To investigate the rate of emergence of subgenome dominance, we examined gene expression, gene methylation, and transposable element (TE) methylation in a natural, <140-year-old allopolyploid (Mimulus peregrinus), a resynthesized interspecies triploid hybrid (M. robertsii), a resynthesized allopolyploid (M. peregrinus), and progenitor species (M. guttatus and M. luteus). We show that subgenome expression dominance occurs instantly following the hybridization of divergent genomes and significantly increases over generations. Additionally, CHH methylation levels are reduced in regions near genes and within TEs in the first-generation hybrid, intermediate in the resynthesized allopolyploid, and are repatterned differently between the dominant and recessive subgenomes in the natural allopolyploid. Subgenome differences in levels of TE methylation mirror the increase in expression bias observed over the generations following hybridization. These findings provide important insights into genomic and epigenomic shock that occurs following hybridization and polyploid events and may also contribute to uncovering the mechanistic basis of heterosis and subgenome dominance.

Transposable elements in Drosophila

Transposable elements (TEs) are mobile genetic elements that can mobilize within host genomes. As TEs comprise more than 40% of the human genome and are linked to numerous diseases, understanding their mechanisms of mobilization and regulation is important. Drosophila melanogaster is an ideal model organism for the study of eukaryotic TEs as its genome contains a diverse array of active TEs. TEs universally impact host genome size via transposition and deletion events, but may also adopt unique functional roles in host organisms. There are 2 main classes of TEs: DNA transposons and retrotransposons. These classes are further divided into subgroups of TEs with unique structural and functional characteristics, demonstrating the significant variability among these elements. Despite this variability, D. melanogaster and other eukaryotic organisms utilize conserved mechanisms to regulate TEs. This review focuses on the transposition mechanisms and regulatory pathways of TEs, and their functional roles in D. melanogaster.

LINE Insertion Polymorphisms are Abundant but at Low Frequencies across Populations of Anolis carolinensis

Vertebrate genomes differ considerably in size and structure. Among the features that show the most variation is the abundance of Long Interspersed Nuclear Elements (LINEs). Mammalian genomes contain 100,000s LINEs that belong to a single clade, L1, and in most species a single family is usually active at a time. In contrast, non-mammalian vertebrates (fish, amphibians and reptiles) contain multiple active families, belonging to several clades, but each of them is represented by a small number of recently inserted copies. It is unclear why vertebrate genomes harbor such drastic differences in LINE composition. To address this issue, we conducted whole genome resequencing to investigate the population genomics of LINEs across 13 genomes of the lizard Anolis carolinensis sampled from two geographically and genetically distinct populations in the Eastern Florida and the Gulf Atlantic regions of the United States. We used the Mobile Element Locator Tool to identify and genotype polymorphic insertions from five major clades of LINEs (CR1, L1, L2, RTE and R4) and the 41 subfamilies that constitute them. Across these groups we found large variation in the frequency of polymorphic insertions and the observed length distributions of these insertions, suggesting these groups vary in their activity and how frequently they successfully generate full-length, potentially active copies. Though we found an abundance of polymorphic insertions (over 45,000) most of these were observed at low frequencies and typically appeared as singletons. Site frequency spectra for most LINEs showed a significant shift toward low frequency alleles compared to the spectra observed for total genomic single nucleotide polymorphisms. Using Tajima's D, FST and the mean number of pairwise differences in LINE insertion polymorphisms, we found evidence that negative selection is acting on LINE families in a length-dependent manner, its effects being stronger in the larger Eastern Florida population. Our results suggest that a large effective population size and negative selection limit the expansion of polymorphic LINE insertions across these populations and that the probability of LINE polymorphisms reaching fixation is extremely low.

Convergence of retrotransposons in oomycetes and plants

Background

Retrotransposons comprise a ubiquitous and abundant class of eukaryotic transposable elements. All members of this class rely on reverse transcriptase activity to produce a DNA copy of the element from the RNA template. However, other activities of the retrotransposon-encoded polyprotein may differ between diverse retrotransposons. The polyprotein domains corresponding to each of these activities may have their own evolutionary history independent from that of the reverse transcriptase, thus underlying the modular view on the evolution of retrotransposons. Furthermore, some transposable elements can independently evolve similar domain architectures by acquiring functionally similar but phylogenetically distinct modules. This convergent evolution of retrotransposons may ultimately suggest similar regulatory pathways underlying the lifecycle of the elements.

Results

Here, we provide new examples of the convergent evolution of retrotransposons of species from two unrelated taxa: green plants and parasitic protozoan oomycetes. In the present study we first analyzed the available genomic sequences of oomycete species and characterized two groups of Ty3/Gypsy long terminal repeat retrotransposons, namely Chronos and Archon, and a subgroup of L1 non-long terminal repeat retrotransposons. The results demonstrated that the retroelements from these three groups each have independently acquired plant-related ribonuclease H domains. This process closely resembles the evolution of retrotransposons in the genomes of green plants. In addition, we showed that Chronos elements captured a chromodomain, mimicking the process of chromodomain acquisition by Chromoviruses, another group of Ty3/Gypsy retrotransposons of plants, fungi, and vertebrates.

Conclusions

Repeated and strikingly similar acquisitions of ribonuclease H domains and chromodomains by different retrotransposon groups from unrelated taxa indicate similar selection pressure acting on these elements. Thus, there are some major trends in the evolution of the structural composition of retrotransposons, and characterizing these trends may enhance the current understanding of the retrotransposon life cycle.

Electronic supplementary material

The online version of this article (doi:10.1186/s13100-017-0087-y) contains supplementary material, which is available to authorized users.

Genome-wide analysis of transposable elements in the coffee berry borer Hypothenemus hampei (Coleoptera: Curculionidae): description of novel families

The coffee berry borer (CBB) Hypothenemus hampei is the most limiting pest of coffee production worldwide. The CBB genome has been recently sequenced; however, information regarding the presence and characteristics of transposable elements (TEs) was not provided. Using systematic searching strategies based on both de novo and homology-based approaches, we present a library of TEs from the draft genome of CBB sequenced by the Colombian Coffee Growers Federation. The library consists of 880 sequences classified as 66% Class I (LTRs: 46%, non-LTRs: 20%) and 34% Class II (DNA transposons: 8%, Helitrons: 16% and MITEs: 10%) elements, including families of the three main LTR (Gypsy, Bel-Pao and Copia) and non-LTR (CR1, Daphne, I/Nimb, Jockey, Kiri, R1, R2 and R4) clades and DNA superfamilies (Tc1-mariner, hAT, Merlin, P, PIF-Harbinger, PiggyBac and Helitron). We propose the existence of novel families: Hypo, belonging to the LTR Gypsy superfamily; Hamp, belonging to non-LTRs; and rosa, belonging to Class II or DNA transposons. Although the rosa clade has been previously described, it was considered to be a basal subfamily of the mariner family. Based on our phylogenetic analysis, including Tc1, mariner, pogo, rosa and Lsra elements from other insects, we propose that rosa and Lsra elements are subfamilies of an independent family of Class II elements termed rosa. The annotations obtained indicate that a low percentage of the assembled CBB genome (approximately 8.2%) consists of TEs. Although these TEs display high diversity, most sequences are degenerate, with few full-length copies of LTR and DNA transposons and several complete and putatively active copies of non-LTR elements. MITEs constitute approximately 50% of the total TEs content, with a high proportion associated with DNA transposons in the Tc1-mariner superfamily.

The Wide Distribution and Change of Target Specificity of R2 Non-LTR Retrotransposons in Animals

Transposons, or transposable elements, are the major components of genomes in most eukaryotes. Some groups of transposons have developed target specificity that limits the integration sites to a specific nonessential sequence or a genomic region to avoid gene disruption caused by insertion into an essential gene. R2 is one of the most intensively investigated groups of sequence-specific non-LTR retrotransposons and is inserted at a specific site inside of 28S ribosomal RNA (rRNA) genes. R2 is known to be distributed among at least six animal phyla even though its occurrence is reported to be patchy. Here, in order to obtain a more detailed picture of the distribution of R2, we surveyed R2 using both in silico screening and degenerate PCR, particularly focusing on actinopterygian fish. We found two families of the R2C lineage from vertebrates, although it has previously only been found in platyhelminthes. We also revealed the apparent movement of insertion sites of a lineage of actinopterygian R2, which was likely concurrent with the acquisition of a 28S rRNA-derived sequence in their 3' UTR. Outside of actinopterygian fish, we revealed the maintenance of a single R2 lineage in birds; the co-existence of four lineages of R2 in the leafcutter bee Megachile rotundata; the first examples of R2 in Ctenophora, Mollusca, and Hemichordata; and two families of R2 showing no target specificity. These findings indicate that R2 is relatively stable and universal, while differences in the distribution and maintenance of R2 lineages probably reflect characteristics of some combination of both R2 lineages and host organisms.

Restricting retrotransposons: a review

Retrotransposons have generated about 40 % of the human genome. This review examines the strategies the cell has evolved to coexist with these genomic "parasites", focussing on the non-long terminal repeat retrotransposons of humans and mice. Some of the restriction factors for retrotransposition, including the APOBECs, MOV10, RNASEL, SAMHD1, TREX1, and ZAP, also limit replication of retroviruses, including HIV, and are part of the intrinsic immune system of the cell. Many of these proteins act in the cytoplasm to degrade retroelement RNA or inhibit its translation. Some factors act in the nucleus and involve DNA repair enzymes or epigenetic processes of DNA methylation and histone modification. RISC and piRNA pathway proteins protect the germline. Retrotransposon control is relaxed in some cell types, such as neurons in the brain, stem cells, and in certain types of disease and cancer, with implications for human health and disease. This review also considers potential pitfalls in interpreting retrotransposon-related data, as well as issues to consider for future research.

[STRUCTURE AND DISTRIBUTION OF THE RETROTRANSPOSON BOV-B LINE]

The classification of mobile elements was discussed. Special attention was devoted to the retroelement of the LINE group: retrotransposon Bov-B LINE. The history of its origin and distribution in the nature was considered. The results of the phenomenon of horizontal transition of the retrotransposon Bov-B LINE between evolutionally distant classes were discussed.

Site-specific non-LTR retrotransposons

Although most of non-long terminal repeat (non-LTR) retrotransposons are incorporated in the host genome almost randomly, some non-LTR retrotransposons are incorporated into specific sequences within a target site. On the basis of structural and phylogenetic features, non-LTR retrotransposons are classified into two large groups, restriction enzyme-like endonuclease (RLE)-encoding elements and apurinic/apyrimidinic endonuclease (APE)-encoding elements. All clades of RLE-encoding non-LTR retrotransposons include site-specific elements. However, only two of more than 20 APE-encoding clades, Tx1 and R1, contain site-specific non-LTR elements. Site-specific non-LTR retrotransposons usually target within multi-copy RNA genes, such as rRNA gene (rDNA) clusters, or repetitive genomic sequences, such as telomeric repeats; this behavior may be a symbiotic strategy to reduce the damage to the host genome. Site- and sequence-specificity are variable even among closely related non-LTR elements and appeared to have changed during evolution. In the APE-encoding elements, the primary determinant of the sequence- specific integration is APE itself, which nicks one strand of the target DNA during the initiation of target primed reverse transcription (TPRT). However, other factors, such as interaction between mRNA and the target DNA, and access to the target region in the nuclei also affect the sequence-specificity. In contrast, in the RLE-encoding elements, DNA-binding motifs appear to affect their sequence-specificity, rather than the RLE domain itself. Highly specific integration properties of these site-specific non-LTR elements make them ideal alternative tools for sequence-specific gene delivery, particularly for therapeutic purposes in human diseases.

Reverse transcriptase genes are highly abundant and transcriptionally active in marine plankton assemblages

Genes encoding reverse transcriptases (RTs) are found in most eukaryotes, often as a component of retrotransposons, as well as in retroviruses and in prokaryotic retroelements. We investigated the abundance, classification and transcriptional status of RTs based on Tara Oceans marine metagenomes and metatranscriptomes encompassing a wide organism size range. Our analyses revealed that RTs predominate large-size fraction metagenomes (>5 μm), where they reached a maximum of 13.5% of the total gene abundance. Metagenomic RTs were widely distributed across the phylogeny of known RTs, but many belonged to previously uncharacterized clades. Metatranscriptomic RTs showed distinct abundance patterns across samples compared with metagenomic RTs. The relative abundances of viral and bacterial RTs among identified RT sequences were higher in metatranscriptomes than in metagenomes and these sequences were detected in all metatranscriptome size fractions. Overall, these observations suggest an active proliferation of various RT-assisted elements, which could be involved in genome evolution or adaptive processes of plankton assemblage.

Evolution of the EKA family of powdery mildew avirulence-effector genes from the ORF 1 of a LINE retrotransposon

Background

The Avrk1 and Avra10 avirulence (AVR) genes encode effectors that increase the pathogenicity of the fungus Blumeria graminis f.sp. hordei (Bgh), the powdery mildew pathogen, in susceptible barley plants. In resistant barley, MLK1 and MLA10 resistance proteins recognize the presence of AVRK1 and AVRA10, eliciting the hypersensitive response typical of gene for gene interactions. Avrk1 and Avra10 have more than 1350 homologues in Bgh genome, forming the EKA (Effectors homologous to Avrk1 and Avra10) gene family.

Results

We tested the hypothesis that the EKA family originated from degenerate copies of Class I LINE retrotransposons by analysing the EKA family in the genome of Bgh isolate DH14 with bioinformatic tools specially developed for the analysis of Transposable Elements (TE) in genomes. The Class I LINE retrotransposon copies homologous to Avrk1 and Avra10 represent 6.5 % of the Bgh annotated genome and, among them, we identified 293 AVR/effector candidate genes. We also experimentally identified peptides that indicated the translation of several predicted proteins from EKA family members, which had higher relative abundance in haustoria than in hyphae.

Conclusions

Our analyses indicate that Avrk1 and Avra10 have evolved from part of the ORF1 gene of Class I LINE retrotransposons. The co-option of Avra10 and Avrk1 as effectors from truncated copies of retrotransposons explains the huge number of homologues in Bgh genome that could act as dynamic reservoirs from which new effector genes may evolve. These data provide further evidence for recruitment of retrotransposons in the evolution of new biological functions.

Electronic supplementary material

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

Ancient Origin of the U2 Small Nuclear RNA Gene-Targeting Non-LTR Retrotransposons Utopia

Most non-long terminal repeat (non-LTR) retrotransposons encoding a restriction-like endonuclease show target-specific integration into repetitive sequences such as ribosomal RNA genes and microsatellites. However, only a few target-specific lineages of non-LTR retrotransposons are distributed widely and no lineage is found across the eukaryotic kingdoms. Here we report the most widely distributed lineage of target sequence-specific non-LTR retrotransposons, designated Utopia. Utopia is found in three supergroups of eukaryotes: Amoebozoa, SAR, and Opisthokonta. Utopia is inserted into a specific site of U2 small nuclear RNA genes with different strength of specificity for each family. Utopia families from oomycetes and wasps show strong target specificity while only a small number of Utopia copies from reptiles are flanked with U2 snRNA genes. Oomycete Utopia families contain an “archaeal” RNase H domain upstream of reverse transcriptase (RT), which likely originated from a plant RNase H gene. Analysis of Utopia from oomycetes indicates that multiple lineages of Utopia have been maintained inside of U2 genes with few copy numbers. Phylogenetic analysis of RT suggests the monophyly of Utopia, and it likely dates back to the early evolution of eukaryotes.

Social evolution. Genomic signatures of evolutionary transitions from solitary to group living

The evolution of eusociality is one of the major transitions in evolution, but the underlying genomic changes are unknown. We compared the genomes of 10 bee species that vary in social complexity, representing multiple independent transitions in social evolution, and report three major findings. First, many important genes show evidence of neutral evolution as a consequence of relaxed selection with increasing social complexity. Second, there is no single road map to eusociality; independent evolutionary transitions in sociality have independent genetic underpinnings. Third, though clearly independent in detail, these transitions do have similar general features, including an increase in constrained protein evolution accompanied by increases in the potential for gene regulation and decreases in diversity and abundance of transposable elements. Eusociality may arise through different mechanisms each time, but would likely always involve an increase in the complexity of gene networks.

Repbase Update, a database of repetitive elements in eukaryotic genomes

Repbase Update (RU) is a database of representative repeat sequences in eukaryotic genomes. Since its first development as a database of human repetitive sequences in 1992, RU has been serving as a well-curated reference database fundamental for almost all eukaryotic genome sequence analyses. Here, we introduce recent updates of RU, focusing on technical issues concerning the submission and updating of Repbase entries and will give short examples of using RU data. RU sincerely invites a broader submission of repeat sequences from the research community.

A New Class of SINEs with snRNA Gene-Derived Heads

Eukaryotic genomes are colonized by various transposons including short interspersed elements (SINEs). The 5' region (head) of the majority of SINEs is derived from one of the three types of RNA genes--7SL RNA, transfer RNA (tRNA), or 5S ribosomal RNA (rRNA)--and the internal promoter inside the head promotes the transcription of the entire SINEs. Here I report a new group of SINEs whose heads originate from either the U1 or U2 small nuclear RNA gene. These SINEs, named SINEU, are distributed among crocodilians and classified into three families. The structures of the SINEU-1 subfamilies indicate the recurrent addition of a U1- or U2-derived sequence onto the 5' end of SINEU-1 elements. SINEU-1 and SINEU-3 are ancient and shared among alligators, crocodiles, and gharials, while SINEU-2 is absent in the alligator genome. SINEU-2 is the only SINE family that was active after the split of crocodiles and gharials. All SINEU families, especially SINEU-3, are preferentially inserted into a family of Mariner DNA transposon, Mariner-N4_AMi. A group of Tx1 non-long terminal repeat retrotransposons designated Tx1-Mar also show target preference for Mariner-N4_AMi, indicating that SINEU was mobilized by Tx1-Mar.

Identification and characterization of expressed retrotransposons in the genome of the Paracoccidioides species complex

BACKGROUND

Species from the Paracoccidioides complex are thermally dimorphic fungi and the causative agents of paracoccidioidomycosis, a deep fungal infection that is the most prevalent systemic mycosis in Latin America and represents the most important cause of death in immunocompetent individuals with systemic mycosis in Brazil. We previously described the identification of eight new families of DNA transposons in Paracoccidioides genomes. In this work, we aimed to identify potentially active retrotransposons in Paracoccidioides genomes.

RESULTS

We identified five different retrotransposon families (four LTR-like and one LINE-like element) in the genomes of three Paracoccidioides isolates. Retrotransposons were present in all of the genomes analyzed. P. brasiliensis and P. lutzii species harbored the same retrotransposon lineages but differed in their copy numbers. In the Pb01, Pb03 and Pb18 genomes, the number of LTR retrotransposons was higher than the number of LINE-like elements, and the LINE-like element RtPc5 was transcribed in Paracoccidioides lutzii (Pb01) but could not be detected in P. brasiliensis (Pb03 and Pb18) by semi-quantitative RT-PCR.

CONCLUSION

Five new potentially active retrotransposons have been identified in the genomic assemblies of the Paracoccidioides species complex using a combined computational and experimental approach. The distribution across the two known species, P. brasiliensis and P. lutzii, and phylogenetics analysis indicate that these elements could have been acquired before speciation occurred. The presence of active retrotransposons in the genome may have implications regarding the evolution and genetic diversification of the Paracoccidioides genus.

The genomes of two key bumblebee species with primitive eusocial organization

Background

The shift from solitary to social behavior is one of the major evolutionary transitions. Primitively eusocial bumblebees are uniquely placed to illuminate the evolution of highly eusocial insect societies. Bumblebees are also invaluable natural and agricultural pollinators, and there is widespread concern over recent population declines in some species. High-quality genomic data will inform key aspects of bumblebee biology, including susceptibility to implicated population viability threats.

Results

We report the high quality draft genome sequences of Bombus terrestris and Bombus impatiens, two ecologically dominant bumblebees and widely utilized study species. Comparing these new genomes to those of the highly eusocial honeybee Apis mellifera and other Hymenoptera, we identify deeply conserved similarities, as well as novelties key to the biology of these organisms. Some honeybee genome features thought to underpin advanced eusociality are also present in bumblebees, indicating an earlier evolution in the bee lineage. Xenobiotic detoxification and immune genes are similarly depauperate in bumblebees and honeybees, and multiple categories of genes linked to social organization, including development and behavior, show high conservation. Key differences identified include a bias in bumblebee chemoreception towards gustation from olfaction, and striking differences in microRNAs, potentially responsible for gene regulation underlying social and other traits.

Conclusions

These two bumblebee genomes provide a foundation for post-genomic research on these key pollinators and insect societies. Overall, gene repertoires suggest that the route to advanced eusociality in bees was mediated by many small changes in many genes and processes, and not by notable expansion or depauperation.

Electronic supplementary material

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

Do larger genomes contain more diverse transposable elements?

Background

The genomes of eukaryotes vary enormously in size, with much of this diversity driven by differences in the abundances of transposable elements (TEs). There is also substantial structural and phylogenetic diversity among TEs, such that they can be classified into distinct classes, superfamilies, and families. Possible relationships between TE diversity (and not just abundance) and genome size have not been investigated to date, though there are reasons to expect either a positive or a negative correlation. This study compares data from 257 species of animals, plants, fungi, and “protists” to determine whether TE diversity at the superfamily level is related to genome size.

Results

No simple relationship was found between TE diversity and genome size. There is no significant correlation across all eukaryotes, but there is a positive correlation for genomes below 500Mbp and a negative correlation among land plants. No relationships were found across animals or within vertebrates. Some TE superfamilies tend to be present across all major groups of eukaryotes, but there is considerable variance in TE diversity in different taxa.

Conclusions

Differences in genome size are thought to arise primarily through accumulation of TEs, but beyond a certain point (~500 Mbp), TE diversity does not increase with genome size. Several possible explanations for these complex patterns are discussed, and recommendations to facilitate future analyses are provided.

Electronic supplementary material

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

A survey of transposable element classification systems--a call for a fundamental update to meet the challenge of their diversity and complexity

The increase of publicly available sequencing data has allowed for rapid progress in our understanding of genome composition. As new information becomes available we should constantly be updating and reanalyzing existing and newly acquired data. In this report we focus on transposable elements (TEs) which make up a significant portion of nearly all sequenced genomes. Our ability to accurately identify and classify these sequences is critical to understanding their impact on host genomes. At the same time, as we demonstrate in this report, problems with existing classification schemes have led to significant misunderstandings of the evolution of both TE sequences and their host genomes. In a pioneering publication Finnegan (1989) proposed classifying all TE sequences into two classes based on transposition mechanisms and structural features: the retrotransposons (class I) and the DNA transposons (class II). We have retraced how ideas regarding TE classification and annotation in both prokaryotic and eukaryotic scientific communities have changed over time. This has led us to observe that: (1) a number of TEs have convergent structural features and/or transposition mechanisms that have led to misleading conclusions regarding their classification, (2) the evolution of TEs is similar to that of viruses by having several unrelated origins, (3) there might be at least 8 classes and 12 orders of TEs including 10 novel orders. In an effort to address these classification issues we propose: (1) the outline of a universal TE classification, (2) a set of methods and classification rules that could be used by all scientific communities involved in the study of TEs, and (3) a 5-year schedule for the establishment of an International Committee for Taxonomy of Transposable Elements (ICTTE).