We present a new algorithm for the identification of bound regions from
ChIP-seq experiments. Our method for identifying statistically significant
peaks from read coverage is inspired by the notion of persistence in
topological data analysis and provides a non-parametric approach that is robust
to noise in experiments. Specifically, our method reduces the peak calling
problem to the study of tree-based statistics derived from the data. We
demonstrate the accuracy of our method on existing datasets, and we show that
it can discover previously missed regions and can more clearly discriminate
between multiple binding events. The software T-PIC (Tree shape Peak
Identification for ChIP-Seq) is available at
http://math.berkeley.edu/~vhower/tpic.htmlComment: 12 pages, 6 figure

A Barski

AA Bhinge

B Wold

EG Wilbanks

ET Wang

G Carlsson

G Robertson

GR Grimmett

J Rozowsky

Lior Pachter

M Lupien

MB Noyes

PJ Park

R Development Core Team

RK Bradley

S Bhamidi

S Evans

S MacArthur

S Pepke

SN Evans

Steven N Evans

T Barrett

T Laajala

Valerie Hower

WJ Kent

Y Benjamini

Y Zhang

BMC Bioinformatics

English

arXiv

Springer - Publisher Connector

Shape-based peak identification for ChIP-Seq

Crossref

BackgroundThe identification of binding targets for proteins using ChIP-Seq has gained popularity as an alternative to ChIP-chip. Sequencing can, in principle, eliminate artifacts associated with microarrays, and cheap sequencing offers the ability to sequence deeply and obtain a comprehensive survey of binding. A number of algorithms have been developed to call "peaks" representing bound regions from mapped reads. Most current algorithms incorporate multiple heuristics, and despite much work it remains difficult to accurately determine individual peaks corresponding to distinct binding events.ResultsOur method for identifying statistically significant peaks from read coverage is inspired by the notion of persistence in topological data analysis and provides a non-parametric approach that is statistically sound and robust to noise in experiments. Specifically, our method reduces the peak calling problem to the study of tree-based statistics derived from the data. We validate our approach using previously published data and show that it can discover previously missed regions.ConclusionsThe difficulty in accurately calling peaks for ChIP-Seq data is partly due to the difficulty in defining peaks, and we demonstrate a novel method that improves on the accuracy of previous methods in resolving peaks. Our introduction of a robust statistical test based on ideas from topological data analysis is also novel. Our methods are implemented in a program called T-PIC (Tree shape Peak Identification for ChIP-Seq) is available at http://bio.math.berkeley.edu/tpic/

Hower, Valerie

Evans, Steven N

Pachter, Lior

eScholarship - University of California

Background: The identification of binding targets for proteins using ChIP-Seq has gained popularity as an alternative to ChIP-chip. Sequencing can, in principle, eliminate artifacts associated with microarrays, and cheap sequencing offers the ability to sequence deeply and obtain a comprehensive survey of binding. A number of algorithms have been developed to call "peaks" representing bound regions from mapped reads. Most current algorithms incorporate multiple heuristics, and despite much work it remains difficult to accurately determine individual peaks corresponding to distinct binding events. 



Results: Our method for identifying statistically significant peaks from read coverage is inspired by the notion of persistence in topological data analysis and provides a non-parametric approach that is statistically sound and robust to noise in experiments. Specifically, our method reduces the peak calling problem to the study of tree-based statistics derived from the data. We validate our approach using previously published data and show that it can discover previously missed regions. 



Conclusions: The difficulty in accurately calling peaks for ChIP-Seq data is partly due to the difficulty in defining peaks, and we demonstrate a novel method that improves on the accuracy of previous methods in resolving peaks. Our introduction of a robust statistical test based on ideas from topological data analysis is also novel. Our methods are implemented in a program called T-PIC (T ree shape P eak I dentification for C hIP-Seq) is available at http://bio.math.berkeley.edu/tpic/

Evans, Steven N.

Caltech Authors - Main

A: Computation for ChIP-seq and RNA-seq studies.

A: Spectra of large random trees

CB: Alternative isoform regulation in human tissue transcriptomes. Nature

ChIP-seq: advantages and challenges of a maturing technology. Nat Rev Genet

Controlling the false discovery rate: a practical and powerful approach to multiple testing.

Core Team: R: A Language and Environment for Statistical Computing R Foundation for Statistical Computing,

Coverage statistics for sequence census methods.

Developmental roles of 21 Drosophila transcription factors are determined by quantitative differences in binding to an overlapping set of thousands of genomic regions. Genome Biology

Elo L: A practical comparison of methods for detecting transcription factor binding sites in ChIP-seq experiments.

Facciotti MT: Evaluation of Algorithm Performance

FoxA1 Translates Epigenetic Signatures into EnhancerDriven Lineage-Specific Transcription. Cell

Haussler D: The Human Genome Browser at UCSC. Genome Research

Liu XS: Model-based Analysis of ChIP-Seq (MACS). Genome Biology

MB: Binding Site Turnover Produces Pervasive Quantitative Changes in Transcription Factor Binding between Closely Related Drosophila Species. PLoS Biol

MB: PeakSeq enables systematic scoring of ChIP-Seq experiments relative to controls. Nat Biotech

MB: Supplemental code and data for PeakSeq: scoring ChIP-seq experiments relative to controls [http://www.gersteinlab.org/ proj/PeakSeq/].

Myers RM: Sequence census methods for functional genomics.

NCBI GEO: archive for high-throughput functional genomic data. Nucl Acids Res 2009, 37(suppl 1):D885-890 [http://nar.

Probability and real trees, Volume

S: Genome-wide profiles of STAT1 DNA association using chromatin immunoprecipitation and massively parallel sequencing.

SA: A systematic characterization of factors that regulate Drosophila segmentation via a bacterial one-hybrid system. Nucleic Acids Research

Stirzaker DR: Probability and random processes. third edition.

Topology and data. Bull Amer Math Soc (N.S.)

Yekutieli D: The control of the false discovery rate in multiple testing under dependency. Ann Statist

Zhao K: Genomic location analysis by ChIP-Seq.

https://authors.library.caltech.edu/74785/3/12859_2010_4332_MOESM2_ESM.pdf

Shape-based peak identification for ChIP-Seq

Abstract

Similar works

Full text

Available Versions

Springer - Publisher Connector

Crossref

eScholarship - University of California

Caltech Authors - Main