Genome-wide Single-Molecule Footprinting Reveals High RNA Polymerase II Turnover at Paused Promoters

Multiple Forms and Functions of Premature Termination by RNA Polymerase II

Eukaryotic genomes are widely transcribed by RNA polymerase II (pol II) both within genes and in intergenic regions. POL II elongation complexes comprising the polymerase, the DNA template and nascent RNA transcript must be extremely processive in order to transcribe the longest genes which are over 1 megabase long and take many hours to traverse. Dedicated termination mechanisms are required to disrupt these highly stable complexes. Transcription termination occurs not only at the 3' ends of genes once a full length transcript has been made, but also within genes and in promiscuously transcribed intergenic regions. Termination at these latter positions is termed "premature" because it is not triggered in response to a specific signal that marks the 3' end of a gene, like a polyA site. One purpose of premature termination is to remove polymerases from intergenic regions where they are "not wanted" because they may interfere with transcription of overlapping genes or the progress of replication forks. Premature termination has recently been appreciated to occur at surprisingly high rates within genes where it is speculated to serve regulatory or quality control functions. In this review I summarize current understanding of the different mechanisms of premature termination and its potential functions.

Defining gene ends: RNA polymerase II CTD threonine 4 phosphorylation marks transcription termination regions genome-wide

Defining the beginning of a eukaryotic protein-coding gene is relatively simple. It corresponds to the first ribonucleotide incorporated by RNA polymerase II (Pol II) into the nascent RNA molecule. This nucleotide is protected by capping and maintained in the mature messenger RNA (mRNA). However, in higher eukaryotes, the end of mRNA is separated from the sites of transcription termination by hundreds to thousands of base pairs. Currently used genomic annotations only take account of the end of the mature transcript - the sites where pre-mRNA cleavage occurs, while the regions in which transcription terminates are unannotated. Here, we describe the evidence for a marker of transcription termination, which could be widely applicable in genomic studies. Pol II termination regions can be determined genome-wide by detecting Pol II phosphorylated on threonine 4 of its C-terminal domain (Pol II CTD-T4ph). Pol II in this state pauses before leaving the DNA template. Up to date this potent mark has been underused because the evidence for its place and role in termination is scattered across multiple publications. We summarize the observations regarding Pol II CTD-T4ph in termination regions and present bioinformatic analyses that further support Pol II CTD-T4ph as a global termination mark in animals.

The PNUTS phosphatase complex controls transcription pause release

Gene expression is regulated by controlling distinct steps of the transcriptional cycle, including initiation, pausing, elongation, and termination. Kinases phosphorylate RNA polymerase II (RNA Pol II) and associated factors to control transitions between these steps and to act as central gene regulatory nodes. Similarly, phosphatases that dephosphorylate these components are emerging as important regulators of transcription, although their roles remain less well understood. Here, we discover that the mouse PNUTS-PP1 phosphatase complex plays an essential role in controlling transcription pause release in addition to its previously described function in transcription termination. Transcription pause release by the PNUTS complex is essential for almost all RNA Pol II-dependent gene transcription, relies on its PP1 phosphatase subunit, and controls the phosphorylation of factors required for pause release and elongation. Together, these observations reveal an essential new role for a phosphatase complex in transcription pause release and show that the PNUTS complex is essential for RNA Pol II-dependent transcription.

Single-molecule states link transcription factor binding to gene expression

The binding of multiple transcription factors (TFs) to genomic enhancers drives gene expression in mammalian cells1. However, the molecular details that link enhancer sequence to TF binding, promoter state and transcription levels remain unclear. Here we applied single-molecule footprinting2,3 to measure the simultaneous occupancy of TFs, nucleosomes and other regulatory proteins on engineered enhancer-promoter constructs with variable numbers of TF binding sites for both a synthetic TF and an endogenous TF involved in the type I interferon response. Although TF binding events on nucleosome-free DNA are independent, activation domains recruit cofactors that destabilize nucleosomes, driving observed TF binding cooperativity. Average TF occupancy linearly determines promoter activity, and we decompose TF strength into separable binding and activation terms. Finally, we develop thermodynamic and kinetic models that quantitatively predict both the enhancer binding microstates and gene expression dynamics. This work provides a template for the quantitative dissection of distinct contributors to gene expression, including TF activation domains, concentration, binding affinity, binding site configuration and recruitment of chromatin regulators.

Systematic screening of enhancer-blocking insulators in Drosophila identifies their DNA sequence determinants

Long-range transcriptional activation of gene promoters by abundant enhancers in animal genomes calls for mechanisms to limit inappropriate regulation. DNA elements called insulators serve this purpose by shielding promoters from an enhancer when interposed. Unlike promoters and enhancers, insulators have not been systematically characterized due to lacking high-throughput screening assays, and questions regarding how insulators are distributed and encoded in the genome remain. Here, we establish "insulator-seq" as a plasmid-based massively parallel reporter assay in Drosophila cultured cells to perform a systematic insulator screen of selected genomic loci. Screening developmental gene loci showed that not all insulator protein binding sites effectively block enhancer-promoter communication. Deep insulator mutagenesis identified sequences flexibly positioned around the CTCF insulator protein binding motif that are critical for functionality. The ability to screen millions of DNA sequences without positional effect has enabled functional mapping of insulators and provided further insights into the determinants of insulators.

Spatially resolved epigenome sequencing via Tn5 transposition and deterministic DNA barcoding in tissue

