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

Retrotransposon targeting to RNA polymerase III-transcribed genes

Retrotransposons are genetic elements that are similar in structure and life cycle to retroviruses by replicating via an RNA intermediate and inserting into a host genome. The Saccharomyces cerevisiae (S. cerevisiae) Ty1-5 elements are long terminal repeat (LTR) retrotransposons that are members of the Ty1-copia (Pseudoviridae) or Ty3-gypsy (Metaviridae) families. Four of the five S. cerevisiae Ty elements are inserted into the genome upstream of RNA Polymerase (Pol) III-transcribed genes such as transfer RNA (tRNA) genes. This particular genomic locus provides a safe environment for Ty element insertion without disruption of the host genome and is a targeting strategy used by retrotransposons that insert into compact genomes of hosts such as S. cerevisiae and the social amoeba Dictyostelium. The mechanism by which Ty1 targeting is achieved has been recently solved due to the discovery of an interaction between Ty1 Integrase (IN) and RNA Pol III subunits. We describe the methods used to identify the Ty1-IN interaction with Pol III and the Ty1 targeting consequences if the interaction is perturbed. The details of Ty1 targeting are just beginning to emerge and many unexplored areas remain including consideration of the 3-dimensional shape of genome. We present a variety of other retrotransposon families that insert adjacent to Pol III-transcribed genes and the mechanism by which the host machinery has been hijacked to accomplish this targeting strategy. Finally, we discuss why retrotransposons selected Pol III-transcribed genes as a target during evolution and how retrotransposons have shaped genome architecture.

Viral proteins as a potential driver of histone depletion in dinoflagellates

Within canonical eukaryotic nuclei, DNA is packaged with highly conserved histone proteins into nucleosomes, which facilitate DNA condensation and contribute to genomic regulation. Yet the dinoflagellates, a group of unicellular algae, are a striking exception to this otherwise universal feature as they have largely abandoned histones and acquired apparently viral-derived substitutes termed DVNPs (dinoflagellate-viral-nucleoproteins). Despite the magnitude of this transition, its evolutionary drivers remain unknown. Here, using Saccharomyces cerevisiae as a model, we show that DVNP impairs growth and antagonizes chromatin by localizing to histone binding sites, displacing nucleosomes, and impairing transcription. Furthermore, DVNP toxicity can be relieved through histone depletion and cells diminish their histones in response to DVNP expression suggesting that histone reduction could have been an adaptive response to these viral proteins. These findings provide insights into eukaryotic chromatin evolution and highlight the potential for horizontal gene transfer to drive the divergence of cellular systems.

Major Determinants of Nucleosome Positioning

The compact structure of the nucleosome limits DNA accessibility and inhibits the binding of most sequence-specific proteins. Nucleosomes are not randomly located on the DNA but positioned with respect to the DNA sequence, suggesting models in which critical binding sites are either exposed in the linker, resulting in activation, or buried inside a nucleosome, resulting in repression. The mechanisms determining nucleosome positioning are therefore of paramount importance for understanding gene regulation and other events that occur in chromatin, such as transcription, replication, and repair. Here, we review our current understanding of the major determinants of nucleosome positioning: DNA sequence, nonhistone DNA-binding proteins, chromatin-remodeling enzymes, and transcription. We outline the major challenges for the future: elucidating the precise mechanisms of chromatin opening and promoter activation, identifying the complexes that occupy promoters, and understanding the multiscale problem of chromatin fiber organization.

Chromatin Higher-Order Folding: A Perspective with Linker DNA Angles

The mechanism by which the "beads-on-a-string" nucleosome chain folds into various higher-order chromatin structures in eukaryotic cell nuclei is still poorly understood. The various models depicting higher-order chromatin as regular helical fibers and the very opposite "polymer melt" theory imply that interactions between nucleosome "beads" make the main contribution to the chromatin compaction. Other models in which the geometry of linker DNA "strings" entering and exiting the nucleosome define the three-dimensional structure predict that small changes in the linker DNA configuration may strongly affect nucleosome chain folding and chromatin higher-order structure. Among those studies, the cross-disciplinary approach pioneered by Jörg Langowski that combines computational modeling with biophysical and biochemical experiments was most instrumental for understanding chromatin higher-order structure in vitro. Strikingly, many recent studies, including genome-wide nucleosome interaction mapping and chromatin imaging, show an excellent agreement with the results of three-dimensional computational modeling based on the primary role of linker DNA geometry in chromatin compaction. This perspective relates nucleosome array models with experimental studies of nucleosome array folding in vitro and in situ. I argue that linker DNA configuration plays a key role in determining nucleosome chain flexibility, topology, and propensity for self-association, thus providing new implications for regulation of chromatin accessibility to DNA binding factors and RNA transcription machinery as well as long-range communications between distant genomic sites.

Relationship between histone modifications and transcription factor binding is protein family specific

The very small fraction of putative binding sites (BSs) that are occupied by transcription factors (TFs) in vivo can be highly variable across different cell types. This observation has been partly attributed to changes in chromatin accessibility and histone modification (HM) patterns surrounding BSs. Previous studies focusing on BSs within DNA regulatory regions found correlations between HM patterns and TF binding specificities. However, a mechanistic understanding of TF-DNA binding specificity determinants is still not available. The ability to predict in vivo TF binding on a genome-wide scale requires the identification of features that determine TF binding based on evolutionary relationships of DNA binding proteins. To reveal protein family-dependent mechanisms of TF binding, we conducted comprehensive comparisons of HM patterns surrounding BSs and non-BSs with exactly matched core motifs for TFs in three cell lines: 33 TFs in GM12878, 37 TFs in K562, and 18 TFs in H1-hESC. These TFs displayed protein family-specific preferences for HM patterns surrounding BSs, with high agreement among cell lines. Moreover, compared to models based on DNA sequence and shape at flanking regions of BSs, HM-augmented quantitative machine-learning methods resulted in increased performance in a TF family-specific manner. Analysis of the relative importance of features in these models indicated that TFs, displaying larger HM pattern differences between BSs and non-BSs, bound DNA in an HM-specific manner on a protein family-specific basis. We propose that TF family-specific HM preferences reveal distinct mechanisms that assist in guiding TFs to their cognate BSs by altering chromatin structure and accessibility.

Precise genome-wide mapping of single nucleosomes and linkers in vivo

We developed a chemical cleavage method that releases single nucleosome dyad-containing fragments, allowing us to precisely map both single nucleosomes and linkers with high accuracy genome-wide in yeast. Our single nucleosome positioning data reveal that nucleosomes occupy preferred positions that differ by integral multiples of the DNA helical repeat. By comparing nucleosome dyad positioning maps to existing genomic and transcriptomic data, we evaluated the contributions of sequence, transcription, and histones H1 and H2A.Z in defining the chromatin landscape. We present a biophysical model that neglects DNA sequence and shows that steric occlusion suffices to explain the salient features of nucleosome positioning.

