CTCF-dependent chromatin boundaries formed by asymmetric nucleosome arrays with decreased linker length

The importance of DNA sequence for nucleosome positioning in transcriptional regulation

Nucleosome positioning is a key factor for transcriptional regulation. Nucleosomes regulate the dynamic accessibility of chromatin and interact with the transcription machinery at every stage. Influences to steer nucleosome positioning are diverse, and the according importance of the DNA sequence in contrast to active chromatin remodeling has been the subject of long discussion. In this study, we evaluate the functional role of DNA sequence for all major elements along the process of transcription. We developed a random forest classifier based on local DNA structure that assesses the sequence-intrinsic support for nucleosome positioning. On this basis, we created a simple data resource that we applied genome-wide to the human genome. In our comprehensive analysis, we found a special role of DNA in mediating the competition of nucleosomes with cis-regulatory elements, in enabling steady transcription, for positioning of stable nucleosomes in exons, and for repelling nucleosomes during transcription termination. In contrast, we relate these findings to concurrent processes that generate strongly positioned nucleosomes in vivo that are not mediated by sequence, such as energy-dependent remodeling of chromatin.

The relationship between nucleosome positioning and higher-order genome folding

Eukaryotic genomes are organized into 3D structures, which range from small-scale nucleosome arrays to large-scale chromatin domains. These structures have an important role in the regulation of transcription and other nuclear processes. Despite advances in our understanding of the properties, functions, and underlying mechanisms of genome structures, there are many open questions about the interplay between these structures across scales. In particular, it is not well understood if and how 1D features of nucleosome arrays influence large-scale 3D genome folding patterns. In this review, we discuss recent studies that address these questions and summarize our current understanding of the relationship between nucleosome positioning and higher-order genome folding.

guidedNOMe-seq quantifies chromatin states at single allele resolution for hundreds of custom regions in parallel

Since the introduction of next generation sequencing technologies, the field of epigenomics has evolved rapidly. However, most commonly used assays are enrichment-based methods and thus only semi-quantitative. Nucleosome occupancy and methylome sequencing (NOMe-seq) allows for quantitative inference of chromatin states with single locus resolution, but this requires high sequencing depth and is therefore prohibitively expensive to routinely apply to organisms with large genomes. To overcome this limitation, we introduce guidedNOMe-seq, where we combine NOMe profiling with large scale sgRNA synthesis and Cas9-mediated region-of-interest (ROI) liberation. To facilitate quantitative comparisons between multiple samples, we additionally develop an R package to standardize differential analysis of any type of NOMe-seq data. We extensively benchmark guidedNOMe-seq in a proof-of-concept study, dissecting the interplay of ChAHP and CTCF on chromatin. In summary we present a cost-effective, scalable, and customizable target enrichment extension to the existing NOMe-seq protocol allowing genome-scale quantification of nucleosome occupancy and transcription factor binding at single allele resolution.

DNA methylation patterns of transcription factor binding regions characterize their functional and evolutionary contexts

BACKGROUND

DNA methylation is an important epigenetic modification which has numerous roles in modulating genome function. Its levels are spatially correlated across the genome, typically high in repressed regions but low in transcription factor (TF) binding sites and active regulatory regions. However, the mechanisms establishing genome-wide and TF binding site methylation patterns are still unclear.

RESULTS

Here we use a comparative approach to investigate the association of DNA methylation to TF binding evolution in mammals. Specifically, we experimentally profile DNA methylation and combine this with published occupancy profiles of five distinct TFs (CTCF, CEBPA, HNF4A, ONECUT1, FOXA1) in the liver of five mammalian species (human, macaque, mouse, rat, dog). TF binding sites are lowly methylated, but they often also have intermediate methylation levels. Furthermore, biding sites are influenced by the methylation status of CpGs in their wider binding regions even when CpGs are absent from the core binding motif. Employing a classification and clustering approach, we extract distinct and species-conserved patterns of DNA methylation levels at TF binding regions. CEBPA, HNF4A, ONECUT1, and FOXA1 share the same methylation patterns, while CTCF's differ. These patterns characterize alternative functions and chromatin landscapes of TF-bound regions. Leveraging our phylogenetic framework, we find DNA methylation gain upon evolutionary loss of TF occupancy, indicating coordinated evolution. Furthermore, each methylation pattern has its own evolutionary trajectory reflecting its genomic contexts.

CONCLUSIONS

Our epigenomic analyses indicate a role for DNA methylation in TF binding changes across species including that specific DNA methylation profiles characterize TF binding and are associated with their regulatory activity, chromatin contexts, and evolutionary trajectories.

Aging clock based on nucleosome reorganisation derived from cell-free DNA

