How well does "Trinity-Differential-Expression" work on normalized raw counts without RSEM abundance of estimation
Entering edit mode
harelarik • 50
Last seen 12 weeks ago


We have RNAseq data for ~120 samples. We are considering of using "Trinity Differential Expression" protocol described by Brian Haas in
How well does it work without the RSEM abundance of estimation?
We were thinking of feeding normalized gene counts (e.g., by TMM logCPM) into the Trinity protocol.

Thank you very much for your help,


Trinity DifferentialExpression RSEM RNAseq edgeR • 90 views
Entering edit mode

Not being the developer, but I would recommend to do these kind of analysis manually rather than in a closed end-to-end pipeline where you essentially have no control over any parameters or details. You are probably interested in differential gene rather than transcript-level expression? I'd start from raw gene-level counts, e.g. obtained via tximport, not normalized counts as these are not compatible with most common DE tools as explicitely mentioned in their vignettes. Be sure to properly perform QC to detect potential batch effects and outlier samples.

Entering edit mode

Thank you.

Regarding: "I'd start from raw gene-level counts, e.g. obtained via tximport" Did you mean without RSEM fix for abundance?

The protocol is recommending to use logcpmTMM normalization, but they do it on rsem counts not. I was wondering how important is that the input will be rsem. We have allready did the DE analysis with dseq2.


Login before adding your answer.

Similar Posts
Loading Similar Posts
Traffic: 216 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.4