The esterase and PHD domains in CR1-like non-LTR retrotransposons

The Chlamydomonas Genome Project, version 6: Reference assemblies for mating-type plus and minus strains reveal extensive structural mutation in the laboratory

Five versions of the Chlamydomonas reinhardtii reference genome have been produced over the last two decades. Here we present version 6, bringing significant advances in assembly quality and structural annotations. PacBio-based chromosome-level assemblies for two laboratory strains, CC-503 and CC-4532, provide resources for the plus and minus mating-type alleles. We corrected major misassemblies in previous versions and validated our assemblies via linkage analyses. Contiguity increased over ten-fold and >80% of filled gaps are within genes. We used Iso-Seq and deep RNA-seq datasets to improve structural annotations, and updated gene symbols and textual annotation of functionally characterized genes via extensive manual curation. We discovered that the cell wall-less classical reference strain CC-503 exhibits genomic instability potentially caused by deletion of the helicase RECQ3, with major structural mutations identified that affect >100 genes. We therefore present the CC-4532 assembly as the primary reference, although this strain also carries unique structural mutations and is experiencing rapid proliferation of a Gypsy retrotransposon. We expect all laboratory strains to harbor gene-disrupting mutations, which should be considered when interpreting and comparing experimental results. Collectively, the resources presented here herald a new era of Chlamydomonas genomics and will provide the foundation for continued research in this important reference organism.

An ancient clade of Penelope-like retroelements with permuted domains is present in the green lineage and protists, and dominates many invertebrate genomes

Penelope-like elements (PLEs) are an enigmatic clade of retrotransposons whose reverse transcriptases (RTs) share a most recent common ancestor with telomerase RTs. The single ORF of canonical endonuclease (EN)+ PLEs encodes RT and a C-terminal GIY–YIG EN that enables intrachromosomal integration, whereas EN− PLEs lack EN and are generally restricted to chromosome termini. EN+ PLEs have only been found in animals, except for one case of horizontal transfer to conifers, whereas EN− PLEs occur in several kingdoms. Here, we report a new, deep-branching PLE clade with a permuted domain order, whereby an N-terminal GIY–YIG EN is linked to a C-terminal RT by a short domain with a characteristic CxC motif. These N-terminal EN+ PLEs share a structural organization, including pseudo-LTRs and complex tandem/inverted insertions, with canonical EN+ PLEs from Penelope/Poseidon, Neptune, and Nematis clades, and show insertion bias for microsatellites, but lack canonical hammerhead ribozyme motifs. However, their phylogenetic distribution is much broader. The Naiads, found in numerous invertebrate phyla, can reach tens of thousands of copies per genome. In spiders and clams, Naiads independently evolved to encode selenoproteins containing multiple selenocysteines. Chlamys, which lack the CCHH motif universal to PLE ENs, occur in green algae, spike mosses (targeting ribosomal DNA), and slime molds. Unlike canonical PLEs, RTs of N-terminal EN+ PLEs contain the insertion-in-fingers domain (IFD), strengthening the link between PLEs and telomerases. Additionally, we describe Hydra, a novel metazoan C-terminal EN+ clade. Overall, we conclude that PLE diversity, taxonomic distribution, and abundance are comparable with non-LTR and LTR-retrotransposons.

Transposons-Based Clonal Diversity in Trematode Involves Parts of CR1 (LINE) in Eu- and Heterochromatin

Trematode parthenitae have long been believed to form clonal populations, but clonal diversity has been discovered in this asexual stage of the lifecycle. Clonal polymorphism in the model species Himasthla elongata has been previously described, but the source of this phenomenon remains unknown. In this work, we traced cercarial clonal diversity using a simplified amplified fragment length polymorphism (SAFLP) method and characterised the nature of fragments in diverse electrophoretic bands. The repetitive elements were identified in both the primary sequence of the H. elongata genome and in the transcriptome data. Long-interspersed nuclear elements (LINEs) and long terminal repeat retrotransposons (LTRs) were found to represent an overwhelming majority of the genome and the transposon transcripts. Most sequenced fragments from SAFLP pattern contained the reverse transcriptase (RT, ORF2) domains of LINEs, and only a few sequences belonged to ORFs of LTRs and ORF1 of LINEs. A fragment corresponding to a CR1-like (LINE) spacer region was discovered and named CR1-renegade (CR1-rng). In addition to RT-containing CR1 transcripts, we found short CR1-rng transcripts in the redia transcriptome and short contigs in the mobilome. Probes against CR1-RT and CR1-rng presented strikingly different pictures in FISH mapping, despite both being fragments of CR1. In silico data and Southern blotting indicated that CR1-rng is not tandemly organised. CR1 involvement in clonal diversity is discussed.

The substrate selectivity of the two homologous SGNH hydrolases from Streptomyces bacteria: Molecular dynamics and experimental study

Two extracellular enzymes of the SGNH hydrolase superfamily reveal highly homologous 3D structures, but act on different substrates; one is a true phospholipase A1 from Streptomyces albidoflavus (SaPLA1, EC: 3.1.1.32, PDB code: 4HYQ), whereas the promiscuous enzyme from Streptomyces rimosus (SrLip, EC: 3.1.1.3, PDB code: 5MAL) exhibits lipase, phospholipase, esterase, thioesterase, and Tweenase activities. To get insight into binding modes of phospholipid and triglyceride substrates in both enzymes and understand their chain-length preferences, we opted for computational approach based on in silico prepared enzyme-substrate complexes. Docking procedure and molecular dynamics simulations at microsecond time scale were applied. The modelled complexes of SaPLA1 and SrLip enzymes revealed substrate accommodation: a) the acyl-chain attached to sn-1 position fits into the hydrophobic pocket, b) the acyl-chain attached to sn-2 position fits in the hydrophobic cleft, whereas c) the sn-3 bound acyl chain of the triglyceride or polar head of the glycerophospholipid fits into the binding groove. Moreover, our results pinpointed subtle amino acid differences in the hydrophobic pockets of these two enzymes which accommodate the acyl chain attached to sn-1 position of glycerol to be responsible for the chain length preference. Slight differences in the binding grooves of SaPLA1 and SrLip, which accommodate the acyl chain attached to sn-3 position are responsible for exclusive phospholipase and both phospholipase/lipase activities of these two enzymes, respectively. The results of modelling correlate with the experimentally obtained kinetic parameters given in the literature and are important for protein engineering that aims to obtain a variant of enzyme, which would preferably act on the substrate of interest.

Successful Invasions of Short Internally Deleted Elements (SIDEs) and Its Partner CR1 in Lepidoptera Insects

Although DNA transposons often generated internal deleted derivatives such as miniature inverted-repeat transposable elements, short internally deleted elements (SIDEs) derived from nonlong terminal-repeat retrotransposons are rare. Here, we found a novel SIDE, named Persaeus, that originated from the chicken repeat 1 (CR1) retrotransposon Zenon and it has been found widespread in Lepidoptera insects. Our findings suggested that Persaeus and the partner Zenon have experienced a transposition burst in their host genomes and the copy number of Persaeus and Zenon in assayed genomes are significantly correlated. Accordingly, the activity though age analysis indicated that the replication wave of Persaeus coincided with that of Zenon. Phylogenetic analyses suggested that Persaeus may have evolved at least four times independently, and that it has been vertically transferred into its host genomes. Together, our results provide new insights into the evolution dynamics of SIDEs and its partner non-LTRs.

Capture of a functionally active methyl-CpG binding domain by an arthropod retrotransposon family

