iterative-peak-merging

Merge MACS2 summit peaks into fixed-width 501bp reference sets.

1|Updated Nov 20, 2025
One-click install
npx skills add https://github.com/tony-zhelonkin/SciAgent-toolkit --skill iterative-peak-merging
Or copy as Structured Prompt for Agent
Please help me install this Agent Skill.
Skill: iterative-peak-merging
Source: https://github.com/tony-zhelonkin/SciAgent-toolkit/tree/main/skills/iterative-peak-merging
Command: npx skills add https://github.com/tony-zhelonkin/SciAgent-toolkit --skill iterative-peak-merging

SYSTEM DOCUMENTATION & REQUIREMENTS

What problem does it solve?

Iterative overlap peak merging resolves the challenge of consolidating per-sample ATAC-seq peak sets into a single, consistent reference by avoiding daisy-chaining and preserving cell-type-specific peaks through a fixed-width peak framework.

Core Features & Use Cases

  • Fixed-width peak construction: Generates 501bp peaks centered on summits for uniform downstream analysis.
  • Tiered merging: Performs within-group then across-group merging to preserve reproducibility and capture group-specific accessibility.
  • Compatibility: Designed to work with MACS2 summit files, reproducibility rules, and downstream tools in ATAC-seq pipelines (e.g., ArchR, chromVAR, DESeq2).

Quick Start

Run the createIterativeOverlapPeakSet.R script with your MACS2 summit files and metadata to generate the final merged peak set

Frequently Asked Questions about iterative-peak-merging

High-intent search queries and answers about installing and using this skill.

FAQPage Schema
How do I merge ATAC-seq peaks from multiple samples into a single reference set?

To merge ATAC-seq peaks from multiple samples, you can use iterative overlap merging to consolidate per-sample MACS2 summit files into a single fixed-width 501bp reference set, avoiding daisy-chaining while preserving cell-type-specific peaks for downstream quantification.

Why does iterative overlap peak merging use fixed-width 501bp peaks?

Fixed-width 501bp peaks are used in iterative overlap merging to center peaks on MACS2 summits, ensuring uniform downstream analysis and preventing daisy-chaining artifacts that occur with variable-width peak consolidation across multiple samples.

What is the best way to consolidate MACS2 summit files for downstream ATAC-seq quantification?

The best way to consolidate MACS2 summit files is through tiered merging, which performs within-group then across-group merging to preserve reproducibility and capture group-specific accessibility before downstream quantification with tools like DESeq2.

Does this iterative peak merging approach work with ArchR and chromVAR pipelines?

Yes, this iterative peak merging approach is designed for compatibility with downstream ATAC-seq pipeline tools including ArchR, chromVAR, and DESeq2, ensuring the consolidated 501bp peak set integrates smoothly with standard quantification workflows.

When do I need iterative within- and across-group merging for ATAC-seq peaks?

You need iterative within- and across-group merging when consolidating ATAC-seq peaks across cell types or conditions, as it applies reproducibility rules to preserve group-specific accessibility and ensure robust peak detection in the final reference set.

Can I use this peak merging method without reproducibility rules for my ATAC-seq data?

Reproducibility rules are integral to this peak merging method, ensuring that the tiered within-group and across-group merging process accurately preserves cell-type-specific accessibility and generates a robust 501bp peak set for downstream analysis.