Nucleosome positions predicted through comparative genomics

Nucleosomal DNA has topological memory

One elusive aspect of the chromosome architecture is how it constrains the DNA topology. Nucleosomes stabilise negative DNA supercoils by restraining a DNA linking number difference (∆Lk) of about -1.26. However, whether this capacity is uniform across the genome is unknown. Here, we calculate the ∆Lk restrained by over 4000 nucleosomes in yeast cells. To achieve this, we insert each nucleosome in a circular minichromosome and perform Topo-seq, a high-throughput procedure to inspect the topology of circular DNA libraries in one gel electrophoresis. We show that nucleosomes inherently restrain distinct ∆Lk values depending on their genomic origin. Nucleosome DNA topologies differ at gene bodies (∆Lk = -1.29), intergenic regions (∆Lk = -1.23), rDNA genes (∆Lk = -1.24) and telomeric regions (∆Lk = -1.07). Nucleosomes near the transcription start and termination sites also exhibit singular DNA topologies. Our findings demonstrate that nucleosome DNA topology is imprinted by its native chromatin context and persists when the nucleosome is relocated.

Emerging Approaches to Profile Accessible Chromatin from Formalin-Fixed Paraffin-Embedded Sections

Nucleosomes are non-uniformly distributed across eukaryotic genomes, with stretches of 'open' chromatin strongly associated with transcriptionally active promoters and enhancers. Understanding chromatin accessibility patterns in normal tissue and how they are altered in pathologies can provide critical insights to development and disease. With the advent of high-throughput sequencing, a variety of strategies have been devised to identify open regions across the genome, including DNase-seq, MNase-seq, FAIRE-seq, ATAC-seq, and NicE-seq. However, the broad application of such methods to FFPE (formalin-fixed paraffin-embedded) tissues has been curtailed by the major technical challenges imposed by highly fixed and often damaged genomic material. Here, we review the most common approaches for mapping open chromatin regions, recent optimizations to overcome the challenges of working with FFPE tissue, and a brief overview of a typical data pipeline with analysis considerations.

Genome organization across scales: mechanistic insights from in vitro reconstitution studies

Eukaryotic genomes are compacted and organized into distinct three-dimensional (3D) structures, which range from small-scale nucleosome arrays to large-scale chromatin domains. These chromatin structures play an important role in the regulation of transcription and other nuclear processes. The molecular mechanisms that drive the formation of chromatin structures across scales and the relationship between chromatin structure and function remain incompletely understood. Because the processes involved are complex and interconnected, it is often challenging to dissect the underlying principles in the nuclear environment. Therefore, in vitro reconstitution systems provide a valuable approach to gain insight into the molecular mechanisms by which chromatin structures are formed and to determine the cause-consequence relationships between the processes involved. In this review, we give an overview of in vitro approaches that have been used to study chromatin structures across scales and how they have increased our understanding of the formation and function of these structures. We start by discussing in vitro studies that have given insight into the mechanisms of nucleosome positioning. Next, we discuss recent efforts to reconstitute larger-scale chromatin domains and loops and the resulting insights into the principles of genome organization. We conclude with an outlook on potential future applications of chromatin reconstitution systems and how they may contribute to answering open questions concerning chromatin architecture.

A single fiber view of the nucleosome organization in eukaryotic chromatin

Eukaryotic cells are thought to arrange nucleosomes into extended arrays with evenly spaced nucleosomes phased at genomic landmarks. Here we tested to what extent this stereotypic organization describes the nucleosome organization in Saccharomyces cerevisiae using Fiber-Seq, a long-read sequencing technique that maps entire nucleosome arrays on individual chromatin fibers in a high throughput manner. With each fiber coming from a different cell, Fiber-Seq uncovers cell-to-cell heterogeneity. The long reads reveal the nucleosome architecture even over repetitive DNA such as the ribosomal DNA repeats. The absolute nucleosome occupancy, a parameter that is difficult to obtain with conventional sequencing approaches, is a direct readout of Fiber-Seq. We document substantial deviations from the stereotypical nucleosome organization with unexpectedly long linker DNAs between nucleosomes, gene bodies missing entire nucleosomes, cell-to-cell heterogeneity in nucleosome occupancy, heterogeneous phasing of arrays and irregular nucleosome spacing. Nucleosome array structures are indistinguishable throughout the gene body and with respect to the direction of transcription arguing against transcription promoting array formation. Acute nucleosome depletion destroyed most of the array organization indicating that nucleosome remodelers cannot efficiently pack nucleosomes under those conditions. Given that nucleosomes are cis-regulatory elements, the cell-to-cell heterogeneity uncovered by Fiber-Seq provides much needed information to understand chromatin structure and function.

Galaxy Dnpatterntools for Computational Analysis of Nucleosome Positioning Sequence Patterns

Nucleosomes are basic units of DNA packing in eukaryotes. Their structure is well conserved from yeast to human and consists of the histone octamer core and 147 bp DNA wrapped around it. Nucleosomes are bound to a majority of the eukaryotic genomic DNA, including its regulatory regions. Hence, they also play a major role in gene regulation. For the latter, their precise positioning on DNA is essential. In the present paper, we describe Galaxy dnpatterntools—software package for nucleosome DNA sequence analysis and mapping. This software will be useful for computational biologists practitioners to conduct more profound studies of gene regulatory mechanisms.

NucPosDB: a database of nucleosome positioning in vivo and nucleosomics of cell-free DNA

Nucleosome positioning is involved in many gene regulatory processes happening in the cell, and it may change as cells differentiate or respond to the changing microenvironment in a healthy or diseased organism. One important implication of nucleosome positioning in clinical epigenetics is its use in the “nucleosomics” analysis of cell-free DNA (cfDNA) for the purpose of patient diagnostics in liquid biopsies. The rationale for this is that the apoptotic nucleases that digest chromatin of the dying cells mostly cut DNA between nucleosomes. Thus, the short pieces of DNA in body fluids reflect the positions of nucleosomes in the cells of origin. Here, we report a systematic nucleosomics database — NucPosDB — curating published nucleosome positioning datasets in vivo as well as datasets of sequenced cell-free DNA (cfDNA) that reflect nucleosome positioning in situ in the cells of origin. Users can select subsets of the database by a number of criteria and then obtain raw or processed data. NucPosDB also reports the originally determined regions with stable nucleosome occupancy across several individuals with a given condition. An additional section provides a catalogue of computational tools for the analysis of nucleosome positioning or cfDNA experiments and theoretical algorithms for the prediction of nucleosome positioning preferences from DNA sequence. We provide an overview of the field, describe the structure of the database in this context, and demonstrate data variability using examples of different medical conditions. NucPosDB is useful both for the analysis of fundamental gene regulation processes and the training of computational models for patient diagnostics based on cfDNA. The database currently curates ~ 400 publications on nucleosome positioning in cell lines and in situ as well as cfDNA from > 10,000 patients and healthy volunteers. For open-access cfDNA datasets as well as key MNase-seq datasets in human cells, NucPosDB allows downloading processed mapped data in addition to the regions with stable nucleosome occupancy. NucPosDB is available at https://generegulation.org/nucposdb/.

Predicting the configuration and energy of DNA in a nucleosome by coarse-grain modelling

We present a novel algorithm which uses a coarse-grained model and an energy minimisation procedure to predict the sequence-dependent DNA configuration in a nucleosome together with its energetic cost.

The biogenesis and function of nucleosome arrays

