<p>I have a RNA seq data which I am trying to identify DEGs. Dataset contains:  </p>

<p>Untreated - 2 replicates (time point 1st day)<br>
Treated - 4 replicates (2 replicates at 3rd day &amp; 2 replicates at 7th day ) </p>

<p>I merged all the samples and tried to identify DEGs untreated vs treated. The 1st day sample is considered as control (untreated) for the study.</p>

<ol>
<li>Is it correct to merge all samples and compare. Or should I compare the two time point treated samples with untreated as 1st day vs 3rd day, 1st day vs 7th day and 3rd day vs 7th day separately.</li>
</ol>

<p>I tried LRT test in DEseq2 to understand the changes in gene expression across time between untreated and treated. The code is given below.</p>

<pre><code>sampleFiles&lt;-grep('SRR',list.files(directory),value=TRUE)
sampleCondition&lt;-c('untreated','untreated','treated','treated','treated','treated')
sampleTable&lt;-data.frame(sampleName=sampleFiles, fileName=sampleFiles, condition=sampleCondition)
ddsHTSeq&lt;-DESeqDataSetFromHTSeqCount(sampleTable=sampleTable, directory=directory, design=~condition)
colData(ddsHTSeq)$condition&lt;-factor(colData(ddsHTSeq)$condition, levels=c('untreated','treated'))
dds &lt;- DESeq(ddsHTSeq, test="LRT", reduced=~1)
resLRT &lt;- results(dds)
resOrdered &lt;- resLRT[order(resLRT$padj),]
</code></pre>

<p>After running the code I am getting huge number of DEGs (~15000 genes) as significant DEGs after applying FDR correction cutoff &lt;0.05.</p>

<ol start="2">
<li><p>Is there any error in mentioning the condition in the code.</p></li>
<li><p>Can you please tell the condition I have specified for LRT test in the code is enough to get the desired output.</p></li>
<li><p>Whether LRT test gives huge no of DEGs after analysis as significant DEGs.</p></li>
</ol>

CC BY-SA 4.0
2018-05-15T10:02:20.327Z
Redo this with all 3 groups indicated in `condition` rather than merging the two latter time-points together.
Devon Ryan
CC BY-SA 4.0
2018-05-15T11:39:24.337Z
So all the analysis for the above mentioned conditions(Untreated, Treated T1, Treated T2) must be performed separately? and How can I infer results based on DEG from 3 conditions?
AISWARYA JAYAPRAKASH
CC BY-SA 4.0
2018-05-15T12:42:57.327Z
All this analysis is really going to tell you is which genes were DE on day 3/7 vs. day 1. In order to determine which gene expression changes are caused by your treatment (as opposed to any number of other factors that might vary between days), you need to compare treated vs. mock treated samples on the same days.
heathobrien
CC BY-SA 4.0
2018-05-15T12:44:58.063Z
You can still perform the exact same LRT with 3 groups, it will be more informative than with 2.
Devon Ryan
CC BY-SA 4.0
2022-03-29T15:01:12.663Z
" LRT with 3 groups" this step dds <- DESeq(ddsHTSeq, test="LRT", reduced=~1)
 i have very fundamental doubt. Here it is single condition which is treated/untreatetd under sample condition. LRT is explained in the manual to see the difference between full and reduced model. How does it work here since it compares across the three condition. How do i infer which is being here used as reference here?
kcm
<r><rna-seq><deseq2>
Likelihood Ratio Test in DEseq2
# Likelihood Ratio Test in DEseq2

(Asked by: username_0 on 2018-05-15)

I have a RNA seq data which I am trying to identify DEGs. Dataset contains:  
Untreated - 2 replicates (time point 1st day)
Treated - 4 replicates (2 replicates at 3rd day & 2 replicates at 7th day ) 
I merged all the samples and tried to identify DEGs untreated vs treated. The 1st day sample is considered as control (untreated) for the study.

Is it correct to merge all samples and compare. Or should I compare the two time point treated samples with untreated as 1st day vs 3rd day, 1st day vs 7th day and 3rd day vs 7th day separately.

I tried LRT test in DEseq2 to understand the changes in gene expression across time between untreated and treated. The code is given below.
sampleFiles<-grep('SRR',list.files(directory),value=TRUE)
sampleCondition<-c('untreated','untreated','treated','treated','treated','treated')
sampleTable<-data.frame(sampleName=sampleFiles, fileName=sampleFiles, condition=sampleCondition)
ddsHTSeq<-DESeqDataSetFromHTSeqCount(sampleTable=sampleTable, directory=directory, design=~condition)
colData(ddsHTSeq)$condition<-factor(colData(ddsHTSeq)$condition, levels=c('untreated','treated'))
dds <- DESeq(ddsHTSeq, test="LRT", reduced=~1)
resLRT <- results(dds)
resOrdered <- resLRT[order(resLRT$padj),]

After running the code I am getting huge number of DEGs (~15000 genes) as significant DEGs after applying FDR correction cutoff <0.05.

Is there any error in mentioning the condition in the code.
Can you please tell the condition I have specified for LRT test in the code is enough to get the desired output.
Whether LRT test gives huge no of DEGs after analysis as significant DEGs.

**Tags:** r, rna-seq, deseq2
> **Comments:**
> **username_1**: Redo this with all 3 groups indicated in `condition` rather than merging the two latter time-points together.  
> **username_0**: So all the analysis for the above mentioned conditions(Untreated, Treated T1, Treated T2) must be performed separately? and How can I infer results based on DEG from 3 conditions?  
> **username_2**: All this analysis is really going to tell you is which genes were DE on day 3/7 vs. day 1. In order to determine which gene expression changes are caused by your treatment (as opposed to any number of other factors that might vary between days), you need to compare treated vs. mock treated samples on the same days.  
> **username_1**: You can still perform the exact same LRT with 3 groups, it will be more informative than with 2.  
> **username_3**: " LRT with 3 groups" this step dds <- DESeq(ddsHTSeq, test="LRT", reduced=~1)  i have very fundamental doubt. Here it is single condition which is treated/untreatetd under sample condition. LRT is explained in the manual to see the difference between full and reduced model. How does it work here since it compares across the three condition. How do i infer which is being here used as reference here?  

---

## 0 answers

