Hierarchies in eukaryotic genome organization: Insights from polymer theory and simulations

Mesoscale Modeling Reveals Hierarchical Looping of Chromatin Fibers Near Gene Regulatory Elements

While it is well-recognized that chromatin loops play an important role in gene regulation, structural details regarding higher order chromatin loops are only emerging. Here we present a systematic study of restrained chromatin loops ranging from 25 to 427 nucleosomes (fibers of 5-80 Kb DNA in length), mimicking gene elements studied by 3C contact data. We find that hierarchical looping represents a stable configuration that can effectively bring distant regions of the GATA-4 gene together, satisfying connections reported by 3C experiments. Additionally, we find that restrained chromatin fibers larger than 100 nucleosomes (∼20Kb) form closed plectonemes, whereas fibers shorter than 100 nucleosomes form simple hairpin loops. By studying the dependence of loop structures on internal parameters, we show that loop features are sensitive to linker histone concentration, loop length, divalent ions, and DNA linker length. Specifically, increasing loop length, linker histone concentration, and divalent ion concentration are associated with increased persistence length (or decreased bending), while varying DNA linker length in a manner similar to experimentally observed "nucleosome free regions" (found near transcription start sites) disrupts intertwining and leads to loop opening and increased persistence length in linker histone depleted (-LH) fibers. Chromatin fiber structure sensitivity to these parameters, all of which vary throughout the cell cycle, tissue type, and species, suggests that caution is warranted when using uniform polymer models to fit chromatin conformation capture genome-wide data. Furthermore, the folding geometry we observe near the transcription initiation site of the GATA-4 gene suggests that hierarchical looping provides a structural mechanism for gene inhibition, and offers tunable parameters for design of gene regulation elements.

Inferring 3D chromatin structure using a multiscale approach based on quaternions

Background

The knowledge of the spatial organisation of the chromatin fibre in cell nuclei helps researchers to understand the nuclear machinery that regulates dna activity. Recent experimental techniques of the type Chromosome Conformation Capture (3c, or similar) provide high-resolution, high-throughput data consisting in the number of times any possible pair of dna fragments is found to be in contact, in a certain population of cells. As these data carry information on the structure of the chromatin fibre, several attempts have been made to use them to obtain high-resolution 3d reconstructions of entire chromosomes, or even an entire genome. The techniques proposed treat the data in different ways, possibly exploiting physical-geometric chromatin models. One popular strategy is to transform contact data into Euclidean distances between pairs of fragments, and then solve a classical distance-to-geometry problem.

Results

We developed and tested a reconstruction technique that does not require translating contacts into distances, thus avoiding a number of related drawbacks. Also, we introduce a geometrical chromatin chain model that allows us to include sound biochemical and biological constraints in the problem. This model can be scaled at different genomic resolutions, where the structures of the coarser models are influenced by the reconstructions at finer resolutions. The search in the solution space is then performed by a classical simulated annealing, where the model is evolved efficiently through quaternion operators. The presence of appropriate constraints permits the less reliable data to be overlooked, so the result is a set of plausible chromatin configurations compatible with both the data and the prior knowledge.

Conclusions

To test our method, we obtained a number of 3d chromatin configurations from hi-c data available in the literature for the long arm of human chromosome 1, and validated their features against known properties of gene density and transcriptional activity. Our results are compatible with biological features not introduced a priori in the problem: structurally different regions in our reconstructions highly correlate with functionally different regions as known from literature and genomic repositories.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0667-0) contains supplementary material, which is available to authorized users.

Chiral symmetry breaking of a double-stranded helical chain through bend-writhe coupling

This paper explores asymmetric elasticity of a double-stranded helical chain, which serves as a minimal model of biopolymers. The model consists of two elastic chains that mutually intertwine in a right-handed manner, forming a double-stranded helix. A simple numerical experiment for structural relaxation, which reduces the total elastic energy of the model monotonically without thermal fluctuations, reveals possible asymmetric elasticity inherent in the helical chain. It is first shown that a short segment of the double-stranded helical chain has a tendency to unwind when it is bent. It is also shown that a short segment of the helical chain has a tendency to writhe in the left direction upon bending. This tendency gives rise to a propensity for a longer segment of the chain to form a left-handed superhelix spontaneously upon bending. Finally, this propensity of the helical chain to form a left-handed superhelix is proposed to be a possible origin of the uniform left-handed wrapping of DNA around nucleosome core particles in nature. The results presented here could provide deeper insights into the roles and significance of helical chirality of biopolymers.

Evaluation of novel design strategies for developing zinc finger nucleases tools for treating human diseases

Zinc finger nucleases (ZFNs) are associated with cell death and apoptosis by binding at countless undesired locations. This cytotoxicity is associated with the binding ability of engineered zinc finger domains to bind dissimilar DNA sequences with high affinity. In general, binding preferences of transcription factors are associated with significant degenerated diversity and complexity which convolutes the design and engineering of precise DNA binding domains. Evolutionary success of natural zinc finger proteins, however, evinces that nature created specific evolutionary traits and strategies, such as modularity and rank-specific recognition to cope with binding complexity that are critical for creating clinical viable tools to precisely modify the human genome. Our findings indicate preservation of general modularity and significant alteration of the rank-specific binding preferences of the three-finger binding domain of transcription factor SP1 when exchanging amino acids in the 2nd finger.