Numerous chromatin remodeling enzymes position nucleosomes in eukaryotic cells. Aside from these factors, transcription, DNA sequence, and statistical positioning of nucleosomes also shape the nucleosome landscape. The precise contributions of these processes remain unclear due to their functional redundancy in vivo. By incisive genome engineering, we radically decreased their redundancy in Saccharomyces cerevisiae. The transcriptional machinery strongly disrupts evenly spaced nucleosomes. Proper nucleosome density and DNA sequence are critical for their biogenesis. The INO80 remodeling complex helps space nucleosomes in vivo and positions the first nucleosome over genes in an H2A.Z-independent fashion. INO80 requires its Arp8 subunit but unexpectedly not the Nhp10 module for spacing. Cells with irregularly spaced nucleosomes suffer from genotoxic stress including DNA damage, recombination and transpositions. We derive a model of the biogenesis of the nucleosome landscape and suggest that it evolved not only to regulate but also to protect the genome.

NucHMM: a method for quantitative modeling of nucleosome organization identifying functional nucleosome states distinctly associated with splicing potentiality

We develop a novel computational method, NucHMM, to identify functional nucleosome states associated with cell type-specific combinatorial histone marks and nucleosome organization features such as phasing, spacing and positioning. We test it on publicly available MNase-seq and ChIP-seq data in MCF7, H1, and IMR90 cells and identify 11 distinct functional nucleosome states. We demonstrate these nucleosome states are distinctly associated with the splicing potentiality of skipping exons. This advances our understanding of the chromatin function at the nucleosome level and offers insights into the interplay between nucleosome organization and splicing processes.

Effects of size, cooperativity, and competitive binding on protein positioning on DNA

Accurate positioning of proteins on chromosomal DNA is crucial for its proper organization as well as gene transcription regulation. Recent experiments revealed existence of periodic patterns of nucleoprotein complexes on DNA, which frequently cannot be explained by sequence-dependent binding of proteins. Previous theoretical studies suggest that such patterns typically emerge as a result of the proteins' volume-exclusion effect. However, the role of other physical factors in patterns' formation, such as the length of DNA, its sequence heterogeneity, and protein binding cooperativity/binding competition to DNA, remains unclear. To address these less understood yet important aspects, we investigated potential effects of these factors on protein positioning on finite-size DNA by using transfer-matrix calculations. It has been found that upon binding to DNA, proteins form oscillatory patterns that span over the length of up to ∼10 times the size of the protein binding site, with the shape of the patterns being strongly dependent on the length of DNA and the proteins' binding cooperativity to DNA. Furthermore, calculations showed that small variations in the proteins' affinity to DNA due to its sequence heterogeneity do not much change the main geometric characteristics of the observed protein patterns. Finally, competition between two different types of proteins for binding to DNA has been found to lead to formation of highly diverse and complex alternating positioning of the two proteins. Altogether, these results provide new insights into the roles of physicochemical properties of proteins, the DNA length, and DNA-binding competition between proteins in formation of protein positioning patterns on DNA.

Computational methods and next-generation sequencing approaches to analyze epigenetics data: Profiling of methods and applications

Epigenetics is mainly comprised of features that regulate genomic interactions thereby playing a crucial role in a vast array of biological processes. Epigenetic mechanisms such as DNA methylation and histone modifications influence gene expression by modulating the packaging of DNA in the nucleus. A plethora of studies have emphasized the importance of analyzing epigenetics data through genome-wide studies and high-throughput approaches, thereby providing key insights towards epigenetics-based diseases such as cancer. Recent advancements have been made towards translating epigenetics research into a high throughput approach such as genome-scale profiling. Amongst all, bioinformatics plays a pivotal role in achieving epigenetics-related computational studies. Despite significant advancements towards epigenomic profiling, it is challenging to understand how various epigenetic modifications such as chromatin modifications and DNA methylation regulate gene expression. Next-generation sequencing (NGS) provides accurate and parallel sequencing thereby allowing researchers to comprehend epigenomic profiling. In this review, we summarize different computational methods such as machine learning and other bioinformatics tools, publicly available databases and resources to identify key modifications associated with epigenetic machinery. Additionally, the review also focuses on understanding recent methodologies related to epigenome profiling using NGS methods ranging from library preparation, different sequencing platforms and analytical techniques to evaluate various epigenetic modifications such as DNA methylation and histone modifications. We also provide detailed information on bioinformatics tools and computational strategies responsible for analyzing large scale data in epigenetics.

Yeast chromatin remodeling complexes and their roles in transcription

The nucleosome is a small unit of chromatin, which is dynamic in eukaryotes. Chromatin conformation and post-translational modifications affect nucleosome dynamics under certain conditions, playing an important role in the epigenetic regulation of transcription, replication and reprogramming. The Snf2 remodeling family is one of the crucial remodeling complexes that tightly regulate chromatin structure and affect nucleosome dynamics. This family alters nucleosome positioning, exchanges histone variants, and assembles and disassembles nucleosomes at certain locations. Moreover, the Snf2 family, in conjunction with other co-factors, regulates gene expression in Saccharomyces cerevisiae. Here we first review recent findings on the Snf2 family remodeling complexes and then use some examples to illustrate the cooperation between different members of Snf2 family, and the cooperation between Snf2 family and other co-factors in gene regulation especially during transcription initiation.

Nucleosome positioning sequence patterns as packing or regulatory

Nucleosome positioning DNA sequence patterns (NPS)-usually distributions of particular dinucleotides or other sequence elements in nucleosomal DNA-at least partially determine chromatin structure and arrangements of nucleosomes that in turn affect gene expression. Statistically, NPS are defined as oscillations of the dinucleotide periodicity of about 10 base pairs (bp) which reflects the double helix period. We compared the nucleosomal DNA patterns in mouse, human and yeast organisms and observed few distinctive patterns that can be termed as packing and regulatory referring to distinctive modes of chromatin function. For the first time the NPS patterns in nucleus accumbens cells (NAC) in mouse brain were characterized and compared to the patterns in human CD4+ and apoptotic lymphocyte cells and well studied patterns in yeast. The NPS patterns in human CD4+ cells and mouse brain cells had very high positive correlation. However, there was no correlation between them and patterns in human apoptotic lymphocyte cells and yeast, but the latter two were highly correlated with each other. By their dinucleotide arrangements the analyzed NPS patterns classified into stable canonical WW/SS (W = A or T and S = C or G dinucleotide) and less stable RR/YY (R = A or G and Y = C or T dinucleotide) patterns and anti-patterns. In the anti-patterns positioning of the dinucleotides is flipped compared to those in the regular patterns. Stable canonical WW/SS patterns and anti-patterns are ubiquitously observed in many organisms and they had high resemblance between yeast and human apoptotic cells. Less stable RR/YY patterns had higher positive correlation between mouse and normal human cells. Our analysis and evidence from scientific literature lead to idea that various distinct patterns in nucleosomal DNA can be related to the two roles of the chromatin: packing (WW/SS) and regulatory (RR/YY and "anti").

Intracellular nucleosomes constrain a DNA linking number difference of -1.26 that reconciles the Lk paradox

