The BET Protein BRD2 Cooperates with CTCF to Enforce Transcriptional and Architectural Boundaries

Boundary sequences flanking the mouse tyrosinase locus ensure faithful pattern of gene expression

Control of gene expression is dictated by cell-type specific regulatory sequences that physically organize the structure of chromatin, including promoters, enhancers and insulators. While promoters and enhancers convey cell-type specific activating signals, insulators prevent the cross-talk of regulatory elements within adjacent loci and safeguard the specificity of action of promoters and enhancers towards their targets in a tissue specific manner. Using the mouse tyrosinase (Tyr) locus as an experimental model, a gene whose mutations are associated with albinism, we described the chromatin structure in cells at two distinct transcriptional states. Guided by chromatin structure, through the use of Chromosome Conformation Capture (3C), we identified sequences at the 5' and 3' boundaries of this mammalian gene that function as enhancers and insulators. By CRISPR/Cas9-mediated chromosomal deletion, we dissected the functions of these two regulatory elements in vivo in the mouse, at the endogenous chromosomal context, and proved their mechanistic role as genomic insulators, shielding the Tyr locus from the expression patterns of adjacent genes.

Erythropoietin Regulates Transcription and YY1 Dynamics in a Pre-established Chromatin Architecture

The three-dimensional architecture of the genome plays an essential role in establishing and maintaining cell identity. However, the magnitude and temporal kinetics of changes in chromatin structure that arise during cell differentiation remain poorly understood. Here, we leverage a murine model of erythropoiesis to study the relationship between chromatin conformation, the epigenome, and transcription in erythroid cells. We discover that acute transcriptional responses induced by erythropoietin (EPO), the hormone necessary for erythroid differentiation, occur within an invariant chromatin topology. Within this pre-established landscape, Yin Yang 1 (YY1) occupancy dynamically redistributes to sites in proximity of EPO-regulated genes. Using HiChIP, we identify chromatin contacts mediated by H3K27ac and YY1 that are enriched for enhancer-promoter interactions of EPO-responsive genes. Taken together, these data are consistent with an emerging model that rapid, signal-dependent transcription occurs in the context of a pre-established chromatin architecture.

•

EPO induces rapid RNA Pol II response at a key subset of genes

•

YY1 is redistributed in the genome following 1 h EPO stimulation

•

CTCF and YY1 bind different locations pre and post 1 h EPO stimulation

•

E-P loops mediated by H3K27ac are largely invariant in response to EPO

GATA-1-dependent histone H3K27 acetylation mediates erythroid cell-specific chromatin interaction between CTCF sites

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

Inhibitors of bromodomain and extra-terminal proteins for treating multiple human diseases

Clinical development of bromodomain and extra‐terminal (BET) protein inhibitors differs from the traditional course of drug development. These drugs are simultaneously being evaluated for treating a wide spectrum of human diseases due to their novel mechanism of action. BET proteins are epigenetic “readers,” which play a primary role in transcription. Here, we briefly describe the BET family of proteins, of which BRD4 has been studied most extensively. We discuss BRD4 activity at latent enhancers as an example of BET protein function. We examine BRD4 redistribution and enhancer reprogramming in embryonic development, cancer, cardiovascular, autoimmune, and metabolic diseases, presenting hallmark studies that highlight BET proteins as attractive targets for therapeutic intervention. We review the currently available approaches to targeting BET proteins, methods of selectively targeting individual bromodomains, and review studies that compare the effects of selective BET inhibition to those of pan‐BET inhibition. Lastly, we examine the current clinical landscape of BET inhibitor development.

Domain-selective targeting of BET proteins in cancer and immunological diseases

Cancer and inflammation are strongly interconnected processes. Chronic inflammatory pathologies can be at the heart of tumor development; similarly, tumor-elicited inflammation is a consequence of many cancers. The mechanistic interdependence between cancer and inflammatory pathologies points toward common protein effectors which represent potential shared targets for pharmacological intervention. Epigenetic mechanisms often drive resistance to cancer therapy and immunomodulatory strategies. The bromodomain and extraterminal domain (BET) proteins are epigenetic adapters which play a major role in controlling cell proliferation and the production of inflammatory mediators. A plethora of small molecules aimed at inhibiting BET protein function to treat cancer and inflammatory diseases have populated academic and industry efforts in the last 10 years. In this review, we will discuss recent pharmacological approaches aimed at targeting a single or a subset of the eight bromodomains within the BET family which have the potential to tease apart clinical efficacy and safety signals of BET inhibitors.

CGGBP1-regulated cytosine methylation at CTCF-binding motifs resists stochasticity

Background

The human CGGBP1 binds to GC-rich regions and interspersed repeats, maintains homeostasis of stochastic cytosine methylation and determines DNA-binding of CTCF. Interdependence between regulation of cytosine methylation and CTCF occupancy by CGGBP1 remains unknown.

Results

By analyzing methylated DNA-sequencing data obtained from CGGBP1-depleted cells, we report that some transcription factor-binding sites, including CTCF, resist stochastic changes in cytosine methylation. By analysing CTCF-binding sites we show that cytosine methylation changes at CTCF motifs caused by CGGBP1 depletion resist stochastic changes. These CTCF-binding sites are positioned at locations where the spread of cytosine methylation in cis depends on the levels of CGGBP1.

Conclusion

Our findings suggest that CTCF occupancy and functions are determined by CGGBP1-regulated cytosine methylation patterns.

Cohesin-Dependent and -Independent Mechanisms Mediate Chromosomal Contacts between Promoters and Enhancers

It is currently assumed that 3D chromosomal organization plays a central role in transcriptional control. However, depletion of cohesin and CTCF affects the steady-state levels of only a minority of transcripts. Here, we use high-resolution Capture Hi-C to interrogate the dynamics of chromosomal contacts of all annotated human gene promoters upon degradation of cohesin and CTCF. We show that a majority of promoter-anchored contacts are lost in these conditions, but many contacts with distinct properties are maintained, and some new ones are gained. The rewiring of contacts between promoters and active enhancers upon cohesin degradation associates with rapid changes in target gene transcription as detected by SLAM sequencing (SLAM-seq). These results provide a mechanistic explanation for the limited, but consistent, effects of cohesin and CTCF depletion on steady-state transcription and suggest the existence of both cohesin-dependent and -independent mechanisms of enhancer-promoter pairing.

