A high-resolution map of nucleosome positioning on a fission yeast centromere

Epigenetic engineering of yeast reveals dynamic molecular adaptation to methylation stress and genetic modulators of specific DNMT3 family members

Cytosine methylation is a ubiquitous modification in mammalian DNA generated and maintained by several DNA methyltransferases (DNMTs) with partially overlapping functions and genomic targets. To systematically dissect the factors specifying each DNMT's activity, we engineered combinatorial knock-in of human DNMT genes in Komagataella phaffii, a yeast species lacking endogenous DNA methylation. Time-course expression measurements captured dynamic network-level adaptation of cells to DNMT3B1-induced DNA methylation stress and showed that coordinately modulating the availability of S-adenosyl methionine (SAM), the essential metabolite for DNMT-catalyzed methylation, is an evolutionarily conserved epigenetic stress response, also implicated in several human diseases. Convolutional neural networks trained on genome-wide CpG-methylation data learned distinct sequence preferences of DNMT3 family members. A simulated annealing interpretation method resolved these preferences into individual flanking nucleotides and periodic poly(A) tracts that rotationally position highly methylated cytosines relative to phased nucleosomes. Furthermore, the nucleosome repeat length defined the spatial unit of methylation spreading. Gene methylation patterns were similar to those in mammals, and hypo- and hypermethylation were predictive of increased and decreased transcription relative to control, respectively, in the absence of mammalian readers of DNA methylation. Introducing controlled epigenetic perturbations in yeast thus enabled characterization of fundamental genomic features directing specific DNMT3 proteins.

Application of MNase-Seq in the Global Mapping of Nucleosome Positioning in Plants

The precise positioning of nucleosomes along the underlying DNA is critical for a variety of biological processes, especially in regulating transcription. The interplay between nucleosomes and transcription factors for accessing the underlying DNA sequences is one of the key determinants that affect transcriptional regulation. Moreover, nucleosomes with various packing statuses confer distinct functions in regulating gene expressions in response to various internal or external signals. Therefore, global mapping of nucleosome positions is one informative way to elucidate the relationship between patterns of nucleosome positioning/occupancy and transcriptional regulations. MNase digestion coupled with high-throughput sequencing (MNase-seq) has been utilized widely for global mapping of nucleosome positioning in eukaryotes that have a sequenced genome. We have developed a robust MNase-seq procedure in plants. It mainly includes plant nuclei isolation, treatment of purified nuclei with MNase, gel recovery of MNase-trimmed mononucleosomal DNA with an approximate size of 150 bp, MNase-seq library preparation followed by Illumina sequencing, and data analysis. MNase-seq has already been successfully applied to identify genome-wide nucleosome positioning in model plants, rice, and Arabidopsis thaliana.

Inner Kinetochore Protein Interactions with Regional Centromeres of Fission Yeast

Centromeres of the fission yeast Schizosaccharomyces pombe lack the highly repetitive sequences that make most other "regional" centromeres refractory to analysis. To map fission yeast centromeres, we applied H4S47C-anchored cleavage mapping and native and cross-linked chromatin immunoprecipitation with paired-end sequencing. H3 nucleosomes are nearly absent from the central domain, which is occupied by centromere-specific H3 (cenH3 or CENP-A) nucleosomes with two H4s per particle that are mostly unpositioned and are more widely spaced than nucleosomes elsewhere. Inner kinetochore proteins CENP-A, CENP-C, CENP-T, CENP-I, and Scm3 are highly enriched throughout the central domain except at tRNA genes, with no evidence for preferred kinetochore assembly sites. These proteins are weakly enriched and less stably incorporated in H3-rich heterochromatin. CENP-A nucleosomes protect less DNA from nuclease digestion than H3 nucleosomes, while CENP-T protects a range of fragment sizes. Our results suggest that CENP-T particles occupy linkers between CENP-A nucleosomes and that classical regional centromeres differ from other centromeres by the absence of CENP-A nucleosome positioning.

The yeast kinetochore - structural insights from optical microscopy

In eukaryotic organisms each chromosome is captured by spindle microtubules. This interaction is mediated by an evolutionarily conserved, multi-protein complex called kinetochore. The main function of the kinetochore is to maintain correct chromosome segregation during cell division by transporting each sister chromatid to the spindle pole bodies localized at the opposite sites of the yeast nucleus. The kinetochore of budding yeast (Saccharomyces cerevisiae) is often used as a model system due to its simple composition compared to higher eukaryotes. This review highlights results obtained using optical imaging that revealed relative positions and stoichiometry of the major components of the budding yeast kinetochore.