The interplay between chromatin structure and DNA topology is a fundamental, yet elusive, regulator of genome activities. A paradigmatic case is the “linking number paradox” of nucleosomal DNA, which refers to the incongruence between the near two left-handed superhelical turns of DNA around the histone octamer and the DNA linking number difference (∆Lk) stabilized by individual nucleosomes, which has been experimentally estimated to be about −1.0. Here, we analyze the DNA topology of a library of mononucleosomes inserted into small circular minichromosomes to determine the average ∆Lk restrained by individual nucleosomes in vivo. Our results indicate that most nucleosomes stabilize about −1.26 units of ∆Lk. This value balances the twist (∆Tw  ≈ + 0.2) and writhe (∆Wr ≈ −1.5) deformations of nucleosomal DNA in terms of the equation ∆Lk = ∆Tw + ∆Wr. Our finding reconciles the existing discrepancy between theoretical and observed measurement of the ΔLk constrained by nucleosomes.

There had been an enduring discrepancy between theoretical and observed measurement of the DNA linking number (∆Lk) constrained by nucleosomes. Here the authors provide measurements of the ∆Lk constrained by individual nucleosomes in native chromatin that reconcile this discrepancy.

Structure and Chromosomal Organization of Yeast Genes Regulated by Topoisomerase II

Cellular DNA topoisomerases (topo I and topo II) are highly conserved enzymes that regulate the topology of DNA during normal genome transactions, such as DNA transcription and replication. In budding yeast, topo I is dispensable whereas topo II is essential, suggesting fundamental and exclusive roles for topo II, which might include the functions of the topo IIa and topo IIb isoforms found in mammalian cells. In this review, we discuss major findings of the structure and chromosomal organization of genes regulated by topo II in budding yeast. Experimental data was derived from short (10 min) and long term (120 min) responses to topo II inactivation in top-2 ts mutants. First, we discuss how short term responses reveal a subset of yeast genes that are regulated by topo II depending on their promoter architecture. These short term responses also uncovered topo II regulation of transcription across multi-gene clusters, plausibly by common DNA topology management. Finally, we examine the effects of deactivated topo II on the elongation of RNA transcripts. Each study provides an insight into the particular chromatin structure that interacts with the activity of topo II. These findings are of notable clinical interest as numerous anti-cancer therapies interfere with topo II activity.

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 molecular basis of metabolic cycles and their relationship to circadian rhythms

Metabolic cycles result from the partitioning of oxidative and reductive metabolism into rhythmic phases of gene expression and oscillating post-translational protein modifications. Relatively little is known about how these switches in gene expression are controlled, although recent studies have suggested that transcription itself may play a central role. This review explores the molecular basis of the metabolic and gene-expression oscillations in the yeast Saccharomyces cerevisiae, as well as how they relate to other biological time-keeping mechanisms, such as circadian rhythms.

Genomes of Multicellular Organisms Have Evolved to Attract Nucleosomes to Promoter Regions

Sequences that influence nucleosome positioning in promoter regions, and their relation to gene regulation, have been the topic of much research over the last decade. In yeast, significant nucleosome-depleted regions are found, which facilitate transcription. With the arrival of nucleosome positioning maps for the human genome, it was discovered that in our genome, unlike in that of yeast, promoters encode for high nucleosome occupancy. In this work, we look at the genomes of a range of different organisms, to provide a catalog of nucleosome positioning signals in promoters across the tree of life. We utilize a computational model of the nucleosome, based on crystallographic analyses of the structure and elasticity of the nucleosome, to predict the nucleosome positioning signals in promoter regions. To be able to apply our model to large genomic datasets, we introduce an approximative scheme that makes use of the limited range of correlations in nucleosomal sequence preferences to create a computationally efficient approximation of the full biophysical model. Our predictions show that a clear distinction between unicellular and multicellular life is visible in the intrinsically encoded nucleosome affinity. Furthermore, the strength of the nucleosome positioning signals correlates with the complexity of the organism. We conclude that encoding for high nucleosome occupancy, as in the human genome, is in fact a universal feature of multicellular life.

Hierarchical regulation of the genome: global changes in nucleosome organization potentiate genome response

Nucleosome occupancy is critically important in regulating access to the eukaryotic genome. Few studies in human cells have measured genome-wide nucleosome distributions at high temporal resolution during a response to a common stimulus. We measured nucleosome distributions at high temporal resolution following Kaposi's-sarcoma-associated herpesvirus (KSHV) reactivation using our newly developed mTSS-seq technology, which maps nucleosome distribution at the transcription start sites (TSS) of all human genes. Nucleosomes underwent widespread changes in organization 24 hours after KSHV reactivation and returned to their basal nucleosomal architecture 48 hours after KSHV reactivation. The widespread changes consisted of an indiscriminate remodeling event resulting in the loss of nucleosome rotational phasing signals. Additionally, one in six TSSs in the human genome possessed nucleosomes that are translationally remodeled. 72% of the loci with translationally remodeled nucleosomes have nucleosomes that moved to positions encoded by the underlying DNA sequence. Finally we demonstrated that these widespread alterations in nucleosomal architecture potentiated regulatory factor binding. These descriptions of nucleosomal architecture changes provide a new framework for understanding the role of chromatin in the genomic response, and have allowed us to propose a hierarchical model for chromatin-based regulation of genome response.

Drosophila H2A and H2A.Z Nucleosome Sequences Reveal Different Nucleosome Positioning Sequence Patterns

Nucleosomes are implicated in transcriptional regulation as well as in packing and stabilizing the DNA. Nucleosome positions affect the transcription by impeding or facilitating the binding of transcription factors. The DNA sequence, especially the periodic occurrences of dinucleotides, is a major factor that affects the nucleosome positioning. We analyzed the Drosophila DNA sequences bound by H2A and H2A.Z nucleosomes. Periodic patterns of dinucleotides (weak-weak/strong-strong or purine-purine/pyrimidine-pyrimidine) were identified as WW/SS and RR/YY nucleosome positioning sequence (NPS) patterns. The WW/SS NPS pattern of the H2A nucleosome has a 10-bp period of weak-weak/strong-strong (W = A or T; S = G or C) dinucleotides. The 10-bp periodicity, however, is disrupted in the middle of the sequence. At the dyad, the SS dinucleotide is preferred. On the other hand, the RR/YY NPS pattern has an 18-bp periodicity of purine-purine/pyrimidine-pyrimidine (R = A or G; Y = T or C) dinucleotides. The NPS patterns from H2A.Z nucleosomes differ from the NPS patterns from H2A nucleosomes. The RR/YY pattern of H2A.Z nucleosomes has major peaks shifted by 10 bp deviated from the H2A nucleosome pattern. The H2A and H2A.Z nucleosomes have different sequence preferences. The shifted peaks coincide with DNA regions interacting with the histone loops.

Genomic Nucleosome Organization Reconstituted with Pure Proteins

Chromatin remodelers regulate genes by organizing nucleosomes around promoters, but their individual contributions are obfuscated by the complex in vivo milieu of factor redundancy and indirect effects. Genome-wide reconstitution of promoter nucleosome organization with purified proteins resolves this problem and is therefore a critical goal. Here, we reconstitute four stages of nucleosome architecture using purified components: yeast genomic DNA, histones, sequence-specific Abf1/Reb1, and remodelers RSC, ISW2, INO80, and ISW1a. We identify direct, specific, and sufficient contributions that in vivo observations validate. First, RSC clears promoters by translating poly(dA:dT) into directional nucleosome removal. Second, partial redundancy is recapitulated where INO80 alone, or ISW2 at Abf1/Reb1sites, positions +1 nucleosomes. Third, INO80 and ISW2 each align downstream nucleosomal arrays. Fourth, ISW1a tightens the spacing to canonical repeat lengths. Such a minimal set of rules and proteins establishes core mechanisms by which promoter chromatin architecture arises through a blend of redundancy and specialization.

