APE-type non-LTR retrotransposons: determinants involved in target site recognition

Structural dissection of sequence recognition and catalytic mechanism of human LINE-1 endonuclease

Long interspersed nuclear element-1 (L1) is an autonomous non-LTR retrotransposon comprising ∼20% of the human genome. L1 self-propagation causes genomic instability and is strongly associated with aging, cancer and other diseases. The endonuclease domain of L1’s ORFp2 protein (L1-EN) initiates de novo L1 integration by nicking the consensus sequence 5′-TTTTT/AA-3′. In contrast, related nucleases including structurally conserved apurinic/apyrimidinic endonuclease 1 (APE1) are non-sequence specific. To investigate mechanisms underlying sequence recognition and catalysis by L1-EN, we solved crystal structures of L1-EN complexed with DNA substrates. This showed that conformational properties of the preferred sequence drive L1-EN’s sequence-specificity and catalysis. Unlike APE1, L1-EN does not bend the DNA helix, but rather causes ‘compression’ near the cleavage site. This provides multiple advantages for L1-EN’s role in retrotransposition including facilitating use of the nicked poly-T DNA strand as a primer for reverse transcription. We also observed two alternative conformations of the scissile bond phosphate, which allowed us to model distinct conformations for a nucleophilic attack and a transition state that are likely applicable to the entire family of nucleases. This work adds to our mechanistic understanding of L1-EN and related nucleases and should facilitate development of L1-EN inhibitors as potential anticancer and antiaging therapeutics.

Chromosomal Locations of a Non-LTR Retrotransposon, , in and , and Its Implication on Genome Evolution of Species

Mobile elements are major regulators of genome evolution through their effects on genome size and chromosome structure in higher organisms. Non-long terminal repeat (non-LTR) retrotransposons, one of the subclasses of transposons, are specifically inserted into repetitive DNA sequences. While studies on the insertion of non-LTR retrotransposons into ribosomal RNA genes and other repetitive DNA sequences have been reported in the animal kingdom, studies in the plant kingdom are limited. Here, using FISH, we confirmed that <i>Menolird18</i>, a member of LINE (long interspersed nuclear element) in non-LTR retrotransposons and found in <i>Cucumis melo</i>, was inserted into ITS and ETS (internal and external transcribed spacers) regions of 18S rDNA in melon and cucumber. Beside the 18S rDNA regions, <i>Menolird18</i> was also detected in all centromeric regions of melon, while it was located at pericentromeric and sub-telomeric regions in cucumber. The fact that FISH signals of <i>Menolird18</i> were found in centromeric and rDNA regions of mitotic chromosomes suggests that <i>Menolird18</i> is a rDNA and centromere-specific non-LTR retrotransposon in melon. Our findings are the first report on a non-LTR retrotransposon that is highly conserved in 2 different plant species, melon and cucumber. The clear distinction of chromosomal localization of <i>Menolird18</i> in melon and cucumber implies that it might have been involved in the evolutionary processes of the melon (<i>C. melo</i>) and cucumber (<i>C. sativus</i>) genomes.

Completion of LINE integration involves an open '4-way' branched DNA intermediate

Long Interspersed Elements (LINEs), also known as non-LTR retrotransposons, encode a multifunctional protein that reverse transcribes its mRNA into DNA at the site of insertion by target primed reverse transcription. The second half of the integration reaction remains very poorly understood. Second-strand DNA cleavage and second-strand DNA synthesis were investigated in vitro using purified components from a site-specific restriction-like endonuclease (RLE) bearing LINE. DNA structure was shown to be a critical component of second-strand DNA cleavage. A hitherto unknown and unexplored integration intermediate, an open ‘4-way’ DNA junction, was recognized by the element protein and cleaved in a Holliday junction resolvase-like reaction. Cleavage of the 4-way junction resulted in a natural primer-template pairing used for second-strand DNA synthesis. A new model for RLE LINE integration is presented.

Cross-Kingdom Commonality of a Novel Insertion Signature of RTE-Related Short Retroposons

