As sequencing depth of chromatin studies continually grows deeper for sensitive profiling of regulatory elements or chromatin spatial structures, aligning and preprocessing of these sequencing data have become the bottleneck for analysis. Here we present Chromap, an ultrafast method for aligning and preprocessing high throughput chromatin profiles. Chromap is comparable to BWA-MEM and Bowtie2 in alignment accuracy and is over 10 times faster than traditional workflows on bulk ChIP-seq/Hi-C profiles and than 10x Genomics’ CellRanger v2.0.0 pipeline on single-cell ATAC-seq profiles.
Chromatin, which consists of DNA and associated proteins, contains genetic information and is a mechanical component of the nucleus. Heterochromatic histone methylation controls nucleus and chromosome stiffness, but the contribution of heterochromatin protein HP1α (CBX5) is unknown. We used a novel HP1α auxin-inducible degron human cell line to rapidly degrade HP1α. Degradation did not alter transcription, local chromatin compaction, or histone methylation, but did decrease chromatin stiffness. Single-nucleus micromanipulation reveals that HP1α is essential to chromatin-based mechanics and maintains nuclear morphology, separate from histone methylation. Further experiments with dimerization-deficient HP1αI165E indicate that chromatin crosslinking via HP1α dimerization is critical, while polymer simulations demonstrate the importance of chromatin-chromatin crosslinkers in mechanics. In mitotic chromosomes, HP1α similarly bolsters stiffness while aiding in mitotic alignment and faithful segregation. HP1α is therefore a critical chromatin-crosslinking protein that provides mechanical strength to chromosomes and the nucleus throughout the cell cycle and supports cellular functions.
The Hi-C technique has been shown to be a promising method to detect structural variations (SVs) in human genomes. However, algorithms that can use Hi-C data for a full-range SV detection have been severely lacking. Current methods can only identify interchromosomal translocations and long-range intrachromosomal SVs (>1 Mb) at less-than-optimal resolution. Therefore, we develop EagleC, a framework that combines deep-learning and ensemble-learning strategies to predict a full range of SVs at high resolution. We show that EagleC can uniquely capture a set of fusion genes that are missed by whole-genome sequencing or nanopore. Furthermore, EagleC also effectively captures SVs in other chromatin interaction platforms, such as HiChIP, Chromatin interaction analysis with paired-end tag sequencing (ChIA-PET), and capture Hi-C. We apply EagleC in more than 100 cancer cell lines and primary tumors and identify a valuable set of high-quality SVs. Last, we demonstrate that EagleC can be applied to single-cell Hi-C and used to study the SV heterogeneity in primary tumors.