Nucleosome positioning, nucleosome spacing and the nucleosome code

Genome-wide mapping of nucleosomes in yeast using paired-end sequencing

The DNA of eukaryotic cells is packaged into chromatin by histone proteins, which play a central role in regulating access to genetic information. The nucleosome core is the basic structural unit of chromatin: it is composed of an octamer of the four major core histones (two molecules each of H2A, H2B, H3, and H4), around which are wrapped ∼1.75 negative superhelical turns of DNA, a total of 145-147bp. Nucleosome cores are regularly spaced along the DNA in vivo, separated by linker DNA. Nucleosomes are compact structures capable of blocking access to the DNA that they contain. For example, they may prevent the binding of transcription factors to their cognate sites. It is therefore very important to obtain quantitative information on the positions of nucleosomes with respect to regulatory regions in vivo. The advent of high-throughput sequencing methods has revolutionized this field. We describe the use and advantages of paired-end sequencing to map nucleosomal DNA obtained by micrococcal nuclease digestion of budding yeast nuclei. This approach provides high-quality genome-wide nucleosome occupancy and position maps.

The histone variant H2A.Z interconverts two stable epigenetic chromatin states

The nucleosomes occupying the chromosomal start sites of transcription contain the histone H2A variant H2A.Z in place of H2A. Upon galactose induction, nucleosomes are evicted from the GAL1 locus in Saccharomyces cerevisiae cells. H2A.Z (which is encoded by the HTZ1 gene in S. cerevisiae) is required for the eviction of the GAL1 promoter nucleosome and for the transcriptional activation of the GAL1 gene; however, histones are also important for transcriptional repression and we asked in the present paper if H2A.Z also plays a role in the glucose repression of the GAL1 promoter. With the help of a fusion of the URA3 ORF (open reading frame) to the GAL1 promoter, we were able to detect two different epigenetic transcription states of the GAL1 promoter in glucose-grown cells lacking H2A.Z: a repressed state that is occupied by a H2A-containing nucleosome and a derepressed state that is nucleosome-free. These two chromatin states are inherited stably through many cell divisions. According to the model described in the present paper, the role of H2A.Z is to facilitate the addition and removal of promoter nucleosomes and to prevent the formation of unfavourable stable epigenetic chromatin structures, which are not in accordance with the environmental conditions.

Finding a balance: how diverse dosage compensation strategies modify histone h4 to regulate transcription

Dosage compensation balances gene expression levels between the sex chromosomes and autosomes and sex-chromosome-linked gene expression levels between the sexes. Different dosage compensation strategies evolved in different lineages, but all involve changes in chromatin. This paper discusses our current understanding of how modifications of the histone H4 tail, particularly changes in levels of H4 lysine 16 acetylation and H4 lysine 20 methylation, can be used in different contexts to either modulate gene expression levels twofold or to completely inhibit transcription.

Predicting sumoylation site by feature selection method

The small ubiquitin-like modifier (SUMO) proteins are a kind of proteins that can be attached to a series of proteins. The sumoylation of protein is an important posttranslational modification. Thus, the prediction of the sumoylation site of a given protein is significant. Here we employed a combined method to perform this task. We predicted the sumoylation site of a protein by a two-staged procedure. At the first stage, whether a protein would be sumoylated was predicted; whereas at the second stage, the sumoylation sites of the protein were predicted if it was determined to be modified by SUMO at the first stage. At the first stage, we encoded a protein with protein families (PFAM) and trained the predictor with nearest network algorithm (NNA); at the second stage, we encoded nonapeptides (peptides that contain nine residues) of the protein containing the lysine residues, with Amino Acid Index, and trained the predictor with NNA. The predictor was tested by the k-fold cross-validation method. The highest accuracy of the second-staged predictor was 99.55% when 12 features were incorporated in the predictor. The corresponding Matthews Correlation Coefficient was 0.7952. These results indicate that the method is a promising tool to predict the sumoylation site of a protein. At last, the features used in the predictor are discussed. The software is available at request.

Activation-induced disruption of nucleosome position clusters on the coding regions of Gcn4-dependent genes extends into neighbouring genes

We have used paired-end sequencing of yeast nucleosomal DNA to obtain accurate genomic maps of nucleosome positions and occupancies in control cells and cells treated with 3-aminotriazole (3AT), an inducer of the transcriptional activator Gcn4. In control cells, 3AT-inducible genes exhibit a series of distinct nucleosome occupancy peaks. However, the underlying position data reveal that each nucleosome peak actually consists of a cluster of mutually exclusive overlapping positions, usually including a dominant position. Thus, each nucleosome occupies one of several possible positions and consequently, different cells have distinct local chromatin structures. Induction results in a major disruption of nucleosome positioning, sometimes with altered spacing and a dramatic loss of occupancy over the entire gene, often extending into a neighbouring gene. Nucleosome-depleted regions are generally unaffected. Genes repressed by 3AT show the same changes, but in reverse. We propose that yeast genes exist in one of several alternative nucleosomal arrays, which are disrupted by activation. We conclude that activation results in gene-wide chromatin remodelling and that this remodelling can even extend into the chromatin of flanking genes.

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.

