Question: Differential Histone Chip-Seq Analysis Without Replicates
gravatar for Diana
5.9 years ago by
Diana770 wrote:

Hi everyone,

I have histone chip-seq data for 3 developmental stages. I want to find out the differential histone acetyl marks between these 3 stages. I've been trying to work with DiffBind but unfortunately I have no replicates so I cant use that. Does anyone have any suggestions as to how I can do this without replicates?Are there any tools out there or packages in R that deal with this problem. I know there is Bedtools but anything else?

Thanks a lot!!!!

R chip-seq • 3.8k views
ADD COMMENTlink modified 4.8 years ago by Ian5.4k • written 5.9 years ago by Diana770
gravatar for shenli.sam
5.9 years ago by
shenli.sam190 wrote:

You can use diffReps. It is an automated pipeline for ChIP-seq differential analysis. Since you don't have replicates, you can choose G-test (recommended) or Chi-square test. diffReps has implemented them all.

ADD COMMENTlink written 5.9 years ago by shenli.sam190
gravatar for sirenbeserker
5.9 years ago by
sirenbeserker70 wrote:

This page: Get the most out of your chip-seq experiments, contains a list of quantitative tools that look like what you are after. I'm doing similar things with my ChIP-seq data, good luck!

ADD COMMENTlink written 5.9 years ago by sirenbeserker70
gravatar for Ian
4.8 years ago by
University of Manchester, UK
Ian5.4k wrote:

I realise this is an old question, but one i have recently been considering. 

A recent article summarised the problem nicely for me ( and highlights the two main strategies for identifying regions to compare in differential binding analysis.

Use of windows:

E.g. diffReps

This does not require peak calling to have been perfomed.  As histone peak profiles can be narrow and broad and everything in between the use of windows seems to be a better option. 

Use of pre-called peaks:

E.g. MAnorm

If you have confidence in the peak calling results then this is a popular option.  To my mind it is not so useful for diffuse histone marks were peak calling is often not as accurate as with transcription factors or narrow histone marks.

ADD COMMENTlink written 4.8 years ago by Ian5.4k
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: 1453 users visited in the last hour