Acute depletion of CTCF directly affects MYC regulation through loss of enhancer-promoter looping

DeepOCR: A multi-species deep-learning framework for accurate identification of open chromatin regions in livestock

A wealth of experimental evidence has suggested that open chromatin regions (OCRs) are involved in many critical biological activities, such as DNA replication, enhancer activity, and gene transcription. Accurately identifying OCRs in livestock species can provide critical insights into the distribution and characteristics of OCRs for disease treatment in livestock, thereby improving animal welfare. However, most current machine-learning methods for OCR prediction were originally designed for a limited number of model organisms, such as humans and some model organisms, and thus their performance on non-model organisms, specifically livestock, is often unsatisfactory. To bridge this gap, we propose DeepOCR, a lightweight depth-separable residual network model for predicting OCRs in livestock, including chicken, cattle, and sheep. DeepOCR integrates a single convolution layer and two improved residue structure blocks to extract and learn important features from the input DNA sequences. A fully connected layer was also employed to further process the extracted features and improve the robustness of the entire network. Our benchmarking experiments demonstrated superior prediction performance of DeepOCR compared to state-of-the-art approaches on testing datasets of the three species. The source code of DeepOCR is freely available for academic purposes at https://github.com/jasonzhao371/DeepOCR/. We anticipate DeepOCR servers as a practical and reliable computational tool for OCR-related studies in livestock species.

Tet-mediated DNA methylation dynamics affect chromosome organization

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

Chromatin and aberrant enhancer activity in KMT2A rearranged acute lymphoblastic leukemia

To make a multicellular organism, genes need to be transcribed at the right developmental stages and in the right tissues. DNA sequences termed 'enhancers' are crucial to achieve this. Despite concerted efforts, the exact mechanisms of enhancer activity remain elusive. Mixed lineage leukemia (MLL or KMT2A) rearrangements (MLLr), commonly observed in cases of acute lymphoblastic leukemia (ALL) and acute myeloid leukemia, produce novel in-frame fusion proteins. Recent work has shown that the MLL-AF4 fusion protein drives aberrant enhancer activity at key oncogenes in ALL, dependent on the continued presence of MLL-AF4 complex components. As well as providing some general insights into enhancer function, these observations may also provide an explanation for transcriptional heterogeneity observed in MLLr patients.

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.

Enhancer selectivity in space and time: from enhancer-promoter interactions to promoter activation

The primary regulators of metazoan gene expression are enhancers, originally functionally defined as DNA sequences that can activate transcription at promoters in an orientation-independent and distance-independent manner. Despite being crucial for gene regulation in animals, what mechanisms underlie enhancer selectivity for promoters, and more fundamentally, how enhancers interact with promoters and activate transcription, remain poorly understood. In this Review, we first discuss current models of enhancer-promoter interactions in space and time and how enhancers affect transcription activation. Next, we discuss different mechanisms that mediate enhancer selectivity, including repression, biochemical compatibility and regulation of 3D genome structure. Through 3D polymer simulations, we illustrate how the ability of 3D genome folding mechanisms to mediate enhancer selectivity strongly varies for different enhancer-promoter interaction mechanisms. Finally, we discuss how recent technical advances may provide new insights into mechanisms of enhancer-promoter interactions and how technical biases in methods such as Hi-C and Micro-C and imaging techniques may affect their interpretation.

Spatial chromatin accessibility sequencing resolves high-order spatial interactions of epigenomic markers

As the genome is organized into a three-dimensional structure in intracellular space, epigenomic information also has a complex spatial arrangement. However, most epigenetic studies describe locations of methylation marks, chromatin accessibility regions, and histone modifications in the horizontal dimension. Proper spatial epigenomic information has rarely been obtained. In this study, we designed spatial chromatin accessibility sequencing (SCA-seq) to resolve the genome conformation by capturing the epigenetic information in single-molecular resolution while simultaneously resolving the genome conformation. Using SCA-seq, we are able to examine the spatial interaction of chromatin accessibility (e.g. enhancer-promoter contacts), CpG island methylation, and spatial insulating functions of the CCCTC-binding factor. We demonstrate that SCA-seq paves the way to explore the mechanism of epigenetic interactions and extends our knowledge in 3D packaging of DNA in the nucleus.

Chuanxiong Rhizoma extracts prevent liver fibrosis via targeting CTCF-c-MYC-H19 pathway

Objective

Hepatic fibrosis has been widely considered as a conjoint consequence of almost all chronic liver diseases. Chuanxiong Rhizoma (Chuanxiong in Chinese, CX) is a traditional Chinese herbal product to prevent cerebrovascular, gynecologic and hepatic diseases. Our previous study found that CX extracts significantly reduced collagen contraction force of hepatic stellate cells (HSCs). Here, this study aimed to compare the protection of different CX extracts on bile duct ligation (BDL)-induced liver fibrosis and investigate plausible underlying mechanisms.

Methods

The active compounds of CX extracts were identified by high performance liquid chromatography (HPLC). Network pharmacology was used to determine potential targets of CX against hepatic fibrosis. Bile duct hyperplasia and liver fibrosis were evaluated by serologic testing and histopathological evaluation. The expression of targets of interest was determined by quantitative real-time PCR (qPCR) and Western blot.

Results

Different CX extracts were identified by tetramethylpyrazine, ferulic acid and senkyunolide A. Based on the network pharmacological analysis, 42 overlap targets were obtained via merging the candidates targets of CX and liver fibrosis. Different aqueous, alkaloid and phthalide extracts of CX (CXAE, CXAL and CXPHL) significantly inhibited diffuse severe bile duct hyperplasia and thus suppressed hepatic fibrosis by decreasing CCCTC binding factor (CTCF)-c-MYC-long non-coding RNA H19 (H19) pathway in the BDL-induced mouse model. Meanwhile, CX extracts, especially CXAL and CXPHL also suppressed CTCF-c-MYC-H19 pathway and inhibited ductular reaction in cholangiocytes stimulated with taurocholate acid (TCA), lithocholic acid (LCA) and transforming growth factor beta (TGF-β), as illustrated by decreased bile duct proliferation markers.

Conclusion

Our data supported that different CX extracts, especially CXAL and CXPHL significantly alleviated hepatic fibrosis and bile duct hyperplasia via inhibiting CTCF-c-MYC-H19 pathway, providing novel insights into the anti-fibrotic mechanism of CX.

An alternative NURF complex sustains acute myeloid leukemia by regulating the accessibility of insulator regions

Efficient treatment of acute myeloid leukemia (AML) patients remains a challenge despite recent therapeutic advances. Here, using a CRISPRi screen targeting chromatin factors, we identified the nucleosome-remodeling factor (NURF) subunit BPTF as an essential regulator of AML cell survival. We demonstrate that BPTF forms an alternative NURF chromatin remodeling complex with SMARCA5 and BAP18, which regulates the accessibility of a large set of insulator regions in leukemic cells. This ensures efficient CTCF binding and boundary formation between topologically associated domains that is essential for maintaining the leukemic transcriptional programs. We also demonstrate that the well-studied PHD2-BROMO chromatin reader domains of BPTF, while contributing to complex recruitment to chromatin, are dispensable for leukemic cell growth. Taken together, our results uncover how the alternative NURF complex contributes to leukemia and provide a rationale for its targeting in AML.

