kangaroo, a mobile element from Volvox carteri, is a member of a newly recognized third class of retrotransposons

Mobilization of retrotransposons as a cause of chromosomal diversification and rapid speciation: the case for the Antarctic teleost genus Trematomus

Background

The importance of transposable elements (TEs) in the genomic remodeling and chromosomal rearrangements that accompany lineage diversification in vertebrates remains the subject of debate. The major impediment to understanding the roles of TEs in genome evolution is the lack of comparative and integrative analyses on complete taxonomic groups. To help overcome this problem, we have focused on the Antarctic teleost genus Trematomus (Notothenioidei: Nototheniidae), as they experienced rapid speciation accompanied by dramatic chromosomal diversity. Here we apply a multi-strategy approach to determine the role of large-scale TE mobilization in chromosomal diversification within Trematomus species.

Results

Despite the extensive chromosomal rearrangements observed in Trematomus species, our measurements revealed strong interspecific genome size conservation. After identifying the DIRS1, Gypsy and Copia retrotransposon superfamilies in genomes of 13 nototheniid species, we evaluated their diversity, abundance (copy numbers) and chromosomal distribution. Four families of DIRS1, nine of Gypsy, and two of Copia were highly conserved in these genomes; DIRS1 being the most represented within Trematomus genomes. Fluorescence in situ hybridization mapping showed preferential accumulation of DIRS1 in centromeric and pericentromeric regions, both in Trematomus and other nototheniid species, but not in outgroups: species of the Sub-Antarctic notothenioid families Bovichtidae and Eleginopsidae, and the non-notothenioid family Percidae.

Conclusions

In contrast to the outgroups, High-Antarctic notothenioid species, including the genus Trematomus, were subjected to strong environmental stresses involving repeated bouts of warming above the freezing point of seawater and cooling to sub-zero temperatures on the Antarctic continental shelf during the past 40 millions of years (My). As a consequence of these repetitive environmental changes, including thermal shocks; a breakdown of epigenetic regulation that normally represses TE activity may have led to sequential waves of TE activation within their genomes. The predominance of DIRS1 in Trematomus species, their transposition mechanism, and their strategic location in “hot spots” of insertion on chromosomes are likely to have facilitated nonhomologous recombination, thereby increasing genomic rearrangements. The resulting centric and tandem fusions and fissions would favor the rapid lineage diversification, characteristic of the nototheniid adaptive radiation.

Electronic supplementary material

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

Retrotransposon Domestication and Control in Dictyostelium discoideum

Transposable elements, identified in all eukaryotes, are mobile genetic units that can change their genomic position. Transposons usually employ an excision and reintegration mechanism, by which they change position, but not copy number. In contrast, retrotransposons amplify via RNA intermediates, increasing their genomic copy number. Hence, they represent a particular threat to the structural and informational integrity of the invaded genome. The social amoeba Dictyostelium discoideum, model organism of the evolutionary Amoebozoa supergroup, features a haploid, gene-dense genome that offers limited space for damage-free transposition. Several of its contemporary retrotransposons display intrinsic integration preferences, for example by inserting next to transfer RNA genes or other retroelements. Likely, any retrotransposons that invaded the genome of the amoeba in a non-directed manner were lost during evolution, as this would result in decreased fitness of the organism. Thus, the positional preference of the Dictyostelium retroelements might represent a domestication of the selfish elements. Likewise, the reduced danger of such domesticated transposable elements led to their accumulation, and they represent about 10% of the current genome of D. discoideum. To prevent the uncontrolled spreading of retrotransposons, the amoeba employs control mechanisms including RNA interference and heterochromatization. Here, we review TRE5-A, DIRS-1 and Skipper-1, as representatives of the three retrotransposon classes in D. discoideum, which make up 5.7% of the Dictyostelium genome. We compile open questions with respect to their mobility and cellular regulation, and suggest strategies, how these questions might be addressed experimentally.

Tyrosine Recombinase Retrotransposons and Transposons

Retrotransposons carrying tyrosine recombinases (YR) are widespread in eukaryotes. The first described tyrosine recombinase mobile element, DIRS1, is a retroelement from the slime mold Dictyostelium discoideum. The YR elements are bordered by terminal repeats related to their replication via free circular dsDNA intermediates. Site-specific recombination is believed to integrate the circle without creating duplications of the target sites. Recently a large number of YR retrotransposons have been described, including elements from fungi (mucorales and basidiomycetes), plants (green algae) and a wide range of animals including nematodes, insects, sea urchins, fish, amphibia and reptiles. YR retrotransposons can be divided into three major groups: the DIRS elements, PAT-like and the Ngaro elements. The three groups form distinct clades on phylogenetic trees based on alignments of reverse transcriptase/ribonuclease H (RT/RH) and YR sequences, and also having some structural distinctions. A group of eukaryote DNA transposons, cryptons, also carry tyrosine recombinases. These DNA transposons do not encode a reverse transcriptase. They have been detected in several pathogenic fungi and oomycetes. Sequence comparisons suggest that the crypton YRs are related to those of the YR retrotransposons. We suggest that the YR retrotransposons arose from the combination of a crypton-like YR DNA transposon and the RT/RH encoding sequence of a retrotransposon. This acquisition must have occurred at a very early point in the evolution of eukaryotes.

The evolution of tyrosine-recombinase elements in Nematoda

Transposable elements can be categorised into DNA and RNA elements based on their mechanism of transposition. Tyrosine recombinase elements (YREs) are relatively rare and poorly understood, despite sharing characteristics with both DNA and RNA elements. Previously, the Nematoda have been reported to have a substantially different diversity of YREs compared to other animal phyla: the Dirs1-like YRE retrotransposon was encountered in most animal phyla but not in Nematoda, and a unique Pat1-like YRE retrotransposon has only been recorded from Nematoda. We explored the diversity of YREs in Nematoda by sampling broadly across the phylum and including 34 genomes representing the three classes within Nematoda. We developed a method to isolate and classify YREs based on both feature organization and phylogenetic relationships in an open and reproducible workflow. We also ensured that our phylogenetic approach to YRE classification identified truncated and degenerate elements, informatively increasing the number of elements sampled. We identified Dirs1-like elements (thought to be absent from Nematoda) in the nematode classes Enoplia and Dorylaimia indicating that nematode model species do not adequately represent the diversity of transposable elements in the phylum. Nematode Pat1-like elements were found to be a derived form of another Pat1-like element that is present more widely in animals. Several sequence features used widely for the classification of YREs were found to be homoplasious, highlighting the need for a phylogenetically-based classification scheme. Nematode model species do not represent the diversity of transposable elements in the phylum.

DIRS1-like retrotransposons are widely distributed among Decapoda and are particularly present in hydrothermal vent organisms

BACKGROUND

Transposable elements are major constituents of eukaryote genomes and have a great impact on genome structure and stability. Considering their mutational abilities, TEs can contribute to the genetic diversity and evolution of organisms. Knowledge of their distribution among several genomes is an essential condition to study their dynamics and to better understand their role in species evolution. DIRS1-like retrotransposons are a particular group of retrotransposons according to their mode of transposition that implies a tyrosine recombinase. To date, they have been described in a restricted number of species in comparison with the LTR retrotransposons. In this paper, we determine the distribution of DIRS1-like elements among 25 decapod species, 10 of them living in hydrothermal vents that correspond to particularly unstable environments.

RESULTS

Using PCR approaches, we have identified 15 new DIRS1-like families in 15 diverse decapod species (shrimps, lobsters, crabs and galatheid crabs). Hydrothermal organisms show a particularly great diversity of DIRS1-like elements with 5 families characterized among Alvinocarididae shrimps and 3 in the galatheid crab Munidopsis recta. Phylogenic analyses show that these elements are divergent toward the DIRS1-like families previously described in other crustaceans and arthropods and form a new clade called AlDIRS1. At larger scale, the distribution of DIRS1-like retrotransposons appears more or less patchy depending on the taxa considered. Indeed, a scattered distribution can be observed in the infraorder Brachyura whereas all the species tested in infraorders Caridea and Astacidea harbor some DIRS1-like elements.

CONCLUSION

Our results lead to nearly double both the number of DIRS1-like elements described to date, and the number of species known to harbor these ones. In this study, we provide the first degenerate primers designed to look specifically for DIRS1-like retrotransposons. They allowed for revealing for the first time a widespread distribution of these elements among a large phylum, here the order Decapoda. They also suggest some peculiar features of these retrotransposons in hydrothermal organisms where a great diversity of elements is already observed. Finally, this paper constitutes the first essential step which allows for considering further studies based on the dynamics of the DIRS1-like retrotransposons among several genomes.

