An integrative epigenome-based strategy for unbiased functional profiling of clinical kinase inhibitors.

More than 500 kinases are implicated in the control of most cellular process in mammals, and deregulation of their activity is linked to cancer and inflammatory disorders. 80 clinical kinase inhibitors (CKIs) have been approved for clinical use and hundreds are in various stages of development. However, CKIs inhibit other kinases in addition to the intended target(s), causing both enhanced clinical effects and undesired side effects that are only partially predictable based on in vitro selectivity profiling. Here, we report an integrative approach grounded on the use of chromatin modifications as unbiased, information-rich readouts of the functional effects of CKIs on macrophage activation. This approach exceeded the performance of transcriptome-based approaches and allowed us to identify similarities and differences among CKIs with identical intended targets, to recognize novel CKI specificities and to pinpoint CKIs that may be repurposed to control inflammation, thus supporting the utility of this strategy to improve selection and use of CKIs in clinical settings.

Intra- and interchromosomal contact mapping reveals the Igh locus has extensive conformational heterogeneity and interacts with B-lineage genes.

To produce a diverse antibody repertoire, immunoglobulin heavy-chain (Igh) loci undergo large-scale alterations in structure to facilitate juxtaposition and recombination of spatially separated variable (VH), diversity (DH), and joining (JH) genes. These chromosomal alterations are poorly understood. Uncovering their patterns shows how chromosome dynamics underpins antibody diversity. Using tiled Capture Hi-C, we produce a comprehensive map of chromatin interactions throughout the 2.8-Mb Igh locus in progenitor B cells. We find that the Igh locus folds into semi-rigid subdomains and undergoes flexible looping of the VH genes to its 3' end, reconciling two views of locus organization. Deconvolution of single Igh locus conformations using polymer simulations identifies thousands of different structures. This heterogeneity may underpin the diversity of V(D)J recombination events. All three immunoglobulin loci also participate in a highly specific, developmentally regulated network of interchromosomal interactions with genes encoding B cell-lineage factors. This suggests a model of interchromosomal coordination of B cell development.

Locality of contacts determines the subdiffusion exponents in polymeric models of chromatin.

Loop extrusion by motor proteins mediates the attractive interactions in chromatin on the length scale of megabases, providing the polymer with a well-defined structure and at the same time determining its dynamics. The mean-square displacement of chromatin loci varies from a Rouse-like scaling to a more constrained subdiffusion, depending on cell type, genomic region, and time scale. With a simple polymeric model, we show that such a Rouse-like dynamics occurs when the parameters of the model are chosen so that contacts are local along the chain, while in the presence of nonlocal contacts we observe subdiffusion at short time scales with exponents smaller than 0.5. Such exponents are independent of the detailed choice of the parameters and build a master curve that depends only on the mean locality of the resulting contacts. We compare the loop-extrusion model with a polymeric model with static links, showing that also in this case only the presence of nonlocal contacts can produce low-exponent subdiffusion. We interpret these results in terms of a simple analytical model.

Analyzing the Genome-Wide Distribution of Histone Marks by CUT&Tag in Drosophila Embryos.

CUT&Tag is a method to map the genome-wide distribution of histone modifications and some chromatin-associated proteins. CUT&Tag relies on antibody-targeted chromatin tagmentation and can easily be scaled up or automatized. This protocol provides clear experimental guidelines and helpful considerations when planning and executing CUT&Tag experiments.

Chromosome Conformation Capture Followed by Genome-Wide Sequencing (Hi-C) in Drosophila Embryos.

This protocol provides specific details on how to perform Hi-C, the genome-wide version of Chromosome Conformation Capture (3C) followed by high-throughput sequencing, in Drosophila embryos. Hi-C provides a genome-wide population-averaged snapshot of the 3D genome organization within nuclei. In Hi-C, formaldehyde-cross-linked chromatin is enzymatically digested using restriction enzymes; digested fragments are biotinylated and subjected to proximity ligation; ligated fragments are purified using streptavidin followed by paired-end sequencing. Hi-C allows the detection of higher order folding structures such as topologically associated domains (TADs) and active/inactive compartments (A/B compartments, respectively). Performing this assay in developing embryos gives the unique opportunity to investigate dynamic chromatin changes when 3D chromatin structure is established in embryogenesis.

