Question: DESeq differential expression
2
gravatar for pmanga
6.4 years ago by
pmanga40
United States
pmanga40 wrote:

Hi,

Upon performing differential expression analysis for my RNA seq data set using DESeq I found that the log2 transformed differential expression values that I get from DESeq are very close to the RT-qPCR values but opposite in sign. I am new in the field and was concerned so as to why this would be happening? Is this something inherent to the way DESeq calculates the diff. values? 

Thanks.

rna-seq • 4.0k views
ADD COMMENTlink modified 6.4 years ago by Martombo2.7k • written 6.4 years ago by pmanga40
3

Have you tried to swap the "reference" transcriptome? Like you did AxB, if you do BxA, maybe?

ADD REPLYlink written 6.4 years ago by eddie.im130

Hi,

Ya, I tried that and it would change the signs but was confused whether it was the right thing to do by swapping the reference as I wasn't sure if DESeq was comparing b to a or a to b.The following reply cleared that too. Thanks :)

ADD REPLYlink written 6.4 years ago by pmanga40
4
gravatar for Martombo
6.4 years ago by
Martombo2.7k
Seville, ES
Martombo2.7k wrote:

when you perform a DESeq analysis, the program is comparing conditionB to conditionA, therefore the fold change will be a ratio of baseMeanB over baseMeanA. for example if you call DESeq like this:

nbinomTest( cds, "untreated", "treated" )

the fold change will be greater than 1 (and the log2 fold change positive), if the gene is more expressed in "treated". otherwise it will be lower than one (and the log2 fold change negative).

so you can simply change the sign to your log2 fold changes, or repeat the analysis changing the order of the conditions:

nbinomTest( cds, "treated", "untreated" )

ADD COMMENTlink modified 6.4 years ago • written 6.4 years ago by Martombo2.7k

Hi,

Thanks a lot, that was really helpful!!

ADD REPLYlink written 6.4 years ago by pmanga40

Also.. like we can write the final differential expression analysis file with the log2, p-values etc. and save it, is there a way to obtain the file with the normalized values before we perform the DE analysis (right after normalization)? I tried a similar write command with the normalization file name instead of the DE analysis file but that gives an error.

ADD REPLYlink written 6.4 years ago by pmanga40

 cdsBlind<- estimateDispersions(cds,method=”blind”)

 vsd<- varianceStabilizingTransformation(cdsBlind)

exprs(vsd)  will be the normalized counts table 

ADD REPLYlink written 6.4 years ago by Ming Tang2.6k

well that will give you the variance stabilized counts, which in general is not what you use in the DE.

to get the normalized counts table you can just use the function counts, setting the option normalized to TRUE:

counts( cds, normalized=TRUE )

(where cds is your CountDataSet object)

ADD REPLYlink modified 6.4 years ago • written 6.4 years ago by Martombo2.7k

It worked..Thanks!

ADD REPLYlink written 6.4 years ago by pmanga40

Hi,

This may be a weird question but, I have this normalized dataset that I want to analyze the same way as DESeq analyzes a data set normalized within it for differential expression.Is there a way to get DESeq to analyze it for differential expression (DE)? Or if not (as I read that it needs raw read data and normalized data set cant be used) how can I use similar test of variance i.e. if I am not getting it wrong DESeq uses a chi square test can I replicate the same conditions(?) it uses elswhere to simulate a similar DE analysis as by DESeq?

ADD REPLYlink written 6.4 years ago by pmanga40
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 2.3.0
Traffic: 1757 users visited in the last hour