Influence of Rotational Nucleosome Positioning on Transcription Start Site Selection in Animal Promoters

The recruitment of RNA-Pol-II to the transcription start site (TSS) is an important step in gene regulation in all organisms. Core promoter elements (CPE) are conserved sequence motifs that guide Pol-II to the TSS by interacting with specific transcription factors (TFs). However, only a minority of animal promoters contains CPEs. It is still unknown how Pol-II selects the TSS in their absence. Here we present a comparative analysis of promoters' sequence composition and chromatin architecture in five eukaryotic model organisms, which shows the presence of common and unique DNA-encoded features used to organize chromatin. Analysis of Pol-II initiation patterns uncovers that, in the absence of certain CPEs, there is a strong correlation between the spread of initiation and the intensity of the 10 bp periodic signal in the nearest downstream nucleosome. Moreover, promoters' primary and secondary initiation sites show a characteristic 10 bp periodicity in the absence of CPEs. We also show that DNA natural variants in the region immediately downstream the TSS are able to affect both the nucleosome-DNA affinity and Pol-II initiation pattern. These findings support the notion that, in addition to CPEs mediated selection, sequence-induced nucleosome positioning could be a common and conserved mechanism of TSS selection in animals.

Prediction of nucleosome positioning by the incorporation of frequencies and distributions of three different nucleotide segment lengths into a general pseudo k-tuple nucleotide composition

Motivation

Nucleosome positioning plays important roles in many eukaryotic intranuclear processes, such as transcriptional regulation and chromatin structure formation. The investigations of nucleosome positioning rules provide a deeper understanding of these intracellular processes.

Results

Nucleosome positioning prediction was performed using a model consisting of three types of variables characterizing a DNA sequence—the number of five-nucleotide sequences, the number of three-nucleotide combinations in one period of a helix, and mono- and di-nucleotide distributions in DNA fragments. Using recently proposed stringent benchmark datasets with low biases for Saccharomyces cerevisiae, Homo sapiens, Caenorhabditis elegans and Drosophila melanogaster, the present model was shown to have a better prediction performance than the recently proposed predictors. This model was able to display the common and organism-dependent factors that affect nucleosome forming and inhibiting sequences as well. Therefore, the predictors developed here can accurately predict nucleosome positioning and help determine the key factors influencing this process.

Supplementary information

Supplementary data are available at Bioinformatics online.

Mesoscale Modeling Reveals Hierarchical Looping of Chromatin Fibers Near Gene Regulatory Elements

While it is well-recognized that chromatin loops play an important role in gene regulation, structural details regarding higher order chromatin loops are only emerging. Here we present a systematic study of restrained chromatin loops ranging from 25 to 427 nucleosomes (fibers of 5-80 Kb DNA in length), mimicking gene elements studied by 3C contact data. We find that hierarchical looping represents a stable configuration that can effectively bring distant regions of the GATA-4 gene together, satisfying connections reported by 3C experiments. Additionally, we find that restrained chromatin fibers larger than 100 nucleosomes (∼20Kb) form closed plectonemes, whereas fibers shorter than 100 nucleosomes form simple hairpin loops. By studying the dependence of loop structures on internal parameters, we show that loop features are sensitive to linker histone concentration, loop length, divalent ions, and DNA linker length. Specifically, increasing loop length, linker histone concentration, and divalent ion concentration are associated with increased persistence length (or decreased bending), while varying DNA linker length in a manner similar to experimentally observed "nucleosome free regions" (found near transcription start sites) disrupts intertwining and leads to loop opening and increased persistence length in linker histone depleted (-LH) fibers. Chromatin fiber structure sensitivity to these parameters, all of which vary throughout the cell cycle, tissue type, and species, suggests that caution is warranted when using uniform polymer models to fit chromatin conformation capture genome-wide data. Furthermore, the folding geometry we observe near the transcription initiation site of the GATA-4 gene suggests that hierarchical looping provides a structural mechanism for gene inhibition, and offers tunable parameters for design of gene regulation elements.

Nucleosome Presence at AML-1 Binding Sites Inversely Correlates with Ly49 Expression: Revelations from an Informatics Analysis of Nucleosomes and Immune Cell Transcription Factors

Beyond its role in genomic organization and compaction, the nucleosome is believed to participate in the regulation of gene transcription. Here, we report a computational method to evaluate the nucleosome sensitivity for a transcription factor over a given stretch of the genome. Sensitive factors are predicted to be those with binding sites preferentially contained within nucleosome boundaries and lacking 10 bp periodicity. Based on these criteria, the Acute Myeloid Leukemia-1a (AML-1a) transcription factor, a regulator of immune gene expression, was identified as potentially sensitive to nucleosomal regulation within the mouse Ly49 gene family. This result was confirmed in RMA, a cell line with natural expression of Ly49, using MNase-Seq to generate a nucleosome map of chromosome 6, where the Ly49 gene family is located. Analysis of this map revealed a specific depletion of nucleosomes at AML-1a binding sites in the expressed Ly49A when compared to the other, silent Ly49 genes. Our data suggest that nucleosome-based regulation contributes to the expression of Ly49 genes, and we propose that this method of predicting nucleosome sensitivity could aid in dissecting the regulatory role of nucleosomes in general.

The nucleosome—a large protein complex with DNA wound around it—is the fundamental unit of genomic organization in the eukaryotic cell. More than just a DNA organizer, however, nucleosomes may control gene expression by interfering with the cell’s ability to access the wound-up DNA, as shown by recent research. In this report, we demonstrate a computational method for predicting which elements of the genome are sensitive to regulation by nucleosomes. As a proof-of-concept, we identify AML-1a binding sites—important sequences in DNA regulation—as being specifically nucleosome sensitive. We then show that AML-1a sites are specifically depleted of nucleosomes when a gene is expressed, indicating the ability for nucleosomes to suppress the expression of that gene. This finding confirms that nucleosomes are likely involved in genome regulation, and provides a method for predicting which areas of the genome are probably affected most by nucleosomes. This paper also highlights the usefulness of the Ly49 gene family in testing computer-derived genomic predictions, and is of interest to anyone studying how gene expression is regulated from cell to cell.

Classification of Promoters Based on the Combination of Core Promoter Elements Exhibits Different Histone Modification Patterns

Four different histones (H2A, H2B, H3, and H4; two subunits each) constitute a histone octamer, around which DNA wraps to form histone-DNA complexes called nucleosomes. Amino acid residues in each histone are occasionally modified, resulting in several biological effects, including differential regulation of transcription. Core promoters that encompass the transcription start site have well-conserved DNA motifs, including the initiator (Inr), TATA box, and DPE, which are collectively called the core promoter elements (CPEs). In this study, we systematically studied the associations between the CPEs and histone modifications by integrating the Drosophila Core Promoter Database and time-series ChIP-seq data for histone modifications (H3K4me3, H3K27ac, and H3K27me3) during development in Drosophila melanogaster via the modENCODE project. We classified 96 core promoters into four groups based on the presence or absence of the TATA box or DPE, calculated the histone modification ratio at the core promoter region, and transcribed region for each core promoter. We found that the histone modifications in TATA-less groups were static during development and that the core promoters could be clearly divided into three types: i) core promoters with continuous active marks (H3K4me3 and H3K27ac), ii) core promoters with a continuous inactive mark (H3K27me3) and occasional active marks, and iii) core promoters with occasional histone modifications. Linear regression analysis and non-linear regression by random forest showed that the TATA-containing groups included core promoters without histone modifications, for which the measured RNA expression values were not predictable accurately from the histone modification status. DPE-containing groups had a higher relative frequency of H3K27me3 in both the core promoter region and transcribed region. In summary, our analysis showed that there was a systematic link between the existence of the CPEs and the dynamics, frequency and influence on transcriptional activity of histone modifications.

