Incorporating spike-ins to rna-seq analysis
1
0
Entering edit mode
cicely.day • 0
@e356d75d
Last seen 4 weeks ago
France

Hi,

I have been analysing RNA-seq data using DESEQ2, but keep being asked the same question. Can you / why haven't you incorporated the spike-ins for normalisation.

From my perspective, this isn't necessary due to the tests DESEQ2 runs, which consider composition bias/sequencing depth, which spike-in normalisation doesn't. Is this right? If not, does anyone know a way spike-ins can be incorporated into DESEQ2 analysis?

For context, I have been using the DESeqDataSetFromHTSeqCount function as I have read count files produced from my alignment script.

Thanks in advance, C

Spike-in RNASeqR DESeq2 Normalization rna-seq • 110 views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 5 hours ago
United States

Yes, use controlGenes, run estimateSizeFactors() before running DESeq() and that gives you spike-in based scaling factors.

If there are large changes in the distribution of abundance across samples, then the methods that assume focal changes (the default library size correction) will not work. So this is why such experiments require additional information on what features can be considered unchanging.

ADD COMMENT

Login before adding your answer.

Traffic: 267 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