Volume 13 Supplement 8

## The International Conference on Intelligent Biology and Medicine (ICIBM) Genomics

# A Bayesian decision fusion approach for microRNA target prediction

- Dong Yue
^{1}, - Maozu Guo
^{2}, - Yidong Chen
^{3}and - Yufei Huang
^{1, 3}Email author

**13(Suppl 8)**:S13

**DOI: **10.1186/1471-2164-13-S8-S13

© Yue et al.; licensee BioMed Central Ltd. 2012

**Published: **17 December 2012

## Abstract

MicroRNAs (miRNAs) are 19-25 nucleotides non-coding RNAs known to have important post-transcriptional regulatory functions. The computational target prediction algorithm is vital to effective experimental testing. However, since different existing algorithms rely on different features and classifiers, there is a poor agreement among the results of different algorithms. To benefit from the advantages of different algorithms, we proposed an algorithm called BCmicrO that combines the prediction of different algorithms with Bayesian Network. BCmicrO was evaluated using the training data and the proteomic data. The results show that BCmicrO improves both the sensitivity and the specificity of each individual algorithm. All the related materials including genome-wide prediction of human targets and a web-based tool are available at http://compgenomics.utsa.edu/gene/gene_1.php.

## Background

Gene regulation in human genome assumes multiple modes including transcriptional regulation by the regulatory proteins or transcription factors (TFs), and post-transcriptional regulation by including most notably microRNA (miRNA). MiRNA is a small non-coding RNA that has been discovered to repress transcription and/or protein translation of hundreds of genes by binding to the 3' Untranslated Region (UTR) of target genes [1, 2]. Understanding the functions and regulatory mechanisms of miRNA comprises one of the most active areas of research; such understanding will greatly advance our knowledge about the complexity of gene regulation and will help us to identify new therapeutic targets for effective treatment of various diseases.

Identifying miRNAs' target genes is an important first step in elucidating its function. Past work produced many target prediction algorithms based on miRNA-target sequence paring including TargetScan [3–5], miRanda [6, 7], PicTar [8], mirTarget [9, 10], PITA [11], DianamicroT [12] and others [13–21]. However, the prediction results of existing algorithms are still of low precision (i.e., low percentage of true targets among the predicted targets) and poor sensitivity (i.e., small percentage of true targets being predicted). In a recent study [22], Bartel et al., validated the prediction results of TargetScan, miRanda, PicTar, and PITA using a mass spectrometry (MS) approach. It was found that two thirds of their predicted targets appeared to be false positives, indicating a precision of only about 30%. As a result, the existing algorithms still cannot be used as target screening for subsequent bench testing.

There seems to be a poor agreement between the results of different algorithms and yet they achieve similar performance; this fact indicates that different algorithms rely on different mechanisms in making prediction, each of which has its own advantages. Indeed, the aforementioned sequence-based algorithms make predictions based on various important features of miRNA and mRNA nucleotide sequence interaction. Although a few important features including "seed region complementary", "binding free energy", and "sequence conservation" are among the most common adopted ones, different algorithms do utilize different sets of features. The differences in features and classifiers contribute to the differences in their prediction results. It is therefore desirable to integrate the predictions of different algorithms in order to combine their different advantages.

To do so, we propose a Bayesian decision fusion algorithm, BCmicrO. The goal of this algorithm is to improve the performance of existing target prediction algorithms. BCmicrO explicitly models the distributions of prediction results for each algorithm based on a training dataset composed of carefully constructed positive and negative miRNA-target pairs. These distributions capture the distinctions among different algorithms and weigh the differences at the decision level. With these distributions, the integration of different decisions is carried out based on Bayesian Network (BN). We tested the performance BCmicrO (combining TargetScan, miRanda, PicTar, mirTarget, PITA, and DianamicroT) with our training data, and validate it on the proteomics data. BCmicrO show clear improvement.

## Methods

### Overview of BCmicrO