Assessing sufficiency and necessity of enhancer activities for gene expression and the mechanisms of transcription activation

Enhancers are important genomic regulatory elements directing cell type-specific transcription. They assume a key role during development and disease, and their identification and functional characterization have long been the focus of scientific interest. The advent of next-generation sequencing and clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9-based genome editing has revolutionized the means by which we study enhancer biology. In this review, we cover recent developments in the prediction of enhancers based on chromatin characteristics and their identification by functional reporter assays and endogenous DNA perturbations. We discuss that the two latter approaches provide different and complementary insights, especially in assessing enhancer sufficiency and necessity for transcription activation. Furthermore, we discuss recent insights into mechanistic aspects of enhancer function, including findings about cofactor requirements and the role of post-translational histone modifications such as monomethylation of histone H3 Lys4 (H3K4me1). Finally, we survey how these approaches advance our understanding of transcription regulation with respect to promoter specificity and transcriptional bursting and provide an outlook covering open questions and promising developments.

A history of genome editing in Saccharomyces cerevisiae

Genome editing is a form of highly precise genetic engineering which produces alterations to an organism's genome as small as a single base pair with no incidental or auxiliary modifications; this technique is crucial to the field of synthetic biology, which requires such precision in the installation of novel genetic circuits into host genomes. While a new methodology for most organisms, genome editing capabilities have been used in the budding yeast Saccharomyces cerevisiae for decades. In this review, I will present a brief history of genome editing in S. cerevisiae, discuss the current gold standard method of Cas9-mediated genome editing, and speculate on future directions of the field.

The Latest Twists in Chromatin Remodeling

In its most restrictive interpretation, the notion of chromatin remodeling refers to the action of chromatin-remodeling enzymes on nucleosomes with the aim of displacing and removing them from the chromatin fiber (the effective polymer formed by a DNA molecule and proteins). This local modification of the fiber structure can have consequences for the initiation and repression of the transcription process, and when the remodeling process spreads along the fiber, it also results in long-range effects essential for fiber condensation. There are three regulatory levels of relevance that can be distinguished for this process: the intrinsic sequence preference of the histone octamer, which rules the positioning of the nucleosome along the DNA, notably in relation to the genetic information coded in DNA; the recognition or selection of nucleosomal substrates by remodeling complexes; and, finally, the motor action on the nucleosome exerted by the chromatin remodeler. Recent work has been able to provide crucial insights at each of these three levels that add new twists to this exciting and unfinished story, which we highlight in this perspective.

Generation of Native Chromatin Immunoprecipitation Sequencing Libraries for Nucleosome Density Analysis

We present a modified native chromatin immunoprecipitation sequencing (ChIP-seq) experimental protocol compatible with a Gaussian mixture distribution based analysis methodology (nucleosome density ChIP-seq; ndChIP-seq) that enables the generation of combined measurements of micrococcal nuclease (MNase) accessibility with histone modification genome-wide. Nucleosome position and local density, and the posttranslational modification of their histone subunits, act in concert to regulate local transcription states. Combinatorial measurements of nucleosome accessibility with histone modification generated by ndChIP-seq allows for the simultaneous interrogation of these features. The ndChIP-seq methodology is applicable to small numbers of primary cells inaccessible to cross-linking based ChIP-seq protocols. Taken together, ndChIP-seq enables the measurement of histone modification in combination with local nucleosome density to obtain new insights into shared mechanisms that regulate RNA transcription within rare primary cell populations.

Subtracting the sequence bias from partially digested MNase-seq data reveals a general contribution of TFIIS to nucleosome positioning

BACKGROUND

TFIIS stimulates RNA cleavage by RNA polymerase II and promotes the resolution of backtracking events. TFIIS acts in the chromatin context, but its contribution to the chromatin landscape has not yet been investigated. Co-transcriptional chromatin alterations include subtle changes in nucleosome positioning, like those expected to be elicited by TFIIS, which are elusive to detect. The most popular method to map nucleosomes involves intensive chromatin digestion by micrococcal nuclease (MNase). Maps based on these exhaustively digested samples miss any MNase-sensitive nucleosomes caused by transcription. In contrast, partial digestion approaches preserve such nucleosomes, but introduce noise due to MNase sequence preferences. A systematic way of correcting this bias for massively parallel sequencing experiments is still missing.

RESULTS

To investigate the contribution of TFIIS to the chromatin landscape, we developed a refined nucleosome-mapping method in Saccharomyces cerevisiae. Based on partial MNase digestion and a sequence-bias correction derived from naked DNA cleavage, the refined method efficiently mapped nucleosomes in promoter regions rich in MNase-sensitive structures. The naked DNA correction was also important for mapping gene body nucleosomes, particularly in those genes whose core promoters contain a canonical TATA element. With this improved method, we analyzed the global nucleosomal changes caused by lack of TFIIS. We detected a general increase in nucleosomal fuzziness and more restricted changes in nucleosome occupancy, which concentrated in some gene categories. The TATA-containing genes were preferentially associated with decreased occupancy in gene bodies, whereas the TATA-like genes did so with increased fuzziness. The detected chromatin alterations correlated with functional defects in nascent transcription, as revealed by genomic run-on experiments.

CONCLUSIONS

The combination of partial MNase digestion and naked DNA correction of the sequence bias is a precise nucleosomal mapping method that does not exclude MNase-sensitive nucleosomes. This method is useful for detecting subtle alterations in nucleosome positioning produced by lack of TFIIS. Their analysis revealed that TFIIS generally contributed to nucleosome positioning in both gene promoters and bodies. The independent effect of lack of TFIIS on nucleosome occupancy and fuzziness supports the existence of alternative chromatin dynamics during transcription elongation.

Regulation of tRNA gene transcription by the chromatin structure and nucleosome dynamics

The short, non-coding genes transcribed by the RNA polymerase (pol) III, necessary for survival of a cell, need to be repressed under the stress conditions in vivo. The pol III-transcribed genes have adopted several novel chromatin-based regulatory mechanisms to their advantage. In the budding yeast, the sub-nucleosomal size tRNA genes are found in the nucleosome-free regions, flanked by positioned nucleosomes at both the ends. With their chromosomes-wide distribution, all tRNA genes have a different chromatin context. A single nucleosome dynamics controls the accessibility of the genes for transcription. This dynamics operates under the influence of several chromatin modifiers in a gene-specific manner, giving the scope for differential regulation of even the isogenes within a tRNA gene family. The chromatin structure around the pol III-transcribed genes provides a context conducive for steady-state transcription as well as gene-specific transcriptional regulation upon signaling from the environmental cues. This article is part of a Special Issue entitled: SI: Regulation of tRNA synthesis and modification in physiological conditions and disease edited by Dr. Boguta Magdalena.

Remarkable Evolutionary Plasticity of Centromeric Chromatin

