Molecular Co-occupancy Identifies Transcription Factor Binding Cooperativity In Vivo

cBAF generates subnucleosomes that expand OCT4 binding and function beyond DNA motifs at enhancers

The canonical BRG/BRM-associated factor (cBAF) complex is essential for chromatin opening at enhancers in mammalian cells. However, the nature of the open chromatin remains unclear. Here, we show that, in addition to producing histone-free DNA, cBAF generates stable hemisome-like subnucleosomal particles containing the four core histones associated with 50-80 bp of DNA. Our genome-wide analysis indicates that cBAF makes these particles by targeting and splitting fragile nucleosomes. In mouse embryonic stem cells, these subnucleosomes become an in vivo binding substrate for the master transcription factor OCT4 independently of the presence of OCT4 DNA motifs. At enhancers, the OCT4-subnucleosome interaction increases OCT4 occupancy and amplifies the genomic interval bound by OCT4 by up to one order of magnitude compared to the region occupied on histone-free DNA. We propose that cBAF-dependent subnucleosomes orchestrate a molecular mechanism that projects OCT4 function in chromatin opening beyond its DNA motifs.

Systematic assessment of ISWI subunits shows that NURF creates local accessibility for CTCF

Catalytic activity of the imitation switch (ISWI) family of remodelers is critical for nucleosomal organization and DNA binding of certain transcription factors, including the insulator protein CTCF. Here we define the contribution of individual subcomplexes by deriving a panel of isogenic mouse stem cell lines, each lacking one of six ISWI accessory subunits. Individual deletions of subunits of either CERF, RSF, ACF, WICH or NoRC subcomplexes only moderately affect the chromatin landscape, while removal of the NURF-specific subunit BPTF leads to a strong reduction in chromatin accessibility and SNF2H ATPase localization around CTCF sites. This affects adjacent nucleosome occupancy and CTCF binding. At a group of sites with reduced chromatin accessibility, CTCF binding persists but cohesin occupancy is reduced, resulting in decreased insulation. These results suggest that CTCF binding can be separated from its function as an insulator in nuclear organization and identify a specific role for NURF in mediating SNF2H localization and chromatin opening at bound CTCF sites.

Transfer learning reveals sequence determinants of the quantitative response to transcription factor dosage

Deep learning approaches have made significant advances in predicting cell type-specific chromatin patterns from the identity and arrangement of transcription factor (TF) binding motifs. However, most models have been applied in unperturbed contexts, precluding a predictive understanding of how chromatin state responds to TF perturbation. Here, we used transfer learning to train and interpret deep learning models that use DNA sequence to predict, with accuracy approaching experimental reproducibility, how the concentration of two dosage-sensitive TFs (TWIST1, SOX9) affects regulatory element (RE) chromatin accessibility in facial progenitor cells. High-affinity motifs that allow for heterotypic TF co-binding and are concentrated at the center of REs buffer against quantitative changes in TF dosage and strongly predict unperturbed accessibility. In contrast, motifs with low-affinity or homotypic binding distributed throughout REs lead to sensitive responses with minimal contributions to unperturbed accessibility. Both buffering and sensitizing features show signatures of purifying selection. We validated these predictive sequence features using reporter assays and showed that a biophysical model of TF-nucleosome competition can explain the sensitizing effect of low-affinity motifs. Our approach of combining transfer learning and quantitative measurements of the chromatin response to TF dosage therefore represents a powerful method to reveal additional layers of the cis-regulatory code.

The Hill function is the universal Hopfield barrier for sharpness of input-output responses

The Hill functions, [Formula: see text], have been widely used in biology for over a century but, with the exception of [Formula: see text], they have had no justification other than as a convenient fit to empirical data. Here, we show that they are the universal limit for the sharpness of any input-output response arising from a Markov process model at thermodynamic equilibrium. Models may represent arbitrary molecular complexity, with multiple ligands, internal states, conformations, coregulators, etc, under core assumptions that are detailed in the paper. The model output may be any linear combination of steady-state probabilities, with components other than the chosen input ligand held constant. This formulation generalizes most of the responses in the literature. We use a coarse-graining method in the graph-theoretic linear framework to show that two sharpness measures for input-output responses fall within an effectively bounded region of the positive quadrant, [Formula: see text], for any equilibrium model with [Formula: see text] input binding sites. [Formula: see text] exhibits a cusp which approaches, but never exceeds, the sharpness of [Formula: see text], but the region and the cusp can be exceeded when models are taken away from thermodynamic equilibrium. Such fundamental thermodynamic limits are called Hopfield barriers, and our results provide a biophysical justification for the Hill functions as the universal Hopfield barriers for sharpness. Our results also introduce an object, [Formula: see text], whose structure may be of mathematical interest, and suggest the importance of characterizing Hopfield barriers for other forms of cellular information processing.

Dynamic modes of Notch transcription hubs conferring memory and stochastic activation revealed by live imaging the co-activator Mastermind

Developmental programming involves the accurate conversion of signalling levels and dynamics to transcriptional outputs. The transcriptional relay in the Notch pathway relies on nuclear complexes containing the co-activator Mastermind (Mam). By tracking these complexes in real time, we reveal that they promote the formation of a dynamic transcription hub in Notch ON nuclei which concentrates key factors including the Mediator CDK module. The composition of the hub is labile and persists after Notch withdrawal conferring a memory that enables rapid reformation. Surprisingly, only a third of Notch ON hubs progress to a state with nascent transcription, which correlates with polymerase II and core Mediator recruitment. This probability is increased by a second signal. The discovery that target-gene transcription is probabilistic has far-reaching implications because it implies that stochastic differences in Notch pathway output can arise downstream of receptor activation.

Epigenetic reprogramming in the transition from pluripotency to totipotency

Mammalian development commences with the zygote, which can differentiate into both embryonic and extraembryonic tissues, a capability known as totipotency. Only the zygote and embryos around zygotic genome activation (ZGA) (two-cell embryo stage in mice and eight-cell embryo in humans) are totipotent cells. Epigenetic modifications undergo extremely extensive changes during the acquisition of totipotency and subsequent development of differentiation. However, the underlying molecular mechanisms remain elusive. Recently, the discovery of mouse two-cell embryo-like cells, human eight-cell embryo-like cells, extended pluripotent stem cells and totipotent-like stem cells with extra-embryonic developmental potential has greatly expanded our understanding of totipotency. Experiments with these in vitro models have led to insights into epigenetic changes in the reprogramming of pluri-to-totipotency, which have informed the exploration of preimplantation development. In this review, we highlight the recent findings in understanding the mechanisms of epigenetic remodeling during totipotency capture, including RNA splicing, DNA methylation, chromatin configuration, histone modifications, and nuclear organization.

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.

The Hill function is the universal Hopfield barrier for sharpness of input-output responses

The Hill functions, ℋ h ( x ) = x h / 1 + x h , have been widely used in biology for over a century but, with the exception of ℋ 1 , they have had no justification other than as a convenient fit to empirical data. Here, we show that they are the universal limit for the sharpness of any input-output response arising from a Markov process model at thermodynamic equilibrium. Models may represent arbitrary molecular complexity, with multiple ligands, internal states, conformations, co-regulators, etc, under core assumptions that are detailed in the paper. The model output may be any linear combination of steady-state probabilities, with components other than the chosen input ligand held constant. This formulation generalises most of the responses in the literature. We use a coarse-graining method in the graph-theoretic linear framework to show that two sharpness measures for input-output responses fall within an effectively bounded region of the positive quadrant, Ω m ⊂ ℝ + 2 , for any equilibrium model with m input binding sites. Ω m exhibits a cusp which approaches, but never exceeds, the sharpness of ℋ m but the region and the cusp can be exceeded when models are taken away from thermodynamic equilibrium. Such fundamental thermodynamic limits are called Hopfield barriers and our results provide a biophysical justification for the Hill functions as the universal Hopfield barriers for sharpness. Our results also introduce an object, Ω m , whose structure may be of mathematical interest, and suggest the importance of characterising Hopfield barriers for other forms of cellular information processing.

Transcription factor exchange enables prolonged transcriptional bursts

Single-molecule imaging inside living cells has revealed that transcription factors (TFs) bind to DNA transiently, but a long-standing question is how this transient binding is related to transcription activation. Here, we devised a microscopy method to simultaneously measure transient TF binding at a single locus and the effect of these binding events on transcription. We show that DNA binding of the yeast TF Gal4 activates transcription of a target gene within a few seconds, with at least ∼20% efficiency and with a high initiation rate of ∼1 RNA/s. Gal4 DNA dissociation decreases transcription rapidly. Moreover, at a gene with multiple binding sites, individual Gal4 molecules only rarely stay bound throughout the entire burst but instead frequently exchange during a burst to increase the transcriptional burst duration. Our results suggest a mechanism for enhancer regulation in more complex eukaryotes, where TF cooperativity and exchange enable robust and responsive transcription regulation.

Single-molecule chromatin configurations link transcription factor binding to expression in human cells

The binding of multiple transcription factors (TFs) to genomic enhancers activates gene expression in mammalian cells. However, the molecular details that link enhancer sequence to TF binding, promoter state, and gene expression levels remain opaque. We applied single-molecule footprinting (SMF) to measure the simultaneous occupancy of TFs, nucleosomes, and components of the transcription machinery on engineered enhancer/promoter constructs with variable numbers of TF binding sites for both a synthetic and an endogenous TF. We find that activation domains enhance a TF's capacity to compete with nucleosomes for binding to DNA in a BAF-dependent manner, TF binding on nucleosome-free DNA is consistent with independent binding between TFs, and average TF occupancy linearly contributes to promoter activation rates. We also decompose TF strength into separable binding and activation terms, which can be tuned and perturbed independently. Finally, we develop thermodynamic and kinetic models that quantitatively predict both the binding microstates observed at the enhancer and subsequent time-dependent gene expression. This work provides a template for quantitative dissection of distinct contributors to gene activation, including the activity of chromatin remodelers, TF activation domains, chromatin acetylation, TF concentration, TF binding affinity, and TF binding site configuration.

Time will tell: comparing timescales to gain insight into transcriptional bursting

Recent imaging studies have captured the dynamics of regulatory events of transcription inside living cells. These events include transcription factor (TF) DNA binding, chromatin remodeling and modification, enhancer-promoter (E-P) proximity, cluster formation, and preinitiation complex (PIC) assembly. Together, these molecular events culminate in stochastic bursts of RNA synthesis, but their kinetic relationship remains largely unclear. In this review, we compare the timescales of upstream regulatory steps (input) with the kinetics of transcriptional bursting (output) to generate mechanistic models of transcription dynamics in single cells. We highlight open questions and potential technical advances to guide future endeavors toward a quantitative and kinetic understanding of transcription regulation.

A single fiber view of the nucleosome organization in eukaryotic chromatin

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

The BAF chromatin remodeler synergizes with RNA polymerase II and transcription factors to evict nucleosomes

Chromatin accessibility is a hallmark of active transcription and entails ATP-dependent nucleosome remodeling, which is carried out by complexes such as Brahma-associated factor (BAF). However, the mechanistic links between transcription, nucleosome remodeling and chromatin accessibility are unclear. Here, we used a chemical-genetic approach coupled with time-resolved chromatin profiling to dissect the interplay between RNA Polymerase II (RNAPII), BAF and DNA-sequence-specific transcription factors in mouse embryonic stem cells. We show that BAF dynamically unwraps and evicts nucleosomes at accessible chromatin regions, while RNAPII promoter-proximal pausing stabilizes BAF chromatin occupancy and enhances ATP-dependent nucleosome eviction by BAF. We find that although RNAPII and BAF dynamically probe both transcriptionally active and Polycomb-repressed genomic regions, pluripotency transcription factor chromatin binding confers locus specificity for productive chromatin remodeling and nucleosome eviction by BAF. Our study suggests a paradigm for how functional synergy between dynamically acting chromatin factors regulates locus-specific nucleosome organization and chromatin accessibility.

Joint epigenome profiling reveals cell-type-specific gene regulatory programmes in human cortical organoids

Gene expression is regulated by multiple epigenetic mechanisms, which are coordinated in development and disease. However, current multiomics methods are frequently limited to one or two modalities at a time, making it challenging to obtain a comprehensive gene regulatory signature. Here, we describe a method-3D genome, RNA, accessibility and methylation sequencing (3DRAM-seq)-that simultaneously interrogates spatial genome organization, chromatin accessibility and DNA methylation genome-wide and at high resolution. We combine 3DRAM-seq with immunoFACS and RNA sequencing in cortical organoids to map the cell-type-specific regulatory landscape of human neural development across multiple epigenetic layers. Finally, we apply a massively parallel reporter assay to profile cell-type-specific enhancer activity in organoids and to functionally assess the role of key transcription factors for human enhancer activation and function. More broadly, 3DRAM-seq can be used to profile the multimodal epigenetic landscape in rare cell types and different tissues.

Abnormal expression of PRKAG2-AS results in dysfunction of cardiomyocytes through regulating PRKAG2 transcription by interacting with PPARG

The role of PRKAG2 in the maintenance of heart function is well established, but little is known about how PRKAG2 is regulated in cardiomyocytes. In this study, we investigated the role of the lncRNA PRKAG2-AS, which is present at the PRKAG2 promoter, in the regulation of PRKAG2 expression. PRKAG2-AS expression was predominantly nuclear, as determined by RNA nucleoplasmic separation and fluorescence in situ hybridization. Knockdown of PRKAG2-AS in the nucleus, but not the cytoplasm, significantly decreased the expression of PRKAG2b and PRKAG2d. Interestingly, we found that PRKAG2-AS and its target genes, PRKAG2b and PRKAG2d, were reduced in the hearts of patients with ischemic cardiomyopathy, suggesting a potential role for PRKAG2-AS in myocardial ischemia. Indeed, knockdown of PRKAG2-AS in the nucleus resulted in apoptosis of cardiomyocytes. We further elucidated the mechanism by which PRKAG2-AS regulates PRKAG2 transcription by identifying 58 PRKAG2-AS interacting proteins. Among them, PPARG was selected for further investigation based on its correlation and potential interaction with PRKAG2-AS in regulating transcription. Overexpression of PPARG, or its activation with rosiglitazone, led to a significant increase in the expression of PRKAG2b and PRKAG2d in cardiomyocytes, which could be attenuated by PRKAG2-AS knockdown. This finding suggests that PRKAG2-AS mediates, at least partially, the protective effects of rosiglitazone on hypoxia-induced apoptosis. However, given the risk of rosiglitazone in heart failure, we also examined the involvement of PRKAG2-AS in this condition and found that PRKAG2-AS, as well as PRKAG2b and PRKAG2d, was elevated in hearts with dilated cardiomyopathy (DCM) and that overexpression of PRKAG2-AS led to a significant increase in PRKAG2b and PRKAG2d expression, indicating that up-regulation of PRKAG2-AS may contribute to the mechanism of heart failure by promoting transcription of PRKAG2. Consequently, proper expression of PRKAG2-AS is essential for maintaining cardiomyocyte function, and aberrant PRKAG2-AS expression induced by hypoxia or other stimuli may cause cardiac dysfunction.

Gene regulatory network inference in the era of single-cell multi-omics

The interplay between chromatin, transcription factors and genes generates complex regulatory circuits that can be represented as gene regulatory networks (GRNs). The study of GRNs is useful to understand how cellular identity is established, maintained and disrupted in disease. GRNs can be inferred from experimental data - historically, bulk omics data - and/or from the literature. The advent of single-cell multi-omics technologies has led to the development of novel computational methods that leverage genomic, transcriptomic and chromatin accessibility information to infer GRNs at an unprecedented resolution. Here, we review the key principles of inferring GRNs that encompass transcription factor-gene interactions from transcriptomics and chromatin accessibility data. We focus on the comparison and classification of methods that use single-cell multimodal data. We highlight challenges in GRN inference, in particular with respect to benchmarking, and potential further developments using additional data modalities.

Paired yeast one-hybrid assays to detect DNA-binding cooperativity and antagonism across transcription factors

Cooperativity and antagonism between transcription factors (TFs) can drastically modify their binding to regulatory DNA elements. While mapping these relationships between TFs is important for understanding their context-specific functions, existing approaches either rely on DNA binding motif predictions, interrogate one TF at a time, or study individual TFs in parallel. Here, we introduce paired yeast one-hybrid (pY1H) assays to detect cooperativity and antagonism across hundreds of TF-pairs at DNA regions of interest. We provide evidence that a wide variety of TFs are subject to modulation by other TFs in a DNA region-specific manner. We also demonstrate that TF-TF relationships are often affected by alternative isoform usage and identify cooperativity and antagonism between human TFs and viral proteins from human papillomaviruses, Epstein-Barr virus, and other viruses. Altogether, pY1H assays provide a broadly applicable framework to study how different functional relationships affect protein occupancy at regulatory DNA regions.

Genome-wide single-molecule analysis of long-read DNA methylation reveals heterogeneous patterns at heterochromatin that reflect nucleosome organisation

High-throughput sequencing technology is central to our current understanding of the human methylome. The vast majority of studies use chemical conversion to analyse bulk-level patterns of DNA methylation across the genome from a population of cells. While this technology has been used to probe single-molecule methylation patterns, such analyses are limited to short reads of a few hundred basepairs. DNA methylation can also be directly detected using Nanopore sequencing which can generate reads measuring megabases in length. However, thus far these analyses have largely focused on bulk-level assessment of DNA methylation. Here, we analyse DNA methylation in single Nanopore reads from human lymphoblastoid cells, to show that bulk-level metrics underestimate large-scale heterogeneity in the methylome. We use the correlation in methylation state between neighbouring sites to quantify single-molecule heterogeneity and find that heterogeneity varies significantly across the human genome, with some regions having heterogeneous methylation patterns at the single-molecule level and others possessing more homogeneous methylation patterns. By comparing the genomic distribution of the correlation to epigenomic annotations, we find that the greatest heterogeneity in single-molecule patterns is observed within heterochromatic partially methylated domains (PMDs). In contrast, reads originating from euchromatic regions and gene bodies have more ordered DNA methylation patterns. By analysing the patterns of single molecules in more detail, we show the existence of a nucleosome-scale periodicity in DNA methylation that accounts for some of the heterogeneity we uncover in long single-molecule DNA methylation patterns. We find that this periodic structure is partially masked in bulk data and correlates with DNA accessibility as measured by nanoNOMe-seq, suggesting that it could be generated by nucleosomes. Our findings demonstrate the power of single-molecule analysis of long-read data to understand the structure of the human methylome.

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.

Research Note: Identification of core promoter region of the polyunsaturated fatty acid synthesis-related gene family in chicken

Chicken is considered an ideal model species to study the synthesis of polyunsaturated fatty acids (PUFAs) due to its appropriate proportions of fatty acids and abundant content of PUFAs, suitable for human consumption. However, the molecular mechanisms regulating poultry PUFA synthesis remain unclear. Here, we systematically explored the transcriptional regulation activity of the gene family related to PUFA synthesis in chicken by carrying out the Dual-Luciferase Reporter Assay. We identified the core promoter regions of members of the chicken PUFA synthesis-related gene family, including ELOVL1, ELOVL2, ELOVL3, ELOVL4, ELOVL5, ELOVL6, ELOVL7, FADS1, FADS2, FADS6, SCD, and SCD5. Additionally, changes in relative fluorescence values of different truncated segments in the upstream regulatory region of these genes indicate the existence of regulatory regions. Furthermore, we predicted the transcription factors that bind to the identified core promoter regions of multiple genes, including Sp1, NF-1, C/EBPalpha, etc. These findings provide a basis for the molecular mechanisms regulating poultry PUFA synthesis and offer new scientific insight into the potential improvement of poultry meat quality in the future.

Beyond assembly: the increasing flexibility of single-molecule sequencing technology

The maturation of high-throughput short-read sequencing technology over the past two decades has shaped the way genomes are studied. Recently, single-molecule, long-read sequencing has emerged as an essential tool in deciphering genome structure and function, including filling gaps in the human reference genome, measuring the epigenome and characterizing splicing variants in the transcriptome. With recent technological developments, these single-molecule technologies have moved beyond genome assembly and are being used in a variety of ways, including to selectively sequence specific loci with long reads, measure chromatin state and protein-DNA binding in order to investigate the dynamics of gene regulation, and rapidly determine copy number variation. These increasingly flexible uses of single-molecule technologies highlight a young and fast-moving part of the field that is leading to a more accessible era of nucleic acid sequencing.

Histone modifications regulate pioneer transcription factor cooperativity

Pioneer transcription factors have the ability to access DNA in compacted chromatin1. Multiple transcription factors can bind together to a regulatory element in a cooperative way, and cooperation between the pioneer transcription factors OCT4 (also known as POU5F1) and SOX2 is important for pluripotency and reprogramming2-4. However, the molecular mechanisms by which pioneer transcription factors function and cooperate on chromatin remain unclear. Here we present cryo-electron microscopy structures of human OCT4 bound to a nucleosome containing human LIN28B or nMATN1 DNA sequences, both of which bear multiple binding sites for OCT4. Our structural and biochemistry data reveal that binding of OCT4 induces changes to the nucleosome structure, repositions the nucleosomal DNA and facilitates cooperative binding of additional OCT4 and of SOX2 to their internal binding sites. The flexible activation domain of OCT4 contacts the N-terminal tail of histone H4, altering its conformation and thus promoting chromatin decompaction. Moreover, the DNA-binding domain of OCT4 engages with the N-terminal tail of histone H3, and post-translational modifications at H3K27 modulate DNA positioning and affect transcription factor cooperativity. Thus, our findings suggest that the epigenetic landscape could regulate OCT4 activity to ensure proper cell programming.

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.

Integrative approaches to study enhancer-promoter communication

The spatiotemporal control of gene expression in complex multicellular organisms relies on noncoding regulatory sequences such as enhancers, which activate transcription of target genes often over large genomic distances. Despite the advances in the identification and characterization of enhancers, the principles and mechanisms by which enhancers select and control their target genes remain largely unknown. Here, we review recent interdisciplinary and quantitative approaches based on emerging techniques that aim to address open questions in the field, notably how regulatory information is encoded in the DNA sequence, how this information is transferred from enhancers to promoters, and how these processes are regulated in time.

Epigenetic Control of Cell Potency and Fate Determination during Mammalian Gastrulation

Pluripotent embryonic stem cells have a unique and characteristic epigenetic profile, which is critical for differentiation to all embryonic germ lineages. When stem cells exit the pluripotent state and commit to lineage-specific identities during the process of gastrulation in early embryogenesis, extensive epigenetic remodelling mediates both the switch in cellular programme and the loss of potential to adopt alternative lineage programmes. However, it remains to be understood how the stem cell epigenetic profile encodes pluripotency, or how dynamic epigenetic regulation helps to direct cell fate specification. Recent advances in stem cell culture techniques, cellular reprogramming, and single-cell technologies that can quantitatively profile epigenetic marks have led to significant insights into these questions, which are important for understanding both embryonic development and cell fate engineering. This review provides an overview of key concepts and highlights exciting new advances in the field.

NOMe-HiC: joint profiling of genetic variant, DNA methylation, chromatin accessibility, and 3D genome in the same DNA molecule

Cis-regulatory elements are coordinated to regulate the expression of their targeted genes. However, the joint measurement of cis-regulatory elements' activities and their interactions in spatial proximity is limited by the current sequencing approaches. We describe a method, NOMe-HiC, which simultaneously captures single-nucleotide polymorphisms, DNA methylation, chromatin accessibility (GpC methyltransferase footprints), and chromosome conformation changes from the same DNA molecule, together with the transcriptome, in a single assay. NOMe-HiC shows high concordance with state-of-the-art mono-omic assays across different molecular measurements and reveals coordinated chromatin accessibility at distal genomic segments in spatial proximity and novel types of long-range allele-specific chromatin accessibility.

Histone modifications regulate pioneer transcription factor binding and cooperativity

Pioneer transcription factors have the ability to access DNA in compacted chromatin. Multiple transcription factors can bind together to a regulatory element in a cooperative way and cooperation between pioneer transcription factors Oct4 and Sox2 is important for pluripotency and reprogramming. However, the molecular mechanisms by which pioneer transcription factors function and cooperate remain unclear. Here we present cryo-EM structures of human Oct4 bound to a nucleosome containing human Lin28B and nMatn1 DNA sequences, which bear multiple binding sites for Oct4. Our structural and biochemistry data reveal that Oct4 binding induces changes to the nucleosome structure, repositions the nucleosomal DNA and facilitates cooperative binding of additional Oct4 and of Sox2 to their internal binding sites. The flexible activation domain of Oct4 contacts the histone H4 N-terminal tail, altering its conformation and thus promoting chromatin decompaction. Moreover, the DNA binding domain of Oct4 engages with histone H3 N-terminal tail, and posttranslational modifications at H3K27 modulate DNA positioning and affect transcription factor cooperativity. Thus, our data show that the epigenetic landscape can regulate Oct4 activity to ensure proper cell reprogramming.

Single-molecule footprinting identifies context-dependent regulation of enhancers by DNA methylation

Enhancers are cis-regulatory elements that control the establishment of cell identities during development. In mammals, enhancer activation is tightly coupled with DNA demethylation. However, whether this epigenetic remodeling is necessary for enhancer activation is unknown. Here, we adapted single-molecule footprinting to measure chromatin accessibility and transcription factor binding as a function of the presence of methylation on the same DNA molecules. We leveraged natural epigenetic heterogeneity at active enhancers to test the impact of DNA methylation on their chromatin accessibility in multiple cell lineages. Although reduction of DNA methylation appears dispensable for the activity of most enhancers, we identify a class of cell-type-specific enhancers where DNA methylation antagonizes the binding of transcription factors. Genetic perturbations reveal that chromatin accessibility and transcription factor binding require active demethylation at these loci. Thus, in addition to safeguarding the genome from spurious activation, DNA methylation directly controls transcription factor occupancy at active enhancers.

Deciphering the multi-scale, quantitative cis-regulatory code

Uncovering the cis-regulatory code that governs when and how much each gene is transcribed in a given genome and cellular state remains a central goal of biology. Here, we discuss major layers of regulation that influence how transcriptional outputs are encoded by DNA sequence and cellular context. We first discuss how transcription factors bind specific DNA sequences in a dosage-dependent and cooperative manner and then proceed to the cofactors that facilitate transcription factor function and mediate the activity of modular cis-regulatory elements such as enhancers, silencers, and promoters. We then consider the complex and poorly understood interplay of these diverse elements within regulatory landscapes and its relationships with chromatin states and nuclear organization. We propose that a mechanistically informed, quantitative model of transcriptional regulation that integrates these multiple regulatory layers will be the key to ultimately cracking the cis-regulatory code.

RNA Polymerase II, the BAF remodeler and transcription factors synergize to evict nucleosomes

Chromatin accessibility is a hallmark of active transcription and requires ATP-dependent nucleosome remodeling by Brahma-Associated Factor (BAF). However, the mechanistic link between transcription, nucleosome remodeling, and chromatin accessibility is unclear. Here, we used a chemical-genetic approach to dissect the interplay between RNA Polymerase II (RNAPII), BAF, and DNA-sequence-specific transcription factors (TFs) in mouse embryonic stem cells. By time-resolved chromatin profiling with acute transcription block at distinct stages, we show that RNAPII promoter-proximal pausing stabilizes BAF chromatin occupancy and enhances nucleosome eviction by BAF. We find that RNAPII and BAF probe both transcriptionally active and Polycomb-repressed genomic regions and provide evidence that TFs capture transient site exposure due to nucleosome unwrapping by BAF to confer locus specificity for persistent chromatin remodeling. Our study reveals the mechanistic basis of cell-type-specific chromatin accessibility. We propose a new paradigm for how functional synergy between dynamically acting chromatin factors regulates nucleosome organization.

Patient-specific identification of genome-wide DNA-methylation differences between intracranial and extracranial melanoma metastases

