Collaborative competition mechanism for gene activation in vivo

Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities

Genome-scale studies have revealed extensive, cell type-specific colocalization of transcription factors, but the mechanisms underlying this phenomenon remain poorly understood. Here, we demonstrate in macrophages and B cells that collaborative interactions of the common factor PU.1 with small sets of macrophage- or B cell lineage-determining transcription factors establish cell-specific binding sites that are associated with the majority of promoter-distal H3K4me1-marked genomic regions. PU.1 binding initiates nucleosome remodeling, followed by H3K4 monomethylation at large numbers of genomic regions associated with both broadly and specifically expressed genes. These locations serve as beacons for additional factors, exemplified by liver X receptors, which drive both cell-specific gene expression and signal-dependent responses. Together with analyses of transcription factor binding and H3K4me1 patterns in other cell types, these studies suggest that simple combinations of lineage-determining transcription factors can specify the genomic sites ultimately responsible for both cell identity and cell type-specific responses to diverse signaling inputs.

A cis-regulatory signature in ascidians and flies, independent of transcription factor binding sites

BACKGROUND

Transcription initiation is controlled by cis-regulatory modules. Although these modules are usually made of clusters of short transcription factor binding sites, a small minority of such clusters in the genome have cis-regulatory activity. This paradox is currently unsolved.

RESULTS

To identify what discriminates active from inactive clusters, we focused our attention on short topologically unconstrained clusters of two ETS and two GATA binding sites, similar to the early neural enhancer of Ciona intestinalis Otx. We first computationally identified 55 such clusters, conserved between the two Ciona genomes. In vivo assay of the activity of 19 hits identified three novel early neural enhancers, all located next to genes coexpressed with Otx. Optimization of ETS and GATA binding sites was not always sufficient to confer activity to inactive clusters. Rather, a dinucleotide sequence code associated to nucleosome depletion showed a robust correlation with enhancer potential. Identification of a large collection of Ciona regulatory regions revealed that predicted nucleosome depletion constitutes a general signature of Ciona enhancers, which is conserved between orthologous loci in the two Ciona genomes and which partitions conserved noncoding sequences into a major nucleosome-bound fraction and a minor nucleosome-free fraction with higher cis-regulatory potential. We also found this signature in a large fraction of short Drosophila cis-regulatory modules.

CONCLUSION

This study indicates that a sequence-based dinucleotide signature, previously associated with nucleosome depletion and independent of transcription factor binding sites, contributes to the definition of a local cis-regulatory potential in two metazoa, Ciona intestinalis and Drosophila melanogaster.

Noise can induce bimodality in positive transcriptional feedback loops without bistability

Transcriptional positive-feedback loops are widely associated with bistability, characterized by two stable expression states that allow cells to respond to analog signals in a digital manner. Using a synthetic system in budding yeast, we show that positive feedback involving a promoter with multiple transcription factor (TF) binding sites can induce a steady-state bimodal response without cooperative binding of the TF. Deterministic models of this system do not predict bistability. Rather, the bimodal response requires a short-lived TF and stochastic fluctuations in the TF's expression. Multiple binding sites provide these fluctuations. Because many promoters possess multiple binding sites and many TFs are unstable, positive-feedback loops in gene regulatory networks may exhibit bimodal responses, but not necessarily because of deterministic bistability, as is commonly thought.

Blurring of high-resolution data shows that the effect of intrinsic nucleosome occupancy on transcription factor binding is mostly regional, not local

Genome wide maps of nucleosome occupancy in yeast have recently been produced through deep sequencing of nuclease-protected DNA. These maps have been obtained from both crosslinked and uncrosslinked chromatin in vivo, and from chromatin assembled from genomic DNA and nucleosomes in vitro. Here, we analyze these maps in combination with existing ChIP-chip data, and with new ChIP-qPCR experiments reported here. We show that the apparent nucleosome density in crosslinked chromatin, when compared to uncrosslinked chromatin, is preferentially increased at transcription factor (TF) binding sites, suggesting a strategy for mapping generic transcription factor binding sites that would not require immunoprecipitation of a particular factor. We also confirm previous conclusions that the intrinsic, sequence dependent binding of nucleosomes helps determine the localization of TF binding sites. However, we find that the association between low nucleosome occupancy and TF binding is typically greater if occupancy at a site is averaged over a 600bp window, rather than using the occupancy at the binding site itself. We have also incorporated intrinsic nucleosome binding occupancies as weights in a computational model for TF binding, and by this measure as well we find better prediction if the high resolution nucleosome occupancy data is averaged over 600bp. We suggest that the intrinsic DNA binding specificity of nucleosomes plays a role in TF binding site selection not so much through the specification of precise nucleosome positions that permit or occlude binding, but rather through the creation of low occupancy regions that can accommodate competition from TFs through rearrangement of nucleosomes.

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

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

Dynamics and function of compact nucleosome arrays

The packaging of eukaryotic DNA into chromatin sterically occludes polymerases, recombinases and repair enzymes. How chromatin structure changes to allow their actions is unknown. We constructed defined fluorescently labeled trinucleosome arrays, allowing analysis of chromatin conformational dynamics via fluorescence resonance energy transfer (FRET). The arrays undergo reversible Mg2+-dependent folding similar to that of longer arrays studied previously. We define two intermediate conformational states in the reversible folding of the nucleosome arrays and characterize the microscopic rate constants. Nucleosome arrays are highly dynamic even when compact, undergoing conformational fluctuations on timescales in the second to microsecond range. Compact states of the arrays allow binding to DNA within the central nucleosome via site exposure. Protein binding can also drive decompaction of the arrays. Thus, our results reveal multiple modes by which spontaneous chromatin fiber dynamics allow for the invasion and action of DNA-processing protein complexes.

An ensemble model of competitive multi-factor binding of the genome

Hundreds of different factors adorn the eukaryotic genome, binding to it in large number. These DNA binding factors (DBFs) include nucleosomes, transcription factors (TFs), and other proteins and protein complexes, such as the origin recognition complex (ORC). DBFs compete with one another for binding along the genome, yet many current models of genome binding do not consider different types of DBFs together simultaneously. Additionally, binding is a stochastic process that results in a continuum of binding probabilities at any position along the genome, but many current models tend to consider positions as being either binding sites or not. Here, we present a model that allows a multitude of DBFs, each at different concentrations, to compete with one another for binding sites along the genome. The result is an "occupancy profile," a probabilistic description of the DNA occupancy of each factor at each position. We implement our model efficiently as the software package COMPETE. We demonstrate genome-wide and at specific loci how modeling nucleosome binding alters TF binding, and vice versa, and illustrate how factor concentration influences binding occupancy. Binding cooperativity between nearby TFs arises implicitly via mutual competition with nucleosomes. Our method applies not only to TFs, but also recapitulates known occupancy profiles of a well-studied replication origin with and without ORC binding. Importantly, the sequence preferences our model takes as input are derived from in vitro experiments. This ensures that the calculated occupancy profiles are the result of the forces of competition represented explicitly in our model and the inherent sequence affinities of the constituent DBFs.

What controls nucleosome positions?

The DNA of eukaryotic genomes is wrapped in nucleosomes, which strongly distort and occlude the DNA from access to most DNA-binding proteins. An understanding of the mechanisms that control nucleosome positioning along the DNA is thus essential to understanding the binding and action of proteins that carry out essential genetic functions. New genome-wide data on in vivo and in vitro nucleosome positioning greatly advance our understanding of several factors that can influence nucleosome positioning, including DNA sequence preferences, DNA methylation, histone variants and post-translational modifications, higher order chromatin structure, and the actions of transcription factors, chromatin remodelers and other DNA-binding proteins. We discuss how these factors function and ways in which they might be integrated into a unified framework that accounts for both the preservation of nucleosome positioning and the dynamic nucleosome repositioning that occur across biological conditions, cell types, developmental processes and disease.

From DNA sequence to transcriptional behaviour: a quantitative approach

Complex transcriptional behaviours are encoded in the DNA sequences of gene regulatory regions. Advances in our understanding of these behaviours have been recently gained through quantitative models that describe how molecules such as transcription factors and nucleosomes interact with genomic sequences. An emerging view is that every regulatory sequence is associated with a unique binding affinity landscape for each molecule and, consequently, with a unique set of molecule-binding configurations and transcriptional outputs. We present a quantitative framework based on existing methods that unifies these ideas. This framework explains many experimental observations regarding the binding patterns of factors and nucleosomes and the dynamics of transcriptional activation. It can also be used to model more complex phenomena such as transcriptional noise and the evolution of transcriptional regulation.

Incorporating nucleosomes into thermodynamic models of transcription regulation

Transcriptional control is central to many cellular processes, and, consequently, much effort has been devoted to understanding its underlying mechanisms. The organization of nucleosomes along promoter regions is important for this process, since most transcription factors cannot bind nucleosomal sequences and thus compete with nucleosomes for DNA access. This competition is governed by the relative concentrations of nucleosomes and transcription factors and by their respective sequence binding preferences. However, despite its importance, a mechanistic understanding of the quantitative effects that the competition between nucleosomes and factors has on transcription is still missing. Here we use a thermodynamic framework based on fundamental principles of statistical mechanics to explore theoretically the effect that different nucleosome organizations along promoters have on the activation dynamics of promoters in response to varying concentrations of the regulating factors. We show that even simple landscapes of nucleosome organization reproduce experimental results regarding the effect of nucleosomes as general repressors and as generators of obligate binding cooperativity between factors. Our modeling framework also allows us to characterize the effects that various sequence elements of promoters have on the induction threshold and on the shape of the promoter activation curves. Finally, we show that using only sequence preferences for nucleosomes and transcription factors, our model can also predict expression behavior of real promoter sequences, thereby underscoring the importance of the interplay between nucleosomes and factors in determining expression kinetics.

Two strategies for gene regulation by promoter nucleosomes

Chromatin structure is central for the regulation of gene expression, but its genome-wide organization is only beginning to be understood. Here, we examine the connection between patterns of nucleosome occupancy and the capacity to modulate gene expression upon changing conditions, i.e., transcriptional plasticity. By analyzing genome-wide data of nucleosome positioning in yeast, we find that the presence of nucleosomes close to the transcription start site is associated with high transcriptional plasticity, while nucleosomes at more distant upstream positions are negatively correlated with transcriptional plasticity. Based on this, we identify two typical promoter structures associated with low or high plasticity, respectively. The first class is characterized by a relatively large nucleosome-free region close to the start site coupled with well-positioned nucleosomes further upstream, whereas the second class displays a more evenly distributed and dynamic nucleosome positioning, with high occupancy close to the start site. The two classes are further distinguished by multiple promoter features, including histone turnover, binding site locations, H2A.Z occupancy, expression noise, and expression diversity. Analysis of nucleosome positioning in human promoters reproduces the main observations. Our results suggest two distinct strategies for gene regulation by chromatin, which are selectively employed by different genes.

Spontaneous access to DNA target sites in folded chromatin fibers

DNA wrapped in nucleosomes is sterically occluded from many protein complexes that must act on it; how such complexes gain access to nucleosomal DNA is not known. In vitro studies on isolated nucleosomes show that they undergo spontaneous partial unwrapping conformational transitions, which make the wrapped nucleosomal DNA transiently accessible. Thus, site exposure might provide a general mechanism allowing access of protein complexes to nucleosomal DNA. However, existing quantitative analyses of site exposure focused on single nucleosomes, while the presence of neighbor nucleosomes and concomitant chromatin folding might significantly influence site exposure. In this work, we carried out quantitative studies on the accessibility of nucleosomal DNA in homogeneous nucleosome arrays. Two striking findings emerged. Organization into chromatin fibers changes the accessibility of nucleosomal DNA only modestly, from approximately 3-fold decreases to approximately 8-fold increases in accessibility. This means that nucleosome arrays are intrinsically dynamic and accessible even when they are visibly condensed. In contrast, chromatin folding decreases the accessibility of linker DNA by as much as approximately 50-fold. Thus, nucleosome positioning dramatically influences the accessibility of target sites located inside nucleosomes, while chromatin folding dramatically regulates access to target sites in linker DNA.

Chromatin decouples promoter threshold from dynamic range

Chromatin influences gene expression by restricting access of DNA binding proteins to their cognate sites in the genome. Large-scale characterization of nucleosome positioning in Saccharomyces cerevisiae has revealed a stereotyped promoter organization in which a nucleosome-free region (NFR) is present within several hundred base pairs upstream of the translation start site. Many transcription factors bind within NFRs and nucleate chromatin remodelling events which then expose other cis-regulatory elements. However, it is not clear how transcription-factor binding and chromatin influence quantitative attributes of gene expression. Here we show that nucleosomes function largely to decouple the threshold of induction from dynamic range. With a series of variants of one promoter, we establish that the affinity of exposed binding sites is a primary determinant of the level of physiological stimulus necessary for substantial gene activation, and sites located within nucleosomal regions serve to scale expression once chromatin is remodelled. Furthermore, we find that the S. cerevisiae phosphate response (PHO) pathway exploits these promoter designs to tailor gene expression to different environmental phosphate levels. Our results suggest that the interplay of chromatin and binding-site affinity provides a mechanism for fine-tuning responses to the same cellular state. Moreover, these findings may be a starting point for more detailed models of eukaryotic transcriptional control.

Ultrasensitive gene regulation by positive feedback loops in nucleosome modification

Eukaryotic transcription involves the synergistic interaction of many different proteins. However, the question remains how eukaryotic promoters achieve ultrasensitive or threshold responses to changes in the concentration or activity of a single transcription factor (TF). We show theoretically that by recruiting a histone-modifying enzyme, a TF binding non-cooperatively to a single site can change the balance between opposing positive feedback loops in histone modification to produce a large change in gene expression in response to a small change in concentration of the TF. This mechanism can also generate bistable promoter responses, allowing a gene to be on in some cells and off in others, despite the cells being in identical conditions. In addition, the system provides a simple means by which the activities of many TFs could be integrated at a promoter.

Topological comparison of methods for predicting transcriptional cooperativity in yeast

Background

The cooperative interaction between transcription factors has a decisive role in the control of the fate of the eukaryotic cell. Computational approaches for characterizing cooperative transcription factors in yeast, however, are based on different rationales and provide a low overlap between their results. Because the wealth of information contained in protein interaction networks and regulatory networks has proven highly effective in elucidating functional relationships between proteins, we compared different sets of cooperative transcription factor pairs (predicted by four different computational methods) within the frame of those networks.

Results

Our results show that the overlap between the sets of cooperative transcription factors predicted by the different methods is low yet significant. Cooperative transcription factors predicted by all methods are closer and more clustered in the protein interaction network than expected by chance. On the other hand, members of a cooperative transcription factor pair neither seemed to regulate each other nor shared similar regulatory inputs, although they do regulate similar groups of target genes.

Conclusion

Despite the different definitions of transcriptional cooperativity and the different computational approaches used to characterize cooperativity between transcription factors, the analysis of their roles in the framework of the protein interaction network and the regulatory network indicates a common denominator for the predictions under study. The knowledge of the shared topological properties of cooperative transcription factor pairs in both networks can be useful not only for designing better prediction methods but also for better understanding the complexities of transcriptional control in eukaryotes.

Stress-dependent dynamics of global chromatin remodeling in yeast: dual role for SWI/SNF in the heat shock stress response