Centromeres were familiar to cell biologists in the late 19th century, but for most eukaryotes the basis for centromere specification has remained enigmatic. Much attention has been focused on the cenH3 (CENP-A) histone variant, which forms the foundation of the centromere. To investigate the DNA sequence requirements for centromere specification, we applied a variety of epigenomic approaches, which have revealed surprising diversity in centromeric chromatin properties. Whereas each point centromere of budding yeast is occupied by a single precisely positioned tetrameric nucleosome with one cenH3 molecule, the "regional" centromeres of fission yeast contain unphased presumably octameric nucleosomes with two cenH3s. In Caenorhabditis elegans, kinetochores assemble all along the chromosome at sites of cenH3 nucleosomes that resemble budding yeast point centromeres, whereas holocentric insects lack cenH3 entirely. The "satellite" centromeres of most animals and plants consist of cenH3-containing particles that are precisely positioned over homogeneous tandem repeats, but in humans, different α-satellite subfamilies are occupied by CENP-A nucleosomes with very different conformations. We suggest that this extraordinary evolutionary diversity of centromeric chromatin architectures can be understood in terms of the simplicity of the task of equal chromosome segregation that is continually subverted by selfish DNA sequences.

Resetting the Yeast Epigenome with Human Nucleosomes

Humans and yeast are separated by a billion years of evolution, yet their conserved histones retain central roles in gene regulation. Here, we "reset" yeast to use core human nucleosomes in lieu of their own (a rare event taking 20 days), which initially only worked with variant H3.1. The cells adapt by acquiring suppressor mutations in cell-division genes or by acquiring certain aneuploid states. Converting five histone residues to their yeast counterparts restored robust growth. We reveal that humanized nucleosomes are positioned according to endogenous yeast DNA sequence and chromatin-remodeling network, as judged by a yeast-like nucleosome repeat length. However, human nucleosomes have higher DNA occupancy, globally reduce RNA content, and slow adaptation to new conditions by delaying chromatin remodeling. These humanized yeasts (including H3.3) pose fundamental new questions about how chromatin is linked to many cell processes and provide a platform to study histone variants via yeast epigenome reprogramming.

A code for transcription elongation speed

The two major steps of gene expression are transcription and translation. While hundreds of studies regarding the effect of sequence features on the translation elongation process have been published, very few connect sequence features to the transcription elongation rate. We suggest, for the first time, that short transcript sub-sequences have a typical effect on RNA polymerase (RNAP) speed: we show that nucleotide 5-mers tend to have typical RNAP speed (or transcription rate), which is consistent along different parts of genes and among different groups of genes with high correlation. We also demonstrate that relative RNAP speed correlates with mRNA levels of endogenous and heterologous genes. Furthermore, we show that the estimated transcription and translation elongation rates correlate in endogenous genes. Finally, we demonstrate that our results are consistent for different high resolution experimental measurements of RNAP densities. These results suggest for the first time that transcription elongation is partly encoded in the transcript, affected by the codon-usage, and optimized by evolution with a significant effect on gene expression and organismal fitness.

Crystal structure of the overlapping dinucleosome composed of hexasome and octasome

Nucleosomes are dynamic entities that are repositioned along DNA by chromatin remodeling processes. A nucleosome repositioned by the switch-sucrose nonfermentable (SWI/SNF) remodeler collides with a neighbor and forms the intermediate "overlapping dinucleosome." Here, we report the crystal structure of the overlapping dinucleosome, in which two nucleosomes are associated, at 3.14-angstrom resolution. In the overlapping dinucleosome structure, the unusual "hexasome" nucleosome, composed of the histone hexamer lacking one H2A-H2B dimer from the conventional histone octamer, contacts the canonical "octasome" nucleosome, and they intimately associate. Consequently, about 250 base pairs of DNA are left-handedly wrapped in three turns, without a linker DNA segment between the hexasome and octasome moieties. The overlapping dinucleosome structure may provide important information to understand how nucleosome repositioning occurs during the chromatin remodeling process.

Site-Specific Disulfide Crosslinked Nucleosomes with Enhanced Stability

We engineered nucleosome core particles (NCPs) with two site-specific cysteine crosslinks that increase the stability of the particle. The first disulfide was introduced between the two copies of H2A via an H2A-N38C point mutation, effectively crosslinking the two H2A/H2B heterodimers together to stabilize the histone octamer against H2A/H2B dimer dissociation. The second crosslink was engineered between an R40C point mutation on the N-terminal tail of H3 and the NCP DNA ends by the introduction of a convertible nucleotide. This crosslink maintains the nucleosome DNA in a fixed translational setting relative to the histone octamer and prevents dilution-driven dissociation. The X-ray crystal structures of NCPs containing the disulfides in isolation and in combination were determined. Both disulfides stabilize the structure of the NCP without disturbing the overall structure. Nucleosomes containing these modifications will be advantageous for biochemical and structural studies as a consequence of their greater resistance to dissociation during high dilution in purification, elevated salt for crystallization and vitrification for cryogenic electron microscopy.

•

Crosslinked nucleosome core particles have increased stability against H2A/H2B dimer loss and DNA dissociation.

•

A site-specific disulfide crosslink was introduced between the two copies of H2A in the histone octamer to stabilize its quaternary structure.

•

Site-specific disulfide crosslinks were introduced between histone H3 and DNA within the nucleosome core particle.

•

Three X-ray crystal structures of crosslinked nucleosome core particles were determined at high resolution.

Maximum entropy methods for extracting the learned features of deep neural networks

New architectures of multilayer artificial neural networks and new methods for training them are rapidly revolutionizing the application of machine learning in diverse fields, including business, social science, physical sciences, and biology. Interpreting deep neural networks, however, currently remains elusive, and a critical challenge lies in understanding which meaningful features a network is actually learning. We present a general method for interpreting deep neural networks and extracting network-learned features from input data. We describe our algorithm in the context of biological sequence analysis. Our approach, based on ideas from statistical physics, samples from the maximum entropy distribution over possible sequences, anchored at an input sequence and subject to constraints implied by the empirical function learned by a network. Using our framework, we demonstrate that local transcription factor binding motifs can be identified from a network trained on ChIP-seq data and that nucleosome positioning signals are indeed learned by a network trained on chemical cleavage nucleosome maps. Imposing a further constraint on the maximum entropy distribution also allows us to probe whether a network is learning global sequence features, such as the high GC content in nucleosome-rich regions. This work thus provides valuable mathematical tools for interpreting and extracting learned features from feed-forward neural networks.

Deep learning is a state-of-the-art reformulation of artificial neural networks that have a long history of development. It can perform superbly well in diverse automated classification and prediction problems, including handwriting recognition, image identification, and biological pattern recognition. Its modern success can be attributed to improved training algorithms, clever network architecture, rapid explosion of available data, and advanced computing power–all of which have allowed the great expansion in the number of unknown parameters to be estimated by the model. These parameters, however, are so intricately connected through highly nonlinear functions that interpreting which essential features of given data are actually used by a deep neural network for its excellent performance has been difficult. We address this problem by using ideas from statistical physics to sample new unseen data that are likely to behave similarly to original data points when passed through the trained network. This synthetic data cloud around each original data point retains informative features while averaging out nonessential ones, ultimately allowing us to extract important network-learned features from the original data set and thus improving the human interpretability of deep learning methods. We demonstrate how our method can be applied to biological sequence analysis.

