Integrative Genomics Viewer (IGV) Exploration & Biological Synthesis of HD Dysregulation
Integrative Genomics Viewer (IGV) Exploration & Biological Synthesis of HD Dysregulation
Loading multi-track BigWig files, narrowPeak calls, and gene models into IGV for visual verification; exploring critical Huntington's Disease loci including BDNF, HTT, and DARPP-32/PPP1R1B; synthesizing the entire multi-day analytical journey into a compelling, publication-ready biological narrative.
1. Visual Inspection in the Integrative Genomics Viewer (IGV)
Statistical metrics and composite heatmaps provide macro-level validation, but inspectable biological confirmation requires examining individual gene loci in IGV.
Step-by-Step IGV Inspection Protocol
- Launch IGV and set the reference genome to Human (hg38).
- Load the four normalized BigWig files:
Control_Rep1.bw,Control_Rep2.bw,HD_Rep1.bw,HD_Rep2.bw. - Load the consensus peak calls:
HD_peaks.narrowPeakandControl_peaks.narrowPeak. - Select all BigWig tracks, right-click, and enable "Group Autoscale". This ensures identical vertical scales ($0 - 45\text{ BPM}$) across all tracks for unbiased visual comparison.
- Navigate to canonical locus coordinates: BDNF (chr11:27,700,000-27,740,000) and PPP1R1B / DARPP-32 (chr17:39,350,000-39,370,000).
2. The Biological Story: Uncovering the Epigenetic Lesion
At the BDNF promoter IV locus, healthy control tracks display robust, sharp H3K4me3 peaks reaching $50\text{ BPM}$, indicating an open, actively transcribed chromatin state that drives continuous synthesis of neurotrophic support for striatal projection neurons. In stark contrast, both Huntington's Disease tracks exhibit a collapse of H3K4me3 signal ($< 8\text{ BPM}$), representing profound promoter de-methylation and transcriptional shutdown.
Combined with our HOMER motif discoveries (SP1 and REST/NRSF enrichment) and DiffBind statistical testing ($\log_2\text{FC} = -2.8$, $\text{FDR} = 4.2 \times 10^{-7}$), we can synthesize the complete causal chain:
- Mutant huntingtin aggregates disrupt the MLL/COMPASS H3K4 methyltransferase complex and recruit HDACs / co-repressors to neuronal promoters.
- Promoter H3K4me3 is progressively lost across post-mitotic striatal and cortical neurons.
- Downstream transcription of survival factors (BDNF, DARPP-32, dopamine receptor $D2$) is silenced.
- Deprived of essential trophic support, striatal medium spiny neurons undergo progressive apoptosis, culminating in the motor and cognitive decline of Huntington's Disease.
Congratulations! You have taken raw, multi-gigabyte FASTQ reads from public archives, performed rigorous quality and alignment filtering, modeled statistical peak differentials, mapped regulatory transcription factor motifs, and translated raw sequencing reads into a publication-grade biological story.
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.