Question: Normalizing counts using TIN metric before DESeq2
0
gravatar for hsbio
4 weeks ago by
hsbio10
hsbio10 wrote:

Hi,

I have some total RNA-seq data that was obtained from RNA that had quite low RIN values (unfortunately due to the source material, this was pretty much unavoidable). I've been looking at different ways people have tried to approach this problem and try and correct for it. I've seen that the transcript integrity number (TIN) has been developed in part to address this problem and one of the ways they did this was to first calculate a median TIN value for each sample and then correct the gene read counts based on this. In the paper they regressed the gene level read counts to the TIN value using a loess method.

And I was basically wondering how would be the most sensible way/if it is even possible to implement this sort of normalization in a Salmon -> tximport > DESeq2 pipeline.

I guess my main question would be, on what counts I would perform this normalization.

Would it be best to just do it on the the counts after importing to gene level with tximport? and then feeding these normalised counts back to the tximport matrix before running DESeq2 (with DESeqDataSetFromTximport)?

Or would running tximport using countsFromAbundance and then using these counts for normalization and then feeding the normalised counts straight to DESeq2 (with DESeqDataSetFromMatrix) be better?

I am just unsure of how the offset matrix and things would play into all this. Are any of these approaches sensible here? Apologies if this is all way off track/not clear.

Any other suggestions on how to best deal with data with low RINs would also be appreciated.

Thanks.

deseq2 salmon tximport • 60 views
ADD COMMENTlink modified 4 weeks ago by Michael Love25k • written 4 weeks ago by hsbio10
Answer: Normalizing counts using TIN metric before DESeq2
1
gravatar for Michael Love
4 weeks ago by
Michael Love25k
United States
Michael Love25k wrote:

A week ago I posted this about RIN and TIN correction. Can you take a look and reply here if you have more questions:

https://support.bioconductor.org/p/124506/#124526

ADD COMMENTlink written 4 weeks ago by Michael Love25k

Hi Michael,

Many thanks for the response.

I have never used RUV or SVA before but I will have a read of those and try them here and will post if I have any further questions.

Thanks again.

ADD REPLYlink written 21 days ago by hsbio10
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 196 users visited in the last hour