Aging induces systematic changes in the distribution of nucleosomes, which affect gene expression programs. Here we reconstructed nucleosome maps based on cell-free DNA (cfDNA) extracted from blood plasma using four cohorts of people of different ages. We show that nucleosomes tend to be separated by larger genomic distances in older people, and age correlates with the nucleosome repeat length (NRL). Furthermore, we developed the first aging clock based on cfDNA nucleosomics. Machine learning based on cfDNA distance distributions allowed predicting person's age with the median absolute error of 3-3.5 years.

Nucleosome reorganisation in breast cancer tissues

BACKGROUND

Nucleosome repositioning in cancer is believed to cause many changes in genome organisation and gene expression. Understanding these changes is important to elucidate fundamental aspects of cancer. It is also important for medical diagnostics based on cell-free DNA (cfDNA), which originates from genomic DNA regions protected from digestion by nucleosomes.

RESULTS

We have generated high-resolution nucleosome maps in paired tumour and normal tissues from the same breast cancer patients using MNase-assisted histone H3 ChIP-seq and compared them with the corresponding cfDNA from blood plasma. This analysis has detected single-nucleosome repositioning at key regulatory regions in a patient-specific manner and common cancer-specific patterns across patients. The nucleosomes gained in tumour versus normal tissue were particularly informative of cancer pathways, with ~ 20-fold enrichment at CpG islands, a large fraction of which marked promoters of genes encoding DNA-binding proteins. The tumour tissues were characterised by a 5-10 bp decrease in the average distance between nucleosomes (nucleosome repeat length, NRL), which is qualitatively similar to the differences between pluripotent and differentiated cells. This effect was correlated with gene activity, differential DNA methylation and changes in local occupancy of linker histone variants H1.4 and H1X.

CONCLUSIONS

Our study offers a novel resource of high-resolution nucleosome maps in breast cancer patients and reports for the first time the effect of systematic decrease of NRL in paired tumour versus normal breast tissues from the same patient. Our findings provide a new mechanistic understanding of nucleosome repositioning in tumour tissues that can be valuable for patient diagnostics, stratification and monitoring.

In vitro reconstitution of chromatin domains shows a role for nucleosome positioning in 3D genome organization

Eukaryotic genomes are organized into chromatin domains. The molecular mechanisms driving the formation of these domains are difficult to dissect in vivo and remain poorly understood. Here we reconstitute Saccharomyces cerevisiae chromatin in vitro and determine its 3D organization at subnucleosome resolution by micrococcal nuclease-based chromosome conformation capture and molecular dynamics simulations. We show that regularly spaced and phased nucleosome arrays form chromatin domains in vitro that resemble domains in vivo. This demonstrates that neither loop extrusion nor transcription is required for basic domain formation in yeast. In addition, we find that the boundaries of reconstituted domains correspond to nucleosome-free regions and that insulation strength scales with their width. Finally, we show that domain compaction depends on nucleosome linker length, with longer linkers forming more compact structures. Together, our results demonstrate that regular nucleosome positioning is important for the formation of chromatin domains and provide a proof-of-principle for bottom-up 3D genome studies.

Corticosteroid-induced chromatin loop dynamics at the FKBP5 gene

FKBP5 is a 115-kb-long glucocorticoid-inducible gene implicated in psychiatric disorders. To investigate the complexities of chromatin interaction frequencies at the FKBP5 topologically associated domain (TAD), we deployed 15 one-to-all chromatin capture viewpoints near gene promoters, enhancers, introns, and CTCF-loop anchors. This revealed a "one-TAD-one-gene" structure encompassing the FKBP5 promoter and its enhancers. The FKBP5 promoter and its two glucocorticoid-stimulated enhancers roam the entire TAD while displaying subtle cell type-specific interactomes. The FKBP5 TAD consists of two nested CTCF loops that are coordinated by one CTCF site in the eighth intron of FKBP5 and another beyond its polyadenylation site, 61 kb further. Loop extension correlates with transcription increases through the intronic CTCF site. This is efficiently compensated for, since the short loop is restored even under high transcription regimes. The boundaries of the FKBP5 TAD consist of divergent CTCF site patterns, harbor multiple smaller genes, and are resilient to glucocorticoid stimulation. Interestingly, both FKBP5 TAD boundaries harbor H3K27me3-marked heterochromatin blocks that may reinforce them. We propose that cis-acting genetic and epigenetic polymorphisms underlying FKBP5 expression variation are likely to reside within a 240-kb region that consists of the FKBP5 TAD, its left sub-TAD, and both its boundaries.

Inferring CTCF-binding patterns and anchored loops across human tissues and cell types