The repressive capacity of cytosine DNA methylation is mediated by recruitment of silencing complexes by methyl-CpG binding domain (MBD) proteins. Despite MBD proteins being associated with silencing, we discovered that a family of arthropod Copia retrotransposons have incorporated a host-derived MBD. We functionally show how retrotransposon-encoded MBDs preferentially bind to CpG-dense methylated regions, which correspond to transposable element regions of the host genome, in the myriapod Strigamia maritima Consistently, young MBD-encoding Copia retrotransposons (CopiaMBD) accumulate in regions with higher CpG densities than other LTR-retrotransposons also present in the genome. This would suggest that retrotransposons use MBDs to integrate into heterochromatic regions in Strigamia, avoiding potentially harmful insertions into host genes. In contrast, CopiaMBD insertions in the spider Stegodyphus dumicola genome disproportionately accumulate in methylated gene bodies compared with other spider LTR-retrotransposons. Given that transposons are not actively targeted by DNA methylation in the spider genome, this distribution bias would also support a role for MBDs in the integration process. Together, these data show that retrotransposons can co-opt host-derived epigenome readers, potentially harnessing the host epigenome landscape to advantageously tune the retrotransposition process.

Distribution of CR1-like transposable element in woodpeckers (Aves Piciformes): Z sex chromosomes can act as a refuge for transposable elements

Birds have relatively few repetitive sequences compared to other groups of vertebrates; however, the members of order Piciformes (woodpeckers) have more of these sequences, composed mainly of transposable elements (TE). The TE most often found in birds is a retrotransposon chicken repeat 1 (CR1). Piciformes lineages were subjected to an expansion of the CR1 elements, carrying a larger fraction of transposable elements. This study compared patterns of chromosome distribution among five bird species, through chromosome mapping of the CR1 sequence and reconstructed their phylogenetic tree. We analyzed several members of Piciformes (Colaptes campestris, Colaptes melanochloros, Melanerpes candidus, and Veniliornis spilogaster), as well as Galliformes (Gallus gallus). Gallus gallus is the species with which most genomic and hence cytogenetic studies have been performed. The results showed that CR1 sequences are a monophyletic group and do not depend on their hosts. All species analyzed showed a hybridization signal by fluorescence in situ hybridization (FISH). In all species, the chromosomal distribution of CR1 was not restricted to heterochromatin regions in the macrochromosomes, principally pair 1 and the Z sex chromosome. Accumulation in the Z sex chromosomes can serve as a refuge for transposable elements. These results highlight the importance of transposable elements in host genomes and karyotype evolution.

Human LINE-1 retrotransposition requires a metastable coiled coil and a positively charged N-terminus in L1ORF1p

LINE-1 (L1) is an autonomous retrotransposon, which acted throughout mammalian evolution and keeps contributing to human genotypic diversity, genetic disease and cancer. L1 encodes two essential proteins: L1ORF1p, a unique RNA-binding protein, and L1ORF2p, an endonuclease and reverse transcriptase. L1ORF1p contains an essential, but rapidly evolving N-terminal portion, homo-trimerizes via a coiled coil and packages L1RNA into large assemblies. Here, we determined crystal structures of the entire coiled coil domain of human L1ORF1p. We show that retrotransposition requires a non-ideal and metastable coiled coil structure, and a strongly basic L1ORF1p amino terminus. Human L1ORF1p therefore emerges as a highly calibrated molecular machine, sensitive to mutation but functional in different hosts. Our analysis rationalizes the locally rapid L1ORF1p sequence evolution and reveals striking mechanistic parallels to coiled coil-containing membrane fusion proteins. It also suggests how trimeric L1ORF1p could form larger meshworks and indicates critical novel steps in L1 retrotransposition.

Almost half of the human genome consists of DNA strings that have been copied and pasted from one part of the genome to another many thousands of times. These strings of DNA are called mobile genetic elements. Mobile elements can disrupt important genes, causing disease and cancer, but they can also drive evolution.

Presently, only one type of mobile element, called LINE-1, is active in the human genome and able to multiply without help from other mobile elements. LINE-1 DNA is ‘transcribed’ to form molecules of LINE-1 RNA, which can then be ‘translated’ into two distinct proteins. These bind to LINE-1 RNA, which then gets back-transcribed into DNA and inserted as a new LINE-1 element in a new region of the genome. One of the two proteins, called L1ORF1p, forms complexes where three copies of the protein come together. These ‘trimers’ cover and protect LINE-1 RNA and are required for LINE-1 mobility.

Different versions of L1ORF1p are found in different animals. Part of the protein is the same across all mammals, and this ‘conserved’ part controls the ability of L1ORF1p to bind to RNA. The non-conserved part of L1ORF1p differs even between humans and their closest animal relatives and little was known about its structure or role. However, this rapidly evolving part of L1ORF1p is essential for LINE-1 mobility.

Using X-ray crystallography, Khazina and Weichenrieder obtained a molecular snapshot of the part of L1ORF1p that interacts with other copies of the protein to form trimers. Combined with earlier snapshots of L1ORF1p’s conserved part, this generated a complete structural model of the L1ORF1p trimer. Additional biophysical characterizations suggest that L1ORF1p trimers form a semi-stable structure that can partially open up, indicating how trimers could form larger assemblies of L1ORF1p on LINE-1 RNA. Indeed, the need to maintain a semi-stable structure could explain why L1ORF1p is evolving so rapidly. A second important finding is that the beginning of L1ORF1p needs to be positively charged – a requirement that warrants further exploration.

The structural and mechanistic insight into L1ORF1p points to critical new steps in LINE-1 mobilization. It will help to design inhibitor molecules with the goal to halt the mobilization process at various points and to dissect such steps in great detail. Understanding how to control LINE-1 mobility could help to improve stem cell therapies and reproduction assistance techniques, due to the fact that LINE-1 mobility is a potential source of mutation in stem cells, egg and sperm cells, and newly formed embryos.

Giant Reverse Transcriptase-Encoding Transposable Elements at Telomeres

Transposable elements are omnipresent in eukaryotic genomes and have a profound impact on chromosome structure, function and evolution. Their structural and functional diversity is thought to be reasonably well-understood, especially in retroelements, which transpose via an RNA intermediate copied into cDNA by the element-encoded reverse transcriptase, and are characterized by a compact structure. Here, we report a novel type of expandable eukaryotic retroelements, which we call Terminons. These elements can attach to G-rich telomeric repeat overhangs at the chromosome ends, in a process apparently facilitated by complementary C-rich repeats at the 3′-end of the RNA template immediately adjacent to a hammerhead ribozyme motif. Terminon units, which can exceed 40 kb in length, display an unusually complex and diverse structure, and can form very long chains, with host genes often captured between units. As the principal polymerizing component, Terminons contain Athena reverse transcriptases previously described in bdelloid rotifers and belonging to the enigmatic group of Penelope-like elements, but can additionally accumulate multiple cooriented ORFs, including DEDDy 3′-exonucleases, GDSL esterases/lipases, GIY-YIG-like endonucleases, rolling-circle replication initiator (Rep) proteins, and putatively structural ORFs with coiled-coil motifs and transmembrane domains. The extraordinary length and complexity of Terminons and the high degree of interfamily variability in their ORF content challenge the current views on the structural organization of eukaryotic retroelements, and highlight their possible connections with the viral world and the implications for the elevated frequency of gene transfer.

A conserved role for the ESCRT membrane budding complex in LINE retrotransposition

Long interspersed nuclear element-1s (LINE-1s, or L1s) are an active family of retrotransposable elements that continue to mutate mammalian genomes. Despite the large contribution of L1 to mammalian genome evolution, we do not know where active L1 particles (particles in the process of retrotransposition) are located in the cell, or how they move towards the nucleus, the site of L1 reverse transcription. Using a yeast model of LINE retrotransposition, we identified ESCRT (endosomal sorting complex required for transport) as a critical complex for LINE retrotransposition, and verified that this interaction is conserved for human L1. ESCRT interacts with L1 via a late domain motif, and this interaction facilitates L1 replication. Loss of the L1/ESCRT interaction does not impair RNP formation or enzymatic activity, but leads to loss of retrotransposition and reduced L1 endonuclease activity in the nucleus. This study highlights the importance of the ESCRT complex in the L1 life cycle and suggests an unusual mode for L1 RNP trafficking.

