A map of nucleosome positions in yeast at base-pair resolution.

The exact positions of nucleosomes along genomic DNA can influence many aspects of chromosome function. However, existing methods for mapping nucleosomes do not provide the necessary single-base-pair accuracy to determine these positions. Here we develop and apply a new approach for direct mapping of nucleosome centres on the basis of chemical modification of engineered histones. The resulting map locates nucleosome positions genome-wide in unprecedented detail and accuracy. It shows new aspects of the in vivo nucleosome organization that are linked to transcription factor binding, RNA polymerase pausing and the higher-order structure of the chromatin fibre.

Single-cell nucleosome mapping reveals the molecular basis of gene expression heterogeneity.

Nucleosomes, the basic unit of chromatin, have a critical role in the control of gene expression. Nucleosome positions have generally been determined by examining bulk populations of cells and then correlated with overall gene expression. Here, we describe a technique to determine nucleosome positioning in single cells by virtue of the ability of the nucleosome to protect DNA from GpC methylation. In the acid phosphatase inducible PHO5 gene, we find that there is significant cell-to-cell variation in nucleosome positions and shifts in nucleosome positioning correlate with changes in gene expression. However, nucleosome positioning is not absolute, and even with major shifts in gene expression, some cells fail to change nucleosome configuration. Mutations of the PHO5 promoter that introduce a poly(dA:dT) tract-stimulated gene expression under nonpermissive conditions led to shifts of positioned nucleosomes similar to induction of PHO5. By contrast, mutations that altered AA/TT/AT periodicity reduced gene expression upon PHO5 induction and stabilized nucleosomes in most cells, suggesting that enhanced nucleosome affinity for DNA antagonizes chromatin remodelers. Finally, we determined nucleosome positioning in two regions described as "fuzzy" or nucleosome-free when examined in a bulk assay. These regions consisted of distinct nucleosomes with a larger footprint for potential location and an increase population of cells lacking a nucleosome altogether. These data indicate an underlying complexity of nucleosome positioning that may contribute to the flexibility and heterogeneity of gene expression.

Chemical map of Schizosaccharomyces pombe reveals species-specific features in nucleosome positioning.

Using a recently developed chemical approach, we have generated a genome-wide map of nucleosomes in vivo in Schizosaccharomyces pombe (S. pombe) at base pair resolution. The shorter linker length previously identified in S. pombe is due to a preponderance of nucleosomes separated by ∼4/5 bp, placing nucleosomes on opposite faces of the DNA. The periodic dinucleotide feature thought to position nucleosomes is equally strong in exons as in introns, demonstrating that nucleosome positioning information can be superimposed on coding information. Unlike the case in Saccharomyces cerevisiae, A/T-rich sequences are enriched in S. pombe nucleosomes, particularly at ±20 bp around the dyad. This difference in nucleosome binding preference gives rise to a major distinction downstream of the transcription start site, where nucleosome phasing is highly predictable by A/T frequency in S. pombe but not in S. cerevisiae, suggesting that the genomes and DNA binding preferences of nucleosomes have coevolved in different species. The poly (dA-dT) tracts affect but do not deplete nucleosomes in S. pombe, and they prefer special rotational positions within the nucleosome, with longer tracts enriched in the 10- to 30-bp region from the dyad. S. pombe does not have a well-defined nucleosome-depleted region immediately upstream of most transcription start sites; instead, the -1 nucleosome is positioned with the expected spacing relative to the +1 nucleosome, and its occupancy is negatively correlated with gene expression. Although there is generally very good agreement between nucleosome maps generated by chemical cleavage and micrococcal nuclease digestion, the chemical map shows consistently higher nucleosome occupancy on DNA with high A/T content.

In vitro selection of DNAs with an increased propensity to form small circles.

