Importing data into the DESeq2 pipeline using tximport
1
0
Entering edit mode
@nikolay-ivanov-23079
Last seen 2.7 years ago
USA/New York City/Weill Cornell Medicine

I have a question regarding importing data into the DESeq2 pipeline using tximport.

I’m reading in salmon quant files using txiport:

tx_salmon = tximport(salmonFiles, type = "salmon", tx2gene = tx2, ignoreTxVersion=T)

After that, I make a DESeqDataSet:

dds = DESeqDataSetFromTximport(tx_salmon, colData = metadata, design = ~ condition)

Since tx_salmon$counts gives you un-normalized counts, the code above is absolutely equivalent to:

dds = DESeqDataSetFromMatrix(countData = tx_salmon$counts, colData = metadata, design = ~ condition)

Is that correct? I just wanted to confirm.

Thank you!

DESeq2 tximport • 4.4k views
ADD COMMENT
1
Entering edit mode
@mikelove
Last seen 9 hours ago
United States

Not equivalent.

https://bioconductor.org/packages/release/bioc/vignettes/DESeq2/inst/doc/DESeq2.html#transcript-abundance-files-and-tximport-tximeta

In particular, this quote from DESeq2 vignette, see bold section:

"Full details on the motivation and methods for importing transcript level abundance and count estimates, summarizing to gene-level count matrices and producing an offset which corrects for potential changes in average transcript length across samples are described in (Soneson, Love, and Robinson 2015)."

Also read here, from tximport vignette:

https://bioconductor.org/packages/release/bioc/vignettes/tximport/inst/doc/tximport.html#Downstream_DGE_in_Bioconductor

In particular:

Note: Do not manually pass the original gene-level counts to downstream methods without an offset.

ADD COMMENT
0
Entering edit mode

Got it. Thank you so much!

ADD REPLY

Login before adding your answer.

Traffic: 849 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6