Spatial epigenetic mapping of tissues enables the study of gene regulation programs and cellular functions with the dependency on their local tissue environment. Here we outline a complete procedure for two spatial epigenomic profiling methods: spatially resolved genome-wide profiling of histone modifications using in situ cleavage under targets and tagmentation (CUT&Tag) chemistry (spatial-CUT&Tag) and transposase-accessible chromatin sequencing (spatial-ATAC-sequencing) for chromatin accessibility. Both assays utilize in-tissue Tn5 transposition to recognize genomic DNA loci followed by microfluidic deterministic barcoding to incorporate spatial address codes. Furthermore, these two methods do not necessitate prior knowledge of the transcription or epigenetic markers for a given tissue or cell type but permit genome-wide unbiased profiling pixel-by-pixel at the 10 μm pixel size level and single-base resolution. To support the widespread adaptation of these methods, details are provided in five general steps: (1) sample preparation; (2) Tn5 transposition in spatial-ATAC-sequencing or antibody-controlled pA-Tn5 tagmentation in CUT&Tag; (3) library preparation; (4) next-generation sequencing; and (5) data analysis using our customed pipelines available at: https://github.com/dyxmvp/Spatial_ATAC-seq and https://github.com/dyxmvp/spatial-CUT-Tag . The whole procedure can be completed on four samples in 2-3 days. Familiarity with basic molecular biology and bioinformatics skills with access to a high-performance computing environment are required. A rudimentary understanding of pathology and specimen sectioning, as well as deterministic barcoding in tissue-specific skills (e.g., design of a multiparameter barcode panel and creation of microfluidic devices), are also advantageous. In this protocol, we mainly focus on spatial profiling of tissue region-specific epigenetic landscapes in mouse embryos and mouse brains using spatial-ATAC-sequencing and spatial-CUT&Tag, but these methods can be used for other species with no need for species-specific probe design.

Nuclear sorting of short RNA polymerase II transcripts

Mammalian genomes produce an abundance of short RNA. This is, to a large extent, due to the genome-wide and spurious activity of RNA polymerase II (RNAPII). However, it is also because the vast majority of initiating RNAPII, regardless of the transcribed DNA unit, terminates within a ∼3-kb early "pausing zone." Given that the resultant RNAs constitute both functional and non-functional species, their proper sorting is critical. One way to think about such quality control (QC) is that transcripts, from their first emergence, are relentlessly targeted by decay factors, which may only be avoided by engaging protective processing pathways. In a molecular materialization of this concept, recent progress has found that both "destructive" and "productive" RNA effectors assemble at the 5' end of capped RNA, orchestrated by the essential arsenite resistance protein 2 (ARS2) protein. Based on this principle, we here discuss early QC mechanisms and how these might sort short RNAs to their final fates.

Genome access is transcription factor-specific and defined by nucleosome position

Mammalian gene expression is controlled by transcription factors (TFs) that engage sequence motifs in a chromatinized genome, where nucleosomes can restrict DNA access. Yet, how nucleosomes affect individual TFs remains unclear. Here, we measure the ability of over one hundred TF motifs to recruit TFs in a defined chromosomal locus in mouse embryonic stem cells. This identifies a set sufficient to enable the binding of TFs with diverse tissue specificities, functions, and DNA-binding domains. These chromatin-competent factors are further classified when challenged to engage motifs within a highly phased nucleosome. The pluripotency factors OCT4-SOX2 preferentially engage non-nucleosomal and entry-exit motifs, but not nucleosome-internal sites, a preference that also guides binding genome wide. By contrast, factors such as BANP, REST, or CTCF engage throughout, causing nucleosomal displacement. This supports that TFs vary widely in their sensitivity to nucleosomes and that genome access is TF specific and influenced by nucleosome position in the cell.

RNA polymerases reshape chromatin architecture and couple transcription on individual fibers

RNA polymerases must initiate and pause within a complex chromatin environment, surrounded by nucleosomes and other transcriptional machinery. This environment creates a spatial arrangement along individual chromatin fibers ripe for both competition and coordination, yet these relationships remain largely unknown owing to the inherent limitations of traditional structural and sequencing methodologies. To address this, we employed long-read chromatin fiber sequencing (Fiber-seq) in Drosophila to visualize RNA polymerase (Pol) within its native chromatin context with single-molecule precision along up to 30 kb fibers. We demonstrate that Fiber-seq enables the identification of individual Pol II, nucleosome, and transcription factor footprints, revealing Pol II pausing-driven destabilization of downstream nucleosomes. Furthermore, we demonstrate pervasive direct distance-dependent transcriptional coupling between nearby Pol II genes, Pol III genes, and transcribed enhancers, modulated by local chromatin architecture. Overall, transcription initiation reshapes surrounding nucleosome architecture and couples nearby transcriptional machinery along individual chromatin fibers.

Defining a chromatin architecture that supports transcription at RNA polymerase II promoters

Mammalian RNA polymerase II preinitiation complexes assemble adjacent to a nucleosome whose proximal edge (NPE) is typically 40 to 50 bp downstream of the transcription start site. At active promoters, that +1 nucleosome is universally modified by trimethylation on lysine 4 of histone H3 (H3K4me3). The Pol II preinitiation complex only extends 35 bp beyond the transcription start site, but nucleosomal templates with an NPE at +51 are nearly inactive in vitro with promoters that lack a TATA element and thus depend on TFIID for promoter recognition. Significantly, this inhibition is relieved when the +1 nucleosome contains H3K4me3, which can interact with TFIID subunits. Here, we show that H3K4me3 templates with both TATA and TATA-less promoters are active with +35 NPEs when transcription is driven by TFIID. Templates with +20 NPE are also active but at reduced levels compared to +35 and +51 NPEs, consistent with a general inhibition of promoter function when the proximal nucleosome encroaches on the preinitiation complex. Remarkably, dinucleosome templates support transcription when H3K4me3 is only present in the distal nucleosome, suggesting that TFIID-H3K4me3 interaction does not require modification of the +1 nucleosome. Transcription reactions performed with an alternative protocol retaining most nuclear factors results primarily in early termination, with a minority of complexes successfully traversing the first nucleosome. In such reactions, the +1 nucleosome does not substantially affect the level of termination even with an NPE of +20, indicating that a nucleosome barrier is not a major driver of early termination by Pol II.

