Collaborative competition mechanism for gene activation in vivo

Deciphering Allosteric Modulation of Cancer-Associated Histone Missense Mutations

Histone mutations have been implicated in various cancers, but their mechanistic effects on chromatin dynamics remain largely unexplored. In this study, we investigated the allosteric modulation effects induced by 40 cancer-associated histone missense mutations. By combining computational approaches with experimental evidence, we assessed the allosteric and functional impact of these mutations. Our results reveal that the allosteric effects of histone mutations are position-specific, with mutations near the H3 and H4 histone N-terminal tails exhibiting the strongest long-range perturbations. Notably, we predicted seven mutations with significant allosteric effects, potentially altering nucleosome interactions. Experimental verification of H2BS64Y and H2BS64F mutations demonstrated that they disrupted normal histone function, altered the H2BK120 ubiquitination levels and genome stability, findings suggestive of their potential cancerogenic effects. Collectively, these results show that allostery may serve as a critical mechanism underlying the oncogenic potential of some histone mutations, highlighting the need for further exploration of allosteric pathways in cancer epigenetics.

Interpreting regulatory mechanisms of Hippo signaling through a deep learning sequence model

Signaling pathway components are well studied, but how they mediate cell-type-specific transcription responses is an unresolved problem. Using the Hippo pathway in mouse trophoblast stem cells as a model, we show that the DNA binding of signaling effectors is driven by cell-type-specific sequence rules that can be learned genome wide by deep learning models. Through model interpretation and experimental validation, we show that motifs for the cell-type-specific transcription factor TFAP2C enhance TEAD4/YAP1 binding in a nucleosome-range and distance-dependent manner, driving synergistic enhancer activation. We also discovered that Tead double motifs are widespread, highly active canonical response elements. Molecular dynamics simulations suggest that TEAD4 binds them cooperatively through surprisingly labile protein-protein interactions that depend on the DNA template. These results show that the response to signaling pathways is encoded in the cis-regulatory sequences and that interpreting the rules reveals insights into the mechanisms by which signaling effectors influence cell-type-specific enhancer activity.

Aberrant expression of histone H2B variants reshape chromatin and alter oncogenic gene expression programs

Chromatin architecture governs DNA accessibility and gene expression. Thus, any perturbations to chromatin can significantly alter gene expression programs and promote disease. Prior studies demonstrate that every amino acid in a histone is functionally significant, and that even a single amino acid substitution can drive specific cancers. We previously observed that naturally occurring H2B variants are dysregulated during the epithelial to mesenchymal transition (EMT) in bronchial epithelial cells. Naturally occurring H2B variants differ from canonical H2B by only a few amino acids, yet single amino acid changes in other histone variants (e.g., H3.3) can drive cancer. We therefore hypothesized that H2B variants might function like oncohistones, and investigated how they modify chromatin architecture, dynamics, and function. We find that H2B variants are frequently dysregulated in many cancers, and correlate with patient prognosis. Despite high sequence similarity, mutations in each H2B variant tend to occur at specific "hotspots" in cancer. Some H2B variants cause tighter DNA wrapping around nucleosomes, leading to more compact chromatin structures and reduced transcription factor accessibility to nucleosomal DNA. They also altered genome-wide accessibility to oncogenic regulatory elements and genes, with concomitant changes in oncogenic gene expression programs. Although we did not observe changes in cell proliferation or migration in vitro , our Gene Ontology (GO) analyses of ATAC-seq peaks and RNA-seq data indicated significant changes in oncogenic pathways. These findings suggest that H2B variants may influence early-stage, cancer-associated regulatory mechanisms, potentially setting the stage for oncogenesis later on. Thus, H2B variant expression could serve as an early cancer biomarker, and H2B variants might be novel therapeutic targets.

Comprehensive analysis of computational approaches in plant transcription factors binding regions discovery

Transcription factors (TFs) are regulatory proteins which bind to a specific DNA region known as the transcription factor binding regions (TFBRs) to regulate the rate of transcription process. The identification of TFBRs has been made possible by a number of experimental and computational techniques established during the past few years. The process of TFBR identification involves peak identification in the binding data, followed by the identification of motif characteristics. Using the same binding data attempts have been made to raise computational models to identify such binding regions which could save time and resources spent for binding experiments. These computational approaches depend a lot on what way they learn and how. These existing computational approaches are skewed heavily around human TFBRs discovery, while plants have drastically different genomic setup for regulation which these approaches have grossly ignored. Here, we provide a comprehensive study of the current state of the matters in plant specific TF discovery algorithms. While doing so, we encountered several software tools' issues rendering the tools not useable to researches. We fixed them and have also provided the corrected scripts for such tools. We expect this study to serve as a guide for better understanding of software tools' approaches for plant specific TFBRs discovery and the care to be taken while applying them, especially during cross-species applications. The corrected scripts of these software tools are made available at https://github.com/SCBB-LAB/Comparative-analysis-of-plant-TFBS-software.

Transcriptional silencing in Saccharomyces cerevisiae: known unknowns

Transcriptional silencing in Saccharomyces cerevisiae is a persistent and highly stable form of gene repression. It involves DNA silencers and repressor proteins that bind nucleosomes. The silenced state is influenced by numerous factors including the concentration of repressors, nature of activators, architecture of regulatory elements, modifying enzymes and the dynamics of chromatin.Silencers function to increase the residence time of repressor Sir proteins at silenced domains while clustering of silenced domains enables increased concentrations of repressors and helps facilitate long-range interactions. The presence of an accessible NDR at the regulatory regions of silenced genes, the cycling of chromatin configurations at regulatory sites, the mobility of Sir proteins, and the non-uniform distribution of the Sir proteins across the silenced domain, all result in silenced chromatin that only stably silences weak promoters and enhancers via changes in transcription burst duration and frequency.These data collectively suggest that silencing is probabilistic and the robustness of silencing is achieved through sub-optimization of many different nodes of action such that a stable expression state is generated and maintained even though individual constituents are in constant flux.

