miRNA annotation tools comparison

In summary: I will show which is the best miRNA mapping tool. I used 4 options for this benchmarking:

I think that these are the most used, and other not used but good to try them. They were clearly developed for other purposes, but as well, they generate the input of many miRNA pipelines. I just wanted to know how well my tool was doing.

The first aim to develop miraligner was to get annotated additions of nucleotides at the end of miRNA sequences, something that is very common in mirna biogenesis: isomirs and often they are missed by short read and fast mappers.

I have a repository for this kind of things, so anybody can reproduce my results, and check if I did something wrong, or comment on it.

In this post I just want to know which tool detects more miRNA, for that I did two main steps:

  1. simulate a bunch of miRNAs (isomirs) with my python script that is part of SeqBuster suit. I generated around 10000 sequences. Normally, one small-RNAseq library produces around half million different sequences.
  2. use only miRBase human precursors as reference genome

I used default parameter for all, so probably there is a set of parameters that would be better for some tools, but I didn’t search for them this time (happy to accept issues to add them).

The results are:

  1. miraligner and STAR map more sequences. miraligner loses sequences shorter than 15 nt, normally (miRNA) are around 21, and those sequences map to repeat elements
  2. STAR is the best mapping, but need some parsing to reduce false positive. I think that pipelines should change to this tool.
  3. GEM has problems with additions and nt substitutions in many cases, same as novoaligner (I will look at this in the future)
  4. bowtie2/bowtie is the second that annotated most (and best)
  5. microrazer has a problem with mismatches, but maybe there is some parameter to trick
  6. miRExpress with default options only will map perfect matches sequences to precursor, so strongly recommended to allow errors to increase sensibility.

I have to say that the advantage of miraligner is that it gives you the sequence annotated as miRNA or precursor, and gives you the exactly modifications that sequences have if they are compared to miRBase database. And you could feed the results to my R package to plot isomers distribution of samples from different groups, and do differential expression analysis with DESeq2, or another tool.

I didn’t add time consumption because all of them run in a couple of minutes.

In my next post I will focus in the the correct annotation of each sequence, and the possible problems with cross-mapping events, when the sequence comes from another regions of the genome but map to miRBase precursor as well [see my previous post for more details].

As well I will use STAR with the full genome and see if the mapping continues being the best. In that case I will add a script to SeqBuster to parse the output of STAR for those who can map with STAR (need up to 32G for human genome)

Advertisements

About lpantano

Post-doctoral fellow at IBB PhD in Biomedicine at UPF Msc in Bioinformatics for health science Bachelor in Biochemistry at UG
This entry was posted in science and tagged , , , . Bookmark the permalink.

2 Responses to miRNA annotation tools comparison

  1. Ido Tamir says:

    Nice blog. You seem quite diligent in your analyses. Could you benchmark http://arn.ugr.es/srnabench/ also? It also detects isomirs.

    When will seqbuster be published?

    What do you use for differential expression of miRNAs? I looked briefly through your github repo but did not find anything.

    • lpantano says:

      Hi,

      thanks! I will add it!
      About DE analysis, actually the first version used a implemented z-score test, but nowadays I just use DESeq, quite robust for miRNA. I always do one filter before DE, just take miRNA expressed in 75% of the samples.
      SeqBuster is already published:

      SeqBuster is a bioinformatic tool for the processing and analysis of small RNAs datasets, reveals ubiquitous miRNA modifications in human embryonic cells. Pantano L, Estivill X, Martí E. Nucleic Acids Res. 2010 Mar;38(5):e34. Epub 2009 Dec 11.

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s