Using salmon on HISAT2 BAM files for transcript quantification
1
1
Entering edit mode
12 weeks ago
skden2020 ▴ 10

Hi,

I am wondering if there is any need to align data at all if the goal is to quantify transcripts with salmon.

Are there any advantages to aligning versus doing direct. quantification using salmon, since both options seem to exist.

I have previously used HISAT2 mapping to ref genome and am wondering if I can perform Salmon quantification for transcripts on these BAM files.

Thanks, Swati

salmon alignment HISAT2 • 709 views
0
Entering edit mode

I am wondering if there is any need to align data at all if the goal is to quantify transcripts with salmon.

No, if the goal is quantification using salmon, there is no need for mapping.

Are there any advantages to aligning versus doing direct?

I think it would depend on your specie. Repeat content, G+C%, and multicopy gene families might affect the results, especially those produced by quasi-mappers (like salmon).

I have previously used HISAT2 mapping to ref genome and am wondering if I can perform Salmon quantification for transcripts on these BAM files.

You can't use bam files as input for salmon, can you? Besides, salmon has its own algorithm for "mapping", using a bam file produced from a different aligner would definitely affect the result. From salmon's description:

Salmon is a tool for quantifying the expression of transcripts using RNA-seq data. Salmon uses new algorithms (specifically, coupling the concept of quasi-mapping with a two-phase inference procedure) to provide accurate expression estimates very quickly (i.e. wicked-fast) and while using little memory.

0
Entering edit mode

Salmon can't use bam files, but it's easy enough to convert bam back to fastq files. I typically do this, but I don't have a fantastic reason why other than being stubborn.

0
Entering edit mode

If you don't have access to the fastq files yes, you can convert them, otherwise, I wouldn't say it is being stubborn.

¯\_(ツ)_/¯

0
Entering edit mode
6 weeks ago
ATpoint 70k

You can only use BAM files with salmon if these are in transcriptome space, so the aligner used the transcriptome rather than genome for mapping. See the manual for this option called 'alignment-based mode'. Still, I never got why people would do that rather than using salmon directly.

Traffic: 1439 users visited in the last hour
FAQ
API
Stats

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