BACKGROUND
rna interference  is a post-transcriptional gene regulatory mechanism by which a double-stranded rna  induces sequence-specific gene silencing
<cit> .
the rnai pathway consists of multiple steps, including the cleavage of the dsrna by dicer to form a 19-nucleotide short interfering rna  with 3' overhangs, the incorporation of the sirna molecule into the rna-induced silencing complex , and the recognition of the target gene transcript by the risc-sirna complex to induce mrna degradation or translational repression.
in mammals and many other organisms, chemically synthesized sirna molecules can be introduced into cells to knockdown the expression of a specific gene.
because of the simplicity and low cost, sirna-based gene silencing has quickly become an important technique in functional genomics
<cit> .
since not all sirnas are equally effective, sirna design is one of the critical steps in gene silencing studies.
earlier experimental data have suggested several sets of empirical rules for designing potent sirnas.
for example, ui-tei et al.
<cit>  proposed several criteria for potent sirnas, including the presence of au-rich 5' terminal region and g/c at the 3' end of the antisense strand, and the absence of long gc stretches .
the thermodynamic properties of sirna duplexes were shown to affect target mrna degradation
<cit> .
in addition, secondary structure in sirnas could reduce the efficacy of gene silencing
<cit> .
although these findings provided important insight into rna interference, the empirical rules were often derived from relatively small datasets, and thus might not cover all the relevant features affecting sirna potency.
with the accumulation of sirna data, machine learning methods have been developed for both classification and regression analysis of sirna potency.
strom
<cit>  reported a genetic programming  method with boosting algorithms for the binary classification of effective and ineffective sirnas.
it was shown that the boosted gp classifier outperformed support vector machines  trained on the same dataset.
ladunga
<cit>  trained svms for sirna classification based on biophysical signatures of free energy, target site accessibility and dinucleotide characteristics.
wang et al.
<cit>  also constructed an svm classifier to identify hyperfunctional sirnas by using simple sequence and structural features such as base composition at each position, gc content, and secondary structure.
artificial neural networks  have been trained for regression analysis of sirna efficacy.
huesken et al.
<cit>  developed the biopredsi model with nucleotide sequence information to predict the inhibitory activity of sirnas, and used the ann model to design a human sirna library.
shabalina et al.
<cit>  performed thermodynamic and correlation analyses on a set of sirnas, and constructed an ann model with three parameters characterizing sirna sequences.
in addition, vert et al.
<cit>  constructed the simple linear model dsir using the lasso procedure for sirna efficacy prediction with basic sequence features.
the above-mentioned previous studies suggest that many sirna features of sequence composition, thermodynamic stability and secondary structure are related to the effectiveness of gene silencing.
however, few advanced methods have been used to select the most relevant features for predicting sirna potency.
most of the previous studies selected some sirna features based on empirical knowledge or simple statistical analysis .
more recently, klingelhoefer et al.
<cit>  used a stochastic logistic regression-based algorithm to identify relevant features associated with sirna potency.
the feature selection method revealed several sequence motifs such as ucu in potent sirnas.
in this study, random forests  were constructed to select important sequence features for predicting sirna potency.
rf-based variable importance measures were previously used in microarray expression data analyses to select a relatively small set of informative genes for disease/sample classification
<cit> .
however, it remained to be demonstrated whether the rf methodology could be applied to sequence feature selection for sirna classification.
the selected sequence features were used to construct svm models for predicting sirna potency.
our results suggest that sirna potency is significantly affected by its nucleotide dimer and trimer compositions.
some nucleotide motifs such as ucc appear to be positively correlated with sirna efficacy, whereas other motifs such as gag may have a negative effect on gene expression inhibition.
the findings will likely be useful for rational sirna design in large-scale functional genomics projects.
methods
data
as described in the previous study
<cit> , a non-redundant set of experimentally evaluated sirnas were collected from several published studies.
for each sirna, the relative level of target gene mrna was the ratio of the remaining mrna level after sirna treatment to the wild-type control level.
the relative mrna level  was used to determine the effectiveness of mrna knockdown.
effective sirnas gave rise to lower levels of remaining gene expression.
the cut-off level of  <dig>  was used to define positive instances  and negative instances .
the dataset used in this study consisted of  <dig> positive instances and  <dig> negative instances.
each sirna instance was a sequence of  <dig> nucleotides  representing the antisense strand of the target gene mrna transcript.
the two-nucleotide overhang at the 3' end of a sirna was not included in the data instance.
although more sirna data have recently become available and will be used in the future study, our approach was first tested on this relatively small dataset so that our findings could be compared with the previously published results.
sequence features
of many potential features, only a few may be relevant for sirna classification.
this study examined  <dig> sequence features belonging to six groups .
the first group has  <dig> features, each of which is the nucleotide identity  of a sequence position in a sirna.
in this study, real numeric values are used to represent the nucleotide identity:  <dig>  for a,  <dig>  for u,  <dig>  for g, and  <dig>  for c. the second feature group indicates the base composition of a sirna .
the third group has  <dig> features representing the frequencies of all possible dinucleotides .
the fourth feature group consists of  <dig> frequencies of all possible trinucleotides .
the fifth feature group indicates the global and local g/c contents in a sirna.
there are  <dig> features, one for the overall g/c content and  <dig> for local g/c contents.
with a sliding window size of five nucleotides, local g/c contents are calculated for all the possible windows along a 19-nucleotide sirna sequence.
the sixth group has the feature of sirna secondary structure stability.
the free energy of secondary structure was calculated by using the rnafold program in the vienna rna package
<cit> .
random forests for feature selection
in this study, important sequence features for sirna classification were identified by using the random forest  algorithm implemented in the software package available at http://www.stat.berkeley.edu/~breiman/randomforests/.
the rf algorithm uses a combination of independent decision trees to model data and measure variable importance
<cit> .
each decision tree in a forest is constructed using a bootstrap sample from the data, and about one-third of the data instances are not used to grow the tree.
these instances are called the out-of-bag  data for the tree.
at each node of the tree, m variables out of all the n input variables  are randomly selected, and the tree node is split using the selected m variables.
the random selection of features at each node decreases the correlation between the trees in the forest.
thus, the rf algorithm can handle many redundant features and avoid model overfitting.
it has been shown that rfs outperform adaboost ensembles on noisy datasets, and can perform well on data with many weak input variables
<cit> .
to evaluate the importance of variable x, its values in the oob instances associated with each tree in the forest are permuted randomly.
the permuted oob instances as well as the original oob instances are then classified using the tree.
the number of correct classifications on the original oob instances is subtracted by the number of predictions for the correct class on the permuted oob instances to calculate a raw score based on the tree.
the importance score of variable x is defined as the average of raw scores over all the trees in the forest.
for a fixed number of trees in the forest, the larger importance score a variable has, the more important it is for classification.
in addition, a z-score can be obtained by dividing the variable importance score by its standard error, and a statistical significance level may be assigned to the z-score assuming normality
<cit> .
support vector machine classifiers
support vector machines  were trained with the selected sequence features to predict sirna potency.
the svmlight software package  was used to construct svm classifiers.
the svm learning algorithm has been applied to a variety of biological problems for pattern classification, and may have superior generalization power with the ability to avoid overfitting
<cit> .
for a given set of binary-labelled training examples, the svm algorithm maps the input space into a higher-dimensional space, and seeks a hyperplane to separate the positive data instances from the negative ones
<cit> .
the optimal hyperplane maximizes the separation margin between the two classes of training data, and is defined by a fraction of the input data instances  close to the hyperplane.
the distance measurement between the data points in the high-dimensional space is defined by the kernel function.
in this study, we used the radial basis function  kernel:
where and are two data vectors, and γ is a training parameter.
a smaller γ value makes the decision boundary smoother.
another parameter for svm training is the regularization factor c, which controls the trade-off between low training error and large margin
<cit> .
different values for the γ and c parameters have been tested in this study to optimize the classifier performance.
classifier performance evaluation
we used a fivefold cross-validation approach to evaluate the performance of svm classifiers.positive and negative instances were distributed randomly into five folds.
in each of the five iterative steps, four of the five folds were used to train a classifier, and then the classifier was evaluated using the holdout fold .
the predictions made for the test instances in all the five iterations were combined and used to compute the following performance measures:
where tp is the number of true positives; tn is the number of true negatives; fp is the number of false positives; and fn is the number of false negatives.
in addition to overall accuracy, sensitivity and specificity, matthews correlation coefficient  is also commonly used as a measure of the quality of binary classifications
<cit> .
mcc measures the correlation between predictions and the actual class labels.
the receiver operating characteristic  curve is probably the most robust approach for classifier evaluation and comparison
<cit> .
the roc curve is drawn by plotting the true positive rate  against the false positive rate, which equals to .
in this work, the roc curve has been generated by varying the output threshold of a classifier and plotting the true positive rate against false positive rate for each threshold value.
the area under the roc curve  can be used as a reliable measure of classifier performance
<cit> .
since the roc plot is a unit square, the maximum value of auc is  <dig>  which is achieved by a perfect classifier.
weak classifiers have auc values close to  <dig> .
RESULTS
random forest-based selection of important features
there are many potential features for sirna classification.
to select the important features, sirna sequences were coded with the  <dig> features shown in table  <dig>  and then used to construct random forests  with different settings of the m parameter .
as suggested by the rf software
<cit> , the ceiling of the square root of the total number of input variables might be used as the default value of m .
in this study,  <dig> rfs were constructed by varying the m parameter setting from  <dig> to  <dig> .
each rf with  <dig> trees selected the top  <dig> features based on the z-score of variable importance.
some of the common features selected by the rfs were then identified for sirna classification.
the use of multiple rfs might increase the reliability for identifying relevant features.
some trinucleotide or dinucleotide features show negative correlation with sirna efficacy.
for example, the frequencies of cag and gag on the antisense strand have average z-scores of  <dig>  and  <dig> , respectively, and are negatively correlated with sirna efficacy .
the pearson's correlation coefficients are - <dig>  for cag% and - <dig>  for gag%.
the other composition features with negative effects on sirna efficacy include gca%, aua%, cug%, ag%, gg%, gga%, and ggc% .
these nucleotide motifs should be avoided in designing potent sirnas.
several features of base composition and g/c content were also selected by the rfs.
as shown in table  <dig>  the frequency of g on the antisense strand  is negatively correlated with sirna efficacy , whereas u% shows a positive correlation with gene expression inhibition .
both the global g/c content and the g/c content in the first  <dig> bases  are negatively correlated with sirna efficacy .
in addition, the nucleotide identity at the third position  of the antisense strand may be an important feature for sirna classification.
a nucleotide c at the third position has a negative effect on sirna efficacy .
however, the nucleotide identities at the other positions as well as sirna secondary structure stability were not selected by most of the rfs.
some of the important features in table  <dig> were previously shown to be related to the efficacy of gene silencing.
very high g/c contents were found to have a negative effect on sirna efficacy  <cit> .
the presence of au-rich 5' terminal region of the antisense strand was proposed to be one of the criteria for designing potent sirnas
<cit> .
consistent with the previous findings, both the overall g/c% and the 5' terminal g/c%  of the antisense strand were selected by the rfs in this study.
it was previously shown that the frequency of u, but not g or gg, was positively correlated with sirna efficacy
<cit> .
this observation has also been confirmed in the present study.
nevertheless, many of the other features in table  <dig>  including the top three features , have not been well documented in the literature.
thus, the findings in this study provide new insights into the rational design of potent sirnas.
support vector machine classifiers of sirna potency
to validate the important features selected by the rfs, support vector machines  were trained with these features for sirna classification.
the sirna instances were labeled with either 'potent'  or 'non-potent' .
classifier performance was evaluated using a fivefold cross-validation approach.
different settings of svm training parameters were tested, and the best classifier was constructed with γ =  <dig>  and c =  <dig> .
as shown in table  <dig>  the svm classifier, named rf_features, achieved  <dig> % overall accuracy with  <dig> % sensitivity and  <dig> % specificity.
the matthews correlation coefficient  was  <dig> , and the area under the roc curve  reached  <dig> .
thus, the features selected by the rfs can be used to construct relatively accurate svm models for predicting sirna potency.
for performance comparison, svm classifiers were also constructed using all the  <dig> features or only the  <dig> sirna sequence features .
as shown in table  <dig>  the classifier rf_features outperformed these two classifiers, namely all_features and seq_features.
the all_features classifier  achieved  <dig> % overall accuracy with  <dig> % sensitivity and  <dig> % specificity, mcc =  <dig> , and auc =  <dig> .
the results suggest that there may be some redundant or correlated information in the full feature set, and the rf-based feature selection can be used to improve classifier performance.
the seq_features classifier  showed significantly worse performance than rf_features , suggesting that the selected composition features contain some important information for sirna classification.
in figure  <dig>  the receiver operating characteristic  curves of the three svm classifiers are compared, and the result is consistent with the performance measures shown in table  <dig>  the classifier rf_features appears to be slightly more accurate than all_features.
both rf_features and all_features are significantly better than seq_features.
to further evaluate the classifier rf_features, we examined the svm output used to predict sirna potency.
the classifier was constructed using sirna instances with binary labels .
the actual magnitude of gene expression inhibition was not included in the training data.
in figure  <dig>  the svm output for each test instance is plotted against the level of gene expression inhibition.
interestingly, the output from rf_features is positively correlated with the level of gene expression inhibition .
the result further suggests that the classifier rf_features has learned some important sirna patterns related to the efficacy of gene silencing.
the classifier rf_features compared favorably with several existing models.
in the previous study by s trom  <cit> , a similar dataset of sirnas was used to construct svm and genetic programming  classifiers.
the most accurate svm classifier achieved an auc of  <dig> , mcc =  <dig> , 61% sensitivity and 68% specificity.
the boosted gp classifier showed better performance with auc =  <dig> , mcc =  <dig> , 50% sensitivity and 73% specificity.
in this study, the classifier rf_features achieved an auc of  <dig> , mcc =  <dig> ,  <dig> % sensitivity and  <dig> % specificity.
however, it is not straightforward to compare rf_features with the other existing models.
in two previous studies
<cit> , svms were also trained for sirna classification, but with different definitions of positive and negative instances.
in several other studies
<cit> , regression models instead of classifiers were constructed for predicting sirna efficacy.
CONCLUSIONS
we have developed a new machine learning approach for predicting sirna potency based on random forests and support vector machines.
since there were many potential features for sirna classification, random forests were used for feature selection based on variable importance scores.
interestingly, most of the selected features were nucleotide dimer and trimer compositions of sirna sequence.
some nucleotide motifs  showed positive correlation with sirna efficacy, whereas other motifs  might have a negative effect on gene silencing.
these important features were used to train support vector machines for predicting sirna potency with relatively high accuracy.
in the future, we will apply our approach to a large, integrated dataset of sirnas, and develop a software system for rational sirna design in functional geneomic studies.
competing interests
the authors declare that they have no competing interests.
authors' contributions
lw initiated and designed the study.
lw and ch conducted the data analysis.
lw drafted the manuscript.
jyy participated in experiment design, result interpretation and manuscript preparation.

 rna interference is a post-transcriptional gene regulatory mechanism by which a double-stranded rna induces sequence-specific gene silencing. the rnai pathway consists of multiple steps, including the cleavage of the dsrna by dicer to form a 19-nucleotide short interfering rna with 3' overhangs, the incorporation of the sirna molecule into the rna-induced silencing complex, and the recognition of the target gene transcript by the risc-sirna complex to induce mrna degradation or translational repression. in mammals and many other organisms, chemically synthesized sirna molecules can be introduced into cells to knockdown the expression of a specific gene. because of the simplicity and low cost, sirna-based gene silencing has quickly become an important technique in functional genomics. since not all sirnas are equally effective, sirna design is one of the critical steps in gene silencing studies. earlier experimental data have suggested several sets of empirical rules for designing potent sirnas. for example, ui-tei et al. proposed several criteria for potent sirnas, including the presence of au-rich 5' terminal region and g/c at the 3' end of the antisense strand, and the absence of long gc stretches. the thermodynamic properties of sirna duplexes were shown to affect target mrna degradation. in addition, secondary structure in sirnas could reduce the efficacy of gene silencing. although these findings provided important insight into rna interference, the empirical rules were often derived from relatively small datasets, and thus might not cover all the relevant features affecting sirna potency. with the accumulation of sirna data, machine learning methods have been developed for both classification and regression analysis of sirna potency. strom reported a genetic programming method with boosting algorithms for the binary classification of effective and ineffective sirnas. it was shown that the boosted gp classifier outperformed support vector machines trained on the same dataset. ladunga trained svms for sirna classification based on biophysical signatures of free energy, target site accessibility and dinucleotide characteristics. wang et al. also constructed an svm classifier to identify hyperfunctional sirnas by using simple sequence and structural features such as base composition at each position, gc content, and secondary structure. artificial neural networks have been trained for regression analysis of sirna efficacy. huesken et al. developed the biopredsi model with nucleotide sequence information to predict the inhibitory activity of sirnas, and used the ann model to design a human sirna library. shabalina et al. performed thermodynamic and correlation analyses on a set of sirnas, and constructed an ann model with three parameters characterizing sirna sequences. in addition, vert et al. constructed the simple linear model dsir using the lasso procedure for sirna efficacy prediction with basic sequence features. the above-mentioned previous studies suggest that many sirna features of sequence composition, thermodynamic stability and secondary structure are related to the effectiveness of gene silencing. however, few advanced methods have been used to select the most relevant features for predicting sirna potency. most of the previous studies selected some sirna features based on empirical knowledge or simple statistical analysis. more recently, klingelhoefer et al. used a stochastic logistic regression-based algorithm to identify relevant features associated with sirna potency. the feature selection method revealed several sequence motifs such as ucu in potent sirnas. in this study, random forests were constructed to select important sequence features for predicting sirna potency. rf-based variable importance measures were previously used in microarray expression data analyses to select a relatively small set of informative genes for disease/sample classification. however, it remained to be demonstrated whether the rf methodology could be applied to sequence feature selection for sirna classification. the selected sequence features were used to construct svm models for predicting sirna potency. our results suggest that sirna potency is significantly affected by its nucleotide dimer and trimer compositions. some nucleotide motifs such as ucc appear to be positively correlated with sirna efficacy, whereas other motifs such as gag may have a negative effect on gene expression inhibition. the findings will likely be useful for rational sirna design in large-scale functional genomics projects. methods data as described in the previous study, a non-redundant set of experimentally evaluated sirnas were collected from several published studies. for each sirna, the relative level of target gene mrna was the ratio of the remaining mrna level after sirna treatment to the wild-type control level. the relative mrna level was used to determine the effectiveness of mrna knockdown. effective sirnas gave rise to lower levels of remaining gene expression. the cut-off level of was used to define positive instances and negative instances. the dataset used in this study consisted of positive instances and negative instances. each sirna instance was a sequence of nucleotides representing the antisense strand of the target gene mrna transcript. the two-nucleotide overhang at the 3' end of a sirna was not included in the data instance. although more sirna data have recently become available and will be used in the future study, our approach was first tested on this relatively small dataset so that our findings could be compared with the previously published results. sequence features of many potential features, only a few may be relevant for sirna classification. this study examined sequence features belonging to six groups. the first group has features, each of which is the nucleotide identity of a sequence position in a sirna. in this study, real numeric values are used to represent the nucleotide identity: for a, for u, for g, and for c. the second feature group indicates the base composition of a sirna. the third group has features representing the frequencies of all possible dinucleotides. the fourth feature group consists of frequencies of all possible trinucleotides. the fifth feature group indicates the global and local g/c contents in a sirna. there are features, one for the overall g/c content and for local g/c contents. with a sliding window size of five nucleotides, local g/c contents are calculated for all the possible windows along a 19-nucleotide sirna sequence. the sixth group has the feature of sirna secondary structure stability. the free energy of secondary structure was calculated by using the rnafold program in the vienna rna package. random forests for feature selection in this study, important sequence features for sirna classification were identified by using the random forest algorithm implemented in the software package available at http://www.stat.berkeley.edu/~breiman/randomforests/. the rf algorithm uses a combination of independent decision trees to model data and measure variable importance. each decision tree in a forest is constructed using a bootstrap sample from the data, and about one-third of the data instances are not used to grow the tree. these instances are called the out-of-bag data for the tree. at each node of the tree, m variables out of all the n input variables are randomly selected, and the tree node is split using the selected m variables. the random selection of features at each node decreases the correlation between the trees in the forest. thus, the rf algorithm can handle many redundant features and avoid model overfitting. it has been shown that rfs outperform adaboost ensembles on noisy datasets, and can perform well on data with many weak input variables. to evaluate the importance of variable x, its values in the oob instances associated with each tree in the forest are permuted randomly. the permuted oob instances as well as the original oob instances are then classified using the tree. the number of correct classifications on the original oob instances is subtracted by the number of predictions for the correct class on the permuted oob instances to calculate a raw score based on the tree. the importance score of variable x is defined as the average of raw scores over all the trees in the forest. for a fixed number of trees in the forest, the larger importance score a variable has, the more important it is for classification. in addition, a z-score can be obtained by dividing the variable importance score by its standard error, and a statistical significance level may be assigned to the z-score assuming normality. support vector machine classifiers support vector machines were trained with the selected sequence features to predict sirna potency. the svmlight software package was used to construct svm classifiers. the svm learning algorithm has been applied to a variety of biological problems for pattern classification, and may have superior generalization power with the ability to avoid overfitting. for a given set of binary-labelled training examples, the svm algorithm maps the input space into a higher-dimensional space, and seeks a hyperplane to separate the positive data instances from the negative ones. the optimal hyperplane maximizes the separation margin between the two classes of training data, and is defined by a fraction of the input data instances close to the hyperplane. the distance measurement between the data points in the high-dimensional space is defined by the kernel function. in this study, we used the radial basis function kernel: where and are two data vectors, and γ is a training parameter. a smaller γ value makes the decision boundary smoother. another parameter for svm training is the regularization factor c, which controls the trade-off between low training error and large margin. different values for the γ and c parameters have been tested in this study to optimize the classifier performance. classifier performance evaluation we used a fivefold cross-validation approach to evaluate the performance of svm classifiers.positive and negative instances were distributed randomly into five folds. in each of the five iterative steps, four of the five folds were used to train a classifier, and then the classifier was evaluated using the holdout fold. the predictions made for the test instances in all the five iterations were combined and used to compute the following performance measures: where tp is the number of true positives; tn is the number of true negatives; fp is the number of false positives; and fn is the number of false negatives. in addition to overall accuracy, sensitivity and specificity, matthews correlation coefficient is also commonly used as a measure of the quality of binary classifications. mcc measures the correlation between predictions and the actual class labels. the receiver operating characteristic curve is probably the most robust approach for classifier evaluation and comparison. the roc curve is drawn by plotting the true positive rate against the false positive rate, which equals to. in this work, the roc curve has been generated by varying the output threshold of a classifier and plotting the true positive rate against false positive rate for each threshold value. the area under the roc curve can be used as a reliable measure of classifier performance. since the roc plot is a unit square, the maximum value of auc is which is achieved by a perfect classifier. weak classifiers have auc values close to.  random forest-based selection of important features there are many potential features for sirna classification. to select the important features, sirna sequences were coded with the features shown in table and then used to construct random forests with different settings of the m parameter. as suggested by the rf software, the ceiling of the square root of the total number of input variables might be used as the default value of m. in this study, rfs were constructed by varying the m parameter setting from to. each rf with trees selected the top features based on the z-score of variable importance. some of the common features selected by the rfs were then identified for sirna classification. the use of multiple rfs might increase the reliability for identifying relevant features. some trinucleotide or dinucleotide features show negative correlation with sirna efficacy. for example, the frequencies of cag and gag on the antisense strand have average z-scores of and, respectively, and are negatively correlated with sirna efficacy. the pearson's correlation coefficients are - for cag% and - for gag%. the other composition features with negative effects on sirna efficacy include gca%, aua%, cug%, ag%, gg%, gga%, and ggc%. these nucleotide motifs should be avoided in designing potent sirnas. several features of base composition and g/c content were also selected by the rfs. as shown in table the frequency of g on the antisense strand is negatively correlated with sirna efficacy, whereas u% shows a positive correlation with gene expression inhibition. both the global g/c content and the g/c content in the first bases are negatively correlated with sirna efficacy. in addition, the nucleotide identity at the third position of the antisense strand may be an important feature for sirna classification. a nucleotide c at the third position has a negative effect on sirna efficacy. however, the nucleotide identities at the other positions as well as sirna secondary structure stability were not selected by most of the rfs. some of the important features in table were previously shown to be related to the efficacy of gene silencing. very high g/c contents were found to have a negative effect on sirna efficacy. the presence of au-rich 5' terminal region of the antisense strand was proposed to be one of the criteria for designing potent sirnas. consistent with the previous findings, both the overall g/c% and the 5' terminal g/c% of the antisense strand were selected by the rfs in this study. it was previously shown that the frequency of u, but not g or gg, was positively correlated with sirna efficacy. this observation has also been confirmed in the present study. nevertheless, many of the other features in table including the top three features, have not been well documented in the literature. thus, the findings in this study provide new insights into the rational design of potent sirnas. support vector machine classifiers of sirna potency to validate the important features selected by the rfs, support vector machines were trained with these features for sirna classification. the sirna instances were labeled with either'potent' or'non-potent'. classifier performance was evaluated using a fivefold cross-validation approach. different settings of svm training parameters were tested, and the best classifier was constructed with γ = and c =. as shown in table the svm classifier, named rf_features, achieved % overall accuracy with % sensitivity and % specificity. the matthews correlation coefficient was, and the area under the roc curve reached. thus, the features selected by the rfs can be used to construct relatively accurate svm models for predicting sirna potency. for performance comparison, svm classifiers were also constructed using all the features or only the sirna sequence features. as shown in table the classifier rf_features outperformed these two classifiers, namely all_features and seq_features. the all_features classifier achieved % overall accuracy with % sensitivity and % specificity, mcc =, and auc =. the results suggest that there may be some redundant or correlated information in the full feature set, and the rf-based feature selection can be used to improve classifier performance. the seq_features classifier showed significantly worse performance than rf_features, suggesting that the selected composition features contain some important information for sirna classification. in figure the receiver operating characteristic curves of the three svm classifiers are compared, and the result is consistent with the performance measures shown in table the classifier rf_features appears to be slightly more accurate than all_features. both rf_features and all_features are significantly better than seq_features. to further evaluate the classifier rf_features, we examined the svm output used to predict sirna potency. the classifier was constructed using sirna instances with binary labels. the actual magnitude of gene expression inhibition was not included in the training data. in figure the svm output for each test instance is plotted against the level of gene expression inhibition. interestingly, the output from rf_features is positively correlated with the level of gene expression inhibition. the result further suggests that the classifier rf_features has learned some important sirna patterns related to the efficacy of gene silencing. the classifier rf_features compared favorably with several existing models. in the previous study by s trom, a similar dataset of sirnas was used to construct svm and genetic programming classifiers. the most accurate svm classifier achieved an auc of, mcc =, 61% sensitivity and 68% specificity. the boosted gp classifier showed better performance with auc =, mcc =, 50% sensitivity and 73% specificity. in this study, the classifier rf_features achieved an auc of, mcc =, % sensitivity and % specificity. however, it is not straightforward to compare rf_features with the other existing models. in two previous studies, svms were also trained for sirna classification, but with different definitions of positive and negative instances. in several other studies, regression models instead of classifiers were constructed for predicting sirna efficacy.  we have developed a new machine learning approach for predicting sirna potency based on random forests and support vector machines. since there were many potential features for sirna classification, random forests were used for feature selection based on variable importance scores. interestingly, most of the selected features were nucleotide dimer and trimer compositions of sirna sequence. some nucleotide motifs showed positive correlation with sirna efficacy, whereas other motifs might have a negative effect on gene silencing. these important features were used to train support vector machines for predicting sirna potency with relatively high accuracy. in the future, we will apply our approach to a large, integrated dataset of sirnas, and develop a software system for rational sirna design in functional geneomic studies. competing interests the authors declare that they have no competing interests. authors' contributions lw initiated and designed the study. lw and ch conducted the data analysis. lw drafted the manuscript. jyy participated in experiment design, result interpretation and manuscript preparation. 