guidedNOMe-seq quantifies chromatin states at single allele resolution for hundreds of custom regions in parallel

Since the introduction of next generation sequencing technologies, the field of epigenomics has evolved rapidly. However, most commonly used assays are enrichment-based methods and thus only semi-quantitative. Nucleosome occupancy and methylome sequencing (NOMe-seq) allows for quantitative inference of chromatin states with single locus resolution, but this requires high sequencing depth and is therefore prohibitively expensive to routinely apply to organisms with large genomes. To overcome this limitation, we introduce guidedNOMe-seq, where we combine NOMe profiling with large scale sgRNA synthesis and Cas9-mediated region-of-interest (ROI) liberation. To facilitate quantitative comparisons between multiple samples, we additionally develop an R package to standardize differential analysis of any type of NOMe-seq data. We extensively benchmark guidedNOMe-seq in a proof-of-concept study, dissecting the interplay of ChAHP and CTCF on chromatin. In summary we present a cost-effective, scalable, and customizable target enrichment extension to the existing NOMe-seq protocol allowing genome-scale quantification of nucleosome occupancy and transcription factor binding at single allele resolution.

Basis of gene-specific transcription regulation by the Integrator complex

The Integrator complex attenuates gene expression via the premature termination of RNA polymerase II (RNAP2) at promoter-proximal pausing sites. It is required for stimulus response, cell differentiation, and neurodevelopment, but how gene-specific and adaptive regulation by Integrator is achieved remains unclear. Here, we identify two sites on human Integrator subunits 13/14 that serve as binding hubs for sequence-specific transcription factors (TFs) and other transcription effector complexes. When Integrator is attached to paused RNAP2, these hubs are positioned upstream of the transcription bubble, consistent with simultaneous TF-promoter tethering. The TFs co-localize with Integrator genome-wide, increase Integrator abundance on target genes, and co-regulate responsive transcriptional programs. For instance, sensory cilia formation induced by glucose starvation depends on Integrator-TF contacts. Our data suggest TF-mediated promoter recruitment of Integrator as a widespread mechanism for targeted transcription regulation.

Global control of RNA polymerase II

RNA polymerase II (Pol II) is the multi-protein complex responsible for transcribing all protein-coding messenger RNA (mRNA). Most research on gene regulation is focused on the mechanisms controlling which genes are transcribed when, or on the mechanics of transcription. How global Pol II activity is determined receives comparatively less attention. Here, we follow the life of a Pol II molecule from 'assembly of the complex' to nuclear import, enzymatic activity, and degradation. We focus on how Pol II spends its time in the nucleus, and on the two-way relationship between Pol II abundance and activity in the context of homeostasis and global transcriptional changes.

Genome-wide distribution of 5-hydroxymethyluracil and chromatin accessibility in the Breviolum minutum genome

BACKGROUND

In dinoflagellates, a unique and extremely divergent genomic and nuclear organization has evolved. The highly unusual features of dinoflagellate nuclei and genomes include permanently condensed liquid crystalline chromosomes, primarily packaged by proteins other than histones, genes organized in very long unidirectional gene arrays, a general absence of transcriptional regulation, high abundance of the otherwise very rare DNA modification 5-hydroxymethyluracil (5-hmU), and many others. While most of these fascinating properties are originally identified in the 1970s and 1980s, they have not yet been investigated using modern genomic tools.

RESULTS

In this work, we address some of the outstanding questions regarding dinoflagellate genome organization by mapping the genome-wide distribution of 5-hmU (using both immunoprecipitation-based and basepair-resolution chemical mapping approaches) and of chromatin accessibility in the genome of the Symbiodiniaceae dinoflagellate Breviolum minutum. We find that the 5-hmU modification is preferentially enriched over certain classes of repetitive elements, often coincides with the boundaries between gene arrays, and is generally correlated with decreased chromatin accessibility, the latter otherwise being largely uniform along the genome. We discuss the potential roles of 5-hmU in the functional organization of dinoflagellate genomes and its relationship to the transcriptional landscape of gene arrays.

CONCLUSIONS

Our results provide the first window into the 5-hmU and chromatin accessibility landscapes in dinoflagellates.

NELF focuses sites of initiation and maintains promoter architecture

Many factors control the elongation phase of transcription by RNA polymerase II (Pol II), a process that plays an essential role in regulating gene expression. We utilized cells expressing degradation tagged subunits of NELFB, PAF1 and RTF1 to probe the effects of depletion of the factors on nascent transcripts using PRO-Seq and on chromatin architecture using DFF-ChIP. Although NELF is involved in promoter proximal pausing, depletion of NELFB had only a minimal effect on the level of paused transcripts and almost no effect on control of productive elongation. Instead, NELF depletion increased the utilization of downstream transcription start sites and caused a dramatic, genome-wide loss of H3K4me3 marked nucleosomes. Depletion of PAF1 and RTF1 both had major effects on productive transcript elongation in gene bodies and also caused initiation site changes like those seen with NELFB depletion. Our study confirmed that the first nucleosome encountered during initiation and early elongation is highly positioned with respect to the major TSS. In contrast, the positions of H3K4me3 marked nucleosomes in promoter regions are heterogeneous and are influenced by transcription. We propose a model defining NELF function and a general role of the H3K4me3 modification in blocking transcription initiation.

The role of O-GlcNAcylation in RNA polymerase II transcription

Eukaryotic RNA polymerase II (RNAPII) is responsible for the transcription of the protein-coding genes in the cell. Enormous progress has been made in discovering the protein activities that are required for transcription to occur, but the effects of post-translational modifications (PTMs) on RNAPII transcriptional regulation are much less understood. Most of our understanding relates to the cyclin-dependent kinases (CDKs), which appear to act relatively early in transcription. However, it is becoming apparent that other PTMs play a crucial role in the transcriptional cycle, and it is doubtful that any sort of complete understanding of this regulation is attainable without understanding the spectra of PTMs that occur on the transcriptional machinery. Among these is O-GlcNAcylation. Recent experiments have shown that the O-GlcNAc PTM likely has a prominent role in transcription. This review will cover the role of the O-GlcNAcylation in RNAPII transcription during initiation, pausing, and elongation, which will hopefully be of interest to both O-GlcNAc and RNAPII transcription researchers.

Power-law behavior of transcriptional bursting regulated by enhancer-promoter communication

Revealing how transcriptional bursting kinetics are genomically encoded is challenging because genome structures are stochastic at the organization level and are suggestively linked to gene transcription. To address this challenge, we develop a generic theoretical framework that integrates chromatin dynamics, enhancer-promoter (E-P) communication, and gene-state switching to study transcriptional bursting. The theory predicts that power law can be a general rule to quantitatively describe bursting modulations by E-P spatial communication. Specifically, burst frequency and burst size are up-regulated by E-P communication strength, following power laws with positive exponents. Analysis of the scaling exponents further reveals that burst frequency is preferentially regulated. Bursting kinetics are down-regulated by E-P genomic distance with negative power-law exponents, and this negative modulation desensitizes at large distances. The mutual information between burst frequency (or burst size) and E-P spatial distance further reveals essential characteristics of the information transfer from E-P communication to transcriptional bursting kinetics. These findings, which are in agreement with experimental observations, not only reveal fundamental principles of E-P communication in transcriptional bursting but also are essential for understanding cellular decision-making.

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.

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.

O-GlcNAcylation and Its Roles in Neurodegenerative Diseases

 As a non-classical post-translational modification, O-linked β-N-acetylglucosamine (O-GlcNAc) modification (O-GlcNAcylation) is widely found in human organ systems, particularly in our brains, and is indispensable for healthy cell biology. With the increasing age of the global population, the incidence of neurodegenerative diseases is increasing, too. The common characteristic of these disorders is the aggregation of abnormal proteins in the brain. Current research has found that O-GlcNAcylation dysregulation is involved in misfolding or aggregation of these abnormal proteins to mediate disease progression, but the specific mechanism has not been defined. This paper reviews recent studies on O-GlcNAcylation's roles in several neurodegenerative disorders such as Alzheimer's disease, Parkinson's disease, amyotrophic lateral sclerosis, Huntington's disease, Machado-Joseph's disease, and giant axonal neuropathy, and shows that O-GlcNAcylation, as glucose metabolism sensor, mediating synaptic function, participating in oxidative stress response and signaling pathway conduction, directly or indirectly regulates characteristic pathological protein toxicity and affects disease progression. The existing results suggest that targeting O-GlcNAcylation will provide new ideas for clinical diagnosis, prevention, and treatment of neurodegenerative diseases.

RNA polymerases reshape chromatin and coordinate transcription on individual fibers

During eukaryotic transcription, RNA polymerases must initiate and pause within a crowded, complex environment, surrounded by nucleosomes and other transcriptional activity. This environment creates a spatial arrangement along individual chromatin fibers ripe for both competition and coordination, yet these relationships remain largely unknown owing to the inherent limitations of traditional structural and sequencing methodologies. To address these limitations, we employed long-read chromatin fiber sequencing (Fiber-seq) to visualize RNA polymerases within their native chromatin context at single-molecule and near single-nucleotide resolution along up to 30 kb fibers. We demonstrate that Fiber-seq enables the identification of single-molecule RNA Polymerase (Pol) II and III transcription associated footprints, which, in aggregate, mirror bulk short-read sequencing-based measurements of transcription. We show that Pol II pausing destabilizes downstream nucleosomes, with frequently paused genes maintaining a short-term memory of these destabilized nucleosomes. Furthermore, we demonstrate pervasive direct coordination and anti-coordination between nearby Pol II genes, Pol III genes, transcribed enhancers, and insulator elements. This coordination is largely limited to spatially organized elements within 5 kb of each other, implicating short-range chromatin environments as a predominant determinant of coordinated polymerase initiation. Overall, transcription initiation reshapes surrounding nucleosome architecture and coordinates nearby transcriptional machinery along individual chromatin fibers.

Single molecule MATAC-seq reveals key determinants of DNA replication origin efficiency

