Evolutionary footprints of nucleosome positions in yeast

Using genome-wide maps of nucleosome positions in yeast, we have analyzed the influence of chromatin structure on the molecular evolution of genomic DNA. We have observed, on average, 10-15% lower substitution rates in linker regions than in nucleosomal DNA. This widespread local rate heterogeneity represents an evolutionary footprint of nucleosome positions and reveals that nucleosome organization is a genomic feature conserved over evolutionary timescales.     

Predicting nucleosome positions on the DNA: combining intrinsic sequence preferences and remodeler activities

Nucleosome positions on the DNA are determined by the intrinsic affinities of histone proteins to a given DNA sequence and by the ATP-dependent activities of chromatin remodeling complexes that can translocate nucleosomes with respect to the DNA. Here, we report a theoretical approach that takes into account both contributions. In the theoretical analysis two types of experiments have been considered: in vitro experiments with a single reconstituted nucleosome and in vivo genome-scale mapping of nucleosome positions. The effect of chromatin remodelers was described by iteratively redistributing the nucleosomes according to certain rules until a new steady state was reached. Three major classes of remodeler activities were identified: (i) the establishment of a regular nucleosome spacing in the vicinity of a strong positioning signal acting as a boundary, (ii) the enrichment/depletion of nucleosomes through amplification of intrinsic DNA-sequence-encoded signals and (iii) the removal of nucleosomes from high-affinity binding sites. From an analysis of data for nucleosome positions in resting and activated human CD4⁺ T cells [Schones et al., Cell 132, p. 887] it was concluded that the redistribution of a nucleosome map to a new state is greatly facilitated if the remodeler complex translocates the nucleosome with a preferred directionality.     

Prediction of nucleosome positions in the yeast genome based on matched mirror position filtering

Nucleosome positioning can affect the accessibility of the underlying DNA to the nuclear environment and as such plays an essential role in the regulation of cellular processes. Specific patterns have been found in the underlying DNA sequences of the nucleosome, and one of the most important patterns includes dinucleotides distributed every 10 to 11 base pairs. Based on this property, we propose to match each dinucleotide in the sequence against its mirror occurrences for 10 to 11 base pairs on both left-hand and right­hand sides. A large number of matches in a local region will then signify the existence of a nucleosome. In this paper, we propose the matched mirror position filters for efficient matching of periodic dinucleotide patterns and computationally predict the nucleosome positions. Experimental results on the Saccharomyces cerevisiae (yeast) genome show that the proposed algorithm can predict nucleosome positions effectively. More than 50% of our predicted nucleosomes are within 35 base pairs of those detected by biological experiments.     

Protein-DNA interactions and nuclease-sensitive regions determine nucleosome positions on yeast plasmid chromatin

To study mechanisms of nucleosome positioning, small circular plasmids were constructed, assembled into chromatin in vivo in Saccharomyces cerevisiae, and their chromatin structures were analysed with respect to positions of nucleosomes and nuclease-sensitive regions. Plasmids used include insertions of the URA3 gene into the TRP1 gene of the TRP1ARS1 circular plasmid in the same (TRURAP) or opposite (TRARUP) orientation. The URA3 gene has six precisely positioned, stable nucleosomes flanked by nuclease-sensitive regions at the 5' and 3' ends of the gene. Three of these nucleosome positions do not depend on the flanking nuclease-sensitive regions, since they are formed at similar positions in a derivative plasmid (TUmidL) that contains the middle of the URA3 sequence but not the 5' and 3' ends. These positions are probably due to protein-DNA interactions. In both TRURAP and TRARUP, the positions of the nucleosomes on the TRP1 gene were, however, shifted compared with the positions on the parental TRP1ARS1 circle and TUmidL. These changes are interpreted to be due to changes in the positions of flanking nuclease-sensitive regions that might act as boundaries to position nucleosomes. Thus, two independent mechanisms for nucleosome positioning have been demonstrated in vivo. The ARS1 region contains the 3' end of the TRP1 gene and the putative origin of replication. Since in TRURAP and TRARUP the TRP1 gene is interrupted, but the ARS1 region remains nuclease sensitive, this non-nucleosomal conformation of the ARS1 region probably reflects a chromatin structure important for replication.     

DNA sequence plays a major role in determining nucleosome positions in yeast CUP1 chromatin.

