Topologically Associating Domains Boundaries Research Articles

TAD border deletion at the Kit locus causes tissue-specific ectopic activation of a neighboring gene

Topologically associated domains (TADs) restrict promoter-enhancer interactions, thereby maintaining the spatiotemporal pattern of gene activity. However, rearrangements of the TADs boundaries do not always lead to significant changes in the activity pattern. Here, we investigated the consequences of the TAD boundaries deletion on the expression of developmentally important genes encoding tyrosine kinase receptors: Kit, Kdr, Pdgfra. We used genome editing in mice to delete the TADs boundaries at the Kit locus and characterized chromatin folding and gene expression in pure cultures of fibroblasts, mast cells, and melanocytes. We found that although Kit is highly active in both mast cells and melanocytes, deletion of the TAD boundary between the Kit and Kdr genes results in ectopic activation only in melanocytes. Thus, the epigenetic landscape, namely the mutual arrangement of enhancers and actively transcribing genes, is important for predicting the consequences of the TAD boundaries removal. We also found that mice without a TAD border between the Kit and Kdr genes have a phenotypic manifestation of the mutation — a lighter coloration. Thus, the data obtained shed light on the principles of interaction between the 3D chromatin organization and epigenetic marks in the regulation of gene activity.

TAD boundary deletion causes PITX2-related cardiac electrical and structural defects

While 3D chromatin organization in topologically associating domains (TADs) and loops mediating regulatory element-promoter interactions is crucial for tissue-specific gene regulation, the extent of their involvement in human Mendelian disease is largely unknown. Here, we identify 7 families presenting a new cardiac entity associated with a heterozygous deletion of 2 CTCF binding sites on 4q25, inducing TAD fusion and chromatin conformation remodeling. The CTCF binding sites are located in a gene desert at 1 Mb from the Paired-like homeodomain transcription factor 2 gene (PITX2). By introducing the ortholog of the human deletion in the mouse genome, we recapitulate the patient phenotype and characterize an opposite dysregulation of PITX2 expression in the sinoatrial node (ectopic activation) and ventricle (reduction), respectively. Chromatin conformation assay performed in human induced pluripotent stem cell-derived cardiomyocytes harboring the minimal deletion identified in family#1 reveals a conformation remodeling and fusion of TADs. We conclude that TAD remodeling mediated by deletion of CTCF binding sites causes a new autosomal dominant Mendelian cardiac disorder.

Abstract 2860: High-resolution mapping of oncogenic structural changes in osteosarcoma

Abstract We have leveraged long-read and spatial sequencing modalities to identify cancer-specific oncogenic alterations for treatment of osteosarcomas. Osteosarcoma is archetypal of cancers driven by aneuploidy and structural rearrangement rather than point mutations. The functional consequences of many SVs stem from their effect on multidimensional genome organization. SVs can enable enhancer hijacking, alter boundaries of topologically associated domains (TADs), and move gene loci to different regulatory compartments to affect gene expression. SVs can also lead to the formation of extrachromosomal DNA amplicons (ecDNA), which can be megabases in size and can incorporate both genes and regulatory elements to result in substantial increases in gene transcription and intratumor heterogeneity. Complex SVs may have profound implications for prognosis and treatment in OS and other cancers yet they remain poorly understood. By incorporating long-read optical genome mapping (OGM) and chromatin conformation capture sequencing (HiC), we studied SVs in osteosarcoma with the goal of elucidating their contribution to cancer development. We used a unique panel of 11 OS patient-derived xenograft (PDX) cell lines which have been previously characterized regarding response to targeted therapies as well as genomically via whole genome (WGS), RNA (RNA-seq), and chromatin accessibility (ATAC-seq) sequencing. To understand the complex genomic reorganization that occurs in osteosarcoma, we integrated OGM and WGS data for all 11 cell lines to describe detailed genome-wide SVs including ecDNA formation. PDX lines assayed using OGM exhibited between 106 and 704 translocations each (versus zero detected in germline controls) and we identified dozens of potential ecDNA amplicons, including one amplicon which may explain exceptionally high expression of YAP1 and BIRC3 transcripts. Further, we used HiC to describe the topological environment resulting from these SVs and observed potential enhancer hijacking events affecting dozens of known oncogenes including MYC and CDK4. HiC interaction maps among six cell lines reveal genomic compartments and TAD boundaries in detail, with 4,167 unique TAD boundaries identified among six cell lines using 50kb genomic bins. Despite the observed complexity of these genomes, we note that 71% of TAD boundaries are consistent in at least four of the six HiC assayed cell lines, and 44% are present in all six lines. Changing TAD boundaries are associated with numerous SVs and/or copy number discrepancies, and are associated with genomic regions that contain genes with dysregulated RNA transcript levels. Overall, OGM and HiC modalities add considerably to our ability to detect oncogenic SVs and the genes they affect. Citation Format: Andrew Scott Clugston, Stan Leung, Eunice Fuentes, Leanne Sayles, Megan Ostrowski, Coco Wu, Christina Curtis, Zhicheng Ma, Yanding Zhao, Georgi K. Marinov, Vijay Ramani, Marcus Breese, Alejandro Sweet-Cordero. High-resolution mapping of oncogenic structural changes in osteosarcoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 2860.