Long interspersed nuclear elements (LINEs) are a class of retrotransposable elements that mutate mammalian genomes. LINEs have been highly successful in the human genome, multiplying to over 800,000 copies. The LINE-encoded replication machinery is also used by other retrotransposons, and in total, has been responsible for the generation of over 1/3 of human DNA sequence. To replicate, a LINE mRNA forms a ribonucleoprotein particle (RNP) with its proteins. This RNP eventually enters the nucleus to integrate a cDNA copy of itself into chromosomes. The events between RNP formation and successful integration are difficult to study and largely unknown. Here we show that the ESCRT complex plays a conserved role in LINE retrotransposition in both yeast and humans. ESCRT is a membrane budding complex involved in cellular trafficking and membrane budding/fusion. Our results imply that membranes play an integral part of LINE replication, and ESCRT may be required for RNP trafficking towards the nucleus.

LTR-Retrotransposons from Bdelloid Rotifers Capture Additional ORFs Shared between Highly Diverse Retroelement Types

Rotifers of the class Bdelloidea, microscopic freshwater invertebrates, possess a highlydiversified repertoire of transposon families, which, however, occupy less than 4% of genomic DNA in the sequenced representative Adineta vaga. We performed a comprehensive analysis of A. vaga retroelements, and found that bdelloid long terminal repeat (LTR)retrotransposons, in addition to conserved open reading frame (ORF) 1 and ORF2 corresponding to gag and pol genes, code for an unusually high variety of ORF3 sequences. Retrovirus-like LTR families in A. vaga belong to four major lineages, three of which are rotiferspecific and encode a dUTPase domain. However only one lineage contains a canonical envlike fusion glycoprotein acquired from paramyxoviruses (non-segmented negative-strand RNA viruses), although smaller ORFs with transmembrane domains may perform similar roles. A different ORF3 type encodes a GDSL esterase/lipase, which was previously identified as ORF1 in several clades of non-LTR retrotransposons, and implicated in membrane targeting. Yet another ORF3 type appears in unrelated LTR-retrotransposon lineages, and displays strong homology to DEDDy-type exonucleases involved in 3'-end processing of RNA and single-stranded DNA. Unexpectedly, each of the enzymatic ORF3s is also associated with different subsets of Penelope-like Athena retroelement families. The unusual association of the same ORF types with retroelements from different classes reflects their modular structure with a high degree of flexibility, and points to gene sharing between different groups of retroelements.

LINEs between Species: Evolutionary Dynamics of LINE-1 Retrotransposons across the Eukaryotic Tree of Life

LINE-1 (L1) retrotransposons are dynamic elements. They have the potential to cause great genomic change because of their ability to ‘jump’ around the genome and amplify themselves, resulting in the duplication and rearrangement of regulatory DNA. Active L1, in particular, are often thought of as tightly constrained, homologous and ubiquitous elements with well-characterized domain organization. For the past 30 years, model organisms have been used to define L1s as 6–8 kb sequences containing a 5′-UTR, two open reading frames working harmoniously in cis, and a 3′-UTR with a polyA tail. In this study, we demonstrate the remarkable and overlooked diversity of L1s via a comprehensive phylogenetic analysis of elements from over 500 species from widely divergent branches of the tree of life. The rapid and recent growth of L1 elements in mammalian species is juxtaposed against the diverse lineages found in other metazoans and plants. In fact, some of these previously unexplored mammalian species (e.g. snub-nosed monkey, minke whale) exhibit L1 retrotranspositional ‘hyperactivity’ far surpassing that of human or mouse. In contrast, non-mammalian L1s have become so varied that the current classification system seems to inadequately capture their structural characteristics. Our findings illustrate how both long-term inherited evolutionary patterns and random bursts of activity in individual species can significantly alter genomes, highlighting the importance of L1 dynamics in eukaryotes.

Evolutionary Dynamics of 5S rDNA and Recurrent Association of Transposable Elements in Electric Fish of the Family Gymnotidae (Gymnotiformes): The Case of Gymnotus mamiraua

Gymnotidae is a family of electric fish endemic to the Neotropics consisting of 2 genera: Electrophorus and Gymnotus. The genus Gymnotus is widely distributed and is found in all of the major Brazilian river systems. Physical and molecular mapping data for the ribosomal DNA (rDNA) in this genus are still scarce, with its chromosomal location known in only 11 species. As other species of Gymnotus with 2n = 54 chromosomes from the Paraná-Paraguay basin, G. mamiraua was found to have a large number of 5S rDNA sites. Isolation and cloning of the 5S rDNA sequences from G. mamiraua identified a fragment of a transposable element similar to the Tc1/mariner transposon associated with a non-transcribed spacer. Double fluorescence in situ hybridization analysis of this element and the 5S rDNA showed that they were colocalized on several chromosomes, in addition to acting as nonsyntenic markers on others. Our data show the association between these sequences and suggest that the Tc1 retrotransposon may be the agent that drives the spread of these 5S rDNA-like sequences in the G. mamiraua genome.

Distribution patterns and impact of transposable elements in genes of green algae

Transposable elements (TEs) are DNA sequences able to transpose in the host genome, a remarkable feature that enables them to influence evolutive trajectories of species. An investigation about the TE distribution and TE impact in different gene regions of the green algae species Chlamydomonas reinhardtii and Volvox carteri was performed. Our results indicate that TEs are very scarce near introns boundaries, suggesting that insertions in this region are negatively selected. This contrasts with previous results showing enrichment of tandem repeats in introns boundaries and suggests that different evolutionary forces are acting in these different classes of repeats. Despite the relatively low abundance of TEs in the genome of green algae when compared to mammals, the proportion of poly(A) sites derived from TEs found in C. reinhardtii was similar to that described in human and mice. This fact, associated with the enrichment of TEs in gene 5' and 3' flanks of C. reinhardtii, opens up the possibility that TEs may have considerably contributed for gene regulatory sequences evolution in this species. Moreover, it was possible identify several instances of TE exonization for C. reinhardtii, with a particularly interesting case from a gene coding for Condensin II, a protein involved in the maintenance of chromosomal structure, where the addition of a transposomal PHD finger may contribute to binding specificity of this protein. Taken together, our results suggest that the low abundance of TEs in green algae genomes is correlated with a strict negative selection process, combined with the retention of copies that contribute positively with gene structures.

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.

Multiple lineages of ancient CR1 retroposons shaped the early genome evolution of amniotes

Chicken repeat 1 (CR1) retroposons are long interspersed elements (LINEs) that are ubiquitous within amniote genomes and constitute the most abundant family of transposed elements in birds, crocodilians, turtles, and snakes. They are also present in mammalian genomes, where they reside as numerous relics of ancient retroposition events. Yet, despite their relevance for understanding amniote genome evolution, the diversity and evolution of CR1 elements has never been studied on an amniote-wide level. We reconstruct the temporal and quantitative activity of CR1 subfamilies via presence/absence analyses across crocodilian phylogeny and comparative analyses of 12 crocodilian genomes, revealing relative genomic stasis of retroposition during genome evolution of extant Crocodylia. Our large-scale phylogenetic analysis of amniote CR1 subfamilies suggests the presence of at least seven ancient CR1 lineages in the amniote ancestor; and amniote-wide analyses of CR1 successions and quantities reveal differential retention (presence of ancient relics or recent activity) of these CR1 lineages across amniote genome evolution. Interestingly, birds and lepidosaurs retained the fewest ancient CR1 lineages among amniotes and also exhibit smaller genome sizes. Our study is the first to analyze CR1 evolution in a genome-wide and amniote-wide context and the data strongly suggest that the ancestral amniote genome contained myriad CR1 elements from multiple ancient lineages, and remnants of these are still detectable in the relatively stable genomes of crocodilians and turtles. Early mammalian genome evolution was thus characterized by a drastic shift from CR1 prevalence to dominance and hyperactivity of L2 LINEs in monotremes and L1 LINEs in therians.

