encode-toolkit

List, search, and inspect ENCODE files by format, type, and assembly

Generate ENCODE citations for publications, grants, and presentations

Check if two ENCODE experiments are compatible for combined analysis

Cross-reference ENCODE data with PubMed, GEO, ClinicalTrials, and bioRxiv

Download ENCODE files (BED, FASTQ, BAM, bigWig) with MD5 verification

Execute ENCODE ATAC-seq pipeline from FASTQ to accessibility peaks with Tn5 correction, Bowtie2, and MACS2

Execute ENCODE ChIP-seq pipeline from FASTQ to peaks and signal tracks using BWA-MEM, MACS2, and IDR

Execute CUT&RUN pipeline from FASTQ to peaks with Bowtie2, SEACR, and spike-in normalization

Execute ENCODE DNase-seq pipeline from FASTQ to hotspots and footprints using BWA, Hotspot2, and HINT-ATAC

Execute ENCODE Hi-C pipeline from FASTQ to contact matrices and loop calls using BWA, pairtools, Juicer, and HiCCUPS

Guide for using JASPAR transcription factor binding profiles with ENCODE ChIP-seq data. Use when users need to find TF binding motifs in ENCODE peaks, validate ChIP-seq targets with known motifs, or scan regulatory regions for TF binding potential. Trigger on: JASPAR, motif database, binding profile, PWM, position weight matrix, TF motif, motif enrichment, motif scanning, binding site prediction.

Convert genomic coordinates between assembly versions (GRCh37/hg19 to GRCh38/hg38, mm9 to mm10). Guides UCSC liftOver for BED files, CrossMap for VCF/bigWig, and handles unmapped regions with provenance logging.

Build comprehensive DNA methylation maps by aggregating WGBS (Whole Genome Bisulfite Sequencing) data across multiple ENCODE experiments, donors, and labs. Use when the user wants to answer "where is DNA methylated/unmethylated in my tissue?" by combining per-CpG methylation data into tissue-level methylation profiles. Handles coverage filtering, identifies hypomethylated regions (HMRs) and partially methylated domains (PMDs), and manages cross-lab variation.

Guide for de novo and known motif enrichment analysis of ENCODE ChIP-seq and ATAC-seq peaks using HOMER and MEME Suite. Use when users need to discover TF binding motifs in peaks, validate ChIP-seq targets, or find co-binding partners. Trigger on: motif analysis, HOMER, MEME, de novo motif, motif enrichment, findMotifsGenome, AME, MEME-ChIP, known motif, TF binding motif, co-factor, motif discovery.

Integrate multiple ENCODE data types (RNA-seq, ATAC-seq, Histone ChIP-seq, TF ChIP-seq) for a tissue/cell type to build a comprehensive regulatory landscape. Use when the user wants to answer "what are the enhancers, promoters, and regulatory elements active in my tissue, and which transcription factors control them?" by layering expression, chromatin accessibility, histone marks, and TF binding data. Follows the Mawla et al. 2023 framework for cross-assay integration of islet cell type-specific data. Handles chromatin state annotation (ChromHMM), enhancer-gene linkage, TF motif enrichment, and cell type-specific regulatory element identification. Use for ANY multi-omic analysis, enhancer discovery, regulatory network construction, or epigenomic characterization using ENCODE data.