Global identification of SWI/SNF targets reveals compensation by EP400.

Mammalian SWI/SNF chromatin remodeling complexes move and evict nucleosomes at gene promoters and enhancers to modulate DNA access. Although SWI/SNF subunits are commonly mutated in disease, therapeutic options are limited by our inability to predict SWI/SNF gene targets and conflicting studies on functional significance. Here, we leverage a fast-acting inhibitor of SWI/SNF remodeling to elucidate direct targets and effects of SWI/SNF. Blocking SWI/SNF activity causes a rapid and global loss of chromatin accessibility and transcription. Whereas repression persists at most enhancers, we uncover a compensatory role for the EP400/TIP60 remodeler, which reestablishes accessibility at most promoters during prolonged loss of SWI/SNF. Indeed, we observe synthetic lethality between EP400 and SWI/SNF in cancer cell lines and human cancer patient data. Our data define a set of molecular genomic features that accurately predict gene sensitivity to SWI/SNF inhibition in diverse cancer cell lines, thereby improving the therapeutic potential of SWI/SNF inhibitors.

Human SMARCA5 is continuously required to maintain nucleosome spacing.

Genetic models suggested that SMARCA5 was required for DNA-templated events including transcription, DNA replication, and DNA repair. We engineered a degron tag into the endogenous alleles of SMARCA5, a catalytic component of the imitation switch complexes in three different human cell lines to define the effects of rapid degradation of this key regulator. Degradation of SMARCA5 was associated with a rapid increase in global nucleosome repeat length, which may allow greater chromatin compaction. However, there were few changes in nascent transcription within the first 6 h of degradation. Nevertheless, we demonstrated a requirement for SMARCA5 to control nucleosome repeat length at G1/S and during the S phase. SMARCA5 co-localized with CTCF and H2A.Z, and we found a rapid loss of CTCF DNA binding and disruption of nucleosomal phasing around CTCF binding sites. This spatiotemporal analysis indicates that SMARCA5 is continuously required for maintaining nucleosomal spacing.

Human histone H1 variants impact splicing outcome by controlling RNA polymerase II elongation.

Histones shape chromatin structure and the epigenetic landscape. H1, the most diverse histone in the human genome, has 11 variants. Due to the high structural similarity between the H1s, their unique functions in transferring information from the chromatin to mRNA-processing machineries have remained elusive. Here, we generated human cell lines lacking up to five H1 subtypes, allowing us to characterize the genomic binding profiles of six H1 variants. Most H1s bind to specific sites, and binding depends on multiple factors, including GC content. The highly expressed H1.2 has a high affinity for exons, whereas H1.3 binds intronic sequences. H1s are major splicing regulators, especially of exon skipping and intron retention events, through their effects on the elongation of RNA polymerase II (RNAPII). Thus, H1 variants determine splicing fate by modulating RNAPII elongation.

PAX3-FOXO1 coordinates enhancer architecture, eRNA transcription, and RNA polymerase pause release at select gene targets.

Transcriptional control is a highly dynamic process that changes rapidly in response to various cellular and extracellular cues, making it difficult to define the mechanism of transcription factor function using slow genetic methods. We used a chemical-genetic approach to rapidly degrade a canonical transcriptional activator, PAX3-FOXO1, to define the mechanism by which it regulates gene expression programs. By coupling rapid protein degradation with the analysis of nascent transcription over short time courses and integrating CUT&RUN, ATAC-seq, and eRNA analysis with deep proteomic analysis, we defined PAX3-FOXO1 function at a small network of direct transcriptional targets. PAX3-FOXO1 degradation impaired RNA polymerase pause release and transcription elongation at most regulated gene targets. Moreover, the activity of PAX3-FOXO1 at enhancers controlling this core network was surprisingly selective, affecting single elements in super-enhancers. This combinatorial analysis indicated that PAX3-FOXO1 was continuously required to maintain chromatin accessibility and enhancer architecture at regulated enhancers.

Definition of a small core transcriptional circuit regulated by AML1-ETO.

Transcription factors regulate gene networks controlling normal hematopoiesis and are frequently deregulated in acute myeloid leukemia (AML). Critical to our understanding of the mechanism of cellular transformation by oncogenic transcription factors is the ability to define their direct gene targets. However, gene network cascades can change within minutes to hours, making it difficult to distinguish direct from secondary or compensatory transcriptional changes by traditional methodologies. To overcome this limitation, we devised cell models in which the AML1-ETO protein could be quickly degraded upon addition of a small molecule. The rapid kinetics of AML1-ETO removal, when combined with analysis of transcriptional output by nascent transcript analysis and genome-wide AML1-ETO binding by CUT&RUN, enabled the identification of direct gene targets that constitute a core AML1-ETO regulatory network. Moreover, derepression of this gene network was associated with RUNX1 DNA binding and triggered a transcription cascade ultimately resulting in myeloid differentiation.

