Skip to main content

The discriminant power of RNA features for pre-miRNA recognition

Abstract

Background

Computational discovery of microRNAs (miRNA) is based on pre-determined sets of features from miRNA precursors (pre-miRNA). Some feature sets are composed of sequence-structure patterns commonly found in pre-miRNAs, while others are a combination of more sophisticated RNA features. In this work, we analyze the discriminant power of seven feature sets, which are used in six pre-miRNA prediction tools. The analysis is based on the classification performance achieved with these feature sets for the training algorithms used in these tools. We also evaluate feature discrimination through the F-score and feature importance in the induction of random forests.

Results

Small or non-significant differences were found among the estimated classification performances of classifiers induced using sets with diversification of features, despite the wide differences in their dimension. Inspired in these results, we obtained a lower-dimensional feature set, which achieved a sensitivity of 90% and a specificity of 95%. These estimates are within 0.1% of the maximal values obtained with any feature set (SELECT, Section “Results and discussion”) while it is 34 times faster to compute. Even compared to another feature set (FS2, see Section “Results and discussion”), which is the computationally least expensive feature set of those from the literature which perform within 0.1% of the maximal values, it is 34 times faster to compute. The results obtained by the tools used as references in the experiments carried out showed that five out of these six tools have lower sensitivity or specificity.

Conclusion

In miRNA discovery the number of putative miRNA loci is in the order of millions. Analysis of putative pre-miRNAs using a computationally expensive feature set would be wasteful or even unfeasible for large genomes. In this work, we propose a relatively inexpensive feature set and explore most of the learning aspects implemented in current ab-initio pre-miRNA prediction tools, which may lead to the development of efficient ab-initio pre-miRNA discovery tools.

The material to reproduce the main results from this paper can be downloaded from http://bioinformatics.rutgers.edu/Static/Software/discriminant.tar.gz.

Background

A microRNA (miRNA) is a small (approx. 17–25 nucleotides) non-coding RNA molecule (ncRNA) that modulates the stability of mRNA targets and their rate of translation into proteins [1]. MiRNAs are present in the genome of vertebrates, plants, algae and even viruses and are involved in diverse and complex biological processes, like development and cell differentiation [2], tumorigenesis [3] and immunity [4]. They can also alter plant gene expression in response to environmental stresses [5].

In animals, maturation of canonical miRNAs occurs in two steps: First, the long primary miRNA transcript is processed within the nucleus into a 60–120 nucleotides (nt) stem-loop hairpin precursor (pre-miRNA) by the enzyme Drosha [6]. Afterwards, within the cytoplasm, the enzyme Dicer cleaves the pre-miRNA into a double stranded RNA duplex (miRNA/miRNA*) and into a loop. The loop is degraded as a by-product [7], whereas the RNA duplex is unwound by helicase activity, releasing the mature miRNA and the star sequence [6]. The last is typically degraded whereas the mature miRNA guides the microribonucleo-protein complex (miRNP) to target messengers RNAs (mRNAs) by partial sequence complementarity [7].

Machine learning in miRNA recognition

Differently from protein-coding genes, ncRNA genes do not contain easily detectable signals in the sequence level [8]. Therefore, computational pipelines for miRNA discovery rely on characteristic features of pre-miRNAs. Coupled with comparative genomic computational pipelines, using RNAseq read libraries or completely ab-initio methods, machine learning (ML) algorithms have played an important role in miRNA discovery [7, 920]. ML algorithms induce models which are able to predict novel miRNAs based on patterns learned from known pre-miRNA sequences and from other RNA hairpin-like sequences, such as pseudo pre-miRNAs, transfer RNA (tRNA) and mRNAs. ProMiR [9, 21], a probabilistic method, searches for pre-miRNA in genomic sequences using sequence and structure features. Naïve Bayes based probabilistic models were proposed to score miRNAs [11] and pre-miRNA [7, 19] candidates. Support vector machines (SVMs) [10, 12, 13, 16, 17], random forest (RF) [14, 20], relaxed variable kernel density estimator (RVKDE) [15] and generalized Gaussian components based density estimation (G 2DE) [18] were used to induce classifiers for pre-miRNA prediction.

A limitation of working directly with the sequence is that the information available may not be sufficient to infer accurate models for novel miRNA prediction. For example, the approximate number of RNAs folding into hairpin-like secondary structures in the human genome, without filtering based on phylogenetic conservation, was estimated as 11 million [22]. Filters derived from other pre-miRNA features reduced the number of pre-miRNA candidates to around 5,300. Therefore, the use of features that consider different aspects may reduce false positive rates in miRNA detection.

Feature sets investigated in the literature

Some of the features commonly extracted from RNA sequences for pre-miRNA recognition may not help to distinguish between positive (true pre-miRNAs) and negative (pseudo pre-miRNAs) classes. Therefore, the feature sets considered may have an important effect in the learning process. Several feature sets have been proposed for pre-miRNA recognition [10, 13, 14, 16, 18, 20, 23]. To comparatively evaluate the effect of these different sets, we investigated seven feature sets proposed in the literature, named here FSi, i{1,..,7}. They were used to induce six classifiers and they contain most of the features employed in computational pipelines for pre-miRNA discovery. Next, we briefly define each of these sets.

