Structural constraints in collaborative competition of transcription factors against the nucleosome

Circadian regulation of stereotypic chromatin conformations at enhancers

Cooperation between the circadian transcription factor (TF) CLOCK:BMAL1 and other TFs at cis-regulatory elements (CREs) is critical to daily rhythms of transcription. Yet, the modalities of this cooperation are unclear. Here, we analyzed the co-binding of multiple TFs on single DNA molecules in mouse liver using single molecule footprinting (SMF). We found that SMF reads clustered in stereotypic chromatin states that reflect distinguishable organization of TFs and nucleosomes, and that were remarkably conserved between all samples. DNA protection at CLOCK:BMAL1 binding motif (E-box) varied between CREs, from E-boxes being solely bound by CLOCK:BMAL1 to situations where other TFs competed with CLOCK:BMAL1 for E-box binding. SMF also uncovered CLOCK:BMAL1 cooperative binding at E-boxes separated by 250 bp, which structurally altered the CLOCK:BMAL1-DNA interface. Importantly, we discovered multiple nucleosomes with E-boxes at entry/exit sites that were removed upon CLOCK:BMAL1 DNA binding, thereby promoting the formation of open chromatin states that facilitate DNA binding of other TFs and that were associated with rhythmic transcription. These results demonstrate the utility of SMF for studying how CLOCK:BMAL1 and other TFs regulate stereotypical chromatin states at CREs to promote transcription.

Physical interaction between MSL2 and CLAMP assures direct cooperativity and prevents competition at composite binding sites

MSL2, the DNA-binding subunit of the Drosophila dosage compensation complex, cooperates with the ubiquitous protein CLAMP to bind MSL recognition elements (MREs) on the X chromosome. We explore the nature of the cooperative binding to these GA-rich, composite sequence elements in reconstituted naïve embryonic chromatin. We found that the cooperativity requires physical interaction between both proteins. Remarkably, disruption of this interaction does not lead to indirect, nucleosome-mediated cooperativity as expected, but to competition. The protein interaction apparently not only increases the affinity for composite binding sites, but also locks both proteins in a defined dimeric state that prevents competition. High Affinity Sites of MSL2 on the X chromosome contain variable numbers of MREs. We find that the cooperation between MSL2/CLAMP is not influenced by MRE clustering or arrangement, but happens largely at the level of individual MREs. The sites where MSL2/CLAMP bind strongly in vitro locate to all chromosomes and show little overlap to an expanded set of X-chromosomal MSL2 in vivo binding sites generated by CUT&RUN. Apparently, the intrinsic MSL2/CLAMP cooperativity is limited to a small selection of potential sites in vivo. This restriction must be due to components missing in our reconstitution, such as roX2 lncRNA.

A spatio-temporally constrained gene regulatory network directed by PBX1/2 acquires limb patterning specificity via HAND2

A lingering question in developmental biology has centered on how transcription factors with widespread distribution in vertebrate embryos can perform tissue-specific functions. Here, using the murine hindlimb as a model, we investigate the elusive mechanisms whereby PBX TALE homeoproteins, viewed primarily as HOX cofactors, attain context-specific developmental roles despite ubiquitous presence in the embryo. We first demonstrate that mesenchymal-specific loss of PBX1/2 or the transcriptional regulator HAND2 generates similar limb phenotypes. By combining tissue-specific and temporally controlled mutagenesis with multi-omics approaches, we reconstruct a gene regulatory network (GRN) at organismal-level resolution that is collaboratively directed by PBX1/2 and HAND2 interactions in subsets of posterior hindlimb mesenchymal cells. Genome-wide profiling of PBX1 binding across multiple embryonic tissues further reveals that HAND2 interacts with subsets of PBX-bound regions to regulate limb-specific GRNs. Our research elucidates fundamental principles by which promiscuous transcription factors cooperate with cofactors that display domain-restricted localization to instruct tissue-specific developmental programs.

Cooperation between bHLH transcription factors and histones for DNA access

The basic helix-loop-helix (bHLH) family of transcription factors recognizes DNA motifs known as E-boxes (CANNTG) and includes 108 members1. Here we investigate how chromatinized E-boxes are engaged by two structurally diverse bHLH proteins: the proto-oncogene MYC-MAX and the circadian transcription factor CLOCK-BMAL1 (refs. 2,3). Both transcription factors bind to E-boxes preferentially near the nucleosomal entry-exit sites. Structural studies with engineered or native nucleosome sequences show that MYC-MAX or CLOCK-BMAL1 triggers the release of DNA from histones to gain access. Atop the H2A-H2B acidic patch4, the CLOCK-BMAL1 Per-Arnt-Sim (PAS) dimerization domains engage the histone octamer disc. Binding of tandem E-boxes5-7 at endogenous DNA sequences occurs through direct interactions between two CLOCK-BMAL1 protomers and histones and is important for circadian cycling. At internal E-boxes, the MYC-MAX leucine zipper can also interact with histones H2B and H3, and its binding is indirectly enhanced by OCT4 elsewhere on the nucleosome. The nucleosomal E-box position and the type of bHLH dimerization domain jointly determine the histone contact, the affinity and the degree of competition and cooperativity with other nucleosome-bound factors.

HAND2 Assists MYCN Enhancer Invasion to Regulate a Noradrenergic Neuroblastoma Phenotype

Noradrenergic neuroblastoma is characterized by a core transcriptional regulatory circuitry (CRC) comprised of transcription factors (TF) such as PHOX2B, HAND2, and GATA3, which form a network with MYCN. At normal physiologic levels, MYCN mainly binds to promoters but when aberrantly upregulated as in neuroblastoma, MYCN also binds to enhancers. Here, we investigated how MYCN invades enhancers and whether CRC TFs play a role in this process. HAND2 was found to regulate chromatin accessibility and to assist MYCN binding to enhancers. Moreover, HAND2 cooperated with MYCN to compete with nucleosomes to regulate global gene transcription. The cooperative interaction between MYCN and HAND2 could be targeted with an Aurora A kinase inhibitor plus a histone deacetylase inhibitor, resulting in potent downregulation of both MYCN and the CRC TFs and suppression of MYCN-amplified neuroblastoma tumor growth. This study identifies cooperation between MYCN and HAND2 in neuroblastoma and demonstrates that simultaneously targeting MYCN and CRC TFs is an effective way to treat this aggressive pediatric tumor.

SIGNIFICANCE

HAND2 and MYCN compete with nucleosomes to regulate global gene transcription and to drive a malignant neuroblastoma phenotype.

Profiling Accessible Chromatin and Nucleosomes in the Mammalian Genome

Genomic DNA wraps around core histones to form nucleosomes, which provides steric constraints on how transcription factors (TFs) can interact with gene regulatory sequences. It is increasingly apparent that well-positioned, accessible nucleosomes are an inherent feature of active enhancers and can facilitate cooperative TF binding, referred to as nucleosome-mediated cooperativity. Thus, profiling chromatin and nucleosome properties (accessibility, positioning, and occupancy) on the genome is crucial to understand cell-type-specific gene regulation. Here we describe a simplified protocol to profile accessible nucleosomes in the mammalian genome using low-level and high-level micrococcal nuclease (MNase) digestion followed by genome-wide sequencing.

Characterization of sequence determinants of enhancer function using natural genetic variation

Sequence variation in enhancers that control cell-type-specific gene transcription contributes significantly to phenotypic variation within human populations. However, it remains difficult to predict precisely the effect of any given sequence variant on enhancer function due to the complexity of DNA sequence motifs that determine transcription factor (TF) binding to enhancers in their native genomic context. Using F₁-hybrid cells derived from crosses between distantly related inbred strains of mice, we identified thousands of enhancers with allele-specific TF binding and/or activity. We find that genetic variants located within the central region of enhancers are most likely to alter TF binding and enhancer activity. We observe that the AP-1 family of TFs (Fos/Jun) are frequently required for binding of TEAD TFs and for enhancer function. However, many sequence variants outside of core motifs for AP-1 and TEAD also impact enhancer function, including sequences flanking core TF motifs and AP-1 half sites. Taken together, these data represent one of the most comprehensive assessments of allele-specific TF binding and enhancer function to date and reveal how sequence changes at enhancers alter their function across evolutionary timescales.

Mapping cis-regulatory elements in human neurons links psychiatric disease heritability and activity-regulated transcriptional programs

Genome-wide association studies (GWASs) have identified hundreds of loci associated with psychiatric diseases, yet there is a lack of understanding of disease pathophysiology. Common risk variants can shed light on the underlying molecular mechanisms; however, identifying causal variants remains challenging. We map cis-regulatory elements in human neurons derived from pluripotent stem cells. This system allows us to determine enhancers that activate the transcription of neuronal activity-regulated gene programs, which are thought to be critical for synaptic plasticity and are not possible to identify from postmortem tissues. Using the activity-by-contact model, we create variant-to-gene maps to interpret the function of GWAS variants. Our work nominates a subset of variants to elucidate the molecular mechanisms involving GWAS-significant loci. It also highlights that in vitro human cellular models are a powerful platform for identifying and mechanistic studies of human trait-associated genetic variants in cell states that are inaccessible from other types of human samples.

Managing the Steady State Chromatin Landscape by Nucleosome Dynamics

Gene regulation arises out of dynamic competition between nucleosomes, transcription factors, and other chromatin proteins for the opportunity to bind genomic DNA. The timescales of nucleosome assembly and binding of factors to DNA determine the outcomes of this competition at any given locus. Here, we review how these properties of chromatin proteins and the interplay between the dynamics of different factors are critical for gene regulation. We discuss how molecular structures of large chromatin-associated complexes, kinetic measurements, and high resolution mapping of protein-DNA complexes in vivo set the boundary conditions for chromatin dynamics, leading to models of how the steady state behaviors of regulatory elements arise. Expected final online publication date for the Annual Review of Biochemistry, Volume 91 is June 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Cooperative binding between distant transcription factors is a hallmark of active enhancers

Enhancers harbor binding motifs that recruit transcription factors (TFs) for gene activation. While cooperative binding of TFs at enhancers is known to be critical for transcriptional activation of a handful of developmental enhancers, the extent of TF cooperativity genome-wide is unknown. Here, we couple high-resolution nuclease footprinting with single-molecule methylation profiling to characterize TF cooperativity at active enhancers in the Drosophila genome. Enrichment of short micrococcal nuclease (MNase)-protected DNA segments indicates that the majority of enhancers harbor two or more TF-binding sites, and we uncover protected fragments that correspond to co-bound sites in thousands of enhancers. From the analysis of co-binding, we find that cooperativity dominates TF binding in vivo at the majority of active enhancers. Cooperativity is highest between sites spaced 50 bp apart, indicating that cooperativity occurs without apparent protein-protein interactions. Our findings suggest nucleosomes promoting cooperativity because co-binding may effectively clear nucleosomes and promote enhancer function.

HOX paralogs selectively convert binding of ubiquitous transcription factors into tissue-specific patterns of enhancer activation

Gene expression programs determine cell fate in embryonic development and their dysregulation results in disease. Transcription factors (TFs) control gene expression by binding to enhancers, but how TFs select and activate their target enhancers is still unclear. HOX TFs share conserved homeodomains with highly similar sequence recognition properties, yet they impart the identity of different animal body parts. To understand how HOX TFs control their specific transcriptional programs in vivo, we compared HOXA2 and HOXA3 binding profiles in the mouse embryo. HOXA2 and HOXA3 directly cooperate with TALE TFs and selectively target different subsets of a broad TALE chromatin platform. Binding of HOX and tissue-specific TFs convert low affinity TALE binding into high confidence, tissue-specific binding events, which bear the mark of active enhancers. We propose that HOX paralogs, alone and in combination with tissue-specific TFs, generate tissue-specific transcriptional outputs by modulating the activity of TALE TFs at selected enhancers.

Super-enhancers in transcriptional regulation and genome organization

Gene expression is precisely controlled in a stage and cell-type-specific manner, largely through the interaction between cis-regulatory elements and their associated trans-acting factors. Where these components aggregate in promoters and enhancers, they are able to cooperate to modulate chromatin structure and support the engagement in long-range 3D superstructures that shape the dynamics of a cell's genomic architecture. Recently, the term 'super-enhancer' has been introduced to describe a hyper-active regulatory domain comprising a complex array of sequence elements that work together to control the key gene networks involved in cell identity. Here, we survey the unique characteristics of super-enhancers compared to other enhancer types and summarize the recent advances in our understanding of their biological role in gene regulation. In particular, we discuss their capacity to attract the formation of phase-separated condensates, and capacity to generate three-dimensional genome structures that precisely activate their target genes. We also propose a multi-stage transition model to explain the evolutionary pressure driving the development of super-enhancers in complex organisms, and highlight the potential for involvement in tumorigenesis. Finally, we discuss more broadly the role of super-enhancers in human health disorders and related potential in therapeutic interventions.

Ensembles of Breathing Nucleosomes: A Computational Study

About three-fourths of the human DNA molecules are wrapped into nucleosomes, protein spools with DNA. Nucleosomes are highly dynamic, transiently exposing their DNA through spontaneous unspooling. Recent experiments allowed to observe the DNA of an ensemble of such breathing nucleosomes through x-ray diffraction with contrast matching between the solvent and the protein core. In this study, we calculate such an ensemble through a Monte Carlo simulation of a coarse-grained nucleosome model with sequence-dependent DNA mechanics. Our analysis gives detailed insights into the sequence dependence of nucleosome breathing observed in the experiment and allows us to determine the adsorption energy of the DNA bound to the protein core as a function of the ionic strength. Moreover, we predict the breathing behavior of other potentially interesting sequences and compare the findings to earlier related experiments.

DNA sequence influences hexasome orientation to regulate DNA accessibility

Nucleosomes, the fundamental organizing units of eukaryotic genomes, contain ∼146 base pairs of DNA wrapped around a histone H3–H4 tetramer and two histone H2A–H2B dimers. Converting nucleosomes into hexasomes by removal of a H2A–H2B dimer is an important regulatory event, but its regulation and functional consequences are not well-understood. To investigate the influence of hexasomes on DNA accessibility, we used the property of the Widom-601 Nucleosome Positioning Sequence (NPS) to form homogeneously oriented hexasomes in vitro. We find that DNA accessibility to transcription factors (TF) on the hexasome H2A–H2B distal side is identical to naked DNA, while the accessibility on the H2A–H2B proximal side is reduced by 2-fold, which is due to a 2-fold reduction in hexasome unwrapping probability. We then determined that a 23 bp region of the Widom-601 NPS is responsible for forming homogeneously oriented hexasomes. Analysis of published ChIP-exo data of hexasome containing genes identified two DNA sequence motifs that correlate with hexasome orientation in vivo, while ExoIII mapping studies of these sequences revealed they generate homogeneously oriented hexasomes in vitro. These results indicate that hexasome orientation, which is influenced by the underlying DNA sequence in vivo, is important for modulating DNA accessibility to regulate transcription.

Quantitative Modeling of Nucleosome Unwrapping from Both Ends

In eukaryotic cells, DNA is packaged into chromatin where nucleosomes are the basic packaging unit. Important cellular processes including gene expression, DNA replication, and DNA repair require nucleosomal DNA to be unwrapped so that functional proteins can access their target sites, which otherwise are sterically occluded. A key question in this process is what the unwrapped conformations individual nucleosomes adopt within chromatin are. Here, we develop a concurrent nucleosome unwrapping model to address this question. We hypothesize that for a given end-to-end distance of the nucleosomal DNA, the nucleosomal DNA stochastically unwraps from the histone core from both ends independently and that this combination of unwrapping from both sides results in a significant increase in the average distance between the DNA extending from both sides of the nucleosomes. We test our model on recently published experiments using a DNA origami nanocaliper that quantifies nucleosome unwrapping and achieve good agreement between experiment and model prediction. We then investigate the DNA origami caliper distribution when attached to a hexasome (a nucleosome lacking an H2A/H2B dimer). A significant shift in the caliper angle distribution caused by the asymmetric structural features of the hexasome seen experimentally is consistent with the model. Our modeling approach may be more broadly useful to the interpretation of other studies of nucleosome dynamics, chromatin dynamics, and regulatory processes involving nucleosome unwrapping, as well as more generally to optimization of future DNA origami designs to probe mechanical properties of biomolecules.

Biophysics of Chromatin Dynamics

Nucleosomes and chromatin control eukaryotic genome accessibility and thereby regulate DNA processes, including transcription, replication, and repair. Conformational dynamics within the nucleosome and chromatin structure play a key role in this regulatory function. Structural fluctuations continuously expose internal DNA sequences and nucleosome surfaces, thereby providing transient access for the nuclear machinery. Progress in structural studies of nucleosomes and chromatin has provided detailed insight into local chromatin organization and has set the stage for recent in-depth investigations of the structural dynamics of nucleosomes and chromatin fibers. Here, we discuss the dynamic processes observed in chromatin over different length scales and timescales and review current knowledge about the biophysics of distinct structural transitions.

Natural Genetic Variation Reveals Key Features of Epigenetic and Transcriptional Memory in Virus-Specific CD8 T Cells

Stable changes in chromatin states and gene expression in cells of the immune system form the basis for memory of infections and other challenges. Here, we used naturally occurring cis-regulatory variation in wild-derived inbred mouse strains to explore the mechanisms underlying long-lasting versus transient gene regulation in CD8 T cells responding to acute viral infection. Stably responsive DNA elements were characterized by dramatic and congruent chromatin remodeling events affecting multiple neighboring sites and required distinct transcription factor (TF) binding motifs for their accessibility. Specifically, we found that cooperative recruitment of T-box and Runx family transcription factors to shared targets mediated stable chromatin remodeling upon T cell activation. Our observations provide insights into the molecular mechanisms driving virus-specific CD8 T cell responses and suggest a general mechanism for the formation of transcriptional and epigenetic memory applicable to other immune and non-immune cells.

Pou5f3, SoxB1, and Nanog remodel chromatin on high nucleosome affinity regions at zygotic genome activation

The zebrafish embryo is transcriptionally mostly quiescent during the first 10 cell cycles, until the main wave of zygotic genome activation (ZGA) occurs, accompanied by fast chromatin remodeling. At ZGA, homologs of the mammalian stem cell transcription factors (TFs) Pou5f3, Nanog, and Sox19b bind to thousands of developmental enhancers to initiate transcription. So far, how these TFs influence chromatin dynamics at ZGA has remained unresolved. To address this question, we analyzed nucleosome positions in wild-type and maternal-zygotic (MZ) mutants for pou5f3 and nanog by MNase-seq. We show that Nanog, Sox19b, and Pou5f3 bind to the high nucleosome affinity regions (HNARs). HNARs are spanning over 600 bp, featuring high in vivo and predicted in vitro nucleosome occupancy and high predicted propeller twist DNA shape value. We suggest a two-step nucleosome destabilization-depletion model, in which the same intrinsic DNA properties of HNAR promote both high nucleosome occupancy and differential binding of TFs. In the first step, already before ZGA, Pou5f3 and Nanog destabilize nucleosomes at HNAR centers genome-wide. In the second step, post-ZGA, Nanog, Pou5f3, and SoxB1 maintain open chromatin state on the subset of HNARs, acting synergistically. Nanog binds to the HNAR center, whereas the Pou5f3 stabilizes the flanks. The HNAR model will provide a useful tool for genome regulatory studies in a variety of biological systems.

AP-1 Transcription Factors and the BAF Complex Mediate Signal-Dependent Enhancer Selection

Enhancer elements are genomic regulatory sequences that direct the selective expression of genes so that genetically identical cells can differentiate and acquire the highly specialized forms and functions required to build a functioning animal. To differentiate, cells must select from among the ∼10⁶ enhancers encoded in the genome the thousands of enhancers that drive the gene programs that impart their distinct features. We used a genetic approach to identify transcription factors (TFs) required for enhancer selection in fibroblasts. This revealed that the broadly expressed, growth-factor-inducible TFs FOS/JUN (AP-1) play a central role in enhancer selection. FOS/JUN selects enhancers together with cell-type-specific TFs by collaboratively binding to nucleosomal enhancers and recruiting the SWI/SNF (BAF) chromatin remodeling complex to establish accessible chromatin. These experiments demonstrate how environmental signals acting via FOS/JUN and BAF coordinate with cell-type-specific TFs to select enhancer repertoires that enable differentiation during development.

Regulation of circadian clock transcriptional output by CLOCK:BMAL1

The mammalian circadian clock relies on the transcription factor CLOCK:BMAL1 to coordinate the rhythmic expression of 15% of the transcriptome and control the daily regulation of biological functions. The recent characterization of CLOCK:BMAL1 cistrome revealed that although CLOCK:BMAL1 binds synchronously to all of its target genes, its transcriptional output is highly heterogeneous. By performing a meta-analysis of several independent genome-wide datasets, we found that the binding of other transcription factors at CLOCK:BMAL1 enhancers likely contribute to the heterogeneity of CLOCK:BMAL1 transcriptional output. While CLOCK:BMAL1 rhythmic DNA binding promotes rhythmic nucleosome removal, it is not sufficient to generate transcriptionally active enhancers as assessed by H3K27ac signal, RNA Polymerase II recruitment, and eRNA expression. Instead, the transcriptional activity of CLOCK:BMAL1 enhancers appears to rely on the activity of ubiquitously expressed transcription factors, and not tissue-specific transcription factors, recruited at nearby binding sites. The contribution of other transcription factors is exemplified by how fasting, which effects several transcription factors but not CLOCK:BMAL1, either decreases or increases the amplitude of many rhythmically expressed CLOCK:BMAL1 target genes. Together, our analysis suggests that CLOCK:BMAL1 promotes a transcriptionally permissive chromatin landscape that primes its target genes for transcription activation rather than directly activating transcription, and provides a new framework to explain how environmental or pathological conditions can reprogram the rhythmic expression of clock-controlled genes.

Multifaceted Targeting of the Chromatin Mediates Gonadotropin-Releasing Hormone Effects on Gene Expression in the Gonadotrope

Gonadotropin-releasing hormone (GnRH) stimulates the expression of multiple genes in the pituitary gonadotropes, most notably to induce synthesis of the gonadotropins, luteinizing hormone (LH), and follicle-stimulating hormone (FSH), but also to ensure the appropriate functioning of these cells at the center of the mammalian reproductive endocrine axis. Aside from the activation of gene-specific transcription factors, GnRH stimulates through its membrane-bound receptor, alterations in the chromatin that facilitate transcription of its target genes. These include changes in the histone and DNA modifications, nucleosome positioning, and chromatin packaging at the regulatory regions of each gene. The requirements for each of these events vary according to the DNA sequence which determines the basal chromatin packaging at the regulatory regions. Despite considerable progress in this field in recent years, we are only beginning to understand some of the complexities involved in the role and regulation of this chromatin structure, including new modifications, extensive cross talk, histone variants, and the actions of distal enhancers and non-coding RNAs. This short review aims to integrate the latest findings on GnRH-induced alterations in the chromatin of its target genes, which indicate multiple and diverse actions. Understanding these processes is illuminating not only in the context of the activation of these hormones during the reproductive life span but may also reveal how aberrant epigenetic regulation of these genes leads to sub-fertility.

The role of DNA sequence in nucleosome breathing

Roughly 3/4 of human genomes are sequestered by nucleosomes, DNA spools with a protein core, dictating a broad range of biological processes, ranging from gene regulation, recombination, and replication, to chromosome condensation. Nucleosomes are dynamical structures and temporarily expose wrapped DNA through spontaneous unspooling from either end, a process called site exposure or nucleosome breathing. Here we ask how this process is influenced by the mechanical properties of the wrapped DNA, which is known to depend on the underlying base pair sequence. Using a coarse-grained nucleosome model we calculate the accessibility profiles for site exposure. We find that the process is very sensitive to sequence effects, so that evolution could potentially tune the accessibility of nucleosomal DNA and would only need a small number of mutations to do so.

MADS-domain transcription factors and the floral quartet model of flower development: linking plant development and evolution

The floral quartet model of floral organ specification poses that different tetramers of MIKC-type MADS-domain transcription factors control gene expression and hence the identity of floral organs during development. Here, we provide a brief history of the floral quartet model and review several lines of recent evidence that support the model. We also describe how the model has been used in contemporary developmental and evolutionary biology to shed light on enigmatic topics such as the origin of land and flowering plants. Finally, we suggest a novel hypothesis describing how floral quartet-like complexes may interact with chromatin during target gene activation and repression.

Combinatorial function of transcription factors and cofactors

Differential gene expression gives rise to the many cell types of complex organisms. Enhancers regulate transcription by binding transcription factors (TFs), which in turn recruit cofactors to activate RNA Polymerase II at core promoters. Transcriptional regulation is typically mediated by distinct combinations of TFs, enabling a relatively small number of TFs to generate a large diversity of cell types. However, how TFs achieve combinatorial enhancer control and how enhancers, enhancer-bound TFs, and the cofactors they recruit regulate RNA Polymerase II activity is not entirely clear. Here, we review how TF synergy is mediated at the level of DNA binding and after binding, the role of cofactors and the post-translational modifications they catalyze, and discuss different models of enhancer-core-promoter communication.

CENP-A and H3 Nucleosomes Display a Similar Stability to Force-Mediated Disassembly

Centromere-specific nucleosomes are a central feature of the kinetochore complex during mitosis, in which microtubules exert pulling and pushing forces upon the centromere. CENP-A nucleosomes have been assumed to be structurally unique, thereby providing resilience under tension relative to their H3 canonical counterparts. Here, we directly test this hypothesis by subjecting CENP-A and H3 octameric nucleosomes, assembled on random or on centromeric DNA sequences, to varying amounts of applied force by using single-molecule magnetic tweezers. We monitor individual disassembly events of CENP-A and H3 nucleosomes. Regardless of the DNA sequence, the force-mediated disassembly experiments for CENP-A and H3 nucleosomes demonstrate similar rupture forces, life time residency and disassembly steps. From these experiments, we conclude that CENP-A does not, by itself, contribute unique structural features to the nucleosome that lead to a significant resistance against force-mediated disruption. The data present insights into the mechanistic basis for how CENP-A nucleosomes might contribute to the structural foundation of the centromere in vivo.

Tension-Dependent Free Energies of Nucleosome Unwrapping

Nucleosomes form the basic unit of compaction within eukaryotic genomes, and their locations represent an important, yet poorly understood, mechanism of genetic regulation. Quantifying the strength of interactions within the nucleosome is a central problem in biophysics and is critical to understanding how nucleosome positions influence gene expression. By comparing to single-molecule experiments, we demonstrate that a coarse-grained molecular model of the nucleosome can reproduce key aspects of nucleosome unwrapping. Using detailed simulations of DNA and histone proteins, we calculate the tension-dependent free energy surface corresponding to the unwrapping process. The model reproduces quantitatively the forces required to unwrap the nucleosome and reveals the role played by electrostatic interactions during this process. We then demonstrate that histone modifications and DNA sequence can have significant effects on the energies of nucleosome formation. Most notably, we show that histone tails contribute asymmetrically to the stability of the outer and inner turn of nucleosomal DNA and that depending on which histone tails are modified, the tension-dependent response is modulated differently.

Nucleosome dynamics: Sequence matters

About three quarter of all eukaryotic DNA is wrapped around protein cylinders, forming nucleosomes. Even though the histone proteins that make up the core of nucleosomes are highly conserved in evolution, nucleosomes can be very different from each other due to posttranslational modifications of the histones. Another crucial factor in making nucleosomes unique has so far been underappreciated: the sequence of their DNA. This review provides an overview of the experimental and theoretical progress that increasingly points to the importance of the nucleosomal base pair sequence. Specifically, we discuss the role of the underlying base pair sequence in nucleosome positioning, sliding, breathing, force-induced unwrapping, dissociation and partial assembly and also how the sequence can influence higher-order structures. A new view emerges: the physical properties of nucleosomes, especially their dynamical properties, are determined to a large extent by the mechanical properties of their DNA, which in turn depends on DNA sequence.