CCCTC-binding factor (CTCF) is a transcription regulator with a complex role in gene regulation. The recognition and effects of CTCF on DNA sequences, chromosome barriers, and enhancer blocking are not well understood. Existing computational tools struggle to assess the regulatory potential of CTCF-binding sites and their impact on chromatin loop formation. Here we have developed a deep-learning model, DeepAnchor, to accurately characterize CTCF binding using high-resolution genomic/epigenomic features. This has revealed distinct chromatin and sequence patterns for CTCF-mediated insulation and looping. An optimized implementation of a previous loop model based on DeepAnchor score excels in predicting CTCF-anchored loops. We have established a compendium of CTCF-anchored loops across 52 human tissue/cell types, and this suggests that genomic disruption of these loops could be a general mechanism of disease pathogenesis. These computational models and resources can help investigate how CTCF-mediated cis-regulatory elements shape context-specific gene regulation in cell development and disease progression.

Predicting scale-dependent chromatin polymer properties from systematic coarse-graining

Simulating chromatin is crucial for predicting genome organization and dynamics. Although coarse-grained bead-spring polymer models are commonly used to describe chromatin, the relevant bead dimensions, elastic properties, and the nature of inter-bead potentials are unknown. Using nucleosome-resolution contact probability (Micro-C) data, we systematically coarse-grain chromatin and predict quantities essential for polymer representation of chromatin. We compute size distributions of chromatin beads for different coarse-graining scales, quantify fluctuations and distributions of bond lengths between neighboring regions, and derive effective spring constant values. Unlike the prevalent notion, our findings argue that coarse-grained chromatin beads must be considered as soft particles that can overlap, and we derive an effective inter-bead soft potential and quantify an overlap parameter. We also compute angle distributions giving insights into intrinsic folding and local bendability of chromatin. While the nucleosome-linker DNA bond angle naturally emerges from our work, we show two populations of local structural states. The bead sizes, bond lengths, and bond angles show different mean behavior at Topologically Associating Domain (TAD) boundaries and TAD interiors. We integrate our findings into a coarse-grained polymer model and provide quantitative estimates of all model parameters, which can serve as a foundational basis for all future coarse-grained chromatin simulations.

Histone exchange sensors reveal variant specific dynamics in mouse embryonic stem cells

Eviction of histones from nucleosomes and their exchange with newly synthesized or alternative variants is a central epigenetic determinant. Here, we define the genome-wide occupancy and exchange pattern of canonical and non-canonical histone variants in mouse embryonic stem cells by genetically encoded exchange sensors. While exchange of all measured variants scales with transcription, we describe variant-specific associations with transcription elongation and Polycomb binding. We found considerable exchange of H3.1 and H2B variants in heterochromatin and repeat elements, contrasting the occupancy and little exchange of H3.3 in these regions. This unexpected association between H3.3 occupancy and exchange of canonical variants is also evident in active promoters and enhancers, and further validated by reduced H3.1 dynamics following depletion of H3.3-specific chaperone, HIRA. Finally, analyzing transgenic mice harboring H3.1 or H3.3 sensors demonstrates the vast potential of this system for studying histone exchange and its impact on gene expression regulation in vivo.

Human SMARCA5 is continuously required to maintain nucleosome spacing

Genetic models suggested that SMARCA5 was required for DNA-templated events including transcription, DNA replication, and DNA repair. We engineered a degron tag into the endogenous alleles of SMARCA5, a catalytic component of the imitation switch complexes in three different human cell lines to define the effects of rapid degradation of this key regulator. Degradation of SMARCA5 was associated with a rapid increase in global nucleosome repeat length, which may allow greater chromatin compaction. However, there were few changes in nascent transcription within the first 6 h of degradation. Nevertheless, we demonstrated a requirement for SMARCA5 to control nucleosome repeat length at G1/S and during the S phase. SMARCA5 co-localized with CTCF and H2A.Z, and we found a rapid loss of CTCF DNA binding and disruption of nucleosomal phasing around CTCF binding sites. This spatiotemporal analysis indicates that SMARCA5 is continuously required for maintaining nucleosomal spacing.

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

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

Cohesin and CTCF complexes mediate contacts in chromatin loops depending on nucleosome positions

