Question: Files to compare ChIP vs Input on genome browser: after alignment or peak calling?
0
gravatar for salamandra
3 months ago by
salamandra130
salamandra130 wrote:

When visualizing ChIP-seq data on genome browser should we use files obtained upon read alignment (bam files) or instead use the files treat_qvalue.bdg and control_lambda.bdg (from macs2) provided upon peak calling?

chip-seq macs2 genome browser • 332 views
ADD COMMENTlink modified 3 months ago by Friederike1.8k • written 3 months ago by salamandra130
1
gravatar for Satyajeet Khare
3 months ago by
Satyajeet Khare1.2k
Pune, India
Satyajeet Khare1.2k wrote:

You can use treat_pileup.bdg.gz files and peaks.bed file. Bam files are not normalized to sequencing depth and are bigger in size.

ADD COMMENTlink modified 3 months ago • written 3 months ago by Satyajeet Khare1.2k

so treat_pileup.bdg is for chip, but then for input do i use peaks.bed? cause I though .bed didn't give 'read intensity'...

ADD REPLYlink written 3 months ago by salamandra130

peak.bed file tells you which peaks are "real". There is only one peaks.bed file and it contains the co-ordinates of the statistically significant peaks.

ADD REPLYlink modified 3 months ago • written 3 months ago by Satyajeet Khare1.2k

thank you, but I wanted chip and input to be different tracks on genomebrowser, which should I use for chip and which for input? do you know?

ADD REPLYlink written 3 months ago by salamandra130
1

You can use treat_pileup.bdg.gz for ChIP and control_lambda.bdg.gz file for input.

ADD REPLYlink modified 3 months ago • written 3 months ago by Satyajeet Khare1.2k

ok then if wanted to show only one track with the difference in intensity of ChIP minus Input I would use one of these two files: treat_pvalue.bdg or treat_qvalue.bdg files? or are those two files the ChIP intensity, but just for the significant peaks?

ADD REPLYlink written 3 months ago by salamandra130
1

for this, you could use deepTools bamCompare

ADD REPLYlink written 3 months ago by Friederike1.8k

but still in that case, bamCompare doesn't select significant peaks, i'd still have to select those significant peaks based on macs2 peak.bed file, as sugested by Satyajeet?

ADD REPLYlink written 3 months ago by salamandra130

also, do we still have to remove duplicates when using bamCompare?

ADD REPLYlink written 3 months ago by salamandra130

yes, deepTools does not do any statistical assessment, that's what the MACS2 output is for.

ADD REPLYlink written 3 months ago by Friederike1.8k

You can show reads only from treat_pileup.bdg.gz file for regions where you observe significant peak. The peak is called at a region because no significant reads were observed in Input sample in that locus.

ADD REPLYlink written 3 months ago by Satyajeet Khare1.2k
1
gravatar for Friederike
3 months ago by
Friederike1.8k
United States
Friederike1.8k wrote:

deepTools' bamCoverage will allow you to generate bigwig files that are normalized for sequencing depth, i.e., they will be smaller in size than the bedgraph files and you can somewhat judge the differences in the peak magnitudes between different samples visually.

You could run bamCoverage for the ChIP and input BAM files separately and you will obtain bigwig files that you can easily load into a genome browser.

ADD COMMENTlink written 3 months ago by Friederike1.8k

But in that case would I have to first remove duplicates? I'm asking because I am not removing sequence duplicates after alignment to genome, as during peak calling MACS2 automatically removes them. But if use bigwigs got directly from bam files they will have duplicates right?

ADD REPLYlink written 3 months ago by salamandra130
1

if the duplicates are marked as such, e.g. after PICARD's MarkDuplicates, they can be ignored using one of the many options of bamCoverage

ADD REPLYlink written 3 months ago by Friederike1.8k
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: 1857 users visited in the last hour