Comprehensive ChIP-seq QC Pipeline

Overview

This skill performs a full ChIP-seq quality control analysis from aligned BAM files and peak files.

Main steps include:

• Refer to the Inputs & Outputs section to check inputs and build the output architecture. All the output file should located in ${proj_dir} in Step 0.
• Perform cross-correlation analysis to calculate NSC and RSC.
• Compute FRiP (Fraction of Reads in Peaks) using peak files and aligned BAMs.

Inputs & Outputs

Inputs

${sample}.bam # filtered bam files
${sample}.narrowPeak # or broadPeak

Outputs

all_chip_qc/
    ${sample}_spp.txt
    ${sample}_crosscorr.pdf
    ${sample}_frip.txt

Step 0: Initialize Project

Call:

• mcpproject-init-toolsproject_init

with:

• sample: all
• task: atac_qc

The tool will:

• Createallchipqc directory.
• Return the full path of the allchipqc directory, which will be used as ${proj_dir}.

Step 1: Calculate Cross-Correlation Metrics (NSC, RSC)

Call:

• mcpqc-toolsrun_phantompeakqualtools

with:

• bam_file: Path to BAM file
• outputdir: ${projdir}/

Output: ${sample}spp.txt, ${sample}crosscorr.pdf

Step 2: Calculate the fraction of reads falling within peak regions.

Call:

• mcpqc-toolscalculate_frip

with: bamfile: Path to BAM file. peakfile: Path to Peak file (BED/narrowPeak/broadPeak). outputdir: ${projdir}/

Output: ${sample}_frip.txt