Parallel mapping with site-directed hydroxyl radicals and micrococcal nuclease reveals structural features of positioned nucleosomes in vivo

Micrococcal nuclease (MNase) has been widely used for analyses of nucleosome locations in many organisms. However, due to its sequence preference, the interpretations of the positions and occupancies of nucleosomes using MNase have remained controversial. Next-generation sequencing (NGS) has also been utilized for analyses of MNase-digests, but some technical biases are commonly present in the NGS experiments. Here, we established a gel-based method to map nucleosome positions in Saccharomyces cerevisiae, using isolated nuclei as the substrate for the histone H4 S47C-site-directed chemical cleavage in parallel with MNase digestion. The parallel mapping allowed us to compare the chemically and enzymatically cleaved sites by indirect end-labeling and primer extension mapping, and thus we could determine the nucleosome positions and the sizes of the nucleosome-free regions (or nucleosome-depleted regions) more accurately, as compared to nucleosome mapping by MNase alone. The analysis also revealed that the structural features of the nucleosomes flanked by the nucleosome-free region were different from those within regularly arrayed nucleosomes, showing that the structures and dynamics of individual nucleosomes strongly depend on their locations. Moreover, we demonstrated that the parallel mapping results were generally consistent with the previous genome-wide chemical mapping and MNase-Seq results. Thus, the gel-based parallel mapping will be useful for the analysis of a specific locus under various conditions.

A Reply to "MNase-Sensitive Complexes in Yeast: Nucleosomes and Non-histone Barriers," by Chereji et al

In this issue of Molecular Cell, Chereji et al. (2017) present new data on MNase-sensitive particles previously identified upstream of transcription start sites at many promoters in budding yeast, and they argue, based upon negative histone-ChIP results, that they are non-nucleosomal signals generated by transcription factors (TFs). We show instead, based upon functional experiments where the relevant TFs are rapidly depleted, that this explanation does not hold, and we argue instead that histone ChIP and chemical cleavage assays have a limited capacity to capture these highly dynamic, MNase-sensitive "fragile" nucleosomes.

Enabling Graded and Large-Scale Multiplex of Desired Genes Using a Dual-Mode dCas9 Activator in Saccharomyces cerevisiae

Standard approaches for dCas9-based modification of gene expression are limited in the ability to multiplex targets, establish streamlined cassettes, and utilize commonly studied Pol II promoters. In this work, we repurpose the dCas9-VPR activator to act as a dual-mode activator/repressor that can be programmed solely on the basis of target position at gene loci. Furthermore, we implement this approach using a streamlined Pol II-ribozyme system that allows expression of many sgRNAs from a single transcript. By "stepping" dCas9-VPR within the promoter region and ORF we create graded activation and repression (respectively) of target genes, allowing precise control over multiplexed gene modulation. Expression from the Pol II system increased the net amount of sgRNA production in cells by 3.88-fold relative to the Pol III SNR52 promoter, leading to a significant improvement in dCas9-VPR repression strength. Finally, we utilize our Pol II system to create galactose-inducible switching of gene expression states and multiplex constructs capable of modulating up to 4 native genes from a single vector. Our approach represents a significant step toward minimizing DNA required to assemble CRISPR systems in eukaryotes while enhancing the efficacy (greater repression strength), scale (more sgRNAs), and scope (inducibility) of dCas9-mediated gene regulation.

Crowding-induced transcriptional bursts dictate polymerase and nucleosome density profiles along genes

During eukaryotic transcription, RNA polymerase (RNAP) translocates along DNA molecules covered with nucleosomes and other DNA binding proteins. Though the interactions between a single nucleosome and RNAP are by now fairly well understood, this understanding has not been synthesized into a description of transcription on crowded genes, where multiple RNAP transcribe through nucleosomes while preserving the nucleosome coverage. We here take a deductive modeling approach to establish the consequences of RNAP–nucleosome interactions for transcription in crowded environments. We show that under physiologically crowded conditions, the interactions of RNAP with nucleosomes induce a strong kinetic attraction between RNAP molecules, causing them to self-organize into stable and moving pelotons. The peloton formation quantitatively explains the observed nucleosome and RNAP depletion close to the initiation site on heavily transcribed genes. Pelotons further translate into short-timescale transcriptional bursts at termination, resulting in burst characteristics consistent with instances of bursty transcription observed in vivo. To facilitate experimental testing of our proposed mechanism, we present several analytic relations that make testable quantitative predictions.

INO80 represses osmostress induced gene expression by resetting promoter proximal nucleosomes

The conserved INO80 chromatin remodeling complex is involved in regulation of DNA damage repair, replication and transcription. It is commonly recruited to the transcription start region and contributes to the establishment of promoter-proximal nucleosomes. We find a substantial influence of INO80 on nucleosome dynamics and gene expression during stress induced transcription. Transcription induced by osmotic stress leads to genome-wide remodeling of promoter proximal nucleosomes. INO80 function is required for timely return of evicted nucleosomes to the 5΄ end of induced genes. Reduced INO80 function in Arp8-deficient cells leads to correlated prolonged transcription and nucleosome eviction. INO80 and the related complex SWR1 regulate incorporation of the H2A.Z isoform at promoter proximal nucleosomes. However, H2A.Z seems not to influence osmotic stress induced gene regulation. Furthermore, we show that high rates of transcription promote INO80 recruitment to promoter regions, suggesting a connection between active transcription and promoter proximal nucleosome remodeling. In addition, we find that absence of INO80 enhances bidirectional promoter activity at highly induced genes and expression of a number of stress induced transcripts. We suggest that INO80 has a direct repressive role via promoter proximal nucleosome remodeling to limit high levels of transcription in yeast.

DNA topology in chromatin is defined by nucleosome spacing

In eukaryotic nucleosomes, DNA makes ~1.7 superhelical turns around histone octamer. However, there is a long-standing discrepancy between the nucleosome core structure determined by x-ray crystallography and measurements of DNA topology in circular minichromosomes, indicating that there is only ~1.0 superhelical turn per nucleosome. Although several theoretical assumptions were put forward to explain this paradox by conformational variability of the nucleosome linker, none was tested experimentally. We analyzed topological properties of DNA in circular nucleosome arrays with precisely positioned nucleosomes. Using topological electrophoretic assays and electron microscopy, we demonstrate that the DNA linking number per nucleosome strongly depends on the nucleosome spacing and varies from -1.4 to -0.9. For the predominant {10n + 5} class of nucleosome repeats found in native chromatin, our results are consistent with the DNA topology observed earlier. Thus, we reconcile the topological properties of nucleosome arrays with nucleosome core structure and provide a simple explanation for the DNA topology in native chromatin with variable DNA linker length. Topological polymorphism of the chromatin fibers described here may reflect a more general tendency of chromosomal domains containing active or repressed genes to acquire different nucleosome spacing to retain topologically distinct higher-order structures.