PU.1 and BCL11B sequentially cooperate with RUNX1 to anchor mSWI/SNF to poise the T cell effector landscape

Adaptive immunity relies on specialized effector functions elicited by lymphocytes, yet how antigen recognition activates appropriate effector responses through nonspecific signaling intermediates is unclear. Here we examined the role of chromatin priming in specifying the functional outputs of effector T cells and found that most of the cis-regulatory landscape active in effector T cells was poised early in development before the expression of the T cell antigen receptor. We identified two principal mechanisms underpinning this poised landscape: the recruitment of the nucleosome remodeler mammalian SWItch/Sucrose Non-Fermentable (mSWI/SNF) by the transcription factors RUNX1 and PU.1 to establish chromatin accessibility at T effector loci; and a 'relay' whereby the transcription factor BCL11B succeeded PU.1 to maintain occupancy of the chromatin remodeling complex mSWI/SNF together with RUNX1, after PU.1 silencing during lineage commitment. These mechanisms define modes by which T cells acquire the potential to elicit specialized effector functions early in their ontogeny and underscore the importance of integrating extrinsic cues to the developmentally specified intrinsic program.

PML::RARA and GATA2 proteins interact via DNA templates to induce aberrant self-renewal in mouse and human hematopoietic cells

The underlying mechanism(s) by which the PML::RARA fusion protein initiates acute promyelocytic leukemia is not yet clear. We defined the genomic binding sites of PML::RARA in primary mouse and human hematopoietic progenitor cells with V5-tagged PML::RARA, using anti-V5-PML::RARA chromatin immunoprecipitation sequencing and CUT&RUN approaches. Most genomic PML::RARA binding sites were found in regions that were already chromatin-accessible (defined by ATAC-seq) in unmanipulated, wild-type promyelocytes, suggesting that these regions are "open" prior to PML::RARA expression. We found that GATA binding motifs, and the direct binding of the chromatin "pioneering factor" GATA2, were significantly enriched near PML::RARA binding sites. Proximity labeling studies revealed that PML::RARA interacts with ~250 proteins in primary mouse hematopoietic cells; GATA2 and 33 others require PML::RARA binding to DNA for the interaction to occur, suggesting that binding to their cognate DNA target motifs may stabilize their interactions. In the absence of PML::RARA, Gata2 overexpression induces many of the same epigenetic and transcriptional changes as PML::RARA. These findings suggested that PML::RARA may indirectly initiate its transcriptional program by activating Gata2 expression: Indeed, we demonstrated that inactivation of Gata2 prior to PML::RARA expression prevented its ability to induce self-renewal. These data suggested that GATA2 binding creates accessible chromatin regions enriched for both GATA and Retinoic Acid Receptor Element motifs, where GATA2 and PML::RARA can potentially bind and interact with each other. In turn, PML::RARA binding to DNA promotes a feed-forward transcriptional program by positively regulating Gata2 expression. Gata2 may therefore be required for PML::RARA to establish its transcriptional program.

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.

Cis-regulatory landscapes in the evolution and development of the mammalian skull

Extensive morphological variation found in mammals reflects the wide spectrum of their ecological adaptations. The highest morphological diversity is present in the craniofacial region, where geometry is mainly dictated by the bony skull. Mammalian craniofacial development represents complex multistep processes governed by numerous conserved genes that require precise spatio-temporal control. A central question in contemporary evolutionary biology is how a defined set of conserved genes can orchestrate formation of fundamentally different structures, and therefore how morphological variability arises. In principle, differential gene expression patterns during development are the source of morphological variation. With the emergence of multicellular organisms, precise regulation of gene expression in time and space is attributed to cis-regulatory elements. These elements contribute to higher-order chromatin structure and together with trans-acting factors control transcriptional landscapes that underlie intricate morphogenetic processes. Consequently, divergence in cis-regulation is believed to rewire existing gene regulatory networks and form the core of morphological evolution. This review outlines the fundamental principles of the genetic code and genomic regulation interplay during development. Recent work that deepened our comprehension of cis-regulatory element origin, divergence and function is presented here to illustrate the state-of-the-art research that uncovered the principles of morphological novelty. This article is part of the theme issue 'The mammalian skull: development, structure and function'.

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.

Understanding the Target Search by Multiple Transcription Factors on Nucleosomal DNA

The association of multiple Transcription Factors (TFs) in the cis-regulatory region is imperative for developmental changes in eukaryotes. The underlying process is exceedingly complex, and it is not at all clear what orchestrates the overall search process by multiple TFs. In this study, by developing a theoretical model based on a discrete-state stochastic approach, we investigated the target search mechanism of multiple TFs on nucleosomal DNA. Experimental kinetic rate constants of different TFs are taken as input to estimate the Mean-First-Passage time to recognize the binding motifs by two TFs on a dynamic nucleosome model. The theory systematically analyzes when the TFs search their binding motifs hierarchically and when simultaneously by proceeding via the formation of a protein-protein complex. Our results, validated by extensive Monte Carlo simulations, elucidate the molecular basis of the complex target search phenomenon of multiple TFs on nucleosomal DNA.

The Role of PARP1 and PAR in ATP-Independent Nucleosome Reorganisation during the DNA Damage Response

The functioning of the eukaryotic cell genome is mediated by sophisticated protein-nucleic-acid complexes, whose minimal structural unit is the nucleosome. After the damage to genomic DNA, repair proteins need to gain access directly to the lesion; therefore, the initiation of the DNA damage response inevitably leads to local chromatin reorganisation. This review focuses on the possible involvement of PARP1, as well as proteins acting nucleosome compaction, linker histone H1 and non-histone chromatin protein HMGB1. The polymer of ADP-ribose is considered the main regulator during the development of the DNA damage response and in the course of assembly of the correct repair complex.

