DEG with several conditions
Dear all,
I have RNA-Seq count data I want to analyze with DESeq2. My data has the following characteristics:
- there are 2 different cell lines
- there are 2 different treatments
- there are 3 different timepoints of measurement after treatment
- For both the cell lines I also have control samples
My main question is: should I create multiple DESeqDataSetFromMatrix objects, or can I just create one object and then perform all the analysis on that object? If yes, how would such a design formula look like?
Furthermore, how can I make generalized conclusions out of all these pairwise comparisons?
Thank you in advance for your inpunt!
• 1,081 views
•
link
updated
by
BioinfGuru
•
written
by
Kylie •
•
written
by
Kylie •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Tools to model DEG/pathways in (fast-dividing) CANCER cell lines across timepoints? (in bulk RNA-se…
written by txema.heredia •
-
DESeq2 for multiple groups
written by banaalkhayrat123 •
-
single cell data integration approach
written by Jay •
-
Can I infer the fraction of replicating cells from bulk RNA-seq data?
written by txema.heredia •
-
DESeq2 design comparing changes through 2 timepoints of 2 different treatments.
written by txema.heredia •
-
Multi factor analysis in DESeq2 and design formula
written by nik.kraemer •
-
Visualize overlap differentially expressed genes of multiple comparisons.
written by lucatucciarone
-
Differential Expression Help
written by mlai2567 •
-
Where to find DNase-seq and RNA-seq assays for various cancer lines?
written by bhchiang •
-
qRT-PCR analysis : average triplicates for statistical testing
written by Nicolas Rosewick
First of all, which are the comparisons you want to test? The question on whether to split or keep together usually comes down to whether you want to compare between cell lines, and whether one line has more spread than the other. Consider checking data using PCA, see vignette, to get an idea how the spread is. The question on "general conclusions" again depends on what you want to compare.
Ah good point. I want to do multiple things:
I already found with PCA that most of the variance (91%) is explained by difference in cell line.
Take a look at this blog, it was the biggest help to me: A guide to designs and contrasts in DESeq2
Section 4 will have the core of what you need but I'd strongly recommend also using the author's contraster() function so you can check what you are comparing (the comparison wording in the code is more intuitive).