Systematic characterization of the HOXA9 downstream targets in MLL-r leukemia by noncoding CRISPR screens

Accumulating evidence indicates that HOXA9 dysregulation is necessary and sufficient for leukemic transformation and maintenance. However, it remains largely unknown how HOXA9, as a homeobox transcriptional factor, binds to noncoding regulatory sequences and controls the downstream genes. Here, we conduct dropout CRISPR screens against 229 HOXA9-bound peaks identified by ChIP-seq. Integrative data analysis identifies reproducible noncoding hits, including those located in the distal enhancer of FLT3 and intron of CDK6. The Cas9-editing and dCas9-KRAB silencing of the HOXA9-bound sites significantly reduce corresponding gene transcription and impair cell proliferation in vitro, and in vivo by transplantation into NSG female mice. In addition, RNA-seq, Q-PCR analysis, chromatin accessibility change, and chromatin conformation evaluation uncover the noncoding regulation mechanism of HOXA9 and its functional downstream genes. In summary, our work improves our understanding of how HOXA9-associated transcription programs reconstruct the regulatory network specifying MLL-r dependency.

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

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

CTCF/cohesin organize the ground state of chromatin-nuclear speckle association

The interchromatin space in the cell nucleus contains various membrane-less nuclear bodies. Recent findings indicate that nuclear speckles, comprising a distinct nuclear body, exhibit interactions with certain chromatin regions in a ground state. Key questions are how this ground state of chromatin-nuclear speckle association is established and what are the gene regulatory roles of this layer of nuclear organization. We report here that chromatin structural factors CTCF and cohesin are required for full ground state association between DNA and nuclear speckles. Disruption of ground state DNA-speckle contacts via either CTCF depletion or cohesin depletion had minor effects on basal level expression of speckle-associated genes, however we show strong negative effects on stimulus-dependent induction of speckle-associated genes. We identified a putative speckle targeting motif (STM) within cohesin subunit RAD21 and demonstrated that the STM is required for chromatin-nuclear speckle association. In contrast to reduction of CTCF or RAD21, depletion of the cohesin releasing factor WAPL stabilized cohesin on chromatin and DNA-speckle contacts, resulting in enhanced inducibility of speckle-associated genes. In addition, we observed disruption of chromatin-nuclear speckle association in patient derived cells with Cornelia de Lange syndrome (CdLS), a congenital neurodevelopmental diagnosis involving defective cohesin pathways, thus revealing nuclear speckles as an avenue for therapeutic inquiry. In summary, our findings reveal a mechanism to establish the ground organizational state of chromatin-speckle association, to promote gene inducibility, and with relevance to human disease.

A truncated ETHYLENE INSENSITIVE3-like protein, GhLYI, regulates senescence in cotton

Numerous endogenous and environmental signals regulate the intricate and highly orchestrated process of plant senescence. Ethylene (ET), which accumulates as senescence progresses, is a major promoter of leaf senescence. The master transcription activator ETHYLENE INSENSITIVE3 (EIN3) activates the expression of a wide range of downstream genes during leaf senescence. Here, we found that a unique EIN3-LIKE 1 (EIL1) gene, cotton LINT YIELD INCREASING (GhLYI), encodes a truncated EIN3 protein in upland cotton (Gossypium hirsutum L.) that functions as an ET signal response factor and a positive regulator of senescence. Ectopic expression or overexpression of GhLYI accelerated leaf senescence in both Arabidopsis (Arabidopsis thaliana) and cotton. Cleavage under targets and tagmentation (CUT&Tag) analyses revealed that SENESCENCE-ASSOCIATED GENE 20 (SAG20) was a target of GhLYI. Electrophoretic mobility shift assay (EMSA), yeast 1-hybrid (Y1H), and dual-luciferase transient expression assay confirmed that GhLYI directly bound the promoter of SAG20 to activate its expression. Transcriptome analysis revealed that transcript levels of a series of senescence-related genes, SAG12, NAC-LIKE, ACTIVATED by APETALA 3/PISTILLATA (NAP/ANAC029), and WRKY53, are substantially induced in GhLYI overexpression plants compared with wild-type (WT) plants. Virus-induced gene silencing (VIGS) preliminarily confirmed that knockdown of GhSAG20 delayed leaf senescence. Collectively, our findings provide a regulatory module involving GhLYI-GhSAG20 in controlling senescence in cotton.

Allelic chromatin structure precedes imprinted expression of Kcnk9 during neurogenesis

Differences in chromatin state inherited from the parental gametes influence the regulation of maternal and paternal alleles in offspring. This phenomenon, known as genomic imprinting, results in genes preferentially transcribed from one parental allele. While local epigenetic factors such as DNA methylation are known to be important for the establishment of imprinted gene expression, less is known about the mechanisms by which differentially methylated regions (DMRs) lead to differences in allelic expression across broad stretches of chromatin. Allele-specific higher-order chromatin structure has been observed at multiple imprinted loci, consistent with the observation of allelic binding of the chromatin-organizing factor CTCF at multiple DMRs. However, whether allelic chromatin structure impacts allelic gene expression is not known for most imprinted loci. Here we characterize the mechanisms underlying brain-specific imprinted expression of the Peg13-Kcnk9 locus, an imprinted region associated with intellectual disability. We performed region capture Hi-C on mouse brains from reciprocal hybrid crosses and found imprinted higher-order chromatin structure caused by the allelic binding of CTCF to the Peg13 DMR. Using an in vitro neuron differentiation system, we showed that imprinted chromatin structure precedes imprinted expression at the locus. Additionally, activation of a distal enhancer induced imprinted expression of Kcnk9 in an allelic chromatin structure-dependent manner. This work provides a high-resolution map of imprinted chromatin structure and demonstrates that chromatin state established in early development can promote imprinted expression upon differentiation.

MLL-AF4 cooperates with PAF1 and FACT to drive high-density enhancer interactions in leukemia

Aberrant enhancer activation is a key mechanism driving oncogene expression in many cancers. While much is known about the regulation of larger chromosome domains in eukaryotes, the details of enhancer-promoter interactions remain poorly understood. Recent work suggests co-activators like BRD4 and Mediator have little impact on enhancer-promoter interactions. In leukemias controlled by the MLL-AF4 fusion protein, we use the ultra-high resolution technique Micro-Capture-C (MCC) to show that MLL-AF4 binding promotes broad, high-density regions of enhancer-promoter interactions at a subset of key targets. These enhancers are enriched for transcription elongation factors like PAF1C and FACT, and the loss of these factors abolishes enhancer-promoter contact. This work not only provides an additional model for how MLL-AF4 is able to drive high levels of transcription at key genes in leukemia but also suggests a more general model linking enhancer-promoter crosstalk and transcription elongation.

Patterns of transcription factor binding and epigenome at promoters allow interpretable predictability of multiple functions of non-coding and coding genes