Plasticity and epigenetic inheritance of centromere-specific histone H3 (CENP-A)-containing nucleosome positioning in the fission yeast

Nucleosomes containing the specific histone H3 variant CENP-A mark the centromere locus on each chromatin and initiate kinetochore assembly. For the common type of regional centromeres, little is known in molecular detail of centromeric chromatin organization, its propagation through cell division, and how distinct organization patterns may facilitate kinetochore assembly. Here, we show that in the fission yeast S. pombe, a relatively small number of CENP-A/Cnp1 nucleosomes are found within the centromeric core and that their positioning relative to underlying DNA varies among genetically homogenous cells. Consistent with the flexible positioning of Cnp1 nucleosomes, a large portion of the endogenous centromere is dispensable for its essential activity in mediating chromosome segregation. We present biochemical evidence that Cnp1 occupancy directly correlates with silencing of the underlying reporter genes. Furthermore, using a newly developed pedigree analysis assay, we demonstrated the epigenetic inheritance of Cnp1 positioning and quantified the rate of occasional repositioning of Cnp1 nucleosomes throughout cell generations. Together, our results reveal the plasticity and the epigenetically inheritable nature of centromeric chromatin organization.

NSeq: a multithreaded Java application for finding positioned nucleosomes from sequencing data

We introduce NSeq, a fast and efficient Java application for finding positioned nucleosomes from the high-throughput sequencing of MNase-digested mononucleosomal DNA. NSeq includes a user-friendly graphical interface, computes false discovery rates (FDRs) for candidate nucleosomes from Monte Carlo simulations, plots nucleosome coverage and centers, and exploits the availability of multiple processor cores by parallelizing its computations. Java binaries and source code are freely available at https://github.com/songlab/NSeq. The software is supported on all major platforms equipped with Java Runtime Environment 6 or later.

Myb-domain protein Teb1 controls histone levels and centromere assembly in fission yeast

The TTAGGG motif is common to two seemingly unrelated dimensions of chromatin function-the vertebrate telomere repeat and the promoter regions of many Schizosaccharomyces pombe genes, including all of those encoding canonical histones. The essential S. pombe protein Teb1 contains two Myb-like DNA binding domains related to those found in telomere proteins and binds the human telomere repeat sequence TTAGGG. Here, we analyse Teb1 binding throughout the genome and the consequences of reduced Teb1 function. Chromatin immunoprecipitation (ChIP)-on-chip analysis reveals robust Teb1 binding at many promoters, notably including all of those controlling canonical histone gene expression. A hypomorphic allele, teb1-1, confers reduced binding and reduced levels of histone transcripts. Prompted by previously suggested connections between histone expression and centromere identity, we examined localization of the centromeric histone H3 variant Cnp1 and found reduced centromeric binding along with reduced centromeric silencing. These data identify Teb1 as a master regulator of histone levels and centromere identity.

Neocentromeres and epigenetically inherited features of centromeres

Neocentromeres are ectopic sites where new functional kinetochores assemble and permit chromosome segregation. Neocentromeres usually form following genomic alterations that remove or disrupt centromere function. The ability to form neocentromeres is conserved in eukaryotes ranging from fungi to mammals. Neocentromeres that rescue chromosome fragments in cells with gross chromosomal rearrangements are found in several types of human cancers, and in patients with developmental disabilities. In this review, we discuss the importance of neocentromeres to human health and evaluate recently developed model systems to study neocentromere formation, maintenance, and function in chromosome segregation. Additionally, studies of neocentromeres provide insight into native centromeres; analysis of neocentromeres found in human clinical samples and induced in model organisms distinguishes features of centromeres that are dependent on centromere DNA from features that are epigenetically inherited together with the formation of a functional kinetochore.

The process of kinetochore assembly in yeasts