A protocol was devised to select for DNA molecules that efficiently form circles from a library of 126 base pair DNAs containing 90 randomized base pairs. After six rounds of selection, individual molecules from the library showed 20- to 100-fold greater j-factors compared with the starting library, validating the selection protocol. High-throughput sequencing revealed a sinusoidal pattern of enrichment and de-enrichment of A/T dinucleotides in the random region with a 10.4 base pair period associated with the helicity of DNA. A similar, but more moderate pattern of C/G dinucleotides was offset by precisely half a helical turn. While C/G dinucleotide enrichments were evenly distributed, A/T dinucleotide enrichments displayed a preference to cluster in individual DNA molecules. The most highly enriched 10 base pair sequences in the random region contained adjacent blocks of A/T and C/G trinucleotides present in some, but not all, rapidly cyclizing molecules. The phased dinucleotide enrichments closely match those present in accurately mapped yeast nucleosomes, confirming the importance of DNA bending in nucleosome formation. However, at certain sites the nucleosomal DNAs show dinucleotide enrichments that differ substantially from the cyclization data. These discrepancies can often be correlated with sequence specific contacts that form between histones and DNA.

From DNA sequence to transcriptional behaviour: a quantitative approach.

Complex transcriptional behaviours are encoded in the DNA sequences of gene regulatory regions. Advances in our understanding of these behaviours have been recently gained through quantitative models that describe how molecules such as transcription factors and nucleosomes interact with genomic sequences. An emerging view is that every regulatory sequence is associated with a unique binding affinity landscape for each molecule and, consequently, with a unique set of molecule-binding configurations and transcriptional outputs. We present a quantitative framework based on existing methods that unifies these ideas. This framework explains many experimental observations regarding the binding patterns of factors and nucleosomes and the dynamics of transcriptional activation. It can also be used to model more complex phenomena such as transcriptional noise and the evolution of transcriptional regulation.

The DNA-encoded nucleosome organization of a eukaryotic genome.

Nucleosome organization is critical for gene regulation. In living cells this organization is determined by multiple factors, including the action of chromatin remodellers, competition with site-specific DNA-binding proteins, and the DNA sequence preferences of the nucleosomes themselves. However, it has been difficult to estimate the relative importance of each of these mechanisms in vivo, because in vivo nucleosome maps reflect the combined action of all influencing factors. Here we determine the importance of nucleosome DNA sequence preferences experimentally by measuring the genome-wide occupancy of nucleosomes assembled on purified yeast genomic DNA. The resulting map, in which nucleosome occupancy is governed only by the intrinsic sequence preferences of nucleosomes, is similar to in vivo nucleosome maps generated in three different growth conditions. In vitro, nucleosome depletion is evident at many transcription factor binding sites and around gene start and end sites, indicating that nucleosome depletion at these sites in vivo is partly encoded in the genome. We confirm these results with a micrococcal nuclease-independent experiment that measures the relative affinity of nucleosomes for approximately 40,000 double-stranded 150-base-pair oligonucleotides. Using our in vitro data, we devise a computational model of nucleosome sequence preferences that is significantly correlated with in vivo nucleosome occupancy in Caenorhabditis elegans. Our results indicate that the intrinsic DNA sequence preferences of nucleosomes have a central role in determining the organization of nucleosomes in vivo.

High-resolution nucleosome mapping of targeted regions using BAC-based enrichment.

We report a target enrichment method to map nucleosomes of large genomes at unprecedented coverage and resolution by deeply sequencing locus-specific mononucleosomal DNA enriched via hybridization with bacterial artificial chromosomes. We achieved ≈ 10 000-fold enrichment of specific loci, which enabled sequencing nucleosomes at up to ≈ 500-fold higher coverage than has been reported in a mammalian genome. We demonstrate the advantages of generating high-sequencing coverage for mapping the center of discrete nucleosomes, and we show the use of the method by mapping nucleosomes during T cell differentiation using nuclei from effector T-cells differentiated from clonal, isogenic, naïve, primary murine CD4 and CD8 T lymphocytes. The analysis reveals that discrete nucleosomes exhibit cell type-specific occupancy and positioning depending on differentiation status and transcription. This method is widely applicable to mapping many features of chromatin and discerning its landscape in large genomes at unprecedented resolution.

