miRNA-seq / Correlate miRNA with gene expression array data


Given a list of miRNAs and matching gene expression array data, this tool identifies which of the target genes correlate positively or negatively with the miRNA expression. Note that the miRNA and gene expression data sets must have the same number of samples.



You need to have miRNA-seq and gene expression array data from the same samples, and indicate the matching samples with numbers in a given phenodata column.

The analysis consists of two parts:

Firstly, predicted target genes for the list of miRNAs are fetched from the TargetScan and PicTar databases, and only those genes that are predicted in both databases are included. Then, another intersection is made with the list of genes present in the matching gene expression data set.

Secondly, the correlation is calculated between the expression of a miRNA and each of its target genes, using Pearson, Spearman or Kendall's measure of correlation. Pearson correlation is parametric, so it can be sensitive to outliers, while the rank-based Spearman method is more robust. In cases when the direction of the change between two-conditions, or time-points, is more important than the absolute values of the expression profile, Kendall's correlation method is recommended. The statistically significant results are reported in separate tables for the positively and negatively correlating genes.


Output consists of two text files, one listing the significantly positively correlating target genes and one listing the significantly negatively correlating ones.


This tool leverages some of the functionality from the RmiR Bioconductor package. More info on the databases of predicted miRNA targets can be found at: