ReMap : Atlas of Regulatory Regions

Description ReMap 2020 Homo Sapiens

This track hub represents the ReMap Atlas of regulatory regions which consists in a large scale integrative analysis of all Public ChIP-seq data for transcriptional regulators from GEO, ArrayExpress and ENCODE.
To enable genome-wide identification of regulatory elements we have collected, curated, analysed a total of 6,498 ChIP-seq and Exo-seq data sets from Public source (GEO, ArrayExpress and ENCODE). After applying our quality filters we retained 6,000 data sets.
Those merged analyses covers a total of 960 DNA-binding protein (transcriptional regulators) such as a variety of transcription factors (TFs), transcription co-activators (TCFs) and chromatin-remodeling factors (CRFs) for 164 million peaks.

Below a schematic diagram of the types of regulatory regions:

Schematic diagram

Individual BED files for specific TFs, or Cells or datasets can be found and downloaded on the ReMap website or

Display Conventions and Configuration


Data sets sources

GEO & ArrayExpress

Public ChIP-seq data sets were extracted from Gene Expression Omnibus (GEO) and ArrayExpress (AE) databases. For GEO, the query '('chip seq' OR 'chipseq' OR 'chip sequencing') AND 'Genome binding/occupancy profiling by high throughput sequencing' AND 'homo sapiens'[organism] AND NOT 'ENCODE'[project]' was used to return a list of all potential data sets to analyse, which were then manually assessed for further analyses. Data sets involving polymerases (i.e. Pol2 and Pol3), and some mutated or fused TFs (e.g. KAP1 N/C terminal mutation, GSE27929) were exckuded.


Available ENCODE ChIP-seq data sets for transcriptional regulators from portal were processed with the our uniform workflow. We retrieved the list of ENCODE data as FASTQ files from the ENCODE portal ( using the following filters: Assay: "ChIP-seq", Organism: "Homo sapiens", Target of assay: "transcription factor", Available data: "fastq" on 2016 June 21st. Metadata information in JSON format and FASTQ files were retrieved using the Python requests module.

ChIP-seq processing

Both Public and ENCODE data were processed similarly. Bowtie 2 (PMC3322381) (version 2.2.9) with options -end-to-end -sensitive was used to align all reads on the human genome (GRCh38/hg38 assembly). Biological and technical replicates for each unique combination of GSE/TF/Cell type or Biological condition were used for peak calling. TFBS were identified using MACS2 peak-calling tool (PMC3120977) (version in order to follow ENCODE ChIP-seq guidelines, with stringent thresholds (MACS2 default thresholds, p-value: 1e-5). An input data set was used when available.

Quality assessment

To assess the quality of public data sets, we computed a score based on the cross-correlation and the FRiP (fraction of reads in peaks) metrics developed by the ENCODE Consortium ( Two thresholds were defined for each of the two cross-correlation ratios (NSC, normalized strand coefficient: 1.05 and 1.10; RSC, relative strand coefficient: 0.8 and 1.0). Detailed descriptions of the ENCODE quality coefficients can be found at We used the phantompeak tools suite ( to compute RSC and NSC.

Full details of our methods can be found in the references below.

Data Availability

The ReMap BED files are available for download at the ReMap website or in the download tab.

Papers to cite

  1. Cheneby J., Gheorghe M., Artufel M., Mathelier A., Ballester, B.
    ReMap 2018: An updated regulatory regions atlas from an integrative analysis of DNA-binding ChIP-seq experiments.
    Nucleic Acids Research (2018) gkx1092
  2. Griffon A., Barbier Q., Dalino J., van Helden J., Spicuglia S., Ballester B.
    Integrative analysis of public ChIP-seq experiments reveals a complex multi-cell regulatory landscape.
    Nucleic Acids Research (2015) 43 (4): e27.


If you have questions or comments, please write to: