Computational methods for analysis of cellular functions and pathways collectively targeted by differentially expressed microRNA.

This report presents computational methods of analysis of cellular processes, functions, and pathways affected by differentially expressed microRNA, a statistical basis of the gene enrichment analysis method, a modification of enrichment analysis method accounting for combinatorial targeting of Gene Ontology categories by multiple miRNAs and examples of the global functional profiling of predicted targets of differentially expressed miRNAs in cancer. We have also summarized an application of Ingenuity Pathway Analysis tools for in depth analysis of microRNA target sets that may be useful for the biological interpretation of microRNA profiling data. To illustrate the utility of these methods, we report the main results of our recent computational analysis of five published datasets of aberrantly expressed microRNAs in five human cancers (pancreatic cancer, breast cancer, colon cancer, lung cancer, and lymphoma). Using a combinatorial target prediction algorithm and statistical enrichment analysis, we have determined Gene Ontology categories as well as biological functions, disease categories, toxicological categories, and signaling pathways that are: targeted by multiple microRNAs; statistically significantly enriched with target genes; and known to be affected in specific cancers. Our recent computational analysis of predicted targets of co-expressed miRNAs in five human cancers suggests that co-expressed miRNAs provide systemic compensatory response to the abnormal phenotypic changes in cancer cells by targeting a broad range of functional categories and signaling pathways reportedly affected in a particular cancer.     

Refining microRNA target predictions: Sorting the wheat from the chaff

microRNAs are short RNAs that reduce gene expression by binding to their targets. The accurate prediction of microRNA targets is essential to understanding the function of microRNAs. Computational predictions indicate that all human genes may be regulated by microRNAs, with each microRNA possibly targeting thousands of genes. Here we discuss computational methods for identifying mammalian microRNA targets and refining them for further experimental validation. We describe microRNA target prediction resources and procedures and how they integrate with various types of experimental techniques that aim to validate them or further explore their function. We also provide a list of target prediction databases and explain how these are curated.     

Prediction of personalized microRNA activity

MicroRNAs have been known to regulate almost all physiological and pathological processes by suppressing their target genes. In humans, more than 1000 microRNAs have been identified, each of which targets dozens or even hundreds of genes. Facing this huge repertoire of microRNA targeting, it is important to identify which microRNAs are active, i.e., down-regulating their targets, in specific physiological or pathological conditions. Predicting active microRNAs is different from predicting microRNA targets because the authentic target genes of a microRNA are often not directly and solely regulated by that microRNA, leading to inconsistent expression changes between the microRNA and its true targets. Several computational programs have been proposed to predict the activity of a microRNA from the expressions of its target genes. These programs performed well when being applied on the expression data obtained from distinct tissue types or from experiments that transfect a microRNA into cells (i.e., non-physiological). But the performance of microRNA activity prediction is not clear on the expression data from the same tissue type in two physiological conditions, e.g., liver tissues from cancer patients and healthy people. In this work, we evaluate the performance of two microRNA activity prediction programs using seven expression data sets, all of which compare samples in two physiological conditions, as well as propose a new approach that predicts microRNA activity with an accuracy of over 80%. Unlike current methods, which predict active microRNAs by comparing two groups of samples, e.g., tumor versus normal, our new approach compares each diseased sample with all the samples in the control group. In other words, it can predict the microRNA activity of a person. In this work, this new application is named to predict “personalized microRNA activity”.     

Assessing the utility of thermodynamic features for microRNA target prediction under relaxed seed and no conservation requirements

Background

Many computational microRNA target prediction tools are focused on several key features, including complementarity to 5′seed of miRNAs and evolutionary conservation. While these features allow for successful target identification, not all miRNA target sites are conserved and adhere to canonical seed complementarity. Several studies have propagated the use of energy features of mRNA:miRNA duplexes as an alternative feature. However, different independent evaluations reported conflicting results on the reliability of energy-based predictions. Here, we reassess the usefulness of energy features for mammalian target prediction, aiming to relax or eliminate the need for perfect seed matches and conservation requirement.

Methodology/Principal Findings