Synthetic Lethal and Resistance Interactions with BET Bromodomain Inhibitors in Triple-Negative Breast Cancer

BET bromodomain inhibitors (BBDIs) are candidate therapeutic agents for triple-negative breast cancer (TNBC) and other cancer types, but inherent and acquired resistance to BBDIs limits their potential clinical use. Using CRISPR and small-molecule inhibitor screens combined with comprehensive molecular profiling of BBDI response and resistance, we identified synthetic lethal interactions with BBDIs and genes that, when deleted, confer resistance. We observed synergy with regulators of cell cycle progression, YAP, AXL, and SRC signaling, and chemotherapeutic agents. We also uncovered functional similarities and differences among BRD2, BRD4, and BRD7. Although deletion of BRD2 enhances sensitivity to BBDIs, BRD7 loss leads to gain of TEAD-YAP chromatin binding and luminal features associated with BBDI resistance. Single-cell RNA-seq, ATAC-seq, and cellular barcoding analysis of BBDI responses in sensitive and resistant cell lines highlight significant heterogeneity among samples and demonstrate that BBDI resistance can be pre-existing or acquired.

A WIZ/Cohesin/CTCF Complex Anchors DNA Loops to Define Gene Expression and Cell Identity

Chromosome structure is a key regulator of gene expression. CTCF and cohesin play critical roles in structuring chromosomes by mediating physical interactions between distant genomic sites. The resulting DNA loops often contain genes and their cis-regulatory elements. Despite the importance of DNA loops in maintaining proper transcriptional regulation and cell identity, there is limited understanding of the molecular mechanisms that regulate their dynamics and function. We report a previously unrecognized role for WIZ (widely interspaced zinc finger-containing protein) in DNA loop architecture and regulation of gene expression. WIZ forms a complex with cohesin and CTCF that occupies enhancers, promoters, insulators, and anchors of DNA loops. Aberrant WIZ function alters cohesin occupancy and increases the number of DNA loop structures in the genome. WIZ is required for proper gene expression and transcriptional insulation. Our results uncover an unexpected role for WIZ in DNA loop architecture, transcriptional control, and maintenance of cell identity.

DeepMILO: a deep learning approach to predict the impact of non-coding sequence variants on 3D chromatin structure

Non-coding variants have been shown to be related to disease by alteration of 3D genome structures. We propose a deep learning method, DeepMILO, to predict the effects of variants on CTCF/cohesin-mediated insulator loops. Application of DeepMILO on variants from whole-genome sequences of 1834 patients of twelve cancer types revealed 672 insulator loops disrupted in at least 10% of patients. Our results show mutations at loop anchors are associated with upregulation of the cancer driver genes BCL2 and MYC in malignant lymphoma thus pointing to a possible new mechanism for their dysregulation via alteration of insulator loops.

Bromodomain and extra-terminal domain (BET) proteins regulate melanocyte differentiation

BACKGROUND

Pharmacologic inhibition of bromodomain and extra-terminal (BET) proteins is currently being explored as a new therapeutic approach in cancer. Some studies have also implicated BET proteins as regulators of cell identity and differentiation through their interactions with lineage-specific factors. However, the role of BET proteins has not yet been investigated in melanocyte differentiation. Melanocyte inducing transcription factor (MITF) is the master regulator of melanocyte differentiation, essential for pigmentation and melanocyte survival. In this study, we tested the hypothesis that BET proteins regulate melanocyte differentiation through interactions with MITF.

RESULTS

Here we show that chemical inhibition of BET proteins prevents differentiation of unpigmented melanoblasts into pigmented melanocytes and results in de-pigmentation of differentiated melanocytes. BET inhibition also slowed cell growth, without causing cell death, increasing the number of cells in G1. Transcriptional profiling revealed that BET inhibition resulted in decreased expression of pigment-specific genes, including many MITF targets. The expression of pigment-specific genes was also down-regulated in melanoma cells, but to a lesser extent. We found that RNAi depletion of the BET family members, bromodomain-containing protein 4 (BRD4) and bromodomain-containing protein 2 (BRD2) inhibited expression of two melanin synthesis enzymes, TYR and TYRP1. Both BRD4 and BRD2 were detected on melanocyte promoters surrounding MITF-binding sites, were associated with open chromatin structure, and promoted MITF binding to these sites. Furthermore, BRD4 and BRD2 physically interacted with MITF.

CONCLUSION

These findings indicate a requirement for BET proteins in the regulation of pigmentation and melanocyte differentiation. We identified changes in pigmentation specific gene expression that occur upon BET inhibition in melanoblasts, melanocytes, and melanoma cells.

An integrative view of the regulatory and transcriptional landscapes in mouse hematopoiesis

Thousands of epigenomic data sets have been generated in the past decade, but it is difficult for researchers to effectively use all the data relevant to their projects. Systematic integrative analysis can help meet this need, and the VISION project was established for validated systematic integration of epigenomic data in hematopoiesis. Here, we systematically integrated extensive data recording epigenetic features and transcriptomes from many sources, including individual laboratories and consortia, to produce a comprehensive view of the regulatory landscape of differentiating hematopoietic cell types in mouse. By using IDEAS as our integrative and discriminative epigenome annotation system, we identified and assigned epigenetic states simultaneously along chromosomes and across cell types, precisely and comprehensively. Combining nuclease accessibility and epigenetic states produced a set of more than 200,000 candidate cis-regulatory elements (cCREs) that efficiently capture enhancers and promoters. The transitions in epigenetic states of these cCREs across cell types provided insights into mechanisms of regulation, including decreases in numbers of active cCREs during differentiation of most lineages, transitions from poised to active or inactive states, and shifts in nuclease accessibility of CTCF-bound elements. Regression modeling of epigenetic states at cCREs and gene expression produced a versatile resource to improve selection of cCREs potentially regulating target genes. These resources are available from our VISION website to aid research in genomics and hematopoiesis.

Comparative structure-function analysis of bromodomain and extraterminal motif (BET) proteins in a gene-complementation system

