Diverse DNA transposons in rotifers of the class Bdelloidea

Genomic Signature of Sexual Reproduction in the Bdelloid Rotifer Macrotrachella quadricornifera

Bdelloid rotifers, common freshwater invertebrates of ancient origin and worldwide distribution have long been thought to be entirely asexual, being the principal exception to the view that in eukaryotes the loss of sex leads to early extinction. That bdelloids are facultatively sexual is shown by a study of allele sharing within a group of closely related bdelloids of the species Macrotrachella quadricornifera, supporting the view that sexual reproduction is essential for long-term success in all eukaryotes.

Diversity and Evolution of pogo and Tc1/mariner Transposons in the Apoidea Genomes

Bees (Apoidea), the largest and most crucial radiation of pollinators, play a vital role in the ecosystem balance. Transposons are widely distributed in nature and are important drivers of species diversity. However, transposons are rarely reported in important pollinators such as bees. Here, we surveyed 37 bee genomesin Apoidea, annotated the pogo and Tc1/mariner transposons in the genome of each species, and performed a phylogenetic analysis and determined their overall distribution. The pogo and Tc1/mariner families showed high diversity and low abundance in the 37 species, and their proportion was significantly higher in solitary bees than in social bees. DD34D/mariner was found to be distributed in almost all species and was found in Apis mellifera, Apis mellifera carnica, Apis mellifera caucasia, and Apis mellifera mellifera, and Euglossa dilemma may still be active. Using horizontal transfer analysis, we found that DD29-30D/Tigger may have experienced horizontal transfer (HT) events. The current study displayed the evolution profiles (including diversity, activity, and abundance) of the pogo and Tc1/mariner transposons across 37 species of Apoidea. Our data revealed their contributions to the genomic variations across these species and facilitated in understanding of the genome evolution of this lineage.

A Survey of Transposon Landscapes in the Putative Ancient Asexual Ostracod Darwinula stevensoni

How asexual reproduction shapes transposable element (TE) content and diversity in eukaryotic genomes remains debated. We performed an initial survey of TE load and diversity in the putative ancient asexual ostracod Darwinula stevensoni. We examined long contiguous stretches of DNA in clones from a genomic fosmid library, totaling about 2.5 Mb, and supplemented these data with results on TE abundance and diversity from an Illumina draft genome. In contrast to other TE studies in putatively ancient asexuals, which revealed relatively low TE content, we found that at least 19% of the fosmid dataset and 26% of the genome assembly corresponded to known transposons. We observed a high diversity of transposon families, including LINE, gypsy, PLE, mariner/Tc, hAT, CMC, Sola2, Ginger, Merlin, Harbinger, MITEs and helitrons, with the prevalence of DNA transposons. The predominantly low levels of sequence diversity indicate that many TEs are or have recently been active. In the fosmid data, no correlation was found between telomeric repeats and non-LTR retrotransposons, which are present near telomeres in other taxa. Most TEs in the fosmid data were located outside of introns and almost none were found in exons. We also report an N-terminal Myb/SANT-like DNA-binding domain in site-specific R4/Dong non-LTR retrotransposons. Although initial results on transposable loads need to be verified with high quality draft genomes, this study provides important first insights into TE dynamics in putative ancient asexual ostracods.

Traveler, a New DD35E Family of Tc1/Mariner Transposons, Invaded Vertebrates Very Recently

The discovery of new members of the Tc1/mariner superfamily of transposons is expected based on the increasing availability of genome sequencing data. Here, we identified a new DD35E family termed Traveler (TR). Phylogenetic analyses of its DDE domain and full-length transposase showed that, although TR formed a monophyletic clade, it exhibited the highest sequence identity and closest phylogenetic relationship with DD34E/Tc1. This family displayed a very restricted taxonomic distribution in the animal kingdom and was only detected in ray-finned fish, anura, and squamata, including 91 vertebrate species. The structural organization of TRs was highly conserved across different classes of animals. Most intact TR transposons had a length of ∼1.5 kb (range 1,072-2,191 bp) and harbored a single open reading frame encoding a transposase of ∼340 aa (range 304-350 aa) flanked by two short-terminal inverted repeats (13-68 bp). Several conserved motifs, including two helix-turn-helix motifs, a GRPR motif, a nuclear localization sequence, and a DDE domain, were also identified in TR transposases. This study also demonstrated the presence of horizontal transfer events of TRs in vertebrates, whereas the average sequence identities and the evolutionary dynamics of TR elements across species and clusters strongly indicated that the TR family invaded the vertebrate lineage very recently and that some of these elements may be currently active, combining the intact TR copies in multiple lineages of vertebrates. These data will contribute to the understanding of the evolutionary history of Tc1/mariner transposons and that of their hosts.

Intruder (DD38E), a recently evolved sibling family of DD34E/Tc1 transposons in animals

BACKGROUND

A family of Tc1/mariner transposons with a characteristic DD38E triad of catalytic amino acid residues, named Intruder (IT), was previously discovered in sturgeon genomes, but their evolutionary landscapes remain largely unknown.

RESULTS

Here, we comprehensively investigated the evolutionary profiles of ITs, and evaluated their cut-and-paste activities in cells. ITs exhibited a narrow taxonomic distribution pattern in the animal kingdom, with invasions into two invertebrate phyla (Arthropoda and Cnidaria) and three vertebrate lineages (Actinopterygii, Agnatha, and Anura): very similar to that of the DD36E/IC family. Some animal orders and species seem to be more hospitable to Tc1/mariner transposons, one order of Amphibia and seven Actinopterygian orders are the most common orders with horizontal transfer events and have been invaded by all four families (DD38E/IT, DD35E/TR, DD36E/IC and DD37E/TRT) of Tc1/mariner transposons, and eight Actinopterygii species were identified as the major hosts of these families. Intact ITs have a total length of 1.5-1.7 kb containing a transposase gene flanked by terminal inverted repeats (TIRs). The phylogenetic tree and sequence identity showed that IT transposases were most closely related to DD34E/Tc1. ITs have been involved in multiple events of horizontal transfer in vertebrates and have invaded most lineages recently (< 5 million years ago) based on insertion age analysis. Accordingly, ITs presented high average sequence identity (86-95%) across most vertebrate species, suggesting that some are putatively active. ITs can transpose in human HeLa cells, and the transposition efficiency of consensus TIRs was higher than that of the TIRs of natural isolates.

