Question: Recommendation for using limma-voom workflow for ChIPseq analysis. Or other tool for capturing global histone mark difference?
gravatar for bioinfouser
8 months ago by
bioinfouser80 wrote:

Hello awesome community!

I just wanted to know your expert opinion, whether limma-voom protocol is suitable for analysing ChIPseq data(with macs2 peaks and raw read counts generated from featureCounts). Does it assume that majority of the expression will not change among groups, just like DESeq2? I am having trouble with my ChIPseq data analysis as I am expecting global change of histone modification, which doesn't follow DESeq2's assumption. DiffBind seemed to be working to some extent to capture the global change with its simple BAM file library size normalization, but my experimental design is too complex to block multiple unwanted factors and get proper differential binding regions.... I was thinking of bypassing DESeq2's default sizeFactor normalization and feed simple BAM library size factors(just like DiffBind) and then doing estimateDispersions and nbinomWaldTest functions to perform differential analysis. I know in this case, running a spike-in control would have been optimal, but it is not possible right now.

I would really appreciate if you can comment on which sets of tool should be optimal to capture global histone modification differences when spike in control is not available.

chip-seq • 267 views
ADD COMMENTlink written 8 months ago by bioinfouser80

Check the csaw package. It manual contains extensive suggestions for analysis of ChIP-seq data.

ADD REPLYlink modified 8 months ago • written 8 months ago by ATpoint44k

Thank you! I will look into it. My current knowledge tells me, it is a window based tool, rather then "analysing predefined peaks(by macs2, for example)" that I am doing. Could you comment if csaw is more advantageous than defined peak based analysis?

ADD REPLYlink written 8 months ago by bioinfouser80
Please log in to add an answer.


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