Evolutionary analysis of nucleosome positioning sequences based on New Symmetric Relative Entropy

New Symmetric Relative Entropy (NSRE) was applied innovatively to analyze the nucleosome sequences in S. cerevisiae, S. pombe and Drosophila. NSRE distributions could well reflect the characteristic differences of nucleosome sequences among three organisms, and the differences indicate a concerted evolution in the sequence usage of nucleosome. Further analysis about the nucleosomes around TSS shows that the constitutive property of +1/-1 nucleosomes in S. cerevisiae is different from that in S. pombe and Drosophila, which indicates that S. cerevisiae has a different transcription regulation mechanism based on nucleosome. However, in either case, the nucleosome dyad region is conserved and always has a higher NSRE. Base composition analysis shows that this conservative property in nucleosome dyad region is mainly determined by base A and T, and the dependence degrees on base A and T are consistent in three organisms.

Genome-wide maps of alkylation damage, repair, and mutagenesis in yeast reveal mechanisms of mutational heterogeneity

DNA base damage is an important contributor to genome instability, but how the formation and repair of these lesions is affected by the genomic landscape and contributes to mutagenesis is unknown. Here, we describe genome-wide maps of DNA base damage, repair, and mutagenesis at single nucleotide resolution in yeast treated with the alkylating agent methyl methanesulfonate (MMS). Analysis of these maps revealed that base excision repair (BER) of alkylation damage is significantly modulated by chromatin, with faster repair in nucleosome-depleted regions, and slower repair and higher mutation density within strongly positioned nucleosomes. Both the translational and rotational settings of lesions within nucleosomes significantly influence BER efficiency; moreover, this effect is asymmetric relative to the nucleosome dyad axis and is regulated by histone modifications. Our data also indicate that MMS-induced mutations at adenine nucleotides are significantly enriched on the nontranscribed strand (NTS) of yeast genes, particularly in BER-deficient strains, due to higher damage formation on the NTS and transcription-coupled repair of the transcribed strand (TS). These findings reveal the influence of chromatin on repair and mutagenesis of base lesions on a genome-wide scale and suggest a novel mechanism for transcription-associated mutation asymmetry, which is frequently observed in human cancers.

MNase-Sensitive Complexes in Yeast: Nucleosomes and Non-histone Barriers

Micrococcal nuclease (MNase) is commonly used to map nucleosomes genome-wide, but nucleosome maps are affected by the degree of digestion. It has been proposed that many yeast promoters are not nucleosome-free but instead occupied by easily digested, unstable, "fragile" nucleosomes. We analyzed the histone content of all MNase-sensitive complexes by MNase-ChIP-seq and sonication-ChIP-seq. We find that yeast promoters are predominantly bound by non-histone protein complexes, with little evidence for fragile nucleosomes. We do detect MNase-sensitive nucleosomes elsewhere in the genome, including at transcription termination sites. However, they have high A/T content, suggesting that MNase sensitivity does not indicate instability, but rather the preference of MNase for A/T-rich DNA, such that A/T-rich nucleosomes are digested faster than G/C-rich nucleosomes. We confirm our observations by analyzing ChIP-exo, chemical mapping, and ATAC-seq data from other laboratories. Thus, histone ChIP-seq experiments are essential to distinguish nucleosomes from other DNA-binding proteins that protect against MNase.

Towards quantitative analysis of gene regulation by enhancers

Enhancers are regulatory DNA sequences that can activate transcription over large distances. Recent studies have revealed the widespread role of distant activation in eukaryotic gene regulation and in the development of various human diseases, including cancer. Here we review recent progress in the field, focusing on new experimental and computational approaches that quantify the role of chromatin structure and dynamics during enhancer-promoter interactions in vitro and in vivo.

Histone Acetylation, Not Stoichiometry, Regulates Linker Histone Binding in Saccharomyces cerevisiae

Linker histones play a fundamental role in shaping chromatin structure, but how their interaction with chromatin is regulated is not well understood. In this study, we used a combination of genetic and genomic approaches to explore the regulation of linker histone binding in the yeast, Saccharomyces cerevisiae We found that increased expression of Hho1, the yeast linker histone, resulted in a severe growth defect, despite only subtle changes in chromatin structure. Further, this growth defect was rescued by mutations that increase histone acetylation. Consistent with this, genome-wide analysis of linker histone occupancy revealed an inverse correlation with histone tail acetylation in both yeast and mouse embryonic stem cells. Collectively, these results suggest that histone acetylation negatively regulates linker histone binding in S. cerevisiae and other organisms and provide important insight into how chromatin structure is regulated and maintained to both facilitate and repress transcription.

Genome-wide Mapping of the Nucleosome Landscape by Micrococcal Nuclease and Chemical Mapping

Nucleosomes regulate the transcription output of the genome by occluding the underlying DNA sequences from DNA-binding proteins that must act on it. Knowledge of the precise locations of nucleosomes in the genome is thus essential towards understanding how transcription is regulated. Current nucleosome-mapping strategies involve digesting chromatin with nucleases or chemical cleavage followed by high-throughput sequencing. In this review, we compare the traditional micrococcal nuclease (MNase)-based approach with a chemical cleavage strategy, with discussion on the important insights each has uncovered about the role of nucleosomes in shaping transcriptional processes.

Genome-wide Nucleosome Occupancy and Organization Modulates the Plasticity of Gene Transcriptional Status in Maize

Nucleosomes are fundamental units of chromatin that play critical roles in gene regulation by modulating DNA accessibility. However, their roles in regulating tissue-specific gene transcription are poorly understood. Here, we present genome-wide nucleosome maps of maize shoot and endosperm generated by sequencing the micrococcal nuclease digested nucleosomal DNA. The changes of gene transcriptional status between shoot and endosperm were accompanied by preferential nucleosome loss from the promoters and shifts in the first nucleosome downstream of the transcriptional start sites (+1 nucleosome) and upstream of transcriptional termination sites (-1 nucleosome). Intrinsically DNA-encoded nucleosome organization was largely associated with the capacity of a gene to alter its transcriptional status among different tissues. Compared with tissue-specific genes, constitutively expressed genes showed more pronounced 5' and 3' nucleosome-depleted regions as well as further +1 nucleosome to transcriptional start sites and -1 nucleosome to transcriptional termination sites. Moreover, nucleosome organization was more highly correlated with the plasticity of gene transcriptional status than with its expression level when examined using in vivo and predicted nucleosome data. In addition, the translational efficiencies of tissue-specific genes appeared to be greater than those of constitutively expressed genes. Taken together, our results indicate that intrinsically DNA-encoded nucleosome organization is important, beyond its role in regulating gene expression levels, in determining the plasticity of gene transcriptional status.

