miRPlant: an integrated tool for identification of plant miRNA from RNA sequencing data (2024)

  • Journal List
  • BMC Bioinformatics
  • v.15(1); 2014
  • PMC4141084

As a library, NLM provides access to scientific literature. Inclusion in an NLM database does not imply endorsem*nt of, or agreement with, the contents by NLM or the National Institutes of Health.
Learn more: PMC Disclaimer | PMC Copyright Notice

miRPlant: an integrated tool for identification of plant miRNA from RNA sequencing data (1)

Link to Publisher's site

BMC Bioinformatics. 2014; 15(1): 275.

Published online 2014 Aug 12. doi:10.1186/1471-2105-15-275

PMCID: PMC4141084

PMID: 25117656

Jiyuan An,miRPlant: an integrated tool for identification of plant miRNA from RNA sequencing data (2) John Lai, Atul Sajjanhar, Melanie L Lehman, and Colleen C NelsonmiRPlant: an integrated tool for identification of plant miRNA from RNA sequencing data (3)

Author information Article notes Copyright and License information PMC Disclaimer

Associated Data

Supplementary Materials



Small RNA sequencing is commonly used to identify novel miRNAs and to determine their expression levels in plants. There are several miRNA identification tools for animals such as miRDeep, miRDeep2 and miRDeep*. miRDeep-P was developed to identify plant miRNA using miRDeep’s probabilistic model of miRNA biogenesis, but it depends on several third party tools and lacks a user-friendly interface. The objective of our miRPlant program is to predict novel plant miRNA, while providing a user-friendly interface with improved accuracy of prediction.


We have developed a user-friendly plant miRNA prediction tool called miRPlant. We show using 16 plant miRNA datasets from four different plant species that miRPlant has at least a 10% improvement in accuracy compared to miRDeep-P, which is the most popular plant miRNA prediction tool. Furthermore, miRPlant uses a Graphical User Interface for data input and output, and identified miRNA are shown with all RNAseq reads in a hairpin diagram.


We have developed miRPlant which extends miRDeep* to various plant species by adopting suitable strategies to identify hairpin excision regions and hairpin structure filtering for plants. miRPlant does not require any third party tools such as mapping or RNA secondary structure prediction tools. miRPlant is also the first plant miRNA prediction tool that dynamically plots miRNA hairpin structure with small reads for identified novel miRNAs. This feature will enable biologists to visualize novel pre-miRNA structure and the location of small RNA reads relative to the hairpin. Moreover, miRPlant can be easily used by biologists with limited bioinformatics skills.

miRPlant and its manual are freely available at http://www.australianprostatecentre.org/research/software/mirplant or http://sourceforge.net/projects/mirplant/.

Keywords: RNA-seq, miRNA, Plant small RNA, RNA secondary structure


miRNA is a class of non-coding endogenous small RNA that post transcriptionally regulates target genes [1]. miRDeep-P [2] is one of the most commonly used computational plant miRNA identification tool, which is based on the miRDeep [3] algorithm.

The most challenging problem in identifying novel plant miRNA is to find a suitable genomic region as a miRNA precursor candidate (to test whether it forms hairpins) because the majority of precursor miRNA in plants are between 100-200bp [4], which is much longer than those in animals. Approaches using a shorter miRNA precursor may result in false negatives if the miRNA is longer and more variable than the predicted precursor region. Conversely, using a longer candidate precursor region to test whether it forms a hairpin structure may result in a non-complimentary match for the mature miRNA within the candidate precursor miRNA. Thus, in miRPlant, after small RNA sequencing reads are mapped to the genome, genomic regions around mapped reads are extended by 200bp to determine whether they form hairpin structures. To ensure detection of short plant miRNA, we also scan 100bp regions to see if we can detect a hairpin. This strategy can detect bona fide miRNAs that would otherwise be missed if only the longer (200bp) precursor candidate length was used.

The strategy for determining the precursor region is different between miRDeep-P and miRPlant. miRDeep-P determines the precursor region based on the genomic region having overlapping reads, while miRPlant determines a precursor region based on the mature miRNA region (or highest expressed read). The latter strategy can reduce the number of false negative results [5, 6], as it guarantees that the mature miRNA is located at the end of one arm of the stem loop.It is important that biologists with basic computer skills can easily use RNAseq tools in order to broaden research within this field. Thus, miRPlant was developed using the platform independent computer language Java. A Graphical User Interface (GUI) is employed whereby a complete pipeline analysis of raw data input is achieved in a few clicks of buttons: (.fastq files) - > mapping (.bam files) - > miRNA identification, expression, and secondary structure display - > mRNA target prediction. To further streamline accessibility of miRPlant, the tool does not require any third party tool. miRPlant also has a detailed but concise data output display that can be exported for publication in different file formats such as eps, pdf and svg (Figure1). miRPlant images are generated dynamically.

