I am new to this RNAseq world
I want to perform differential gene expression analysis using DESeq2. I read the DESeq2 vignette but got confused in the Contrasts and Interaction section.
My metadata consists of 22 samples of adipose tissue, 12 from high diet-fed mice (HFD) and 10 from chow-fed mice (CD) each co-cultured with tumor (T) and control (C).
I have seen that adipose tissue from HFD already differs significantly from adipose tissue from CD (without tumor) but I want to test and filter out only those differences that resulted due to tumor co-culture in HFD compared to CD.
I used the following code to create the DESeq object
dds <- DESeqDataSetFromMatrix(countData = main_data, colData = meta_data, design = ~treatment + tumor_status + treatment:tumor_status)
dds <- DESeq(dds)
resultsNames(dds)
[1] "Intercept" "treatment_HFD_vs_CD" "tumor_status_T_vs_C" "treatmentHFD.tumor_statusT"
Is my approach of creating the DESeq object correct ? Also which contrast will give me the desired result?