Tet-mediated DNA methylation dynamics affect chromosome organization

DNA Methylation is a significant epigenetic modification that can modulate chromosome states, but its role in orchestrating chromosome organization has not been well elucidated. Here we systematically assessed the effects of DNA Methylation on chromosome organization with a multi-omics strategy to capture DNA Methylation and high-order chromosome interaction simultaneously on mouse embryonic stem cells with DNA methylation dioxygenase Tet triple knock-out (Tet-TKO). Globally, upon Tet-TKO, we observed weakened compartmentalization, corresponding to decreased methylation differences between CpG island (CGI) rich and poor domains. Tet-TKO could also induce hypermethylation for the CTCF binding peaks in TAD boundaries and chromatin loop anchors. Accordingly, CTCF peak generally weakened upon Tet-TKO, which results in weakened TAD structure and depletion of long-range chromatin loops. Genes that lost enhancer-promoter looping upon Tet-TKO showed DNA hypermethylation in their gene bodies, which may compensate for the disruption of gene expression. We also observed distinct effects of Tet1 and Tet2 on chromatin organization and increased DNA methylation correlation on spatially interacted fragments upon Tet inactivation. Our work showed the broad effects of Tet inactivation and DNA methylation dynamics on chromosome organization.

DNA fragility at topologically associated domain boundaries is promoted by alternative DNA secondary structure and topoisomerase II activity

CCCTC-binding factor (CTCF) binding sites are hotspots of genome instability. Although many factors have been associated with CTCF binding site fragility, no study has integrated all fragility-related factors to understand the mechanism(s) of how they work together. Using an unbiased, genome-wide approach, we found that DNA double-strand breaks (DSBs) are enriched at strong, but not weak, CTCF binding sites in five human cell types. Energetically favorable alternative DNA secondary structures underlie strong CTCF binding sites. These structures coincided with the location of topoisomerase II (TOP2) cleavage complex, suggesting that DNA secondary structure acts as a recognition sequence for TOP2 binding and cleavage at CTCF binding sites. Furthermore, CTCF knockdown significantly increased DSBs at strong CTCF binding sites and at CTCF sites that are located at topologically associated domain (TAD) boundaries. TAD boundary-associated CTCF sites that lost CTCF upon knockdown displayed increased DSBs when compared to the gained sites, and those lost sites are overrepresented with G-quadruplexes, suggesting that the structures act as boundary insulators in the absence of CTCF, and contribute to increased DSBs. These results model how alternative DNA secondary structures facilitate recruitment of TOP2 to CTCF binding sites, providing mechanistic insight into DNA fragility at CTCF binding sites.

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.

Dissection of a CTCF topological boundary uncovers principles of enhancer-oncogene regulation

Enhancer-gene communication is dependent on topologically associating domains (TADs) and boundaries enforced by the CCCTC-binding factor (CTCF) insulator, but the underlying structures and mechanisms remain controversial. Here, we investigate a boundary that typically insulates fibroblast growth factor (FGF) oncogenes but is disrupted by DNA hypermethylation in gastrointestinal stromal tumors (GISTs). The boundary contains an array of CTCF sites that enforce adjacent TADs, one containing FGF genes and the other containing ANO1 and its putative enhancers, which are specifically active in GIST and its likely cell of origin. We show that coordinate disruption of four CTCF motifs in the boundary fuses the adjacent TADs, allows the ANO1 enhancer to contact FGF3, and causes its robust induction. High-resolution micro-C maps reveal specific contact between transcription initiation sites in the ANO1 enhancer and FGF3 promoter that quantitatively scales with FGF3 induction such that modest changes in contact frequency result in strong changes in expression, consistent with a causal relationship.

RPL35A drives ovarian cancer progression by promoting the binding of YY1 to CTCF promoter

Ovarian cancer is one of the most common gynaecological malignancies with poor prognosis and lack of effective treatment. The improvement of the situation of ovarian cancer urgently requires the exploration of its molecular mechanism to develop more effective molecular targeted drugs. In this study, the role of human ribosomal protein l35a (RPL35A) in ovarian cancer was explored in vitro and in vivo. Our data identified that RPL35A expression was abnormally elevated in ovarian cancer. Clinically, high expression of RPL35A predicted short survival and poor TNM staging in patients with ovarian cancer. Functionally, RPL35A knock down inhibited ovarian cancer cell proliferation and migration, enhanced apoptosis, while overexpression had the opposite effect. Mechanically, RPL35A promoted the direct binding of transcription factor YY1 to CTCF in ovarian cancer cells. Consistently, RPL35A regulated ovarian cancer progression depending on CTCF in vitro and in vivo. Furthermore, RPL35A affected the proliferation and apoptosis of ovarian cancer cells through PPAR signalling pathway. In conclusion, RPL35A drove ovarian cancer progression by promoting the binding of YY1 and CTCF promoter, and inhibiting this process may be an effective strategy for targeted therapy of this disease.

DARDN: A Deep-Learning Approach for CTCF Binding Sequence Classification and Oncogenic Regulatory Feature Discovery

Characterization of gene regulatory mechanisms in cancer is a key task in cancer genomics. CCCTC-binding factor (CTCF), a DNA binding protein, exhibits specific binding patterns in the genome of cancer cells and has a non-canonical function to facilitate oncogenic transcription programs by cooperating with transcription factors bound at flanking distal regions. Identification of DNA sequence features from a broad genomic region that distinguish cancer-specific CTCF binding sites from regular CTCF binding sites can help find oncogenic transcription factors in a cancer type. However, the presence of long DNA sequences without localization information makes it difficult to perform conventional motif analysis. Here, we present DNAResDualNet (DARDN), a computational method that utilizes convolutional neural networks (CNNs) for predicting cancer-specific CTCF binding sites from long DNA sequences and employs DeepLIFT, a method for interpretability of deep learning models that explains the model's output in terms of the contributions of its input features. The method is used for identifying DNA sequence features associated with cancer-specific CTCF binding. Evaluation on DNA sequences associated with CTCF binding sites in T-cell acute lymphoblastic leukemia (T-ALL) and other cancer types demonstrates DARDN's ability in classifying DNA sequences surrounding cancer-specific CTCF binding from control constitutive CTCF binding and identifying sequence motifs for transcription factors potentially active in each specific cancer type. We identify potential oncogenic transcription factors in T-ALL, acute myeloid leukemia (AML), breast cancer (BRCA), colorectal cancer (CRC), lung adenocarcinoma (LUAD), and prostate cancer (PRAD). Our work demonstrates the power of advanced machine learning and feature discovery approach in finding biologically meaningful information from complex high-throughput sequencing data.

STAG2: Computational Analysis of Missense Variants Involved in Disease

The human STAG2 protein is an essential component of the cohesin complex involved in cellular processes of gene expression, DNA repair, and genomic integrity. Somatic mutations in the STAG2 sequence have been associated with various types of cancer, while congenital variants have been linked to developmental disorders such as Mullegama-Klein-Martinez syndrome, X-linked holoprosencephaly-13, and Cornelia de Lange syndrome. In the cohesin complex, the direct interaction of STAG2 with DNA and with NIPBL, RAD21, and CTCF proteins has been described. The function of STAG2 within the complex is still unknown, but it is related to its DNA binding capacity and is modulated by its binding to the other three proteins. Every missense variant described for STAG2 is located in regions involved in one of these interactions. In the present work, we model the structure of 12 missense variants described for STAG2, as well as two other variants of NIPBl and two of RAD21 located at STAG2 interaction zone, and then analyze their behavior through molecular dynamic simulations, comparing them with the same simulation of the wild-type protein. This will allow the effects of variants to be rationalized at the atomic level and provide clues as to how STAG2 functions in the cohesin complex.

Differential DNA methylation and CTCF binding between the ESR1 promoter a of MCF-7 and MDA-MB-231 breast cancer cells

BACKGROUND

ESR1 is expressed by 60-70% of breast tumours. it's a good prognosis factor and the target of hormone therapy. Optimization of ESR1 reactivation therapy is currently ongoing. Here we probe if the transcription factor CTCF plays a role in the differential expression of ESR1 in the breast cancer cell lines MCF-7 (ESR1+) and MDA-MB-231 (ESR1-).

METHODS AND RESULTS

Knockdown of CTCF in MCF-7 resulted in decreased ESR1 gene expression. CTCF binds to the promoter of ESR1 in MCF-7 but not in MDA-MB-231 cells. CTCF ESR1 binding sites are unmethylated in MCF7 but methylated in MDA-MB-231 cells.

CONCLUSION

ESR1 expression in MCF7 cells is dependent on CTCF expression. CTCF can bind to specific regions of the promotor of ESR1 gene in MCF-7 cells but not in MDA-MB-231 cells, this correlates with the methylation status of these regions and could be involved in the transcriptional regulation of ESR1.

SENP1 knockdown-mediated CTCF SUMOylation enhanced its stability and alleviated lipopolysaccharide-evoked inflammatory injury in human lung fibroblasts via regulation of FOXA2 transcription

BACKGROUND

Excessive inflammation is the main cause of treatment failure in neonatal pneumonia (NP). CCCTC-binding factor (CTCF) represents an important node in various inflammatory diseases. In the present study, we tried to clarify the function and underlying molecular mechanism of CTCF on an in vitro cellular model of NP, which was generated by simulating the human lung fibroblast cell line WI-38 with lipopolysaccharide (LPS).

METHODS

The SUMOylation level and protein interaction were verified by Co-immunoprecipitation assay. Cell viability was measured by Cell Counting Kit-8 assay. Inflammatory factors were examined by Enzyme-linked immunosorbent assay. Cell apoptosis was evaluated by TUNEL assay. The binding activity of CTCF to target promoter was tested by chromatin immunoprecipitation and luciferase reporter assay.

RESULTS

LPS treatment restrained cell viability, promoted the production of inflammatory factors, and enhanced cell apoptosis. CTCF overexpression played anti-inflammatory and anti-apoptotic roles. Furthermore, CTCF was modified by SUMOylation with small ubiquitin-like modifier protein 1 (SUMO1). Interfering with sumo-specific protease 1 (SENP1) facilitated CTCF SUMOylation and protein stability, thus suppressing LPS-evoked inflammatory and apoptotic injuries. Moreover, CTCF could bind to the forkhead box protein A2 (FOXA2) promoter region to promote FOXA2 expression. The anti-inflammatory and anti-apoptotic roles of CTCF are associated with FOXA2 activation. In addition, SENP1 knockdown increased FOXA2 expression by enhancing the abundance and binding ability of CTCF.

CONCLUSIONS

SUMOylation of CTCF by SENP1 knockdown enhanced its protein stability and binding ability and it further alleviated LPS-evoked inflammatory injury in human lung fibroblasts by positively regulating FOXA2 transcription.

Transcriptional regulation of FACT involves Coordination of chromatin accessibility and CTCF binding

Histone chaperone FACT (facilitates chromatin transcription) is well known to promote chromatin recovery during transcription. However, the mechanism how FACT regulates genome-wide chromatin accessibility and transcription factor binding has not been fully elucidated. Through loss-of-function studies, we show here that FACT component Ssrp1 is required for DNA replication and DNA damage repair and is also essential for progression of cell phase transition and cell proliferation in mouse embryonic fibroblast cells. On the molecular level, absence of the Ssrp1 leads to increased chromatin accessibility, enhanced CTCF binding, and a remarkable change in dynamic range of gene expression. Our study thus unequivocally uncovers a unique mechanism by which FACT complex regulates transcription by coordinating genome-wide chromatin accessibility and CTCF binding.

CTCF regulates hepatitis B virus cccDNA chromatin topology

Hepatitis B Virus (HBV) is a small DNA virus that replicates via an episomal covalently closed circular DNA (cccDNA) that serves as the transcriptional template for viral mRNAs. The host protein, CCCTC-binding factor (CTCF), is a key regulator of cellular transcription by maintaining epigenetic boundaries, nucleosome phasing, stabilisation of long-range chromatin loops and directing alternative exon splicing. We previously reported that CTCF binds two conserved motifs within Enhancer I of the HBV genome and represses viral transcription, however, the underlying mechanisms were not identified. We show that CTCF depletion in cells harbouring cccDNA-like HBV molecules and in de novo infected cells resulted in an increase in spliced transcripts, which was most notable in the abundant SP1 spliced transcript. In contrast, depletion of CTCF in cell lines with integrated HBV DNA had no effect on the abundance of viral transcripts and in line with this observation there was limited evidence for CTCF binding to viral integrants, suggesting that CTCF-regulation of HBV transcription is specific to episomal cccDNA. Analysis of HBV chromatin topology by Assay for Transposase Accessible Chromatin Sequencing (ATAC-Seq) revealed an accessible region spanning Enhancers I and II and the basal core promoter (BCP). Mutating the CTCF binding sites within Enhancer I resulted in a dramatic rearrangement of chromatin accessibility where the open chromatin region was no longer detected, indicating loss of the phased nucleosome up- and down-stream of the HBV enhancer/BCP. These data demonstrate that CTCF functions to regulate HBV chromatin conformation and nucleosomal positioning in episomal maintained cccDNA, which has important consequences for HBV transcription regulation.

Antigen exposure reshapes chromatin architecture in central memory CD8+ T cells and imprints enhanced recall capacity

CD62L+ central memory CD8+ T (TCM) cells provide enhanced protection than naive cells; however, the underlying mechanism, especially the contribution of higher-order genomic organization, remains unclear. Systematic Hi-C analyses reveal that antigen-experienced CD8+ T cells undergo extensive rewiring of chromatin interactions (ChrInt), with TCM cells harboring specific interaction hubs compared with naive CD8+ T cells, as observed at cytotoxic effector genes such as Ifng and Tbx21. TCM cells also acquire de novo CTCF (CCCTC-binding factor) binding sites, which are not only strongly associated with TCM-specific hubs but also linked to increased activities of local gene promoters and enhancers. Specific ablation of CTCF in TCM cells impairs rapid induction of genes in cytotoxic program, energy supplies, transcription, and translation by recall stimulation. Therefore, acquisition of CTCF binding and ChrInt hubs by TCM cells serves as a chromatin architectural basis for their transcriptomic dynamics in primary response and for imprinting the code of "recall readiness" against secondary challenge.

Outward-oriented sites within clustered CTCF boundaries are key for intra-TAD chromatin interactions and gene regulation

CTCF plays an important role in 3D genome organization by adjusting the strength of chromatin insulation at TAD boundaries, where clustered CBS (CTCF-binding site) elements are often arranged in a tandem array with a complex divergent or convergent orientation. Here, using Pcdh and HOXD loci as a paradigm, we look into the clustered CTCF TAD boundaries and find that, counterintuitively, outward-oriented CBS elements are crucial for inward enhancer-promoter interactions as well as for gene regulation. Specifically, by combinatorial deletions of a series of putative enhancer elements in mice in vivo or CBS elements in cultured cells in vitro, in conjunction with chromosome conformation capture and RNA-seq analyses, we show that deletions of outward-oriented CBS elements weaken the strength of long-distance intra-TAD promoter-enhancer interactions and enhancer activation of target genes. Our data highlight the crucial role of outward-oriented CBS elements within the clustered CTCF TAD boundaries in developmental gene regulation and have interesting implications on the organization principles of clustered CTCF sites within TAD boundaries.

Involvement of CCCTC-binding factor in epigenetic regulation of cancer

A major global health burden continues to be borne by the complex and multifaceted disease of cancer. Epigenetic changes, which are essential for the emergence and spread of cancer, have drawn a huge amount of attention recently. The CCCTC-binding factor (CTCF), which takes part in a wide range of cellular processes including genomic imprinting, X chromosome inactivation, 3D chromatin architecture, local modifications of histone, and RNA polymerase II-mediated gene transcription, stands out among the diverse array of epigenetic regulators. CTCF not only functions as an architectural protein but also modulates DNA methylation and histone modifications. Epigenetic regulation of cancer has already been the focus of plenty of studies. Understanding the role of CTCF in the cancer epigenetic landscape may lead to the development of novel targeted therapeutic strategies for cancer. CTCF has already earned its status as a tumor suppressor gene by acting like a homeostatic regulator of genome integrity and function. Moreover, CTCF has a direct effect on many important transcriptional regulators that control the cell cycle, apoptosis, senescence, and differentiation. As we learn more about CTCF-mediated epigenetic modifications and transcriptional regulations, the possibility of utilizing CTCF as a diagnostic marker and therapeutic target for cancer will also increase. Thus, the current review intends to promote personalized and precision-based therapeutics for cancer patients by shedding light on the complex interplay between CTCF and epigenetic processes.

Transcription induces context-dependent remodeling of chromatin architecture during differentiation

Metazoan chromosomes are organized into discrete spatial domains (TADs), believed to contribute to the regulation of transcriptional programs. Despite extensive correlation between domain organization and gene activity, a direct mechanistic link is unclear, with perturbation studies often showing little effect. To follow chromatin architecture changes during development, we used Capture Hi-C to interrogate the domains around key differentially expressed genes during mouse thymocyte maturation, uncovering specific remodeling events. Notably, one TAD boundary was broadened to accommodate RNA polymerase elongation past the border, and subdomains were formed around some activated genes without changes in CTCF binding. The ectopic induction of some genes was sufficient to recapitulate domain formation in embryonic stem cells, providing strong evidence that transcription can directly remodel chromatin structure. These results suggest that transcriptional processes drive complex chromosome folding patterns that can be important in certain genomic contexts.

XPF interacts with TOP2B for R-loop processing and DNA looping on actively transcribed genes

Co-transcriptional RNA-DNA hybrids can not only cause DNA damage threatening genome integrity but also regulate gene activity in a mechanism that remains unclear. Here, we show that the nucleotide excision repair factor XPF interacts with the insulator binding protein CTCF and the cohesin subunits SMC1A and SMC3, leading to R-loop-dependent DNA looping upon transcription activation. To facilitate R-loop processing, XPF interacts and recruits with TOP2B on active gene promoters, leading to double-strand break accumulation and the activation of a DNA damage response. Abrogation of TOP2B leads to the diminished recruitment of XPF, CTCF, and the cohesin subunits to promoters of actively transcribed genes and R-loops and the concurrent impairment of CTCF-mediated DNA looping. Together, our findings disclose an essential role for XPF with TOP2B and the CTCF/cohesin complex in R-loop processing for transcription activation with important ramifications for DNA repair-deficient syndromes associated with transcription-associated DNA damage.

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, 61 kb 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.

Long non-coding RNA MIR22HG suppresses the chondrogenic differentiation of human adipose-derived stem cells by interacting with CTCF to upregulate CRLF1

Improved chondrogenic differentiation of mesenchymal stem cells (MSCs) by genetic regulation is a potential method for regenerating articular cartilage. LncRNA MIR22HG has been proven to accelerate osteogenic differentiation, but the regulation mechanism of chondrogenic differentiation is still unclear. Human adipose-derived stem cells (hADSCs) have been widely utilised for bone tissue engineering applications. The present study aimed to examine the effect of MIR22HG on the chondrogenic differentiation of hADSCs. The results confirmed that MIR22HG was downregulated in the process of chondrogenic differentiation. Subsequently, gain- and loss-of-function of MIR22HG experiments showed that the overexpression of MIR22HG suppressed the deposition of cartilage matrix proteoglycans and decreased the expression of cartilage-related markers (e.g. Sox9, ACAN and Col2A1), whereas the knockdown of MIR22HG had the opposite effect. MIR22HG could bind to CTCF (CCCTC-binding factor), and CTCF could bind to the CRLF1 (cytokine receptor-like factor 1) promoter and upregulate CRLF1 gene expression. Besides, inhibition of CRLF1 can reverse the effect of MIR22HG on cell chondrogenic differentiation of hADSCs. Taken together, our outcomes reveal that MIR22HG suppressed chondrogenic differentiation by interaction with CTCF to stabilise CRLF1.

CTCF coordinates cell fate specification via orchestrating regulatory hubs with pioneer transcription factors

CCCTC-binding factor (CTCF), a ubiquitously expressed architectural protein, has emerged as a key regulator of cell identity gene transcription. However, the precise molecular mechanism underlying specialized functions of CTCF remains elusive. Here, we investigate the mechanism through integrative analyses of primary hepatocytes, myocytes, and B cells from mouse and human. We demonstrate that CTCF cooperates with lineage-specific pioneer transcription factors (TFs), including MyoD, FOXA, and PU.1, to control cell identity at 1D and 3D levels. At the 1D level, pioneer TFs facilitate lineage-specific CTCF occupancy via opening chromatin. At the 3D level, CTCF and pioneer TFs form regulatory hubs to govern the expression of cell identity genes. This mechanism is validated using MyoD-null mice, CTCF knockout mice, and CRISPR editing during myogenic differentiation. Collectively, these findings uncover a general mechanism whereby CTCF acts as a cell identity cofactor to control cell identity genes via orchestrating regulatory hubs with pioneer TFs.

Decitabine disrupts EBV genomic epiallele DNA methylation patterns around CTCF binding sites to increase chromatin accessibility and lytic transcription in gastric cancer

IMPORTANCE

