deeptools

deepTools - An NGS Data Analysis and Visualization Toolkit

Skill Overview

deepTools is a collection of Python command-line tools specifically designed for high-throughput sequencing data analysis. It helps researchers process, perform quality control on, and visualize data from experiments such as ChIP-seq, RNA-seq, and ATAC-seq. With deepTools, BAM alignment files can be converted into normalized coverage tracks (bigWig), sample correlation analyses and signal enrichment assessments can be performed, and publication-quality heatmaps and signal profile plots can be generated.

Applicable Scenarios

1. Complete ChIP-seq Data Analysis Workflow

A complete workflow from raw BAM files to publication-quality figures, including quality control, normalized coverage calculation, treatment-versus-control comparisons, and signal visualization around TSSs and peak regions. The skill provides prebuilt workflow templates and supports one-click generation of ChIP-seq QC scripts for rapidly assessing replicate correlation, enrichment strength, and sequencing depth.

2. Strand-Specific Coverage Analysis for RNA-seq

Provides strand-specific coverage track generation tailored to the characteristics of RNA-seq data. Supports separate visualization of the positive and negative strands, RPKM/CPM normalization, and proper handling of splice regions without incorrectly extending reads. It is particularly suitable for analyzing data generated using common strand-specific library construction methods such as the dUTP method.

3. Tn5 Offset Correction and Signal Visualization for ATAC-seq

When processing ATAC-seq data, automatically performs Tn5 cleavage-site offset correction (--ATACshift) to generate accurate chromatin accessibility signal profiles. Supports nucleosome positioning pattern analysis (nucleosome ladder plots), fragment length distribution assessment, and signal enrichment analysis in peak regions.

Core Functions

1. BAM/bigWig File Processing and Normalization

Provides tools such as bamCoverage and bamCompare to convert BAM alignment files into bigWig coverage tracks. It includes multiple normalization methods (RPGC, CPM, RPKM, and BPM), automatically handles genome size parameters, and supports common operations such as read extension, duplicate filtering, and mapping-quality filtering to ensure signal comparability between samples.

2. Quality Control and Sample Assessment

Uses tools such as plotFingerprint, plotCorrelation, plotPCA, and plotCoverage to comprehensively assess sequencing data quality. These tools can detect ChIP enrichment strength, evaluate replicate correlation, determine whether sequencing depth is sufficient, and identify abnormal samples. They provide intuitive visual outputs and guidance for interpreting the results.

3. Genome Signal Visualization

Uses computeMatrix to calculate signal matrices over genomic feature regions, and plotHeatmap and plotProfile to generate publication-quality heatmaps and signal profile plots. It supports reference-point mode (such as TSSs) and scale-regions mode (such as gene bodies), along with extensive visualization options including clustering, color mapping, and custom legends.

Frequently Asked Questions

What types of sequencing data analysis does deepTools support?

deepTools primarily supports DNA sequencing-based analyses, including ChIP-seq, ATAC-seq, MNase-seq, and DNA-seq, as well as coverage visualization for RNA-seq. For ChIP-seq, it provides complete QC, normalization, and visualization workflows; for ATAC-seq, it includes Tn5 offset correction; and for RNA-seq, it supports strand-specific coverage analysis. The skill is not intended for variant calling, transcript quantification, or other tasks that require different tools.

How should I choose an appropriate normalization method?

The choice of normalization method depends on the experiment type and analysis objective. For ChIP-seq coverage analysis, RPGC (which requires the effective genome size) or CPM is recommended. For comparisons between ChIP-seq samples, use bamCompare together with a log2 transformation. For RNA-seq fixed-region analysis, use CPM; for gene-level analysis, use RPKM to account for gene length. For ATAC-seq, use RPGC or CPM. The skill includes a table of effective genome sizes for common species, along with detailed guidance on selecting normalization methods.

What special requirements apply to ATAC-seq analysis?

ATAC-seq data processing requires particular attention to the offset caused by Tn5 cleavage. The --ATACshift parameter of alignmentSieve should be used for correction (equivalent to --shift 4 -5 5 -4, while retaining only properly paired fragments). Using --extendReads for ATAC-seq data is not recommended. During analysis, the fragment length distribution should be examined; ideally, it should display a nucleosome ladder pattern with approximately 200 bp periodicity. The skill provides dedicated ATAC-seq workflow templates covering offset correction, coverage calculation, and visualization steps.