A comparison of nucleosome organization in Drosophila cell lines

Changes in the distribution of nucleosomes along the genome influence chromatin structure and impact gene expression by modulating the accessibility of DNA to transcriptional machinery. However, the role of genome-wide nucleosome positioning in gene expression and in maintaining differentiated cell states remains poorly understood. Drosophila melanogaster cell lines represent distinct tissue types and exhibit cell-type specific gene expression profiles. They thus could provide a useful tool for investigating cell-type specific nucleosome organization of an organism's genome. To evaluate this possibility, we compared genome-wide nucleosome positioning and occupancy in five different Drosophila tissue-specific cell lines, and in reconstituted chromatin, and then tested for correlations between nucleosome positioning, transcription factor binding motifs, and gene expression. Nucleosomes in all cell lines were positioned in accordance with previously known DNA-nucleosome interactions, with helically repeating A/T di-nucleotide pairs arranged within nucleosomal DNAs and AT-rich pentamers generally excluded from nucleosomal DNA. Nucleosome organization in all cell lines differed markedly from in vitro reconstituted chromatin, with highly expressed genes showing strong nucleosome organization around transcriptional start sites. Importantly, comparative analysis identified genomic regions that exhibited cell line-specific nucleosome enrichment or depletion. Further analysis of these regions identified 91 out of 16,384 possible heptamer sequences that showed differential nucleosomal occupation between cell lines, and 49 of the heptamers matched one or more known transcription factor binding sites. These results demonstrate that there is differential nucleosome positioning between these Drosophila cell lines and therefore identify a system that could be used to investigate the functional significance of differential nucleosomal positioning in cell type specification.

Understanding nucleosome dynamics and their links to gene expression and DNA replication

Advances in genomics technology have provided the means to probe myriad chromatin interactions at unprecedented spatial and temporal resolution. This has led to a profound understanding of nucleosome organization within the genome, revealing that nucleosomes are highly dynamic. Nucleosome dynamics are governed by a complex interplay of histone composition, histone post-translational modifications, nucleosome occupancy and positioning within chromatin, which are influenced by numerous regulatory factors, including general regulatory factors, chromatin remodellers, chaperones and polymerases. It is now known that these dynamics regulate diverse cellular processes ranging from gene transcription to DNA replication and repair.

The proto-chromatosome: A fundamental subunit of chromatin?

Eukaryotic DNA is packaged into regularly spaced nucleosomes, resembling beads on a string. Each bead contains ∼147 bp wrapped around a core histone octamer. Linker histone (H1) binds to the linker DNA to drive chromatin folding. Micrococcal nuclease (MNase) digestion studies reveal 2 mono-nucleosomal intermediates: the core particle (∼147 bp) and the chromatosome (∼160 bp; a core particle with additional DNA protected by H1). We have recently developed an improved method for mapping nucleosomes, using exonuclease III to remove residual linker (MNase-Exo-seq). (1) We discovered 2 new intermediate particles corresponding to core particles with ∼7 bp of linker protruding from one side (∼154 bp) or both sides (∼161 bp), which are formed in the absence of H1. We propose that these "proto-chromatosomes" are stabilized by core histone-DNA contacts in the linker, ∼7 bp from the nucleosome boundaries. These contacts may determine the topography of the H1 binding site.

Designing nucleosomal force sensors

About three quarters of our DNA is wrapped into nucleosomes: DNA spools with a protein core. It is well known that the affinity of a given DNA stretch to be incorporated into a nucleosome depends on the geometry and elasticity of the basepair sequence involved, causing the positioning of nucleosomes. Here we show that DNA elasticity can have a much deeper effect on nucleosomes than just their positioning: it affects their "identities". Employing a recently developed computational algorithm, the mutation Monte Carlo method, we design nucleosomes with surprising physical characteristics. Unlike any other nucleosomes studied so far, these nucleosomes are short-lived when put under mechanical tension whereas other physical properties are largely unaffected. This suggests that the nucleosome, the most abundant DNA-protein complex in our cells, might more properly be considered a class of complexes with a wide array of physical properties, and raises the possibility that evolution has shaped various nucleosome species according to their genomic context.

Kilobase Pair Chromatin Fiber Contacts Promoted by Living-System-Like DNA Linker Length Distributions and Nucleosome Depletion

Nucleosome placement, or DNA linker length patterns, are believed to yield specific spatial features in chromatin fibers, but details are unknown. Here we examine by mesoscale modeling how kilobase (kb) range contacts and fiber looping depend on linker lengths ranging from 18 to 45 bp, with values modeled after living systems, including nucleosome free regions (NFRs) and gene encoding segments. We also compare artificial constructs with alternating versus randomly distributed linker lengths in the range of 18-72 bp. We show that nonuniform distributions with NFRs enhance flexibility and encourage kb-range contacts. NFRs between neighboring gene segments diminish short-range contacts between flanking nucleosomes, while enhancing kb-range contacts via hierarchical looping. We also demonstrate that variances in linker lengths enhance such contacts. In particular, moderate sized variations in fiber linker lengths (∼27 bp) encourage long-range contacts in randomly distributed linker length fibers. Our work underscores the importance of linker length patterns, alongside bound proteins, in biological regulation. Contacts formed by kb-range chromatin folding are crucial to gene activity. Because we find that special linker length distributions in living systems promote kb contacts, our work suggests ways to manipulate these patterns for regulation of gene activity.

Establishing nucleosome architecture and stability at promoters: Roles of pioneer transcription factors and the RSC chromatin remodeler

Improvements in deep sequencing, together with methods to rapidly deplete essential transcription factors (TFs) and chromatin remodelers, have recently led to a more detailed picture of promoter nucleosome architecture in yeast and its relationship to transcriptional regulation. These studies revealed that ∼40% of all budding yeast protein-coding genes possess a unique promoter structure, where we propose that an unusually unstable nucleosome forms immediately upstream of the transcription start site (TSS). This "fragile" nucleosome (FN) promoter architecture relies on the combined action of the essential RSC (Remodels Structure of Chromatin) nucleosome remodeler and pioneer transcription factors (PTFs). FNs are associated with genes whose expression is high, coupled to cell growth, and characterized by low cell-to-cell variability (noise), suggesting that they may promote these features. Recent studies in metazoans suggest that the presence of dynamic nucleosomes upstream of the TSS at highly expressed genes may be conserved throughout evolution.

Well-positioned nucleosomes punctuate polycistronic pol II transcription units and flank silent VSG gene arrays in Trypanosoma brucei

Background

The compaction of DNA in chromatin in eukaryotes allowed the expansion of genome size and coincided with significant evolutionary diversification. However, chromatin generally represses DNA function, and mechanisms coevolved to regulate chromatin structure and its impact on DNA. This included the selection of specific nucleosome positions to modulate accessibility to the DNA molecule. Trypanosoma brucei, a member of the Excavates supergroup, falls in an ancient evolutionary branch of eukaryotes and provides valuable insight into the organization of chromatin in early genomes.

