Evolutionary divergence of intrinsic and trans-regulated nucleosome positioning sequences reveals plastic rules for chromatin organization

Determinants of nucleosome positioning

Nucleosome positioning is critical for gene expression and most DNA-related processes. Here we review the dominant patterns of nucleosome positioning that have been observed and summarize the current understanding of their underlying determinants. The genome-wide pattern of nucleosome positioning is determined by the combination of DNA sequence, ATP-dependent nucleosome remodeling enzymes and transcription factors that include activators, components of the preinitiation complex and elongating RNA polymerase II. These determinants influence each other such that the resulting nucleosome positioning patterns are likely to differ among genes and among cells in a population, with consequent effects on gene expression.

Back to the origin: reconsidering replication, transcription, epigenetics, and cell cycle control

In bacteria, replication is a carefully orchestrated event that unfolds the same way for each bacterium and each cell division. The process of DNA replication in bacteria optimizes cell growth and coordinates high levels of simultaneous replication and transcription. In metazoans, the organization of replication is more enigmatic. The lack of a specific sequence that defines origins of replication has, until recently, severely limited our ability to define the organizing principles of DNA replication. This question is of particular importance as emerging data suggest that replication stress is an important contributor to inherited genetic damage and the genomic instability in tumors. We consider here the replication program in several different organisms including recent genome-wide analyses of replication origins in humans. We review recent studies on the role of cytosine methylation in replication origins, the role of transcriptional looping and gene gating in DNA replication, and the role of chromatin's 3-dimensional structure in DNA replication. We use these new findings to consider several questions surrounding DNA replication in metazoans: How are origins selected? What is the relationship between replication and transcription? How do checkpoints inhibit origin firing? Why are there early and late firing origins? We then discuss whether oncogenes promote cancer through a role in DNA replication and whether errors in DNA replication are important contributors to the genomic alterations and gene fusion events observed in cancer. We conclude with some important areas for future experimentation.

Genetic and epigenetic determinants of DNA replication origins, position and activation

In the genome of eukaryotic cells, DNA synthesis is initiated at multiple sites called origins of DNA replication. Origins must fire only once per cell cycle and how this is achieved is now well understood. However, little is known about the mechanisms that determine when and where replication initiates in a given cell. A large body of evidence indicates that origins are not equal in terms of efficiency and timing of activation. Origin usage also changes concomitantly with the different cell differentiation programs. As DNA replication occurs in the context of chromatin, initiation could be influenced by multiple parameters, such as nucleosome positioning, histone modifications, and three-dimensional (3D) organization of the nucleus. This view is supported by recent genome-wide studies showing that DNA replication profiles are shaped by genetic and epigenetic processes that act both at the local and global levels to regulate origin function in eukaryotic cells.

Ubiquitous heterogeneity and asymmetry of the chromatin environment at regulatory elements

Gene regulation at functional elements (e.g., enhancers, promoters, insulators) is governed by an interplay of nucleosome remodeling, histone modifications, and transcription factor binding. To enhance our understanding of gene regulation, the ENCODE Consortium has generated a wealth of ChIP-seq data on DNA-binding proteins and histone modifications. We additionally generated nucleosome positioning data on two cell lines, K562 and GM12878, by MNase digestion and high-depth sequencing. Here we relate 14 chromatin signals (12 histone marks, DNase, and nucleosome positioning) to the binding sites of 119 DNA-binding proteins across a large number of cell lines. We developed a new method for unsupervised pattern discovery, the Clustered AGgregation Tool (CAGT), which accounts for the inherent heterogeneity in signal magnitude, shape, and implicit strand orientation of chromatin marks. We applied CAGT on a total of 5084 data set pairs to obtain an exhaustive catalog of high-resolution patterns of histone modifications and nucleosome positioning signals around bound transcription factors. Our analyses reveal extensive heterogeneity in how histone modifications are deposited, and how nucleosomes are positioned around binding sites. With the exception of the CTCF/cohesin complex, asymmetry of nucleosome positioning is predominant. Asymmetry of histone modifications is also widespread, for all types of chromatin marks examined, including promoter, enhancer, elongation, and repressive marks. The fine-resolution signal shapes discovered by CAGT unveiled novel correlation patterns between chromatin marks, nucleosome positioning, and sequence content. Meta-analyses of the signal profiles revealed a common vocabulary of chromatin signals shared across multiple cell lines and binding proteins.

