GenoGAM: genome-wide generalized additive models for ChIP-Seq analysis

Stricker, Georg; Engelhardt, Alexander; Schulz, Daniel; Schmid, Matthias; Tresch, Achim und Gagneur, Julien (2017): GenoGAM: genome-wide generalized additive models for ChIP-Seq analysis. In: Bioinformatics, Bd. 33, Nr. 15: S. 2258-2265

Abstract

Motivation: Chromatin immunoprecipitation followed by deep sequencing (ChIP-Seq) is a widely used approach to study protein-DNA interactions. Often, the quantities of interest are the differential occupancies relative to controls, between genetic backgrounds, treatments, or combinations thereof. Current methods for differential occupancy of ChIP-Seq data rely however on binning or sliding window techniques, for which the choice of the window and bin sizes are subjective. Results: Here, we present GenoGAM (Genome-wide Generalized Additive Model), which brings the well-established and flexible generalized additive models framework to genomic applications using a data parallelism strategy. We model ChIP-Seq read count frequencies as products of smooth functions along chromosomes. Smoothing parameters are objectively estimated from the data by cross-validation, eliminating ad hoc binning and windowing needed by current approaches. GenoGAM provides base-level and region-level significance testing for full factorial designs. Application to a ChIP-Seq dataset in yeast showed increased sensitivity over existing differential occupancy methods while controlling for type I error rate. By analyzing a set of DNA methylation data and illustrating an extension to a peak caller, we further demonstrate the potential of GenoGAM as a generic statistical modeling tool for genome-wide assays.

• Dokumententyp: Zeitschriftenartikel
• Fakultät: Chemie und Pharmazie > Department Biochemie
• Themengebiete: 500 Naturwissenschaften und Mathematik > 540 Chemie
• ISSN: 1367-4803
• Sprache: Englisch