The diversity of retrotransposons and the properties of their reverse transcriptases

A number of abundant mobile genetic elements called retrotransposons reverse transcribe RNA to generate DNA for insertion into eukaryotic genomes. Four major classes of retrotransposons are described here. First, the long-terminal-repeat (LTR) retrotransposons have similar structures and mechanisms to those of the vertebrate retroviruses. Genes that may enable these retrotransposons to leave a cell have been acquired by these elements in a number of animal and plant lineages. Second, the tyrosine recombinase retrotransposons are similar to the LTR retrotransposons except that they have substituted a recombinase for the integrase and recombine into the host chromosomes. Third, the non-LTR retrotransposons use a cleaved chromosomal target site generated by an encoded endonuclease to prime reverse transcription. Finally, the Penelope-like retrotransposons are not well understood but appear to also use cleaved DNA or the ends of chromosomes as primer for reverse transcription. Described in the second part of this review are the enzymatic properties of the reverse transcriptases (RTs) encoded by retrotransposons. The RTs of the LTR retrotransposons are highly divergent in sequence but have similar enzymatic activities to those of retroviruses. The RTs of the non-LTR retrotransposons have several unique properties reflecting their adaptation to a different mechanism of retrotransposition.

Repetitive sequences in complex genomes: structure and evolution

Eukaryotic genomes contain vast amounts of repetitive DNA derived from transposable elements (TEs). Large-scale sequencing of these genomes has produced an unprecedented wealth of information about the origin, diversity, and genomic impact of what was once thought to be "junk DNA." This has also led to the identification of two new classes of DNA transposons, Helitrons and Polintons, as well as several new superfamilies and thousands of new families. TEs are evolutionary precursors of many genes, including RAG1, which plays a role in the vertebrate immune system. They are also the driving force in the evolution of epigenetic regulation and have a long-term impact on genomic stability and evolution. Remnants of TEs appear to be overrepresented in transcription regulatory modules and other regions conserved among distantly related species, which may have implications for our understanding of their impact on speciation.

Maintenance DNA methyltransferase (Met1) and silencing of CpG-methylated foreign DNA in Volvox carteri

DNA methylation plays an important role in the gene-silencing network of higher eukaryotes. We have analyzed the 21.5-kb maintenance methyltransferase (M-MTase) gene, met1, of the multicellular green alga Volvox carteri. The met1 transcript was detected only during the period when DNA replication and cell division are taking place. It encodes a 238 kDa protein containing eight C-terminal activity domains typical of M-MTases, plus upstream DNA-binding domains including the ProDom domain PD003757, which experimental analyses in animal systems have indicated is required for targeting the enzyme to DNA-replication foci. Several insertions of unknown function make Volvox Met1 the largest known member of the Met1/Dnmt1 family. Here we also show that several endogenous transposon families are CpG-methylated in Volvox, which we think causes them to be inactive. This view is supported by the observation that an in vitro CpG-methylated gene introduced into Volvox was maintained in the methylated and silent state over >100 generations. Thus, we believe that Met1 recognizes and perpetuates the in vitro methylation signal, and that the silencing machinery is then able to transduce such a methylation-only signal into a stable heterochromatic (and silent) state.

Quantitative analysis of cell-type specific gene expression in the green alga Volvox carteri

Background

The multicellular alga Volvox carteri possesses only two cell types: mortal, motile somatic cells and potentially immortal, immotile reproductive cells. It is therefore an attractive model system for studying how cell-autonomous cytodifferentiation is programmed within a genome. Moreover, there are ongoing genome projects both in Volvox carteri and in the closely related unicellular alga Chlamydomonas reinhardtii. However, gene sequencing is only the beginning. To identify cell-type specific expression and to determine relative expression rates, we evaluate the potential of real-time RT-PCR for quantifying gene transcript levels.

Results

Here we analyze a diversified pool of 39 target genes by real-time RT-PCR for each cell type. This gene pool contains previously known genes with unknown localization of cellular expression, 28 novel genes which are described in this study for the first time, and a few known, cell-type specific genes as a control. The respective gene products are, for instance, part of photosynthesis, cellular regulation, stress response, or transport processes. We provide expression data for all these genes.

Conclusion

