behaviour of bowtie compared to junction aligners?
1
0
Entering edit mode
6.8 years ago
popayekid55 ▴ 110

Hi,

I have human RNASeq data generated using lexogen kit. I have followed their analysis protocol (to process raw data). Then used tophat2 for alignment, observed only 53%. check unalinged using blast they are going human rRNA sequences end2end. Then i used hisat2, the % increased to 75. Then tried bowtie2 end2end parameter (which does not align junction read i assume) has given 79%.

Unaligned reads from different aligners checked against NCBI NT using BLAST they are hitting to human rRNAs/ncRNAs. Questions,

  1. How bowtie is preforming better than junction aligners?

  2. How do i improve the overall alignment?

Reference used was GRCh38 from ensemble

I have seen the same trend for fungi sample as well.

Any in site what could be the reason?

RNA-seq alignment • 1.7k views
ADD COMMENT
0
Entering edit mode

Did you use equivalent settings for each aligner? Could you post the command-line for each aligner?

ADD REPLY
0
Entering edit mode

used default parameters, bowtie2 end-end, same in tophat enabled b2-sensitive.

ADD REPLY
1
Entering edit mode
6.8 years ago
  1. It shouldn't, since it'll miss everything that's spliced
  2. Use STAR or another mapper that allows local alignment. You should get 90+% alignment.
ADD COMMENT
0
Entering edit mode

I think you missed this part.

Unaligned reads from different aligners checked against NCBI NT using BLAST they are hitting to human rRNAs/ncRNAs.

popayekid55 : Not much you can do if above is true.

ADD REPLY
0
Entering edit mode

I saw that, my presumption is that the edit distance is high enough that the reads aren't aligned by the aligners used.

ADD REPLY
0
Entering edit mode

Or the reference did not contain rDNA repeat.

ADD REPLY
0
Entering edit mode

That's usually part of the human reference, though if this were mouse I'd agree.

ADD REPLY
0
Entering edit mode

Unaligned reads from different aligners checked against NCBI NT using BLAST they are hitting to human rRNAs/ncRNAs.

Reason I thought that the reference here must not have rDNA.

ADD REPLY
0
Entering edit mode

Not much you can do if above is true.

how bowtie is performing better than splice aligners?

ADD REPLY
0
Entering edit mode

TopHat2 generally uses bowtie2. Are you sure you were using bowtie v.1?

ADD REPLY
0
Entering edit mode

i have used bowtie2-2.2.5 version. Same version was imported to tophat as well

ADD REPLY

Login before adding your answer.

Traffic: 1570 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6