Nucleosome positioning by genomic excluding-energy barriers

Nucleosome sliding can influence the spreading of histone modifications

Nucleosomes are the fundamental building blocks of chromatin that not only help in the folding of chromatin, but also in carrying epigenetic information. It is known that nucleosome sliding is responsible for dynamically organizing chromatin structure and the resulting gene regulation. Since sliding can move two neighboring nucleosomes physically close or away, can it play a role in the spreading of histone modifications? We investigate this by simulating a stochastic model that couples nucleosome dynamics with the kinetics of histone modifications. We show that the sliding of nucleosomes can affect the modification pattern as well as the time it takes to modify a given region of chromatin. Exploring different nucleosome densities and modification kinetic parameters, we show that nucleosome sliding can be important for creating histone modification domains. Our model predicts that nucleosome density coupled with sliding dynamics can create an asymmetric histone modification profile around regulatory regions. We also compute the probability distribution of modified nucleosomes and relaxation kinetics of modifications. Our predictions are comparable with known experimental results.

Coupling between Sequence-Mediated Nucleosome Organization and Genome Evolution

The nucleosome is a major modulator of DNA accessibility to other cellular factors. Nucleosome positioning has a critical importance in regulating cell processes such as transcription, replication, recombination or DNA repair. The DNA sequence has an influence on the position of nucleosomes on genomes, although other factors are also implicated, such as ATP-dependent remodelers or competition of the nucleosome with DNA binding proteins. Different sequence motifs can promote or inhibit the nucleosome formation, thus influencing the accessibility to the DNA. Sequence-encoded nucleosome positioning having functional consequences on cell processes can then be selected or counter-selected during evolution. We review the interplay between sequence evolution and nucleosome positioning evolution. We first focus on the different ways to encode nucleosome positions in the DNA sequence, and to which extent these mechanisms are responsible of genome-wide nucleosome positioning in vivo. Then, we discuss the findings about selection of sequences for their nucleosomal properties. Finally, we illustrate how the nucleosome can directly influence sequence evolution through its interactions with DNA damage and repair mechanisms. This review aims to provide an overview of the mutual influence of sequence evolution and nucleosome positioning evolution, possibly leading to complex evolutionary dynamics.

Effects of size, cooperativity, and competitive binding on protein positioning on DNA

Accurate positioning of proteins on chromosomal DNA is crucial for its proper organization as well as gene transcription regulation. Recent experiments revealed existence of periodic patterns of nucleoprotein complexes on DNA, which frequently cannot be explained by sequence-dependent binding of proteins. Previous theoretical studies suggest that such patterns typically emerge as a result of the proteins' volume-exclusion effect. However, the role of other physical factors in patterns' formation, such as the length of DNA, its sequence heterogeneity, and protein binding cooperativity/binding competition to DNA, remains unclear. To address these less understood yet important aspects, we investigated potential effects of these factors on protein positioning on finite-size DNA by using transfer-matrix calculations. It has been found that upon binding to DNA, proteins form oscillatory patterns that span over the length of up to ∼10 times the size of the protein binding site, with the shape of the patterns being strongly dependent on the length of DNA and the proteins' binding cooperativity to DNA. Furthermore, calculations showed that small variations in the proteins' affinity to DNA due to its sequence heterogeneity do not much change the main geometric characteristics of the observed protein patterns. Finally, competition between two different types of proteins for binding to DNA has been found to lead to formation of highly diverse and complex alternating positioning of the two proteins. Altogether, these results provide new insights into the roles of physicochemical properties of proteins, the DNA length, and DNA-binding competition between proteins in formation of protein positioning patterns on DNA.

The implication of DNA bending energy for nucleosome positioning and sliding

Nucleosome not only directly affects cellular processes, such as DNA replication, recombination, and transcription, but also severs as a fundamentally important target of epigenetic modifications. Our previous study indicated that the bending property of DNA is important in nucleosome formation, particularly in predicting the dyad positions of nucleosomes on a DNA segment. Here, we investigated the role of bending energy in nucleosome positioning and sliding in depth to decipher sequence-directed mechanism. The results show that bending energy is a good physical index to predict the free energy in the process of nucleosome reconstitution in vitro. Our data also imply that there are at least 20% of the nucleosomes in budding yeast do not adopt canonical positioning, in which underlying sequences wrapped around histones are structurally symmetric. We also revealed distinct patterns of bending energy profile for distinctly organized chromatin structures, such as well-positioned nucleosomes, fuzzy nucleosomes, and linker regions and discussed nucleosome sliding in terms of bending energy. We proposed that the stability of a nucleosome is positively correlated with the strength of the bending anisotropy of DNA segment, and both accessibility and directionality of nucleosome sliding is likely to be modulated by diverse patterns of DNA bending energy profile.

Coupling of replisome movement with nucleosome dynamics can contribute to the parent-daughter information transfer

Positioning of nucleosomes along the genomic DNA is crucial for many cellular processes that include gene regulation and higher order packaging of chromatin. The question of how nucleosome-positioning information from a parent chromatin gets transferred to the daughter chromatin is highly intriguing. Accounting for experimentally known coupling between replisome movement and nucleosome dynamics, we propose a model that can obtain de novo nucleosome assembly similar to what is observed in recent experiments. Simulating nucleosome dynamics during replication, we argue that short pausing of the replication fork, associated with nucleosome disassembly, can be a event crucial for communicating nucleosome positioning information from parent to daughter. We show that the interplay of timescales between nucleosome disassembly (τ_p) at the replication fork and nucleosome sliding behind the fork (τ_s) can give rise to a rich ‘phase diagram’ having different inherited patterns of nucleosome organization. Our model predicts that only when τ_p ≥ τ_s the daughter chromatin can inherit nucleosome positioning of the parent.

Evidence for DNA Sequence Encoding of an Accessible Nucleosomal Array across Vertebrates

Nucleosome-depleted regions around which nucleosomes order following the "statistical" positioning scenario were recently shown to be encoded in the DNA sequence in human. This intrinsic nucleosomal ordering strongly correlates with oscillations in the local GC content as well as with the interspecies and intraspecies mutation profiles, revealing the existence of both positive and negative selection. In this letter, we show that these predicted nucleosome inhibitory energy barriers (NIEBs) with compacted neighboring nucleosomes are indeed ubiquitous to all vertebrates tested. These 1 kb-sized chromatin patterns are widely distributed along vertebrate chromosomes, overall covering more than a third of the genome. We have previously observed in human deviations from neutral evolution at these genome-wide distributed regions, which we interpreted as a possible indication of the selection of an open, accessible, and dynamic nucleosomal array to constitutively facilitate the epigenetic regulation of nuclear functions in a cell-type-specific manner. As a first, very appealing observation supporting this hypothesis, we report evidence of a strong association between NIEB borders and the poly(A) tails of Alu sequences in human. These results suggest that NIEBs provide adequate chromatin patterns favorable to the integration of Alu retrotransposons and, more generally to various transposable elements in the genomes of primates and other vertebrates.

Evidence of selection for an accessible nucleosomal array in human

BACKGROUND

Recently, a physical model of nucleosome formation based on sequence-dependent bending properties of the DNA double-helix has been used to reveal some enrichment of nucleosome-inhibiting energy barriers (NIEBs) nearby ubiquitous human "master" replication origins. Here we use this model to predict the existence of about 1.6 millions NIEBs over the 22 human autosomes.

RESULTS

We show that these high energy barriers of mean size 153 bp correspond to nucleosome-depleted regions (NDRs) in vitro, as expected, but also in vivo. On either side of these NIEBs, we observe, in vivo and in vitro, a similar compacted nucleosome ordering, suggesting an absence of chromatin remodeling. This nucleosomal ordering strongly correlates with oscillations of the GC content as well as with the interspecies and intraspecies mutation profiles along these regions. Comparison of these divergence rates reveals the existence of both positive and negative selections linked to nucleosome positioning around these intrinsic NDRs. Overall, these NIEBs and neighboring nucleosomes cover 37.5 % of the human genome where nucleosome occupancy is stably encoded in the DNA sequence. These 1 kb-sized regions of intrinsic nucleosome positioning are equally found in GC-rich and GC-poor isochores, in early and late replicating regions, in intergenic and genic regions but not at gene promoters.

CONCLUSION

The source of selection pressure on the NIEBs has yet to be resolved in future work. One possible scenario is that these widely distributed chromatin patterns have been selected in human to impair the condensation of the nucleosomal array into the 30 nm chromatin fiber, so as to facilitate the epigenetic regulation of nuclear functions in a cell-type-specific manner.

A deformation energy-based model for predicting nucleosome dyads and occupancy

Nucleosome plays an essential role in various cellular processes, such as DNA replication, recombination, and transcription. Hence, it is important to decode the mechanism of nucleosome positioning and identify nucleosome positions in the genome. In this paper, we present a model for predicting nucleosome positioning based on DNA deformation, in which both bending and shearing of the nucleosomal DNA are considered. The model successfully predicted the dyad positions of nucleosomes assembled in vitro and the in vitro map of nucleosomes in Saccharomyces cerevisiae. Applying the model to Caenorhabditis elegans and Drosophila melanogaster, we achieved satisfactory results. Our data also show that shearing energy of nucleosomal DNA outperforms bending energy in nucleosome occupancy prediction and the ability to predict nucleosome dyad positions is attributed to bending energy that is associated with rotational positioning of nucleosomes.

Theoretical estimates of exposure timescales of protein binding sites on DNA regulated by nucleosome kinetics

It is being increasingly realized that nucleosome organization on DNA crucially regulates DNA-protein interactions and the resulting gene expression. While the spatial character of the nucleosome positioning on DNA has been experimentally and theoretically studied extensively, the temporal character is poorly understood. Accounting for ATPase activity and DNA-sequence effects on nucleosome kinetics, we develop a theoretical method to estimate the time of continuous exposure of binding sites of non-histone proteins (e.g. transcription factors and TATA binding proteins) along any genome. Applying the method to Saccharomyces cerevisiae, we show that the exposure timescales are determined by cooperative dynamics of multiple nucleosomes, and their behavior is often different from expectations based on static nucleosome occupancy. Examining exposure times in the promoters of GAL1 and PHO5, we show that our theoretical predictions are consistent with known experiments. We apply our method genome-wide and discover huge gene-to-gene variability of mean exposure times of TATA boxes and patches adjacent to TSS (+1 nucleosome region); the resulting timescale distributions have non-exponential tails.

Genome-wide profiling of nucleosome sensitivity and chromatin accessibility in Drosophila melanogaster

Nucleosomal DNA is thought to be generally inaccessible to DNA-binding factors, such as micrococcal nuclease (MNase). Here, we digest Drosophila chromatin with high and low concentrations of MNase to reveal two distinct nucleosome types: MNase-sensitive and MNase-resistant. MNase-resistant nucleosomes assemble on sequences depleted of A/T and enriched in G/C-containing dinucleotides, whereas MNase-sensitive nucleosomes form on A/T-rich sequences found at transcription start and termination sites, enhancers and DNase I hypersensitive sites. Estimates of nucleosome formation energies indicate that MNase-sensitive nucleosomes tend to be less stable than MNase-resistant ones. Strikingly, a decrease in cell growth temperature of about 10°C makes MNase-sensitive nucleosomes less accessible, suggesting that observed variations in MNase sensitivity are related to either thermal fluctuations of chromatin fibers or the activity of enzymatic machinery. In the vicinity of active genes and DNase I hypersensitive sites nucleosomes are organized into periodic arrays, likely due to 'phasing' off potential barriers formed by DNA-bound factors or by nucleosomes anchored to their positions through external interactions. The latter idea is substantiated by our biophysical model of nucleosome positioning and energetics, which predicts that nucleosomes immediately downstream of transcription start sites are anchored and recapitulates nucleosome phasing at active genes significantly better than sequence-dependent models.

DNA Physical Properties and Nucleosome Positions Are Major Determinants of HIV-1 Integrase Selectivity