In multicellular organisms, such as vertebrates and flowering plants, horizontal transfer (HT) of genetic information is thought to be a rare event. However, recent findings unveiled unexpectedly frequent HT of RTE-clade LINEs. To elucidate the molecular footprints of the genomic integration machinery of RTE-related retroposons, the sequence patterns surrounding the insertion sites of plant Au-like SINE families were analyzed in the genomes of a wide variety of flowering plants. A novel and remarkable finding regarding target site duplications (TSDs) for SINEs was they start with thymine approximately one helical pitch (ten nucleotides) downstream of a thymine stretch. This TSD pattern was found in RTE-clade LINEs, which share the 3'-end sequence of these SINEs, in the genome of leguminous plants. These results demonstrably show that Au-like SINEs were mobilized by the enzymatic machinery of RTE-clade LINEs. Further, we discovered the same TSD pattern in animal SINEs from lizard and mammals, in which the RTE-clade LINEs sharing the 3'-end sequence with these animal SINEs showed a distinct TSD pattern. Moreover, a significant correlation was observed between the first nucleotide of TSDs and microsatellite-like sequences found at the 3'-ends of SINEs and LINEs. We propose that RTE-encoded protein could preferentially bind to a DNA region that contains a thymine stretch to cleave a phosphodiester bond downstream of the stretch. Further, determination of cleavage sites and/or efficiency of primer sites for reverse transcription may depend on microsatellite-like repeats in the RNA template. Such a unique mechanism may have enabled retroposons to successfully expand in frontier genomes after HT.

Globular domain structure and function of restriction-like-endonuclease LINEs: similarities to eukaryotic splicing factor Prp8

Background

R2 elements are a clade of early branching Long Interspersed Elements (LINEs). LINEs are retrotransposable elements whose replication can have profound effects on the genomes in which they reside. No crystal or EM structures exist for the reverse transcriptase (RT) and linker regions of LINEs.

Results

Using limited proteolysis as a probe for globular domain structure, we show that the protein encoded by the Bombyx mori R2 element has two major globular domains: (1) a small globular domain consisting of the N-terminal zinc finger and Myb motifs, and (2) a large globular domain consisting of the RT, linker, and type II restriction-like endonuclease (RLE). Further digestion of the large globular domain occurred within the RT. Mapping these RT cleavages onto an updated model of the R2Bm RT indicated that the thumb of the RT was largely protected from proteolytic cleavage. The crystal structure of the large globular domain of Prp8, a eukaryotic splicing factor, was a major template used in building the R2Bm RT model, particularly the thumb region. The large fragment of Prp8 consists not only of a RT similar to R2Bm, but also an RLE and a linker connecting the two regions. The linker sequences adjacent to the RLE in LINEs and Prp8 share a set of two important α-helices and a (presumptive) knuckle/ββα structural motif that are closely associated with the thumb. The RLEs of LINEs and Prp8 share a unique catalytic core residue spacing as well as other key residues.

Conclusions

The protein encoded by RLE LINEs consists of two major globular domains. The larger of the two globular domain contains the RT, linker, and RLE and is similar to the large fragment of the spliceosomal protein Prp8. The similarities are suggestive of possible common ancestry.

Electronic supplementary material

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

Evaluating different DNA binding domains to modulate L1 ORF2p-driven site-specific retrotransposition events in human cells

DNA binding domains (DBDs) have been used with great success to impart targeting capabilities to a variety of proteins creating highly useful genomic tools. We evaluated the ability of five types of DBDs and strategies (AAV Rep proteins, Cre, TAL effectors, zinc finger proteins, and Cas9/gRNA system) to target the L1 ORF2 protein to drive retrotransposition of Alu inserts to specific sequences in the human genome. First, we find that the L1 ORF2 protein tolerates the addition of protein domains both at the amino- and carboxy-terminus. Although in some instances retrotransposition efficiencies slightly diminished, all fusion proteins containing an intact ORF2 were capable of driving retrotransposition. Second, the stability of individual ORF2 fusion proteins varies and difficult to predict. Third, DBDs that require the formation of multimers for target recognition are unlikely to modify targeting of ORF2p-driven insertions. Fourth, the more components needed to assemble into a complex to drive targeted retrotransposition, the less likely the strategy will increase targeted insertions. Fifth, abundance of target sequences present in the genome will likely dictate the effectiveness and efficiency of targeted insertions. Lastly, the cleavage capabilities of Cas9 (or a Cas9 nickase variant) are unable to substitute for the L1 ORF2 endonuclease domain functions, suggestive that the endonuclease domain has alternate functions needed for retrotransposition. From these studies, we conclude that the most critical component for the modification of the human L1 ORF2 protein to drive targeted insertions is the selection of the DBD due to the varying functional requirements and impacts on protein stability.

