如何进行DAP-seq的数据挖掘,筛选验证位点

从样本准备到寄送公司,每一天都在“祈祷”有个心仪的分析结果,终于在这天随着邮件提示音的响起,收到了分析结果......

分析前工作

爱基在进行数据分析之前,会有两次质控报告反馈给老师们。第一个,基因组DNA的提取质控报告(图1):保证DNA的完整性以及足够的量进行后续的富集亲和纯化;第二个,富集建库报告:构建DNA文库,利用磁珠富集与加完halo Tag标签表达的目的蛋白结合DNA片段,并纯化获得IP文库。这个过程中,为了检测蛋白表达的正常,爱基利用抗体对富集产物进行 WB 检测,同样对于文库也会进行质检(图2)。

图片

图1 DNA提取质控报告

图2 WB结果显示目的蛋白表达正常

分析思路

1. 使用MEME进行motif分析 MEME是一种常用的motif分析工具,可以用于从序列数据中识别出潜在的motif。利用dap-seq输出的peak序列,可以使用MEME进行motif分析。 首先,需要将dap-seq输出的peak序列转换为fasta格式的文件。可以使用bedtools将peak序列提取出来,并将其转换为fasta格式: ``` bedtools getfasta -fi genome.fa -bed peaks.bed -fo peaks.fa ``` 其中,genome.fa为参考基因组序列文件,peaks.bed为dap-seq输出的peak文件,peaks.fa为输出的fasta格式的peak序列文件。 然后,可以使用MEME对peaks.fa进行motif分析: ``` meme peaks.fa -oc meme_output -nmotifs 10 -minw 6 -maxw 20 ``` 其中,meme_output为输出结果的文件夹,-nmotifs指定需要识别的motif数量,-minw和-maxw分别指定motif的最小和最大长度。 2. 使用Homer进行motif分析 Homer是另一个常用的motif分析工具,也可以用于从序列数据中识别出潜在的motif。类似地,利用dap-seq输出的peak序列,可以使用Homer进行motif分析。 首先,需要将dap-seq输出的peak序列转换为bed格式的文件。可以使用bedtools将peak序列提取出来,并将其转换为bed格式: ``` bedtools sort -i peaks.bed > sorted_peaks.bed bedtools merge -i sorted_peaks.bed > merged_peaks.bed awk 'BEGIN{OFS="\t"}{print $1,$2,$3,"peak_"NR,".",$6}' merged_peaks.bed > peaks_homer.bed ``` 其中,peaks.bed为dap-seq输出的peak文件,sorted_peaks.bed和merged_peaks.bed为中间文件,peaks_homer.bed为转换后的bed格式的peak文件。 然后,可以使用Homer对peaks_homer.bed进行motif分析: ``` findMotifsGenome.pl peaks_homer.bed genome_dir homer_output -size 200 -p 8 ``` 其中,genome_dir为参考基因组序列文件夹,homer_output为输出结果的文件夹,-size指定motif的长度,-p指定使用的线程数。 3. 对比分析motif 对于使用不同的工具进行motif分析得到的结果,可以使用Tomtom进行对比分析。Tomtom是一个用于motif比对和聚类的工具,可以帮助用户在已知的motif数据库中搜索相似的motif,并将它们聚类为同一个motif家族。 首先,需要将使用不同工具得到的motif结果转换为meme格式的文件,并将其放入同一个文件夹中,如motif_dir。 然后,可以使用Tomtom进行对比分析: ``` tomtom -o tomtom_output -verbosity 1 -thresh 0.1 -eps -text -min-overlap 5 -dist pearson -no-ssc motif_dir/motif1.meme motif_dir/motif2.meme ``` 其中,tomtom_output为输出结果的文件夹,-thresh指定使用的阈值,-dist指定使用的距离度量方式,-no-ssc表示不使用自身比对,motif1.meme和motif2.meme为需要比对的motif文件。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值