Melanomas frequently metastasize to distant organs and especially intracranial metastases still represent a major clinical challenge. Epigenetic reprogramming of intracranial metastases is thought to be involved in therapy failure, but so far only little is known about patient-specific DNA-methylation differences between intra- and extracranial melanoma metastases. Hierarchical clustering of the methylomes of 24 patient-matched intra- and extracranial melanoma metastases pairs revealed that intra- and extracranial metastases of individual patients were more similar to each other than to metastases in the same tissue from other patients. Therefore, a personalized analysis of each metastases pair was done by a Hidden Markov Model to classify methylation levels of individual CpGs as decreased, unchanged or increased in the intra- compared to the extracranial metastasis. The predicted DNA-methylation alterations were highly patient-specific differing in the number and methylation states of altered CpGs. Nevertheless, four important general observations were made: (i) intracranial metastases of most patients mainly showed a reduction of DNA-methylation, (ii) cytokine signaling was most frequently affected by differential methylation in individual metastases pairs, but also MAPK, PI3K/Akt and ECM signaling were often altered, (iii) frequently affected genes were mainly involved in signaling, growth, adhesion or apoptosis, and (iv) an enrichment of functional terms related to channel and transporter activities supports previous findings for a brain-like phenotype. In addition, the derived set of 17 signaling pathway genes that distinguished intra- from extracranial metastases in more than 50% of patients included well-known oncogenes (e.g. PRKCA, DUSP6, BMP4) and several other genes known from neuronal disorders (e.g. EIF4B, SGK1, CACNG8). Moreover, associations of gene body methylation alterations with corresponding gene expression changes revealed that especially the three signaling pathway genes JAK3, MECOM, and TNXB differ strongly in their expression between patient-matched intra- and extracranial metastases. Our analysis contributes to an in-depth characterization of DNA-methylation differences between patient-matched intra- and extracranial melanoma metastases and may provide a basis for future experimental studies to identify targets for new therapeutic approaches.

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.

Chromatin accessibility: methods, mechanisms, and biological insights

Access to DNA is a prerequisite to the execution of essential cellular processes that include transcription, replication, chromosomal segregation, and DNA repair. How the proteins that regulate these processes function in the context of chromatin and its dynamic architectures is an intensive field of study. Over the past decade, genome-wide assays and new imaging approaches have enabled a greater understanding of how access to the genome is regulated by nucleosomes and associated proteins. Additional mechanisms that may control DNA accessibility in vivo include chromatin compaction and phase separation - processes that are beginning to be understood. Here, we review the ongoing development of accessibility measurements, we summarize the different molecular and structural mechanisms that shape the accessibility landscape, and we detail the many important biological functions that are linked to chromatin accessibility.

Long-range phasing of dynamic, tissue-specific and allele-specific regulatory elements

Epigenomic maps identify gene regulatory elements by their chromatin state. However, prevailing short-read sequencing methods cannot effectively distinguish alleles, evaluate the interdependence of elements in a locus or capture single-molecule dynamics. Here, we apply targeted nanopore sequencing to profile chromatin accessibility and DNA methylation on contiguous ~100-kb DNA molecules that span loci relevant to development, immunity and imprinting. We detect promoters, enhancers, insulators and transcription factor footprints on single molecules based on exogenous GpC methylation. We infer relationships among dynamic elements within immune loci, and order successive remodeling events during T cell stimulation. Finally, we phase primary sequence and regulatory elements across the H19/IGF2 locus, uncovering primate-specific features. These include a segmental duplication that stabilizes the imprinting control region and a noncanonical enhancer that drives biallelic IGF2 expression in specific contexts. Our study advances emerging strategies for phasing gene regulatory landscapes and reveals a mechanism that overrides IGF2 imprinting in human cells.

Structures and consequences of pioneer factor binding to nucleosomes

Pioneer transcription factors are able to bind a partially exposed motif on the surface of a nucleosome, enabling the proteins to target sites in silent regions of chromatin that have been compacted by linker histone. The targeting of nucleosomal DNA by pioneer factors has been observed in vitro and in vivo, where binding can promote local nucleosome exposure that allows other transcription factors, nucleosome remodelers, and histone modifiers to engage the chromatin and elicit gene activation or further repression. Pioneer factors thereby establish new gene expression programs during cell fate changes that occur during embryonic development, regeneration, and cancer. Here, we review recent biophysical studies that reveal the structural features and strategies used by pioneer factors to accomplish nucleosome binding and the consequential changes to nucleosomes that can lead to DNA accessibility.

Elucidating the biology of transcription factor-DNA interaction for accurate identification of cis-regulatory elements

Transcription factors (TFs) play a critical role in determining cell fate decisions by integrating developmental and environmental signals through binding to specific cis-regulatory modules and regulating spatio-temporal specificity of gene expression patterns. Precise identification of functional TF binding sites in time and space not only will revolutionize our understanding of regulatory networks governing cell fate decisions but is also instrumental to uncover how genetic variations cause morphological diversity or disease. In this review, we discuss recent advances in mapping TF binding sites and characterizing the various parameters underlying the complexity of binding site recognition by TFs.

A landscape of gene regulation in the parasitic amoebozoa Entamoeba spp

Entamoeba are amoeboid extracellular parasites that represent an important group of organisms for which the regulatory networks must be examined to better understand how genes and functional processes are interrelated. In this work, we inferred the gene regulatory networks (GRNs) in four Entamoeba species, E. histolytica, E. dispar, E. nuttalli, and E. invadens, and the GRN topological properties and the corresponding biological functions were evaluated. From these analyses, we determined that transcription factors (TFs) of E. histolytica, E. dispar, and E. nuttalli are associated mainly with the LIM family, while the TFs in E. invadens are associated with the RRM_1 family. In addition, we identified that EHI_044890 regulates 121 genes in E. histolytica, EDI_297980 regulates 284 genes in E. dispar, ENU1_120230 regulates 195 genes in E. nuttalli, and EIN_249270 regulates 257 genes in E. invadens. Finally, we identified that three types of processes, Macromolecule metabolic process, Cellular macromolecule metabolic process, and Cellular nitrogen compound metabolic process, are the main biological processes for each network. The results described in this work can be used as a basis for the study of gene regulation in these organisms.

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.

TF-COMB - discovering grammar of transcription factor binding sites

Cooperativity between transcription factors is important to regulate target gene expression. In particular, the binding grammar of TFs in relation to each other, as well as in the context of other genomic elements, is crucial for TF functionality. However, tools to easily uncover co-occurrence between DNA-binding proteins, and investigate the regulatory modules of TFs, are limited. Here we present TF-COMB (Transcription Factor Co-Occurrence using Market Basket analysis) - a tool to investigate co-occurring TFs and binding grammar within regulatory regions. We found that TF-COMB can accurately identify known co-occurring TFs from ChIP-seq data, as well as uncover preferential localization to other genomic elements. With the use of ATAC-seq footprinting and TF motif locations, we found that TFs exhibit both preferred orientation and distance in relation to each other, and that these are biologically significant. Finally, we extended the analysis to not only investigate individual TF pairs, but also TF pairs in the context of networks, which enabled the investigation of TF complexes and TF hubs. In conclusion, TF-COMB is a flexible tool to investigate various aspects of TF binding grammar.