Site-specific non-LTR retrotransposons

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

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.

Structural features and mechanism of translocation of non-LTR retrotransposons in Candida albicans

A number of abundant mobile genetic elements called retrotransposons reverse transcribe RNA to generate DNA for insertion into eukaryotic genomes. Non-long-terminal repeat (non-LTR) retrotransposons represent a major class of retrotransposons, and transposons that move by target-primed reverse transcription lack LTRs characteristic of retroviruses and retroviral-like transposons. Yeast model systems in Candida albicans and Saccharomyces cerevisiae have been developed for the study of non-LTR retrotransposons. Non-LTR retrotransposons are divided into LINEs (long interspersed nuclear elements), SINEs (short interspersed nuclear elements), and SVA (SINE, VNTR, and Alu). LINE-1 elements have been described in fungi, and several families called Zorro elements have been detected from C. albicans. They are all members of L1 clades. Through a mechanism named target-primed reverse transcription (TPRT), LINEs translocate the new copy into the target site to initiate DNA synthesis primed by the 3′ OH of the broken strand. In this article, we describe some advances in the research on structural features and origin of non-LTR retrotransposons in C. albicans, and discuss mechanisms underlying their reverse transcription and integration of the donor copy into the target site.

Non-LTR retrotransposons and microsatellites: Partners in genomic variation

The human genome is laden with both non-LTR (long-terminal repeat) retrotransposons and microsatellite repeats. Both types of sequences are able to, either actively or passively, mutagenize the genomes of human individuals and are therefore poised to dynamically alter the human genomic landscape across generations. Non-LTR retrotransposons, such as L1 and Alu, are a major source of new microsatellites, which are born both concurrently and subsequently to L1 and Alu integration into the genome. Likewise, the mutation dynamics of microsatellite repeats have a direct impact on the fitness of their non-LTR retrotransposon parent owing to microsatellite expansion and contraction. This review explores the interactions and dynamics between non-LTR retrotransposons and microsatellites in the context of genomic variation and evolution.

Expression of a LINE-1 endonuclease variant in gastric cancer: its association with clinicopathological parameters

Background

Long interspersed nuclear element-1 (LINE-1 or L1), the most abundant and only autonomously active family of non-LTR retrotransposons in the human genome, expressed not only in the germ lines but also in somatic tissues. It contributes to genetic instability, aging, and age-related diseases, such as cancer. Our previous study identified in human gastric adenocarcinoma an upregulated transcript GCRG213, which shared 88% homology with human L1 sequence and contained a putative conserved apurinic/apyrimidinic endonucleas1 domain.

Methods

Immunohistochemistry was carried out by using a monoclonal mouse anti-human GCRG213 protein (GCRG213p) antibody produced in our laboratory, on tissue microarray constructed with specimens from 175 gastric adenocarcinoma patients. The correlation between GCRG213p expression and patient clinicopathological parameters was evaluated. GCRG213p expression in gastric cancer cell lines were studied using Western blotting analysis. L1 promoter methylation status of gastric cancer cells was tested using methylation-specific PCR. BLASTP was used at the NCBI Blast server to identify GCRG213p sequence to any alignments in the Protein Data Bank databases.

Results

