Skip to main content
  • Discovery notes
  • Open access
  • Published:

Transduplication resulted in the incorporation of two protein-coding sequences into the Turmoil-1 transposable element of C. elegans

Abstract

Transposable elements may acquire unrelated gene fragments into their sequences in a process called transduplication. Transduplication of protein-coding genes is common in plants, but is unknown of in animals. Here, we report that the Turmoil-1 transposable element in C. elegans has incorporated two protein-coding sequences into its inverted terminal repeat (ITR) sequences. The ITRs of Turmoil-1 contain a conserved RNA recognition motif (RRM) that originated from the rsp-2 gene and a fragment from the protein-coding region of the cpg-3 gene. We further report that an open reading frame specific to C. elegans may have been created as a result of a Turmoil-1 insertion. Mutations at the 5' splice site of this open reading frame may have reactivated the transduplicated RRM motif.

Reviewers

This article was reviewed by Dan Graur and William Martin. For the full reviews, please go to the Reviewers' Reports section.

Findings

The possible contribution of transposable elements to the proteome has been discussed in several publications [1–7] and has provoked much debate [8]. Many mechanisms are known to increase the protein-domain repertoire, e.g., domain duplication, substitution mutations, insertions, deletions, and domain rearrangements [9]. In metazoans, a transposable element may result in transduction, in which a DNA segment downstream of transposable elements is mobilized as part of an aberrant transposition. This may result in gene duplication or exon shuffling, subsequently enriching the protein repertoire [10–13]. However, in the process of transduction, the transposable element does not acquire gene fragments as part of its sequence.

In plants, on the other hand, thousands of transposable elements contain duplicated gene fragments, captured in a process termed transduplication. Transduplication is a potentially rich source of novel coding sequences within rice and Arabidopsis thaliana [14–16]. Recently, transduplications of small nucleolar (sno) RNA by retroposon-like non-LTR transposable elements were found in the C. elegans [17] and platypus genomes [18].

The Harbinger superfamily of "cut-and-paste" DNA transposons was discovered through in silico studies [19]. This superfamily is characterized by Harbinger-specific transposases that are distantly related to the transposases encoded by the IS5-like group of bacterial transposons, such as IS5, IS112, and ISL2. Harbinger transposons are not as widespread as the eukaryotic hAT and mariner/Tc1 transposons; they are found in plants and nematodes [19–23] but not in mammals. Usually, Harbinger transposons are flanked by 3-bp target site duplications and 25- to 50-bp inverted terminal repeats (ITRs).

Turmoil-1 is a 5,024-bp long DNA transposon with 760-bp long ITRs and a Harbinger-specific transposase (Figure 1A). These ITRs are unique to Turmoil-1 and are not found in other members of the Harbinger superfamily of DNA transposons [24]. One complete copy of the Turmoil-1 was found on chromosome II of C. elegans; eight Turmoil-1 fragments exist in the genome (for detailed information, see Table 1). An analysis of C. elegans transposable elements [25, 26] revealed that a 205-bp ITR sequence within Turmoil-1 is highly similar to a region of two exons separated by an intron of the rsp-2 gene (see pairwise alignment using bl2seq [27] Figure 1B). These two exons encode the RNA Recognition Motif (RRM), which is found in many eukaryal and bacterial proteins. Specifically, the type of RRM domain present in the rsp-2 (called RRM1) gene is highly conserved evolutionarily [28]. The high similarity between the ITR sequence and the fragment of the rsp-2 gene implies that one originated from the other. The antiquity of this domain and a phylogenetic analysis (Figure 2) indicate that Turmoil-1 has recently acquired a portion of the rsp-2 gene sequence into its ITR. Tree reconstruction was performed with the PhyML program version 2.4.5 [29] using among-site rate variation with four discrete rate categories, and the JTT model [30] of sequence evolution.

Table 1 Turmoil-1 elements within the C. elegans genome
Figure 1
figure 1

Evidence for two transduplication events within the ITR sequence of the Turmoil -1 transposon. (A) A schematic illustration of the full sequence of Turmoil-1 on chromosome II (see table 1, copy number 6) and the structure and sites of inclusion of fragments of the rsp-2 and cpg-3 genes within the ITR sequence (B) Pairwise alignment (using bl2seq [27]) between the Turmoil-1 sequence and the rsp-2 sequence. The start codon, which is also the first amino acid of the RRM1 domain in the rsp-2 gene, is boxed in red. The intron sequence is indicated by lower-case letters; the protein-coding region is indicated by upper-case letters. The positions of the 5' and 3' splice sites (5'ss and 3'ss, respectively) are marked with arrows. The nucleotides in the 5' splice site of the rsp-2 gene that were mutated in the sequence of the transposable element, thereby abolishing the original 5' splice site, are shown in red. (C) Pairwise alignment between Turmoil-1 and the cpg-3 sequence.