CONCLUSIONS

We conclude that DD38E/IT originated from DD34E/Tc1 and can be detected in two invertebrate phyla (Arthropoda and Cnidaria), and in three vertebrate lineages (Actinopterygii, Agnatha and Anura). IT has experienced multiple HT events in animals, dominated by recent amplifications in most species and has high identity among vertebrate taxa. Our reconstructed IT transposon vector designed according to the sequence from the "cat" genome showed high cut-and-paste activity. The data suggest that IT has been acquired recently and is active in many species. This study is meaningful for understanding the evolution of the Tc1/mariner superfamily members and their hosts.

The genome, transcriptome, and proteome of the fish parasite Pomphorhynchus laevis (Acanthocephala)

Thorny-headed worms (Acanthocephala) are endoparasites exploiting Mandibulata (Arthropoda) and Gnathostomata (Vertebrata). Despite their world-wide occurrence and economic relevance as a pest, genome and transcriptome assemblies have not been published before. However, such data might hold clues for a sustainable control of acanthocephalans in animal production. For this reason, we present the first draft of an acanthocephalan nuclear genome, besides the mitochondrial one, using the fish parasite Pomphorhynchus laevis (Palaeacanthocephala) as a model. Additionally, we have assembled and annotated the transcriptome of this species and the proteins encoded. A hybrid assembly of long and short reads resulted in a near-complete P. laevis draft genome of ca. 260 Mb, comprising a large repetitive portion of ca. 63%. Numbers of transcripts and translated proteins (35,683) were within the range of other members of the Rotifera-Acanthocephala clade. Our data additionally demonstrate a significant reorganization of the acanthocephalan gene repertoire. Thus, more than 20% of the usually conserved metazoan genes were lacking in P. laevis. Ontology analysis of the retained genes revealed many connections to the incorporation of carotinoids. These are probably taken up via the surface together with lipids, thus accounting for the orange coloration of P. laevis. Furthermore, we found transcripts and protein sequences to be more derived in P. laevis than in rotifers from Monogononta and Bdelloidea. This was especially the case in genes involved in energy metabolism, which might reflect the acanthocephalan ability to use the scarce oxygen in the host intestine for respiration and simultaneously carry out fermentation. Increased plasticity of the gene repertoire through the integration of foreign DNA into the nuclear genome seems to be another underpinning factor of the evolutionary success of acanthocephalans. In any case, energy-related genes and their proteins may be considered as candidate targets for the acanthocephalan control.

Genome Evolution in Outcrossing vs. Selfing vs. Asexual Species

A major current molecular evolution challenge is to link comparative genomic patterns to species' biology and ecology. Breeding systems are pivotal because they affect many population genetic processes and thus genome evolution. We review theoretical predictions and empirical evidence about molecular evolutionary processes under three distinct breeding systems-outcrossing, selfing, and asexuality. Breeding systems may have a profound impact on genome evolution, including molecular evolutionary rates, base composition, genomic conflict, and possibly genome size. We present and discuss the similarities and differences between the effects of selfing and clonality. In reverse, comparative and population genomic data and approaches help revisiting old questions on the long-term evolution of breeding systems.

Comparative genomics of bdelloid rotifers: Insights from desiccating and nondesiccating species

Bdelloid rotifers are a class of microscopic invertebrates that have existed for millions of years apparently without sex or meiosis. They inhabit a variety of temporary and permanent freshwater habitats globally, and many species are remarkably tolerant of desiccation. Bdelloids offer an opportunity to better understand the evolution of sex and recombination, but previous work has emphasised desiccation as the cause of several unusual genomic features in this group. Here, we present high-quality whole-genome sequences of 3 bdelloid species: Rotaria macrura and R. magnacalcarata, which are both desiccation intolerant, and Adineta ricciae, which is desiccation tolerant. In combination with the published assembly of A. vaga, which is also desiccation tolerant, we apply a comparative genomics approach to evaluate the potential effects of desiccation tolerance and asexuality on genome evolution in bdelloids. We find that ancestral tetraploidy is conserved among all 4 bdelloid species, but homologous divergence in obligately aquatic Rotaria genomes is unexpectedly low. This finding is contrary to current models regarding the role of desiccation in shaping bdelloid genomes. In addition, we find that homologous regions in A. ricciae are largely collinear and do not form palindromic repeats as observed in the published A. vaga assembly. Consequently, several features interpreted as genomic evidence for long-term ameiotic evolution are not general to all bdelloid species, even within the same genus. Finally, we substantiate previous findings of high levels of horizontally transferred nonmetazoan genes in both desiccating and nondesiccating bdelloid species and show that this unusual feature is not shared by other animal phyla, even those with desiccation-tolerant representatives. These comparisons call into question the proposed role of desiccation in mediating horizontal genetic transfer.

Molecular Evolution of piggyBac Superfamily: From Selfishness to Domestication

The piggyBac transposable element was originally isolated from the cabbage looper moth, Trichoplusia ni, in the 1980s. Despite its early discovery and specificity compared to the other Class II elements, the diversity and evolution of this superfamily have been only partially analyzed. Two main types of elements can be distinguished: the piggyBac-like elements (PBLE) with terminal inverted repeats, untranslated region, and an open reading frame encoding a transposase, and the piggyBac-derived sequences (PGBD), containing a sequence derived from a piggyBac transposase, and which correspond to domesticated elements. To define the distribution, their structural diversity and phylogenetic relationships, analyses were conducted using known PBLE and PGBD sequences to scan databases. From this data mining, numerous new sequences were characterized (50 for PBLE and 396 for PGBD). Structural analyses suggest that four groups of PBLE can be defined according to the presence/absence of sub-terminal repeats. The transposase is characterized by highly variable catalytic domain and C-terminal region. There is no relationship between the structural groups and the phylogeny of these PBLE elements. The PGBD are clearly structured into nine main groups. A new group of domesticated elements is suspected in Neopterygii and the remaining eight previously described elements have been investigated in more detail. In all cases, these sequences are no longer transposable elements, the catalytic domain of the ancestral transposase is not always conserved, but they are under strong purifying selection. The phylogeny of both PBLE and PGBD suggests multiple and independent domestication events of PGBD from different PBLE ancestors.