The first set, named FS1 [16], has 48 sequence and structural features. The second feature set FS2, which corresponds to a subset of 21 features of FS1, was used by [16] to induce a classifier, named microPred. The third feature set FS3, was used to induce a classifier called G 2DE [18] and is composed by seven features also present in FS1. FS4 is a set of 32 sequence-structural features used by the triplet-SVM [10]. FS5 is a set of 1,300 sequence-structure motifs used by another classifier, mirident [23]. FS6 is the feature set used by the MiPred [14]. This feature set merged FS4, the minimum free energy of folding (MFE) and a stability measure (randfold). Finally, FS7 merged FS2 with four features of FS4 and three other features: percentage of low complexity regions detected in the sequence, maximal length of the amino acid string without stop codons and cumulative size of internal loops. It was used to induce the HuntMi [20]. The classification performances reported by the mentioned tools are among the highest in previous works. Their approximate specificities and sensitivities are: triplet-SVM (90%, 93%), MiPred (93%, 89%), microPred (97%, 90%), G 2DE (98%, 87%), mirident (99%, 98%) and HuntMi (97%, 95%).

Proposal and key findings

In this study, we investigated the discriminative power of seven RNA feature sets, previously used in six tools developed for pre-miRNA prediction. Among them are two sets composed of sequence-structure features (FS4 and FS5) and five sets are a miscellany of RNA features (FS1-FS3 and FS6-FS7). The investigation of a specific feature set, using a particular training data and learning algorithm, may insert learning biases. As a consequence, the predictive performance for other training sets and learning algorithms could be different. To deal with this problem, we evaluated each feature set using the learning algorithms, SVMs, RF and G 2DE, which were used in the publications proposing those feature sets. According to the experimental results, the miscellaneous feature sets produced more accurate predictive models than features sets composed from only sequence-structure patterns. However, the differences in accuracy among miscellaneous feature sets are small or insignificant, despite their large differences in composition and dimensionality. Inspired by these results, we selected a subset of 13 features, of lower computational cost, but with a similar classification performance, when compared with FS1-FS3 and FS6-FS7 feature sets. The classes of positive and negative test sets used in the experiments presented in this paper were predicted by the tools that we used as reference. Except for one tool, higher sensitivity was tied to lower specificity and vice versa.

Methods

Our goal was to investigate the predictive performance of RNA features in distinguishing pre-miRNAs from pseudo hairpins. As such, we adopted seven feature sets and three learning algorithms. The feature sets were used to induce classifiers for pre-miRNAs in triplet-SVM (FS4), MiPred (FS6), microPred (FS2 and FS1), G 2DE (FS3), mirident (FS5) and HuntMi (FS7). SVMs were used in triplet-SVM, microPred and mirident, whereas RF was used in MiPred and HuntMi. Generalized Gaussian density estimator (G 2DE) [24] is not a tool for pre-miRNA prediction in the sense that the features have to be computed by the user in his/her own pipeline. Nevertheless, we included G 2DE because of its predictive performance and class distribution interpretability. The subsections below provide details of our experiments.

Data sets

Human pre-miRNA sequences were downloaded from mirbase 19 [25] as the primary source of positive examples. In an attempt to avoid overfitting, we removed redundant sequences. As such, we clustered the available 1,600 human pre-miRNAs sequences using dnaclust [26] such that sequences within a cluster shared 80% similarity. Then, one sequence of each cluster was randomly picked. This yielded the set of positives composed of 1,377 non-redundant pre-miRNAs sequences.

The negative examples were the 8,494 pseudo hairpins from human RefSeq genes, originally obtained by [10] and subsequently used by [13, 14, 16, 18, 20, 23]. These sequences were obtained in order to keep basic features such as length distribution and minimal free energy of folding (MFE), similar to those observed in human pre-miRNAs. Moreover, this set has no redundant sequences. But, in order to adopt a uniform criterion for redundancy removal, we applied the same procedure adopted in the positive set. Only singleton clusters were formed. In practice, it is expected that high similarity between positive and negative examples leads to higher specificity [10].

Experiments

The predictive performance of any model is dependent on the training set representativeness, which usually increases with the increase in the training set size. Typically, density based algorithms, such as G 2DE, are more sensitive to the course of dimensionality and larger training sets are more likely to provide higher predictive performances. We determined experimentally the training set size which would be suitable for any algorithm and feature set. Each experiment was repeated 10 times, in order to provide standard deviations of each classification performance estimation. One repetition consisted of a test set, named here as GEN and 13 training sets. For a given repetition, GEN was composed by 459 sequences of each class, which corresponded to 1/3 of the 1,377 non-redundant pre-miRNA sequences. The remaining positive and negative sequences were used to sample increments of 67 sequences of each class to compose the training sets of 134,268,…,1,742 instances. Each feature set was computed from the same training and test sets. In total, we worked with 10 test sets and 13×10 training sets. The classification performances of the three algorithms converged to a threshold for training set sizes equal to 1,608, for all feature sets. Therefore, we presented the results for the largest training set, which contained 1,742 sequences.

Classification performance measures

Classification performance was measured as accuracy (Acc), sensitivity (Se), specificity (Sp), F-measure (Fm) and Mathew correlation coefficient (Mcc); see below. This measures can be computed as given below, such that TP, FN, TN and FP are the numbers of true positives, false negatives, true negatives and false positives, respectively.

Acc = 100 × TP + TN TP + FN + TN + FP Se = 100 × TP TP + FN Sp = 100 × TP TN + FP Fm = 100 × 2 × TP 2 × TP + FN + FP Mcc = 100 × TP × TN FP × FN ( ( TP + FP ) × ( TN + FN ) × ( TP + FN ) × ( TN + FP ) )