Understanding the biological roles of all genes only through experimental methods is challenging. A computational approach with reliable interpretability is needed to infer the function of genes, particularly for non-coding RNAs. We have analyzed genomic features that are present across both coding and non-coding genes like transcription factor (TF) and cofactor ChIP-seq (823), histone modifications ChIP-seq (n = 621), cap analysis gene expression (CAGE) tags (n = 255), and DNase hypersensitivity profiles (n = 255) to predict ontology-based functions of genes. Our approach for gene function prediction was reliable (>90% balanced accuracy) for 486 gene-sets. PubMed abstract mining and CRISPR screens supported the inferred association of genes with biological functions, for which our method had high accuracy. Further analysis revealed that TF-binding patterns at promoters have high predictive strength for multiple functions. TF-binding patterns at the promoter add an unexplored dimension of explainable regulatory aspects of genes and their functions. Therefore, we performed a comprehensive analysis for the functional-specificity of TF-binding patterns at promoters and used them for clustering functions to reveal many latent groups of gene-sets involved in common major cellular processes. We also showed how our approach could be used to infer the functions of non-coding genes using the CRISPR screens of coding genes, which were validated using a long non-coding RNA CRISPR screen. Thus our results demonstrated the generality of our approach by using gene-sets from CRISPR screens. Overall, our approach opens an avenue for predicting the involvement of non-coding genes in various functions.

Dynamic nucleosome remodeling mediated by YY1 underlies early mouse development

Nucleosome positioning can alter the accessibility of DNA-binding proteins to their cognate DNA elements, and thus its precise control is essential for cell identity and function. Mammalian preimplantation embryos undergo temporal changes in gene expression and cell potency, suggesting the involvement of dynamic epigenetic control during this developmental phase. However, the dynamics of nucleosome organization during early development are poorly understood. In this study, using a low-input MNase-seq method, we show that nucleosome positioning is globally obscure in zygotes but becomes well defined during subsequent development. Down-regulation of the chromatin assembly in embryonic stem cells can partially reverse nucleosome organization into a zygote-like pattern, suggesting a possible link between the chromatin assembly pathway and fuzzy nucleosomes in zygotes. We also reveal that YY1, a zinc finger-containing transcription factor expressed upon zygotic genome activation, regulates the de novo formation of well-positioned nucleosome arrays at the regulatory elements through identifying YY1-binding sites in eight-cell embryos. The YY1-binding regions acquire H3K27ac enrichment around the eight-cell and morula stages, and YY1 depletion impairs the morula-to-blastocyst transition. Thus, our study delineates the remodeling of nucleosome organization and its underlying mechanism during early mouse development.

Allelic chromatin structure primes imprinted expression of Kcnk9 during neurogenesis

Differences in chromatin state inherited from the parental gametes influence the regulation of maternal and paternal alleles in offspring. This phenomenon, known as genomic imprinting, results in genes preferentially transcribed from one parental allele. While local epigenetic factors such as DNA methylation are known to be important for the establishment of imprinted gene expression, less is known about the mechanisms by which differentially methylated regions (DMRs) lead to differences in allelic expression across broad stretches of chromatin. Allele-specific higher-order chromatin structure has been observed at multiple imprinted loci, consistent with the observation of allelic binding of the chromatin-organizing factor CTCF at multiple DMRs. However, whether allelic chromatin structure impacts allelic gene expression is not known for most imprinted loci. Here we characterize the mechanisms underlying brain-specific imprinted expression of the Peg13-Kcnk9 locus, an imprinted region associated with intellectual disability. We performed region capture Hi-C on mouse brain from reciprocal hybrid crosses and found imprinted higher-order chromatin structure caused by the allelic binding of CTCF to the Peg13 DMR. Using an in vitro neuron differentiation system, we show that on the maternal allele enhancer-promoter contacts formed early in development prime the brain-specific potassium leak channel Kcnk9 for maternal expression prior to neurogenesis. In contrast, these enhancer-promoter contacts are blocked by CTCF on the paternal allele, preventing paternal Kcnk9 activation. This work provides a high-resolution map of imprinted chromatin structure and demonstrates that chromatin state established in early development can promote imprinted expression upon differentiation.

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.

Hematopoietic/erythroid enhancers activate nearby target genes by extending histone H3K27ac and transcribing intergenic RNA

Enhancers activate gene transcription remotely, which requires tissue specific transcription factors binding to them. GATA1 and TAL1 are hematopoietic/erythroid-specific factors and often bind together to enhancers, activating target genes. Interestingly, we found that some hematopoietic/erythroid genes are transcribed in a GATA1-dependent but TAL1-independnet manner. They appear to have enhancers within a relatively short distance. In this study, we paired highly transcribed hematopoietic/erythroid genes with the nearest GATA1/TAL1-binding enhancers and analyzed these putative enhancer-gene pairs depending on distance between them. Enhancers located at various distances from genes in the pairs, which was not related to transcription level of the genes. However, genes with enhancers at short distances away tended to be transcriptionally unaffected by TAL1 depletion. Histone H3K27ac extended from the enhancers to target genes. The H3K27ac extension was maintained without TAL1, even though it disappeared owing to the loss of GATA1. Intergenic RNA was highly transcribed from the enhancers to nearby target genes, independent of TAL1. Taken together, TAL1-independent transcription of hematopoietic/erythroid genes appears to be promoted by enhancers present in a short distance. These enhancers are likely to activate nearby target genes by tracking the intervening regions.

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.

CTCF controls three-dimensional enhancer network underlying the inflammatory response of bone marrow-derived dendritic cells

Dendritic cells are antigen-presenting cells orchestrating innate and adaptive immunity. The crucial role of transcription factors and histone modifications in the transcriptional regulation of dendritic cells has been extensively studied. However, it is not been well understood whether and how three-dimensional chromatin folding controls gene expression in dendritic cells. Here we demonstrate that activation of bone marrow-derived dendritic cells induces extensive reprogramming of chromatin looping as well as enhancer activity, both of which are implicated in the dynamic changes in gene expression. Interestingly, depletion of CTCF attenuates GM-CSF-mediated JAK2/STAT5 signaling, resulting in defective NF-κB activation. Moreover, CTCF is necessary for establishing NF-κB-dependent chromatin interactions and maximal expression of pro-inflammatory cytokines, which prime Th1 and Th17 cell differentiation. Collectively, our study provides mechanistic insights into how three-dimensional enhancer networks control gene expression during bone marrow-derived dendritic cells activation, and offers an integrative view of the complex activities of CTCF in the inflammatory response of bone marrow-derived dendritic cells.

3D genomics and its applications in precision medicine

