High-throughput chromosome conformation capture assays, such as Hi-C, have shown that the genome is organized into organizational units such as topologically associating domains (TADs), which can impact gene regulatory processes. The sparsity of Hi-C matrices poses a challenge for reliable detection of these units. We present GRiNCH, a constrained matrix-factorization-based approach for simultaneous smoothing and discovery of TADs from sparse contact count matrices. GRiNCH shows superior performance against seven TAD-calling methods and three smoothing methods. GRiNCH is applicable to multiple platforms including SPRITE and HiChIP and can predict novel boundary factors with potential roles in genome organization.

The three-dimensional (3D) organization of the genome plays a critical role in gene regulation for diverse normal and disease processes. High-throughput chromosome conformation capture (3C) assays, such as Hi-C, SPRITE, GAM, and HiChIP, have revealed higher-order organizational units such as topologically associating domains (TADs), which can shape the regulatory landscape governing downstream phenotypes. Analysis of high-throughput 3C data depends on the sequencing depth, which directly affects the resolution and the sparsity of the generated 3D contact count map. Identification of TADs remains a significant challenge due to the sensitivity of existing methods to resolution and sparsity. Here we present GRiNCH, a novel matrix-factorization-based approach for simultaneous TAD discovery and smoothing of contact count matrices from high-throughput 3C data. GRiNCH TADs are enriched in known architectural proteins and chromatin modification signals and are stable to the resolution, and sparsity of the input data. GRiNCH smoothing improves the recovery of structure and significant interactions from low-depth datasets. Furthermore, enrichment analysis of 746 transcription factor motifs in GRiNCH TADs from developmental time-course and cell-line Hi-C datasets predicted transcription factors with potentially novel genome organization roles. GRiNCH is a broadly applicable tool for the analysis of high throughput 3C datasets from a variety of platforms including SPRITE and HiChIP to understand 3D genome organization in diverse biological contexts.

Cancer risk by environmental exposure is modulated by an individual’s genetics and age at exposure. This age-specific period of susceptibility is referred to as the “Window of Susceptibility” (WOS). Rats have a similar WOS for developing breast cancer. A previous study in rat identified an age-specific long-range regulatory interaction for the cancer gene, Pappa, that is associated with breast cancer susceptibility. However, the global role of three-dimensional genome organization and downstream gene expression programs in the WOS is not known. Therefore, we generated Hi-C and RNA-seq data in rat mammary epithelial cells within and outside the WOS. To systematically identify higher-order changes in 3D genome organization, we developed NE-MVNMF that combines network enhancement followed by multitask non-negative matrix factorization. We examined three-dimensional genome organization dynamics at the level of individual loops as well as higher-order domains. Differential chromatin interactions tend to be associated with differentially up-regulated genes with the WOS and recapitulate several human SNP-gene interactions associated with breast cancer susceptibility. Our approach identified genomic blocks of regions with greater overall differences in contact count between the two time points when the cluster assignments change and identified genes and pathways implicated in early carcinogenesis and cancer treatment. Our results suggest that WOS-specific changes in 3D genome organization are linked to transcriptional changes that may influence susceptibility to breast cancer.

Changes in the three-dimensional (3D) structure of the genome are an emerging hallmark of cancer. Cancer-associated copy number variants and single nucleotide polymorphisms promote rewiring of chromatin loops, disruption of topologically associating domains (TADs), active/inactive chromatin state switching, leading to oncogene expression and silencing of tumor suppressors. However, little is known about 3D changes during cancer progression to a chemotherapy-resistant state. We integrated chromatin conformation capture (Hi-C), RNA-seq, and whole-genome sequencing obtained from triple-negative breast cancer patient-derived xenograft primary tumors (UCD52) and carboplatin-resistant samples and found increased short-range (< 2Mb) interactions, chromatin looping, formation of topologically associating domains (TAD), chromatin state switching into a more active state, and amplification of ATP-binding cassette (ABC) transporters. Transcriptome changes suggested the role of long-noncoding RNAs in carboplatin resistance. Rewiring of the 3D genome was associated with TP53, TP63, BATF, FOS-JUN family of transcription factors and led to activation of aggressiveness-, metastasis- and other cancer-related pathways. Integrative analysis highlighted increased ribosome biogenesis and oxidative phosphorylation, suggesting the role of mitochondrial energy metabolism. Our results suggest that 3D genome remodeling may be a key mechanism underlying carboplatin resistance. Citation Format: Mikhail Dozmorov, Maggie Marshall, Narmeen Rashid, Jacqueline Grible, Aaron D. Valentine, Amy Olex, Kavita Murthy, Abhijit Chakraborty, Joaquin Reyna, Daniela Salgado Figueroa, Da-Inn Lee, Brittany Baur, Sushmita Roy, Ferhat Ay, Chuck Harrell. Carboplatin resistance-associated changes in the 3D chromatin landscape of a triple-negative breast cancer Patient-Derived Xenograft [abstract]. In: Proceedings of the 2022 San Antonio Breast Cancer Symposium; 2022 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2023;83(5 Suppl):Abstract nr P1-13-23.