Study of Transposable Elements and Their Genomic Impact

Transposable elements (TEs) have been considered traditionally as junk DNA, i.e., DNA sequences that despite representing a high proportion of genomes had no evident cellular functions. However, over the last decades, it has become undeniable that not only TE-derived DNA sequences have (and had) a fundamental role during genome evolution, but also TEs have important implications in the origin and evolution of many genomic disorders. This concise review provides a brief overview of the different types of TEs that can be found in genomes, as well as a list of techniques and methods used to study their impact and mobilization. Some of these techniques will be covered in detail in this Method Book.

TRT, a Vertebrate and Protozoan Tc1-Like Transposon: Current Activity and Horizontal Transfer

We report a Danio rerio transposon named DrTRT, for D. rerio Transposon Related to Tc1. The complete sequence of the DrTRT transposon is 1,563 base pairs (bp) in length, and its transposase putatively encodes a 338-amino acid protein that harbors a DD37E motif in its catalytic domain. We present evidence based on searches of publicly available genomes that TRT elements commonly occur in vertebrates and protozoa. Phylogenetic and functional domain comparisons confirm that TRT constitutes a new subfamily within the Tc1 family. Hallmark features of having no premature termination codons within the transposase, the presence of all expected functional domains, and its occurrence in the bony fish transcriptome suggest that TRT might have current or recent activity in these species. Further analysis showed that the activity of TRT elements in these species might have arisen about between 4 and 19 Ma. Interestingly, our results also implied that the widespread distribution of TRT among fishes, frog, and snakes is the result of multiple independent HT events, probably from bony fishes to snakes or frog. Finally, the mechanisms underlying horizontal transfer of TRT elements are discussed.

Multitasking of the piRNA Silencing Machinery: Targeting Transposable Elements and Foreign Genes in the Bdelloid Rotifer Adineta vaga

RNA-mediated silencing processes play a key role in silencing of transposable elements, especially in the germ line, where piwi-interacting RNAs (piRNAs) are responsible for suppressing transposon mobility and maintaining genome integrity. We previously reported that the genome of Adineta vaga, the first sequenced representative of the phylum Rotifera (class Bdelloidea), is characterized by massive levels of horizontal gene transfer, by unusually low transposon content, and by highly diversified RNA-mediated silencing machinery. Here, we investigate genome-wide distribution of pi-like small RNAs, which in A. vaga are 25-31 nucleotides in length and have a strong 5'-uridine bias, while lacking ping-pong amplification signatures. In agreement with expectations, 71% of mapped reads corresponded to annotated transposons, with 93% of these reads being in the antisense orientation. Unexpectedly, a significant fraction of piRNAs originate from predicted coding regions corresponding to genes of putatively foreign origin. The distribution of piRNAs across foreign genes is not biased toward 3'-UTRs, instead resembling transposons in uniform distribution pattern throughout the gene body, and in predominantly antisense orientation. We also find that genes with small RNA coverage, including a number of genes of metazoan origin, are characterized by higher occurrence of telomeric repeats in the surrounding genomic regions, and by higher density of transposons in the vicinity, which have the potential to promote antisense transcription. Our findings highlight the complex interplay between RNA-based silencing processes and acquisition of genes at the genome periphery, which can result either in their loss or eventual domestication and integration into the host genome.

piggyBac Transposon

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

A Novel Type Pathway-Specific Regulator and Dynamic Genome Environments of a Solanapyrone Biosynthesis Gene Cluster in the Fungus Ascochyta rabiei

Secondary metabolite genes are often clustered together and situated in particular genomic regions, like the subtelomere, that can facilitate niche adaptation in fungi. Solanapyrones are toxic secondary metabolites produced by fungi occupying different ecological niches. Full-genome sequencing of the ascomycete Ascochyta rabiei revealed a solanapyrone biosynthesis gene cluster embedded in an AT-rich region proximal to a telomere end and surrounded by Tc1/Mariner-type transposable elements. The highly AT-rich environment of the solanapyrone cluster is likely the product of repeat-induced point mutations. Several secondary metabolism-related genes were found in the flanking regions of the solanapyrone cluster. Although the solanapyrone cluster appears to be resistant to repeat-induced point mutations, a P450 monooxygenase gene adjacent to the cluster has been degraded by such mutations. Among the six solanapyrone cluster genes (sol1 to sol6), sol4 encodes a novel type of Zn(II)2Cys6 zinc cluster transcription factor. Deletion of sol4 resulted in the complete loss of solanapyrone production but did not compromise growth, sporulation, or virulence. Gene expression studies with the sol4 deletion and sol4-overexpressing mutants delimited the boundaries of the solanapyrone gene cluster and revealed that sol4 is likely a specific regulator of solanapyrone biosynthesis and appears to be necessary and sufficient for induction of the solanapyrone cluster genes. Despite the dynamic surrounding genomic regions, the solanapyrone gene cluster has maintained its integrity, suggesting important roles of solanapyrones in fungal biology.

Do Species Exist in Asexuals? Theory and Evidence from Bdelloid Rotifers

The possibility for independently evolving entities to form and persist in the absence of sexual recombination in eukaryotes has been questioned; nevertheless, there are organisms that are known to be asexual and that have apparently diversified into multiple species as recognized by taxonomists. These organisms have therefore been identified as an evolutionary paradox. We explore three alternative hypotheses attempting to solve the apparent paradox, focusing on bdelloid rotifers, the most studied group of organisms in which all species are considered asexual: (1) they may have some hidden form of sex; (2) species do not represent biological entities but simply convenient names; and (3) sex may not be a necessary requirement for speciation. We provide ample evidence against the first two hypotheses, reporting several studies supporting (1) bdelloids asexuality from different approaches, and (2) the existence of species from genetics, jaw morphology, ecology, and physiology. Thus, we (3) explore the role of sex in speciation comparing bdelloid and monogonont rotifers, and conclude with some caveats that could still change our understanding of bdelloid species.

Genome evolution in outcrossing versus selfing versus asexual species

