Beyond A and B Compartments: how major nuclear locales define nuclear genome organization and function.

Models of nuclear genome organization often propose a binary division into active versus inactive compartments, yet they overlook nuclear bodies. Here we integrated analysis of sequencing and image-based data to compare genome organization in four human cell types relative to three different nuclear locales: the nuclear lamina, nuclear speckles, and nucleoli. Whereas gene expression correlates mostly with nuclear speckle proximity, DNA replication timing correlates with proximity to multiple nuclear locales. Speckle attachment regions emerge as DNA replication initiation zones whose replication timing and gene composition vary with their attachment frequency. Most facultative LADs retain a partially repressed state as iLADs, despite their positioning in the nuclear interior. Knock out of two lamina proteins, Lamin A and LBR, causes a shift of H3K9me3-enriched LADs from lamina to nucleolus, and a reciprocal relocation of H3K27me3-enriched partially repressed iLADs from nucleolus to lamina. Thus, these partially repressed iLADs appear to compete with LADs for nuclear lamina attachment with consequences for replication timing. The nuclear organization in adherent cells is polarized with nuclear bodies and genomic regions segregating both radially and relative to the equatorial plane. Together, our results underscore the importance of considering genome organization relative to nuclear locales for a more complete understanding of the spatial and functional organization of the human genome.

Nucleolus and centromere TSA-Seq reveals variable localization of heterochromatin in different cell types.

Genome differential positioning within interphase nuclei remains poorly explored. We extended and validated TSA-seq to map genomic regions near nucleoli and pericentric heterochromatin in four human cell lines. Our study confirmed that smaller chromosomes localize closer to nucleoli but further deconvolved this by revealing a preference for chromosome arms below 36-46 Mbp in length. We identified two lamina associated domain subsets through their differential nuclear lamina versus nucleolar positioning in different cell lines which showed distinctive patterns of DNA replication timing and gene expression across all cell lines. Unexpectedly, active, nuclear speckle-associated genomic regions were found near typically repressive nuclear compartments, which is attributable to the close proximity of nuclear speckles and nucleoli in some cell types, and association of centromeres with nuclear speckles in hESCs. Our study points to a more complex and variable nuclear genome organization than suggested by current models, as revealed by our TSA-seq methodology.

Perturbations in 3D genome organization can promote acquired drug resistance.

Acquired drug resistance is a major problem in the treatment of cancer. hTERT-immortalized, untransformed RPE-1 cells can acquire resistance to Taxol by derepressing the ABCB1 gene, encoding for the multidrug transporter P-gP. Here, we investigate how the ABCB1 gene is derepressed. ABCB1 activation is associated with reduced H3K9 trimethylation, increased H3K27 acetylation, and ABCB1 displacement from the nuclear lamina. While altering DNA methylation and H3K27 methylation had no major impact on ABCB1 expression, nor did it promote resistance, disrupting the nuclear lamina component Lamin B Receptor did promote the acquisition of a Taxol-resistant phenotype in a subset of cells. CRISPRa-mediated gene activation supported the notion that lamina dissociation influences ABCB1 derepression. We propose a model in which nuclear lamina dissociation of a repressed gene allows for its activation, implying that deregulation of the 3D genome topology could play an important role in tumor evolution and the acquisition of drug resistance.

Optimisation of TP53 reporters by systematic dissection of synthetic TP53 response elements.

TP53 is a transcription factor that controls multiple cellular processes, including cell cycle arrest, DNA repair and apoptosis. The relation between TP53 binding site architecture and transcriptional output is still not fully understood. Here, we systematically examined in three different cell lines the effects of binding site affinity and copy number on TP53-dependent transcriptional output, and also probed the impact of spacer length and sequence between adjacent binding sites, and of core promoter identity. Paradoxically, we found that high-affinity TP53 binding sites are less potent than medium-affinity sites. TP53 achieves supra-additive transcriptional activation through optimally spaced adjacent binding sites, suggesting a cooperative mechanism. Optimally spaced adjacent binding sites have a ∼10-bp periodicity, suggesting a role for spatial orientation along the DNA double helix. We leveraged these insights to construct a log-linear model that explains activity from sequence features, and to identify new highly active and sensitive TP53 reporters.

Heritable transcriptional defects from aberrations of nuclear architecture.

