Peak Calling For Non-Chip Data?

I want to call peaks for a nucleosome density signal from the ENCODE dataset. Is it ok to use a package like MACS to do so, or are there underlying assumptions to such packages that I am overlooking?

Thanks!

EDIT: Upon further inspection, the data I want to call peaks for is in the .bigWig file format. Any way to call peaks on such a format? (i.e. converting it to a .bed file?)

macs peak-calling encode

• 3,931 views • link

updated 13.1 years ago by Ming Tommy Tang

• written 13.1 years ago by vedtopkar •

Why don't you want to use ENCODE's peak calls ?

• 0 views • link

written 13.1 years ago by dario.garvan

As far as I can tell, peak calls are not published for the Nucleosome position dataset. (I'd love to be proven wrong, however!)

• 0 views • link

written 13.1 years ago by vedtopkar •

2 answers

I'd take a look at this very recent publication (http://www.nature.com/nbt/journal/vaop/ncurrent/full/nbt.2596.html) which I haven't tried yet but am about to do.

• 115 views • link

written 13.1 years ago by Mikael Huss

This is perfect, thanks!

• 0 views • link

written 13.1 years ago by vedtopkar •

check if the ones below are what you are looking for

http://code.google.com/p/danpos/ http://code.google.com/p/dantools/

Bioconductor packages http://bioconductor.org/packages/release/bioc/html/nucleR.html http://bioconductor.org/packages/release/bioc/html/PING.html

• 116 views • link

written 13.1 years ago by Ming Tommy Tang

Very helpful, thanks!

• 0 views • link

written 13.1 years ago by vedtopkar •

More posts like this

Define enhancer regions
written 2.6 years ago by Faiza •

Hi! I need to get enhancer regions from publicly available data e.g ENCODE. Encode contain Histone marks bigwig file which is Signal p value and …
Encode project signal p-value
written 3.2 years ago by jomo018

In TF ChIP-Seq experiments, ENCODE project provides two types of signal BigWig files: fold change over control and signal p-value. I'm unable to locate the …
Using DANPOS2 to determine nucleosome positions
written 3.5 years ago by elisheva

Hi everyone. I want to determine the nucleosome positions from ENCODE data ([ENCFF000VME][1]). The file is a signal file in Bigwig format - so there …
ATACseq: Extract Mono-nucleosome information from di- and tri-nucleosome data
written 4.4 years ago by Deep •

Hi, My objective is to split di- and tri-nucleosome ATAC-seq data points into mono-nucleosome and then consolidate the total ATAC-seq (mono-, di- and tri-) nucleosome …
Differential Peak Analysis from ENCODE Pipeline normalized outputs
written 6.8 years ago by sjquon •

I processed my ChIP-seq and ATAC samples using the ENCODE pipeline (https://github.com/ENCODE-DCC/chip-seq-pipeline2 and https://github.com/ENCODE-DCC/atac-seq-pipeline). These pipelines gives outputs such as aligned bam files, bigwig files …
Bam to BigWig loss of data at MACS peaks (nothing in the BigWig why?)
written 8.3 years ago by contreras.xavier •

Hi all, We have some Chip-Seq data (from Human HeLa cells, 30 million reads) for our factor of interest. If I check my data on …
ENCODE DNAse 1 HyperSensitivity Data - Signal of unique reads vs Base Overlap Signal vs Raw Signal
written 9.4 years ago by morovatunc

Hi, I would like to know the difference between base overlap signal, raw signal, signal of unique peaks output types which belong to Encode DNAse …
tag density v.s. binding affinity v.s. signal value
written 11.7 years ago by liu4gre

I am wondering if it is OK to use signal value stored in narrowPeak format (ENCODE) to represent the binding affinity of a specific TF? …
bigwig to peak
written 12.1 years ago by Nicolas Rosewick

Hi, I've a bigwig containing a chip-seq signal for a histone mark (H3K4me1) and I want to produce a bed file containing the peaks (the …
Teaching Yourself Bioinformatics/Encode?
written 13.1 years ago by vedtopkar •

Hi all, I am an undergraduate with a fair bit of programming experience and am beginning to dip my toes into bioinformatics and the ENCODE …