Why Do Nucleosomes Unwrap Asymmetrically?

Nucleosomes, DNA spools with a protein core, engage about three-quarters of eukaryotic DNA and play a critical role in chromosomal processes, ranging from gene regulation, recombination, and replication to chromosome condensation. For more than a decade, micromanipulation experiments where nucleosomes are put under tension, as well as the theoretical interpretations of these experiments, have deepened our understanding of the stability and dynamics of nucleosomes. Here we give a theoretical explanation for a surprising new experimental finding: nucleosomes wrapped onto the 601 positioning sequence (the sequence used in most laboratories) respond highly asymmetrically to external forces by always unwrapping from the same end. Using a computational nucleosome model, we show that this asymmetry can be explained by differences in the DNA mechanics of two very short stretches on the wrapped DNA portion. Our finding suggests that the physical properties of nucleosomes, here the response to forces, can be tuned locally by the choice of the underlying base-pair sequence. This leads to a new view of nucleosomes: a physically highly varied set of DNA-protein complexes whose properties can be tuned on evolutionary time scales to their specific function in the genomic context.

Complementary Activities of TELOMERE REPEAT BINDING Proteins and Polycomb Group Complexes in Transcriptional Regulation of Target Genes

In multicellular organisms, Polycomb Repressive Complex 1 (PRC1) and PRC2 repress target genes through histone modification and chromatin compaction. Arabidopsis thaliana mutants strongly compromised in the pathway cannot develop differentiated organs. LIKE HETEROCHROMATIN PROTEIN1 (LHP1) is so far the only known plant PRC1 component that directly binds to H3K27me3, the histone modification set by PRC2, and also associates genome-wide with trimethylation of lysine 27 of histone H3 (H3K27me3). Surprisingly, lhp1 mutants show relatively mild phenotypic alterations. To explain this paradox, we screened for genetic enhancers of lhp1 mutants to identify novel components repressing target genes together with, or in parallel to, LHP1. Two enhancing mutations were mapped to TELOMERE REPEAT BINDING PROTEIN1 (TRB1) and its paralog TRB3. We show that TRB1 binds to thousands of genomic sites containing telobox or related cis-elements with a significant increase of sites and strength of binding in the lhp1 background. Furthermore, in combination with lhp1, but not alone, trb1 mutants show increased transcription of LHP1 targets, such as floral meristem identity genes, which are more likely to be bound by TRB1 in the lhp1 background. By contrast, expression of a subset of LHP1-independent TRB1 target genes, many involved in primary metabolism, is decreased in the absence of TRB1 alone. Thus, TRB1 is a bivalent transcriptional modulator that maintains downregulation of Polycomb Group (PcG) target genes in lhp1 mutants, while it sustains high expression of targets that are regulated independently of PcG.

Dynamic regulation of transcription factors by nucleosome remodeling

The chromatin landscape and promoter architecture are dominated by the interplay of nucleosome and transcription factor (TF) binding to crucial DNA sequence elements. However, it remains unclear whether nucleosomes mobilized by chromatin remodelers can influence TFs that are already present on the DNA template. In this study, we investigated the interplay between nucleosome remodeling, by either yeast ISW1a or SWI/SNF, and a bound TF. We found that a TF serves as a major barrier to ISW1a remodeling, and acts as a boundary for nucleosome repositioning. In contrast, SWI/SNF was able to slide a nucleosome past a TF, with concurrent eviction of the TF from the DNA, and the TF did not significantly impact the nucleosome positioning. Our results provide direct evidence for a novel mechanism for both nucleosome positioning regulation by bound TFs and TF regulation via dynamic repositioning of nucleosomes.

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

Cells contain thousands of genes that are encoded by molecules of DNA. In yeast and other eukaryotic organisms, this DNA is wrapped around proteins called histones to make structures called nucleosomes. This compacts the DNA and allows it to fit inside the tiny nucleus within the cell. The positioning of the nucleosomes influences how tightly packed the DNA is, which in turn influences the activity of genes. Less active genes tend to be found within regions of DNA that are tightly packed, while more active genes are found in less tightly packed regions.

To activate a gene, proteins called transcription factors bind to a section of DNA within the gene called the promoter. Enzymes known as ‘chromatin remodelers’ can alter the locations of nucleosomes on DNA to allow the transcription factors access to the promoters of particular genes. In yeast, the SWI/SNF family of chromatin remodelers can disassemble nucleosomes to promote gene activity, while the ISW1 family organises nucleosomes into closely spaced groups to repress gene activity. However, it is not clear if, or how, chromatin remodelers can influence transcription factors that are already bound to DNA.

Here, Li et al. studied the interactions between a transcription factor and the chromatin remodelers in yeast. The experiment used a piece of DNA that contained a bound transcription factor and a single nucleosome. Li et al. used a technique called ‘single molecule DNA unzipping’, which enabled them to precisely locate the position of the nucleosome and transcription factor before and after the nucleosome was remodeled. The experiments found that a chromatin remodeler called ISW1a moved the nucleosome away from the transcription factor, while a SWI/SNF chromatin remodeler moved the nucleosome towards it.

Significantly, Li et al. also found that a transcription factor is a major barrier to ISW1a's remodeling activity, suggesting that ISW1a may use transcription factors as reference points to position nucleosomes. In contrast, SWI/SNF was able to slide a nucleosome past the transcription factor, which led to the transcription factor falling off the DNA. Therefore, SWI/SNF is able to move transcription factors out of the way to deactivate genes.

Li et al. propose a new model for how chromatin remodelers can move nucleosomes and regulate transcription factors to alter gene activity. A future challenge will be to observe these types of activities in living cells.

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

Differential Nucleosome Occupancies across Oct4-Sox2 Binding Sites in Murine Embryonic Stem Cells

The binding sequence for any transcription factor can be found millions of times within a genome, yet only a small fraction of these sequences encode functional transcription factor binding sites. One of the reasons for this dichotomy is that many other factors, such as nucleosomes, compete for binding. To study how the competition between nucleosomes and transcription factors helps determine a functional transcription factor site from a predicted transcription factor site, we compared experimentally-generated in vitro nucleosome occupancy with in vivo nucleosome occupancy and transcription factor binding in murine embryonic stem cells. Using a solution hybridization enrichment technique, we generated a high-resolution nucleosome map from targeted regions of the genome containing predicted sites and functional sites of Oct4/Sox2 regulation. We found that at Pax6 and Nes, which are bivalently poised in stem cells, functional Oct4 and Sox2 sites show high amounts of in vivo nucleosome displacement compared to in vitro. Oct4 and Sox2, which are active, show no significant displacement of in vivo nucleosomes at functional sites, similar to nonfunctional Oct4/Sox2 binding. This study highlights a complex interplay between Oct4 and Sox2 transcription factors and nucleosomes among different target genes, which may result in distinct patterns of stem cell gene regulation.

Reorganization of enhancer patterns in transition from naive to primed pluripotency

Naive and primed pluripotency is characterized by distinct signaling requirements, transcriptomes, and developmental properties, but both cellular states share key transcriptional regulators: Oct4, Sox2, and Nanog. Here, we demonstrate that transition between these two pluripotent states is associated with widespread Oct4 relocalization, mirrored by global rearrangement of enhancer chromatin landscapes. Our genomic and biochemical analyses identified candidate mediators of primed state-specific Oct4 binding, including Otx2 and Zic2/3. Even when differentiation cues are blocked, premature Otx2 overexpression is sufficient to exit the naive state, induce transcription of a substantial subset of primed pluripotency-associated genes, and redirect Oct4 to previously inaccessible enhancer sites. However, the ability of Otx2 to engage new enhancer regions is determined by its levels, cis-encoded properties of the sites, and the signaling environment. Our results illuminate regulatory mechanisms underlying pluripotency and suggest that the capacity of transcription factors such as Otx2 and Oct4 to pioneer new enhancer sites is highly context dependent.

The influence of DNA shape fluctuations on fluorescence resonance energy transfer efficiency measurements in nucleosomes

Fluorescence resonance energy transfer (FRET) measurements allow one to observe site exposure in nucleosomes, i.e. the transient unwrapping of a part of the wrapped DNA from the histone octamer. In such experiments one can typically distinguish between a closed state and an open state but in principle one might hope to detect several states, each corresponding to a certain number of open binding sites. Here we show that even in an ideal FRET setup it would be hard to detect unwrapping states with intermediate levels of FRET efficiencies. As the unwrapped DNA molecule, modelled here as a wormlike chain, has a finite stiffness, shape fluctuations smear out FRET signals completely from such intermediate states.

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.

Functional roles of nucleosome stability and dynamics

Nucleosome is a histone-DNA complex known as the fundamental repeating unit of chromatin. Up to 90% of eukaryotic DNA is wrapped around consecutive octamers made of the core histones H2A, H2B, H3 and H4. Nucleosome positioning affects numerous cellular processes that require robust and timely access to genomic DNA, which is packaged into the tight confines of the cell nucleus. In living cells, nucleosome positions are determined by intrinsic histone-DNA sequence preferences, competition between histones and other DNA-binding proteins for genomic sequence, and ATP-dependent chromatin remodelers. We discuss the major energetic contributions to nucleosome formation and remodeling, focusing especially on partial DNA unwrapping off the histone octamer surface. DNA unwrapping enables efficient access to nucleosome-buried binding sites and mediates rapid nucleosome removal through concerted action of two or more DNA-binding factors. High-resolution, genome-scale maps of distances between neighboring nucleosomes have shown that DNA unwrapping and nucleosome crowding (mutual invasion of nucleosome territories) are much more common than previously thought. Ultimately, constraints imposed by nucleosome energetics on the rates of ATP-dependent and spontaneous chromatin remodeling determine nucleosome occupancy genome-wide, and shape pathways of cellular response to environmental stresses.

Arabidopsis SEPALLATA proteins differ in cooperative DNA-binding during the formation of floral quartet-like complexes

The SEPALLATA (SEP) genes of Arabidopsis thaliana encode MADS-domain transcription factors that specify the identity of all floral organs. The four Arabidopsis SEP genes function in a largely yet not completely redundant manner. Here, we analysed interactions of the SEP proteins with DNA. All of the proteins were capable of forming tetrameric quartet-like complexes on DNA fragments carrying two sequence elements termed CArG-boxes. Distances between the CArG-boxes for strong cooperative DNA-binding were in the range of 4-6 helical turns. However, SEP1 also bound strongly to CArG-box pairs separated by smaller or larger distances, whereas SEP2 preferred large and SEP4 preferred small inter-site distances for binding. Cooperative binding of SEP3 was comparatively weak for most of the inter-site distances tested. All SEP proteins constituted floral quartet-like complexes together with the floral homeotic proteins APETALA3 (AP3) and PISTILLATA (PI) on the target genes AP3 and SEP3. Our results suggest an important part of an explanation for why the different SEP proteins have largely, but not completely redundant functions in determining floral organ identity: they may bind to largely overlapping, but not identical sets of target genes that differ in the arrangement and spacing of the CArG-boxes in their cis-regulatory regions.

The sterol-based transcriptional control of human 7-dehydrocholesterol reductase (DHCR7): Evidence of a cooperative regulatory program in cholesterol synthesis

The enzyme 7-dehydrocholesterol reductase (DHCR7) catalyzes the final step of cholesterol synthesis via the Kandutsch-Russell pathway, and is crucial in maintaining cellular cholesterol levels. Its absence leads to the devastating fetal developmental disorder Smith-Lemli-Opitz Syndrome (SLOS). How this enzyme is regulated has implications in controlling not only cholesterol synthesis, but also the synthesis of Vitamin D - another product of 7-dehydrocholesterol. In this study, we look specifically at how DHCR7 is regulated by the sterol regulatory element-binding protein-2 (SREBP-2) transcription factor. Sterol regulation has previously been studied in the rat DHCR7 promoter, but we have found that its regulatory elements are not all conserved in humans. Rather, the human promoter contains two binding sites for SREBP-2 (at -155 and -55) and a binding site for the nuclear factor-Y (NF-Y) cofactor (at -136). The -155 site is a particularly responsive sterol regulatory element (SRE) which is well conserved in mammals, and was possibly overlooked in the rat promoter study. The exact location of the weaker -55 site (close to the known rat SRE) may have shifted during evolution. Furthermore, we established that the two SREs that bind SREBP-2 work in cooperation to synergistically activate DHCR7. We have previously characterized the SREs in DHCR24, the final enzyme in the alternate Bloch pathway of cholesterol synthesis. Here, comparison of the sterol regulation of these terminal enzymes demonstrates the unique cooperative system that helps to control cholesterol synthesis.

Ubiquitous nucleosome crowding in the yeast genome

Nucleosomes may undergo a conformational change in which a stretch of DNA peels off the histone octamer surface as a result of thermal fluctuations or interactions with chromatin remodelers. Thus, neighboring nucleosomes may invade each other's territories by DNA unwrapping and translocation, or through initial assembly in partially wrapped states. A recent high-resolution map of distances between dyads of neighboring nucleosomes in Saccharomyces cerevisiae reveals that nucleosomes frequently overlap DNA territories of their neighbors. This conclusion is supported by lower-resolution maps of S. cerevisiae nucleosome lengths based on micrococcal nuclease digestion and paired-end sequencing. The average length of wrapped DNA follows a stereotypical pattern in genes and promoters, correlated with the well-known distribution of nucleosome occupancy: nucleosomal DNA tends to be shorter in promoters and longer in coding regions. To explain these observations, we have developed a biophysical model that uses a 10-11-bp periodic histone-DNA binding energy profile. The profile is based on the pattern of histone-DNA contacts in nucleosome crystal structures, as well as the idea of linker length discretization caused by higher-order chromatin structure. Our model is in agreement with the observed genome-wide distributions of interdyad distances, wrapped DNA lengths, and nucleosome occupancies. Furthermore, our approach explains in vitro measurements of the accessibility of nucleosome-covered target sites and nucleosome-induced cooperativity between DNA-binding factors. We rule out several alternative scenarios of histone-DNA interactions as inconsistent with the genomic data.

