PeakSeq

From GersteinInfo

(Difference between revisions)
Jump to: navigation, search
Line 1: Line 1:
-
PeakSeq
 
-
PeakSeq: Systematic Scoring of ChIP-Seq Experiments Relative to Controls
+
 
 +
'''PeakSeq: Systematic Scoring of ChIP-Seq Experiments Relative to Controls'''
 +
 
Rozowsky J, Euskirchen G, Auerbach R, Zhang Z, Gibson T, Bjornson R, Carriero N, Snyder M, Gerstein M
Rozowsky J, Euskirchen G, Auerbach R, Zhang Z, Gibson T, Bjornson R, Carriero N, Snyder M, Gerstein M

Revision as of 09:21, 16 September 2010


PeakSeq: Systematic Scoring of ChIP-Seq Experiments Relative to Controls

Rozowsky J, Euskirchen G, Auerbach R, Zhang Z, Gibson T, Bjornson R, Carriero N, Snyder M, Gerstein M

Chromatin Immunoprecipitation followed by tag sequencing (ChIP-Seq) using high-throughput next-generation instrumentation is replacing ChIP-chip for genome-wide mapping of sites of transcription factor binding and chromatin modification, especially for mammalian genomes. Here we develop a methodology for identifying punctate binding sites in ChIP-Seq experiments based on their characteristics. In particular, we produce two deeply sequenced datasets for human RNA polymerase II and STAT1 with matching input DNA controls. In these sets, we observe that signal peaks, corresponding to sites of potential binding, are strongly correlated with peaks in input DNA that likely reveal features of open chromatin. Based on these observations we develop a two-pass approach for scoring ChIP-Seq data relative to controls. The first pass identifies putative binding sites and compensates for variation in the mappability of sequences across the genome. The second pass filters out sites that are not significantly enriched compared to the normalized input DNA and computes a precise enrichment and significance. Using our scoring approach we investigate the design of an optimal ChIP-Seq experiment. We examine the number of identified binding sites as a function of sequencing depth (i.e. saturation) and the value of multiple replicas. In particular, we find that little additional biological information is gained from more than two replicas.


Mappability Map

Mappability Map Code

Maps:

Map for C. elegans

Map for D. melanogaster

Map for H. sapiens"

Map for M. musculus


ChIP-Seq Scoring

Preprocessing Code

C Code for ChIP-Seq Scoring

Perl Code for ChIP-Seq Scoring


Scored results:

Pol II Results

STAT1 Results


Raw ChIP-Seq Sequence Data:

Pol II Sequence Data

STAT1 Sequence Data


Contact Joel Rozowsky at joel DOT rozowsky AT yale DOT edu.

Personal tools