Transcriptional heterogeneity due to plasticity of the epigenetic state of chromatin contributes to tumour evolution, metastasis and drug resistance1-3. However, the mechanisms that cause this epigenetic variation are incompletely understood. Here we identify micronuclei and chromosome bridges, aberrations in the nucleus common in cancer4,5, as sources of heritable transcriptional suppression. Using a combination of approaches, including long-term live-cell imaging and same-cell single-cell RNA sequencing (Look-Seq2), we identified reductions in gene expression in chromosomes from micronuclei. With heterogeneous penetrance, these changes in gene expression can be heritable even after the chromosome from the micronucleus has been re-incorporated into a normal daughter cell nucleus. Concomitantly, micronuclear chromosomes acquire aberrant epigenetic chromatin marks. These defects may persist as variably reduced chromatin accessibility and reduced gene expression after clonal expansion from single cells. Persistent transcriptional repression is strongly associated with, and may be explained by, markedly long-lived DNA damage. Epigenetic alterations in transcription may therefore be inherently coupled to chromosomal instability and aberrations in nuclear architecture.

Defining the fine structure of promoter activity on a genome-wide scale with CISSECTOR.

Classic promoter mutagenesis strategies can be used to study how proximal promoter regions regulate the expression of particular genes of interest. This is a laborious process, in which the smallest sub-region of the promoter still capable of recapitulating expression in an ectopic setting is first identified, followed by targeted mutation of putative transcription factor binding sites. Massively parallel reporter assays such as survey of regulatory elements (SuRE) provide an alternative way to study millions of promoter fragments in parallel. Here we show how a generalized linear model (GLM) can be used to transform genome-scale SuRE data into a high-resolution genomic track that quantifies the contribution of local sequence to promoter activity. This coefficient track helps identify regulatory elements and can be used to predict promoter activity of any sub-region in the genome. It thus allows in silico dissection of any promoter in the human genome to be performed. We developed a web application, available at cissector.nki.nl, that lets researchers easily perform this analysis as a starting point for their research into any promoter of interest.

Dynamic chromosomal interactions and control of heterochromatin positioning by Ki-67.

Ki-67 is a chromatin-associated protein with a dynamic distribution pattern throughout the cell cycle and is thought to be involved in chromatin organization. The lack of genomic interaction maps has hampered a detailed understanding of its roles, particularly during interphase. By pA-DamID mapping in human cell lines, we find that Ki-67 associates with large genomic domains that overlap mostly with late-replicating regions. Early in interphase, when Ki-67 is present in pre-nucleolar bodies, it interacts with these domains on all chromosomes. However, later in interphase, when Ki-67 is confined to nucleoli, it shows a striking shift toward small chromosomes. Nucleolar perturbations indicate that these cell cycle dynamics correspond to nucleolar maturation during interphase, and suggest that nucleolar sequestration of Ki-67 limits its interactions with larger chromosomes. Furthermore, we demonstrate that Ki-67 does not detectably control chromatin-chromatin interactions during interphase, but it competes with the nuclear lamina for interaction with late-replicating DNA, and it controls replication timing of (peri)centromeric regions. Together, these results reveal a highly dynamic choreography of genome interactions and roles for Ki-67 in heterochromatin organization.

CTCF and cohesin promote focal detachment of DNA from the nuclear lamina.

BACKGROUND: Lamina-associated domains (LADs) are large genomic regions that are positioned at the nuclear lamina. It has remained largely unclear what drives the positioning and demarcation of LADs. Because the insulator protein CTCF is enriched at LAD borders, it was postulated that CTCF binding could position some LAD boundaries, possibly through its function in stalling cohesin and hence preventing cohesin invading into the LAD. To test this, we mapped genome-nuclear lamina interactions in mouse embryonic stem cells after rapid depletion of CTCF and other perturbations of cohesin dynamics. RESULTS: CTCF and cohesin contribute to a sharp transition in lamina interactions at LAD borders, while LADs are maintained after depletion of these proteins, also at borders marked by CTCF. CTCF and cohesin may thus reinforce LAD borders, but do not position these. CTCF binding sites within LADs are locally detached from the lamina and enriched for accessible DNA and active histone modifications. Remarkably, despite lamina positioning being strongly correlated with genome inactivity, this DNA remains accessible after the local detachment is lost following CTCF depletion. At a chromosomal scale, cohesin depletion and cohesin stabilization by depletion of the unloading factor WAPL quantitatively affect lamina interactions, indicative of perturbed chromosomal positioning in the nucleus. Finally, while H3K27me3 is locally enriched at CTCF-marked LAD borders, we find no evidence for an interplay between CTCF and H3K27me3 on lamina interactions. CONCLUSIONS: These findings illustrate that CTCF and cohesin are not primary determinants of LAD patterns. Rather, these proteins locally modulate NL interactions.

From fluorescent foci to sequence: Illuminating DNA double strand break repair by high-throughput sequencing technologies.