The role of DNA sequence in determining nucleosome positions in vivo was investigated by comparing the positions adopted by nucleosomes reconstituted on a yeast plasmid in vitro using purified core histones with those in native chromatin containing the same DNA, described previously. Nucleosomes were reconstituted on a 2.5 kilobase pair DNA sequence containing the yeast TRP1ARS1 plasmid with CUP1 as an insert (TAC-DNA). Multiple, alternative, overlapping nucleosome positions were mapped on TAC-DNA. For the 58 positioned nucleosomes identified, the relative positioning strengths and the stabilities to salt and temperature were determined. These positions were, with a few exceptions, identical to those observed in native, remodeled TAC chromatin containing an activated CUP1 gene. Only some of these positions are utilized in native, unremodeled chromatin. These observations suggest that DNA sequence is likely to play a very important role in positioning nucleosomes in vivo. We suggest that events occurring in yeast CUP1 chromatin determine which positions are occupied in vivo and when they are occupied.     

Characteristics of nucleosome core DNA and their applications in predicting nucleosome positions

By analyzing dinucleotide position-frequency data of yeast nucleosome-bound DNA sequences, dinucleotide periodicities of core DNA sequences were investigated. Within frequency domains, weakly bound dinucleotides (AA, AT, and the combinations AA-TT-TA and AA-TT-TA-AT) present doublet peaks in a periodicity range of 10-11 bp, and strongly bound dinucleotides present a single peak. A time-frequency analysis, based on wavelet transformation, indicated that weakly bound dinucleotides of core DNA sequences were spaced smaller (∼10.3 bp) at the two ends, with larger (∼11.1 bp) spacing in the middle section. The finding was supported by DNA curvature and was prevalent in all core DNA sequences. Therefore, three approaches were developed to predict nucleosome positions. After analyzing a 2200-bp DNA sequence, results indicated that the predictions were feasible; areas near protein-DNA binding sites resulted in periodicity profiles with irregular signals. The effects of five dinucleotide patterns were evaluated, indicating that the AA-TT pattern exhibited better performance. A chromosome-scale prediction demonstrated that periodicity profiles perform better than previously described, with up to 59% accuracy. Based on predictions, nucleosome distributions near the beginning and end of open reading frames were analyzed. Results indicated that the majority of open reading frames’ start and end sites were occupied by nucleosomes.     

Global nucleosome occupancy in yeast

Background  

Although eukaryotic genomes are generally thought to be entirely chromatin-associated, the activated PHO5 promoter in yeast is largely devoid of nucleosomes. We systematically evaluated nucleosome occupancy in yeast promoters by immunoprecipitating nucleosomal DNA and quantifying enrichment by microarrays.     

Tension-dependent nucleosome remodeling at the pericentromere in yeast

Nucleosome positioning is important for the structural integrity of chromosomes. During metaphase the mitotic spindle exerts physical force on pericentromeric chromatin. The cell must adjust the pericentromeric chromatin to accommodate the changing tension resulting from microtubule dynamics to maintain a stable metaphase spindle. Here we examine the effects of spindle-based tension on nucleosome dynamics by measuring the histone turnover of the chromosome arm and the pericentromere during metaphase in the budding yeast Saccharomyces cerevisiae. We find that both histones H2B and H4 exhibit greater turnover in the pericentromere during metaphase. Loss of spindle-based tension by treatment with the microtubule-depolymerizing drug nocodazole or compromising kinetochore function results in reduced histone turnover in the pericentromere. Pericentromeric histone dynamics are influenced by the chromatin-remodeling activities of STH1/NPS1 and ISW2. Sth1p is the ATPase component of the Remodels the Structure of Chromatin (RSC) complex, and Isw2p is an ATP-dependent DNA translocase member of the Imitation Switch (ISWI) subfamily of chromatin-remodeling factors. The balance between displacement and insertion of pericentromeric histones provides a mechanism to accommodate spindle-based tension while maintaining proper chromatin packaging during mitosis.     

Predicting nucleosome positioning in genomes: physical and bioinformatic approaches

