用途:从MeRIP-seq(甲基化RNA免疫沉淀测序)数据中,通过IP样本与input对照的比较,调用m6A修饰峰值,识别m6A修饰位点。
主要能力与工具:
- exomePeak2(推荐):支持生物学重复的峰值调用。用法为调用exomePeak2()函数,传入bam_ip、bam_input、gff、genome、paired_end等参数,随后用exportResults(result, format = 'BED')导出峰值。
- MACS3替代方案:将input视为对照进行峰值调用,命令为macs3 callpeak -t IP_rep1.bam IP_rep2.bam -c Input_rep1.bam Input_rep2.bam -f BAMPE -g hs -n m6a_peaks --nomodel --extsize 150 -q 0.05。
- MeTPeak:支持GTF注释的峰值调用,调用metpeak()函数,传入IP_BAM、INPUT_BAM、GENE_ANNO_GTF、OUTPUT_DIR参数。
峰值过滤:可按倍数富集与q值过滤,典型阈值为FC > 2、q < 0.05,示例命令为awk '$7 > 2 && $9 < 0.05' peaks.xls > filtered_peaks.bed。
相关技能:merip-preprocessing用于峰值调用前的数据准备;m6a-differential用于比较不同条件下的峰值;chip-seq/peak-calling概念类似。
限制:资料未提供各工具的具体版本要求、运行环境依赖、输入文件格式细节及性能说明,暂未提供。