Question: the bam file used for Diffbind
0
gravatar for mikysyc2016
4 months ago by
mikysyc201620
mikysyc201620 wrote:

Hi, when use diffbind to check ChIP-seq data we need bam file. Do I need to use the bam file which already removed unmapped reads /duplicates/blacklist? And what will be the right range of FRiP for a successful experiment? Thanks.

sequencing chip-seq R • 279 views
ADD COMMENTlink modified 4 months ago by jared.andrews071.4k • written 4 months ago by mikysyc201620
0
gravatar for jared.andrews07
4 months ago by
jared.andrews071.4k
St. Louis, MO
jared.andrews071.4k wrote:

Yes, it is a good idea to remove blacklisted reads before using DiffBind. I usually don't bother removing duplicates unless I have an unusually high percentage of them - particularly for factors with sharp peaks where some duplication is expected.

The FRiP really depends on your sample and how heavily your factor binds. Histone marks get up into the 20-30+% range sometimes if you have a really good enrichment - TFs I usually expect more in the ~15% range, but again, it really depends on your factor and how well your ChIP worked. A general rule I like myself is that I want at least 2x enrichment over my inputs. So if my input has 6% of reads in peaks, I'd want to see at least 12% of reads in peaks in my actual ChIP. But ideally it should be higher, it just depends on how lenient you can afford to be with your samples.

The DiffBind authors have another R package called ChIPQC that is quite good and can tell you the percentage of reads in blacklisted regions, duplication rate, in peaks, general enrichment rate, range of signal, etc. It also makes some nice figures and is easy to use.

ADD COMMENTlink written 4 months ago by jared.andrews071.4k

how about maped and unmaped reads? did you remove use reads mapped bam file?

ADD REPLYlink written 4 months ago by mikysyc201620

Typically I leave them in, as many programs ignore them.

ADD REPLYlink written 4 months ago by jared.andrews071.4k
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: 1423 users visited in the last hour