Epstein-Barr virus (EBV) latency is controlled by epigenetic silencing by DNA methylation [5-methyl cytosine (5mC)], histone modifications, and chromatin looping. However, how they dictate the transcriptional program in EBV-associated gastric cancers remains incompletely understood. EBV-associated gastric cancer displays a 5mC hypermethylated phenotype. A potential treatment for this cancer subtype is the DNA hypomethylating agent, which induces EBV lytic reactivation and targets hypermethylation of the cellular DNA. In this study, we identified a heterogeneous pool of EBV epialleles within two tumor-derived gastric cancer cell lines that are disrupted with a hypomethylating agent. Stochastic DNA methylation patterning at critical regulatory regions may be an underlying mechanism for spontaneous reactivation. Our results highlight the critical role of epigenetic modulation on EBV latency and life cycle, which is maintained through the interaction between 5mC and the host protein CCCTC-binding factor.

Abdominal surgery induces long-lasting changes in expression and binding of CTCF with impact on Major Histocompatibility Complex II transcription in circulating human monocytes

BACKGROUND

Postoperative immunosuppression has been recognized as an important driver of surgery-related morbidity and mortality. It is characterized by lymphocyte depression and impaired monocyte capability to present foreign antigens to T-cells via Major Histocompatibility Complex, Class II (MHC-II) molecules. In patients with postoperative abdominal sepsis, we previously detected a persisting differential binding of the CCCTC-Binding Factor (CTCF), a superordinate regulator of transcription, inside the MHC-II region with specific impact on human leucocyte antigen (HLA) gene expression. In this prospective exploratory study, we investigated to which extent major surgery affects the MHC-II region of circulating CD14+-monocytes.

RESULTS

In non-immunocompromised patients undergoing elective major abdominal surgery, a postoperative loss of monocyte HLA-DR surface receptor density was accompanied by a decline in the transcription levels of the classical MHC-II genes HLA-DRA, HLA-DRB1, HLA-DPA1 and HLA-DPB1. The surgical event decreased the expression of the transcriptional MHC-II regulators CIITA and CTCF and led to a lower CTCF enrichment at an intergenic sequence within the HLA-DR subregion. During the observation period, we found a slow and only incomplete restoration of monocyte HLA-DR surface receptor density as well as a partial recovery of CIITA, HLA-DRA and HLA-DRB1 expression. In contrast, transcription of HLA-DPA1, HLA-DPB1, CTCF and binding of CTCF within the MHC-II remained altered.

CONCLUSION

In circulating monocytes, major surgery does not globally affect MHC-II transcription but rather induces specific changes in the expression of selected HLA genes, followed by differential recovery patterns and accompanied by a prolonged reduction of CTCF expression and binding within the MHC-II region. Our results hint toward a long-lasting impact of a major surgical intervention on monocyte functionality, possibly mediated by epigenetic changes that endure the life span of the individual cell.

Regulation of CTCF loop formation during pancreatic cell differentiation

Transcription reprogramming during cell differentiation involves targeting enhancers to genes responsible for establishment of cell fates. To understand the contribution of CTCF-mediated chromatin organization to cell lineage commitment, we analyzed 3D chromatin architecture during the differentiation of human embryonic stem cells into pancreatic islet organoids. We find that CTCF loops are formed and disassembled at different stages of the differentiation process by either recruitment of CTCF to new anchor sites or use of pre-existing sites not previously involved in loop formation. Recruitment of CTCF to new sites in the genome involves demethylation of H3K9me3 to H3K9me2, demethylation of DNA, recruitment of pioneer factors, and positioning of nucleosomes flanking the new CTCF sites. Existing CTCF sites not involved in loop formation become functional loop anchors via the establishment of new cohesin loading sites containing NIPBL and YY1 at sites between the new anchors. In both cases, formation of new CTCF loops leads to strengthening of enhancer promoter interactions and increased transcription of genes adjacent to loop anchors. These results suggest an important role for CTCF and cohesin in controlling gene expression during cell differentiation.

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.

Structures of CTCF-DNA complexes including all 11 zinc fingers

The CCCTC-binding factor (CTCF) binds tens of thousands of enhancers and promoters on mammalian chromosomes by means of its 11 tandem zinc finger (ZF) DNA-binding domain. In addition to the 12-15-bp CORE sequence, some of the CTCF binding sites contain 5' upstream and/or 3' downstream motifs. Here, we describe two structures for overlapping portions of human CTCF, respectively, including ZF1-ZF7 and ZF3-ZF11 in complex with DNA that incorporates the CORE sequence together with either 3' downstream or 5' upstream motifs. Like conventional tandem ZF array proteins, ZF1-ZF7 follow the right-handed twist of the DNA, with each finger occupying and recognizing one triplet of three base pairs in the DNA major groove. ZF8 plays a unique role, acting as a spacer across the DNA minor groove and positioning ZF9-ZF11 to make cross-strand contacts with DNA. We ascribe the difference between the two subgroups of ZF1-ZF7 and ZF8-ZF11 to residues at the two positions -6 and -5 within each finger, with small residues for ZF1-ZF7 and bulkier and polar/charged residues for ZF8-ZF11. ZF8 is also uniquely rich in basic amino acids, which allows salt bridges to DNA phosphates in the minor groove. Highly specific arginine-guanine and glutamine-adenine interactions, used to recognize G:C or A:T base pairs at conventional base-interacting positions of ZFs, also apply to the cross-strand interactions adopted by ZF9-ZF11. The differences between ZF1-ZF7 and ZF8-ZF11 can be rationalized structurally and may contribute to recognition of high-affinity CTCF binding sites.

The crucial role of CTCF in mitotic progression during early development of sea urchin

CCCTC-binding factor (CTCF), an insulator protein with 11 zinc fingers, is enriched at the boundaries of topologically associated domains (TADs) in eukaryotic genomes. In this study, we isolated and analyzed the cDNAs encoding HpCTCF, the CTCF homolog in the sea urchin Hemicentrotus pulcherrimus, to investigate its expression patterns and functions during the early development of sea urchin. HpCTCF contains nine zinc fingers corresponding to fingers 2-10 of the vertebrate CTCF. Expression pattern analysis revealed that HpCTCF mRNA was detected at all developmental stages and in the entire embryo. Upon expressing the HpCTCF-GFP fusion protein in early embryos, we observed its uniform distribution within interphase nuclei. However, during mitosis, it disappeared from the chromosomes and subsequently reassembled on the chromosome during telophase. Moreover, the morpholino-mediated knockdown of HpCTCF resulted in mitotic arrest during the morula to blastula stage. Most of the arrested chromosomes were not phospholylated at serine 10 of histone H3, indicating that mitosis was arrested at the telophase by HpCTCF depletion. Furthermore, impaired sister chromatid segregation was observed using time-lapse imaging of HpCTCF-knockdown embryos. Thus, HpCTCF is essential for mitotic progression during the early development of sea urchins, especially during the telophase-to-interphase transition. However, the normal development of pluteus larvae in CRISPR-mediated HpCTCF-knockout embryos suggests that disruption of zygotic HpCTCF expression has little effect on embryonic and larval development.

CTCF and R-loops are boundaries of cohesin-mediated DNA looping