Cohesin and CTCF control the dynamics of chromosome folding.

In mammals, interactions between sequences within topologically associating domains enable control of gene expression across large genomic distances. Yet it is unknown how frequently such contacts occur, how long they last and how they depend on the dynamics of chromosome folding and loop extrusion activity of cohesin. By imaging chromosomal locations at high spatial and temporal resolution in living cells, we show that interactions within topologically associating domains are transient and occur frequently during the course of a cell cycle. Interactions become more frequent and longer in the presence of convergent CTCF sites, resulting in suppression of variability in chromosome folding across time. Supported by physical models of chromosome dynamics, our data suggest that CTCF-anchored loops last around 10 min. Our results show that long-range transcriptional regulation might rely on transient physical proximity, and that cohesin and CTCF stabilize highly dynamic chromosome structures, facilitating selected subsets of chromosomal interactions.

Inversion of a topological domain leads to restricted changes in its gene expression and affects interdomain communication.

The interplay between the topological organization of the genome and the regulation of gene expression remains unclear. Depletion of molecular factors (e.g. CTCF) underlying topologically associating domains (TADs) leads to modest alterations in gene expression, whereas genomic rearrangements involving TAD boundaries disrupt normal gene expression and can lead to pathological phenotypes. Here, we targeted the TAD neighboring that of the noncoding transcript Xist, which controls X-chromosome inactivation. Inverting 245 kb within the TAD led to expected rearrangement of CTCF-based contacts but revealed heterogeneity in the 'contact' potential of different CTCF sites. Expression of most genes therein remained unaffected in mouse embryonic stem cells and during differentiation. Interestingly, expression of Xist was ectopically upregulated. The same inversion in mouse embryos led to biased Xist expression. Smaller inversions and deletions of CTCF clusters led to similar results: rearrangement of contacts and limited changes in local gene expression, but significant changes in Xist expression in embryos. Our study suggests that the wiring of regulatory interactions within a TAD can influence the expression of genes in neighboring TADs, highlighting the existence of mechanisms of inter-TAD communication.

Nonlinear control of transcription through enhancer-promoter interactions.

Chromosome structure in mammals is thought to regulate transcription by modulating three-dimensional interactions between enhancers and promoters, notably through CTCF-mediated loops and topologically associating domains (TADs)1-4. However, how chromosome interactions are actually translated into transcriptional outputs remains unclear. Here, to address this question, we use an assay to position an enhancer at large numbers of densely spaced chromosomal locations relative to a fixed promoter, and measure promoter output and interactions within a genomic region with minimal regulatory and structural complexity. A quantitative analysis of hundreds of cell lines reveals that the transcriptional effect of an enhancer depends on its contact probabilities with the promoter through a nonlinear relationship. Mathematical modelling suggests that nonlinearity might arise from transient enhancer-promoter interactions being translated into slower promoter bursting dynamics in individual cells, therefore uncoupling the temporal dynamics of interactions from those of transcription. This uncovers a potential mechanism of how distal enhancers act from large genomic distances, and of how topologically associating domain boundaries block distal enhancers. Finally, we show that enhancer strength also determines absolute transcription levels as well as the sensitivity of a promoter to CTCF-mediated transcriptional insulation. Our measurements establish general principles for the context-dependent role of chromosome structure in long-range transcriptional regulation.

Polymer Folding Simulations from Hi-C Data.

In the absence of a clear molecular understanding of the mechanism that stabilizes specific contacts in interphasic chromatin, we resort to the principle of maximum entropy to build a polymeric model based on the Hi-C data of the specific system one wants to study. The interactions are set by an iterative Monte Carlo algorithm to reproduce the average contacts summarized by the Hi-C map. The study of the ensemble of conformations generated by the algorithm can report a much richer set of information than the experimental map alone, including colocalization of multiple sites, fluctuations of the contacts, and kinetical properties.

deepBlink: threshold-independent detection and localization of diffraction-limited spots.