Nucleosome mapping across the CFTR locus identifies novel regulatory factors.

Nucleosome positioning on the chromatin strand plays a critical role in regulating accessibility of DNA to transcription factors and chromatin modifying enzymes. Hence, detailed information on nucleosome depletion or movement at cis-acting regulatory elements has the potential to identify predicted binding sites for trans-acting factors. Using a novel method based on enrichment of mononucleosomal DNA by bacterial artificial chromosome hybridization, we mapped nucleosome positions by deep sequencing across 250 kb, encompassing the cystic fibrosis transmembrane conductance regulator (CFTR) gene. CFTR shows tight tissue-specific regulation of expression, which is largely determined by cis-regulatory elements that lie outside the gene promoter. Although multiple elements are known, the repertoire of transcription factors that interact with these sites to activate or repress CFTR expression remains incomplete. Here, we show that specific nucleosome depletion corresponds to well-characterized binding sites for known trans-acting factors, including hepatocyte nuclear factor 1, Forkhead box A1 and CCCTC-binding factor. Moreover, the cell-type selective nucleosome positioning is effective in predicting binding sites for novel interacting factors, such as BAF155. Finally, we identify transcription factor binding sites that are overrepresented in regions where nucleosomes are depleted in a cell-specific manner. This approach recognizes the glucocorticoid receptor as a novel trans-acting factor that regulates CFTR expression in vivo.

Controls of nucleosome positioning in the human genome.

Nucleosomes are important for gene regulation because their arrangement on the genome can control which proteins bind to DNA. Currently, few human nucleosomes are thought to be consistently positioned across cells; however, this has been difficult to assess due to the limited resolution of existing data. We performed paired-end sequencing of micrococcal nuclease-digested chromatin (MNase-seq) from seven lymphoblastoid cell lines and mapped over 3.6 billion MNase-seq fragments to the human genome to create the highest-resolution map of nucleosome occupancy to date in a human cell type. In contrast to previous results, we find that most nucleosomes have more consistent positioning than expected by chance and a substantial fraction (8.7%) of nucleosomes have moderate to strong positioning. In aggregate, nucleosome sequences have 10 bp periodic patterns in dinucleotide frequency and DNase I sensitivity; and, across cells, nucleosomes frequently have translational offsets that are multiples of 10 bp. We estimate that almost half of the genome contains regularly spaced arrays of nucleosomes, which are enriched in active chromatin domains. Single nucleotide polymorphisms that reduce DNase I sensitivity can disrupt the phasing of nucleosome arrays, which indicates that they often result from positioning against a barrier formed by other proteins. However, nucleosome arrays can also be created by DNA sequence alone. The most striking example is an array of over 400 nucleosomes on chromosome 12 that is created by tandem repetition of sequences with strong positioning properties. In summary, a large fraction of nucleosomes are consistently positioned--in some regions because they adopt favored sequence positions, and in other regions because they are forced into specific arrangements by chromatin remodeling or DNA binding proteins.

Two-step mechanism for modifier of transcription 1 (Mot1) enzyme-catalyzed displacement of TATA-binding protein (TBP) from DNA.

The TATA box binding protein (TBP) is a central component of the transcription preinitiation complex, and its occupancy at a promoter is correlated with transcription levels. The TBP-promoter DNA complex contains sharply bent DNA and its interaction lifetime is limited by the ATP-dependent TBP displacement activity of the Snf2/Swi2 ATPase Mot1. Several mechanisms for Mot1 action have been proposed, but how it catalyzes TBP removal from DNA is unknown. To better understand the Mot1 mechanism, native gel electrophoresis and FRET were used to determine how Mot1 affects the trajectory of DNA in the TBP-DNA complex. Strikingly, in the absence of ATP, Mot1 acts to unbend DNA, whereas TBP remains closely associated with the DNA in a stable Mot1-TBP-DNA ternary complex. Interestingly, and in contrast to full-length Mot1, the isolated Mot1 ATPase domain binds DNA, and its affinity for DNA is nucleotide-dependent, suggesting parallels between the Mot1 mechanism and DNA translocation-based mechanisms of chromatin remodeling enzymes. Based on these findings, a model is presented for Mot1 that links a DNA conformational change with ATP-induced DNA translocation.