A quantitative metric of pioneer activity reveals that HNF4A has stronger in vivo pioneer activity than FOXA1

BACKGROUND

We and others have suggested that pioneer activity - a transcription factor's (TF's) ability to bind and open inaccessible loci - is not a qualitative trait limited to a select class of pioneer TFs. We hypothesize that most TFs display pioneering activity that depends on the TF concentration and the motif content at their target loci.

RESULTS

Here, we present a quantitative in vivo measure of pioneer activity that captures the relative difference in a TF's ability to bind accessible versus inaccessible DNA. The metric is based on experiments that use CUT&Tag to measure the binding of doxycycline-inducible TFs. For each location across the genome, we determine the concentration of doxycycline required for a TF to reach half-maximal occupancy; lower concentrations reflect higher affinity. We propose that the relative difference in a TF's affinity between ATAC-seq labeled accessible and inaccessible binding sites is a measure of its pioneer activity. We estimate binding affinities at tens of thousands of genomic loci for the endodermal TFs FOXA1 and HNF4A and show that HNF4A has stronger pioneer activity than FOXA1. We show that both FOXA1 and HNF4A display higher binding affinity at inaccessible sites with more copies of their respective motifs. The quantitative analysis of binding suggests different modes of binding for FOXA1, including an anti-cooperative mode of binding at certain accessible loci.

CONCLUSIONS

Our results suggest that relative binding affinities are reasonable measures of pioneer activity and support the model wherein most TFs have some degree of context-dependent pioneer activity.

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.

Enhancer RNAs step forward: new insights into enhancer function

Enhancers confer precise spatiotemporal patterns of gene expression in response to developmental and environmental stimuli. Over the last decade, the transcription of enhancer RNAs (eRNAs) – nascent RNAs transcribed from active enhancers – has emerged as a key factor regulating enhancer activity. eRNAs are relatively short-lived RNA species that are transcribed at very high rates but also quickly degraded. Nevertheless, eRNAs are deeply intertwined within enhancer regulatory networks and are implicated in a number of transcriptional control mechanisms. Enhancers show changes in function and sequence over evolutionary time, raising questions about the relationship between enhancer sequences and eRNA function. Moreover, the vast majority of single nucleotide polymorphisms associated with human complex diseases map to the non-coding genome, with causal disease variants enriched within enhancers. In this Primer, we survey the diverse roles played by eRNAs in enhancer-dependent gene expression, evaluating different models for eRNA function. We also explore questions surrounding the genetic conservation of enhancers and how this relates to eRNA function and dysfunction.

Summary: This Primer evaluates the ideas that underpin developing models for eRNA function, exploring cases in which perturbed eRNA function contributes to disease.

Stochastic models of nucleosome dynamics reveal regulatory rules of stimulus-induced epigenome remodeling

The genomic positions of nucleosomes are a defining feature of the cell's epigenomic state, but signal-dependent transcription factors (SDTFs), upon activation, bind to specific genomic locations and modify nucleosome positioning. Here we leverage SDTFs as perturbation probes to learn about nucleosome dynamics in living cells. We develop Markov models of nucleosome dynamics and fit them to time course sequencing data of DNA accessibility. We find that (1) the dynamics of DNA unwrapping are significantly slower in cells than reported from cell-free experiments, (2) only models with cooperativity in wrapping and unwrapping fit the available data, (3) SDTF activity produces the highest eviction probability when its binding site is adjacent to but not on the nucleosome dyad, and (4) oscillatory SDTF activity results in high location variability. Our work uncovers the regulatory rules governing SDTF-induced nucleosome dynamics in live cells, which can predict chromatin accessibility alterations during inflammation at single-nucleosome resolution.

Generating specificity in genome regulation through transcription factor sensitivity to chromatin

Cell type-specific gene expression relies on transcription factors (TFs) binding DNA sequence motifs embedded in chromatin. Understanding how motifs are accessed in chromatin is crucial to comprehend differential transcriptional responses and the phenotypic impact of sequence variation. Chromatin obstacles to TF binding range from DNA methylation to restriction of DNA access by nucleosomes depending on their position, composition and modification. In vivo and in vitro approaches now enable the study of TF binding in chromatin at unprecedented resolution. Emerging insights suggest that TFs vary in their ability to navigate chromatin states. However, it remains challenging to link binding and transcriptional outcomes to molecular characteristics of TFs or the local chromatin substrate. Here, we discuss our current understanding of how TFs access DNA in chromatin and novel techniques and directions towards a better understanding of this critical step in genome regulation.

Synthetic multistability in mammalian cells

In multicellular organisms, gene regulatory circuits generate thousands of molecularly distinct, mitotically heritable states through the property of multistability. Designing synthetic multistable circuits would provide insight into natural cell fate control circuit architectures and would allow engineering of multicellular programs that require interactions among distinct cell types. We created MultiFate, a naturally inspired, synthetic circuit that supports long-term, controllable, and expandable multistability in mammalian cells. MultiFate uses engineered zinc finger transcription factors that transcriptionally self-activate as homodimers and mutually inhibit one another through heterodimerization. Using a model-based design, we engineered MultiFate circuits that generate as many as seven states, each stable for at least 18 days. MultiFate permits controlled state switching and modulation of state stability through external inputs and can be expanded with additional transcription factors. These results provide a foundation for engineering multicellular behaviors in mammalian cells.

A test of the pioneer factor hypothesis using ectopic liver gene activation