CHD1 remodelers regulate nucleosome spacing in vitro and align nucleosomal arrays over gene coding regions in S. pombe

Nucleosome positioning governs access to eukaryotic genomes. Many genes show a stereotypic organisation at their 5'end: a nucleosome free region just upstream of the transcription start site (TSS) followed by a regular nucleosomal array over the coding region. The determinants for this pattern are unclear, but nucleosome remodelers are likely critical. Here we study the role of remodelers in global nucleosome positioning in S. pombe and the corresponding changes in expression. We find a striking evolutionary shift in remodeler usage between budding and fission yeast. The S. pombe RSC complex does not seem to be involved in nucleosome positioning, despite its prominent role in S. cerevisiae. While S. pombe lacks ISWI-type remodelers, it has two CHD1-type ATPases, Hrp1 and Hrp3. We demonstrate nucleosome spacing activity for Hrp1 and Hrp3 in vitro, and that together they are essential for linking regular genic arrays to most TSSs in vivo. Impaired arrays in the absence of either or both remodelers may lead to increased cryptic antisense transcription, but overall gene expression levels are only mildly affected.

CpG islands and GC content dictate nucleosome depletion in a transcription-independent manner at mammalian promoters

One clear hallmark of mammalian promoters is the presence of CpG islands (CGIs) at more than two-thirds of genes, whereas TATA boxes are only present at a minority of promoters. Using genome-wide approaches, we show that GC content and CGIs are major promoter elements in mammalian cells, able to govern open chromatin conformation and support paused transcription. First, we define three classes of promoters with distinct transcriptional directionality and pausing properties that correlate with their GC content. We further analyze the direct influence of GC content on nucleosome positioning and depletion and show that CpG content and CGI width correlate with nucleosome depletion both in vivo and in vitro. We also show that transcription is not essential for nucleosome exclusion but influences both a weak +1 and a well-positioned nucleosome at CGI borders. Altogether our data support the idea that CGIs have become an essential feature of promoter structure defining novel regulatory properties in mammals.

Nucleosome Positioning

Nucleosome positioning is not only related to genomic DNA compaction but also to other biological functions. After the chromatin is digested by micrococcal nuclease, nucleosomal (nucleosome-bound) DNA fragments can be sequenced and mapped on the genomic DNA sequence. Due to the development of modern DNA sequencing technology, genome-wide nucleosome mapping has been performed in a wide range of eukaryotic species. Comparative analyses of the nucleosome positions have revealed that the nucleosome is more frequently formed in exonic than intronic regions, and that most of transcription start and translation (or transcription) end sites are located in nucleosome linker DNA regions, indicating that nucleosome positioning influences transcription initiation, transcription termination, and gene splicing. In addition, nucleosomal DNA contains guanine and cytosine (G + C)-rich sequences and a high level of cytosine methylation. Thus, the nucleosome positioning system has been conserved during eukaryotic evolution.

A functional evolutionary approach to identify determinants of nucleosome positioning: a unifying model for establishing the genome-wide pattern

Although the genomic pattern of nucleosome positioning is broadly conserved, quantitative aspects vary over evolutionary timescales. We identify the cis and trans determinants of nucleosome positioning using a functional evolutionary approach involving S. cerevisiae strains containing large genomic regions from other yeast species. In a foreign species, nucleosome depletion at promoters is maintained over poly(dA:dT) tracts, whereas internucleosome spacing and all other aspects of nucleosome positioning tested are not. Interestingly, the locations of the +1 nucleosome and RNA start sites shift in concert. Strikingly, in a foreign species, nucleosome-depleted regions occur fortuitously in coding regions, and they often act as promoters that are associated with a positioned nucleosome array linked to the length of the transcription unit. We suggest a three-step model in which nucleosome remodelers, general transcription factors, and the transcriptional elongation machinery are primarily involved in generating the nucleosome positioning pattern in vivo.

Active nucleosome positioning beyond intrinsic biophysics is revealed by in vitro reconstitution