Theoretical estimates of exposure timescales of protein binding sites on DNA regulated by nucleosome kinetics

It is being increasingly realized that nucleosome organization on DNA crucially regulates DNA-protein interactions and the resulting gene expression. While the spatial character of the nucleosome positioning on DNA has been experimentally and theoretically studied extensively, the temporal character is poorly understood. Accounting for ATPase activity and DNA-sequence effects on nucleosome kinetics, we develop a theoretical method to estimate the time of continuous exposure of binding sites of non-histone proteins (e.g. transcription factors and TATA binding proteins) along any genome. Applying the method to Saccharomyces cerevisiae, we show that the exposure timescales are determined by cooperative dynamics of multiple nucleosomes, and their behavior is often different from expectations based on static nucleosome occupancy. Examining exposure times in the promoters of GAL1 and PHO5, we show that our theoretical predictions are consistent with known experiments. We apply our method genome-wide and discover huge gene-to-gene variability of mean exposure times of TATA boxes and patches adjacent to TSS (+1 nucleosome region); the resulting timescale distributions have non-exponential tails.

Genome-wide profiling of nucleosome sensitivity and chromatin accessibility in Drosophila melanogaster

Nucleosomal DNA is thought to be generally inaccessible to DNA-binding factors, such as micrococcal nuclease (MNase). Here, we digest Drosophila chromatin with high and low concentrations of MNase to reveal two distinct nucleosome types: MNase-sensitive and MNase-resistant. MNase-resistant nucleosomes assemble on sequences depleted of A/T and enriched in G/C-containing dinucleotides, whereas MNase-sensitive nucleosomes form on A/T-rich sequences found at transcription start and termination sites, enhancers and DNase I hypersensitive sites. Estimates of nucleosome formation energies indicate that MNase-sensitive nucleosomes tend to be less stable than MNase-resistant ones. Strikingly, a decrease in cell growth temperature of about 10°C makes MNase-sensitive nucleosomes less accessible, suggesting that observed variations in MNase sensitivity are related to either thermal fluctuations of chromatin fibers or the activity of enzymatic machinery. In the vicinity of active genes and DNase I hypersensitive sites nucleosomes are organized into periodic arrays, likely due to 'phasing' off potential barriers formed by DNA-bound factors or by nucleosomes anchored to their positions through external interactions. The latter idea is substantiated by our biophysical model of nucleosome positioning and energetics, which predicts that nucleosomes immediately downstream of transcription start sites are anchored and recapitulates nucleosome phasing at active genes significantly better than sequence-dependent models.

Categorical spectral analysis of periodicity in nucleosomal DNA

DNA helical twist imposes geometric constraints on the location of histone–DNA interaction sites along nucleosomal DNA. Certain 10.5-bp periodic nucleotides in phase with these geometric constraints have been suggested to facilitate nucleosome positioning. However, the extent of nucleotide periodicity in nucleosomal DNA and its significance in directing nucleosome positioning still remain unclear. We clarify these issues by applying categorical spectral analysis to high-resolution nucleosome maps in two yeast species. We find that only a small fraction of nucleosomal sequences contain significant 10.5-bp periodicity. We further develop a spectral decomposition method to show that the previously observed periodicity in aligned nucleosomal sequences mainly results from proper phasing among nucleosomal sequences, and not from a preponderant occurrence of periodicity within individual sequences. Importantly, we show that this phasing may arise from the histones’ proclivity for putting preferred nucleotides at some of the evenly spaced histone–DNA contact points with respect to the dyad axis. We demonstrate that 10.5-bp periodicity, when present, significantly facilitates rotational, but not translational, nucleosome positioning. Finally, although periodicity only moderately affects nucleosome occupancy genome wide, reduced periodicity is an evolutionarily conserved signature of nucleosome-depleted regions around transcription start/termination sites.

A species-specific nucleosomal signature defines a periodic distribution of amino acids in proteins

Nucleosomes are the basic structural units of chromatin. Most of the yeast genome is organized in a pattern of positioned nucleosomes that is stably maintained under a wide range of physiological conditions. In this work, we have searched for sequence determinants associated with positioned nucleosomes in four species of fission and budding yeasts. We show that mononucleosomal DNA follows a highly structured base composition pattern, which differs among species despite the high degree of histone conservation. These nucleosomal signatures are present in transcribed and non-transcribed regions across the genome. In the case of open reading frames, they correctly predict the relative distribution of codons on mononucleosomal DNA, and they also determine a periodicity in the average distribution of amino acids along the proteins. These results establish a direct and species-specific connection between the position of each codon around the histone octamer and protein composition.

Nucleosome positioning: resources and tools online

Nucleosome positioning is an important process required for proper genome packing and its accessibility to execute the genetic program in a cell-specific, timely manner. In the recent years hundreds of papers have been devoted to the bioinformatics, physics and biology of nucleosome positioning. The purpose of this review is to cover a practical aspect of this field, namely, to provide a guide to the multitude of nucleosome positioning resources available online. These include almost 300 experimental datasets of genome-wide nucleosome occupancy profiles determined in different cell types and more than 40 computational tools for the analysis of experimental nucleosome positioning data and prediction of intrinsic nucleosome formation probabilities from the DNA sequence. A manually curated, up to date list of these resources will be maintained at http://generegulation.info.

The chromatin remodelers RSC and ISW1 display functional and chromatin-based promoter antagonism

ISWI family chromatin remodelers typically organize nucleosome arrays, while SWI/SNF family remodelers (RSC) typically disorganize and eject nucleosomes, implying an antagonism that is largely unexplored in vivo. Here, we describe two independent genetic screens for rsc suppressors that yielded mutations in the promoter-focused ISW1a complex or mutations in the ‘basic patch’ of histone H4 (an epitope that regulates ISWI activity), strongly supporting RSC-ISW1a antagonism in vivo. RSC and ISW1a largely co-localize, and genomic nucleosome studies using rsc isw1 mutant combinations revealed opposing functions: promoters classified with a nucleosome-deficient region (NDR) gain nucleosome occupancy in rsc mutants, but this gain is attenuated in rsc isw1 double mutants. Furthermore, promoters lacking NDRs have the highest occupancy of both remodelers, consistent with regulation by nucleosome occupancy, and decreased transcription in rsc mutants. Taken together, we provide the first genetic and genomic evidence for RSC-ISW1a antagonism and reveal different mechanisms at two different promoter architectures.

DOI:http://dx.doi.org/10.7554/eLife.06073.001

The genome of an organism can contain hundreds to thousands of genes. However, these genes are not all active at the same time. Instead, mechanisms exist that control when genes are switched off or on. One such mechanism alters how tightly DNA is packaged into a structure called chromatin. To form chromatin, DNA is wrapped around histone proteins at different points along its length; these structures are known as nucleosomes. Once formed, chromatin can either adopt a tightly packed form that represses gene activity or a less compact form associated with gene activation.