Stochastic origin activation gives rise to significant cell-to-cell variability in the pattern of genome replication. The molecular basis for heterogeneity in efficiency and timing of individual origins is a long-standing question. Here, we developed Methylation Accessibility of TArgeted Chromatin domain Sequencing (MATAC-Seq) to determine single-molecule chromatin accessibility of four specific genomic loci. MATAC-Seq relies on preferential modification of accessible DNA by methyltransferases combined with Nanopore-Sequencing for direct readout of methylated DNA-bases. Applying MATAC-Seq to selected early-efficient and late-inefficient yeast replication origins revealed large heterogeneity of chromatin states. Disruption of INO80 or ISW2 chromatin remodeling complexes leads to changes at individual nucleosomal positions that correlate with changes in their replication efficiency. We found a chromatin state with an accessible nucleosome-free region in combination with well-positioned +1 and +2 nucleosomes as a strong predictor for efficient origin activation. Thus, MATAC-Seq identifies the large spectrum of alternative chromatin states that co-exist on a given locus previously masked in population-based experiments and provides a mechanistic basis for origin activation heterogeneity during eukaryotic DNA replication. Consequently, our single-molecule chromatin accessibility assay will be ideal to define single-molecule heterogeneity across many fundamental biological processes such as transcription, replication, or DNA repair in vitro and ex vivo.

Model-based characterization of the equilibrium dynamics of transcription initiation and promoter-proximal pausing in human cells

In metazoans, both transcription initiation and the escape of RNA polymerase (RNAP) from promoter-proximal pausing are key rate-limiting steps in gene expression. These processes play out at physically proximal sites on the DNA template and appear to influence one another through steric interactions. Here, we examine the dynamics of these processes using a combination of statistical modeling, simulation, and analysis of real nascent RNA sequencing data. We develop a simple probabilistic model that jointly describes the kinetics of transcription initiation, pause-escape, and elongation, and the generation of nascent RNA sequencing read counts under steady-state conditions. We then extend this initial model to allow for variability across cells in promoter-proximal pause site locations and steric hindrance of transcription initiation from paused RNAPs. In an extensive series of simulations, we show that this model enables accurate estimation of initiation and pause-escape rates. Furthermore, we show by simulation and analysis of real data that pause-escape is often strongly rate-limiting and that steric hindrance can dramatically reduce initiation rates. Our modeling framework is applicable to a variety of inference problems, and our software for estimation and simulation is freely available.

RNA polymerase II dynamics and mRNA stability feedback scale mRNA amounts with cell size

A fundamental feature of cellular growth is that total protein and RNA amounts increase with cell size to keep concentrations approximately constant. A key component of this is that global transcription rates increase in larger cells. Here, we identify RNA polymerase II (RNAPII) as the limiting factor scaling mRNA transcription with cell size in budding yeast, as transcription is highly sensitive to the dosage of RNAPII but not to other components of the transcriptional machinery. Our experiments support a dynamic equilibrium model where global RNAPII transcription at a given size is set by the mass action recruitment kinetics of unengaged nucleoplasmic RNAPII to the genome. However, this only drives a sub-linear increase in transcription with size, which is then partially compensated for by a decrease in mRNA decay rates as cells enlarge. Thus, limiting RNAPII and feedback on mRNA stability work in concert to scale mRNA amounts with cell size.

The chromatin landscape of the euryarchaeon Haloferax volcanii

BACKGROUND

Archaea, together with Bacteria, represent the two main divisions of life on Earth, with many of the defining characteristics of the more complex eukaryotes tracing their origin to evolutionary innovations first made in their archaeal ancestors. One of the most notable such features is nucleosomal chromatin, although archaeal histones and chromatin differ significantly from those of eukaryotes, not all archaea possess histones and it is not clear if histones are a main packaging component for all that do. Despite increased interest in archaeal chromatin in recent years, its properties have been little studied using genomic tools.

RESULTS

Here, we adapt the ATAC-seq assay to archaea and use it to map the accessible landscape of the genome of the euryarchaeote Haloferax volcanii. We integrate the resulting datasets with genome-wide maps of active transcription and single-stranded DNA (ssDNA) and find that while H. volcanii promoters exist in a preferentially accessible state, unlike most eukaryotes, modulation of transcriptional activity is not associated with changes in promoter accessibility. Applying orthogonal single-molecule footprinting methods, we quantify the absolute levels of physical protection of H. volcanii and find that Haloferax chromatin is similarly or only slightly more accessible, in aggregate, than that of eukaryotes. We also evaluate the degree of coordination of transcription within archaeal operons and make the unexpected observation that some CRISPR arrays are associated with highly prevalent ssDNA structures.

CONCLUSIONS

Our results provide the first comprehensive maps of chromatin accessibility and active transcription in Haloferax across conditions and thus a foundation for future functional studies of archaeal chromatin.

The landscape of the histone-organized chromatin of Bdellovibrionota bacteria

Histone proteins have traditionally been thought to be restricted to eukaryotes and most archaea, with eukaryotic nucleosomal histones deriving from their archaeal ancestors. In contrast, bacteria lack histones as a rule. However, histone proteins have recently been identified in a few bacterial clades, most notably the phylum Bdellovibrionota, and these histones have been proposed to exhibit a range of divergent features compared to histones in archaea and eukaryotes. However, no functional genomic studies of the properties of Bdellovibrionota chromatin have been carried out. In this work, we map the landscape of chromatin accessibility, active transcription and three-dimensional genome organization in a member of Bdellovibrionota (a Bacteriovorax strain). We find that, similar to what is observed in some archaea and in eukaryotes with compact genomes such as yeast, Bacteriovorax chromatin is characterized by preferential accessibility around promoter regions. Similar to eukaryotes, chromatin accessibility in Bacteriovorax positively correlates with gene expression. Mapping active transcription through single-strand DNA (ssDNA) profiling revealed that unlike in yeast, but similar to the state of mammalian and fly promoters, Bacteriovorax promoters exhibit very strong polymerase pausing. Finally, similar to that of other bacteria without histones, the Bacteriovorax genome exists in a three-dimensional (3D) configuration organized by the parABS system along the axis defined by replication origin and termination regions. These results provide a foundation for understanding the chromatin biology of the unique Bdellovibrionota bacteria and the functional diversity in chromatin organization across the tree of life.

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.