Archaeal nucleosome positioning in vivo and in vitro is directed by primary sequence motifs.

BACKGROUND: Histone wrapping of DNA into nucleosomes almost certainly evolved in the Archaea, and predates Eukaryotes. In Eukaryotes, nucleosome positioning plays a central role in regulating gene expression and is directed by primary sequence motifs that together form a nucleosome positioning code. The experiments reported were undertaken to determine if archaeal histone assembly conforms to the nucleosome positioning code. RESULTS: Eukaryotic nucleosome positioning is favored and directed by phased helical repeats of AA/TT/AT/TA and CC/GG/CG/GC dinucleotides, and disfavored by longer AT-rich oligonucleotides. Deep sequencing of genomic DNA protected from micrococcal nuclease digestion by assembly into archaeal nucleosomes has established that archaeal nucleosome assembly is also directed and positioned by these sequence motifs, both in vivo in Methanothermobacter thermautotrophicus and Thermococcus kodakarensis and in vitro in reaction mixtures containing only one purified archaeal histone and genomic DNA. Archaeal nucleosomes assembled at the same locations in vivo and in vitro, with much reduced assembly immediately upstream of open reading frames and throughout the ribosomal rDNA operons. Providing further support for a common positioning code, archaeal histones assembled into nucleosomes on eukaryotic DNA and eukaryotic histones into nucleosomes on archaeal DNA at the same locations. T. kodakarensis has two histones, designated HTkA and HTkB, and strains with either but not both histones deleted grow normally but do exhibit transcriptome differences. Comparisons of the archaeal nucleosome profiles in the intergenic regions immediately upstream of genes that exhibited increased or decreased transcription in the absence of HTkA or HTkB revealed substantial differences but no consistent pattern of changes that would correlate directly with archaeal nucleosome positioning inhibiting or stimulating transcription. CONCLUSIONS: The results obtained establish that an archaeal histone and a genome sequence together are sufficient to determine where archaeal nucleosomes preferentially assemble and where they avoid assembly. We confirm that the same nucleosome positioning code operates in Archaea as in Eukaryotes and presumably therefore evolved with the histone-fold mechanism of DNA binding and compaction early in the archaeal lineage, before the divergence of Eukaryotes.

p53 binds preferentially to genomic regions with high DNA-encoded nucleosome occupancy.

The human transcription factor TP53 is a pivotal roadblock against cancer. A key unresolved question is how the p53 protein selects its genomic binding sites in vivo out of a large pool of potential consensus sites. We hypothesized that chromatin may play a significant role in this site-selection process. To test this, we used a custom DNA microarray to measure p53 binding at approximately 2000 sites predicted to possess high-sequence specificity, and identified both strongly bound and weakly bound sites. When placed within a plasmid, weakly bound sites become p53 responsive and regain p53 binding when stably integrated into random genomic locations. Notably, strongly bound sites reside preferentially within genomic regions whose DNA sequence is predicted to encode relatively high intrinsic nucleosome occupancy. Using in vivo nucleosome occupancy measurements under conditions where p53 is inactive, we experimentally confirmed this prediction. Furthermore, upon p53 activation, nucleosomes are partially displaced from a relatively broad region surrounding the bound p53 sites, and this displacement is rapidly reversed upon inactivation of p53. Thus, in contrast to the general assumption that transcription-factor binding is preferred in sites that have low nucleosome occupancy prior to factor activation, we find that p53 binding occurs preferentially within a chromatin context of high intrinsic nucleosome occupancy.

What controls nucleosome positions?