H3K27ac bookmarking promotes rapid post-mitotic activation of the pluripotent stem cell program without impacting 3D chromatin reorganization.

During self-renewal, cell-type-defining features are drastically perturbed in mitosis and must be faithfully reestablished upon G1 entry, a process that remains largely elusive. Here, we characterized at a genome-wide scale the dynamic transcriptional and architectural resetting of mouse pluripotent stem cells (PSCs) upon mitotic exit. We captured distinct waves of transcriptional reactivation with rapid induction of stem cell genes and transient activation of lineage-specific genes. Topological reorganization at different hierarchical levels also occurred in an asynchronous manner and showed partial coordination with transcriptional resetting. Globally, rapid transcriptional and architectural resetting associated with mitotic retention of H3K27 acetylation, supporting a bookmarking function. Indeed, mitotic depletion of H3K27ac impaired the early reactivation of bookmarked, stem-cell-associated genes. However, 3D chromatin reorganization remained largely unaffected, suggesting that these processes are driven by distinct forces upon mitotic exit. This study uncovers principles and mediators of PSC molecular resetting during self-renewal.

TFIID Enables RNA Polymerase II Promoter-Proximal Pausing.

RNA polymerase II (RNAPII) transcription is governed by the pre-initiation complex (PIC), which contains TFIIA, TFIIB, TFIID, TFIIE, TFIIF, TFIIH, RNAPII, and Mediator. After initiation, RNAPII enzymes pause after transcribing less than 100 bases; precisely how RNAPII pausing is enforced and regulated remains unclear. To address specific mechanistic questions, we reconstituted human RNAPII promoter-proximal pausing in vitro, entirely with purified factors (no extracts). As expected, NELF and DSIF increased pausing, and P-TEFb promoted pause release. Unexpectedly, the PIC alone was sufficient to reconstitute pausing, suggesting RNAPII pausing is an inherent PIC function. In agreement, pausing was lost upon replacement of the TFIID complex with TATA-binding protein (TBP), and PRO-seq experiments revealed widespread disruption of RNAPII pausing upon acute depletion (t = 60 min) of TFIID subunits in human or Drosophila cells. These results establish a TFIID requirement for RNAPII pausing and suggest pause regulatory factors may function directly or indirectly through TFIID.

NELF Regulates a Promoter-Proximal Step Distinct from RNA Pol II Pause-Release.

RNA polymerase II (RNA Pol II) is generally paused at promoter-proximal regions in most metazoans, and based on in vitro studies, this function has been attributed to the negative elongation factor (NELF). Here, we show that upon rapid depletion of NELF, RNA Pol II fails to be released into gene bodies, stopping instead around the +1 nucleosomal dyad-associated region. The transition to the 2nd pause region is independent of positive transcription elongation factor P-TEFb. During the heat shock response, RNA Pol II is rapidly released from pausing at heat shock-induced genes, while most genes are paused and transcriptionally downregulated. Both of these aspects of the heat shock response remain intact upon NELF loss. We find that NELF depletion results in global loss of cap-binding complex from chromatin without global reduction of nascent transcript 5' cap stability. Thus, our studies implicate NELF functioning in early elongation complexes distinct from RNA Pol II pause-release.

Protocol variations in run-on transcription dataset preparation produce detectable signatures in sequencing libraries.

BACKGROUND: A variety of protocols exist for producing whole genome run-on transcription datasets. However, little is known about how differences between these protocols affect the signal within the resulting libraries. RESULTS: Using run-on transcription datasets generated from the same biological system, we show that a variety of GRO- and PRO-seq preparation methods leave identifiable signatures within each library. Specifically we show that the library preparation method results in differences in quality control metrics, as well as differences in the signal distribution at the 5 ' end of transcribed regions. These shifts lead to disparities in eRNA identification, but do not impact analyses aimed at inferring the key regulators involved in changes to transcription. CONCLUSIONS: Run-on sequencing protocol variations result in technical signatures that can be used to identify both the enrichment and library preparation method of a particular data set. These technical signatures are batch effects that limit detailed comparisons of pausing ratios and eRNAs identified across protocols. However, these batch effects have only limited impact on our ability to infer which regulators underlie the observed transcriptional changes.