Most primary gastric cancer, lymph node metastases and gastric intestinal metaplasia glands showed positive GCRG213p immunoreactivity. High GCRG213p immunostaining score in the primary gastric cancer was positively correlated with tumor differentiation (well differentiated, p = 0.001), Lauren’s classification (intestinal type, p < 0.05) and a late age onset of gastric adenocarcinoma (≥65 yrs; p < 0.05). GCRG213p expression has no association with other clinicopathological parameters, including survival. Western blotting analysis of GCRG213p expression in gastric cancer cells indicated that GCRG213p level was higher in gastric cancer cell lines than in human normal gastric epithelium immortalized cell line GES-1. Partial methylation of L1 in gastric cancer cells was confirmed by methylation-specific PCR. BLASTP program analysis revealed that GCRG213p peptide shared 83.0% alignment with the C-terminal region of L1 endonuclease (L1-EN). GCRG213p sequence possesses the important residues that compose the conserved features of L1-EN.

Conclusions

GCRG213p could be a variant of L1-EN, a functional member of L1-EN family. Overexpression of GCRG213p is common in both primary gastric cancer and lymph node metastasis. These findings provide evidence of somatic L1 expression in gastric cancer, and its potential consequences in the form of tumor.

Endonuclease domain of the Drosophila melanogaster R2 non-LTR retrotransposon and related retroelements: a new model for transposition

The molecular mechanisms of the transposition of non-long terminal repeat (non-LTR) retrotransposons are not well understood; the key questions of how the 3′-ends of cDNA copies integrate and how site-specific integration occurs remain unresolved. Integration depends on properties of the endonuclease (EN) domain of retrotransposons. Using the EN domain of the Drosophila R2 retrotransposon as a model for other, closely related non-LTR retrotransposons, we investigated the EN domain and found that it resembles archaeal Holliday-junction resolvases. We suggest that these non-LTR retrotransposons are co-transcribed with the host transcript. Combined with the proposed resolvase activity of the EN domain, this model yields a novel mechanism for site-specific retrotransposition within this class of retrotransposons, with resolution proceeding via a Holliday junction intermediate.

Independently derived targeting of 28S rDNA by A- and D-clade R2 retrotransposons: Plasticity of integration mechanism

Restriction-like endonuclease (RLE) bearing non-LTR retrotransposons are site-specific elements that integrate into the genome through a target primed reverse transcription mechanism (TPRT). R2 elements have been used as a model system for investigating non-LTR retrotransposon integration. We previously demonstrated that R2 retrotransposons require two subunits of the element-encoded multifunctional protein to integrate-one subunit bound upstream of the insertion site and one bound downstream. R2 elements have been phylogenetically categorized into four clades: R2-A, B, C and D, that diverged from a common ancestor more than 850 million years ago. All R2 elements target the same sequence within 28S rDNA. The amino-terminal domain of R2Bm, an R2-D clade element, contains a single zinc finger and a Myb motif that are responsible for binding R2 protein downstream of the insertion site. Target site recognition is of interest as it is the first step in the integration reaction and may help elucidate evolutionary history and integration mechanism. The amino-terminal domain of R2-A clade members contains three zinc fingers and a Myb motif. We show here that R2Lp, an R2-A clade member, uses its amino-terminal DNA binding motifs to bind upstream of the insertion site. Because the R2-A and R2-D clade elements recognize 28S rDNA differently, we conclude the A- and D-clades represent independent targeting events to the 28S site. Our results also indicate a certain plasticity of insertional mechanics exists between the two clades.

Rapid decay of unstable Leishmania mRNAs bearing a conserved retroposon signature 3'-UTR motif is initiated by a site-specific endonucleolytic cleavage without prior deadenylation

