Volume 9 Supplement 12
Seventh International Conference on Bioinformatics (InCoB2008)
Fast splice site detection using information content and feature reduction
- AKMA Baten^{1}Email author,
- SK Halgamuge^{1} and
- BCH Chang^{2}
DOI: 10.1186/1471-2105-9-S12-S8
© Baten et al; licensee BioMed Central Ltd. 2008
Published: 12 December 2008
Abstract
Background
Accurate identification of splice sites in DNA sequences plays a key role in the prediction of gene structure in eukaryotes. Already many computational methods have been proposed for the detection of splice sites and some of them showed high prediction accuracy. However, most of these methods are limited in terms of their long computation time when applied to whole genome sequence data.
Results
In this paper we propose a hybrid algorithm which combines several effective and informative input features with the state of the art support vector machine (SVM). To obtain the input features we employ information content method based on Shannon's information theory, Shapiro's score scheme, and Markovian probabilities. We also use a feature elimination scheme to reduce the less informative features from the input data.
Conclusion
In this study we propose a new feature based splice site detection method that shows improved acceptor and donor splice site detection in DNA sequences when the performance is compared with various state of the art and well known methods.
Background
Over the past decades, the scientific community has experienced a major growth in numbers of sequence data. With the emergence of novel and efficient sequencing technology, DNA sequencing is now much faster. Sequencing of several genomes including the human genome have been completed successfully. This massive amount of sequence data demands sophisticated tools for the analysis of data.
Identifying genes accurately is one of the most important and challenging tasks in bioinformatics and it requires the prediction of the complete gene structure. Identification of splice sites is the core component of eukaryotic gene finding algorithms. Their success depends on the precise identification of the exon-intron structure and the splice sites. Most of the eukaryotic protein coding genes are characterized by exons and introns. Exons are the protein coding portion of a gene and they are segmented with intervening sequences of introns. The border between an exon and an intron is known as the splice site. The splice site upstream of an intron is called the donor splice site (in the direction 5' to 3') and one that is downstream of an intron is the acceptor splice site (in the direction 3' to 5'). The consensus sequence refers to the nucleotides, which are conserved or most frequently observed in a particular position. The acceptor and donor splice sites with consensus AG (corresponding to the end of an intron) and GT (corresponding to the beginning of an intron) dinucleotides respectively are known as canonical splice sites. Approximately 99% of the splice sites are canonical [1]. As AG and GT represent possible acceptor and donor splice sites, every AG and GT in a DNA sequence is a candidate acceptor or donor splice site and they need to be classified as either a real (true) splice site or a pseudo (false) splice site.
Over the years many computational methods have been proposed for the identification of splice sites. Most of those methods are designed to identify the apparent consensus AG and GT in the splicing junction. These methods can be largely classified into probabilistic methods [2–8], neural network and support vector machine methods [9–19], and methods based on discriminant analysis [20, 21]. Neural networks and support vector machines (SVM) learn the complex features of neighbourhoods surrounding the consensus AG/GT dinucleotides by a complex non-linear transformation. Probabilistic models estimate position specific probabilities of splice sites by computing likelihoods of candidate signal sequences. The discriminant analysis uses several statistical measures to evaluate the presence of specific nucleotides, recognizing the splice sites without explicitly determining the probability distributions [18].
In DNA sequences, true consensus AG/GT dinucleotides are outnumbered by many false AG/GTs. However, nucleotides surrounding true AG/GTs show a certain nucleotide dependency and sequential relationship compared to those surrounding false AG/GTs. There are several methods which are particularly designed to capture this relationship and to identify true splice sites among numerous false ones. Weight matrix methods (WMM) and methods based on Markov models are popular methods of this category. WMM was successfully adopted in methods like NetPlantGene [22] and NNSplice [10]. Salzberg et al. and Zhang et al. [2, 6], used a linear first order Markov model (MM1) also known as the weight array method (WAM) and they have achieved a good splice site prediction accuracy. MM1 only utilizes first order sequential relationship. It is desirable to use a higher order Markov model to capture the higher order and extended sequential relationship. However, the computational complexity increases polynomialy with the increase of the order of the Markov model, and also higher order Markov models require a large number of training samples. The maximal dependence decomposition (MDD) algorithm was proposed by Burge et al. [23] to overcome these limitations. MDD is a decision tree process and models the dependency between adjacent nucleotides. To take the advantages of both MDD and Markov models, Pertea et al. [4] proposed the GeneSplicer method which combines MDD and second order Markov models (MM2). GeneSplicer showed an improved splice site detection performance. More recently, Rajapakse et al. [17] proposed a complex splice site detection method by combining mostly second order Markov models with backpropagation neural networks (BPNN). This method showed an improved performance over GeneSplicer, however, BPNN is already computationally expensive and this method requires a larger sequence window. In contrast, a machine learning technique such as SVM has the advantage of inferring an optimal classifier from the training data. SVM has been used to classify splice site data with limited success [9, 12, 14–16].
Most of the existing splice site detection methods focused on the improvement of classification performance. However other studies suggest that, considering the increasing growth of sequence data, the focus of new methods should be towards developing faster methods [24–27]. In our previous work we showed an improved splice site detection performance by using several preprocessing methods including WMM0/MM0, WMM1, MM1 with SVM [18]. However, the training time and the number of input features to SVM is a major concern. SVM performs better when it is trained with the most important and meaningful features. So, the reduction of less important features may improve both the classification performance and training time of SVM. In this paper, we propose a feature selection strategy which reduces the less important features from the input data. We also combine the well studied information content method based on Shannon's information theory [28–30] and the Shapiro's score method [31] to extract meaningful information from sequence that can potentially identify splice sites. Our method showed an improved splice site detection performance when compared to the existing methods in terms of classification accuracy and training time.
Results
Classification performance comparison
Our hybrid algorithm combines several effective and informative input features with the state of the art support vector machine (SVM). To obtain the input features we employ information content method based on Shannon's information theory, Shapiro's score scheme, and Markovian probabilities. We also use the F-score feature elimination scheme to reduce the less informative features from the input data. We use the publicly available NN269 [10] splice site dataset to evaluate the performance of our method. The MM1 parameters are calculated from the dataset and F-score method (refer to the method section) is applied to reduce the number of MM1 parameters, which is referred as Reduced MM1 SVM method. We also calculate the information content and Shapiro's score from the dataset and use the proposed the IC Shapiro SVM method, which is a linear combination of information content and Shapiro's score. We compare the performance of our methods with MM1 SVM method as proposed [18]. To evaluate the classification performance we use several performance evaluation methods such as the sensitivity, specificity, receiver operating characteristics curve (ROC), and the area under ROC (AUC) as described in the method section.
AUC and training time for different models for NN269 acceptor splice sites.
Model | SVM kernel | AUC | Training time until convergence (hh.mm.ss) |
---|---|---|---|
Reduced MM1 SVM (Best in terms of accuracy) | GRBF | 0.9741389 | 00.22.17 |
Reduced MM1 SVM | Polynomial | 0.9695822 | 00.10.48 |
MM1 SVM [18] | Polynomial | 0.9674048 | 00.11.02 |
IC Shapiro SVM (Best In terms of Time) | Polynomial | 0.96287 | 00:01:18 |
AUC and training time improvement for different models compared to MM1-SVM method for NN269 acceptor splice sites.
Model | SVM kernel | AUC | Training time until convergence (mm.ss) | Performance Improvement | Time Improvement |
---|---|---|---|---|---|
MM1 SVM [18] | Polynomial | 0.9674 | 11.02 | - | - |
Reduced MM1 SVM (Best in terms of accuracy) | GRBF | 0.9741 | 22.17 | 0.69% | -101.96% |
Reduced MM1 SVM | Polynomial | 0.9695 | 10.48 | 0.2171% | 2.11% |
IC Shapiro SVM (Best In terms of Time) | Polynomial | 0.9628 | 01:18 | -0.4755% | 88.21% |
As shown in Table 2, the best acceptor splice site detection performance is produced by Reduced MM1 SVM with GRBF kernel which is 0.69% superior then MM1 SVM with polynomial kernel. However, Reduced MM1 SVM GRBF requires much longer training time (more than 100%) than MM1 SVM Polynomial. Reduced MM1 SVM Polynomial improves the performance by 0.21% and it also 2.11% faster than MM1 SVM Polynomial. Finally, IC Shapiro SVM Polynomial is just 0.47% worse then MM1 SVM Polynomial, however, it shows a significant improvement in the training time and is 88.21% faster than MM1 SVM Polynomial.
AUC and training time for different models for NN269 donor splice sites.
Model | SVM kernel | AUC | Training time until convergence (hh.mm.ss) |
---|---|---|---|
Reduced MM1 SVM (Best in terms of accuracy) | GRBF | 0.9790232 | 00:20:04 |
Reduced MM1 SVM | Polynomial | 0.9764903 | 00:09:30 |
MM1 SVM | Polynomial | 0.9761952 | 00:10:02 |
IC Shapiro SVM (Best In terms of Time) | Polynomial | 0.9665982 | 00:02:59 |
AUC and training time improvement for different models compared to MM1-SVM method for NN269 donor splice sites.
Model | SVM kernel | AUC | Training time until convergence (mm.ss) | Performance Improvement | Time Improvement |
---|---|---|---|---|---|
MM1 SVM | Polynomial | 0.9761 | 10:02 | - | - |
Reduced MM1 SVM (Best in terms of accuracy) | GRBF | 0.9790 | 20:04 | 0.297% | -100% |
Reduced MM1 SVM | Polynomial | 0.9764 | 09:30 | 0.0102% | 5.31% |
IC Shapiro SVM (Best In terms of Time) | Polynomial | 0.9665 | 02:59 | -0.9835% | 70.26% |
As shown in Table 4, Reduced MM1 SVM GRBF marginally produces the best acceptor splice site detection performance which is 0.29% superior then MM1 SVM with polynomial kernel. However, it requires much longer training time (more than 100%) than MM1 SVM Polynomial. Reduced MM1 SVM Polynomial performs almost equally as well as MM1 SVM Polynomial, though it is 5.31% faster. Finally, IC Shapiro SVM Polynomial is almost 1% worse then MM1 SVM Polynomial, however, it shows a significant improvement in the training time and is 70.26% faster than MM1 SVM Polynomial. All the parameters regarding the SVM implementations with GRBF and Polynomial kernels are provided in Additional file 1: Table S1.
Discussions
One of the biological machineries involved in the splicing process is known as the Spliceosome, which binds in a splice site after determining the information available in that site. Information content and Shapiro's score are two well known methods to determine the information in splice sites [29, 30]. Previously we have used other methods such as first order Markov model (MM1), first order weight matrix model (WMM1) and zero order Markov/Weight matrix model (MM0/WMM0) to capture such information [18]. However, our results in this paper show that information content and Shapiro's score are more capable of capturing more meaningful information than those methods we have used previously, which justifies the use of these methods to extract features. Our method based on information content and Shapiro's score is also proved to be much faster.
We also use the F-score feature ranking measure to select the most meaningful features and use it to eliminate less important features. We use the F-score measure to reduce MM1 parameters and when compared with MM1 SVM method [18], we find that Reduced MM1 SVM method based on reduced MM1 parameters performs better in terms of classification accuracy as shown in Figures 1, 2, 3, and 4. The performance of IC Shapiro SVM is marginally worse then the best performing methods. However, it shows much faster training time than others as IC Shapiro SVM uses much less number of features. This IC and Shapiro's score scheme and their integration as a set of features is an important step towards faster splice site identification and can be effectively used in the splice site detection for the whole genome where vast amount of sequences data is available. However, it is worthwhile to further investigate this method to improve its classification accuracy.
Conclusion
Modern sequencing techniques can produce a massive amount of data in short time and the number of sequence data is almost exponentially increasing. Fast splice site detection is very useful when we consider the very large volumes of available data for the training and testing of a method. To cope with such a large volume of data we also need faster methods. The fast splice site detection method we propose in this paper can also be applied to identify other signals in the sequence such as promoters and translation initiation sites.
Methods
Proposed models
We propose several models for the identification of acceptor and donor splice sites. Corresponding to the two types of splice sites, the splice site classification problem is subdivided into two classification problems: acceptor splice site classification and donor splice site classification. Two separate models are constructed for the identification of acceptor splice sites and donor splice sites respectively.
Proposed models and their description.
Model | Description |
---|---|
Reduced MM1 SVM Polynomial | Only reduced MM1 parameters and SVM with polynomial kernel |
Reduced MM1 SVM GRBF | Only reduced MM1 parameters and SVM with GRBF kernel |
IC Shapiro SVM Polynomial | Information content, Shapiro's score and SVM with polynomial kernel |
Markov model
Markov model can be regarded as a finite state machine with Markov property. Let us consider a sequence of random variables X_{1}, X_{2},..... X_{ n }which takes on values from a finite state space A = {A_{1}, A_{2},... A_{ n }}. If the probability of transition from state A_{ i }at time n to state A_{ j }at time n +1 depends only on A_{ i }, and not any previous history of process, then the process is said to have the Markov property or to be a Markov model or Markov chain.
Where, the Markovian probability P_{ i }(s_{ i }) = P(s_{ i }|s_{i-1}, s_{i-2},...., s_{i-k}) denotes the conditional probability of a nucleotide at location i given the k predecessors. In the current application we use a first order Markov model to model the sequences and hence, k = 1.
Information content
where, f(b, l) is the probability of base b at position l.
We first generated an individual information weight matrix from the frequencies of each nucleotide at each position to calculate the information content (R_{ i }, bits) of each splice site sequence. The individual information weight matrix can be calculated by the following equation [29]:
R_{ iw }(b, l) = 2 + log_{2} f(b, l) (3)
The information content of each splice site was calculated by summing up R_{ iw }(b, l) at each position of the splice site sequences. The relationship between R_{ iw }(b, l) and R_{ sequence }(l) is provided in Additional file 2.
Shapiro's score
Shapiro et al. [31] proposed a method to score the strength of splice sites based on percentage of each nucleotide at each position. First they create a frequency matrix of nucleotides in each of the positions of the splice site sequence. Shapiro's score for acceptor splice site is given by the equation [31]:
SS_{ acceptor }= 100*((t 1-l 1)/(h 1-l 1) + (t 2-l 2)/(h 2-l 2))/2 (4)
where, t 1 is the sum of best 8 of 10 nucleotide percentages at position -13 to -4
l 1 is the sum of lowest 8 of 10 nucleotide percentages at position -13 to -4
h 1 is the sum of highest 8 of 10 percentages at position -13 to -4
t 2 is the sum of best nucleotide percentages at position -3 to + 1
l 2 is the sum of lowest nucleotide percentages at position -3 to + 1
h 2 is the sum of highest nucleotide percentages at position -3 to +1
Similarly, Shapiro's score for donor splice site is given by the equation [31]:
SS_{ donor }= 100*(t - min)/(max - min) (5)
where, t is the sum of percentages at position -3 to + 7
min is the sum of lowest percentages at position -3 to + 7
max is the sum of highest percentages at position -3 to + 7
Sequence feature elimination based on F-score
Sequence feature elimination is an important step towards the classification task. Classifiers like neural networks, SVM's etc. perform better when they are trained with meaningful input data. Redundant data often causes misclassification and hence, the reduction of classification performance. So it is desirable to eliminate the less important features from the input data and to select those features that can potentially discriminate between true and false class. According to Dror et al. [32], there are three potential benefits of feature selection: improving the performance of the classifier, producing a cost-effective classifier, and providing a better understanding of the problem.
Support vector machine
0 ≤ α_{ i }≤ C, i = 1,..., l, (9)
where v is the set of support vectors.
Dataset
To evaluate the performance of the proposed models, we used publicly available NN269 [10] splice site dataset. The dataset is divided into two groups namely: the acceptor splice sites and the donor splice sites. It contains 1324 confirmed true acceptor splice sites, 5552 false acceptor sites, 1324 confirmed true donor sites, and 4922 false donor sites collected from 269 human genes. The pseudo or false acceptor/donor splice sites are those having AG/GT in the splicing junction but not a real acceptor or donor splice site according to the annotation. Acceptor splice sites have a window of 90 nucleotides (-70 to +20) with the consensus nucleotides AG at positions -69 and -70. This window includes the last 70 nucleotides of an intron and the first 20 nucleotides of the succeeding exon. On the other hand, donor splice sites have a window of 15 nucleotides (-7 to +8) with the consensus nucleotides GT at positions +1 and +2. This window includes the last 9 nucleotides of an exon and the first 6 nucleotides of the succeeding intron. The acceptor and donor splice site datasets are divided into a unique training and test dataset. The test datasets do not contain any sequence which is in training dataset. The training dataset contains 1116 true acceptor, 1116 true donor, 4672 false acceptor, and 4140 false donor sites. The test data set contains 208 true acceptor sites, 208 true donor sites, 881 false acceptor sites, and 782 false donor sites.
Model learning
The learning of the model is designed in two phases. Phase one consists of estimation of Markov parameters, scoring information content of sequences, and calculation of Shapiro's score. In phase two SVM is trained with polynomial and Gaussian kernels.
A position specific nucleotide background matrix is required for the calculation of information content and Shapiro's scores. A generalized frequency matrix of the splice site regions of the whole genome is preferable as it gives the most reasonable statistics of the occurrence of nucleotides in the splice site. However, we only use the NN269 true training data to construct the frequency matrix. As we compared the performance of MM1 SVM [18] with that of IC Shapiro SVM, it is required that the same training data be used to create MM1 parameters, information content and Shapiro's scores. To calculate the information content score the individual information content weight matrix R_{ iw }(b, l) is created from the nucleotide background matrix following equation (2) (refer to the method section). Then the information content is calculated by summing up R_{ iw }(b, l) of the specified positions. Similarly the nucleotide background matrix is used to calculate the Shapiro's score for acceptor and donor splice sites following equations (3) and (4) respectively.
We used the leave one out cross validation technique is applied to determine the splice site prediction accuracy and to compare the predictive accuracy with other methods. The cross validation is performed by randomly partitioning the data into five independent subsets. Each of the subsets does not share any repeating sequences. Each model was trained by selecting four of the subsets (training data) and was tested on the remaining one. Finally, we took the average of the five prediction accuracies as the final prediction measure of the model.
Performance measures
Definition of TP, TN, FP and FN
Predicted positive | Predicted negative | |
---|---|---|
Real positive | true positives, TP | false negatives, FN |
Real negative | true negatives, TN | false positives, FP |
Also receiver operator curve (ROC) is drawn using the sensitivity and specificity values. ROC analysis is an effective and widely used method for assessing the classification performance [40]. When a ROC is created from the sensitivity (the y axis) and specificity (the x axis) of a model, the closer a curve follows the left-hand border and then the top of the border of the ROC plot, the more accurate the model is (refer to Figure 3, 4, 5 and 6). We also calculate the area under ROC curve (AUC), as classification performance of some of the models are very close and may not clearly distinguish performance of two models when we view them in the ROC curve. However, AUC accurately measure the total ROC area covered by a model.
Declarations
Acknowledgements
We gratefully acknowledge the helpful discussions and comments provided by Dr. Arthur Hsu.
This article has been published as part of BMC Bioinformatics Volume 9 Supplement 12, 2008: Asia Pacific Bioinformatics Network (APBioNet) Seventh International Conference on Bioinformatics (InCoB2008). The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2105/9?issue=S12.
Authors’ Affiliations
References
- Burset M, Seledtsov A, Solovyeva VV: Analysis of canonical and non-canonical splice sites in mammalian genomes. Nucleic Acids Research. 2000, 28 (21): 4364-4375.PubMed CentralView ArticlePubMedGoogle Scholar
- Salzberg S: A method for identifying splice sites and translation start site in eucaryotic mRNA. Computer Applications in the Biosciences. 1997, 13 (4): 384-390.Google Scholar
- Chen T-M, Chung-Chin Lu, Wen-Hsiung Li: Prediction of splice sites with dependency graphs and their expanded bayesian networks. Bioinformatics. 2005, 21 (4): 471-482.View ArticlePubMedGoogle Scholar
- Pertea M, Xiao Ying L, Salzberg SL: GeneSplicer: a new computational method for splice site detection. Nucleic Acids Research. 2001, 29 (5): 1185-1190.PubMed CentralView ArticlePubMedGoogle Scholar
- Marashi SA, Changiz Eslahchi, Pezeshk H, Sadeghi M: Impact of RNA structure on the prediction of donor and acceptor splice sites. BMC Bioinformatics. 2006, 7: 297-PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang M, Marr T: A weight array method for splicing signal analysis. Comput Appl Biosci. 1993, 9: 499-509.PubMedGoogle Scholar
- Castelo R, Guigo R: Splice site identification by idlBNs. Bioinformatics. 2004, 20 (1): 69-76.View ArticleGoogle Scholar
- Cai D, Delcher A, Kao B, Kasif S: Modeling splice sites with Bayes networks. Bioinformatics. 2000, 16 (2): 152-158.View ArticlePubMedGoogle Scholar
- Ratsch G, Sonnenburg S, Schafer C: Learning Interpretable SVMs for Biological Sequence Classification. BMC Bioinformatics. 2006, 7 (suppl 1): S9-PubMed CentralView ArticlePubMedGoogle Scholar
- Reese MG, Eeckman F, Kupl D, Haussler D: Improved splice site detection in Genie. Journal of Computational Biology. 1997, 4 (3): 311-324.View ArticlePubMedGoogle Scholar
- Brunak S, Engelbrecht J, Knudsen S: Prediction of mRNA donor and acceptor sites from the DNA sequence. Journal of Molecular Biology. 1991, 220: 49-65.View ArticlePubMedGoogle Scholar
- Sun YF, Fan XD, Li YD: Identifying splicing sites in eukaryotic RNA: Support vector machine approach. Computers in biology and medicine. 2003, 33: 17-29.View ArticlePubMedGoogle Scholar
- Reese MG: Application of a time-delay neural network to promoter annotation in the Drosophila melanogaster. Computer chem. 2001, 26 (1): 51-56.View ArticleGoogle Scholar
- Zhang XH-F, Katherine AH, Ilana H, Christina SL, Lawrence AC: Sequence information for the splicing of human pre-mRNA identified by support vector machine classification. Genome Research. 2003, 13: 2637-2650.PubMed CentralView ArticlePubMedGoogle Scholar
- Sonnenburg S: New methods for detecting splice junction sites in DNA sequence. Master's Thesis. 2002, Humbold University, GermanyGoogle Scholar
- Degroeve S, Saeys Y, Baets BD, Rouze P, Peer YVD: SpliceMachine: predicting splice sites from high-dimensional local context representations. Bioinformatics. 2005, 21 (8): 1332-1338.View ArticlePubMedGoogle Scholar
- Rajapakse JCaHLS: Markov encoding for detecting signals in genomic sequences. IEEE/ACM Transactions on Computational Biology and Bioinformatics. 2005, 2 (2): 131-142.View ArticlePubMedGoogle Scholar
- Baten AKMA, Chang BCH, Halgamuge SK, Li J: Splice site identification using probabilistic parameters and SVM classification. BMC Bioinformatics. 2006, 7 (Suppl 5):Google Scholar
- Sonnenburg S, Schweikert G, Philips P, Behr J, Rätsch G: Accurate splice site prediction using support vector machines. BMC Bioinformatics. 2007, 8 (10):Google Scholar
- Chuang JSaRD: Splice site prediction using a sparse network of winnows. Technical Report. 2001, University of Illinois, Urbana-ChampaignGoogle Scholar
- Zhang L: aLL: Splice site prediction with quadratic discriminant analysis using diversity measure. Nucleic Acids Research. 2003, 31 (21): 6214-6220.PubMed CentralView ArticlePubMedGoogle Scholar
- Hebsgaard SM, korning PG, Tolstrup N, Engelbrecht J, Rouze P, Brunak S: Splice site prediction in Arabidopsis Thaliana pre-mRNA by combining local and global sequence information. Nucleic Acids Research. 1996, 24: 3439-3452.PubMed CentralView ArticlePubMedGoogle Scholar
- Burge C: Modeling dependencies in pre-mRNA splicing signals. Computational methods in Molecular Biology. 1998, Elsevier press, chapter 8: 129-163.View ArticleGoogle Scholar
- Saeys Y, Degroeve S, Aeyels D, Peer Van de Y, Rouze P: Fast feature selection using a simple estimation of distribution algorithm: a case study on splice site prediction. Bioinformatics. 2003, Oxford Univ Press, 19: 179-188.Google Scholar
- Waddell P, Kishino H, Ota R: Very fast algorithms for evaluating the stability of ML and Bayesian phylogenetic trees from sequence data. Genome Informatics. 2002, 13: 82-92.PubMedGoogle Scholar
- Buckingham S: Bioinformatics: Programmed for success. Nature. 2003, 425 (6954): 209-215.Google Scholar
- Washietl S, Hofacker I, Stadler P: From The Cover: Fast and reliable prediction of noncoding RNAs. Proceedings of the National Academy of Sciences. 2005, 102 (7): 2454-View ArticleGoogle Scholar
- Shannon CE: A mathematical theory of communication. Bell System Tech J. 1948, 27: 379-423. 623–656.View ArticleGoogle Scholar
- Schneider TD: Information content of individual genetic sequences. Journal of Theoretical Biology. 1997, 189: 427-441.View ArticlePubMedGoogle Scholar
- Itoh H, Washio T, Masaru Tomita: Computational comparative analyses of alternative splicing regulation using full-length cDNA of various eukaryotes. RNA. 2004, 10: 1005-1018.PubMed CentralView ArticlePubMedGoogle Scholar
- Shapiro MB, Senapathy P: RNA splice junctions of different classes of eukaryotes: Sequence statistics and functional implications in gene expression. Nucleic Acids Research. 1987, 15: 7155-7174.PubMed CentralView ArticlePubMedGoogle Scholar
- Dror G, Sorek R, Shamir R: Accurate identification of alternatively spliced exons using support vector machine. Bioinformatics. 2004, 21 (7): 897-901.View ArticlePubMedGoogle Scholar
- Golub T, Slomin D, Tamayo P, Huard C, Gaasenbeek M, Mesirov J, Coller H, Loh M, Downing J, Caliguiri M, Bloomfield C, Lander E: Molecular classification of cancer: class discovery and class prediction by gene expression monitoring. Science. 1999, 286: 531-537.View ArticlePubMedGoogle Scholar
- Cortes C, Vapnik V: Support vector network, Machine Learning. Machine Learning. 1995, 20: 273-293.Google Scholar
- Vapnik v: The nature of statistical learning theory. 1995, Springer, New YorkView ArticleGoogle Scholar
- Cristianini N, Shawe-Taylor J: An introduction to support vector machine and kernel based learning methods. 2000, Cambridge University press, CambridgeView ArticleGoogle Scholar
- Drucker H, Wu D, Vapnik V: Support vector machines for spam categorization. IEEE transaction on Neural Networks. 1995, 10: 1054-1084.Google Scholar
- Burge C: A tutorial on support vector machines for pattern recognition. Data Mining and Knowledge Discovery. 1998, 2: 121-167.View ArticleGoogle Scholar
- Liu H, Wong L: Data mining tools for biological sequences. Journal of bioinformatics and computational biology. 2003, 1: 139-160.View ArticlePubMedGoogle Scholar
- Yeo G, Burge CB: Maximum Entropy Modeling of Short Sequence Motifs with Applications to RNA Splicing Signals. Journal of Computational Biology. 2004, 11 (2–3): 377-394.View ArticlePubMedGoogle Scholar
Copyright
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.