Characterizing RNA stability genome-wide through combined analysis of PRO-seq and RNA-seq data.

BACKGROUND: The concentrations of distinct types of RNA in cells result from a dynamic equilibrium between RNA synthesis and decay. Despite the critical importance of RNA decay rates, current approaches for measuring them are generally labor-intensive, limited in sensitivity, and/or disruptive to normal cellular processes. Here, we introduce a simple method for estimating relative RNA half-lives that is based on two standard and widely available high-throughput assays: Precision Run-On sequencing (PRO-seq) and RNA sequencing (RNA-seq). RESULTS: Our method treats PRO-seq as a measure of transcription rate and RNA-seq as a measure of RNA concentration, and estimates the rate of RNA decay required for a steady-state equilibrium. We show that this approach can be used to assay relative RNA half-lives genome-wide, with good accuracy and sensitivity for both coding and noncoding transcription units. Using a structural equation model (SEM), we test several features of transcription units, nearby DNA sequences, and nearby epigenomic marks for associations with RNA stability after controlling for their effects on transcription. We find that RNA splicing-related features are positively correlated with RNA stability, whereas features related to miRNA binding and DNA methylation are negatively correlated with RNA stability. Furthermore, we find that a measure based on U1 binding and polyadenylation sites distinguishes between unstable noncoding and stable coding transcripts but is not predictive of relative stability within the mRNA or lincRNA classes. We also identify several histone modifications that are associated with RNA stability. CONCLUSION: We introduce an approach for estimating the relative half-lives of individual RNAs. Together, our estimation method and systematic analysis shed light on the pervasive impacts of RNA stability on cellular RNA concentrations.

Use of the nuclear walk-on methodology to determine sites of RNA polymerase II initiation and pausing and quantify nascent RNAs in cells.

Transcription by RNA polymerase II (Pol II) is controlled during initiation, elongation, and termination by a large variety of transcription factors, the state of chromatin modifications, and environmental conditions. Herein we describe experimental approaches for the examination of Pol II transcription at semi-global and genome-wide scales through analysis of nascent Pol II transcripts. We begin with a description of the nuclear walk-on (NWO) assay, which involves rapid isolation of nuclei in the presence of EDTA, followed by extension of about a quarter of the nascent transcripts with 32P-CTP. Labeled nascent transcripts are then analyzed by denaturing PAGE and phosphorimaging followed by densitometry analysis to quantify the signal on the gel. A parallel reaction containing α-amanitin to inhibit Pol II reveals transcription due to Pol I and Pol III, which can be subtracted to yield a profile of Pol II transcription. We then describe how to use the NWO as a front end for PRO-Seq and PRO-Cap methods, which permit the genome-wide characterization of Pol II transcription at nucleotide resolution and provide precise information about sites of transcription initiation and pausing. We discuss strategies for optimizing sequencing methods that capture nascent Pol II transcripts, methods of bias reduction, and approaches for normalizing these and other sequencing datasets using spike-in controls.

Mutant FOXO1 controls an oncogenic network via enhancer accessibility.

Transcriptional dysregulation is a hallmark of diffuse large B cell lymphoma (DLBCL), as transcriptional regulators are frequently mutated. However, our mechanistic understanding of how normal transcriptional programs are co-opted in DLBCL has been hindered by a lack of methodologies that provide the temporal resolution required to separate direct and indirect effects on transcriptional control. We applied a chemical-genetic approach to engineer the inducible degradation of the transcription factor FOXO1, which is recurrently mutated (mFOXO1) in DLBCL. The combination of rapid degradation of mFOXO1, nascent transcript detection, and assessment of chromatin accessibility allowed us to identify the direct targets of mFOXO1. mFOXO1 was required to maintain accessibility at specific enhancers associated with multiple oncogenes, and mFOXO1 degradation impaired RNA polymerase pause-release at some targets. Wild-type FOXO1 appeared to weakly regulate many of the same targets as mFOXO1 and was able to complement the degradation of mFOXO1 in the context of AKT inhibition.

Mapping Transcription Regulation with Run-on and Sequencing Data Using the Web-Based tfTarget Gateway.

