peak-atlas-unpaired

Builds consensus scATAC peak atlases from unpaired RNA and ATAC data via label transfer.

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

SYSTEM DOCUMENTATION & REQUIREMENTS

💡 This Skill requires Signac, Seurat, GenomicRanges, MACS3, and includes scripts (resource) and references (resource) and assets (resource) components.

What problem does it solve? When scRNA-seq and scATAC-seq come from separate cell pools, there are no shared barcodes, so cell identity cannot be read directly off ATAC cells and peak calling cannot be grouped by cell type. This Skill bridges the modalities with a gene-activity matrix and staged CCA label transfer, then calls, merges, tiers, and validates a consensus peak atlas. ## Core Features & Use Cases - Staged label transfer: Transfers coarse lineage labels (R1) then refined subtype labels (R2) from an annotated RNA reference onto ATAC cells via gene-activity space, excluding low-confidence cells from label-resolved peak calling. - Three complementary peak-calling strategies: Calls peaks per coarse cell type (A), per refined subtype × condition (B), and label-free per condition (C) as a transfer-error-robust hedge, then merges them with support voting. - External consensus and tier stratification: Folds published ATAC datasets into a consensus peak source with a 0.25 support threshold and bias gate, and stratifies peaks into priority tiers including protected rare-biology Tier 0 (0a/0b/0c). - Use Case: Given separate scRNA and scATAC experiments on the same immune system, transfer cell-type labels onto ATAC cells, call peaks three ways, gate rare single-strategy peaks on an external-correlation test (r >= 0.4), and produce a validated consensus atlas. ## Quick Start Use the peak-atlas-unpaired skill to transfer labels from my annotated scRNA reference onto my scATAC object and build a tiered consensus peak atlas across conditions.

Frequently Asked Questions about peak-atlas-unpaired

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

FAQPage Schema
How do I call peaks on scATAC data when RNA and ATAC are from different cells?▼

Transfer cell identity from an annotated RNA reference onto ATAC cells using a gene-activity matrix and CCA anchors, then call peaks per transferred label group with Signac's CallPeaks. A label-free per-condition strategy hedges against transfer error.

What is the difference between unpaired and multiome peak atlas workflows?▼

Unpaired data has separate RNA and ATAC cell pools, so identity must be transferred via gene activity before peak calling. Paired 10x Multiome shares barcodes, so cell identity is observed directly and a multiome-specific workflow applies instead.

Why does Signac CallPeaks return identical peaks for every group?▼

CallPeaks with only the idents argument does not subset fragments, so MACS sees all fragments regardless of group. Always pass group.by with the grouping column so Signac filters fragments per group before calling peaks.

How do I validate rare cell-type peaks in an scATAC atlas?▼

Correlate pseudobulk accessibility profiles per cell type against binary accessibility of external datasets, and gate single-strategy rare peaks (Tier 0c) on a named hypothesis with Pearson r >= 0.4. Peaks failing the test are dropped as likely noise.

When should I not use the unpaired peak atlas approach?▼

Do not use it for paired 10x Multiome data where RNA and ATAC come from the same cells, or when no annotated RNA reference exists at all. Without a reference, degrade to ATAC clusters plus external consensus instead.