The pioneer factor hypothesis (PFH) states that pioneer factors (PFs) are a subclass of transcription factors (TFs) that bind to and open inaccessible sites and then recruit non-pioneer factors (non-PFs) that activate batteries of silent genes. The PFH predicts that ectopic gene activation requires the sequential activity of qualitatively different TFs. We tested the PFH by expressing the endodermal PF FOXA1 and non-PF HNF4A in K562 lymphoblast cells. While co-expression of FOXA1 and HNF4A activated a burst of endoderm-specific gene expression, we found no evidence for a functional distinction between these two TFs. When expressed independently, both TFs bound and opened inaccessible sites, activated endodermal genes, and ‘pioneered’ for each other, although FOXA1 required fewer copies of its motif for binding. A subset of targets required both TFs, but the predominant mode of action at these targets did not conform to the sequential activity predicted by the PFH. From these results, we hypothesize an alternative to the PFH where ‘pioneer activity’ depends not on categorically different TFs but rather on the affinity of interaction between TF and DNA.

Cells only use a fraction of their genetic information to make the proteins they need. The rest is carefully packaged away and tightly bundled in structures called nucleosomes. This physically shields the DNA from being accessed by transcription factors – the molecular actors that can read genes and kickstart the protein production process. Effectively, the genetic sequences inside nucleosomes are being silenced.

However, during development, transcription factors must overcome this nucleosome barrier and activate silent genes to program cells. The pioneer factor hypothesis describes how this may be possible: first, ‘pioneer’ transcription factors can bind to and ‘open up’ nucleosomes to make target genes accessible. Then, non-pioneer factors can access the genetic sequence and recruit cofactors that begin copying the now-exposed genetic information.

The widely accepted theory is based on studies of two proteins – FOXA1, an archetypal pioneer factor, and HNF4A, a non-pioneer factor – but the predictions of the pioneer factor hypothesis have yet to be explicitly tested. To do so, Hansen et al. expressed FOXA1 and HNF4A, separately and together, in cells which do not usually make these proteins. They then assessed how the proteins could bind to DNA and impact gene accessibility and transcription.

The experiments demonstrate that FOXA1 and HNF4A do not necessarily follow the two-step activation predicted by the pioneer factor hypothesis. When expressed independently, both transcription factors bound and opened inaccessible sites, activated target genes, and ‘pioneered’ for each other. Similar patterns were observed across the genome. The only notable distinction between the two factors was that FOXA1, the archetypal pioneering factor, required fewer copies of its target sequence to bind DNA than HNF4A.

These findings led Hansen et al. to propose an alternative theory to the pioneer factor hypothesis which eliminates the categorical distinction between pioneer and non-pioneer factors. Overall, this work has implications for how biologists understand the way that transcription factors activate silent genes during development.

Extracting multi-way chromatin contacts from Hi-C data

There is a growing realization that multi-way chromatin contacts formed in chromosome structures are fundamental units of gene regulation. However, due to the paucity and complexity of such contacts, it is challenging to detect and identify them using experiments. Based on an assumption that chromosome structures can be mapped onto a network of Gaussian polymer, here we derive analytic expressions for n-body contact probabilities (n > 2) among chromatin loci based on pairwise genomic contact frequencies available in Hi-C, and show that multi-way contact probability maps can in principle be extracted from Hi-C. The three-body (triplet) contact probabilities, calculated from our theory, are in good correlation with those from measurements including Tri-C, MC-4C and SPRITE. Maps of multi-way chromatin contacts calculated from our analytic expressions can not only complement experimental measurements, but also can offer better understanding of the related issues, such as cell-line dependent assemblies of multiple genes and enhancers to chromatin hubs, competition between long-range and short-range multi-way contacts, and condensates of multiple CTCF anchors.

Role of Enhancers in Development and Diseases

Enhancers are cis-regulatory elements containing short DNA sequences that serve as binding sites for pioneer/regulatory transcription factors, thus orchestrating the regulation of genes critical for lineage determination. The activity of enhancer elements is believed to be determined by transcription factor binding, thus determining the cell state identity during development. Precise spatio-temporal control of the transcriptome during lineage specification requires the coordinated binding of lineage-specific transcription factors to enhancers. Thus, enhancers are the primary determinants of cell identity. Numerous studies have explored the role and mechanism of enhancers during development and disease, and various basic questions related to the functions and mechanisms of enhancers have not yet been fully answered. In this review, we discuss the recently published literature regarding the roles of enhancers, which are critical for various biological processes governing development. Furthermore, we also highlight that altered enhancer landscapes provide an essential context to understand the etiologies and mechanisms behind numerous complex human diseases, providing new avenues for effective enhancer-based therapeutic interventions.

The transcription factor code in iPSC reprogramming

Transcription factor (TF)-induced reprogramming of somatic cells across lineages and to induced pluripotent stem cells (iPSCs) has revealed a remarkable plasticity of differentiated cells and presents great opportunities for generating clinically relevant cell types for disease modeling and regenerative medicine. The understanding of iPSC reprogramming provides insights into the mechanisms that safeguard somatic cell identity, drive epigenetic reprogramming, and underlie cell fate specification in vivo. The combinatorial action of TFs has emerged as the key mechanism for the direct and indirect effects of reprogramming factors that induce the remodelling of the enhancer landscape. The interplay of TFs in iPSC reprogramming also yields trophectoderm- and extraembryonic endoderm-like cell populations, uncovering an intriguing plasticity of cell states and opening new avenues for exploring cell fate decisions during early embryogenesis.

NF-κB dynamics determine the stimulus specificity of epigenomic reprogramming in macrophages