Mechanisms that regulate the activities of TET proteins

The ten-eleven translocation (TET) family of dioxygenases consists of three members, TET1, TET2, and TET3. All three TET enzymes have Fe⁺² and α-ketoglutarate (α-KG)-dependent dioxygenase activities, catalyzing the 1st step of DNA demethylation by converting 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC), and further oxidize 5hmC to 5-formylcytosine (5fC) and 5-carboxylcytosine (5caC). Gene knockout studies demonstrated that all three TET proteins are involved in the regulation of fetal organ generation during embryonic development and normal tissue generation postnatally. TET proteins play such roles by regulating the expression of key differentiation and fate-determining genes via (1) enzymatic activity-dependent DNA methylation of the promoters and enhancers of target genes; and (2) enzymatic activity-independent regulation of histone modification. Interacting partner proteins and post-translational regulatory mechanisms regulate the activities of TET proteins. Mutations and dysregulation of TET proteins are involved in the pathogenesis of human diseases, specifically cancers. Here, we summarize the research on the interaction partners and post-translational modifications of TET proteins. We also discuss the molecular mechanisms by which these partner proteins and modifications regulate TET functioning and target gene expression. Such information will help in the design of medications useful for targeted therapy of TET-mutant-related diseases.

Hormone-controlled cooperative binding of transcription factors drives synergistic induction of fasting-regulated genes

During fasting, hepatocytes produce glucose in response to hormonal signals. Glucagon and glucocorticoids are principal fasting hormones that cooperate in regulating glucose production via gluconeogenesis. However, how these hormone signals are integrated and interpreted to a biological output is unknown. Here, we use genome-wide profiling of gene expression, enhancer dynamics and transcription factor (TF) binding in primary mouse hepatocytes to uncover the mode of cooperation between glucagon and glucocorticoids. We found that compared to a single treatment with each hormone, a dual treatment directs hepatocytes to a pro-gluconeogenic gene program by synergistically inducing gluconeogenic genes. The cooperative mechanism driving synergistic gene expression is based on ‘assisted loading’ whereby a glucagon-activated TF (cAMP responsive element binding protein; CREB) leads to enhancer activation which facilitates binding of the glucocorticoid receptor (GR) upon glucocorticoid stimulation. Glucagon does not only activate single enhancers but also activates enhancer clusters, thereby assisting the loading of GR also across enhancer units within the cluster. In summary, we show that cells integrate extracellular signals by an enhancer-specific mechanism: one hormone-activated TF activates enhancers, thereby assisting the loading of a TF stimulated by a second hormone, leading to synergistic gene induction and a tailored transcriptional response to fasting.

Transcription factor dynamics in plants: Insights and technologies for in vivo imaging

Biochemical and genetic approaches have been extensively used to study transcription factor (TF) functions, but their dynamic behaviors and the complex ways in which they regulate transcription in plant cells remain unexplored, particularly behaviors such as translocation and binding to DNA. Recent developments in labeling and imaging techniques provide the necessary sensitivity and resolution to study these behaviors in living cells. In this review, we present an up-to-date portrait of the dynamics and regulation of TFs under physiologically relevant conditions and then summarize recent advances in fluorescent labeling strategies and imaging techniques. We then discuss future prospects and challenges associated with the application of these techniques to examine TFs' intricate dance in living plants.

Different responses to risperidone treatment in Schizophrenia: a multicenter genome-wide association and whole exome sequencing joint study

Risperidone is routinely used in the clinical management of schizophrenia, but the treatment response is highly variable among different patients. The genetic underpinnings of the treatment response are not well understood. We performed a pharmacogenomic study of the treatment response to risperidone in patients with schizophrenia by using a SNP microarray -based genome-wide association study (GWAS) and whole exome sequencing (WES)-based GWAS. DNA samples were collected from 189 patients for the GWAS and from 222 patients for the WES after quality control in multiple centers of China. Antipsychotic response phenotypes of patients who received eight weeks of risperidone treatment were quantified with percentage change on the Positive and Negative Syndrome Scale (PANSS). The GWAS revealed a significant association between several SNPs and treatment response, such as three GRM7 SNPs (rs141134664, rs57521140, and rs73809055). Gene-based analysis in WES revealed 13 genes that were associated with antipsychotic response, such as GPR12 and MAP2K3. We did not identify shared loci or genes between GWAS and WES, but association signals tended to cluster into the GPCR gene family and GPCR signaling pathway, which may play an important role in the treatment response etiology. This study may provide a research paradigm for pharmacogenomic research, and these data provide a promising illustration of our potential to identify genetic variants underlying antipsychotic responses and may ultimately facilitate precision medicine in schizophrenia.

HNF4A modulates glucocorticoid action in the liver

The glucocorticoid receptor (GR) is a nuclear receptor critical to the regulation of energy metabolism and inflammation. The actions of GR are dependent on cell type and context. Here, we demonstrate the role of liver lineage-determining factor hepatocyte nuclear factor 4A (HNF4A) in defining liver specificity of GR action. In mouse liver, the HNF4A motif lies adjacent to the glucocorticoid response element (GRE) at GR binding sites within regions of open chromatin. In the absence of HNF4A, the liver GR cistrome is remodeled, with loss and gain of GR recruitment evident. Loss of chromatin accessibility at HNF4A-marked sites associates with loss of GR binding at weak GRE motifs. GR binding and chromatin accessibility are gained at sites characterized by strong GRE motifs, which show GR recruitment in non-liver tissues. The functional importance of these HNF4A-regulated GR sites is indicated by an altered transcriptional response to glucocorticoid treatment in the Hnf4a-null liver.