The goal of BCmicrO is to generate the probability of an mRNA to be the target of a mRNA by integrating the predictions of different existing algorithms. In this paper, we focus on integrating TargetScan, miRanda, PicTar, mirTarget, PITA and Diana-microT's prediction scores. It should be noted that predictions from additional algorithms can be included in a similar fashion. TargetScan utilizes mainly seed region complementary and sequence conservation features for identifying potential binding sites and also applies a linear regression model to combine UTR features including 3' pairing score, local AU content, and distance from nearest 3'UTR terminus to produce a prediction context score for a UTR. On the other hand, miRanda relies on nucleotide complementariness and binding free energy in making the prediction. In contract, PicTar assumes a Hidden Markov Model (HMM) for seed region complementary and binding free energy to predict the potential binding sites. MirTarget is a SVM based algorithm with 113 features defined for a miRNA and target pairs. The key of PITA is a novel miRNA-target interaction model, based on the experimental observation - a strong secondary structure formed by 3'UTR itself will prevent the binding of miRNA. Diana-microT is a rule based miRNA target prediction algorithm applying a modified dynamic programming algorithm to determine the minimum free energy for each segment with a miRNA.

### Model formulation

Genome-wide predictions of TargetScan, miRanda, PicTar, mirTarget, PITA and Diana-microT are all reported in terms of scores. Particularly, TargetScan predicts miRNA's potential binding sites in the mRNA's 3' UTR, a context score is calculated for each site and the total context score is computed to represent the confidence of an mRNA to be a target. MiRanda indentifies all possible target sites for an mRNA and the highest target site score is selected to represent the confidence of the corresponding mRNA being a target. PicTar and other algorithms also compute a score reflecting the likelihood that the mRNA is a target.

To integrate these scores, BCmicrO adopts a BN model. BN is also known as directed graphical models, where the links of the graphs have a particular directionality indicated by arrows. The unique feature of BN is that the joint distribution over all of the random variables can be decomposed into a product of factors, each depending only on a subset of the variables [23].

*x*

_{1},

*x*

_{2},...

*x*

_{6}denote the scores of a miRNA-mRNA pair by TargetScan, miRanda, PicTar, mirTarget, PITA and Diana-microT, respectively. Also, set

*y*as an indicator variable such that

*y*= 1, if the mRNA is a real miRNA target, and

*y*= 0, otherwise. The goal of BCmicrO is to calculate

*P*(

*y*= 1|

*x*

_{1},

*x*

_{2}, ...,

*x*

_{6}), the posterior probability of the mRNA to be the miRNA's target given the TargetScan, miRanda, PicTar, mirTarget, PITA and Diana-microT scores. In reality, not all six scores are available for a miRNA-mRNA pair. Commonly, each algorithm only provides the prediction scores meeting a cutoff threshold. Therefore, we introduce the score indicators

*s*

_{1},

*s*

_{2},...,

*s*

_{6}to denote whether TargetScan, miRanda, PicTar mirTarget, PITA and Diana-microT report scores, or

*s*

_{ i }= 1(

*i*∈ {1, 2, ..., 6}) if the algorithm i reports a score, and

*s*

_{ i }= 0 otherwise. Also note that

*x*

_{ i }may be a score or no score (NaN) because of the cutoff value that mentioned before. The posterior probability can be expressed based on the BN model as

*p*(

*y*= 1) is the prior probability of an mRNA being a target and

*p*(

*x*

_{1},

*x*

_{2}, ...,

*x*

_{6}|

*y*) is the likelihood function. It is noted from the graphical model that given y,

*x*

_{1},

*x*

_{2},...,

*x*

_{6}are conditional independent and thus

It becomes clear that *p*(*y* = 1|*x*_{1}, *x*_{2}, ..., *x*_{6}) can be calculated from (1)-(3) once we have the conditional distributions *p*(*x*_{
i
}|*s*_{
i
}*, y*) and *p*(*s*_{
i
}|*y*), for *s*_{
i
} ∈ {0, 1}, *y* ∈ {0, 1}, and *i* ∈ {1, 2, ...,6}. In addition, based on Graphical model (Figure 2), given y, the conditional distributions of different algorithms are independent, such as *p*(*x*_{1}|*s*_{1}, *y*) is conditional independent of *p*(*x*_{2}|*s*_{2}, *y*). In the later section, we will discuss the process of acquiring all the above conditional probabilities in detail.