Genome-wide distribution of 5-hydroxymethyluracil and chromatin accessibility in the Breviolum minutum genome

In dinoflagellates, a unique and extremely divergent genomic and nuclear organization has evolved. The highly unusual features of dinoflagellate nuclei and genomes include permanently condensed liquid crystalline chromosomes, primarily packaged by proteins other than histones, genes organized in very long unidirectional gene arrays, a general absence of transcriptional regulation, high abundance of the otherwise very rare DNA modification 5-hydroxymethyluracil (5-hmU), and many others. While most of these fascinating properties were originally identified in the 1970s and 1980s, they have not yet been investigated using modern genomic tools. In this work, we address some of the outstanding questions regarding dinoflagellate genome organization by mapping the genome-wide distribution of 5-hmU (using both immunoprecipitation-based and basepair-resolution chemical mapping approaches) and of chromatin accessibility in the genome of the Symbiodiniaceae dinoflagellate Breviolum minutum. We find that the 5-hmU modification is preferentially enriched over certain classes of repetitive elements, often coincides with the boundaries between gene arrays, and is generally correlated with decreased chromatin accessibility, the latter otherwise being largely uniform along the genome. We discuss the potential roles of 5-hmU in the functional organization of dinoflagellate genomes and its relationship to the transcriptional landscape of gene arrays.

BurstDECONV: a signal deconvolution method to uncover mechanisms of transcriptional bursting in live cells

Monitoring transcription in living cells gives access to the dynamics of this complex fundamental process. It reveals that transcription is discontinuous, whereby active periods (bursts) are separated by one or several types of inactive periods of distinct lifetimes. However, decoding temporal fluctuations arising from live imaging and inferring the distinct transcriptional steps eliciting them is a challenge. We present BurstDECONV, a novel statistical inference method that deconvolves signal traces into individual transcription initiation events. We use the distribution of waiting times between successive polymerase initiation events to identify mechanistic features of transcription such as the number of rate-limiting steps and their kinetics. Comparison of our method to alternative methods emphasizes its advantages in terms of precision and flexibility. Unique features such as the direct determination of the number of promoter states and the simultaneous analysis of several potential transcription models make BurstDECONV an ideal analytic framework for live cell transcription imaging experiments. Using simulated realistic data, we found that our method is robust with regards to noise or suboptimal experimental designs. To show its generality, we applied it to different biological contexts such as Drosophila embryos or human cells.

Regulation of Pol II Pausing during Daily Gene Transcription in Mouse Liver

Cell autonomous circadian oscillation is present in central and various peripheral tissues. The intrinsic tissue clock and various extrinsic cues drive gene expression rhythms. Transcription regulation is thought to be the main driving force for gene rhythms. However, how transcription rhythms arise remains to be fully characterized due to the fact that transcription is regulated at multiple steps. In particular, Pol II recruitment, pause release, and premature transcription termination are critical regulatory steps that determine the status of Pol II pausing and transcription output near the transcription start site (TSS) of the promoter. Recently, we showed that Pol II pausing exhibits genome-wide changes during daily transcription in mouse liver. In this article, we review historical as well as recent findings on the regulation of transcription rhythms by the circadian clock and other transcription factors, and the potential limitations of those results in explaining rhythmic transcription at the TSS. We then discuss our results on the genome-wide characteristics of daily changes in Pol II pausing, the possible regulatory mechanisms involved, and their relevance to future research on circadian transcription regulation.

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.

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.

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.

Recycling of modified H2A-H2B provides short-term memory of chromatin states

Chromatin landscapes are disrupted during DNA replication and must be restored faithfully to maintain genome regulation and cell identity. The histone H3-H4 modification landscape is restored by parental histone recycling and modification of new histones. How DNA replication impacts on histone H2A-H2B is currently unknown. Here, we measure H2A-H2B modifications and H2A.Z during DNA replication and across the cell cycle using quantitative genomics. We show that H2AK119ub1, H2BK120ub1, and H2A.Z are recycled accurately during DNA replication. Modified H2A-H2B are segregated symmetrically to daughter strands via POLA1 on the lagging strand, but independent of H3-H4 recycling. Post-replication, H2A-H2B modification and variant landscapes are quickly restored, and H2AK119ub1 guides accurate restoration of H3K27me3. This work reveals epigenetic transmission of parental H2A-H2B during DNA replication and identifies cross talk between H3-H4 and H2A-H2B modifications in epigenome propagation. We propose that rapid short-term memory of recycled H2A-H2B modifications facilitates restoration of stable H3-H4 chromatin states.

A CpG island-encoded mechanism protects genes from premature transcription termination

Transcription must be tightly controlled to regulate gene expression and development. However, our understanding of the molecular mechanisms that influence transcription and how these are coordinated in cells to ensure normal gene expression remains rudimentary. Here, by dissecting the function of the SET1 chromatin-modifying complexes that bind to CpG island-associated gene promoters, we discover that they play a specific and essential role in enabling the expression of low to moderately transcribed genes. Counterintuitively, this effect can occur independently of SET1 complex histone-modifying activity and instead relies on an interaction with the RNA Polymerase II-binding protein WDR82. Unexpectedly, we discover that SET1 complexes enable gene expression by antagonising premature transcription termination by the ZC3H4/WDR82 complex at CpG island-associated genes. In contrast, at extragenic sites of transcription, which typically lack CpG islands and SET1 complex occupancy, we show that the activity of ZC3H4/WDR82 is unopposed. Therefore, we reveal a gene regulatory mechanism whereby CpG islands are bound by a protein complex that specifically protects genic transcripts from premature termination, effectively distinguishing genic from extragenic transcription and enabling normal gene expression.

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.

