PeakSnake/config.yaml

#######################
# globe  
#######################
PEAKSNAKE_HOME: /home/zhxd/workspace/20231225_ChIP_ATAC/PeakSnake
genome: ref/genome.fa
gtf: ref/genes_longest.gtf 
tss_tes_bed: ref/genes_longest.bed6  
tss_tes_shuffle_bed: ref/genes_shuffle.bed6
seq_type: PE  # 或者 "PE"
peak_type: broad # narrow | broad
# 有效基因组大小, 参考 https://deeptools.readthedocs.io/en/develop/content/feature/effectiveGenomeSize.html#effective-genome-size
gsize: 119481543 
# idr: 只支持 2 重复,适合 narrowPeak) ; intersect (2 + 重复 ) ; norep: 没有重复 
peak_selection: intersect

#######################
# samples  
#######################
samples:
  bZIP11_CHXDEX_rep1: 
  bZIP11_CHXDEX_rep2:
  bZIP11_CHXMOCK_rep1:
  bZIP11_CHXMOCK_rep2:
  WT_CHXDEX_rep1:
  WT_CHXDEX_rep2:

contrasts:
  - bZIP11_CHXDEX_vs_bZIP11_CHXMOCK
  - bZIP11_CHXDEX_vs_WT_CHXDEX

#######################
# fastp  
#######################
fastp: "--adapter_sequence AGATCGGAAGAGCACACGTCTGAACTCCAGTCA --adapter_sequence_r2 AGATCGGAAGAGCGTCGTGTAGGGAAAGAGTGT --qualified_quality_phred 20 --length_required 50"

#######################
# bowtie2 
# PE 建议设置 --no-mixed --no-discordant, 只保留成对比对
# ATAC-seq 必须设置 --maxins 1000 ~ 2000
#######################
bowtie2_index: ref/genome  # 不设置则从genome.fa 构建
bowtie2: --no-mixed --no-discordant --maxins 1000

#######################
# alignmentSieve 
# blackregion: 黑名单区域 线粒体 叶绿体 大肠杆菌等
# ATAC-seq 在 extra 中添加 --ATACshift 
#######################
alignmentSieve:
  minMappingQuality: 25
  blacklist: ref/blacklist.bed
  extra:

#######################
# bam to bw
#######################
bamCoverage: "--binSize 50 --normalizeUsing RPKM"

#######################
# MACS3
# fold: 用于过滤 peaks
# qscore: 用于过滤 peaks
# macs3 callpeak ChIP-seq 一般不需要设置特殊参数
# ATAC-seq 推荐: --nomodel --shift -100 --extsize 200 或 --nomodel --shift -75 --extsize 150
#######################
macs3: 
  qscore: 3
  extra: --nomodel --shift -100 --extsize 200 

#######################
# replicate intersect 
#######################
intersect:
  min_overlap: 0.5

#######################
# IDR 
#######################
idr:
  idr_threshold: 0.05

#######################
# 提取序列用于 motif 分析
#######################
motif:
  top_n: 1000
  summit_flank: 50

#######################
# Peank 注释
# 只关注 TSS 附近 用 ChIPseeker
# 否则UROPA 
#######################
uropa:
  do: yes
  feature: ['gene']
  feature_anchor: ['start', 'center', 'end']
  relative_location:  ['Upstream', 'Downstream', 'OverlapStart', 'OverlapEnd','PeakInsideFeature', 'FeatureInsidePeak']
  distance: [100000, 100000]

#######################
# diff analysis
#######################
diff_peaks:
  do: yes
  method: edgeR # edgeR | DESeq2
  padj: 0.05

#######################
# cross correlation
# ChIP-seq 需要进行, ATAC-seq 不需要 
#######################
cross_correlation:
  do: no
  spp: "-s=0:5:500"