I have bulk RNAseq data for patients with a disease, and for control patients.
Cell deconvolution (using a single-cell RNA-seq reference - a separate study) shows that cell type proportions are changing a lot in cases vs controls. In addition, differential gene expression analysis (DESeq2) highlights many genes which vary in cases vs controls.
Is there a way to 'deconfound' the differential gene expression measures by the known changes in the cell type proportions? One option I thought of is to include the proportion of each cell type as covariates in DESeq2 but I wasn't 100% sure if that was valid.
The ultimate aim is to be able to say whether gene expression changes in disease vs controls are happening in specific cell types. Do you think this is possible ?
Thank you very much!
Hey Peter, I was just thinking about doing the same in my analysis, but unsure whether it would work.
Did we end up going forward with this? Any advice you could share?