The first three measures are commonly used whereas Fm is prefered when a compromise between sensitivity and precision is desirable. Mcc measures the correlation between real and predicted classes and it is considered less biased towards class imbalance. We presented the predictive performances by the mean and the standard deviation (mean ± SD), over the 10 repetitions.

Features

Features used in this work are presented in Table 1, with references for the detailed descriptions. The prediction of the secondary structure in this work considered the energy model, as implemented in RNAfold [27] and UNAFold [28]. They predict the structure which gives the minimum free energy of folding (MFE). We kept the same parameters used in the original publications.

Table 1 Features used in each feature set

The sequence-structure features combined sequence nucleotide information and its predicted state at the secondary structure. In FS4, each feature represents the relative frequency of three contiguous nucleotides states at the secondary structure, fixing the middle character ({X s s s }, X{A,C,G,U} and s{paired, unpaired }). Because the triplet-SVM script excludes sequences with multiple loops, we implemented a Python script to compute FS4 in any sequence. The motifs in FS5 give the counts of its occurrence in the sequence-structure string. This string is obtained by padding the nucleotide sequence with its respective predicted state at the secondary structure ({left-paired, right-paired, unpaired }). This set was obtained using the Python script provided in the authors’ website. To compute FS1, we implemented a Python script, based on the microPred Perl pipeline. FS1 contains the largest diversity of features and depends on several independent scripts. We used the same scripts and options used in microPred. However, we used RNAfold from ViennaRNA-2.0 [37] and UNAFold v3.8 [28], instead of the older versions used in microPred. Initially, we obtained implausible values for features based on pair-probabilities. They were linked to the function get _pr() from the RNA Perl package of ViennaRNA-2.0. We bypassed this problem by restarting Perl for each new sequence. FS2 and FS3 were obtained from FS1. A customized Python script computed the randfold (p) and the MFE with the RANDFold [31] package and merged these two features with FS4 to obtain FS6. FS7 was obtained merging FS2 and the seven additional features obtained using the Python script obtained from the authors’ website.

Algorithms

The algorithms we adopted have different learning biases. This is important for the present work, since learning biases can play in favor of a feature set over others. SVM and RF are the two most applied algorithms for pre-miRNA classification, whereas G 2DE offered class distribution interpretability. Similar interpretation would be obtained using RVKDE but [18] showed that RVKDE produced accuracies similar to G 2DE and slightly lower than SVM, even though the number of kernels constructed by each algorithm were on average 920 (RVKDE), 361 (SVM) and six (G 2DE).

Support vector machines

SVMs deal with classification tasks by finding a hyperplane that separates training instances from two different classes with the maximum margin. The examples used to determine the hyperplane are the support vectors. Because many problems are not linearly separable, for these problems, the original feature space is mapped into a higher-dimensional space, where linear separation becomes feasible. Points from the original space are mapped to the new space by a kernel function. The RBF (radial basis function) kernel is a reasonable choice as it performs well for a wide range of problems [38]. For the training of SVMs, we used a Python interface for the library libsvm 3.12 [38]. This interface implements the C-SVM algorithm using the RBF kernel. The kernel parameters γ and C were tuned by 5-fold cross validation (CV) over the grid 2−5,2−3,…,215×2−15,2−13,…,23. The pair (C,γ) that led to the highest CV accuracy was used to train the SVMs using the whole training set. The induced model was then applied to the corresponding GEN test set.

Random forest

RF is an ensemble learning algorithm that induces a set of decision trees based on the concepts of “bagging” and random feature selection. Bagging is an important approach to estimate generalization error, whereas the latter is important to generate tree diversity. It was shown [39] that the strength of the ensemble depends on the strength of individual trees and the correlation between any two trees in the forest [40]. The number of features affects the strength of individual trees as well as the tree diversity, while the number of trees affects the generalization error. In order to obtain an ensemble with lower generalization error, a sufficiently large number of trees shall be chosen, taking into consideration two facts: RFs do not overfit, but limit the generalization error. This means that the number of trees has to be large enough to ensure lower generalization error, but after a certain value it does not have any effect on the generalization error estimate. For our experiments, we adopted the R package randomForest[40]. Each ensemble was generated over the grid (30,40,50,60,70,80,90,100,150,250,350,450)× ( 0.5 , 0.75 , 1 , 1.25 , 1.5 ) d , representing respectively the number of trees and the number of features. The is the default number of features tried in each node split and d is the dimension of the feature space. We chose the ensemble with the lowest generalization error over the grid and applied it to the corresponding GEN test set.

Generalized Gaussian density estimator

G 2DE [24] was designed to predict an instance class based on the probability density functions (pdf) of both positive and negative classes. Each pdf is fitted as mixture of generalized Gaussian components, using a limited user-defined number of components. One important feature of G 2DE is to provide the coefficients and parameters associated with these generalized components [24].

The learning process of G 2DE involves the estimation of the pdf parameters of each class, in addition to the weights of each component. If k is the maximum number of components and the feature space has dimension d, the number of parameters will be k(d+2)(d+1)/2. An evolutionary optimization algorithm finds the solution by maximizing the number of instances correctly classified in the training set plus the likelihood of class distributions. It requires two user-defined parameters: the number of Gaussian components (k) and the number of individuals for the initial population in the genetic algorithm (N). The first was kept six as in [18], and N was set to 100k, instead of 10k. High values of N implicate in more running time. On the other hand, it is expected that high N increases chances of findings an optimal solution. Since the solution is not deterministic, we ran G 2DE five times and chose the solution which gave the highest CV accuracy. The number five was determined by us through computational experiments.

