MBE Advance Access originally published online on June 22, 2005
Molecular Biology and Evolution 2005 22(10):2048-2054; doi:10.1093/molbev/msi200
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Research Article |
Spliced-Leader trans-Splicing in Freshwater Planarians
Department of Cell and Developmental Biology, University of Illinois at Urbana, Champaign
E-mail: pnewmark{at}life.uiuc.edu.
| Abstract |
|---|
|
|
|---|
trans-Splicing, in which a spliced-leader (SL) RNA is appended to the most 5' exon of independently transcribed pre-mRNAs, has been described in a wide range of eukaryotes, from protozoans to chordates. Here we describe trans-splicing in the freshwater planarian Schmidtea mediterranea, a free-living member of the phylum Platyhelminthes. Analysis of an expressed sequence tag (EST) collection from this organism showed that over 300 transcripts shared one of two
35-base sequences (Smed SL-1 and SL-2) at their 5' ends. Examination of genomic sequences encoding representatives of these transcripts revealed that these shared sequences were transcribed elsewhere in the genome. RNA blot analysis, 5' and 3' rapid amplification of cDNA ends, as well as genomic sequence data showed that 42-nt SL sequences were derived from small RNAs of
110 nt. Similar sequences were also found at the 5' ends of ESTs from the planarian Dugesia japonica. trans-Splicing has already been described in numerous representatives of the phylum Platyhelminthes (trematodes, cestodes, and polyclads); its presence in two representatives of the triclads supports the hypothesis that this mode of RNA processing is ancestral within this group. The upcoming complete genome sequence of S. mediterranea, combined with this animal's experimental accessibility and susceptibility to RNAi, provide another model organism in which to study the function of the still-enigmatic trans-splicing.
Key Words: spliced-leader RNA trans-splicing planarian Platyhelminthes Schmidtea mediterranea Dugesia japonica
| Introduction |
|---|
|
|
|---|
Spliced-leader trans-splicing is an unconventional form of mRNA processing in which an independently transcribed leader sequence is appended to the 5' end of pre-mRNAs (Davis 1996
While analyzing a collection of expressed sequence tags (ESTs) from a free-living flatworm, the freshwater planarian Schmidtea mediterranea (R. M. Zayas, A. Hernández, B. Habermann, Y. Wang, and P. A. Newmark, in preparation) we noticed a stretch of
35 nt that was shared between the 5' ends of hundreds of different EST clones. Here we present evidence that these sequences are derived from trans-splicing of SL RNAs to a population of transcripts in S. mediterranea. We also show that similar SL sequences are found in another planarian, Dugesia japonica, providing further evidence that trans-splicing may be an ancestral feature within the Platyhelminthes. The planarian's experimental accessibility, susceptibility to RNAi, ease of culture, and the forthcoming genome sequence of S. mediterranea provide another model organism in which to study the function(s) of trans-splicing.
| Materials and Methods |
|---|
|
|
|---|
Informatics
Potential trans-spliced transcripts were identified by BlastN (Stand-alone BLAST 2.2.8 [Altschul et al. 1997
trans-Spliced cDNAs were examined by BlastX against the nonredundant protein database at NCBI. Those with significant matches were classified further by Gene Ontology (GO), using terms in the Biological Process category based on the closest annotated homolog (Ashburner et al. 2000
; Gene Ontology Consortium 2001
). RPS-Blast (Marchler-Bauer et al. 2005
) was used to query the conserved domain database for conserved domains encoded by the predicted products of trans-spliced RNAs. All ESTs containing the Smed SL sequences have been deposited in GenBank; accession numbers for these sequences are available as a supplementary file (Supplementary Material online).
RNA Isolation
Total RNA was isolated from a clonal line (CIWsx2) of the hermaphroditic strain of S. mediterranea. Worms were placed on ice in RNAlater (Ambion, Austin, Tex.), and then Totally RNA (Ambion) was used to extract RNA as recommended by the supplier. Total RNA obtained in this manner was LiCl precipitated prior to use.
Rapid Amplification of cDNA Ends
Rapid amplification of cDNA ends (RACE) experiments were performed using FirstChoice RLM-RACE (Ambion). For 5' RACE gene-specific primers were designed from S. mediterranea genomic DNA using the predicted SL RNA sequence. The reverse primers used were 5'-CAAGTGACTGTCAAAAATTAACC-3' and 5'-CTAATGTTGGATAACGGTCC-3' for the outer and inner amplification reactions, respectively. For SL 3' RACE, poly(A) tails were added to 10 µg S. mediterranea total RNA using yeast poly(A) polymerase (Ambion); first-strand cDNA was then synthesized from the polyadenylated RNA. The gene-specific forward primer used for 3' RACE was 5'-GACGGTCTTATCGAAATC-3'.
An
70-bp fragment obtained from the second round of 5' RACE amplification and an
100-bp fragment obtained from 3' RACE were gel purified using Qiaex II resin (Qiagen, Valencia, Calif.), cloned into pCR II T-A cloning vector, and transformed into Oneshot TOP10F' competent cells (Invitrogen, Carlsbad, Calif.). Clones were checked for inserts by colony polymerase chain reaction; plasmid DNA was then purified by miniprep (Wizard Plus SV minipreps, Promega, Madison, Wisc.). Inserts were sequenced using the standard Big Dye 3.1 sequencing reaction (ABI, Foster City, Calif.) with the addition of 0.7 M betaine (Sigma, St. Louis, Mo.) and were analyzed using Sequencher 4.2.2 (Gene Codes Co., Ann Arbor, Mich.).
Northern Blot Analysis
Total RNA (10 µg) was separated in a formaldehyde/agarose gel. After electrophoresis, RNA was capillary transferred to Hybond N+ nylon membrane (Amersham, Piscataway, N.J.) and UV cross-linked to the membrane using a Spectrolinker XL-1500 (Spectronics Corporation, Westbury, N.Y.). The sequence of the oligonucleotide probe was 5'-GACGGTCTTATCGAAATCTATATAAATC-3'; 100
mol were 3'-end labeled with digoxigenin-ddUTP (Roche, Indianapolis, Ind.) using terminal transferase (Roche). Hybridizations were carried out in DIG Easy Hyb (Roche) for 18 h at room temperature. The blot was washed twice in 2 x saline sodium citrate (SSC)/0.1% sodium dodecyl sulfate (SDS) and twice in 0.5 x SSC/0.1% SDS at room temperature. CDP-Star (Roche), a chemiluminescent substrate for alkaline phosphatase, was used to detect hybridized probe. The blot was incubated in 1% Blocking Reagent (Roche) in maleic acid buffer (MAB; 0.1 M maleic acid, 0.15 NaCl, pH 7.5) containing 1:20,000 Anti-Digoxigenin-AP (Roche), washed in MAB/0.3% Tween 20, equilibrated in detection buffer (0.1 M Tris-HCl, 0.1 M NaCl, pH 7.5), and then incubated in CDP-Star working solution for 5 min before exposure to X-ray film.
RNA Secondary Structure Modeling and Sequence Alignments
The SL sequences derived from RACE experiments were modeled using web-based mfold 2.3 (Zuker 2003
). Folding was performed at 18°C and 22°C (the temperatures at which the planarians are reared in the laboratory) as well as 37°C with the constraint that the putative Sm-binding site (AAUUUUUGA) remained single stranded. To generate graphical output, the mfold predictions were imported into RnaViz 2.0 (De Rijk, Wuyts, and De Wachter 2003
). ClustalW 1.8 (Thompson, Higgins, and Gibson 1994
; Chenna et al. 2003
) was used to produce multiple sequence alignments, and graphical display of the output was performed using Boxshade 3.21.
| Results and Discussion |
|---|
|
|
|---|
Identification of Putative SL Sequences in S. mediterranea ESTs
To study the process of epigenetic germ cell specification, we have generated an EST collection from the hermaphroditic strain of the planarian S. mediterranea. This collection of
27,000 ESTs assembles into
10,000 putative transcripts,
6,500 of which are present in contiguous sequences (contigs) of two or more ESTs (R. M. Zayas, A. Hernández, B. Habermann, Y. Wang, and P. A. Newmark, in preparation). While comparing these sequences to an EST collection generated from the asexual strain of the same species (Sánchez Alvarado et al. 2002
35 nt with the 5' ends of
25 ESTs from the asexual planarian strain. The simplest explanation for this shared 5' sequence was that it resulted from a trans-splicing event in which an SL sequence was appended to numerous different transcripts.
To identify the potential targets of SL addition within our EST collection we used the shared 35-nt sequence to query the assembled ESTs using BlastN. This search identified 320 putative targets of trans-splicing or
3.2% of the transcripts represented in the EST collection; this number likely underestimates the frequency of the occurrence of the putative SL sequence because its identification requires full-length or nearly full-length cDNAs. Closer inspection of the putative SL sequences revealed two distinct types: one present on 304 of the assembled sequences and another present on 16 of the sequences. Because of its abundance the former will be referred to as Smed SL-1 (fig. 1A), whereas the less prevalent form will be called Smed SL-2 (fig. 1B); they differ from each other by 2 nt (fig. 1A and B). Upon inspection of the contigs that contained SL-2, we found four in which either SL-1 or SL-2 was observed at the 5' end of the aligned ESTs; an additional eight SL-2containing ESTs were found aligned in contigs with either SL-1 or putative 5' cis-sequences (see examples in fig. 1C in which SL-1, SL-2, and cis-sequences are all represented by ESTs). Examination of genomic DNA sequences encoding representatives of the potentially trans-spliced mRNAs revealed that the putative SL sequences were not present in the vicinity of the transcribed sequences; rather, potential splice acceptor sites were identified at the presumptive site of SL addition (fig. 1C). These results are all consistent with trans-splicing appending SL sequences to a sizable population of planarian transcripts; furthermore, it appears that mRNAs encoded by the same gene may receive either SL sequence. It will be of interest to determine if the utilization of different SL sequences has any functional significance, for example, whether they play roles in distinct cell types or at different stages of development.
|
Characterization of Planarian SL RNAs
RNA blot analysis of total RNA hybridized with an antisense oligonucleotide to the predicted SL sequence revealed a prominent band of
100 bases as well as a smear ranging in size 0.38 kb (fig. 2A). This result is consistent with a heterogeneous population of RNAs being the recipient of SL trans-splicing, with the SL derived from a small RNA. 5' and 3' RACE were performed to characterize the termini of the
100-base SL RNA (see Materials and Methods). These sequences were then compared with genomic SL sequences identified by BlastN analysis of S. mediterranea whole-genome shotgun reads available via the NCBI Trace Archives. Two distinct sequences were identified that corresponded to SL-1 and SL-2 and that matched the sequences obtained by RACE (fig. 2B): SL-1 RNA is 107 bases in length and SL-2 is 106 bases long; both produce a 42-base spliced leader. Additional genomic sequences encoding potential SL-1-type RNAs were identified from the whole-genome shotgun reads. In total, we identified seven distinct genomic contigs that contain SL sequences, one of which contains three distinct SL repeats; an alignment of these sequences is shown in figure 2C. None of these sequences appears to be located within 5S RNA gene repeats (not shown), consistent with observations of SL genes in other flatworms (Davis 1997
|
Secondary structure predictions for the Smed SL sequences were generated using mfold version 2.3 (Zuker 2003
|
Conservation of SL Sequences in D. japonica
To address the question of whether trans-splicing was more widespread within the freshwater planarians, we used the SL sequence to query Paludicola (freshwater planarian) sequences in GenBank. In addition to identifying 30 SL-containing ESTs from the asexual strain of S. mediterranea, we found 72 ESTs from the planarian D. japonica (Mineta et al. 2003
|
Annotation of trans-Spliced mRNAs
To examine potential functional relationships between the products encoded by trans-spliced mRNAs in planarians, we used these sequences as queries in BlastX searches against the nonredundant protein database (nr, NCBI). Of 320 trans-spliced sequences from the hermaphrodite S. mediterranea, 242 had significant (E < 1 x 104) hits to the nr database. Based on the closest annotated homolog, ESTs were assigned GO terms in the Biological Process category; 110 ESTs could be assigned to a Biological Process (Supplementary Spreadsheet 1, Supplementary Material online), and the major categories are summarized in table 1. We found that there were no obvious trends in the predicted biological functions of the products of trans-spliced mRNAs, similar to what has been described in other flatworms (Davis et al. 1995
149 different domains represented (Supplementary Spreadsheet 1, Supplementary Material online). Similar to the GO annotation, the distribution of domains suggests that trans-splicing is widespread and not specific to particular gene families or functional categories.
|
In this report we have provided evidence for trans-splicing of SL RNAs in freshwater planarians. Planarians are best known for their remarkable regenerative abilities (Reddien and Sánchez Alvarado 2004
| Supplementary Material |
|---|
|
|
|---|
The complete Gene Ontology term annotation and conserved protein domain assignment results are available as a supplementary spreadsheet (Zayas et al-S1.xls). The accession numbers for all of the trans-spliced ESTs reported here are also available as a supplementary file (SmedSL_accession.txt). The supplementary spreadsheet 1 and the supplementary file are available at Molecular Biology and Evolution online (http://www.mbe.oxfordjournals.org/).
| Acknowledgements |
|---|
|
|
|---|
We would like to thank Francesc Cebrià and Richard Davis for helpful comments on the manuscript, Bianca Habermann for annotation of planarian EST sequences, the Washington University Genome Sequencing Center and the National Human Genome Research Institute for shotgun sequencing the S. mediterranea genome, as well as Jeffrey Haas and Phil Anders for help with computer analysis. R.M.Z. is a fellow of the Jane Coffin Childs Memorial Fund for Medical Research. This work was supported by an NSF CAREER Award (IBN-0237825) and National Institutes of Health R01 HD-43403 to P.A.N. P.A.N. is a Damon Runyon Scholar supported by the Damon Runyon Cancer Research Foundation (DRS 33-03).
| Footnotes |
|---|
1 These authors contributed equally to this work.
Billie Swalla, Associate Editor
| References |
|---|
|
|
|---|
Altschul, S. F., T. L. Madden, A. A. Schaffer, J. Zhang, Z. Zhang, W. Miller, and D. J. Lipman. 1997. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25:33893402.
Ashburner, M., C. A. Ball, J. A. Blake et al. (20 co-authors). 2000. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat. Genet. 25:2529.[CrossRef][Web of Science][Medline]
Baguñà, J., and M. Riutort. 2004. Molecular phylogeny of the Platyhelminthes. Can. J. Zool. 82:168193.[CrossRef]
Blumenthal, T. 1995. Trans-splicing and polycistronic transcription in Caenorhabditis elegans. Trends Genet. 11:132136.[CrossRef][Web of Science][Medline]
Brehm, K., K. Hubert, E. Sciutto, T. Garate, and M. Frosch. 2002. Characterization of a spliced leader gene and of trans-spliced mRNAs from Taenia solium. Mol. Biochem. Parasitol. 122:105110.[CrossRef][Web of Science][Medline]
Brehm, K., K. Jensen, and M. Frosch. 2000. mRNA trans-splicing in the human parasitic cestode Echinococcus multilocularis. J. Biol. Chem. 275:3831138318.
Bruzik, J. P., K. Van Doren, D. Hirsh, and J. A. Steitz. 1988. Trans splicing involves a novel form of small nuclear ribonucleoprotein particles. Nature 335:559562.[CrossRef][Medline]
Campbell, D. A., D. A. Thornton, and J. C. Boothroyd. 1984. Apparent discontinuous transcription of Trypanosoma brucei variant surface antigen genes. Nature 311:350355.[CrossRef][Medline]
Chenna, R., H. Sugawara, T. Koike, R. Lopez, T. J. Gibson, D. G. Higgins, and J. D. Thompson. 2003. Multiple sequence alignment with the Clustal series of programs. Nucleic Acids Res. 31:34973500.
Davis, R. E. 1996. Spliced leader RNA trans-splicing in metazoa. Parasitol. Today 12:3340.[CrossRef][Web of Science][Medline]
. 1997. Surprising diversity and distribution of spliced leader RNAs in flatworms. Mol. Biochem. Parasitol. 87:2948.[CrossRef][Web of Science][Medline]
Davis, R. E., C. Hardwick, P. Tavernier, S. Hodgson, and H. Singh. 1995. RNA trans-splicing in flatworms. Analysis of trans-spliced mRNAs and genes in the human parasite, Schistosoma mansoni. J. Biol. Chem. 270:2181321819.
De Rijk, P., J. Wuyts, and R. De Wachter. 2003. RnaViz 2: an improved representation of RNA secondary structure. Bioinformatics 19:299300.
Ganot, P., T. Kallesoe, R. Reinhardt, D. Chourrout, and E. M. Thompson. 2004. Spliced-leader RNA trans splicing in a chordate, Oikopleura dioica, with a compact genome. Mol. Cell. Biol. 24:77957805.
Gene Ontology Consortium. 2001. Creating the gene ontology resource: design and implementation. Genome Res. 11:14251433.
Kooter, J. M., T. De Lange, and P. Borst. 1984. Discontinuous synthesis of mRNA in trypanosomes. EMBO J. 3:23872392.[Web of Science][Medline]
Krause, M., and D. Hirsh. 1987. A trans-spliced leader sequence on actin mRNA in C. elegans. Cell 49:753761.[CrossRef][Web of Science][Medline]
Lee, K. Z., and R. J. Sommer. 2003. Operon structure and trans-splicing in the nematode Pristionchus pacificus. Mol. Biol. Evol. 20:20972103.
Marchler-Bauer, A., J. B. Anderson, P. F. Cherukuri et al. (24 co-authors). 2005. CDD: a Conserved Domain Database for protein classification. Nucleic Acids Res. 33(Database Issue):D192D196.
Maroney, P. A., J. A. Denker, E. Darzynkiewicz, R. Laneve, and T. W. Nilsen. 1995. Most mRNAs in the nematode Ascaris lumbricoides are trans-spliced: a role for spliced leader addition in translational efficiency. RNA 1:714723.[Abstract]
Milhausen, M., R. G. Nelson, S. Sather, M. Selkirk, and N. Agabian. 1984. Identification of a small RNA containing the trypanosome spliced leader: a donor of shared 5' sequences of trypanosomatid mRNAs? Cell 38:721729.[CrossRef][Web of Science][Medline]
Mineta, K., M. Nakazawa, F. Cebrià, K. Ikeo, K. Agata, and T. Gojobori. 2003. Origin and evolutionary process of the CNS elucidated by comparative genomics analysis of planarian ESTs. Proc. Natl. Acad. Sci. USA 100:76667671.
Newmark, P. A., and A. Sánchez Alvarado. 2002. Not your father's planarian: a classic model enters the era of functional genomics. Nat. Rev. Genet. 3:210219.[CrossRef][Web of Science][Medline]
Nilsen, T. W. 2001. Evolutionary origin of SL-addition trans-splicing: still an enigma. Trends Genet. 17:678680.[CrossRef][Web of Science][Medline]
Pouchkina-Stantcheva, N. N., and A. Tunnacliffe. 2005. Spliced leader RNA-mediated trans-splicing in Phylum Rotifera. Mol. Biol. Evol. 22:14821489.
Rajkovic, A., R. E. Davis, J. N. Simonsen, and F. M. Rottman. 1990. A spliced leader is present on a subset of mRNAs from the human parasite Schistosoma mansoni. Proc. Natl. Acad. Sci. USA 87:88798883.
Reddien, P. W., A. L. Bermange, K. J. Murfitt, J. R. Jennings, and A. Sánchez Alvarado. 2005. Identification of genes needed for regeneration, stem cell function, and tissue homeostasis by systematic gene perturbation in planaria. Dev. Cell 8:635649.[CrossRef][Web of Science][Medline]
Reddien, P. W., and A. Sánchez Alvarado. 2004. Fundamentals of planarian regeneration. Annu. Rev. Cell Dev. Biol. 20:725757.[CrossRef][Web of Science][Medline]
Sánchez Alvarado, A., and P. A. Newmark. 1999. Double-stranded RNA specifically disrupts gene expression during planarian regeneration. Proc. Natl. Acad. Sci. USA 96:50495054.
Sánchez Alvarado, A., P. A. Newmark, S. M. Robb, and R. Juste. 2002. The Schmidtea mediterranea database as a molecular resource for studying platyhelminthes, stem cells and regeneration. Development 129:56595665.
Stover, N. A., and R. E. Steele. 2001. Trans-spliced leader addition to mRNAs in a cnidarian. Proc. Natl. Acad. Sci. USA 98:56935698.
Tessier, L. H., M. Keller, R. L. Chan, R. Fournier, J. H. Weil, and P. Imbault. 1991. Short leader sequences may be transferred from small RNAs to pre-mature mRNAs by trans-splicing in Euglena. EMBO J. 10:26212625.[Web of Science][Medline]
Thompson, J. D., D. G. Higgins, and T. J. Gibson. 1994. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22:46734680.
Umesono, Y., K. Watanabe, and K. Agata. 1997. A planarian orthopedia homolog is specifically expressed in the branch region of both the mature and regenerating brain. Dev. Growth Differ. 39:723727.[CrossRef][Web of Science][Medline]
Vandenberghe, A. E., T. H. Meedel, and K. E. Hastings. 2001. mRNA 5'-leader trans-splicing in the chordates. Genes Dev. 15:294303.
Zeiner, G. M., N. R. Sturm, and D. A. Campbell. 2003. The Leishmania tarentolae spliced leader contains determinants for association with polysomes. J. Biol. Chem. 278:3826938275.
Zuker, M. 2003. Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 31:34063415.
![]()
CiteULike
Connotea
Del.icio.us What's this?
This article has been cited by other articles:
![]() |
H. Zhang, D. A. Campbell, N. R. Sturm, and S. Lin Dinoflagellate Spliced Leader RNA Genes Display a Variety of Sequences and Genomic Arrangements Mol. Biol. Evol., August 1, 2009; 26(8): 1757 - 1771. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Pettitt, B. Muller, I. Stansfield, and B. Connolly Spliced leader trans-splicing in the nematode Trichinella spiralis uses highly polymorphic, noncanonical spliced leaders RNA, April 1, 2008; 14(4): 760 - 770. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Zhang, Y. Hou, L. Miranda, D. A. Campbell, N. R. Sturm, T. Gaasterland, and S. Lin Spliced leader RNA trans-splicing in dinoflagellates PNAS, March 13, 2007; 104(11): 4618 - 4623. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Cheng, L. Cohen, D. Ndegwa, and R. E. Davis The Flatworm Spliced Leader 3'-Terminal AUG as a Translation Initiator Methionine J. Biol. Chem., January 13, 2006; 281(2): 733 - 743. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||