A major current molecular evolution challenge is to link comparative genomic patterns to species' biology and ecology. Breeding systems are pivotal because they affect many population genetic processes, and thus genome evolution. We review theoretical predictions and empirical evidence about molecular evolutionary processes under three distinct breeding systems-outcrossing, selfing, and asexuality. Breeding systems may have a profound impact on genome evolution, including molecular evolutionary rates, base composition, genomic conflict, and possibly genome size. However, while asexual species essentially conform to theoretical predictions, the situation is less simple in selfing species. We discuss the possible reasons to potentially explain this paradox. In reverse, comparative and population genomic data and approaches help revisiting old questions on the long-term evolution of breeding systems.

Genome walking in eukaryotes

Genome walking is a molecular procedure for the direct identification of nucleotide sequences from purified genomes. The only requirement is the availability of a known nucleotide sequence from which to start. Several genome walking methods have been developed in the last 20 years, with continuous improvements added to the first basic strategies, including the recent coupling with next generation sequencing technologies. This review focuses on the use of genome walking strategies in several aspects of the study of eukaryotic genomes. In a first part, the analysis of the numerous strategies available is reported. The technical aspects involved in genome walking are particularly intriguing, also because they represent the synthesis of the talent, the fantasy and the intelligence of several scientists. Applications in which genome walking can be employed are systematically examined in the second part of the review, showing the large potentiality of this technique, including not only the simple identification of nucleotide sequences but also the analysis of large collections of mutants obtained from the insertion of DNA of viral origin, transposons and transfer DNA (T-DNA) constructs. The enormous amount of data obtained indicates that genome walking, with its large range of applicability, multiplicity of strategies and recent developments, will continue to have much to offer for the rapid identification of unknown sequences in several fields of genomic research.

The struggle for life of the genome's selfish architects

Transposable elements (TEs) were first discovered more than 50 years ago, but were totally ignored for a long time. Over the last few decades they have gradually attracted increasing interest from research scientists. Initially they were viewed as totally marginal and anecdotic, but TEs have been revealed as potentially harmful parasitic entities, ubiquitous in genomes, and finally as unavoidable actors in the diversity, structure, and evolution of the genome. Since Darwin's theory of evolution, and the progress of molecular biology, transposable elements may be the discovery that has most influenced our vision of (genome) evolution. In this review, we provide a synopsis of what is known about the complex interactions that exist between transposable elements and the host genome. Numerous examples of these interactions are provided, first from the standpoint of the genome, and then from that of the transposable elements. We also explore the evolutionary aspects of TEs in the light of post-Darwinian theories of evolution.

Genome structure of bdelloid rotifers: shaped by asexuality or desiccation?

Bdelloid rotifers are microscopic invertebrate animals best known for their ancient asexuality and the ability to survive desiccation at any life stage. Both factors are expected to have a profound influence on their genome structure. Recent molecular studies demonstrated that, although the gene-rich regions of bdelloid genomes are organized as colinear pairs of closely related sequences and depleted in repetitive DNA, subtelomeric regions harbor diverse transposable elements and horizontally acquired genes of foreign origin. Although asexuality is expected to result in depletion of deleterious transposons, only desiccation appears to have the power to produce all the uncovered genomic peculiarities. Repair of desiccation-induced DNA damage would require the presence of a homologous template, maintaining colinear pairs in gene-rich regions and selecting against insertion of repetitive DNA that might cause chromosomal rearrangements. Desiccation may also induce a transient state of competence in recovering animals, allowing them to acquire environmental DNA. Even if bdelloids engage in rare or obscure forms of sexual reproduction, all these features could still be present. The relative contribution of asexuality and desiccation to genome organization may be clarified by analyzing whole-genome sequences and comparing foreign gene and transposon content in species which lost the ability to survive desiccation.

A subtelomeric non-LTR retrotransposon Hebe in the bdelloid rotifer Adineta vaga is subject to inactivation by deletions but not 5' truncations

BACKGROUND

Rotifers of the class Bdelloidea are microscopic freshwater invertebrates best known for: their capacity for anhydrobiosis; the lack of males and meiosis; and for the ability to capture genes from other non-metazoan species. Although genetic exchange between these animals might take place by non-canonical means, the overall lack of meiosis and syngamy should greatly impair the ability of transposable elements (TEs) to spread in bdelloid populations. Previous studies demonstrated that bdelloid chromosome ends, in contrast to gene-rich regions, harbour various kinds of TEs, including specialized telomere-associated retroelements, as well as DNA TEs and retrovirus-like retrotransposons which are prone to horizontal transmission. Vertically-transmitted retrotransposons have not previously been reported in bdelloids and their identification and studies of the patterns of their distribution and evolution could help in the understanding of the high degree of TE compartmentalization within bdelloid genomes.

RESULTS

We identified and characterized a non-long terminal repeat (LTR) retrotransposon residing primarily in subtelomeric regions of the genome in the bdelloid rotifer Adineta vaga. Contrary to the currently prevailing views on the mode of proliferation of non-LTR retrotransposons, which results in frequent formation of 5'-truncated ('dead-on-arrival') copies due to the premature disengagement of the element-encoded reverse transcriptase from its template, this non-LTR element, Hebe, is represented only by non-5'-truncated copies. Most of these copies, however, were subject to internal deletions associated with microhomologies, a hallmark of non-homologous end-joining events.

CONCLUSIONS

The non-LTR retrotransposon Hebe from the bdelloid rotifer A. vaga was found to undergo frequent microhomology-associated deletions, rather than 5'-terminal truncations characteristic of this class of retrotransposons, and to exhibit preference for telomeric localization. These findings represent the first example of a vertically transmitted putatively deleterious TE in bdelloids, and may indicate the involvement of microhomology-mediated non-homologous end-joining in desiccation-induced double-strand break repair at the genome periphery.

DNA transposons: nature and applications in genomics

Repeated DNA makes up a large fraction of a typical mammalian genome, and some repetitive elements are able to move within the genome (transposons and retrotransposons). DNA transposons move from one genomic location to another by a cut-and-paste mechanism. They are powerful forces of genetic change and have played a significant role in the evolution of many genomes. As genetic tools, DNA transposons can be used to introduce a piece of foreign DNA into a genome. Indeed, they have been used for transgenesis and insertional mutagenesis in different organisms, since these elements are not generally dependent on host factors to mediate their mobility. Thus, DNA transposons are useful tools to analyze the regulatory genome, study embryonic development, identify genes and pathways implicated in disease or pathogenesis of pathogens, and even contribute to gene therapy. In this review, we will describe the nature of these elements and discuss recent advances in this field of research, as well as our evolving knowledge of the DNA transposons most widely used in these studies.

