Abstract:BackgroundSeveral methods exist for the prediction of precursor miRNAs (pre-miRNAs) in genomic or sRNA-seq (small RNA sequences) data produced by NGS (Next Generation Sequencing). One key information used for this task is the characteristic hairpin structure adopted by pre-miRNAs, that in general are identified using RNA folders whose complexity is cubic in the size of the input. The vast majority of pre-miRNA predictors then rely on further information learned from previously validated miRNAs from the same or a closely related genome for the final prediction of new miRNAs. With this paper, we wished to address three main issues. The first was methodological and aimed at obtaining a more time-efficient predictor, however without losing in accuracy which represented a second issue. We indeed aimed at better predicting miRNAs at a genome scale, but also from sRNAseq data where in some cases, notably of plants, the current folding methods often infer the wrong structure. The third issue is related to the fact that it is important to rely as little as possible on previously recorded examples of miRNAs. We therefore also sought a method that is less dependent on previous miRNA records.ResultsAs concerns the first and second issues, we present a novel alternative to a classical folder based on a thermodynamic Nearest-Neighbour (NN) model for computing the free energy and predicting the classical hairpin structure of a pre-miRNA. We show that the free energies thus computed correlate well with those of RNAfold. This novel method, called Mirinho, has quadratic instead of cubic complexity and is much more efficient also in practice. When applied to sRNAseq data of plants, it gives in general better results than classical folders. On the third issue, we show that Mirinho, which uses as only knowledge the length of the loops and stem-arms and the free energy of the pre-miRNA hairpin, compares well with algorithms that require more information. The results, obtained with different datasets, are indeed similar to those of other approaches with which such a comparison was possible. These needed to be publicly available softwares that could be used on a large input. In some cases, Mirinho is even better in terms of sensitivity or precision.ConclusionWe provide a simpler and much faster method with very reasonable sensitivity and precision, which can be applied without special adaptation to the prediction of both animal and plant pre-miRNAs, using as input either genomic sequences or sRNA-seq data.

mirtronDB: a mirtron knowledge base

Pmirkb: A Plant Microrna Knowledge Base

Mirtrans: a Resource of Transcriptional Regulation on Micrornas for Human Cell Lines.

The Effect of Alternative Splicing Sites on Mirtron Formation and Arm Selection of Precursor microRNAs

TRmir: A Comprehensive Resource for Human Transcriptional Regulatory Information of MiRNAs

Mirfans: an Integrated Database for Arabidopsis Thaliana Microrna Function Annotations

Mirinho: An efficient and general plant and animal pre-miRNA predictor for genomic and deep sequencing data

Mir+pathway: the Integration and Visualization of Mirna and KEGG Pathways.

Mirmat: Mature Microrna Sequence Prediction

Bioinformatics Analysis on Structural Features of Microrna Precursors in Insects

Miras: a Data Processing System for Mirna Expression Profiling Study

miRBase: tools for microRNA genomics

miRBase: integrating microRNA annotation and deep-sequencing data

Mirsm: an R Package to Infer and Analyse Mirna Sponge Modules in Heterogeneous Data

MirGeneDB 3.0: Improved taxonomic sampling, uniform nomenclature of novel conserved microRNA families, and updated covariance models.

miRDB: a microRNA target prediction and functional annotation database with a wiki interface.

Simira: A Tool To Identify Coregulation Between Micrornas And Rna-Binding Proteins

mirPRo–a novel standalone program for differential expression and variation analysis of miRNAs

Mitalos V2: Analyzing Tissue Specific Microrna Function.

Genome-wide hairpins datasets of animals and plants for novel miRNA prediction

CG-content log-ratio distributions of Caenorhabditis elegans and Drosophila melanogaster mirtrons