Although chromatin structure is known to affect transcriptional activity, it is not clear how broadly patterns of changes in histone modifications and nucleosome occupancy affect the dynamic regulation of transcription in response to perturbations. The identity and role of chromatin remodelers that mediate some of these changes are also unclear. Here, we performed temporal genome-wide analyses of gene expression, nucleosome occupancy, and histone H4 acetylation during the response of yeast (Saccharomyces cerevisiae) to different stresses and report several findings. First, a large class of predominantly ribosomal protein genes, whose transcription was repressed during both heat shock and stationary phase, showed strikingly contrasting histone acetylation patterns. Second, the SWI/SNF complex was required for normal activation as well as repression of genes during heat shock, and loss of SWI/SNF delayed chromatin remodeling at the promoters of activated genes. Third, Snf2 was recruited to ribosomal protein genes and Hsf1 target genes, and its occupancy of this large set of genes was altered during heat shock. Our results suggest a broad and direct dual role for SWI/SNF in chromatin remodeling, during heat shock activation as well as repression, at promoters and coding regions.

General transfer matrix formalism to calculate DNA-protein-drug binding in gene regulation: application to OR operator of phage lambda

The transfer matrix methodology is proposed as a systematic tool for the statistical-mechanical description of DNA-protein-drug binding involved in gene regulation. We show that a genetic system of several cis-regulatory modules is calculable using this method, considering explicitly the site-overlapping, competitive, cooperative binding of regulatory proteins, their multilayer assembly and DNA looping. In the methodological section, the matrix models are solved for the basic types of short- and long-range interactions between DNA-bound proteins, drugs and nucleosomes. We apply the matrix method to gene regulation at the O(R) operator of phage lambda. The transfer matrix formalism allowed the description of the lambda-switch at a single-nucleotide resolution, taking into account the effects of a range of inter-protein distances. Our calculations confirm previously established roles of the contact CI-Cro-RNAP interactions. Concerning long-range interactions, we show that while the DNA loop between the O(R) and O(L) operators is important at the lysogenic CI concentrations, the interference between the adjacent promoters P(R) and P(RM) becomes more important at small CI concentrations. A large change in the expression pattern may arise in this regime due to anticooperative interactions between DNA-bound RNA polymerases. The applicability of the matrix method to more complex systems is discussed.

Stopped-flow fluorescence resonance energy transfer for analysis of nucleosome dynamics

Macromolecular assemblies and machines undergo large-scale conformational changes as essential features of their normal function. Modern stopped-flow instrumentation and biotechnology combine to provide a powerful tool for characterizing the rates and natures of these conformational changes. Standard commercially available instruments provide extraordinary sensitivity and speed, allowing analysis of millisecond or longer timescale processes, with concentrations as low as a few nanomolar and volumes of just a few hundred microliters. One can now place specific dyes anywhere desired on a nucleic acid, and often on a protein as well. This ability allows the use of fluorescence resonance energy transfer experiments for detailed conformational analyses, even as the system is evolving rapidly over time following the initiation of a reaction. This approach is ideally suited for analysis of intrinsic properties of chromatin and of the machines that control chromatin assembly, disassembly, and function.

Chromatin-dependent cooperativity between site-specific transcription factors in vivo

Accessing binding sites in DNA wrapped around histones in condensed chromatin is an obstacle that transcription factors must overcome to regulate gene expression. Here we demonstrate cooperativity between two transcription factors, the glucocorticoid receptor (GR) and nuclear factor 1 (NF1) to bind the mouse mammary tumor virus promoter organized as regular chromatin in vivo. This cooperativity is not observed when the promoter is introduced transiently into cells. Using RNA interference to deplete NF1 protein levels in the cells, we confirmed that NF1 promotes binding of GR to the promoter. Furthermore, we observed a similar synergism between GR and NF1 binding on the endogenous 11beta-hydroxysteroid dehydrogenase promoter, also regulated by GR and NF1. Our results suggest that the chromatin architecture of the promoters does not permit strong association of GR in the absence of NF1. Therefore we propose that cooperativity among DNA binding factors in binding to their cognate recognition sites in chromatin may be an important feature in the regulation of gene expression.

The features of Drosophila core promoters revealed by statistical analysis

BACKGROUND

Experimental investigation of transcription is still a very labor- and time-consuming process. Only a few transcription initiation scenarios have been studied in detail. The mechanism of interaction between basal machinery and promoter, in particular core promoter elements, is not known for the majority of identified promoters. In this study, we reveal various transcription initiation mechanisms by statistical analysis of 3393 nonredundant Drosophila promoters.

RESULTS

Using Drosophila-specific position-weight matrices, we identified promoters containing TATA box, Initiator, Downstream Promoter Element (DPE), and Motif Ten Element (MTE), as well as core elements discovered in Human (TFIIB Recognition Element (BRE) and Downstream Core Element (DCE)). Promoters utilizing known synergetic combinations of two core elements (TATA_Inr, Inr_MTE, Inr_DPE, and DPE_MTE) were identified. We also establish the existence of promoters with potentially novel synergetic combinations: TATA_DPE and TATA_MTE. Our analysis revealed several motifs with the features of promoter elements, including possible novel core promoter element(s). Comparison of Human and Drosophila showed consistent percentages of promoters with TATA, Inr, DPE, and synergetic combinations thereof, as well as most of the same functional and mutual positions of the core elements. No statistical evidence of MTE utilization in Human was found. Distinct nucleosome positioning in particular promoter classes was revealed.

CONCLUSION

We present lists of promoters that potentially utilize the aforementioned elements/combinations. The number of these promoters is two orders of magnitude larger than the number of promoters in which transcription initiation was experimentally studied. The sequences are ready to be experimentally tested or used for further statistical analysis. The developed approach may be utilized for other species.