Technologies to study DNA double-strand break (DSB) repair have traditionally mostly relied on fluorescence read-outs, either by microscopy or flow cytometry. The advent of high throughput sequencing (HTS) has created fundamentally new opportunities to study the mechanisms underlying DSB repair. Here, we review the suite of HTS-based assays that are used to study three different aspects of DNA repair: detection of broken ends, protein recruitment and pathway usage. We highlight new opportunities that HTS technology offers towards a better understanding of the DSB repair process.

Systematic analysis of intrinsic enhancer-promoter compatibility in the mouse genome.

Gene expression is in part controlled by cis-regulatory elements (CREs) such as enhancers and repressive elements. Anecdotal evidence has indicated that a CRE and a promoter need to be biochemically compatible for promoter regulation to occur, but this compatibility has remained poorly characterized in mammalian cells. We used high-throughput combinatorial reporter assays to test thousands of CRE-promoter pairs from three Mb-sized genomic regions in mouse cells. This revealed that CREs vary substantially in their promoter compatibility, ranging from striking specificity to broad promiscuity. More than half of the tested CREs exhibit significant promoter selectivity. Housekeeping promoters tend to have similar CRE preferences, but other promoters exhibit a wide diversity of compatibilities. Higher-order transcription factors (TF) motif combinations may account for compatibility. CRE-promoter selectivity does not correlate with looping interactions in the native genomic context, suggesting that chromatin folding and compatibility are two orthogonal mechanisms that confer specificity to gene regulation.

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.

Lamina-associated domains: Tethers and looseners.

Lamina-associated domains (LADs) are large heterochromatic regions that are positioned at the nuclear lamina (NL). A major question is how LAD-NL interactions are mediated and controlled. Here, we review recent progress in the search for molecular tethers and looseners of LADs and we discuss the link between LAD-NL tethering, transcription regulation, and genome replication. We also provide a brief summary of technological advances that may uncover new aspects of LAD biology.

Genome-Wide Mapping and Microscopy Visualization of Protein-DNA Interactions by pA-DamID.

Several methods have been developed to map protein-DNA interactions genome-wide in the last decades. Protein A-DamID (pA-DamID) is a recent addition to this list with distinct advantages. pA-DamID relies on antibody-based targeting of the bacterial Dam enzyme, resulting in adenine methylation of DNA in contact with the protein of interest. This m6A can then be visualized by microscopy, or mapped genome-wide. The main advantages of pA-DamID are an easy and direct visualization of DNA that is in contact with the protein of interest, unbiased mapping of protein-DNA interactions, and the possibility to select specific subpopulations of cells by flow cytometry before further sample processing. pA-DamID is particularly suited to study proteins that form large chromatin domains or that are part of distinct nuclear structures such as the nuclear lamina. This chapter describes the pA-DamID procedure from cell harvesting to the preparation of microscopy slides and high-throughput sequencing libraries.

Impact of chromatin context on Cas9-induced DNA double-strand break repair pathway balance.

DNA double-strand break (DSB) repair is mediated by multiple pathways. It is thought that the local chromatin context affects the pathway choice, but the underlying principles are poorly understood. Using a multiplexed reporter assay in combination with Cas9 cutting, we systematically measure the relative activities of three DSB repair pathways as a function of chromatin context in >1,000 genomic locations. This reveals that non-homologous end-joining (NHEJ) is broadly biased toward euchromatin, while the contribution of microhomology-mediated end-joining (MMEJ) is higher in specific heterochromatin contexts. In H3K27me3-marked heterochromatin, inhibition of the H3K27 methyltransferase EZH2 reverts the balance toward NHEJ. Single-stranded template repair (SSTR), often used for precise CRISPR editing, competes with MMEJ and is moderately linked to chromatin context. These results provide insight into the impact of chromatin on DSB repair pathway balance and guidance for the design of Cas9-mediated genome editing experiments.

m6A methylation potentiates cytosolic dsDNA recognition in a sequence-specific manner.

Nucleic acid sensing through pattern recognition receptors is critical for immune recognition of microbial infections. Microbial DNA is frequently methylated at the N6 position of adenines (m6A), a modification that is rare in mammalian host DNA. We show here how that m6A methylation of 5'-GATC-3' motifs augments the immunogenicity of synthetic double-stranded (ds)DNA in murine macrophages and dendritic cells. Transfection with m6A-methylated DNA increased the expression of the activation markers CD69 and CD86, and of Ifnß, iNos and Cxcl10 mRNA. Similar to unmethylated cytosolic dsDNA, recognition of m6A DNA occurs independently of TLR and RIG-I signalling, but requires the two key mediators of cytosolic DNA sensing, STING and cGAS. Intriguingly, the response to m6A DNA is sequence-specific. m6A is immunostimulatory in some motifs, but immunosuppressive in others, a feature that is conserved between mouse and human macrophages. In conclusion, epigenetic alterations of DNA depend on the context of the sequence and are differentially perceived by innate cells, a feature that could potentially be used for the design of immune-modulating therapeutics.

