Discover and Characterize Regulatory Elements with ENCODE
When to Use
- User wants to find enhancers, promoters, silencers, or insulators in a specific tissue using ENCODE
- User asks about "regulatory elements", "cCREs", "enhancer discovery", "ChromHMM", or "super-enhancers"
- User needs to classify chromatin states or identify active regulatory regions from histone mark data
- User wants to use ENCODE's 926,535 cCRE catalog or run functional validation (CRISPR/MPRA/reporter)
- Example queries: "find active enhancers in liver", "classify chromatin states for my tissue", "identify super-enhancers from H3K27ac data"
Identify, classify, and functionally characterize regulatory elements using ENCODE's catalog of 926,535 human candidate cis-regulatory elements (cCREs) and layered functional genomics data.
Scientific Rationale
The question: "What regulatory elements are active in my tissue of interest, and what are they doing?"
The human genome contains an estimated 1–2 million regulatory elements — far outnumbering the ~20,000 protein-coding genes. These elements (enhancers, promoters, silencers, insulators) control when, where, and how much each gene is expressed. No single biochemical assay can definitively identify a regulatory element; instead, combinatorial patterns of chromatin marks, accessibility, and TF binding are used to classify candidate elements.
The ENCODE cCRE Registry
The ENCODE Phase 3 project (ENCODE Project Consortium 2020) established a registry of 926,535 human and 339,815 mouse cCREs covering 7.9% and 3.4% of their respective genomes. These are classified using combinations of DNase-seq, H3K4me3, H3K27ac, and CTCF ChIP-seq signals across hundreds of biosamples. The registry is accessible via the SCREEN web server and represents the most comprehensive catalog of candidate regulatory elements in any organism.
An expanded registry (Moore et al. 2024, bioRxiv preprint) extends this to 2.35 million human cCREs with functional characterization from STARR-seq, MPRA, and CRISPR perturbation covering >90% of human cCREs.
Key Distinction: Candidate vs. Validated
ENCODE cCREs are candidate regulatory elements identified by biochemical signatures. Biochemical activity (histone marks, accessibility) is necessary but not sufficient for function. A region marked by H3K27ac is likely regulatory, but functional validation (perturbation, reporter assays) is required to confirm that it actually regulates a target gene. The gap between biochemical annotation and validated function is the central challenge.
Literature Support
- ENCODE Project Consortium 2020 (Nature, ~1,656 citations): Registry of 926,535 human cCREs. Introduces the SCREEN web server. DOI
- Kundaje et al. 2015 (Nature, ~5,810 citations): Roadmap Epigenomics — integrative analysis of 111 reference human epigenomes. Chromatin state maps across tissues. Disease variants enriched in tissue-specific epigenomic marks. DOI
- Ernst & Kellis 2012 (Nature Methods, ~2,294 citations): ChromHMM — multivariate hidden Markov model for chromatin state discovery from combinatorial histone modification patterns. DOI
- Hnisz et al. 2013 (Cell, ~3,215 citations): Super-enhancer catalog across human cell types. Disease-associated variation enriched in super-enhancers of disease-relevant cells. DOI
- Whyte et al. 2013 (Cell, ~2,500 citations): Defined super-enhancers as large enhancer clusters occupied by master TFs and Mediator. Introduced the ROSE algorithm. DOI
- Shlyueva et al. 2014 (Nature Reviews Genetics, ~1,200 citations): Authoritative review of enhancer sequence properties, chromatin signatures, genome-wide prediction, and high-throughput activity assays. DOI
- Schoenfelder & Fraser 2019 (Nature Reviews Genetics, ~869 citations): How enhancer-promoter interactions are established through 3D genome architecture (TADs, CTCF loops). DOI
- Visel et al. 2007 (Nucleic Acids Research, ~1,079 citations): VISTA Enhancer Browser — in vivo transgenic mouse validation of enhancers. 4,500+ experiments. DOI
- Gasperini et al. 2019 (Cell, ~465 citations): CRISPRi screen of 5,920 candidate enhancers with scRNA-seq readout. Identified 664 enhancer-gene pairs. Established the "crisprQTL" framework. DOI
- Yao et al. 2024 (Nature Methods, ~26 citations): ENCODE4 Functional Characterization Centers — 108 CRISPRi screens, >540,000 perturbations. Pre-designed sgRNAs targeting 3.27M ENCODE SCREEN cCREs. DOI
- Nasser et al. 2021 (Nature, ~468 citations): ABC model enhancer-gene maps in 131 cell types. DOI
- Heintzman et al. 2007 (Nature Genetics, ~2,300 citations): Discovered that H3K4me1 marks enhancers while H3K4me3 marks promoters — the foundational chromatin signature for distinguishing regulatory element classes. DOI
- Rada-Iglesias et al. 2011 (Nature, ~1,200 citations): Identified "poised enhancers" marked by H3K4me1+H3K27me3 (without H3K27ac) in hESCs. These activate during differentiation by gaining H3K27ac. DOI
- Amemiya et al. 2019 (Scientific Reports, ~1,372 citations): ENCODE Blacklist — regions producing artifactual signal across ChIP-seq, ATAC-seq, and DNase-seq. Must be filtered before any regulatory element classification. DOI
Step 1: Define the Element Type and Tissue Context
ENCODE cCRE Classification System
| cCRE Class | Abbreviation | Biochemical Signature | Genomic Context | Example |
|---|---|---|---|---|
| Promoter-like | PLS | DNase+ H3K4me3+ (±H3K27ac) | Within 200bp of annotated TSS | Gene promoter |
| Proximal enhancer-like | pELS | DNase+ H3K27ac+ (H3K4me3-) | Within 2kb of TSS | Proximal enhancer |
| Distal enhancer-like | dELS | DNase+ H3K27ac+ (H3K4me3-) | >2kb from TSS | Distal enhancer |
| CTCF-only | CTCF-only | DNase+ CTCF+ (no H3K4me3/H3K27ac) | Any | Insulator/boundary |
| DNase-H3K4me3 | DNase-H3K4me3 | DNase+ H3K4me3+ | >200bp from TSS | Unannotated promoter-like |
Extended Element Types (Beyond cCRE Classification)
| Element Type | Key Signatures | ENCODE Assays | Notes |
|---|---|---|---|
| Active promoter | H3K4me3+ H3K27ac+ accessible | Histone ChIP-seq, ATAC/DNase | Corresponds to PLS cCREs |
| Active enhancer | H3K4me1+ H3K27ac+ H3K4me3- accessible | Histone ChIP-seq, ATAC/DNase | Corresponds to pELS/dELS |
| Poised enhancer | H3K4me1+ H3K27me3+ H3K27ac- | Histone ChIP-seq | Bivalent; may activate upon differentiation |
| Primed enhancer | H3K4me1+ only (no H3K27ac, no H3K27me3) | Histone ChIP-seq | Ready for activation but not currently active |
| Super-enhancer | Broad H3K27ac, multiple TFs, high Mediator | ChIP-seq, TF ChIP-seq | ROSE algorithm (Whyte 2013) |
| Silencer | H3K27me3+ (Polycomb) or H3K9me3+ (heterochromatin) | Histone ChIP-seq | Two distinct repressive mechanisms |
| Insulator | CTCF binding at TAD boundary | TF ChIP-seq (CTCF), Hi-C | Blocks enhancer-promoter communication |
| Stretch enhancer | >3kb H3K27ac domain (not classified as super-enhancer) | Histone ChIP-seq | Parker et al. 2013; enriched for disease variants |
Check data availability for the target tissue:
encode_get_facets(organ="...", biosample_type="tissue")
Step 2: Search for Tissue-Specific Functional Data
Minimal Data Requirements
For basic regulatory element identification, you need