The tale beyond the tail: histone core domain modifications and the regulation of chromatin structure

Histone post-translational modifications occur, not only in the N-terminal tail domains, but also in the core domains. While modifications in the N-terminal tail function largely through the regulation of the binding of non-histone proteins to chromatin, based on their location in the nucleosome, core domain modifications may also function through distinct mechanisms involving structural alterations to the nucleosome. This article reviews the recent developments in regards to these novel histone modifications and discusses their important role in the regulation of chromatin structure.

Gene expression from random libraries of yeast promoters

Genomewide techniques to assay gene expression and transcription factor binding are in widespread use, but are far from providing predictive rules for the function of regulatory DNA. To investigate more intensively the grammar rules for active regulatory sequence, we made libraries from random ligations of a very restricted set of sequences. Working with the yeast Saccharomyces cerevisiae, we developed a novel screen based on the sensitivity of ascospores lacking dityrosine to treatment with lytic enzymes. We tested two separate libraries built by random ligation of a single type of activator site either for a well-characterized sporulation factor, Ndt80, or for a new sporulation-specific regulatory site that we identified and several neutral spacer elements. This selective system achieved up to 1:10(4) enrichment of the artificial sequences that were active during sporulation, allowing a high-throughput analysis of large libraries of synthetic promoters. This is not practical with methods involving direct screening for expression, such as those based on fluorescent reporters. There were very few false positives, since active promoters always passed the screen when retested. The survival rate of our libraries containing roughly equal numbers of spacers and activators was a few percent that of libraries made from activators alone. The sequences of approximately 100 examples of active and inactive promoters could not be distinguished by simple binary rules; instead, the best model for the data was a linear regression fit of a quantitative measure of gene activity to multiple features of the regulatory sequence.

Explicit equilibrium modeling of transcription-factor binding and gene regulation

A computational model, GOMER, is presented that predicts transcription-factor binding and incorporates effects of cooperativity and competition.

We have developed a computational model that predicts the probability of transcription factor binding to any site in the genome. GOMER (generalizable occupancy model of expression regulation) calculates binding probabilities on the basis of position weight matrices, and incorporates the effects of cooperativity and competition by explicit calculation of coupled binding equilibria. GOMER can be used to test hypotheses regarding gene regulation that build upon this physically principled prediction of protein-DNA binding.

Transcriptional regulation by the numbers: applications

With the increasing amount of experimental data on gene expression and regulation, there is a growing need for quantitative models to describe the data and relate them to their respective context. Thermodynamic models provide a useful framework for the quantitative analysis of bacterial transcription regulation. This framework can facilitate the quantification of vastly different forms of gene expression from several well-characterized bacterial promoters that are regulated by one or two species of transcription factors; it is useful because it requires only a few parameters. As such, it provides a compact description useful for higher-level studies (e.g. of genetic networks) without the need to invoke the biochemical details of every component. Moreover, it can be used to generate hypotheses on the likely mechanisms of transcriptional control.

Chromatin remodeling complexes: strength in diversity, precision through specialization

Chromatin remodeling complexes all have an ATPase subunit important for nucleosome remodeling, but are diversified and specialized by additional associated proteins. Remodelers have well-established roles in a wide range of chromosomal processes, including transcriptional regulation and chromatin assembly. Recent work, however, has revealed remarkable new functions for remodelers, such as histone variant deposition, cohesin function, and RNA transcript elongation and termination. Remodeler complexes are tailored both compositionally and mechanistically to perform particular chromatin functions.

Low-Dose Exposure to γ Rays Induces Specific Gene Regulations in Normal Human Keratinocytes

Skin is the organ most exposed to various environmental aggressors, including ionizing radiation. Low-dose and low-dose-rate exposures to gamma rays account for most occupational, medical or environmental irradiations. To examine whether this type of exposure triggers specific molecular responses, cultured primary keratinocytes isolated from adult normal skin were irradiated with single acute doses of 1 cGy or 2 Gy. DNA microarrays containing 10,500 probes were used to assess transcriptional changes over a time course between 3 and 72 h postirradiation. Keratinocytes were studied at a differentiated stage to mimic the response of cells from the suprabasal layers of the epidermis. A major finding of this study was the identification of an important number of low-dose-specific genes (140), most of which were modulated at 48 h. Clustering analysis also revealed low-dose-specific profiles. One of these clusters (17 known genes) was further analyzed using Gibbs sampling algorithm, which led to the identification of 7 putative promoter sequences. These results show for the first time that low-dose ionizing radiation is able to induce specific transcriptional responses in human keratinocytes. Our findings support the potential usefulness of microarrays in biological dosimetry studies after low-dose exposures.

Intrinsic Histone-DNA Interactions and Low Nucleosome Density Are Important for Preferential Accessibility of Promoter Regions in Yeast

In yeast cells, preferential accessibility of the HIS3-PET56 promoter region is determined by a general property of the DNA sequence, not by defined sequence elements. In vivo, this region is largely devoid of nucleosomes, and accessibility is directly related to reduced histone density. The HIS3-PET56 and DED1 promoter regions associate poorly with histones in vitro, indicating that intrinsic nucleosome stability is a major determinant of preferential accessibility. Specific and genome-wide analyses indicate that low nucleosome density is a very common feature of yeast promoter regions that correlates poorly with transcriptional activation. Thus, the yeast genome is organized into structurally distinct promoter and nonpromoter regions whose DNA sequences inherently differ with respect to nucleosome formation. This organization ensures that transcription factors bind preferentially to appropriate sites in promoters, rather than to the excess of irrelevant sites in nonpromoter regions.