Results

We have mapped nucleosome positions in T. brucei and identified important differences compared to other eukaryotes: The RNA polymerase II initiation regions in T. brucei do not exhibit pronounced nucleosome depletion, and show little evidence for defined −1 and +1 nucleosomes. In contrast, a well-positioned nucleosome is present directly on the splice acceptor sites within the polycistronic transcription units. The RNA polyadenylation sites were depleted of nucleosomes, with a single well-positioned nucleosome present immediately downstream of the predicted sites. The regions flanking the silent variant surface glycoprotein (VSG) gene cassettes showed extensive arrays of well-positioned nucleosomes, which may repress cryptic transcription initiation. The silent VSG genes themselves exhibited a less regular nucleosomal pattern in both bloodstream and procyclic form trypanosomes. The DNA replication origins, when present within silent VSG gene cassettes, displayed a defined nucleosomal organization compared with replication origins in other chromosomal core regions.

Conclusions

Our results indicate that some organizational features of chromatin are evolutionarily ancient, and may already have been present in the last eukaryotic common ancestor.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-017-0121-9) contains supplementary material, which is available to authorized users.

Sequence-specific thermodynamic properties of nucleic acids influence both transcriptional pausing and backtracking in yeast

RNA Polymerase II pauses and backtracks during transcription, with many consequences for gene expression and cellular physiology. Here, we show that the energy required to melt double-stranded nucleic acids in the transcription bubble predicts pausing in Saccharomyces cerevisiae far more accurately than nucleosome roadblocks do. In addition, the same energy difference also determines when the RNA polymerase backtracks instead of continuing to move forward. This data-driven model corroborates—in a genome wide and quantitative manner—previous evidence that sequence-dependent thermodynamic features of nucleic acids influence both transcriptional pausing and backtracking.

NucTools: analysis of chromatin feature occupancy profiles from high-throughput sequencing data

Background

Biomedical applications of high-throughput sequencing methods generate a vast amount of data in which numerous chromatin features are mapped along the genome. The results are frequently analysed by creating binary data sets that link the presence/absence of a given feature to specific genomic loci. However, the nucleosome occupancy or chromatin accessibility landscape is essentially continuous. It is currently a challenge in the field to cope with continuous distributions of deep sequencing chromatin readouts and to integrate the different types of discrete chromatin features to reveal linkages between them.

Results

Here we introduce the NucTools suite of Perl scripts as well as MATLAB- and R-based visualization programs for a nucleosome-centred downstream analysis of deep sequencing data. NucTools accounts for the continuous distribution of nucleosome occupancy. It allows calculations of nucleosome occupancy profiles averaged over several replicates, comparisons of nucleosome occupancy landscapes between different experimental conditions, and the estimation of the changes of integral chromatin properties such as the nucleosome repeat length. Furthermore, NucTools facilitates the annotation of nucleosome occupancy with other chromatin features like binding of transcription factors or architectural proteins, and epigenetic marks like histone modifications or DNA methylation. The applications of NucTools are demonstrated for the comparison of several datasets for nucleosome occupancy in mouse embryonic stem cells (ESCs) and mouse embryonic fibroblasts (MEFs).

Conclusions

The typical workflows of data processing and integrative analysis with NucTools reveal information on the interplay of nucleosome positioning with other features such as for example binding of a transcription factor CTCF, regions with stable and unstable nucleosomes, and domains of large organized chromatin K9me2 modifications (LOCKs). As potential limitations and problems we discuss how inter-replicate variability of MNase-seq experiments can be addressed.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3580-2) contains supplementary material, which is available to authorized users.

Evolutionary mechanism and biological functions of 8-mers containing CG dinucleotide in yeast

The rules of k-mer non-random usage and the biological functions are worthy of special attention. Firstly, the article studied human 8-mer spectra and found that only the spectra of cytosine-guanine (CG) dinucleotide classification formed independent unimodal distributions when the 8-mers were classified into three subsets under 16 dinucleotide classifications. Secondly, the distribution rules were reproduced by other seven species including yeast, which showed that the evolution phenomenon had species universality. It followed that we proposed two theoretical conjectures: (1) CG₁ motifs (8-mers including 1 CG) are the nucleosome-binding motifs. (2) CG₂ motifs (8-mers including two or more than two CG) are the modular units of CpG islands. Our conjectures were confirmed in yeast by the following results: a maximum of average area under the receiver operating characteristic (AUC) resulted from CG₁ information during nucleosome core sequences, and linker sequences were distinguished by three CG subsets; there was a one-to-one relationship between abundant CG₁ signal regions and histone positions; the sequence changing of squeezed nucleosomes was relevant with the strength of CG₁ signals; and the AUC value of 0.986 was based on CG₂ information when CpG islands and non-CpG islands were distinguished by the three CG subsets.

Paf1 Has Distinct Roles in Transcription Elongation and Differential Transcript Fate

RNA polymerase II (Pol2) movement through chromatin and the co-transcriptional processing and fate of nascent transcripts is coordinated by transcription elongation factors (TEFs) such as polymerase-associated factor 1 (Paf1), but it is not known whether TEFs have gene-specific functions. Using strand-specific nucleotide resolution techniques, we show that levels of Paf1 on Pol2 vary between genes, are controlled dynamically by environmental factors via promoters, and reflect levels of processing and export factors on the encoded transcript. High levels of Paf1 on Pol2 promote transcript nuclear export, whereas low levels reflect nuclear retention. Strains lacking Paf1 show marked elongation defects, although low levels of Paf1 on Pol2 are sufficient for transcription elongation. Our findings support distinct Paf1 functions: a core general function in transcription elongation, satisfied by the lowest Paf1 levels, and a regulatory function in determining differential transcript fate by varying the level of Paf1 on Pol2.

Linking Chromatin Fibers to Gene Folding by Hierarchical Looping

While much is known about DNA structure on the basepair level, this scale represents only a fraction of the structural levels involved in folding the genomic material. With recent advances in experimental and theoretical techniques, a variety of structures have been observed on the fiber, gene, and chromosome levels of genome organization. Here we view chromatin architecture from nucleosomes and fibers to genes and chromosomes, highlighting the rich structural diversity and fiber fluidity emerging from both experimental and theoretical techniques. In this context, we discuss our recently proposed folding mechanism, which we call "hierarchical looping", similar to rope flaking used in mountain climbing, where 10-nm zigzag chromatin fibers are compacted laterally into self-associating loops which then stack and fold in space. We propose that hierarchical looping may act as a bridge between fibers and genes as well as provide a mechanism to relate key features of interphase and metaphase chromosome architecture to genome structural changes. This motif emerged by analysis of ultrastructural internucleosome contact data by electron microscopy-assisted nucleosome interaction capture cross-linking experiments, in combination with mesoscale modeling. We suggest that while the local folding of chromatin can be regulated at the fiber level by adjustment of internal factors such as linker-histone binding affinities, linker DNA lengths, and divalent ion levels, hierarchical looping on the gene level can additionally be controlled by posttranslational modifications and external factors such as polycomb group proteins. From a combination of 3C data and mesoscale modeling, we suggest that hierarchical looping could also play a role in epigenetic gene silencing, as stacked loops may occlude access to transcription start sites. With advances in crystallography, single-molecule in vitro biochemistry, in vivo imaging techniques, and genome-wide contact data experiments, various modeling approaches are allowing for previously unavailable structural interpretation of these data at multiple spatial and temporal scales. An unprecedented level of productivity and opportunity is on the horizon for the chromatin structure field.

