Splicing regulatory elements (SREs) are short, degenerate sequences on pre-mRNA molecules that enhance or inhibit the splicing process via the binding of splicing factors, proteins that regulate the functioning of the spliceosome. Existing methods for identifying SREs in a genome are either experimental or computational. Here, we propose a formalism based on de Bruijn graphs that combines genomic structure, word count enrichment analysis, and experimental evidence to identify SREs found in exons. In our approach, SREs are not restricted to a fixed length (i.e., k-mers, for a fixed k). As a result, we identify 2001 putative exonic enhancers and 3080 putative exonic silencers for human genes, with lengths varying from 6 to 15 nucleotides. Many of the predicted SREs overlap with experimentally verified binding sites. Our model provides a novel method to predict variable length putative regulatory elements computationally for further experimental investigation.
Get full access to this article
View all access options for this article.
References
1.
AlexandreP., GalanteF., SakabeN.J., et al. 2004. Detection and evaluation of intron retention events in the human transcriptome. RNA, 10, 757–765.
2.
BarashY., BlencoweB.J., and FreyB.J.2010a. Model-based detection of alternative splicing signals. Bioinformatics, 26, i325–i333.
3.
BarashY., CalarcoJ.a., GaoW., et al. 2010b. Deciphering the splicing code. Nature, 465, 53–59.
4.
BarisO., DelettreC., Amati-BonneauP., et al. 2003. Fourteen novel OPA1 mutations in autosomal dominant optic atrophy including two de novo mutations in sporadic optic atrophy. Hum. Mutat., 21, 656–656.
5.
BatageljV., and ZaversnikM.2003. An O(m) algorithm for cores decomposition of networks. arXiv Prepr. cs/0310049, 1–10.
6.
BauerS., GrossmannS., VingronM., et al. 2008. Ontologizer 2.0–a multifunctional tool for GO term enrichment analysis and data exploration. Bioinformatics, 24, 1650–1651.
7.
BellareP., SmallE.C., HuangX., et al. 2008. A role for ubiquitin in the spliceosome assembly pathway. Nat. Struct. Mol. Biol., 15, 444–451.
8.
BenjaminiY., and HochbergY.1995. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B, 57, 289–300.
9.
BuendiaP., TyreeJ., LoredoR., et al. 2012. Identification of conserved splicing motifs in mutually exclusive exons of 15 insect species. BMC Genomics, 13, S1.
10.
BurnettB.G., MunozE., TandonA., et al. 2008. Regulation of SMN protein stability. Mol. Cell. Biol., 29, 1107–1115.
DjordjevicM.2007. SELEX experiments: New prospects, applications and data analysis in inferring regulatory pathways. Biomol. Eng., 24, 179–189.
13.
EZ., WangL., and ZhouJ.2013. Splicing and alternative splicing in rice and humans. BMB Rep., 46, 439–447.
14.
EichnerJ., ZellerG., LaubingerS., et al. 2011. Support vector machines-based identification of alternative splicing in Arabidopsis thaliana from whole-genome tiling arrays. BMC Boinformatics, 12, 55–55.
15.
EngleE.C.2010. Human genetic disorders of axon guidance. Cold Spring Harb. Perspect. Biol., 2, a001784.
16.
FairbrotherW.G., YehR.F., SharpP.A., et al. 2002. Predictive identification of exonic splicing enhancers in human genes. Science, 297, 1007–1013.
17.
FedorovA., SaxonovS., FedorovaL., et al. 2001. Comparison of intron-containing and intron-lacking human genes elucidates putative exonic splicing enhancers. Nucleic Acids Res., 29, 1464–1469.
18.
FerreiraE.N., GalanteP.A.F., CarraroD.M., et al. 2007. Alternative splicing: A bioinformatics perspective. Mol. Biosyst., 3, 473–477.
19.
Garcia-BlancoM.A., BaraniakA.P., and LasdaE.L.2004. Alternative splicing in disease and therapy. Nat. Biotechnol., 22, 535–546.
20.
GiuliettiM., PivaF., D'AntonioM., et al. 2013. SpliceAid-F: a database of human splicing factors and their RNA-binding sites. Nucleic Acids Res., 41, D125–D131.
21.
GorenA., RamO., AmitM., et al. 2006. Comparative analysis identifies exonic splicing regulatory sequences–The complex definition of enhancers and silencers. Mol. Cell, 22, 769–781.
22.
GuilS., and CáceresJ.F.2007. The multifunctional RNA-binding protein hnRNP A1 is required for processing of miR-18a. Nat. Struct. Mol. Biol., 14, 591–596.
HopcroftJ.E., and UllmanJ.D.1979. Introduction to Automata Theory, Languages, and Computation. Addison-Wesley Publishing Company, Reading, Massachusetts.
25.
KarolchikD., HinrichsA.S., FureyT.S., et al. 2004. The UCSC table browser data retrieval tool. Nucleic Acids Res., 32, D493–D496.
26.
KeS., and ChasinL.A.2010. Intronic motif pairs cooperate across exons to promote pre-mRNA splicing. Genome Biol., 11, R84.
27.
KeS., ShangS., KalachikovS.M., et al. 2011. Quantitative evaluation of all hexamers as exonic splicing elements. Genome Res., 21, 1360–1374.
28.
KerenH., Lev-MaorG., and AstG.2010. Alternative splicing and evolution: Diversification, exon definition and function. Nat. Rev. Genet., 11, 345–355.
29.
KimJ., ZhaoS., HowardB.E., et al. 2009. Mining of cis-regulatory motifs associated with tissue-specific alternative splicing. Lecture Notes in Computer Science., 5542, 260–271.
30.
KitsakM., GallosL.K., HavlinS., et al. 2010. Identification of influential spreaders in complex networks. Nat. Phys., 6, 888–893.
31.
LvY., ZuoZ., and XuX.2013. Global detection and identification of developmental stage specific transcripts in mouse brain using subtractive cross-screening algorithm. Genomics, 102, 229–236.
32.
MatlinA.J., ClarkF., and SmithC.W.J.2005. Understanding alternative splicing: towards a cellular code. Nat. Rev. Mol. Cell Biol., 6, 386–398.
33.
MillevoiS., BernatS., TellyD., et al. 2010. The c.5242C>A BRCA1 missense variant induces exon skipping by increasing splicing repressors binding. Breast Cancer Res. Treat., 120, 391–399.
34.
NugentA.a., KolpakA.L., and EngleE.C.2012. Human disorders of axon guidance. Curr. Opin. Neurobiol., 22, 837–843.
35.
PemmarajuS., and SkienaS.2003. Computational discrete mathematics: combinatorics and graph theory with mathematica. Cambridge University Press, Cambridge.
36.
PerteaM., MountS.M., and SalzbergS.L.2007. A computational survey of candidate exonic splicing enhancer motifs in the model plant Arabidopsis thaliana. BMC Bioinformatics, 8, 159–159.
37.
RamalhoR.F., GelfmanS., De SouzaJ.E., et al. 2013. Testing for natural selection in human exonic splicing regulators associated with evolutionary rate shifts. J. Mol. Evol., 76, 228–239.
38.
RosenbergA.L., and HeathL.S.2000. Graph Separators, With Applications. Kluwer Academic/Plenum Publishers.
39.
SakabeN.J., and De SouzaS.J.2007. Sequence features responsible for intron retention in human. BMC Genomics, 8, 59–59.
40.
SanfordJ.R., CoutinhoP., HackettJ.a., et al. 2008. Identification of nuclear and cytoplasmic mRNA targets for the shuttling protein SF2/ASF. PLoS One, 3, e3369.
41.
SelvakumarM., and HelfmanD.M.1999. Exonic splicing enhancers contribute to the use of both 3′ and 5′ splice site usage of rat beta-tropomyosin pre-mRNA. RNA, 5, 378–394.
42.
StammS., RiethovenJ.-J., Le TexierV., et al. 2006. ASD: A bioinformatics resource on alternative splicing. Nucleic Acids Res., 34, D46–D55.
43.
SzcześniakM.W., KabzaM., PokrzywaR., et al. 2013. ERISdb: a database of plant splice sites and splicing signals. Plant Cell Physiol., 54, e10.
44.
TackeR., and ManleyJ.L.1995. The human splicing factors ASF/SF2 and SC35 possess distinct, functionally significant RNA binding specificities. EMBO J., 14, 3540–3551.
45.
TollerveyJ.R., CurkT., RogeljB., et al. 2011. Characterizing the RNA targets and position-dependent splicing regulation by TDP-43. Nat. Neurosci., 14, 452–458.
46.
UleJ., JensenK.B., RuggiuM., et al. 2003. CLIP identifies Nova-regulated RNA networks in the brain. Science, 302, 1212–1215.
47.
WangZ., and BurgeC.B.2008. Splicing regulation: From a parts list of regulatory elements to an integrated splicing code. RNA, 14, 802–813.
48.
WangZ., RolishM.E., YeoG., et al. 2004. Systematic identification and analysis of exonic splicing silencers. Cell, 119, 831–845.
49.
WarnerL.E., and ChamberlainJ.S.2006. Minigenes. Encycl. Life Sci., 1–6.
WenJ., ChenZ., and CaiX.2013. A biophysical model for identifying splicing regulatory elements and their interactions. PLoS One, 8, e54885.
52.
WenJ., ChibaA., and CaiX.2010. Computational identification of tissue-specific alternative splicing elements in mouse genes from RNA-Seq. Nucleic Acids Res., 38, 7895–7907.
53.
ZhangJ., KuoC.C.J., and ChenL.2011. GC content around splice sites affects splicing through pre-mRNA secondary structures. BMC Genomics, 12, 90–90.
54.
ZhangJ., KuoC.-C.J., and ChenL.2012. VERSE: a varying effect regression for splicing elements discovery. J. Comput. Biol., 19, 855–865.
55.
ZhangX.H., KangsamaksinT., MannS.P., et al. 2005. Exon inclusion is dependent on predictable exonic splicing enhancers. Mol. Cell. Biol., 25, 7323–7332.
56.
ZhangX.H.F., and ChasinL.A.2004. Computational definition of sequence motifs governing constitutive exon splicing. Genes Dev., 18, 1241–1250.
57.
ZhangX.H.F., HellerK.A., HefterI., et al. 2003. Sequence information for the splicing of human pre-mRNA identified by support vector machine classification. Genome Res., 13, 2637–2650.
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