Abstract 1648: RB represses cohesin-dependent loop formation and activates E2F-independent transcription

Abstract The retinoblastoma protein (RB) suppresses the activity of the E2F transcription factor family, controlling cellular proliferation. Recent studies, however, indicate that RB's role in chromatin organization, which is not yet fully understood, might be distinct from this E2F-dependent regulation. Employing chromosome conformation capture at a single nucleosome resolution, we identified RB as a cell cycle-regulated repressor of cohesin-dependent loop formation at topologically associating domain (TAD) boundaries. RB depletion increased the number and size of cohesin-dependent loops and strengthened topologically associating domains (TADs). This phenomenon was specific to the G1 phase and was not observed in the S phase, indicating that this novel RB function is regulated by the cell cycle. Mechanistically, RB showed extensive colocalization with cohesin in the human genome, and it impacted cohesin’s distribution on the chromatin. Active RB reduced cohesin from RB-bound TAD boundaries and decreased cohesin activity therein, as assessed by the increased K105/106 acetylation of the cohesin subunit SMC3. This led to reduced insulation in chromosome conformation capture assays. Importantly, by weakening the insulation activity of the adjacent insulators, RB non-canonically enhanced the expression of non-E2F target genes salient for cell adhesion and extracellular matrix organization. When RB was lost, cells showed a more rapid cellular detachment rate and an elevated migration rate, indicating that this novel RB function controls specific transcriptional programs rather than arbitrary genes to regulate the process of cell adhesion and migration. Overall, we conclude that RB has a central role in the interplay between cell cycle and chromatin organization, by repressing the cohesin-dependent loop formation at TAD boundaries. This RB function safeguards E2F-independent transcriptional programs driven by active enhancers and helps maintain cellular adhesion to the extracellular matrix. Citation Format: Hanjun Lee, Ioanna-Maria Gkotinakou, Badri Krishnan, Nicholas J. Dyson, Michael S. Lawrence, Ioannis Sanidas. RB represses cohesin-dependent loop formation and activates E2F-independent transcription [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 1648.

Identifying TAD-like domains on single-cell Hi-C data by graph embedding and changepoint detection.

Topologically associating domains (TADs) are fundamental building blocks of 3D genome. TAD-like domains in single cells are regarded as the underlying genesis of TADs discovered in bulk cells. Understanding the organization of TAD-like domains helps to get deeper insights into their regulatory functions. Unfortunately, it remains a challenge to identify TAD-like domains on single-cell Hi-C data due to its ultra-sparsity. We propose scKTLD, an in silico tool for the identification of TAD-like domains on single-cell Hi-C data. It takes Hi-C contact matrix as the adjacency matrix for a graph, embeds the graph structures into a low-dimensional space with the help of sparse matrix factorization followed by spectral propagation, and the TAD-like domains can be identified using a kernel-based changepoint detection in the embedding space. The results tell that our scKTLD is superior to the other methods on the sparse contact matrices, including downsampled bulk Hi-C data as well as simulated and experimental single-cell Hi-C data. Besides, we demonstrated the conservation of TAD-like domain boundaries at single-cell level apart from heterogeneity within and across cell types, and found that the boundaries with higher frequency across single cells are more enriched for architectural proteins and chromatin marks, and they preferentially occur at TAD boundaries in bulk cells, especially at those with higher hierarchical levels. scKTLD is freely available at https://github.com/lhqxinghun/scKTLD.