In eukaryotic genomes, nucleosomes are responsible for packaging DNA and controlling gene expression. For this reason, an increasing interest is arising on computational methods capable of predicting the nucleosome positioning along genomes. In this review we describe and compare bioinformatic and physical approaches adopted to predict nucleosome occupancy along genomes. Computational analyses attempt at decoding the experimental nucleosome maps of genomes in terms of certain dinucleotide step periodicity observed along DNA. Such investigations show that highly significant information about the occurrence of a nucleosome along DNA is intrinsic in certain features of the sequence suggesting that DNA of eukaryotic genomes encodes nucleosome organization. Besides the bioinformatic approaches, physical models were proposed based on the sequence dependent conformational features of the DNA chain, which govern the free energy needed to transform recurrent DNA tracts along the genome into the nucleosomal shape.     

Self-association of the yeast nucleosome assembly protein 1

The self-association properties of the yeast nucleosome assembly protein 1 (yNAP1) have been investigated using biochemical and biophysical methods. Protein cross-linking and calibrated gel filtration chromatography of yNAP1 indicate the protein exists as a complex mixture of species at physiologic ionic strength (75-150 mM). Sedimentation velocity reveals a distribution of species of 4.5-12 Svedbergs (S) over a 50-fold range of concentrations. The solution-state complexity is reduced at higher ionic strength, allowing for examination of the fundamental oligomer. Sedimentation equilibrium of a homogeneous 4.5 S population at 500 mM sodium chloride reveals these species to be yNAP1 dimers. These dimers self-associate to form higher order oligomers at more moderate ionic strength. Titration of guanidine hydrochloride converts the higher order oligomers to the homogeneous 4.5 S dimer and then converts the 4.5 S dimers to 2.5 S monomers. Circular dichroism shows that guanidine-mediated dissociation of higher order oligomers into yNAP1 dimers is accompanied by only slight changes in secondary structure. Dissociation of the dimer requires a nearly complete denaturation event.     

酵母基因组中核小体偏好展列的识别

应用多样性增量结合二次判别分析（Increment of Diversity with Quadratic Discriminant analysis, IDQD）方法,对酵母基因组中的核小体强/弱偏好序列进行了识别。10交叉检验的预测成功率超过了97％,受试者操作特性（receiver operating characteristic,ROC）曲线下面积达到了0．99,预测成功率高于现有SVM算法。最后利用构建好的分类器对酵母基因组中三类包含TATA盒基因的起始密码子ATC上游400nt下游100nt区域进行了分析。结果表明,IDQD算法有能力应用于基因组中核小体序列的识别。     

Chromatin folding modulates nucleosome positioning in yeast minichromosomes

Based on the chromatin structures of the yeast URA3 gene and the TRP1ARS1 circle, we have designed circular minichromosomes of different sizes that should each form a tight tetranucleosome. This structure was assumed to be stiff and bulky and therefore likely to be sensitive to packaging into a three-dimensional structure. The structures of the minichromosomes were determined using micrococcal nuclease. Only one of the minichromosomes showed a protected region of about 570 bp, compatible with the predicted tight tetranucleosome, while all other constructs showed alternative structures. A comparison of the structures revealed that neither histone-DNA interactions nor influences from flanking boundaries are sufficient determinants of nucleosome positions. The data strongly suggest that chromatin folding modulates the nucleosome arrangement along the DNA.     

酵母基因组中核小体偏好序列的识别

应用多样性增量结合二次判别分析(Increment of Diversity with Quadratic Discriminant analysis,IDQD)方法,对酵母基因组中的核小体强/弱偏好序列进行了识别.10交叉检验的预测成功率超过了97%,受试者操作特性(receiver operating characteristic,ROC)曲线下面积达到了0.99,预测成功率高于现有SVM算法.最后利用构建好的分类器对酵母基因组中三类包含TATA盒基因的起始密码子ATG上游400nt下游100nt区域进行了分析.结果表明,IDQD算法有能力应用于基因组中核小体序列的识别.     

Biophysical characterization of the centromere-specific nucleosome from budding yeast

The centromeric DNA of all eukaryotes is assembled upon a specialized nucleosome containing a histone H3 variant known as CenH3. Despite the importance and conserved nature of this protein, the characteristics of the centromeric nucleosome are still poorly understood. In particular, the stoichiometry and DNA-binding properties of the CenH3 nucleosome have been the subject of some debate. We have characterized the budding yeast centromeric nucleosome by biochemical and biophysical methods and show that it forms a stable octamer containing two copies of the Cse4 protein and wraps DNA in a left-handed supercoil, similar to the canonical H3 nucleosome. The DNA-binding properties of the recombinant nucleosome are identical to those observed in vivo demonstrating that the octameric structure is physiologically relevant.     

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.