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.

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.

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 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.

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.

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.

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.

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.

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.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-022-08352-8).

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.

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.

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.

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.

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.

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.

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.

A Role for FACT in RNA Polymerase II Promoter-Proximal Pausing

FACT (facilitates chromatin transcription) is an evolutionarily conserved histone chaperone that was initially identified as an activity capable of promoting RNA polymerase II (Pol II) transcription through nucleosomes in vitro. In this report, we describe a global analysis of FACT function in Pol II transcription in Drosophila. We present evidence that loss of FACT has a dramatic impact on Pol II elongation-coupled processes including histone H3 lysine 4 (H3K4) and H3K36 methylation, consistent with a role for FACT in coordinating histone modification and chromatin architecture during Pol II transcription. Importantly, we identify a role for FACT in the maintenance of promoter-proximal Pol II pausing, a key step in transcription activation in higher eukaryotes. These findings bring to light a broader role for FACT in the regulation of Pol II transcription.

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.

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.

Discovering Transcriptional Regulatory Elements From Run-On and Sequencing Data Using the Web-Based dREG Gateway

Transcription is a chromatin mark that can be used effectively to identify the location of active enhancers and promoters, collectively known as transcriptional regulatory elements (TREs). We recently introduced dREG, a tool for the identification of TREs using run-on and sequencing (RO-seq) assays, including global run-on and sequencing (GRO-seq), precision run-on and sequencing (PRO-seq), and chromatin run-on and sequencing (ChRO-seq). In this protocol, we present step-by-step instructions for running dREG on an arbitrary run-on and sequencing dataset. Users provide dREG with bigWig files (in which each read is represented by a single base) representing the location of RNA polymerase in a cell or tissue sample of interest, and dREG returns a list of genomic regions that are predicted to be active TREs. Finally, we demonstrate the use of dREG regions in discovering transcription factors controlling response to a stimulus and predicting their target genes. Together, this protocol provides detailed instructions for running dREG on arbitrary run-on and sequencing data. © 2018 by John Wiley & Sons, Inc.

R1 retrotransposons in the nucleolar organizers of Drosophila melanogaster are transcribed by RNA polymerase I upon heat shock

The ribosomal RNA genes (rDNA) of Drosophila melanogaster reside within centromere-proximal nucleolar organizers on both the X and Y chromosomes. Each locus contains between 200-300 tandem repeat rDNA units that encode 18S, 5.8S, 2S, and 28S ribosomal RNAs (rRNAs) necessary for ribosome biogenesis. In arthropods like Drosophila, about 60% of the rDNA genes have R1 and/or R2 retrotransposons inserted at specific sites within their 28S regions; these units likely fail to produce functional 28S rRNA. We showed earlier that R2 expression increases upon nucleolar stress caused by the loss of the ribosome assembly factor, Nucleolar Phosphoprotein of 140 kDa (Nopp140). Here we show that R1 expression is selectively induced by heat shock. Actinomycin D, but not α-amanitin, blocked R1 expression in S2 cells upon heat shock, indicating that R1 elements are transcribed by Pol I. A series of RT-PCRs established read-through transcription by Pol I from the 28S gene region into R1. Sequencing the RT-PCR products confirmed the 28S-R1 RNA junction and the expression of R1 elements within nucleolar rDNA rather than R1 elements known to reside in centromeric heterochromatin. Using a genome-wide precision run-on sequencing (PRO-seq) data set available at NCBI-GEO, we show that Pol I activity on R1 elements is negligible under normal non-heat shock conditions but increases upon heat shock. We propose that prior to heat shock Pol I pauses within the 5' end of R1 where we find a consensus "pause button", and that heat shock releases Pol I for read-through transcription farther into R1.

Use of conditioned media is critical for studies of regulation in response to rapid heat shock

Heat shock response (HSR) maintains and restores protein homeostasis when cells are exposed to proteotoxic heat stress. Heat shock (HS) triggers a rapid and robust change in genome-wide transcription, protein synthesis, and chaperone activity; and therefore, the HSR has been widely used as a model system in these studies. The conventional method of performing instantaneous HS in the laboratory uses heated fresh media to induce HSR when added to cells. However, addition of fresh media to cells may evoke additional cellular responses and signaling pathways. Here, we compared the change in global transcription profile when HS is performed with either heated fresh media or heated conditioned media. We found that the use of heated fresh media induces transcription of hundreds of genes that HS alone does not induce, and masks or partially masks HS-mediated downregulation of thousands of genes. The fresh-media-dependent upregulated genes encode ribosomal subunit proteins involved in translation and RNA processing factors. More importantly, fresh media also induce transcription of several heat shock protein genes (Hsps) in a heat shock factor 1 (HSF1)-independent manner. Thus, we conclude that a conventional method of HS with heated fresh media causes changes in transcription regulation that confound the actual change caused solely by elevated temperature of cells.