SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
DESeq2 model design + contrasts gamb Bioinformatics 39 06-25-2015 02:06 PM
DESeq2 multifactorial design Erinyes Bioinformatics 5 09-19-2014 08:32 AM
Deseq2 design Wei-HD Bioinformatics 11 09-04-2014 11:13 AM
Paired design versus unpaired design in DESeq2 KristenC RNA Sequencing 1 05-29-2014 11:05 AM
non balanced design in DESeq2 libro.s Bioinformatics 1 01-10-2014 05:10 AM

Reply
 
Thread Tools
Old 11-07-2014, 07:48 AM   #1
sisterdot
Junior Member
 
Location: Europe

Join Date: Apr 2013
Posts: 6
Default DESeq2 multifactorial design

I hope someone can help me out in defining the right design formulas and contrasts to use with DESeq2 for a (subjectively) rather complex experimental design.

The two infection types infA and infB were considered to be equivalent by the experimental scientists.

The experimental design was as follows:
Code:
	genotype	sex	inf_type	inf
m_KO_iA	KO	m	infA	inf
m_KO_1	KO	m	NOinf	NOinf
m_KO_iB	KO	m	infB	inf
m_KO_2	KO	m	NOinf	NOinf
m_WT_iA	WT	m	infA	inf
m_WT_1	WT	m	NOinf	NOinf
m_WT_iB	WT	m	infB	inf
m_WT_2	WT	m	NOinf	NOinf
f_KO_A	KO	f	infA	inf
f_KO_1	KO	f	NOinf	NOinf
f_KO_iB	KO	f	infB	inf
f_KO_2	KO	f	NOinf	NOinf
f_WT_iA	WT	f	infA	inf
f_WT_1	WT	f	NOinf	NOinf
f_WT_iA	WT	f	infA	inf
f_WT_2	WT	f	NOinf	NOinf
And the questions i would like to pose with DESeq2 are:
- how do the KO cells differ in their response to any treatment compared to WT cells.
- how do the KO cells differ in their response to treatment A (infA) compared to WT cells.
- how do the KO cells differ in their response to treatment B (infB) compared to WT cells.

any suggestions on design formulas and contrasts to use are very much appreciated.
i am really sorry for not providing a suggestion myself, but i am afraid that might get embarrassing

THANKS
sisterdot is offline   Reply With Quote
Old 11-07-2014, 08:25 AM   #2
Wallysb01
Senior Member
 
Location: San Francisco, CA

Join Date: Feb 2011
Posts: 286
Default

You should read up on the interaction terms and changing the comparisons when you call the results function.

The design probably needs to look something like:

design ~ genotype + sex + inf_type + inf + sex:inf + sex:inf_type + genotype:inf + genotype:inf_type

By default the results function will test the last term in your design function, so you can change that or just use the comparisons option with results.
Wallysb01 is offline   Reply With Quote
Old 11-09-2014, 11:33 AM   #3
sisterdot
Junior Member
 
Location: Europe

Join Date: Apr 2013
Posts: 6
Default

Thanks for your reply...

i am going for
~ genotype + sex + inf + sex:inf + genotype:inf
for now...

gave a larger set of significant genes than not accounting for sex
~ genotype + inf + genotype:inf

is it possible to check in advance which main effects and interactions are best to include ?
sisterdot is offline   Reply With Quote
Reply

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 04:31 AM.


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