Histone H3K4 and H3K36 Methylation Independently Recruit the NuA3 Histone Acetyltransferase in Saccharomyces cerevisiae

Histone post-translational modifications (PTMs) alter chromatin structure by promoting the interaction of chromatin-modifying complexes with nucleosomes. The majority of chromatin-modifying complexes contain multiple domains that preferentially interact with modified histones, leading to speculation that these domains function in concert to target nucleosomes with distinct combinations of histone PTMs. In Saccharomyces cerevisiae, the NuA3 histone acetyltransferase complex contains three domains, the PHD finger in Yng1, the PWWP domain in Pdp3, and the YEATS domain in Taf14; which in vitro bind to H3K4 methylation, H3K36 methylation, and acetylated and crotonylated H3K9, respectively. While the in vitro binding has been well characterized, the relative in vivo contributions of these histone PTMs in targeting NuA3 is unknown. Here, through genome-wide colocalization and by mutational interrogation, we demonstrate that the PHD finger of Yng1, and the PWWP domain of Pdp3 independently target NuA3 to H3K4 and H3K36 methylated chromatin, respectively. In contrast, we find no evidence to support the YEATS domain of Taf14 functioning in NuA3 recruitment. Collectively our results suggest that the presence of multiple histone PTM binding domains within NuA3, rather than restricting it to nucleosomes containing distinct combinations of histone PTMs, can serve to increase the range of nucleosomes bound by the complex. Interestingly, however, the simple presence of NuA3 is insufficient to ensure acetylation of the associated nucleosomes, suggesting a secondary level of acetylation regulation that does not involve control of HAT-nucleosome interactions.

Plasmodium falciparum Nucleosomes Exhibit Reduced Stability and Lost Sequence Dependent Nucleosome Positioning

The packaging and organization of genomic DNA into chromatin represents an additional regulatory layer of gene expression, with specific nucleosome positions that restrict the accessibility of regulatory DNA elements. The mechanisms that position nucleosomes in vivo are thought to depend on the biophysical properties of the histones, sequence patterns, like phased di-nucleotide repeats and the architecture of the histone octamer that folds DNA in 1.65 tight turns. Comparative studies of human and P. falciparum histones reveal that the latter have a strongly reduced ability to recognize internal sequence dependent nucleosome positioning signals. In contrast, the nucleosomes are positioned by AT-repeat sequences flanking nucleosomes in vivo and in vitro. Further, the strong sequence variations in the plasmodium histones, compared to other mammalian histones, do not present adaptations to its AT-rich genome. Human and parasite histones bind with higher affinity to GC-rich DNA and with lower affinity to AT-rich DNA. However, the plasmodium nucleosomes are overall less stable, with increased temperature induced mobility, decreased salt stability of the histones H2A and H2B and considerable reduced binding affinity to GC-rich DNA, as compared with the human nucleosomes. In addition, we show that plasmodium histone octamers form the shortest known nucleosome repeat length (155bp) in vitro and in vivo. Our data suggest that the biochemical properties of the parasite histones are distinct from the typical characteristics of other eukaryotic histones and these properties reflect the increased accessibility of the P. falciparum genome.

Nucleosomes are not positioned randomly on DNA but on preferential sites with respect to the underlying DNA sequence. Histones belong to the most conserved eukaryotic proteins, as sequence dependent nucleosome positioning is an essential regulatory feature of nucleosomes, determining the accessibility of regulatory factors to DNA. We determined the biochemical properties of plasmodium histones and show that they are distinct from human forms, explaining the accessible chromatin structure of P. falciparum. Amino acid exchanges in the histones do not present an adaption to the AT-rich genome, but rather reduce the binding affinity to GC-rich DNA sequences, resulting in rather unstable nucleosomes with labile H2A and H2B, requiring extra-nucleosomal positioning signals to keep them on place. Plasmodium chromatin exhibits the shortest nucleosome spacing known to date potentially inhibiting the formation of higher order structures and maintaining chromatin accessible.

Contributions of Sequence to the Higher-Order Structures of DNA

One of the critical unanswered questions in genome biophysics is how the primary sequence of DNA bases influences the global properties of very-long-chain molecules. The local sequence-dependent features of DNA found in high-resolution structures introduce irregularities in the disposition of adjacent residues that facilitate the specific binding of proteins and modulate the global folding and interactions of double helices with hundreds of basepairs. These features also determine the positions of nucleosomes on DNA and the lengths of the interspersed DNA linkers. Like the patterns of basepair association within DNA, the arrangements of nucleosomes in chromatin modulate the properties of longer polymers. The intrachromosomal loops detected in genomic studies contain hundreds of nucleosomes, and given that the simulated configurations of chromatin depend on the lengths of linker DNA, the formation of these loops may reflect sequence-dependent information encoded within the positioning of the nucleosomes. With knowledge of the positions of nucleosomes on a given genome, methods are now at hand to estimate the looping propensities of chromatin in terms of the spacing of nucleosomes and to make a direct connection between the DNA base sequence and larger-scale chromatin folding.

Dissecting relative contributions of cis- and trans-determinants to nucleosome distribution by comparing Tetrahymena macronuclear and micronuclear chromatin

The ciliate protozoan Tetrahymena thermophila contains two types of structurally and functionally differentiated nuclei: the transcriptionally active somatic macronucleus (MAC) and the transcriptionally silent germ-line micronucleus (MIC). Here, we demonstrate that MAC features well-positioned nucleosomes downstream of transcription start sites and flanking splice sites. Transcription-associated trans-determinants promote nucleosome positioning in MAC. By contrast, nucleosomes in MIC are dramatically delocalized. Nucleosome occupancy in MAC and MIC are nonetheless highly correlated with each other, as well as with in vitro reconstitution and predictions based upon DNA sequence features, revealing unexpectedly strong contributions from cis-determinants. In particular, well-positioned nucleosomes are often matched with GC content oscillations. As many nucleosomes are coordinately accommodated by both cis- and trans-determinants, we propose that their distribution is shaped by the impact of these nucleosomes on the mutational and transcriptional landscape, and driven by evolutionary selection.