Detection of diffraction-limited spots in single-molecule microscopy images is traditionally performed with mathematical operators designed for idealized spots. This process requires manual tuning of parameters that is time-consuming and not always reliable. We have developed deepBlink, a neural network-based method to detect and localize spots automatically. We demonstrate that deepBlink outperforms other state-of-the-art methods across six publicly available datasets containing synthetic and experimental data.

HP1 drives de novo 3D genome reorganization in early Drosophila embryos.

Fundamental features of 3D genome organization are established de novo in the early embryo, including clustering of pericentromeric regions, the folding of chromosome arms and the segregation of chromosomes into active (A-) and inactive (B-) compartments. However, the molecular mechanisms that drive de novo organization remain unknown1,2. Here, by combining chromosome conformation capture (Hi-C), chromatin immunoprecipitation with high-throughput sequencing (ChIP-seq), 3D DNA fluorescence in situ hybridization (3D DNA FISH) and polymer simulations, we show that heterochromatin protein 1a (HP1a) is essential for de novo 3D genome organization during Drosophila early development. The binding of HP1a at pericentromeric heterochromatin is required to establish clustering of pericentromeric regions. Moreover, HP1a binding within chromosome arms is responsible for overall chromosome folding and has an important role in the formation of B-compartment regions. However, depletion of HP1a does not affect the A-compartment, which suggests that a different molecular mechanism segregates active chromosome regions. Our work identifies HP1a as an epigenetic regulator that is involved in establishing the global structure of the genome in the early embryo.

Effective model of loop extrusion predicts chromosomal domains.

An active loop-extrusion mechanism is regarded as the main out-of-equilibrium mechanism responsible for the structuring of megabase-sized domains in chromosomes. We developed a model to study the dynamics of the chromosome fiber by solving the kinetic equations associated with the motion of the extruder. By averaging out the position of the extruder along the chain, we build an effective equilibrium model capable of reproducing experimental contact maps based solely on the positions of extrusion-blocking proteins. We assessed the quality of the effective model using numerical simulations of chromosomal segments and comparing the results with explicit-extruder models and experimental data.

A Conserved Noncoding Locus Regulates Random Monoallelic Xist Expression across a Topological Boundary.

cis-Regulatory communication is crucial in mammalian development and is thought to be restricted by the spatial partitioning of the genome in topologically associating domains (TADs). Here, we discovered that the Xist locus is regulated by sequences in the neighboring TAD. In particular, the promoter of the noncoding RNA Linx (LinxP) acts as a long-range silencer and influences the choice of X chromosome to be inactivated. This is independent of Linx transcription and independent of any effect on Tsix, the antisense regulator of Xist that shares the same TAD as Linx. Unlike Tsix, LinxP is well conserved across mammals, suggesting an ancestral mechanism for random monoallelic Xist regulation. When introduced in the same TAD as Xist, LinxP switches from a silencer to an enhancer. Our study uncovers an unsuspected regulatory axis for X chromosome inactivation and a class of cis-regulatory effects that may exploit TAD partitioning to modulate developmental decisions.

DamC reveals principles of chromatin folding in vivo without crosslinking and ligation.

Current understanding of chromosome folding is largely reliant on chromosome conformation capture (3C)-based experiments, where chromosomal interactions are detected as ligation products after chromatin crosslinking. To measure chromosome structure in vivo, quantitatively and without crosslinking and ligation, we implemented a modified version of DNA adenine methyltransferase identification (DamID) named DamC, which combines DNA methylation-based detection of chromosomal interactions with next-generation sequencing and biophysical modeling of methylation kinetics. DamC performed in mouse embryonic stem cells provides the first in vivo validation of the existence of topologically associating domains (TADs), CTCF loops and confirms 3C-based measurements of the scaling of contact probabilities. Combining DamC with transposon-mediated genomic engineering shows that new loops can be formed between ectopic and endogenous CTCF sites, which redistributes physical interactions within TADs. DamC provides the first crosslinking- and ligation-free demonstration of the existence of key structural features of chromosomes and provides novel insights into how chromosome structure within TADs can be manipulated.

