DESeq usage of transcript TPM after tximport
1
0
Entering edit mode
KIR_F • 0
@kir_f
Last seen 7 months ago

Hello,

I have seen tximport will import TPM, estimated counts and effective lengths from transcript quantification tools into a list object.

As far as I know, DESeqDataSetFromTximport() can directly read this list ('txi' in the vignette) to create a DESeq object.

I understand TPM can be used by tximport() to recalculate the original estimated counts using 'lengthScaledTPM' or 'scaledTPM'.

Since DESeq will work on transcript counts, I wanted to know if calculated TPM are going to be somehow used later by DESeq or otherwise they will be completely ignored (I would say no, but wanted to be sure).

tximport • 315 views
2
Entering edit mode
@mikelove
Last seen 51 minutes ago
United States

When DESeq2 is used with tximport via DESeqDataSetFromTximport, the TPM data (abundance matrix) is used to calculate a per-gene-and-sample offset. After that only the count and the offset are used. Or, if countsFromAbundance is used when calling tximport, then only the count is used.