We detect significant differences of energy features at experimentally supported human miRNA target sites and at genome-wide sites of AGO protein interaction. This trend is confirmed on datasets that assay the effect of miRNAs on mRNA and protein expression changes, and a simple linear regression model leads to significant correlation of predicted versus observed expression change. Compared to 6-mer seed matches as baseline, application of our energy-based model leads to ∼3–5-fold enrichment on highly down-regulated targets, and allows for prediction of strictly imperfect targets with enrichment above baseline.

Conclusions/Significance

In conclusion, our results indicate significant promise for energy-based miRNA target prediction that includes a broader range of targets without having to use conservation or impose stringent seed match rules.     

A scoring matrix approach to detecting miRNA target sites

Background  

Experimental identification of microRNA (miRNA) targets is a difficult and time consuming process. As a consequence several computational prediction methods have been devised in order to predict targets for follow up experimental validation. Current computational target prediction methods use only the miRNA sequence as input. With an increasing number of experimentally validated targets becoming available, utilising this additional information in the search for further targets may help to improve the specificity of computational methods for target site prediction.     

Flanking region sequence information to refine microRNA target predictions

The non-coding elements of a genome, with many of them considered as junk earlier, have now started gaining long due respectability, with microRNAs as the best current example. MicroRNAs bind preferentially to the 3′ untranslated regions (UTRs) of the target genes and negatively regulate their expression most of the time. Several microRNA:target prediction softwares have been developed based upon various assumptions and the majority of them consider the free energy of binding of a target to its microRNA and seed conservation. However, the average concordance between the predictions made by these softwares is limited and compounded by a large number of false-positive results. In this study, we describe a methodology developed by us to refine microRNA:target prediction by target prediction softwares through observations made from a comprehensive study. We incorporated the information obtained from dinucleotide content variation patterns recorded for flanking regions around the target sites using support vector machines (SVMs) trained over two different major sources of experimental data, besides other sources. We assessed the performance of our methodology with rigorous tests over four different dataset models and also compared it with a recently published refinement tool, MirTif. Our methodology attained a higher average accuracy of 0.88, average sensitivity and specificity of 0.81 and 0.94, respectively, and areas under the curves (AUCs) for all the four models scored above 0.9, suggesting better performance by our methodology and a possible role of flanking regions in microRNA targeting control. We used our methodology over genes of three different pathways — toll-like receptor (TLR), apoptosis and insulin — to finally predict the most probable targets. We also investigated their possible regulatory associations, and identified a hsa-miR-23a regulatory module.     

Systematic identification of microRNA functions by combining target prediction and expression profiling

Target predictions and validations are major obstacles facing microRNA (miRNA) researchers. Animal miRNA target prediction is challenging because of limited miRNA sequence complementarity to the targets. In addition, only a small number of predicted targets have been experimentally validated and the miRNA mechanism is poorly understood. Here we present a novel algorithm for animal miRNA target prediction. The algorithm combines relevant parameters for miRNA target recognition and heuristically assigns different weights to these parameters according to their relative importance. A score calculation scheme is introduced to reflect the strength of each parameter. We also performed microarray time course experiments to identify downregulated genes due to miRNA overexpression. The computational target prediction is combined with the miRNA transfection experiment to systematically identify the gene targets of human miR-124. miR-124 overexpression led to a significant downregulation of many cell cycle related genes. This may be the result of direct suppression of a few cell growth inhibitors at the early stage of miRNA overexpression, and these targeted genes were continuously suppressed over a long period of time. Our high-throughput approach can be generalized to globally identify the targets and functions of other miRNAs.     

Interplay between HIV-1 infection and host microRNAs

Using microRNA array analyses of in vitro HIV-1-infected CD4(+) cells, we find that several host microRNAs are significantly up- or downregulated around the time HIV-1 infection peaks in vitro. While microRNA-223 levels were significantly enriched in HIV-1-infected CD4(+)CD8(-) PBMCs, microRNA-29a/b, microRNA-155 and microRNA-21 levels were significantly reduced. Based on the potential for microRNA binding sites in a conserved sequence of the Nef-3'-LTR, several host microRNAs potentially could affect HIV-1 gene expression. Among those microRNAs, the microRNA-29 family has seed complementarity in the HIV-1 3'-UTR, but the potential suppressive effect of microRNA-29 on HIV-1 is severely blocked by the secondary structure of the target region. Our data support a possible regulatory circuit at the peak of HIV-1 replication which involves downregulation of microRNA-29, expression of Nef, the apoptosis of host CD4 cells and upregulation of microRNA-223.     

