Identifying clusters of cis-regulatory elements underpinning TAD structures and lineage-specific regulatory networks

TULIPs decorate the three-dimensional genome of PFA ependymoma

Posterior fossa group A (PFA) ependymoma is a lethal brain cancer diagnosed in infants and young children. The lack of driver events in the PFA linear genome led us to search its 3D genome for characteristic features. Here, we reconstructed 3D genomes from diverse childhood tumor types and uncovered a global topology in PFA that is highly reminiscent of stem and progenitor cells in a variety of human tissues. A remarkable feature exclusively present in PFA are type B ultra long-range interactions in PFAs (TULIPs), regions separated by great distances along the linear genome that interact with each other in the 3D nuclear space with surprising strength. TULIPs occur in all PFA samples and recur at predictable genomic coordinates, and their formation is induced by expression of EZHIP. The universality of TULIPs across PFA samples suggests a conservation of molecular principles that could be exploited therapeutically.

In vivo CRISPR screens identify a dual function of MEN1 in regulating tumor-microenvironment interactions

Functional genomic screens in two-dimensional cell culture models are limited in identifying therapeutic targets that influence the tumor microenvironment. By comparing targeted CRISPR-Cas9 screens in a two-dimensional culture with xenografts derived from the same cell line, we identified MEN1 as the top hit that confers differential dropout effects in vitro and in vivo. MEN1 knockout in multiple solid cancer types does not impact cell proliferation in vitro but significantly promotes or inhibits tumor growth in immunodeficient or immunocompetent mice, respectively. Mechanistically, MEN1 knockout redistributes MLL1 chromatin occupancy, increasing H3K4me3 at repetitive genomic regions, activating double-stranded RNA expression and increasing neutrophil and CD8+ T cell infiltration in immunodeficient and immunocompetent mice, respectively. Pharmacological inhibition of the menin-MLL interaction reduces tumor growth in a CD8+ T cell-dependent manner. These findings reveal tumor microenvironment-dependent oncogenic and tumor-suppressive functions of MEN1 and provide a rationale for targeting MEN1 in solid cancers.

Pushing the TAD boundary: Decoding insulator codes of clustered CTCF sites in 3D genomes

Topologically associating domain (TAD) boundaries are the flanking edges of TADs, also known as insulated neighborhoods, within the 3D structure of genomes. A prominent feature of TAD boundaries in mammalian genomes is the enrichment of clustered CTCF sites often with mixed orientations, which can either block or facilitate enhancer-promoter (E-P) interactions within or across distinct TADs, respectively. We will discuss recent progress in the understanding of fundamental organizing principles of the clustered CTCF insulator codes at TAD boundaries. Specifically, both inward- and outward-oriented CTCF sites function as topological chromatin insulators by asymmetrically blocking improper TAD-boundary-crossing cohesin loop extrusion. In addition, boundary stacking and enhancer clustering facilitate long-distance E-P interactions across multiple TADs. Finally, we provide a unified mechanism for RNA-mediated TAD boundary function via R-loop formation for both insulation and facilitation. This mechanism of TAD boundary formation and insulation has interesting implications not only on how the 3D genome folds in the Euclidean nuclear space but also on how the specificity of E-P interactions is developmentally regulated.

Permeable TAD boundaries and their impact on genome-associated functions

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

The developmental and evolutionary characteristics of transcription factor binding site clustered regions based on an explainable machine learning model

Gene expression is temporally and spatially regulated by the interaction of transcription factors (TFs) and cis-regulatory elements (CREs). The uneven distribution of TF binding sites across the genome poses challenges in understanding how this distribution evolves to regulate spatio-temporal gene expression and consequent heritable phenotypic variation. In this study, chromatin accessibility profiles and gene expression profiles were collected from several species including mammals (human, mouse, bovine), fish (zebrafish and medaka), and chicken. Transcription factor binding sites clustered regions (TFCRs) at different embryonic stages were characterized to investigate regulatory evolution. The study revealed dynamic changes in TFCR distribution during embryonic development and species evolution. The synchronization between TFCR complexity and gene expression was assessed across species using RegulatoryScore. Additionally, an explainable machine learning model highlighted the importance of the distance between TFCR and promoter in the coordinated regulation of TFCRs on gene expression. Our results revealed the developmental and evolutionary dynamics of TFCRs during embryonic development from fish, chicken to mammals. These data provide valuable resources for exploring the relationship between transcriptional regulation and phenotypic differences during embryonic development.