Origin of the cell nucleus, mitosis and sex: roles of intracellular coevolution

BACKGROUND

The transition from prokaryotes to eukaryotes was the most radical change in cell organisation since life began, with the largest ever burst of gene duplication and novelty. According to the coevolutionary theory of eukaryote origins, the fundamental innovations were the concerted origins of the endomembrane system and cytoskeleton, subsequently recruited to form the cell nucleus and coevolving mitotic apparatus, with numerous genetic eukaryotic novelties inevitable consequences of this compartmentation and novel DNA segregation mechanism. Physical and mutational mechanisms of origin of the nucleus are seldom considered beyond the long-standing assumption that it involved wrapping pre-existing endomembranes around chromatin. Discussions on the origin of sex typically overlook its association with protozoan entry into dormant walled cysts and the likely simultaneous coevolutionary, not sequential, origin of mitosis and meiosis.

RESULTS

I elucidate nuclear and mitotic coevolution, explaining the origins of dicer and small centromeric RNAs for positionally controlling centromeric heterochromatin, and how 27 major features of the cell nucleus evolved in four logical stages, making both mechanisms and selective advantages explicit: two initial stages (origin of 30 nm chromatin fibres, enabling DNA compaction; and firmer attachment of endomembranes to heterochromatin) protected DNA and nascent RNA from shearing by novel molecular motors mediating vesicle transport, division, and cytoplasmic motility. Then octagonal nuclear pore complexes (NPCs) arguably evolved from COPII coated vesicle proteins trapped in clumps by Ran GTPase-mediated cisternal fusion that generated the fenestrated nuclear envelope, preventing lethal complete cisternal fusion, and allowing passive protein and RNA exchange. Finally, plugging NPC lumens by an FG-nucleoporin meshwork and adopting karyopherins for nucleocytoplasmic exchange conferred compartmentation advantages. These successive changes took place in naked growing cells, probably as indirect consequences of the origin of phagotrophy. The first eukaryote had 1-2 cilia and also walled resting cysts; I outline how encystation may have promoted the origin of meiotic sex. I also explain why many alternative ideas are inadequate.

CONCLUSION

Nuclear pore complexes are evolutionary chimaeras of endomembrane- and mitosis-related chromatin-associated proteins. The keys to understanding eukaryogenesis are a proper phylogenetic context and understanding organelle coevolution: how innovations in one cell component caused repercussions on others.

Transposable elements and factors influencing their success in eukaryotes

Recent advances in genome sequencing have led to a vast accumulation of transposable element data. Consideration of the genome sequencing projects in a phylogenetic context reveals that despite the hundreds of eukaryotic genomes that have been sequenced, a strong bias in sampling exists. There is a general under-representation of unicellular eukaryotes and a dearth of genome projects in many branches of the eukaryotic phylogeny. Among sequenced genomes, great variation in genome size exists, however, little difference in the total number of cellular genes is observed. For many eukaryotes, the remaining genomic space is extremely dynamic and predominantly composed of a menagerie of populations of transposable elements. Given the dynamic nature of the genomic niche filled by transposable elements, it is evident that these elements have played an important role in genome evolution. The contribution of transposable elements to genome architecture and to the advent of genetic novelty is likely to be dependent, at least in part, on the transposition mechanism, diversity, number, and rate of turnover of transposable elements in the genome at any given time. The focus of this review is the discussion of some of the forces that act to shape transposable element diversity within and between genomes.

Structure-based discovery and description of plant and animal Helitrons

Helitrons are recently discovered eukaryotic transposons that are predicted to amplify by a rolling-circle mechanism. They are present in most plant and animal species investigated, but were previously overlooked partly because they lack terminal repeats and do not create target site duplications. Helitrons are particularly abundant in flowering plants, where they frequently acquire, and sometimes express, 1 or more gene fragments. A structure-based search protocol was developed to find Helitrons and was used to analyze several plant and animal genomes, leading to the discovery of hundreds of new Helitrons. Analysis of these Helitrons has uncovered mechanisms of element evolution, including end creation and sequence acquisition. Preferential accumulation in gene-poor regions and target site specificities were also identified. Overall, these studies provide insights into the transposition and evolution of Helitrons and their contributions to evolved gene content and genome structure.

A single-copy IS5-like transposon in the genome of a bdelloid rotifer

In the course of sequencing telomeric chromosomal regions of the bdelloid rotifer Adineta vaga, we encountered an unusual DNA transposon. Unlike other bdelloid and, more generally, eukaryotic transposable elements (TEs), it exhibits similarity to prokaryotic insertion sequences (ISs). Phylogenetic analysis indicates that this transposon, named IS5_Av, is related to the ISL2 group of the IS5 family of bacterial IS elements. Despite the apparent intactness of the single open reading frame coding for a DDE transposase and the perfect identity of its 213-bp terminal inverted repeats (TIRs), the element is present in only one copy per diploid genome. It does not exhibit any detectable levels of transcription, so that its transposase gene appears to be silent in the bdelloid host. Although horizontal transfers of TEs between kingdoms are not known to happen in nature, it appears likely that IS5_Av underwent integration into the A. vaga genome relatively recently, but was not successful in adapting to the new host and failed to increase in copy number. Alternatively, it might be the only known member of a novel eukaryotic DNA TE superfamily which is so rare that its other members, if any, have not yet been identified in eukaryotic genomes sequenced to date.

Automated paleontology of repetitive DNA with REANNOTATE

BACKGROUND

Dispersed repeats are a major component of eukaryotic genomes and drivers of genome evolution. Annotation of DNA sequences homologous to known repetitive elements has been mainly performed with the program REPEATMASKER. Sequences annotated by REPEATMASKER often correspond to fragments of repetitive elements resulting from the insertion of younger elements or other rearrangements. Although REPEATMASKER annotation is indispensable for studying genome biology, this annotation does not contain much information on the common origin of fossil fragments that share an insertion event, especially where clusters of nested insertions of repetitive elements have occurred.

RESULTS

