SuperAPI AI 工具集

bio-epitranscriptomics-m6a-peak-calling

数据与分析

查看项目

更新于 10-15 21:55

技能介绍

用途:从MeRIP-seq(甲基化RNA免疫沉淀测序)数据中,通过IP样本与input对照的比较,调用m6A修饰峰值,识别m6A修饰位点。

主要能力与工具:

  1. exomePeak2(推荐):支持生物学重复的峰值调用。用法为调用exomePeak2()函数,传入bam_ip、bam_input、gff、genome、paired_end等参数,随后用exportResults(result, format = 'BED')导出峰值。
  2. 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。
  3. 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概念类似。

限制:资料未提供各工具的具体版本要求、运行环境依赖、输入文件格式细节及性能说明,暂未提供。