Integrating Dynamic 3D Chromatin Architecture and Gene Expression Alterations Reveal Heterosis in Brassica rapa.

Heterosis plays a significant role in enhancing variety, boosting yield, and raising economic value in crops, but the molecular mechanism is still unclear. We analyzed the transcriptomes and 3D genomes of a hybrid (F1) and its parents (w30 and 082). The analysis of the expression revealed a total of 485 specially expressed genes (SEGs), 173 differentially expressed genes (DEGs) above the parental expression level, more actively expressed genes, and up-regulated DEGs in the F1. Further study revealed that the DEGs detected in the F1 and its parents were mainly involved in the response to auxin, plant hormone signal transduction, DNA metabolic process, purine metabolism, starch, and sucrose metabolism, which suggested that these biological processes may play a crucial role in the heterosis of Brassica rapa. The analysis of 3D genome data revealed that hybrid F1 plants tend to contain more transcriptionally active A chromatin compartments after hybridization. Supplementaryly, the F1 had a smaller TAD (topologically associated domain) genome length, but the number was the highest, and the expression change in activated TAD was higher than that of repressed TAD. More specific TAD boundaries were detected between the parents and F1. Subsequently, 140 DEGs with genomic structural variants were selected as potential candidate genes. We found two DEGs with consistent expression changes in A/B compartments and TADs. Our findings suggested that genomic structural variants, such as TADs and A/B chromatin compartments, may affect gene expression and contribute to heterosis in Brassica rapa. This study provides further insight into the molecular mechanism of heterosis in Brassica rapa.

Enhancer-promoter interactions can form independently of genomic distance and be functional across TAD boundaries.

Topologically Associating Domains (TADs) have been suggested to facilitate and constrain enhancer-promoter interactions. However, the role of TAD boundaries in effectively restricting these interactions remains unclear. Here, we show that a significant proportion of enhancer-promoter interactions are established across TAD boundaries in Drosophila embryos, but that developmental genes are strikingly enriched in intra- but not inter-TAD interactions. We pursued this observation using the twist locus, a master regulator of mesoderm development, and systematically relocated one of its enhancers to various genomic locations. While this developmental gene can establish inter-TAD interactions with its enhancer, the functionality of these interactions remains limited, highlighting the existence of topological constraints. Furthermore, contrary to intra-TAD interactions, the formation of inter-TAD enhancer-promoter interactions is not solely driven by genomic distance, with distal interactions sometimes favored over proximal ones. These observations suggest that other general mechanisms must exist to establish and maintain specific enhancer-promoter interactions across large distances.

The 3D Chromatin Landscape of B-Cell Precursor Childhood Acute Lymphoblastic Leukemia