Cohesin and CTCF control the dynamics of chromosome folding

In mammals, interactions between sequences within topologically associating domains enable control of gene expression across large genomic distances. Yet it is unknown how frequently such contacts occur, how long they last and how they depend on the dynamics of chromosome folding and loop extrusion activity of cohesin. By imaging chromosomal locations at high spatial and temporal resolution in living cells, we show that interactions within topologically associating domains are transient and occur frequently during the course of a cell cycle. Interactions become more frequent and longer in the presence of convergent CTCF sites, resulting in suppression of variability in chromosome folding across time. Supported by physical models of chromosome dynamics, our data suggest that CTCF-anchored loops last around 10 min. Our results show that long-range transcriptional regulation might rely on transient physical proximity, and that cohesin and CTCF stabilize highly dynamic chromosome structures, facilitating selected subsets of chromosomal interactions.

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.

Modification mapping by nanopore sequencing

Next generation sequencing (NGS) has provided biologists with an unprecedented view into biological processes and their regulation over the past 2 decades, fueling a wave of development of high throughput methods based on short read DNA and RNA sequencing. For nucleic acid modifications, NGS has been coupled with immunoprecipitation, chemical treatment, enzymatic treatment, and/or the use of reverse transcriptase enzymes with fortuitous activities to enrich for and to identify covalent modifications of RNA and DNA. However, the majority of nucleic acid modifications lack commercial monoclonal antibodies, and mapping techniques that rely on chemical or enzymatic treatments to manipulate modification signatures add additional technical complexities to library preparation. Moreover, such approaches tend to be specific to a single class of RNA or DNA modification, and generate only indirect readouts of modification status. Third generation sequencing technologies such as the commercially available "long read" platforms from Pacific Biosciences and Oxford Nanopore Technologies are an attractive alternative for high throughput detection of nucleic acid modifications. While the former can indirectly sense modified nucleotides through changes in the kinetics of reverse transcription reactions, nanopore sequencing can in principle directly detect any nucleic acid modification that produces a signal distortion as the nucleic acid passes through a nanopore sensor embedded within a charged membrane. To date, more than a dozen endogenous DNA and RNA modifications have been interrogated by nanopore sequencing, as well as a number of synthetic nucleic acid modifications used in metabolic labeling, structure probing, and other emerging applications. This review is intended to introduce the reader to nanopore sequencing and key principles underlying its use in direct detection of nucleic acid modifications in unamplified DNA or RNA samples, and outline current approaches for detecting and quantifying nucleic acid modifications by nanopore sequencing. As this technology matures, we anticipate advances in both sequencing chemistry and analysis methods will lead to rapid improvements in the identification and quantification of these epigenetic marks.

Xrn2 substrate mapping identifies torpedo loading sites and extensive premature termination of RNA pol II transcription

The exonuclease torpedo Xrn2 loads onto nascent RNA 5'-PO₄ ends and chases down pol II to promote termination downstream from polyA sites. We report that Xrn2 is recruited to preinitiation complexes and "travels" to 3' ends of genes. Mapping of 5'-PO₄ ends in nascent RNA identified Xrn2 loading sites stabilized by an active site mutant, Xrn2(D235A). Xrn2 loading sites are approximately two to 20 bases downstream from where CPSF73 cleaves at polyA sites and histone 3' ends. We propose that processing of all mRNA 3' ends comprises cleavage and limited 5'-3' trimming by CPSF73, followed by handoff to Xrn2. A similar handoff occurs at tRNA 3' ends, where cotranscriptional RNase Z cleavage generates novel Xrn2 substrates. Exonuclease-dead Xrn2 increased transcription in 3' flanking regions by inhibiting polyA site-dependent termination. Surprisingly, the mutant Xrn2 also rescued transcription in promoter-proximal regions to the same extent as in 3' flanking regions. eNET-seq revealed Xrn2-mediated degradation of sense and antisense nascent RNA within a few bases of the TSS, where 5'-PO₄ ends may be generated by decapping or endonucleolytic cleavage. These results suggest that a major fraction of pol II complexes terminates prematurely close to the start site under normal conditions by an Xrn2-mediated torpedo mechanism.

BRD4: a general regulator of transcription elongation

Transcription elongation by RNA polymerase II (Pol II) has emerged as a regulatory hub in gene expression. A key control point occurs during early transcription elongation when Pol II pauses in the promoter-proximal region at the majority of genes in mammalian cells and at a large set of genes in Drosophila. An increasing number of trans-acting factors have been linked to promoter-proximal pausing. Some factors help to establish the pause, whereas others are required for the release of Pol II into productive elongation. A dysfunction of this elongation control point leads to aberrant gene expression and can contribute to disease development. The BET bromodomain protein BRD4 has been implicated in elongation control. However, only recently direct BRD4-specific functions in Pol II transcription elongation have been uncovered. This mainly became possible with technological advances that allow selective and rapid ablation of BRD4 in cells along with the availability of approaches that capture the immediate consequences on nascent transcription. This review sheds light on the experimental breakthroughs that led to the emerging view of BRD4 as a general regulator of transcription elongation.

RNA Pol II pausing facilitates phased pluripotency transitions by buffering transcription

