Question: DESeq2 for Cut&Run and differential enrichment testing
0
gravatar for MPerch
5 months ago by
MPerch0
MPerch0 wrote:

Hi there,

I wish to compare differential enrichment of a histone modification (H3K9me2) across 5 replicates of control and treated samples, for each of two cell types. This modification is very broad and not amenable to peak analysis. Therefore I have used featurecounts to count up reads falling in gene-bodies for each gene and sample, and now have a counts matrix of K9me2 enrichment at each gene.

I would like to use this counts matrix in DESeq2 for differential enrichment analysis, as if treating it as RNA-seq in essence. However one issue I have is that for Cut&Run and for analysis of global changes in enrichment, I first normalize the data to the number of mapped reads generated from S.cerevisiae spike-in, within each sample. This means my data table now contains spike-in normalized values and not raw counts, which I know DESeq2 normally requires for DEseq(). If I round my values so that no error is thrown, is it OK to continue to use DEseq() to find genes with differential enrichment of the modification (and what would you recommend to set for the betaPrior?) If not advisable, is there a better way to allow normalization to spike-in values plus differential testing with DESeq2?

Many thanks for any help.

 

deseq2 cut&run • 174 views
ADD COMMENTlink modified 5 months ago by Michael Love23k • written 5 months ago by MPerch0

Hi, I have a question which is not related with your question. I wonder what do you mean "normalize the data to the number of mapped reads generated from S.cerevisiae spike-in"?   Did you divide the ChIP read count by the corresponding spike in mapped reads?

ADD REPLYlink written 5 months ago by Shuang He0
Answer: DESeq2 for Cut&Run and differential enrichment testing
0
gravatar for alex.gos90
5 months ago by
alex.gos9010
alex.gos9010 wrote:

Hi,

I was recently very happy to find RUVSeq (https://bioconductor.org/packages/release/bioc/html/RUVSeq.html) for normalisation prior differential analysis. They have a method specifically developed for spike-in normalisation (RUVg), and provide further advice how to proceed with DEseq2.

Maybe this helps you. 

Best,

Alex

ADD COMMENTlink written 5 months ago by alex.gos9010

I should add that I have used it for ATACseq data, but I do not know how it performs with broad ChIPseq. 

ADD REPLYlink written 5 months ago by alex.gos9010
Answer: DESeq2 for Cut&Run and differential enrichment testing
0
gravatar for Michael Love
5 months ago by
Michael Love23k
United States
Michael Love23k wrote:

Do not normalize the data for sequencing depth and provide to DESeq2. This point is made frequently in our documentation.

You can instead do:

dds <- estimateSizeFactors(dds, controlGenes=...)
dds <- DESeq(dds)

 

ADD COMMENTlink written 5 months ago by Michael Love23k
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 16.09
Traffic: 109 users visited in the last hour