The epigenome of macrophages can be reprogrammed by extracellular cues, but the extent to which different stimuli achieve this is unclear. Nuclear factor κB (NF-κB) is a transcription factor that is activated by all pathogen-associated stimuli and can reprogram the epigenome by activating latent enhancers. However, we show that NF-κB does so only in response to a subset of stimuli. This stimulus specificity depends on the temporal dynamics of NF-κB activity, in particular whether it is oscillatory or non-oscillatory. Non-oscillatory NF-κB opens chromatin by sustained disruption of nucleosomal histone-DNA interactions, enabling activation of latent enhancers that modulate expression of immune response genes. Thus, temporal dynamics can determine a transcription factor's capacity to reprogram the epigenome in a stimulus-specific manner.

Allosteric conformational ensembles have unlimited capacity for integrating information

Integration of binding information by macromolecular entities is fundamental to cellular functionality. Recent work has shown that such integration cannot be explained by pairwise cooperativities, in which binding is modulated by binding at another site. Higher-order cooperativities (HOCs), in which binding is collectively modulated by multiple other binding events, appear to be necessary but an appropriate mechanism has been lacking. We show here that HOCs arise through allostery, in which effective cooperativity emerges indirectly from an ensemble of dynamically interchanging conformations. Conformational ensembles play important roles in many cellular processes but their integrative capabilities remain poorly understood. We show that sufficiently complex ensembles can implement any form of information integration achievable without energy expenditure, including all patterns of HOCs. Our results provide a rigorous biophysical foundation for analysing the integration of binding information through allostery. We discuss the implications for eukaryotic gene regulation, where complex conformational dynamics accompanies widespread information integration.

Gene Regulation in and out of Equilibrium

Determining whether and how a gene is transcribed are two of the central processes of life. The conceptual basis for understanding such gene regulation arose from pioneering biophysical studies in eubacteria. However, eukaryotic genomes exhibit vastly greater complexity, which raises questions not addressed by this bacterial paradigm. First, how is information integrated from many widely separated binding sites to determine how a gene is transcribed? Second, does the presence of multiple energy-expending mechanisms, which are absent from eubacterial genomes, indicate that eukaryotes are capable of improved forms of genetic information processing? An updated biophysical foundation is needed to answer such questions. We describe the linear framework, a graph-based approach to Markov processes, and show that it can accommodate many previous studies in the field. Under the assumption of thermodynamic equilibrium, we introduce a language of higher-order cooperativities and show how it can rigorously quantify gene regulatory properties suggested by experiment. We point out that fundamental limits to information processing arise at thermodynamic equilibrium and can only be bypassed through energy expenditure. Finally, we outline some of the mathematical challenges that must be overcome to construct an improved biophysical understanding of gene regulation.

Molecular Co-occupancy Identifies Transcription Factor Binding Cooperativity In Vivo

Gene activation requires the cooperative activity of multiple transcription factors at cis-regulatory elements (CREs). Yet, most transcription factors have short residence time, questioning the requirement of their physical co-occupancy on DNA to achieve cooperativity. Here, we present a DNA footprinting method that detects individual molecular interactions of transcription factors and nucleosomes with DNA in vivo. We apply this strategy to quantify the simultaneous binding of multiple transcription factors on single DNA molecules at mouse CREs. Analysis of the binary occupancy patterns at thousands of motif combinations reveals that high DNA co-occupancy occurs for most types of transcription factors, in the absence of direct physical interaction, at sites of competition with nucleosomes. Perturbation of pairwise interactions demonstrates the function of molecular co-occupancy in binding cooperativity. Our results reveal the interactions regulating CREs at molecular resolution and identify DNA co-occupancy as a widespread cooperativity mechanism used by transcription factors to remodel chromatin.

Quantitative dissection of transcription in development yields evidence for transcription-factor-driven chromatin accessibility

Thermodynamic models of gene regulation can predict transcriptional regulation in bacteria, but in eukaryotes, chromatin accessibility and energy expenditure may call for a different framework. Here, we systematically tested the predictive power of models of DNA accessibility based on the Monod-Wyman-Changeux (MWC) model of allostery, which posits that chromatin fluctuates between accessible and inaccessible states. We dissected the regulatory dynamics of hunchback by the activator Bicoid and the pioneer-like transcription factor Zelda in living Drosophila embryos and showed that no thermodynamic or non-equilibrium MWC model can recapitulate hunchback transcription. Therefore, we explored a model where DNA accessibility is not the result of thermal fluctuations but is catalyzed by Bicoid and Zelda, possibly through histone acetylation, and found that this model can predict hunchback dynamics. Thus, our theory-experiment dialogue uncovered potential molecular mechanisms of transcriptional regulatory dynamics, a key step toward reaching a predictive understanding of developmental decision-making.

Genome-wide off-rates reveal how DNA binding dynamics shape transcription factor function

Protein-DNA interactions are dynamic, and these dynamics are an important aspect of chromatin-associated processes such as transcription or replication. Due to a lack of methods to study on- and off-rates across entire genomes, protein-DNA interaction dynamics have not been studied extensively. Here, we determine in vivo off-rates for the Saccharomyces cerevisiae chromatin organizing factor Abf1, at 191 sites simultaneously across the yeast genome. Average Abf1 residence times span a wide range, varying between 4.2 and 33 min. Sites with different off-rates are associated with different functional characteristics. This includes their transcriptional dependency on Abf1, nucleosome positioning and the size of the nucleosome-free region, as well as the ability to roadblock RNA polymerase II for termination. The results show how off-rates contribute to transcription factor function and that DIVORSEQ (Determining In Vivo Off-Rates by SEQuencing) is a meaningful way of investigating protein-DNA binding dynamics genome-wide.

Single-molecule regulatory architectures captured by chromatin fiber sequencing