Interspecies insertion polymorphism analysis reveals recent activity of transposable elements in extant coelacanths

Coelacanths are lobe-finned fish represented by two extant species, Latimeria chalumnae in South Africa and Comoros and L. menadoensis in Indonesia. Due to their intermediate phylogenetic position between ray-finned fish and tetrapods in the vertebrate lineage, they are of great interest from an evolutionary point of view. In addition, extant specimens look similar to 300 million-year-old fossils; because of their apparent slowly evolving morphology, coelacanths have been often described as « living fossils ». As an underlying cause of such a morphological stasis, several authors have proposed a slow evolution of the coelacanth genome. Accordingly, sequencing of the L. chalumnae genome has revealed a globally low substitution rate for protein-coding regions compared to other vertebrates. However, genome and gene evolution can also be influenced by transposable elements, which form a major and dynamic part of vertebrate genomes through their ability to move, duplicate and recombine. In this work, we have searched for evidence of transposition activity in coelacanth genomes through the comparative analysis of orthologous genomic regions from both Latimeria species. Comparison of 5.7 Mb (0.2%) of the L. chalumnae genome with orthologous Bacterial Artificial Chromosome clones from L. menadoensis allowed the identification of 27 species-specific transposable element insertions, with a strong relative contribution of CR1 non-LTR retrotransposons. Species-specific homologous recombination between the long terminal repeats of a new coelacanth endogenous retrovirus was also detected. Our analysis suggests that transposon activity is responsible for at least 0.6% of genome divergence between both Latimeria species. Taken together, this study demonstrates that coelacanth genomes are not evolutionary inert: they contain recently active transposable elements, which have significantly contributed to post-speciation genome divergence in Latimeria.

Integrated mechanism for the generation of the 5' junctions of LINE inserts

To elucidate the molecular mechanism of the integration of long interspersed elements (LINEs), we characterized the 5′ ends of more than 200 LINE de novo retrotransposition events into chicken DT40 or human HeLa cells. Human L1 inserts produced 15-bp target-site duplications (TSDs) and zebrafish ZfL2-1 inserts produced 5-bp TSDs in DT40 cells, suggesting that TSD length depends on the LINE species. Further analysis of 5′ junctions revealed that the 5′-end-joining pathways of LINEs can be divided into two fundamental types—annealing or direct. We also found that the generation of 5′ inversions depends on host and LINE species. These results led us to propose a new model for 5′-end joining, the type of which is determined by the extent of exposure of 3′ overhangs generated after the second-strand cleavage and by the involvement of host factors.

Modular organization and reticulate evolution of the ORF1 of Jockey superfamily transposable elements

Background

Long interspersed nuclear elements (LINES) are the most common transposable element (TE) in almost all metazoan genomes examined. In most LINE superfamilies there are two open reading frames (ORFs), and both are required for transposition. The ORF2 is well characterized, while the structure and function of the ORF1 is less well understood. ORF1s have been classified into five types based on structural organization and the domains identified. Here we perform a large scale analysis of ORF1 domains of 448 elements from the Jockey superfamily using multiple alignments and Hidden Markov Model (HMM)-HMM comparisons.

Results

Three major lineages, Chicken repeat 1 (CR1), LINE2 (L2) and Jockey, were identified. All Jockey lineage elements have the same type of ORF1. In contrast, in the L2 and CR1 lineage elements, all five ORF1 types are found, with no one type of ORF1 predominating. A plant homeodomain (PHD) is much more prevalent than previously suspected. ORF1 type variations involving the PHD domain were found in many subgroups of the L2 and CR1 lineages. A Jockey lineage-like ORF1 with a PHD domain was found in both lineages. A phylogenetic analysis of this ORF1 suggests that it has been horizontally transferred. Likewise, an esterase containing ORF1 type was only found in two exclusively vertebrate L2 and CR1 groups, indicating that it may have been acquired in a vertebrate common ancestor and then transferred between the lineages.

Conclusions

The ORF1 of the CR1 and L2 lineages is very structurally diverse. The presence of a PHD domain in many ORF1s of the L2 and CR1 lineages is suggestive of domain shuffling. There is also evidence of possible horizontal transfer of entire ORF1s between lineages. In conclusion, while the structure of the ORF2 appears to be highly constrained and its evolution tree-like, the structure of the ORF1 within the CR1 and L2 lineages is much more variable and its evolution reticulate.

A multicopy Y-chromosomal SGNH hydrolase gene expressed in the testis of the platyfish has been captured and mobilized by a Helitron transposon

Background

Teleost fish present a high diversity of sex determination systems, with possible frequent evolutionary turnover of sex chromosomes and sex-determining genes. In order to identify genes involved in male sex determination and differentiation in the platyfish Xiphophorus maculatus, bacterial artificial chromosome contigs from the sex-determining region differentiating the Y from the X chromosome have been assembled and analyzed.

Results

A novel three-copy gene called teximY (for testis-expressed in Xiphophorus maculatus on the Y) was identified on the Y but not on the X chromosome. A highly related sequence called texim1, probably at the origin of the Y-linked genes, as well as three more divergent texim genes were detected in (pseudo)autosomal regions of the platyfish genome. Texim genes, for which no functional data are available so far in any organism, encode predicted esterases/lipases with a SGNH hydrolase domain. Texim proteins are related to proteins from very different origins, including proteins encoded by animal CR1 retrotransposons, animal platelet-activating factor acetylhydrolases (PAFah) and bacterial hydrolases. Texim gene distribution is patchy in animals. Texim sequences were detected in several fish species including killifish, medaka, pufferfish, sea bass, cod and gar, but not in zebrafish. Texim-like genes are also present in Oikopleura (urochordate), Amphioxus (cephalochordate) and sea urchin (echinoderm) but absent from mammals and other tetrapods. Interestingly, texim genes are associated with a Helitron transposon in different fish species but not in urochordates, cephalochordates and echinoderms, suggesting capture and mobilization of an ancestral texim gene in the bony fish lineage. RT-qPCR analyses showed that Y-linked teximY genes are preferentially expressed in testis, with expression at late stages of spermatogenesis (late spermatids and spermatozeugmata).

Conclusions

These observations suggest either that TeximY proteins play a role in Helitron transposition in the male germ line in fish, or that texim genes are spermatogenesis genes mobilized and spread by transposable elements in fish genomes.

Structure and properties of the esterase from non-LTR retrotransposons suggest a role for lipids in retrotransposition