The spatial organization of the eukaryotic genome plays an important role in regulating transcriptional activity. In the nucleus, chromatin forms loops that assemble into fundamental units called topologically associating domains that facilitate or inhibit long-range contacts. These loops are formed and held together by the ring-shaped cohesin protein complex, and this can involve binding of CCCTC-binding factor (CTCF). High-resolution conformation capture experiments provide the frequency at which two DNA fragments physically associate in three-dimensional space. However, technical limitations of this approach, such as low throughput, low resolution, or noise in contact maps, make data interpretation and identification of chromatin intraloop contacts, e.g., between distal regulatory elements and their target genes, challenging. Herein, an existing coarse-grained model of chromatin at single-nucleosome resolution was extended by integrating potentials describing CTCF and cohesin. We performed replica-exchange Monte Carlo simulations with regularly spaced nucleosomes and experimentally determined nucleosome positions in the presence of cohesin-CTCF, as well as depleted systems as controls. In fully extruded loops caused by the presence of cohesin and CTCF, the number of contacts within the formed loops was increased. The number and types of these contacts were impacted by the nucleosome distribution and loop size. Microloops were observed within cohesin-mediated loops due to thermal fluctuations without additional influence of other factors, and the number, size, and shape of microloops were determined by nucleosome distribution and loop size. Nucleosome positions directly affect the spatial structure and contact probability within a loop, with presumed consequences for transcriptional activity.

Deciphering the mechanical code of the genome and epigenome

Diverse DNA-deforming processes are impacted by the local mechanical and structural properties of DNA, which in turn depend on local sequence and epigenetic modifications. Deciphering this mechanical code (that is, this dependence) has been challenging due to the lack of high-throughput experimental methods. Here we present a comprehensive characterization of the mechanical code. Utilizing high-throughput measurements of DNA bendability via loop-seq, we quantitatively established how the occurrence and spatial distribution of dinucleotides, tetranucleotides and methylated CpG impact DNA bendability. We used our measurements to develop a physical model for the sequence and methylation dependence of DNA bendability. We validated the model by performing loop-seq on mouse genomic sequences around transcription start sites and CTCF-binding sites. We applied our model to test the predictions of all-atom molecular dynamics simulations and to demonstrate that sequence and epigenetic modifications can mechanically encode regulatory information in diverse contexts.

Heterogeneous interactions and polymer entropy decide organization and dynamics of chromatin domains

Chromatin is known to be organized into multiple domains of varying sizes and compaction. While these domains are often imagined as static structures, they are highly dynamic and show cell-to-cell variability. Since processes such as gene regulation and DNA replication occur in the context of these domains, it is important to understand their organization, fluctuation, and dynamics. To simulate chromatin domains, one requires knowledge of interaction strengths among chromatin segments. Here, we derive interaction-strength parameters from experimentally known contact maps and use them to predict chromatin organization and dynamics. Taking two domains on the human chromosome as examples, we investigate its three-dimensional organization, size/shape fluctuations, and dynamics of different segments within a domain, accounting for hydrodynamic effects. Considering different cell types, we quantify changes in interaction strengths and chromatin shape fluctuations in different epigenetic states. Perturbing the interaction strengths systematically, we further investigate how epigenetic-like changes can alter the spatio-temporal nature of the domains. Our results show that heterogeneous weak interactions are crucial in determining the organization of the domains. Computing effective stiffness and relaxation times, we investigate how perturbations in interactions affect the solid- and liquid-like nature of chromatin domains. Quantifying dynamics of chromatin segments within a domain, we show how the competition between polymer entropy and interaction energy influence the timescales of loop formation and maintenance of stable loops.

KSHV episome tethering sites on host chromosomes and regulation of latency-lytic switch by CHD4

Kaposi sarcoma-associated herpesvirus (KSHV) establishes a latent infection in the cell nucleus, but where KSHV episomal genomes are tethered and the mechanisms underlying KSHV lytic reactivation are unclear. Here, we study the nuclear microenvironment of KSHV episomes and show that the KSHV latency-lytic replication switch is regulated via viral long non-coding (lnc)RNA-CHD4 (chromodomain helicase DNA binding protein 4) interaction. KSHV episomes localize with CHD4 and ADNP proteins, components of the cellular ChAHP complex. The CHD4 and ADNP proteins occupy the 5'-region of the highly inducible lncRNAs and terminal repeats of the KSHV genome together with latency-associated nuclear antigen (LANA). Viral lncRNA binding competes with CHD4 DNA binding, and KSHV reactivation sequesters CHD4 from the KSHV genome, which is also accompanied by detachment of KSHV episomes from host chromosome docking sites. We propose a model in which robust KSHV lncRNA expression determines the latency-lytic decision by regulating LANA/CHD4 binding to KSHV episomes.

Analysis of sub-kilobase chromatin topology reveals nano-scale regulatory interactions with variable dependence on cohesin and CTCF