Here I present REANNOTATE, a computational tool to process REPEATMASKER annotation for automated i) defragmentation of dispersed repetitive elements, ii) resolution of the temporal order of insertions in clusters of nested elements, and iii) estimating the age of the elements, if they have long terminal repeats. I have re-annotated the repetitive content of human chromosomes, providing evidence for a recent expansion of satellite repeats on the Y chromosome and, from the retroviral age distribution, for a higher rate of evolution on the Y relative to autosomes.

CONCLUSION

REANNOTATE is ready to process existing annotation for automated evolutionary analysis of all types of complex repeats in any genome. The tool is freely available under the GPL at http://www.bioinformatics.org/reannotate.

Degenerate tetraploidy was established before bdelloid rotifer families diverged

Rotifers of Class Bdelloidea are abundant freshwater invertebrates known for their remarkable ability to survive desiccation and their lack of males and meiosis. Sequencing and annotation of approximately 50-kb regions containing the four hsp82 heat shock genes of the bdelloid Philodina roseola, each located on a separate chromosome, have suggested that its genome is that of a degenerate tetraploid. In order to determine whether a similar structure exists in a bdelloid distantly related to P. roseola and if degenerate tetraploidy was established before the two species separated, we sequenced regions containing the hsp82 genes of a bdelloid belonging to a different family, Adineta vaga, and the histone gene clusters of P. roseola and A. vaga. Our findings are entirely consistent with degenerate tetraploidy and show that it was established before the two bdelloid families diverged and therefore probably before the bdelloid radiation.

Transposition of Mboumar-9: identification of a new naturally active mariner-family transposon

Although mariner transposons are widespread in animal genomes, the vast majority harbor multiple inactivating mutations and only two naturally occurring elements are known to be active. Previously, we discovered a mariner-family transposon, Mboumar, in the satellite DNA of the ant Messor bouvieri. Several copies of the transposon contain a full-length open reading frame, including Mboumar-9, which has 64% nucleotide identity to Mos1 of Drosophila mauritiana. To determine whether Mboumar is currently active, we expressed and purified the Mboumar-9 transposase and demonstrate that it is able to catalyze the movement of a transposon from one plasmid to another in a genetic in vitro hop assay. The efficiency is comparable to that of the well-characterized mariner transposon Mos1. Transposon insertions were precise and were flanked by TA duplications, a hallmark of mariner transposition. Mboumar has been proposed to have a role in the evolution and maintenance of satellite DNA in M. bouvieri and its activity provides a means to examine the involvement of the transposon in the genome dynamics of this organism.

piggyBac can bypass DNA synthesis during cut and paste transposition

DNA synthesis is considered a defining feature in the movement of transposable elements. In determining the mechanism of piggyBac transposition, an insect transposon that is being increasingly used for genome manipulation in a variety of systems including mammalian cells, we have found that DNA synthesis can be avoided during piggyBac transposition, both at the donor site following transposon excision and at the insertion site following transposon integration. We demonstrate that piggyBac transposon excision occurs through the formation of transient hairpins on the transposon ends and that piggyBac target joining occurs by the direct attack of the 3'OH transposon ends on to the target DNA. This is the same strategy for target joining used by the members of DDE superfamily of transposases and retroviral integrases. Analysis of mutant piggyBac transposases in vitro and in vivo using a piggyBac transposition system we have established in Saccharomyces cerevisiae suggests that piggyBac transposase is a member of the DDE superfamily of recombinases, an unanticipated result because of the lack of sequence similarity between piggyBac and DDE family of recombinases.

Extreme resistance of bdelloid rotifers to ionizing radiation

Rotifers of class Bdelloidea are common invertebrate animals with highly unusual characteristics, including apparently obligate asexuality, the ability to resume reproduction after desiccation at any life stage, and a paucity of transposable genetic elements of types not prone to horizontal transmission. We find that bdelloids are also extraordinarily resistant to ionizing radiation (IR). Reproduction of the bdelloids Adineta vaga and Philodina roseola is much more resistant to IR than that of Euchlanis dilatata, a rotifer belonging to the desiccation-intolerant and facultatively sexual class Monogononta, and all other animals for which we have found relevant data. By analogy with the desiccation- and radiation-resistant bacterium Deinococcus radiodurans, we suggest that the extraordinary radiation resistance of bdelloid rotifers is a consequence of their evolutionary adaptation to survive episodes of desiccation encountered in their characteristic habitats and that the damage incurred in such episodes includes DNA breakage that is repaired upon rehydration. Such breakage and repair may have maintained bdelloid chromosomes as colinear pairs and kept the load of transposable genetic elements low and may also have contributed to the success of bdelloid rotifers in avoiding the early extinction suffered by most asexuals.

Evidence for degenerate tetraploidy in bdelloid rotifers

Rotifers of class Bdelloidea have evolved for millions of years apparently without sexual reproduction. We have sequenced 45- to 70-kb regions surrounding the four copies of the hsp82 gene of the bdelloid rotifer Philodina roseola, each of which is on a separate chromosome. The four regions comprise two colinear gene-rich pairs with gene content, order, and orientation conserved within each pair. Only a minority of genes are common to both pairs, also in the same orientation and order, but separated by gene-rich segments present in only one or the other pair. The pattern is consistent with degenerate tetraploidy with numerous segmental deletions, some in one pair of colinear chromosomes and some in the other. Divergence in 1,000-bp windows varies along an alignment of a colinear pair, from zero to as much as 20% in a pattern consistent with gene conversion associated with recombinational repair of DNA double-strand breaks. Although pairs of colinear chromosomes are a characteristic of sexually reproducing diploids and polyploids, a quite different explanation for their presence in bdelloids is suggested by the recent finding that bdelloid rotifers can recover and resume reproduction after suffering hundreds of radiation-induced DNA double-strand breaks per oocyte nucleus. Because bdelloid primary oocytes are in G(1) and therefore lack sister chromatids, we propose that bdelloid colinear chromosome pairs are maintained as templates for the repair of DNA double-strand breaks caused by the frequent desiccation and rehydration characteristic of bdelloid habitats.

Transposable elements in Coffea (Gentianales: Rubiacea) transcripts and their role in the origin of protein diversity in flowering plants