Gene regulation is chiefly determined at the level of individual linear chromatin molecules, yet our current understanding of cis-regulatory architectures derives from fragmented sampling of large numbers of disparate molecules. We developed an approach for precisely stenciling the structure of individual chromatin fibers onto their composite DNA templates using nonspecific DNA N⁶-adenine methyltransferases. Single-molecule long-read sequencing of chromatin stencils enabled nucleotide-resolution readout of the primary architecture of multikilobase chromatin fibers (Fiber-seq). Fiber-seq exposed widespread plasticity in the linear organization of individual chromatin fibers and illuminated principles guiding regulatory DNA actuation, the coordinated actuation of neighboring regulatory elements, single-molecule nucleosome positioning, and single-molecule transcription factor occupancy. Our approach and results open new vistas on the primary architecture of gene regulation.

Sequence and chromatin determinants of transcription factor binding and the establishment of cell type-specific binding patterns

Transcription factors (TFs) selectively bind distinct sets of sites in different cell types. Such cell type-specific binding specificity is expected to result from interplay between the TF's intrinsic sequence preferences, cooperative interactions with other regulatory proteins, and cell type-specific chromatin landscapes. Cell type-specific TF binding events are highly correlated with patterns of chromatin accessibility and active histone modifications in the same cell type. However, since concurrent chromatin may itself be a consequence of TF binding, chromatin landscapes measured prior to TF activation provide more useful insights into how cell type-specific TF binding events became established in the first place. Here, we review the various sequence and chromatin determinants of cell type-specific TF binding specificity. We identify the current challenges and opportunities associated with computational approaches to characterizing, imputing, and predicting cell type-specific TF binding patterns. We further focus on studies that characterize TF binding in dynamic regulatory settings, and we discuss how these studies are leading to a more complex and nuanced understanding of dynamic protein-DNA binding activities. We propose that TF binding activities at individual sites can be viewed along a two-dimensional continuum of local sequence and chromatin context. Under this view, cell type-specific TF binding activities may result from either strongly favorable sequence features or strongly favorable chromatin context.

Developmental Transcriptional Enhancers: A Subtle Interplay between Accessibility and Activity: Considering Quantitative Accessibility Changes between Different Regulatory States of an Enhancer Deconvolutes the Complex Relationship between Accessibility and Activity

Measurements of open chromatin in specific cell types are widely used to infer the spatiotemporal activity of transcriptional enhancers. How reliable are these predictions? In this review, it is argued that the relationship between the accessibility and activity of an enhancer is insufficiently described by simply considering open versus closed chromatin, or active versus inactive enhancers. Instead, recent studies focusing on the quantitative nature of accessibility signal reveal subtle differences between active enhancers and their different inactive counterparts: the closed silenced state and the accessible primed and repressed states. While the open structure as such is not a specific indicator of enhancer activity, active enhancers display a higher degree of accessibility than the primed and repressed states. Molecular mechanisms that may account for these quantitative differences are discussed. A model that relates molecular events at an enhancer to changes in its activity and accessibility in a developing tissue is also proposed.

Improved recovery of cell-cycle gene expression in Saccharomyces cerevisiae from regulatory interactions in multiple omics data

BACKGROUND

Gene expression is regulated by DNA-binding transcription factors (TFs). Together with their target genes, these factors and their interactions collectively form a gene regulatory network (GRN), which is responsible for producing patterns of transcription, including cyclical processes such as genome replication and cell division. However, identifying how this network regulates the timing of these patterns, including important interactions and regulatory motifs, remains a challenging task.

RESULTS

We employed four in vivo and in vitro regulatory data sets to investigate the regulatory basis of expression timing and phase-specific patterns cell-cycle expression in Saccharomyces cerevisiae. Specifically, we considered interactions based on direct binding between TF and target gene, indirect effects of TF deletion on gene expression, and computational inference. We found that the source of regulatory information significantly impacts the accuracy and completeness of recovering known cell-cycle expressed genes. The best approach involved combining TF-target and TF-TF interactions features from multiple datasets in a single model. In addition, TFs important to multiple phases of cell-cycle expression also have the greatest impact on individual phases. Important TFs regulating a cell-cycle phase also tend to form modules in the GRN, including two sub-modules composed entirely of unannotated cell-cycle regulators (STE12-TEC1 and RAP1-HAP1-MSN4).

CONCLUSION

Our findings illustrate the importance of integrating both multiple omics data and regulatory motifs in order to understand the significance regulatory interactions involved in timing gene expression. This integrated approached allowed us to recover both known cell-cycles interactions and the overall pattern of phase-specific expression across the cell-cycle better than any single data set. Likewise, by looking at regulatory motifs in the form of TF-TF interactions, we identified sets of TFs whose co-regulation of target genes was important for cell-cycle expression, even when regulation by individual TFs was not. Overall, this demonstrates the power of integrating multiple data sets and models of interaction in order to understand the regulatory basis of established biological processes and their associated gene regulatory networks.

CTCF confers local nucleosome resiliency after DNA replication and during mitosis

The access of Transcription Factors (TFs) to their cognate DNA binding motifs requires a precise control over nucleosome positioning. This is especially important following DNA replication and during mitosis, both resulting in profound changes in nucleosome organization over TF binding regions. Using mouse Embryonic Stem (ES) cells, we show that the TF CTCF displaces nucleosomes from its binding site and locally organizes large and phased nucleosomal arrays, not only in interphase steady-state but also immediately after replication and during mitosis. Correlative analyses suggest this is associated with fast gene reactivation following replication and mitosis. While regions bound by other TFs (Oct4/Sox2), display major rearrangement, the post-replication and mitotic nucleosome positioning activity of CTCF is not unique: Esrrb binding regions are also characterized by persistent nucleosome positioning. Therefore, selected TFs such as CTCF and Esrrb act as resilient TFs governing the inheritance of nucleosome positioning at regulatory regions throughout the cell-cycle.

Identification of DNA motifs that regulate DNA methylation

