Question: Differential expressed genes in DESeq
0
gravatar for Lila M
3.0 years ago by
Lila M 820
UK
Lila M 820 wrote:

Hi guys, I have an issue... I'm doing some analysis using DESEq2. By default I know that when I create the object dds using

dds <- DESeqDataSetFromMatrix(countData = count.table,
  colData = data.frame(treatment=treatment),
  design = ~ treatment)

an as treatment (control vs recurrence) by default R R will choose as reference level the factors based on alphabetical order, so in this case, the reference level would be control. But I have a question about this, for the downstream analysis, when I run for example:

dds <- DESeq(dds)
res <- results(dds)
res2Order <-res2[order(res2$pvalue), ]
write.table(res2Order, file="my _table")

I assume that the genes reported in my table with p-value < 0.05 are deferentially expressed in cotrol vs. recurred, but does it means that they are more or less expressed in control than recurred, but how can I know if this difference is based on a higher or lower expression? For example,can I confirm that a gene with p-value < 0.05 is more expressed in control than recurrence?

Many thanks in advance!

rna-seq deseq rna expression • 1.4k views
ADD COMMENTlink modified 3.0 years ago by Devon Ryan96k • written 3.0 years ago by Lila M 820
2
gravatar for Devon Ryan
3.0 years ago by
Devon Ryan96k
Freiburg, Germany
Devon Ryan96k wrote:

You'll note that there's a column with the fold change. Its sign dictates the direction of change.

ADD COMMENTlink modified 3.0 years ago • written 3.0 years ago by Devon Ryan96k

Thank you for the information, as is the first time that I have to deal with it, could you please help me to understand the results based on log2FoldChange?

ID                      log2FoldChange          pvalue
ENSG00000132002.6       -1.3096108594    6.38592315283661E-15
ENSG00000188517.13      2.568933441      1.05420723121764E-13
ENSG00000204388.6       1.4787338766     1.99686292186552E-12
ENSG00000240747.6       0.1054737261     1.53129979101136E-11

So in this specific case, does it means that the gene ENSG00000132002.6 is differential expressed in control vs recurred samples and indeed is less expressed in controls? Is the same for the log2fold change values ~ 0?

Many thanks!

ADD REPLYlink modified 3.0 years ago • written 3.0 years ago by Lila M 820

Yes , as Devon explain the sign of the log2FC tells you the direction of change. In the case of ENSG00000132002.6 the log2FC is -1.3096108594, thus foldchange = 2^-1.3096108594 = 0.4034297 so 2.478747 (1/0.4034297) times less expressed.

ADD REPLYlink written 3.0 years ago by Nicolas Rosewick9.0k

But is less expressed in recurred? I understand that if I've included as treatment: control and recurred, and by default R will choose as reference level the factors based on alphabetical order, control should be my reference level and in that case, recurred is less expressed than control. Is this assumption ok?

ADD REPLYlink written 3.0 years ago by Lila M 820
1

Yes. You can check the normalized read counts for ENSG00000132002 to be sure.

norm_count <- counts(dds,normalized=T)
ADD REPLYlink written 3.0 years ago by Nicolas Rosewick9.0k

Ok, many thanks for this tip!!!

ADD REPLYlink written 3.0 years ago by Lila M 820
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: 1291 users visited in the last hour