High fidelity chromosome segregation is essential for efficient transfer of the genetic material from the mother to daughter cells. The kinetochore (KT), which connects the centromere DNA to the spindle apparatus, plays a pivotal role in this process. In spite of considerable divergence in the centromere DNA sequence, basic architecture of a KT is evolutionarily conserved from yeast to humans. However, the identification of a large number of KT proteins paved the way of understanding conserved and diverged regulatory steps that lead to the formation of a multiprotein KT super-complex on the centromere DNA in different organisms. Because it is a daunting task to summarize the entire spectrum of information in a minireview, we focus here on the recent understanding in the process of KT assembly in three yeasts: Saccharomyces cerevisiae, Schizosaccharomyces pombe and Candida albicans. Studies in these unicellular organisms suggest that although the basic process of KT assembly remains the same, the dependence of a conserved protein for its KT localization may vary in these organisms.

Chromatin architectures at fission yeast transcriptional promoters and replication origins

We have used micrococcal nuclease (MNase) digestion followed by deep sequencing in order to obtain a higher resolution map than previously available of nucleosome positions in the fission yeast, Schizosaccharomyces pombe. Our data confirm an unusually short average nucleosome repeat length, ∼152 bp, in fission yeast and that transcriptional start sites (TSSs) are associated with nucleosome-depleted regions (NDRs), ordered nucleosome arrays downstream and less regularly spaced upstream nucleosomes. In addition, we found enrichments for associated function in four of eight groups of genes clustered according to chromatin configurations near TSSs. At replication origins, our data revealed asymmetric localization of pre-replication complex (pre-RC) proteins within large NDRs—a feature that is conserved in fission and budding yeast and is therefore likely to be conserved in other eukaryotic organisms.

Integrity of chromatin and replicating DNA in nuclei released from fission yeast by semi-automated grinding in liquid nitrogen

Background

Studies of nuclear function in many organisms, especially those with tough cell walls, are limited by lack of availability of simple, economical methods for large-scale preparation of clean, undamaged nuclei.

Findings

Here we present a useful method for nuclear isolation from the important model organism, the fission yeast, Schizosaccharomyces pombe. To preserve in vivo molecular configurations, we flash-froze the yeast cells in liquid nitrogen. Then we broke their tough cell walls, without damaging their nuclei, by grinding in a precision-controlled motorized mortar-and-pestle apparatus. The cryo-ground cells were resuspended and thawed in a buffer designed to preserve nuclear morphology, and the nuclei were enriched by differential centrifugation. The washed nuclei were free from contaminating nucleases and have proven well-suited as starting material for genome-wide chromatin analysis and for preparation of fragile DNA replication intermediates.

Conclusions

We have developed a simple, reproducible, economical procedure for large-scale preparation of endogenous-nuclease-free, morphologically intact nuclei from fission yeast. With appropriate modifications, this procedure may well prove useful for isolation of nuclei from other organisms with, or without, tough cell walls.

Diversity in requirement of genetic and epigenetic factors for centromere function in fungi

A centromere is a chromosomal region on which several proteins assemble to form the kinetochore. The centromere-kinetochore complex helps in the attachment of chromosomes to spindle microtubules to mediate segregation of chromosomes to daughter cells during mitosis and meiosis. In several budding yeast species, the centromere forms in a DNA sequence-dependent manner, whereas in most other fungi, factors other than the DNA sequence also determine the centromere location, as centromeres were able to form on nonnative sequences (neocentromeres) when native centromeres were deleted in engineered strains. Thus, in the absence of a common DNA sequence, the cues that have facilitated centromere formation on a specific DNA sequence for millions of years remain a mystery. Kinetochore formation is facilitated by binding of a centromere-specific histone protein member of the centromeric protein A (CENP-A) family that replaces a canonical histone H3 to form a specialized centromeric chromatin structure. However, the process of kinetochore formation on the rapidly evolving and seemingly diverse centromere DNAs in different fungal species is largely unknown. More interestingly, studies in various yeasts suggest that the factors required for de novo centromere formation (establishment) may be different from those required for maintenance (propagation) of an already established centromere. Apart from the DNA sequence and CENP-A, many other factors, such as posttranslational modification (PTM) of histones at centric and pericentric chromatin, RNA interference, and DNA methylation, are also involved in centromere formation, albeit in a species-specific manner. In this review, we discuss how several genetic and epigenetic factors influence the evolution of structure and function of centromeres in fungal species.

Absence of positive selection on centromeric histones in Tetrahymena suggests unsuppressed centromere: drive in lineages lacking male meiosis