### Training data construction

Since the desired conditional distributions depend on y, i.e. the true target status of the mRNA, we need to collect high confidence positive and negative miRNA-target pairs as training data.

**The positive miRNA-target pairs** are collected from miRecords, which stores high-quality experimentally verified miRNA targets [24]. Only mammalian - human, mouse and rat records (852 records) are of our interest, since our goal is to identify mammalian miRNA target. Moreover, Karginov et al. [25] found the mRNAs whose association with Argonaute 2 (Ago2) increased upon miRNA over-expression were more likely targeted by miRNA, and 293 mRNAs were obtained as miR-124's targets. Lastly, 22 experimentally validated miR-124 targets are also collected from paper [25]. We combined the above positive miRNA-target pairs, removed the duplicate records, and ended up with a set of 929 positive miRNA-target pairs.

**The negative miRNA-target pairs** are currently unavailable in any annotated database. We constructed our negative database from two sources. First, it is known that negative targets are mostly up-regulated under miRNA over-expression. Therefore, first of all, negative targets were extracted as the up-regulated genes in 20 microarray data due to miRNA over-expression from NCBI Gene Expression Omnibus (GEO). To assure the high quality of negative data, we only chose the most confident up-regulated genes by restricting the differential expression p value, the fold change and consistency of the samples over time whenever available. To be more specific, the differential expression p value of the negative target must be less than 0.001 to ensure it is differentially expressed and the fold change (FC) of the negative target must be greater than 1.5 to ensure it is not down-regulated. In this process, 3542 negative miRNA-target pairs were gained. This is a high confident negative set compared to those miRNA-gene pairs with un-changed expression or random sampling.

Second, we focus on the existing results of miR-124 using immunoprecipitation (IP) of Ago2, since this technology has both higher sensitivity and specificity than other technologies including microarray and proteomics. Therefore, we obtained 19780 negative miR-124 targets by excluding 22 luciferase validated targets validated and 293 miR-124 target genes predicted in [25]. In the end, 23319 negative miRNA-target pairs were acquired.