Enhancers and promoters predominantly interact within large-scale topologically associating domains (TADs), which are formed by loop extrusion mediated by cohesin and CTCF. However, it is unclear whether complex chromatin structures exist at sub-kilobase-scale and to what extent fine-scale regulatory interactions depend on loop extrusion. To address these questions, we present an MNase-based chromosome conformation capture (3C) approach, which has enabled us to generate the most detailed local interaction data to date (20 bp resolution) and precisely investigate the effects of cohesin and CTCF depletion on chromatin architecture. Our data reveal that cis-regulatory elements have distinct internal nano-scale structures, within which local insulation is dependent on CTCF, but which are independent of cohesin. In contrast, we find that depletion of cohesin causes a subtle reduction in longer-range enhancer-promoter interactions and that CTCF depletion can cause rewiring of regulatory contacts. Together, our data show that loop extrusion is not essential for enhancer-promoter interactions, but contributes to their robustness and specificity and to precise regulation of gene expression.

DNA sequence-dependent formation of heterochromatin nanodomains

The mammalian epigenome contains thousands of heterochromatin nanodomains (HNDs) marked by di- and trimethylation of histone H3 at lysine 9 (H3K9me2/3), which have a typical size of 3–10 nucleosomes. However, what governs HND location and extension is only partly understood. Here, we address this issue by introducing the chromatin hierarchical lattice framework (ChromHL) that predicts chromatin state patterns with single-nucleotide resolution. ChromHL is applied to analyse four HND types in mouse embryonic stem cells that are defined by histone methylases SUV39H1/2 or GLP, transcription factor ADNP or chromatin remodeller ATRX. We find that HND patterns can be computed from PAX3/9, ADNP and LINE1 sequence motifs as nucleation sites and boundaries that are determined by DNA sequence (e.g. CTCF binding sites), cooperative interactions between nucleosomes as well as nucleosome-HP1 interactions. Thus, ChromHL rationalizes how patterns of H3K9me2/3 are established and changed via the activity of protein factors in processes like cell differentiation.

The ability to predict epigenetic regulation is an important challenge in biology. Here the authors describe heterochromatin nanodomains (HNDs) and compare four different types of H3K9me2/3-marked HNDs in mouse embryonic stem cells. They further develop a computational framework to predict genome-wide HND maps from DNA sequence and protein concentrations, at single-nucleotide resolution.

Prediction and comparative analysis of CTCF binding sites based on a first principle approach

We calculated the patterns for the CCCTC transcription factor (CTCF) binding sites across many genomes on a first principle approach. The validation of the first principle method was done on the human as well as on the mouse genome. The predicted human CTCF binding sites are consistent with the consensus sequence, ChIP-seq data for the K562 cell, nucleosome positions for IMR90 cell as well as the CTCF binding sites in the mouse HOXA gene. The analysis of Homo sapiens, Mus musculus, Sus scrofa, Capra hircus and Drosophila melanogaster whole genomes shows: binding sites are organized in cluster-like groups, where two consecutive sites obey a power-law with coefficient ranging from to 0.3292 0.0068 to 0.5409 0.0064; the distance between these groups varies from 18.08 0.52kbp to 42.1 2.0kbp. The genome of Aedes aegypti does not show a power law, but 19.9% of binding sites are 144 4 and 287 5bp distant of each other. We run negative tests, confirming the under-representation of CTCF binding sites in Caenorhabditis elegans, Plasmodium falciparum and Arabidopsis thaliana complete genomes.

NucPosDB: a database of nucleosome positioning in vivo and nucleosomics of cell-free DNA

Nucleosome positioning is involved in many gene regulatory processes happening in the cell, and it may change as cells differentiate or respond to the changing microenvironment in a healthy or diseased organism. One important implication of nucleosome positioning in clinical epigenetics is its use in the “nucleosomics” analysis of cell-free DNA (cfDNA) for the purpose of patient diagnostics in liquid biopsies. The rationale for this is that the apoptotic nucleases that digest chromatin of the dying cells mostly cut DNA between nucleosomes. Thus, the short pieces of DNA in body fluids reflect the positions of nucleosomes in the cells of origin. Here, we report a systematic nucleosomics database — NucPosDB — curating published nucleosome positioning datasets in vivo as well as datasets of sequenced cell-free DNA (cfDNA) that reflect nucleosome positioning in situ in the cells of origin. Users can select subsets of the database by a number of criteria and then obtain raw or processed data. NucPosDB also reports the originally determined regions with stable nucleosome occupancy across several individuals with a given condition. An additional section provides a catalogue of computational tools for the analysis of nucleosome positioning or cfDNA experiments and theoretical algorithms for the prediction of nucleosome positioning preferences from DNA sequence. We provide an overview of the field, describe the structure of the database in this context, and demonstrate data variability using examples of different medical conditions. NucPosDB is useful both for the analysis of fundamental gene regulation processes and the training of computational models for patient diagnostics based on cfDNA. The database currently curates ~ 400 publications on nucleosome positioning in cell lines and in situ as well as cfDNA from > 10,000 patients and healthy volunteers. For open-access cfDNA datasets as well as key MNase-seq datasets in human cells, NucPosDB allows downloading processed mapped data in addition to the regions with stable nucleosome occupancy. NucPosDB is available at https://generegulation.org/nucposdb/.

Epigenetic Regulation of Human T-Cell Leukemia Virus Gene Expression

Viral and cellular gene expression are regulated by epigenetic alterations, including DNA methylation, histone modifications, nucleosome positioning, and chromatin looping. Human T-cell leukemia virus type 1 (HTLV-1) is a pathogenic retrovirus associated with inflammatory disorders and T-cell lymphoproliferative malignancy. The transforming activity of HTLV-1 is driven by the viral oncoprotein Tax, which acts as a transcriptional activator of the cAMP response element-binding protein (CREB) and nuclear factor kappa B (NFκB) pathways. The epigenetic effects of Tax and the induction of lymphoproliferative malignancy include alterations in DNA methylation and histone modifications. In addition, alterations in nucleosome positioning and DNA looping also occur in HTLV-1-induced malignant cells. A mechanistic definition of these effects will pave the way to new therapies for HTLV-1-associated disorders.

Nucleosome-induced homology recognition in chromatin

One of the least understood properties of chromatin is the ability of its similar regions to recognize each other through weak interactions. Theories based on electrostatic interactions between helical macromolecules suggest that the ability to recognize sequence homology is an innate property of the non-ideal helical structure of DNA. However, this theory does not account for the nucleosomal packing of DNA. Can homologous DNA sequences recognize each other while wrapped up in the nucleosomes? Can structural homology arise at the level of nucleosome arrays? Here, we present a theoretical model for the recognition potential well between chromatin fibres sliding against each other. This well is different from the one predicted for bare DNA; the minima in energy do not correspond to literal juxtaposition, but are shifted by approximately half the nucleosome repeat length. The presence of this potential well suggests that nucleosome positioning may induce mutual sequence recognition between chromatin fibres and facilitate the formation of chromatin nanodomains. This has implications for nucleosome arrays enclosed between CTCF-cohesin boundaries, which may form stiffer stem-like structures instead of flexible entropically favourable loops. We also consider switches between chromatin states, e.g. through acetylation/deacetylation of histones, and discuss nucleosome-induced recognition as a precursory stage of genetic recombination.

Chromatin stability generated by stochastic binding and unbinding of cross-linkers at looping sites revealed by Markov models

Chromatin loops inside the nucleus can be stable for a very long time, which remains poorly understood. Such a time is crucial for chromatin organization maintenance and stability. We explore here several physical scenarios, where loop maintenance is due to diffusing cross-linkers (cohesin stabilized by two CTCF molecules) that can bind and unbind at the base of chromatin loops. Using a Markov chain approach to coarse-grain the binding and unbinding, we consider that a stable loop disappears when the last cross-linker is unbound. We derive expressions for this last passage time that we use to quantify the loop stability for various parameters, such as the chemical rate constant or the number of cross-linkers. The present analysis suggests that the balance between binding and unbinding events regulates the number of cross-linkers in place, based on a positive feed-back mechanism that stabilizes the loop over long-time. To conclude, we found that short- and long-lasting stable loops can vary from minutes to the entire cell cycle lifetime, when the number of cross-linkers increases from 1 to 10. This result suggests that a large spectrum of loop time scales is expected with such a few numbers of cross-linkers per local binding sites.

At the Crossroad of Gene Regulation and Genome Organization: Potential Roles for ATP-Dependent Chromatin Remodelers in the Regulation of CTCF-Mediated 3D Architecture

Simple Summary

The way DNA is packaged in the nucleus of a cell is important for when and how genes are expressed. There are many levels of packaging, and new techniques have revealed that long-range interactions are important for both promoting and restricting the transcription of genes. Some long-range interactions are mediated by physical loops in the genome where, like a rubber band, the ring-shaped cohesin complex loops sections of DNA bound by CCCTC-binding factor (CTCF). Both cohesin and CTCF act on DNA, and increasing evidence indicates that their function is inhibited by nucleosomes bound to the DNA. In this review, we summarize the current knowledge of how individual chromatin remodelers, which utilize ATP to move nucleosomes on DNA, facilitate or inhibit cohesin/CTCF-dependent looping interactions.

Abstract

