Question: Remove batch effect for differential expression
1
gravatar for Ron
2.8 years ago by
Ron970
United States
Ron970 wrote:

Hi,

I want to do differential expression between two groups(two batches).For clustering analysis of these samples from two groups,I did batch normalization using COMBAT (sva package)on expression values so as to remove batch effects.

But for differential expression,should I do batch normalization as well before proceeding, as DEseq expects raw counts ? Or is there any way DEseq bioconductor takes into account for different batches and does normalization,before doing differential expression?

Any suggestions,please?

Thanks,

Ron

rna-seq deseq next-gen R • 3.1k views
ADD COMMENTlink modified 2.8 years ago by WouterDeCoster40k • written 2.8 years ago by Ron970
2
gravatar for WouterDeCoster
2.8 years ago by
Belgium
WouterDeCoster40k wrote:

Assuming you did a proper design and treatments are distributed over batches (not all treatment A in batch 1 and treatment B in batch 2 but mixed) you can specify the batch in the design formula and this will be taken into account by DESeq2 and similar tools.

An example of a discussion about a design like that can be found at correcting for a batch in DESeq2 and https://support.bioconductor.org/p/76099/ (among many more)

ADD COMMENTlink written 2.8 years ago by WouterDeCoster40k

Hi WouterDeCoster,

I read this post,that I can use Combat batch normalization on Counts and then used the normalized values for differential expression with limma (not EdgeR and DESeq). What do you think on this?

ADD REPLYlink written 2.7 years ago by Ron970

Don't use Combat. Go by Wouter's advice, i.e., to include batch in the design model / formula.

If you need even more advice, then take that of one of the core developers of Limma himself:

ADD REPLYlink modified 16 months ago • written 16 months ago by Kevin Blighe47k

Hi WouterDeCoster,

How to deal with batch effects in this situation (all treatment A in batch 1 and treatment B in batch 2)? Specifying batch in the design formula gives error. Is it even possible to remove in that case?

ADD REPLYlink written 10 weeks ago by kunal3680
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 750 users visited in the last hour