A dynamic interplay of nucleosome and Msn2 binding regulates kinetics of gene activation and repression following stress

The transcription factor Msn2 mediates a significant proportion of the environmental stress response, in which a common cohort of genes changes expression in a stereotypic fashion upon exposure to any of a wide variety of stresses. We have applied genome-wide chromatin immunoprecipitation and nucleosome profiling to determine where Msn2 binds under stressful conditions and how that binding affects, and is affected by, nucleosome positioning. We concurrently determined the effect of Msn2 activity on gene expression following stress and demonstrated that Msn2 stimulates both activation and repression. We found that some genes responded to both intermittent and continuous Msn2 nuclear occupancy while others responded only to continuous occupancy. Finally, these studies document a dynamic interplay between nucleosomes and Msn2 such that nucleosomes can restrict access of Msn2 to its canonical binding sites while Msn2 can promote reposition, expulsion and recruitment of nucleosomes to alter gene expression. This interplay may allow the cell to discriminate between different types of stress signaling.

Development of a novel prediction method of cis-elements to hypothesize collaborative functions of cis-element pairs in iron-deficient rice

BACKGROUND

Cis-acting elements are essential genomic sequences that control gene expression. In higher eukaryotes, a series of cis-elements function cooperatively. However, further studies are required to examine the co-regulation of multiple cis-elements on a promoter. The aim of this study was to propose a model of cis-element networks that cooperatively regulate gene expression in rice under iron (Fe) deficiency.

RESULTS

We developed a novel clustering-free method, microarray-associated motif analyzer (MAMA), to predict novel cis-acting elements based on weighted sequence similarities and gene expression profiles in microarray analyses. Simulation of gene expression was performed using a support vector machine and based on the presence of predicted motifs and motif pairs. The accuracy of simulated gene expression was used to evaluate the quality of prediction and to optimize the parameters used in this method. Based on sequences of Oryza sativa genes upregulated by Fe deficiency, MAMA returned experimentally identified cis-elements responsible for Fe deficiency in O. sativa. When this method was applied to O. sativa subjected to zinc deficiency and Arabidopsis thaliana subjected to salt stress, several novel candidate cis-acting elements that overlap with known cis-acting elements, such as ZDRE, ABRE, and DRE, were identified. After optimization, MAMA accurately simulated more than 87% of gene expression. Predicted motifs strongly co-localized in the upstream regions of regulated genes and sequences around transcription start sites. Furthermore, in many cases, the separation (in bp) between co-localized motifs was conserved, suggesting that predicted motifs and the separation between them were important in the co-regulation of gene expression.

CONCLUSIONS

Our results are suggestive of a typical sequence model for Fe deficiency-responsive promoters and some strong candidate cis-elements that function cooperatively with known cis-elements.

Computational identification of diverse mechanisms underlying transcription factor-DNA occupancy

ChIP-based genome-wide assays of transcription factor (TF) occupancy have emerged as a powerful, high-throughput method to understand transcriptional regulation, especially on a global scale. This has led to great interest in the underlying biochemical mechanisms that direct TF-DNA binding, with the ultimate goal of computationally predicting a TF's occupancy profile in any cellular condition. In this study, we examined the influence of various potential determinants of TF-DNA binding on a much larger scale than previously undertaken. We used a thermodynamics-based model of TF-DNA binding, called “STAP,” to analyze 45 TF-ChIP data sets from Drosophila embryonic development. We built a cross-validation framework that compares a baseline model, based on the ChIP'ed (“primary”) TF's motif, to more complex models where binding by secondary TFs is hypothesized to influence the primary TF's occupancy. Candidates interacting TFs were chosen based on RNA-SEQ expression data from the time point of the ChIP experiment. We found widespread evidence of both cooperative and antagonistic effects by secondary TFs, and explicitly quantified these effects. We were able to identify multiple classes of interactions, including (1) long-range interactions between primary and secondary motifs (separated by ≤150 bp), suggestive of indirect effects such as chromatin remodeling, (2) short-range interactions with specific inter-site spacing biases, suggestive of direct physical interactions, and (3) overlapping binding sites suggesting competitive binding. Furthermore, by factoring out the previously reported strong correlation between TF occupancy and DNA accessibility, we were able to categorize the effects into those that are likely to be mediated by the secondary TF's effect on local accessibility and those that utilize accessibility-independent mechanisms. Finally, we conducted in vitro pull-down assays to test model-based predictions of short-range cooperative interactions, and found that seven of the eight TF pairs tested physically interact and that some of these interactions mediate cooperative binding to DNA.

Chromatin Immunoprecipitation (ChIP)-based genome-wide assays of transcription factor (TF) occupancy have emerged as a powerful, high throughput method to understand transcriptional regulation, especially on a global scale. Here, we utilize 45 ChIP-chip and ChIP-SEQ data sets from Drosophila to explore the underlying mechanisms of TF-DNA binding. For this, we employ a biophysically motivated computational model, in conjunction with over 300 TF motifs (binding specificities) as well as gene expression and DNA accessibility data from different developmental stages in Drosophila embryos. Our findings provide robust statistical evidence of the role played by TF-TF interactions in shaping genome-wide TF-DNA binding profiles, and thus in directing gene regulation. Our method allows us to go beyond simply recognizing the existence of such interactions, to quantifying their effects on TF occupancy. We are able to categorize the probable mechanisms of these effects as involving direct physical interactions versus accessibility-mediated indirect interactions, long-range versus short-range interactions, and cooperative versus antagonistic interactions. Our analysis reveals widespread evidence of combinatorial regulation present in recently generated ChIP data sets, and sets the stage for rich integrative models of the future that will predict cell type-specific TF occupancy values from sequence and expression data.

Widespread evidence of cooperative DNA binding by transcription factors in Drosophila development

Regulation of eukaryotic gene transcription is often combinatorial in nature, with multiple transcription factors (TFs) regulating common target genes, often through direct or indirect mutual interactions. Many individual examples of cooperative binding by directly interacting TFs have been identified, but it remains unclear how pervasive this mechanism is during animal development. Cooperative TF binding should be manifest in genomic sequences as biased arrangements of TF-binding sites. Here, we explore the extent and diversity of such arrangements related to gene regulation during Drosophila embryogenesis. We used the DNA-binding specificities of 322 TFs along with chromatin accessibility information to identify enriched spacing and orientation patterns of TF-binding site pairs. We developed a new statistical approach for this task, specifically designed to accurately assess inter-site spacing biases while accounting for the phenomenon of homotypic site clustering commonly observed in developmental regulatory regions. We observed a large number of short-range distance preferences between TF-binding site pairs, including examples where the preference depends on the relative orientation of the binding sites. To test whether these binding site patterns reflect physical interactions between the corresponding TFs, we analyzed 27 TF pairs whose binding sites exhibited short distance preferences. In vitro protein–protein binding experiments revealed that >65% of these TF pairs can directly interact with each other. For five pairs, we further demonstrate that they bind cooperatively to DNA if both sites are present with the preferred spacing. This study demonstrates how DNA-binding motifs can be used to produce a comprehensive map of sequence signatures for different mechanisms of combinatorial TF action.