Whilst the molecular pathogenesis of childhood B-cell precursor (BCP) acute lymphoblastic leukemia (ALL) has been studied extensively, its 3D chromatin landscape remains poorly explored. Genome-wide chromosome conformation capture methods have provided the tools to investigate the different units of chromatin organization, such as transcriptionally active (A) and inactive (B) compartments, topologically associating domains (TADs), and fine-scale chromatin loops and enhancer-promoter interactions. The aim of this study is to elucidate the chromatin architecture and topological gene regulation in childhood BCP ALL. To date, 29 primary patient samples were included, comprising the high hyperdiploid (HeH) (n=11), ETV6:: RUNX1-positive (n=8), BCR:: ABL1-positive (n=2), TCF3:: PBX1-positive (n=2), DUX4-rearranged (n=2), intrachromosomal amplification of chromosome 21 (iAMP21) (n=1), KMT2A-rearranged (n=1), near-haploid (n=1) and near-triploid (n=1) genetic subtypes. Leukemic blast cells obtained at diagnosis were analyzed using Micro-C, a high-resolution variation of Hi-C (average number of total reads = 1.4 billion, highest resolution = 5 kb) combined with pair-end sequencing. Chromatin contact heatmaps were generated for each case using Juicer and Cooler. A/B compartments were identified using FanC at 500 kb resolution and visualized in the software IGV, while TAD calling was performed by Juicer, Domaincaller and Insulation Score. Differential chromatin interaction loop calling was made using Pareidolia and Mustache, and structural variant (SV) calling was carried out by EagleC. Preliminary principal component analysis of the 29 cases, based on the first two eigenvectors of the contact matrix (500 kb resolution), showed that HeH, ETV6:: RUNX1-positive, TCF3:: PBX1-positive and DUX4-rearranged cases each clustered based on their chromatin 3D organization. Furthermore, the A/B compartments of 11 HeH and 8 ETV6:: RUNX1-positive cases were analyzed at 500 kb resolution and compartment shifts among the two subtypes were annotated. A total of 390 shifts were detected, where activating shifts (from B to A compartment) happened more often in HeH (263 shifts) than in ETV6:: RUNX1-positive cases (127 shifts). Analysis of TAD boundary strength at 25 kb resolution revealed that HeH cases displayed significantly weaker boundaries compared to ETV6:: RUNX1-positive cases. TAD boundary strength showed no bias towards the frequently gained or non-gained chromosomes in HeH ALL. By merging individual heatmaps of all HeH and ETV6:: RUNX1-positive cases using Cooler, we created subtype-specific profiles and compared the intensity of chromatin interactions between the two genetic subtypes. Chromatin interaction intensity analysis was then combined with previously published RNA-sequencing data to identify transcriptional dysregulation events that could be associated with chromatin interaction changes. Preliminary results show that there was a chromatin loop missing close to the well-known leukemia-related gene IKZF1 in HeH compared to ETV6:: RUNX1-positive cases; this gene also showed lower expression in the RNA-sequencing data. FLT3 was associated with weakened chromatin interactions and down-regulated in ETV6:: RUNX1-positive cases compared to HeH, in agreement with its known high expression in HeH. Finally, we performed screening of SVs using EagleC and Micro-C heatmaps in HeH and ETV6:: RUNX1-positive samples. Out of the 19 included cases, previous whole-genome sequencing (WGS) data were available for 16. We detected 75 SVs, of which 50 were intrachromosomal rearrangements and 25 were translocations. Micro-C heatmaps allowed visual detection of SVs smaller than 1 Mb and permitted identification of the type of SVs. WGS detected 61% of the SVs found in the HeH samples with Micro-C and 51% of those in the ETV6:: RUNX1-positive cases. In summary, we present the first high-resolution genome-wide map of chromatin 3D organization in pediatric ALL. Our results indicate that different subtypes of childhood BCP ALL have distinct 3D chromatin landscapes and that abnormal chromatin architectures affect the regulation of leukemia-related genes.

Three-dimensional chromatin landscapes in hepatocellular carcinoma associated with hepatitis B virus.

Three-dimensional (3D) chromatin architecture frequently altered in cancer. However, its changes during the pathogenesis of hepatocellular carcinoma (HCC) remained elusive. Hi-C and RNA-seq were applied to study the 3D chromatin landscapes and gene expression of HCC and ANHT. Hi-C Pro was used to generate genome-wide raw interaction matrices, which were normalized via iterative correction (ICE). Moreover, the chromosomes were divided into different compartments according to the first principal component (E1). Furthermore, topologically associated domains (TADs) were visualized via WashU Epigenome Browser. Furthermore, differential expression analysis of ANHT and HCC was performed using the DESeq2 R package. Additionally, dysregulated genes associated with 3D genome architecture altered were confirmed using TCGA, qRT-PCR, immunohistochemistry (IHC), etc. RESULTS: First, the intrachromosomal interactions of chr1, chr2, chr5, and chr11 were significantly different, and the interchromosomal interactions of chr4-chr10, chr13-chr21, chr15-chr22, and chr16-chr19 are remarkably different between ANHT and HCC, which resulted in the up-regulation of TP53I3 and ZNF738 and the down-regulation of APOC3 and APOA5 in HCC. Second, 49 compartment regions on 18 chromosomes have significantly switched (A-B or B-A) during HCC tumorigenesis, contributing to up-regulation of RAP2A. Finally, a tumor-specific TAD boundary located on chr5: 6271000-6478000 and enhancer hijacking were identified in HCC tissues, potentially associated with the elevated expression of MED10, whose expression were associated with poor prognosis of HCC patients. This study demonstrates the crucial role of chromosomal structure variation in HCC oncogenesis and potential novel biomarkers of HCC, laying a foundation for cancer precision medicine development.

Corticosteroid-induced chromatin loop dynamics at the FKBP5 gene.