Three-dimensional (3D) genomics is an emerging discipline that studies the three-dimensional structure of chromatin and the three-dimensional and functions of genomes. It mainly studies the three-dimensional conformation and functional regulation of intranuclear genomes, such as DNA replication, DNA recombination, genome folding, gene expression regulation, transcription factor regulation mechanism, and the maintenance of three-dimensional conformation of genomes. Self-chromosomal conformation capture (3C) technology has been developed, and 3D genomics and related fields have developed rapidly. In addition, chromatin interaction analysis techniques developed by 3C technologies, such as paired-end tag sequencing (ChIA-PET) and whole-genome chromosome conformation capture (Hi-C), enable scientists to further study the relationship between chromatin conformation and gene regulation in different species. Thus, the spatial conformation of plant, animal, and microbial genomes, transcriptional regulation mechanisms, interaction patterns of chromosomes, and the formation mechanism of spatiotemporal specificity of genomes are revealed. With the help of new experimental technologies, the identification of key genes and signal pathways related to life activities and diseases is sustaining the rapid development of life science, agriculture, and medicine. In this paper, the concept and development of 3D genomics and its application in agricultural science, life science, and medicine are introduced, which provides a theoretical basis for the study of biological life processes.

Low-affinity CTCF binding drives transcriptional regulation whereas high-affinity binding encompasses architectural functions

CTCF is a DNA-binding protein which plays critical roles in chromatin structure organization and transcriptional regulation; however, little is known about the functional determinants of different CTCF-binding sites (CBS). Using a conditional mouse model, we have identified one set of CBSs that are lost upon CTCF depletion (lost CBSs) and another set that persists (retained CBSs). Retained CBSs are more similar to the consensus CTCF-binding sequence and usually span tandem CTCF peaks. Lost CBSs are enriched at enhancers and promoters and associate with active chromatin marks and higher transcriptional activity. In contrast, retained CBSs are enriched at TAD and loop boundaries. Integration of ChIP-seq and RNA-seq data has revealed that retained CBSs are located at the boundaries between distinct chromatin states, acting as chromatin barriers. Our results provide evidence that transient, lost CBSs are involved in transcriptional regulation, whereas retained CBSs are critical for establishing higher-order chromatin architecture.

Auxin-inducible degron 2 system deciphers functions of CTCF domains in transcriptional regulation

BACKGROUND

CTCF is a well-established chromatin architectural protein that also plays various roles in transcriptional regulation. While CTCF biology has been extensively studied, how the domains of CTCF function to regulate transcription remains unknown. Additionally, the original auxin-inducible degron 1 (AID1) system has limitations in investigating the function of CTCF.

RESULTS

We employ an improved auxin-inducible degron technology, AID2, to facilitate the study of acute depletion of CTCF while overcoming the limitations of the previous AID system. As previously observed through the AID1 system and steady-state RNA analysis, the new AID2 system combined with SLAM-seq confirms that CTCF depletion leads to modest nascent and steady-state transcript changes. A CTCF domain sgRNA library screening identifies the zinc finger (ZF) domain as the region within CTCF with the most functional relevance, including ZFs 1 and 10. Removal of ZFs 1 and 10 reveals genomic regions that independently require these ZFs for DNA binding and transcriptional regulation. Notably, loci regulated by either ZF1 or ZF10 exhibit unique CTCF binding motifs specific to each ZF.

CONCLUSIONS

By extensively comparing the AID1 and AID2 systems for CTCF degradation in SEM cells, we confirm that AID2 degradation is superior for achieving miniAID-tagged protein degradation without the limitations of the AID1 system. The model we create that combines AID2 depletion of CTCF with exogenous overexpression of CTCF mutants allows us to demonstrate how peripheral ZFs intricately orchestrate transcriptional regulation in a cellular context for the first time.

Isocitrate dehydrogenase 1 mutation drives leukemogenesis by PDGFRA activation due to insulator disruption in acute myeloid leukemia (AML)

Acute myeloid leukemia (AML) is characterized by complex molecular alterations and driver mutations. Elderly patients show increased frequencies of IDH mutations with high chemoresistance and relapse rates despite recent therapeutic advances. Besides being associated with global promoter hypermethylation, IDH1 mutation facilitated changes in 3D DNA-conformation by CTCF-anchor methylation and upregulated oncogene expression in glioma, correlating with poor prognosis. Here, we investigated the role of IDH1 p.R132H mutation in altering 3D DNA-architecture and subsequent oncogene activation in AML. Using public RNA-Seq data, we identified upregulation of tyrosine kinase PDGFRA in IDH1-mutant patients, correlating with poor prognosis. DNA methylation analysis identified CpG hypermethylation within a CTCF-anchor upstream of PDGFRA in IDH1-mutant patients. Increased PDGFRA expression, PDGFRA-CTCF methylation and decreased CTCF binding were confirmed in AML CRISPR cells with heterozygous IDH1 p.R132H mutation and upon exogenous 2-HG treatment. IDH1-mutant cells showed higher sensitivity to tyrosine kinase inhibitor dasatinib, which was supported by reduced blast count in a patient with refractory IDH1-mutant AML after dasatinib treatment. Our data illustrate that IDH1 p.R132H mutation leads to CTCF hypermethylation, disrupting DNA-looping and insulation of PDGFRA, resulting in PDGFRA upregulation in IDH1-mutant AML. Treatment with dasatinib may offer a novel treatment strategy for IDH1-mutant AML.

A global high-density chromatin interaction network reveals functional long-range and trans-chromosomal relationships

BACKGROUND

Chromatin contacts are essential for gene-expression regulation; however, obtaining a high-resolution genome-wide chromatin contact map is still prohibitively expensive owing to large genome sizes and the quadratic scale of pairwise data. Chromosome conformation capture (3C)-based methods such as Hi-C have been extensively used to obtain chromatin contacts. However, since the sparsity of these maps increases with an increase in genomic distance between contacts, long-range or trans-chromatin contacts are especially challenging to sample.

RESULTS

Here, we create a high-density reference genome-wide chromatin contact map using a meta-analytic approach. We integrate 3600 human, 6700 mouse, and 500 fly Hi-C experiments to create species-specific meta-Hi-C chromatin contact maps with 304 billion, 193 billion, and 19 billion contacts in respective species. We validate that meta-Hi-C contact maps are uniquely powered to capture functional chromatin contacts in both cis and trans. We find that while individual dataset Hi-C networks are largely unable to predict any long-range coexpression (median 0.54 AUC), meta-Hi-C networks perform comparably in both cis and trans (0.65 AUC vs 0.64 AUC). Similarly, for long-range expression quantitative trait loci (eQTL), meta-Hi-C contacts outperform all individual Hi-C experiments, providing an improvement over the conventionally used linear genomic distance-based association. Assessing between species, we find patterns of chromatin contact conservation in both cis and trans and strong associations with coexpression even in species for which Hi-C data is lacking.

CONCLUSIONS

We have generated an integrated chromatin interaction network which complements a large number of methodological and analytic approaches focused on improved specificity or interpretation. This high-depth "super-experiment" is surprisingly powerful in capturing long-range functional relationships of chromatin interactions, which are now able to predict coexpression, eQTLs, and cross-species relationships. The meta-Hi-C networks are available at https://labshare.cshl.edu/shares/gillislab/resource/HiC/ .

CTCF blocks antisense transcription initiation at divergent promoters