We have previously shown that the Leishmania genome possess two widespread families of extinct retroposons termed Short Interspersed DEgenerated Retroposons (SIDER1/2) that play a role in post-transcriptional regulation. Moreover, we have demonstrated that SIDER2 retroposons promote mRNA degradation. Here we provide new insights into the mechanism by which unstable Leishmania mRNAs harboring a SIDER2 retroposon in their 3′-untranslated region are degraded. We show that, unlike most eukaryotic transcripts, SIDER2-bearing mRNAs do not undergo poly(A) tail shortening prior to rapid turnover, but instead, they are targeted for degradation by a site-specific endonucleolytic cleavage. The main cleavage site was mapped in two randomly selected SIDER2-containing mRNAs in vivo between an AU dinucleotide at the 5′-end of the second 79-nt signature (signature II), which represents the most conserved sequence amongst SIDER2 retroposons. Deletion of signature II abolished endonucleolytic cleavage and deadenylation-independent decay and increased mRNA stability. Interestingly, we show that overexpression of SIDER2 anti-sense RNA can increase sense transcript abundance and stability, and that complementarity to the cleavage region is required for protecting SIDER2-containing transcripts from degradation. These results establish a new paradigm for how unstable mRNAs are degraded in Leishmania and could serve as the basis for a better understanding of mRNA decay pathways in general.

Insights into the DNA cleavage mechanism of human LINE-1 retrotransposon endonuclease

The human LINE-1 endonuclease (L1-EN) contributes in defining the genomic integration sites of the abundant human L1 and Alu retrotransposons. LINEs have been considered as possible vehicles for gene delivery and understanding the mechanism of L1-EN could help engineering them as genetic tools. We tested the in vitro activity of point mutants in three L1-EN residues--Asp145, Arg155, Ile204--that are key for DNA cleavage, and determined their crystal structures. The L1-EN structure remains overall unaffected by the mutations, which change the enzyme activity but leave DNA cleavage sequence specificity mostly unaffected. To better understand the mechanism of L1-EN, we performed molecular dynamics simulations using as model the structures of wild type EN-L1, of two betaB6-betaB5 loop exchange mutants we have described previously to be important for DNA recognition, of the R155A mutant from this study, and of the homologous TRAS1 endonuclease: all confirm a rigid scaffold. The simulations crucially indicate that the betaB6-betaB5 loop shows an anticorrelated motion with the surface loops betaA6-betaA5 and betaB3-alphaB1. The latter loop harbors N118, a residue that alters DNA cleavage specificity in homologous endonucleases, and implies that the plasticity and correlated motion of these loops has a functional importance in DNA recognition and binding. To further explore how these loops are possibly involved in DNA binding, we docked computationally two DNA substrates to our structure, one involving a flipped-out nucleotide downstream the scissile phosphodiester; and one not. The models for both scenarios are feasible and agree with the hypotheses derived from the dynamic simulations. The reduced cleavage activity we have observed for the I204Y mutant above however, favors the flipped out nucleotide model.

Multiple duplication of the bucentaur gene family, which recruits the APE-like domain of retrotransposon: Identification of a novel homolog and distinct cellular expression

The p97bcnt/cfdp2 is a ruminant-specific gene created by a combination of gene duplication of ancestral bcnt (bucentaur) or cfdp1 (craniofacial developmental protein 1), bcnt/cfdp1, and the insertion of a retrotransposable element-1 (RTE). As a result, p97Bcnt recruits the whole apurinic/apyrimidinic endonuclease (APE)-like domain of RTE in the middle of the molecule (RTE domain) as a region encoded by an exon. In addition, p97Bcnt contains two intramolecular repeats (IRs) of 40 amino acids each in the C-terminal region, whereas Bcnt/Cfdp1 contains one IR. We have identified an additional bovine homolog with a structure highly similar to p97Bcnt, designated p97Bcnt2, which contains three IRs. p97bcnt2 is located in tandem with bcnt/cfdp1 and p97bcnt within a 177-kb range on bovine chromosome 18, a syntenic region of human chromosome 16. The gene product is expressed as a protein with an apparent molecular mass of 102 kDa. The phylogenetic tree strongly suggests that p97bcnt-2 forms a third clade of the bcnt family and that the first duplication of the IR unit occurred prior to the divergence of p97bcnt and p97bcnt-2. To address the question of whether these bcnt members have distinct functions, we first examined the expression and localization of the p97Bcnt family members. p97Bcnt is substantially expressed in many tissues involved in responses to external and internal stress. In the testis, p97Bcnt localizes preferentially in the nuclei of spermatozoa, while Bcnt/Cfdp1 localizes predominantly in the cytosol of Leydig cells and some spermatogenic cells, implying that at least these two molecules of the Bcnt family play different functional roles. These results provide evidence for the direct contribution of RTE to gene diversity to form gene families that may support cellular function.

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.