FKBP5 is a 115-kb-long glucocorticoid-inducible gene implicated in psychiatric disorders. To investigate the complexities of chromatin interaction frequencies at the FKBP5 topologically associated domain (TAD), we deployed 15 one-to-all chromatin capture viewpoints near gene promoters, enhancers, introns, and CTCF-loop anchors. This revealed a "one-TAD-one-gene" structure encompassing the FKBP5 promoter and its enhancers. The FKBP5 promoter and its two glucocorticoid-stimulated enhancers roam the entire TAD while displaying subtle cell type-specific interactomes. The FKBP5 TAD consists of two nested CTCF loops that are coordinated by one CTCF site in the eighth intron of FKBP5 and another beyond its polyadenylation site, 61kb further. Loop extension correlates with transcription increases through the intronic CTCF site. This is efficiently compensated for, since the short loop is restored even under high transcription regimes. The boundaries of the FKBP5 TAD consist of divergent CTCF site patterns, harbor multiple smaller genes, and are resilient to glucocorticoid stimulation. Interestingly, both FKBP5 TAD boundaries harbor H3K27me3-marked heterochromatin blocks that may reinforce them. We propose that cis-acting genetic and epigenetic polymorphisms underlying FKBP5 expression variation are likely to reside within a 240-kb region that consists of the FKBP5 TAD, its left sub-TAD, and both its boundaries.

Topological reorganization and functional alteration of distinct genomic components in gallbladder cancer.

Altered three-dimensional architecture of chromatin influences various genomic regulators and subsequent gene expression in human cancer. However, knowledge of the topological rearrangement of genomic hierarchical layers in cancer is largely limited. Here, by taking advantage of in situ Hi-C, RNA-sequencing, and chromatin immunoprecipitation sequencing (ChIP-seq), we investigated structural reorganization and functional changes in chromosomal compartments, topologically associated domains (TADs), and CCCTC binding factor (CTCF)-mediated loops in gallbladder cancer (GBC) tissues and cell lines. We observed that the chromosomal compartment A/B switch was correlated with CTCF binding levels and gene expression changes. Increased inter-TAD interactions with weaker TAD boundaries were identified in cancer cell lines relative to normal controls. Furthermore, the chromatin short loops and cancer unique loops associated with chromatin remodeling and epithelial-mesenchymal transition activation were enriched in cancer compared with their control counterparts. Cancer-specific enhancer-promoter loops, which contain multiple transcription factor binding motifs, acted as a central element to regulate aberrant gene expression. Depletion of individual enhancers in each loop anchor that connects with promoters led to the inhibition of their corresponding gene expressions. Collectively, our data offer the landscape of hierarchical layers of cancer genome and functional alterations that contribute to the development of GBC.

Multi-feature clustering of CTCF binding creates robustness for loop extrusion blocking and Topologically Associating Domain boundaries

Topologically Associating Domains (TADs) separate vertebrate genomes into insulated regulatory neighborhoods that focus genome-associated processes. TADs are formed by Cohesin-mediated loop extrusion, with many TAD boundaries consisting of clustered binding sites of the CTCF insulator protein. Here we determine how this clustering of CTCF binding contributes to the blocking of loop extrusion and the insulation between TADs. We identify enrichment of three features of CTCF binding at strong TAD boundaries, consisting of strongly bound and closely spaced CTCF binding peaks, with a further enrichment of DNA-binding motifs within these peaks. Using multi-contact Nano-C analysis in cells with normal and perturbed CTCF binding, we establish that individual CTCF binding sites contribute to the blocking of loop extrusion, but in an incomplete manner. When clustered, individual CTCF binding sites thus create a stepwise insulation between neighboring TADs. Based on these results, we propose a model whereby multiple instances of temporal loop extrusion blocking create strong insulation between TADs.

Epigenetic germline variants predict cancer prognosis and risk and distribute uniquely in topologically associating domains