Run-on and sequencing assays like GRO-seq, PRO-seq, and ChRO-seq allow for joint profiling of transcription activity of transcriptional regulatory elements (TREs), i.e., promoters and active enhancers, and target genes. Variation in biological conditions, such as treated vs. control, results in changes in the activity of transcription factors (TFs), which induces concerted changes in TREs and target genes. By modeling the differences between two biological conditions, we developed the computational pipeline known as tfTarget that predicts a set of putative TREs and target genes responding to each TF under the biological condition of interest. In this chapter, we demonstrate the use of the new web-based tfTarget in mapping transcription regulation using run-on sequencing data.

Nascent transcription upon interferon-α2 stimulation on human and rhesus macaque lymphoblastoid cell lines.

OBJECTIVES: The interferon-triggered innate immune response has been observed to be under strong diversifying selection to counteract the many pathogens hosts have to defend against. In particular, rewiring of gene transcription regulation allows organisms to rapidly acquire new phenotypes by removing and adding genes into the innate immune gene network. Dissecting the molecular processes by which this rewiring takes place, either by changing the DNA regulatory elements or by changing the activity of the regulators across species, is key to better understand this evolutionary process. DATA DESCRIPTION: To better comprehend the evolutionary dynamics that have occurred in the initial transcriptional response to interferon in primates, we present Precision Run-On (PRO-seq) datasets made after 1 h of interferon-α2 stimulation on human and rhesus macaque lymphoblastoid cell lines. Further, we tested the difference between using either species' cognate interferon versus using the other orthologous interferon to account for any potential impacts in the interaction of the orthologous interferons with their cellular membrane receptors. This data provides insights into the regulatory mechanisms that drive species-specific responses to environmental perturbations, such as the one driven by the interactions of pathogens and their hosts.

PRO-seq: Precise Mapping of Engaged RNA Pol II at Single-Nucleotide Resolution.

Gene regulation is dependent on the production of mRNAs and a repertoire of non-coding RNAs by RNA polymerase II (RNAPII). Precision run-on sequencing (PRO-seq) maps the position of engaged RNAPII complexes at single-nucleotide resolution and can reveal direct targets of regulation, locations of enhancers, and transcription mechanisms that are difficult or impossible to measure by analysis of total cellular RNA. Briefly, this method first involves permeabilizing cells with mild detergents to remove intracellular NTPs and halt transcription. Transcription is then resumed in the presence of biotin-NTPs and sarkosyl to allow transcriptional incorporation of a single biotinylated NTP by RNAPII. The biotin moiety is then bound to streptavidin beads to stringently enrich for nascent RNAs. Sequencing libraries are then generated such that the first base read corresponds to the 3' end of the nascent transcript. Here, we describe our current protocol for generating PRO-seq libraries from metazoan cells, including adaptations of previously published protocols to incorporate unique molecular identifiers, reduce ligation bias, and improve library yields. Additional commentary describes quality control and processing of PRO-seq data and references for more advanced downstream analysis such as gene and enhancer identification. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Cell permeabilization for PRO-seq Basic Protocol 2: Construction of PRO-seq libraries Support Protocol: Adenylation of 3' adapter.

Acute depletion of human core nucleoporin reveals direct roles in transcription control but dispensability for 3D genome organization.

The nuclear pore complex (NPC) comprises more than 30 nucleoporins (NUPs) and is a hallmark of eukaryotes. NUPs have been suggested to be important in regulating gene transcription and 3D genome organization. However, evidence in support of their direct roles remains limited. Here, by Cut&Run, we find that core NUPs display broad but also cell-type-specific association with active promoters and enhancers in human cells. Auxin-mediated rapid depletion of two NUPs demonstrates that NUP93, but not NUP35, directly and specifically controls gene transcription. NUP93 directly activates genes with high levels of RNA polymerase II loading and transcriptional elongation by facilitating full BRD4 recruitment to their active enhancers. dCas9-based tethering confirms a direct and causal role of NUP93 in gene transcriptional activation. Unexpectedly, in situ Hi-C and H3K27ac or H3K4me1 HiChIP results upon acute NUP93 depletion show negligible changesS2211-1247(22)01437-1 of 3D genome organization ranging from A/B compartments and topologically associating domains (TADs) to enhancer-promoter contacts.

RNA polymerase mapping in plants identifies intergenic regulatory elements enriched in causal variants.