Target site analysis of RTE1_LA and its AfroSINE partner in the elephant genome

SINEs retrotranspose using their partner LINE's enzymatic machinery. It has recently been proposed that AfroSINEs ending with GGTTT 3' tandem repeats were mobilized by RTE elements ending with CAA 3' tandem repeats in the Afrotherian genome. Using sequences from the elephant genome, we show that AfroSINEs derive from RTE ending with GGTTT-like 3' tandem repeats, a subgroup of RTE1_LA that only reached low copy number, and confirm that they were most likely mobilized by RTE ending with CAA(n) tandem repeats (RTE1_LA-CAA(n)). This partnership is supported by sequence similarity between two regions of the elements, overlap in the timing of their activity, common features of their target site consensus that are not shared by other members of the RTE family, and their high copy number. Detailed analyses of pre-insertion loci reveal that like many other apurinic/apyrimidinic endonuclease encoding elements, RTE1_LA-CAA(n) shows loose target site specificity. In addition, the RTE1_LA-CAA(n) target site consensus shares several structural and primary sequence features with that of LINE1, suggesting that these two elements share close functional similarity in the target primed reverse transcription (TPRT) reaction. Interestingly, although globally similar, the target site consensus of AfroSINE(Anc) and RTE1_LA-CAA(n) differ in several aspects. These differences, not observed among all SINE/LINE pairs so far examined, are most likely due to the fact that AfroSINEs and RTE1_LA-CAA(n) are terminated by a different tandem repeat motif. We propose that these differences reflect constraints imposed by base pairing interactions between the mRNA 3' terminal tandem repeats and the target DNA at the onset of TPRT. So in addition to the endonuclease nicking preference, the mRNA of these elements appears to play an important role in integration site choice through a passive, post-nicking, selective process.

Determinants for DNA target structure selectivity of the human LINE-1 retrotransposon endonuclease

The human LINE-1 endonuclease (L1-EN) is the targeting endonuclease encoded by the human LINE-1 (L1) retrotransposon. L1-EN guides the genomic integration of new L1 and Alu elements that presently account for ∼28% of the human genome. L1-EN bears considerable technological interest, because its target selectivity may ultimately be engineered to allow the site-specific integration of DNA into defined genomic locations. Based on the crystal structure, we generated L1-EN mutants to analyze and manipulate DNA target site recognition. Crystal structures and their dynamic and functional analysis show entire loop grafts to be feasible, resulting in altered specificity, while individual point mutations do not change the nicking pattern of L1-EN. Structural parameters of the DNA target seem more important for recognition than the nucleotide sequence, and nicking profiles on DNA oligonucleotides in vitro are less well defined than the respective integration site consensus in vivo. This suggests that additional factors other than the DNA nicking specificity of L1-EN contribute to the targeted integration of non-LTR retrotransposons.

Getting closer to a pre-vertebrate genome: the non-LTR retrotransposons of Branchiostoma floridae

Non-LTR retrotransposons are common in vertebrate genomes and although present in invertebrates they appear at a much lower frequency. The cephalochordate amphioxus is the closest living relative to vertebrates and has been considered a good model for comparative analyses of genome expansions during vertebrate evolution. With the aim to assess the involvement of transposable elements in these events, we have analysed the non-LTR retrotransposons of Branchiostoma floridae. In silico searches have allowed to reconstruct non-LTR elements of six different clades (CR1, I, L1, L2, NeSL and RTE) and assess their structural features. According to the estimated copy number of these elements they account for less than 1% of the haploid genome, which reminds of the low abundance also encountered in the urochordate Ciona intestinalis. Amphioxus (B. floridae) and Ciona share a pre-vertebrate-like organization for the non-LTR retrotransposons (<150 copies, < 1% of the genome) versus the complexity associated to higher vertebrates (Homo sapiens >1.3·10⁶ copies, > 20% of the genome).