Question: Very low aligbnment rate for miRNA
0
gravatar for prp491
3.5 years ago by
prp49110
prp49110 wrote:

I know this kind of question have been asked before but I am not able to solve my problem. I have fastq file containg reads from small RNAseq of non-model plant. I cleaved my adapter (AGATCGGAAGAGCACA) and found peak at 24 nt. I search out the conserved miRNA from miRbase database and predicted the novel miRNA by mirPREFFER. But I got very low percentage (>1% )of alignment when I aligned the my fastq file the predicted miRNAs of my plant. I tried two aligner Bowtie2 and novoalign

I have uploaded the the raw, adapter-trimmed and length sorted fastq files and their respectice fastqc reports HERE

Any help will be much appreciated Thanks

mirna rna-seq • 1.4k views
ADD COMMENTlink modified 3.5 years ago by igor8.6k • written 3.5 years ago by prp49110

How long are your sequences? If it is less than 25, it would be the limit for most mappers. Since it is RNAseq could you try segemehl or STAR and see if there are more valid alignments?

ADD REPLYlink written 3.5 years ago by Rohit1.4k

I am trying to map 18-24 nt long reads only. Since in literature people have successfully used the Bowtie2 for miRNAseq that's why I used it. I will try to use the START as well as segemehl . Thanks

ADD REPLYlink written 3.5 years ago by prp49110
1
gravatar for harold.smith.tarheel
3.5 years ago by
United States
harold.smith.tarheel4.4k wrote:

Also, check your miRNA reference database to ensure that it contains only a single copy of each unique miRNA sequence. Most miRNAs are encoded by multiple loci and/or are part of highly homologous families, and will be flagged as multi-mappers unless you use the proper reference. You should also specify a shorter seed length for your aligner (per Rohit), and it may also be useful to increase the stringency of alignment (zero or one mismatch).

ADD COMMENTlink written 3.5 years ago by harold.smith.tarheel4.4k

Here I am confused little bit. Should I used the mature or precursor sequences as reference? My Bowtie2 command looks like this bowtie2 --local -L 15 -N 1

Is it OK? Thanks

ADD REPLYlink written 3.5 years ago by prp49110

+1 For suggesting to tweak the aligner's settings. With bwa mem alignments with score <30 are not returned by default (-T option), this means that even if a match is perfect but too short (say 22bp) you get nothing in output. Again in bwa mem, the minimum seed length is 19 (-k option) which might be too long for miRNA. (I assume bowtie has similar settings to play with).

ADD REPLYlink written 3.5 years ago by dariober10k
1
gravatar for igor
3.5 years ago by
igor8.6k
United States
igor8.6k wrote:

Most RNA-seq tools aren't really designed for very short sequences that you get from small RNA-seq. I would suggest using a miRNA-specific tool, such as:

ADD COMMENTlink written 3.5 years ago by igor8.6k

Thanks for the advice. miRge is using the Bowtie for alignment that is what I am using for my analysis.

ADD REPLYlink written 3.5 years ago by prp49110

Right, but they are also adjusting Bowtie settings.

ADD REPLYlink written 3.5 years ago by igor8.6k

I+ just tricked the miRge to run with Bowtie2 but I found the same low rate of alignment :(

ADD REPLYlink written 3.5 years ago by prp49110
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2623 users visited in the last hour