This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Building RNAseq pipeline STAR-RSEM-tximport-DESeq2

Hello,

I did some googling about mRNAseq pipeline and built that pipeline. So I would like to know if my pipeline is okay.

Reference genome : Human gencode GRCh37
Aligner : STAR
Expression quantification : RSEM
Gene count table convert? : tximport
Statistical test : DESeq2

it seems like people use htseq or featureCount for gene count table and put it into edgeR or DESeq2. I am not sure, if I go for RSEM to run DESeq2 instead of htseq or featureCount. I tried kallisto and slueth, but I want to stick to this traditional alignment method.

Any advice would be apprieciated.
Thank you,

rna-seq pipeline
• 6,769 views • link
updated by Matteo Schiavinato Matteo Schiavinato • written by jkim • jkim

Unless you have a very specific reason to use GRCh37 you should consider using GRCh38. 38 has been available for a number of years and is now mature.

• 0 views • link
written by GenoMax GenoMax

Thank you for your comment. I have some old data based on GRCh37, So I go with GRCh37. Maybe next project I will go for GRCh38.

• 0 views • link
written by jkim • jkim

1 answer

I would like to know if my pipeline is okay.

RSEM can use STAR within its own script, so I would do that instead of running STAR yourself, unless you want to sophisticate your mapping more.

I am not sure, if I go for RSEM to run DESeq2 instead of htseq or featureCount.

The output is comparable, it's a matter of choice, perhaps RSEM is more sophisticated. But in the end it's counting mapped reads per gene/transcript.

I tried kallisto and slueth

Read this: http://www.rna-seqblog.com/limitation-of-alignment-free-tools-in-total-rna-seq-quantification/

• 0 views • link
written by Matteo Schiavinato Matteo Schiavinato

Log in to answer this question.

More posts like this