Taking into account nucleosomes for predicting gene expression

The eukaryotic genome is organized in a chain of nucleosomes that consist of 145-147 bp of DNA wrapped around a histone octamer protein core. Binding of transcription factors (TF) to nucleosomal DNA is frequently impeded, which makes it a challenging task to calculate TF occupancy at a given regulatory genomic site for predicting gene expression. Here, we review methods to calculate TF binding to DNA in the presence of nucleosomes. The main theoretical problems are (i) the computation speed that is becoming a bottleneck when partial unwrapping of DNA from the nucleosome is considered, (ii) the perturbation of the binding equilibrium by the activity of ATP-dependent chromatin remodelers, which translocate nucleosomes along the DNA, and (iii) the model parameterization from high-throughput sequencing data and fluorescence microscopy experiments in living cells. We discuss strategies that address these issues to efficiently compute transcription factor binding in chromatin.

Nucleosome mapping across the CFTR locus identifies novel regulatory factors

Nucleosome positioning on the chromatin strand plays a critical role in regulating accessibility of DNA to transcription factors and chromatin modifying enzymes. Hence, detailed information on nucleosome depletion or movement at cis-acting regulatory elements has the potential to identify predicted binding sites for trans-acting factors. Using a novel method based on enrichment of mononucleosomal DNA by bacterial artificial chromosome hybridization, we mapped nucleosome positions by deep sequencing across 250 kb, encompassing the cystic fibrosis transmembrane conductance regulator (CFTR) gene. CFTR shows tight tissue-specific regulation of expression, which is largely determined by cis-regulatory elements that lie outside the gene promoter. Although multiple elements are known, the repertoire of transcription factors that interact with these sites to activate or repress CFTR expression remains incomplete. Here, we show that specific nucleosome depletion corresponds to well-characterized binding sites for known trans-acting factors, including hepatocyte nuclear factor 1, Forkhead box A1 and CCCTC-binding factor. Moreover, the cell-type selective nucleosome positioning is effective in predicting binding sites for novel interacting factors, such as BAF155. Finally, we identify transcription factor binding sites that are overrepresented in regions where nucleosomes are depleted in a cell-specific manner. This approach recognizes the glucocorticoid receptor as a novel trans-acting factor that regulates CFTR expression in vivo.

Codons support the maintenance of intrinsic DNA polymer flexibility over evolutionary timescales

Despite our long familiarity with how the genetic code specifies the amino acid sequence, we still know little about why it is organized in the way that it is. Contrary to the view that the organization of the genetic code is a "frozen accident" of evolution, recent studies have demonstrated that it is highly nonrandom, with implications for both codon assignment and usage. We hypothesize that this inherent nonrandomness may facilitate the coexistence of both sequence and structural information in DNA. Here, we take advantage of a simple metric of intrinsic DNA flexibility to analyze mutational effects on the four phosphate linkages present in any given codon. Application of a simple evolutionary neutral model of substitution to random sequences, translated with alternative genetic codes, reveals that the standard code is highly optimized to favor synonymous substitutions that maximize DNA polymer flexibility, potentially counteracting neutral evolutionary drift toward stiffer DNA caused by spontaneous deamination. Comparison to existing mutational patterns in yeast also demonstrates evidence of strong selective constraint on DNA flexibility, especially at so-called "silent" sites. We also report a fundamental relationship between DNA flexibility, codon usage bias, and several important evolutionary descriptors of comparative genomics (e.g., base composition, transition/transversion ratio, and nonsynonymous vs. synonymous substitution rate). Recent advances in structural genomics have emphasized the role of the DNA polymer's flexibility in both gene function and whole genome folding, thereby implicating possible reasons for codons to facilitate the multiplexing of both genetic and structural information within the same molecular context.

Disentangling the many layers of eukaryotic transcriptional regulation

Regulation of gene expression in eukaryotes is an extremely complex process. In this review, we break down several critical steps, emphasizing new data and techniques that have expanded current gene regulatory models. We begin at the level of DNA sequence where cis-regulatory modules (CRMs) provide important regulatory information in the form of transcription factor (TF) binding sites. In this respect, CRMs function as instructional platforms for the assembly of gene regulatory complexes. We discuss multiple mechanisms controlling complex assembly, including cooperative DNA binding, combinatorial codes, and CRM architecture. The second section of this review places CRM assembly in the context of nucleosomes and condensed chromatin. We discuss how DNA accessibility and histone modifications contribute to TF function. Lastly, new advances in chromosomal mapping techniques have provided increased understanding of intra- and interchromosomal interactions. We discuss how these topological maps influence gene regulatory models.

Sterols regulate 3β-hydroxysterol Δ24-reductase (DHCR24) via dual sterol regulatory elements: cooperative induction of key enzymes in lipid synthesis by Sterol Regulatory Element Binding Proteins

3β-Hydroxysterol Δ24-reductase (DHCR24) catalyzes a final step in cholesterol synthesis, and has been ascribed diverse functions, such as being anti-apoptotic and anti-inflammatory. How this enzyme is regulated transcriptionally by sterols is currently unclear. Some studies have suggested that its expression is regulated by Sterol Regulatory Element Binding Proteins (SREBPs) while another suggests it is through the Liver X Receptor (LXR). However, these transcription factors have opposing effects on cellular sterol levels, so it is likely that one predominates. Here we establish that sterol regulation of DHCR24 occurs predominantly through SREBP-2, and identify the particular region of the DHCR24 promoter to which SREBP-2 binds. We demonstrate that sterol regulation is mediated by two sterol regulatory elements (SREs) in the promoter of the gene, assisted by two nearby NF-Y binding sites. Moreover, we present evidence that the dual SREs work cooperatively to regulate DHCR24 expression by comparison to two known SREBP target genes, the LDL receptor with one SRE, and farnesyl-diphosphate farnesyltransferase 1, with two SREs.

Histone chaperone FACT coordinates nucleosome interaction through multiple synergistic binding events

In eukaryotic cells, DNA maintenance requires ordered disassembly and re-assembly of chromatin templates. These processes are highly regulated and require extrinsic factors such as chromatin remodelers and histone chaperones. The histone chaperone FACT (facilitates chromatin transcription) is a large heterodimeric complex with roles in transcription, replication, and repair. FACT promotes and subsequently restricts access to DNA as a result of dynamic nucleosome reorganization. However, until now, there lacked a truly quantitative assessment of the critical contacts mediating FACT function. Here, we demonstrate that FACT binds histones, DNA, and intact nucleosomes at nanomolar concentrations. We also determine roles for the histone tails in free histone and nucleosome binding by FACT. Furthermore, we propose that the conserved acidic C-terminal domain of the FACT subunit Spt16 actively displaces nucleosomal DNA to provide access to the histone octamer. Experiments with tri-nucleosome arrays indicate a possible mode for FACT binding within chromatin. Together, the data reveal that specific FACT subunits synchronize interactions with various target sites on individual nucleosomes to generate a high affinity binding event and promote reorganization.