Figure 2
figure 2

Maximum likelihood tree of the RRM domain. The RRM domain sequence, which is part of the ITR of Turmoil-1, is indicated in red. The tree shows that the RRM domain within Turmoil-1 is derived from the ancestral RRM domain rather than vice versa.

A comparative analysis of the rsp-2 gene and the 205-bp region of the gene found in the ITR sequence, revealed that the Turmoil-1 sequence has accumulated several point mutations within the 5' splice site that make it non-functional, whereas the 3' splice site of the intron remains intact (the mutations in the 5' splice site region are marked in red in Figure 1B). Since the RRM domain within the Turmoil-1 DNA transposon is not under purifying selection to maintain the reading frame or the functionality of the splice sites, these mutations are not unexpected.

Within the same ITR domain of Turmoil-1, and very close to the site of insertion of the RRM domain of rsp-2 gene, there is evidence of another "DNA kidnapping" event. A 131-bp fragment from the coding region of C. elegans gene cpg-3, which is unique to nematodes, was inserted into the ITR (Figure 1C). No sequences homologous to Turmoil-1 flank the cpg-3 gene. Thus, similar to the rsp-2 case, a gene fragment from the cpg-3 most likely was acquired by Turmoil-1, and not vice-versa. As the gene fragments are present on both sides of the ITR, capture may have occurred through non-homologous recombination.

One of the Turmoil-1 copies (number 1 in Table 1) contains within it an open reading frame (ORF) with the accession number Y48G1BL.4. It contains two putative exons and an intron (Figure 3), which are similar to the RRM domain. The 5' splice site that corresponds to that in the rsp-2 transcript has been mutated. A novel 5' splice site is most likely located nine nucleotides downstream from the original one. At this site, a point mutation changed an AT into a GT dinucleotide (marked in red in Figure 3). Usage of this 5' splice site maintains the ORF equivalent to that of the RRM domain of rsp-2 with the exception of the addition of three amino acids. These additional residues should have a negligible effect on the three-dimensional structure of the RRM domain (Figure 3C). This ORF, however, may not be transcriptionally active as its sequence is only found in the UNIPROT database (accession number Q9N3P9), and there is no EST or cDNA supporting evidence. If this is the case, it would be consistent with reports that indicate that all known transduplicates in rice, in spite of their genomic abundance, are pseudogenes [16].

Figure 3
figure 3

Y48G1BL.4 gene with predicted RRM domain. (A) A schematic illustration of the Y48G1BL.4 gene with predicted exons and an intron, and the position of the Turmoil-1 sequence within the genome. (B) Pairwise alignment (using bl2seq) between Turmoil-1 and the Y48G1BL.4 DNA sequence. The intron sequence is indicated by lower-case letters and the protein-coding region by upper-case letters. The potential new 5' splice site is marked in red and the 5' splice site position within the original RRM domain is marked with an arrow. (C) The predicted three-dimensional structure of the RRM domain generated from Y48G1BL.4.

Our analysis indicates that Turmoil-1 of C. elegans has captured two unrelated coding sequences within its ITRs at proximate locations. The presence of a transduplication "hotspot" in this region may be tentatively inferred. This analysis reveals a transduplication of protein-coding regions in C. elegans and strengthens the hypothesis that protein domains may be mobilized by transposable elements.

Reviewers' comments

Reviewer's report 1: Dan Graur, Department of Biology & Biochemistry University of Houston, Texas, USA

A very straightforward report – I have no other comments.

Reviewer's report 2: William Martin, Institut fuer Botanik III, Heinrich-Heine Universitaet Duesseldorf, Germany

This is an interesting and straightforward paper reporting the presence of transduplication in Caenorhabditis. The report of transduplication in animals would appear to be novel and certainly of sufficient interest to warrant publication. It might be the seed of a larger transduplication avalanche in animals, we'll see. I think the paper is fine for publication with the exception of "open read [ing] frame" in the abstract.

Author's response

Thanks for your comment – the typo was corrected.

Abbreviations

RRM:

RNA Recognition motif

ITR:

inverted terminal repeat

ORF:

open reading frame.

References

  1. Lev-Maor G, Sorek R, Shomron N, Ast G: The birth of an alternatively spliced exon: 3' splice-site selection in Alu exons. Science. 2003, 300 (5623): 1288-1291. 10.1126/science.1082588.

    Article  PubMed  CAS  Google Scholar 

  2. Sorek R, Ast G, Graur D: Alu-containing exons are alternatively spliced. Genome Res. 2002, 12 (7): 1060-1067. 10.1101/gr.229302.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  3. Krull M, Brosius J, Schmitz J: Alu-SINE exonization: en route to protein-coding function. Mol Biol Evol. 2005, 22 (8): 1702-1711. 10.1093/molbev/msi164.

    Article  PubMed  CAS  Google Scholar 

  4. Krull M, Petrusma M, Makalowski W, Brosius J, Schmitz J: Functional persistence of exonized mammalian-wide interspersed repeat elements (MIRs). Genome Res. 2007, 17 (8): 1139-1145. 10.1101/gr.6320607.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  5. Nekrutenko A, Li WH: Transposable elements are found in a large number of human protein-coding genes. Trends Genet. 2001, 17 (11): 619-621. 10.1016/S0168-9525(01)02445-3.

    Article  PubMed  CAS  Google Scholar 

  6. Sela N, Mersch B, Gal-Mark N, Lev-Maor G, Hotz-Wagenblatt A, Ast G: Comparative analysis of transposed elements' insertion within human and mouse genomes reveals Alu's unique role in shaping the human transcriptome. Genome Biol. 2007, 8 (6): R127-10.1186/gb-2007-8-6-r127.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Lorenc A, Makalowski W: Transposable elements and vertebrate protein diversity. Genetica. 2003, 118 (2–3): 183-191. 10.1023/A:1024105726123.

    Article  PubMed  CAS  Google Scholar 

  8. Gotea V, Makalowski W: Do transposable elements really contribute to proteomes?. Trends Genet. 2006, 22 (5): 260-267. 10.1016/j.tig.2006.03.006.

    Article  PubMed  CAS  Google Scholar 

  9. Chothia C, Gough J, Vogel C, Teichmann SA: Evolution of the protein repertoire. Science. 2003, 300 (5626): 1701-1703. 10.1126/science.1085371.

    Article  PubMed  CAS  Google Scholar 

  10. Pickeral OK, Makalowski W, Boguski MS, Boeke JD: Frequent human genomic DNA transduction driven by LINE-1 retrotransposition. Genome Res. 2000, 10 (4): 411-415. 10.1101/gr.10.4.411.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  11. Szak ST, Pickeral OK, Makalowski W, Boguski MS, Landsman D, Boeke JD: Molecular archeology of L1 insertions in the human genome. Genome Biol. 2002, 3 (10): research0052-10.1186/gb-2002-3-10-research0052.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Xing J, Wang H, Belancio VP, Cordaux R, Deininger PL, Batzer MA: Emergence of primate genes by retrotransposon-mediated sequence transduction. Proc Natl Acad Sci USA. 2006, 103 (47): 17608-17613. 10.1073/pnas.0603224103.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  13. Goodier JL, Ostertag EM, Kazazian HH: Transduction of 3'-flanking sequences is common in L1 retrotransposition. Hum Mol Genet. 2000, 9 (4): 653-657. 10.1093/hmg/9.4.653.

    Article  PubMed  CAS  Google Scholar 

  14. Hoen DR, Park KC, Elrouby N, Yu Z, Mohabir N, Cowan RK, Bureau TE: Transposon-mediated expansion and diversification of a family of ULP-like genes. Mol Biol Evol. 2006, 23 (6): 1254-1268. 10.1093/molbev/msk015.

    Article  PubMed  CAS  Google Scholar 

  15. Jiang N, Bao Z, Zhang X, Eddy SR, Wessler SR: Pack-MULE transposable elements mediate gene evolution in plants. Nature. 2004, 431 (7008): 569-573. 10.1038/nature02953.

    Article  PubMed  CAS  Google Scholar 

  16. Juretic N, Hoen DR, Huynh ML, Harrison PM, Bureau TE: The evolutionary fate of MULE-mediated duplications of host gene fragments in rice. Genome Res. 2005, 15 (9): 1292-1297. 10.1101/gr.4064205.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  17. Zemann A, op de Bekke A, Kiefmann M, Brosius J, Schmitz J: Evolution of small nucleolar RNAs in nematodes. Nucleic Acids Res. 2006, 34 (9): 2676-2685. 10.1093/nar/gkl359.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  18. Schmitz J, Zemann A, Churakov G, Kuhl H, Grutzner F, Reinhardt R, Brosius J: Retroposed SNOfall – a mammalian-wide comparison of platypus snoRNAs. Genome Res. 2008, 18 (6): 1005-1010. 10.1101/gr.7177908.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  19. Kapitonov VV, Jurka J: Molecular paleontology of transposable elements from Arabidopsis thaliana. Genetica. 1999, 107 (1–3): 27-37. 10.1023/A:1004030922447.

    Article  PubMed  CAS  Google Scholar 

  20. Jurka J, Kapitonov VV: PIFs meet Tourists and Harbingers: a superfamily reunion. Proc Natl Acad Sci USA. 2001, 98 (22): 12315-12316. 10.1073/pnas.231490598.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  21. Zhang X, Feschotte C, Zhang Q, Jiang N, Eggleston WB, Wessler SR: P instability factor: an active maize transposon system associated with the amplification of Tourist-like MITEs and a new superfamily of transposases. Proc Natl Acad Sci USA. 2001, 98 (22): 12572-12577. 10.1073/pnas.211442198.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  22. Jiang N, Bao Z, Zhang X, Hirochika H, Eddy SR, McCouch SR, Wessler SR: An active DNA transposon family in rice. Nature. 2003, 421 (6919): 163-167. 10.1038/nature01214.

    Article  PubMed  CAS  Google Scholar 

  23. Kikuchi K, Terauchi K, Wada M, Hirano HY: The plant MITE mPing is mobilized in anther culture. Nature. 2003, 421 (6919): 167-170. 10.1038/nature01218.

    Article  PubMed  CAS  Google Scholar 

  24. Kapitonov VV, Jurka J: Harbinger transposons and an ancient HARBI1 gene derived from a transposase. DNA Cell Biol. 2004, 23 (5): 311-324. 10.1089/104454904323090949.

    Article  PubMed  CAS  Google Scholar 

  25. Levy A, Sela N, Ast G: TranspoGene and microTranspoGene: transposed elements influence on the transcriptome of seven vertebrates and invertebrates. Nucleic Acids Res. 2008, D47-52. 36 Database

  26. Transpogene Database. [http://transpogene.tau.ac.il/]

  27. Tatusova TA, Madden TL: BLAST 2 Sequences, a new tool for comparing protein and nucleotide sequences. FEMS Microbiol Lett. 1999, 174 (2): 247-250. 10.1111/j.1574-6968.1999.tb13575.x.

    Article  PubMed  CAS  Google Scholar 

  28. Longman D, Johnstone IL, Caceres JF: Functional characterization of SR and SR-related genes in Caenorhabditis elegans. Embo J. 2000, 19 (7): 1625-1637. 10.1093/emboj/19.7.1625.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  29. Guindon S, Lethiec F, Duroux P, Gascuel O: PHYML Online – a web server for fast maximum likelihood-based phylogenetic inference. Nucleic Acids Res. 2005, W557-559. 10.1093/nar/gki352. 33 Web Server

  30. Jones DT, Taylor WR, Thornton JM: The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci. 1992, 8 (3): 275-282.

    PubMed  CAS  Google Scholar 

Download references

Acknowledgements

We thank Prof. Jerzy Jurka for critical reading of the manuscript. This work was supported by the Israeli Ministry of Science and Technology (MOST) and by grants from the Israel Science Foundation (1449/04 and 40/05), MOP Germany-Israel, GIF, ICA through the Ber-Lehmsdorf Memorial Fund, and DIP and EURASNET. AS is a fellow of the Complexity Science Scholarship program and is supported by a fellowship from the Israeli Ministry of Science.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Gil Ast.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

NS and AS did the experiments and analysis. WM interpreted the results. GA and TP supervised the study. NS, AS, WM, TP and GA wrote the paper.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Sela, N., Stern, A., Makalowski, W. et al. Transduplication resulted in the incorporation of two protein-coding sequences into the Turmoil-1 transposable element of C. elegans. Biol Direct 3, 41 (2008). https://doi.org/10.1186/1745-6150-3-41

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1745-6150-3-41

Keywords