The DNA of eukaryotic genomes is wrapped in nucleosomes, which strongly distort and occlude the DNA from access to most DNA-binding proteins. An understanding of the mechanisms that control nucleosome positioning along the DNA is thus essential to understanding the binding and action of proteins that carry out essential genetic functions. New genome-wide data on in vivo and in vitro nucleosome positioning greatly advance our understanding of several factors that can influence nucleosome positioning, including DNA sequence preferences, DNA methylation, histone variants and post-translational modifications, higher order chromatin structure, and the actions of transcription factors, chromatin remodelers and other DNA-binding proteins. We discuss how these factors function and ways in which they might be integrated into a unified framework that accounts for both the preservation of nucleosome positioning and the dynamic nucleosome repositioning that occur across biological conditions, cell types, developmental processes and disease.

A genomic code for nucleosome positioning.

Eukaryotic genomes are packaged into nucleosome particles that occlude the DNA from interacting with most DNA binding proteins. Nucleosomes have higher affinity for particular DNA sequences, reflecting the ability of the sequence to bend sharply, as required by the nucleosome structure. However, it is not known whether these sequence preferences have a significant influence on nucleosome position in vivo, and thus regulate the access of other proteins to DNA. Here we isolated nucleosome-bound sequences at high resolution from yeast and used these sequences in a new computational approach to construct and validate experimentally a nucleosome-DNA interaction model, and to predict the genome-wide organization of nucleosomes. Our results demonstrate that genomes encode an intrinsic nucleosome organization and that this intrinsic organization can explain approximately 50% of the in vivo nucleosome positions. This nucleosome positioning code may facilitate specific chromosome functions including transcription factor binding, transcription initiation, and even remodelling of the nucleosomes themselves.

Using DNA mechanics to predict in vitro nucleosome positions and formation energies.

In eukaryotic genomes, nucleosomes function to compact DNA and to regulate access to it both by simple physical occlusion and by providing the substrate for numerous covalent epigenetic tags. While competition with other DNA-binding factors and action of chromatin remodeling enzymes significantly affect nucleosome formation in vivo, nucleosome positions in vitro are determined by steric exclusion and sequence alone. We have developed a biophysical model, DNABEND, for the sequence dependence of DNA bending energies, and validated it against a collection of in vitro free energies of nucleosome formation and a set of in vitro nucleosome positions mapped at high resolution. We have also made a first ab initio prediction of nucleosomal DNA geometries, and checked its accuracy against the nucleosome crystal structure. We have used DNABEND to design both strong and weak histone- binding sequences, and measured the corresponding free energies of nucleosome formation. We find that DNABEND can successfully predict in vitro nucleosome positions and free energies, providing a physical explanation for the intrinsic sequence dependence of histone-DNA interactions.

Predicting nucleosome positioning using a duration Hidden Markov Model.

BACKGROUND: The nucleosome is the fundamental packing unit of DNAs in eukaryotic cells. Its detailed positioning on the genome is closely related to chromosome functions. Increasing evidence has shown that genomic DNA sequence itself is highly predictive of nucleosome positioning genome-wide. Therefore a fast software tool for predicting nucleosome positioning can help understanding how a genome's nucleosome organization may facilitate genome function. RESULTS: We present a duration Hidden Markov model for nucleosome positioning prediction by explicitly modeling the linker DNA length. The nucleosome and linker models trained from yeast data are re-scaled when making predictions for other species to adjust for differences in base composition. A software tool named NuPoP is developed in three formats for free download. CONCLUSIONS: Simulation studies show that modeling the linker length distribution and utilizing a base composition re-scaling method both improve the prediction of nucleosome positioning regarding sensitivity and false discovery rate. NuPoP provides a user-friendly software tool for predicting the nucleosome occupancy and the most probable nucleosome positioning map for genomic sequences of any size. When compared with two existing methods, NuPoP shows improved performance in sensitivity.

H2A.Z-mediated localization of genes at the nuclear periphery confers epigenetic memory of previous transcriptional state.

