accessibility-aggregat…
Build comprehensive chromatin accessibility maps by aggregating ATAC-seq and DNase-seq narrowPeak data across multiple ENCODE…
Compare ENCODE experiments across different biosamples, tissues, or cell lines to identify tissue-specific regulatory patterns. Use when the user wants cross-tissue comparison, cell-type comparison, tissue-specific elements, differential chromatin, biosample matching, disease vs
Claude Code
$ npx -y skills add ammawla/encode-toolkit --skill compare-biosamples --agent claude-code
Codex
$ npx skills add ammawla/encode-toolkit --skill compare-biosamples -a codex -y
Cursor
$ npx skills add ammawla/encode-toolkit --skill compare-biosamples -a cursor -y
Gemini CLI
$ npx skills add ammawla/encode-toolkit --skill compare-biosamples -a gemini-cli -y
GitHub Copilot
$ npx skills add ammawla/encode-toolkit --skill compare-biosamples -a github-copilot -y
OpenCode
$ npx skills add ammawla/encode-toolkit --skill compare-biosamples -a opencode -y
Windsurf
$ npx skills add ammawla/encode-toolkit --skill compare-biosamples -a windsurf -y
name: compare-biosamples description: Compare ENCODE experiments across different biosamples, tissues, or cell lines to identify tissue-specific regulatory patterns. Use when the user wants cross-tissue comparison, cell-type comparison, tissue-specific elements, differential chromatin, biosample matching, disease vs normal comparison, developmental time course, constitutive vs variable regulation, or multi-tissue data availability mapping. Handles batch effect detection, biosample hierarchy, and comparison design.
Help the user systematically compare data availability and experiments across different biosamples to identify tissue-specific regulatory patterns, constitutive elements, and cross-tissue differences.
Cross-biosample comparison is the foundation of understanding tissue-specific gene regulation. Regulatory elements -- particularly enhancers -- are the primary drivers of cell-type identity, with promoters being largely shared across tissues. Comparing the same assay across multiple biosamples reveals which regulatory elements are constitutive (shared) versus tissue-specific (unique to one or few cell types).
**The core question**: "Which regulatory features distinguish tissue A from tissue B, and which are shared?"
This requires careful matching of datasets, awareness of batch effects, and understanding of the biosample hierarchy to avoid confounding biological signal with technical variation.
| # | Reference | Key Contribution | |---|-----------|-----------------| | 1 | Roadmap Epigenomics Consortium 2015, Nature, DOI:10.1038/nature14248 (~5,810 cit) | Generated 111 reference epigenomes across tissues/cell types; established the framework for cross-tissue epigenomic comparison. Showed that enhancer chromatin states are the most tissue-variable elements. | | 2 | ENCODE Phase 3 2020, Nature, DOI:10.1038/s41586-020-2493-4 (~1,656 cit) | Expanded functional annotations to 1.3M candidate cis-regulatory elements (cCREs) across hundreds of biosamples; defined tissue-activity indices for regulatory elements. | | 3 | Andersson et al. 2014, Nature, DOI:10.1038/nature12787 (~1,500 cit) | FANTOM5 atlas of active enhancers across 808 samples; demonstrated that only ~5% of enhancers are active across all tissues, with the majority being highly tissue-specific. | | 4 | Heintzman et al. 2009, Nature, DOI:10.1038/nature07917 (~2,200 cit) | Showed histone modifications distinguish cell types: H3K4me1/H3K27ac at enhancers are the most discriminating tissue-specific marks, while H3K4me3 at promoters is largely shared. | | 5 | Thurman et al. 2012, Nature, DOI:10.1038/nature11232 (~2,000 cit) | Mapped accessible chromatin across 125 cell types; demonstrated that DNase I hypersensitive sites define cell-type identity and that accessibility patterns cluster by tissue of origin. | | 6 | Leek et al. 2010, Nat Rev Genet, DOI:10.1038/nrg2825 (~1,200 cit) | Comprehensive review of batch effects in genomic data; showed that lab, platform, and processing date can dominate biological variation if not properly controlled. | | 7 | Forrest et al. 2014, Nature, DOI:10.1038/nature13182 (~1,100 cit) | FANTOM5 promoter-level expression atlas across 975 samples; demonstrated that promoter usage (not just gene expression) is tissue-specific and defines cell identity. |
Understanding what varies across tissues and what does not is essential before designing a comparison.
| Feature | Cross-Tissue Behavior | Implication for Comparison | |---------|----------------------|---------------------------| | **Promoters (H3K4me3)** | Largely shared (~70% active in most tissues) | Poor discriminators between tissues | | **Enhancers (H3K27ac + H3K4me1)** | Highly tissue-specific (~5% shared across all tissues) | Best discriminators; focus comparison here | | **Chromatin accessibility (ATAC/DNase)** | Moderate tissue-specificity (~20-30% shared) | Good secondary discriminator; clusters by tissue of origin | | **Polycomb repression (H3K27me3)** | Tissue-specific (marks silenced developmental genes) | Useful for identifying repressed lineage programs | | **Gene expression (RNA-seq)** | Moderate tissue-specificity | Housekeeping genes shared; tissue-specific TFs are key | | **CTCF binding** | Largely constitutive (~70% conserved) | Defines structural boundaries; less tissue-variable | | **DNA methylation** | Bimodal; enhancers show tissue-variable methylation | Hypomethylation at active enhancers is tissue-specific |
H3K27ac at enhancers is the single most informative mark for distinguishing tissues (Heintzman et al. 2009, Roadmap 2015). If the user can only compare one mark across tissues, H3K27ac should be the first choice, followed by chromatin accessibility (ATAC-seq or DNase-seq).
| Level | Description | Biological Relevance | Reproducibility | Caveats | |-------|-------------|---------------------|-----------------|---------| | **Tissue** | Primary tissue from donor (e.g., pancreas, liver) | Highest -- in vivo biology preserved | Lower -- donor variation, cell-type heterogeneity | Mixed cell populations; composition varies by donor age/sex/health | | **Primary cell** | Cells isolated from tissue (e.g., hepatocytes, islets) | High --
Search ENCODE, cross-reference 14 databases, run 7 analysis pipelines, and generate publication-ready methods — all from natural language in Claude Code.
Build comprehensive chromatin accessibility maps by aggregating ATAC-seq and DNase-seq narrowPeak data across multiple ENCODE…
Guide for multi-experiment batch operations: QC screening, batch download, comparison, and report generation across many ENCODE…
Install bioinformatics tools for ENCODE data analysis. Covers CLI tools (BWA, STAR, samtools, MACS2), R/Bioconductor packages…
Guide for integrating CellxGene Census single-cell data with ENCODE bulk experiments. Use when users need cell-type-specific…
Generate proper ENCODE citations for publications, grants, and presentations. Use when the user needs to cite ENCODE data, create…
Guide for annotating ENCODE regulatory variants with ClinVar clinical significance. Use when users need to check if variants in…
Use this skill to extract structured Markdown/JSON from PDFs and document images—tables with cell-level precision, formulas as…
Estimates intracellular metabolic fluxes from steady-state carbon-13 isotope-tracing measurements using validated atom maps,…