Non-LTR retrotransposons are mobile genetic elements and play a major role in eukaryotic genome evolution and disease. Similar to retroviruses they encode a reverse transcriptase, but their genomic integration mechanism is fundamentally different, and they lack homologs of the retroviral nucleocapsid-forming protein Gag. Instead, their first open reading frames encode distinct multi-domain proteins (ORF1ps) presumed to package the retrotransposon-encoded RNA into ribonucleoprotein particles (RNPs). The mechanistic roles of ORF1ps are poorly understood, particularly of ORF1ps that appear to harbor an enzymatic function in the form of an SGNH-type lipolytic acetylesterase. We determined the crystal structures of the coiled coil and esterase domains of the ORF1p from the Danio rerio ZfL2-1 element. We demonstrate a dimerization of the coiled coil and a hydrolytic activity of the esterase. Furthermore, the esterase binds negatively charged phospholipids and liposomes, but not oligo-(A) RNA. Unexpectedly, the esterase can split into two dynamic half-domains, suited to engulf long fatty acid substrates extending from the active site. These properties indicate a role for lipids and membranes in non-LTR retrotransposition. We speculate that Gag-like membrane targeting properties of ORF1ps could play a role in RNP assembly and in membrane-dependent transport or localization processes.

RNA-Mediated Gene Duplication and Retroposons: Retrogenes, LINEs, SINEs, and Sequence Specificity

A substantial number of “retrogenes” that are derived from the mRNA of various intron-containing genes have been reported. A class of mammalian retroposons, long interspersed element-1 (LINE1, L1), has been shown to be involved in the reverse transcription of retrogenes (or processed pseudogenes) and non-autonomous short interspersed elements (SINEs). The 3′-end sequences of various SINEs originated from a corresponding LINE. As the 3′-untranslated regions of several LINEs are essential for retroposition, these LINEs presumably require “stringent” recognition of the 3′-end sequence of the RNA template. However, the 3′-ends of mammalian L1s do not exhibit any similarity to SINEs, except for the presence of 3′-poly(A) repeats. Since the 3′-poly(A) repeats of L1 and Alu SINE are critical for their retroposition, L1 probably recognizes the poly(A) repeats, thereby mobilizing not only Alu SINE but also cytosolic mRNA. Many flowering plants only harbor L1-clade LINEs and a significant number of SINEs with poly(A) repeats, but no homology to the LINEs. Moreover, processed pseudogenes have also been found in flowering plants. I propose that the ancestral L1-clade LINE in the common ancestor of green plants may have recognized a specific RNA template, with stringent recognition then becoming relaxed during the course of plant evolution.

Structural and functional characterisation of TesA - a novel lysophospholipase A from Pseudomonas aeruginosa

TesA from Pseudomonas aeruginosa belongs to the GDSL hydrolase family of serine esterases and lipases that possess a broad substrate- and regiospecificity. It shows high sequence homology to TAP, a multifunctional enzyme from Escherichia coli exhibiting thioesterase, lysophospholipase A, protease and arylesterase activities. Recently, we demonstrated high arylesterase activity for TesA, but only minor thioesterase and no protease activity. Here, we present a comparative analysis of TesA and TAP at the structural, biochemical and physiological levels. The crystal structure of TesA was determined at 1.9 Å and structural differences were identified, providing a possible explanation for the differences in substrate specificities. The comparison of TesA with other GDSL-hydrolase structures revealed that the flexibility of active-site loops significantly affects their substrate specificity. This assumption was tested using a rational approach: we have engineered the putative coenzyme A thioester binding site of E. coli TAP into TesA of P. aeruginosa by introducing mutations D17S and L162R. This TesA variant showed increased thioesterase activity comparable to that of TAP. TesA is the first lysophospholipase A described for the opportunistic human pathogen P. aeruginosa. The enzyme is localized in the periplasm and may exert important functions in the homeostasis of phospholipids or detoxification of lysophospholipids.

Evolution of the Australian lungfish (Neoceratodus forsteri) genome: a major role for CR1 and L2 LINE elements

Haploid genomes greater than 25,000 Mb are rare, within the animals only the lungfish and some of the salamanders and crustaceans are known to have genomes this large. There is very little data on the structure of genomes this size. It is known, however, that for animal genomes up to 3,000 Mb, there is in general a good correlation between genome size and the percent of the genome composed of repetitive sequence and that this repetitive component is highly dynamic. In this study, we sampled the Australian lungfish genome using three mini-genomic libraries and found that with very little sequence, the results converged on an estimate of 40% of the genome being composed of recognizable transposable elements (TEs), chiefly from the CR1 and L2 long interspersed nuclear element clades. We further characterized the CR1 and L2 elements in the lungfish genome and show that although most CR1 elements probably represent recent amplifications, the L2 elements are more diverse and are more likely the result of a series of amplifications. We suggest that our sampling method has probably underestimated the recognizable TE content. However, on the basis of the most likely sources of error, we suggest that this very large genome is not largely composed of recently amplified, undetected TEs but may instead include a large component of older degenerate TEs. Based on these estimates, and on Thomson's (Thomson K. 1972. An attempt to reconstruct evolutionary changes in the cellular DNA content of lungfish. J Exp Zool. 180:363-372) inference that in the lineage leading to the extant Australian lungfish, there was massive increase in genome size between 350 and 200 mya, after which the size of the genome changed little, we speculate that the very large Australian lungfish genome may be the result of a massive amplification of TEs followed by a long period with a very low rate of sequence removal and some ongoing TE activity.

CR1 retroposons provide a new insight into the phylogeny of Phasianidae species (Aves: Galliformes)

Chicken repeat 1 (CR1) elements, a class of retroposons belonging to non-long-terminal repeats, have been recognized as powerful tools for phylogenetic studies. Here we examine the phylogenetic relationships of 11 Phasianidae species based on CR1 retroposons. Together with 19 loci reported previously, a total of 99 CR1 loci were identified from chicken genome and turkey BAC clone sequences. 75 insertion events were used to address the branching order of 11 species in Phasianidae. The topology of our tree suggests that: 1) Gallus gallus possessed a basal phylogenetic position within Phasianidae and was related to Bambusicola thoracica (BSP=100%); 2) After the split of G. gallus and B. thoracica, Arborophila rufipectus diverged from Phasianidae (BSP=100%). Nine unambiguous insertion events supported a phylogenetic position of A. rufipectus different to previous mitochondrial data suggesting a hybrid origin or an ancient introgression of A. rufipectus; and 3) 22 CR1 insertion events strongly supported the eight phasianids under investigation sharing a common ancestor. Our study has revisited the phylogenetic position of G. gallus and A. rufipectus and provided a new insight into the phylogeny of Phasianidae birds. It showed that a CR1-based methodology has a great potential to be informative within Phasianidae in resolving relationships of closely related species whose radiation and speciation have occurred very recently.

Low dependency of retrotransposition on the ORF1 protein of the zebrafish LINE, ZfL2-1

The zebrafish long interspersed element (LINE), ZfL2-1, which belongs to the L2 clade, contains two open reading frames, ORF1 and ORF2. ORF1 encodes a protein containing a coiled-coil motif and an esterase domain, whereas ORF2 encodes a protein containing an endonuclease and a reverse transcriptase domain. To elucidate the functional significance of ORF1 in retrotransposition, we constructed many variants of ZfL2-1 and examined their retrotransposition ability. We concluded: 1) the ORF1 protein is not essential for ZfL2-1 retrotransposition in cultured cells; 2) the translation of ORF1 is required for the translation of ORF2; and 3) ORF2 translation probably occurs via suppression of the ORF1 stop codon, the efficiency of which is influenced by the context of the sequence juxtaposed to the 3' side of the stop codon. These results offer a new perspective on the evolution of the L2 clade LINEs.

A new mechanism to ensure integration during LINE retrotransposition: a suggestion from analyses of the 5' extra nucleotides