Feature selection

Once the prediction model is induced, the highest computational cost in the evaluation of putative pre-miRNAs is the feature extraction from the sequences to be classified. Since this procedure is performed on millions of sequences, we performed feature selection on different data sets, excluding features which depend on shuffled sequences, which have a higher extraction cost. We also analyzed the importance of the whole set of 85 features obtained by combining all features from FS1 and FS6, and three features from FS7. The analyses of feature importance was performed using the feature importance estimated by randomForest [40], and the F-score, as described in [41]. Briefly, randomForest estimates the importance of a feature x i by computing the difference between the number of correctly classified out-of-bag (OOB) vectors before and after the permutation of x i in those vectors, during the training phase. For example, in the data used for this study, the number of OOB vectors was approximately 580. If an ensemble correctly classifies on average 500 and 36 OOB vectors before and after the permutation of x i , the estimated importance is approximately 464. This value indicates that x i was crucial for the correct classification. However, the interpretation must consider that the importance is conditioned to the induced ensemble. Thus, another feature x j with importance 200 in the ensemble including x i could obtain much higher importance in another ensemble excluding x i . Nevertheless, this measure provides a criterion to evaluate the relevance of each feature given the whole set. Differently, the F-score estimates the ratio of between and within classes distances and is computed before the learning step. Features with higher F-scores are more likely to be more discriminative, even though there is no objective criterion to decide on a specific score cut-off. In our experiments, we trained SVMs eliminating features with F-score below different score thresholds.

Results and discussion

Effect of feature sets and training algorithm

Data dimensionality may affect the learning process, particularly for parametric models. In our experiments, G 2DE only converged to a predictive model for the feature set FS3, which has only seven features. This algorithm uses a genetic algorithm to estimate the parameters of Gaussian components and their corresponding weights. In order to obtain the individuals for the initial population, the genetic algorithm uses another algorithm which generates random covariance matrices. In our experiments, this algorithm generated non-positive definite matrices, which caused the non convergence of G 2DE for higher dimensions.

Table 2 shows that feature sets composed by a miscellany of RNA features produced higher classification performances than feature sets composed by sequence-structure patterns. However, the small differences in classification performance produced by FS3 and FS6, compared to FS1, FS2 and FS7, shows that either the diversity or the dimensionality may affect classification. Indeed, FS6 has a larger dimension than FS2, but it is not composed of the same level of feature diversity than FS2. On the other hand, FS3 is more diverse than FS6, but does not contain enough features to produce sensitivity comparable to FS2.

Table 2 Predicted performances of classifiers trained with 1,742 examples, presented as the mean and standard deviation (Mean ± SD)

Nevertheless, it must be stressed that the highest predictive performances by the classifiers were not significantly different when the feature sets FS1, FS2, and FS7 were used (Table 2). As previously mentioned, FS2 and FS3 are subsets of FS1, while FS7 merged FS2 with seven additional features. These characteristics together with the very similar results obtained when using FS1, FS2, and FS7 suggest that the increase in the number of features leads to a limited increase of the predictive performance, even though the additional features were shown to be distinct features of pre-miRNAs [16, 20, 29, 31].

Feature discrimination and feature selection

Initially, we analyzed the importance of each one of the 85 features, obtained by merging FS1, FS6 and three features from FS7, when they were all used to induce classification ensembles by RF. In parallel, we also computed the F-scores. The Pearson correlation coefficient between the averages of these two measures was 75%, showing a relatively high correlation between importance and discrimination. Figure 1 shows the features whose importance was considered higher than 6. It can be seen that only 5 features obtained average importances higher than 40, corresponding to approx. 6% of the 580 OOB training instances. This results suggested that most of the 85 features may be redundant or irrelevant.

Figure 1
figure 1

Average feature importance estimated during the induction of RF ensembles. Features with importance lower than five were omitted. The average feature importance drops-off quickly after the 10th feature, indicating that for each ensemble there are few distinguishing features.

Interestingly, the features depeding on shuffled sequences appeared among those with the hightest importance or F-score. However, these features were not included in our feature selection step, due to their high computational cost and redundancy to the selected features. Moreover, since the features sets FS1 and FS7 share the 21 features of FS2 and they all produced classifiers with the highest predictive performances, we assumed that the relevant features for pre-miRNA classification were among these 21 common features. Therefore, the feature selection was performed using FS2, eliminating zD, which depends on shuffled sequences. The features selected (SELECT) from this set, in order of relevance, were: MFEI1, MFEI2, dG, dQ, dF, NEFE, Diff, dS, dS/L, |G-C|/L, |G-U|/L, %(G-U)/stems and MFEI3. Interestingly, six features of this set are energy-based measures (MFEI1, MFEI2, dG,NEFE, Diff, dS, MFEI3). The other relevant features are: entropy (dQ), compactness of the tree graph representation (dF), two thermodynamical features (dS and dS/L), normalized frequencies of G-C and G-U pairing (|G-C|/L, |G-U|/L, %(G-U)/stems).

The computational cost of each feature set was estimated by the computation time for a data set composed of 100 pre-miRNAs sequences randomly sampled from mirbase 20. Among the feature sets that produced the highest classification performance, FS1 had the highest cost (3:41 h), followed by FS7 (1:18 h), FS2 (1:17 h) and FS6 (39:03 min). Because SELECT and FS3 do not contain any stability measure, their costs are significantly lower than the cost for FS2. They took 2:17 min and 10 s to be computed. Among the sequence-structure based feature sets, FS5 took 3:24 min to be computed, whereas FS4 took 2 s.

