Transcription Factor Motif Enrichment & Regulatory Factor Discovery (HOMER)
Transcription Factor Motif Enrichment & Regulatory Factor Discovery (HOMER)
Uncovering transcription factor binding motifs within differential H3K4me3 peak summits. Running HOMER findMotifsGenome.pl with GC-content matched background, identifying known and de novo motifs (SP1, CREB, REST/NRSF), and mapping transcriptional repressor occupancy in Huntington's Disease.
1. Regulatory Factor & Motif Discovery with HOMER
Histone modification changes are often driven by upstream transcription factor binding alterations or recruitment of chromatin remodelers. By scanning peak summits for over-represented sequence motifs, we identify the exact transcription factor families whose binding sites are disrupted in Huntington's Disease.
2. Running HOMER findMotifsGenome.pl
# Convert differential peaks to HOMER BED format (chr, start, end, name, score, strand)
awk '{print $1"\t"$2"\t"$3"\tPeak_"NR"\t"$5"\t+"}' \
05_diffbind/hd_lost_peaks.bed > 05_diffbind/hd_lost_peaks_homer.bed
# Execute HOMER motif discovery (+/- 200 bp around peak center)
findMotifsGenome.pl \
05_diffbind/hd_lost_peaks_homer.bed \
hg38 \
05_diffbind/homer_motifs/ \
-size 200 \
-mask \
-p 83. Key Transcription Factor Hits in Huntington's Disease
HOMER compares peak sequences against randomly generated background genomic intervals matched for GC content. In our Huntington's Disease analysis, top enriched motifs include:
- SP1 (Specificity Protein 1): Essential GC-box transcription factor known to interact directly with mutant huntingtin; sequestration by $mHTT$ leads to promoter inactivation.
- CREB / ATF: cAMP response element-binding protein, critical for neuroprotection and BDNF transcription.
- REST/NRSF: RE1-Silencing Transcription factor; in healthy cells, wild-type huntingtin sequesters REST in the cytoplasm. In HD, mutant huntingtin allows REST to translocate into the nucleus and repress neuronal genes.
Common Pitfalls & Troubleshooting Checklist
- Java Heap Exhaustion (GATK / Picard): If processes terminate with
java.lang.OutOfMemoryError, set your environment variable explicitly:export _JAVA_OPTIONS="-Xmx16g -Djava.io.tmpdir=/tmp". - Index Timestamp Mismatch: If tools complain that a BAM index is missing or out of date, regenerate with
samtools index -@ 4 sample.bam. - Conda Channel Conflicts: Avoid ABI incompatibilities between bioconda and conda-forge by setting strict priority:
conda config --set channel_priority strict. - Chromosome Naming Contamination: Always verify whether reference FASTA headers use
chr1or1; divergent naming causes silent alignment drops.