CTCF mediates dosage- and sequence-context-dependent transcriptional insulation by forming local chromatin domains

Permeable TAD boundaries and their impact on genome-associated functions

TAD boundaries are genomic elements that separate biological processes in neighboring domains by blocking DNA loops that are formed through Cohesin-mediated loop extrusion. Most TAD boundaries consist of arrays of binding sites for the CTCF protein, whose interaction with the Cohesin complex blocks loop extrusion. TAD boundaries are not fully impermeable though and allow a limited amount of inter-TAD loop formation. Based on the reanalysis of Nano-C data, a multicontact Chromosome Conformation Capture assay, we propose a model whereby clustered CTCF binding sites promote the successive stalling of Cohesin and subsequent dissociation from the chromatin. A fraction of Cohesin nonetheless achieves boundary read-through. Due to a constant rate of Cohesin dissociation elsewhere in the genome, the maximum length of inter-TAD loops is restricted though. We speculate that the DNA-encoded organization of stalling sites regulates TAD boundary permeability and discuss implications for enhancer-promoter loop formation and other genomic processes.

Polymer physics models reveal structural folding features of single-molecule gene chromatin conformations

Here, we employ polymer physics models of chromatin to investigate the 3D folding of a 2Mb wide genomic region encompassing the human LTN1 gene, a crucial DNA locus involved in key cellular functions. Through extensive Molecular Dynamics simulations, we reconstruct in-silico the ensemble of single-molecule LTN1 3D structures, which we benchmark against recent in-situ Hi-C 2.0 data. The model-derived single molecules are then used to predict structural folding features at the single-cell level, providing testable predictions for super-resolution microscopy experiments.

SnapFISH-IMPUTE: an imputation method for multiplexed DNA FISH data

Chromatin spatial organization plays a crucial role in gene regulation. Recently developed and prospering multiplexed DNA FISH technologies enable direct visualization of chromatin conformation in the nucleus. However, incomplete data caused by limited detection efficiency can substantially complicate and impair downstream analysis. Here, we present SnapFISH-IMPUTE that imputes missing values in multiplexed DNA FISH data. Analysis on multiple published datasets shows that the proposed method preserves the distribution of pairwise distances between imaging loci, and the imputed chromatin conformations are indistinguishable from the observed conformations. Additionally, imputation greatly improves downstream analyses such as identifying enhancer-promoter loops and clustering cells into distinct cell types. SnapFISH-IMPUTE is freely available at https://github.com/hyuyu104/SnapFISH-IMPUTE .

CTCF mutation at R567 causes developmental disorders via 3D genome rearrangement and abnormal neurodevelopment

The three-dimensional genome structure organized by CTCF is required for development. Clinically identified mutations in CTCF have been linked to adverse developmental outcomes. Nevertheless, the underlying mechanism remains elusive. In this investigation, we explore the regulatory roles of a clinically relevant R567W point mutation, located within the 11th zinc finger of CTCF, by introducing this mutation into both murine models and human embryonic stem cell-derived cortical organoid models. Mice with homozygous CTCFR567W mutation exhibit growth impediments, resulting in postnatal mortality, and deviations in brain, heart, and lung development at the pathological and single-cell transcriptome levels. This mutation induces premature stem-like cell exhaustion, accelerates the maturation of GABAergic neurons, and disrupts neurodevelopmental and synaptic pathways. Additionally, it specifically hinders CTCF binding to peripheral motifs upstream to the core consensus site, causing alterations in local chromatin structure and gene expression, particularly at the clustered protocadherin locus. Comparative analysis using human cortical organoids mirrors the consequences induced by this mutation. In summary, this study elucidates the influence of the CTCFR567W mutation on human neurodevelopmental disorders, paving the way for potential therapeutic interventions.

Senescence of human pancreatic beta cells enhances functional maturation through chromatin reorganization and promotes interferon responsiveness

Senescent cells can influence the function of tissues in which they reside, and their propensity for disease. A portion of adult human pancreatic beta cells express the senescence marker p16, yet it is unclear whether they are in a senescent state, and how this affects insulin secretion. We analyzed single-cell transcriptome datasets of adult human beta cells, and found that p16-positive cells express senescence gene signatures, as well as elevated levels of beta-cell maturation genes, consistent with enhanced functionality. Senescent human beta-like cells in culture undergo chromatin reorganization that leads to activation of enhancers regulating functional maturation genes and acquisition of glucose-stimulated insulin secretion capacity. Strikingly, Interferon-stimulated genes are elevated in senescent human beta cells, but genes encoding senescence-associated secretory phenotype (SASP) cytokines are not. Senescent beta cells in culture and in human tissue show elevated levels of cytoplasmic DNA, contributing to their increased interferon responsiveness. Human beta-cell senescence thus involves chromatin-driven upregulation of a functional-maturation program, and increased responsiveness of interferon-stimulated genes, changes that could increase both insulin secretion and immune reactivity.

RUNX1 C-terminal mutations impair blood cell differentiation by perturbing specific enhancer-promoter networks

ABSTRACT

The transcription factor RUNX1 is a master regulator of hematopoiesis and is frequently mutated in myeloid malignancies. Mutations in its runt homology domain (RHD) frequently disrupt DNA binding and result in loss of RUNX1 function. However, it is not clearly understood how other RUNX1 mutations contribute to disease development. Here, we characterized RUNX1 mutations outside of the RHD. Our analysis of the patient data sets revealed that mutations within the C-terminus frequently occur in hematopoietic disorders. Remarkably, most of these mutations were nonsense or frameshift mutations and were predicted to be exempt from nonsense-mediated messenger RNA decay. Therefore, this class of mutation is projected to produce DNA-binding proteins that contribute to the pathogenesis in a distinct manner. To model this, we introduced the RUNX1R320∗ mutation into the endogenous gene locus and demonstrated the production of RUNX1R320∗ protein. Expression of RUNX1R320∗ resulted in the disruption of RUNX1 regulated processes such as megakaryocytic differentiation, through a transcriptional signature different from RUNX1 depletion. To understand the underlying mechanisms, we used Global RNA Interactions with DNA by deep sequencing (GRID-seq) to examine enhancer-promoter connections. We identified widespread alterations in the enhancer-promoter networks within RUNX1 mutant cells. Additionally, we uncovered enrichment of RUNX1R320∗ and FOXK2 binding at the MYC super enhancer locus, significantly upregulating MYC transcription and signaling pathways. Together, our study demonstrated that most RUNX1 mutations outside the DNA-binding domain are not subject to nonsense-mediated decay, producing protein products that act in concert with additional cofactors to dysregulate hematopoiesis through mechanisms distinct from those induced by RUNX1 depletion.

Longitudinal analysis of epigenome-wide DNA methylation reveals novel loci associated with BMI change in East Asians

BACKGROUND

Obesity is a global public health concern linked to chronic diseases such as cardiovascular disease and type 2 diabetes (T2D). Emerging evidence suggests that epigenetic modifications, particularly DNA methylation, may contribute to obesity. However, the molecular mechanism underlying the longitudinal change of BMI has not been well-explored, especially in East Asian populations.

METHODS

This study performed a longitudinal epigenome-wide association analysis of DNA methylation to uncover novel loci associated with BMI change in 533 individuals across two Chinese cohorts with repeated DNA methylation and BMI measurements over four years.

RESULTS

We identified three novel CpG sites (cg14671384, cg25540824, and cg10848724) significantly associated with BMI change. Two of the identified CpG sites were located in regions previously associated with body shape and basal metabolic rate. Annotation of the top 20 BMI change-associated CpGs revealed strong connections to obesity and T2D. Notably, these CpGs exhibited active regulatory roles and located in genes with high expression in the liver and digestive tract, suggesting a potential regulatory pathway from genome to phenotypes of energy metabolism and absorption via DNA methylation. Cross-sectional and longitudinal EWAS comparisons indicated different mechanisms between CpGs related to BMI and BMI change.

CONCLUSION

This study enhances our understanding of the epigenetic dynamics underlying BMI change and emphasizes the value of longitudinal analyses in deciphering the complex interplay between epigenetics and obesity.

Genome-wide Cas9-mediated screening of essential non-coding regulatory elements via libraries of paired single-guide RNAs

The functions of non-coding regulatory elements (NCREs), which constitute a major fraction of the human genome, have not been systematically studied. Here we report a method involving libraries of paired single-guide RNAs targeting both ends of an NCRE as a screening system for the Cas9-mediated deletion of thousands of NCREs genome-wide to study their functions in distinct biological contexts. By using K562 and 293T cell lines and human embryonic stem cells, we show that NCREs can have redundant functions, and that many ultra-conserved elements have silencer activity and play essential roles in cell growth and in cellular responses to drugs (notably, the ultra-conserved element PAX6_Tarzan may be critical for heart development, as removing it from human embryonic stem cells led to defects in cardiomyocyte differentiation). The high-throughput screen, which is compatible with single-cell sequencing, may allow for the identification of druggable NCREs.

Recruitment of CTCF to the SIRT1 promoter after Oxidative Stress mediates Cardioprotective Transcription

Because most DNA-binding transcription factors (dbTFs), including the architectural regulator CTCF, bind RNA and exhibit di-/multimerization, a central conundrum is whether these distinct properties are regulated post-transcriptionally to modulate transcriptional programs. Here, investigating stress-dependent activation of SIRT1, encoding an evolutionarily-conserved protein deacetylase, we show that induced phosphorylation of CTCF acts as a rheostat to permit CTCF occupancy of low-affinity promoter DNA sites to precisely the levels necessary. This CTCF recruitment to the SIRT1 promoter is eliciting a cardioprotective cardiomyocyte transcriptional activation program and provides resilience against the stress of the beating heart in vivo . Mice harboring a mutation in the conserved low-affinity CTCF promoter binding site exhibit an altered, cardiomyocyte-specific transcriptional program and a systolic heart failure phenotype. This transcriptional role for CTCF reveals that a covalent dbTF modification regulating signal-dependent transcription serves as a previously unsuspected component of the oxidative stress response.

Members of an array of zinc finger proteins specify distinct Hox chromatin boundaries

The partitioning of repressive from actively transcribed chromatin domains in mammalian cells fosters cell-type specific gene expression patterns. During differentiation, this partitioning is reconstructed, reflecting gene expression profiles appropriate to new cellular identities. Yet, the chromatin occupancy of the key chromatin insulator, CTCF, at the developmentally important Hox clusters remains unchanged during differentiation. Thus, dynamic changes in chromatin boundaries must entail other activities, such as the previously identified MAZ insulator. Given its requirement for chromatin loop formation, we examined cohesin-based chromatin occupancy without the known insulators, CTCF and MAZ, and identified a novel family of zinc finger proteins (ZNFs), some of which exhibit tissue-specific expression. Thus far, two of these novel ZNFs facilitate the formation of chromatin boundaries at the Hox clusters that are distinct from each other and from that of MAZ. PATZ1 was critical to the thoracolumbar boundary of Hox clusters in differentiating motor neurons in vitro , skeletal development in vivo, and to looping interactions within the genome. On the other hand, ZNF263 contributed to cervicothoracic boundaries in motor neurons. We propose that these novel insulating activities act in concert with cohesin, alone or combinatorially, with or without CTCF, to implement precise positional identity and cell fate during development.

Super-enhancer interactomes from single cells link clustering and transcription

Regulation of gene expression hinges on the interplay between enhancers and promoters, traditionally explored through pairwise analyses. Recent advancements in mapping genome folding, like GAM, SPRITE, and multi-contact Hi-C, have uncovered multi-way interactions among super-enhancers (SEs), spanning megabases, yet have not measured their frequency in single cells or the relationship between clustering and transcription. To close this gap, here we used multiplexed imaging to map the 3D positions of 376 SEs across thousands of mammalian nuclei. Notably, our single-cell images reveal that while SE-SE contacts are rare, SEs often form looser associations we termed "communities". These communities, averaging 4-5 SEs, assemble cooperatively under the combined effects of genomic tethers, Pol2 clustering, and nuclear compartmentalization. Larger communities are associated with more frequent and larger transcriptional bursts. Our work provides insights about the SE interactome in single cells that challenge existing hypotheses on SE clustering in the context of transcriptional regulation.

Extensive folding variability between homologous chromosomes in mammalian cells

Genetic variation and 3D chromatin structure have major roles in gene regulation. Due to challenges in mapping chromatin conformation with haplotype-specific resolution, the effects of genetic sequence variation on 3D genome structure and gene expression imbalance remain understudied. Here, we applied Genome Architecture Mapping (GAM) to a hybrid mouse embryonic stem cell (mESC) line with high density of single nucleotide polymorphisms (SNPs). GAM resolved haplotype-specific 3D genome structures with high sensitivity, revealing extensive allelic differences in chromatin compartments, topologically associating domains (TADs), long-range enhancer-promoter contacts, and CTCF loops. Architectural differences often coincide with allele-specific differences in gene expression, mediated by Polycomb repression. We show that histone genes are expressed with allelic imbalance in mESCs, are involved in haplotype-specific chromatin contact marked by H3K27me3, and are targets of Polycomb repression through conditional knockouts of Ezh2 or Ring1b. Our work reveals highly distinct 3D folding structures between homologous chromosomes, and highlights their intricate connections with allelic gene expression.

Bento: a toolkit for subcellular analysis of spatial transcriptomics data

The spatial organization of molecules in a cell is essential for their functions. While current methods focus on discerning tissue architecture, cell-cell interactions, and spatial expression patterns, they are limited to the multicellular scale. We present Bento, a Python toolkit that takes advantage of single-molecule information to enable spatial analysis at the subcellular scale. Bento ingests molecular coordinates and segmentation boundaries to perform three analyses: defining subcellular domains, annotating localization patterns, and quantifying gene-gene colocalization. We demonstrate MERFISH, seqFISH + , Molecular Cartography, and Xenium datasets. Bento is part of the open-source Scverse ecosystem, enabling integration with other single-cell analysis tools.

Toward a comprehensive view of gene architecture during transcription

The activation of genes within the nucleus of eukaryotic cells is a tightly regulated process, orchestrated by a complex interplay of various physical properties and interacting factors. Studying the multitude of components and features that collectively contribute to gene activation has proven challenging due to the complexities of simultaneously visualizing the dynamic and transiently interacting elements that coalesce within the small space occupied by each individual gene. However, various labeling and imaging advances are now starting to overcome this challenge, enabling visualization of gene activation at different lengths and timescales. In this review, we aim to highlight these microscopy-based advances and suggest how they can be combined to provide a comprehensive view of the mechanisms regulating gene activation.

3C methods in cancer research: recent advances and future prospects

In recent years, Hi-C technology has revolutionized cancer research by elucidating the mystery of three-dimensional chromatin organization and its role in gene regulation. This paper explored the impact of Hi-C advancements on cancer research by delving into high-resolution techniques, such as chromatin loops, structural variants, haplotype phasing, and extrachromosomal DNA (ecDNA). Distant regulatory elements interact with their target genes through chromatin loops. Structural variants contribute to the development and progression of cancer. Haplotype phasing is crucial for understanding allele-specific genomic rearrangements and somatic clonal evolution in cancer. The role of ecDNA in driving oncogene amplification and drug resistance in cancer cells has also been revealed. These innovations offer a deeper understanding of cancer biology and the potential for personalized therapies. Despite these advancements, challenges, such as the accurate mapping of repetitive sequences and precise identification of structural variants, persist. Integrating Hi-C with multiomics data is key to overcoming these challenges and comprehensively understanding complex cancer genomes. Thus, Hi-C is a powerful tool for guiding precision medicine in cancer research and treatment.

Enhancer-promoter specificity in gene transcription: molecular mechanisms and disease associations

Although often located at a distance from their target gene promoters, enhancers are the primary genomic determinants of temporal and spatial transcriptional specificity in metazoans. Since the discovery of the first enhancer element in simian virus 40, there has been substantial interest in unraveling the mechanism(s) by which enhancers communicate with their partner promoters to ensure proper gene expression. These research efforts have benefited considerably from the application of increasingly sophisticated sequencing- and imaging-based approaches in conjunction with innovative (epi)genome-editing technologies; however, despite various proposed models, the principles of enhancer-promoter interaction have still not been fully elucidated. In this review, we provide an overview of recent progress in the eukaryotic gene transcription field pertaining to enhancer-promoter specificity. A better understanding of the mechanistic basis of lineage- and context-dependent enhancer-promoter engagement, along with the continued identification of functional enhancers, will provide key insights into the spatiotemporal control of gene expression that can reveal therapeutic opportunities for a range of enhancer-related diseases.

The impact of DNA methylation on CTCF-mediated 3D genome organization

Cytosine DNA methylation is a highly conserved epigenetic mark in eukaryotes. Although the role of DNA methylation at gene promoters and repetitive elements has been extensively studied, the function of DNA methylation in other genomic contexts remains less clear. In the nucleus of mammalian cells, the genome is spatially organized at different levels, and strongly influences myriad genomic processes. There are a number of factors that regulate the three-dimensional (3D) organization of the genome, with the CTCF insulator protein being among the most well-characterized. Pertinently, CTCF binding has been reported as being DNA methylation-sensitive in certain contexts, perhaps most notably in the process of genomic imprinting. Therefore, it stands to reason that DNA methylation may play a broader role in the regulation of chromatin architecture. Here we summarize the current understanding that is relevant to both the mammalian DNA methylation and chromatin architecture fields and attempt to assess the extent to which DNA methylation impacts the folding of the genome. The focus is in early embryonic development and cellular transitions when the epigenome is in flux, but we also describe insights from pathological contexts, such as cancer, in which the epigenome and 3D genome organization are misregulated.

Spatially organized cellular communities form the developing human heart

The heart, which is the first organ to develop, is highly dependent on its form to function1,2. However, how diverse cardiac cell types spatially coordinate to create the complex morphological structures that are crucial for heart function remains unclear. Here we integrated single-cell RNA-sequencing with high-resolution multiplexed error-robust fluorescence in situ hybridization to resolve the identity of the cardiac cell types that develop the human heart. This approach also provided a spatial mapping of individual cells that enables illumination of their organization into cellular communities that form distinct cardiac structures. We discovered that many of these cardiac cell types further specified into subpopulations exclusive to specific communities, which support their specialization according to the cellular ecosystem and anatomical region. In particular, ventricular cardiomyocyte subpopulations displayed an unexpected complex laminar organization across the ventricular wall and formed, with other cell subpopulations, several cellular communities. Interrogating cell-cell interactions within these communities using in vivo conditional genetic mouse models and in vitro human pluripotent stem cell systems revealed multicellular signalling pathways that orchestrate the spatial organization of cardiac cell subpopulations during ventricular wall morphogenesis. These detailed findings into the cellular social interactions and specialization of cardiac cell types constructing and remodelling the human heart offer new insights into structural heart diseases and the engineering of complex multicellular tissues for human heart repair.

Role of H3K4 monomethylation in gene regulation

Methylation of histone H3 on the lysine-4 residue (H3K4me) is found throughout the eukaryotic domain, and its initial discovery as a conserved epigenetic mark of active transcription from yeast to mammalian cells has contributed to the histone code hypothesis. However, recent studies have raised questions on whether the different forms of H3K4me play a direct role in gene regulation or are simply by-products of the transcription process. Here, we review the often-conflicting experimental evidence, focusing on the monomethylation of lysine 4 on histone H3 that has been linked to the transcriptional state of enhancers in metazoans. We suggest that this epigenetic mark acts in a context-dependent manner to directly facilitate the transcriptional output of the genome and the establishment of cellular identity.

Boundary stacking interactions enable cross-TAD enhancer-promoter communication during limb development

Although promoters and their enhancers are frequently contained within a topologically associating domain (TAD), some developmentally important genes have their promoter and enhancers within different TADs. Hypotheses about molecular mechanisms enabling cross-TAD interactions remain to be assessed. To test these hypotheses, we used optical reconstruction of chromatin architecture to characterize the conformations of the Pitx1 locus on single chromosomes in developing mouse limbs. Our data support a model in which neighboring boundaries are stacked as a result of loop extrusion, bringing boundary-proximal cis-elements into contact. This stacking interaction also contributes to the appearance of architectural stripes in the population average maps. Through molecular dynamics simulations, we found that increasing boundary strengths facilitates the formation of the stacked boundary conformation, counter-intuitively facilitating border bypass. This work provides a revised view of the TAD borders' function, both facilitating and preventing cis-regulatory interactions, and introduces a framework to distinguish border-crossing from border-respecting enhancer-promoter pairs.

Epigenetic regulation in cancer

Epigenetic modifications are defined as heritable changes in gene activity that do not involve changes in the underlying DNA sequence. The oncogenic process is driven by the accumulation of alterations that impact genome's structure and function. Genetic mutations, which directly disrupt the DNA sequence, are complemented by epigenetic modifications that modulate gene expression, thereby facilitating the acquisition of malignant characteristics. Principals among these epigenetic changes are shifts in DNA methylation and histone mark patterns, which promote tumor development and metastasis. Notably, the reversible nature of epigenetic alterations, as opposed to the permanence of genetic changes, positions the epigenetic machinery as a prime target in the discovery of novel therapeutics. Our review delves into the complexities of epigenetic regulation, exploring its profound effects on tumor initiation, metastatic behavior, metabolic pathways, and the tumor microenvironment. We place a particular emphasis on the dysregulation at each level of epigenetic modulation, including but not limited to, the aberrations in enzymes responsible for DNA methylation and histone modification, subunit loss or fusions in chromatin remodeling complexes, and the disturbances in higher-order chromatin structure. Finally, we also evaluate therapeutic approaches that leverage the growing understanding of chromatin dysregulation, offering new avenues for cancer treatment.

DEAD-Box Helicase 17 exacerbates non-alcoholic steatohepatitis via transcriptional repression of cyp2c29, inducing hepatic lipid metabolism disorder and eliciting the activation of M1 macrophages

OBJECTIVE

Our study was to elucidate the role of RNA helicase DEAD-Box Helicase 17 (DDX17) in NAFLD and to explore its underlying mechanisms.

METHODS

We created hepatocyte-specific Ddx17-deficient mice aim to investigate the impact of Ddx17 on NAFLD induced by a high-fat diet (HFD) as well as methionine and choline-deficient l-amino acid diet (MCD) in adult male mice. RNA-seq and lipidomic analyses were conducted to depict the metabolic landscape, and CUT&Tag combined with chromatin immunoprecipitation (ChIP) and luciferase reporter assays were conducted.

RESULTS

In this work, we observed a notable increase in DDX17 expression in the livers of patients with NASH and in murine models of NASH induced by HFD or MCD. After introducing lentiviruses into hepatocyte L02 for DDX17 knockdown or overexpression, we found that lipid accumulation induced by palmitic acid/oleic acid (PAOA) in L02 cells was noticeably weakened by DDX17 knockdown but augmented by DDX17 overexpression. Furthermore, hepatocyte-specific DDX17 knockout significantly alleviated hepatic steatosis, inflammatory response and fibrosis in mice after the administration of MCD and HFD. Mechanistically, our analysis of RNA-seq and CUT&Tag results combined with ChIP and luciferase reporter assays indicated that DDX17 transcriptionally represses Cyp2c29 gene expression by cooperating with CCCTC binding factor (CTCF) and DEAD-Box Helicase 5 (DDX5). Using absolute quantitative lipidomics analysis, we identified a hepatocyte-specific DDX17 deficiency that decreased lipid accumulation and altered lipid composition in the livers of mice after MCD administration. Based on the RNA-seq analysis, our findings suggest that DDX17 could potentially have an impact on the modulation of lipid metabolism and the activation of M1 macrophages in murine NASH models.

CONCLUSION

These results imply that DDX17 is involved in NASH development by promoting lipid accumulation in hepatocytes, inducing the activation of M1 macrophages, subsequent inflammatory responses and fibrosis through the transcriptional repression of Cyp2c29 in mice. Therefore, DDX17 holds promise as a potential drug target for the treatment of NASH.

ZNF143 deletion alters enhancer/promoter looping and CTCF/cohesin geometry

The transcription factor ZNF143 contains a central domain of seven zinc fingers in a tandem array and is involved in 3D genome construction. However, the mechanism by which ZNF143 functions in chromatin looping remains unclear. Here, we show that ZNF143 directionally recognizes a diverse range of genomic sites directly within enhancers and promoters and is required for chromatin looping between these sites. In addition, ZNF143 is located between CTCF and cohesin at numerous CTCF sites, and ZNF143 removal narrows the space between CTCF and cohesin. Moreover, genetic deletion of ZNF143, in conjunction with acute CTCF degradation, reveals that ZNF143 and CTCF collaborate to regulate higher-order topological chromatin organization. Finally, CTCF depletion enlarges direct ZNF143 chromatin looping. Thus, ZNF143 is recruited by CTCF to the CTCF sites to regulate CTCF/cohesin configuration and TAD (topologically associating domain) formation, whereas directional recognition of genomic DNA motifs directly by ZNF143 itself regulates promoter activity via chromatin looping.

SnapFISH-IMPUTE: an imputation method for multiplexed DNA FISH data

Chromatin spatial organization plays a crucial role in gene regulation. Recently developed and prospering multiplexed DNA FISH technologies enable direct visualization of chromatin conformation in nucleus. However, incomplete data caused by limited detection efficiency can substantially complicate and impair downstream analysis. Here, we present SnapFISH-IMPUTE that imputes missing values in multiplexed DNA FISH data. Analysis on multiple published datasets shows that the proposed method preserves the distribution of pairwise distances between imaging loci, and the imputed chromatin conformations are indistinguishable from the observed conformations. Additionally, imputation greatly improves downstream analyses such as identifying enhancer-promoter loops and clustering cells into distinct cell types. SnapFISH-IMPUTE is freely available at https://github.com/hyuyu104/SnapFISH-IMPUTE.

Cancer Evolution: A Multifaceted Affair

Cancer cells adapt and survive through the acquisition and selection of molecular modifications. This process defines cancer evolution. Building on a theoretical framework based on heritable genetic changes has provided insights into the mechanisms supporting cancer evolution. However, cancer hallmarks also emerge via heritable nongenetic mechanisms, including epigenetic and chromatin topological changes, and interactions between tumor cells and the tumor microenvironment. Recent findings on tumor evolutionary mechanisms draw a multifaceted picture where heterogeneous forces interact and influence each other while shaping tumor progression. A comprehensive characterization of the cancer evolutionary toolkit is required to improve personalized medicine and biomarker discovery.

SIGNIFICANCE

Tumor evolution is fueled by multiple enabling mechanisms. Importantly, genetic instability, epigenetic reprogramming, and interactions with the tumor microenvironment are neither alternative nor independent evolutionary mechanisms. As demonstrated by findings highlighted in this perspective, experimental and theoretical approaches must account for multiple evolutionary mechanisms and their interactions to ultimately understand, predict, and steer tumor evolution.

Rapid and Quantitative Functional Interrogation of Human Enhancer Variant Activity in Live Mice

Functional analysis of non-coding variants associated with human congenital disorders remains challenging due to the lack of efficient in vivo models. Here we introduce dual-enSERT, a robust Cas9-based two-color fluorescent reporter system which enables rapid, quantitative comparison of enhancer allele activities in live mice of any genetic background. We use this new technology to examine and measure the gain- and loss-of-function effects of enhancer variants linked to limb polydactyly, autism, and craniofacial malformation. By combining dual-enSERT with single-cell transcriptomics, we characterize variant enhancer alleles at cellular resolution, thereby implicating candidate molecular pathways in pathogenic enhancer misregulation. We further show that independent, polydactyly-linked enhancer variants lead to ectopic expression in the same cell populations, indicating shared genetic mechanisms underlying non-coding variant pathogenesis. Finally, we streamline dual-enSERT for analysis in F0 animals by placing both reporters on the same transgene separated by a synthetic insulator. Dual-enSERT allows researchers to go from identifying candidate enhancer variants to analysis of comparative enhancer activity in live embryos in under two weeks.

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.

Conserved and divergent gene regulatory programs of the mammalian neocortex

Divergence of cis-regulatory elements drives species-specific traits1, but how this manifests in the evolution of the neocortex at the molecular and cellular level remains unclear. Here we investigated the gene regulatory programs in the primary motor cortex of human, macaque, marmoset and mouse using single-cell multiomics assays, generating gene expression, chromatin accessibility, DNA methylome and chromosomal conformation profiles from a total of over 200,000 cells. From these data, we show evidence that divergence of transcription factor expression corresponds to species-specific epigenome landscapes. We find that conserved and divergent gene regulatory features are reflected in the evolution of the three-dimensional genome. Transposable elements contribute to nearly 80% of the human-specific candidate cis-regulatory elements in cortical cells. Through machine learning, we develop sequence-based predictors of candidate cis-regulatory elements in different species and demonstrate that the genomic regulatory syntax is highly preserved from rodents to primates. Finally, we show that epigenetic conservation combined with sequence similarity helps to uncover functional cis-regulatory elements and enhances our ability to interpret genetic variants contributing to neurological disease and traits.

Topology regulatory elements: From shaping genome architecture to gene regulation

The importance of 3D genome topology in the control of gene expression is becoming increasingly apparent, while regulatory mechanisms remain incompletely understood. Several recent studies have identified architectural elements that influence developmental gene expression by shaping locus topology. We refer to these elements as topological regulatory elements (TopoREs) to reflect their dual roles in genome organisation and gene expression. Importantly, these elements do not harbour autonomous transcriptional activation capacity, and instead appear to facilitate enhancer-promoter interactions, contributing to robust and precise timing of transcription. We discuss examples of TopoREs from two classes that are either dependent or independent of CTCF binding. Importantly, identification and interpretation of TopoRE function may shed light on multiple aspects of gene regulation, including the relationship between enhancer-promoter proximity and transcription, and enhancer-promoter specificity. Ultimately, understanding TopoRE diversity and function will aid in the interpretation of how human sequence variation can impact transcription and contribute to disease phenotypes.

Droplet-based single-cell joint profiling of histone modifications and transcriptomes

We previously reported Paired-Tag, a combinatorial indexing-based method that can simultaneously map histone modifications and gene expression at single-cell resolution at scale. However, the lengthy procedure of Paired-Tag has hindered its general adoption in the community. To address this bottleneck, we developed a droplet-based Paired-Tag protocol that is faster and more accessible than the previous method. Using cultured mammalian cells and primary brain tissues, we demonstrate its superior performance at identifying candidate cis-regulatory elements and associating their dynamic chromatin state to target gene expression in each constituent cell type in a complex tissue.

Context-dependent perturbations in chromatin folding and the transcriptome by cohesin and related factors

Cohesin regulates gene expression through context-specific chromatin folding mechanisms such as enhancer-promoter looping and topologically associating domain (TAD) formation by cooperating with factors such as cohesin loaders and the insulation factor CTCF. We developed a computational workflow to explore how three-dimensional (3D) structure and gene expression are regulated collectively or individually by cohesin and related factors. The main component is CustardPy, by which multi-omics datasets are compared systematically. To validate our methodology, we generated 3D genome, transcriptome, and epigenome data before and after depletion of cohesin and related factors and compared the effects of depletion. We observed diverse effects on the 3D genome and transcriptome, and gene expression changes were correlated with the splitting of TADs caused by cohesin loss. We also observed variations in long-range interactions across TADs, which correlated with their epigenomic states. These computational tools and datasets will be valuable for 3D genome and epigenome studies.

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.

Whole exome sequencing highlights rare variants in CTCF, DNMT1, DNMT3A, EZH2 and SUV39H1 as associated with FSHD

Introduction: Despite the progress made in the study of Facioscapulohumeral Dystrophy (FSHD), the wide heterogeneity of disease complicates its diagnosis and the genotype-phenotype correlation among patients and within families. In this context, the present work employed Whole Exome Sequencing (WES) to investigate known and unknown genetic contributors that may be involved in FSHD and may represent potential disease modifiers, even in presence of a D4Z4 Reduced Allele (DRA). Methods: A cohort of 126 patients with clinical signs of FSHD were included in the study, which were characterized by D4Z4 sizing, methylation analysis and WES. Specific protocols were employed for D4Z4 sizing and methylation analysis, whereas the Illumina® Next-Seq 550 system was utilized for WES. The study included both patients with a DRA compatible with FSHD diagnosis and patients with longer D4Z4 alleles. In case of patients harboring relevant variants from WES, the molecular analysis was extended to the family members. Results: The WES data analysis highlighted 20 relevant variants, among which 14 were located in known genetic modifiers (SMCHD1, DNMT3B and LRIF1) and 6 in candidate genes (CTCF, DNMT1, DNMT3A, EZH2 and SUV39H1). Most of them were found together with a permissive short (4-7 RU) or borderline/long DRA (8-20 RU), supporting the possibility that different genes can contribute to disease heterogeneity in presence of a FSHD permissive background. The segregation and methylation analysis among family members, together with clinical findings, provided a more comprehensive picture of patients. Discussion: Our results support FSHD pathomechanism being complex with a multigenic contribution by several known (SMCHD1, DNMT3B, LRIF1) and possibly other candidate genes (CTCF, DNMT1, DNMT3A, EZH2, SUV39H1) to disease penetrance and expressivity. Our results further emphasize the importance of extending the analysis of molecular findings within the proband's family, with the purpose of providing a broader framework for understanding single cases and allowing finer genotype-phenotype correlations in FSHD-affected families.

Competition between transcription and loop extrusion modulates promoter and enhancer dynamics

The spatiotemporal configuration of genes with distal regulatory elements, and the impact of chromatin mobility on transcription, remain unclear. Loop extrusion is an attractive model for bringing genetic elements together, but how this functionally interacts with transcription is also largely unknown. We combine live tracking of genomic loci and nascent transcripts with molecular dynamics simulations to assess the spatiotemporal arrangement of the Sox2 gene and its enhancer, in response to a battery of perturbations. We find a close link between chromatin mobility and transcriptional status: active elements display more constrained mobility, consistent with confinement within specialized nuclear sites, and alterations in enhancer mobility distinguish poised from transcribing alleles. Strikingly, we find that whereas loop extrusion and transcription factor-mediated clustering contribute to promoter-enhancer proximity, they have antagonistic effects on chromatin dynamics. This provides an experimental framework for the underappreciated role of chromatin dynamics in genome regulation.

SnapFISH: a computational pipeline to identify chromatin loops from multiplexed DNA FISH data

Multiplexed DNA fluorescence in situ hybridization (FISH) imaging technologies have been developed to map the folding of chromatin fibers at tens of nanometers and up to several kilobases in resolution in single cells. However, computational methods to reliably identify chromatin loops from such imaging datasets are still lacking. Here we present a Single-Nucleus Analysis Pipeline for multiplexed DNA FISH (SnapFISH), to process the multiplexed DNA FISH data and identify chromatin loops. SnapFISH can identify known chromatin loops from mouse embryonic stem cells with high sensitivity and accuracy. In addition, SnapFISH obtains comparable results of chromatin loops across datasets generated from diverse imaging technologies. SnapFISH is freely available at https://github.com/HuMingLab/SnapFISH .

Inferring CTCF-binding patterns and anchored loops across human tissues and cell types

CCCTC-binding factor (CTCF) is a transcription regulator with a complex role in gene regulation. The recognition and effects of CTCF on DNA sequences, chromosome barriers, and enhancer blocking are not well understood. Existing computational tools struggle to assess the regulatory potential of CTCF-binding sites and their impact on chromatin loop formation. Here we have developed a deep-learning model, DeepAnchor, to accurately characterize CTCF binding using high-resolution genomic/epigenomic features. This has revealed distinct chromatin and sequence patterns for CTCF-mediated insulation and looping. An optimized implementation of a previous loop model based on DeepAnchor score excels in predicting CTCF-anchored loops. We have established a compendium of CTCF-anchored loops across 52 human tissue/cell types, and this suggests that genomic disruption of these loops could be a general mechanism of disease pathogenesis. These computational models and resources can help investigate how CTCF-mediated cis-regulatory elements shape context-specific gene regulation in cell development and disease progression.

Spatial and temporal organization of the genome: Current state and future aims of the 4D nucleome project

The four-dimensional nucleome (4DN) consortium studies the architecture of the genome and the nucleus in space and time. We summarize progress by the consortium and highlight the development of technologies for (1) mapping genome folding and identifying roles of nuclear components and bodies, proteins, and RNA, (2) characterizing nuclear organization with time or single-cell resolution, and (3) imaging of nuclear organization. With these tools, the consortium has provided over 2,000 public datasets. Integrative computational models based on these data are starting to reveal connections between genome structure and function. We then present a forward-looking perspective and outline current aims to (1) delineate dynamics of nuclear architecture at different timescales, from minutes to weeks as cells differentiate, in populations and in single cells, (2) characterize cis-determinants and trans-modulators of genome organization, (3) test functional consequences of changes in cis- and trans-regulators, and (4) develop predictive models of genome structure and function.

TADs: Dynamic structures to create stable regulatory functions

Mammalian chromosomes are organized at different length scales within the cell nucleus. Topologically Associating Domains (TADs) are structural units of 3D genome organization with functions in gene regulation, DNA replication, recombination and repair. Whereas TADs were initially interpreted as insulated domains, recent studies are revealing that these domains should be interpreted as dynamic collections of actively extruding loops. This process of loop extrusion is subsequently blocked at dedicated TAD boundaries, thereby promoting intra-domain interactions over their surroundings. In this review, we discuss how mammalian TAD structure can emerge from this dynamic process and we discuss recent evidence that TAD boundaries can have regulatory functions.

Optical genome mapping in an atypical Pelizaeus-Merzbacher prenatal challenge

Pathogenic genetic variants represent a challenge in prenatal counseling, especially when clinical presentation in familial carriers is atypical. We describe a prenatal case involving a microarray-detected duplication of PLP1 which causes X-linked Pelizaeus-Merzbacher disease, a progressive hypomyelinating leukodystrophy. Because of atypical clinical presentation in an older male child, the duplication was examined using a novel technology, optical genome mapping, and was found to be an inverted duplication, which has not been previously described. Simultaneously, segregation analysis identified another healthy adult male carrier of this unique structural rearrangement. The novel PLP1 structural variant was reclassified, and a healthy boy was delivered. In conclusion, we suggest that examining structural variants with novel methods is warranted especially in cases with atypical clinical presentation and may in these cases lead to improved prenatal and postnatal genetic counseling.

A spatial genome aligner for resolving chromatin architectures from multiplexed DNA FISH

Multiplexed fluorescence in situ hybridization (FISH) is a widely used approach for analyzing three-dimensional genome organization, but it is challenging to derive chromosomal conformations from noisy fluorescence signals, and tracing chromatin is not straightforward. Here we report a spatial genome aligner that parses true chromatin signal from noise by aligning signals to a DNA polymer model. Using genomic distances separating imaged loci, our aligner estimates spatial distances expected to separate loci on a polymer in three-dimensional space. Our aligner then evaluates the physical probability observed signals belonging to these loci are connected, thereby tracing chromatin structures. We demonstrate that this spatial genome aligner can efficiently model chromosome architectures from DNA FISH data across multiple scales and be used to predict chromosome ploidies de novo in interphase cells. Reprocessing of previous whole-genome chromosome tracing data with this method indicates the spatial aggregation of sister chromatids in S/G2 phase cells in asynchronous mouse embryonic stem cells and provides evidence for extranumerary chromosomes that remain tightly paired in postmitotic neurons of the adult mouse cortex.

Integrative approaches to study enhancer-promoter communication

The spatiotemporal control of gene expression in complex multicellular organisms relies on noncoding regulatory sequences such as enhancers, which activate transcription of target genes often over large genomic distances. Despite the advances in the identification and characterization of enhancers, the principles and mechanisms by which enhancers select and control their target genes remain largely unknown. Here, we review recent interdisciplinary and quantitative approaches based on emerging techniques that aim to address open questions in the field, notably how regulatory information is encoded in the DNA sequence, how this information is transferred from enhancers to promoters, and how these processes are regulated in time.

Genome-wide tracing to decipher nuclear organization

Nuclear organization impacts gene expression activity and cell phenotype. Our current understanding is mainly derived from ensemble-level sequencing studies that reflect the 3D genome structure of millions of cells. These approaches have provided invaluable details on the 3D organizations of the genome and their relation to other nuclear landmarks. However, they mostly lack the ability to provide multimodal information simultaneously at the single-cell level. In recent years, cutting-edge imaging technologies have risen to the challenge of simultaneously describing multiple components of the nuclear space at the single-cell level, paving the way for a deeper understanding of the genome structure-function relationship. This review will focus on the development and utilization of such technologies to gain a multi-component view of the nucleus at single-cell resolution, dissecting the complexity and heterogeneity of nuclear organization.

Loop stacking organizes genome folding from TADs to chromosomes

Although population-level analyses revealed significant roles for CTCF and cohesin in mammalian genome organization, their contributions at the single-cell level remain incompletely understood. Here, we used a super-resolution microscopy approach to measure the effects of removal of CTCF or cohesin in mouse embryonic stem cells. Single-chromosome traces revealed cohesin-dependent loops, frequently stacked at their loop anchors forming multi-way contacts (hubs), bridging across TAD boundaries. Despite these bridging interactions, chromatin in intervening TADs was not intermixed, remaining separated in distinct loops around the hub. At the multi-TAD scale, steric effects from loop stacking insulated local chromatin from ultra-long range (>4 Mb) contacts. Upon cohesin removal, the chromosomes were more disordered and increased cell-cell variability in gene expression. Our data revise the TAD-centric understanding of CTCF and cohesin and provide a multi-scale, structural picture of how they organize the genome on the single-cell level through distinct contributions to loop stacking.

Competition between transcription and loop extrusion modulates promoter and enhancer dynamics

The spatiotemporal configuration of genes with distal regulatory elements, and the impact of chromatin mobility on transcription, remain unclear. Loop extrusion is an attractive model for bringing genetic elements together, but how this functionally interacts with transcription is also largely unknown. We combine live tracking of genomic loci and nascent transcripts with molecular dynamics simulations to assess the 4D arrangement of the Sox2 gene and its enhancer, in response to a battery of perturbations. We find that alterations in chromatin mobility, not promoter-enhancer distance, is more informative about transcriptional status. Active elements display more constrained mobility, consistent with confinement within specialized nuclear sites, and alterations in enhancer mobility distinguish poised from transcribing alleles. Strikingly, we find that whereas loop extrusion and transcription factor-mediated clustering contribute to promoter-enhancer proximity, they have antagonistic effects on chromatin dynamics. This provides an experimental framework for the underappreciated role of chromatin dynamics in genome regulation.

SATB1 regulates 3D genome architecture in T cells by constraining chromatin interactions surrounding CTCF-binding sites

Special AT-rich sequence binding protein 1 (SATB1) has long been proposed to act as a global chromatin loop organizer in T cells. However, the exact functions of SATB1 in spatial genome organization remain elusive. Here we show that the depletion of SATB1 in human and murine T cells leads to transcriptional dysregulation for genes involved in T cell activation, as well as alterations of 3D genome architecture at multiple levels, including compartments, topologically associating domains, and loops. Importantly, SATB1 extensively colocalizes with CTCF throughout the genome. Depletion of SATB1 leads to increased chromatin contacts among and across the SATB1/CTCF co-occupied sites, thereby affecting the transcription of critical regulators of T cell activation. The loss of SATB1 does not affect CTCF occupancy but significantly reduces the retention of CTCF in the nuclear matrix. Collectively, our data show that SATB1 contributes to 3D genome organization by constraining chromatin topology surrounding CTCF-binding sites.

Topologically associating domain boundaries are required for normal genome function

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

Synthetic regulatory genomics uncovers enhancer context dependence at the Sox2 locus

Sox2 expression in mouse embryonic stem cells (mESCs) depends on a distal cluster of DNase I hypersensitive sites (DHSs), but their individual contributions and degree of interdependence remain a mystery. We analyzed the endogenous Sox2 locus using Big-IN to scarlessly integrate large DNA payloads incorporating deletions, rearrangements, and inversions affecting single or multiple DHSs, as well as surgical alterations to transcription factor (TF) recognition sequences. Multiple mESC clones were derived for each payload, sequence-verified, and analyzed for Sox2 expression. We found that two DHSs comprising a handful of key TF recognition sequences were each sufficient for long-range activation of Sox2 expression. By contrast, three nearby DHSs were entirely context dependent, showing no activity alone but dramatically augmenting the activity of the autonomous DHSs. Our results highlight the role of context in modulating genomic regulatory element function, and our synthetic regulatory genomics approach provides a roadmap for the dissection of other genomic loci.

Prediction of CTCF loop anchor based on machine learning

Introduction: Various activities in biological cells are affected by three-dimensional genome structure. The insulators play an important role in the organization of higher-order structure. CTCF is a representative of mammalian insulators, which can produce barriers to prevent the continuous extrusion of chromatin loop. As a multifunctional protein, CTCF has tens of thousands of binding sites in the genome, but only a portion of them can be used as anchors of chromatin loops. It is still unclear how cells select the anchor in the process of chromatin looping.

Methods: In this paper, a comparative analysis is performed to investigate the sequence preference and binding strength of anchor and non-anchor CTCF binding sites. Furthermore, a machine learning model based on the CTCF binding intensity and DNA sequence is proposed to predict which CTCF sites can form chromatin loop anchors.

Results: The accuracy of the machine learning model that we constructed for predicting the anchor of the chromatin loop mediated by CTCF reached 0.8646. And we find that the formation of loop anchor is mainly influenced by the CTCF binding strength and binding pattern (which can be interpreted as the binding of different zinc fingers).

Discussion: In conclusion, our results suggest that The CTCF core motif and it’s flanking sequence may be responsible for the binding specificity. This work contributes to understanding the mechanism of loop anchor selection and provides a reference for the prediction of CTCF-mediated chromatin loops.

Recent progress and challenges in single-cell imaging of enhancer-promoter interaction

In the past two years, approaches relying on high-resolution microscopy and live-cell imaging have increasingly contributed to our understanding of the 3D genome organization and its importance for transcriptional control. Here, we describe recent progress that has highlighted how flexible and heterogeneous 3D chromatin structure is, on the length scales relevant to transcriptional control. We describe work that has investigated how robust transcriptional outcomes may be derived from such flexible organization without the need for clearly distinct structures in active and silent cells. We survey the latest state of the art in directly observing the dynamics of chromatin interactions, and suggest how some recent, apparently contradictory conclusions may be reconciled.