DNA methylation is an important epigenetic mark but how its locus-specificity is decided in relation to DNA sequence is not fully understood. Here, we have analyzed 34 diverse whole-genome bisulfite sequencing datasets in human and identified 313 motifs, including 92 and 221 associated with methylation (methylation motifs, MMs) and unmethylation (unmethylation motifs, UMs), respectively. The functionality of these motifs is supported by multiple lines of evidence. First, the methylation levels at the MM and UM motifs are respectively higher and lower than the genomic background. Second, these motifs are enriched at the binding sites of methylation modifying enzymes including DNMT3A and TET1, indicating their possible roles of recruiting these enzymes. Third, these motifs significantly overlap with "somatic QTLs" (quantitative trait loci) of methylation and expression. Fourth, disruption of these motifs by mutation is associated with significantly altered methylation level of the CpGs in the neighbor regions. Furthermore, these motifs together with somatic mutations are predictive of cancer subtypes and patient survival. We revealed some of these motifs were also associated with histone modifications, suggesting a possible interplay between the two types of epigenetic modifications. We also found some motifs form feed forward loops to contribute to DNA methylation dynamics.

Chromatin accessibility and the regulatory epigenome

Physical access to DNA is a highly dynamic property of chromatin that plays an essential role in establishing and maintaining cellular identity. The organization of accessible chromatin across the genome reflects a network of permissible physical interactions through which enhancers, promoters, insulators and chromatin-binding factors cooperatively regulate gene expression. This landscape of accessibility changes dynamically in response to both external stimuli and developmental cues, and emerging evidence suggests that homeostatic maintenance of accessibility is itself dynamically regulated through a competitive interplay between chromatin-binding factors and nucleosomes. In this Review, we examine how the accessible genome is measured and explore the role of transcription factors in initiating accessibility remodelling; our goal is to illustrate how chromatin accessibility defines regulatory elements within the genome and how these epigenetic features are dynamically established to control gene expression.

An Optogenetic Platform for Real-Time, Single-Cell Interrogation of Stochastic Transcriptional Regulation

Transcription is a highly regulated and inherently stochastic process. The complexity of signal transduction and gene regulation makes it challenging to analyze how the dynamic activity of transcriptional regulators affects stochastic transcription. By combining a fast-acting, photo-regulatable transcription factor with nascent RNA quantification in live cells and an experimental setup for precise spatiotemporal delivery of light inputs, we constructed a platform for the real-time, single-cell interrogation of transcription in Saccharomyces cerevisiae. We show that transcriptional activation and deactivation are fast and memoryless. By analyzing the temporal activity of individual cells, we found that transcription occurs in bursts, whose duration and timing are modulated by transcription factor activity. Using our platform, we regulated transcription via light-driven feedback loops at the single-cell level. Feedback markedly reduced cell-to-cell variability and led to qualitative differences in cellular transcriptional dynamics. Our platform establishes a flexible method for studying transcriptional dynamics in single cells.

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.

Enhancers and chromatin structures: regulatory hubs in gene expression and diseases

Gene expression requires successful communication between enhancer and promoter regions, whose activities are regulated by a variety of factors and associated with distinct chromatin structures; in addition, functionally related genes and their regulatory repertoire tend to be arranged in the same subchromosomal regulatory domains. In this review, we discuss the importance of enhancers, especially clusters of enhancers (such as super-enhancers), as key regulatory hubs to integrate environmental cues and encode spatiotemporal instructions for genome expression, which are critical for a variety of biological processes governing mammalian development. Furthermore, we emphasize that the enhancer–promoter interaction landscape provides a critical context to understand the aetiologies and mechanisms behind numerous complex human diseases and provides new avenues for effective transcription-based interventions.

Competition between histone and transcription factor binding regulates the onset of transcription in zebrafish embryos

Upon fertilization, the genome of animal embryos remains transcriptionally inactive until the maternal-to-zygotic transition. At this time, the embryo takes control of its development and transcription begins. How the onset of zygotic transcription is regulated remains unclear. Here, we show that a dynamic competition for DNA binding between nucleosome-forming histones and transcription factors regulates zebrafish genome activation. Taking a quantitative approach, we found that the concentration of non-DNA-bound core histones sets the time for the onset of transcription. The reduction in nuclear histone concentration that coincides with genome activation does not affect nucleosome density on DNA, but allows transcription factors to compete successfully for DNA binding. In agreement with this, transcription factor binding is sensitive to histone levels and the concentration of transcription factors also affects the time of transcription. Our results demonstrate that the relative levels of histones and transcription factors regulate the onset of transcription in the embryo.

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

The DNA in a fertilized egg contains all the information required to form an animal’s body. In order for the animal to develop properly, particular genes encoded in the DNA are only active at specific times. The DNA is wrapped around proteins called histones, which allows the DNA to be tightly packed inside the cell. However, histones can block other proteins called transcription factors from binding to the DNA to activate the genes. Young embryos initially develop with all of their genes switched off, relying on the nutrients and other molecules provided by their mother. After some time, the embryo starts to switch on its own genes to take control of its own development, but it was not clear how this happens.

Joseph et al. investigated how genes are activated in zebrafish embryos, which are often used as models to study how animals develop. The experiments show that competition between histones and transcription factors for binding to DNA controls when genes are switched on. In young fish embryos, there are so many histones present that transcription factors have no opportunity to bind to DNA. Over time, however, the numbers of histones decrease, allowing transcription factors to bind to DNA and switch on genes.

Histones and transcription factors regulate the activity of genes throughout the life of the animal. Therefore, competition between these two types of protein may also control gene activity in other situations. A better understanding of how gene activity is controlled could allow researchers to more easily grow different types of cell in the laboratory or to reprogram specific cells in the body. As such, these new findings may aid the development of therapies to regenerate organs or tissues that have been damaged by injury or disease.

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