The centromeric nucleosome of budding yeast is perfectly positioned and covers the entire centromere

The centromeres of budding yeast are ~120 bp in size and contain three functional elements: an AT-rich region flanked by binding sites for Cbf1 and CBF3. A specialized nucleosome containing the H3 variant Cse4 (CenH3) is formed at the centromere. Our genome-wide paired-end sequencing of nucleosomal DNA reveals that the centromeric nucleosome contains a micrococcal nuclease-resistant kernel of 123-135 bp, depending on the centromere, and is therefore significantly shorter than the canonical nucleosome. Unlike canonical nucleosomes, the centromeric nucleosome is essentially perfectly positioned. The entire centromere is included, together with at least 1 bp of DNA upstream of the Cbf1 site and at least 4 bp downstream of the CBF3 site. The fact that the binding sites for Cbf1 and CBF3 are included within the centromeric nucleosome has important implications for models of the centromeric nucleosome and for kinetochore function.

Nucleosome positioning pattern derived from oligonucleotide compositions of genomic sequences

Availability of nucleosome positioning pattern(s) is crucial for chromatin studies. The matrix form of the pattern has been recently derived (I. Gabdank, D. Barash, E. N. Trifonov. J Biomol Struct Dyn 26, 403-412 (2009), and E. N. Trifonov. J Biomol Struct Dyn 27, 741-746 (2010)). In its simplified linear form it is described by the motif CGRAAATTTYCG. Oligonucleotide components of the motif (say, triplets GRA, RAA, AAA, etc.) would be expected to appear in eukaryotic sequences more frequently. In this work we attempted the reconstruction of the bendability patterns for 13 genomes by a novel approach-extension of highest frequency trinucleotides. The consensus of the patterns reconstructed on the basis of trinucleotide frequencies in 13 eukaryotic genomes is derived: CRAAAATTTTYG. It conforms to the earlier established sequence motif. The reconstruction, thus, attests to the universality of the nucleosome DNA bendability pattern.

A polar barrier to transcription can be circumvented by remodeler-induced nucleosome translocation

Many eukaryotic genes are regulated at the level of transcript elongation. Nucleosomes are likely targets for this regulation. Previously, we have shown that nucleosomes formed on very strong positioning sequences (601 and 603), present a high, orientation-dependent barrier to transcription by RNA polymerase II in vitro. The existence of this polar barrier correlates with the interaction of a 16-bp polar barrier signal (PBS) with the promoter-distal histone H3-H4 dimer. Here, we show that the polar barrier is relieved by ISW2, an ATP-dependent chromatin remodeler, which translocates the nucleosome over a short distance, such that the PBS no longer interacts with the distal H3-H4 dimer, although it remains within the nucleosome. In vivo, insertion of the 603 positioning sequence into the yeast CUP1 gene results in a modest reduction in transcription, but this reduction is orientation-independent, indicating that the polar barrier can be circumvented. However, the 603-nucleosome is present at the expected position in only a small fraction of cells. Thus, the polar barrier is probably non-functional in vivo because the nucleosome is not positioned appropriately, presumably due to nucleosome sliding activities. We suggest that interactions between PBSs and chromatin remodelers might have significant regulatory potential.

High nucleosome occupancy is encoded at X-linked gene promoters in C. elegans

We mapped nucleosome occupancy by paired-end Illumina sequencing in C. elegans embryonic cells, adult somatic cells, and a mix of adult somatic and germ cells. In all three samples, the nucleosome occupancy of gene promoters on the X chromosome differed from autosomal promoters. While both X and autosomal promoters exhibit a typical nucleosome-depleted region upstream of transcript start sites and a well-positioned +1 nucleosome, X-linked gene promoters on average exhibit higher nucleosome occupancy relative to autosomal promoters. We show that the difference between X and autosomes does not depend on the somatic dosage compensation machinery. Instead, the chromatin difference at promoters is partly encoded by DNA sequence, because a model trained on nucleosome sequence preferences from S. cerevisiae in vitro data recapitulate nearly completely the experimentally observed difference between X and autosomal promoters. The model predictions also correlate very well with experimentally determined occupancy values genome-wide. The nucleosome occupancy differences observed on X promoters may bear on mechanisms of X chromosome dosage compensation in the soma, and chromosome-wide repression of X in the germline.