Transposable elements are major components of plant genomes and they influence their evolution, acting as recombination hot spots, acquiring specific cell functions or becoming part of protein-coding regions. The latter is the subject of the present analysis. This study is a report on the annotation of transposable elements (TEs) in expressed sequences of Coffea arabica, Coffea canephora and Coffea racemosa, showing the occurrence of 383 ESTs and 142 unigenes with TE fragments in these three Coffea species. Based on selected unigenes, it was possible to suggest 26 putative proteins with TE-cassette insertions, demonstrating a likely contribution to protein variability. The genes for two of those proteins, the fertility restorer (FR) and the pyrophosphate-dependent phosphofructokinase (PPi-PFKs) genes, were selected for evaluating the impact of TE-cassettes on host gene evolution of other plant genomes (Arabidopsis thaliana, Oryza sativa and Populus trichocarpa). This survey allowed identifying a FR gene in O. sativa harboring multiple insertions of LTR retrotransposons that originated new exons, which however does not necessarily mean a case of molecular domestication. A possible transduction event of a fragment of the PPi-PFK beta-subunit gene mediated by Helitron ATREPX1 in Arabidopsis thaliana was also highlighted.

Evolution and horizontal transfer of a DD37E DNA transposon in mosquitoes

ITmD37E, a unique class II transposable element (TE) with an ancient origin, appears to have been involved in multiple horizontal transfers in mosquitoes as ITmD37E sequences from 10 mosquito species of five genera share high nucleotide (nt) identities. For example, ITmD37E sequences from Aedes aegypti and Anopheles gambiae, which have an estimated common ancestor of 145-200 million years ago, display 92% nt identity. The comparison of ITmD37E and host mosquito phylogenies shows a lack of congruence. The wide distribution of conserved ITmD37Es in mosquitoes and the presence of intact copies suggest that this element may have been recently active.

Highly similar piggyBac transposase-like sequences from various Bactrocera (Diptera, Tephritidae) species

The piggyBac transposable element is currently the vector of choice for transgenesis, enhancer trapping, gene discovery and gene function determination in both insects and mammals. However, the recent discovery of sequences with similarity to piggyBac in a wide diversity of organisms suggests that piggyBac may be horizontally transferred to distantly related species. This has raised concern on the wide-range application of piggyBac-based transformation vectors and their stability. In this paper, the presence of sequences homologous to the piggyBac transposase was investigated in 17 species belonging to six genera within the Tephritidae family, including many pest species for which transformation has already been achieved. piggyBac-like sequences, with a high degree of similarity to the original Trichoplusia ni transposase sequence were identified only in six species of the Bactrocera genus.

Helitrons on a roll: eukaryotic rolling-circle transposons

Rolling-circle eukaryotic transposons, known as Helitron transposons, were first discovered in plants (Arabidopsis thaliana and Oryza sativa) and in the nematode Caenorhabditis elegans. To date, Helitrons have been identified in a diverse range of species, from protists to mammals. They represent a major class of eukaryotic transposons and are fundamentally different from classical transposons in terms of their structure and mechanism of transposition. Helitrons seem to have a major role in the evolution of host genomes. They frequently capture diverse host genes, some of which can evolve into novel host genes or become essential for helitron transposition.

Characterization of EamaT1, a member of maT family of transposable elements from the earthworm Eisenia andrei (Annelida, Oligochaeta)

The maT family is a unique clade within the Tc1-mariner superfamily, and their distribution is to date known as being limited to invertebrates. A novel transposon named EamaT1 is described from the genome of the earthworm Eisenia andrei. The full sized EamaT1 was obtained by degenerate and inverse PCR-based amplification. Sequence analysis of multiple copies of the EamaT1, which consisted of 0.9 and 1.4 kb elements, showed that the consensual EamaT1 with inverted terminal repeats (ITRs) of 69 bp was 1,422 bp long and flanked by a duplicated TA dinucleotide. The EamaT1 is present in approximately 120-250 copies per diploid genome but undergoes an inactivation process as a result of accumulating multiple mutations and is nonfunctional. The open reading frame (ORF) of the EamaT1 consensus encoding 356 amino acid sequences of transposase contained a DD37D signature and a conserved paired-like DNA binding motif for the transposition mechanism. The result of ITRs comparison confirmed their consensus terminal sequences (5'-CAGGGTG-3') and AT-rich region on the internal bases for ITRs-transposase interaction.

A modified universal fast walking method for single-tube transposon mapping

An enhanced universal fast walking (UFW) method adapted for the mapping of transposons is described. This protocol combines the original UFW method with the use of agarase to unravel composite nucleotide sequence, thereby forgoing molecular cloning steps and the use of restriction enzymes and ligases necessary in other available genome walking methods such as the prominent inverse PCR. The minuscule automatable chemistry of UFW is completed within one reaction vessel using a constant enzyme buffer, and the intrinsic DNA fingerprints, from which amplicons may be quantitatively recovered, offer quality assurance. The core steps of the protocol, spanning half a day or less, comprise first-strand synthesis, primer destruction, random-ended-primer annealing, distal branched-end repair, second-primer destruction, lariat formation and final amplification. Distinctively, no starting or intermediate templates are wasted during the reaction series, thus achieving yields comparable to direct PCR. Ultimate per-reaction walk-lengths are schematically illimitable and sequence-ready amplicons can be produced immediately from prevalent single-copy genomic walk origins. The core UFW protocol may be applied, as described here, to expedited transposon boundary retrieval, but is also applicable to general genome walking and cDNA walking, as well as viral and other insertional element mapping.

Multiple subfamilies of mariner transposable elements are present in stalk-eyed flies (Diptera: Diopsidae)

The Diopsid stalk-eyed flies are an increasingly well-studied group. Presented here is evidence of the first known transposable elements discovered in these flies. The vertumnana mariner subfamily was identified in the Diopsini tribe, but could not be amplified in species of the Sphyracephalini tribe. PCR screening with degenerate primers revealed that multiple mariner subfamilies are present within the Diopsidae. Most of the sequenced elements appear to be pseudogenes; however two subfamilies are shown to be evolving under purifying selection, raising the possibility that mariner is active in some Diopsid species. Evidence is presented of a possible horizontal transfer event involving an unknown Teleopsis species and the Tephritid fly Bactrocera neohumeralis.