Transcription at most promoters is divergent, initiating at closely spaced oppositely oriented core promoters to produce sense transcripts along with often unstable upstream antisense transcripts (uasTrx). How antisense transcription is regulated and to what extent it is coordinated with sense transcription is not well understood. Here, by combining acute degradation of the multi-functional transcription factor CTCF and nascent transcription measurements, we find that CTCF specifically suppresses antisense but not sense transcription at hundreds of divergent promoters. Primary transcript RNA-FISH shows that CTCF lowers burst fraction but not burst intensity of uasTrx and that co-bursting of sense and antisense transcripts is disfavored. Genome editing, chromatin conformation studies and high-resolution transcript mapping revealed that precisely positioned CTCF directly suppresses the initiation of uasTrx, in a manner independent of its architectural function. In sum, CTCF shapes the transcriptional landscape in part by suppressing upstream antisense transcription.

CTCF blocks antisense transcription initiation at divergent promoters

Transcription at most promoters is divergent, initiating at closely spaced oppositely oriented core promoters to produce sense transcripts along with often unstable upstream antisense transcripts (uasTrx). How antisense transcription is regulated and to what extent it is coordinated with sense transcription is not well understood. Here, by combining acute degradation of the multi-functional transcription factor CTCF and nascent transcription measurements, we find that CTCF specifically suppresses antisense but not sense transcription at hundreds of divergent promoters. Primary transcript RNA-FISH shows that CTCF lowers burst fraction but not burst intensity of uasTrx and that co-bursting of sense and antisense transcripts is disfavored. Genome editing, chromatin conformation studies and high-resolution transcript mapping revealed that precisely positioned CTCF directly suppresses the initiation of uasTrx, in a manner independent of its architectural function. In sum, CTCF shapes the transcriptional landscape in part by suppressing upstream antisense transcription.

Induced hepatic stem cells maintain self-renewal through the high expression of Myc coregulated by TET1 and CTCF

Background

Induced hepatic stem cells (iHepSCs) with the capacities of self-renewal and bidifferentiation into hepatocytes and cholangiocytes were generated from mouse embryonic fibroblasts (MEFs) by lineage reprogramming in our previous research. However, the mechanism of iHepSC self-renewal has not been elucidated. Active demethylation regulated by Tet1 plays an important role in the self-renewal of stem cells, including pluripotent stem cells and adult stem cells. Here, we investigated the role and mechanism of Tet1-regulated demethylation in the self-renewal of iHepSCs.

Methods

The methylation levels and the expression of Tet1 in iHepSCs and MEFs were analyzed by immunofluorescent staining, quantitative reverse transcription PCR and western blotting. Then, the effects of Tet1 knockdown on the proliferation and self-renewal of iHepSCs were analyzed by CCK8, colony formation, and sphere formation assays. The mechanism by which Tet1 regulates the self-renewal of iHepSCs was investigated by chromatin immunoprecipitation, bisulfite sequence PCR, and methylation-sensitive restriction endonuclease-PCR.

Results

The high level of 5hmC and the low level of 5mC in iHepSCs were accompanied by high expression of Tet1. After Tet1 expression was knocked down by shRNA in iHepSCs, the proliferation and self-renewal capacities were inhibited, and the expression of Myc was also decreased. The higher expression level of Myc in iHepSCs maintained its self-renewal and was regulated by Tet1, which directly binds to CBS-1 and site A regions of the Myc promoter and demethylates the CpG cytosine. In addition, CTCF also binds to the CBS-1 and site A regions of the Myc promoter and regulates Myc expression along with TET1.

Conclusion

The self-renewal of iHepSCs was maintained by the higher expression of Myc, which was coregulated by TET1 and CTCF. This study may provide new insights into the self-renewal of stem cells, which can promote the research and application of ‘reprogrammed’ stem cells.

Structural variations in cancer and the 3D genome

Structural variations (SVs) affect more of the cancer genome than any other type of somatic genetic alteration but difficulties in detecting and interpreting them have limited our understanding. Clinical cancer sequencing also increasingly aims to detect SVs, leading to a widespread necessity to interpret their biological and clinical relevance. Recently, analyses of large whole-genome sequencing data sets revealed features that impact rates of SVs across the genome in different cancers. A striking feature has been the extent to which, in both their generation and their influence on the selective fitness of cancer cells, SVs are more specific to individual cancer types than other genetic alterations such as single-nucleotide variants. This Perspective discusses how the folding of the 3D genome, and differences in its folding across cell types, affect observed SV rates in different cancer types as well as how SVs can impact cancer cell fitness.

RNA helicase-dependent gene looping impacts messenger RNA processing

DDX5 and DDX17 are DEAD-box RNA helicase paralogs which regulate several aspects of gene expression, especially transcription and splicing, through incompletely understood mechanisms. A transcriptome analysis of DDX5/DDX17-depleted human cells confirmed the large impact of these RNA helicases on splicing and revealed a widespread deregulation of 3' end processing. In silico analyses and experiments in cultured cells showed the binding and functional contribution of the genome organizing factor CTCF to chromatin sites at or near a subset of DDX5/DDX17-dependent exons that are characterized by a high GC content and a high density of RNA Polymerase II. We propose the existence of an RNA helicase-dependent relationship between CTCF and the dynamics of transcription across DNA and/or RNA structured regions, that contributes to the processing of internal and terminal exons. Moreover, local DDX5/DDX17-dependent chromatin loops spatially connect RNA helicase-regulated exons with their cognate promoter, and we provide the first direct evidence that de novo gene looping modifies alternative splicing and polyadenylation. Overall our findings uncover the impact of DDX5/DDX17-dependent chromatin folding on pre-messenger RNA processing.

Super-Enhancers, Phase-Separated Condensates, and 3D Genome Organization in Cancer

3D chromatin organization plays an important role in transcription regulation and gene expression. The 3D genome is highly maintained by several architectural proteins, such as CTCF, Yin Yang 1, and cohesin complex. This structural organization brings regulatory DNA elements in close proximity to their target promoters. In this review, we discuss the 3D chromatin organization of super-enhancers and their relationship to phase-separated condensates. Super-enhancers are large clusters of DNA elements. They can physically contact with their target promoters by chromatin looping during transcription. Multiple transcription factors can bind to enhancer and promoter sequences and recruit a complex array of transcriptional co-activators and RNA polymerase II to effect transcriptional activation. Phase-separated condensates of transcription factors and transcriptional co-activators have been implicated in assembling the transcription machinery at particular enhancers. Cancer cells can hijack super-enhancers to drive oncogenic transcription to promote cell survival and proliferation. These dysregulated transcriptional programs can cause cancer cells to become highly dependent on transcriptional regulators, such as Mediator and BRD4. Moreover, the expression of oncogenes that are driven by super-enhancers is sensitive to transcriptional perturbation and often occurs in phase-separated condensates, supporting therapeutic rationales of targeting SE components, 3D genome organization, or dysregulated condensates in cancer.

Coming full circle: on the origin and evolution of the looping model for enhancer-promoter communication