The widely expressed bromodomain and extraterminal motif (BET) proteins bromodomain-containing protein 2 (BRD2), BRD3, and BRD4 are multifunctional transcriptional regulators that bind acetylated chromatin via their conserved tandem bromodomains. Small molecules that target BET bromodomains are being tested for various diseases but typically do not discern between BET family members. Genomic distributions and protein partners of BET proteins have been described, but the basis for differences in BET protein function within a given lineage remains unclear. By establishing a gene knockout-rescue system in a Brd2-null erythroblast cell line, here we compared a series of mutant and chimeric BET proteins for their ability to modulate cell growth, differentiation, and gene expression. We found that the BET N-terminal halves bearing the bromodomains convey marked differences in protein stability but do not account for specificity in BET protein function. Instead, when BET proteins were expressed at comparable levels, their specificity was largely determined by the C-terminal half. Remarkably, a chimeric BET protein comprising the N-terminal half of the structurally similar short BRD4 isoform (BRD4S) and the C-terminal half of BRD2 functioned similarly to intact BRD2. We traced part of the BRD2-specific activity to a previously uncharacterized short segment predicted to harbor a coiled-coil (CC) domain. Deleting the CC segment impaired BRD2's ability to restore growth and differentiation, and the CC region functioned in conjunction with the adjacent ET domain to impart BRD2-like activity onto BRD4S. In summary, our results identify distinct BET protein domains that regulate protein turnover and biological activities.

Upon microbial challenge, human neutrophils undergo rapid changes in nuclear architecture and chromatin folding to orchestrate an immediate inflammatory gene program

Differentiating neutrophils undergo large-scale changes in nuclear morphology. How such alterations in structure are established and modulated upon exposure to microbial agents is largely unknown. Here, we found that prior to encounter with bacteria, an armamentarium of inflammatory genes was positioned in a transcriptionally passive environment suppressing premature transcriptional activation. Upon microbial exposure, however, human neutrophils rapidly (<3 h) repositioned the ensemble of proinflammatory genes toward the transcriptionally permissive compartment. We show that the repositioning of genes was closely associated with the swift recruitment of cohesin across the inflammatory enhancer landscape, permitting an immediate transcriptional response upon bacterial exposure. We found that activated enhancers, marked by increased deposition of H3K27Ac, were highly enriched for cistromic elements associated with PU.1, CEBPB, TFE3, JUN, and FOSL2 occupancy. These data reveal how upon microbial challenge the cohesin machinery is recruited to an activated enhancer repertoire to instruct changes in chromatin folding, nuclear architecture, and to activate an inflammatory gene program.

Systematic integration of GATA transcription factors and epigenomes via IDEAS paints the regulatory landscape of hematopoietic cells

Members of the GATA family of transcription factors play key roles in the differentiation of specific cell lineages by regulating the expression of target genes. Three GATA factors play distinct roles in hematopoietic differentiation. In order to better understand how these GATA factors function to regulate genes throughout the genome, we are studying the epigenomic and transcriptional landscapes of hematopoietic cells in a model-driven, integrative fashion. We have formed the collaborative multi-lab VISION project to conduct ValIdated Systematic IntegratiON of epigenomic data in mouse and human hematopoiesis. The epigenomic data included nuclease accessibility in chromatin, CTCF occupancy, and histone H3 modifications for 20 cell types covering hematopoietic stem cells, multilineage progenitor cells, and mature cells across the blood cell lineages of mouse. The analysis used the Integrative and Discriminative Epigenome Annotation System (IDEAS), which learns all common combinations of features (epigenetic states) simultaneously in two dimensions-along chromosomes and across cell types. The result is a segmentation that effectively paints the regulatory landscape in readily interpretable views, revealing constitutively active or silent loci as well as the loci specifically induced or repressed in each stage and lineage. Nuclease accessible DNA segments in active chromatin states were designated candidate cis-regulatory elements in each cell type, providing one of the most comprehensive registries of candidate hematopoietic regulatory elements to date. Applications of VISION resources are illustrated for the regulation of genes encoding GATA1, GATA2, GATA3, and Ikaros. VISION resources are freely available from our website http://usevision.org.

OnTAD: hierarchical domain structure reveals the divergence of activity among TADs and boundaries

The spatial organization of chromatin in the nucleus has been implicated in regulating gene expression. Maps of high-frequency interactions between different segments of chromatin have revealed topologically associating domains (TADs), within which most of the regulatory interactions are thought to occur. TADs are not homogeneous structural units but appear to be organized into a hierarchy. We present OnTAD, an optimized nested TAD caller from Hi-C data, to identify hierarchical TADs. OnTAD reveals new biological insights into the role of different TAD levels, boundary usage in gene regulation, the loop extrusion model, and compartmental domains. OnTAD is available at https://github.com/anlin00007/OnTAD.

Chromatin structure dynamics during the mitosis-to-G1 phase transition

Features of higher-order chromatin organization-such as A/B compartments, topologically associating domains and chromatin loops-are temporarily disrupted during mitosis1,2. Because these structures are thought to influence gene regulation, it is important to understand how they are re-established after mitosis. Here we examine the dynamics of chromosome reorganization by Hi-C after mitosis in highly purified, synchronous mouse erythroid cell populations. We observed rapid establishment of A/B compartments, followed by their gradual intensification and expansion. Contact domains form from the 'bottom up'-smaller subTADs are formed initially, followed by convergence into multi-domain TAD structures. CTCF is partially retained on mitotic chromosomes and immediately resumes full binding in ana/telophase. By contrast, cohesin is completely evicted from mitotic chromosomes and regains focal binding at a slower rate. The formation of CTCF/cohesin co-anchored structural loops follows the kinetics of cohesin positioning. Stripe-shaped contact patterns-anchored by CTCF-grow in length, which is consistent with a loop-extrusion process after mitosis. Interactions between cis-regulatory elements can form rapidly, with rates exceeding those of CTCF/cohesin-anchored contacts. Notably, we identified a group of rapidly emerging transient contacts between cis-regulatory elements in ana/telophase that are dissolved upon G1 entry, co-incident with the establishment of inner boundaries or nearby interfering chromatin loops. We also describe the relationship between transcription reactivation and architectural features. Our findings indicate that distinct but mutually influential forces drive post-mitotic chromatin reconfiguration.

