DESeq2 comparison two cell types and two genotypes
1
0
Entering edit mode
5 months ago

Dear All,

I have the samples of two cell types (Embryo and Ovule) with two genotypes (WT and MT), which indicate as condition and genotype below. Could you please suggest how to extract information from these four combinations?

I would like to compare:

  1. EmbryoWT vs EmbryoMT
  2. OvuleWT vs OvuleMT.
  3. EmbryoWT vs OvuleWT
  4. EmbryoMT vs OvuleMT

enter image description here

I used the code below with model design = ~ genotype + condition + genotype:condition to compare my samples and try to look for solutions. However, It seems to be the analysis will include the same condition or genotype together which is not our interest.

`dds <- DESeqDataSetFromMatrix(countData = clean_data, colData = meta, design = ~ genotype + condition + genotype:condition)

I highly appreciate your help.

Thanks in advance,

DESEQ2 • 534 views
ADD COMMENT
0
Entering edit mode
ADD REPLY
0
Entering edit mode
5 months ago

You can use interactions, but don't. Make a new column of tissue_genotype, make that your design.

ADD COMMENT
0
Entering edit mode

Hi swbarbes2

Thank you very much for your comments. Do you mean the interaction of genotype:condition from this model?

    design = ~ genotype + condition + genotype:condition

Could you please suggest how to extract this condition?

  1. EmbryoWT vs EmbryoMT
  2. OvuleWT vs OvuleMT.
  3. EmbryoWT vs OvuleWT
  4. EmbryoMT vs OvuleMT

    Best regards,

ADD REPLY
0
Entering edit mode

I think what swbarbes2 meant was to create third column concatenating the two other columns. You'll than have this column to make the design from

EmbryoMT
EmbryoMT
EmbryoMT
EmbryoWT
EmbryoWT
EmbryoWT
OvuleMT
OvuleMT
OvuleMT
OvuleWT
OvuleWT
OvuleWT

This you can use to make your analysis with.

To make your pair0wise comparison you can look at the resultNames() command or manually create your pairs, if what you need is not there by looking at the contrasts option

ADD REPLY
0
Entering edit mode

Hi Assa,

Thank you very much for your suggestion. In this case, which model would you suggest?

design = ~ genotype + condition + genotype:condition or only design = ~ concatenation

Best regards,

ADD REPLY
0
Entering edit mode

The second. With a complex call to results, you can get the first design to answer your questions, but using the simpler design is much harder to mess up, and much easier for you or someone else to follow down the road.

ADD REPLY
0
Entering edit mode

Dear swbarnes2,

Thank you very much for pointing this out. Based on this model, I can find EmbryoWT vs EmbryoMT and OvuleWT vs OvuleMT.

Would it be possible to find ALL WT (EmbryoWT+OvuleWT) vs ALL MT (EmbryoMT+OvuleMT)**?

Thank you very much for your help

ADD REPLY
1
Entering edit mode

Just use the original column in your design with these two factors WT and MT for that.

ADD REPLY

Login before adding your answer.

Traffic: 1368 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6