The proteins that control how chromatin is packed are called ‘chromatin remodelers’. These proteins work in complexes that either disassemble chromatin—for example, by repositioning nucleosomes or removing histones—or organize chromatin by replacing nucleosomes and making it more compact.

Studies in many organisms have identified two key families of chromatin remodelers. In yeast, the ISWI family of complexes assembles chromatin and a protein complex called RSC disassembles chromatin. Parnell, Schlichter et al. used a range of genetic techniques to investigate whether these two chromatin-remodeling complexes work against each other in a species of yeast called Saccharomyces cerevisiae. The results suggest that this is indeed the case.

Both the ISWI complex and the RSC complex bind to regions of DNA called promoters, which are found at the start of a gene and help to regulate its activity. Parnell, Schlichter et al. found that the RSC complex helps to activate genes by establishing or maintaining regions of nucleosome-poor chromatin at a promoter. The chromatin is relaxed in these regions, which allows the proteins that activate genes to access the DNA. This effect is partially counteracted by the ISWI complex, which repositions nucleosomes across the promoters to fill the gaps created by the RSC complex.

In comparison, Parnell, Schlichter et al. found that promoters that do not have regions of nucleosome-poor chromatin contain a larger number of both of the remodeling complexes and have a high turnover of histone proteins. This suggests that at these sites, the RSC proteins are needed to overcome the assembly of nucleosomes by the ISWI complex in order to activate the gene. Thus, these two chromatin remodelers, ISWI and RSC, compete at promoters to determine whether they contain or lack nucleosomes, which helps determine whether the gene is silent or active, respectively. Future work will look further at how the ‘winner’ is determined: how transcription factors or signaling systems within the cell bias the recruitment or activity of RSC or ISWI at particular promoters, to determine gene activity.

DOI:http://dx.doi.org/10.7554/eLife.06073.002

Mutation bias, rather than binding preference, underlies the nucleosome-associated G+C% variation in eukaryotes

The effects of genetic content on epigenetic status have been extensively studied, but how epigenetic status affects genetic content is not well understood. As a key epigenetic factor the nucleosome structure is highly correlated with local G+C% in eukaryotic genomes. The prevailing explanation to the pattern is that nucleosome occupancy favors higher G+C% sequences more than lower G+C% sequences. However, recent observation of a biased mutation spectrum caused by nucleosome occupancy suggests that the higher G+C% of nucleosomal DNA might be the evolutionary consequence of nucleosome occupancy. To distinguish the two explanations, we examined data from an in vitro nucleosome reconstitution experiment in which histones are incubated with yeast Saccharomyces cerevisiae and Escherichia coli genomic DNA, the former has been shaped by nucleosome structure while the latter has not. There is a strong positive correlation between nucleosome density and G+C% for the yeast DNA, an observation consistent with in vivo data, and such a pattern nearly vanishes for E. coli genomic DNA, suggesting that biased mutation, rather than biased occupancy, explains the most nucleosome-associated G+C% variation in eukaryotic genomes.

Pioneer transcription factors in cell reprogramming

Biochemical and genomic studies have shown that transcription factors with the highest reprogramming activity often have the special ability to engage their target sites on nucleosomal DNA, thus behaving as “pioneer factors” to initiate events in closed chromatin. This review by Iwafuchi-Doi and Zaret focuses on the most recent studies of pioneer factors in cell programming and reprogramming, how pioneer factors have special chromatin-binding properties, and facilitators and impediments to chromatin binding.

A subset of eukaryotic transcription factors possesses the remarkable ability to reprogram one type of cell into another. The transcription factors that reprogram cell fate are invariably those that are crucial for the initial cell programming in embryonic development. To elicit cell programming or reprogramming, transcription factors must be able to engage genes that are developmentally silenced and inappropriate for expression in the original cell. Developmentally silenced genes are typically embedded in “closed” chromatin that is covered by nucleosomes and not hypersensitive to nuclease probes such as DNase I. Biochemical and genomic studies have shown that transcription factors with the highest reprogramming activity often have the special ability to engage their target sites on nucleosomal DNA, thus behaving as “pioneer factors” to initiate events in closed chromatin. Other reprogramming factors appear dependent on pioneer factors for engaging nucleosomes and closed chromatin. However, certain genomic domains in which nucleosomes are occluded by higher-order chromatin structures, such as in heterochromatin, are resistant to pioneer factor binding. Understanding the means by which pioneer factors can engage closed chromatin and how heterochromatin can prevent such binding promises to advance our ability to reprogram cell fates at will and is the topic of this review.

Autonomous actions of the human growth hormone long-range enhancer

The human growth hormone (hGH) gene is controlled by a long-range enhancer, HSI, located 14.5 kb 5′ to the hGH promoter. HSI establishes a domain of noncoding transcription that is ‘looped’ to the hGH promoter as an essential step in initiating hGH gene expression. Thus, defining how HSI generates its domain of noncoding transcription is central to understanding its long-range function. Here, we demonstrate that activation of noncoding transcription reflects an HSI-autonomous activity fully independent of interactions with linked gene promoters and occurring in spatial and temporal synchrony with initiation of GH expression in the embryonic pituitary. HSI establishes its noncoding transcription start sites (TSS) over a defined distance from its core determinants and in a manner independent of local primary sequences. The interval between HSI and it TSS co-maps with a domain of disordered and/or highly mobile nucleosomes specific to the pituitary locus. Thus, a localized chromatin reconfiguration by HSI and consequent establishment of an adjacent domain of noncoding transcription constitute initiating events in long-range enhancer function within the hGH locus.

Global nucleosome positioning regulates salicylic acid mediated transcription in Arabidopsis thaliana

BACKGROUND

The nucleosome positioning regulates the gene expression and many other DNA-related processes in eukaryotes. Genome-wide mapping of nucleosome positions and correlation of genome-wide nucleosomal remodeling with the changes in the gene expression can help us understanding gene regulation on genome level.

RESULTS

In the present study, we correlate the gene expression and the genomic nucleosomal remodeling in response to salicylic acid (SA) treatment in A. thaliana. We have mapped genome-wide nucleosomes by performing tiling microarray using 146 bp mononucleosomal template DNA. The average nucleosomal coverage is approximately 346 bp per nucleosome both under the control and the SA-treated conditions. The nucleosomal coverage is more in the coding region than in the 5' regulatory regions. We observe approximately 50% nucleosomal remodeling on SA treatment where significant nucleosomal depletion and nucleosomal enrichment around the transcription start site (TSS) occur in SA induced genes and SA repressed genes respectively in response to SA treatment. Especially in the case of the SA-induced group, the nucleosomal remodeling over the minimal promoter in response to SA is especially significant in the Non-expresser of PR1 (NPR1)-dependent genes. A detailed investigation of npr1-1 mutant confirms a distinct role of NPR1 in the nucleosome remodeling over the core promoter. We have also identified several motifs for various hormonal responses; including ABRE elements in the remodeled nucleosomal regions around the promoter region in the SA regulated genes. We have further identified that the W-box and TGACG/C motif, reported to play an important role in SA-mediated induction, are enriched in nucleosome free regions (NFRs) of the promoter region of the SA induced genes.

CONCLUSIONS

This is the first study reporting genome-wide effects of SA treatment on the chromatin architecture of A. thaliana. It also reports significant role of NPR1 in genome-wide nucleosomal remodeling in response to SA.

Nucleosome Positioning with Set of Key Positions and Nucleosome Affinity