Retroviral integrases (INs) catalyse the integration of the reverse transcribed viral DNA into the host cell genome. This process is selective, and chromatin has been proposed to be a major factor regulating this step in the viral life cycle. However, the precise underlying mechanisms are still under investigation. We have developed a new in vitro integration assay using physiologically-relevant, reconstituted genomic acceptor chromatin and high-throughput determination of nucleosome positions and integration sites, in parallel. A quantitative analysis of the resulting data reveals a chromatin-dependent redistribution of the integration sites and establishes a link between integration sites and nucleosome positions. The co-activator LEDGF/p75 enhanced integration but did not modify the integration sites under these conditions. We also conducted an in cellulo genome-wide comparative study of nucleosome positions and human immunodeficiency virus type-1 (HIV-1) integration sites identified experimentally in vivo. These studies confirm a preferential integration in nucleosome-covered regions. Using a DNA mechanical energy model, we show that the physical properties of DNA probed by IN binding are important in determining IN selectivity. These novel in vitro and in vivo approaches confirm that IN has a preference for integration into a nucleosome, and suggest the existence of two levels of IN selectivity. The first depends on the physical properties of the target DNA and notably, the energy required to fit DNA into the IN catalytic pocket. The second depends on the DNA deformation associated with DNA wrapping around a nucleosome. Taken together, these results indicate that HIV-1 IN is a shape-readout DNA binding protein.

Global mapping of DNA conformational flexibility on Saccharomyces cerevisiae

In this study we provide the first comprehensive map of DNA conformational flexibility in Saccharomyces cerevisiae complete genome. Flexibility plays a key role in DNA supercoiling and DNA/protein binding, regulating DNA transcription, replication or repair. Specific interest in flexibility analysis concerns its relationship with human genome instability. Enrichment in flexible sequences has been detected in unstable regions of human genome defined fragile sites, where genes map and carry frequent deletions and rearrangements in cancer. Flexible sequences have been suggested to be the determinants of fragile gene proneness to breakage; however, their actual role and properties remain elusive. Our in silico analysis carried out genome-wide via the StabFlex algorithm, shows the conserved presence of highly flexible regions in budding yeast genome as well as in genomes of other Saccharomyces sensu stricto species. Flexibile peaks in S. cerevisiae identify 175 ORFs mapping on their 3'UTR, a region affecting mRNA translation, localization and stability. (TA)n repeats of different extension shape the central structure of peaks and co-localize with polyadenylation efficiency element (EE) signals. ORFs with flexible peaks share common features. Transcripts are characterized by decreased half-life: this is considered peculiar of genes involved in regulatory systems with high turnover; consistently, their function affects biological processes such as cell cycle regulation or stress response. Our findings support the functional importance of flexibility peaks, suggesting that the flexible sequence may be derived by an expansion of canonical TAYRTA polyadenylation efficiency element. The flexible (TA)n repeat amplification could be the outcome of an evolutionary neofunctionalization leading to a differential 3'-end processing and expression regulation in genes with peculiar function. Our study provides a new support to the functional role of flexibility in genomes and a strategy for its characterization inside human fragile sites.

Ubiquitous human 'master' origins of replication are encoded in the DNA sequence via a local enrichment in nucleosome excluding energy barriers

As the elementary building block of eukaryotic chromatin, the nucleosome is at the heart of the compromise between the necessity of compacting DNA in the cell nucleus and the required accessibility to regulatory proteins. The recent availability of genome-wide experimental maps of nucleosome positions for many different organisms and cell types has provided an unprecedented opportunity to elucidate to what extent the DNA sequence conditions the primary structure of chromatin and in turn participates in the chromatin-mediated regulation of nuclear functions, such as gene expression and DNA replication. In this study, we use in vivo and in vitro genome-wide nucleosome occupancy data together with the set of nucleosome-free regions (NFRs) predicted by a physical model of nucleosome formation based on sequence-dependent bending properties of the DNA double-helix, to investigate the role of intrinsic nucleosome occupancy in the regulation of the replication spatio-temporal programme in human. We focus our analysis on the so-called replication U/N-domains that were shown to cover about half of the human genome in the germline (skew-N domains) as well as in embryonic stem cells, somatic and HeLa cells (mean replication timing U-domains). The 'master' origins of replication (MaOris) that border these megabase-sized U/N-domains were found to be specified by a few hundred kb wide regions that are hyper-sensitive to DNase I cleavage, hypomethylated, and enriched in epigenetic marks involved in transcription regulation, the hallmarks of localized open chromatin structures. Here we show that replication U/N-domain borders that are conserved in all considered cell lines have an environment highly enriched in nucleosome-excluding-energy barriers, suggesting that these ubiquitous MaOris have been selected during evolution. In contrast, MaOris that are cell-type-specific are mainly regulated epigenetically and are no longer favoured by a local abundance of intrinsic NFRs encoded in the DNA sequence. At the smaller few hundred bp scale of gene promoters, CpG-rich promoters of housekeeping genes found nearby ubiquitous MaOris as well as CpG-poor promoters of tissue-specific genes found nearby cell-type-specific MaOris, both correspond to in vivo NFRs that are not coded as nucleosome-excluding-energy barriers. Whereas the former promoters are likely to correspond to high occupancy transcription factor binding regions, the latter are an illustration that gene regulation in human is typically cell-type-specific.

Intasome architecture and chromatin density modulate retroviral integration into nucleosome

Background

Retroviral integration depends on the interaction between intasomes, host chromatin and cellular targeting cofactors as LEDGF/p75 or BET proteins. Previous studies indicated that the retroviral integrase, by itself, may play a role in the local integration site selection within nucleosomal target DNA. We focused our study on this local association by analyzing the intrinsic properties of various retroviral intasomes to functionally accommodate different chromatin structures in the lack of other cofactors.

Results

Using in vitro conditions allowing the efficient catalysis of full site integration without these cofactors, we show that distinct retroviral integrases are not equally affected by chromatin compactness. Indeed, while PFV and MLV integration reactions are favored into dense and stable nucleosomes, HIV-1 and ASV concerted integration reactions are preferred into poorly dense chromatin regions of our nucleosomal acceptor templates. Predicted nucleosome occupancy around integration sites identified in infected cells suggests the presence of a nucleosome at the MLV and HIV-1 integration sites surrounded by differently dense chromatin. Further analyses of the relationships between the in vitro integration site selectivity and the structure of the inserted DNA indicate that structural constraints within intasomes could account for their ability to accommodate nucleosomal DNA and could dictate their capability to bind nucleosomes functionally in these specific chromatin contexts.

Conclusions

Thus, both intasome architecture and compactness of the chromatin surrounding the targeted nucleosome appear important determinants of the retroviral integration site selectivity. This supports a mechanism involving a global targeting of the intasomes toward suitable chromatin regions followed by a local integration site selection modulated by the intrinsic structural constraints of the intasomes governing the target DNA bending and dictating their sensitivity toward suitable specific nucleosomal structures and density.

Electronic supplementary material

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

Embryonic stem cell specific "master" replication origins at the heart of the loss of pluripotency

Epigenetic regulation of the replication program during mammalian cell differentiation remains poorly understood. We performed an integrative analysis of eleven genome-wide epigenetic profiles at 100 kb resolution of Mean Replication Timing (MRT) data in six human cell lines. Compared to the organization in four chromatin states shared by the five somatic cell lines, embryonic stem cell (ESC) line H1 displays (i) a gene-poor but highly dynamic chromatin state (EC4) associated to histone variant H2AZ rather than a HP1-associated heterochromatin state (C4) and (ii) a mid-S accessible chromatin state with bivalent gene marks instead of a polycomb-repressed heterochromatin state. Plastic MRT regions (≲ 20% of the genome) are predominantly localized at the borders of U-shaped timing domains. Whereas somatic-specific U-domain borders are gene-dense GC-rich regions, 31.6% of H1-specific U-domain borders are early EC4 regions enriched in pluripotency transcription factors NANOG and OCT4 despite being GC poor and gene deserts. Silencing of these ESC-specific “master” replication initiation zones during differentiation corresponds to a loss of H2AZ and an enrichment in H3K9me3 mark characteristic of late replicating C4 heterochromatin. These results shed a new light on the epigenetically regulated global chromatin reorganization that underlies the loss of pluripotency and lineage commitment.

During development, embryonic stem cell (ESC) enter a program of cell differentiation eventually leading to all the necessary differentiated cell types. Understanding the mechanisms responsible for the underlying modifications of the gene expression program is of fundamental importance, as it will likely have strong impact on the development of regenerative medicine. We show that besides some epigenetic regulation, ubiquitous master replication origins at replication timing U-domain borders shared by 6 human cell types are transcriptionally active open chromatin regions specified by a local enrichment in nucleosome free regions encoded in the DNA sequence suggesting that they have been selected during evolution. In contrast, ESC specific master replication origins bear a unique epigenetic signature (enrichment in CTCF, H2AZ, NANOG, OCT4, …) likely contributing to maintain ESC chromatin in a highly dynamic and accessible state that is refractory to polycomb and HP1 heterochromatin spreading. These ESC specific master origins thus appear as key genomic regions where epigenetic control of chromatin organization is at play to maintain pluripotency of stem cell lineages and to guide lineage commitment to somatic cell types.

Large replication skew domains delimit GC-poor gene deserts in human

Besides their large-scale organization in isochores, mammalian genomes display megabase-sized regions, spanning both genes and intergenes, where the strand nucleotide composition asymmetry decreases linearly, possibly due to replication activity. These so-called skew-N domains cover about a third of the human genome and are bordered by two skew upward jumps that were hypothesized to compose a subset of "master" replication origins active in the germline. Skew-N domains were shown to exhibit a particular gene organization. Genes with CpG-rich promoters likely expressed in the germline are over represented near the master replication origins, with large genes being co-oriented with replication fork progression, which suggests some coordination of replication and transcription. In this study, we describe another skew structure that covers ∼13% of the human genome and that is bordered by putative master replication origins similar to the ones flanking skew-N domains. These skew-split-N domains have a shape reminiscent of a N, but split in half, leaving in the center a region of null skew whose length increases with domain size. These central regions (median size ∼860 kb) have a homogeneous composition, i.e. both a null and constant skew and a constant and low GC content. They correspond to heterochromatin gene deserts found in low-GC isochores with an average gene density of 0.81 promoters/Mb as compared to 7.73 promoters/Mb genome wide. The analysis of epigenetic marks and replication timing data confirms that, in these late replicating heterochomatic regions, the initiation of replication is likely to be random. This contrasts with the transcriptionally active euchromatin state found around the bordering well positioned master replication origins. Altogether skew-N domains and skew-split-N domains cover about 50% of the human genome.

Quantifying the role of steric constraints in nucleosome positioning

Statistical positioning, the localization of nucleosomes packed against a fixed barrier, is conjectured to explain the array of well-positioned nucleosomes at the 5′ end of genes, but the extent and precise implications of statistical positioning in vivo are unclear. We examine this hypothesis quantitatively and generalize the idea to include moving barriers as well as nucleosomes actively packed against a barrier. Early experiments noted a similarity between the nucleosome profile aligned and averaged across genes and that predicted by statistical positioning; however, we demonstrate that aligning random nucleosomes also generates the same profile, calling the previous interpretation into question. New rigorous results reformulate statistical positioning as predictions on the variance structure of nucleosome locations in individual genes. In particular, a quantity termed the variance gradient, describing the change in variance between adjacent nucleosomes, is tested against recent high-throughput nucleosome sequencing data. Constant variance gradients provide support for generalized statistical positioning in ∼50% of long genes. Genes that deviate from predictions have high nucleosome turnover and cell-to-cell gene expression variability. The observed variance gradient suggests an effective nucleosome size of 158 bp, instead of the commonly perceived 147 bp. Our analyses thus clarify the role of statistical positioning in vivo.

Time to be versatile: regulation of the replication timing program in budding yeast

Eukaryotic replication origins are activated at different times during the S phase of the cell cycle, following a temporal program that is stably transmitted to daughter cells. Although the mechanisms that control initiation at the level of individual origins are now well understood, much less is known on how cells coordinate replication at hundreds of origins distributed on the chromosomes. In this review, we discuss recent advances shedding new light on how this complex process is regulated in the budding yeast Saccharomyces cerevisiae. The picture that emerges from these studies is that replication timing is regulated in cis by mechanisms modulating the chromatin structure and the subnuclear organization of origins. These mechanisms do not affect the licensing of replication origins but determine their ability to compete for limiting initiation factors, which are recycled from early to late origins throughout the length of the S phase.

Nucleosome positioning and kinetics near transcription-start-site barriers are controlled by interplay between active remodeling and DNA sequence