Control of gene expression is fundamental at every level of cell function. Promoter-proximal pausing and divergent transcription at promoters and enhancers, which are prominent features in animals, have only been studied in a handful of research experiments in plants. PRO-Seq analysis in cassava (Manihot esculenta) identified peaks of transcriptionally engaged RNA polymerase at both the 5' and 3' end of genes, consistent with paused or slowly moving Polymerase. In addition, we identified divergent transcription at intergenic sites. A full genome search for bi-directional transcription using an algorithm for enhancer detection developed in mammals (dREG) identified many intergenic regulatory element (IRE) candidates. These sites showed distinct patterns of methylation and nucleotide conservation based on genomic evolutionary rate profiling (GERP). SNPs within these IRE candidates explained significantly more variation in fitness and root composition than SNPs in chromosomal segments randomly ascertained from the same intergenic distribution, strongly suggesting a functional importance of these sites. Maize GRO-Seq data showed RNA polymerase occupancy at IREs consistent with patterns in cassava. Furthermore, these IREs in maize significantly overlapped with sites previously identified on the basis of open chromatin, histone marks, and methylation, and were enriched for reported eQTL. Our results suggest that bidirectional transcription can identify intergenic genomic regions in plants that play an important role in transcription regulation and whose identification has the potential to aid crop improvement.

The chromatin remodeler Ino80 mediates RNAPII pausing site determination.

BACKGROUND: Promoter-proximal pausing of RNA polymerase II (RNAPII) is a critical step for the precise regulation of gene expression. Despite the apparent close relationship between promoter-proximal pausing and nucleosome, the role of chromatin remodeler governing this step has mainly remained elusive. RESULTS: Here, we report highly confined RNAPII enrichments downstream of the transcriptional start site in Saccharomyces cerevisiae using PRO-seq experiments. This non-uniform distribution of RNAPII exhibits both similar and different characteristics with promoter-proximal pausing in Schizosaccharomyces pombe and metazoans. Interestingly, we find that Ino80p knockdown causes a significant upstream transition of promoter-proximal RNAPII for a subset of genes, relocating RNAPII from the main pausing site to the alternative pausing site. The proper positioning of RNAPII is largely dependent on nucleosome context. We reveal that the alternative pausing site is closely associated with the + 1 nucleosome, and nucleosome architecture around the main pausing site of these genes is highly phased. In addition, Ino80p knockdown results in an increase in fuzziness and a decrease in stability of the + 1 nucleosome. Furthermore, the loss of INO80 also leads to the shift of promoter-proximal RNAPII toward the alternative pausing site in mouse embryonic stem cells. CONCLUSIONS: Based on our collective results, we hypothesize that the highly conserved chromatin remodeler Ino80p is essential in establishing intact RNAPII pausing during early transcription elongation in various organisms, from budding yeast to mouse.

Distinct properties and functions of CTCF revealed by a rapidly inducible degron system.

CCCTC-binding factor (CTCF) is a conserved zinc finger transcription factor implicated in a wide range of functions, including genome organization, transcription activation, and elongation. To explore the basis for CTCF functional diversity, we coupled an auxin-induced degron system with precision nuclear run-on. Unexpectedly, oriented CTCF motifs in gene bodies are associated with transcriptional stalling in a manner independent of bound CTCF. Moreover, CTCF at different binding sites (CBSs) displays highly variable resistance to degradation. Motif sequence does not significantly predict degradation behavior, but location at chromatin boundaries and chromatin loop anchors, as well as co-occupancy with cohesin, are associated with delayed degradation. Single-molecule tracking experiments link chromatin residence time to CTCF degradation kinetics, which has ramifications regarding architectural CTCF functions. Our study highlights the heterogeneity of CBSs, uncovers properties specific to architecturally important CBSs, and provides insights into the basic processes of genome organization and transcription regulation.

Timing RNA polymerase pausing with TV-PRO-seq.

Transcription of many genes in metazoans is subject to polymerase pausing, which is the transient stop of transcriptionally engaged polymerases. This is known to mainly occur in promoter-proximal regions but it is not well understood. In particular, a genome-wide measurement of pausing times at high resolution has been lacking. We present here the time-variant precision nuclear run-on and sequencing (TV-PRO-seq) assay, an extension of the standard PRO-seq that allows us to estimate genome-wide pausing times at single-base resolution. Its application to human cells demonstrates that, proximal to promoters, polymerases pause more frequently but for shorter times than in other genomic regions. Comparison with single-cell gene expression data reveals that the polymerase pausing times are longer in highly expressed genes, while transcriptionally noisier genes have higher pausing frequencies and slightly longer pausing times. Analyses of histone modifications suggest that the marker H3K36me3 is related to the polymerase pausing.