Bioinformatics Advance Access originally published online on January 29, 2009
Bioinformatics 2009 25(5):674-675; doi:10.1093/bioinformatics/btp020
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
SECISaln, a web-based tool for the creation of structure-based alignments of eukaryotic SECIS elements
1Institut Municipal d'Investigació Médica, 2Centre de Regulació Genómica, Universitat Pompeu Fabra and Parc de Recerca Biomedica de Barcelona, Carrer del Doctor Aiguader 88, 08003, Barcelona, Catalonia, Spain and 3Unité Architecture et Réactivité de l'ARN, Université Louis Pasteur de Strasbourg, CNRS, 15 rue René Descartes, F-67084 Strasbourg, France
*To whom correspondence should be addressed.
 |
ABSTRACT |
|---|
 TOP ABSTRACT ACKNOWLEDGEMENTSREFERENCES |
|---|
Summary: Selenoproteins contain the 21st amino acid selenocysteine which is encoded by an inframe UGA codon, usually read as a stop. In eukaryotes, its co-translational recoding requires the presence of an RNA stem–loop structure, the SECIS element in the 3 untranslated region of (UTR) selenoprotein mRNAs. Despite little sequence conservation, SECIS elements share the same overall secondary structure. Until recently, the lack of a significantly high number of selenoprotein mRNA sequences hampered the identification of other potential sequence conservation. In this work, the web-based tool SECISaln provides for the first time an extensive structure-based sequence alignment of SECIS elements resulting from the well-defined secondary structure of the SECIS RNA and the increased size of the eukaryotic selenoproteome. We have used SECISaln to improve our knowledge of SECIS secondary structure and to discover novel, conserved nucleotide positions and we believe it will be a useful tool for the selenoprotein and RNA scientific communities.
Availability: SECISaln is freely available as a web-based tool at http://genome.crg.es/software/secisaln/.
Contact: charles.chapple{at}crg.es
Supplementary information: Supplementary data are available at Bioinformatics online.
Selenoproteins are a diverse family of proteins characterized by the presence of the 21st amino acid, selenocysteine (Sec or U). Selenocysteine is co-translationally inserted into the growing polypeptide chain in response to UGA, otherwise read as a stop codon. The correct recoding of UGA to Sec requires the presence of a stem-loop structure, the SECIS element in the 3 untranslated region (UTR) of selenoprotein gene transcripts. Accordingly, the presence of a suitable SECIS element has been used in many studies as a tool for the computational prediction of novel selenoproteins (Castellano et al., 2001; Kryukov et al., 1999; Lescure et al.; 1999) and a specialized tool for SECIS prediction, SECISearch (Kryukov et al., 2003), has already been described and has been widely used.
There are two types of eukaryotic SECISes, type I and type II differing at the apex by the presence of the additional helix 3 in type II (Fagegaltier et al., 2000; Grundner-Culemann et al., 1999; Walczak et al., 1996 see Fig. 1. Although the SECIS structure is conserved, there is little sequence conservation beyond the consecutive non-Watson-Crick base pairs UGAN/KGAW constituting the quartet, an unpaired A 5 to UGAN and a run of As in the apical loop/internal loop 2 (Fagegaltier et al., 2000; Walczak et al., 1996). Of these only the UGA/GA of the quartet is invariable1 (e.g. Buettner et al., 1996; Lobanov et al., 2007). Here, we describe SECISaln, a web-based tool that creates structure-based alignments of an extensive dataset of eukaryotic SECIS sequences. Its implementation led us to uncover novel, conserved sequence elements.
|
SECISaln will predict a SECIS element in the query sequence, split it into its constituent parts and align these against a precompiled database of eukaryotic SECIS elements. The user can choose whether the database sequences are sorted by protein family or by species, thereby offering the possibility of comparing the submitted sequence to other, known SECISes. In addition, SECISaln returns a graphical image of the predicted structure of the user-submitted sequence as well as a multiple structural alignment of all SECIS elements of that type already present in the database. SECISaln uses SECISearch for the SECIS prediction step, described in detail in (Kryukov et al. (2003) and is not intended as a replacement for SECISearch. Our patterns and free-energy cutoffs are not stringent and will result in a high false positive rate if used to identify novel SECIS elements. Ideally, SECISaln should be used on sequences which are known to contain a SECIS element, and its main application is the detailed characterization of structural features in the identified SECIS elements, through the multiple structural comparison to other known SECIS elements.
In addition to being the first structural alignment tool for SECIS elements, SECISaln also provides the largest available, manually curated collection of eukaryotic SECISes. Our SECIS collection was built by searching for homologs of all known eukaryotic selenoproteins in NCBIs Refseq mRNA and TIGRs EGO databases. We ran TBLASTN searches using the human (when available, other species when not) selenoproteins as queries. We then extracted the relevant mRNA sequence from the database and identified its SECIS element. We also manually added insect SECIS sequences that had been previously identified (Chapple and Guigó, 2008), but which are not yet present in mRNA databases. This process resulted in a collection of 62 type I and 224 type II SECISes, a clear indication that type II constitute the major part of SECIS elements. Interestingly, although all selenoprotein families had a type II SECIS in at least one species, SelO, SelT, MsrA, DI2, SelS, 15kDa, TR3, SelI, Gpx3 and TR2 had type II SECISes in all species investigated. GPx1 and DI1 had type I SECISes in all species except Danio rerio.
Analyzing the structural alignments produced by SECISaln provided a more detailed picture of SECIS structural features. For instance the length of helix 2, which was previously set to 14 bp, is less constrained and ranges in fact from 11 bp to 16 bp. SECISaln also highlighted previously unknown conserved residues in eukaryotic SECIS elements (see Supplementary Table 1), which can be summarized as a new consensus core sequence for eukaryotic SECIS elements as shown in Figure 1. Most striking of these is an overrepresentation of G at position 1 (3 to abcd) and a corresponding overrepresentation of Y (C or U) at position 1. We also observed a clear overrepresentation of U in type I elements, and Y in type II at position –4. This is particularly surprising since no cross-species sequence conservation has ever been observed five to the quartet, with the exception of the conserved R, and may be connected to the SBP2-SECIS contacts observed in this area (Cléry et al., 2007; Fletcher et al., 2001)
In conclusion, we believe that SECISaln, as has already been demonstrated by the analyses presented here, will be a very useful tool for the analysis and understanding of SECIS elements.
|
ACKNOWLEDGEMENTS |
|---|
|
TOPABSTRACTACKNOWLEDGEMENTSREFERENCES |
|---|
The authors would like to thank David Martin for his help with CGI scripting and Marco Mariotti for beta testing the software.
Funding: Spanish Ministry of Education and Science (to R.G.); BioSapiens European Network of Excellence (to R.G.); National Institute for Bioinformatics (www.inab.org) a platform of ‘Genoma Espa na’ (to R.G.); ToxNuc-E program (to A.K.); ACI BCMS of the French Ministry of Research (to A.K.); Pre-doctoral Fellowship from the Spanish Ministry of Education and Science (to C.E.C.).
Conflict of Interest: none declared.
|
FOOTNOTES |
|---|
Associate Editor: Ivo Hafacker
1 With one exception, the SelT genes of Toxoplasma gondii and Neospora canine have a non-canonical GGA/GA sequence instead (Novoselov et al., 2007). 
Received on October 7, 2008; revised on December 18, 2008; accepted on January 7, 2009
|
REFERENCES |
|---|
|
TOPABSTRACTACKNOWLEDGEMENTSREFERENCES |
|---|
Buettner C, et al. The Caenorhabditis elegans homologue of thioredoxin reductase contains a selenocysteine insertion sequence (secis) element that differs from mammalian secis elements but directs selenocysteine incorporation. J. Biol. Chem. (1999) 274:21598–21602.
Castellano S, et al. In silico identification of novel selenoproteins in the Drosophila melanogaster genome. EMBO Rep. (2001) 2:697–702.[CrossRef][Web of Science][Medline]
Chapple CE, Guigó R. Relaxation of selective constraints causes independent selenoprotein extinction in insect genomes. PLoS ONE (2008) 3:e2968.[CrossRef][Medline]
Cléry A, et al. An improved definition of the rna-binding specificity of secis-binding protein 2, an essential component of the selenocysteine incorporation machinery. Nucleic Acids Res. (2007) 35:1868–1884.
Fagegaltier D, et al. Structural analysis of new local features in secis RNA hairpins. Nucleic Acids Res. (2000) 28:2679–2689.
Fletcher JE, et al. The selenocysteine incorporation machinery: interactions between the secis RNA and the secis-binding protein sbp2. RNA (2001) 7:1442–1453.[Abstract]
Grundner-Culemann E, et al. Two distinct secis structures capable of directing selenocysteine incorporation in eukaryotes. RNA (1999) 5:625–635.[Abstract]
Kryukov GV, et al. New mammalian selenocysteine-containing proteins identified with an algorithm that searches for selenocysteine insertion sequence elements. J. Biol. Chem. (1999) 274:33888–33897.
Kryukov GV, et al. Characterization of mammalian selenoproteomes. Science (2003) 300:1439–1443.
Lescure A, et al. Novel selenoproteins identified in silico and in vivo by using a conserved rna structural motif. J. Biol. Chem. (1999) 274:38147–38154.
Lobanov AV, et al. Evolutionary dynamics of eukaryotic selenoproteomes: large selenoproteomes may associate with aquatic life and small with terrestrial life. Genome Biol. (2007) 8:R198.[CrossRef][Medline]
Novoselov SV, et al. A highly efficient form of the selenocysteine insertion sequence element in protozoan parasites and its use in mammalian cells. Proc. Natl Acad. Sci. USA (2007) 104:7857–7862.
Walczak R, et al. A novel RNA structural motif in the selenocysteine insertion element of eukaryotic selenoprotein mRNAs. RNA (1996) 2:367–379.[Abstract]
CiteULike 
Connotea 
Del.icio.us What's this?
This article has been cited by other articles:
|
|
 |
|
  L. Latreche, O. Jean-Jean, D. M. Driscoll, and L. Chavatte Novel structural determinants in human SECIS elements modulate the translational recoding of UGA as selenocysteine Nucleic Acids Res., September 1, 2009; 37(17): 5868 - 5880. [Abstract] [Full Text] [PDF]
|
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||