用途:从 ATAC-seq 片段长度模式中定位核小体中心、占据度与模糊度,并分析调控元件周围的核小体排布。
主要能力:利用 Tn5 切割特征区分片段类别(无核小体 <100 bp、单核小体 180-247 bp、双核小体 315-473 bp);生成 V-plot 诊断定位信息是否可恢复;识别 +1 核小体;估计核小体重复长度(NRL);比较条件间核小体位移、增减与模糊度变化;单细胞场景可用 scprinter。
工具与用法:NucleoATAC 提供单碱基占据度与核小体调用,命令为 nucleoatac run --bed regions.bed --bam sample.bam --fasta genome.fa;ATACseqQC 用 splitGAlignmentsByCut() 分片段类、factorFootprints() 做 TF 侧翼分析、vPlot 画 V-plot;DANPOS3 用 python danpos.py dpos sample.bam 做差异定位,ATAC 需调参如 --smooth_width 80 --width 145;scprinter 用于多尺度与单细胞推断。
适用场景:核小体占据度轨迹、TSS 或 motif 中心 V-plot、条件间差异定位、全基因 +1 核小体调用、单细胞定位、片段大小 QC、NRL 估计。
限制:NucleoATAC 自 2018 年起不再维护,依赖旧版 Python/NumPy,短区域或低深度(<30M 核 reads)易失败;ATACseqQC 无单碱基核小体调用;DANPOS3 默认参数面向 MNase,需 ATAC 调参;scprinter 较新、基准仍在演进且建议 GPU。ATAC 片段大小分析不能确证 H2A.Z,需 ChIP-seq 验证。使用前应核对各工具版本与参数。