We investigate how DNA sequence, ATP-dependent chromatin remodeling and nucleosome-depleted ‘barriers’ co-operate to determine the kinetics of nucleosome organization, in a stochastic model of nucleosome positioning and dynamics. We find that ‘statistical’ positioning of nucleosomes against ‘barriers’, hypothesized to control chromatin structure near transcription start sites, requires active remodeling and therefore cannot be described using equilibrium statistical mechanics. We show that, unlike steady-state occupancy, DNA site exposure kinetics near a barrier is dominated by DNA sequence rather than by proximity to the barrier itself. The timescale for formation of positioning patterns near barriers is proportional to the timescale for active nucleosome eviction. We also show that there are strong gene-to-gene variations in nucleosome positioning near barriers, which are eliminated by averaging over many genes. Our results suggest that measurement of nucleosome kinetics can reveal information about sequence-dependent regulation that is not apparent in steady-state nucleosome occupancy.

Fitness landscape for nucleosome positioning

Histone-DNA complexes, so-called nucleosomes, are the building blocks of DNA packaging in eukaryotic cells. The histone-binding affinity of a local DNA segment depends on its elastic properties and determines its accessibility within the nucleus, which plays an important role in the regulation of gene expression. Here, we derive a fitness landscape for intergenic DNA segments in yeast as a function of two molecular phenotypes: their elasticity-dependent histone affinity and their coverage with transcription factor binding sites. This landscape reveals substantial selection against nucleosome formation over a wide range of both phenotypes. We use it as the core component of a quantitative evolutionary model for intergenic DNA segments. This model consistently predicts the observed diversity of histone affinities within wild Saccharomyces paradoxus populations, as well as the affinity divergence between neighboring Saccharomyces species. Our analysis establishes histone binding and transcription factor binding as two separable modes of sequence evolution, each of which is a direct target of natural selection.

Gene promoters show chromosome-specificity and reveal chromosome territories in humans

Background

Gene promoters have guided evolution processes for millions of years. It seems that they were the main engine responsible for the integration of different mutations favorable for the environmental conditions. In cooperation with different transcription factors and other biochemical components, these regulatory regions dictate the synthesis frequency of RNA molecules. Predominantly in the last decade, it has become clear that nuclear organization impacts upon gene regulation. To fully understand the connections between Homo sapiens chromosomes and their gene promoters, we analyzed 1200 promoter sequences using our Kappa Index of Coincidence method.

Results

In order to measure the structural similarity of gene promoters, we used two-dimensional image-based patterns obtained through Kappa Index of Coincidence (Kappa IC) and (C+G)% values. The center of weight of each promoter pattern indicated a structure similarity between promoters of each chromosome. Furthermore, the proximity of chromosomes seems to be in accordance to the structural similarity of their gene promoters. The arrangement of chromosomes according to Kappa IC values of promoters, shows a striking symmetry between the chromosome length and the structure of promoters located on them. High Kappa IC and (C+G)% values of gene promoters were also directly associated with the most frequent genetic diseases. Taking into consideration these observations, a general hypothesis for the evolutionary dynamics of the genome has been proposed. In this hypothesis, heterochromatin and euchromatin domains exchange DNA sequences according to a difference in the rate of Slipped Strand Mispairing and point mutations.

Conclusions

In this paper we showed that gene promoters appear to be specific to each chromosome. Furthermore, the proximity between chromosomes seems to be in accordance to the structural similarity of their gene promoters. Our findings are based on comprehensive data from Transcriptional Regulatory Element Database and a new computer model whose core is using Kappa index of coincidence.

iNuc-PhysChem: a sequence-based predictor for identifying nucleosomes via physicochemical properties

Nucleosome positioning has important roles in key cellular processes. Although intensive efforts have been made in this area, the rules defining nucleosome positioning is still elusive and debated. In this study, we carried out a systematic comparison among the profiles of twelve DNA physicochemical features between the nucleosomal and linker sequences in the Saccharomyces cerevisiae genome. We found that nucleosomal sequences have some position-specific physicochemical features, which can be used for in-depth studying nucleosomes. Meanwhile, a new predictor, called iNuc-PhysChem, was developed for identification of nucleosomal sequences by incorporating these physicochemical properties into a 1788-D (dimensional) feature vector, which was further reduced to a 884-D vector via the IFS (incremental feature selection) procedure to optimize the feature set. It was observed by a cross-validation test on a benchmark dataset that the overall success rate achieved by iNuc-PhysChem was over 96% in identifying nucleosomal or linker sequences. As a web-server, iNuc-PhysChem is freely accessible to the public at http://lin.uestc.edu.cn/server/iNuc-PhysChem. For the convenience of the vast majority of experimental scientists, a step-by-step guide is provided on how to use the web-server to get the desired results without the need to follow the complicated mathematics that were presented just for the integrity in developing the predictor. Meanwhile, for those who prefer to run predictions in their own computers, the predictor's code can be easily downloaded from the web-server. It is anticipated that iNuc-PhysChem may become a useful high throughput tool for both basic research and drug design.

Sequence-based prediction of single nucleosome positioning and genome-wide nucleosome occupancy

Nucleosome positioning dictates eukaryotic DNA compaction and access. To predict nucleosome positions in a statistical mechanics model, we exploited the knowledge that nucleosomes favor DNA sequences with specific periodically occurring dinucleotides. Our model is the first to capture both dyad position within a few base pairs, and free binding energy within 2 k(B)T, for all the known nucleosome positioning sequences. By applying Percus's equation to the derived energy landscape, we isolate sequence effects on genome-wide nucleosome occupancy from other factors that may influence nucleosome positioning. For both in vitro and in vivo systems, three parameters suffice to predict nucleosome occupancy with correlation coefficients of respectively 0.74 and 0.66. As predicted, we find the largest deviations in vivo around transcription start sites. This relatively simple algorithm can be used to guide future studies on the influence of DNA sequence on chromatin organization.

DNA structure, nucleosome placement and chromatin remodelling: a perspective

A major question in chromatin biology is to what extent the sequence of DNA directly determines the genetic and chromatin organization of a eukaryotic genome? We consider two aspects to this question: the DNA sequence-specified positioning of nucleosomes and the determination of NDRs (nucleosome-depleted regions) or barriers. We argue that, in budding yeast, while DNA sequence-specified nucleosome positioning may contribute to positions flanking the regions lacking nucleosomes, DNA thermodynamic stability is a major component determinant of the genetic organization of this organism.

