Search

Clustered regularly interspaced short palindromic repeats (CRISPR) technology emerges a remarkable potential for cure of refractory cancer like metastatic breast cancer. However, how to efficiently deliver the CRISPR system with non-viral carrier remains a major issue to be solved. Here, we report a kind of targeted core-shell nanoparticles (NPs) carrying dual plasmids (pHR-pCas9) for precise CCCTC-binding factor (CTCF) gene insert to circumvent metastatic breast cancer. The targeted core-shell NPs carrying pHR-pCas9 can accomplish γGTP-mediated cellular uptake and endosomal escape, facilitate the precise insert and stable expression of CTCF gene, inhibit the migration, metastasis, and colonization of metastatic breast cancer cells. Besides, the finding further reveals that the inhibitory mechanism of metastasis could be associated with up-regulating CTCF protein, followed by down-regulating stomatin (STOM) protein. The study offers a universal nanostrategy enabling the robust non-viral delivery of gene-editing system for treatment of severe illness.

The three-dimensional structure of the human genome has been proven to have a significant functional impact on gene expression. The high-order spatial chromatin is organised first by looping mediated by multiple protein factors, and then it is further formed into larger structures of topologically associated domains (TADs) or chromatin contact domains (CCDs), followed by A/B compartments and finally the chromosomal territories (CTs). The genetic variation observed in human population influences the multi-scale structures, posing a question regarding the functional impact of structural variants reflected by the variability of the genes expression patterns. The current methods of evaluating the functional effect include eQTLs analysis which uses statistical testing of influence of variants on spatially close genes. Rarely, non-coding DNA sequence changes are evaluated by their impact on the biomolecular interaction network (BIN) reflecting the cellular interactome that can be analysed by the classical graph-theoretic algorithms. Therefore, in the second part of the review, we introduce the concept of BIN, i.e. a meta-network model of the complete molecular interactome developed by integrating various biological networks. The BIN meta-network model includes DNA-protein binding by the plethora of protein factors as well as chromatin interactions, therefore allowing connection of genomics with the downstream biomolecular processes present in a cell. As an illustration, we scrutinise the chromatin interactions mediated by the CTCF protein detected in a ChIA-PET experiment in the human lymphoblastoid cell line GM12878. In the corresponding BIN meta-network the DNA spatial proximity is represented as a graph model, combined with the Proteins-Interaction Network (PIN) of human proteome using the Gene Association Network (GAN). Furthermore, we enriched the BIN with the signalling and metabolic pathways and Gene Ontology (GO) terms to assert its functional context. Finally, we mapped the Single Nucleotide Polymorphisms (SNPs) from the GWAS studies and identified the chromatin mutational hot-spots associated with a significant enrichment of SNPs related to autoimmune diseases. Afterwards, we mapped Structural Variants (SVs) from healthy individuals of 1000 Genomes Project and identified an interesting example of the missing protein complex associated with protein Q6GYQ0 due to a deletion on chromosome 14. Such an analysis using the meta-network BIN model is therefore helpful in evaluating the influence of genetic variation on spatial organisation of the genome and its functional effect in a cell.

Insulators characterized in Drosophila and mammals have been shown to play a key role in the restriction of promiscuous enhancer-promoter interactions, as well as reshaping the topological landscape of chromosomes. Yet the role of insulators in plants remains poorly understood, in large part because of a lack of well-characterized insulators and binding factor(s). In this study, we isolated a 1.2-kb RS2-9 insulator from the Oryza sativa (rice) genome that can, when interposed between an enhancer and promoter, efficiently block the activation function of both constitutive and floral organ-specific enhancers in transgenic Arabidopsis and Nicotiana tabacum (tobacco). In the rice genome, the genes flanking RS2-9 exhibit an absence of mutual transcriptional interactions, as well as a lack of histone modification spread. We further determined that O. sativa Homeobox 1 (OSH1) bound two regions of RS2-9, as well as over 50 000 additional sites in the rice genome, the majority of which resided in intergenic regions. Mutation of one of the two OSH1-binding sites in RS2-9 impaired insulation activity by up to 60%, whereas the mutation of both binding sites virtually abolished insulator function. We also demonstrated that OSH1 binding sites were associated with 72% of the boundaries of topologically associated domains (TADs) identified in the rice genome, which is comparable to the 77% of TAD boundaries bound by the insulator CCCTC-binding factor (CTCF) in mammals. Taken together, our findings indicate that OSH1-RS2-9 acts as a true insulator in plants, and highlight a potential role for OSH1 in gene insulation and topological organization in plant genomes.

Objective Genetic variants spanning the ubiquitin-conjugating enzyme E2 L3 (UBE2L3) gene are associated with increased expression of the UBE2L3-encoded E2 ubiquitin-conjugating enzyme, UbcH7, that facilitates activation of proinflammatory NF-κB signaling, and susceptibility to autoimmune diseases. This study aims to delineate how genetic variants carried on the UBE2L3-YDJC autoimmune risk haplotype function to drive hypermorphic UBE2L3 expression. Methods We used bioinformatic analyses, electrophoretic mobility shift assays, and luciferase reporter assays to identify and functionally characterize allele-specific effects of risk variants positioned in chromatin accessible regions of immune cells. Chromatin conformation capture (3C)-qPCR, ChIP-qPCR, and siRNA knockdown assays were performed on patient-derived EBV-transformed B cells homozygous for the UBE2L3-YDJC non-risk or risk haplotype to determine if the risk haplotype increases UBE2L3 expression by altering the regulatory chromatin architecture in the region. Results Five of the seven prioritized variants demonstrated allele-specific increases in nuclear protein binding affinity and regulatory activity. HiChIP and 3C-qPCR uncovered a long-range interaction between the UBE2L3 promoter (rs140490, rs140491, rs11089620) and downstream YDJC promoter (rs3747093) that was strengthened in the presence of the UBE2L3-YDJC risk haplotype, and correlated with the loss of CTCF binding and gain of YY1 binding at the risk alleles. Depleting YY1 by siRNA disrupted the long-range interaction between the two promoters and reduced UBE2L3 expression. Conclusion The UBE2L3-YDJC autoimmune risk haplotype increases UBE2L3 expression through strengthening a YY1-mediated interaction between the UBE2L3 and YDJC promoters.

Cisplatin-based chemotherapy remains the standard care for non-small cell lung cancer (NSCLC) patients. Relapse after chemotherapy-induced dormancy affects the overall survival of patients. The evolution of cancer cells under chemotherapy stress is regulated by transcription factors (TFs) with binding sites initially buried deep within inaccessible chromatin. The transcription machinery and dynamic epigenetic alterations during the process of dormancy-reactivation of lung cancer cells after chemotherapy need to be investigated. Here, we investigated the chromatin accessibility of lung cancer cells after cisplatin treatment, using an assay for transposase-accessible chromatin sequencing (ATAC-seq). We observed that global chromatin accessibility was extensively improved. Transcriptional Regulatory Relationships Unraveled by Sentence-based Text mining (TRRUST) v.2 was used to elucidate TF-target interaction during the process of dormancy and reactivation. Enhancer regions and motifs specific to key TFs including JUN, MYC, SMAD3, E2F1, SP1, CTCF, SMAD4, STAT3, NFKB1, and KLF4 were enriched in differential loci ATAC-seq peaks of dormant and reactivated cancer cells induced by chemotherapy. The findings suggest that these key TFs regulated gene expressions during the process of dormancy and reactivation of cancer cells through altering promoter accessibility of target genes. Our study helps advance understanding of how cancer cells adapt to the stress induced by chemotherapy through TF binding motif accessibility., Graphical abstract, In this study, we profiled chromatin accessibility and transcriptome in cisplatin-responsive dormant and reactivated cells. We observed features of the chromatin landscape of lung cancer cells after chemotherapy. Our integrative analyses revealed that adaptive cellular response is often achieved by transcription factor binding motif accessibility.

Although previous studies identified numerous single nucleotide polymorphisms (SNPs) and their target genes predisposed to prostate cancer (PrCa) risks, SNP-related splicing associations are rarely reported. In this study, we applied distance-based sQTL analysis (sQTLseekeR) using RNA-seq and SNP genotype data from benign prostate tissue (n = 467) and identified significant associations in 3344 SNP-transcript pairs (P ≤ 0.05) at PrCa risk loci. We characterized a common SNP (rs7247241) and its target gene (PPP1R14A) located in chr19q13, an sQTL with risk allele T associated with upregulation of long isoform (P = 9.99E−7). We confirmed the associations in both TCGA (P = 2.42E−24) and GTEX prostate cohorts (P = 9.08E−78). To functionally characterize this SNP, we performed chromatin immunoprecipitation qPCR and confirmed stronger CTCF and PLAGL2 binding in rs7247241 C than T allele. We found that CTCF binding enrichment was negatively associated with methylation level at the SNP site in human cell lines (r = −0.58). Bisulfite sequencing showed consistent association of rs7247241-T allele with nearby sequence CpG hypermethylation in prostate cell lines and tissues. Moreover, the methylation level at CpG sites nearest to the CTCF binding and first exon splice-in (ψ) of PPP1R14A was significantly associated with aggressive phenotype in the TCGA PrCa cohort. Meanwhile, the long isoform of the gene also promoted cell proliferation. Taken together, with the most updated gene annotations, we reported a set of sQTL associated with multiple traits related to human prostate diseases and revealed a unique role of PrCa risk SNP rs7247241 on PPP1R14A isoform transition.

Deregulations of DNA-methylation-related genes are common in cancers, but frameshift mutation status in colon cancer (CC) is unknown. Our study aims to assess whether CTCF, ZFP57, and ATF7IP genes in this category are mutated in CC. CTCF, ZFP57, and ATF7IP genes have repeat coding sequences, which are frequently deleted or duplicated in CC, harboring the phenotype of unstable or high microsatellite instability (MSI-H). We studied 140 CCs [95 MSI-H CCs and 45 stable MSI (MSS) CCs], and found 7 CCs with MSI-H (6/95: 6.3%) harbored frameshift mutations within the repeats, whereas those with MSS did not. Of note, the CTCF frameshift mutations showed the regional difference in the 2 (12.5%) of 16 MSI-H CCs, indicating there was intratumoral heterogeneity. In the immunohistochemistry for ATF7IP, the MSI-H CC showed low intensity compared to MSS CC. Together, CTCF, ZFP57, and ATF7IP genes, despite the low incidence of the mutations, are altered in several ways (mutation, expression, and intratumoral heterogeneity) and could contribute to MSI-H CC development.

Transcription is regulated and mediated by multiprotein complexes in a chromatin context. Transcription causes changes in DNA topology which is modulated by DNA topoisomerases, enzymes that catalyse changes in DNA topology via transient breaking and re-joining of one or both strands of the phosphodiester backbone. Mammals have six DNA topoisomerases, this review focuses on one, DNA topoisomerase II beta (TOP2B). In the absence of TOP2B transcription of many developmentally regulated genes is altered. Long genes seem particularly susceptible to the lack of TOP2B. Biochemical studies of the role of TOP2B in transcription regulated by ligands such as nuclear hormones, growth factors and insulin has revealed PARP1 associated with TOP2B and also PRKDC, XRCC5 and XRCC6. Analysis of publicly available databases of protein interactions confirms these interactions and illustrates interactions with other key transcriptional regulators including TRIM28. TOP2B has been shown to interact with proteins involved in chromosome organisation including CTCF and RAD21. Comparison of publicly available Chip-seq datasets reveals the location at which these proteins interact with genes. The availability of resources such as large datasets of protein–protein interactions, e.g. BioGrid and IntAct and protein–DNA interactions such as Chip-seq in GEO enables scientists to extend models and propose new hypotheses.

Background Cardiotoxicity is a common complication following anthracycline chemotherapy and represents one of the serious adverse reactions affecting life, which severely limits the effective use of anthracyclines in cancer therapy. Although some genes have been investigated by individual studies, the comprehensive analysis of key genes and molecular regulatory network in anthracyclines-induced cardiotoxicity (AIC) is lacking but urgently needed. Methods The present study integrating several transcription profiling datasets aimed to identify key genes associated with AIC by weighted correlation network analysis (WGCNA) and differentially expressed analysis (DEA) and also constructed miRNA-transcription factor-gene regulatory network. A total of three transcription profiling datasets involving 47 samples comprising 41 rat heart tissues and 6 human induced pluripotent stem cell-derived cardiomyocytes (hiPSCMs) samples were enrolled. Results The WGCNA and DEA with E-MTAB-1168 identified 14 common genes affected by doxorubicin administrated by 4 weeks or 6 weeks. Functional and signal enrichment analyses revealed that these genes were mainly enriched in the regulation of heart contraction, muscle contraction, heart process, and oxytocin signaling pathway. Ten (Ryr2, Casq1, Fcgr2b, Postn, Tceal5, Ccn2, Tnfrsf12a, Mybpc2, Ankrd23, Scn3b) of the 14 genes were verified by another gene expression profile GSE154603. Importantly, three key genes (Ryr2, Tnfrsf12a, Scn3b) were further validated in a hiPSCMs-based in-vitro model. Additionally, the miRNA-transcription factor-gene regulatory revealed several top-ranked transcription factors including Tcf12, Ctcf, Spdef, Ebf1, Sp1, Rcor1 and miRNAs including miR-124-3p, miR-195-5p, miR-146a-5p, miR-17-5p, miR-15b-5p, miR-424-5p which may be involved in the regulation of genes associated with AIC. Conclusions Collectively, the current study suggested the important role of the key genes, oxytocin signaling pathway, and the miRNA-transcription factor-gene regulatory network in elucidating the molecular mechanism of AIC.

ObjectiveIn a recent genome-wide association study, a significant genetic association between rs34330 of CDKN1B and risk of systemic lupus erythematosus (SLE) in Han Chinese was identified. This study was undertaken to validate the reported association and elucidate the biochemical mechanisms underlying the effect of the variant.MethodsWe performed an allelic association analysis in patients with SLE, followed by a meta-analysis assessing genome-wide association data across 11 independent cohorts (n = 28,872). In silico bioinformatics analysis and experimental validation in SLE-relevant cell lines were applied to determine the functional consequences of rs34330.ResultsWe replicated a genetic association between SLE and rs34330 (meta-analysis P = 5.29 × 10-22 , odds ratio 0.84 [95% confidence interval 0.81-0.87]). Follow-up bioinformatics and expression quantitative trait locus analysis suggested that rs34330 is located in active chromatin and potentially regulates several target genes. Using luciferase and chromatin immunoprecipitation-real-time quantitative polymerase chain reaction, we demonstrated substantial allele-specific promoter and enhancer activity, and allele-specific binding of 3 histone marks (H3K27ac, H3K4me3, and H3K4me1), RNA polymerase II (Pol II), CCCTC-binding factor, and a critical immune transcription factor (interferon regulatory factor 1 [IRF-1]). Chromosome conformation capture revealed long-range chromatin interactions between rs34330 and the promoters of neighboring genes APOLD1 and DDX47, and effects on CDKN1B and the other target genes were directly validated by clustered regularly interspaced short palindromic repeat (CRISPR)-based genome editing. Finally, CRISPR/dead CRISPR-associated protein 9-based epigenetic activation/silencing confirmed these results. Gene-edited cell lines also showed higher levels of proliferation and apoptosis.ConclusionCollectively, these findings suggest a mechanism whereby the rs34330 risk allele (C) influences the presence of histone marks, RNA Pol II, and IRF-1 transcription factor to regulate expression of several target genes linked to proliferation and apoptosis. This process could potentially underlie the association of rs34330 with SLE.

CCCTC-binding factor (CTCF) plays fundamental roles in transcriptional regulation and chromatin architecture maintenance. CTCF is also a tumour suppressor frequently mutated in cancer, however, the structural and functional impact of mutations have not been examined. We performed molecular and structural characterisation of five cancer-specific CTCF missense zinc finger (ZF) mutations occurring within key intra- and inter-ZF residues. Functional characterisation of CTCF ZF mutations revealed a complete (L309P, R339W, R377H) or intermediate (R339Q) abrogation as well as an enhancement (G420D) of the anti-proliferative effects of CTCF. DNA binding at select sites was disrupted and transcriptional regulatory activities abrogated. Molecular docking and molecular dynamics confirmed that mutations in residues specifically contacting DNA bases or backbone exhibited loss of DNA binding. However, R339Q and G420D were stabilised by the formation of new primary DNA bonds, contributing to gain-of-function. Our data confirm that a spectrum of loss-, change- and gain-of-function impacts on CTCF zinc fingers are observed in cell growth regulation and gene regulatory activities. Hence, diverse cellular phenotypes of mutant CTCF are clearly explained by examining structure–function relationships. Supplementary Information The online version contains supplementary material available at 10.1007/s00018-021-03946-z.

ETO2 is a nuclear co-repressor, which plays a critical role in the regulation of the cell cycle, self-renewal capacity, and differentiation of hematopoietic progenitor cells. We identified novel fusion transcripts involving ETO2 and CTCF by RNA-seq in a multiple relapsed AML case. The CTCF-ETO2 and ETO2-CTCF chimeric genes were validated by RT-PCR and Sanger sequencing. In addition, both transcripts apparently promoted cell proliferation via JAK/STAT3 pathway that is sensitive to STAT3 inhibitors. The novel fusions may have prognostic value and pathogenic mechanisms in acute myeloid leukemia.

Binding sites of the chromatin regulator protein CTCF function as important landmarks in the human genome. The recently characterized CTCF-binding sites at LINE-1 repeats depend on another repeat-regulatory protein CGGBP1. These CGGBP1-dependent CTCF-binding sites serve as potential barrier elements for epigenetic marks such as H3K9me3. Such CTCF-binding sites are associated with asymmetric H3K9me3 levels as well as RNA levels in their flanks. The functions of these CGGBP1-dependent CTCF-binding sites remain unknown. By performing targeted studies on candidate CGGBP1-dependent CTCF-binding sites cloned in an SV40 promoter-enhancer episomal system we show that these regions act as inhibitors of ectopic transcription from the SV40 promoter. CGGBP1-dependent CTCF-binding sites that recapitulate their genomic function of loss of CTCF binding upon CGGBP1 depletion and H3K9me3 asymmetry in immediate flanks are also the ones that show the strongest inhibition of ectopic transcription. By performing a series of strand-specific reverse transcription PCRs we demonstrate that this ectopic transcription results in the synthesis of RNA from the SV40 promoter in a direction opposite to the downstream reporter gene in a strand-specific manner. The unleashing of the bidirectionality of the SV40 promoter activity and a breach of the transcription barrier seems to depend on depletion of CGGBP1 and loss of CTCF binding proximal to the SV40 promoter. RNA-sequencing reveals that CGGBP1-regulated CTCF-binding sites act as barriers to transcription at multiple locations genome-wide. These findings suggest a role of CGGBP1-dependent binding sites in restricting ectopic transcription.

Haoming Mai,* Haisheng Xie,* Jia Hou,* Haitao Chen, Bin Zhou, Jinlin Hou, Deke Jiang State Key Laboratory of Organ Failure Research, Guangdong Key Laboratory of Viral Hepatitis Research, Guangdong Institute of Liver Diseases, Department of Infectious Diseases and Hepatology Unit, Nanfang Hospital, Southern Medical University, Guangzhou, Guangdong, 510515, People’s Republic of China*These authors contributed equally to this workCorrespondence: Deke Jiang; Jinlin HouState Key Laboratory of Organ Failure Research, Guangdong Key Laboratory of Viral Hepatitis Research, Guangdong Institute of Liver Diseases, Department of Infectious Diseases and Hepatology Unit, Nanfang Hospital, Southern Medical University, Guangzhou, Guangdong, 510515, People’s Republic of ChinaTel +86-20-62786533; +86-20-61641941Email dekejiang17@smu.edu.cn; jlhousmu@163.comPurpose: Activation of actin cytoskeleton remodeling is an important stage preceding cancer cell metastasis. Previous genome-wide association studies (GWAS) have identified multiple hepatitis B virus (HBV)-related hepatocellular carcinoma (HCC)-associated risk loci. However, limited sample size or strict significance threshold of GWAS may cause HBV-related HCC risk-associated genetic loci to be undetected. We aimed to investigate the performance of the SNP rs13025377 in PPP1CB in HCC.Patients and Methods: We performed a case–control study including 1161 cases and 1353 controls to evaluate associations between single nucleotide polymorphisms (SNPs) from 98 actin-cytoskeleton regulatory genes and risk of HBV-related HCC. The effects of SNPs on HBV-related HCC risk were assessed under logistic regression model and corrected by false discovery rate (FDR).Results: We found that rs13025377 in PPP1CB was significantly associated with HBV-related HCC risk [odds ratio (OR) = 0.81, 95% confidence interval (CI) = 0.72∼ 0.91, P = 4.88× 10– 4]. The risk allele A of rs13025377 increased PPP1CB expression levels in normal liver tissue. SNP rs4665434 was tagged by rs13025377 (r2 = 0.9) and its protective allele disrupted CTCF and Cohesin motifs. According to public datasets, PPP1CB, CTCF and Cohesin expression levels are increased in tumor tissues. Kaplan–Meier plots demonstrated that higher PPP1CB expression was significantly associated with shorter overall survival (OS). Moreover, we observed strong correlation between CTCF, Cohesin, and PPP1CB in various liver tissues. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis confirmed that PPP1CB plays a role in HCC through actin-cytoskeleton regulation.Conclusion: Thus, these findings indicated that PPP1CB may be a key gene in actin-cytoskeleton regulation and rs13025377 contributes to the risk of HBV-related HCC by regulating PPP1CB expression.Keywords: genome-wide association studies, PPP1CB, CTCF, cohesin, SNP, HBV-related hepatocellular carcinoma

Coordinated chromatin interactions between enhancers and promoters are critical for gene regulation. The architectural protein CTCF mediates chromatin looping and is enriched at the boundaries of topologically associating domains (TADs), which are sub-megabase chromatin structures. In vitro CTCF depletion leads to a loss of TADs but has only limited effects over gene expression, challenging the concept that CTCF-mediated chromatin structures are a fundamental requirement for gene regulation. However, how CTCF and a perturbed chromatin structure impacts gene expression during development remains poorly understood. Here we link the loss of CTCF and gene regulation during patterning and organogenesis in a ctcf knockout zebrafish model. CTCF absence leads to loss of chromatin structure and affects the expression of thousands of genes, including many developmental regulators. Our results demonstrate the essential role of CTCF in providing the structural context for enhancer-promoter interactions, thus regulating developmental genes., J.L.G.-S. received funding from the ERC (Grant Agreement No. 740041), the Spanish Ministerio de Economía y Competitividad (Grant No. BFU2016-74961-P) and the institutional grant Unidad de Excelencia María de Maeztu (MDM-2016-0687). J.T. was funded by a 2019 Leonardo Grant for Researchers and Cultural Creators, BBVA Foundation. I.I.-A. acknowledges support from the Federation of European Biochemical Societies (FEBS Long-Term Fellowship). M.F. was funded by the European Union’s Horizon 2020 research and innovation program under the Marie Skłodowska-Curie grant agreement [#800396] and a Juan de la Cierva-Formación fellow from the Spanish Ministry of Science and Innovation (FJC2018-038233-I). JMS-P was funded by a postdoctoral fellowship from Junta de Andalucía (DOC_00512).

Chromatin interactions play important roles in regulating gene expression. However, the availability of genome-wide chromatin interaction data is limited. We develop a computational method, chromatin interaction neural network (ChINN), to predict chromatin interactions between open chromatin regions using only DNA sequences. ChINN predicts CTCF- and RNA polymerase II-associated and Hi-C chromatin interactions. ChINN shows good across-sample performances and captures various sequence features for chromatin interaction prediction. We apply ChINN to 6 chronic lymphocytic leukemia (CLL) patient samples and a published cohort of 84 CLL open chromatin samples. Our results demonstrate extensive heterogeneity in chromatin interactions among CLL patient samples. Ministry of Education (MOE) National Research Foundation (NRF) Published version This research is supported by the National Research Foundation (NRF) Singapore through an NRF Fellowship awarded to M.J.F (NRF-NRFF2012-054) and NTU start-up funds awarded to M.J.F. This research is supported by the RNA Biology Center at the Cancer Science Institute of Singapore, NUS, as part of funding under the Singapore Ministry of Education Academic Research Fund Tier 3 awarded to Daniel Tenen as lead PI with M.J.F as co-investigator (MOE2014-T3-1-006). This research is supported by a National Research Foundation Competitive Research Programme grant awarded to V.T. as lead PI and M.J.F. as co-PI (NRF-CRP17-2017-02). This research is supported by the National Research Foundation Singapore and the Singapore Ministry of Education under its Research Centres of Excellence initiative. This research is supported by a Ministry of Education Tier II grant awarded to M.J.F (T2EP30120-0020).

The 11 zinc finger (ZF) protein CTCF regulates topologically associating domain formation and transcription through selective binding to thousands of genomic sites. Here, we replaced endogenous CTCF in mouse embryonic stem cells with green-fluorescent-protein-tagged wild-type or mutant proteins lacking individual ZFs to identify additional determinants of CTCF positioning and function. While ZF1 and ZF8–ZF11 are not essential for cell survival, ZF8 deletion strikingly increases the DNA binding off-rate of mutant CTCF, resulting in reduced CTCF chromatin residence time. Loss of ZF8 results in widespread weakening of topologically associating domains, aberrant gene expression and increased genome-wide DNA methylation. Thus, important chromatin-templated processes rely on accurate CTCF chromatin residence time, which we propose depends on local sequence and chromatin context as well as global CTCF protein concentration.

High-throughput sequencing-based assays measure different biochemical activities pertaining to gene regulation, genome-wide. These activities include transcription factor (TF)–DNA binding, enhancer activity, open chromatin, and more. A major goal is to understand underlying sequence components, or motifs, that can explain the measured activity. It is usually not one motif but a combination of motifs bound by cooperatively acting proteins that confers activity to such regions. Furthermore, regions can be diverse, governed by different combinations of TFs/motifs. Current approaches do not take into account this issue of combinatorial diversity. We present a new statistical framework, cisDIVERSITY, which models regions as diverse modules characterized by combinations of motifs while simultaneously learning the motifs themselves. Because cisDIVERSITY does not rely on knowledge of motifs, modules, cell type, or organism, it is general enough to be applied to regions reported by most high-throughput assays. For example, in enhancer predictions resulting from different assays—GRO-cap, STARR-seq, and those measuring chromatin structure—cisDIVERSITY discovers distinct modules and combinations of TF binding sites, some specific to the assay. From protein–DNA binding data, cisDIVERSITY identifies potential cofactors of the profiled TF, whereas from ATAC-seq data, it identifies tissue-specific regulatory modules. Finally, analysis of single-cell ATAC-seq data suggests that regions open in one cell-state encode information about future states, with certain modules staying open and others closing down in the next time point.

Retinal Ganglion Cells (RGCs) lose their ability to grow axons during development. Adult RGCs thus fail to regenerate their axons after injury, leading to vision loss. To uncover mechanisms that promote regeneration of RGC axons, we identified transcription factors (TF) and open chromatin regions that are enriched in rat embryonic RGCs (high axon growth capacity) compared to postnatal RGCs (low axon growth capacity). We found that developmental stage-specific gene expression changes correlated with changes in promoter chromatin accessibility. Binding motifs for TFs such as CREB, CTCF, JUN and YY1 were enriched in the regions of the chromatin that were more accessible in embryonic RGCs. Proteomic analysis of purified rat RGC nuclei confirmed the expression of TFs with potential role in axon growth such as CREB, CTCF, YY1, and JUND. The CREB/ATF binding motif was widespread at the open chromatin region of known pro-regenerative TFs, supporting a role of CREB in regulating axon regeneration. Consistently, overexpression of CREB fused to the VP64 transactivation domain in mouse RGCs promoted axon regeneration after optic nerve injury. Our study provides a map of the chromatin accessibility during RGC development and highlights that TF associated with developmental axon growth can stimulate axon regeneration in mature RGC.

Cohesin is a key organizer of chromatin folding in eukaryotic cells. The two main activities of this ring-shaped protein complex are the maintenance of sister chromatid cohesion and the establishment of long-range DNA–DNA interactions through the process of loop extrusion. Although the basic principles of both cohesion and loop extrusion have been described, we still do not understand several crucial mechanistic details. One of such unresolved issues is the question of whether a cohesin ring topologically embraces DNA string(s) during loop extrusion. Here, we show that cohesin complexes residing on CTCF-occupied genomic sites in mammalian cells do not interact with DNA topologically. We assessed the stability of cohesin-dependent loops and cohesin association with chromatin in high-ionic-strength conditions in G1-synchronized HeLa cells. We found that increased salt concentration completely displaces cohesin from those genomic regions that correspond to CTCF-defined loop anchors. Unsurprisingly, CTCF-anchored cohesin loops also dissipate in these conditions. Because topologically engaged cohesin is considered to be salt resistant, our data corroborate a non-topological model of loop extrusion. We also propose a model of cohesin activity throughout the interphase, which essentially equates the termination of non-topological loop extrusion with topological loading of cohesin. This theoretical framework enables a parsimonious explanation of various seemingly contradictory experimental findings. Supplementary Information The online version contains supplementary material available at 10.1186/s13072-021-00411-w.

Three decades of research have established the CCCTC-binding factor (CTCF) as a ubiquitously expressed chromatin organizing factor and master regulator of gene expression. A new role for CTCF as a regulator of alternative splicing (AS) has now emerged. CTCF has been directly and indirectly linked to the modulation of AS at the individual transcript and at the transcriptome-wide level. The emerging role of CTCF-mediated regulation of AS involves diverse mechanisms; including transcriptional elongation, DNA methylation, chromatin architecture, histone modifications, and regulation of splicing factor expression and assembly. CTCF thereby appears to not only co-ordinate gene expression regulation but contributes to the modulation of transcriptomic complexity. In this review, we highlight previous discoveries regarding the role of CTCF in AS. In addition, we summarize detailed mechanisms by which CTCF mediates AS regulation. We propose opportunities for further research designed to examine the possible fate of CTCF-mediated alternatively spliced genes and associated biological consequences. CTCF has been widely acknowledged as the ‘master weaver of the genome’. Given its multiple connections, further characterization of CTCF’s emerging role in splicing regulation might extend its functional repertoire towards a ‘conductor of the splicing orchestra’.

DNA methylation, chromatin-binding proteins, and DNA looping are common components regulating genomic imprinting which leads to parent-specific monoallelic gene expression. Loss of methylation (LOM) at the human imprinting center 2 (IC2) on chromosome 11p15 is the most common cause of the imprinting overgrowth disorder Beckwith-Wiedemann Syndrome (BWS). Here, we report a familial transmission of a 7.6 kB deletion that ablates the core promoter of KCNQ1. This structural alteration leads to IC2 LOM and causes recurrent BWS. We find that occupancy of the chromatin organizer CTCF is disrupted proximal to the deletion, which causes chromatin architecture changes both in cis and in trans. We also profile the chromatin architecture of IC2 in patients with sporadic BWS caused by isolated LOM to identify conserved features of IC2 regulatory disruption. A strong interaction between CTCF sites around KCNQ1 and CDKN1C likely drive their expression on the maternal allele, while a weaker interaction involving the imprinting control region element may impede this connection and mediate gene silencing on the paternal allele. We present an imprinting model in which KCNQ1 transcription is necessary for appropriate CTCF binding and a novel chromatin conformation to drive allele-specific gene expression.

Animal interphase chromosomes are organized into topologically associating domains (TADs). How TADs are formed is not fully understood. Here, we combined high-throughput chromosome conformation capture and gene silencing to obtain insights into TAD dynamics in Xenopus tropicalis embryos. First, TAD establishment in X. tropicalis is similar to that in mice and flies and does not depend on zygotic genome transcriptional activation. This process is followed by further refinements in active and repressive chromatin compartments and the appearance of loops and stripes. Second, within TADs, higher self-interaction frequencies at one end of the boundary are associated with higher DNA occupancy of the architectural proteins CTCF and Rad21. Third, the chromatin remodeling factor ISWI is required for de novo TAD formation. Finally, TAD structures are variable in different tissues. Our work shows that X. tropicalis is a powerful model for chromosome architecture analysis and suggests that chromatin remodeling plays an essential role in de novo TAD establishment., Hi-C and single-molecule sequencing analysis provide an improved assembly of the Xenopus tropicalis genome and insights into three-dimensional genome dynamics throughout embryogenesis.

CCCTC-binding factor (CTCF), a zinc finger protein, is a transcription factor and regulator of chromatin structure. Forebrain excitatory neuron-specific CTCF deficiency contributes to inflammation via enhanced transcription of inflammation-related genes in the cortex and hippocampus. However, little is known about the long-term effect of CTCF deficiency on postnatal neurons, astrocytes, or microglia in the hippocampus of adult mice. To address this, we knocked out the Ctcf gene in forebrain glutamatergic neurons (Ctcf cKO) by crossing Ctcf-floxed mice with Camk2a-Cre mice and examined the hippocampi of 7.5-10-month-old male mice using immunofluorescence microscopy. We found obvious neuronal cell death and reactive gliosis in the hippocampal cornu ammonis (CA)1 in 7.5-10-month-old cKO mice. Prominent rod-shaped microglia that participate in immune surveillance were observed in the stratum pyramidale and radiatum layer, indicating a potential increase in inflammatory mediators released by hippocampal neurons. Although neuronal loss was not observed in CA3, and dentate gyrus (DG) CTCF depletion induced a significant increase in the number of microglia in the stratum oriens of CA3 and reactive microgliosis and astrogliosis in the molecular layer and hilus of the DG in 7.5-10-month-old cKO mice. These results suggest that long-term Ctcf deletion from forebrain excitatory neurons may contribute to reactive gliosis induced by neuronal damage and consequent neuronal loss in the hippocampal CA1, DG, and CA3 in sequence over 7 months of age. [BMB Reports 2021; 54(6): 317-322].

Mammalian genomes are organized and regulated through long-range chromatin interactions. Structural loops formed by CCCTC-binding factor (CTCF) and cohesin fold the genome into domains, while enhancers interact with promoters across vast genomic distances to regulate gene expression. Although genomics and fixed-cell imaging approaches help illuminate many aspects of chromatin interactions, temporal information is usually lost. Here, we discuss how 3D super-resolution live-cell imaging (SRLCI) can resolve open questions on the dynamic formation and dissolution of chromatin interactions. We discuss SRLCI experimental design, implementation strategies, and data interpretation and highlight associated pitfalls. We conclude that, while technically demanding, SRLCI approaches will likely emerge as a critical tool to dynamically probe 3D genome structure and function and to study enhancer–promoter interactions and chromatin looping.

Lay abstract Chronic kidney disease (CKD) is a growing public health problem, as its prevalence and incidence have almost doubled over the last three decades. Renal interstitial fibrosis is a final common pathological manifestation that contributes to the progression of CKD. CTCF can arrest renal interstitial fibrosis. In this study, we aimed to elucidate the therapeutic mechanism of CTCF in renal interstitial fibrosis in CKD. The results suggest that CTCF could inhibit miR-185-5p transcription, upregulate NPHS2 and thus activate the downstream cGMP-PKG signaling pathway, consequently alleviating renal interstitial fibrosis in CKD. By examining this relationship, our data support that CTCF overexpression, miR-185-5p inhibition or NPHS2 overexpression could suppress renal interstitial fibrosis and thereby impede CKD progression in mice with CKD. In this sense, the present study deepened our understanding of the etiology of CKD, and, more importantly, provided at least three novel therapeutic targets for the development of effective targeted therapy for CKD.

The functional engagement between an enhancer and its target promoter ensures precise gene transcription1. Understanding the basis of promoter choice by enhancers has important implications for health and disease. Here we report that functional loss of a preferred promoter can release its partner enhancer to loop to and activate an alternative promoter (or alternative promoters) in the neighbourhood. We refer to this target-switching process as 'enhancer release and retargeting'. Genetic deletion, motif perturbation or mutation, and dCas9-mediated CTCF tethering reveal that promoter choice by an enhancer can be determined by the binding of CTCF at promoters, in a cohesin-dependent manner-consistent with a model of 'enhancer scanning' inside the contact domain. Promoter-associated CTCF shows a lower affinity than that at chromatin domain boundaries and often lacks a preferred motif orientation or a partnering CTCF at the cognate enhancer, suggesting properties distinct from boundary CTCF. Analyses of cancer mutations, data from the GTEx project and risk loci from genome-wide association studies, together with a focused CRISPR interference screen, reveal that enhancer release and retargeting represents an overlooked mechanism that underlies the activation of disease-susceptibility genes, as exemplified by a risk locus for Parkinson's disease (NUCKS1-RAB7L1) and three loci associated with cancer (CLPTM1L-TERT, ZCCHC7-PAX5 and PVT1-MYC).

In both humans and mice, CTCF-binding elements form a series of interacting loops across the MHC class II (MHC-II) locus, and CTCF is required for maximal MHC-II gene expression. In humans, a CTCF-bound chromatin insulator termed XL9 and a super enhancer (SE) DR/DQ-SE situated in the intergenic region between HLA-DRB1 and HLA-DQA1 play critical roles in regulating MHC-II expression. In this study, we identify a similar SE, termed IA/IE-SE, located between H2-Eb1 and H2-Aa of the mouse that contains a CTCF site (C15) and a novel region of high histone H3K27 acetylation. A genetic knockout of C15 was created and its role on MHC-II expression tested on immune cells. We found that C15 deletion did not alter MHC-II expression in B cells, macrophages, and macrophages treated with IFN-γ because of functional redundancy of the remaining MHC-II CTCF sites. Surprisingly, embryonic fibroblasts derived from C15-deleted mice failed to induce MHC-II gene expression in response to IFN-γ, suggesting that at least in this developmental lineage, C15 was required. Examination of the three-dimensional interactions with C15 and the H2-Eb1 and H2-Aa promoters identified interactions within the novel region of high histone acetylation within the IA/IE-SE (termed N1) that contains a PU.1 binding site. CRISPR/Cas9 deletion of N1 altered chromatin interactions across the locus and resulted in reduced MHC-II expression. Together, these data demonstrate the functional redundancy of the MHC-II CTCF elements and identify a functionally conserved SE that is critical for maximal expression of MHC-II genes.

Studies of nuclear architecture using chromosome conformation capture methods have provided a detailed view of how chromatin folds in the 3D nuclear space. New variants of this technology now afford unprecedented resolution and allow the identification of ever smaller folding domains that offer new insights into the mechanisms by which this organization is established and maintained. Here we review recent results in this rapidly evolving field with an emphasis on CTCF function, with the goal of gaining a mechanistic understanding of the principles by which chromatin is folded in the eukaryotic nucleus.

While constitutive CCCTC-binding factor (CTCF)–binding sites are needed to maintain relatively invariant chromatin structures, such as topologically associating domains, the precise roles of CTCF to control cell-type–specific transcriptional regulation remain poorly explored. We examined CTCF occupancy in different types of primary blood cells derived from the same donor to elucidate a new role for CTCF in gene regulation during blood cell development. We identified dynamic, cell-type–specific binding sites for CTCF that colocalize with lineage-specific transcription factors. These dynamic sites are enriched for single-nucleotide polymorphisms that are associated with blood cell traits in different linages, and they coincide with the key regulatory elements governing hematopoiesis. CRISPR-Cas9–based perturbation experiments demonstrated that these dynamic CTCF-binding sites play a critical role in red blood cell development. Furthermore, precise deletion of CTCF-binding motifs in dynamic sites abolished interactions of erythroid genes, such as RBM38, with their associated enhancers and led to abnormal erythropoiesis. These results suggest a novel, cell-type–specific function for CTCF in which it may serve to facilitate interaction of distal regulatory emblements with target promoters. Our study of the dynamic, cell-type–specific binding and function of CTCF provides new insights into transcriptional regulation during hematopoiesis.

Human CTCF (hCTCF) is a major architectural protein in mammals. In Drosophila, the CTCF homologue (dCTCF) interacts with the BTB domain of the CP190 protein, which is involved in the establishment of open chromatin and activity of insulators. Previously, it was shown that the BTB protein Kaiso interacts with hCTCF and regulates its activity. We have carried out a detailed study of the interaction between these proteins in the yeast two-hybrid assay. Surprisingly, Kaiso did not interact with hCTCF and its Drosophila homologue. On the other hand, CP190 interacted with the C-terminus of hCTCF. The results obtained demonstrate that the interaction between CTCF and CP190 proteins is highly conserved. It is likely that humans have other BTB proteins that perform the functions described for the Drosophila CP190.

Numerous studies of relationship between epigenomic features have focused on their strong correlation across the genome, likely because such relationship can be easily identified by many established methods for correlation analysis. However, two features with little correlation may still colocalize at many genomic sites to implement important functions. There is no bioinformatic tool for researchers to specifically identify such feature pairs. Here, we develop a method to identify feature pairs in which two features have maximal colocalization minimal correlation (MACMIC) across the genome. By MACMIC analysis of 3306 feature pairs in 16 human cell types, we reveal a dual role of CCCTC-binding factor (CTCF) in epigenetic regulation of cell identity genes. Although super-enhancers are associated with activation of target genes, only a subset of super-enhancers colocalized with CTCF regulate cell identity genes. At super-enhancers colocalized with CTCF, CTCF is required for the active marker H3K27ac in cell types requiring the activation, and also required for the repressive marker H3K27me3 in other cell types requiring repression. Our work demonstrates the biological utility of the MACMIC analysis and reveals a key role for CTCF in epigenetic regulation of cell identity. The code for MACMIC is available at https://github.com/bxia888/MACMIC.

Recent advancements in sequencing and imaging technologies are providing new perspectives in solving the mystery of three-dimensional folding of genome in a nucleus. Chromosome conformation capture sequencing has discovered new chromatin structures such as topologically associated domains and loops in hundreds of kilobases. Super-resolution fluorescence microscopy with nanometer resolutions, in particular multiplexed approaches with sequence-specificity, has visualized chromatin structures from the rough folds of whole chromosomes to the fine loops of cis-regulatory elements in intact individual nuclei. Here, recent advancements in genome visualization tools with highly multiplexed labeling and reading are introduced. These imaging technologies have found ensemble behavior consistent to sequencing results, while unveiling single-cell variations. But, they also generated contradictory results on the roles of architectural proteins (like cohesion and CTCF) and enhancer-promoter interactions. Live-cell labeling methods for imaging specific genomic loci, especially the CRISPR/dCas9 system, are reviewed in order to give perspectives in the emergence of tools for visualizing genome structural dynamics.

Switch/sucrose-nonfermentable (SWI/SNF) chromatin-remodeling complexes are critical regulators of chromatin dynamics during transcription, DNA replication, and DNA repair. A recently identified SWI/SNF subcomplex termed GLTSCR1/1L-BAF (GBAF; or “noncanonical BAF”, ncBAF) uniquely contains bromodomain-containing protein BRD9 and glioma tumor suppressor candidate region 1 (GLTSCR1) or its paralog GLTSCR1-like (GLTSCR1L). Recent studies have identified a unique dependency on GBAF (ncBAF) complexes in synovial sarcoma and malignant rhabdoid tumors, both of which possess aberrations in canonical BAF (cBAF) and Polybromo-BAF (PBAF) complexes. Dependencies on GBAF in malignancies without SWI/SNF aberrations, however, are less defined. Here, we show that GBAF, particularly its BRD9 subunit, is required for the viability of prostate cancer cell lines in vitro and for optimal xenograft tumor growth in vivo. BRD9 interacts with androgen receptor (AR) and CCCTC-binding factor (CTCF), and modulates AR-dependent gene expression. The GBAF complex exhibits overlapping genome localization and transcriptional targets as bromodomain and extraterminal domain–containing (BET) proteins, which are established AR coregulators. Our results demonstrate that GBAF is critical for coordinating SWI/SNF–BET cooperation and uncover a new druggable target for AR-positive prostate cancers, including those resistant to androgen deprivation or antiandrogen therapies. Significance: Advanced prostate cancers resistant to androgen receptor antagonists are still susceptible to nontoxic BRD9 inhibitors, making them a promising alternative for halting AR signaling in progressed disease.

Recent technological advancements in the field of chromatin biology have rewritten the textbook on nuclear organization. We now appreciate that the folding of chromatin in the three-dimensional space (i.e. its 3D "architecture") is non-random, hierarchical, and highly complex. While 3D chromatin structure is partially encoded in the primary sequence and thereby broadly conserved across cell types and states, a substantial portion of the genome seems to be dynamic during development or in disease. Moreover, there is growing evidence that at least some of the 3D structure of chromatin is functionally linked to gene regulation, both being modulated by and impacting on multiple nuclear processes (including DNA replication, transcription, and RNA splicing). In recent years, these new concepts have nourished several investigations about the functional role of 3D chromatin topology dynamics in the heart during development and disease. This review aims to provide a comprehensive overview of our current understanding in this field, and to discuss how this knowledge can inform further research as well as clinical practice.

The three-dimensional (3D) organization of the human genome is of crucial importance for gene regulation, and the CCCTC-binding factor (CTCF) plays an important role in chromatin interactions. However, it is still unclear what sequence patterns in addition to CTCF motif pairs determine chromatin loop formation. To discover the underlying sequence patterns, we have developed a deep learning model, called DeepCTCFLoop, to predict whether a chromatin loop can be formed between a pair of convergent or tandem CTCF motifs using only the DNA sequences of the motifs and their flanking regions. Our results suggest that DeepCTCFLoop can accurately distinguish the CTCF motif pairs forming chromatin loops from the ones not forming loops. It significantly outperforms CTCF-MP, a machine learning model based on word2vec and boosted trees, when using DNA sequences only. Furthermore, we show that DNA motifs binding to several transcription factors, including ZNF384, ZNF263, ASCL1, SP1, and ZEB1, may constitute the complex sequence patterns for CTCF-mediated chromatin loop formation. DeepCTCFLoop has also been applied to disease-associated sequence variants to identify candidates that may disrupt chromatin loop formation. Therefore, our results provide useful information for understanding the mechanism of 3D genome organization and may also help annotate and prioritize the noncoding sequence variants associated with human diseases.

During human spermatogenesis, germ cells undergo dynamic changes in chromatin organization/re-packaging and in transcriptomes. In order to better understand the underlying mechanism(s), scATAC-Seq of 5376 testicular cells from 3 normal men were performed. Data were analyzed in parallel with the scRNA-Seq data of human testicular cells. In all, 10 germ cell types associated with spermatogenesis and 6 testicular somatic cell types were identified, along with 142 024 peaks located in promoter, genebody and CpG Island. We had examined chromatin accessibility of all chromosomes, with chromosomes 19 and 17 emerged as the leading chromosomes that displayed high chromatin accessibility. In accessible chromatin regions, transcription factor-binding sites were identified and specific motifs with high frequencies at different spermatogenesis stages were detected, including CTCF, BORIS, NFY, DMRT6, EN1, ISL1 and GLI3. Two most remarkable observations were noted. First, TLE3 was specifically expressed in differentiating spermatogonia. Second, PFN4 was found to be involved in actin cytoskeletal organization during meiosis. More important, unique regions upstream of PFN4 and TLE3 were shown to display high accessibility, illustrating their significance in supporting human spermatogenesis.

Tumor metastasis remains an obstacle in cancer treatment and is responsible for most cancer-related deaths. Nm23-H1 is one of the first metastasis suppressor proteins discovered with the ability to inhibit metastasis of many cancers including breast, colon, and liver cancer. Although loss of Nm23-H1 is observed in aggressive cancers and correlated with metastatic potential, little is known regarding the mechanisms that regulate its cellular level. Here, we examined the mechanisms that control Nm23-H1 expression in breast cancer cells. Initial studies in aggressive MDA-MB-231 cells (expressing low Nm23-H1) and less invasive MCF-7 cells (expressing high Nm23-H1) revealed that mRNA levels correlated with protein expression, suggesting that transcriptional mechanisms may control Nm23-H1 expression. Truncational analysis of the Nm23-H1 promoter revealed a proximal and minimal promoter that harbor putative binding sites for transcription factors including CTCF and EGR1. CTCF and EGR1 induced Nm23-H1 expression and reduced cell migration of MDA-MB-231 cells. Moreover, CTCF and EGR1 were recruited to the Nm23-H1 promoter in MCF-7 cells and their expression correlated with Nm23-H1 levels. This study indicates that loss of Nm23-H1 in aggressive breast cancer is apparently caused by downregulation of CTCF and EGR1, which potentially drive Nm23-H1 expression to promote a less invasive phenotype.

The human genome is partitioned into a collection of genomic features, inclusive of genes, transposable elements, lamina interacting regions, early replicating control elements and cis-regulatory elements, such as promoters, enhancers, and anchors of chromatin interactions. Uneven distribution of these features within chromosomes gives rise to clusters, such as topologically associating domains (TADs), lamina-associated domains, clusters of cis-regulatory elements or large organized chromatin lysine (K) domains (LOCKs). Here we show that LOCKs from diverse histone modifications discriminate primitive from differentiated cell types. Active LOCKs (H3K4me1, H3K4me3 and H3K27ac) cover a higher fraction of the genome in primitive compared to differentiated cell types while repressive LOCKs (H3K9me3, H3K27me3 and H3K36me3) do not. Active LOCKs in differentiated cells lie proximal to highly expressed genes while active LOCKs in primitive cells tend to be bivalent. Genes proximal to bivalent LOCKs are minimally expressed in primitive cells. Furthermore, bivalent LOCKs populate TAD boundaries and are preferentially bound by regulators of chromatin interactions, including CTCF, RAD21 and ZNF143. Together, our results argue that LOCKs discriminate primitive from differentiated cell populations., Large organized chromatin lysine (K) domains (LOCKs) are a genomic feature derived from the cluster of modified nucleosomes. Here the authors provide evidence that LOCKs discriminate primitive from differentiated cell populations and that they relate to transcription regulating events.

Graphical abstract, Topologically associating domains, or TADs, play important roles in genome organization and gene regulation; however, they are often altered in diseases. High-throughput chromatin conformation capturing assays, such as Hi-C, can capture domains of increased interactions, and TADs and boundaries can be identified using well-established analytical tools. However, generating Hi-C data is expensive. In our study, we addressed the relationship between multi-omics data and higher-order chromatin structures using a newly developed machine-learning model called PredTAD. Our tool uses already-available and cost-effective datatypes such as transcription factor and histone modification ChIPseq data. Specifically, PredTAD utilizes both epigenetic and genetic features as well as neighboring information to classify the entire human genome as boundary or non-boundary regions. Our tool can predict boundary changes between normal and breast cancer genomes. Among the most important features for predicting boundary alterations were CTCF, subunits of cohesin (RAD21 and SMC3), and chromosome number, suggesting their roles in conserved and dynamic boundaries formation. Upon further analysis, we observed that genes near altered TAD boundaries were found to be involved in several important breast cancer signaling pathways such as Ras, Jak-STAT, and estrogen signaling pathways. We also discovered a TAD boundary alteration that contributes to RET oncogene overexpression. PredTAD can also successfully predict TAD boundary changes in other conditions and diseases. In conclusion, our newly developed machine learning tool allowed for a more complete understanding of the dynamic 3D chromatin structures involved in signaling pathway activation, altered gene expression, and disease state in breast cancer cells.

Genome-wide association studies have identified 3p21.31 as the main risk locus for severe COVID-19, although underlying mechanisms remain elusive. We perform an epigenomic dissection of 3p21.31, identifying a CTCF-dependent tissue-specific 3D regulatory chromatin hub that controls the activity of several chemokine receptor genes. Risk SNPs colocalize with regulatory elements and are linked to increased expression of CCR1, CCR2 and CCR5 in monocytes and macrophages. As excessive organ infiltration of inflammatory monocytes and macrophages is a hallmark of severe COVID-19, our findings provide a rationale for the genetic association of 3p21.31 variants with elevated risk of hospitalization upon SARS-CoV-2 infection. B.S. and R.W.H. are supported by Dutch Lung Foundation grant 4.1.18.226. G.S. was supported by the ‘Fundación Científica de la Asociación Española Contra el Cáncer’. R.S. is supported by an Erasmus MC Fellowship, a Dutch Lung Foundation Junior Investigator grant (4.2.19.041JO) and a VIDI grant (09150172010068) from the Dutch Research Council (NWO). A.F.v.O. was supported by a postdoctoral fellowship from the Foundation Recherche Médicale. S.S. was supported by an Agence National pour la Recherche” grant (ANR-18-CE12-0019)

The HoxD gene cluster is critical for proper limb formation in tetrapods. In the emerging limb buds, different sub-groups of Hoxd genes respond first to a proximal regulatory signal, then to a distal signal that organizes digits. These two regulations are exclusive from one another and emanate from two distinct TADs flanking HoxD, both containing a range of appropriate enhancer sequences. The telomeric TAD (T-DOM) contains several enhancers active in presumptive forearm cells and is divided into two sub-TADs separated by a CTCF-rich boundary, which defines two regulatory sub-modules. To understand the importance of this particular regulatory topology to control Hoxd gene transcription in time and space, we either deleted or inverted this sub-TAD boundary, eliminated the CTCF binding sites or inverted the entire T-DOM to exchange the respective positions of the two sub-TADs. The effects of such perturbations on the transcriptional regulation of Hoxd genes illustrate the requirement of this regulatory topology for the precise timing of gene activation. However, the spatial distribution of transcripts was eventually resumed, showing that the presence of enhancers sequences, rather than either their exact topology or a particular chromatin architecture, is the key factor. We also show that the affinity of enhancers to find their natural target genes can overcome the presence of both a strong TAD border and an unfavourable orientation of CTCF sites.SIGNIFICANCE STATEMENTMany genes important for vertebrate development are surrounded by series of remote enhancer sequences. Such regulatory landscapes and their target genes are usually located within the same chromatin domains, which appears to constrain the action of these regulatory sequences and hence to facilitate enhancer-promoter recognition and gene expression. We used the HoxD locus to assess the impact of modifying the regulatory topology upon gene activation in space and time. A series of chromosomal re-arrangements involving deletions and inversions reveals that the enhancer topology plays a role in the timing of gene activation. However, gene expression was often recovered, subsequently, illustrating the intrinsic capacity of some enhancers to find their target promoters despite an apparently adverse chromatin topology.

Gene transcription is regulated by distant regulatory elements via combinatorial binding of transcription factors. It is increasingly recognized that alterations in chromatin state and transcription factor binding in these distant regulatory elements may have key roles in cancer development. Here we focused on the first stages of oncogene-induced carcinogenic transformation, and characterized the regulatory network underlying transcriptional changes associated with this process. Using Hi-C data, we observe spatial coupling between differentially expressed genes and their differentially accessible regulatory elements and reveal two candidate transcription factors, p53 and CTCF, as determinants of transcriptional alterations at the early stages of oncogenic HRas-induced transformation in human mammary epithelial cells. Strikingly, the malignant transcriptional reprograming is promoted by redistribution of chromatin binding of these factors without major variation in their expression level. Our results demonstrate that alterations in the regulatory landscape have a major role in driving oncogene-induced transcriptional reprogramming., By integrating transcriptome, regulatory element and chromosome topology profiles, Schwartz et al identify redistribution of p53 and CTCF genomic binding to be associated with transcriptional and phenotypical changes during Ras-induced transformation of mammary epithelial cells.

Testicular germ cell tumours (TGCTs) respond well to cisplatin-based therapy. However, cisplatin resistance and poor outcomes do occur. It has been suggested that a shift towards DNA hypermethylation mediates cisplatin resistance in TGCT cells, although there is little direct evidence to support this claim. Here we utilized a series of isogenic cisplatin-resistant cell models and observed a strong association between cisplatin resistance in TGCT cells and a net increase in global CpG and non-CpG DNA methylation spanning regulatory, intergenic, genic and repeat elements. Hypermethylated loci were significantly enriched for repressive DNA segments, CTCF and RAD21 sites and lamina associated domains, suggesting that global nuclear reorganization of chromatin structure occurred in resistant cells. Hypomethylated CpG loci were significantly enriched for EZH2 and SUZ12 binding and H3K27me3 sites. Integrative transcriptome and methylome analyses showed a strong negative correlation between gene promoter and CpG island methylation and gene expression in resistant cells and a weaker positive correlation between gene body methylation and gene expression. A bidirectional shift between gene promoter and gene body DNA methylation occurred within multiple genes that was associated with upregulation of polycomb targets and downregulation of tumour suppressor genes. These data support the hypothesis that global remodelling of DNA methylation is a key factor in mediating cisplatin hypersensitivity and chemoresistance of TGCTs and furthers the rationale for hypomethylation therapy for refractory TGCT patients.

Both HOX gene expression and CTCF regulation have been well demonstrated to play a critical role in regulating maintenance of leukemic stem cells (LSCs) that are known to be resistant to BET inhibitor (BETi). To investigate the regulatory role of CTCF boundary in aberrant HOX gene expression and the therapeutic sensitivity of BETi in AML, we employed CRISPR-Cas9 genome editing technology to delete 47 base pairs of the CTCF binding motif which is located between HOXA7 and HOXA9 genes (CBS7/9) in different subtypes of AML with either MLL-rearrangement or NPM1 mutation. Our results revealed that HOXA9 is significantly downregulated in response to the CBS7/9 deletion. Moreover, CBS7/9 boundary deletion sensitized the BETi treatment reaction in both MOLM-13 and OCI-AML3 cells. To further examine whether BETi therapeutic sensitivity in AML is depended on the expression level of the HOXA9 gene, we overexpressed the HOXA9 in the CBS7/9 deleted AML cell lines, which can rescue and restore the resistance to BETi treatment of the CBS7/9 KO cells by activating MAPK signaling pathway. Deletion of CBS7/9 specifically decreased the recruitment of BRD4 and RNA pol II to the posterior HOXA genes, in which, a transcription elongation factor ELL3 is the key factor in regulating HOXA gene transcription monitored by CBS7/9 chromatin boundary. Thus, disruption of CBS7/9 boundary perturbs HOXA9 transcription and regulates BETi sensitivity in AML treatment. Moreover, alteration of CTCF boundaries in the oncogene loci may provide a novel strategy to overcome the drug resistance of LSCs. Graphical abstract.

CCCTC-binding factor (CTCF) sites interact with each other in the chromatin environment, establishing chromatin domains. Our previous study showed that interaction between CTCF sites is cell type-specific around the β-globin locus and is dependent on erythroid-specific activator GATA-1. To find out molecular mechanisms of the cell type-specific interaction, we directly inhibited GATA-1 binding to the β-globin enhancers by deleting its binding motifs and found that histone H3K27 acetylation (H3K27ac) was decreased at CTCF sites surrounding the β-globin locus, even though CTCF binding itself was maintained at the sites. Forced H3K27ac by Trichostatin A treatment or CBP/p300 KD affected the interactions between CTCF sites around the β-globin locus without changes in CTCF binding. Analysis of public ChIA-PET data revealed that H3K27ac is higher at CTCF sites forming short interactions than long interactions. GATA-1 was identified as a representative transcription factor that relates with genes present inside the short interactions in erythroid K562 cells. Depletion of GATA-1-reduced H3K27ac at CTCF sites near erythroid-specific enhancers. These results indicate that H3K27ac at CTCF sites is required for cell type-specific chromatin interactions between them. Tissue-specific activator GATA-1 appears to play a role in H3K27ac at CTCF sites in erythroid cells.

Background: The cohesin complex plays an essential role in genome organisation and cell division. A full complement of the cohesin complex and its regulators is important for normal development, since heterozygous mutations in genes encoding these components can be sufficient to produce a disease phenotype. The implication that genes encoding the cohesin subunits or cohesin regulators must be tightly controlled and resistant to variability in expression has not yet been formally tested.Methods: Here, we identify spatial-regulatory connections with potential to regulate expression of cohesin loci (Mitotic: SMC1A, SMC3, STAG1, STAG2, RAD21/RAD21-AS; Meiotic: SMC1B, STAG3, REC8, RAD21L1), cohesin-ring support genes (NIPBL, MAU2, WAPL, PDS5A, PDS5B) and CTCF, including linking their expression to that of other genes. We searched the genome-wide association studies (GWAS) catalogue for SNPs mapped or attributed to cohesin genes by GWAS (GWAS-attributed) and the GTEx catalogue for SNPs mapped to cohesin genes by cis-regulatory variants in one or more of 44 tissues across the human body (expression quantitative trail locus-attributed).Results: Connections that centre on the cohesin ring subunits provide evidence of coordinated regulation that has little tolerance for perturbation. We used the CoDeS3D SNP-gene attribution methodology to identify transcriptional changes across a set of genes coregulated with the cohesin loci that include biological pathways such as extracellular matrix production and proteasome-mediated protein degradation. Remarkably, many of the genes that are coregulated with cohesin loci are themselves intolerant to loss-of-function.Conclusions: The results highlight the importance of robust regulation of cohesin genes and implicate novel pathways that may be important in the human cohesinopathy disorders.

Background Lung tumors and normal lung tissues show large differences in epigenetic modification which can affect the chromosome structure and expression of genes. However, the epigenetic reprogramming in lung adenocarcinoma remains unclear. Methods and results With the bioinformatics analysis, we found that some activated super-enhancers (SEs) only appear in lung adenocarcinoma cells, and 781 abnormal activated super-enhancers (AASEs) were found. Not only are the traditional oncogenes found to be activated by AASEs, such as MET and SLC2A1, but also some new genes were activated by AASEs, which probably contributes to the carcinogenic process in lung cancer. The enrichment analysis of the genes activated by AASEs shows that the glycolysis process and cell proliferation were enhanced and the apoptotic process was negatively regulated. Two AASEs were separately knockout by CRISPR/Cas9 in A549, PC-9, and H1299 cell lines and the expression of target genes decreased. The motif of CTCF, SMARCA1, SOX4, FOXM1, IRF3, IRF7, and STAT2 was enriched in AASEs, supporting that the chromosome structure changed and these transcription factors would be the master regulators on the formation of AASEs. Conclusion This study provided comprehensive insight into the mechanisms of SEs, as well as a potential therapeutic target for lung cancer.

Although CCCTC binding factor (CTCF) has been demonstrated to play a variety of often contradictory roles in tumorigenesis, little is known about its function in the tumorigenesis of ocular melanoma. Here, we generated two artificial CTCF peptides (Decoy-CTCFs) combining the zinc finger domain of wild-type CTCF and artificial marker region. This Decoy-CTCF retained the DNA binding region but lost the functional regions of wild-type CTCF. Transferring artificial CTCF into ocular melanoma cells suppressed proliferation and migration in the tumor cells, while no effect was observed in normal cells. Intriguingly, we first showed that decoy-CTCF inhibited tumorigenesis by preventing the histone acetyltransferase EP300 from binding to the promoter of SELL. Thus SELL was a novel oncogene in the tumorigenesis of ocular melanoma. These studies provide efficient decoy CTCF-based therapeutic concept in malignant ocular melanoma and reveal the potential mechanism underlying decoy-based tumor therapy., Graphical Abstract, An artificial CTCF peptide (Decoy-CTCFs) prevents the histone acetyltransferase EP300 to interact the promoter of SELL for inactivating SELL expression and triggering efficient therapeutic efficacy in ocular tumors. These studies delineate a novel decoy CTCF-guided mechanism in tumor therapy and provide a promising therapeutic concept for ocular melanoma.

Colorectal cancer (CRC) mostly arises from progressive accumulation of somatic mutations within cells. Most commonly mutated genes like TP53, APC and KRAS can promote survival and proliferation of cancer cells. Although the molecular alterations and landscape of some specific mutations in CRC are well known, the presence of a somatic mutation signature related to genomic regions and epigenetic markers remain unclear. To find the signatures from a random distribution of somatic mutations in CRCs, we carried out enrichment analysis in different genomic regions and identified peaks of epigenetic markers. We validated that the mutation frequency in miRNA is dramatically higher than in flanking genomic regions. Moreover, we observed that somatic mutations in CRC and colon cancer cell lines are significantly enriched in CTCF binding sites. We also found these mutations are enriched for H3K27me3 in both normal sigmoid colon and colon cancer cell lines. Taken together, our findings suggest that there are some common somatic mutations signatures which provide new directions to study CRC.