Centromere-drive is a process where centromeres compete for transmission through asymmetric "female" meiosis for inclusion into the oocyte. In symmetric "male" meiosis, all meiotic products form viable germ cells. Therefore, the primary incentive for centromere-drive, a potential transmission bias, is believed to be missing from male meiosis. In this article, we consider whether male meiosis also bears the primary cost of centromere-drive. Because different taxa carry out different combinations of meiotic programs (symmetric + asymmetric, symmetric only, asymmetric only), it is possible to consider the evolutionary consequences of centromere-drive in the context of these differing systems. Groups with both types of meiosis have large, rapidly evolving centromeric regions, and their centromeric histones (CenH3s) have been shown to evolve under positive selection, suggesting roles as suppressors of centromere-drive. In contrast, taxa with only symmetric male meiosis have shown no evidence of positive selection in their centromeric histones. In this article, we present the first evolutionary analysis of centromeric histones in ciliated protozoans, a group that only undergoes asymmetric "female" meiosis. We find no evidence of positive selection acting on CNA1, the CenH3 of Tetrahymena species. Cytological observations of a panel of Tetrahymena species are consistent with dynamic karyotype evolution in this lineage. Our findings suggest that defects in male meiosis, and not mitosis or female meiosis, are the primary selective force behind centromere-drive suppression. Our study raises the possibility that taxa like ciliates, with only female meiosis, may therefore undergo unsuppressed centromere drive.

Heterochromatin is required for normal distribution of Neurospora crassa CenH3

Centromeres serve as platforms for the assembly of kinetochores and are essential for nuclear division. Here we identified Neurospora crassa centromeric DNA by chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) of DNA associated with tagged versions of the centromere foundation proteins CenH3 (CENP-A) and CEN-C (CENP-C) and the kinetochore protein CEN-T (CENP-T). On each chromosome we found an ∼150- to 300-kbp region of enrichment for all three proteins. These regions correspond to intervals predicted to be centromeric DNA by genetic mapping and DNA sequence analyses. By ChIP-seq we found extensive colocalization of CenH3, CEN-C, CEN-T, and histone H3K9 trimethylation (H3K9me3). In contrast, H3K4me2, which has been found at the cores of plant, fission yeast, Drosophila, and mammalian centromeres, was not enriched in Neurospora centromeric DNA. DNA methylation was most pronounced at the periphery of centromeric DNA. Mutation of dim-5, which encodes an H3K9 methyltransferase responsible for nearly all H3K9me3, resulted in altered distribution of CenH3-green fluorescent protein (GFP). Similarly, CenH3-GFP distribution was altered in the absence of HP1, the chromodomain protein that binds to H3K9me3. We conclude that eukaryotes with regional centromeres make use of different strategies for maintenance of CenH3 at centromeres, and we suggest a model in which centromere proteins nucleate at the core but require DIM-5 and HP1 for spreading.

Combinatorial, site-specific requirement for heterochromatic silencing factors in the elimination of nucleosome-free regions

High-resolution nucleosome occupancy maps of heterochromatic regions of wild-type and silencing-defective mutants of the fission yeast Schizosaccharomyces pombe revealed that heterochromatin induces the elimination of nucleosome-free regions (NFRs). NFRs associated with transcription initiation sites as well as those not associated with promoters are affected. We dissected the roles of the histone H3K9 methyltransferase Clr4 and the HP1 proteins Swi6 and Chp2, as well as the two catalytic activities of the SHREC histone deacetylase (HDAC)/ATPase effector complex. Strikingly, different DNA sites have distinct combinatorial requirements for these factors: Five classes of NFRs were identified that are eliminated by silencing factors through a mechanistic hierarchy governed by Clr4. The SHREC HDAC activity plays a major role in the elimination of class I-IV NFRs by antagonizing the action of RSC, a remodeling complex implicated in NFR formation. We propose that heterochromatin formation involves the deployment in several sequence-specific mechanisms to eliminate gaps between nucleosomes, thereby blocking access to the DNA.

Genome-wide maps of mononucleosomes and dinucleosomes containing hyperacetylated histones of Aspergillus fumigatus