Specificity of miR-378a-5p targeting rodent fibronectin

One criterion for microRNA identification is based on their conservation across species, and prediction of miRNA targets by empirical approaches using computational analysis relies on the presence of conservative mRNA 3′UTR. Because most miRNA target sites identified are highly conserved across different species, it is not clear whether miRNA targeting is species-specific. To predict miRNA targeting, we aligned all available fibronectin 3′UTRs and observed significant conservation of all 20 species. Twelve miRNAs were predicted to target most fibronectin 3′UTRs, but rodent fibronectin showed potential binding sites specific for five different miRNAs. One of them, the miR-378a-5p, contained a complete matching seed-region for all rodent fibronectin, which could not be found in any other species. We designed experiments to test whether the species-specific targeting possessed biological function and found that expression of miR-378a-5p decreased cancer cell proliferation, migration, and invasion, resulting in inhibition of tumor growth. Silencing fibronectin expression produced similar effects as miR-378a-5p, while transfection with a construct targeting miR-378-5p produced opposite results. Tumor formation assay showed that enhanced expression of fibronectin in the stromal tissues as a background environment suppressed tumor growth, while increased fibronectin expression inside the tumor cells promoted tumor growth. This was likely due to the different signaling direction, either inside-out or outside-in signal. Our results demonstrated that species-specific targeting by miRNA could also exert functional effects. Thus, one layer of regulation has been added to the complex network of miRNA signaling.     

Naïve Bayes for microRNA target predictions--machine learning for microRNA targets

MOTIVATION: Most computational methodologies for miRNA:mRNA target gene prediction use the seed segment of the miRNA and require cross-species sequence conservation in this region of the mRNA target. Methods that do not rely on conservation generate numbers of predictions, which are too large to validate. We describe a target prediction method (NBmiRTar) that does not require sequence conservation, using instead, machine learning by a na?ve Bayes classifier. It generates a model from sequence and miRNA:mRNA duplex information from validated targets and artificially generated negative examples. Both the 'seed' and 'out-seed' segments of the miRNA:mRNA duplex are used for target identification. RESULTS: The application of machine-learning techniques to the features we have used is a useful and general approach for microRNA target gene prediction. Our technique produces fewer false positive predictions and fewer target candidates to be tested. It exhibits higher sensitivity and specificity than algorithms that rely on conserved genomic regions to decrease false positive predictions.     

Predicting new molecular targets for rhein using network pharmacology

Background

Drugs can influence the whole biological system by targeting interaction reactions. The existence of interactions between drugs and network reactions suggests a potential way to discover targets. The in silico prediction of potential interactions between drugs and target proteins is of core importance for the identification of new drugs or novel targets for existing drugs. However, only a tiny portion of drug-targets in current datasets are validated interactions. This motivates the need for developing computational methods that predict true interaction pairs with high accuracy. Currently, network pharmacology has used in identifying potential drug targets to predicting the spread of drug activity and greatly contributed toward the analysis of biological systems on a much larger scale than ever before.

Methods

In this article, we present a computational method to predict targets for rhein by exploring drug-reaction interactions. We have implemented a computational platform that integrates pathway, protein-protein interaction, differentially expressed genome and literature mining data to result in comprehensive networks for drug-target interaction. We used Cytoscape software for prediction rhein-target interactions, to facilitate the drug discovery pipeline.

Results

Results showed that 3 differentially expressed genes confirmed by Cytoscape as the central nodes of the complicated interaction network (99 nodes, 153 edges). Of note, we further observed that the identified targets were found to encompass a variety of biological processes related to immunity, cellular apoptosis, transport, signal transduction, cell growth and proliferation and metabolism.

Conclusions

Our findings demonstrate that network pharmacology can not only speed the wide identification of drug targets but also find new applications for the existing drugs. It also implies the significant contribution of network pharmacology to predict drug targets.     

Weighted sequence motifs as an improved seeding step in microRNA target prediction algorithms

