<p>STAR quantMode (GeneCounts) essentially provides the same output as HTSeq-Count would, ie. number of reads that cover a given gene. This is the most simple measure of expression you could get from RNA-seq data. Kallisto and Salmon utilize pseudo-alignment to determine expression measures of transcripts (as opposed to genes). RSEM uses some algorithm to determine isoform fractions, but still based on alignments. (Note that I don't have any personal experience using RSEM.)</p>

<p>It ultimately depends on the research question which method will be most appropriate for you.</p>

<h1>read counts</h1>

<p>pros:</p>

<ul>
<li>relatively easy to understand concept-wise</li>
<li>can be used to generate many different expression measures (FPKMs, TPM, CPM, etc.) and is, thus, more versatile in down-stream analyses</li>
</ul>

<p>cons:</p>

<ul>
<li>you can't distinguish isoforms (from the counts themselves)</li>
<li>requires alignment and, thus, more time to obtain</li>
</ul>

<h1>salmon/kallisto:</h1>

<p>pros:</p>

<ul>
<li>faster (due to lack of alignment step)</li>
<li>quantifies (known) isoforms</li>
</ul>

<p>cons:</p>

<ul>
<li>expression values are more abstract, less easy to explain</li>
<li>gene-level analyses will require some extra work</li>
</ul>

<h1>RSEM</h1>

<p>Likely lies somewhere in between, still requiring alignments will make it slower. But it does give you expression values for isoforms. </p>

CC BY-SA 4.0
2020-02-21T09:01:16.057Z
2020-02-21T09:01:16.057Z
DavyCats
<p>A shorter version...STAR's is convenient, in that it gets done right alongside the alignment, but RSEM and Kallisto are smarter about dealing with ambiguous reads.  So I'd say your results from either of those are superior.</p>

<p>Personally, I use STAR for alignment, then give the bam to RSEM.  STAR's gene count is good for a reality check to make sure that I got the strandedness right, and RSEM's gene level expected counts are good for DESeq once they are rounded.</p>

CC BY-SA 4.0
2020-02-22T01:52:11.633Z
You can also feed the transcriptome coordinate BAM file generated by STAR directly to salmon to quantify abundances if you want to rely on the STAR alignments.  As an added bonus, if you process the STAR alignments with salmon, you need not disallow indels (gaps) or soft clips in the alignments (--quantTranscriptomeBan Singleend instead of IndelSoftclipSingleend).  This can be useful if your transcriptome has variation with respect to the reference.
nomad
CC BY-SA 4.0
2020-02-21T19:53:43.243Z
2020-02-21T19:53:43.243Z
swbarnes2
<p>I recently discovered <a href="https://github.com/snakemake-workflows/rna-seq-star-deseq2" rel="nofollow noreferrer">this Snakemake pipeline</a> for RNASeq that uses STAR's quantMode to quantify gene expression for DESeq2 differential expression analysis.</p>

<p>In the past I've always seen the workflow as:</p>

<ol>
<li>STAR alignment</li>
<li>RSEM/HTSeqCount quantification</li>
<li>DESeq2 analysis</li>
</ol>

<p>More recently I've also seen:</p>

<ol>
<li>Kallisto/Salmon pseudo-alignment quantification</li>
<li>DESeq2/Sleuth analysis</li>
</ol>

<p>I tried searching online for differences between STAR quantMode vs other quantification algorithms but couldn't find many details. What are the benefits and drawbacks of using STAR quantMode vs RSEM/Kallisto/Salmon?</p>

<rna-seq><differential-expression><deseq2><rsem><kallisto>
STAR quantMode vs RSEM vs Kallisto
# STAR quantMode vs RSEM vs Kallisto

(Asked by: username_0 on 2020-02-20)

I recently discovered [this Snakemake pipeline](https://github.com/snakemake-workflows/rna-seq-star-deseq2) for RNASeq that uses STAR's quantMode to quantify gene expression for DESeq2 differential expression analysis.

In the past I've always seen the workflow as:

1. STAR alignment
2. RSEM/HTSeqCount quantification
3. DESeq2 analysis

More recently I've also seen:

1. Kallisto/Salmon pseudo-alignment quantification
2. DESeq2/Sleuth analysis

I tried searching online for differences between STAR quantMode vs other quantification algorithms but couldn't find many details. What are the benefits and drawbacks of using STAR quantMode vs RSEM/Kallisto/Salmon?

**Tags:** rna-seq, differential-expression, deseq2, rsem, kallisto

---

## 2 answers

### Answer by username_1 on 2020-02-21. Score: 3 | Accepted answer

STAR quantMode (GeneCounts) essentially provides the same output as HTSeq-Count would, ie. number of reads that cover a given gene. This is the most simple measure of expression you could get from RNA-seq data. Kallisto and Salmon utilize pseudo-alignment to determine expression measures of transcripts (as opposed to genes). RSEM uses some algorithm to determine isoform fractions, but still based on alignments. (Note that I don't have any personal experience using RSEM.)

It ultimately depends on the research question which method will be most appropriate for you.

# read counts

pros:

* relatively easy to understand concept-wise
* can be used to generate many different expression measures (FPKMs, TPM, CPM, etc.) and is, thus, more versatile in down-stream analyses

cons:

* you can't distinguish isoforms (from the counts themselves)
* requires alignment and, thus, more time to obtain

# salmon/kallisto:

pros:

* faster (due to lack of alignment step)
* quantifies (known) isoforms

cons:

* expression values are more abstract, less easy to explain
* gene-level analyses will require some extra work

# RSEM

Likely lies somewhere in between, still requiring alignments will make it slower. But it does give you expression values for isoforms.


---

### Answer by username_2 on 2020-02-21. Score: 0

A shorter version...STAR's is convenient, in that it gets done right alongside the alignment, but RSEM and Kallisto are smarter about dealing with ambiguous reads. So I'd say your results from either of those are superior.

Personally, I use STAR for alignment, then give the bam to RSEM. STAR's gene count is good for a reality check to make sure that I got the strandedness right, and RSEM's gene level expected counts are good for DESeq once they are rounded.

> **Comments:**
> **username_3**: You can also feed the transcriptome coordinate BAM file generated by STAR directly to salmon to quantify abundances if you want to rely on the STAR alignments.  As an added bonus, if you process the STAR alignments with salmon, you need not disallow indels (gaps) or soft clips in the alignments (--quantTranscriptomeBan Singleend instead of IndelSoftclipSingleend).  This can be useful if your transcriptome has variation with respect to the reference.  

---