The formation and precise positioning of nucleosome in chromatin occupies a very important role in studying life process. Today, there are many researchers who discovered that the positioning where the location of a DNA sequence fragment wraps around a histone octamer in genome is not random but regular. However, the positioning is closely relevant to the concrete sequence of core DNA. So in this paper, we analyzed the relation between the affinity and sequence structure of core DNA, and extracted the set of key positions. In these positions, the nucleotide sequences probably occupy mainly action in the binding. First, we simplified and formatted the experimental data with the affinity. Then, to find the key positions in the wrapping, we used neural network to analyze the positive and negative effects of nucleosome generation for each position in core DNA sequences. However, we reached a class of weights with every position to describe this effect. Finally, based on the positions with high weights, we analyzed the reason why the chosen positions are key positions, and used these positions to construct a model for nucleosome positioning prediction. Experimental results show the effectiveness of our method.

An integrated model of transcription factor diffusion shows the importance of intersegmental transfer and quaternary protein structure for target site finding

We present a computational model of transcription factor motion that explains both the observed rapid target finding of transcription factors, and how this motion influences protein and genome structure. Using the Smoldyn software, we modelled transcription factor motion arising from a combination of unrestricted 3D diffusion in the nucleoplasm, sliding along the DNA filament, and transferring directly between filament sections by intersegmental transfer. This presents a fine-grain picture of the way in which transcription factors find their targets two orders of magnitude faster than 3D diffusion alone allows. Eukaryotic genomes contain sections of nucleosome free regions (NFRs) around the promoters; our model shows that the presence and size of these NFRs can be explained as their acting as antennas on which transcription factors slide to reach their targets. Additionally, our model shows that intersegmental transfer may have shaped the quaternary structure of transcription factors: sequence specific DNA binding proteins are unusually enriched in dimers and tetramers, perhaps because these allow intersegmental transfer, which accelerates target site finding. Finally, our model shows that a ‘hopping’ motion can emerge from 3D diffusion on small scales. This explains the apparently long sliding lengths that have been observed for some DNA binding proteins observed in vitro. Together, these results suggest that transcription factor diffusion dynamics help drive the evolution of protein and genome structure.

Prediction of nucleosome rotational positioning in yeast and human genomes based on sequence-dependent DNA anisotropy

BACKGROUND

An organism's DNA sequence is one of the key factors guiding the positioning of nucleosomes within a cell's nucleus. Sequence-dependent bending anisotropy dictates how DNA is wrapped around a histone octamer. One of the best established sequence patterns consistent with this anisotropy is the periodic occurrence of AT-containing dinucleotides (WW) and GC-containing dinucleotides (SS) in the nucleosomal locations where DNA is bent in the minor and major grooves, respectively. Although this simple pattern has been observed in nucleosomes across eukaryotic genomes, its use for prediction of nucleosome positioning was not systematically tested.

RESULTS

We present a simple computational model, termed the W/S scheme, implementing this pattern, without using any training data. This model accurately predicts the rotational positioning of nucleosomes both in vitro and in vivo, in yeast and human genomes. About 65 - 75% of the experimentally observed nucleosome positions are predicted with the precision of one to two base pairs. The program is freely available at http://people.rit.edu/fxcsbi/WS_scheme/. We also introduce a simple and efficient way to compare the performance of different models predicting the rotational positioning of nucleosomes.

CONCLUSIONS

This paper presents the W/S scheme to achieve accurate prediction of rotational positioning of nucleosomes, solely based on the sequence-dependent anisotropic bending of nucleosomal DNA. This method successfully captures DNA features critical for the rotational positioning of nucleosomes, and can be further improved by incorporating additional terms related to the translational positioning of nucleosomes in a species-specific manner.

Apoptotic lymphocytes of H. sapiens lose nucleosomes in GC-rich promoters

We analyzed two sets of human CD4⁺ nucleosomal DNA directly sequenced by Illumina (Solexa) high throughput sequencing method. The first set has ∼40 M sequences and was produced from the normal CD4+ T lymphocytes by micrococcal nuclease. The second set has ∼44 M sequences and was obtained from peripheral blood lymphocytes by apoptotic nucleases. The different nucleosome sets showed similar dinucleotide positioning AA/TT, GG/CC, and RR/YY (R is purine, Y - pyrimidine) patterns with periods of 10–10.4 bp. Peaks of GG/CC and AA/TT patterns were shifted by 5 bp from each other. Two types of promoters in H. sapiens: AT and GC-rich were identified. AT-rich promoters in apoptotic cell had +1 nucleosome shifts 50–60 bp downstream from those in normal lymphocytes. GC-rich promoters in apoptotic cells lost 80% of nucleosomes around transcription start sites as well as in total DNA. Nucleosome positioning was predicted by combination of {AA, TT}, {GG, CC}, {WW, SS} and {RR, YY} patterns. In our study we found that the combinations of {AA, TT} and {GG, CC} provide the best results and successfully mapped 33% of nucleosomes 147 bp long with precision ±15 bp (only 31/147 or 21% is expected).

We analyzed nucleosomal DNA of human CD4+ T normal and apoptotic lymphocytes. Dinucleotide positions (pattern) of AA/TT, GG/CC, WW/SS (W is adenine or thymine, S is guanine or cytosine) and RR/YY (R is purine, Y - pyrimidine) of nucleosome sequences in both cell conditions are similar and have period 10–10.4 bp. We successfully mapped 33% of nucleosomes with precision ±15 bp by combination of {AA, TT}, {GG, CC}, {WW, SS} and {RR, YY} patterns. We identified two types of promoters in H. sapience: AT and GC-rich. AT-rich promoters keep nucleosomes around transcription start site when GC-rich promoters lost 80% of nucleosomes during apoptosis at the same region.

Drosophila Brahma complex remodels nucleosome organizations in multiple aspects

ATP-dependent chromatin remodeling complexes regulate nucleosome organizations. In Drosophila, gene Brm encodes the core Brahma complex, the ATPase subunit of SWI/SNF class of chromatin remodelers. Its role in modulating the nucleosome landscape in vivo is unclear. In this study, we knocked down Brm in Drosophila third instar larvae to explore the changes in nucleosome profiles and global gene transcription. The results show that Brm knockdown leads to nucleosome occupancy changes throughout the entire genome with a bias in occupancy decrease. In contrast, the knockdown has limited impacts on nucleosome position shift. The knockdown also alters another important physical property of nucleosome positioning, fuzziness. Nucleosome position shift, gain or loss and fuzziness changes are all enriched in promoter regions. Nucleosome arrays around the 5' ends of genes are reorganized in five patterns as a result of Brm knockdown. Intriguingly, the concomitant changes in the genes adjacent to the Brahma-dependent remodeling regions have important roles in development and morphogenesis. Further analyses reveal abundance of AT-rich motifs for transcription factors in the remodeling regions.

Computational discovery of feature patterns in nucleosomal DNA sequences

The identification of important factors that affect nucleosome formation is critical to clarify nucleosome-forming mechanisms and the role of the nucleosome in gene regulation. Various features reported in the literature led to our hypothesis that multiple features can together contribute to nucleosome formation. Therefore, we compiled 779 features and developed a pattern discovery and scoring algorithm FFN (Finding Features for Nucleosomes) to identify feature patterns that are differentially enriched in nucleosome-forming sequences and nucleosome-depletion sequences. Applying FFN to genome-wide nucleosome occupancy data in yeast and human, we identified statistically significant feature patterns that may influence nucleosome formation, many of which are common to the two species. We found that both sequence and structural features are important in nucleosome occupancy prediction. We discovered that, even for the same feature combinations, variations in feature values may lead to differences in predictive power. We demonstrated that the identified feature patterns could be used to assist nucleosomal sequence prediction.