Cohesin and CCCTC-binding factor (CTCF) are key regulatory proteins of three-dimensional (3D) genome organization. Cohesin extrudes DNA loops that are anchored by CTCF in a polar orientation. Here, we present direct evidence that CTCF binding polarity controls cohesin-mediated DNA looping. Using single-molecule imaging, we demonstrate that a critical N-terminal motif of CTCF blocks cohesin translocation and DNA looping. The cryo-EM structure of the cohesin-CTCF complex reveals that this CTCF motif ahead of zinc fingers can only reach its binding site on the STAG1 cohesin subunit when the N terminus of CTCF faces cohesin. Remarkably, a C-terminally oriented CTCF accelerates DNA compaction by cohesin. DNA-bound Cas9 and Cas12a ribonucleoproteins are also polar cohesin barriers, indicating that stalling may be intrinsic to cohesin itself. Finally, we show that RNA-DNA hybrids (R-loops) block cohesin-mediated DNA compaction in vitro and are enriched with cohesin subunits in vivo, likely forming TAD boundaries.

Doxorubicin Changes the Spatial Organization of the Genome around Active Promoters

In this study, we delve into the impact of genotoxic anticancer drug treatment on the chromatin structure of human cells, with a particular focus on the effects of doxorubicin. Using Hi-C, ChIP-seq, and RNA-seq, we explore the changes in chromatin architecture brought about by doxorubicin and ICRF193. Our results indicate that physiologically relevant doses of doxorubicin lead to a local reduction in Hi-C interactions in certain genomic regions that contain active promoters, with changes in chromatin architecture occurring independently of Top2 inhibition, cell cycle arrest, and differential gene expression. Inside the regions with decreased interactions, we detected redistribution of RAD21 around the peaks of H3K27 acetylation. Our study also revealed a common structural pattern in the regions with altered architecture, characterized by two large domains separated from each other. Additionally, doxorubicin was found to increase CTCF binding in H3K27 acetylated regions. Furthermore, we discovered that Top2-dependent chemotherapy causes changes in the distance decay of Hi-C contacts, which are driven by direct and indirect inhibitors. Our proposed model suggests that doxorubicin-induced DSBs cause cohesin redistribution, which leads to increased insulation on actively transcribed TAD boundaries. Our findings underscore the significant impact of genotoxic anticancer treatment on the chromatin structure of the human genome.

Chromatin loop dynamics during cellular differentiation are associated with changes to both anchor and internal regulatory features

Three-dimensional (3D) chromatin structure has been shown to play a role in regulating gene transcription during biological transitions. Although our understanding of loop formation and maintenance is rapidly improving, much less is known about the mechanisms driving changes in looping and the impact of differential looping on gene transcription. One limitation has been a lack of well-powered differential looping data sets. To address this, we conducted a deeply sequenced Hi-C time course of megakaryocyte development comprising four biological replicates and 6 billion reads per time point. Statistical analysis revealed 1503 differential loops. Gained loop anchors were enriched for AP-1 occupancy and were characterized by large increases in histone H3K27ac (over 11-fold) but relatively small increases in CTCF and RAD21 binding (1.26- and 1.23-fold, respectively). Linear modeling revealed that changes in histone H3K27ac, chromatin accessibility, and JUN binding were better correlated with changes in looping than RAD21 and almost as well correlated as CTCF. Changes to epigenetic features between-rather than at-boundaries were highly predictive of changes in looping. Together these data suggest that although CTCF and RAD21 may be the core machinery dictating where loops form, other features (both at the anchors and within the loop boundaries) may play a larger role than previously anticipated in determining the relative loop strength across cell types and conditions.

Epigenomic analysis of Alzheimer's disease brains reveals diminished CTCF binding on genes involved in synaptic organization

Epigenetic aberrations are suggested to play an important role in transcriptional alterations in Alzheimer's disease (AD). One of the key mechanisms of epigenetic regulation of gene expression is through the dynamic organization of chromatin structure via the master genome architecture protein, CCCTC-binding factor (CTCF). By forming chromatin loops, CTCF can influence gene transcription in a complex manner. To find out whether genome-wide DNA binding sites for CTCF are altered in AD, we compared CTCF chromatin immunoprecipitation sequencing (ChIP-Seq) data from frontal cortex of human AD patients and normal controls (n = 9 pairs, all females). We have revealed that CTCF-binding affinity on many genes is significantly reduced in AD patients, and these genes are enriched in synaptic organization, cell adhesion, and actin cytoskeleton, including synaptic scaffolding molecules and receptors, such as SHANK2, HOMER1, NRXN1, CNTNAP2 and GRIN2A, and protocadherin (PCDH) and cadherin (CDH) family members. By comparing transcriptomic data from AD patients, we have discovered that many of the synaptic and adhesion genes with reduced CTCF binding in AD are significantly reduced in their mRNA expression. Moreover, a significant overlap of genes with the diminished CTCF binding and the reduced H3K27ac is identified in AD, with the common genes enriched in synaptic organization. These data suggest that the CTCF-controlled 3D chromatin organization is perturbed in AD, which may be linked to the diminished expression of target genes, probably through changes in histone modification.

Multiomic analysis of cohesin reveals that ZBTB transcription factors contribute to chromatin interactions

One bottleneck in understanding the principles of 3D chromatin structures is caused by the paucity of known regulators. Cohesin is essential for 3D chromatin organization, and its interacting partners are candidate regulators. Here, we performed proteomic profiling of the cohesin in chromatin and identified transcription factors, RNA-binding proteins and chromatin regulators associated with cohesin. Acute protein degradation followed by time-series genomic binding quantitation and BAT Hi-C analysis were conducted, and the results showed that the transcription factor ZBTB21 contributes to cohesin chromatin binding, 3D chromatin interactions and transcriptional repression. Strikingly, multiomic analyses revealed that the other four ZBTB factors interacted with cohesin, and double degradation of ZBTB21 and ZBTB7B led to a further decrease in cohesin chromatin occupancy. We propose that multiple ZBTB transcription factors orchestrate the chromatin binding of cohesin to regulate chromatin interactions, and we provide a catalog of many additional proteins associated with cohesin that warrant further investigation.

Sequential and directional insulation by conserved CTCF sites underlies the Hox timer in stembryos

During development, Hox genes are temporally activated according to their relative positions on their clusters, contributing to the proper identities of structures along the rostrocaudal axis. To understand the mechanism underlying this Hox timer, we used mouse embryonic stem cell-derived stembryos. Following Wnt signaling, the process involves transcriptional initiation at the anterior part of the cluster and a concomitant loading of cohesin complexes enriched on the transcribed DNA segments, that is, with an asymmetric distribution favoring the anterior part of the cluster. Chromatin extrusion then occurs with successively more posterior CTCF sites acting as transient insulators, thus generating a progressive time delay in the activation of more posterior-located genes due to long-range contacts with a flanking topologically associating domain. Mutant stembryos support this model and reveal that the presence of evolutionary conserved and regularly spaced intergenic CTCF sites controls the precision and the pace of this temporal mechanism.

Chromatin alternates between A and B compartments at kilobase scale for subgenic organization

Nuclear compartments are prominent features of 3D chromatin organization, but sequencing depth limitations have impeded investigation at ultra fine-scale. CTCF loops are generally studied at a finer scale, but the impact of looping on proximal interactions remains enigmatic. Here, we critically examine nuclear compartments and CTCF loop-proximal interactions using a combination of in situ Hi-C at unparalleled depth, algorithm development, and biophysical modeling. Producing a large Hi-C map with 33 billion contacts in conjunction with an algorithm for performing principal component analysis on sparse, super massive matrices (POSSUMM), we resolve compartments to 500 bp. Our results demonstrate that essentially all active promoters and distal enhancers localize in the A compartment, even when flanking sequences do not. Furthermore, we find that the TSS and TTS of paused genes are often segregated into separate compartments. We then identify diffuse interactions that radiate from CTCF loop anchors, which correlate with strong enhancer-promoter interactions and proximal transcription. We also find that these diffuse interactions depend on CTCF's RNA binding domains. In this work, we demonstrate features of fine-scale chromatin organization consistent with a revised model in which compartments are more precise than commonly thought while CTCF loops are more protracted.

Structural elements promote architectural stripe formation and facilitate ultra-long-range gene regulation at a human disease locus

Enhancer clusters overlapping disease-associated mutations in Pierre Robin sequence (PRS) patients regulate SOX9 expression at genomic distances over 1.25 Mb. We applied optical reconstruction of chromatin architecture (ORCA) imaging to trace 3D locus topology during PRS-enhancer activation. We observed pronounced changes in locus topology between cell types. Subsequent analysis of single-chromatin fiber traces revealed that these ensemble-average differences arise through changes in the frequency of commonly sampled topologies. We further identified two CTCF-bound elements, internal to the SOX9 topologically associating domain, which promote stripe formation, are positioned near the domain's 3D geometric center, and bridge enhancer-promoter contacts in a series of chromatin loops. Ablation of these elements results in diminished SOX9 expression and altered domain-wide contacts. Polymer models with uniform loading across the domain and frequent cohesin collisions recapitulate this multi-loop, centrally clustered geometry. Together, we provide mechanistic insights into architectural stripe formation and gene regulation over ultra-long genomic ranges.

CTCF variant begets to short stature by down-regulation of IGF1

Pathogenic variants in the transcription factor CCCTC-binding factor (CTCF) are associated with mental retardation, autosomal dominant 21 (MRD21, MIM#615502). Current studies supported the strong relationship between CTCF variants and growth, yet the mechanism of CTCF mutation leading to short stature is not known. Clinical information, treatment regimens, and follow-up outcomes of a patient with MRD21 were collected. The possible pathogenic mechanisms of CTCF variants leading to short stature were investigated using immortalized lymphocyte cell lines (LCLs), HEK-293T, and immortalized normal human liver cell lines (LO2). This patient received long-term treatment with recombinant human growth hormone (rhGH) which resulted in an increased height of 1.0 SDS. She had low serum insulin-like growth factor 1 (IGF1) before the treatment and the IGF1 level was not significantly increased during the treatment (-1.38 ± 0.61 SDS). The finding suggested that the CTCF R567W variant could have impaired IGF1 production pathway. We further demonstrated that the mutant CTCF had a reduced ability to bind to the promoter region of IGF1, consequently significantly reducing the transcriptional activation and expression of IGF1. Our novel results demonstrated a direct positive regulation of CTCF on the transcription of the IGF1 promoter. The impaired IGF1 expression due to CTCF mutation may explain the substandard effect of rhGH treatment on MRD21 patients. This study provided novel insights into the molecular basis of CTCF-associated disorder.

Enhancer-promoter contact formation requires RNAPII and antagonizes loop extrusion

Homotypic chromatin interactions and loop extrusion are thought to be the two main drivers of mammalian chromosome folding. Here we tested the role of RNA polymerase II (RNAPII) across different scales of interphase chromatin organization in a cellular system allowing for its rapid, auxin-mediated degradation. We combined Micro-C and computational modeling to characterize subsets of loops differentially gained or lost upon RNAPII depletion. Gained loops, extrusion of which was antagonized by RNAPII, almost invariably formed by engaging new or rewired CTCF anchors. Lost loops selectively affected contacts between enhancers and promoters anchored by RNAPII, explaining the repression of most genes. Surprisingly, promoter-promoter interactions remained essentially unaffected by polymerase depletion, and cohesin occupancy was sustained. Together, our findings reconcile the role of RNAPII in transcription with its direct involvement in setting-up regulatory three-dimensional chromatin contacts genome wide, while also revealing an impact on cohesin loop extrusion.

Elevated enhancer-oncogene contacts and higher oncogene expression levels by recurrent CTCF inactivating mutations in acute T cell leukemia

Monoallelic inactivation of CCCTC-binding factor (CTCF) in human cancer drives altered methylated genomic states, altered CTCF occupancy at promoter and enhancer regions, and deregulated global gene expression. In patients with T cell acute lymphoblastic leukemia (T-ALL), we find that acquired monoallelic CTCF-inactivating events drive subtle and local genomic effects in nearly half of t(5; 14) (q35; q32.2) rearranged patients, especially when CTCF-binding sites are preserved in between the BCL11B enhancer and the TLX3 oncogene. These solitary intervening sites insulate TLX3 from the enhancer by inducing competitive looping to multiple binding sites near the TLX3 promoter. Reduced CTCF levels or deletion of the intervening CTCF site abrogates enhancer insulation by weakening competitive looping while favoring TLX3 promoter to BCL11B enhancer looping, which elevates oncogene expression levels and leukemia burden.

CTCF mediates CD8+ effector differentiation through dynamic redistribution and genomic reorganization

Differentiation of effector CD8+ T cells is instructed by stably and dynamically expressed transcription regulators. Here we show that naive-to-effector differentiation was accompanied by dynamic CTCF redistribution and extensive chromatin architectural changes. Upon CD8+ T cell activation, CTCF acquired de novo binding sites and anchored novel chromatin interactions, and these changes were associated with increased chromatin accessibility and elevated expression of cytotoxic program genes including Tbx21, Ifng, and Klrg1. CTCF was also evicted from its ex-binding sites in naive state, with concomitantly reduced chromatin interactions in effector cells, as observed at memory precursor-associated genes including Il7r, Sell, and Tcf7. Genetic ablation of CTCF indeed diminished cytotoxic gene expression, but paradoxically elevated expression of memory precursor genes. Comparative Hi-C analysis revealed that key memory precursor genes were harbored within insulated neighborhoods demarcated by constitutive CTCF binding, and their induction was likely due to disrupted CTCF-dependent insulation. CTCF thus promotes cytotoxic effector differentiation by integrating local chromatin accessibility control and higher-order genomic reorganization.

Associations of CTCF and FOXA1 with androgen and IGF pathways in men with localized prostate cancer

AIMS

To examine associations between the transcription factors CCCTC-binding factor (CTCF) and forkhead box protein A1 (FOXA1) and the androgen receptor (AR) and their association with components of the insulin-like growth factor (IGF)-pathway in a cohort of men with localized prostate cancer.

METHODS

Using prostate tissue samples collected during the Prostate cancer: Evidence of Exercise and Nutrition Trial (PrEvENT) trial (N = 70 to 92, depending on section availability), we assessed the abundance of CTCF, FOXA1, AR, IGFIR, p-mTOR, PTEN and IGFBP-2 proteins using a modified version of the Allred scoring system. Validation studies were performed using large, publicly available datasets (TCGA) (N = 489).

RESULTS

We identified a strong correlation between CTCF and AR staining with benign prostate tissue. CTCF also strongly associated with the IGFIR, with PTEN and with phospho-mTOR. FOXA1 was also correlated with staining for the IGF-IR, with IGFBP-2 and with staining for activated phosphor-mTOR. The staining for the IGF-IR was strongly correlated with the AR.

CONCLUSION

Our findings emphasise the close and complex links between the endocrine controls, well known to play an important role in prostate cancer, and the transcription factors implicated by the recent genetic evidence.

Transcription shapes 3D chromatin organization by interacting with loop extrusion

Cohesin folds mammalian interphase chromosomes by extruding the chromatin fiber into numerous loops. "Loop extrusion" can be impeded by chromatin-bound factors, such as CTCF, which generates characteristic and functional chromatin organization patterns. It has been proposed that transcription relocalizes or interferes with cohesin and that active promoters are cohesin loading sites. However, the effects of transcription on cohesin have not been reconciled with observations of active extrusion by cohesin. To determine how transcription modulates extrusion, we studied mouse cells in which we could alter cohesin abundance, dynamics, and localization by genetic "knockouts" of the cohesin regulators CTCF and Wapl. Through Hi-C experiments, we discovered intricate, cohesin-dependent contact patterns near active genes. Chromatin organization around active genes exhibited hallmarks of interactions between transcribing RNA polymerases (RNAPs) and extruding cohesins. These observations could be reproduced by polymer simulations in which RNAPs were moving barriers to extrusion that obstructed, slowed, and pushed cohesins. The simulations predicted that preferential loading of cohesin at promoters is inconsistent with our experimental data. Additional ChIP-seq experiments showed that the putative cohesin loader Nipbl is not predominantly enriched at promoters. Therefore, we propose that cohesin is not preferentially loaded at promoters and that the barrier function of RNAP accounts for cohesin accumulation at active promoters. Altogether, we find that RNAP is an extrusion barrier that is not stationary, but rather, translocates and relocalizes cohesin. Loop extrusion and transcription might interact to dynamically generate and maintain gene interactions with regulatory elements and shape functional genomic organization.

Deciphering Hierarchical Chromatin Domains and Preference of Genomic Position Forming Boundaries in Single Mouse Embryonic Stem Cells

The exploration of single-cell 3D genome maps reveals that chromatin domains are indeed physical structures presenting in single cells, and domain boundaries vary from cell to cell. However, systematic analysis of the association between regulatory factor binding and elements and the formation of chromatin domains in single cells has not yet emerged. To this end, a hierarchical chromatin domain structure identification algorithm (named as HiCS) is first developed from individual single-cell Hi-C maps, with superior performance in both accuracy and efficiency. The results suggest that in addition to the known CTCF-cohesin complex, Polycomb, TrxG, pluripotent protein families, and other multiple factors also contribute to shaping chromatin domain boundaries in single embryonic stem cells. Different cooperation patterns of these regulatory factors drive genomic position categories with differential preferences forming boundaries, and the most extensive six types of retrotransposons are differentially distributed in these genomic position categories with preferential localization. The above results suggest that these different retrotransposons within genomic regions interplay with regulatory factors navigating the preference of genomic positions forming boundaries, driving the formation of higher-order chromatin structures, and thus regulating cell functions in single mouse embryonic stem cells.

CTCF, BEAF-32, and CP190 are not required for the establishment of TADs in early Drosophila embryos but have locus-specific roles

The boundaries of topologically associating domains (TADs) are delimited by insulators and/or active promoters; however, how they are initially established during embryogenesis remains unclear. Here, we examined this during the first hours of Drosophila embryogenesis. DNA-FISH confirms that intra-TAD pairwise proximity is established during zygotic genome activation (ZGA) but with extensive cell-to-cell heterogeneity. Most newly formed boundaries are occupied by combinations of CTCF, BEAF-32, and/or CP190. Depleting each insulator individually from chromatin revealed that TADs can still establish, although with lower insulation, with a subset of boundaries (~10%) being more dependent on specific insulators. Some weakened boundaries have aberrant gene expression due to unconstrained enhancer activity. However, the majority of misexpressed genes have no obvious direct relationship to changes in domain-boundary insulation. Deletion of an active promoter (thereby blocking transcription) at one boundary had a greater impact than deleting the insulator-bound region itself. This suggests that cross-talk between insulators and active promoters and/or transcription might reinforce domain boundary insulation during embryogenesis.

Topological screen identifies hundreds of Cp190- and CTCF-dependent Drosophila chromatin insulator elements

Drosophila insulators were the first DNA elements found to regulate gene expression by delimiting chromatin contacts. We still do not know how many of them exist and what impact they have on the Drosophila genome folding. Contrary to vertebrates, there is no evidence that fly insulators block cohesin-mediated chromatin loop extrusion. Therefore, their mechanism of action remains uncertain. To bridge these gaps, we mapped chromatin contacts in Drosophila cells lacking the key insulator proteins CTCF and Cp190. With this approach, we found hundreds of insulator elements. Their study indicates that Drosophila insulators play a minor role in the overall genome folding but affect chromatin contacts locally at many loci. Our observations argue that Cp190 promotes cobinding of other insulator proteins and that the model, where Drosophila insulators block chromatin contacts by forming loops, needs revision. Our insulator catalog provides an important resource to study mechanisms of genome folding.

Enhancer-promoter interactions can bypass CTCF-mediated boundaries and contribute to phenotypic robustness

How enhancers activate their distal target promoters remains incompletely understood. Here we dissect how CTCF-mediated loops facilitate and restrict such regulatory interactions. Using an allelic series of mouse mutants, we show that CTCF is neither required for the interaction of the Sox2 gene with distal enhancers, nor for its expression. Insertion of various combinations of CTCF motifs, between Sox2 and its distal enhancers, generated boundaries with varying degrees of insulation that directly correlated with reduced transcriptional output. However, in both epiblast and neural tissues, enhancer contacts and transcriptional induction could not be fully abolished, and insertions failed to disrupt implantation and neurogenesis. In contrast, Sox2 expression was undetectable in the anterior foregut of mutants carrying the strongest boundaries, and these animals fully phenocopied loss of SOX2 in this tissue. We propose that enhancer clusters with a high density of regulatory activity can better overcome physical barriers to maintain faithful gene expression and phenotypic robustness.

Characterization of a new CCCTC-binding factor binding site as a dual regulator of Epstein-Barr virus latent infection

Distinct viral gene expression characterizes Epstein-Barr virus (EBV) infection in EBV-producing marmoset B-cell (B95-8) and EBV-associated gastric carcinoma (SNU719) cell lines. CCCTC-binding factor (CTCF) is a structural chromatin factor that coordinates chromatin interactions in the EBV genome. Chromatin immunoprecipitation followed by sequencing against CTCF revealed 16 CTCF binding sites in the B95-8 and SNU719 EBV genomes. The biological function of one CTCF binding site (S13 locus) located on the BamHI A right transcript (BART) miRNA promoter was elucidated experimentally. Microscale thermophoresis assay showed that CTCF binds more readily to the stable form than the mutant form of the S13 locus. EBV BART miRNA clusters encode 22 miRNAs, whose roles are implicated in EBV-related cancer pathogenesis. The B95-8 EBV genome lacks a 11.8-kb EcoRI C fragment, whereas the SNU719 EBV genome is full-length. ChIP-PCR assay revealed that CTCF, RNA polymerase II, H3K4me3 histone, and H3K9me3 histone were more enriched at S13 and S16 (167-kb) loci in B95-8 than in the SNU719 EBV genome. 4C-Seq and 3C-PCR assays using B95-8 and SNU719 cells showed that the S13 locus was associated with overall EBV genomic loci including 3-kb and 167-kb region in both EBV genomes. We generated mutations in the S13 locus in bacmids with or without the 11.8-kb BART transcript unit (BART(+/-)). The S13 mutation upregulated BART miRNA expression, weakened EBV latency, and reduced EBV infectivity in the presence of EcoRI C fragment. Another 3C-PCR assay using four types of BART(+/-)·S13(wild-type(Wt)/mutant(Mt)) HEK293-EBV cells revealed that the S13 mutation decreased DNA associations between the 167-kb region and 3-kb in the EBV genome. Based on these results, CTCF bound to the S13 locus along with the 11.8-kb EcoRI C fragment is suggested to form an EBV 3-dimensional DNA loop for coordinated EBV BART miRNA expression and infectivity.

Cohesin and CTCF complexes mediate contacts in chromatin loops depending on nucleosome positions

The spatial organization of the eukaryotic genome plays an important role in regulating transcriptional activity. In the nucleus, chromatin forms loops that assemble into fundamental units called topologically associating domains that facilitate or inhibit long-range contacts. These loops are formed and held together by the ring-shaped cohesin protein complex, and this can involve binding of CCCTC-binding factor (CTCF). High-resolution conformation capture experiments provide the frequency at which two DNA fragments physically associate in three-dimensional space. However, technical limitations of this approach, such as low throughput, low resolution, or noise in contact maps, make data interpretation and identification of chromatin intraloop contacts, e.g., between distal regulatory elements and their target genes, challenging. Herein, an existing coarse-grained model of chromatin at single-nucleosome resolution was extended by integrating potentials describing CTCF and cohesin. We performed replica-exchange Monte Carlo simulations with regularly spaced nucleosomes and experimentally determined nucleosome positions in the presence of cohesin-CTCF, as well as depleted systems as controls. In fully extruded loops caused by the presence of cohesin and CTCF, the number of contacts within the formed loops was increased. The number and types of these contacts were impacted by the nucleosome distribution and loop size. Microloops were observed within cohesin-mediated loops due to thermal fluctuations without additional influence of other factors, and the number, size, and shape of microloops were determined by nucleosome distribution and loop size. Nucleosome positions directly affect the spatial structure and contact probability within a loop, with presumed consequences for transcriptional activity.

The Interplay of Transcription and Genome Topology Programs T Cell Development and Differentiation

T cells are essential for mounting defense against various pathogens and malignantly transformed cells. Thymic development and peripheral T cell differentiation are highly orchestrated biological processes that require precise gene regulation. Higher-order genome organization on multiple scales, in the form of chromatin loops, topologically associating domains and compartments, provides pivotal control of T cell gene expression. CTCF and the cohesin machinery are ubiquitously expressed architectural proteins responsible for establishing chromatin structures. Recent studies indicate that transcription factors, such as T lineage-defining Tcf1 and TCR-induced Batf, may have intrinsic ability and/or engage CTCF to shape chromatin architecture. In this article, we summarize current knowledge on the dynamic changes in genome topology that underlie normal or leukemic T cell development, CD4+ helper T cell differentiation, and CD8+ cytotoxic T cell functions. The knowledge lays a solid foundation for elucidating the causative link of spatial chromatin configuration to transcriptional and functional output in T cells.

A cohesin traffic pattern genetically linked to gene regulation

Cohesin-mediated loop extrusion has been shown to be blocked at specific cis-elements, including CTCF sites, producing patterns of loops and domain boundaries along chromosomes. Here we explore such cis-elements, and their role in gene regulation. We find that transcription termination sites of active genes form cohesin- and RNA polymerase II-dependent domain boundaries that do not accumulate cohesin. At these sites, cohesin is first stalled and then rapidly unloaded. Start sites of transcriptionally active genes form cohesin-bound boundaries, as shown before, but are cohesin-independent. Together with cohesin loading, possibly at enhancers, these sites create a pattern of cohesin traffic that guides enhancer-promoter interactions. Disrupting this traffic pattern, by removing CTCF, renders cells sensitive to knockout of genes involved in transcription initiation, such as the SAGA complexes, and RNA processing such DEAD/H-Box RNA helicases. Without CTCF, these factors are less efficiently recruited to active promoters.

Co-depletion of NIPBL and WAPL balance cohesin activity to correct gene misexpression

The relationship between cohesin-mediated chromatin looping and gene expression remains unclear. NIPBL and WAPL are two opposing regulators of cohesin activity; depletion of either is associated with changes in both chromatin folding and transcription across a wide range of cell types. However, a direct comparison of their individual and combined effects on gene expression in the same cell type is lacking. We find that NIPBL or WAPL depletion in human HCT116 cells each alter the expression of ~2,000 genes, with only ~30% of the genes shared between the conditions. We find that clusters of differentially expressed genes within the same topologically associated domain (TAD) show coordinated misexpression, suggesting some genomic domains are especially sensitive to both more or less cohesin. Finally, co-depletion of NIPBL and WAPL restores the majority of gene misexpression as compared to either knockdown alone. A similar set of NIPBL-sensitive genes are rescued following CTCF co-depletion. Together, this indicates that altered transcription due to reduced cohesin activity can be functionally offset by removal of either its negative regulator (WAPL) or the physical barriers (CTCF) that restrict loop-extrusion events.

Co-depletion of NIPBL and WAPL balance cohesin activity to correct gene misexpression

The relationship between cohesin-mediated chromatin looping and gene expression remains unclear. NIPBL and WAPL are two opposing regulators of cohesin activity; depletion of either is associated with changes in both chromatin folding and transcription across a wide range of cell types. However, a direct comparison of their individual and combined effects on gene expression in the same cell type is lacking. We find that NIPBL or WAPL depletion in human HCT116 cells each alter the expression of ~2,000 genes, with only ~30% of the genes shared between the conditions. We find that clusters of differentially expressed genes within the same topologically associated domain (TAD) show coordinated misexpression, suggesting some genomic domains are especially sensitive to both more or less cohesin. Finally, co-depletion of NIPBL and WAPL restores the majority of gene misexpression as compared to either knockdown alone. A similar set of NIPBL-sensitive genes are rescued following CTCF co-depletion. Together, this indicates that altered transcription due to reduced cohesin activity can be functionally offset by removal of either its negative regulator (WAPL) or the physical barriers (CTCF) that restrict loop-extrusion events.

Chromatin jets define the properties of cohesin-driven in vivo loop extrusion

Complex genomes show intricate organization in three-dimensional (3D) nuclear space. Current models posit that cohesin extrudes loops to form self-interacting domains delimited by the DNA binding protein CTCF. Here, we describe and quantitatively characterize cohesin-propelled, jet-like chromatin contacts as landmarks of loop extrusion in quiescent mammalian lymphocytes. Experimental observations and polymer simulations indicate that narrow origins of loop extrusion favor jet formation. Unless constrained by CTCF, jets propagate symmetrically for 1-2 Mb, providing an estimate for the range of in vivo loop extrusion. Asymmetric CTCF binding deflects the angle of jet propagation as experimental evidence that cohesin-mediated loop extrusion can switch from bi- to unidirectional and is controlled independently in both directions. These data offer new insights into the physiological behavior of in vivo cohesin-mediated loop extrusion and further our understanding of the principles that underlie genome organization.