peaksat: an R package for ChIP-seq peak saturation analysis.

Area of Science:

• Genomics
• Epigenetics
• Bioinformatics

Background:

• Epigenomic profiling assays like ChIP-seq are crucial for mapping chromatin-associated proteins and histone modifications.
• Sequencing depth is a critical factor in experimental design and quality control for these assays.
• Determining optimal sequencing depth can be challenging due to variable requirements across different experimental conditions, targets, and cell types.

Purpose of the Study:

• To develop a tool for estimating target read depth in epigenomic experiments.
• To address the need for informed decisions regarding sequencing data adequacy and peak calling.
• To provide guidance on required read depth for achieving high-quality peak calls.

Main Methods:

• Development of the peaksat R package.
• Analysis of peak saturation curves to estimate target read depth.
• Application of peaksat to ChIP-seq experiments with different histone modifications and cell lines.

Main Results:

• peaksat provides target read depth estimates for epigenomic experiments.
• Distinctive read depth requirements were established for ChIP-seq studies of histone modifications.
• The package successfully estimated target read depths for high-quality peak calls in various cell lines.
• peaksat demonstrated applicability to CUT&RUN and ATAC-seq experiments.

Conclusions:

• peaksat enables researchers to assess if sequencing data depth is adequate for meaningful peak identification.
• The tool quantifies additional reads needed per library when depth is insufficient.
• peaksat is versatile and applicable to various sequencing-based methods involving peak calling.