The bipartite TAD organization of the X-inactivation center ensures opposing developmental regulation of Tsix and Xist.

The mouse X-inactivation center (Xic) locus represents a powerful model for understanding the links between genome architecture and gene regulation, with the non-coding genes Xist and Tsix showing opposite developmental expression patterns while being organized as an overlapping sense/antisense unit. The Xic is organized into two topologically associating domains (TADs) but the role of this architecture in orchestrating cis-regulatory information remains elusive. To explore this, we generated genomic inversions that swap the Xist/Tsix transcriptional unit and place their promoters in each other's TAD. We found that this led to a switch in their expression dynamics: Xist became precociously and ectopically upregulated, both in male and female pluripotent cells, while Tsix expression aberrantly persisted during differentiation. The topological partitioning of the Xic is thus critical to ensure proper developmental timing of X inactivation. Our study illustrates how the genomic architecture of cis-regulatory landscapes can affect the regulation of mammalian developmental processes.

Single-Molecule Imaging of mRNA Localization and Regulation during the Integrated Stress Response.

Biological phase transitions form membrane-less organelles that generate distinct cellular environments. How molecules are partitioned between these compartments and the surrounding cellular space and the functional consequence of this localization is not well understood. Here, we report the localization of mRNA to stress granules (SGs) and processing bodies (PBs) and its effect on translation and degradation during the integrated stress response. Using single mRNA imaging in living human cells, we find that the interactions of mRNAs with SGs and PBs have different dynamics, very few mRNAs directly move between SGs and PBs, and that specific RNA-binding proteins can anchor mRNAs within these compartments. During recovery from stress, we show that mRNAs that were within SGs and PBs are translated and degraded at similar rates as their cytosolic counterparts. Our work provides a framework for using single-molecule measurements to directly investigate the molecular mechanisms of phase-separated compartments within their cellular environment.

The Dynamics of mRNA Turnover Revealed by Single-Molecule Imaging in Single Cells.

RNA degradation plays a fundamental role in regulating gene expression. In order to characterize the spatiotemporal dynamics of RNA turnover in single cells, we developed a fluorescent biosensor based on dual-color, single-molecule RNA imaging that allows intact transcripts to be distinguished from stabilized degradation intermediates. Using this method, we measured mRNA decay in single cells and found that individual degradation events occur independently within the cytosol and are not enriched within processing bodies. We show that slicing of an mRNA targeted for endonucleolytic cleavage by the RNA-induced silencing complex can be observed in real time in living cells. This methodology provides a framework for investigating the entire life history of individual mRNAs from birth to death in single cells.

Reciprocal insulation analysis of Hi-C data shows that TADs represent a functionally but not structurally privileged scale in the hierarchical folding of chromosomes.

Understanding how regulatory sequences interact in the context of chromosomal architecture is a central challenge in biology. Chromosome conformation capture revealed that mammalian chromosomes possess a rich hierarchy of structural layers, from multi-megabase compartments to sub-megabase topologically associating domains (TADs) and sub-TAD contact domains. TADs appear to act as regulatory microenvironments by constraining and segregating regulatory interactions across discrete chromosomal regions. However, it is unclear whether other (or all) folding layers share similar properties, or rather TADs constitute a privileged folding scale with maximal impact on the organization of regulatory interactions. Here, we present a novel algorithm named CaTCH that identifies hierarchical trees of chromosomal domains in Hi-C maps, stratified through their reciprocal physical insulation, which is a single and biologically relevant parameter. By applying CaTCH to published Hi-C data sets, we show that previously reported folding layers appear at different insulation levels. We demonstrate that although no structurally privileged folding level exists, TADs emerge as a functionally privileged scale defined by maximal boundary enrichment in CTCF and maximal cell-type conservation. By measuring transcriptional output in embryonic stem cells and neural precursor cells, we show that the likelihood that genes in a domain are coregulated during differentiation is also maximized at the scale of TADs. Finally, we observe that regulatory sequences occur at genomic locations corresponding to optimized mutual interactions at the same scale. Our analysis suggests that the architectural functionality of TADs arises from the interplay between their ability to partition interactions and the specific genomic position of regulatory sequences.