Long interspersed elements (LINEs) are transposable elements that exist in the chromosomal DNA of most eukaryotes; as such, they have a large impact on the genome evolution of their hosts. LINEs mobilize by a mechanism called retrotransposition in which the LINE RNA is reverse-transcribed into DNA and then integrated into the host chromosome. The integration of the 3' end of the LINE element simultaneously occurs with the initiation of reverse transcription; this process is called target-primed reverse transcription and is one of the important characteristics of LINEs. However, the molecular mechanism of the integration of the 5' end is not well understood. Here, we show that, in cultured cells, the integrants of the zebrafish ZfL2-2 LINE produce extra nucleotides at their 5' ends, and the extra nucleotides originate from their flanking sequences. We also found that, in cultured cells, some integrants of the human L1 LINE and, in their native hosts, some endogenous elements of two other LINEs also contain 5' extra nucleotides of similar origin, suggesting that the mechanism for generation of the 5' extra nucleotides is universal among various LINEs. From these data, we propose a general mechanism for 5' integration in LINE retrotransposition.

Self-interaction, nucleic acid binding, and nucleic acid chaperone activities are unexpectedly retained in the unique ORF1p of zebrafish LINE

Long interspersed elements (LINEs) are mobile elements that comprise a large proportion of many eukaryotic genomes. Although some LINE-encoded open reading frame 1 proteins (ORF1ps) were suggested to be required for LINE mobilization through binding to their RNA, their general role is not known. The ZfL2-1 ORF1p, which belongs to the esterase-type ORF1p, is especially interesting because it has no known RNA-binding domain. Here we demonstrate that ZfL2-1 ORF1p has all the canonical activities associated with known ORF1ps, including self-interaction, nucleic acid binding, and nucleic acid chaperone activities. In particular, we showed that its chaperone activity is reversible, suggesting that the chaperone activities of many other ORF1ps are also reversible. From this discovery, we propose that LINE ORF1ps play a general role in LINE integration by forming a complex with LINE RNA and rearranging its conformation.

Insertion of an intact CR1 retrotransposon in a cadherin gene linked with Bt resistance in the pink bollworm, Pectinophora gossypiella

Three mutations in the Pectinophora gossypiella cadherin gene PgCad1 are linked with resistance to Bacillus thuringiensis (Bt) toxin Cry1Ac. Here we show that the r3 mutation entails recent insertion into PgCad1 of an active chicken repeat (CR1) retrotransposon, designated CR1-1_Pg. Unlike most other CR1 elements, CR1-1_Pg is intact, transcribed by a flanking promoter, contains target site duplications and has a relatively low number of copies. Examination of transcripts from the PgCad1 locus revealed that CR1-1_Pg disrupts both the cadherin protein and a long noncoding RNA of unknown function. Together with previously reported data, these findings show that transposable elements disrupt eight of 12 cadherin alleles linked with resistance to Cry1Ac in three lepidopteran species, indicating that the cadherin locus is a common target for disruption by transposable elements.

Trimeric structure and flexibility of the L1ORF1 protein in human L1 retrotransposition

The LINE-1 (L1) retrotransposon emerges as a major source of human interindividual genetic variation, with important implications for evolution and disease. L1 retrotransposition is poorly understood at the molecular level, and the mechanistic details and evolutionary origin of the L1-encoded L1ORF1 protein (L1ORF1p) are particularly obscure. Here three crystal structures of trimeric L1ORF1p and NMR solution structures of individual domains reveal a sophisticated and highly structured, yet remarkably flexible, RNA-packaging protein. It trimerizes via an N-terminal, ion-containing coiled coil that serves as scaffold for the flexible attachment of the central RRM and the C-terminal CTD domains. The structures explain the specificity for single-stranded RNA substrates, and a mutational analysis indicates that the precise control of domain flexibility is critical for retrotransposition. Although the evolutionary origin of L1ORF1p remains unclear, our data reveal previously undetected structural and functional parallels to viral proteins.

Recent expansion of a new Ingi-related clade of Vingi non-LTR retrotransposons in hedgehogs

Autonomous non-long terminal repeat (non-LTR) retrotransposons and their repetitive remnants are ubiquitous components of mammalian genomes. Recently, we identified non-LTR retrotransposon families, Ingi-1_AAl and Ingi-1_EE, in two hedgehog genomes. Here we rename them to Vingi-1_AAl and Vingi-1_EE and report a new clade "Vingi," which is a sister clade of Ingi that lacks the ribonuclease H domain. In the European hedgehog genome, there are 11 non-autonomous families of elements derived from Vingi-1_EE by internal deletions. No retrotransposons related to Vingi elements were found in any of the remaining 33 mammalian genomes nearly completely sequenced to date, but we identified several new families of Vingi and Ingi retrotransposons outside mammals. Our data suggest the horizontal transfer of Vingi elements to hedgehog, although the vertical transfer cannot be ruled out. The compact structure and trans-mobilization of nonautonomous derivatives of Vingi can make them useful for in vivo retrotransposition assay system.

Endangered species hold clues to human evolution

We report that 18 conserved, and by extension functional, elements in the human genome are the result of retroposon insertions that are evolving under purifying selection in mammals. We show evidence that 1 of the 18 elements regulates the expression of ASXL3 during development by encoding an alternatively spliced exon that causes nonsense-mediated decay of the transcript. The retroposon that gave rise to these functional elements was quickly inactivated in the mammalian ancestor, and all traces of it have been lost due to neutral decay. However, the tuatara has maintained a near-ancestral version of this retroposon in its extant genome, which allows us to connect the 18 human elements to the evolutionary events that created them. We propose that conservation efforts over more than 100 years may not have only prevented the tuatara from going extinct but could have preserved our ability to understand the evolutionary history of functional elements in the human genome. Through simulations, we argue that species with historically low population sizes are more likely to harbor ancient mobile elements for long periods of time and in near-ancestral states, making these species indispensable in understanding the evolutionary origin of functional elements in the human genome.

Evolution of serum albumin intron-1 is shaped by a 5' truncated non-long terminal repeat retrotransposon in western Palearctic water frogs (Neobatrachia)