Multi-feature clustering of CTCF binding creates robustness for loop extrusion blocking and Topologically Associating Domain boundaries

Topologically Associating Domains (TADs) separate vertebrate genomes into insulated regulatory neighborhoods that focus genome-associated processes. TADs are formed by Cohesin-mediated loop extrusion, with many TAD boundaries consisting of clustered binding sites of the CTCF insulator protein. Here we determine how this clustering of CTCF binding contributes to the blocking of loop extrusion and the insulation between TADs. We identify enrichment of three features of CTCF binding at strong TAD boundaries, consisting of strongly bound and closely spaced CTCF binding peaks, with a further enrichment of DNA-binding motifs within these peaks. Using multi-contact Nano-C analysis in cells with normal and perturbed CTCF binding, we establish that individual CTCF binding sites contribute to the blocking of loop extrusion, but in an incomplete manner. When clustered, individual CTCF binding sites thus create a stepwise insulation between neighboring TADs. Based on these results, we propose a model whereby multiple instances of temporal loop extrusion blocking create strong insulation between TADs.

Genomic profiling of HIV-1 integration in microglia cells links viral integration to the topologically associated domains

HIV-1 encounters the hierarchically organized host chromatin to stably integrate and persist in anatomically distinct latent reservoirs. The contribution of genome organization in HIV-1 infection has been largely understudied across different HIV-1 targets. Here, we determine HIV-1 integration sites (ISs), associate them with chromatin and expression signatures at different genomic scales in a microglia cell model, and profile them together with the primary T cell reservoir. HIV-1 insertions into introns of actively transcribed genes with IS hotspots in genic and super-enhancers, characteristic of blood cells, are maintained in the microglia cell model. Genome organization analysis reveals dynamic CCCTC-binding factor (CTCF) clusters in cells with active and repressed HIV-1 transcription, whereas CTCF removal impairs viral integration. We identify CTCF-enriched topologically associated domain (TAD) boundaries with signatures of transcriptionally active chromatin as HIV-1 integration determinants in microglia and CD4+ T cells, highlighting the importance of host genome organization in HIV-1 infection.

Virtual ChIP-seq: predicting transcription factor binding by learning from the transcriptome

Existing methods for computational prediction of transcription factor (TF) binding sites evaluate genomic regions with similarity to known TF sequence preferences. Most TF binding sites, however, do not resemble known TF sequence motifs, and many TFs are not sequence-specific. We developed Virtual ChIP-seq, which predicts binding of individual TFs in new cell types, integrating learned associations with gene expression and binding, TF binding sites from other cell types, and chromatin accessibility data in the new cell type. This approach outperforms methods that predict TF binding solely based on sequence preference, predicting binding for 36 TFs (MCC>0.3).

Polycomb contraction differentially regulates terminal human hematopoietic differentiation programs

Background

Lifelong production of the many types of mature blood cells from less differentiated progenitors is a hierarchically ordered process that spans multiple cell divisions. The nature and timing of the molecular events required to integrate the environmental signals, transcription factor activity, epigenetic modifications, and changes in gene expression involved are thus complex and still poorly understood. To address this gap, we generated comprehensive reference epigenomes of 8 phenotypically defined subsets of normal human cord blood.

Results

We describe a striking contraction of H3K27me3 density in differentiated myelo-erythroid cells that resembles a punctate pattern previously ascribed to pluripotent embryonic stem cells. Phenotypically distinct progenitor cell types display a nearly identical repressive H3K27me3 signature characterized by large organized chromatin K27-modification domains that are retained by mature lymphoid cells but lost in terminally differentiated monocytes and erythroblasts. We demonstrate that inhibition of polycomb group members predicted to control large organized chromatin K27-modification domains influences lymphoid and myeloid fate decisions of primary neonatal hematopoietic progenitors in vitro. We further show that a majority of active enhancers appear in early progenitors, a subset of which are DNA hypermethylated and become hypomethylated and induced during terminal differentiation.

Conclusion