We present a new microRNA target prediction algorithm called TargetBoost, and show that the algorithm is stable and identifies more true targets than do existing algorithms. TargetBoost uses machine learning on a set of validated microRNA targets in lower organisms to create weighted sequence motifs that capture the binding characteristics between microRNAs and their targets. Existing algorithms require candidates to have (1) near-perfect complementarity between microRNAs' 5' end and their targets; (2) relatively high thermodynamic duplex stability; (3) multiple target sites in the target's 3' UTR; and (4) evolutionary conservation of the target between species. Most algorithms use one of the two first requirements in a seeding step, and use the three others as filters to improve the method's specificity. The initial seeding step determines an algorithm's sensitivity and also influences its specificity. As all algorithms may add filters to increase the specificity, we propose that methods should be compared before such filtering. We show that TargetBoost's weighted sequence motif approach is favorable to using both the duplex stability and the sequence complementarity steps. (TargetBoost is available as a Web tool from http://www.interagon.com/demo/.).     

MicroRNA-377 Regulates Mesenchymal Stem Cell-Induced Angiogenesis in Ischemic Hearts by Targeting VEGF

MicroRNAs have been appreciated in various cellular functions, including the regulation of angiogenesis. Mesenchymal-stem-cells (MSCs) transplanted to the MI heart improve cardiac function through paracrine-mediated angiogenesis. However, whether microRNAs regulate MSC induced angiogenesis remains to be clarified. Using microRNA microarray analysis, we identified a microRNA expression profile in hypoxia-treated MSCs and observed that among all dysregulated microRNAs, microRNA-377 was decreased the most significantly. We also validated that vascular endothelial growth factor (VEGF) is a target of microRNA-377 using dual-luciferase reporter assay and Western-blotting. Knockdown of endogenous microRNA-377 promoted tube formation in human umbilical vein endothelial cells. We then engineered rat MSCs with lentiviral vectors to either overexpress microRNA-377 (MSC^miR-377) or knockdown microRNA-377 (MSC^Anti-377) to investigate whether microRNA-377 regulated MSC-induced myocardial angiogenesis, using MSCs infected with lentiviral empty vector to serve as controls (MSC^Null). Four weeks after implantation of the microRNA-engineered MSCs into the infarcted rat hearts, the vessel density was significantly increased in MSC^Anti-377-hearts, and this was accompanied by reduced fibrosis and improved myocardial function as compared to controls. Adverse effects were observed in MSC^miR-377-treated hearts, including reduced vessel density, impaired myocardial function, and increased fibrosis in comparison with MSC^Null-group. These findings indicate that hypoxia-responsive microRNA-377 directly targets VEGF in MSCs, and knockdown of endogenous microRNA-377 promotes MSC-induced angiogenesis in the infarcted myocardium. Thus, microRNA-377 may serve as a novel therapeutic target for stem cell-based treatment of ischemic heart disease.     

miRTar Hunter: A prediction system for identifying human microRNA target sites

MicroRNAs (miRNAs) are important regulators of gene expression and play crucial roles in many biological processes including apoptosis, differentiation, development, and tumorigenesis. Recent estimates suggest that more than 50% of human protein coding genes may be regulated by miRNAs and that each miRNA may bind to 300–400 target genes. Approximately 1,000 human miRNAs have been identified so far with each having up to hundreds of unique target mRNAs. However, the targets for a majority of these miRNAs have not been identified due to the lack of large-scale experimental detection techniques. Experimental detection of miRNA target sites is a costly and time-consuming process, even though identification of miRNA targets is critical to unraveling their functions in various biological processes. To identify miRNA targets, we developed miRTar Hunter, a novel computational approach for predicting target sites regardless of the presence or absence of a seed match or evolutionary sequence conservation. Our approach is based on a dynamic programming algorithm that incorporates more sequence-specific features and reflects the properties of various types of target sites that determine diverse aspects of complementarities between miRNAs and their targets. We evaluated the performance of our algorithm on 532 known human miRNA:target pairs and 59 experimentally-verified negative miRNA:target pairs, and also compared our method with three popular programs for 481 miRNA:target pairs. miRTar Hunter outperformed three popular existing algorithms in terms of recall and precision, indicating that our unique scheme to quantify the determinants of complementary sites is effective at detecting miRNA targets. miRTar Hunter is now available at http://203.230.194.162/~kbkim.