The prediction scores of the positive and negative pairs for the three algorithms were subsequently obtained. The TargetScan (v5.1) scores were downloaded from web site (http://www.TargetScan.org/) [4, 5]. miRanda (2008 Sept) scores were downloaded from web site (http://www.microrna.org) [6, 7]; PicTar (2006) target predictions were downloaded from web site (http://PicTar.mdc-berlin.de/) [8]; mirTarget prediction results were downloaded from website (http://mirdb.org/miRDB/download.html) [9, 10]; PITA scores were downloaded from (http://genie.weizmann.ac.il/pubs/mir07/index.html) [11]. Diana-microT scores are downloaded from (http://diana.cslab.ece.ntua.gr/microT/) [12].

### Training of the conditional distributions

1. *p*(*x*_{
i
} = *score*|*s*_{
i
} = 1, *y* = 1), *i* ∈ {1, 2, ...6}

*p*(

*x*

_{ i }=

*score*|

*s*

_{ i }= 1,

*y*= 1) is the probability of a miRNA-target pair's score of algorithm

*i*given this pair is a positive pair and has a score. Note that "

*x*

_{ i }=

*score*" means

*x*

_{ i }has score. To find this conditional distribution, we should obtain the positive miRNA-target pairs with scores. To this end, we searched the positive miRNA-target pairs in each algorithm prediction results. Specifically, for TargetScan prediction results, 199 scores were obtained for the positive miRNA-target pairs and the histogram is shown in Figure 3. Upon flipping the histogram horizontally at its maximum score, a Gamma distribution was fitted. Maximum likelihood estimator (MLE) was used to estimate the parameters of the Gamma distribution.

*p*(

*x*

_{ i }=

*score*|

*s*

_{ i }= 1,

*y*= 1) (

*i*= 2 to 6), can be determined similarly, with a total of 278,175, 214, 631 and 396 positive pairs with scores for miRanda (Additional file 1), PicTar (Additional file 2), mirTarget (Additional file 3), PITA (Additional file 4), and Diana-microT (Additional file 5), respectively.

2. *p*(*x*_{
i
} = *score*|*s*_{
i
} = 1, *y* = 0), *i* ∈ {1, 2, ... 6}

*p*(

*x*

_{ i }=

*score*|

*s*

_{ i }= 1,

*y*= 0) is the probability of a miRNA-target pair's score

*x*

_{ i }from algorithm

*i*, given this pair is a negative pair and has a reported score. Similar as

*p*(

*x*

_{ i }=

*score*|

*s*

_{ i }= 1,

*y*= 1), we searched the negative miRNA-target pairs prediction results. 1928 negative miRNA-target pairs with TargetScan scores were acquired. The Gamma distribution was fitted to the scores (Figure 4). For miRanda (Additional file 6), PicTar (Additional file 7), mirTarget (Additional file 8), PITA (Additional file 9), and Diana-microT (Additional file 10), 1230, 613,436, 8831 and 3254 scores for the negative pairs were obtained, respectively.

3. *p*(*s*_{
i
}|*y*), *i* ∈ {1, 2, ... 6}, *s*_{
i
} ∈ {0, 1}, *y* ∈ {0, 1}

*p*(

*s*

_{ i }= 0|

*y*= 0),

*p*(

*s*

_{ i }= 0|

*y*= 1),

*p*(

*s*

_{ i }= 1|

*y*= 0), and

*p*(

*s*

_{ i }= 1|

*y*= 1) are the true negative rate (TNR), false negative rate (FNR), false positive rate (FPR), and true positive rate (TPR) for each algorithm. Since the prediction is carried out genome-wide, they should be assessed for a data set of a similar composition of positive and negative targets for human genome. The real miR-124 targets are retrieved from [25] including 22 targets validated by luciferase and 256 Net IP enrichment identified target genes, and the rest of genes are considered as the negative miR-124 targets (19780). Table 1 enlisted the estimated probabilities of each algorithm. We assume that the performance of each algorithm is consistent for all miRNAs. As a byproduct, the prior

*p*(

*y*= 1) is estimated from the Net IP data of miR-124 as 0.0133.

TP, FP, TN and FN rate of the 3 algorithms

TP rate | FP rate | TN rate | FN rate | |
---|---|---|---|---|

TargetScan | 0.4082 | 0.0877 | 0.9123 | 0.5918 |

miRanda | 0.3371 | 0.0783 | 0.9217 | 0.6629 |

PicTar | 0.1798 | 0.0390 | 0.9610 | 0.8202 |

mirTarget | 0.2285 | 0.0302 | 0.9698 | 0.7715 |

PITA | 0.7603 | 0.3942 | 0.6058 | 0.2397 |

Diana-microT | 0.4045 | 0.1474 | 0.8526 | 0.5955 |

4. Other Conditional distributions

It is apparent that *p*(*x*_{
i
} = *score*|*s*_{
i
} = 0, *y*) = 0 for all *i*. Similarly, we have *p*(*x*_{
i
} = *NaN* | *s*_{
i
} = 0, *y*) = 1,*p*(*x*_{
i
} = *NaN* | *s*_{
i
} = 1, *y*) = 0

## Result

### Test of BCmicrO on training data

### Test of BCmicrO on proteomics data

*n*is the rank of the predictions and

*c*(

*t*) denotes the cumulative sum at rank

*t*. In table 2, 3, 4, 5, we calculate

*A*(

*n*) for

*n*∈ {100, 200, 300, 400, 500} for each algorithm. The lower

*A*(

*n*), the better the performance. In sum, BCmicrO has a clear advantage over the rest algorithm in top 400 and top 500 (Table 2 - Table 6). In addition, BCmcirO has the lowest cumulative sum at top 300 of miR-155 (Table 4). To quantify the consistency of the algorithm, we calculate the average cumulative sum of protein fold change in all miRNAs tests:

Cumulative protein down-fold for miR-let-7b

100 | 200 | 300 | 400 | 500 | |
---|---|---|---|---|---|

BCmicrO | -1010 | -2931 | -5369 | -8387 | -12164 |

PicTar | -626.6 | -626.6 | -626.6 | -626.6 | -626.6 |

mirTarget | -79.2 | -79.2 | -79.2 | -79.2 | -79.2 |

miRanda | -595.7 | -1439 | -1439 | -1439 | -1439 |

PITA | -289.9 | -1281 | -2690 | -4394 | -6221 |

Diana-microT | -1021 | -3393 | -5956 | -5956 | -5956 |

TargetScan | -6.9 | -6.9 | -6.9 | -6.9 | -6.9 |

Cumulative protein down-fold for miR-16

100 | 200 | 300 | 400 | 500 | |
---|---|---|---|---|---|

BCmicrO | -1644 | -5935 | -11682 | -17803 | -23623 |

PicTar | -1564 | -1720 | -1720 | -1720 | -1720 |

mirTarget | -1402 | -1402 | -1402 | -1402 | -1402 |

miRanda | -1170 | -4237 | -7143 | -7143 | -7143 |

PITA | -225.2 | -1194 | -2714 | -4745 | -6813 |

Diana-microT | -2058 | -6435 | -12499 | -18696 | -23338 |

TargetScan | -1062 | -1102 | -1102 | -1102 | -1102 |

Cumulative protein down-fold for miR-155

100 | 200 | 300 | 400 | 500 | |
---|---|---|---|---|---|

BCmicrO | -3001 | -9657 | -18582 | -28624 | -39434 |

PicTar | -217 | -217 | -217 | -217 | -217 |

mirTarget | -1061 | -1061 | -1061 | -1061 | -1061 |

miRanda | -2193 | -4683 | -4683 | -4683 | -4683 |

PITA | -868 | -3627 | -7210 | -12089 | -18034 |

Diana-microT | -3050 | -9657 | -15804 | -15804 | -15804 |

TargetScan | -3059 | -10008 | -15808 | -15808 | -15808 |

Cumulative protein down-fold for miR-30a

100 | 200 | 300 | 400 | 500 | |
---|---|---|---|---|---|

BCmicrO | -1298 | -3816 | -6982 | -10523 | -14403 |

PicTar | -1159 | -1237 | -1237 | -1237 | -1237 |

mirTarget | -1186 | -1432 | -1432 | -1432 | -1432 |

miRanda | -1006 | -3014 | -3183 | -3183 | -3183 |

PITA | -218 | -900 | -1713 | -2654 | -4112 |

Diana-microT | -1521 | -4331 | -7423 | -9550 | -9550 |

TargetScan | -99.8 | -99.8 | -99.8 | -99.8 | -99.8 |

Average cumulative protein down-fold for all miRNAs

100 | 200 | 300 | 400 | 500 | |
---|---|---|---|---|---|

BCmicrO | -869.3 | -2792 | -5327 | -8167 | -11203 |

PicTar | -445.9 | -475.2 | -475.2 | -475.2 | -475.2 |

mirTarget | -466.1 | -496.9 | -496.9 | -496.9 | -496.9 |

miRanda | -620.7 | -1671 | -2056 | -2056 | -2056 |

PITA | -200.2 | -875.3 | -1791 | -2985 | -4397 |

Diana-microT | -956 | -2977 | -5210 | -6250 | -6831 |

TargetScan | -528.5 | -1402 | -2127 | -2127 | -2127 |

Here, *A*_{
i
}(*n*) is the score function in (4), and *i* means the ith miRNA. F(n) is the final score at top *n* in all miRNAs tests. The scores of F(x) is shown in Table 6. BCmicrO ranks number 1 at top 300 - 500, which clearly showed that BCmicrO can consistently provide the best prediction when compared with individual algorithms in most of the cases.

## Conclusion

We proposed a new miRNA target prediction algorithm, BCmicrO, which combines the prediction result of 6 algorithms -PicTar, mirTarget, PITA, miRanda, DianamicroT, and TargetScan, using Bayesian Network.

Performance of BCmicrO was first validated based on the training data. It shows that BCmicrO has better AUC than the other 6 algorithms and also has higher sensitivity, given the same specificity. BCmicrO was also tested on proteomic data for miR-16, let-7b, miR-155, and miR-30a. BCmicrO achieved the lowest cumulative sum of protein fold change and proven to consistently deliver the best performance. BCmicrO is of low complexity and can be easily upgraded as each constituent algorithm improves itself. Additional algorithms can be also integrated into BCmicrO in a similar fashion.

## Declarations

### Acknowledgements

Y. Huang is supported by National Institute of Health (R01 CA096512, 5G12RR013646-12), and Qatar National Research Fund (09-874-3-235). M. Guo is supported by Natural Science Foundation of China (60932008 and 61172098) and Fundamental Research Funds for the Central Universities (HIT.ICRST.2010 022). The authors wish to acknowledge computational support provided by the UTSA Computational Systems Biology Core Facility (NIH RCMI grant 5G12RR013646-12).

This article has been published as part of *BMC Bioinformatics* Volume 13 Supplement 17, 2012: Eleventh International Conference on Bioinformatics (InCoB2012): Bioinformatics. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcbioinformatics/supplements/13/S17

## Authors’ Affiliations

## References

- Bushati N, Cohen S: microRNA functions. Annu Rev Cell Dev Biol. 2007, 23: 175-205. 10.1146/annurev.cellbio.23.090506.123406.View ArticlePubMed
- Ying S, Chang D, Lin S: The microRNA (miRNA): overview of the RNA genes that modulate gene function. Molecular biotechnology. 2008, 38 (3): 257-268. 10.1007/s12033-007-9013-8.View ArticlePubMed
- Lewis B, Shih I, Jones-Rhoades M, Bartel D, Burge C, et al: Prediction of mammalian microRNA targets. Cell. 2003, 115 (7): 787-798. 10.1016/S0092-8674(03)01018-3.View ArticlePubMed
- Lewis B, Burge C, Bartel D: Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005, 120: 15-20. 10.1016/j.cell.2004.12.035.View ArticlePubMed
- Grimson A, Farh K, Johnston W, Garrett-Engele P, Lim L, Bartel D: MicroRNA targeting specificity in mammals: determinants beyond seed pairing. Molecular cell. 2007, 27: 91-105. 10.1016/j.molcel.2007.06.017.PubMed CentralView ArticlePubMed
- Enright A, John B, Gaul U, Tuschl T, Sander C, Marks D, et al: MicroRNA targets in Drosophila. Genome biology. 2003, 5: 1-1.View Article
- John B, Enright A, Aravin A, Tuschl T, Sander C, Marks D: Human microRNA targets. PLoS biology. 2004, 2 (11): e363-10.1371/journal.pbio.0020363.PubMed CentralView ArticlePubMed
- Krek A, Grün D, Poy M, Wolf R, Rosenberg L, Epstein E, MacMenamin P, Da Piedade I, Gunsalus K, Stoffel M, et al: Combinatorial microRNA target predictions. Nature genetics. 2005, 37 (5): 495-500. 10.1038/ng1536.View ArticlePubMed
- Wang X: miRDB: a microRNA target prediction and functional annotation database with a wiki interface. Rna. 2008, 14 (6): 1012-1017. 10.1261/rna.965408.PubMed CentralView ArticlePubMed
- Wang X, El Naqa I: Prediction of both conserved and nonconserved microRNA targets in animals. Bioinformatics. 2008, 24 (3): 325-10.1093/bioinformatics/btm595.View ArticlePubMed
- Kertesz M, Iovino N, Unnerstall U, Gaul U, Segal E: The role of site accessibility in microRNA target recognition. Nature genetics. 2007, 39 (10): 1278-1284. 10.1038/ng2135.View ArticlePubMed
- Maragkakis M, Reczko M, Simossis V, Alexiou P, Papadopoulos G, Dalamagas T, Giannopoulos G, Goumas G, Koukis E, Kourtis K, et al: DIANA-microT web server: elucidating microRNA functions through target prediction. Nucleic acids research. 2009, 37 (suppl 2): W273-W276.PubMed CentralView ArticlePubMed
- Kiriakidou M, Nelson P, Kouranov A, Fitziev P, Bouyioukos C, Mourelatos Z, Hatzigeorgiou A: A combined computational-experimental approach predicts human microRNA targets. Genes & development. 2004, 18 (10): 1165-10.1101/gad.1184704.View Article
- Rehmsmeier M, Steffen P, Höchsmann M, Giegerich R: Fast and effective prediction of microRNA/target duplexes. Rna. 2004, 10 (10): 1507-1517. 10.1261/rna.5248604.PubMed CentralView ArticlePubMed
- Rusinov V, Baev V, Minkov I, Tabler M: MicroInspector: a web tool for detection of miRNA binding sites in an RNA sequence. Nucleic acids research. 2005, 33 (suppl 2): W696-W700.PubMed CentralView ArticlePubMed
- Burgler C, Macdonald P: Prediction and verification of microRNA targets by MovingTargets, a highly adaptable prediction method. BMC genomics. 2005, 6: 88-10.1186/1471-2164-6-88.PubMed CentralView ArticlePubMed
- Rajewsky N, Socci N: Computational identification of microRNA targets. Genome Biology. 2004, 5 (2): 5-5. 10.1186/gb-2004-5-2-p5.View Article
- Kim S, Nam J, Rhee J, Lee W, Zhang B: miTarget: microRNA target gene prediction using a support vector machine. BMC bioinformatics. 2006, 7: 411-10.1186/1471-2105-7-411.PubMed CentralView ArticlePubMed
- Miranda K, Huynh T, Tay Y, Ang Y, Tam W, Thomson A, Lim B, Rigoutsos I: A pattern-based method for the identification of MicroRNA binding sites and their corresponding heteroduplexes. Cell. 2006, 126 (6): 1203-1217. 10.1016/j.cell.2006.07.031.View ArticlePubMed
- SaeTrom O, Snøve O, SÆTROM P: Weighted sequence motifs as an improved seeding step in microRNA target prediction algorithms. Rna. 2005, 11 (7): 995-1003. 10.1261/rna.7290705.PubMed CentralView ArticlePubMed
- Huang J, Frey B, Morris Q: Comparing sequence and expression for predicting microRNA targets using GenMiR3. Pacific Symposium on Biocomputing. 2008, 13: 52-63.
- Baek D, VillÃ J, et al: The impact of microRNAs on protein output. Nature. 2008, 455 (7209): 64-71. 10.1038/nature07242.PubMed CentralView ArticlePubMed
- Bishop C: Pattern recognition and machine learning. 2006, springer New York, 4: en ligne) SS
- Xiao F, Zuo Z, Cai G, Kang S, Gao X, Li T: miRecords: an integrated resource for microRNA-target interactions. Nucleic acids research. 2009, 37 (suppl 1): D105-D110.PubMed CentralView ArticlePubMed
- Karginov F, Conaco C, Xuan Z, Schmidt B, Parker J, Mandel G, Hannon G: A biochemical approach to identifying microRNA targets. Proceedings of the National Academy of Sciences. 2007, 104 (49): 19291-10.1073/pnas.0709971104.View Article
- Selbach M, Schwanhäusser B, Thierfelder N, Fang Z, Khanin R, Rajewsky N: Widespread changes in protein synthesis induced by microRNAs. Nature. 2008, 455 (7209): 58-63. 10.1038/nature07228.View ArticlePubMed

## 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.