DiMeLo-seq: a long-read, single-molecule method for mapping protein-DNA interactions genome wide

Studies of genome regulation routinely use high-throughput DNA sequencing approaches to determine where specific proteins interact with DNA, and they rely on DNA amplification and short-read sequencing, limiting their quantitative application in complex genomic regions. To address these limitations, we developed Directed Methylation with Long-read sequencing (DiMeLo-seq), which uses antibody-tethered enzymes to methylate DNA near a target protein’s binding sites in situ. These exogenous methylation marks are then detected simultaneously with endogenous CpG methylation on unamplified DNA using long-read, single-molecule sequencing technologies. We optimized and benchmarked DiMeLo-seq by mapping chromatin-binding proteins and histone modifications across the human genome. Furthermore, we identified where centromere protein A (CENP-A) localizes within highly repetitive regions that we re unmappable with short sequencing reads, and we estimated the density of CENP-A molecules along single chromatin fibers. DiMeLo-seq is a versatile method that provides multimodal, genome-wide information for investigating protein-DNA interactions.

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.

Systematic analysis of naturally occurring insertions and deletions that alter transcription factor spacing identifies tolerant and sensitive transcription factor pairs

Regulation of gene expression requires the combinatorial binding of sequence-specific transcription factors (TFs) at promoters and enhancers. Prior studies showed that alterations in the spacing between TF binding sites can influence promoter and enhancer activity. However, the relative importance of TF spacing alterations resulting from naturally occurring insertions and deletions (InDels) has not been systematically analyzed. To address this question, we first characterized the genome-wide spacing relationships of 73 TFs in human K562 cells as determined by ChIP-seq (chromatin immunoprecipitation sequencing). We found a dominant pattern of a relaxed range of spacing between collaborative factors, including 45 TFs exclusively exhibiting relaxed spacing with their binding partners. Next, we exploited millions of InDels provided by genetically diverse mouse strains and human individuals to investigate the effects of altered spacing on TF binding and local histone acetylation. These analyses suggested that spacing alterations resulting from naturally occurring InDels are generally tolerated in comparison to genetic variants directly affecting TF binding sites. To experimentally validate this prediction, we introduced synthetic spacing alterations between PU.1 and C/EBPβ binding sites at six endogenous genomic loci in a macrophage cell line. Remarkably, collaborative binding of PU.1 and C/EBPβ at these locations tolerated changes in spacing ranging from 5 bp increase to >30 bp decrease. Collectively, these findings have implications for understanding mechanisms underlying enhancer selection and for the interpretation of non-coding genetic variation.

Genome-wide quantification of transcription factor binding at single-DNA-molecule resolution using methyl-transferase footprinting

Precise control of gene expression requires the coordinated action of multiple factors at cis-regulatory elements. We recently developed single-molecule footprinting to simultaneously resolve the occupancy of multiple proteins including transcription factors, RNA polymerase II and nucleosomes on single DNA molecules genome-wide. The technique combines the use of cytosine methyltransferases to footprint the genome with bisulfite sequencing to resolve transcription factor binding patterns at cis-regulatory elements. DNA footprinting is performed by incubating permeabilized nuclei with recombinant methyltransferases. Upon DNA extraction, whole-genome or targeted bisulfite libraries are prepared and loaded on Illumina sequencers. The protocol can be completed in 4-5 d in any laboratory with access to high-throughput sequencing. Analysis can be performed in 2 d using a dedicated R package and requires access to a high-performance computing system. Our method can be used to analyze how transcription factors cooperate and antagonize to regulate transcription.

Regulation of transcription reactivation dynamics exiting mitosis

Proliferating cells experience a global reduction of transcription during mitosis, yet their cell identity is maintained and regulatory information is propagated from mother to daughter cells. Mitotic bookmarking by transcription factors has been proposed as a potential mechanism to ensure the reactivation of transcription at the proper set of genes exiting mitosis. Recently, mitotic transcription and waves of transcription reactivation have been observed in synchronized populations of human hepatoma cells. However, the study did not consider that mitotic-arrested cell populations progressively desynchronize leading to measurements of gene expression on a mixture of cells at different internal cell-cycle times. Moreover, it is not well understood yet what is the precise role of mitotic bookmarking on mitotic transcription as well as on the transcription reactivation waves. Ultimately, the core gene regulatory network driving the precise transcription reactivation dynamics remains to be identified. To address these questions, we developed a mathematical model to correct for the progressive desynchronization of cells and estimate gene expression dynamics with respect to a cell-cycle pseudotime. Furthermore, we used a multiple linear regression model to infer transcription factor activity dynamics. Our analysis allows us to characterize waves of transcription factor activities exiting mitosis and predict a core gene regulatory network responsible of the transcription reactivation dynamics. Moreover, we identified more than 60 transcription factors that are highly active during mitosis and represent new candidates of mitotic bookmarking factors which could be relevant therapeutic targets to control cell proliferation.

Specific gene expression patterns confer particular identities to cells. During proliferation, cells undergo mitosis when chromosomes are formed and segregated into two new cells leading to a global downregulation of gene expression. Yet, cell identity is propagated from mother to daughter cells by the reactivation of gene expression at the appropriate set of genes once mitosis is completed. Mitotic bookmarking has been proposed as a mechanism to regulate this process. Indeed certain regulatory factors tag genes during mitosis to promote gene reactivation in the next cycle. Here we analyze gene expression over time measured on synchronized cell populations by using a new generation sequencing technique. To do so, we proposed a mathematical model to obtain the exact gene expression dynamics with respect to the cell-cycle progression and identified waves of genes reactivation during mitosis and the transition to the next cycle. Also, we developed a computational method that allowed us to predict key regulatory factors that drive this process and predict new candidates that could be involved in mitotic bookmarking. These regulatory factors could be relevant therapeutic targets to control cell proliferation.