Pipeline Overview

The ChIP-seq transcription factor pipeline was developed as a part of the ENCODE Uniform Processing Pipelines series. The ENCODE Consortium has developed two analysis pipelines to study the different classes of protein-chromatin interactions. Both ChIP-seq piplines share the same mapping steps, but differ in the methods for signal and peak calling and in the subsequent statistical treatment of replicates. The full ChIP-seq pipeline code is available on Github and can be run on DNAnexus (link requires account creation) at their current pricing.

Pipeline for Mapping of FASTQs

View the current mapping pipeline instance.

Pipeline for Transcription Factor ChIP-seq

Transcription factor ChIP-seq (TF ChIP-seq) specifically looks at proteins which are thought to associate with specific DNA sequences to influence the rate of transcription.

Replicated experiments

View the current replicated TF pipeline instance.

Unreplicated experiments

View the current unreplicated TF pipeline instance.

Back to the top

Inputs:

Inputs for mapping

File format	Information contained in file	File description	Notes
fastq	reads	G-zipped reads, paired-ended or single ended, stranded or unstranded	Multiple fastqs from a single biological replicate or library are concatenated before mapping. Reads must meet the criteria outlined under the Uniform Processing Pipeline Restrictions.
fasta	genome indices	Indices are dependent on the assembly being used for mapping

Inputs for peak calling

File format	Information contained in file	File description	Notes
bam	alignments (filtered)	Produced by mapping reads to the genome
bam	alignments	Filtered bam file from the control experiment	See the Uniform Processing Pipeline Restrictions for more on controls.

Outputs:

Transcription factor pipeline (replicated)

File format	Information contained in file	File description	Notes
bigWig	fold change over control, signal p-value	Two versions of nucleotide resolution signal coverage tracks.	The signal is expressed in two ways: as fold-over control at each position, and as a p-value to reject the null hypothesis that the signal at that location is present in the control.
bed and bigBed (narrowPeak)	peaks	Relaxed peak calls for each replicate individually and for both replicates' reads pooled together.	These peaks are thresholded to sample enough noise in the experiment for efficient statistical comparison of replicates in subsequent steps; as such, many false positives are expected to be present. They are not meant to be interpreted as definitive binding events, but are rather intended to be used as input for subsequent statistical comparison of replicates.
bed and bigBed (narrowPeak)	conservative IDR peaks	Peaks derived from IDR analysis of biological replicates.	IDR, or Irreproducible Discovery Rate, measures the reproducibility of the assay in question. The optimal set of peaks is more sensitive, especially when one replicate has substantially lower data quality than the other.
bed and bigBed (narrowPeak)	optimal IDR peaks	The largest set of peaks derived from IDR analysis of biological replicates and pseudoreplicates.
Quality control metrics are collected to determine library complexity, read depth, FRiP score, and reproducibility.

Transcription factor pipeline (unreplicated)

File format	Information contained in file	File description	Notes
bigWig	fold change over control, signal p-value	Two versions of nucleotide resolution signal coverage tracks.	The signal is expressed in two ways: as fold-over control at each position, and as a p-value to reject the null hypothesis that the signal at that location is present in the control.
bed and bigBed (narrowPeak)	peaks	Relaxed peak calls from all reads in the single replicate	These peaks are thresholded to sample enough noise in the experiment for efficient statistical comparison of pseudoreplicates in subsequent steps; as such, many false positives are expected to be present. They are not meant to be interpreted as definitive binding events, but are rather intended to be used as input for subsequent statistical comparison of replicates.
bed and bigBed (narrowPeak)	pseudoreplicated IDR thresholded peaks	The set of peak calls from two partitions, or "pseudoreplicates"	In the "unreplicated IDR" step, an IDR-based strategy is used to identify stable peaks across pseudoreplicates (instead of across true replicates). Peaks from the relaxed set (from all reads) as well as from pseudoreplicates are ranked, and those peaks that are well-supported in both pseudoreplicates (i.e. cross the same IDR threshold as for replicated experiments) are reported. There are no corresponding optimal and conservative peak sets in the unreplicated case because there are no true replicates.
Quality control metrics are collected to determine library complexity, read depth, FRiP score, and reproducibility.

References

Genomic References

View the mapping assembly and genome annotation reference files used in this pipeline

Links and Publications

Find data generated by the transcription factor pipeline for replicated experiments
Find data generated by the transcription factor pipeline for unreplicated experiments
Explore all ChIP-seq related publications on the ENCODE portal

Back to the top

Pipeline schematic