Background: Methylation quantitative trait loci (meQTLs) associate with different levels of local DNA methylation in cancers. Here, we investigated whether the distribution of cancer meQTLs reflected functional organization of the genome in the form of chromatin topologically associated domains (TADs) and evaluated whether cancer meQTLs near known driver genes have the potential to influence cancer risk or progression. Methods: Published cancer meQTLs were analyzed according to their location in transcriptionally active or inactive TADs and TAD boundary regions. Cancer meQTLs near known cancer genes were analyzed for association with cancer risk in the UKBioBank and prognosis in The Cancer Genome Atlas (TCGA). Results: In TAD boundary regions, the density of cancer meQTLs was higher near inactive TADs. Furthermore, we observed an enrichment of cancer meQTLs in active TADs near tumor suppressors, whereas there was a depletion of such meQTLs near oncogenes. Several meQTLs were associated with cancer risk in the UKBioBank, and we were able to reproduce breast cancer risk associations in the DRIVE cohort. Survival analysis in TCGA implicated a number of meQTLs in 13 tumor types. In 10 of these, polygenic cancer meQTL scores were associated with increased hazard in a CoxPH analysis. Risk and survival-associated meQTLs tended to affect cancer genes involved in DNA damage repair and cellular adhesion and reproduced cancer-specific associations reported in prior literature. Conclusions: This study provides evidence that genetic variants that influence local DNA methylation are affected by chromatin structure and can impact tumor evolution.

CATAD: exploring topologically associating domains from an insight of core-attachment structure.

Identifying topologically associating domains (TADs), which are considered as the basic units of chromosome structure and function, can facilitate the exploration of the 3D-structure of chromosomes. Methods have been proposed to identify TADs by detecting the boundaries of TADs or identifying the closely interacted regions as TADs, while the possible inner structure of TADs is seldom investigated. In this study, we assume that a TAD is composed of a core and its surrounding attachments, and propose a method, named CATAD, to identify TADs based on the core-attachment structure model. In CATAD, the cores of TADs are identified based on the local density and cosine similarity, and the surrounding attachments are determined based on boundary insulation. CATAD was applied to the Hi-C data of two human cell lines and two mouse cell lines, and the results show that the boundaries of TADs identified by CATAD are significantly enriched by structural proteins, histone modifications, transcription start sites and enzymes. Furthermore, CATAD outperforms other methods in many cases, in terms of the average peak, boundary tagged ratio and fold change. In addition, CATAD is robust and rarely affected by the different resolutions of Hi-C matrices. Conclusively, identifying TADs based on the core-attachment structure is useful, which may inspire researchers to explore TADs from the angles of possible spatial structures and formation process.

Abstract A052: Chromatin conformational changes in prostate cancer progression

Abstract Introduction. Prostate cancer (PCa) is the second most common cancer in men worldwide. Globally, it is among the top three leading causes of cancer-related mortality in men. Although many patients carry an indolent form of PCa, a subset of patients develops the aggressive form of PCa which is often treated by androgen depletion therapy. Despite treatment, the disease, in a subgroup of patients, progresses to develop castration resistant prostate cancer (CRPC). Despite the ongoing efforts, the underlying progression mechanisms are not yet fully understood. A more comprehensive understanding may offer translational benefits. Thus, using patient sample-derived high-throughput sequencing data from multiple epigenetic data layers together with expression levels of genes and proteins, we aim to uncover aberrant regulatory mechanisms and their role in PCa progression. Methods. We characterize the PCa progression-related conformational chromatin changes by applying Hi-C to clinical samples from a cohort of 18 patients, divided into three evenly-sized groups: benign prostatic hyperplasia (BPH), untreated prostate cancer (PC), and CRPC. We further characterize our findings by integrating chromatin accessibility, DNA methylation, gene and protein expression data from matching samples. Results. Looking at active and inactive (A/B) compartments, as inferred by Hi-C data, we observe that they are in high agreement across the sample groups. Integrating these A/B compartments with DNA-methylation data, we observe an elevated and decreased genome-wide methylation pattern respectively. Looking at topologically associating domain (TAD) boundaries, we observe a high agreement across sample groups as 70% of the TAD boundaries are shared by all groups. Focusing in on the chromatin loops, we observe more heterogeneity compared to the overall chromatin conformation as 38% of the loops are sample-specific. Detecting differential loops followed by their integration with gene expression data for genes overlapping the differential loop anchors, we observe an association between differential loop strength and gene expression as PC develops. Further integration of the differential loops with accessible chromatin regions, as inferred by ATAC-seq data, gene expression, and publicly available transcription factor (TF) binding data, we observe that distal enhancers of up-regulated genes in PC are enriched for PCa-relevant TFs such as FOXA1, HOXB13, and GATA2, but we see no enrichment for up-regulated loops in PC. In CRPC, these same TFs are enriched for up-regulated loops, but this does not seem to affect the gene expression. Conclusion. Our results suggest while PC samples undergo a more uniform set of chromatin conformational changes, CRPC samples indicate elevated heterogeneity in their acquired conformational changes potentially contributing to a more aberrant disease progression trajectory. Citation Format: Tuomo Virtanen, Ebrahim Afyounian, Riikka Nurminen, Teuvo Tammela, Alfonso Urbanucci, Juha Kesseli, Kirsi Rautajoki, Tapio Visakorpi, Matti Nykter. Chromatin conformational changes in prostate cancer progression [abstract]. In: Proceedings of the AACR Special Conference: Advances in Prostate Cancer Research; 2023 Mar 15-18; Denver, Colorado. Philadelphia (PA): AACR; Cancer Res 2023;83(11 Suppl):Abstract nr A052.