In mammalian organisms, enhancers can regulate transcription from great genomic distances. How enhancers affect distal gene expression has been a major question in the field of gene regulation. One model to explain how enhancers communicate with their target promoters, the chromatin looping model, posits that enhancers and promoters come in close spatial proximity to mediate communication. Chromatin looping has been broadly accepted as a means for enhancer–promoter communication, driven by accumulating in vitro and in vivo evidence. The genome is now known to be folded into a complex 3D arrangement, created and maintained in part by the interplay of the Cohesin complex and the DNA-binding protein CTCF. In the last few years, however, doubt over the relationship between looping and transcriptional activation has emerged, driven by studies finding that only a modest number of genes are perturbed with acute degradation of looping machinery components. In parallel, newer models describing distal enhancer action have also come to prominence. In this article, we explore the emergence and development of the looping model as a means for enhancer–promoter communication and review the contrasting evidence between historical gene-specific and current global data for the role of chromatin looping in transcriptional regulation. We also discuss evidence for alternative models to chromatin looping and their support in the literature. We suggest that, while there is abundant evidence for chromatin looping as a major mechanism for enhancer function, enhancer–promoter communication is likely mediated by more than one mechanism in an enhancer- and context-dependent manner.

Implications of Dosage Deficiencies in CTCF and Cohesin on Genome Organization, Gene Expression, and Human Neurodevelopment

Properly organizing DNA within the nucleus is critical to ensure normal downstream nuclear functions. CTCF and cohesin act as major architectural proteins, working in concert to generate thousands of high-intensity chromatin loops. Due to their central role in loop formation, a massive research effort has been dedicated to investigating the mechanism by which CTCF and cohesin create these loops. Recent results lead to questioning the direct impact of CTCF loops on gene expression. Additionally, results of controlled depletion experiments in cell lines has indicated that genome architecture may be somewhat resistant to incomplete deficiencies in CTCF or cohesin. However, heterozygous human genetic deficiencies in CTCF and cohesin have illustrated the importance of their dosage in genome architecture, cellular processes, animal behavior, and disease phenotypes. Thus, the importance of considering CTCF or cohesin levels is especially made clear by these heterozygous germline variants that characterize genetic syndromes, which are increasingly recognized in clinical practice. Defined primarily by developmental delay and intellectual disability, the phenotypes of CTCF and cohesin deficiency illustrate the importance of architectural proteins particularly in neurodevelopment. We discuss the distinct roles of CTCF and cohesin in forming chromatin loops, highlight the major role that dosage of each protein plays in the amplitude of observed effects on gene expression, and contrast these results to heterozygous mutation phenotypes in murine models and clinical patients. Insights highlighted by this comparison have implications for future research into these newly emerging genetic syndromes.

MYC overexpression leads to increased chromatin interactions at superenhancers and MYC binding sites

The MYC oncogene encodes for the MYC protein and is frequently dysregulated across multiple cancer cell types, making it an attractive target for cancer therapy. MYC overexpression leads to MYC binding at active enhancers, resulting in a global transcriptional amplification of active genes. Because super-enhancers are frequently dysregulated in cancer, we hypothesized that MYC preferentially invades into super-enhancers and alters the cancer genome organization. To that end, we performed ChIP-seq, RNA-seq, circular chromosome conformation capture (4C-seq), and Spike-in Quantitative Hi-C (SIQHiC) on the U2OS osteosarcoma cell line with tetracycline-inducible MYC. MYC overexpression in U2OS cells modulated histone acetylation and increased MYC binding at super-enhancers. SIQHiC analysis revealed increased global chromatin contact frequency, particularly at chromatin interactions connecting MYC binding sites at promoters and enhancers. Immunofluorescence staining showed that MYC molecules formed punctate foci at these transcriptionally active domains after MYC overexpression. These results demonstrate the accumulation of overexpressed MYC at promoter–enhancer hubs and suggest that MYC invades into enhancers through spatial proximity. At the same time, the increased protein–protein interactions may strengthen these chromatin interactions to increase chromatin contact frequency. CTCF siRNA knockdown in MYC-overexpressed U2OS cells demonstrated that removal of architectural proteins can disperse MYC and abrogate the increase in chromatin contacts. By elucidating the chromatin landscape of MYC-driven cancers, we can potentially target MYC-associated chromatin interactions for cancer therapy.

Capture-C: a modular and flexible approach for high-resolution chromosome conformation capture

Chromosome conformation capture (3C) methods measure the spatial proximity between DNA elements in the cell nucleus. Many methods have been developed to sample 3C material, including the Capture-C family of protocols. Capture-C methods use oligonucleotides to enrich for interactions of interest from sequencing-ready 3C libraries. This approach is modular and has been adapted and optimized to work for sampling of disperse DNA elements (NuTi Capture-C), including from low cell inputs (LI Capture-C), as well as to generate Hi-C like maps for specific regions of interest (Tiled-C) and to interrogate multiway interactions (Tri-C). We present the design, experimental protocol and analysis pipeline for NuTi Capture-C in addition to the variations for generation of LI Capture-C, Tiled-C and Tri-C data. The entire procedure can be performed in 3 weeks and requires standard molecular biology skills and equipment, access to a next-generation sequencing platform, and basic bioinformatic skills. Implemented with other sequencing technologies, these methods can be used to identify regulatory interactions and to compare the structural organization of the genome in different cell types and genetic models.

A CTCF-Binding Element and Histone Deacetylation Cooperatively Maintain Chromatin Loops, Linking to Long-Range Gene Regulation in Cancer Genomes

Background

Genes spanning long chromosomal domains are coordinately regulated in human genome, which contribute to global gene dysregulation and carcinogenesis in cancer. It has been noticed that epigenetic modification and chromatin architecture may participate in the regulation process. However, the regulation patterns and functional elements of long-range gene regulation are unclear.

Methods

Based on the clinical transcriptome data from different tumor sets, a novel expressional correlation analysis pipeline was performed to classify the co-regulated regions and subsets of intercorrelated regions. The GLAM2 program was used to predict conserved DNA elements that enriched in regions. Two conserved elements were selected to delete in Ishikawa and HeLa cells by CRISPR-Cas9. SAHA treatment and HDAC knockdown were used to change the histone acetylation status. Using qPCR, MTT, and scratch healing assay, we evaluate the effect on gene expression and cancer cell phenotype. By DNA pull-down and ChIP, the element-binding proteins were testified. 3C and 3D-FISH were performed to depict the alteration in chromatin architecture.

Results

In multiple cancer genomes, we classified subsets of coordinately regulated regions (sub-CRRs) that possibly shared the same regulatory mechanisms and exhibited similar expression patterns. A new conserved DNA element (CRE30) was enriched in sub-CRRs and associated with cancer patient survival. CRE30 could restrict gene regulation in sub-CRRs and affect cancer cell phenotypes. DNA pull-down showed that multiple proteins including CTCF were recruited on the CRE30 locus, and ChIP assay confirmed the CTCF-binding signals. Subsequent results uncovered that as an essential element, CRE30 maintained chromatin loops and mediated a compact chromatin architecture. Moreover, we found that blocking global histone deacetylation induced chromatin loop disruption and CTCF dropping in the region containing CRE30, linked to promoted gene regulation. Additionally, similar effects were observed with CRE30 deletion in another locus of chromosome 8.

Conclusions

