- Research
- Open Access
HuMiTar: A sequence-based method for prediction of human microRNA targets
https://doi.org/10.1186/1748-7188-3-16
© Ruan et al; licensee BioMed Central Ltd. 2008
- Received: 03 May 2008
- Accepted: 22 December 2008
- Published: 22 December 2008
Abstract
Background
MicroRNAs (miRs) are small noncoding RNAs that bind to complementary/partially complementary sites in the 3' untranslated regions of target genes to regulate protein production of the target transcript and to induce mRNA degradation or mRNA cleavage. The ability to perform accurate, high-throughput identification of physiologically active miR targets would enable functional characterization of individual miRs. Current target prediction methods include traditional approaches that are based on specific base-pairing rules in the miR's seed region and implementation of cross-species conservation of the target site, and machine learning (ML) methods that explore patterns that contrast true and false miR-mRNA duplexes. However, in the case of the traditional methods research shows that some seed region matches that are conserved are false positives and that some of the experimentally validated target sites are not conserved.
Results
We present HuMiTar, a computational method for identifying common targets of miRs, which is based on a scoring function that considers base-pairing for both seed and non-seed positions for human miR-mRNA duplexes. Our design shows that certain non-seed miR nucleotides, such as 14, 18, 13, 11, and 17, are characterized by a strong bias towards formation of Watson-Crick pairing. We contrasted HuMiTar with several representative competing methods on two sets of human miR targets and a set of ten glioblastoma oncogenes. Comparison with the two best performing traditional methods, PicTar and TargetScanS, and a representative ML method that considers the non-seed positions, NBmiRTar, shows that HuMiTar predictions include majority of the predictions of the other three methods. At the same time, the proposed method is also capable of finding more true positive targets as a trade-off for an increased number of predictions. Genome-wide predictions show that the proposed method is characterized by 1.99 signal-to-noise ratio and linear, with respect to the length of the mRNA sequence, computational complexity. The ROC analysis shows that HuMiTar obtains results comparable with PicTar, which are characterized by high true positive rates that are coupled with moderate values of false positive rates.
Conclusion
The proposed HuMiTar method constitutes a step towards providing an efficient model for studying translational gene regulation by miRs.
Keywords
- Seed Region
- Additional Prediction
- High True Positive Rate
- True Positive Prediction
- Traditional Predictor
Background
MicroRNAs (miRs) are endogenously expressed non-coding RNAs, which downregulate expression of their target mRNAs by inhibiting translational initiation or by inducing degradation of mRNA [1]. They are associated with numerous gene families in multi-cellular species and their regulatory functions in various biological processes are widespread [2–14]. The ability to perform accurate, high-throughput identification of physiologically active miR targets is one of the enabling factors for functional characterization of individual miRs. This is also true in case on human miRs, for which only a handful have been experimentally linked to specific functions. The methods for the prediction of miR targets can be subdivided into two classes, traditional approaches, which combine several factors such as sequence complementarity, minimization of free energy, and cross-species conservation, and machine learning (ML) methods that exploit statistical patterns that differentiate between true and false miR-mRNA duplexes. The former methods aim at finding target sites for a given miR by scanning 3' untranslated region (UTR) of the mRNA, while the latter methods classify a given duplex as true or false.
Current traditional sequence-based target predictors are based on the presence of a conserved 'seed region' (nucleotides 2–7) of exact Watson-Crick complementary base-pairing between the 3' UTR of the mRNA and the 5' end of the miR [15, 16]. They are based on two principles: (1) identification of potential miR binding sites according to specific base-pairing rules in the seed region, and (2) implementation of cross-species conservation [17]. Recent survey by Sethupathy and colleagues [18] compared five widely used traditional tools for mammalian target prediction which include DIANA-microT [7], miRanda [19], TargetScan [3], TargetScanS [11], and PicTar [10]. They observed that the earlier methods, i.e., TargetScan and DIANA-microT, achieve a relatively low sensitivity and predict a small number of targets. The miRanda was shown to provide a substantially better sensitivity as a trade-off for large increase in the total number of predictions. The two more recent programs, TargetScanS and PicTar, have almost identical sensitivity when compared with miRanda but they predict several thousand fewer miR-mRNA interactions. Another survey that investigated several traditional predictors including PicTar, TargetScanS, miRanda, and RNAhybrid [8], concludes that miRanda and RNAhybrid obtain lower accuracy and sensitivity when compared with TargetScanS and PicTar [17]. These conclusions were also confirmed in a recent study by Huang and colleagues [16]. They show that the highest quality predictions are obtained by TargetScanS, closely followed by PicTar, while miRanda and DIANA-microT were ranked lower. Most recently, Kuhn and colleagues suggest use of PictTar, TargetScanS, and PicTar to perform computational prediction of miR targets [20]. Based on the above, our experimental section includes three representative miR target prediction methods, TargetScanS, PicTar, and Diana-MicroT. The first two were selected based on their favorable performance, while predictions of Diana-MicroT were used as a point of reference, i.e., representative early generation program characterized by a relatively low sensitivity.
Recent research resulted in development of several ML methods. These methods usually filter predictions provided by the traditional predictors. Their main drawback is that they filter targets by using a predefined and relatively small number of false targets, i.e., they do not scan the mRNA sequence but instead they simulate that by using a small set of negatives (false targets). For instance, a method by Yan and colleagues filters miRanda's predictions based on 48 positive and 16 negative sites [21]. A more recent, NBmiRTar method, which also filters predictions of miRanda, applies 225 true miR targets, 38 confirmed negative sites, and up to 5000 of artificially generated negative sites [22]. The most recent method is based on binding matrix technique, in which the information concerning both the miRNA sequence and a set of experimentally validated targets is used to perform predictions [23]. The main drawback of this approach is the necessity of providing a set of validated targets which is not required in case of the proposed and the abovementioned sequence-based prediction methods. At the same time, we note that the ML methods establish the prediction model based on information concerning both the seed and the non-seed positions, which is also exploited in our research. To this end, we include NBmiRTar method in our experimental section.
Datasets.
Dataset name | Dataset details | Dataset goal |
---|---|---|
Design set | 66 human miR-mRNA duplexes published in TarBase before 2006, see Table 1 [see Additional file 1]; this set includes 29 miRs and 36 genes. | Design of the proposed prediction method. Evaluation and comparison of sensitivity – # predictions trade-off and overlap between the predictions of different methods. |
Independent set | 39 human miRs that were published in TarBase between January 2006 and June 2007; this set includes 20 miRs and 26 genes. | Evaluation and comparison of sensitivity – # predictions trade-off and overlap between the predictions of different methods. |
Interactions set | 190 miR-mRNA interactions pairs experimentally tested in Drophilia. The dataset was taken from [30]. | Evaluation and comparison of the specificity/sensitivity trade-off. The ROC curves and AUC values were compared with results of five competing methods reported in [30]. |
GO (glioblastoma oncogenes) set | Ten glioblastoma oncogenes. The choice is motivated by our expertise in profiling glioblastoma. Although our goal was to compare predictions on 17 glioblastoma oncogenes, only 10 of them could be found in PicTar database. The oncogenes and the associated 328 miRs are given in Tables 2 and 3 [see Additional file 1], respectively. | Comparison of the sensitivity, number of predicted targets and overlap between the predictions of different methods. |
Results and discussion
Datasets and experimental setup
Dataset used to validate and compare the proposed prediction method are summarized in Table 1. The following empirical tests were performed:
1. Comparison of sensitivity – number of predicted targets trade-off. HuMiTar, PicTar, DIANA-MicroT, TargetScanS, and NBmiRTar were compared on the design set and the independent set.
2. Comparison of the overlap of predictions. HuMiTar is compared with the best-performing traditional method PicTar and TargetScanS and ML method NBmiRTar on the GO set. We also include Western blots which are used to verify correctness of some of the HuMiTar predictions.
3. Evaluation and comparison of sensitivity/specificity trade-off based on ROC (receiver operating characteristic) analysis. The predictions of HuMiTar on the interactions set were compared with results of five competing predictions methods reported in [30].
4. Predictions on p53. HuMiTar predicted a set of miRs that target p53, some of which were independently verified in [31].
5. Genome-wide target prediction. HuMiTar was applied to perform genome wide predictions for 16 different species.
We also estimate the signal-to-noise ratio based on predictions for human 3' UTRs and using the procedure introduced to validate PicTar [10]. This ratio and the analysis presented in test 1 are performed to estimate specificity of the proposed method; similar evaluation for the traditional methods was done in [3, 10, 11, 17, 18]. Finally, we also estimate the computational complexity of the proposed method.
Comparison of sensitivity – number of predicted targets trade-off
Summary of prediction results of HuMiTar (HT), PicTar (PT), DIANA-MicroT (DM), TargetScanS (TS), and NBmiRTar (NT). Panel A gives results for the design set of 66 miR-mRNA duplexes. Panel B gives results for the independent set of 39 miR-mRNA duplexes. The hollow circles show results of individual methods, hollow triangles show results of union between HuMiTar and one competing method, and cross corresponds to union of the four competing methods.
Following, we concentrate on the results on the independent set since this set was not used to design the proposed method and thus it allows for an unbiased analysis. A recent study by Nielsen and colleagues reveals several miR targeting determinants [32]. They concern patterns outside of the seed and include presence of adenosine opposite miR base 1 and of adenosine or uridine opposite miR base 9. We applied both of these determinants on the set of 39 duplexes, and found 10 matches, i.e., 10 duplexes satisfy both of the determinants. All of these 10 duplexes were correctly predicted by HuMiTar, while 8 were predicted by TargetScanS, 6 by PicTar, 5 by NBmiRTar, and none by DIANA-MicroT. When considering 13 out of 39 cases for which the adenosine or uridine was opposite miR base 9, all of them were correctly classified by HuMiTar, and 10, 8, 6, and 0 by TargerScanS, PicTar, NBmiRTar, and DIANA-MicroT, respectively. Finally, for 19 duplexes in which the adenosine was opposite miR base 1, 16 of them were found by HuMiTar, 12 by TargerScanS, 10 by PicTar, 7 by NBmiRTar, and none by DIANA-MicroT. This provides an independent validation of the improvements provided by the HuMiTar, which uses scoring function that considers base-pairing outside of the seed, in contrast to the traditional methods that are based on the base-pairing rules only in the seed region.
Comparison of the overlap of predictions
Comparison of the number of predicted targets and their overlap for predictions on the GO set. Panel A shows the number of targets predicted by the three competing methods including PicTar, TargetScanS, and NBmiRTar. The white area represents targets that overlap with predictions of HuMiTar and the black area shows the remaining targets. Panel B shows the number of targets predicted by HuMiTar. The white area shows overlap with predictions of a competing method indicated at the x -axis, and the black area shows predictions specific to HuMiTar. In the case of PicTar, the predictions are reduced to a set of miRs that are included in the PicTar's database http://pictar.mdc-berlin.de/.
HuMiTar predicts 97% of the targets that were predicted by PicTar while only 4 targets were predicted exclusively by PicTar. At the same time, HuMiTar finds numerous extra targets that were not predicted by PicTar. Among them, 646 and 442 extra targets were predicted for miRs that are and that are not included in the PicTar's database, respectively. Although the results show that HiMiTar generates larger number of predictions that cover virtually all predictions made by PicTar, these additional prediction could constitute either true or false positives. Since it would infeasible to verify correctness of the entire set of 1088 additional predictions provided by HuMiTar, we concentrate our efforts on a specific target, Septin7, due to our existing wet-lab expertise. The test considers three sets of miRs:
Set 1
5 miRs from the 18 targets that were predicted by both HuMiTar and PicTar, i.e. miR-19a, miR-127, miR-141, miR-182, and miR183. This set of used to investigate whether the common results in fact concern true positives.
Set 2
5 miRs from the 23 targets that were predicted only by HuMiTar and which are not included in the PicTar's database, i.e. miR-202, miR-248, miR-412, miR-453, and miR-450. This set concerns miRs that have not been predicted by PicTar, i.e. they were predicted only with the use of HuMiTar.
Set 3
11 miRs from the 34 extra targets on Septin7, which were found only by HuMiTar although these miRs are included in PicTar's database, i.e. miR-148, miR-106b, miR-134, miR-106, miR-144, miR-151, miR-384, miR-101, miR-142, miR-129, and miR-126. This set is of particular interest (and thus it is larger), since it concerns targets that were not predicted by PicTar, but which were predicted by HuMiTar.
We run a Western blot according to the following procedure. The human glioblastoma cell line U251 was obtained from China Academia Sinica cell repository in Shanghai, China. All cell lines were grown in Dulbecco's modified Eagle's medium (DMEM) (Gibco, USA) supplemented with 10% fetal bovine serum (Gibco, USA), 2 mM glutamine (Sigma, USA), 100 units of penicillin/ml (Sigma, USA), and 100 μg of streptomycin/ml (Sigma, USA), incubated at 37°C with 5% CO2, and sub-cultured every 2~3 days. The antisense oligonucleotides of the pre-scanned miRNAs were chemically synthesized by GenePharma (Shanghai, China) and were transfected into U251 cells by Oligofectamine (Invitrogen, USA) according to the manufactures' protocol. Parental and transfected cells were washed with ice-cold phosphate-buffered saline (PBS) three times. The cells were then solubilized in 1% Nonidet P-40 lysis buffer (20 mM Tris, pH 8.0, 137 mM NaCl, 1% Nonidet P-40, 10% glycerol, 1 mM CaCl2, 1 mM MgCl2, 1 mM phenylmethylsulfonyl fluoride, 1 mM sodium fluoride, 1 mM sodium orthovanadate, and a protease inhibitor mixture). Homogenates were clarified by centrifugation at 20,000 ×g for 15 minutes at 4°C and protein concentrations were determined by a bicinchoninic acid protein assay kit (Pierce, USA). Equal amounts of lysates were subjected to SDS-PAGE on 8% SDS-acrylamide gel. Separate proteins were transferred to PVDF membranes (Millipore, USA) and incubated with primary antibody against Septin-7 (Santa Cruz, USA), followed by incubation with HRP-conjugated secondary antibody (Zymed, USA). The specific protein was detected by using a SuperSignal protein detection kit (Pierce, USA). The membrane was stripped and re-probed with antibody against β-actin (Santa Cruz, USA).
Western blots for selected 21 miRs and Septin7. The Septin7 expression levels were measured (left to right) for (1) control sample, (2) miR-127, (3) miR-182, (4) miR-412, (5) miR-19a, (6) miR-453, (7) miR-448, (8) miR-450, (9) miR-183, (10) miR-141, (11) miR-202, (12) miR-148, (13) miR-106b, (14) miR-134, (15) miR-106, (16) miR-144, (17) miR-151, (18) miR-384, (19) miR-101, (20) miR-142, (21) miR-129 and (22) miR-126. The position 1 is the control sample; positions 2 to 11 inclusive concern 10 miRs for which predictions were obtained either by both PicTar and HiMiTar (positions 2, 3, 5, 9, and 10) or only by HuMiTar while these MiR were not included in the PicTar's database (positions 4, 6, 7, 8, and 11); positions 12 to 22 concern 11 miRs which were predicted by HuMiTar and which were included in the PicTar's database. We note that our analysis lacks results on the mutant targets that would strengthen the claim that the activation results of up-regulation of the predicted miRs. Due to limited resources and since the goal of this work is to present a new in-silico prediction method rather than to investigate whether miRs can up-regulate translation, we note that our conclusions concerning the up-regulation should not be considered as the primary outcome of this work.
-
In Set 1, up-regulation is shown for miR-19a (position 5), miR-183 (position 9), and miR-141 (position 10); we also observe that miR-182 (position 3) is likely to be a true positive. The experiment shows that miR-127 does not affect the expression levels of Septin7, which may suggest that this is a false positive prediction. To summarize, 4(3) out of 5 of the predictions in Set 1 are shown to be true positives.
-
In Set 2, miR-448 (position 7), miR-450 (position 8), and miR-202 (position 11) show up-regulation; miR-453 (position 6) is a borderline case, although the Western blot suggests that it could be classified as a true positive. Finally, miR-412 has no impact on the Septin7 expression, and thus it should be considered as a false positive. As a result, 4(3) out of 5 predictions in this set are true positives.
-
In Set 3, the Western blots indicate that all 11 miRs (positions 12 to 22 in Figure 3) target Septin7 and thus they constitute true positive predictions.
Overall, the experiment indicates that for the considered set of miRs, HuMiTar obtains about 80% sensitivity when predicting targets on Septin7. The reported up-regulation is consistent with recent research that also indicates that miRs can up-regulate the translation [33, 34]. Although the above results cannot be generalized to other targets, they indicate that predictions generated by the proposed method are characterized by favorable sensitivity when compared with PicTar.
Among the TargetScanS predictions, 91% were also predicted by HuMiTar and the remaining 9%, i.e., 68 targets, were not included in the output of HuMiTar. At the same time, the proposed method provides 562 additional predictions. Similarly as in case of comparison with PicTar, we probe the sensitivity of both prediction methods based on targets predicted for Septin7. Analysis of 39 miRs that were predicted exclusively by HuMiTar shows that 11 of them (miR-101, miR-126, miR-129, miR-134, miR-144, miR-151, miR-202, miR-384, miR-412, miR-450, miR-453) are included in the Western blot on Figure 3. Among them, nine are true positives, miR-453 is a borderline case, and miR-412 is a false positive. Although our limited resources prohibit more extensive experimental analysis, the above analysis suggests that additional predictions provided by HuMiTar include true positives.
Finally, although the overlap between the predictions of HuMiTar and NBmiRTar is the smallest among the three competing methods, it still constitutes over a half (56%) of the NBmiRTar's predictions. The HuMiTar provides 949 predictions which are not included in the output of NBmiRTar, while 213 predictions are exclusive to NBmiRTar.
Overall, the test on the GO set shows that predictions of HuMiTar overlap with the predictions of the competing methods. In particular, we note that the HuMiTar's outputs cover almost all predictions of PicTar and majority of predictions of the other two methods. At the same time, our predictions also include novel targets that could correspond to biologically meaningful sites.
Evaluation and comparison of sensitivity/specificity trade-off
ROC curves and the corresponding AUC scores that quantify sensitivity and specificity of different miR target predictors on the interactions set. The results include five existing prediction methods, PicTar [10], miRanda [19], method by Stark and colleagues [35], STarMir [36], and PITA [30]. The latter method includes two versions, one requiring unpairing of only target-site nucleotides (PITA no flank) and another that also requires unpairing of 3 and 15 flanking nucleotides upstream and downstream of the target site (PITA 3/15 flank), respectively. For each prediction method, the targets were sorted by score and the FP rates (x axis) and TP rates (y axis) were plotted for each possible score prediction threshold. The area under the curve (AUC) for each method is shown in the figure legend. The AUC is computed by extending each plot to the upper right corner as in [30]. The results obtained by a random sorting of the targets are shown using a thin dashed line. The ROC curves and AUC values of PITA, STARK, PicTar, miRanda, and STarMir were taken from [30].
Predictions on p53
List of 15 miRs that were predicted by HuMiTar to target p53 and which were confirmed by Xi and colleagues (Xi et al., 2006).
hsa-let-7a | hsa-miR-296 | hsa-miR-125b | hsa-miR-183 | hsa-miR-19b |
---|---|---|---|---|
hsa-miR-30b | hsa-miR-30c | hsa-miR-30a-5p | hsa-miR-30d | hsa-miR-27a |
hsa-miR-103 | hsa-miR-107 | hsa-miR-92 | hsa-miR-10a | hsa-miR-326 |
Genome-wide target prediction
Summary of genome-wide predictions with HuMiTar.
Species type | # of miRs used | Total # of targets predicted with HuMitar | Average number of targets per miR |
---|---|---|---|
Anopheles gambiae | 37 | 324050 | 8758 |
Bos taurus | 125 | 1234929 | 9879 |
Caenorhabditis briggsae | 93 | 759339 | 8165 |
Caenorhabditis elegan | 134 | 1035742 | 7729 |
Drosophila melanogaster | 78 | 728121 | 9335 |
Drosophila pseudoobscura | 68 | 653553 | 9611 |
Fugu rubripes | 109 | 1083787 | 9943 |
Gallus gallus | 133 | 1302995 | 9797 |
Homo sapiens | 471 | 4591244 | 9748 |
Macaca mulatta | 70 | 738781 | 10554 |
Monodelphis domestica | 100 | 1038582 | 10386 |
Mus musculus | 380 | 3470973 | 9134 |
Pan troglodytes | 80 | 843686 | 10546 |
Rattus norvegicus | 238 | 2379080 | 9996 |
Tetraodon nigroviridis | 109 | 1090377 | 10003 |
Xenopus tropicalis | 160 | 1634721 | 10217 |
average | 149 | 1431873 | 9613 |
One of the accepted ways of assessing the statistical significance of predicted targets, which was performed in [3, 10, 11], is based on using random miR sequences ('mock' miRs) as controls [17]. The motivation is that the mock sequences are unlikely to be biologically relevant and thus observing the ratio between the number of predictions for real miRs and for the mock miRs would indicate how many of the predictions for real miRs are indeed biologically relevant. The ratio of 'real' versus 'mock' predictions is provided as an estimate of the signal-to-noise ratio (SNR) of the target predictions [17]. The HuMiTar's SNR was evaluated using the set of 58 miRs and the randomization procedure to generate the mock miRs that were originally used to estimate SNR for PicTar [10]. The SNR of HuMiTar for the 58 real/mock miRs in the entire human 3'UTR set equals 1.99. PicTar's SNR was estimated to equal 1.8 when considering conservation using human, chimpanzee and mouse genomes, and 2.3 and 3.6 when dog and chicken genomes were added, respectively [10]. In case of TargetScanS the ratio was estimated to be 2.4 and 3.8 when considering all predictions and when considering only the positions conserved in all five genomes, respectively [11]. We note that the latter improved SNR was also accompanied by a 51% loss in sensitivity [11]. The SNR of the proposed method, which does not implement cross-species conservation, is comparable with the ratio of PicTar that was computed when the cross-species conservation was limited to human, chimpanzee and mouse genomes, and to the SNR of TargetScanS when conservation was not included. We anticipate that the SNR would increase if we would incorporate the cross-species conservation to filter our predictions.
Computational complexity
The asymptotic computational complexity of HuMiTar is O(m2n) where m is the length of miR sequence and n is the length of the mRNA. This dominant factor contributing to the overall complexity is computation of 2D-coding that is used to perform initial screening of the mRNA. We note that since m is a small constant, the proposed method is characterized by a linear complexity with respect to the length of the mRNA sequence. We also performed experimental evaluation of the execution time. Using the set of 39,215 human 3'UTR sequences (the total lengths of these sequences equals 3.62e+07) the search for miR-21 targets takes 1,488 seconds using a desktop computer. We also computed execution times for ten randomly drawn miRs shown in Table 9 [see Additional file 1]. The targets were predicted on average in 1,451 seconds (about 24 minutes) per one miR, which shows that the proposed method can be applied on the genomic scale.
Conclusion
The prediction of animal miR targets is an open and difficult problem in spite of several years of the existing research. HuMiTar, which is a prediction method designed based on human miR targets, is shown to provide predictions characterized by favorable sensitivity, which comes at a price of an increased number of predictions. The HuMiTar's predictions cover majority of the predictions of the best-performing competing methods such as PicTar, TargetScanS, and NBmiRTar. HuMiTar has good computational efficiency and comparable signal-to-noise ratio when compared with TargetScanS and PicTar. ROC analysis shows that HuMiTar provides predictions of quality that is comparable with the quality of PicTar, while our predictions are characterized by high true positive rates and moderate values of false positive rates. Our prediction method constitutes a step towards providing an efficient computational model for studying translational gene regulation by microRNAs. Our future work will concentrate of relaxation of the base-pairing requirements in the seed region to accommodate for miR-mRNA with the imperfect pairing and inclusion of information based on the stacked pairs and unpaired regions.
Methods
HuMiTar works in two steps: (1) a 2D-coding method finds candidate targets by scanning 3'UTR of a given mRNA; and (2) the selected candidate targets are filtered using a composite scoring function.
Scoring function
The scoring function was designed based on statistical analysis of 66 miR-mRNAs duplexes from the design set. Based on analysis of several recent studies [2–14, 24], the miR sequence was divided into four regions: (1) position 1; (2) seed region (positions 2 to 8); (3) region 1 (positions 9 to 13); and (4) region 2 (positions 14 to 20). We computed conditional (assuming that they concern only the actual sites) and unconditional (assuming any other position in the mRNA using a sliding window) frequencies of the nucleotide pairs formed between miR's seed region, region 1, and region 2 and the corresponding 3' UTR of the mRNA. These values were used to compute affinity of each pair to form a bond between miR and mRNA; the affinity values are shown in Table 11 [see Additional file 1]. The differences in the obtained affinity values of the same pairs for different regions show that the occurrence of the corresponding pairs in a considered candidate complex should be weighed accordingly. We apply the principles of balance of moments to compute the weight values. The underlying interpretation is that the high affinity to bind in the seed region between a given miR and mRNA fragment should be balanced by sufficiently large affinity to bind in the non-seed regions (regions 1 and 2). We assume that the sum of moments generated by positions in the seed region should be greater than the sum of moment of the positions in non-seed regions. This problem is formulated and solved, i.e., the corresponding scoring function that optimizes the balance between binding in the seed and the non-seed regions is parameterized, using a standard linear programming model. The parameterization shows that the formation of complementary pairs for positions 9, 10, 12, 15, 16, 19 and 20 is less "important" than for the positions 11, 13, 14, 17 and 18. This result can be validated against a recent study that investigated Watson-Crick pairing for contiguous nucleotides. This study concluded that when excluding the seed, positions 13–16 have the strongest preference for the complementary pairing [37]. Although we consider each position individually, while the other study analyzed multimers, we also found that positions 13 and 14 have strong tendency to form complementary pairs. Our prediction model considers existence of two opposing forces: positive effect corresponding to the formation of complementary base pairs (which is quantified by the reward function), and negative effect due to the existence of non-complementary base pairs (quantified by the penalty function). The scoring function is defined as a difference between the values produced by the reward and the cost functions.
2D-coding method
The 3' UTR of mRNA is scanned using a sliding window. The basic principle of the 2D-coding is to scan an mRNA segment starting at a given position (using the sliding window) by finding stretches (segments) of complementary base pairs, which are denoted by A i where i = 1, 2, ..., 5. We start with finding the first segment, denoted by A1, in the miR's seed region, and then continue along the miR's sequence to find the subsequent segments. Given that a sufficient number of such segments can be found, which depends on several parameters explained in the [Additional file 1], the 2D-coding will output a candidate complex. The procedure will stop after finding A5 since no more then five complementary segments were found when scanning duplexes in the design set.
Declarations
Acknowledgements
J.R. thanks NSFC (No 10671100), Liuhui Center for applied mathematics, the joint program of Tianjin and Nankai Universities, and China-Canada interchange project administered by MITACS. This work was supported in part by the Discovery grant to L.K. from NSERC Canada. K.C. was supported by a scholarship sponsored by Alberta Ingenuity and iCORE. C.K. and P.P. would like also to thank NSFC (No 10671100).
The authors would like to thank Michael Kertesz for providing dataset and experimental results for the ROC analysis.
Authors’ Affiliations
References
- Engels BM, Hutvagner G: Principles and effects of microRNA-mediated post-transcriptional gene regulation. Oncogene. 2006, 25: 6163-6169.PubMedView ArticleGoogle Scholar
- Enright AJ, John B, Gaul U, Tuschl T, Sander C, Marks DS: MicroRNA targets in Drosophila. Genome Biol. 2003, 5 (1): R1-PubMedPubMed CentralView ArticleGoogle Scholar
- Lewis BP, Shih I, Jones-Rhoades MW, Bartel DP, Burge CB: Prediction of mammalian microRNA targets. Cell. 2003, 115: 787-798.PubMedView ArticleGoogle Scholar
- Stark A, Brennecke J, Russell RB, Cohen SM: Identification of Drosophila MicroRNA targets. PLoS Biol. 2003, 1 (3): e397-10.1371/journal.pbio.0000060.View ArticleGoogle Scholar
- Doench JG, Sharp PA: Specificity of microRNA target selection in translational repression. Genes Dev. 2004, 18 (5): 504-511.PubMedPubMed CentralView ArticleGoogle Scholar
- John B, Enright AJ, Aravin A, Tuschl T, Sander C, Marks DS: Human microRNA targets. PLoS Biol. 2004, 2 (11): e363-PubMedPubMed CentralView ArticleGoogle Scholar
- Kiriakidou M, Nelson P, Kouranov A, Fitziev P, Bouyioukos C, Mourelatos Z, Hatzigeorgiou AG: A combined computational- experimental approach predicts human miR targets. Genes & Dev. 2004, 18: 1165-1178.View ArticleGoogle Scholar
- Rehmsmeier M, Steffen P, Hochsmann M, Giegerich R: Fast and effective prediction of microRNA/target duplexes. RNA. 2004, 10: 1507-1517.PubMedPubMed CentralView ArticleGoogle Scholar
- Vella MC, Reinert K, Slack FJ: Architecture of a validated microRNA: target interaction. Chem Biol. 2004, 11: 1619-1623.PubMedView ArticleGoogle Scholar
- Krek A, Grun D, Poy MN, Wolf R, Rosenberg L, Epstein EJ, MacMenamin P, Piedade I, Gunsalus KC, Stoffel M: Combinatorial microRNA target predictions. Nat Genet. 2005, 37: 495-500.PubMedView ArticleGoogle Scholar
- Lewis BP, Burge CB, Bartel DP: Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005, 120: 15-20.PubMedView ArticleGoogle Scholar
- Saetrom O, Snove O, Saetrom P: Weighted sequence motifs as an improved seeding step in microRNA target prediction algorithms. RNA. 2005, 11 (7): 995-1003.PubMedPubMed CentralView ArticleGoogle Scholar
- Xie X, Lu J, Kulbokas EJ, Golub TR, Mootha V, Lindblad-Toh K, Lander ES, Kellis M: Systematic discovery of regulatory motifs in human promoters and 3'UTRs by comparison of several mammals. Nature. 2005, 434: 338-345.PubMedPubMed CentralView ArticleGoogle Scholar
- Watanabe Y, Yachie N, Numata K, Saito R, Kanai A, Tomita M: Computational analysis of microRNA targets in Caenorhabditis elegans. Gene. 2006, 365: 2-10.PubMedView ArticleGoogle Scholar
- Brennecke J, Stark A, Russell RB, Cohen SM: Principles of MicroRNA-Target Recognition. PLoS Biol. 2005, 3 (3): e85-PubMedPubMed CentralView ArticleGoogle Scholar
- Huang JC, Babak T, Corson TW, Chua G, Khan S, Gallie BL, Hughes TR, Blencowe BJ, Frey BJ, Morris QD: Using expression profiling data to identify human microRNA targets. Nat Methods. 2007, 4 (12): 1045-9.PubMedView ArticleGoogle Scholar
- Rajewsky N: microRNA target predictions in animals. Nat Genet. 2006, 38 (Suppl): S8-S13.PubMedView ArticleGoogle Scholar
- Sethupathy P, Megraw M, Hatzigeorgiou AG: A guide through present computational approaches for the identification of mammalian microRNA targets. Nat Methods. 2006, 3: 881-886.PubMedView ArticleGoogle Scholar
- Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ: miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 2006, 34: D140-D144.PubMedPubMed CentralView ArticleGoogle Scholar
- Kuhn DE, Martin MM, Feldman DS, Terry AV, Nuovo GJ, Elton TS: Experimental validation of miRNA targets. Methods. 2008, 44 (1): 47-54.PubMedPubMed CentralView ArticleGoogle Scholar
- Yan X, Chaoa T, Tub K, Zhanga Y, Xieb L, Gonga Y, Yuana J, Qianga B, Peng X: Improving the prediction of human microRNA target genes by using ensemble algorithm. FEBS Lett. 2007, 581 (8): 1587-1593.PubMedView ArticleGoogle Scholar
- Yousef M, Jung S, Kossenkov AV, Showe LC, Showe MK: Naïve Bayes for microRNA target predictions – machine learning for microRNA targets. Bioinformatics. 2007, 23 (22): 2987-2992.PubMedView ArticleGoogle Scholar
- Moxon S, Moulton V, Kim JT: A scoring matrix approach to detecting miRNA target sites. Alg Mol Biol. 2008, 3: 3-10.1186/1748-7188-3-3.View ArticleGoogle Scholar
- Didiano D, Hobert O: Perfect seed pairing is not a generally reliable predictor for miRNA-target interactions. Nat Struct Mol Biol. 2006, 13 (9): 849-51.PubMedView ArticleGoogle Scholar
- Sethupathy P, Corda B, Hatziegeorgiou AG: TarBase: A comprehensive database of experimentally supported animal microRNA targets. RNA. 2006, 12: 192-197.PubMedPubMed CentralView ArticleGoogle Scholar
- Rhoades MW, Reinhart BJ, Lim LP, Burge CB, Bartel B: Prediction of plant microRNA targets. Cell. 2002, 110: 513-520.PubMedView ArticleGoogle Scholar
- Zhang Y: miRU: an automated plant miRNA target prediction server. Nucleic Acids Res. 2005, 33: W701-W704.PubMedPubMed CentralView ArticleGoogle Scholar
- Gusev Y: Computational methods for analysis of cellular functions and pathways collectively targeted by differentially expressed microRNA. Methods. 2008, 44 (1): 61-72.PubMedView ArticleGoogle Scholar
- Adams BD, Furneaux H, White B: The Micro-Ribonucleic Acid (miRNA) miR-206 Targets the Human Estrogen Receptor-α (ERα) and Represses ERα Messenger RNA and Protein Expression in Breast Cancer Cell Lines. Mol Endocrinol. 2007, 21 (5): 1132-1147.PubMedView ArticleGoogle Scholar
- Kertesz M, Iovino N, Unnerstall U, Gaul U, Segal E: The role of site accessibility in microRNA target recognition. Nat Genetics. 2007, 39: 1278-84. 10.1038/ng2135.PubMedView ArticleGoogle Scholar
- Xi Y, Shalgi R, Fodstad O, Pilpel Y, Ju J: Differentially Regulated Micro-RNAs and Actively Translated Messenger RNA Transcripts by Tumor suppressor p53 in Colon Cancer. Clin Cancer Res. 2006, 12 (7 Pt 1): 2014-2024.PubMedView ArticleGoogle Scholar
- Nielsen CB, Shomron N, Sandberg R, Hornstein E, Kitzman J, Burge CB: Determinants of targeting by endogenous and exogenous microRNAs and siRNAs. RNA. 2007, 13: 1894-1910.PubMedPubMed CentralView ArticleGoogle Scholar
- Vasudevan S, Tong Y, Steitz JA: Switching from repression to activation: microRNAs can up-regulate translation. Science. 2007, 318 (5858): 1931-4.PubMedView ArticleGoogle Scholar
- Rusk N: When microRNAs activate translation. Nature Methods. 2008, 5: 1223-Google Scholar
- Stark A, Brennecke J, Bushati N, Russell RB, Cohen SM: Animal MicroRNAs confer robustness to gene expression and have a significant impact on 3'UTR evolution. Cell. 2005, 123: 1133-46.PubMedView ArticleGoogle Scholar
- Long D, Lee R, Williams P, Chan CY, Ambros V, Ding Y: Potent effect of target structure on microRNA function. Nat Struct Mol Biol. 2007, 14: 287-94.PubMedView ArticleGoogle Scholar
- Grimson A, Kai-How Farth K, Johnston WK, Garrnet-Engele P, Lim LP, Bartel DP: MicroRNA targeting specificity in mammals: determinants beyond seed pairing. Molecular Cell. 2007, 27: 91-105.PubMedPubMed CentralView ArticleGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.