Differential Gene Expression analysis in heterogeneous tissues
Entering edit mode
12 months ago
txema.heredia ▴ 130


I must run DGE analysis on bulk RNA-seq samples from abdominal aorta. Unfortunately, during the initial QC steps I detected a large variability in two lipid processing genes pervasive across all experimental conditions. This suggests that some of the sequenced samples had a larger amount of adipose tissue still attached to them. Thus, some samples' expression profile might be dominated by adipose tissue instead of vascular.

How should I perform the DEG analysis correcting by this? Can I simply add to the metadata the DESeq2-normalized expression levels of these 2 genes and run DESeq2 using as my design ~ adipose_score * grouping_variable ? where adipose_score is a continuous variable and grouping_variable is a factor with condition+genotype.

I have read about using RUVseq here. However, this dataset contains 3 experiments KOing a certain gene.

  • WT vs whole-body KO
  • WT vs KO only in celltype 1
  • WT vs KO only in celltype 2

The proportions of celltype1 & 2 do also vary from sample to sample (even within the same experiment), although we cannot detect a clear pattern or marker for that. I worry that, by running RUVseq to remove unwanted variation, the software is going to remove also the celltype-of-interest variation instead of only the adipose-specific one.

What would you recommend me to do here? And how should I interpret the results/model?


RNA-seq deseq2 • 388 views

Login before adding your answer.

Traffic: 2569 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6