Our research clarified a new functional element that recruits CTCF and collaborates with histone deacetylation to maintain high-order chromatin organizations, linking to long-range gene regulation in cancer genomes. The findings highlight a close relationship among conserved DNA element, epigenetic modification, and chromatin architecture in long-range gene regulation process.

Canonical WNT signaling-dependent gating of MYC requires a noncanonical CTCF function at a distal binding site

Abnormal WNT signaling increases MYC expression in colon cancer cells in part via oncogenic super-enhancer-(OSE)-mediated gating of the active MYC to the nuclear pore in a poorly understood process. We show here that the principal tenet of the WNT-regulated MYC gating, facilitating nuclear export of the MYC mRNA, is regulated by a CTCF binding site (CTCFBS) within the OSE to confer growth advantage in HCT-116 cells. To achieve this, the CTCFBS directs the WNT-dependent trafficking of the OSE to the nuclear pore from intra-nucleoplasmic positions in a stepwise manner. Once the OSE reaches a peripheral position, which is triggered by a CTCFBS-mediated CCAT1 eRNA activation, its final stretch (≤0.7 μm) to the nuclear pore requires the recruitment of AHCTF1, a key nucleoporin, to the CTCFBS. Thus, a WNT/ß-catenin-AHCTF1-CTCF-eRNA circuit enables the OSE to promote pathological cell growth by coordinating the trafficking of the active MYC gene within the 3D nuclear architecture.

Gene-gating of a MYC oncogenic super-enhancer (OSE) increases its expression in colon cancer cells in a poorly understood process. Here the authors show that MYC gating requires a CTCF binding site (CTCFBS) within the OSE that directs the stepwise trafficking of the OSE to the nuclear pore to facilitate increased nuclear export of MYC mRNA, which results in a growth advantage.

OUP accepted manuscript

Spatial genome organization is tightly controlled by several regulatory mechanisms and is essential for gene expression control. Nuclear receptors are ligand-activated transcription factors that modulate physiological and pathophysiological processes and are primary pharmacological targets. DNA binding of the important loop-forming insulator protein CCCTC-binding factor (CTCF) was modulated by 1α,25-dihydroxyvitamin D₃ (1,25(OH)₂D₃). We performed CTCF HiChIP assays to produce the first genome-wide dataset of CTCF long-range interactions in 1,25(OH)₂D₃-treated cells, and to determine whether dynamic changes of spatial chromatin interactions are essential for fine-tuning of nuclear receptor signaling. We detected changes in 3D chromatin organization upon vitamin D receptor (VDR) activation at 3.1% of all observed CTCF interactions. VDR binding was enriched at both differential loop anchors and within differential loops. Differential loops were observed in several putative functional roles including TAD border formation, promoter-enhancer looping, and establishment of VDR-responsive insulated neighborhoods. Vitamin D target genes were enriched in differential loops and at their anchors. Secondary vitamin D effects related to dynamic chromatin domain changes were linked to location of downstream transcription factors in differential loops. CRISPR interference and loop anchor deletion experiments confirmed the functional relevance of nuclear receptor ligand-induced adjustments of the chromatin 3D structure for gene expression regulation.

CTCF-mediated chromatin looping provides a topological framework for the formation of phase-separated transcriptional condensates

CTCF is crucial to the organization of mammalian genomes into loop structures. According to recent studies, the transcription apparatus is compartmentalized and concentrated at super-enhancers to form phase-separated condensates and drive the expression of cell-identity genes. However, it remains unclear whether and how transcriptional condensates are coupled to higher-order chromatin organization. Here, we show that CTCF is essential for RNA polymerase II (Pol II)-mediated chromatin interactions, which occur as hyperconnected spatial clusters at super-enhancers. We also demonstrate that CTCF clustering, unlike Pol II clustering, is independent of liquid-liquid phase-separation and resistant to perturbation of transcription. Interestingly, clusters of Pol II, BRD4, and MED1 were found to dissolve upon CTCF depletion, but were reinstated upon restoration of CTCF, suggesting a potent instructive function for CTCF in the formation of transcriptional condensates. Overall, we provide evidence suggesting that CTCF-mediated chromatin looping acts as an architectural prerequisite for the assembly of phase-separated transcriptional condensates.

Identification of chromatin states during zebrafish gastrulation using CUT&RUN and CUT&Tag

BACKGROUND

Cell fate decisions are governed by interactions between sequence-specific transcription factors and a dynamic chromatin landscape. Zebrafish offer a powerful system for probing the mechanisms that drive these cell fate choices, especially in the context of early embryogenesis. However, technical challenges associated with conventional methods for chromatin profiling have slowed progress toward understanding the exact relationships between chromatin changes, transcription factor binding, and cellular differentiation during zebrafish embryogenesis.

RESULTS

To overcome these challenges, we adapted the chromatin profiling methods Cleavage Under Targets and Release Using Nuclease (CUT&RUN) and CUT&Tag for use in zebrafish and applied these methods to generate high-resolution enrichment maps for H3K4me3, H3K27me3, H3K9me3, RNA polymerase II, and the histone variant H2A.Z using tissue isolated from whole, mid-gastrula stage embryos. Using this data, we identify a subset of genes that may be bivalently regulated during both zebrafish and mouse gastrulation, provide evidence for an evolving H2A.Z landscape during embryo development, and demonstrate the effectiveness of CUT&RUN for detecting H3K9me3 enrichment at repetitive sequences.

CONCLUSIONS

Our results demonstrate the power of combining CUT&RUN and CUT&Tag methods with the strengths of the zebrafish system to define emerging chromatin landscapes in the context of vertebrate embryogenesis.

CTCF and transcription influence chromatin structure re-configuration after mitosis

During mitosis, transcription is globally attenuated and chromatin architecture is dramatically reconfigured. We exploited the M- to G1-phase progression to interrogate the contributions of the architectural factor CTCF and the process of transcription to genome re-sculpting in newborn nuclei. Depletion of CTCF during the M- to G1-phase transition alters short-range compartmentalization after mitosis. Chromatin domain boundary re-formation is impaired upon CTCF loss, but a subset of boundaries, characterized by transitions in chromatin states, is established normally. Without CTCF, structural loops fail to form, leading to illegitimate contacts between cis-regulatory elements (CREs). Transient CRE contacts that are normally resolved after telophase persist deeply into G1-phase in CTCF-depleted cells. CTCF loss-associated gains in transcription are often linked to increased, normally illegitimate enhancer-promoter contacts. In contrast, at genes whose expression declines upon CTCF loss, CTCF seems to function as a conventional transcription activator, independent of its architectural role. CTCF-anchored structural loops facilitate formation of CRE loops nested within them, especially those involving weak CREs. Transcription inhibition does not significantly affect global architecture or transcription start site-associated boundaries. However, ongoing transcription contributes considerably to the formation of gene domains, regions of enriched contacts along gene bodies. Notably, gene domains emerge in ana/telophase prior to completion of the first round of transcription, suggesting that epigenetic features in gene bodies contribute to genome reconfiguration prior to transcription. The focus on the de novo formation of nuclear architecture during G1 entry yields insights into the contributions of CTCF and transcription to chromatin architecture dynamics during the mitosis to G1-phase progression.