Open in a separate window

Figure 1

Output display of predicted miRNA. The read location and number of reads are shown relative to the precursor hairpin structure. The red sequence represents the mature miRNA.


miRPlant operations can be divided into the following stages:

  • i.

    filter out reads if their length is out of the 10-23bp range, or which have a read-quality below the criteria that is set by user.

  • ii.

    aggregate exact reads into one.

  • iii.

    map aggregated reads to the genome reference without mismatch. miRPlant uses the Java-coded bowtie [7] alignment algorithm. BAM format is used to store mapped reads. Please note that the attribute “XS” in the BAM file is used to record the copy number of the read as introduced by miRDeep*.

  • iv.

    gather sequences in the reference genome flanking the RNAseq read (precursor miRNA region) to determine whether the genomic region forms a hairpin structure using the RNA secondary structure algorithm [8].

  • v.

    use the miRDeep model to calculate the score for each predicted miRNA to measure the strength of the prediction. A higher score equates to a higher probability that the predicted miRNA is true.

The miRPlant interface enables users to customize parameters since different plant species may have different miRNA biogenesis [2] (Figure2). The default precursor miRNA length is set to 200bp. Here the precursor length represents the length between the mature miRNA and the mature star miRNA; the two flanking sequences are excluded. miRPlant generates six output files similar to miRDeep*. Since the precursor length of plant miRNA is much longer than that of animals, the distance between the mature miRNA and mature star miRNA may be very long, which may result in the formation of an internal loop. Therefore, miRPlant allows for internal loops. The default minimum loop (including the distance from loop ends to the mature or star mature miRNA) size is 25bp. In predicting mature miRNA, miRPlant requires less than 10% (max inconRead Ratio option in GUI) of reads falling out of the predicted miRNA and star mature miRNA sequence. In miRDeep, RNAseq reads in the loop are counted as being consistent, but plant miRNA have very long loops. Thus, we exclude reads located within the loop region. The other parameters are the same as with miRDeep*.

Open in a separate window

Figure 2

Parameter settings for miRPlant. Adapter sequences need to be replaced as appropriate. Data processing by miRPlant depends on the extension of the input file. Mapping and identification is performed if the input file extension is “.fastq” or “.fa”. Only identification is performed if the file extension is “.bam”. Output “.result” files are shown after clicking “submit”.

Results and discussion

miRPlant has been tested on two rice datasets [9]. Both miRPlant and miRDeep-P employ the miRDeep score calculation, with miRPlant having better performance than miRDeep-P (Table1), largely because miRPlant uses a flexible method to form the precursor candidates from the genomic region surrounding RNAseq reads. We set a minimum score of four when using miRPlant. A detailed summary of results can be found in Additional file 1 using GEO access number GSM278571 and GSM278572 for the RNAseq datasets.

Table 1

Comparison table

Rice (GSM278571)Rice (GSM278572)
Precision0.82(31/38)0.95(36/38)0.7 (44/63)0.83 (52/63)
Recall0.22 (31/144)0.25 (36/144)0.24 (44/181)0.29 (52/181)

Open in a separate window

Precision = known MiR/predicted MiR Recall = known MiR/total known MiR.

To further confirm the advantaged of miRPlant, we have extended this analysis to three more species (Arabidopsis thaliana, Medicago truncatula and Prunus persica) comprising 16 small RNA sequencing datasets (Detailed information in Additional file 2). To compare the two tools, we rank the predicted miRNAs in descending order of score for each tool, and then take the top 100 miRNAs from miRPlant and miRDeep-P for our comparison. We show that miRPlant consistently outperforms these other tools in all samples (Table2, Additional files 3 and 4).

Table 2

Comparison table (ATH, MTR, PPE)

A. thaliana (Number of known miRNA: 121)M. truncatula (Number of known miRNA: 196)P. persica (Number of known miRNAs: 75)

Open in a separate window

Precision = known MiR/predicted MiR Recall = known MiR/total known MiR.


miRPlant is modelled off miRDeep* [5] for use with plant small RNA sequencing data. We have integrated all third party tools such as genomic mapping and RNA secondary structure prediction [8] into a Java library, which is seamlessly integrated into miRPlant.

Availability and requirements

Project name: miRPlant.

Project home page:http://www.australianprostatecentre.org/research/software/mirplant.

Operating system (s): Windows, Linux, Mac OS.

Programming language: Java.

Other requirements: JRE.

License: GNU General Public License.

Any restrictions to use by non-academics: None.

Electronic supplementary material

Additional file 1:List of all identified miRNAs from two rice small RNAseq data.(XLSX 337 KB)(337K, xlsx)

Additional file 2:Small RNA sequencing data details.(DOCX 44 KB)(44K, docx)

Additional file 3:Detailed result of miRPlant prediction.(XLSX 13 MB)(13M, xlsx)

Additional file 4:Detailed result of miRDeep-P prediction.(XLSX 541 KB)(541K, xlsx)


Colleen C Nelson is supported by the Queensland Smart Futures Premier’s Fellowship. This project was funded by the Prostate Cancer Foundation of Australia (PCFA PG25).


Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

This work was supported by The Commonwealth Government of Australia, Department of Health and Queensland State Government, Smart Futures Premier Fellowship - Colleen C Nelson.

Contributor Information

Jiyuan An, Email: ua.ude.tuq@na.j.

John Lai, Email: ua.ude.tuq@ial.2j.

Atul Sajjanhar, Email: ua.ude.nikaed@rahnajjas.luta.

Melanie L Lehman, Email: ua.ude.tuq@namhel.einalem.

Colleen C Nelson, Email: ua.ude.tuq@noslen.neelloc.


1. Pritchard CC, Cheng HH, Tewari M. MicroRNA profiling: approaches and considerations. Nat Rev Genet. 2012;13(5):358–369. doi:10.1038/nrg3198. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

2. Yang X, Li L. miRDeep-P: a computational tool for analyzing the microRNA transcriptome in plants. Bioinformatics. 2011;27(18):2614–2615. [PubMed] [Google Scholar]

3. Friedlander MR, Chen W, Adamidi C, Maaskola J, Einspanier R, Knespel S, Rajewsky N. Discovering microRNAs from deep sequencing data using miRDeep. Nat Biotechnol. 2008;26(4):407–415. doi:10.1038/nbt1394. [PubMed] [CrossRef] [Google Scholar]

4. Meyers BC, Axtell MJ, Bartel B, Bartel DP, Baulcombe D, Bowman JL, Cao X, Carrington JC, Chen X, Green PJ, Griffithsnes S, Jacobsen SE, Mallory AC, Martienssen RA, Poethig RS, Qi Y, Vaucheret H, Voinnet O, Watanabe Y, Weigel D, Zhu JK. Criteria for annotation of plant MicroRNAs. Plant cell. 2008;20(12):3186–3190. doi:10.1105/tpc.108.064311. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

5. An J, Lai J, Lehman ML, Nelson CC. miRDeep*: an integrated application tool for miRNA identification from RNA sequencing data. Nucleic Acids Res. 2013;41(2):727–737. doi:10.1093/nar/gks1187. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

6. Friedlander MR, Mackowiak SD, Li N, Chen W, Rajewsky N. miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 2012;40(1):37–52. doi:10.1093/nar/gkr688. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

7. Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10(3):R25. doi:10.1186/gb-2009-10-3-r25. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

8. Hofacker IL. Vienna RNA secondary structure server. Nucleic Acids Res. 2003;31(13):3429–3431. doi:10.1093/nar/gkg599. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

9. Zhu QH, Spriggs A, Matthew L, Fan L, Kennedy G, Gubler F, Helliwell C. A diverse set of microRNAs and microRNA-like small RNAs in developing rice grains. Genome Res. 2008;18(9):1456–1465. doi:10.1101/gr.075572.107. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

Articles from BMC Bioinformatics are provided here courtesy of BMC

miRPlant: an integrated tool for identification of plant miRNA from RNA sequencing data (2024)


Top Articles
Latest Posts
Article information

Author: Gregorio Kreiger

Last Updated:

Views: 6358

Rating: 4.7 / 5 (77 voted)

Reviews: 92% of readers found this page helpful

Author information

Name: Gregorio Kreiger

Birthday: 1994-12-18

Address: 89212 Tracey Ramp, Sunside, MT 08453-0951

Phone: +9014805370218

Job: Customer Designer

Hobby: Mountain biking, Orienteering, Hiking, Sewing, Backpacking, Mushroom hunting, Backpacking

Introduction: My name is Gregorio Kreiger, I am a tender, brainy, enthusiastic, combative, agreeable, gentle, gentle person who loves writing and wants to share my knowledge and understanding with you.