RSC-dependent constructive and destructive interference between opposing arrays of phased nucleosomes in yeast

RSC and SWI/SNF are related ATP-dependent chromatin remodeling machines that move nucleosomes, regulating access to DNA. We addressed their roles in nucleosome phasing relative to transcription start sites in yeast. SWI/SNF has no effect on phasing at the global level. In contrast, RSC depletion results in global nucleosome repositioning: Both upstream and downstream nucleosomal arrays shift toward the nucleosome-depleted region (NDR), with no change in spacing, resulting in a narrower and partly filled NDR. The global picture of RSC-depleted chromatin represents the average of a range of chromatin structures, with most genes showing a shift of the +1 or the -1 nucleosome into the NDR. Using RSC ChIP data reported by others, we show that RSC occupancy is highest on the coding regions of heavily transcribed genes, though not at their NDRs. We propose that RSC has a role in restoring chromatin structure after transcription. Analysis of gene pairs in different orientations demonstrates that phasing patterns reflect competition between phasing signals emanating from neighboring NDRs. These signals may be in phase, resulting in constructive interference and a regular array, or out of phase, resulting in destructive interference and fuzzy positioning. We propose a modified barrier model, in which a stable complex located at the NDR acts as a bidirectional phasing barrier. In RSC-depleted cells, this barrier has a smaller footprint, resulting in narrower NDRs. Thus, RSC plays a critical role in organizing yeast chromatin.

Nucleosomes shape DNA polymorphism and divergence

An estimated 80% of genomic DNA in eukaryotes is packaged as nucleosomes, which, together with the remaining interstitial linker regions, generate higher order chromatin structures [1]. Nucleosome sequences isolated from diverse organisms exhibit ∼10 bp periodic variations in AA, TT and GC dinucleotide frequencies. These sequence elements generate intrinsically curved DNA and help establish the histone-DNA interface. We investigated an important unanswered question concerning the interplay between chromatin organization and genome evolution: do the DNA sequence preferences inherent to the highly conserved histone core exert detectable natural selection on genomic divergence and polymorphism? To address this hypothesis, we isolated nucleosomal DNA sequences from Drosophila melanogaster embryos and examined the underlying genomic variation within and between species. We found that divergence along the D. melanogaster lineage is periodic across nucleosome regions with base changes following preferred nucleotides, providing new evidence for systematic evolutionary forces in the generation and maintenance of nucleosome-associated dinucleotide periodicities. Further, Single Nucleotide Polymorphism (SNP) frequency spectra show striking periodicities across nucleosomal regions, paralleling divergence patterns. Preferred alleles occur at higher frequencies in natural populations, consistent with a central role for natural selection. These patterns are stronger for nucleosomes in introns than in intergenic regions, suggesting selection is stronger in transcribed regions where nucleosomes undergo more displacement, remodeling and functional modification. In addition, we observe a large-scale (∼180 bp) periodic enrichment of AA/TT dinucleotides associated with nucleosome occupancy, while GC dinucleotide frequency peaks in linker regions. Divergence and polymorphism data also support a role for natural selection in the generation and maintenance of these super-nucleosomal patterns. Our results demonstrate that nucleosome-associated sequence periodicities are under selective pressure, implying that structural interactions between nucleosomes and DNA sequence shape sequence evolution, particularly in introns.

In eukaryotic cells, the majority of DNA is packaged in nucleosomes comprised of ∼147 bp of DNA wound tightly around the highly conserved histone octamer. Nucleosomal DNA from diverse organisms shows an anti-correlated ∼10 bp periodicity of AT-rich and GC-rich dinucleotides. These sequence features influence DNA bending and shape, facilitating structural interactions. We asked whether natural selection mediated through the periodic sequence preferences of nucleosomes shapes the evolution of non-protein-coding regions of D. melanogaster by examining the inter- and intra-species genomic variation relative to these fundamental chromatin building blocks. The sequence changes across nucleosome-bound regions on the melanogaster lineage mirror the observed nucleosome dinucleotide periodicities. Importantly, we show that the frequencies of polymorphisms in natural populations vary across these regions, paralleling divergence, with higher frequencies of preferred alleles. These patterns are most evident for intronic regions and indicate that non-protein coding regions are evolving toward sequences that facilitate the canonical association with the histone core. This result is consistent with the hypothesis that interactions between DNA and the core have systematic impacts on function that are subject to natural selection and are not solely due to mutational bias. These ubiquitous interactions with the histone core partially account for the evolutionary constraint observed in unannotated genomic regions, and may drive broad changes in base composition.

Reconfiguration of nucleosome-depleted regions at distal regulatory elements accompanies DNA methylation of enhancers and insulators in cancer

It is well established that cancer-associated epigenetic repression occurs concomitant with CpG island hypermethylation and loss of nucleosomes at promoters, but the role of nucleosome occupancy and epigenetic reprogramming at distal regulatory elements in cancer is still poorly understood. Here, we evaluate the scope of global epigenetic alterations at enhancers and insulator elements in prostate and breast cancer cells using simultaneous genome-wide mapping of DNA methylation and nucleosome occupancy (NOMe-seq). We find that the genomic location of nucleosome-depleted regions (NDRs) is mostly cell type specific and preferentially found at enhancers in normal cells. In cancer cells, however, we observe a global reconfiguration of NDRs at distal regulatory elements coupled with a substantial reorganization of the cancer methylome. Aberrant acquisition of nucleosomes at enhancer-associated NDRs is associated with hypermethylation and epigenetic silencing marks, and conversely, loss of nucleosomes with demethylation and epigenetic activation. Remarkably, we show that nucleosomes remain strongly organized and phased at many facultative distal regulatory elements, even in the absence of a NDR as an anchor. Finally, we find that key transcription factor (TF) binding sites also show extensive peripheral nucleosome phasing, suggesting the potential for TFs to organize NDRs genome-wide and contribute to deregulation of cancer epigenomes. Together, our findings suggest that “decommissioning” of NDRs and TFs at distal regulatory elements in cancer cells is accompanied by DNA hypermethylation susceptibility of enhancers and insulator elements, which in turn may contribute to an altered genome-wide architecture and epigenetic deregulation in malignancy.

Coregulation of transcription factor binding and nucleosome occupancy through DNA features of mammalian enhancers

Transcription factors (TFs) preferentially bind sites contained in regions of computationally predicted high nucleosomal occupancy, suggesting that nucleosomes are gatekeepers of TF binding sites. However, because of their complexity mammalian genomes contain millions of randomly occurring, unbound TF consensus binding sites. We hypothesized that the information controlling nucleosome assembly may coincide with the information that enables TFs to bind cis-regulatory elements while ignoring randomly occurring sites. Hence, nucleosomes would selectively mask genomic sites that can be contacted by TFs and thus be potentially functional. The hematopoietic pioneer TF Pu.1 maintained nucleosome depletion at macrophage-specific enhancers that displayed a broad range of nucleosome occupancy in other cell types and in reconstituted chromatin. We identified a minimal set of DNA sequence and shape features that accurately predicted both Pu.1 binding and nucleosome occupancy genome-wide. These data reveal a basic organizational principle of mammalian cis-regulatory elements whereby TF recruitment and nucleosome deposition are controlled by overlapping DNA sequence features.