現(xiàn)在ChIP-seq的數(shù)據(jù)基本是最常見的測序數(shù)據(jù)類型之一活合，主要有Transcription factor ChIP-seq和Histone ChIP-seq作谭。前者是看轉(zhuǎn)錄因子的結(jié)合位置橱乱，后者是組蛋白修飾發(fā)生的位置。下面分享一下一般流程我抠。
- 質(zhì)控 (quality control)
首先要看一下ChIP-seq數(shù)據(jù)的質(zhì)量拙已，數(shù)據(jù)的信號最好比background要強(qiáng)很很多决记。一般要有control，這樣call peaks更準(zhǔn)確可信倍踪， control主要有Input DNA 和 IgG兩種系宫，前一種更常用。
檢測質(zhì)量的一些方式：
1). peaks中reads的數(shù)量建车，如果peaks的reads普遍較少扩借，則質(zhì)量一般。
2). peaks信號高癞志，背景低往枷。
3). 測序深度深 框产。
4). Diverse library (與重復(fù)duplications有關(guān)凄杯，如下圖)
4). 有重復(fù)并且與重復(fù)之間相似性較高…
……
做質(zhì)控的軟件方法：
1). ChIPQC (T Carroll, Front Genet, 2014.)
2). SPP package - Unix/Linux (PV Karchenko, Nature Biotechnol, 2008.)
3). ENCODE中的標(biāo)準(zhǔn)流程
- 序列比對 (mapping of fastq)
序列比對一般用BWA或者Bowtie2，兩者效果差不多秉宿。BWA的bwa samse（單端數(shù)據(jù)）和bwa sampe (雙端數(shù)據(jù)) 跑的速度比較慢戒突，但是效果很不錯(cuò)，用法如下：
bwa index reference.fa # 建立索引 -p可設(shè)置前綴描睦，不設(shè)置前綴就是reference.fa膊存。
# 單端數(shù)據(jù)
bwa aln -t 8 reference.fa test.fq.gz > test.sai
bwa samse -n 10 reference.fa test.sai test.fq.gz > test_se.sam
# 雙端數(shù)據(jù)：
bwa aln reference.fa test_reads1.fq > test1.sai
bwa aln reference.fa test_reads2.fq > test2.sai
bwa sampe reference.fa test1.sai test2.sai test_reads1.fq test_reads2.fq > test_pe.sam
BWA的mem，速度很快：
bwa mem reference.fa reads.fq > test_se.sam # 單端
bwa mem reference.fa read1.fq read2.fq > test_pe.sam # 雙端
bowtie2的用法：
bowtie2-build reference.fa index # 創(chuàng)建index
bowtie2 -p 8 -x index -U test_read.fq -S test_se.sam # 單端比對
bowtie2 -p 8 -x index -1 test_read1.fq -2 test_read2.fq -S test_pe.sam # 雙端比對
效果個(gè)人感覺差不多。
- 去除重復(fù) (remove duplicates)
由于PCR實(shí)驗(yàn)存在不可避免的實(shí)驗(yàn)誤差隔崎，所以會存在重復(fù) (duplicates)今艺。比如兩條不同的reads，起止位置完全一致爵卒。比如：
其中第二條已經(jīng)被picard標(biāo)注出來了虚缎。被標(biāo)注的第二列flag會加1024。
去重的軟件中samtools rmdup （基本已不用）钓株，samtools markdup（更新后的）和picard最常用实牡。rmdup效果不怎么好，而且如果有遇到相同位置的reads轴合， 會優(yōu)先選擇質(zhì)量高的那一條read创坞。picard與samtools markdup效果相似（仿佛調(diào)用的同一個(gè)？并不確定）受葛。都可以標(biāo)記重復(fù)题涨，也可以選擇直接去掉。以下是用法：
samtools markdup -@ 8 -r test.bam filter_test.bam # -r是直接去掉重復(fù)总滩，不加是直接標(biāo)記
picard去重有三種方式可選携栋，在DUPLICATE_SCORING_STRATEGY參數(shù)中，分別是SUM_OF_BASE_QUALITIES, TOTAL_MAPPED_REFERENCE_LENGTH和RANDOM咳秉。即當(dāng)有重復(fù)時(shí)分別選擇留下總堿基質(zhì)量最高的婉支、匹配上參考基因組最長的和隨機(jī)。
picard MarkDuplicates I=test.bam O= filter_test.bam M=dup_metrics.txt REMOVE_DUPLICATES=true
在call peak之前需要去除blacklisted regions澜建，這些區(qū)域可能是有問題的向挖，詳解及下載可參考http://www.reibang.com/p/76edbc772500
#Remove alignments in Encode blacklisted regions :
intersectBed -v -abam in.bam -b ENCFF001TDO.bed > out.bam
- peak calling
peaks是reads信號比較強(qiáng)的區(qū)域，也就是我們找到的轉(zhuǎn)錄因子或者組蛋白修飾最有可能結(jié)合的地方炕舵。call peaks仍然有不少軟件何之，比較常用的是MACS2和Hotspot2。
示例：
macs2 callpeak -t test.bam -c control.bam -f BAM -g hs -n test -B -q 0.01
針對不同的數(shù)據(jù)考慮用不同的參數(shù)咽筋。
- 下游分析 (downstream analysis)
分析完之后下游可以做的事情很多溶推，視情況而定〖楣ィ可以同時(shí)分析DNase-seq或者ATAC-seq的數(shù)據(jù)蒜危，看轉(zhuǎn)錄因子與染色質(zhì)開放區(qū)的關(guān)系；或者Homer等工具注釋peaks睹耐，看不同轉(zhuǎn)錄因子/組蛋白修飾之間的關(guān)系辐赞，或者分析TF的target gene。也可以用MEME等做motif分析硝训。
歡迎關(guān)注响委！