A 5' truncated non-LTR CR1-like retrotransposon, named RanaCR1, was identified in the serum albumin intron-1 (SAI-1) of at least seven species of western Palearctic water frogs (WPWF). Based on sequence similarity of the carboxy-terminal region (CTR) of ORF2 and/or the highly conserved 3' untranslated region (3' UTR), RanaCR1-like elements occur also in the genome of Xenopus tropicalis and Rana temporaria. Unlike other CR1 elements, RanaCR1 contains a CA microsatellite in its 3' UTR. The low nucleotide diversity of the 3' UTR compared to the CTR and to SAI-1 suggests that this region still plays a role in WPWF, either as a structure-stabilizing element, or within a species-specific transcriptional network. Length variation of water frog SAI-1 sequences is caused by deletions that extend in some cases beyond the 5' or 3' ends of RanaCR1, probably a result of selection for structural and functional stability of the primary transcript. The impact of RanaCR1 on SAI-1 evolution is also indicated by the significant negative correlation between the length of both SAI-1 and RanaCR1 and the percentage GC content of RanaCR1. Both SAI-1 and RanaCR1 sequences support the sister group relationship of R. perezi and R. saharica, which are placed in the phylogenetic tree at a basal position, the sister clade to other water frog taxa. It also supports the monophyly of the R. lessonae group; of Anatolian water frogs (R. cf. bedriagae), which are not conspecific with R. bedriagae, and of the European ridibunda group. Within the ridibunda clade, Greek frogs are clearly separated, supporting the hypothesis that Balkan water frogs represent a distinct species. Frogs from Atyrau (Kazakhstan), the type locality of R. ridibunda, were heterozygous for a ridibunda and a cf. bedriagae specific allele.

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

Rapidly growing number of sequenced genomes requires fast and accurate computational tools for analysis of different transposable elements (TEs). In this paper we focus on a rapid and reliable procedure for classification of autonomous non-LTR retrotransposons based on alignment and clustering of their reverse transcriptase (RT) domains. Typically, the RT domain protein sequences encoded by different non-LTR retrotransposons are similar to each other in terms of significant BLASTP E-values. Therefore, they can be easily detected by the routine BLASTP searches of genomic DNA sequences coding for proteins similar to the RT domains of known non-LTR retrotransposons. However, detailed classification of non-LTR retrotransposons, i.e. their assignment to specific clades, is a slow and complex procedure that is not formalized or integrated as a standard set of computational methods and data. Here we describe a tool (RTclass1) designed for the fast and accurate automated assignment of novel non-LTR retrotransposons to known or novel clades using phylogenetic analysis of the RT domain protein sequences. RTclass1 classifies a particular non-LTR retrotransposon based on its RT domain in less than 10 min on a standard desktop computer and achieves 99.5% accuracy. RT1class1 works either as a stand-alone program installed locally or as a web-server that can be accessed distantly by uploading sequence data through the internet (http://www.girinst.org/RTphylogeny/RTclass1).

Non-LTR retrotransposons encode noncanonical RRM domains in their first open reading frame

Non-LTR retrotransposons (NLRs) are a unique class of mobile genetic elements that have significant impact on the evolution of eukaryotic genomes. However, the molecular details and functions of their encoded proteins, in particular of the accessory ORF1p proteins, are poorly understood. Here, we identify noncanonical RNA-recognition-motifs (RRMs) in several phylogenetically unrelated NLR ORF1p proteins. This provides an explanation for their RNA-binding properties and clearly shows that they are not related to the retroviral nucleocapsid protein Gag, despite the frequent presence of CCHC zinc knuckles. In particular, we characterize the ORF1p protein of the human long interspersed nuclear element 1 (LINE-1 or L1). We show that L1ORF1p is a multidomain protein, consisting of a coiled coil (cc), RRM, and C-terminal domain (CTD). Most importantly, we solved the crystal structure of the RRM domain, which is characterized by extended loops stabilized by unique salt bridges. Furthermore, we demonstrate that L1ORF1p trimerizes via its N-terminal cc domain, and we suggest that this property is functionally important for all homologues. The formation of distinct complexes with single-stranded nucleic acids requires the presence of the RRM and CTD domains on the same polypeptide chain as well as their close cooperation. Finally, the phylogenetic analysis of mammalian L1ORF1p shows an ancient origin of the RRM domain and supports a modular evolution of NLRs.

Mobility pathways for vertebrate L1, L2, CR1, and RTE clade retrotransposons

Autonomous non-long terminal repeat retrotransposons (NLRs) are ubiquitous mobile genetic elements that insert their DNA copies at new locations by retrotransposition. In vertebrates, there are 4 NLR clades, L1, L2, CR1, and RTE, which diverged in the Precambrian era. It has been demonstrated that retrotransposition of L1 and L2 members proceeds via coordinated reactions of targeted DNA cleavage and reverse transcription catalyzed by the NLR-encoded proteins, which are followed by the joining of the 5' (upstream) junction. However, the study on the mobility pathways for vertebrate NLRs is so far limited to L1 and L2. In this report, using target analysis of nested transposons for genomic copies, we studied retrotransposition pathways for a variety of vertebrate NLRs, including those of the L1, L2, CR1, and RTE clades in the human, cow, opossum, chicken, and zebrafish genomes. Thus, this study constitutes the first comprehensive analysis of NLR retrotransposition products in vertebrates. Our data revealed that these elements share similar mechanisms for the cleavages of the 2 target DNA strands and for the initiation of reverse transcription. Possible endonuclease-independent insertions were also identified. Overall, our results suggest the existence of multiple retrotransposition pathways that are conserved among the diverse NLR clades in various vertebrate hosts.

Waves of genomic hitchhikers shed light on the evolution of gamebirds (Aves: Galliformes)

BACKGROUND

The phylogenetic tree of Galliformes (gamebirds, including megapodes, currassows, guinea fowl, New and Old World quails, chicken, pheasants, grouse, and turkeys) has been considerably remodeled over the last decades as new data and analytical methods became available. Analyzing presence/absence patterns of retroposed elements avoids the problems of homoplastic characters inherent in other methodologies. In gamebirds, chicken repeats 1 (CR1) are the most prevalent retroposed elements, but little is known about the activity of their various subtypes over time. Ascertaining the fixation patterns of CR1 elements would help unravel the phylogeny of gamebirds and other poorly resolved avian clades.

RESULTS

We analyzed 1,978 nested CR1 elements and developed a multidimensional approach taking advantage of their transposition in transposition character (TinT) to characterize the fixation patterns of all 22 known chicken CR1 subtypes. The presence/absence patterns of those elements that were active at different periods of gamebird evolution provided evidence for a clade (Cracidae + (Numididae + (Odontophoridae + Phasianidae))) not including Megapodiidae; and for Rollulus as the sister taxon of the other analyzed Phasianidae. Genomic trace sequences of the turkey genome further demonstrated that the endangered African Congo Peafowl (Afropavo congensis) is the sister taxon of the Asian Peafowl (Pavo), rejecting other predominantly morphology-based groupings, and that phasianids are monophyletic, including the sister taxa Tetraoninae and Meleagridinae.

CONCLUSION

The TinT information concerning relative fixation times of CR1 subtypes enabled us to efficiently investigate gamebird phylogeny and to reconstruct an unambiguous tree topology. This method should provide a useful tool for investigations in other taxonomic groups as well.

The Juan non-LTR retrotransposon in mosquitoes: genomic impact, vertical transmission and indications of recent and widespread activity

BACKGROUND

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

RESULTS

Systematic analysis of the ~1.3 Gbp Aedes aegypti (Ae. aegypti) genome sequence suggests that Juan-A is the only Juan-type non-LTR in Aedes aegypti. Juan-A is highly reiterated and comprises approximately 3% of the genome. Using minimum cutoffs of 90% length and 70% nucleotide (nt) identity, 663 copies were found by BLAST using the published Juan-A sequence as the query. All 663 copies are at least 95% identical to Juan-A, while 378 of these copies are 99% identical to Juan-A, indicating that the Juan-A family has been transposing recently in evolutionary history. Using the 0.34 Kb 5' UTR as the query, over 2000 copies were identified that may contain internal promoters, leading to questions on the genomic impact of Juan-A. Juan sequences were obtained by PCR, library screening, and database searches for 18 mosquito species of six genera including Aedes, Ochlerotatus, Psorophora, Culex, Deinocerites, and Wyeomyia. Comparison of host and Juan phylogenies shows overall congruence with few exceptions.

CONCLUSION

Juan-A is a major genomic component in Ae. aegypti and it has been retrotransposing recently in evolutionary history. There are also indications that Juan has been recently active in a wide range of mosquito species. Furthermore, our research demonstrates that a Jockey clade non-LTR without target site-specificity has been sustained by vertical transmission in the mosquito family. These results strengthen the argument that non-LTRs tend to be genomic elements capable of persistence by vertical descent over a long evolutionary time.

New Trypanosoma cruzi repeated element that shows site specificity for insertion

A new family of site-specific repeated elements identified in Trypanosoma cruzi, which we named TcTREZO, is described here. TcTREZO appears to be a composite repeated element, since three subregions may be defined within it on the basis of sequence similarities with other T. cruzi sequences. Analysis of the distribution of TcTREZO in the genome clearly indicates that it displays site specificity for insertion. Most TcTREZO elements are flanked by conserved sequences. There is a highly conserved 68-bp sequence at the 5' end of the element and a sequence domain of approximately 500 bp without a well-defined borderline at the 3' end. Northern blot hybridization and reverse transcriptase PCR analyses showed that TcTREZO transcripts are expressed as oligo(A)-terminated transcripts whose length corresponds to the unit size of the element (1.6 kb). Transcripts of approximately 0.2 kb derived from a small part of TcTREZO are also detected in steady-state RNA. TcTREZO transcripts are unspliced and not translated. The copy number of TcTREZO sequences was estimated to be approximately 173 copies per haploid genome. TcTREZO appears to have been assembled by insertions of sequences into a progenitor element. Once associated with each other, these subunits were amplified as a new transposable element. TcTREZO shows site specificity for insertion, suggesting that a sequence-specific endonuclease could be responsible for its insertion at a unique site.

A new system for analyzing LINE retrotransposition in the chicken DT40 cell line widely used for reverse genetics

Long interspersed elements (LINEs) are autonomous transposable elements that proliferate via retrotransposition, which involves reverse transcription of LINE RNAs. It is anticipated that LINE retrotransposition requires both LINE-encoded proteins and host-encoded proteins. However, identification of the host factors, their roles, and the steps at which they act on retrotransposition are poorly understood because of the lack of an appropriate genetic system to study LINE retrotransposition in a series of mutant hosts. To construct such a genetic system, we applied the retrotransposition-indicative cassette method to DT40 cells, a chicken cell line for which a variety of isogenic mutants have been established by gene targeting. Because DT40 cells are non-adherent, we utilized a selective soft agarose medium to allow the formation of colonies of cells that had undergone LINE retrotransposition. Colony formation was completely dependent on the activities of the LINE-encoded proteins and on the presence of the essential 3' region of the LINE RNA. Moreover, the selected colonies indeed carried retrotransposed LINE copies in their chromosomes, with integration features similar to those of genomic (native) LINE copies. This method thus allows the authentic selection of LINE-retrotransposed cells and the approximate recapitulation of retrotransposition events that occur in nature. Therefore, the DT40 cell system established here provides a powerful tool for the elucidation of LINE retrotransposition pathways, the host factors involved, and their roles.

Novel retrotransposon analysis reveals multiple mobility pathways dictated by hosts

Autonomous non-long-terminal-repeat retrotransposons (NLRs) proliferate by retrotransposition via coordinated reactions of target DNA cleavage and reverse transcription by a mechanism called target-primed reverse transcription (TPRT). Whereas this mechanism guarantees the covalent attachment of the NLR and its target site at the 3' junction, mechanisms for the joining at the 5' junction have been conjectural. To better understand the retrotransposition pathways, we analyzed target-NLR junctions of zebrafish NLRs with a new method of identifying genomic copies that reside within other transposons, termed "target analysis of nested transposons" (TANT). Application of the TANT method revealed various features of the zebrafish NLR integrants; for example, half of the integrants carry extra nucleotides at the 5' junction, which is in stark contrast to the major human NLR, LINE-1. Interestingly, in a cell culture assay, retrotransposition of the zebrafish NLR in heterologous human cells did not bear extra 5' nucleotides, indicating that the choice of the 5' joining pathway is affected by the host. Our results suggest that several pathways exist for NLR retrotransposition and argue in favor of host protein involvement. With genomic sequence information accumulating exponentially, our data demonstrate the general applicability of the TANT method for the analysis of a wide variety of retrotransposons.

Two families of non-LTR retrotransposons, Syrinx and Daphne, from the Darwinulid ostracod, Darwinula stevensoni

Two novel families of non-LTR retrotransposons, named Syrinx and Daphne, were cloned and characterized in a putative ancient asexual ostracod Darwinula stevensoni. Phylogenetic analysis reveals that Daphne is the founding member of a novel clade of non-LTR retroelements, which also contains retrotransposon families from the sea urchin and the silkworm and forms a sister clade to L2-like elements. The Syrinx family of non-LTR retrotransposons exhibits evidence of relatively recent activity, manifested in high levels of sequence similarity between individual copies and a three- to ten-fold excess of synonymous substitutions, which is indicative of purifying selection. The Daphne family may have very few copies with intact open reading frames, and exhibits neutral within-family ratio of non-synonymous to synonymous substitutions. It can additionally be characterized by formation of inverted truncated head-to-head structures. All of these features make recent activity less likely than in the Syrinx family. Our results are discussed in light of the evolutionary consequences of long-term asexuality in general and in D. stevensoni in particular.

REM1, a new type of long terminal repeat retrotransposon in Chlamydomonas reinhardtii

A new long terminal repeat (LTR) retrotransposon, named REM1, has been identified in the green alga Chlamydomonas reinhardtii. It was found in low copy number, highly methylated, and with an inducible transpositional activity. This retrotransposon is phylogenetically related to Ty3-gypsy LTR retrotransposons and possesses new and unusual structural features. A regulatory module, ORF3p, is present in an inverse transcriptional orientation to that of the polyprotein and contains PHD-finger and chromodomains, which might confer specificity of the target site and are highly conserved in proteins involved in transcriptional regulation by chromatin remodeling. By using different wild-type and mutant strains, we show that CrREM1 was active with a strong transcriptional activity and amplified its copy number in strains that underwent foreign DNA integration and/or genetic crosses. However, integration of CrREM1 was restricted to these events even though the expression of its full-length transcripts remained highly activated. A regulatory mechanism of CrREM1 retrotransposition which would help to minimize its deleterious effects in the host genome is proposed.

Isolation and characterization of retrotransposition-competent LINEs from zebrafish

Long interspersed elements (LINEs) are a type of retroposon and are widely distributed in most eukaryotic genomes. LINEs are classified into two groups, the stringent type and relaxed type, based on the recognition of the 3' tail of their own RNA by reverse transcriptase (RT) during retrotransposition. Although most LINEs are thought to belong to the stringent type, retrotransposition studies of the stringent type LINEs are relatively limited compared with those of the relaxed type. We have now isolated two retrotransposition-competent LINEs (ZfL2-1 and ZfL2-2) from the zebrafish genome. Both ZfL2-1 and ZfL2-2 are members of the L2 clade; ZfL2-1 encodes two open reading frames (ORFs) and ZfL2-2 encodes one ORF, and each of the ORFs is required for retrotransposition. Using a retrotransposition assay in HeLa cells, we established that both ZfL2-1 and Zfl2-2 belong to the stringent type. We also demonstrated that an esterase (ES) domain encoded by ZfL2-1 ORF1 strongly enhances its own retrotransposition. The ES domain is encoded only in ORF1 of LINEs classified in the CR1 and L2 clades, although its function or significance in retrotransposition has not been elucidated. Thus, this is the first experimental evidence that the ES domain has an enhancing function during retrotransposition. These zebrafish LINEs will be useful for determining the function of ORF1 and the retrotransposition mechanism of stringent-type LINEs.

Repbase Update, a database of eukaryotic repetitive elements

Repbase Update is a comprehensive database of repetitive elements from diverse eukaryotic organisms. Currently, it contains over 3600 annotated sequences representing different families and subfamilies of repeats, many of which are unreported anywhere else. Each sequence is accompanied by a short description and references to the original contributors. Repbase Update includes Repbase Reports, an electronic journal publishing newly discovered transposable elements, and the Transposon Pub, a web-based browser of selected chromosomal maps of transposable elements. Sequences from Repbase Update are used to screen and annotate repetitive elements using programs such as Censor and RepeatMasker. Repbase Update is available on the worldwide web at http://www.girinst.org/Repbase_Update.html.

SINEs and LINEs: symbionts of eukaryotic genomes with a common tail

Many SINEs and LINEs have been characterized to date, and examples of the SINE and LINE pair that have the same 3' end sequence have also increased. We report the phylogenetic relationships of nearly all known LINEs from which SINEs are derived, including a new example of a SINE/LINE pair identified in the salmon genome. We also use several biological examples to discuss the impact and significance of SINEs and LINEs in the evolution of vertebrate genomes.