Morphine Re-arranges Chromatin Spatial Architecture of Primate Cortical Neurons

The expression of linear DNA sequence is precisely regulated by the three-dimensional (3D) architecture of chromatin. Morphine-induced aberrant gene networks of neurons have been extensively investigated; however, how morphine impacts the 3D genomic architecture of neurons is still unknown. Here, we applied digestion-ligation-only high-throughput chromosome conformation capture (DLO Hi-C) technology to investigate the effects of morphine on the 3D chromatin architecture of primate cortical neurons. After receiving continuous morphine administration for 90 days on rhesus monkeys, we discovered that morphine re-arranged chromosome territories, with a total of 391 segmented compartments being switched. Morphine altered over half of the detected topologically associated domains (TADs), most of which exhibited a variety of shifts, followed by separating and fusing types. Analysis of the looping events at kilobase-scale resolution revealed that morphine increased not only the number but also the length of differential loops. Moreover, all identified differentially expressed genes from the RNA sequencing data were mapped to the specific TAD boundaries or differential loops, and were further validated for changed expression. Collectively, an altered 3D genomic architecture of cortical neurons may regulate the gene networks associated with morphine effects. Our finding provides critical hubs connecting chromosome spatial organization and gene networks associated with the morphine effects in humans.

LPAD: using network construction and label propagation to detect topologically associating domains from Hi-C data.

With the development of chromosome conformation capture technique, the study of spatial conformation of a genome based on Hi-C technique has made a quantum leap. Previous studies reveal that genomes are folded into hierarchy of three-dimensional (3D) structures associated with topologically associating domains (TADs), and detecting TAD boundaries is of great significance in the chromosome-level analysis of 3D genome architecture. In this paper, we propose a novel TAD identification method, LPAD, which first extracts node correlations from global interactions of chromosomes based on the random walk with restart and then builds an undirected graph from Hi-C contact matrix. Next, LPAD designs a label propagation-based approach to discover communities and generates TADs. Experimental results verify the effectiveness and quality of TAD detections compared with existing methods. Furthermore, experimental evaluation of chromatin immunoprecipitation sequencing data shows that LPAD performs high enrichment of histone modifications remarkably nearby the TAD boundaries, and these results demonstrate LPAD's advantages on TAD identification accuracy.

Topologically associating domain boundaries are required for normal genome function

Topologically associating domain (TAD) boundaries partition the genome into distinct regulatory territories. Anecdotal evidence suggests that their disruption may interfere with normal gene expression and cause disease phenotypes1–3, but the overall extent to which this occurs remains unknown. Here we demonstrate that targeted deletions of TAD boundaries cause a range of disruptions to normal in vivo genome function and organismal development. We used CRISPR genome editing in mice to individually delete eight TAD boundaries (11–80 kb in size) from the genome. All deletions examined resulted in detectable molecular or organismal phenotypes, which included altered chromatin interactions or gene expression, reduced viability, and anatomical phenotypes. We observed changes in local 3D chromatin architecture in 7 of 8 (88%) cases, including the merging of TADs and altered contact frequencies within TADs adjacent to the deleted boundary. For 5 of 8 (63%) loci examined, boundary deletions were associated with increased embryonic lethality or other developmental phenotypes. For example, a TAD boundary deletion near Smad3/Smad6 caused complete embryonic lethality, while a deletion near Tbx5/Lhx5 resulted in a severe lung malformation. Our findings demonstrate the importance of TAD boundary sequences for in vivo genome function and reinforce the critical need to carefully consider the potential pathogenicity of noncoding deletions affecting TAD boundaries in clinical genetics screening.