Genome-wide nucleosome maps revealed well-positioned nucleosomes as a major theme in eukaryotic genome organization. Promoter regions often show a conserved pattern with an NDR (nucleosome-depleted region) from which regular nucleosomal arrays emanate. Three mechanistic contributions to such NDR-array-organization and nucleosome positioning in general are discussed: DNA sequence, DNA binders and DNA-templated processes. Especially, intrinsic biophysics of DNA sequence preferences for nucleosome formation was prominently suggested to explain the majority of nucleosome positions ('genomic code for nucleosome positioning'). Nonetheless, non-histone factors that bind DNA with high or low specificity, such as transcription factors or remodelling enzymes respectively and processes such as replication, transcription and the so-called 'statistical positioning' may be involved too. Recently, these models were tested for yeast by genome-wide reconstitution. DNA sequence preferences as probed by SGD (salt gradient dialysis) reconstitution generated many NDRs, but only few individual nucleosomes, at their proper positions, and no arrays. Addition of a yeast extract and ATP led to dramatically more in vivo-like nucleosome positioning, including regular arrays for the first time. This improvement depended essentially on the extract and ATP but not on transcription or replication. Nucleosome occupancy and close spacing were maintained around promoters, even at lower histone density, arguing for active packing of nucleosomes against the 5' ends of genes rather than statistical positioning. A first extract fractionation identified a direct, specific, necessary, but not sufficient role for the RSC (remodels the structure of chromatin) remodelling enzyme. Collectively, nucleosome positioning in yeast is actively determined by factors beyond intrinsic biophysics, and in steady-state rather than at equilibrium.

Comparative studies of gene expression and the evolution of gene regulation

The hypothesis that differences in gene regulation have an important role in speciation and adaptation is more than 40 years old. With the advent of new sequencing technologies, we are able to characterize and study gene expression levels and associated regulatory mechanisms in a large number of individuals and species at an unprecedented resolution and scale. We have thus gained new insights into the evolutionary pressures that shape gene expression levels and have developed an appreciation for the relative importance of evolutionary changes in different regulatory genetic and epigenetic mechanisms. The current challenge is to link gene regulatory changes to adaptive evolution of complex phenotypes. Here we mainly focus on comparative studies in primates and how they are complemented by studies in model organisms.

Genome-wide identification and characterization of replication origins by deep sequencing

Background

DNA replication initiates at distinct origins in eukaryotic genomes, but the genomic features that define these sites are not well understood.

Results

We have taken a combined experimental and bioinformatic approach to identify and characterize origins of replication in three distantly related fission yeasts: Schizosaccharomyces pombe, Schizosaccharomyces octosporus and Schizosaccharomyces japonicus. Using single-molecule deep sequencing to construct amplification-free high-resolution replication profiles, we located origins and identified sequence motifs that predict origin function. We then mapped nucleosome occupancy by deep sequencing of mononucleosomal DNA from the corresponding species, finding that origins tend to occupy nucleosome-depleted regions.

Conclusions

The sequences that specify origins are evolutionarily plastic, with low complexity nucleosome-excluding sequences functioning in S. pombe and S. octosporus, and binding sites for trans-acting nucleosome-excluding proteins functioning in S. japonicus. Furthermore, chromosome-scale variation in replication timing is conserved independently of origin location and via a mechanism distinct from known heterochromatic effects on origin function. These results are consistent with a model in which origins are simply the nucleosome-depleted regions of the genome with the highest affinity for the origin recognition complex. This approach provides a general strategy for understanding the mechanisms that define DNA replication origins in eukaryotes.

Nucleosome positioning: bringing order to the eukaryotic genome

Nucleosomes are an essential component of eukaryotic chromosomes. The impact of nucleosomes is seen not just on processes that directly access the genome, such as transcription, but also on an evolutionary timescale. Recent studies in various organisms have provided high-resolution maps of nucleosomes throughout the genome. Computational analysis, in conjunction with many other kinds of data, has shed light on several aspects of nucleosome biology. Nucleosomes are positioned by several means, including intrinsic sequence biases, by stacking against a fixed barrier, by DNA-binding proteins and by chromatin remodelers. These studies underscore the important organizational role of nucleosomes in all eukaryotic genomes. This paper reviews recent genomic studies that have shed light on the determinants of nucleosome positioning and their impact on the genome.