In higher order organisms, the genome is assembled into a protein-dense structure called chromatin. Chromatin is spatially organized in the nucleus through hierarchical folding, which is tightly regulated both in cycling cells and quiescent cells. Assembly and folding are not one-time events in a cell’s lifetime; rather, they are subject to dynamic shifts to allow changes in transcription, DNA replication, or DNA damage repair. Chromatin is regulated at many levels, and recent tools have permitted the elucidation of specific factors involved in the maintenance and regulation of the three-dimensional (3D) genome organization. In this review/perspective, we aim to cover the potential, but relatively unelucidated, crosstalk between 3D genome architecture and the ATP-dependent chromatin remodelers with a specific focus on how the architectural proteins CTCF and cohesin are regulated by chromatin remodeling.

Nucleosome Positioning and Spacing: From Mechanism to Function

Eukaryotes associate their genomes with histone proteins, forming nucleosome particles. Nucleosomes regulate and protect the genetic information. They often assemble into evenly spaced arrays of nucleosomes. These regular nucleosome arrays cover significant portions of the genome, in particular over genes. The presence of these evenly spaced nucleosome arrays is highly conserved throughout the entire eukaryotic domain. Here, we review the mechanisms behind the establishment of this primary structure of chromatin with special emphasis on the biogenesis of evenly spaced nucleosome arrays. We highlight the roles that transcription, nucleosome remodelers, DNA sequence, and histone density play towards the formation of evenly spaced nucleosome arrays and summarize our current understanding of their cellular functions. We end with key unanswered questions that remain to be explored to obtain an in-depth understanding of the biogenesis and function of the nucleosome landscape.

CTCF as a boundary factor for cohesin-mediated loop extrusion: evidence for a multi-step mechanism

Mammalian genome structure is closely linked to function. At the scale of kilobases to megabases, CTCF and cohesin organize the genome into chromatin loops. Mechanistically, cohesin is proposed to extrude chromatin loops bidirectionally until it encounters occupied CTCF DNA-binding sites. Curiously, loops form predominantly between CTCF binding sites in a convergent orientation. How CTCF interacts with and blocks cohesin extrusion in an orientation-specific manner has remained a mechanistic mystery. Here, we review recent papers that have shed light on these processes and suggest a multi-step interaction between CTCF and cohesin. This interaction may first involve a pausing step, where CTCF halts cohesin extrusion, followed by a stabilization step of the CTCF-cohesin complex, resulting in a chromatin loop. Finally, we discuss our own recent studies on an internal RNA-Binding Region (RBRi) in CTCF to elucidate its role in regulating CTCF clustering, target search mechanisms and chromatin loop formation and future challenges.

The Accidental Ally: Nucleosome Barriers Can Accelerate Cohesin-Mediated Loop Formation in Chromatin

An important question in the context of the three-dimensional organization of chromosomes is the mechanism of formation of large loops between distant basepairs. Recent experiments suggest that the formation of loops might be mediated by loop extrusion factor proteins such as cohesin. Experiments on cohesin have shown that cohesins walk diffusively on the DNA and that nucleosomes act as obstacles to the diffusion, lowering the permeability and hence reducing the effective diffusion constant. An estimation of the times required to form the loops of typical sizes seen in Hi-C experiments using these low-effective-diffusion constants leads to times that are unphysically large. The puzzle then is the following: how does a cohesin molecule diffusing on the DNA backbone achieve speeds necessary to form the large loops seen in experiments? We propose a simple answer to this puzzle and show that although at low densities, nucleosomes act as barriers to cohesin diffusion, beyond a certain concentration they can reduce loop formation times because of a subtle interplay between the nucleosome size and the mean linker length. This effect is further enhanced on considering stochastic binding kinetics of nucleosomes on the DNA backbone and leads to predictions of lower loop formation times than might be expected from a naive obstacle picture of nucleosomes.

Chromatin Remodelers in the 3D Nuclear Compartment

Chromatin remodeling complexes (CRCs) use ATP hydrolysis to maintain correct expression profiles, chromatin stability, and inherited epigenetic states. More than 20 CRCs have been described to date, which encompass four large families defined by their ATPase subunits. These complexes and their subunits are conserved from yeast to humans through evolution. Their activities depend on their catalytic subunits which through ATP hydrolysis provide the energy necessary to fulfill cellular functions such as gene transcription, DNA repair, and transposon silencing. These activities take place at the first levels of chromatin compaction, and CRCs have been recognized as essential elements of chromatin dynamics. Recent studies have demonstrated an important role for these complexes in the maintenance of higher order chromatin structure. In this review, we present an overview of the organization of the genome within the cell nucleus, the different levels of chromatin compaction, and importance of the architectural proteins, and discuss the role of CRCs and how their functions contribute to the dynamics of the 3D genome organization.