Mechanism of Transcriptional Silencing in Yeast

Transcriptional silencing is a phenomenon in which the transcription of a gene by RNA polymerase II or III is repressed or not, dependent only on the gene's chromosomal location. Two prevailing models exist for silencing: (1) steric hindrance in silenced chromatin inhibits the binding of upstream activator proteins or polymerase or (2) silencing primarily blocks steps downstream of transcription preinitiation complex formation. Here, we test these models quantitatively for the case of SIR2-dependent silencing in budding yeast, using foreign and endogenous reporter proteins, at transgenic and endogenous loci. Our results contradict both models and show instead that transcriptional silencing at several URA3 transgenes, and at the naturally silenced endogenous HMRa and HMLalpha mating type genes, acts downstream of gene activator protein binding to strongly reduce the occupancy of TFIIB, RNA polymerase II, and TFIIE at the silenced promoters.

Comparison of ABF1 and RAP1 in chromatin opening and transactivator potentiation in the budding yeast Saccharomyces cerevisiae

Autonomously replicating sequence binding factor 1 (ABF1) and repressor/activator protein 1 (RAP1) from budding yeast are multifunctional, site-specific DNA-binding proteins, with roles in gene activation and repression, replication, and telomere structure and function. Previously we have shown that RAP1 can prevent nucleosome positioning in the vicinity of its binding site and have provided evidence that this ability to create a local region of "open" chromatin contributes to RAP1 function at the HIS4 promoter by facilitating binding and activation by GCN4. Here we examine and directly compare to that of RAP1 the ability of ABF1 to create a region of open chromatin near its binding site and to contribute to activated transcription at the HIS4, ADE5,7, and HIS7 promoters. ABF1 behaves similarly to RAP1 in these assays, but it shows some subtle differences from RAP1 in the character of the open chromatin region near its binding site. Furthermore, although the two factors can similarly enhance activated transcription at the promoters tested, RAP1 binding is continuously required for this enhancement, but ABF1 binding is not. These results indicate that ABF1 and RAP1 achieve functional similarity in part via mechanistically distinct pathways.

Global nucleosome occupancy in yeast

A genome-wide study of nucleosome occupancy at yeast promoters shows that promoters that regulate active genes, contain multiple conserved motifs, or contain Rap1 binding sites tend to be depleted of nucleosomes.

Background

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

Results

Nucleosome depletion is observed in promoters that regulate active genes and/or contain multiple evolutionarily conserved motifs that recruit transcription factors. The Rap1 consensus was the only binding motif identified in a completely unbiased search of nucleosome-depleted promoters. Nucleosome depletion in the vicinity of Rap1 consensus sites in ribosomal protein gene promoters was also observed by real-time PCR and micrococcal nuclease digestion. Nucleosome occupancy in these regions was increased by the small molecule rapamycin or, in the case of the RPS11B promoter, by removing the Rap1 consensus sites.

Conclusions

The presence of transcription factor-binding motifs is an important determinant of nucleosome depletion. Most motifs are associated with marked depletion only when they appear in combination, consistent with a model in which transcription factors act collaboratively to exclude nucleosomes and gain access to target sites in the DNA. In contrast, Rap1-binding sites cause marked depletion under steady-state conditions. We speculate that nucleosome depletion enables Rap1 to define chromatin domains and alter them in response to environmental cues.

Opening the chromatin for transcription

Eukaryotic genomes are packaged into a dynamic hierarchy chromatin structure. In such a particular context, the transition from a repressed compacted chromatin to a rather extended fiber is necessary for transcription. The chromatin opening includes three events, the initial factor getting access to nucleosome DNA, local chromatin opening mediated by activator/coactivator, and transcription associated with extensive chromatin opening. Chromatin dynamics, which is DNA sequence dependent, and also occurs in condensed fiber, provides the opportunity for activators binding to DNA. Coactivators recruited by the activator open the chromatin locally. However, it appears that genes adopt distinct chromatin opening mechanisms according to whether the gene is induced expression, developmental and tissue-specific expression, or constitutive expression. In contrast to transcription initiation-related local chromatin opening, large scale of chromatin opening is associated with a functional enhancer as well as high transcription rate. How the transcription initiated from an enhancer or enhancer like modules, i.e. intergenic transcription, conducts the extensive chromatin opening is discussed. A model for long-range interaction that non-coding transcripts from enhancers may promote efficient communication with promoters is proposed.

The histone octamer is invisible when NF-kappaB binds to the nucleosome

The transcription factor NF-kappaB is involved in the transcriptional control of more than 150 genes, but the way it acts at the level of nucleosomal templates is not known. Here we report on a study examining the interaction of NF-kappaB p50 with its DNA recognition sequence in a positioned nucleosome. We demonstrate that NF-kappaB p50 was able to bind to the nucleosome with an apparent association constant close to that for free DNA. In agreement with this, the affinity of NF-kappaB p50 binding does not depend on the localization of its recognition sequence relative to the nucleosome dyad axis. In addition, the binding of NF-kappaB p50 does not induce eviction of histones and does not perturb the overall structure of the nucleosome. The NF-kappaB p50-nucleosome complex exhibits, however, local structural alterations within the NF-kappaB p50 recognition site. Importantly, these alterations were very similar to those found in the NF-kappaB p50-DNA complex. Our data suggest that NF-kappaB p50 can accommodate the distorted, bent DNA within the nucleosome. This peculiar property of NF-kappaB p50 might have evolved to meet the requirements for its function as a central switch for stress responses.

Nucleosomes facilitate their own invasion

DNA wrapped in nucleosomes is sterically occluded, creating obstacles for polymerase, regulatory, remodeling, repair and recombination complexes, which require access to the wrapped DNA. How such complexes recognize and gain access to their DNA target sites is not known. Here we report the direct detection of a dynamic equilibrium conformational transition in nucleosomes that greatly increases the distance between the end of the nucleosomal DNA and the histone core. We quantified the equilibrium constant for this transition under physiological conditions. As predicted by these findings, addition of LexA protein to nucleosomes containing the LexA target site drives this conformational equilibrium toward the unwrapped, accessible state, simultaneously allowing stable LexA binding. This inherent property of nucleosomes allows any protein, whether an energy-dependent machine or a passive binder, to gain access even to buried stretches of nucleosomal DNA.

Identifying cooperativity among transcription factors controlling the cell cycle in yeast

Transcription regulation in eukaryotes is known to occur through the coordinated action of multiple transcription factors (TFs). Recently, a few genome-wide transcription studies have begun to explore the combinatorial nature of TF interactions. We propose a novel approach that reveals how multiple TFs cooperate to regulate transcription in the yeast cell cycle. Our method integrates genome-wide gene expression data and chromatin immunoprecipitation (ChIP-chip) data to discover more biologically relevant synergistic interactions between different TFs and their target genes than previous studies. Given any pair of TFs A and B, we define a novel measure of cooperativity between the two TFs based on the expression patterns of sets of target genes of only A, only B, and both A and B. If the cooperativity measure is significant then there is reason to postulate that the presence of both TFs is needed to influence gene expression. Our results indicate that many cooperative TFs that were previously characterized experimentally indeed have high values of cooperativity measures in our analysis. In addition, we propose several novel, experimentally testable predictions of cooperative TFs that play a role in the cell cycle and other biological processes. Many of them hold interesting clues for cross talk between the cell cycle and other processes including metabolism, stress response and pseudohyphal differentiation. Finally, we have created a web tool where researchers can explore the exhaustive list of cooperative TFs and survey the graphical representation of the target genes' expression profiles. The interface includes a tool to dynamically draw a TF cooperativity network of 113 TFs with user-defined significance levels. This study is an example of how systematic combination of diverse data types along with new functional genomic approaches can provide a rigorous platform to map TF interactions more efficiently.

Formation of boundaries of transcriptionally silent chromatin by nucleosome-excluding structures

The eukaryotic genome is divided into chromosomal domains of distinct gene activities. Transcriptionally silent chromatin tends to encroach upon active chromatin. Barrier elements that can block the spread of silent chromatin have been documented, but the mechanisms of their function are not resolved. We show that the prokaryotic LexA protein can function as a barrier to the propagation of transcriptionally silent chromatin in yeast. The barrier function of LexA correlates with its ability to disrupt local chromatin structure. In accord with this, (CCGNN)(n) and poly(dA-dT), both of which do not favor nucleosome formation, can also act as efficient boundaries of silent chromatin. Moreover, we show that a Rap1p-binding barrier element also disrupts chromatin structure. These results demonstrate that nucleosome exclusion is one of the mechanisms for the establishment of boundaries of silent chromatin domains.

Structural features of transcription factor IIIA bound to a nucleosome in solution

Assembly of a DNA fragment containing a Xenopus borealis somatic-type 5S RNA gene into a nucleosome greatly restricts binding of the 5S gene-specific transcription factor IIIA (TFIIIA) to the 5S internal promoter. However, TFIIIA binds with high affinity to 5S nucleosomes lacking the N-terminal tail domains of the core histones or to nucleosomes in which these domains are hyperacetylated. The degree to which tail acetylation or removal improves TFIIIA binding cannot be simply explained by a commensurate change in the general accessibility of nucleosomal DNA. In order to investigate the molecular basis of how TFIIIA binds to the nucleosome and to ascertain if binding involves all nine zinc fingers and/or displacement of histone-DNA interactions, we examined the TFIIIA-nucleosome complex by hydroxyl radical footprinting and site-directed protein-DNA cross-linking. Our data reveal that the first six fingers of TFIIIA bind and displace approximately 20 bp of histone-DNA interactions at the periphery of the nucleosome, while binding of fingers 7 to 9 appears to overlap with histone-DNA interactions. Molecular modeling based on these results and the crystal structures of a nucleosome core and a TFIIIA-DNA cocomplex yields a precise picture of the ternary complex and a potentially important intermediate in the transition from naïve chromatin structure to productive polymerase III transcription complex.

Sin mutations alter inherent nucleosome mobility

