RNASeq fold change analysis after subtraction of contaminating MEF signal
1
1
Entering edit mode
gpalidwor ▴ 10
@gpalidwor-13703
Last seen 7.3 years ago

I have a set of Mouse RNASeq data for various points of a differentiation time series, with replicates. A known fraction of the cells in each sample are Mouse Embryonic Fibroblasts (MEFs), that fraction varies quite a bit per sample. I have expression data for a pure MEF sample grown under similar conditions.

I'd like to do a fold change analysis between time points, subtracting the MEF expression contamination in such a way that the resulting increased variance per gene is factored into the fold change analysis.

It seems it may be possible to do it within DESeq2 or using svaseq but I can't figure out how. Can anyone recommend a strategy for doing this; I'm really not clear on how to approach it.

Thanks!

(note this is a crosspost from Biostars, as I wasn't able to get a solution there)

deseq2 svaseq • 1.9k views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 6 hours ago
United States

Hi, when you say you have known fraction, do you mean you can give a measured, numeric value for each sample? Or just that you know that some amount will be MEF but you don't know that amount.

ADD COMMENT
0
Entering edit mode

By "known fraction" I mean I can give a measured numeric value for each sample based on the number of cells. 

 

ADD REPLY
0
Entering edit mode

So you could put it in the model as a numeric covariate (you don't do anything special just put it in the design). This however assumes the relationship with expression is log linear (so linear with log expression). You probably want linear with expression though. You can try transforming the MEF variable before putting it in the design, if you expect a certain relationship.

ADD REPLY
0
Entering edit mode

Thanks! I'll give it a try and post this reply on Biostars.

ADD REPLY
0
Entering edit mode

Hello gpalidwor, I am very interested in how you went with this - I am looking to do a very similar thing to 'subtract' RBC contamination from an expression experiment. Can you tell me whether the approach worked? Did you apply a transformation to the MEF variable?

ADD REPLY
0
Entering edit mode

Hello gpalidwor, I am very interested in how you went with this - I am looking to do a very similar thing to 'subtract' RBC contamination from an expression experiment. Can you tell me whether the approach worked? Did you apply a transformation to the MEF variable?

ADD REPLY

Login before adding your answer.

Traffic: 682 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