Question: Mapping to transcriptome with bowtie2 vs mapping to genome using STAR
1
gravatar for geek_y
9 days ago by
geek_y9.3k
Barcelona/CRG/London/Imperial
geek_y9.3k wrote:

Hi,

Does any one know any paper or personally compared the results of the following approaches ?

  1. Mapping RNA-Seq reads to transcriptome fasta using bowtie2 and quantifying using RSEM or other tools.

  2. Mapping RNA-Seq data to genome fasta using STAR and quantifying using featureCounts.

In first case, there is less chances of finding multi-mapped reads because we focus only on transcribed regions and ignore rest of the genome. Don't know how often it's a problem or meaningful.

Some of the data I am looking at shows differences in quantifications coming from two methods for some important genes, so just thought if some one has more insights into it.

rna-seq star bowtie2 • 97 views
ADD COMMENTlink modified 9 days ago • written 9 days ago by geek_y9.3k

Have you looked at one or two of those genes in IGV? Are the reads to them multimappers (presumably due to a similar pseudogene)?

ADD REPLYlink written 9 days ago by Devon Ryan88k

I STAR re-aligned using very low stringent options, no luck so far. I am going to try exact method the authors described in their paper.

ADD REPLYlink written 9 days ago by geek_y9.3k
1
gravatar for geek_y
9 days ago by
geek_y9.3k
Barcelona/CRG/London/Imperial
geek_y9.3k wrote:

figured it out. Sometimes annotations includes 3' UTRs. So if the annotations we are using includes 3' UTRs, we get counts for genes. If not, we dont quantify them. Its not problem of mapping but the annotations we are using.

enter image description here

Here Refseq annotation did not include 3'UTR but the Gencode. So the quantifications differ based on this.

ADD COMMENTlink modified 9 days ago • written 9 days ago by geek_y9.3k
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: 874 users visited in the last hour