Telomere-associated endonuclease-deficient Penelope-like retroelements in diverse eukaryotes

The evolutionary origin of telomerases, enzymes that maintain the ends of linear chromosomes in most eukaryotes, is a subject of debate. Penelope-like elements (PLEs) are a recently described class of eukaryotic retroelements characterized by a GIY-YIG endonuclease domain and by a reverse transcriptase domain with similarity to telomerases and group II introns. Here we report that a subset of PLEs found in bdelloid rotifers, basidiomycete fungi, stramenopiles, and plants, representing four different eukaryotic kingdoms, lack the endonuclease domain and are located at telomeres. The 5' truncated ends of these elements are telomere-oriented and typically capped by species-specific telomeric repeats. Most of them also carry several shorter stretches of telomeric repeats at or near their 3' ends, which could facilitate utilization of the telomeric G-rich 3' overhangs to prime reverse transcription. Many of these telomere-associated PLEs occupy a basal phylogenetic position close to the point of divergence from the telomerase-PLE common ancestor and may descend from the missing link between early eukaryotic retroelements and present-day telomerases.

Spontaneous mutations caused by a Helitron transposon, Hel-It1, in morning glory, Ipomoea tricolor

Helitrons are newcomers among eukaryotic DNA transposons and have originally been identified by computational analysis in the genomes of Arabidopsis, rice and nematode. They are distinguished from other transposons in their structural features, and their proposed transposition mechanisms are involved in rolling circle replication. Computer-predicted autonomous Helitrons with conserved terminal sequences 5'-TC and CTRR-3' are presumed to encode a putative transposase, Rep/Hel-TPase, which contains a characteristic nuclease/ligase domain for the replication-initiation protein (Rep) and a DNA helicase domain (Hel). Plant Helitrons are thought to encode an additional transposase, RPA-TPase, which is related to the largest subunit of the replication protein A (RPA70). Although Helitrons are found in diverse genomes, neither an autonomous element nor a transposition event has been reported. Here we show that a spontaneous pearly-s mutant of Ipomoea tricolor cv. Pearly Gates, exhibiting white flowers and isolated in approximately 1940, has an 11.5-kbp novel Helitron, named Hel-It1, integrated into the DFR-B gene for anthocyanin pigmentation. Hel-It1 shows the predicted plant Helitron structure for an autonomous element with the conserved termini and carrying the two putative transposase genes, Rep/Hel-TPase and RPA-TPase, which contain a nonsense and a frameshift mutation, respectively. Hel-It1-related elements are scattered in the Ipomoea genome, and only a fraction of the pearly-s plants were found to carry Hel-It1 at another insertion site. The pearly-s mutant appears to bear an autonomous element and to express the wild-type RPA-TPase transcripts. The structures of a putative autonomous element and its transposase genes are discussed.

A deep-branching clade of retrovirus-like retrotransposons in bdelloid rotifers

Rotifers of class Bdelloidea, a group of aquatic invertebrates in which males and meiosis have never been documented, are also unusual in their lack of multicopy LINE-like and gypsy-like retrotransposons, groups inhabiting the genomes of nearly all other metazoans. Bdelloids do contain numerous DNA transposons, both intact and decayed, and domesticated Penelope-like retroelements Athena, concentrated at telomeric regions. Here we describe two LTR retrotransposons, each found at low copy number in a different bdelloid species, which define a clade different from previously known clades of LTR retrotransposons. Like bdelloid DNA transposons and Athena, these elements are found preferentially in telomeric regions. Unlike bdelloid DNA transposons, many of which are decayed, the newly described elements, named Vesta and Juno, inhabiting the genomes of Philodina roseola and Adineta vaga, respectively, appear to be intact and represent recent insertions, possibly from an exogenous source. We describe the retrovirus-like structure of the new elements, containing gag, pol, and env-like open reading frames, and discuss their possible origins, transmission, and behavior in bdelloid genomes.

Mutator phenotype of Caenorhabditis elegans DNA damage checkpoint mutants

DNA damage response proteins identify sites of DNA damage and signal to downstream effectors that orchestrate either apoptosis or arrest of the cell cycle and DNA repair. The C. elegans DNA damage response mutants mrt-2, hus-1, and clk-2(mn159) displayed 8- to 15-fold increases in the frequency of spontaneous mutation in their germlines. Many of these mutations were small- to medium-sized deletions, some of which had unusual sequences at their breakpoints such as purine-rich tracts or direct or inverted repeats. Although DNA-damage-induced apoptosis is abrogated in the mrt-2, hus-1, and clk-2 mutant backgrounds, lack of the apoptotic branch of the DNA damage response pathway in cep-1/p53, ced-3, and ced-4 mutants did not result in a Mutator phenotype. Thus, DNA damage checkpoint proteins suppress the frequency of mutation by ensuring that spontaneous DNA damage is accurately repaired in C. elegans germ cells. Although DNA damage response defects that predispose humans to cancer are known to result in large-scale chromosome aberrations, our results suggest that small- to medium-sized deletions may also play roles in the development of cancer.

Assembly of the Tc1 and mariner transposition initiation complexes depends on the origins of their transposase DNA binding domains

In this review, we focus on the assembly of DNA/protein complexes that trigger transposition in eukaryotic members of the IS630-Tc1-mariner (ITm) super-family, the Tc1- and mariner-like elements (TLEs and MLEs). Elements belonging to this super-family encode transposases with DNA binding domains of different origins, and recent data indicate that the chimerization of functional domains has been an important evolutionary aspect in the generation of new transposons within the ITm super-family. These data also reveal that the inverted terminal repeats (ITRs) at the ends of transposons contain three kinds of motif within their sequences. The first two are well known and correspond to the cleavage site on the outer ITR extremities, and the transposase DNA binding site. The organization of ITRs and of the transposase DNA binding domains implies that differing pathways are used by MLEs and TLEs to regulate transposition initiation. These differences imply that the ways ITRs are recognized also differ leading to the formation of differently organized synaptic complexes. The third kind of motif is the transposition enhancers, which have been found in almost all the functional MLEs and TLEs analyzed to date. Finally, in vitro and in vivo assays of various elements all suggest that the transposition initiation complex is not formed randomly, but involves a mechanism of oriented transposon scanning.