ChromHMM Chromatin State Inference

Overview

This skill enables comprehensive chromatin state analysis using chromHMM for histone modification ChIP-seq data. ChromHMM uses a multivariate Hidden Markov Model to segment the genome into discrete chromatin states based on combinatorial patterns of histone modifications.

Main steps include:

• Refer to Inputs & Outputs to verify necessary files.
• Always prompt user if required files are missing.
• Always prompt user for genome assembly used.
• Always prompt user for the bin size for generating binarized files.
• Always prompt user for the bin size for the number of states the ChromHMM target.
• Run chromHMM workflow: Binarization → Learning.

---

When to use this skill

Use this skill when you need to infer chromatin states from histone modification ChIP-seq data using chromHMM.

---

Inputs & Outputs

Inputs

(1) Option 1: BED files of aligned reads

<mark1>.bed
<mark2>.bed
... # Other marks

(1) Option 2: BAM files of aligned reads

<mark1>.bam
<mark2>.bam
... # Other marks

Outputs

chromhmm_output/
  binarized/
    *.txt 
  model/
    *.txt
    ... # other files output by the ChromHMM

---

Decision Tree

Step 0: Initialize Project

Call:

• mcp__project-init-tools__project_init

with:

• sample

  : all

• task

  : chromhmm

Step 1: Prepare the

cellmarkfile

(skip this step if signal files are provided)

• Prepare a .txt file (without header) containing following three columns:

  • sample name
  • marker name
  • name of the BED/BAM file
  • control file of the sample (only provided if the input/control file is available)

• example of the cellmark.txt file

cell1    mark1    cell1_mark2.bam    cell1_control.bam
cell1   mark2    cell1_mark2.bam    cell1/control.bam

Step 2: Data Binarization

• For BAM inputs:
  Call:

  • mcp__chromhmm-tools__binarize_bam

    with:

  • path_chrom_sized

    : Provide by user or detect from the working directory

  • input_dir

    : Directory containing BAM files

  • cellmarkfile

    : Cell mark file defining histone modifications

  • output_dir

    : (e.g.

    binarized/

    )

  • bin_size

    : Provided by user

• For BED inputs:
  Call

  mcp__chromhmm-tools__binarize_bed

  instead.

• For Signal inputs:
  Call:

  mcp__chromhmm-tools__binarize_signal

  with:

  • input_dir

    : Directory of signals

  • output_dir

    : (e.g.

    binarized/

    )

Step 3: Model Learning

Call

• mcp__chromhmm-tools__learn_model

with:

• binarized_dir

  : Directory binarized file located in

• num_states

  : Provide by user (e.g. 15)

• output_model_dir

  : (e.g.

  model_15_states/

  )

• genome

  : Provide by user (e.g.

  hg38

  )

• threads

  : Provide by user (e.g. 16)

Parameter Optimization

Number of States

• 8 states: Basic chromatin states
• 15 states: Standard comprehensive states
• 25 states: High-resolution states
• Optimization: Use Bayesian Information Criterion (BIC)

Bin Size

• 200bp: Standard resolution
• 100bp: High resolution (requires more memory)
• 500bp: Low resolution (faster computation)

State Interpretation

Common Chromatin States

1. Active Promoter: H3K4me3, H3K27ac
2. Weak Promoter: H3K4me3
3. Poised Promoter: H3K4me3, H3K27me3
4. Strong Enhancer: H3K27ac, H3K4me1
5. Weak Enhancer: H3K4me1
6. Insulator: CTCF
7. Transcribed: H3K36me3
8. Repressed: H3K27me3
9. Heterochromatin: Low signal across marks

Troubleshooting

• Memory errors: Reduce bin size or number of states
• Convergence problems: Increase iterations or adjust learning rate
• Uninterpretable states: Check input data quality and mark combinations
• Missing chromosomes: Verify chromosome naming consistency