It is suggested that histone modifications and/or histone variants influence the nucleosomal DNA length. We sequenced both ends of mononucleosomal and dinucleosomal DNA fragments of the filamentous fungus Aspergillus fumigatus, after treatment with the histone deacetylase inhibitor trichostatin A (TSA). After mapping the DNA fragments to the genome, we identified >7 million mononucleosome positions and >7 million dinucleosome positions. We showed that the distributions of the lengths of the mononucleosomal DNA fragments after 15-min and 30-min treatments with micrococcal nuclease (MNase) showed a single peak at 168 nt and 160 nt, respectively. The distributions of the lengths of the dinucleosomal DNA fragments after 15-min- and 30-min-treatment with MNase showed a single peak at 321 nt and 306 nt, respectively. The nucleosomal DNA fragments obtained from the TSA-treated cells were significantly longer than those obtained from the untreated cells. On the other hand, most of the genes did not undergo any change after treatment. Between the TSA-treated and untreated cells, only 77 genes had >or=2-fold change in expression levels. In addition, our results showed that the locations where mononucleosomes were frequently detected were conserved between the TSA-treated cells and untreated cells in the gene promoters (lower density of the nucleosomes). However, these locations were less conserved in the bodies (higher density of the nucleosomes) of genes with >or=2-fold changes. Our findings indicate that TSA influences the nucleosome positions, especially of the regions with high density of the nucleosomes by elongation of the nucleosomal DNA. However, most of the nucleosome positions are conserved in the gene promoters, even after treatment with TSA, because of the low density of nucleosomes in the gene promoters.

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.

The coexistence of the nucleosome positioning code with the genetic code on eukaryotic genomes

It is known that there are several codes residing simultaneously on the DNA double helix. The two best-characterized codes are the genetic code--the code for protein production, and the code for DNA packaging into nucleosomes. Since these codes have to coexist simultaneously on the same DNA region, both must be degenerate to allow this coexistence. A-tracts are homopolymeric stretches of several adjacent deoxyadenosines on one strand of the double helix, having unusual structural properties, which were shown to exclude nucleosomes and as such are instrumental in setting the translational positioning of DNA within nucleosomes. We observe, cross-kingdoms, a strong codon bias toward the avoidance of long A-tracts in exon regions, which enables the formation of high density of nucleosomes in these regions. Moreover, long A-tract avoidance is restricted exclusively to nucleosome-occupied exon regions. We show that this bias in codon usage is sufficient for enabling DNA organization within nucleosomes without constraints on the actual code for proteins. Thus, there is inter-dependency of the two major codes within DNA to allow their coexistence. Furthermore, we show that modulation of A-tract occurrences in exon versus non-exon regions may result in a unique alternation of the diameter of the '30-nm' fiber model.

Genome-wide maps of mono- and di-nucleosomes of Aspergillus fumigatus

We identified 5 975 567 mono- and 6 995 122 di-nucleosome positions of the fungus Aspergillus fumigatus, which was detected at high resolution based on the DNA sequence data obtained from both mono- and di-nucleosomal DNA fragments. We show that the distribution of lengths of the mono-nucleosomal DNA fragments has two peaks at 135 and 150 nt, whereas the distribution of di-nucleosomal DNA fragment lengths has a single peak at 285 nt. Although the gene bodies of the active and inactive genes and the inactive gene promoters had the two peaks of the mono-nucleosomal DNA fragment lengths, the active gene promoter lost the longer peak at 150 nt. Our findings strongly suggest that the nucleosomes protecting longer DNA fragments against MNase at the promoters, thereby inhibiting high gene expression.

The unique structure of A-tracts and intrinsic DNA bending

Short runs of adenines are a ubiquitous DNA element in regulatory regions of many organisms. When runs of 4–6 adenine base pairs (‘A-tracts’) are repeated with the helical periodicity, they give rise to global curvature of the DNA double helix, which can be macroscopically characterized by anomalously slow migration on polyacrylamide gels. The molecular structure of these DNA tracts is unusual and distinct from that of canonical B-DNA. We review here our current knowledge about the molecular details of A-tract structure and its interaction with sequences flanking them of either side and with the environment. Various molecular models were proposed to describe A-tract structure and how it causes global deflection of the DNA helical axis. We review old and recent findings that enable us to amalgamate the various findings to one model that conforms to the experimental data. Sequences containing phased repeats of A-tracts have from the very beginning been synonymous with global intrinsic DNA bending. In this review, we show that very often it is the unique structure of A-tracts that is at the basis of their widespread occurrence in regulatory regions of many organisms. Thus, the biological importance of A-tracts may often be residing in their distinct structure rather than in the global curvature that they induce on sequences containing them.