The next comparison evaluates how the predictive performance associated with each feature set is affected by the use of different classifiers. For FS1, FS2, and FS7, the maximum difference in sensitivity and specificity between SVM and RF was 1.9% and 1.4%, respectively. The predictive performances of the three classifiers using FS3 were very similar. Thus, the learning biases of the three learning algorithms did not seem to have a significant effect on the predictive performance. This small effect of the learning biases is explained by the use of a sufficiently large training set, since most learning algorithms present a clear difference in their predictive performance only when small training sets are used.

Comparison with tools using the same algorithms and feature sets

In order to compare our results with tools used as references for our experiments, we predicted GEN test sets with those tools. Their main characteristics are summarized in Table 3. In Table 4 we show the predictive performance on the GEN sets obtained by the triplet-SVM, MiPred, microPred, G 2DE, mirident and HuntMi classifiers in our experiments. According to Table 4, except for the G 2DE tool, which uses the G 2DE algorithm, the predictive performance values obtained were much lower than the published values or the sensitivity was compromised by the specificity, or vice versa. The comparison in Table 4 used test sets obtained from mirbase 19, whereas the classifiers in those were induced with sequences of older releases. As the representativeness of the pre-miRNA population increases in newer releases, it is likely that the underlying distribution of the positive class would also changes. Therefore, the low sensitivities obtained by tools trained with old releases of mirbase, such as triplet-SVM, microPred, MiPred and mirident are not surprising. However, the low specificity values obtained by microPred and HuntMi, when compared to other older tools, were not expected. The specificity obtained by HuntMi in [20] was 72%, while we obtained 94% in our experiments, using the same algorithm and feature set. In contrast, the corresponding sensitivities were 99% and 88%. Likewise, the specificity of microPred [16] was 68%, while we found 95% using the same algorithm and feature set. Different results are usually obtained for experiments ran by different research groups, but not so different.

Table 3 Main characteristics of tools used as references in this work
Table 4 Predicted performance of tools used as references in our GEN test sets

The largest loss in specificity is observed for microPred and HuntMi tools, which both correct for class imbalance. That is, they attempt to correct the bias due to training sets formed of imbalanced number of examples in each class. This imbalance may cause a bias towards the majority class in the learning algorithm. Ideally, the class imbalance correction would increase the sensitivity without dropping the specificity. MicroPred increased sensitivity from 83% to 90%, while the specificity dropped slightly from 99% to 97%. The imbalance rate, the ratio of positive to negative examples, in the data set was 1:13. Likewise, HuntMi reported sensitivity and specificity values of 94% and 95%, working under an imbalance rate of 1:57. Since the ideal imbalance rate is determined experimentally, it is plausible that the class imbalance correction methods applied by microPred and HuntMi caused generalization problems. A contributing effect, or alternative explanation, might be that HuntMi uses negative sequence which differ greatly from positive examples, whereas the negative sequences used in microPred and in our experiments were selected to be similar to positive sequences. The observed loss of specificity might be countered with modifications to the training procedures. However, the lack of generalization of microPred was also mentioned in [7, 43].

G+C content effect

G+C content is an important feature during the folding of hairpin-like RNA sequences. Because G+C-rich sequences have more alternative high-energy stable binding-pairs, the prediction of the corresponding secondary structure is more complex. We drew the slopes of sensitivity and specificity correspondents to 12.5%-G+C content quantiles, to have a picture of the predictive performance of the feature sets and the algorithms in predicting G+C-rich sequences. Figure 2 shows that the variation in specificity throughout the intervals is random. Nevertheless, the sensitivities depended on the feature set and on the algorithm. All feature sets dropped the sensitivities of RF classifiers in G+C-rich pre-miRNAs. However, when FS1, FS2, FS7 and SELECT (not shown) were used to train SVM classifiers, only random variations in sensitivity along the 12.5%-G+C content quantile intervals were obtained. These four feature sets have %G+C related features in common, such as MFEI1 and normalized frequencies of G-C and G-U pairing (|G-C|/L, |G-U|/L, %(G-U)/stems). As Figure 1 shows, except for MFEI1, the other features appear with relatively low importance in the induction of ensembles by RF. On the other hand, the support vectors from the SVM model contain all the features used. These results confirmed the importance of including %G+C related features to detect G+C-rich pre-miRNA.

Figure 2
figure 2

Predictive performance of classifiers throughout 12 . 5 % -quantile distribution of G+C content. The prediction of the secondary structure of G+C-rich sequences is more challenging. This figure shows that the classification of G+C-rich pre-miRNA sequences is also more complex. As the G+C content increased, the sensitivity dropped, except when SVM was trained with feature sets including % G+C-based features (FS1, FS2 and FS7).

Scope of the investigation

The research reported in this paper was carried out using human sequences, one of the species with the highest abundance of positive sequences. Assuming that the larger the amount of positive sequences, the larger the amount of information about the human pre-miRNA population, the investigation performed with human sequences allowed a more fair comparison of the features sets and learning algorithms. Nevertheless, as it has been indicated that the rising of novel miRNAs is highly correlated with morphological complexity [4449], our results may vary for more distantly related species.

Conclusion