Structural features based genome-wide characterization and prediction of nucleosome organization

Background

Nucleosome distribution along chromatin dictates genomic DNA accessibility and thus profoundly influences gene expression. However, the underlying mechanism of nucleosome formation remains elusive. Here, taking a structural perspective, we systematically explored nucleosome formation potential of genomic sequences and the effect on chromatin organization and gene expression in S. cerevisiae.

Results

We analyzed twelve structural features related to flexibility, curvature and energy of DNA sequences. The results showed that some structural features such as DNA denaturation, DNA-bending stiffness, Stacking energy, Z-DNA, Propeller twist and free energy, were highly correlated with in vitro and in vivo nucleosome occupancy. Specifically, they can be classified into two classes, one positively and the other negatively correlated with nucleosome occupancy. These two kinds of structural features facilitated nucleosome binding in centromere regions and repressed nucleosome formation in the promoter regions of protein-coding genes to mediate transcriptional regulation. Based on these analyses, we integrated all twelve structural features in a model to predict more accurately nucleosome occupancy in vivo than the existing methods that mainly depend on sequence compositional features. Furthermore, we developed a novel approach, named DLaNe, that located nucleosomes by detecting peaks of structural profiles, and built a meta predictor to integrate information from different structural features. As a comparison, we also constructed a hidden Markov model (HMM) to locate nucleosomes based on the profiles of these structural features. The result showed that the meta DLaNe and HMM-based method performed better than the existing methods, demonstrating the power of these structural features in predicting nucleosome positions.

Conclusions

Our analysis revealed that DNA structures significantly contribute to nucleosome organization and influence chromatin structure and gene expression regulation. The results indicated that our proposed methods are effective in predicting nucleosome occupancy and positions and that these structural features are highly predictive of nucleosome organization.

The implementation of our DLaNe method based on structural features is available online.

Distal chromatin structure influences local nucleosome positions and gene expression

The positions of nucleosomes across the genome influence several cellular processes, including gene transcription. However, our understanding of the factors dictating where nucleosomes are located and how this affects gene regulation is still limited. Here, we perform an extensive in vivo study to investigate the influence of the neighboring chromatin structure on local nucleosome positioning and gene expression. Using truncated versions of the Saccharomyces cerevisiae URA3 gene, we show that nucleosome positions in the URA3 promoter are at least partly determined by the local DNA sequence, with so-called ‘antinucleosomal elements’ like poly(dA:dT) tracts being key determinants of nucleosome positions. In addition, we show that changes in the nucleosome positions in the URA3 promoter strongly affect the promoter activity. Most interestingly, in addition to demonstrating the effect of the local DNA sequence, our study provides novel in vivo evidence that nucleosome positions are also affected by the position of neighboring nucleosomes. Nucleosome structure may therefore be an important selective force for conservation of gene order on a chromosome, because relocating a gene to another genomic position (where the positions of neighboring nucleosomes are different from the original locus) can have dramatic consequences for the gene's nucleosome structure and thus its expression.

A comparison study on feature selection of DNA structural properties for promoter prediction

Background

Promoter prediction is an integrant step for understanding gene regulation and annotating genomes. Traditional promoter analysis is mainly based on sequence compositional features. Recently, many kinds of structural features have been employed in promoter prediction. However, considering the high-dimensionality and overfitting problems, it is unfeasible to utilize all available features for promoter prediction. Thus it is necessary to choose some appropriate features for the prediction task.

Results

This paper conducts an extensive comparison study on feature selection of DNA structural properties for promoter prediction. Firstly, to examine whether promoters possess some special structures, we carry out a systematical comparison among the profiles of thirteen structural features on promoter and non-promoter sequences. Secondly, we investigate the correlations between these structural features and promoter sequences. Thirdly, both filter and wrapper methods are utilized to select appropriate feature subsets from thirteen different kinds of structural features for promoter prediction, and the predictive power of the selected feature subsets is evaluated. Finally, we compare the prediction performance of the feature subsets selected in this paper with nine existing promoter prediction approaches.

Conclusions

Experimental results show that the structural features are differentially correlated to promoters. Specifically, DNA-bending stiffness, DNA denaturation and energy-related features are highly correlated with promoters. The predictive power for promoter sequences differentiates greatly among different structural features. Selecting the relevant features can significantly improve the accuracy of promoter prediction.

Revisiting polymer statistical physics to account for the presence of long-range-correlated structural disorder in 2D DNA chains