Previous studies have identified sin mutations that alleviate the requirement for the yeast SWI/SNF chromatin remodelling complex, which include point changes in the yeast genes encoding core histones. Here we characterise the biochemical properties of nucleosomes bearing these mutations. We find that sin mutant nucleosomes have a high inherent thermal mobility. As the SWI/SNF complex can alter nucleosome positioning, the higher mobility of sin mutant nucleosomes provides a means by which sin mutations may substitute for SWI/SNF function. The location of sin mutations also provides a new opportunity for insights into the mechanism for nucleosome mobilisation. We find that both mutations altering histone DNA contacts at the nucleosome dyad and mutations in the dimer-tetramer interface influence nucleosome mobility. Furthermore, incorporation of H2A.Z into nucleosomes, which also alters dimer-tetramer interactions, affects nucleosome mobility. Thus, variation of histone sequence or subtype provides a means by which eukaryotes may regulate access to chromatin through alterations to nucleosome mobility.

HMG-D and histone H1 alter the local accessibility of nucleosomal DNA

There is evidence that HMGB proteins facilitate, while linker histones inhibit chromatin remodelling, respectively. We have examined the effects of HMG-D and histone H1/H5 on accessibility of nucleosomal DNA. Using the 601.2 nucleosome positioning sequence designed by Widom and colleagues we assembled nucleosomes in vitro and probed DNA accessibility with restriction enzymes in the presence or absence of HMG-D and histone H1/H5. For HMG-D our results show increased digestion at two spatially adjacent sites, the dyad and one terminus of nucleosomal DNA. Elsewhere varying degrees of protection from digestion were observed. The C-terminal acidic tail of HMG-D is essential for this pattern of accessibility. Neither the HMG domain by itself nor in combination with the adjacent basic region is sufficient. Histone H1/H5 binding produces two sites of increased digestion on opposite faces of the nucleosome and decreased digestion at all other sites. Our results provide the first evidence of local changes in the accessibility of nucleosomal DNA upon separate interaction with two linker binding proteins.

Combinatorial control of DNase I-hypersensitive site formation and erasure by immunoglobulin heavy chain enhancer-binding proteins

DNase I-hypersensitive sites in cellular chromatin are usually believed to be nucleosome-free regions generated by transcription factor binding. Using a cell-free system we show that hypersensitivity does not simply correlate with the number of DNA-bound proteins. Specifically, the leucine zipper containing basic helix-loop-helix protein TFE3 was sufficient to induce a DNase I-hypersensitive site at the immunoglobulin heavy chain micro enhancer in vitro. TFE3 enhanced binding of an ETS protein PU.1 to the enhancer. However, PU.1 binding erased the DNase I-hypersensitive site without abolishing TFE3 binding. Furthermore, TFE3 binding enhanced transcription in the presence and absence of a hypersensitive site, whereas endonuclease accessibility correlated strictly with DNase I hypersensitivity. We infer that chromatin constraints for transcription and nuclease sensitivity can differ.

Phylogenetically and spatially conserved word pairs associated with gene-expression changes in yeasts

BACKGROUND

Transcriptional regulation in eukaryotes often involves multiple transcription factors binding to the same transcription control region, and to understand the regulatory content of eukaryotic genomes it is necessary to consider the co-occurrence and spatial relationships of individual binding sites. The determination of conserved sequences (often known as phylogenetic footprinting) has identified individual transcription factor binding sites. We extend this concept of functional conservation to higher-order features of transcription control regions.

RESULTS

We used the genome sequences of four yeast species of the genus Saccharomyces to identify sequences potentially involved in multifactorial control of gene expression. We found 989 potential regulatory 'templates': pairs of hexameric sequences that are jointly conserved in transcription regulatory regions and also exhibit non-random relative spacing. Many of the individual sequences in these templates correspond to known transcription factor binding sites, and the sets of genes containing a particular template in their transcription control regions tend to be differentially expressed in conditions where the corresponding transcription factors are known to be active. The incorporation of word pairs to define sequence features yields more specific predictions of average expression profiles and more informative regression models for genome-wide expression data than considering sequence conservation alone.

CONCLUSIONS

The incorporation of both joint conservation and spacing constraints of sequence pairs predicts groups of target genes that are specific for common patterns of gene expression. Our work suggests that positional information, especially the relative spacing between transcription factor binding sites, may represent a common organizing principle of transcription control regions.

On schemes of combinatorial transcription logic

Cells receive a wide variety of cellular and environmental signals, which are often processed combinatorially to generate specific genetic responses. Here we explore theoretically the potentials and limitations of combinatorial signal integration at the level of cis-regulatory transcription control. Our analysis suggests that many complex transcription-control functions of the type encountered in higher eukaryotes are already implementable within the much simpler bacterial transcription system. Using a quantitative model of bacterial transcription and invoking only specific protein-DNA interaction and weak glue-like interaction between regulatory proteins, we show explicit schemes to implement regulatory logic functions of increasing complexity by appropriately selecting the strengths and arranging the relative positions of the relevant protein-binding DNA sequences in the cis-regulatory region. The architectures that emerge are naturally modular and evolvable. Our results suggest that the transcription regulatory apparatus is a "programmable" computing machine, belonging formally to the class of Boltzmann machines. Crucial to our results is the ability to regulate gene expression at a distance. In bacteria, this can be achieved for isolated genes via DNA looping controlled by the dimerization of DNA-bound proteins. However, if adopted extensively in the genome, long-distance interaction can cause unintentional intergenic cross talk, a detrimental side effect difficult to overcome by the known bacterial transcription-regulation systems. This may be a key factor limiting the genome-wide adoption of complex transcription control in bacteria. Implications of our findings for combinatorial transcription control in eukaryotes are discussed.