MACS2、SEACR 与 Genrich
MACS2 是 ChIP-seq 和 ATAC-seq 最常用的峰识别工具,支持窄峰和宽峰;SEACR 为背景低的 CUT&Tag 类数据设计;Genrich 常用于 ATAC-seq。H3K27me3 一般按宽峰处理。
28个分析 Skill · Codar 生信插件
表观基因组学研究不改变 DNA 序列却影响基因表达的调控层面,包括染色质开放、组蛋白修饰、DNA 甲基化和三维基因组结构,典型数据有 ATAC-seq、ChIP-seq、CUT&Tag、亚硫酸氢盐测序和 Hi-C。Codar 的这个插件覆盖峰识别与差异分析、基序与足迹、差异甲基化、增强子靶基因关联和 TAD 分析。
检查测序质量后比对到参考基因组,去除 PCR 重复和低质量读段,ATAC-seq 还要看线粒体读段比例。
评估片段长度分布、TSS 富集和 FRiP 等指标,先判断实验本身是否成功。
按数据类型选择窄峰或宽峰模式,有 input 或 IgG 对照时一并纳入。
在一致峰集上统计计数,按分组做差异可及性、差异结合或差异甲基化分析。
做基序富集与足迹分析,把峰注释到邻近基因,并推断增强子与靶基因的关联。
输出 bigWig 信号轨迹、峰文件、差异区域表、热图和分析脚本。
MACS2 是 ChIP-seq 和 ATAC-seq 最常用的峰识别工具,支持窄峰和宽峰;SEACR 为背景低的 CUT&Tag 类数据设计;Genrich 常用于 ATAC-seq。H3K27me3 一般按宽峰处理。
DiffBind 在一致峰集上计数,再调用 DESeq2 或 edgeR 做差异结合分析,是常用选择;csaw 采用滑动窗口,不依赖预先识别的峰,更适合宽峰信号。两者都需要生物学重复,单样本对比只能做描述。
HOMER 常用于峰区域的已知基序和从头基序富集;chromVAR 在样本或单细胞层面估计各基序的可及性偏差;TOBIAS 利用 Tn5 切割模式做足迹分析,推测转录因子是否实际结合,对测序深度要求较高。
Bismark 是亚硫酸氢盐测序常用的比对和甲基化提取工具;DSS、methylKit 用于识别差异甲基化区域;450K、EPIC 芯片通常用 minfi 或 SeSAMe 预处理。血液和组织样本还要考虑细胞组成的影响。
共 28 项。启用插件后,Codar 会按任务自动选用,不需要逐个记住名称。
在 Codar 左上角选择 Codar 生信。
在插件页找到「表观基因组学」,安装或启用。
用一句话写清数据和目标,确认分析方案后开始执行。
帮我分析这组 ATAC-seq 数据,比较处理组和对照组的差异可及性区域,并做转录因子基序富集。
CUT&Tag 和 CUT&RUN 在原位切割目标区域,背景低、所需细胞少,适合样本有限或研究组蛋白修饰;ChIP-seq 方法成熟,公开数据和抗体资料更多。部分转录因子在 CUT&Tag 中信号不稳定,选择前最好查阅已有数据。
需要。差异可及性分析和 ChIP-seq 一样依赖重复来估计组内变异,通常建议每组至少 2 到 3 个生物学重复,ENCODE 的数据标准也要求重复间有一致性。没有重复时只能做描述性比较,结论需要谨慎。
通常从三个尺度解读:A/B 区室反映开放与紧缩染色质的大尺度分隔,TAD 是内部互作频繁的结构域,染色质环对应增强子与启动子等特定位点之间的接触。能分析到哪个尺度取决于测序深度,染色质环的检测对深度要求最高。