工欲善其事必先利其器 1MACS3 MACS3 是由Liu Tao 主导开发的一款ChIP-Seq 数据分析工具,作为 MACS 系列软件的第三代版本,MACS3 继承并扩展了前两代软件的功能,成为 ChIP-Seq 自动片段大小估计:MACS3 能自动估计 DNA 片段的大小和位移大小,这对于优化峰值检测至关重要,有助于更准确地定位蛋白质-DNA 结合区域。 高效的数据处理 利用对照样本:MACS3 能够利用对照样本,通过比较实验样本和对照样本,减少背景噪音,提高峰值检测的特异性。 批量处理能力:MACS3 能处理大规模数据集,适用于高通量测序数据的分析需求。 灵活的参数设置 自定义参数:MACS3 提供了丰富的参数选项,用户可以根据具体需求调整参数设置。 conda actiavte chipseq conda install bioconda::macs3 4功能简述 MACS3功能 5最小化使用 callpeak MACS3 中的主要函数,它用于处理各种格式的对齐文件
/aligned/sample1.bam & 使用MACS3 进行样本间比较(后续可视化使用) ###重复一起输入找共有peaks macs3 callpeak -t sample1.bam sample2 .bam -c input1.bam input2.bam -f BAMPE -n sample --outdir ./ -g 2.7e9 -B ###重复分别进行比对 macs3 callpeak -t sample1.bam -c input1.bam -f BAMPE -n sample1 --outdir ./ -g 2.7e9 -B macs3 callpeak -t sample2
质控 BWA 比对 MACS3 callpeak Step2: 获取组蛋白修饰 peaks 这部分的数据是:GSE100259 同样的,我们进行老三步:质控、BWA 比对、MACS3 callpeak。
calling at the cluster-levelsnap.tl.macs3(data, groupby='cell_type')peaks = snap.tl.merge_peaks(data.uns['macs3 callingsnap.tl.macs3(data, groupby='leiden', replicate='sample')merged_peaks = snap.tl.merge_peaks(data.uns['macs3
0.环境配置 0.1.安装IGV 从这里下载进行安装:https://www.broadinstitute.org/igv/ 0.2.安装MACS2 MACS2以及新版的MACS3信息可以在这里获取:https