Joint inference and alignment of genome structures enables characterization of compartment-independent reorganization across cell types

BACKGROUND

Comparisons of Hi-C data sets between cell types and conditions have revealed differences in topologically associated domains (TADs) and A/B compartmentalization, which are correlated with differences in gene regulation. However, previous comparisons have focused on known forms of 3D organization while potentially neglecting other functionally relevant differences. We aimed to create a method to quantify all locus-specific differences between two Hi-C data sets.

RESULTS

We developed MultiMDS to jointly infer and align 3D chromosomal structures from two Hi-C data sets, thereby enabling a new way to comprehensively quantify relocalization of genomic loci between cell types. We demonstrate this approach by comparing Hi-C data across a variety of cell types. We consistently find relocalization of loci with minimal difference in A/B compartment score. For example, we identify compartment-independent relocalizations between GM12878 and K562 cells that involve loci displaying enhancer-associated histone marks in one cell type and polycomb-associated histone marks in the other.

CONCLUSIONS

MultiMDS is the first tool to identify all loci that relocalize between two Hi-C data sets. Our method can identify 3D localization differences that are correlated with cell-type-specific regulatory activities and which cannot be identified using other methods.

CGGBP1 regulates CTCF occupancy at repeats

Background

CGGBP1 is a repeat-binding protein with diverse functions in the regulation of gene expression, cytosine methylation, repeat silencing and genomic integrity. CGGBP1 has also been identified as a cooperator of histone-modifying enzymes and as a component of CTCF-containing complexes that regulate the enhancer–promoter looping. CGGBP1–CTCF cross talk in chromatin regulation has been hitherto unknown.

Results

Here, we report that the occupancy of CTCF at repeats depends on CGGBP1. Using ChIP-sequencing for CTCF, we describe its occupancy at repetitive DNA. Our results show that endogenous level of CGGBP1 ensures CTCF occupancy preferentially on repeats over canonical CTCF motifs. By combining CTCF ChIP-sequencing results with ChIP sequencing for three different kinds of histone modifications (H3K4me3, H3K9me3 and H3K27me3), we show that the CGGBP1-dependent repeat-rich CTCF-binding sites regulate histone marks in flanking regions.

Conclusion

CGGBP1 affects the pattern of CTCF occupancy. Our results posit CGGBP1 as a regulator of CTCF and its binding sites in interspersed repeats.

The Cornelia de Lange Syndrome-associated factor NIPBL interacts with BRD4 ET domain for transcription control of a common set of genes

Mutations in NIPBL are the major cause of Cornelia de Lange Syndrome (CdLS). NIPBL is the cohesin-loading factor and has recently been associated with the BET (bromodomains and extra-terminal (ET) domain) proteins BRD2 and BRD4. Related to this, a CdLS-like phenotype has been described associated to BRD4 mutations. Here, we show direct interaction of NIPBL with different BET members in yeast, and selective interaction with BRD4 in cells, being the ET domain involved in the interaction. To understand the relationship between NIPBL and BET proteins, we have performed RNA-Seq expression analysis following depletion of the different proteins. Results indicate that genes regulated by NIPBL largely overlap with those regulated by BRD4 but not with those regulated by BRD2. ChIP-Seq analysis indicates preferential NIPBL occupancy at promoters, and knockdown experiments show mutual stabilization of NIPBL and BRD4 on co-regulated promoters. Moreover, human fibroblasts from CdLS probands with mutations in NIPBL show reduced BRD4 at co-occupied promoters. Functional analysis in vivo, using mutants of Drosophila melanogaster, confirmed the genetic interaction between Nipped-B and fs(1)h, the orthologs of human NIPBL and BRD4, respectively. Thus, we provide evidence for NIPBL and BRD4 cooperation in transcriptional regulation, which should contribute to explain the recently observed CdLS-like phenotype associated with BRD4 mutations.

LADL: light-activated dynamic looping for endogenous gene expression control

Mammalian genomes are folded into tens of thousands of long-range looping interactions. The cause-and-effect relationship between looping and genome function is poorly understood, and the extent to which loops are dynamic on short time scales remains an unanswered question. Here, we engineer a new class of synthetic architectural proteins for directed rearrangement of the three-dimensional genome using blue light. We target our light-activated-dynamic-looping (LADL) system to two genomic anchors with CRISPR guide RNAs and induce their spatial colocalization via light-induced heterodimerization of cryptochrome 2 and a dCas9-CIBN fusion protein. We apply LADL to redirect a stretch enhancer (SE) away from its endogenous Klf4 target gene and to the Zfp462 promoter. Using single-molecule RNA-FISH, we demonstrate that de novo formation of the Zfp462-SE loop correlates with a modest increase in Zfp462 expression. LADL facilitates colocalization of genomic loci without exogenous chemical cofactors and will enable future efforts to engineer reversible and oscillatory loops on short time scales.

Pharmacological Modulation of Transcriptional Coregulators in Cancer

Upon binding of transcription factors to cis-regulatory DNA sequences, transcriptional coregulators are required for the activation or suppression of chromatin-dependent transcriptional signaling. These coregulators are frequently implicated in oncogenesis via causal roles in dysregulated, malignant transcriptional control and represent one of the fastest-growing target classes in small-molecule drug discovery. However, challenges in targeting coregulators include identifying evidence of cancer-specific genetic dependency, matching the pharmacologically addressable protein fold to a functional role in disease pathology, and achieving the necessary selectivity to exploit a given genetic dependency. We discuss here how recent trends in cancer pharmacology have confronted these challenges, positioning coregulators as tractable targets in the development of new cancer therapies.

Identification of Cancer Drivers at CTCF Insulators in 1,962 Whole Genomes