We elaborate on a generalization of the 2D wormlike chain (WLC) model that accounts for the presence of long-range correlations (LRC) in the intrinsic curvature distribution of eukaryotic DNA. This model predicts some decrease of the DNA persistence length resulting from some large-scale intrinsic curvature induced by sequence-dependent persistent random distribution of local bending sites. When assisting exact analytical calculations by numerical DNA simulations, we show that the conjugated contributions of i) the thermal curvature fluctuations characterized by the "dynamic" persistence length ℓ(p)(d) = 2A, where A is the elastic bending modulus, and ii) the intrinsic LRC curvature disorder of amplitude σ(o) and Hurst exponent H > 1/2, characterized by a "static" persistence length ℓ(p)(H) = A(1/2H)σ(o)(-1/H) Γ(1/2H + 1), can be described by a continuum of generalized WLC (GWLC) models parametrized by the LRC exponent H. We use perturbation analysis to investigate the two limiting cases of weak static disorder (w(H) << 1 and weak dynamical fluctuations (1/w (H) << 1), where w(H) = l(p)(d)/l(p)(H) is a dimensionless parameter. From a quantitative point of view, our study demonstrates that even for a small value of the LRC (H approximately equal 0.6-0.8) static disorder amplitude σ(o) ~ 10(-2), as previously reported for genomic DNA, the decrease of the persistence length from the WLC prediction l(p)(d) can be very significant, up to twofold. The implications of these results on the first steps of compaction of DNA in eukaryotic cells are discussed.

HIV-1 integration into chromatin: new insights and future perspectives

Human AIDS, owing to HIV-1, remains a major health issue despite the numerous antiviral compounds that are now available. One of the main obstacles for these antiviral strategies is the appearance of highly resistant viral strains, making the search for new antiviral targets a crucial issue. The insertion of the viral genome into the cell chromosomal DNA, catalyzed by the retroviral integrase enzyme, is required for stable infection and thus constitutes a suitable target for the antiretroviral molecules currently used in therapy. The precise mechanism of interaction between the integration complex (intasome) and the target chromatin remains poorly understood, thus restricting the development of new therapeutic strategies and original methods to control gene transfer in gene therapy. In this article, we summarize the data obtained in this field and underline recent results highlighting this process and paving the way for new and more detailed understanding of this mechanism.

Nucleosome positioning in Saccharomyces cerevisiae

The DNA of eukaryotic cells is spooled around large histone protein complexes, forming nucleosomes that make up the basis for a high-order packaging structure called chromatin. Compared to naked DNA, nucleosomal DNA is less accessible to regulatory proteins and regulatory processes. The exact positions of nucleosomes therefore influence several cellular processes, including gene expression, chromosome segregation, recombination, replication, and DNA repair. Here, we review recent technological advances enabling the genome-wide mapping of nucleosome positions in the model eukaryote Saccharomyces cerevisiae. We discuss the various parameters that determine nucleosome positioning in vivo, including cis factors like AT content, variable tandem repeats, and poly(dA:dT) tracts that function as chromatin barriers and trans factors such as chromatin remodeling complexes, transcription factors, histone-modifying enzymes, and RNA polymerases. In the last section, we review the biological role of chromatin in gene transcription, the evolution of gene regulation, and epigenetic phenomena.

Mechanics of the IL2RA gene activation revealed by modeling and atomic force microscopy

Transcription implies recruitment of RNA polymerase II and transcription factors (TFs) by DNA melting near transcription start site (TSS). Combining atomic force microscopy and computer modeling, we investigate the structural and dynamical properties of the IL2RA promoter and identify an intrinsically negative supercoil in the PRRII region (containing Elf-1 and HMGA1 binding sites), located upstream of a curved DNA region encompassing TSS. Conformational changes, evidenced by time-lapse studies, result in the progressive positioning of curvature apex towards the TSS, likely facilitating local DNA melting. In vitro assays confirm specific binding of the General Transcription Factors (GTFs) TBP and TFIIB over TATA-TSS position, where an inhibitory nucleosome prevented preinitiation complex (PIC) formation and uncontrolled DNA melting. These findings represent a substantial advance showing, first, that the structural properties of the IL2RA promoter are encoded in the DNA sequence and second, that during the initiation process DNA conformation is dynamic and not static.

Functional coupling between HIV-1 integrase and the SWI/SNF chromatin remodeling complex for efficient in vitro integration into stable nucleosomes

Establishment of stable HIV-1 infection requires the efficient integration of the retroviral genome into the host DNA. The molecular mechanism underlying the control of this process by the chromatin structure has not yet been elucidated. We show here that stably associated nucleosomes strongly inhibit in vitro two viral-end integration by decreasing the accessibility of DNA to integrase. Remodeling of the chromatinized template by the SWI/SNF complex, whose INI1 major component interacts with IN, restores and redirects the full-site integration into the stable nucleosome region. These effects are not observed after remodeling by other human remodeling factors such as SNF2H or BRG1 lacking the integrase binding protein INI1. This suggests that the restoration process depends on the direct interaction between IN and the whole SWI/SNF complex, supporting a functional coupling between the remodeling and integration complexes. Furthermore, in silico comparison between more than 40,000 non-redundant cellular integration sites selected from literature and nucleosome occupancy predictions also supports that HIV-1 integration is promoted in the genomic region of weaker intrinsic nucleosome density in the infected cell. Our data indicate that some chromatin structures can be refractory for integration and that coupling between nucleosome remodeling and HIV-1 integration is required to overcome this natural barrier.

A lattice model for transcription factor access to nucleosomal DNA

Nucleosomes, the basic repeating unit of chromatin, consist of 147 basepairs of DNA that are wrapped in almost two turns around a histone protein octamer core. Because ∼3/4 of the human genomic DNA is found within nucleosomes, their position and DNA interaction is an essential determinant for the DNA access of gene-specific transcription factors and other proteins. Here, a DNA lattice model was developed for describing ligand binding in the presence of a nucleosome. The model takes into account intermediate states, in which DNA is partially unwrapped from the histone octamer. This facilitates access of transcription factors to up to 60 DNA basepairs located in the outer turn of nucleosomal DNA, while the inner DNA turn was found to be more resistant to competitive ligand binding. As deduced from quantitative comparisons with recently published experimental data, our model provides a better description than the previously used all-or-none lattice-binding model. Importantly, nucleosome-occupancy maps predicted by the nucleosome-unwrapping model also differed significantly when partial unwrapping of nucleosomal DNA was considered. In addition, large effects on the cooperative binding of transcription factors to multiple binding sites occluded by the nucleosome were apparent. These findings indicate that partial unwrapping of DNA from the histone octamer needs to be taken into account in quantitative models of gene regulation in chromatin.

Contribution of histone sequence preferences to nucleosome organization: proposed definitions and methodology

We propose definitions and procedures for comparing nucleosome maps and discuss current agreement and disagreement on the effect of histone sequence preferences on nucleosome organization in vivo.

Structural constraints revealed in consistent nucleosome positions in the genome of S. cerevisiae

Background

Recent advances in the field of high-throughput genomics have rendered possible the performance of genome-scale studies to define the nucleosomal landscapes of eukaryote genomes. Such analyses are aimed towards providing a better understanding of the process of nucleosome positioning, for which several models have been suggested. Nevertheless, questions regarding the sequence constraints of nucleosomal DNA and how they may have been shaped through evolution remain open. In this paper, we analyze in detail different experimental nucleosome datasets with the aim of providing a hypothesis for the emergence of nucleosome-forming sequences.

Results

We compared the complete sets of nucleosome positions for the budding yeast (Saccharomyces cerevisiae) as defined in the output of two independent experiments with the use of two different experimental techniques. We found that < 10% of the experimentally defined nucleosome positions were consistently positioned in both datasets. This subset of well-positioned nucleosomes, when compared with the bulk, was shown to have particular properties at both sequence and structural levels. Consistently positioned nucleosomes were also shown to occur preferentially in pairs of dinucleosomes, and to be surprisingly less conserved compared with their adjacent nucleosome-free linkers.

Conclusion

Our findings may be combined into a hypothesis for the emergence of a weak nucleosome-positioning code. According to this hypothesis, consistent nucleosomes may be partly guided by nearby nucleosome-free regions through statistical positioning. Once established, a set of well-positioned consistent nucleosomes may impose secondary constraints that further shape the structure of the underlying DNA. We were able to capture these constraints through the application of a recently introduced structural property that is related to the symmetry of DNA curvature. Furthermore, we found that both consistently positioned nucleosomes and their adjacent nucleosome-free regions show an increased tendency towards conservation of this structural feature.

Statistical-mechanical lattice models for protein-DNA binding in chromatin

Statistical-mechanical lattice models for protein-DNA binding are well established as a method to describe complex ligand binding equilibria measured in vitro with purified DNA and protein components. Recently, a new field of applications has opened up for this approach since it has become possible to experimentally quantify genome-wide protein occupancies in relation to the DNA sequence. In particular, the organization of the eukaryotic genome by histone proteins into a nucleoprotein complex termed chromatin has been recognized as a key parameter that controls the access of transcription factors to the DNA sequence. New approaches have to be developed to derive statistical-mechanical lattice descriptions of chromatin-associated protein-DNA interactions. Here, we present the theoretical framework for lattice models of histone-DNA interactions in chromatin and investigate the (competitive) DNA binding of other chromosomal proteins and transcription factors. The results have a number of applications for quantitative models for the regulation of gene expression.

Prediction of nucleosome DNA formation potential and nucleosome positioning using increment of diversity combined with quadratic discriminant analysis

In this work, a novel method was developed to distinguish nucleosome DNA and linker DNA based on increment of diversity combined with quadratic discriminant analysis (IDQD), using k-mer frequency of nucleotides in genome. When used to predict DNA potential for forming nucleosomes, the model achieved a high accuracy of 94.94%, 77.60%, and 86.81%, respectively, for Saccharomyces cerevisiae, Homo sapiens, and Drosophila melanogaster. The area under the receiver operator characteristics curve of our classifier was 0.982 for S. cerevisiae. Our results indicate that DNA sequence preference is critical for nucleosome formation potential and is likely conserved across eukaryotes. The model successfully identified nucleosome-enriched or nucleosome-depleted regions in S. cerevisiae genome, suggesting nucleosome positioning depends on DNA sequence preference. Thus, IDQD classifier is useful for predicting nucleosome positioning.

Quantitative test of the barrier nucleosome model for statistical positioning of nucleosomes up- and downstream of transcription start sites

The positions of nucleosomes in eukaryotic genomes determine which parts of the DNA sequence are readily accessible for regulatory proteins and which are not. Genome-wide maps of nucleosome positions have revealed a salient pattern around transcription start sites, involving a nucleosome-free region (NFR) flanked by a pronounced periodic pattern in the average nucleosome density. While the periodic pattern clearly reflects well-positioned nucleosomes, the positioning mechanism is less clear. A recent experimental study by Mavrich et al. argued that the pattern observed in Saccharomyces cerevisiae is qualitatively consistent with a “barrier nucleosome model,” in which the oscillatory pattern is created by the statistical positioning mechanism of Kornberg and Stryer. On the other hand, there is clear evidence for intrinsic sequence preferences of nucleosomes, and it is unclear to what extent these sequence preferences affect the observed pattern. To test the barrier nucleosome model, we quantitatively analyze yeast nucleosome positioning data both up- and downstream from NFRs. Our analysis is based on the Tonks model of statistical physics which quantifies the interplay between the excluded-volume interaction of nucleosomes and their positional entropy. We find that although the typical patterns on the two sides of the NFR are different, they are both quantitatively described by the same physical model with the same parameters, but different boundary conditions. The inferred boundary conditions suggest that the first nucleosome downstream from the NFR (the +1 nucleosome) is typically directly positioned while the first nucleosome upstream is statistically positioned via a nucleosome-repelling DNA region. These boundary conditions, which can be locally encoded into the genome sequence, significantly shape the statistical distribution of nucleosomes over a range of up to ∼1,000 bp to each side.

Within the last five years, knowledge about nucleosome organization on the genome has grown dramatically. To a large extent, this has been achieved by an increasing number of experimental studies determining nucleosome positions at high resolution over entire genomes. Particular attention has been paid to promoter regions, where a canonical pattern has been established: a nucleosome-free region with pronounced adjacent oscillations in the nucleosome density. Here we tested to what extent this pattern may be quantitatively described by a minimal physical model, a one-dimensional gas of impenetrable particles, commonly referred to as the “Tonks gas.” In this model, density oscillations occur close to a boundary at dense packing. Our systematic quantitative analysis reveals that, in an average over many promoters, a Tonks gas model can indeed account for the nucleosome organization to both sides of the nucleosome-free region, if one allows for different boundary conditions at the two edges. On the downstream side, a single nucleosome is typically directly positioned such that it forms an obstacle for the neighboring nucleosomes, while such a barrier nucleosome is typically missing on the upstream side.

The DNA sequence-dependence of nucleosome positioning in vivo and in vitro

The contribution of histone-DNA interactions to nucleosome positioning in vivo is currently a matter of debate. We argue here that certain nucleosome positions, often in promoter regions, in yeast may be, at least in part, specified by the DNA sequence. In contrast other positions may be poorly specified. Positioning thus has both statistical and DNA-determined components. We further argue that the relative affinity of the octamer for different DNA sequences can vary and therefore the interaction of histones with the DNA is a 'tunable' property.

Prestressed nuclear organization in living cells

The nucleus is maintained in a prestressed state within eukaryotic cells, stabilized mechanically by chromatin structure and other nuclear components on its inside, and cytoskeletal components on its outside. Nuclear architecture is emerging to be critical to the governance of chromatin assembly, regulation of genome function and cellular homeostasis. Elucidating the prestressed organization of the nucleus is thus important to understand how the nuclear architecture impinges on its function. In this chapter, various chemical and mechanical methods have been described to probe the prestressed organization of the nucleus.