rocco 1.4.1

Robust ATAC-seq Peak Calling for Many Samples via Convex Optimization

ROCCO: [R]obust [O]pen [C]hromatin Detection via [C]onvex [O]ptimization

What

ROCCO is an efficient algorithm for detection of "consensus peaks" in large datasets with multiple HTS data samples (namely, ATAC-seq), where an enrichment in read counts/densities is observed in a nontrivial subset of samples.

Input/Output

• Input: Samples' BAM alignments or BigWig tracks
• Output: BED file of consensus peak regions

Note, if BigWig input is used, no preprocessing options can be applied at the alignment level.

How

ROCCO models consensus peak calling as a constrained optimization problem with an upper-bound on the total proportion of the genome selected as open/accessible and a fragmentation penalty to promote spatial consistency in active regions and sparsity elsewhere.

Why

ROCCO offers several attractive features:

1. Consideration of enrichment and spatial characteristics of open chromatin signals
2. Scaling to large sample sizes (100+) with an asymptotic time complexity independent of sample size
3. No required training data or a heuristically determined set of initial candidate peak regions
4. No rigid thresholds on the minimum number/width of supporting samples/replicates
5. Mathematically tractable model permitting worst-case analysis of runtime and performance

Example Behavior

Input

• ENCODE lymphoblastoid data (BEST5, WORST5): 10 real ATAC-seq alignments of varying TSS enrichment (SNR-like)
• Synthetic noisy data (NOISY5)

We run twice under two conditions -- with noisy samples and without

rocco -i *.BEST5.bam *.WORST5.bam -g hg38 -o rocco_output_without_noise.bed
rocco -i *.BEST5.bam *.WORST5.bam *.NOISY5.bam -g hg38 -o rocco_output_with_noise.bed

Output

Comparing each output file:

• ROCCO effectively separates true signal from noise across multiple samples
• ROCCO is robust to noisy samples (e.g., output unaffected by inclusion of NOISY5 inputs)
• ROCCO offers high resolution separation of enriched regions

Paper/Citation

If using ROCCO in your research, please cite the original paper in Bioinformatics (DOI: btad725)

 Nolan H Hamilton, Terrence S Furey, ROCCO: a robust method for detection of open chromatin via convex optimization,
 Bioinformatics, Volume 39, Issue 12, December 2023

Documentation

For additional details, usage examples, etc. please see ROCCO's documentation: https://nolan-h-hamilton.github.io/ROCCO/

Installation

PyPI (pip)

pip install rocco --upgrade

Build from Source

If preferred, ROCCO can easily be built from source:

• Clone or download this repository

  git clone https://github.com/nolan-h-hamilton/ROCCO.git
  cd ROCCO
  python setup.py sdist bdist_wheel
  pip install -e .