Many genes are recruited to the nuclear periphery upon transcriptional activation. The mechanism and functional significance of this recruitment is unclear. We find that recruitment of the yeast INO1 and GAL1 genes to the nuclear periphery is rapid and independent of transcription. Surprisingly, these genes remain at the periphery for generations after they are repressed. Localization at the nuclear periphery serves as a form of memory of recent transcriptional activation, promoting reactivation. Previously expressed GAL1 at the nuclear periphery is activated much more rapidly than long-term repressed GAL1 in the nucleoplasm, even after six generations of repression. Localization of INO1 at the nuclear periphery is necessary and sufficient to promote more rapid activation. This form of transcriptional memory is chromatin based; the histone variant H2A.Z is incorporated into nucleosomes within the recently repressed INO1 promoter and is specifically required for rapid reactivation of both INO1 and GAL1. Furthermore, H2A.Z is required to retain INO1 at the nuclear periphery after repression. Therefore, H2A.Z-mediated localization of recently repressed genes at the nuclear periphery represents an epigenetic state that confers memory of transcriptional activation and promotes reactivation.

Rapid spontaneous accessibility of nucleosomal DNA.

DNA wrapped in nucleosomes is sterically occluded, creating obstacles for proteins that must bind it. How proteins gain access to DNA buried inside nucleosomes is not known. Here we report measurements of the rates of spontaneous nucleosome conformational changes in which a stretch of DNA transiently unwraps off the histone surface, starting from one end of the nucleosome, and then rewraps. The rates are rapid. Nucleosomal DNA remains fully wrapped for only approximately 250 ms before spontaneously unwrapping; unwrapped DNA rewraps within approximately 10-50 ms. Spontaneous unwrapping of nucleosomal DNA allows any protein rapid access even to buried stretches of the DNA. Our results explain how remodeling factors can be recruited to particular nucleosomes on a biologically relevant timescale, and they imply that the major impediment to entry of RNA polymerase into a nucleosome is rewrapping of nucleosomal DNA, not unwrapping.

Distinct modes of regulation by chromatin encoded through nucleosome positioning signals.

The detailed positions of nucleosomes profoundly impact gene regulation and are partly encoded by the genomic DNA sequence. However, less is known about the functional consequences of this encoding. Here, we address this question using a genome-wide map of approximately 380,000 yeast nucleosomes that we sequenced in their entirety. Utilizing the high resolution of our map, we refine our understanding of how nucleosome organizations are encoded by the DNA sequence and demonstrate that the genomic sequence is highly predictive of the in vivo nucleosome organization, even across new nucleosome-bound sequences that we isolated from fly and human. We find that Poly(dA:dT) tracts are an important component of these nucleosome positioning signals and that their nucleosome-disfavoring action results in large nucleosome depletion over them and over their flanking regions and enhances the accessibility of transcription factors to their cognate sites. Our results suggest that the yeast genome may utilize these nucleosome positioning signals to regulate gene expression with different transcriptional noise and activation kinetics and DNA replication with different origin efficiency. These distinct functions may be achieved by encoding both relatively closed (nucleosome-covered) chromatin organizations over some factor binding sites, where factors must compete with nucleosomes for DNA access, and relatively open (nucleosome-depleted) organizations over other factor sites, where factors bind without competition.

Preferentially quantized linker DNA lengths in Saccharomyces cerevisiae.

The exact lengths of linker DNAs connecting adjacent nucleosomes specify the intrinsic three-dimensional structures of eukaryotic chromatin fibers. Some studies suggest that linker DNA lengths preferentially occur at certain quantized values, differing one from another by integral multiples of the DNA helical repeat, approximately 10 bp; however, studies in the literature are inconsistent. Here, we investigate linker DNA length distributions in the yeast Saccharomyces cerevisiae genome, using two novel methods: a Fourier analysis of genomic dinucleotide periodicities adjacent to experimentally mapped nucleosomes and a duration hidden Markov model applied to experimentally defined dinucleosomes. Both methods reveal that linker DNA lengths in yeast are preferentially periodic at the DNA helical repeat ( approximately 10 bp), obeying the forms 10n+5 bp (integer n). This 10 bp periodicity implies an ordered superhelical intrinsic structure for the average chromatin fiber in yeast.