The results show that quantitative real-time RT-PCR is a favorable approach to analyze cell-type specific gene expression in Volvox, which can be extended to a much larger number of genes or to developmental or metabolic mutants. Our expression data also provide a basis for a detailed analysis of individual, previously unknown, cell-type specifically expressed genes.

Distribution and Phylogeny of Penelope-Like Elements in Eukaryotes

Penelope-like elements (PLEs) are a relatively little studied class of eukaryotic retroelements, distinguished by the presence of the GIY-YIG endonuclease domain, the ability of some representatives to retain introns, and the similarity of PLE-encoded reverse transcriptases to telomerases. Although these retrotransposons are abundant in many animal genomes, the reverse transcriptase moiety can also be found in several protists, fungi, and plants, indicating its ancient origin. A comprehensive phylogenetic analysis of PLEs was conducted, based on extended sequence alignments and a considerably expanded data set. PLEs exhibit the pattern of evolution similar to that of non-LTR retrotransposons, which form deep-branching clades dating back to the Precambrian era. However, PLEs seem to have experienced a much higher degree of lineage losses than non-LTR retrotransposons. It is suggested that PLEs and non-LTR retrotransposons are included into a larger eTPRT (eukaryotic target-primed) group of retroelements, characterized by 5' truncation, variable target-site duplication, and the potential of the 3' end to participate in formation of non-autonomous derivatives.

Identification of multiple transcription initiation, polyadenylation, and splice sites in the Drosophila melanogaster TART family of telomeric retrotransposons

The Drosophila non-long terminal repeat (non-LTR) retrotransposons TART and HeT-A specifically retrotranspose to chromosome ends to maintain Drosophila telomeric DNA. Relatively little is known, though, about the regulation of their expression and their retrotransposition to telomeres. We have used rapid amplification of cDNA ends (RACE) to identify multiple transcription initiation and polyadenylation sites for sense and antisense transcripts of three subfamilies of TART elements in Drosophila melanogaster. These results are consistent with the production of an array of TART transcripts. In contrast to other Drosophila non-LTR elements, a major initiation site for sense transcripts was mapped near the 3′ end of the TART 5′-untranslated region (5′-UTR), rather than at the start of the 5′-UTR. A sequence overlapping this sense start site contains a good match to an initiator consensus for the transcription start sites of Drosophila LTR retrotransposons. Interestingly, analysis of 5′ RACE products for antisense transcripts and the GenBank EST database revealed that TART antisense transcripts contain multiple introns. Our results highlight differences between transcription of TART and of other Drosophila non-LTR elements and they provide a foundation for testing the relationship between exceptional aspects of TART transcription and TART's specialized role at telomeres.

Natural history of transposition in the green alga Chlamydomonas reinhardtii: use of the AMT4 locus as an experimental system

The AMT4 locus of the green alga Chlamydomonas reinhardtii, which we mapped to the long arm of chromosome 8, provides a good experimental system for the study of transposition. Most mutations that confer resistance to the toxic ammonium analog methylammonium are in AMT4 and a high proportion of spontaneous mutations are caused by transposon-related events. Among the 15 such events that we have characterized at the molecular level, 9 were associated with insertions of the retrotransposon TOC1, 2 with a small Gulliver-related transposon, and 1 with the Tcr1 transposon. We found that Tcr1 is apparently a foldback transposon with terminal inverted repeats that are much longer and more complex than previously realized. A duplication of Tcr1 yielded a configuration thought to be important for chromosomal evolution. Other mutations in AMT4 were caused by two mobile elements that have not been described before. The sequence of one, which we propose to call the Bill element, indicates that it probably transposes by way of a DNA intermediate and requires functions that it does not encode. The sequence of the other and bioinformatic analysis indicates that it derives from a miniature retrotransposon or TRIM, which we propose to call MRC1 (miniature retrotransposon of Chlamydomonas).

The evolutionary origin of an altruistic gene