Promoter-proximal RNA Pol II pausing is a critical step in transcriptional control. Pol II pausing has been predominantly studied in tissue culture systems. While Pol II pausing has been shown to be required for mammalian development, the phenotypic and mechanistic details of this requirement are unknown. Here, we found that loss of Pol II pausing stalls pluripotent state transitions within the epiblast of the early mouse embryo. Using Nelfb -/- mice and a NELFB degron mouse pluripotent stem cell model, we show that embryonic stem cells (ESCs) representing the naïve state of pluripotency successfully initiate a transition program but fail to balance levels of induced and repressed genes and enhancers in the absence of NELF. We found an increase in chromatin-associated NELF during transition from the naïve to later pluripotent states. Overall, our work defines the acute and long-term molecular consequences of NELF loss and reveals a role for Pol II pausing in the pluripotency continuum as a modulator of cell state transitions.

Structural basis of mRNA maturation: Time to put it together

In eukaryotes, the expression of genetic information begins in the cell nucleus with precursor messenger RNA (pre-mRNA) transcription and processing into mature mRNA. The mRNA is subsequently recognized and packaged by proteins into an mRNA ribonucleoprotein complex (mRNP) and exported to the cytoplasm for translation. Each of the nuclear mRNA maturation steps is carried out by a dedicated molecular machine. Here, we highlight recent structural and mechanistic insights into how these machines function, including the capping enzyme, the spliceosome, the 3'-end processing machinery, and the transcription-export complex. While we increasingly understand individual steps of nuclear gene expression, many questions remain. For example, we are only beginning to reveal how mature mRNAs are recognized and packaged for nuclear export and how mRNA maturation events are coupled to transcription and to each other. Advances in the preparation of recombinant and endogenous protein-nucleic acid complexes, cryo-electron microscopy, and machine learning promise exciting insights into the mechanisms of nuclear gene expression and its spatial organization.

Transcriptional lockdown during acute proteotoxic stress

Cells experiencing proteotoxic stress downregulate the expression of thousands of active genes and upregulate a few stress-response genes. The strategy of downregulating gene expression has conceptual parallels with general lockdown in the global response to the coronavirus disease 2019 (COVID-19) pandemic. The mechanistic details of global transcriptional downregulation of genes, termed stress-induced transcriptional attenuation (SITA), are only beginning to emerge. The reduction in RNA and protein production during stress may spare proteostasis capacity, allowing cells to divert resources to control stress-induced damage. Given the relevance of translational downregulation in a broad variety of diseases, the role of SITA in diseases caused by proteotoxicity should be investigated in future, paving the way for potential novel therapeutics.

DNA damage-induced transcription stress triggers the genome-wide degradation of promoter-bound Pol II

The precise regulation of RNA Polymerase II (Pol II) transcription after genotoxic stress is crucial for proper execution of the DNA damage-induced stress response. While stalling of Pol II on transcription-blocking lesions (TBLs) blocks transcript elongation and initiates DNA repair in cis, TBLs additionally elicit a response in trans that regulates transcription genome-wide. Here we uncover that, after an initial elongation block in cis, TBLs trigger the genome-wide VCP-mediated proteasomal degradation of promoter-bound, P-Ser5-modified Pol II in trans. This degradation is mechanistically distinct from processing of TBL-stalled Pol II, is signaled via GSK3, and contributes to the TBL-induced transcription block, even in transcription-coupled repair-deficient cells. Thus, our data reveal the targeted degradation of promoter-bound Pol II as a critical pathway that allows cells to cope with DNA damage-induced transcription stress and enables the genome-wide adaptation of transcription to genotoxic stress.

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.

A computational pipeline to visualize DNA-protein binding states using dSMF data

Here, we present a pipeline to map states of protein-binding DNA in vivo. Our pipeline infers as well as quantifies cooperative binding. Using dual-enzyme single-molecule footprinting (dSMF) data, we show how our workflow identifies binding states at an enhancer in Drosophila S2 cells. Data from cells lacking endogenous DNA methylation are a prerequisite for this pipeline.

For complete details on the use and execution of this protocol, please refer to Rao et al. (2021) and Krebs et al. (2017).

•

Pipeline that uses dSMF data to quantify DNA-protein binding states

•

This pipeline can reliably map cobinding events on chomatinized DNA

•

Unlike MNase- and DNase-seq, dSMF maps the unbound state of genomic DNA

Nonlinear control of transcription through enhancer-promoter interactions

Chromosome structure in mammals is thought to regulate transcription by modulating three-dimensional interactions between enhancers and promoters, notably through CTCF-mediated loops and topologically associating domains (TADs)1-4. However, how chromosome interactions are actually translated into transcriptional outputs remains unclear. Here, to address this question, we use an assay to position an enhancer at large numbers of densely spaced chromosomal locations relative to a fixed promoter, and measure promoter output and interactions within a genomic region with minimal regulatory and structural complexity. A quantitative analysis of hundreds of cell lines reveals that the transcriptional effect of an enhancer depends on its contact probabilities with the promoter through a nonlinear relationship. Mathematical modelling suggests that nonlinearity might arise from transient enhancer-promoter interactions being translated into slower promoter bursting dynamics in individual cells, therefore uncoupling the temporal dynamics of interactions from those of transcription. This uncovers a potential mechanism of how distal enhancers act from large genomic distances, and of how topologically associating domain boundaries block distal enhancers. Finally, we show that enhancer strength also determines absolute transcription levels as well as the sensitivity of a promoter to CTCF-mediated transcriptional insulation. Our measurements establish general principles for the context-dependent role of chromosome structure in long-range transcriptional regulation.