The ubiquitous transcriptional protein ZNF143 activates a diversity of genes while assisting to organize chromatin structure

Zinc Finger Protein 143 (ZNF143) is a pervasive C2H2 zinc-finger transcriptional activator protein regulating the efficiency of eukaryotic promoter regions. ZNF143 is able to activate transcription at both protein coding genes and small RNA genes transcribed by either RNA polymerase II or RNA polymerase III. Target genes regulated by ZNF143 are involved in an array of different cellular processes including both cancer and development. Although a key player in regulating eukaryotic genes, the molecular mechanism by with ZNF143 binds and activates genes transcribed by two different polymerases is still relatively unknown. In addition to its role as a transcriptional regulator, recent genomics experiments have implicated ZNF143 as a potential co-factor involved in chromatin looping and establishing higher order structure within the genome. This review focuses primarily on possible activation mechanisms of promoters by ZNF143, with less emphasis on the role of ZNF143 in cancer and development, and its function in establishing higher order chromatin contacts within the genome.

Nucleosome positioning and chromatin organization

In our cells, DNA is folded and packed with the help of many proteins into chromatin whose basic unit is a nucleosome-DNA wrapped around octamer of histone proteins. The chain of nucleosomes is further folded and arranged into many layers and has a dynamic organization. How does the complex chromatin organization emerge from interactions among DNA, histones, and non-histone proteins have been a question of great interest. Here we review recent literature that investigated how nucleosome positioning and nucleosome-mediated interactions drive chromatin organization. Unlike our earlier understanding, chromatin is organized into 3D domains of various sizes having irregularly organized nucleosomes. These domains emerge due to heterogeneous nucleosome positioning and diverse inter-nucleosome interactions that vary in space and time.

Linker histone epitopes are hidden by in situ higher-order chromatin structure

BACKGROUND

Histone H1 is the most mobile histone in the cell nucleus. Defining the positions of H1 on chromatin in situ, therefore, represents a challenge. Immunoprecipitation of formaldehyde-fixed and sonicated chromatin, followed by DNA sequencing (xChIP-seq), is traditionally the method for mapping histones onto DNA elements. But since sonication fragmentation precedes ChIP, there is a consequent loss of information about chromatin higher-order structure. Here, we present a new method, xxChIP-seq, employing antibody binding to fixed intact in situ chromatin, followed by extensive washing, a second fixation, sonication and immunoprecipitation. The second fixation is intended to prevent the loss of specifically bound antibody during washing and subsequent sonication and to prevent antibody shifting to epitopes revealed by the sonication process. In many respects, xxChIP-seq is comparable to immunostaining microscopy, which also involves interaction of the primary antibody with fixed and permeabilized intact cells. The only epitopes displayed after immunostaining are the "exposed" epitopes, not "hidden" by the fixation of chromatin higher-order structure. Comparison of immunoprecipitated fragments between xChIP-seq versus xxChIP-seq should indicate which epitopes become inaccessible with fixation and identify their associated DNA elements.

RESULTS

We determined the genomic distribution of histone variants H1.2 and H1.5 in human myeloid leukemia cells HL-60/S4 and compared their epitope exposure by both xChIP-seq and xxChIP-seq, as well as high-resolution microscopy, illustrating the influences of preserved chromatin higher-order structure in situ. We found that xChIP and xxChIP H1 signals are in general negatively correlated, with differences being more pronounced near active regulatory regions. Among the intriguing observations, we find that transcription-related regions and histone PTMs (i.e., enhancers, promoters, CpG islands, H3K4me1, H3K4me3, H3K9ac, H3K27ac and H3K36me3) exhibit significant deficiencies (depletions) in H1.2 and H1.5 xxChIP-seq reads, compared to xChIP-seq. These observations suggest the existence of in situ transcription-related chromatin higher-order structures stabilized by formaldehyde.

CONCLUSION

Comparison of H1 xxChIP-seq to H1 xChIP-seq allows the development of hypotheses on the chromosomal localization of (stabilized) higher-order structure, indicated by the generation of "hidden" H1 epitopes following formaldehyde crosslinking. Changes in H1 epitope exposure surrounding averaged chromosomal binding sites or epigenetic modifications can also indicate whether these sites have chromatin higher-order structure. For example, comparison between averaged active or inactive promoter regions suggests that both regions can acquire stabilized higher-order structure with hidden H1 epitopes. However, the H1 xChIP-seq comparison cannot define their differences. Application of the xxChIP-seq versus H1 xChIP-seq method is particularly relevant to chromatin-associated proteins, such as linker histones, that play dynamic roles in establishing chromatin higher-order structure.