Structural perspective of cooperative transcription factor binding

In prokaryotes, individual transcription factors (TFs) can recognize long DNA motifs that are alone sufficient to define the genes that they induce or repress. In contrast, in higher organisms that have larger genomes, TFs recognize sequences that are too short to define unique genomic positions. In addition, development of multicellular organisms requires molecular systems that are capable of executing combinatorial logical operations. Co-operative recognition of DNA by multiple TFs allows both definition of unique genomic positions in large genomes, and complex information processing at the level of individual regulatory elements. The TFs can co-operate in multiple different ways, and the precise mechanism used for co-operation determines important features of the regulatory interactions. Here, we present an overview of the structural basis of the different mechanisms by which TFs can cooperate, focusing on insight from recent functional studies and structural analyses of specific TF-TF-DNA complexes.

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.

Epigenetic regulation of cellular and cytomegalovirus genes during myeloid cell development

Myeloid cells are important cell types that carry human cytomegalovirus. Latent viral DNA is present in CD34+ progenitor cells and their derived monocytes. However, differentiation of latently infected monocytes to mature macrophages or dendritic cells causes reactivation of latent viruses. During hematopoietic development, pluripotent genes are repressed, and lineage specific genes are activated in a step-wise manner. This process is governed by cell-type specific chromatin states. Enhancers in the hematopoietic system are highly dynamic and established by pioneer (first tier) transcription factors (TFs), which set the stage for second and third tier TF binding. In this review, we examine the epigenetic mechanisms that regulate myeloid cell development, cell identity, and activation with a special focus on factors that regulate viral gene expression and the status of viral infection in myeloid cells.

Selecting optimal combinations of transcription factors to promote axon regeneration: Why mechanisms matter

Recovery from injuries to the central nervous system, including spinal cord injury, is constrained in part by the intrinsically low ability of many CNS neurons to mount an effective regenerative growth response. To improve outcomes, it is essential to understand and ultimately reverse these neuron-intrinsic constraints. Genetic manipulation of key transcription factors (TFs), which act to orchestrate production of multiple regeneration-associated genes, has emerged as a promising strategy. It is likely that no single TF will be sufficient to fully restore neuron-intrinsic growth potential, and that multiple, functionally interacting factors will be needed. An extensive literature, mostly from non-neural cell types, has identified potential mechanisms by which TFs can functionally synergize. Here we examine four potential mechanisms of TF/TF interaction; physical interaction, transcriptional cross-regulation, signaling-based cross regulation, and co-occupancy of regulatory DNA. For each mechanism, we consider how existing knowledge can be used to guide the discovery and effective use of TF combinations in the context of regenerative neuroscience. This mechanistic insight into TF interactions is needed to accelerate the design of effective TF-based interventions to relieve neuron-intrinsic constraints to regeneration and to foster recovery from CNS injury.

Ever-Changing Landscapes: Transcriptional Enhancers in Development and Evolution

A class of cis-regulatory elements, called enhancers, play a central role in orchestrating spatiotemporally precise gene-expression programs during development. Consequently, divergence in enhancer sequence and activity is thought to be an important mediator of inter- and intra-species phenotypic variation. Here, we give an overview of emerging principles of enhancer function, current models of enhancer architecture, genomic substrates from which enhancers emerge during evolution, and the influence of three-dimensional genome organization on long-range gene regulation. We discuss intricate relationships between distinct elements within complex regulatory landscapes and consider their potential impact on specificity and robustness of transcriptional regulation.

Nucleosome dynamics and maintenance of epigenetic states of CpG islands

Methylation of mammalian DNA occurs primarily at CG dinucleotides. These CpG sites are located nonrandomly in the genome, tending to occur within high density clusters of CpGs (islands) or within large regions of low CpG density. Cluster methylation tends to be bimodal, being dominantly unmethylated or mostly methylated. For CpG clusters near promoters, low methylation is associated with transcriptional activity, while high methylation is associated with gene silencing. Alternative CpG methylation states are thought to be stable and heritable, conferring localized epigenetic memory that allows transient signals to create long-lived gene expression states. Positive feedback where methylated CpG sites recruit enzymes that methylate nearby CpGs, can produce heritable bistability but does not easily explain that as clusters increase in size or density they change from being primarily methylated to primarily unmethylated. Here, we show that an interaction between the methylation state of a cluster and its occupancy by nucleosomes provides a mechanism to generate these features and explain genome wide systematics of CpG islands.

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

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

Roles of cofactors and chromatin accessibility in Hox protein target specificity

Background

The regulation of specific target genes by transcription factors is central to our understanding of gene network control in developmental and physiological processes yet how target specificity is achieved is still poorly understood. This is well illustrated by the Hox family of transcription factors as their limited in vitro DNA-binding specificity contrasts with their clear in vivo functional specificity.

Results

We generated genome-wide binding profiles for three Hox proteins, Ubx, Abd-A and Abd-B, following transient expression in Drosophila Kc167 cells, revealing clear target specificity and a striking influence of chromatin accessibility. In the absence of the TALE class homeodomain cofactors Exd and Hth, Ubx and Abd-A bind at a very similar set of target sites in accessible chromatin, whereas Abd-B binds at an additional specific set of targets. Provision of Hox cofactors Exd and Hth considerably modifies the Ubx genome-wide binding profile enabling Ubx to bind at an additional novel set of targets. Both the Abd-B specific targets and the cofactor-dependent Ubx targets are in chromatin that is relatively DNase1 inaccessible prior to the expression of Hox proteins/Hox cofactors.

Conclusions

Our experiments demonstrate a strong role for chromatin accessibility in Hox protein binding and suggest that Hox protein competition with nucleosomes has a major role in Hox protein target specificity in vivo.

Electronic supplementary material

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