Go Back   SEQanswers > Bioinformatics > Bioinformatics
Similar Threads
Thread Thread Starter Forum Replies Last Post
DESeq2 question regarding contrasts disulfidebond Bioinformatics 0 10-15-2015 05:18 PM
DESeq2 Contrasts cacti Bioinformatics 4 08-04-2015 11:20 AM
DESeq2 model design + contrasts gamb Bioinformatics 39 06-25-2015 02:06 PM
Contrasts with DESeq2 mistrm RNA Sequencing 3 07-24-2014 09:03 AM
Differential expression analysis on RPKMs - Contrasts and contrasts of contrasts giorgifm Bioinformatics 6 08-16-2013 10:50 AM

Thread Tools
Old 05-21-2021, 10:43 AM   #1
Junior Member
Location: us

Join Date: May 2021
Posts: 1
Default Need Help Doing Complex Contrasts in DESeq2

I recently completed a 22 day time course with a resistant and susceptible plant with five timepoints (0hr - control), 1dpi, 7dpi, 14dpi, and 22dpi. I was able to do three replicates and sequenced my smaples t get RNAseq libraries (30 libraries total). I now have begun my RNA seq analysis in DESeq2. I have set up my samples object for the comparisons like so:


         timepoint treatment genotype
S1A_rep1         0        SA        S
S1A_rep2         0        SA        S
S1A_rep3         0        SA        S
S1B_rep1         1        SB        S
S1B_rep2         1        SB        S
S1B_rep3         1        SB        S
S1C_rep1         7        SC        S
S1C_rep2         7        SC        S
S1C_rep3         7        SC        S
S1D_rep1        14        SD        S
S1D_rep2        14        SD        S
S1D_rep3        14        SD        S
S1E_rep1        22        SE        S
S1E_rep2        22        SE        S
S1E_rep3        22        SE        S
R1A_rep1         0        RA        R
R1A_rep2         0        RA        R
R1A_rep3         0        RA        R
R1B_rep1         1        RB        R
R1B_rep2         1        RB        R
R1B_rep3         1        RB        R
R1C_rep1         7        RC        R
R1C_rep2         7        RC        R
R1C_rep3         7        RC        R
R1D_rep1        14        RD        R
R1D_rep2        14        RD        R
R1D_rep3        14        RD        R
R1E_rep1        22        RE        R
R1E_rep2        22        RE        R
R1E_rep3        22        RE        R
I am using my treatment column to design my comparisons. I have done the following comparisons so far:


My DESEq2 comparison is done withe the follwing code:

dds <- DESeqDataSetFromTximport(txi, sample, ~treatment)
And my results are retrieved like so:
results(deseqData, contrast=c("treatment","RB","RA"))
I now want to do the following comparisons:

(RB - SB) - (RA - SA) (compare RB and SB, then compare that to RA and SA, compared)
(RE + SE) - (RA + SA)
(RE + RD) - (RC + RB)
[(RE + RD) - (RC + RB)] - [(SE + SD) - (SC - SB)] (compare the early vs late response of the R plants to that of the susceptible)
I am not sure of what code to use. Would anybody be able to point me in the right direction? Any help appreciated.
pthom010 is offline   Reply With Quote
Old 09-22-2021, 04:37 AM   #2
Junior Member
Location: Switzerland

Join Date: Jul 2021
Posts: 1

You have probably solved your issue by now, but if you still struggle with your transcriptomic data, you could try using our Omics Playground platform. We have an interactive GUI that allows you to generate pairwise comparisons with a high degree of flexibility when you input your data into the platform. You can also intersect the results of pairwise comparisons and represent them as Venn diagrams and scatterplots.

If you want to give it a shot, you can try our free access demo version:

Hope that addresses your need!
BigOmics is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 03:18 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2022, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO