CRISPR and biochemical screens identify MAZ as a cofactor in CTCF-mediated insulation at Hox clusters

A continuum of zinc finger transcription factor retention on native chromatin underlies dynamic genome organization

Transcription factor (TF) residence on chromatin translates into quantitative transcriptional or structural outcomes on genome. Commonly used formaldehyde crosslinking fixes TF-DNA interactions cumulatively and compromises the measured occupancy level. Here we mapped the occupancy level of global or individual zinc finger TFs like CTCF and MAZ, in the form of highly resolved footprints, on native chromatin. By incorporating reinforcing perturbation conditions, we established S-score, a quantitative metric to proxy the continuum of CTCF or MAZ retention across different motifs on native chromatin. The native chromatin-retained CTCF sites harbor sequence features within CTCF motifs better explained by S-score than the metrics obtained from other crosslinking or native assays. CTCF retention on native chromatin correlates with local SUMOylation level, and anti-correlates with transcriptional activity. The S-score successfully delineates the otherwise-masked differential stability of chromatin structures mediated by CTCF, or by MAZ independent of CTCF. Overall, our study established a paradigm continuum of TF retention across binding sites on native chromatin, explaining the dynamic genome organization.

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.

Loop Catalog: a comprehensive HiChIP database of human and mouse samples

HiChIP enables cost-effective and high-resolution profiling of regulatory and structural loops. To leverage the increasing number of publicly available HiChIP datasets from diverse cell lines and primary cells, we developed the Loop Catalog (https://loopcatalog.lji.org), a web-based database featuring HiChIP loop calls for 1319 samples across 133 studies and 44 high-resolution Hi-C loop calls. We demonstrate its utility in interpreting fine-mapped GWAS variants (SNP-to-gene linking), in identifying enriched sequence motifs and motif pairs at loop anchors, and in network-level analysis of loops connecting regulatory elements (community detection). Our comprehensive catalog, spanning over 4M unique 5kb loops, along with the accompanying analysis modalities constitutes an important resource for studies in gene regulation and genome organization.

Xenopus tropicalis osteoblast-specific open chromatin regions reveal promoters and enhancers involved in human skeletal phenotypes and shed light on early vertebrate evolution

While understanding the genetic underpinnings of osteogenesis has far-reaching implications for skeletal diseases and evolution, a comprehensive characterization of the osteoblastic regulatory landscape in non-mammalian vertebrates is still lacking. Here, we compared the ATAC-Seq profile of Xenopus tropicalis (Xt) osteoblasts to a variety of non mineralizing control tissues, and identified osteoblast-specific nucleosome free regions (NFRs) at 527 promoters and 6747 distal regions. Sequence analyses, Gene Ontology, RNA-Seq and ChIP-Seq against four key histone marks confirmed that the distal regions correspond to bona fide osteogenic transcriptional enhancers exhibiting a shared regulatory logic with mammals. We report 425 regulatory regions conserved with human and globally associated to skeletogenic genes. Of these, 35 regions have been shown to impact human skeletal phenotypes by GWAS, including one trps1 enhancer and the runx2 promoter, two genes which are respectively involved in trichorhinophalangeal syndrome type I and cleidocranial dysplasia. Intriguingly, 60 osteoblastic NFRs also align to the genome of the elephant shark, a species lacking osteoblasts and bone tissue. To tackle this paradox, we chose to focus on dlx5 because its conserved promoter, known to integrate regulatory inputs during mammalian osteogenesis, harbours an osteoblast-specific NFR in both frog and human. Hence, we show that dlx5 is expressed in Xt and elephant shark odontoblasts, supporting a common cellular and genetic origin of bone and dentine. Taken together, our work (i) unravels the Xt osteogenic regulatory landscape, (ii) illustrates how cross-species comparisons harvest data relevant to human biology and (iii) reveals that a set of genes including bnc2, dlx5, ebf3, mir199a, nfia, runx2 and zfhx4 drove the development of a primitive form of mineralized skeletal tissue deep in the vertebrate lineage.

The scales, mechanisms, and dynamics of the genome architecture

Even when split into several chromosomes, DNA molecules that make up our genome are too long to fit into the cell nuclei unless massively folded. Such folding must accommodate the need for timely access to selected parts of the genome by transcription factors, RNA polymerases, and DNA replication machinery. Here, we review our current understanding of the genome folding inside the interphase nuclei. We consider the resulting genome architecture at three scales with a particular focus on the intermediate (meso) scale and summarize the insights gained from recent experimental observations and diverse computational models.

The N-terminal dimerization domains of human and Drosophila CTCF have similar functionality

BACKGROUND

CTCF is highly likely to be the ancestor of proteins that contain large clusters of C2H2 zinc finger domains, and its conservation is observed across most bilaterian organisms. In mammals, CTCF is the primary architectural protein involved in organizing chromosome topology and mediating enhancer-promoter interactions over long distances. In Drosophila, CTCF (dCTCF) cooperates with other architectural proteins to establish long-range interactions and chromatin boundaries. CTCFs of various organisms contain an unstructured N-terminal dimerization domain (DD) and clusters comprising eleven zinc-finger domains of the C2H2 type. The Drosophila (dCTCF) and human (hCTCF) CTCFs share sequence homology in only five C2H2 domains that specifically bind to a conserved 15 bp motif.

RESULTS

Previously, we demonstrated that CTCFs from different organisms carry unstructured N-terminal dimerization domains (DDs) that lack sequence homology. Here we used the CTCFattP(mCh) platform to introduce desired changes in the Drosophila CTCF gene and generated a series of transgenic lines expressing dCTCF with different variants of the N-terminal domain. Our findings revealed that the functionality of dCTCF is significantly affected by the deletion of the N-terminal DD. Additionally, we observed a strong impact on the binding of the dCTCF mutant to chromatin upon deletion of the DD. However, chromatin binding was restored in transgenic flies expressing a chimeric CTCF protein with the DD of hCTCF. Although the chimeric protein exhibited lower expression levels than those of the dCTCF variants, it efficiently bound to chromatin similarly to the wild type (wt) protein.

CONCLUSIONS

Our findings suggest that one of the evolutionarily conserved functions of the unstructured N-terminal dimerization domain is to recruit dCTCF to its genomic sites in vivo.

A CTCF-dependent mechanism underlies the Hox timer: relation to a segmented body plan

During gastrulation, Hox genes are activated in a time-sequence that follows the order of the genes along their clusters. This property, which is observed in all animals that develop following a progressive rostral-to-caudal morphogenesis, is associated with changes in the chromatin structure and epigenetic profiles of Hox clusters, suggesting a process at least partly based on sequential gene accessibility. Here, we discuss recent work on this issue, as well as a possible mechanism based on the surprising conservation in both the distribution and orientation of CTCF sites inside vertebrate Hox clusters.

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.

Establishing and maintaining Hox profiles during spinal cord development

The chromosomally-arrayed Hox gene family plays central roles in embryonic patterning and the specification of cell identities throughout the animal kingdom. In vertebrates, the relatively large number of Hox genes and pervasive expression throughout the body has hindered understanding of their biological roles during differentiation. Studies on the subtype diversification of spinal motor neurons (MNs) have provided a tractable system to explore the function of Hox genes during differentiation, and have provided an entry point to explore how neuronal fate determinants contribute to motor circuit assembly. Recent work, using both in vitro and in vivo models of MN subtype differentiation, have revealed how patterning morphogens and regulation of chromatin structure determine cell-type specific programs of gene expression. These studies have not only shed light on basic mechanisms of rostrocaudal patterning in vertebrates, but also have illuminated mechanistic principles of gene regulation that likely operate in the development and maintenance of terminal fates in other systems.

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.

Multiple Roles of dXNP and dADD1-Drosophila Orthologs of ATRX Chromatin Remodeler

The Drosophila melanogaster dADD1 and dXNP proteins are orthologues of the ADD and SNF2 domains of the vertebrate ATRX (Alpha-Thalassemia with mental Retardation X-related) protein. ATRX plays a role in general molecular processes, such as regulating chromatin status and gene expression, while dADD1 and dXNP have similar functions in the Drosophila genome. Both ATRX and dADD1/dXNP interact with various protein partners and participate in various regulatory complexes. Disruption of ATRX expression in humans leads to the development of α-thalassemia and cancer, especially glioma. However, the mechanisms that allow ATRX to regulate various cellular processes are poorly understood. Studying the functioning of dADD1/dXNP in the Drosophila model may contribute to understanding the mechanisms underlying the multifunctional action of ATRX and its connection with various cellular processes. This review provides a brief overview of the currently available information in mammals and Drosophila regarding the roles of ATRX, dXNP, and dADD1. It discusses possible mechanisms of action of complexes involving these proteins.

Chromatin in 3D distinguishes dMes-4/NSD and Hypb/dSet2 in protecting genes from H3K27me3 silencing

Cell type-specific barcoding of genomes requires the establishment of hundreds of heterochromatin domains where heterochromatin-associated repressive complexes hinder chromatin accessibility thereby silencing genes. At heterochromatin-euchromatin borders, regulation of accessibility not only depends on the delimitation of heterochromatin but may also involve interplays with nearby genes and their transcriptional activity, or alternatively on histone modifiers, chromatin barrier insulators, and more global demarcation of chromosomes into 3D compartmentalized domains and topological-associating domain (TADs). Here, we show that depletion of H3K36 di- or tri-methyl histone methyltransferases dMes-4/NSD or Hypb/dSet2 induces reproducible increasing levels of H3K27me3 at heterochromatin borders including in nearby promoters, thereby repressing hundreds of genes. Furthermore, dMes-4/NSD influences genes demarcated by insulators and TAD borders, within chromatin hubs, unlike transcription-coupled action of Hypb/dSet2 that protects genes independently of TADs. Insulator mutants recapitulate the increase of H3K27me3 upon dMes-4/NSD depletion unlike Hypb/dSet2. Hi-C data demonstrate how dMes-4/NSD blocks propagation of long-range interactions onto active regions. Our data highlight distinct mechanisms protecting genes from H3K27me3 silencing, highlighting a direct influence of H3K36me on repressive TADs.

The CLAMP GA-binding transcription factor regulates heat stress-induced transcriptional repression by associating with 3D loop anchors

In order to survive when exposed to heat stress (HS), organisms activate stress response genes and repress constitutive gene expression to prevent the accumulation of potentially toxic RNA and protein products. Although many studies have elucidated the mechanisms that drive HS-induced activation of stress response genes across species, little is known about repression mechanisms or how genes are targeted for activation versus repression context-specifically. The mechanisms of heat stress-regulated activation have been well-studied in Drosophila, in which the GA-binding transcription factor GAF is important for activating genes upon heat stress. Here, we show that a functionally distinct GA-binding transcription factor (TF) protein, CLAMP (Chromatin-linked adaptor for MSL complex proteins), is essential for repressing constitutive genes upon heat stress but not activation of the canonical heat stress pathway. HS induces loss of CLAMP-associated 3D chromatin loop anchors associated with different combinations of GA-binding TFs prior to HS if a gene becomes repressed versus activated. Overall, we demonstrate that CLAMP promotes repression of constitutive genes upon HS, and repression and activation are associated with the loss of CLAMP-associated 3D chromatin loops bound by different combinations of GA-binding TFs.

Regulation of the three-dimensional chromatin organization by transposable elements in pig spleen

Like other mammalian species, the pig genome is abundant with transposable elements (TEs). The importance of TEs for three-dimensional (3D) chromatin organization has been observed in species like human and mouse, yet current understanding about pig TEs is absent. Here, we investigated the contribution of TEs for the 3D chromatin organization in three pig tissues, focusing on spleen which is crucial for both adaptive and innate immunity. We identified dozens of TE families overrepresented with CTCF binding sites, including LTR22_SS, LTR15_SS and LTR16_SSc which are pig-specific families of endogenous retroviruses (ERVs). Interestingly, LTR22_SS elements harbor a CTCF motif and create hundreds of CTCF binding sites that are associated with adaptive immunity. We further applied Hi-C to profile the 3D chromatin structure in spleen and found that TE-derived CTCF binding sites correlate with chromatin insulation and frequently overlap TAD borders and loop anchors. Notably, one LTR22_SS-derived CTCF binding site demarcate a TAD boundary upstream of XCL1, which is a spleen-enriched chemokine gene important for lymphocyte trafficking and inflammation. Overall, this study represents a first step toward understanding the function of TEs on 3D chromatin organization regulation in pigs and expands our understanding about the functional importance of TEs in mammals.

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.

Cell-type-specific prediction of 3D chromatin organization enables high-throughput in silico genetic screening

Investigating how chromatin organization determines cell-type-specific gene expression remains challenging. Experimental methods for measuring three-dimensional chromatin organization, such as Hi-C, are costly and have technical limitations, restricting their broad application particularly in high-throughput genetic perturbations. We present C.Origami, a multimodal deep neural network that performs de novo prediction of cell-type-specific chromatin organization using DNA sequence and two cell-type-specific genomic features-CTCF binding and chromatin accessibility. C.Origami enables in silico experiments to examine the impact of genetic changes on chromatin interactions. We further developed an in silico genetic screening approach to assess how individual DNA elements may contribute to chromatin organization and to identify putative cell-type-specific trans-acting regulators that collectively determine chromatin architecture. Applying this approach to leukemia cells and normal T cells, we demonstrate that cell-type-specific in silico genetic screening, enabled by C.Origami, can be used to systematically discover novel chromatin regulation circuits in both normal and disease-related biological systems.

MethNet: a robust approach to identify regulatory hubs and their distal targets in cancer

Aberrations in the capacity of DNA/chromatin modifiers and transcription factors to bind non-coding regions can lead to changes in gene regulation and impact disease phenotypes. However, identifying distal regulatory elements and connecting them with their target genes remains challenging. Here, we present MethNet, a pipeline that integrates large-scale DNA methylation and gene expression data across multiple cancers, to uncover novel cis regulatory elements (CREs) in a 1Mb region around every promoter in the genome. MethNet identifies clusters of highly ranked CREs, referred to as 'hubs', which contribute to the regulation of multiple genes and significantly affect patient survival. Promoter-capture Hi-C confirmed that highly ranked associations involve physical interactions between CREs and their gene targets, and CRISPRi based scRNA Perturb-seq validated the functional impact of CREs. Thus, MethNet-identified CREs represent a valuable resource for unraveling complex mechanisms underlying gene expression, and for prioritizing the verification of predicted non-coding disease hotspots.

Multidimensional profiling reveals GATA1-modulated stage-specific chromatin states and functional associations during human erythropoiesis

Mammalian erythroid development can be divided into three stages: hematopoietic stem and progenitor cell (HSPC), erythroid progenitor (Ery-Pro), and erythroid precursor (Ery-Pre). However, the mechanisms by which the 3D genome changes to establish the stage-specific transcription programs that are critical for erythropoiesis remain unclear. Here, we analyze the chromatin landscape at multiple levels in defined populations from primary human erythroid culture. While compartments and topologically associating domains remain largely unchanged, ∼50% of H3K27Ac-marked enhancers are dynamic in HSPC versus Ery-Pre. The enhancer anchors of enhancer-promoter loops are enriched for occupancy of respective stage-specific transcription factors (TFs), indicating these TFs orchestrate the enhancer connectome rewiring. The master TF of erythropoiesis, GATA1, is found to occupy most erythroid gene promoters at the Ery-Pro stage, and mediate conspicuous local rewiring through acquiring binding at the distal regions in Ery-Pre, promoting productive erythroid transcription output. Knocking out GATA1 binding sites precisely abrogates local rewiring and corresponding gene expression. Interestingly, knocking down GATA1 can transiently revert the cell state to an earlier stage and prolong the window of progenitor state. This study reveals mechanistic insights underlying chromatin rearrangements during development by integrating multidimensional chromatin landscape analyses to associate with transcription output and cellular states.

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

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

Probabilistic establishment of speckle-associated inter-chromosomal interactions

Inter-chromosomal interactions play a crucial role in genome organization, yet the organizational principles remain elusive. Here, we introduce a novel computational method to systematically characterize inter-chromosomal interactions using in situ Hi-C results from various cell types. Our method successfully identifies two apparently hub-like inter-chromosomal contacts associated with nuclear speckles and nucleoli, respectively. Interestingly, we discover that nuclear speckle-associated inter-chromosomal interactions are highly cell-type invariant with a marked enrichment of cell-type common super-enhancers (CSEs). Validation using DNA Oligopaint fluorescence in situ hybridization (FISH) shows a strong but probabilistic interaction behavior between nuclear speckles and CSE-harboring genomic regions. Strikingly, we find that the likelihood of speckle-CSE associations can accurately predict two experimentally measured inter-chromosomal contacts from Hi-C and Oligopaint DNA FISH. Our probabilistic establishment model well describes the hub-like structure observed at the population level as a cumulative effect of summing individual stochastic chromatin-speckle interactions. Lastly, we observe that CSEs are highly co-occupied by MAZ binding and MAZ depletion leads to significant disorganization of speckle-associated inter-chromosomal contacts. Taken together, our results propose a simple organizational principle of inter-chromosomal interactions mediated by MAZ-occupied CSEs.

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

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

Partial erosion on under-methylated regions and chromatin reprogramming contribute to oncogene activation in IDH mutant gliomas

BACKGROUND

IDH1/2 hotspot mutations are well known to drive oncogenic mutations in gliomas and are well-defined in the WHO 2021 classification of central nervous system tumors. Specifically, IDH mutations lead to aberrant hypermethylation of under-methylated regions (UMRs) in normal tissues through the disruption of TET enzymes. However, the chromatin reprogramming and transcriptional changes induced by IDH-related hypermethylation in gliomas remain unclear.

RESULTS

Here, we have developed a precise computational framework based on Hidden Markov Model to identify altered methylation states of UMRs at single-base resolution. By applying this framework to whole-genome bisulfite sequencing data from 75 normal brain tissues and 15 IDH mutant glioma tissues, we identified two distinct types of hypermethylated UMRs in IDH mutant gliomas. We named them partially hypermethylated UMRs (phUMRs) and fully hypermethylated UMRs (fhUMRs), respectively. We found that the phUMRs and fhUMRs exhibit distinct genomic features and chromatin states. Genes related to fhUMRs were more likely to be repressed in IDH mutant gliomas. In contrast, genes related to phUMRs were prone to be up-regulated in IDH mutant gliomas. Such activation of phUMR genes is associated with the accumulation of active H3K4me3 and the loss of H3K27me3, as well as H3K36me3 accumulation in gene bodies to maintain gene expression stability. In summary, partial erosion on UMRs was accompanied by locus-specific changes in key chromatin marks, which may contribute to oncogene activation.

CONCLUSIONS

Our study provides a computational strategy for precise decoding of methylation encroachment patterns in IDH mutant gliomas, revealing potential mechanistic insights into chromatin reprogramming that contribute to oncogenesis.

MYC reshapes CTCF-mediated chromatin architecture in prostate cancer

MYC is a well characterized oncogenic transcription factor in prostate cancer, and CTCF is the main architectural protein of three-dimensional genome organization. However, the functional link between the two master regulators has not been reported. In this study, we find that MYC rewires prostate cancer chromatin architecture by interacting with CTCF protein. Through combining the H3K27ac, AR and CTCF HiChIP profiles with CRISPR deletion of a CTCF site upstream of MYC gene, we show that MYC activation leads to profound changes of CTCF-mediated chromatin looping. Mechanistically, MYC colocalizes with CTCF at a subset of genomic sites, and enhances CTCF occupancy at these loci. Consequently, the CTCF-mediated chromatin looping is potentiated by MYC activation, resulting in the disruption of enhancer-promoter looping at neuroendocrine lineage plasticity genes. Collectively, our findings define the function of MYC as a CTCF co-factor in three-dimensional genome organization.

Mechanisms of Interaction between Enhancers and Promoters in Three Drosophila Model Systems

In higher eukaryotes, the regulation of developmental gene expression is determined by enhancers, which are often located at a large distance from the promoters they regulate. Therefore, the architecture of chromosomes and the mechanisms that determine the functional interaction between enhancers and promoters are of decisive importance in the development of organisms. Mammals and the model animal Drosophila have homologous key architectural proteins and similar mechanisms in the organization of chromosome architecture. This review describes the current progress in understanding the mechanisms of the formation and regulation of long-range interactions between enhancers and promoters at three well-studied key regulatory loci in Drosophila.

Dynamic chromatin organization and regulatory interactions in human endothelial cell differentiation

Vascular endothelial cells are a mesoderm-derived lineage with many essential functions, including angiogenesis and coagulation. The gene-regulatory mechanisms underpinning endothelial specialization are largely unknown, as are the roles of chromatin organization in regulating endothelial cell transcription. To investigate the relationships between chromatin organization and gene expression, we induced endothelial cell differentiation from human pluripotent stem cells and performed Hi-C and RNA-sequencing assays at specific time points. Long-range intrachromosomal contacts increase over the course of differentiation, accompanied by widespread heteroeuchromatic compartment transitions that are tightly associated with transcription. Dynamic topologically associating domain boundaries strengthen and converge on an endothelial cell state, and function to regulate gene expression. Chromatin pairwise point interactions (DNA loops) increase in frequency during differentiation and are linked to the expression of genes essential to vascular biology. Chromatin dynamics guide transcription in endothelial cell development and promote the divergence of endothelial cells from cardiomyocytes.

Mechanisms of enhancer-promoter communication and chromosomal architecture in mammals and Drosophila

The spatial organization of chromosomes is involved in regulating the majority of intranuclear processes in higher eukaryotes, including gene expression. Drosophila was used as a model to discover many transcription factors whose homologs play a key role in regulation of gene expression in mammals. According to modern views, a cohesin complex mostly determines the architecture of mammalian chromosomes by forming chromatin loops on anchors created by the CTCF DNA-binding architectural protein. The role of the cohesin complex in chromosome architecture is poorly understood in Drosophila, and CTCF is merely one of many Drosophila architectural proteins with a proven potential to organize specific long-range interactions between regulatory elements in the genome. The review compares the mechanisms responsible for long-range interactions and chromosome architecture between mammals and Drosophila.

Factors and Mechanisms That Influence Chromatin-Mediated Enhancer–Promoter Interactions and Transcriptional Regulation

Simple Summary

The physical interactions between enhancers and promoters create chromatin conformations involved in gene regulation. In cancer cells, the chromatin conformations can be altered with uncontrolled deposition of histone marks resulting in varied gene expression. Although it is not entirely comprehensive how chromatin-mediated enhancer–promoter (E–P) interactions with various histone marks can affect gene expression, this proximity has been observed in multiple systems at multiple loci and is thought to be essential to control gene expression. In this review, we focus on emerging views of chromatin conformations associated with the E–P interactions and factors that establish or maintain such interactions, which may regulate gene expression.

Abstract

Eukaryotic gene expression is regulated through chromatin conformation, in which enhancers and promoters physically interact (E–P interactions). How such chromatin-mediated E–P interactions affect gene expression is not yet fully understood, but the roles of histone acetylation and methylation, pioneer transcription factors, and architectural proteins such as CCCTC binding factor (CTCF) and cohesin have recently attracted attention. Moreover, accumulated data suggest that E–P interactions are mechanistically involved in biophysical events, including liquid–liquid phase separation, and in biological events, including cancers. In this review, we discuss various mechanisms that regulate eukaryotic gene expression, focusing on emerging views regarding chromatin conformations that are involved in E–P interactions and factors that establish and maintain them.

ExsgRNA: reduce off-target efficiency by on-target mismatched sgRNA

Clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 gene editing technology has been widely used to facilitate efficient genome editing. Current popular sgRNA design tools only consider the sgRNA perfectly matched to the target site and provide the results without any on-target mismatch. We suppose taking on-target gRNA-DNA mismatches into consideration might provide better sgRNA with similar binding activity and reduced off-target sites. Here, we trained a seq2seq-attention model with feedback-loop architecture, to automatically generate sgRNAs with on-target mismatches. Dual-luciferase reporter experiment showed that multiple sgRNAs with three mismatches could achieve the 80% of the relative activity of the perfect matched sgRNA. Meanwhile, it could reduce the number of off-target sites using sgRNAs with on-target mismatches. Finally, we provided a freely accessible web server sgRNA design tool named ExsgRNA. Users could submit their target sequence to this server and get optimal sgRNAs with less off-targets and similar on-target activity compared with the perfect-matched sgRNA.