ChiP-seq小教程收藏chip-seq

ChIP-seq之检峰

2018-07-17  本文已影响1434人  浩渺予怀

ChIP-seq的核心是检峰(Peak Calling),当然现成的工具有很多,其中MACS比较常用,本文以MACS v2为例进行检峰的原理和使用。

一、背景

染色质免疫共沉淀技术(Chromatin Immunoprecipitation,ChIP)也称结合位点分析法,是研究体内蛋白质与DNA相互作用的有力工具,通常用于转录因子结合位点或组蛋白特异性修饰位点的研究。

多种组学手段

二、实验

ChIP实验模式图

一、软件下载

著名代码托管网站github上有最新的开发版本: https://github.com/taoliu/MACS

二、软件使用

(1)双峰模型

检峰过程

MACS2是基于模型(泊松分布)的方法进行检峰的,上述示意图中的模型是双峰模型,目的是为了将比对上的Reads朝3`端偏移(shift),以更准确地得到蛋白-DNA互作的位置。

这里需要注意对于单端数据(SE)需要进行shift以评估平均插入片段大小(Insert Size);若使用的数据为双端(PE)则直接使用真实的Fragment Size进行检峰。

(2)检峰

单链双峰模型

2.1 建立shift模型:

MACS2以固定大小窗口扫描基因组;取前1000个富集程度最好的峰;区分正/负链比对reads,评估d的长度;建立正/负链模型。

2.2 检测显著峰

Reads向中心Shift 1/2d的距离;以长度2d的窗口扫描全基因组得到显著富集的候选峰(基于泊松分布); 合并重叠的峰。

# d:插入片段长度

三、峰型

(1) 窄峰 -转录因子等

专转录因子

(2) 混合峰 - 聚合酶II等

聚合酶II

(3)宽峰 -组蛋白等

组蛋白

四、如何使用

MACS2非常适合转录因子和组蛋白修饰类型的检峰。

使用MACS2时需要根据不同峰型进行参数的调整,软件说明文档参数如下:

常用参数:macs2 callpeak -t ChIP.bam -c Control.bam -f BAM -g hs -n test  -q 0.01

宽峰参数:macs2 callpeak -t ChIP.bam -c Control.bam --broad -g hs --broad -cutoff1

其中,-t为ChIP样本;-c为control样本即input;-f为输入文件的格式,一般使用BWA、Bowtie2比对后的BAM文件;-g有效基因组大小,MACS2内置了一些物种如人类是hs;-n为样本名称;-q为检峰显著性阈值; --broad为检测宽峰的阈值; -cutoff为筛选fold-enrichment的阈值。

针对转录因子类型为了得到结合位置的峰顶可以添加 --call-summits,另外--extsize具体数值可以使用macs2 predictd进行评估(例如评估值为200),检峰参数如下:

macs2 callpeak -t ChIP.bam -c Control.bam -f BAM -g hs -n test --call-summits --nomodel --extsize 200 -q 0.01

针对组蛋白类型,需要合并邻近的峰形成宽峰,检峰参数如下:

macs2 callpeak -t ChIP.bam -c Control.bam -f BAM -g hs -n test --broad --nomodel --extsize 173  -cutoff 2 -q 0.01

对于双端数据,-f参数为BAMPE,这时检峰所使用的插入片段长度为实际的平均插入片段长度,不需要建模(使用--nomodel,同时--extsize失效)

检峰之后就是检测基序(motif)了,链接:检测基序

五、其他软件

参考如下文章:

检峰软件比较 检峰精确性比较(NRSF基序)

六、参   考

[1] https://www.encodeproject.org/

[2] https://en.wikipedia.org/wiki/ChIP-sequencing

[3] https://www.illumina.com/techniques/sequencing/dna-sequencing/chip-seq.html

上一篇下一篇

猜你喜欢

热点阅读