CTCF is a DNA-tension-dependent barrier to cohesin-mediated loop extrusion

In eukaryotes, genomic DNA is extruded into loops by cohesin1. By restraining this process, the DNA-binding protein CCCTC-binding factor (CTCF) generates topologically associating domains (TADs)2,3 that have important roles in gene regulation and recombination during development and disease1,4–7. How CTCF establishes TAD boundaries and to what extent these are permeable to cohesin is unclear8. Here, to address these questions, we visualize interactions of single CTCF and cohesin molecules on DNA in vitro. We show that CTCF is sufficient to block diffusing cohesin, possibly reflecting how cohesive cohesin accumulates at TAD boundaries, and is also sufficient to block loop-extruding cohesin, reflecting how CTCF establishes TAD boundaries. CTCF functions asymmetrically, as predicted; however, CTCF is dependent on DNA tension. Moreover, CTCF regulates cohesin’s loop-extrusion activity by changing its direction and by inducing loop shrinkage. Our data indicate that CTCF is not, as previously assumed, simply a barrier to cohesin-mediated loop extrusion but is an active regulator of this process, whereby the permeability of TAD boundaries can be modulated by DNA tension. These results reveal mechanistic principles of how CTCF controls loop extrusion and genome architecture.

Abstract 214: Patient-specific mapping of oncogenic structural changes in pediatric osteosarcoma

Abstract We have leveraged long-read and spatial sequencing modalities to identify patient-specific oncogenic alterations for treatment of pediatric osteosarcomas. Osteosarcoma (OS) is the most common primary bone malignancy in children and young adults, with over 500 cases identified each year in the United States. OS is archetypal of cancers driven by aneuploidy and structural rearrangement rather than point mutations. These structural variants (SVs) can alter gene expression stoichiometrically, but the functional consequences of many SVs stem from their effect on multidimensional genome organization. SVs can enable enhancer hijacking, alter boundaries of topologically associated domains (TADs), and move gene loci to different regulatory compartments affecting transcription. SVs can also form extrachromosomal DNA amplicons (ecDNA), which can be megabases in size and can incorporate both genes and regulatory elements to result in substantial increases in gene transcription and intratumor heterogeneity. Complex SVs may have profound implications for prognosis and treatment in OS and other cancers. By incorporating long-read optical genome mapping (OGM) and chromatin conformation capture sequencing (HiC), our goal is to detect complex tumor-specific SVs that will inform and refine pediatric OS subtypes, revealing new tumor-specific chemotherapy treatment options. We have generated a unique panel of 11 OS patient-derived xenograft (PDX) cell lines which have been previously characterized regarding response to targeted therapies using whole genome (WGS), RNA (RNA-seq), and chromatin accessibility (ATAC-seq) sequencing. To understand the complex genomic reorganization that occurs in OS, we have integrated OGM and WGS data for six cell lines, and OGM, WGS, and HiC data for two. SVs were detected using OGM, HiC, and WGS modalities, and we observed 59-61% concordance between SVs called using OGM and HiC, versus 22-45% when comparing with WGS calls. PDX lines assayed using OGM exhibit between 5,135 and 7,514 SVs each, and translocation numbers ranged from 212 to 1,422. These lines exhibit between dozens of potential ecDNA amplicons, including one amplicon which may explain that line’s specific and exceptionally high expression of YAP1 and BIRC2 transcripts. HiC interaction maps for two PDX lines reveal genomic compartments and TAD boundaries in detail, with 3,713 unique TAD boundaries identified between them using 50kb genomic bins. We observe that 60% of these boundaries are consistent between these lines, and 27% exhibit differences that suggest boundary loss or gain. These changing boundaries are associated with numerous SVs and/or copy number discrepancies, and with genomic regions containing genes with dysregulated RNA transcript levels. Incorporating OGM and HiC modalities adds substantially to our ability to detect oncogenic SVs and the genes they affect. Citation Format: Andrew S. Clugston, Megan Ostrowski, Coco Wu, Leanne Sayles, Stanley G. Leung, Vijay Ramani, E. Alejandro Sweet-Cordero, Marcus R. Breese. Patient-specific mapping of oncogenic structural changes in pediatric osteosarcoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 214.