Recent studies have shown that mutations at non-coding elements, such as promoters and enhancers, can act as cancer drivers. However, an important class of non-coding elements, namely CTCF insulators, has been overlooked in the previous driver analyses. We used insulator annotations from CTCF and cohesin ChIA-PET and analyzed somatic mutations in 1,962 whole genomes from 21 cancer types. Using the heterogeneous patterns of transcription-factor-motif disruption, functional impact, and recurrence of mutations, we developed a computational method that revealed 21 insulators showing signals of positive selection. In particular, mutations in an insulator in multiple cancer types, including 16% of melanoma samples, are associated with TGFB1 up-regulation. Using CRISPR-Cas9, we find that alterations at two of the most frequently mutated regions in this insulator increase cell growth by 40%-50%, supporting the role of this boundary element as a cancer driver. Thus, our study reveals several CTCF insulators as putative cancer drivers.

Neuronal Transdifferentiation Potential of Human Mesenchymal Stem Cells from Neonatal and Adult Sources by a Small Molecule Cocktail

Human mesenchymal stem cells (MSCs) are good candidates for brain cell replacement strategies and have already been used as adjuvant treatments in neurological disorders. MSCs can be obtained from many different sources, and the present study compares the potential of neuronal transdifferentiation in MSCs from adult and neonatal sources (Wharton's jelly (WhJ), dental pulp (DP), periodontal ligament (PDL), gingival tissue (GT), dermis (SK), placenta (PLAC), and umbilical cord blood (UCB)) with a protocol previously tested in bone marrow- (BM-) MSCs consisting of a cocktail of six small molecules: I-BET151, CHIR99021, forskolin, RepSox, Y-27632, and dbcAMP (ICFRYA). Neuronal morphology and the presence of cells positive for neuronal markers (TUJ1 and MAP2) were considered attributes of neuronal induction. The ICFRYA cocktail did not induce neuronal features in WhJ-MSCs, and these features were only partial in the MSCs from dental tissues, SK-MSCs, and PLAC-MSCs. The best response was found in UCB-MSCs, which was comparable to the response of BM-MSCs. The addition of neurotrophic factors to the ICFRYA cocktail significantly increased the number of cells with complex neuron-like morphology and increased the number of cells positive for mature neuronal markers in BM- and UCB-MSCs. The neuronal cells generated from UCB-MSCs and BM-MSCs showed increased reactivity of the neuronal genes TUJ1, MAP2, NF-H, NCAM, ND1, TAU, ENO2, GABA, and NeuN as well as down- and upregulation of MSC and neuronal genes, respectively. The present study showed marked differences between the MSCs from different sources in response to the transdifferentiation protocol used here. These results may contribute to identifying the best source of MSCs for potential cell replacement therapies.

Epigenetics, DNA Organization, and Inflammatory Bowel Disease

Inflammatory bowel diseases (IBDs) are chronic inflammatory disorders affecting the gastrointestinal tract. The incidence of IBD is increasing, with more cases occurring in developed countries. Multiple factors such as genetics, environmental changes, gut microbiota, and immune abnormalities have been associated with development of IBD. In recent years, it has become increasingly apparent that epigenetic modifications of chromatin and the manner in which chromatin is organized in the nucleus are additionally important elements that can influence responses induced by the factors described above, and may therefore contribute to the onset and pathogenesis of IBD. Epigenetics and chromatin organization regulate diverse functions that include maintenance of homeostasis in the intestinal epithelium, the development and differentiation of immune cells, and modulation of responses generated by the immune system to defend against potential pathogens. Furthermore, changes in epigenetic chromatin marks and in chromatin organization have now been linked to differential gene expression in IBD patient cells. Although direct evidence for a role of histone modifications in IBD is currently very limited, in this review, we summarize the links between various epigenetic modifications, the proteins that catalyze or recognize these modifications, and the development or progression of IBD in human and experimental IBD. We also discuss how epigenetics influence the organization of DNA contacts to regulate gene expression and the implications this may have for diagnosing and treating IBD.

Generative modeling of multi-mapping reads with mHi-C advances analysis of Hi-C studies

Current Hi-C analysis approaches are unable to account for reads that align to multiple locations, and hence underestimate biological signal from repetitive regions of genomes. We developed and validated mHi-C, a multi-read mapping strategy to probabilistically allocate Hi-C multi-reads. mHi-C exhibited superior performance over utilizing only uni-reads and heuristic approaches aimed at rescuing multi-reads on benchmarks. Specifically, mHi-C increased the sequencing depth by an average of 20% resulting in higher reproducibility of contact matrices and detected interactions across biological replicates. The impact of the multi-reads on the detection of significant interactions is influenced marginally by the relative contribution of multi-reads to the sequencing depth compared to uni-reads, cis-to-trans ratio of contacts, and the broad data quality as reflected by the proportion of mappable reads of datasets. Computational experiments highlighted that in Hi-C studies with short read lengths, mHi-C rescued multi-reads can emulate the effect of longer reads. mHi-C also revealed biologically supported bona fide promoter-enhancer interactions and topologically associating domains involving repetitive genomic regions, thereby unlocking a previously masked portion of the genome for conformation capture studies.

The interdependence of gene-regulatory elements and the 3D genome

Imaging studies, high-resolution chromatin conformation maps, and genome-wide occupancy data of architectural proteins have revealed that genome topology is tightly intertwined with gene expression. Cross-talk between gene-regulatory elements is often organized within insulated neighborhoods, and regulatory cues that induce transcriptional changes can reshape chromatin folding patterns and gene positioning within the nucleus. The cause-consequence relationship of genome architecture and gene expression is intricate, and its molecular mechanisms are under intense investigation. Here, we review the interdependency of transcription and genome organization with emphasis on enhancer-promoter contacts in gene regulation.

Integrative view on how erythropoietin signaling controls transcription patterns in erythroid cells

PURPOSE OF REVIEW

Erythropoietin (EPO) is necessary and sufficient to trigger dynamic transcriptional patterns that drive the differentiation of erythroid precursor cells into mature, enucleated red cells. Because the molecular cloning and Food and Drug Administration approval for the therapeutic use of EPO over 30 years ago, a detailed understanding of how EPO works has advanced substantially. Yet, the precise epigenetic and transcriptional mechanisms by which EPO signaling controls erythroid expression patterns remains poorly understood. This review focuses on the current state of erythroid biology in regards to EPO signaling from human genetics and functional genomics perspectives.

RECENT FINDINGS

The goal of this review is to provide an integrative view of the gene regulatory underpinnings for erythroid expression patterns that are dynamically shaped during erythroid differentiation. Here, we highlight vignettes connecting recent insights into a genome-wide association study linking an EPO mutation to anemia, a study linking EPO-signaling to signal transducer and activator of transcription 5 (STAT5) chromatin occupancy and enhancers, and studies that examine the molecular mechanisms driving topological chromatin organization in erythroid cells.

SUMMARY

The genetic, epigenetic, and gene regulatory mechanisms underlying how hormone signal transduction influences erythroid gene expression remains only partly understood. A detailed understanding of these molecular pathways and how they intersect with one another will provide the basis for novel strategies to treat anemia and potentially other hematological diseases. As new regulators and signal transducers of EPO-signaling continue to emerge, new clinically relevant targets may be identified that improve the specificity and effectiveness of EPO therapy.

Disruption of CTCF-YY1-dependent looping of the human papillomavirus genome activates differentiation-induced viral oncogene transcription

The complex life cycle of oncogenic human papillomavirus (HPV) initiates in undifferentiated basal epithelial keratinocytes where expression of the E6 and E7 oncogenes is restricted. Upon epithelial differentiation, E6/E7 transcription is increased through unknown mechanisms to drive cellular proliferation required to support virus replication. We report that the chromatin-organising CCCTC-binding factor (CTCF) promotes the formation of a chromatin loop in the HPV genome that epigenetically represses viral enhancer activity controlling E6/E7 expression. CTCF-dependent looping is dependent on the expression of the CTCF-associated Yin Yang 1 (YY1) transcription factor and polycomb repressor complex (PRC) recruitment, resulting in trimethylation of histone H3 at lysine 27. We show that viral oncogene up-regulation during cellular differentiation results from YY1 down-regulation, disruption of viral genome looping, and a loss of epigenetic repression of viral enhancer activity. Our data therefore reveal a key role for CTCF-YY1-dependent looping in the HPV life cycle and identify a regulatory mechanism that could be disrupted in HPV carcinogenesis.

JQ1 affects BRD2-dependent and independent transcription regulation without disrupting H4-hyperacetylated chromatin states

The bromodomain and extra-terminal domain (BET) proteins are promising drug targets for cancer and immune diseases. However, BET inhibition effects have been studied more in the context of bromodomain-containing protein 4 (BRD4) than BRD2, and the BET protein association to histone H4-hyperacetylated chromatin is not understood at the genome-wide level. Here, we report transcription start site (TSS)-resolution integrative analyses of ChIP-seq and transcriptome profiles in human non-small cell lung cancer (NSCLC) cell line H23. We show that di-acetylation at K5 and K8 of histone H4 (H4K5acK8ac) co-localizes with H3K27ac and BRD2 in the majority of active enhancers and promoters, where BRD2 has a stronger association with H4K5acK8ac than H3K27ac. Although BET inhibition by JQ1 led to complete reduction of BRD2 binding to chromatin, only local changes of H4K5acK8ac levels were observed, suggesting that recruitment of BRD2 does not influence global histone H4 hyperacetylation levels. This finding supports a model in which recruitment of BET proteins via histone H4 hyperacetylation is predominant over hyperacetylation of histone H4 by BET protein-associated acetyltransferases. In addition, we found that a remarkable number of BRD2-bound genes, including MYC and its downstream target genes, were transcriptionally upregulated upon JQ1 treatment. Using BRD2-enriched sites and transcriptional activity analysis, we identified candidate transcription factors potentially involved in the JQ1 response in BRD2-dependent and -independent manner.

Dynamic Nature of CTCF Tandem 11 Zinc Fingers in Multivalent Recognition of DNA As Revealed by NMR Spectroscopy

The 11 zinc fingers (ZFs) of the transcription factor CTCF play a versatile role in the regulation of gene expression. CTCF binds to numerous genomic sites to form chromatin loops and topologically associated domains and thus mediates the 3D architecture of chromatin. Although CTCF inter-ZF plasticity is essential for the recognition of multiple genomic sites, the dynamic nature of its 11 ZFs remains unknown. We assigned the chemical shifts of the CTCF ZFs 1-11 and solved the solution structures of each ZF. NMR backbone dynamics, residual dipolar couplings, and small-angle X-ray scattering experiments suggest a high inter-ZF plasticity of the free-form ZFs 1-11. As exemplified by two different protocadherin DNA sequences, the titration of DNAs to ¹⁵N-labeled CTCF ZFs 1-11 enabled systematic mapping of binding of CTCF ZFs to various chromatin sites. Our work paves the way for illustrating the molecular basis of the versatile DNA recognized by CTCF and has interesting implications for its conformational transition during DNA binding.

Forces driving the three-dimensional folding of eukaryotic genomes

The last decade has radically renewed our understanding of higher order chromatin folding in the eukaryotic nucleus. As a result, most current models are in support of a mostly hierarchical and relatively stable folding of chromosomes dividing chromosomal territories into A‐ (active) and B‐ (inactive) compartments, which are then further partitioned into topologically associating domains (TADs), each of which is made up from multiple loops stabilized mainly by the CTCF and cohesin chromatin‐binding complexes. Nonetheless, the structure‐to‐function relationship of eukaryotic genomes is still not well understood. Here, we focus on recent work highlighting the biophysical and regulatory forces that contribute to the spatial organization of genomes, and we propose that the various conformations that chromatin assumes are not so much the result of a linear hierarchy, but rather of both converging and conflicting dynamic forces that act on it.

Functional roles of CTCF in breast cancer

CTCF, Zinc-finger protein, has been identified as a multifunctional transcription factor that regulates gene expression through various mechanisms, including recruitment of other co-activators and binding to promoter regions of target genes. Furthermore, it has been proposed to be an insulator protein that contributes to the establishment of functional three-dimensional chromatin structures. It can disrupt transcription through blocking the connection between an enhancer and a promoter. Previous studies revealed that the onset of various diseases, including breast cancer, could be attributed to the aberrant expression of CTCF itself or one or more of its target genes. In this review, we will describe molecular dysfunction involving CTCF that induces tumorigenesis and summarize the functional roles of CTCF in breast cancer.

The BRD3 ET domain recognizes a short peptide motif through a mechanism that is conserved across chromatin remodelers and transcriptional regulators

Members of the bromodomain and extra-terminal domain (BET) family of proteins (bromodomain-containing (BRD) 2, 3, 4, and T) are widely expressed and highly conserved regulators of gene expression in eukaryotes. These proteins have been intimately linked to human disease, and more than a dozen clinical trials are currently underway to test BET-protein inhibitors as modulators of cancer. However, although it is clear that these proteins use their bromodomains to bind both histones and transcription factors bearing acetylated lysine residues, the molecular mechanisms by which BET family proteins regulate gene expression are not well defined. In particular, the functions of the other domains such as the ET domain have been less extensively studied. Here, we examine the properties of the ET domain of BRD3 as a protein/protein interaction module. Using a combination of pulldown and biophysical assays, we demonstrate that BRD3 binds to a range of chromatin-remodeling complexes, including the NuRD, BAF, and INO80 complexes, via a short linear "KIKL" motif in one of the complex subunits. NMR-based structural analysis revealed that, surprisingly, this mode of interaction is shared by the AF9 and ENL transcriptional coregulators that contain an acetyl-lysine-binding YEATS domain and regulate transcriptional elongation. This observation establishes a functional commonality between these two families of cancer-related transcriptional regulators. In summary, our data provide insight into the mechanisms by which BET family proteins might link chromatin acetylation to transcriptional outcomes and uncover an unexpected functional similarity between BET and YEATS family proteins.

Emerging Evidence of Chromosome Folding by Loop Extrusion

Chromosome organization poses a remarkable physical problem with many biological consequences: How can molecular interactions between proteins at the nanometer scale organize micron-long chromatinized DNA molecules, insulating or facilitating interactions between specific genomic elements? The mechanism of active loop extrusion holds great promise for explaining interphase and mitotic chromosome folding, yet remains difficult to assay directly. We discuss predictions from our polymer models of loop extrusion with barrier elements and review recent experimental studies that provide strong support for loop extrusion, focusing on perturbations to CTCF and cohesin assayed via Hi-C in interphase. Finally, we discuss a likely molecular mechanism of loop extrusion by structural maintenance of chromosomes complexes.

Enhancer dysfunction in leukemia

Hematopoietic cancers are often initiated by deregulation of the transcriptional machinery. Prominent among such regulators are the sequence-specific DNA-binding transcription factors (TFs), which bind to enhancer and promoter elements in the genome to control gene expression through the recruitment of cofactors. Remarkably, perturbing the function of even a single TF or cofactor can modulate the active enhancer landscape of a cell; conversely, knowledge of the enhancer configuration can be used to discover functionally important TFs in a given cellular process. Our expanding insight into enhancer function can be attributed to the emergence of genome-scale measurements of enhancer activity, which can be applied to virtually any cell type to expose regulatory mechanisms. Such approaches are beginning to reveal the abnormal enhancer configurations present in cancer cells, thereby providing a framework for understanding how transcriptional dysregulation can lead to malignancy. Here, we review the evidence for alterations in enhancer landscapes contributing to the pathogenesis of leukemia, a malignancy in which enhancer-binding proteins and enhancer DNA itself are altered via genetic mutation. We will also highlight examples of small molecules that reprogram the enhancer landscape of leukemia cells in association with therapeutic benefit.

CRISPR-based strategies for studying regulatory elements and chromatin structure in mammalian gene control

The development of high-throughput methods has enabled the genome-wide identification of putative regulatory elements in a wide variety of mammalian cells at an unprecedented resolution. Extensive genomic studies have revealed the important role of regulatory elements and genetic variation therein in disease formation and risk. In most cases, there is only correlative evidence for the roles of these elements and non-coding changes within these elements in pathogenesis. With the advent of genome- and epigenome-editing tools based on the CRISPR technology, it is now possible to test the functional relevance of the regulatory elements and alterations on a genomic scale. Here, we review the various CRISPR-based strategies that have been developed to functionally validate the candidate regulatory elements in mammals as well as the non-coding genetic variants found to be associated with human disease. We also discuss how these synthetic biology tools have helped to elucidate the role of three-dimensional nuclear architecture and higher-order chromatin organization in shaping functional genome and controlling gene expression.

Lyar-Mediated Recruitment of Brd2 to the Chromatin Attenuates Nanog Downregulation Following Induction of Differentiation

During development, cellular differentiation programs need tight regulation for proper display of the activity of multiple factors in time and space. Chromatin adaptors of the BET family (Brd2, Brd3, Brd4 and Brdt in vertebrates) are transcription co-regulators tightly associated with the progression of the cell cycle. A key question regarding their function is whether they work as part of the general transcription machinery or, on the contrary, they are precisely recruited to the chromatin through specific transcription factors. Here, we report the selective recruitment of Brd2 to the chromatin by the transcription factor Lyar. We show that Lyar downregulation results in Brd2 dissociation from a number of promoters studied. On the contrary, dissociation of BET proteins from the chromatin has no effect on Lyar occupancy. Under differentiation conditions, the absence of Lyar leads to impaired downregulation of the pluripotency gene Nanog, with concomitant reduction in the upregulation of differentiation markers. Interestingly, following the induction of differentiation, Brd2 depletion exhibits the same effects as expressing a truncated Lyar molecule lacking the Brd2 interacting domain. Both approaches result in stronger Nanog repression, indicating that Lyar-mediated recruitment of Brd2 moderates Nanog downregulation when differentiation is triggered. Moreover, expression of truncated Lyar leads to impaired differentiation and increased apoptosis. Thus, Lyar-mediated recruitment of Brd2 would participate in preserving a proper timing for Nanog silencing ensuring the appropriate establishment of the differentiation program.

Selective degradation of BET proteins with dBET1, a proteolysis-targeting chimera, potently reduces pro-inflammatory responses in lipopolysaccharide-activated microglia

Bromodomain and extraterminal (BET) proteins are essential to pro-inflammatory gene transcription. The BET family proteins, BRD2, BRD3, BRD4, and testis-specific BRDT, couple chromatin remodeling to gene transcription, acting as histone acetyltransferases, scaffolds for transcription complexes, and markers of histone acetylation. To initiate an inflammatory response, cells undergo de novo gene transcription requiring histone-modifying proteins to make DNA wrapped around histones more or less readily available to transcription complexes. Because BET proteins are the gatekeepers of nuclear factor-κB (NF-κB)-dependent gene transcription, we hypothesized that degradation of BET proteins, particularly BRD2 and BRD4, with the proteolysis-targeting chimera (PROTAC) dBET1 would dampen the pro-inflammatory response in microglia subjected to lipopolysaccharide (LPS) challenge. Degradation of BRD2 and BRD4 was associated with significantly reduced expression of several pro-inflammatory genes: inducible nitric oxide synthase (iNOS), cyclooxygenase-2 (COX-2), interleukin (IL)-1β, tumor necrosis factor-a (TNF-α), IL-6, chemokine (C-C motif) ligand 2 (CCL2), and matrix metalloproteinase-9 (MMP-9). This is the first study showing that dBET1-mediated targeted degradation of BET proteins robustly dampens pro-inflammatory responses in LPS-stimulated microglia. These data suggest that BET degradation with dBET1 will likely reduce expression of pro-inflammatory genes in in vivo neuroinflammatory models associated with microglial/immune cell activation.

Single-cell systems biology: probing the basic unit of information flow

Gene expression varies across cells in a population or a tissue. This heterogeneity has come into sharp focus in recent years through developments in new imaging and sequencing technologies. However, our ability to measure variation has outpaced our ability to interpret it. Much of the variability may arise from random effects occurring in the processes of gene expression (transcription, RNA processing and decay, translation). The molecular basis of these effects is largely unknown. Likewise, a functional role of this variability in growth, differentiation and disease has only been elucidated in a few cases. In this review, we highlight recent experimental and theoretical advances for measuring and analyzing stochastic variation.

The Role of Bromodomain and Extraterminal Motif (BET) Proteins in Chromatin Structure

Bromodomain and extraterminal motif (BET) proteins have been widely investigated for their roles in gene regulation and their potential as therapeutic targets in cancer. Pharmacologic BET inhibitors target the conserved bromodomain-acetyllysine interaction and do not distinguish between BRD2, BRD3, and BRD4. Thus, comparatively little is known regarding the distinct roles played by individual family members, as well as the underlying mechanisms that drive the transcriptional effects of BET inhibitors. Here we review studies regarding the contributions of BET proteins to genome structure and function, including recent work identifying a role for BRD2 as a component of functional and physical chromatin domain boundaries. We also discuss directions of future studies aimed at providing insights into broader architectural functions of BET proteins and their roles in chromatin domain boundary formation.

Harnessing BET Inhibitor Sensitivity Reveals AMIGO2 as a Melanoma Survival Gene

Bromodomain and extraterminal domain inhibitors (BETi) represent promising therapeutic agents for metastatic melanoma, yet their mechanism of action remains unclear. Here we interrogated the transcriptional effects of BETi and identified AMIGO2, a transmembrane molecule, as a BET target gene essential for melanoma cell survival. AMIGO2 is upregulated in melanoma cells and tissues compared to human melanocytes and nevi, and AMIGO2 silencing in melanoma cells induces G1/S arrest followed by apoptosis. We identified the pseudokinase PTK7 as an AMIGO2 interactor whose function is regulated by AMIGO2. Epigenomic profiling and genome editing revealed that AMIGO2 is regulated by a melanoma-specific BRD2/4-bound promoter and super-enhancer configuration. Upon BETi treatment, BETs are evicted from these regulatory elements, resulting in AMIGO2 silencing and changes in PTK7 proteolytic processing. Collectively, this study uncovers mechanisms underlying the therapeutic effects of BETi in melanoma and reveals the AMIGO2-PTK7 axis as a targetable pathway for metastatic melanoma.

Between form and function: the complexity of genome folding

It has been known for over a century that chromatin is not randomly distributed within the nucleus. However, the question of how DNA is folded and the influence of such folding on nuclear processes remain topics of intensive current research. A longstanding, unanswered question is whether nuclear organization is simply a reflection of nuclear processes such as transcription and replication, or whether chromatin is folded by independent mechanisms and this per se encodes function? Evidence is emerging that both may be true. Here, using the α-globin gene cluster as an illustrative model, we provide an overview of the most recent insights into the layers of genome organization across different scales and how this relates to gene activity.

Comparative analysis of three-dimensional chromosomal architecture identifies a novel fetal hemoglobin regulatory element

In this study, Huang et al. compared the chromosomal architectures of fetal and adult human erythroblasts and found that, globally, chromatin structures and compartments A/B are highly similar at both developmental stages. Their results uncover a new critical regulatory region as a potential target for therapeutic genome editing for hemoglobinopathies and highlight the power of chromosome conformation analysis in discovering new cis control elements.

Chromatin structure is tightly intertwined with transcription regulation. Here we compared the chromosomal architectures of fetal and adult human erythroblasts and found that, globally, chromatin structures and compartments A/B are highly similar at both developmental stages. At a finer scale, we detected distinct folding patterns at the developmentally controlled β-globin locus. Specifically, new fetal stage-specific contacts were uncovered between a region separating the fetal (γ) and adult (δ and β) globin genes (encompassing the HBBP1 and BGLT3 noncoding genes) and two distal chromosomal sites (HS5 and 3′HS1) that flank the locus. In contrast, in adult cells, the HBBP1–BGLT3 region contacts the embryonic ε-globin gene, physically separating the fetal globin genes from the enhancer (locus control region [LCR]). Deletion of the HBBP1 region in adult cells alters contact landscapes in ways more closely resembling those of fetal cells, including increased LCR–γ-globin contacts. These changes are accompanied by strong increases in γ-globin transcription. Notably, the effects of HBBP1 removal on chromatin architecture and gene expression closely mimic those of deleting the fetal globin repressor BCL11A, implicating BCL11A in the function of the HBBP1 region. Our results uncover a new critical regulatory region as a potential target for therapeutic genome editing for hemoglobinopathies and highlight the power of chromosome conformation analysis in discovering new cis control elements.