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

  • Noa Sela1,

    Affiliated with

    • Adi Stern2,

      Affiliated with

      • Wojciech Makalowski3,

        Affiliated with

        • Tal Pupko2 and

          Affiliated with

          • Gil Ast1Email author

            Affiliated with

            Biology Direct20083:41

            DOI: 10.1186/1745-6150-3-41

            Received: 03 October 2008

            Accepted: 08 October 2008

            Published: 08 October 2008

            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 [17] 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 [1013]. 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 [1416]. 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 [1923] 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

            #

            Chr (*)

            Start (**)

            End (**)

            orientation

            Containing RRM domain (***)

            Containing cpg-3 fragment

            Length of the element (****)

            1

            I

            216923

            217647

            -

            + (1–827)

            +

            PL

            2

            I

            5979173

            5980185

            +

            + (1–1177)

            +

            PL

            3

            I

            5980183

            5982782

            -

            + (1–2771)

            +

            PL

            4

            I

            7245338

            7245520

            +

            - (125–317)

            +

            PL

            5

            I

            13973017

            13973196

            +

            - (4691–4882)

            +

            PL

            6

            II

            3351574

            3356982

            -

            + (1–5024)

            +

            FL

            7

            II

            10883165

            10884063

            -

            - (519–1090)

            -

            PL

            8

            II

            11452621

            11452867

            -

            + (323–568)

            -

            PL

            9

            IV

            2900974

            2901153

            +

            - (125–316)

            +

            PL

            (*) chromosome number

            (**) Start and End according to the coordinates of C. elegans genome version ce6 (May 2008).

            (***) in brackets the start and end relative to the full sequence of the Turmoil-1 found on chromosome II 3351574 – 3356982 (#6).

            (****) PL – partial length of Turmoil-1/FL – full length of Turmoil-1.

            http://static-content.springer.com/image/art%3A10.1186%2F1745-6150-3-41/MediaObjects/13062_2008_Article_119_Fig1_HTML.jpg
            Figure 1

            Evidence for two transduplication events within the ITR sequence of theTurmoil-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.

            http://static-content.springer.com/image/art%3A10.1186%2F1745-6150-3-41/MediaObjects/13062_2008_Article_119_Fig2_HTML.jpg
            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].
            http://static-content.springer.com/image/art%3A10.1186%2F1745-6150-3-41/MediaObjects/13062_2008_Article_119_Fig3_HTML.jpg
            Figure 3

            Y48G1BL.4gene 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.

            Declarations

            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.

            Authors’ Affiliations

            (1)
            Department of Human Molecular Genetics and Biochemistry, Sackler Faculty of Medicine, Tel Aviv University
            (2)
            Department of Cell Research and Immunology, George S. Wise Faculty of Life Sciences, Tel Aviv University
            (3)
            Institute of Bioinformatics, Faculty of Medicine, University of Muenster

            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.View ArticlePubMed
            2. Sorek R, Ast G, Graur D: Alu-containing exons are alternatively spliced. Genome Res 2002,12(7):1060–1067.View ArticlePubMed
            3. Krull M, Brosius J, Schmitz J: Alu-SINE exonization: en route to protein-coding function. Mol Biol Evol 2005,22(8):1702–1711.View ArticlePubMed
            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.View ArticlePubMed
            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.View ArticlePubMed
            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.View ArticlePubMed
            7. Lorenc A, Makalowski W: Transposable elements and vertebrate protein diversity. Genetica 2003,118(2–3):183–191.View ArticlePubMed
            8. Gotea V, Makalowski W: Do transposable elements really contribute to proteomes? Trends Genet 2006,22(5):260–267.View ArticlePubMed
            9. Chothia C, Gough J, Vogel C, Teichmann SA: Evolution of the protein repertoire. Science 2003,300(5626):1701–1703.View ArticlePubMed
            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.View ArticlePubMed
            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.View ArticlePubMed
            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.View ArticlePubMed
            13. Goodier JL, Ostertag EM, Kazazian HH Jr: Transduction of 3'-flanking sequences is common in L1 retrotransposition. Hum Mol Genet 2000,9(4):653–657.View ArticlePubMed
            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.View ArticlePubMed
            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.View ArticlePubMed
            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.View ArticlePubMed
            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.View ArticlePubMed
            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.View ArticlePubMed
            19. Kapitonov VV, Jurka J: Molecular paleontology of transposable elements from Arabidopsis thaliana. Genetica 1999,107(1–3):27–37.View ArticlePubMed
            20. Jurka J, Kapitonov VV: PIFs meet Tourists and Harbingers: a superfamily reunion. Proc Natl Acad Sci USA 2001,98(22):12315–12316.View ArticlePubMed
            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.View ArticlePubMed
            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.View ArticlePubMed
            23. Kikuchi K, Terauchi K, Wada M, Hirano HY: The plant MITE mPing is mobilized in anther culture. Nature 2003,421(6919):167–170.View ArticlePubMed
            24. Kapitonov VV, Jurka J: Harbinger transposons and an ancient HARBI1 gene derived from a transposase. DNA Cell Biol 2004,23(5):311–324.View ArticlePubMed
            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, (36 Database):D47–52.
            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.View ArticlePubMed
            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.View ArticlePubMed
            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, (33 Web Server):W557–559.
            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

            Copyright

            © Sela et al. 2008

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

            Advertisement