Although the conditions favoring altruism are being increasingly understood, the evolutionary origins of the genetic basis for this behavior remain elusive. Here, we show that reproductive altruism (i.e., a sterile soma) in the multicellular green alga, Volvox carteri, evolved via the co-option of a life-history gene whose expression in the unicellular ancestor was conditioned on an environmental cue (as an adaptive strategy to enhance survival at an immediate cost to reproduction) through shifting its expression from a temporal (environmentally induced) into a spatial (developmental) context. The gene belongs to a diverged and structurally heterogeneous multigene family sharing a SAND-like domain (a DNA-binding module involved in gene transcription regulation). To our knowledge, this is the first example of a social gene specifically associated with reproductive altruism, whose origin can be traced back to a solitary ancestor. These findings complement recent proposals that the differentiation of sterile castes in social insects involved the co-option of regulatory networks that control sequential shifts between phases in the life cycle of solitary insects.

Orthologs and paralogs of regA, a master cell-type regulatory gene in Volvox carteri

The multicellular green alga Volvox carteri forma nagariensis has only two cell types: terminally differentiated somatic cells and reproductive cells. The regA gene maintains the terminally differentiated state of the somatic cells, apparently by repressing transcription of genes required for chloroplast biogenesis and thereby preventing cell growth. Because the RegA protein sequence bore no obvious motifs, we are attempting to identify regions of functional importance by searching for strongly conserved domains in RegA orthologs. Here we report the cloning and characterization of regA from the most closely related known taxon, V. carteri f. kawasakiensis. Given the closeness of the relationship between these two formas, their regA genes are surprisingly different: they differ in the number of introns and by several lengthy indels, and they encode proteins that are only 80% identical. We also serendipitously discovered a paralogous gene immediately upstream of each regA locus. The two regA genes, both upstream paralogs and several genes in Chlamydomonas (the closest unicellular relative of Volvox) encode a conserved region (the VARL domain) that contains what appears to be a DNA-binding SAND domain. This discovery has opened up a new avenue for exploring how regA and the terminally differentiated state that it controls evolved.

Transfer RNA gene-targeted integration: an adaptation of retrotransposable elements to survive in the compact Dictyostelium discoideum genome

Almost every organism carries along a multitude of molecular parasites known as transposable elements (TEs). TEs influence their host genomes in many ways by expanding genome size and complexity, rearranging genomic DNA, mutagenizing host genes, and altering transcription levels of nearby genes. The eukaryotic microorganism Dictyostelium discoideum is attractive for the study of fundamental biological phenomena such as intercellular communication, formation of multicellularity, cell differentiation, and morphogenesis. D. discoideum has a highly compacted, haploid genome with less than 1 kb of genomic DNA separating coding regions. Nevertheless, the D. discoideum genome is loaded with 10% of TEs that managed to settle and survive in this inhospitable environment. In depth analysis of D. discoideum genome project data has provided intriguing insights into the evolutionary challenges that mobile elements face when they invade compact genomes. Two different mechanisms are used by D. discoideum TEs to avoid disruption of host genes upon retrotransposition. Several TEs have invented the specific targeting of tRNA gene-flanking regions as a means to avoid integration into coding regions. These elements have been dispersed on all chromosomes, closely following the distribution of tRNA genes. By contrast, TEs that lack bona fide integration specificities show a strong bias to nested integration, thus forming large TE clusters at certain chromosomal loci that are hardly resolved by bioinformatics approaches. We summarize our current view of D. discoideum TEs and present new data from the analysis of the complete sequences of D. discoideum chromosomes 1 and 2, which comprise more than one third of the total genome.

DIRS-1 and the other tyrosine recombinase retrotransposons