SPIN reveals genome-wide landscape of nuclear compartmentalization.

We report SPIN, an integrative computational method to reveal genome-wide intranuclear chromosome positioning and nuclear compartmentalization relative to multiple nuclear structures, which are pivotal for modulating genome function. As a proof-of-principle, we use SPIN to integrate nuclear compartment mapping (TSA-seq and DamID) and chromatin interaction data (Hi-C) from K562 cells to identify 10 spatial compartmentalization states genome-wide relative to nuclear speckles, lamina, and putative associations with nucleoli. These SPIN states show novel patterns of genome spatial organization and their relation to other 3D genome features and genome function (transcription and replication timing). SPIN provides critical insights into nuclear spatial and functional compartmentalization.

Protocol: A Multiplexed Reporter Assay to Study Effects of Chromatin Context on DNA Double-Strand Break Repair.

DNA double-strand breaks (DSBs) can be repaired through various pathways. Understanding how these pathways are regulated is of great interest for cancer research and optimization of gene editing. The local chromatin environment can affect the balance between repair pathways, but this is still poorly understood. Here we provide a detailed protocol for DSB-TRIP, a technique that utilizes the specific DNA scars left by DSB repair pathways to study pathway usage throughout the genome. DSB-TRIP randomly integrates a repair reporter into many genomic locations, followed by the induction of DSBs in the reporter. Multiplexed sequencing of the resulting scars at all integration sites then reveals the balance between several repair pathways, which can be linked to the local chromatin state of the integration sites. Here we present a step-by-step protocol to perform DSB-TRIP in K562 cells and to analyse the data by a dedicated computational pipeline. We discuss strengths and limitations of the technique, as well as potential additional applications to study DNA repair.

Cell cycle dynamics of lamina-associated DNA.

In mammalian interphase nuclei, more than one thousand large genomic regions are positioned at the nuclear lamina (NL). These lamina-associated domains (LADs) are involved in gene regulation and may provide a backbone for the folding of interphase chromosomes. Little is known about the dynamics of LADs during interphase, in particular at the onset of G1 phase and during DNA replication. We developed an antibody-based variant of the DamID technology (named pA-DamID) that allows us to map and visualize genome-NL interactions with high temporal resolution. Application of pA-DamID combined with synchronization and cell sorting experiments reveals that LAD-NL contacts are generally rapidly established early in G1 phase. However, LADs on the distal ~25 Mb of most chromosomes tend to contact the NL first and then gradually detach, while centromere-proximal LADs accumulate gradually at the NL. Furthermore, our data indicate that S-phase chromatin shows transiently increased lamin interactions. These findings highlight a dynamic choreography of LAD-NL contacts during interphase progression and illustrate the usefulness of pA-DamID to study the dynamics of genome compartmentalization.

Phosphorylated Lamins in Euchromatin: New Clues to Progeria.

Lamin proteins not only form the nuclear lamina, but some are also found in the nuclear interior. In this issue of Developmental Cell, Ikegami et al. describe that phosphorylated Lamin C in the nuclear interior interacts with enhancer-like elements and link this to deregulated transcription in progeria.

Local rewiring of genome-nuclear lamina interactions by transcription.

Transcriptionally inactive genes are often positioned at the nuclear lamina (NL), as part of large lamina-associated domains (LADs). Activation of such genes is often accompanied by repositioning toward the nuclear interior. How this process works and how it impacts flanking chromosomal regions are poorly understood. We addressed these questions by systematic activation or inactivation of individual genes, followed by detailed genome-wide analysis of NL interactions, replication timing, and transcription patterns. Gene activation inside LADs typically causes NL detachment of the entire transcription unit, but rarely more than 50-100 kb of flanking DNA, even when multiple neighboring genes are activated. The degree of detachment depends on the expression level and the length of the activated gene. Loss of NL interactions coincides with a switch from late to early replication timing, but the latter can involve longer stretches of DNA. Inactivation of active genes can lead to increased NL contacts. These extensive datasets are a resource for the analysis of LAD rewiring by transcription and reveal a remarkable flexibility of interphase chromosomes.