Genome-wide profiling of chromosome interactions in Plasmodium falciparum characterizes nuclear architecture and reconfigurations associated with antigenic variation

Spatial relationships within the eukaryotic nucleus are essential for proper nuclear function. In Plasmodium falciparum, the repositioning of chromosomes has been implicated in the regulation of the expression of genes responsible for antigenic variation, and the formation of a single, peri-nuclear nucleolus results in the clustering of rDNA. Nevertheless, the precise spatial relationships between chromosomes remain poorly understood, because, until recently, techniques with sufficient resolution have been lacking. Here we have used chromosome conformation capture and second-generation sequencing to study changes in chromosome folding and spatial positioning that occur during switches in var gene expression. We have generated maps of chromosomal spatial affinities within the P. falciparum nucleus at 25 Kb resolution, revealing a structured nucleolus, an absence of chromosome territories, and confirming previously identified clustering of heterochromatin foci. We show that switches in var gene expression do not appear to involve interaction with a distant enhancer, but do result in local changes at the active locus. These maps reveal the folding properties of malaria chromosomes, validate known physical associations, and characterize the global landscape of spatial interactions. Collectively, our data provide critical information for a better understanding of gene expression regulation and antigenic variation in malaria parasites.

Recovering ensembles of chromatin conformations from contact probabilities

The 3D higher order organization of chromatin within the nucleus of eukaryotic cells has so far remained elusive. A wealth of relevant information, however, is increasingly becoming available from chromosome conformation capture (3C) and related experimental techniques, which measure the probabilities of contact between large numbers of genomic sites in fixed cells. Such contact probabilities (CPs) can in principle be used to deduce the 3D spatial organization of chromatin. Here, we propose a computational method to recover an ensemble of chromatin conformations consistent with a set of given CPs. Compared with existing alternatives, this method does not require conversion of CPs to mean spatial distances. Instead, we estimate CPs by simulating a physically realistic, bead-chain polymer model of the 30-nm chromatin fiber. We then use an approach from adaptive filter theory to iteratively adjust the parameters of this polymer model until the estimated CPs match the given CPs. We have validated this method against reference data sets obtained from simulations of test systems with up to 45 beads and 4 loops. With additional testing against experiments and with further algorithmic refinements, our approach could become a valuable tool for researchers examining the higher order organization of chromatin.

Lattice animal model of chromosome organization

Polymer models tied together by constraints of looping and confinement have been used to explain many of the observed organizational characteristics of interphase chromosomes. Here we introduce a simple lattice animal representation of interphase chromosomes that combines the features of looping and confinement constraints into a single framework. We show through Monte Carlo simulations that this model qualitatively captures both the leveling off in the spatial distance between genomic markers observed in fluorescent in situ hybridization experiments and the inverse decay in the looping probability as a function of genomic separation observed in chromosome conformation capture experiments. The model also suggests that the collapsed state of chromosomes and their segregation into territories with distinct looping activities might be a natural consequence of confinement.

Nuclear organization and chromatin dynamics in yeast: biophysical models or biologically driven interactions?

Over the past decade, tremendous progress has been made in understanding the spatial organization of genes and chromosomes. Nuclear organization can be thought of as information that is not encoded in DNA, but which nevertheless impacts gene expression. Nuclear organizational influences can be cell-specific and are potentially heritable. Thus, nuclear organization fulfills all the criteria necessary for it to be considered an authentic level of epigenetic information. Chromosomal nuclear organization is primarily dictated by the biophysical properties of chromatin. Diffusion models of polymers confined in the crowded nuclear space accurately recapitulate experimental observation. Diffusion is a Brownian process, which implies that the positions of chromosomes and genes are not defined deterministically but are likely to be dictated by the laws of probability. Despite the small size of their nuclei, budding yeast have been instrumental in discovering how epigenetic information is encoded in the spatial organization of the genome. The relatively simple organization of the yeast nucleus and the very high number of genetically identical cells that can be observed under fluorescent microscopy allow statistically robust definitions of the gene and chromosome positions in the nuclear space to be constructed. In this review, we will focus on how the spatial organization of the chromatin in the yeast nucleus might impact transcription. This article is part of a Special Issue entitled: Nuclear Transport and RNA Processing.

Chromatin topology and the regulation of antigen receptor assembly

During an organism's ontogeny and in the adult, each B and T lymphocyte generates a unique antigen receptor, thereby creating the organism's ability to respond to a vast number of different antigens. The antigen receptor loci are organized into distinct regions that contain multiple variable (V), diversity (D), and/or joining (J) and constant (C) coding elements that are scattered across large genomic regions. In this review, we discuss the epigenetic modifications that take place in the different antigen receptor loci, the chromatin structure adopted by the antigen receptor loci to allow recombination of elements separated by large genomic distances, and the relationship between epigenetics and chromatin structure and how they relate to the generation of antigen receptor diversity.