Primitive human hematopoietic cells display a unique repressive H3K27me3 signature that is retained by mature lymphoid cells but is lost in monocytes and erythroblasts. Intervention data implicate that control of this chromatin state change is a requisite part of the process whereby normal human hematopoietic progenitor cells make lymphoid and myeloid fate decisions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01315-1.

Notch-dependent DNA cis-regulatory elements and their dose-dependent control of C. elegans stem cell self-renewal

A long-standing biological question is how DNA cis-regulatory elements shape transcriptional patterns during metazoan development. Reporter constructs, cell culture assays and computational modeling have made major contributions to answering this question, but analysis of elements in their natural context is an important complement. Here, we mutate Notch-dependent LAG-1 binding sites (LBSs) in the endogenous Caenorhabditis elegans sygl-1 gene, which encodes a key stem cell regulator, and analyze the consequences on sygl-1 expression (nascent transcripts, mRNA, protein) and stem cell maintenance. Mutation of one LBS in a three-element cluster approximately halved both expression and stem cell pool size, whereas mutation of two LBSs essentially abolished them. Heterozygous LBS mutant clusters provided intermediate values. Our results lead to two major conclusions. First, both LBS number and configuration impact cluster activity: LBSs act additively in trans and synergistically in cis. Second, the SYGL-1 gradient promotes self-renewal above its functional threshold and triggers differentiation below the threshold. Our approach of coupling CRISPR/Cas9 LBS mutations with effects on both molecular and biological readouts establishes a powerful model for in vivo analyses of DNA cis-regulatory elements.

Chromatin domain alterations linked to 3D genome organization in a large cohort of schizophrenia and bipolar disorder brains

Chromosomal organization, scaling from the 147-base pair (bp) nucleosome to megabase-ranging domains encompassing multiple transcriptional units, including heritability loci for psychiatric traits, remains largely unexplored in the human brain. In this study, we constructed promoter- and enhancer-enriched nucleosomal histone modification landscapes for adult prefrontal cortex from H3-lysine 27 acetylation and H3-lysine 4 trimethylation profiles, generated from 388 controls and 351 individuals diagnosed with schizophrenia (SCZ) or bipolar disorder (BD) (n = 739). We mapped thousands of cis-regulatory domains (CRDs), revealing fine-grained, 104-106-bp chromosomal organization, firmly integrated into Hi-C topologically associating domain stratification by open/repressive chromosomal environments and nuclear topography. Large clusters of hyper-acetylated CRDs were enriched for SCZ heritability, with prominent representation of regulatory sequences governing fetal development and glutamatergic neuron signaling. Therefore, SCZ and BD brains show coordinated dysregulation of risk-associated regulatory sequences assembled into kilobase- to megabase-scaling chromosomal domains.

A primary hierarchically organized patient-derived model enables in depth interrogation of stemness driven by the coding and non-coding genome

Many cancers are organized as cellular hierarchies sustained by cancer stem cells (CSC), whose eradication is crucial for achieving long-term remission. Difficulties to isolate and undertake in vitro and in vivo experimental studies of rare CSC under conditions that preserve their original properties currently constitute a bottleneck for identifying molecular mechanisms involving coding and non-coding genomic regions that govern stemness. We focussed on acute myeloid leukemia (AML) as a paradigm of the CSC model and developed a patient-derived system termed OCI-AML22 that recapitulates the cellular hierarchy driven by leukemia stem cells (LSC). Through classical flow sorting and functional analyses, we established that a single phenotypic population is highly enriched for LSC. The LSC fraction can be easily isolated and serially expanded in culture or in xenografts while faithfully recapitulating functional, transcriptional and epigenetic features of primary LSCs. A novel non-coding regulatory element was identified with a new computational approach using functionally validated primary AML LSC fractions and its role in LSC stemness validated through efficient CRISPR editing using methods optimized for OCI-AML22 LSC. Collectively, OCI-AML22 constitutes a valuable resource to uncover mechanisms governing CSC driven malignancies.

UniBind: maps of high-confidence direct TF-DNA interactions across nine species

BACKGROUND

Transcription factors (TFs) bind specifically to TF binding sites (TFBSs) at cis-regulatory regions to control transcription. It is critical to locate these TF-DNA interactions to understand transcriptional regulation. Efforts to predict bona fide TFBSs benefit from the availability of experimental data mapping DNA binding regions of TFs (chromatin immunoprecipitation followed by sequencing - ChIP-seq).

RESULTS

In this study, we processed ~ 10,000 public ChIP-seq datasets from nine species to provide high-quality TFBS predictions. After quality control, it culminated with the prediction of ~ 56 million TFBSs with experimental and computational support for direct TF-DNA interactions for 644 TFs in > 1000 cell lines and tissues. These TFBSs were used to predict > 197,000 cis-regulatory modules representing clusters of binding events in the corresponding genomes. The high-quality of the TFBSs was reinforced by their evolutionary conservation, enrichment at active cis-regulatory regions, and capacity to predict combinatorial binding of TFs. Further, we confirmed that the cell type and tissue specificity of enhancer activity was correlated with the number of TFs with binding sites predicted in these regions. All the data is provided to the community through the UniBind database that can be accessed through its web-interface ( https://unibind.uio.no/ ), a dedicated RESTful API, and as genomic tracks. Finally, we provide an enrichment tool, available as a web-service and an R package, for users to find TFs with enriched TFBSs in a set of provided genomic regions.

CONCLUSIONS

UniBind is the first resource of its kind, providing the largest collection of high-confidence direct TF-DNA interactions in nine species.

Biological and therapeutic implications of a unique subtype of NPM1 mutated AML

In acute myeloid leukemia (AML), molecular heterogeneity across patients constitutes a major challenge for prognosis and therapy. AML with NPM1 mutation is a distinct genetic entity in the revised World Health Organization classification. However, differing patterns of co-mutation and response to therapy within this group necessitate further stratification. Here we report two distinct subtypes within NPM1 mutated AML patients, which we label as primitive and committed based on the respective presence or absence of a stem cell signature. Using gene expression (RNA-seq), epigenomic (ATAC-seq) and immunophenotyping (CyToF) analysis, we associate each subtype with specific molecular characteristics, disease differentiation state and patient survival. Using ex vivo drug sensitivity profiling, we show a differential drug response of the subtypes to specific kinase inhibitors, irrespective of the FLT3-ITD status. Differential drug responses of the primitive and committed subtype are validated in an independent AML cohort. Our results highlight heterogeneity among NPM1 mutated AML patient samples based on stemness and suggest that the addition of kinase inhibitors to the treatment of cases with the primitive signature, lacking FLT3-ITD, could have therapeutic benefit.

Super enhancers-Functional cores under the 3D genome

Complex biochemical reactions take place in the nucleus all the time. Transcription machines must follow the rules. The chromatin state, especially the three-dimensional structure of the genome, plays an important role in gene regulation and expression. The super enhancers are important for defining cell identity in mammalian developmental processes and human diseases. It has been shown that the major components of transcriptional activation complexes are recruited by super enhancer to form phase-separated condensates. We summarize the current knowledge about super enhancer in the 3D genome. Furthermore, a new related transcriptional regulation model from super enhancer is outlined to explain its role in the mammalian cell progress.

Large organized chromatin lysine domains help distinguish primitive from differentiated cell populations

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

TADs and Their Borders: Free Movement or Building a Wall?

The tridimensional (3D) organization of mammalian genomes combines structures from different length scales. Within this organization, Topologically Associating Domains (TADs) are visible in Hi-C heat maps at the sub-megabase scale. The integrity of TADs is important for correct gene expression, but in a context-dependent and variable manner. The correct structure and function of TADs require the binding of the CTCF protein at both borders, which appears to block an active and dynamic mechanism of "Cohesin-mediated loop extrusion." As a result, mammalian TADs appear as so-called "loop domains" in Hi-C data, which are the focus of this review. Here, we present a reanalysis of TADs from three "golden-standard" mammalian Hi-C data sets. Despite the prominent presence of TADs in Hi-C heat maps from all studies, we find consistently that regions within these domains are only moderately insulated from their surroundings. Moreover, single-cell Hi-C and superresolution microscopy have revealed that the structure of TADs and the position of their borders can vary from cell to cell. The function of TADs as units of gene regulation may thus require additional aspects, potentially incorporating the mechanism of loop extrusion as well. Recent developments in single-cell and multi-contact genomics and superresolution microscopy assays will be instrumental to link TAD formation and structure to their function in transcriptional regulation.