Take the 2-minute tour ×
Biology Stack Exchange is a question and answer site for biology researchers, academics, and students. It's 100% free, no registration required.

So in short, I did a blastn with standalone blast using all known miRNAs against the EST sequences in a genome (word length 7 and e-value 0.01). Now to confirm one of the next steps would be obtain the pre-miRNA sequence of each hit and do an mFOLD for structural confirmation. Each paper I read always mentions to get a sliding window (of 70 - 100 nt length) but how do you get this sliding window from the BLASTN output.. Do you have to code something yourself or is there an easier way to do it. Thanks.

share|improve this question
    
miRNAs in EST- I am a little skeptical on this.. –  WYSIWYG Apr 30 at 12:42
    
sorry for the late answer.. I am not sure if you still need it now. –  WYSIWYG Apr 30 at 12:44

2 Answers 2

up vote 1 down vote accepted

DISCLAIMER: I have never worked with miRNAs.

However, from what little I know, the mature miRNA is produced by modifying the pre-miRNA. This means that you should be able to infer the pre-miRNA by mapping your miRNA to the genome using a tool like BLAT. Since you already have an EST which represents a transcribed sequence, you can map your EST back to the genome. You can simplify this by using a tool like exonerate that also models splicing.

share|improve this answer

A BLASTN would work for mapping the mature miRNA but you should correct some parameters.

E-value is fine but increase the word length. Keep it close to 17 because you are interested in perfect matches. With miRNA small mismatches can mean a different miRNA altogether. Also run for only plus/plus alignments. (You can see this post on biostar for details on how to do it)

What you can also do is to use bowtie (sometimes BLAST is very slow and though BLAST allows multiple threads there is no option for utilize multiple cores i.e true parallelization). Again use the option --norc to avoid plus-minus alignments.

Use bowtie-build on your EST sequences to make an index.

However, there is one more problem; it is unlikely that you'll get pre-miRNA sequences in your EST collection. This is because the stable form is mature miRNA and pre is generally too shortlived to be sequenced properly without appropriate selective measures.

Capturing small RNAs by sequencing is also not that straightforward. It needs size fractionation of RNA and subsequent library preparation from the small RNA fraction. So I guess you will not discover miRNAs from ESTs.

If you have small RNA sequences then you can align it to the genome and then select the $\pm 70$ nt to find the hairpin. In fact there is an algorithm called mirdeep which does that; you can use it and avoid writing your own codes. If you are good at perl then you can modify it to suit your specific needs if any.

share|improve this answer
    
I switched to working with supercontigs –  The Last Word May 2 at 8:51
    
never used bowtie.. I am trying to work on blastn with the number of mismatches set to 2 at the moment. What would you suggest is a good clustering program to remove duplicate sequences –  The Last Word May 27 at 4:26
    
fastx_collapser will collapse redundant sequences. BLAST can get very slow for these things. It is better to use a short read aligner like bowtie. –  WYSIWYG May 27 at 4:46
    
thx.. will have to try bowtie maybe next time. Too far with the current project to retrack. –  The Last Word May 27 at 4:55

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.