Structural and temporal regulation of centromeric chromatin

Normal inheritance of genetic material requires that chromosomes segregate faithfully during mitosis and meiosis. The kinetochore is a unique structure that attaches chromosomes to the microtubule spindle, monitors proper chromosome attachment to the spindle through the mitotic checkpoint, and couples spindle and motor protein forces to move chromosomes during prometaphase and anaphase. The centromere is a specialized chromosomal site that is the structural and functional foundation for kinetochore formation, and is characterized by a unique type of chromatin that needs to be reconstituted after each replication cycle. In this review, recent progress in understanding the structural nature of this chromatin and how it is specifically maintained through cell division are discussed.

High-resolution mapping of sequence-directed nucleosome positioning on genomic DNA

We have mapped in vitro nucleosome positioning on the sheep beta-lactoglobulin gene using high-throughput sequencing to characterise the DNA sequences recovered from reconstituted nucleosomes. This methodology surpasses previous approaches for coverage, accuracy and resolution and, most importantly, offers a simple yet rapid and relatively inexpensive method to characterise genomic DNA sequences in terms of nucleosome positioning capacity. We demonstrate an unambiguous correspondence between in vitro and in vivo nucleosome positioning around the promoter of the gene; identify discrete, sequence-specific nucleosomal structures above the level of the canonical core particle-a feature that has implications for regulatory protein access and higher-order chromatin packing; and reveal new insights into the involvement of periodically organised dinucleotide sequence motifs of the type GG and CC and not AA and TT, as determinants of nucleosome positioning-an observation that supports the idea that the core histone octamer can exploit different patterns of sequence organisation, or structural potential, in the DNA to bring about nucleosome positioning.

Fission yeast Scm3: A CENP-A receptor required for integrity of subkinetochore chromatin

The mechanisms ensuring specific incorporation of CENP-A at centromeres are poorly understood. Mis16 and Mis18 are required for CENP-A localization at centromeres and form a complex that is conserved from fission yeast to human. Fission yeast sim1 mutants that alleviate kinetochore domain silencing are defective in Scm3^Sp, the ortholog of budding yeast Scm3^Sc. Scm3^Sp depends on Mis16/18 for its centromere localization and like them is recruited to centromeres in late anaphase. Importantly, Scm3^Sp coaffinity purifies with CENP-A^Cnp1 and associates with CENP-A^Cnp1 in vitro, yet localizes independently of intact CENP-A^Cnp1 chromatin and is differentially released from chromatin. While Scm3^Sc has been proposed to form a unique hexameric nucleosome with CENP-A^Cse4 and histone H4 at budding yeast point centromeres, we favor a model in which Scm3^Sp acts as a CENP-A^Cnp1 receptor/assembly factor, cooperating with Mis16 and Mis18 to receive CENP-A^Cnp1 from the Sim3 escort and mediate assembly of CENP-A^Cnp1 into subkinetochore chromatin.

The impact of the nucleosome code on protein-coding sequence evolution in yeast

Coding sequence evolution was once thought to be the result of selection on optimal protein function alone. Selection can, however, also act at the RNA level, for example, to facilitate rapid translation or ensure correct splicing. Here, we ask whether the way DNA works also imposes constraints on coding sequence evolution. We identify nucleosome positioning as a likely candidate to set up such a DNA-level selective regime and use high-resolution microarray data in yeast to compare the evolution of coding sequence bound to or free from nucleosomes. Controlling for gene expression and intra-gene location, we find a nucleosome-free "linker" sequence to evolve on average 5-6% slower at synonymous sites. A reduced rate of evolution in linker is especially evident at the 5' end of genes, where the effect extends to non-synonymous substitution rates. This is consistent with regular nucleosome architecture in this region being important in the context of gene expression control. As predicted, codons likely to generate a sequence unfavourable to nucleosome formation are enriched in linker sequence. Amino acid content is likewise skewed as a function of nucleosome occupancy. We conclude that selection operating on DNA to maintain correct positioning of nucleosomes impacts codon choice, amino acid choice, and synonymous and non-synonymous rates of evolution in coding sequence. The results support the exclusion model for nucleosome positioning and provide an alternative interpretation for runs of rare codons. As the intimate association of histones and DNA is a universal characteristic of genic sequence in eukaryotes, selection on coding sequence composition imposed by nucleosome positioning should be phylogenetically widespread.