A considerable part of the computional cost involved in pre-miRNA prediction is due to the feature extraction from candidate sequences. Aiming to recommend effective and less costly sets of features, we investigated the discriminant power of seven RNA feature sets, under controlled sources of variation. Throughout extensive computational experiments, we showed that feature diversity is an important requirement in pre-miRNA recognition. Nevertheless, despite the discriminant power of individual features, higher dimensional sets did not produce higher classification performance classifiers. Based on these results, we proposed a smaller and less costly to compute subset of features, which produced classification performances as high as the produced by higher dimensional and more expensive sets. Because we attempted to avoid all possible sources of bias, we believe that the maximum classification performances reported here are the state-of-the-art for pre-miRNA prediction. Since these maximum classification performances are below experimentally feasible rates, other approaches to increase classification performance are welcome. As our tests showed, the tools used as references in our work either obtained low accuracies or the sensitivities or specificities were compromised.

References

  1. Khorshid M, Hausser J, Zavolan M, van Nimwegen E: A biophysical miRNA-mRNA interaction model infers canonical and noncanonical targets. Nat Methods. 2013, 10 (3): 253-255. 10.1038/nmeth.2341. [http://dx.doi.org/10.1038/nmeth.2341]

    Article  PubMed  CAS  Google Scholar 

  2. Letzen BS, Liu C, Thakor NV, Gearhart JD, All AH, Kerr CL: MicroRNA expression profiling of oligodendrocyte differentiation from human embryonic stem cells. PLoS One. 2010, 5 (5): e10480-10.1371/journal.pone.0010480. [http://dx.plos.org/10.1371/journal.pone.0010480]

    Article  PubMed Central  PubMed  Google Scholar 

  3. Cho WCS: MicroRNAs in cancer - from research to therapy. Biochimica et Biophysica Acta. 2010, 1805 (2): 209-217. [http://dx.doi.org/10.1016/j.bbcan.2009.11.003]

    PubMed  CAS  Google Scholar 

  4. Taganov KD, Boldin MP, Baltimore D: MicroRNAs and immunity: tiny players in a big field. Immunity. 2007, 26 (2): 133-137. 10.1016/j.immuni.2007.02.005. [http://dx.doi.org/10.1016/j.immuni.2007.02.005]

    Article  PubMed  CAS  Google Scholar 

  5. Burklew CE, Ashlock J, Winfrey WB, Zhang B: Effects of aluminum oxide nanoparticles on the growth, development, and microRNA expression of tobacco (Nicotiana tabacum). PloS One. 2012, 7 (5): e34783-10.1371/journal.pone.0034783. [http://dx.plos.org/10.1371/journal.pone.0034783]

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  6. Ritchie W, Gao D, Rasko JEJ: Defining and providing robust controls for microRNA prediction. Bioinformatics. 2012, 28 (8): 1058-1061. 10.1093/bioinformatics/bts114. [http://www.ncbi.nlm.nih.gov/pubmed/22408193]

    Article  PubMed  CAS  Google Scholar 

  7. Friedländer MR, Chen W, Adamidi C, Maaskola J, Einspanier R, Knespel S, Rajewsky N: Discovering microRNAs from deep sequencing data using miRDeep. Nat Biotechnol. 2008, 26 (4): 407-415. 10.1038/nbt1394. [http://dx.doi.org/10.1038/nbt1394]

    Article  PubMed  Google Scholar 

  8. Tran TT, Zhou F, Marshburn S, Stead M, Kushner SR, Xu Y: De novo computational prediction of non-coding RNA genes in prokaryotic genomes. Bioinformatics. 2009, 25 (22): 2897-2905. 10.1093/bioinformatics/btp537. [http://www.ncbi.nlm.nih.gov/pubmed/19744996]

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  9. Nam JW, Shin KR, Han J, Lee Y, Kim NV, Zhang BT: Human microRNA, prediction through a probabilistic co-learning model of sequence and structure. Nucleic Acids Res. 2005, 33 (11): 3570-3581. 10.1093/nar/gki668.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  10. Xue C, Li F, He T, Liu GP, Li Y, Zhang X: Classification of real and pseudo microRNA precursors using local structure-sequence features and support vector machine. BMC Bioinformatics. 2005, 6: 310-10.1186/1471-2105-6-310.

    Article  PubMed Central  PubMed  Google Scholar 

  11. Yousef M, Nebozhyn M, Shatkay H, Kanterakis S, Showe LCC, Showe MKK: Combining multi-species genomic data for microRNA identification using a Naive Bayes classifier machine learning for identification of microRNA genes. Bioinformatics. 2006, 22: 1325-1334. 10.1093/bioinformatics/btl094.

    Article  PubMed  CAS  Google Scholar 

  12. Hertel J, Stadler PF: Hairpins in a Haystack: recognizing microRNA precursors in comparative genomics data. Bioinformatics (Oxford, England). 2006, 22 (14): e197-e202. 10.1093/bioinformatics/btl257. [http://www.ncbi.nlm.nih.gov/pubmed/16873472]

    Article  CAS  Google Scholar 

  13. Ng KL, Mishra SK: De novo SVM classification of precursor microRNAs from genomic pseudo hairpins using global and intrinsic folding measures. Bioinformatics. 2007, 23 (11): 1321-1330. 10.1093/bioinformatics/btm026.

    Article  PubMed  CAS  Google Scholar 

  14. Jiang P, Wu H, Wang W, Ma W, Sun X, Lu Z: MiPred: classification of real and pseudo microRNA precursors using random forest prediction model with combined features. Nuc Ac Res. 2007, 35 (suppl 2): W339-W344.

    Article  Google Scholar 

  15. Chang D, Wang CC, Chen JW: Using a kernel density estimation based classifier to predict species-specific microRNA precursors. BMC Bioinformatics. 2008, 9 (Suppl 12): 241-

    Google Scholar 

  16. Batuwita R, Palade V: microPred: effective classification of pre-miRNAs for human miRNA gene prediction. Bioinformatics. 2009, 25 (8): 989-995. 10.1093/bioinformatics/btp107.

    Article  PubMed  CAS  Google Scholar 

  17. Gerlach D, Kriventseva EV, Rahman N, Vejnar CE, Zdobnov EM: miROrtho computational survey of microRNA genes. Nucleic Acids Res. 2009, 37 (Database issue): gkn707+-

    Google Scholar 

  18. Hsieh CH, Chang DTH, Hsueh CH, Wu CY, Oyang YJ: Predicting microRNA precursors with a generalized Gaussian components based density estimation algorithm. BMC Bioinformatics. 2010, 11 (Suppl 1): S52-10.1186/1471-2105-11-S1-S52. [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3009525&tool=pmcentrez&rendertype=abstract]

    Article  PubMed Central  PubMed  Google Scholar 

  19. Li N, You X, Chen T, Mackowiak SD, Friedländer MR, Weigt M, Du H, Gogol-Döring A, Chang Z, Dieterich C, Hu Y, Chen W: Global profiling of miRNAs and the hairpin precursors: insights into miRNA processing and novel miRNA discovery. Nucleic Acids Res. 2013, 41 (6): 3619-3634. 10.1093/nar/gkt072. [http://nar.oxfordjournals.org/content/41/6/3619.full]

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  20. Szcześniak MW, Sikora M, Makalowska I, Gudyś A: HuntMi: an efficient and taxon-specific approach in pre-miRNA identification. BMC Bioinformatics. 2013, 14: 83-10.1186/1471-2105-14-83. [http://www.biomedcentral.com/1471-2105/14/83]

    Article  PubMed Central  PubMed  Google Scholar 

  21. Nam JW, Kim J, Kim SK, Zhang BT: ProMiR II: a web server for the probabilistic prediction of clustered, nonclustered, conserved and nonconserved microRNAs. Nucleic Acids Res. 2006, 34 (Web-Server-Issue): 455-458.

    Article  Google Scholar 

  22. Bentwich I, Avniel A, Karov Y, Aharonov R, Gilad S, Barad O, Barzilai A, Einat P, Einav U, Meiri E, Sharon E, Spector Y, Bentwich Z: Identification of hundreds of conserved and nonconserved human microRNAs. Nat Genet. 2005, 37 (7): 766-770. 10.1038/ng1590. [http://dx.doi.org/10.1038/ng1590]

    Article  PubMed  CAS  Google Scholar 

  23. Liu X, He S, Skogerbø G, Gong F, Chen R: Integrated sequence-structure motifs suffice to identify microRNA precursors. PloS One. 2012, 7 (3): e32797-10.1371/journal.pone.0032797. [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3305290&tool=pmcentrez&rendertype=abstract]

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  24. Hsieh CH, Chang DTH, Oyang YJ: Data classification with a generalized Gaussian components based density estimation algorithm. Proceedings of the 2009 International Joint Conference on Neural Networks, IJCNN’09. 2009, Piscataway: IEEE Press, 2910-2917. [http://dl.acm.org/citation.cfm?id=1704555.1704693]

    Google Scholar 

  25. Kozomara A, Griffiths-Jones S: miRBase: integrating microRNA annotation and deep-sequencing data. Nucleic Acids Res. 2011, 39 (suppl 1): D152-D157.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  26. Ghodsi M, Liu B, Pop M: DNACLUST: accurate and efficient clustering of phylogenetic marker genes. BMC Bioinformatics. 2011, 12: 271+-10.1186/1471-2105-12-271. [http://dx.doi.org/10.1186/1471-2105-12-271]

    Article  PubMed Central  PubMed  Google Scholar 

  27. Hofacker IL: Vienna RNA secondary structure server. Nucleic Acids Res. 2003, 31 (13): 3429-3431. 10.1093/nar/gkg599.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  28. Markham NR, Zuker M: UNAFold: software for nucleic acid folding and hybridization. Methods Mol Biol. 2008, 453: 3-31. 10.1007/978-1-60327-429-6_1. [http://www.ncbi.nlm.nih.gov/pubmed/18712296]

    Article  PubMed  CAS  Google Scholar 

  29. NG Kwang Loong S, Mishra SK: Unique folding of precursor microRNAs quantitative evidence and implications for de novo identification. RNA. 2007, 13 (2): 170-187.

    Article  PubMed Central  PubMed  Google Scholar 

  30. Morgulis A, Gertz EM, Schäffer AA, Agarwala R: A fast and symmetric DUST implementation to mask low-complexity DNA sequences. J Comput Biol. 2006, 13 (5): 1028-1040. 10.1089/cmb.2006.13.1028. [http://online.liebertpub.com/doi/abs/10.1089/cmb.2006.13.1028]

    Article  PubMed  CAS  Google Scholar 

  31. Bonnet E, Wuyts J, Van de Peer Y, Rouzé P: Evidence that microRNA precursors, unlike other non-coding RNAs, have lower folding free energies than random sequences. Bioinformatics. 2004, 20 (17): 2911-2917. 10.1093/bioinformatics/bth374.

    Article  PubMed  CAS  Google Scholar 

  32. Mathews DH: Using an RNA secondary structure partition function to determine confidence in base pairs predicted by free energy minimization. RNA. 2004, 10 (8): 1178-1190. 10.1261/rna.7650904. [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=1370608&tool=pmcentrez&rendertype=abstract]

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  33. Freyhult E, Gardner PP, Moulton V: A comparison of RNA folding measures. BMC Bioinformatics. 2005, 6: 241-10.1186/1471-2105-6-241.

    Article  PubMed Central  PubMed  Google Scholar 

  34. Markham NR, Zuker M: DINAMelt web server for nucleic acid melting prediction. Nucleic Acids Res. 2005, 33 (Web Server issue): W577-W581. [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=1160267&tool=pmcentrez&rendertype=abstract]

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  35. Gan HH, Fera D, Zorn J, Shiffeldrim N, Tang M, Laserson U, Kim N, Schlick T: RAG: RNA-As-Graphs database–concepts, analysis, and features. Bioinformatics. 2004, 20: 1285-1291. 10.1093/bioinformatics/bth084.

    Article  PubMed  CAS  Google Scholar 

  36. Le SY, Chen JH, Maizel J: Thermodynamic stability and statistical significance of potential stem-loop structures situated at frameshift sites of retroviruses. Nucleic Acids Res. 1989, 17: 6143-6152. 10.1093/nar/17.15.6143.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  37. Lorenz R, Bernhart SH, Höner Zu, Siederdissen C, Tafer H, Flamm C, Stadler PF, Hofacker IL: ViennaRNA Package 2.0. orithm Mol Biol. 2011, 6: 26-[http://www.almob.org/content/6/1/26]

    Google Scholar 

  38. Chang CC, Lin CJ: LIBSVM: A library for support vector machines. ACM Trans Intell Syst Technol. 2011, 2: 27:1-27:27. [Software available at http://www.csie.ntu.edu.tw/~cjlin/libsvm]

    Article  Google Scholar 

  39. Breiman L, Schapire E: Random forests. Machine Learning. 2001, 45: 5-32. 10.1023/A:1010933404324. http://dx.doi.org/10.1023/A:1010933404324.

    Article  Google Scholar 

  40. Liaw A, Wiener M: Classification and regression by random Forest. R News. 2002, 2 (3): 18-22. [http://CRAN.R-project.org/doc/Rnews/]

    Google Scholar 

  41. Chen YW, Lin CJ: Combining, SVMs with various feature selection strategies. Feature Extraction, Volume 207 of Studies in Fuzziness and Soft Computing. Edited by: Guyon I, Nikravesh M, Gunn S, Zadeh LA. 2006, Berlin, Heidelberg: Springer Berlin Heidelberg, 315-324. [http://dx.doi.org/10.1007/978-3-540-35488-8_13]

    Google Scholar 

  42. Scott AJ, Knott M: A cluster analysis method for grouping means in the analysis of variance. Biometrics. 1974, 30 (3): 507-512. 10.2307/2529204. [http://dx.doi.org/10.2307/2529204]

    Article  Google Scholar 

  43. Mathelier A, Carbone A: MIReNA. Bioinformatics. 2010, 26 (18): 2226-2234. 10.1093/bioinformatics/btq329. [http://dx.doi.org/10.1093/bioinformatics/btq329]

    Article  PubMed  CAS  Google Scholar 

  44. Mattick JS: Non-coding RNAs: the architects of eukaryotic complexity. EMBO Reports. 2001, 2 (11): 986-991. 10.1093/embo-reports/kve230. [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=1084129&tool=pmcentrez&rendertype=abstract]

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  45. Mattick JS, Makunin IV: Non-coding RNA. Hum Mol Genet. 2006, 15 Spec No: R17-R29. [http://www.ncbi.nlm.nih.gov/pubmed/16651366]

    Article  Google Scholar 

  46. Heimberg AM, Sempere LF, Moy VN, Donoghue PC, Peterson KJ: MicroRNAs and the advent of vertebrate morphological complexity. Proc Natl Acad Sci USA. 2008, 105: 2946-2950. 10.1073/pnas.0712259105.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  47. Koerner MV, Pauler FM, Huang R, Barlow DP: The function of non-coding RNAs in genomic imprinting. Development. 2009, 136 (11): 1771-1783. 10.1242/dev.030403. [http://dx.doi.org/10.1242/dev.030403]

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  48. Matrajt M: Non-coding RNA in apicomplexan parasites. Mol Biochem Parasitol. 2010, 174: 1-7. 10.1016/j.molbiopara.2010.06.002. [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=2926168&tool=pmcentrez&rendertype=abstract]

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  49. De Mulder K, Berezikov E: Tracing the evolution of tissue identity with microRNAs. Genome Biol. 2010, 11 (3): 111-10.1186/gb-2010-11-3-111. [http://genomebiology.com/2010/11/3/111]

    Article  PubMed Central  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Empresa Brasileira de Pesquisa Agropecuária (Embrapa Soybean) for the continuum financial support to the first author and Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) for the grant awarded to the first author between May 2010 and February 2012. We also thank Dr. Hsieh and Dr. Gudy for their valuable time taken to provide details of G 2DE algorithm and HuntMi classifier, respectively.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Ivani de ON Lopes.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

AS and AC conceived and supervised the study. IL assembled the data, implemented the scripts, ran the experiments and summarized the results. The three authors wrote and approved the final manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( https://creativecommons.org/publicdomain/zero/1.0/ ) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

de ON Lopes, I., Schliep, A. & de LF de Carvalho, A.C. The discriminant power of RNA features for pre-miRNA recognition. BMC Bioinformatics 15, 124 (2014). https://doi.org/10.1186/1471-2105-15-124

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2105-15-124

Keywords