ChIP-seq / ATAC-seq Peak Calling: Unraveling Gene Regulation

What is Peak Calling?

Peak calling is a computational process used in genomics to identify regions of the genome that are bound by proteins or modified by chemical marks. This technique is crucial for understanding gene regulation, as it helps pinpoint where and when genes are turned on or off.

In ChIP-seq (Chromatin Immunoprecipitation sequencing) and ATAC-seq (Assay for Transposase-Accessible Chromatin with high-throughput sequencing), the process involves capturing DNA fragments that have been bound by proteins of interest, such as transcription factors, or are accessible due to histone modifications. These fragments are then sequenced and analyzed to identify enriched regions.

The Poisson Read-Pileup Model

The core of peak calling algorithms like MACS2 is the Poisson read-pileup model, which predicts how reads should be distributed along the genome if there were no binding sites. This model assumes that reads are randomly distributed according to a Poisson distribution, where the mean number of reads per base pair (the lambda value) represents the background noise.

By comparing the observed read counts at each genomic position with the expected counts from the Poisson model, peak callers can identify regions with significantly higher read counts, indicating potential binding sites or accessible chromatin.

Significance Testing and Background Subtraction

To determine if a region is truly enriched for protein binding or histone modifications, peak callers use statistical tests to compare the observed read counts with the expected background. The local-background significance test used in MACS2 calculates the enrichment of reads over a sliding window compared to the surrounding regions.

This approach helps filter out false positives by ensuring that only regions with significantly higher read counts than their neighbors are considered peaks.

Real-World Applications

ChIP-seq and ATAC-seq peak calling have numerous applications in genomics, including identifying transcription factor binding sites, enhancers, promoters, and other regulatory elements. These insights are vital for understanding gene regulation in health and disease.

For example, researchers can use these data to study how changes in chromatin accessibility or protein binding affect gene expression patterns in different cell types or conditions.

Frequently asked questions

Try it live

Everything above runs in your browser — open ChIP-seq / ATAC-seq Peak Caller (2D) and change the parameters while it is running. Nothing is installed, nothing is uploaded, the whole model lives in one tab.