support vector machine classifiers were then constructed using the selected sequence features for predicting sirna potency. interestingly, gene expression inhibition is significantly affected by nucleotide dimer and trimer compositions of sirna sequence. the findings in this study should help design potent sirnas for functional genomics, and might also provide further insights into the molecular mechanism of rna interference. to facilitate gene functional studies, we have developed a new machine learning method to predict sirna potency based on random forests and support vector machines. since there were many potential sequence features, random forests were used to select the most relevant features affecting gene expression inhibition. for a target gene of interest, many sirna molecules may be designed, whereas their efficiency of expression inhibition often varies.
 short interfering rnas can be used to knockdown gene expression in functional genomics. for a target gene of interest, many sirna molecules may be designed, whereas their efficiency of expression inhibition often varies.  to facilitate gene functional studies, we have developed a new machine learning method to predict sirna potency based on random forests and support vector machines. since there were many potential sequence features, random forests were used to select the most relevant features affecting gene expression inhibition. support vector machine classifiers were then constructed using the selected sequence features for predicting sirna potency. interestingly, gene expression inhibition is significantly affected by nucleotide dimer and trimer compositions of sirna sequence.  the findings in this study should help design potent sirnas for functional genomics, and might also provide further insights into the molecular mechanism of rna interference. 3– august the isibm international joint conference on bioinformatics, systems biology and intelligent computing shanghai, china 