Cohesin mutations in myeloid malignancies

Cohesin is a multisubunit protein complex that forms a ring-like structure around DNA. It is essential for sister chromatid cohesion, chromatin organization, transcriptional regulation, and DNA damage repair and plays a major role in dynamically shaping the genome architecture and maintaining DNA integrity. The core complex subunits STAG2, RAD21, SMC1, and SMC3, as well as its modulators PDS5A/B, WAPL, and NIPBL, have been found to be recurrently mutated in hematologic and solid malignancies. These mutations are found across the full spectrum of myeloid neoplasia, including pediatric Down syndrome-associated acute megakaryoblastic leukemia, myelodysplastic syndromes, chronic myelomonocytic leukemia, and de novo and secondary acute myeloid leukemias. The mechanisms by which cohesin mutations act as drivers of clonal expansion and disease progression are still poorly understood. Recent studies have described the impact of cohesin alterations on self-renewal and differentiation of hematopoietic stem and progenitor cells, which are associated with changes in chromatin and epigenetic state directing lineage commitment, as well as genomic integrity. Herein, we review the role of the cohesin complex in healthy and malignant hematopoiesis. We discuss clinical implications of cohesin mutations in myeloid malignancies and discuss opportunities for therapeutic targeting.

Acute depletion of CTCF rewires genome-wide chromatin accessibility

BACKGROUND

The transcription factor CTCF appears indispensable in defining topologically associated domain boundaries and maintaining chromatin loop structures within these domains, supported by numerous functional studies. However, acute depletion of CTCF globally reduces chromatin interactions but does not significantly alter transcription.

RESULTS

Here, we systematically integrate multi-omics data including ATAC-seq, RNA-seq, WGBS, Hi-C, Cut&Run, and CRISPR-Cas9 survival dropout screens, and time-solved deep proteomic and phosphoproteomic analyses in cells carrying auxin-induced degron at endogenous CTCF locus. Acute CTCF protein degradation markedly rewires genome-wide chromatin accessibility. Increased accessible chromatin regions are frequently located adjacent to CTCF-binding sites at promoter regions and insulator sites associated with enhanced transcription of nearby genes. In addition, we use CTCF-associated multi-omics data to establish a combinatorial data analysis pipeline to discover CTCF co-regulatory partners. We successfully identify 40 candidates, including multiple established partners. Interestingly, many CTCF co-regulators that have alterations of their respective downstream gene expression do not show changes of their own expression levels across the multi-omics measurements upon acute CTCF loss, highlighting the strength of our system to discover hidden co-regulatory partners associated with CTCF-mediated transcription.

CONCLUSIONS

This study highlights that CTCF loss rewires genome-wide chromatin accessibility, which plays a critical role in transcriptional regulation.

CTCF as a regulator of alternative splicing: new tricks for an old player

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

Chromatin Alterations in Neurological Disorders and Strategies of (Epi)Genome Rescue

Neurological disorders (NDs) comprise a heterogeneous group of conditions that affect the function of the nervous system. Often incurable, NDs have profound and detrimental consequences on the affected individuals' lives. NDs have complex etiologies but commonly feature altered gene expression and dysfunctions of the essential chromatin-modifying factors. Hence, compounds that target DNA and histone modification pathways, the so-called epidrugs, constitute promising tools to treat NDs. Yet, targeting the entire epigenome might reveal insufficient to modify a chosen gene expression or even unnecessary and detrimental to the patients' health. New technologies hold a promise to expand the clinical toolkit in the fight against NDs. (Epi)genome engineering using designer nucleases, including CRISPR-Cas9 and TALENs, can potentially help restore the correct gene expression patterns by targeting a defined gene or pathway, both genetically and epigenetically, with minimal off-target activity. Here, we review the implication of epigenetic machinery in NDs. We outline syndromes caused by mutations in chromatin-modifying enzymes and discuss the functional consequences of mutations in regulatory DNA in NDs. We review the approaches that allow modifying the (epi)genome, including tools based on TALENs and CRISPR-Cas9 technologies, and we highlight how these new strategies could potentially change clinical practices in the treatment of NDs.

Multiple CTCF sites cooperate with each other to maintain a TAD for enhancer-promoter interaction in the β-globin locus

Insulators are cis-regulatory elements that block enhancer activity and prevent heterochromatin spreading. The binding of CCCTC-binding factor (CTCF) protein is essential for insulators to play the roles in a chromatin context. The β-globin locus, consisting of multiple genes and enhancers, is flanked by two insulators 3'HS1 and HS5. However, it has been reported that the absence of these insulators did not affect the β-globin transcription. To explain the unexpected finding, we have deleted a CTCF motif at 3'HS1 or HS5 in the human β-globin locus and analyzed chromatin interactions around the locus. It was found that a topologically associating domain (TAD) containing the β-globin locus is maintained by neighboring CTCF sites in the CTCF motif-deleted loci. The additional deletions of neighboring CTCF motifs disrupted the β-globin TAD, resulting in decrease of the β-globin transcription. Chromatin interactions of the β-globin enhancers with gene promoter were weakened in the multiple CTCF motifs-deleted loci, even though the enhancers have still active chromatin features such as histone H3K27ac and histone H3 depletion. Genome-wide analysis using public CTCF ChIA-PET and ChIP-seq data showed that chromatin domains possessing multiple CTCF binding sites tend to contain super-enhancers like the β-globin enhancers. Taken together, our results show that multiple CTCF sites surrounding the β-globin locus cooperate with each other to maintain a TAD. The β-globin TAD appears to provide a compact spatial environment that enables enhancers to interact with promoter.

Germline RUNX1 variation and predisposition to childhood acute lymphoblastic leukemia

Genetic alterations in the RUNX1 gene are associated with benign and malignant blood disorders, particularly of megakaryocyte and myeloid lineages. The role of RUNX1 in acute lymphoblastic leukemia (ALL) is less clear, particularly how germline genetic variation influences the predisposition to this type of leukemia. Sequencing 4,836 children with B-ALL and 1,354 cases of T-ALL, we identified 31 and 18 germline RUNX1 variants, respectively. RUNX1 variants in B-ALL consistently showed minimal damaging effects. By contrast, 6 T-ALL-related variants result in drastic loss of RUNX1 activity as a transcription activator in vitro. Ectopic expression of dominant-negative RUNX1 variants in human CD34+ cells repressed differentiation into erythroid, megakaryocytes, and T cells, while promoting myeloid cell development. Chromatin immunoprecipitation sequencing of T-ALL models showed distinctive patterns of RUNX1 binding by variant proteins. Further whole genome sequencing identified JAK3 mutation as the most frequent somatic genomic abnormality in T-ALL with germline RUNX1 variants. Co-introduction of RUNX1 variant and JAK3 mutation in hematopoietic stem and progenitor cells in mice gave rise to T-ALL with early T-cell precursor phenotype. Taken together, these results indicated that RUNX1 is an important predisposition gene for T-ALL and pointed to novel biology of RUNX1-mediated leukemogenesis in the lymphoid lineages.