DIRS-1 is a retroelement from the slime mold Dictyostelium discoideum. Until recently only two related retrotransposons had been described: PAT from the nematode Panagrellus redivivus and Prt1 from the zygomycete fungus Phycomyces blakesleeanus. Analyses of the reverse transcriptase sequences encoded by these three elements suggested that they were closely related to each other and more distantly related to the Ty3/gypsy Long Terminal Repeat (LTR) retroelements. They have several unusual structural features that distinguish them from typical LTR elements. For instance, they each encode a tyrosine recombinase (YR), but not a DDE-type integrase or an aspartic protease. Although the DIRS-1-related elements are bordered by terminal repeats these differ from typical LTRs in a number of ways. In DIRS-1, for example, the terminal repeats are inverted (complementary), non-identical in sequence, and the outer edges of the terminal sequences are repeated (adjacent to each other) in the internal region. PAT has so-called "split" direct repeats in which the unrelated terminal sequences appear as direct repeats adjacent to each other in the internal region. The only repetition displayed by Prt1 is the presence of short inverted terminal repeats, but the sequenced copy of this element is believed to be a truncated version of an element with a structure resembling DIRS-1. The unusual structure of the terminal repeats of the DIRS1-like elements appears to be related to their replication via free circular intermediates. Site-specific recombination is believed to integrate the circle without creating duplications of the target sites. In recognition of these important distinctions it is proposed that the retrotransposons that encode tyrosine recombinases be called the tyrosine recombinase (or YR) retrotransposons. Recently a large number of additional YR retrotransposons have been described, including elements from fungi (zygomycetes and basidiomycetes), plants (green algae) and a wide range of animals including nematodes, insects, sea urchins, fish and amphibia, while remnants of elements related to DIRS-1 occur in the human genome. The complete set of YR retrotransposons can be divided into two major groups, the DIRS elements and the Ngaro elements, the two groups forming distinct clades on phylogenetic trees based on alignments of RT/RH and recombinase sequences, and also having some structural distinctions. A third group of transposable elements, which we call Cryptons, also carry tyrosine recombinases. These elements do not encode a reverse transcriptase and so are believed to be DNA transposons not retrotransposons. They have been detected in several pathogenic fungi, including the basidiomycete Cryptococcus neoformans, and the ascomycetes Coccidioides posadasii and Histoplasma capsulatum. Sequence comparisons suggest that the Crypton YRs are related to those of the YR retrotransposons. We suggest that the YR retrotransposons arose from the combination of a Crypton-like YR DNA transposon and the RT/RH encoding sequence of a retrotransposon.

DIRS retroelements in arthropods: identification of the recently active TcDirs1 element in the red flour beetle Tribolium castaneum

Members of the DIRS family of retrotransposons differ from most other known retrotransposons in that they encode a tyrosine recombinase (YR), a type of enzyme frequently involved in site-specific recombination. This enzyme is believed to insert the extrachromosomal DNA intermediate of DIRS element retrotransposition into the host genome. DIRS elements have been found in plants, a slime mold, fungi, and a variety of animals including vertebrates, echinoderms and nematodes. They have a somewhat patchy distribution, however, apparently being absent from a number of model organisms such as Saccharomyces cerevisiae, Arabidopsis thaliana and Drosophila melanogaster. In this report we describe the first DIRS retroelement to be identified in an arthropod. This element, TcDirs1, was found in the red flour beetle Tribolium castaneum (Coleoptera). It is generally similar in sequence and structure to several previously described members of the DIRS group: it is bordered by inverted terminal repeats and it has a similar set of protein-coding domains (Gag, reverse transcriptase/ribonuclease H, and the YR), although these are arranged in a novel fashion. TcDirs1 elements exhibit several features indicative of recent activity, such as intact coding regions, a high level of sequence similarity between distinct elements and polymorphic insertion sites. Given their presence in an experimentally tractable host, these potentially active elements might serve as useful models for the study of DIRS element retrotransposition. An element closely related to TcDirs1 was also detected in sequences from a second arthropod, the honey bee Apis mellifera (Hymenoptera), suggesting that these retrotransposons are long-term residents of arthropod genomes.

A New Group of Tyrosine Recombinase-Encoding Retrotransposons

A wide variety of novel tyrosine recombinase (YR)-encoding retrotransposons were identified using data emerging from the various eukaryotic genome sequencing projects. Although many of these elements are clearly members of the previously described DIRS group of YR retrotransposons, a substantial number, including elements from a variety of fungi and animals, belong to a distinct and previously unrecognized group. We refer to these latter elements as the Ngaro group after a representative from zebrafish. Like the members of the DIRS group, Ngaro elements encode proteins bearing reverse transcriptase (RT) and ribonuclease H (RH) domains similar to those of long terminal repeat (LTR) retrotransposons. Phylogenetic analyses based on alignments of RT/RH and YR domains, however, indicate that Ngaro and DIRS are anciently diverged groups. Differences in coding capacity also support the distinction between the two groups. For instance, we found that DIRS elements all encode a protein domain which is similar in sequence to the DNA methyltransferases of certain bacteriophages, whereas this domain is absent from all Ngaro elements. Together, the Ngaro and DIRS groups of YR retrotransposons contain elements with an astonishing diversity in structures, with variations in the nature of the associated repeat sequences and in the arrangement and complement of coding regions. In addition they contain elements with some surprising features, such as spliceosomal introns and long overlapping open reading frames.