Halfpipe: a tool for analyzing metabolic labeling RNA-seq data to quantify RNA half-lives

We introduce Halfpipe, a tool for analyzing RNA-seq data from metabolic RNA labeling experiments. Its main features are the absolute quantification of 4-thiouridine-labeling-induced T>C conversions in the data as generated by SLAM-seq, calculating the proportion of newly synthesized transcripts, and estimating subcellular RNA half-lives. Halfpipe excels at correcting critical biases caused by typically low labeling efficiency. We measure and compare the RNA metabolism in the G1 phase and during the mitosis of synchronized human cells. We find that RNA half-lives of constantly expressed RNAs are similar in mitosis and G1 phase, suggesting that RNA stability of those genes is constant throughout the cell cycle. Our estimates correlate well with literature values and with known RNA sequence features. Halfpipe is freely available at https://github.com/IMSBCompBio/Halfpipe.

From computational models of the splicing code to regulatory mechanisms and therapeutic implications

Since the discovery of RNA splicing and its role in gene expression, researchers have sought a set of rules, an algorithm or a computational model that could predict the splice isoforms, and their frequencies, produced from any transcribed gene in a specific cellular context. Over the past 30 years, these models have evolved from simple position weight matrices to deep-learning models capable of integrating sequence data across vast genomic distances. Most recently, new model architectures are moving the field closer to context-specific alternative splicing predictions, and advances in sequencing technologies are expanding the type of data that can be used to inform and interpret such models. Together, these developments are driving improved understanding of splicing regulatory mechanisms and emerging applications of the splicing code to the rational design of RNA- and splicing-based therapeutics.

Pathogen-derived glyoxylate inhibits Tet2 DNA dioxygenase to facilitate bacterial persister formation

Pathogenic bacterial persistence enables survival during antibiotic treatment, leading to treatment failure and recurrent infections, yet its underlying mechanisms remain unclear. Here, we reveal that glyoxylate, a metabolite originally evolved for alternative carbon utilization, functions as a signaling molecule to reprogram the host transcriptome and promote persister formation. Glyoxylate inhibits the DNA dioxygenase TET2, suppressing pro-inflammatory gene expression and attenuating host immune defense. Notably, stimulating TET2 activity with vitamin C or blocking glyoxylate production by Salmonella reduces bacterial antibiotic resistance and improves infection treatment outcomes. Beyond its metabolic role, glyoxylate emerges as a regulator of host-pathogen interactions, while TET2 plays a critical role in preventing bacterial persistence. Our findings suggest that targeting glyoxylate production or enhancing TET2 activity offers promising therapeutic strategies to combat bacterial persistence and enhance the efficacy of antibiotic treatments.

Investigating Gene Expression Noise Reduction by MicroRNAs and MiRISC Reinforcement by Self-Feedback Regulation of mRNA Degradation

MicroRNA (miRNA) induced silencing complex (miRISC) is the targeting apparatus and arguably rate-limiting step of the miRNA-mediated regulatory subsystem - the major noise reducing though metabolically wasteful mechanism. Recently, we reported that miRISC channels miRNA-mediated regulatory activity back onto their own mRNAs to form negative self-feedback loops, a noise-reduction technique in engineering and synthetic/systems biology. Here, we describe mathematical modeling that predicts mRNA expression noise to correlate negatively with degradation rate (K deg ) and noise reduction by self-feedback control ofK deg . We also calculatedK deg and expression noise of mRNAs detected in a cutting-edge total-RNA single-cell RNA-seq (scRNA-seq) dataset. As predicted, miRNA-targeted mRNAs exhibited higherK deg values in conjunction with lower inter-cell expression noise. Moreover, as predicted by our self-feedback loop model, miRISC mRNAs (AGO1/2/3 and TNRC6A/B/C) exhibited further reduced expression noise. In short, mathematical-modeling and total-RNA scRNA-seq data-analysis shed insight into operational trade-off between noise reduction and metabolic/energetic expenditure in producing miRNA-targeted mRNAs destined for enhanced degradation and translational inhibition, as well as negative self-feedback loop reinforcement of miRISC - the core of miRNA-mediated noise-reduction subsystem. To our knowledge, this is the first report of concurrent mRNA degradation and expression noise analyses and of noise reduction by self-feedback control of mRNA degradation.

Promoter-proximal RNA polymerase II termination regulates transcription during human cell type transition

Metazoan gene transcription by RNA polymerase II (Pol II) is regulated in the promoter-proximal region. Pol II can undergo termination in the promoter-proximal region but whether this can contribute to transcription regulation in cells remains unclear. Here we extend our previous multiomics analysis to quantify changes in transcription kinetics during a human cell type transition event. We observe that upregulation of transcription involves an increase in initiation frequency and, at a set of genes, a decrease in promoter-proximal termination. In turn, downregulation of transcription involves a decrease in initiation frequency and an increase in promoter-proximal termination. Thus, promoter-proximal termination of Pol II contributes to the regulation of human gene transcription.

SLAM-seq reveals independent contributions of RNA processing and stability to gene expression in African trypanosomes

Gene expression is a multi-step process that converts DNA-encoded information into proteins, involving RNA transcription, maturation, degradation, and translation. While transcriptional control is a major regulator of protein levels, the role of post-transcriptional processes such as RNA processing and degradation is less well understood due to the challenge of measuring their contributions individually. To address this challenge, we investigated the control of gene expression in Trypanosoma brucei, a unicellular parasite assumed to lack transcriptional control. Instead, mRNA levels in T. brucei are controlled by post-transcriptional processes, which enabled us to disentangle the contribution of both processes to total mRNA levels. In this study, we developed an efficient metabolic RNA labeling approach and combined ultra-short metabolic labeling with transient transcriptome sequencing (TT-seq) to confirm the long-standing assumption that RNA polymerase II transcription is unregulated in T. brucei. In addition, we established thiol (SH)-linked alkylation for metabolic sequencing of RNA (SLAM-seq) to globally quantify RNA processing rates and half-lives. Our data, combined with scRNA-seq data, indicate that RNA processing and stability independently affect total mRNA levels and contribute to the variability seen between individual cells in African trypanosomes.

LDB1 establishes multi-enhancer networks to regulate gene expression

How specific enhancer-promoter pairing is established remains mostly unclear. Besides the CTCF/cohesin machinery, few nuclear factors have been studied for a direct role in physically connecting regulatory elements. Using a murine erythroid cell model, we show via acute degradation experiments that LDB1 directly and broadly promotes connectivity among regulatory elements. Most LDB1-mediated contacts, even those spanning hundreds of kb, can form in the absence of CTCF, cohesin, or YY1 as determined using multiple degron systems. Moreover, an engineered LDB1-driven chromatin loop is cohesin independent. Cohesin-driven loop extrusion does not stall at LDB1-occupied sites but aids the formation of a subset of LDB1-anchored loops. Leveraging the dynamic reorganization of nuclear architecture during the transition from mitosis to G1 phase, we observe that loop formation and de novo LDB1 occupancy correlate and can occur independently of structural loops. Tri-C and Region Capture Micro-C reveal that LDB1 organizes multi-enhancer networks to activate transcription. These findings establish LDB1 as a driver of spatial connectivity.

Recurrent oncogenic ZC3H18 mutations stabilize endogenous retroviral RNA

Endogenous retroviral (ERV) RNA is highly expressed in cancer, although the molecular causes and consequences remain unknown. We found that ZC3H18 (Z18), a component of multiple nuclear RNA surveillance complexes, has recurrent truncating mutations in cancer. We show that Z18trunc mutations are oncogenic and that Z18 plays an evolutionarily conserved role in nuclear RNA surveillance of ERV RNA. In zebrafish, Z18trunc expedited melanoma onset and promoted a specific accumulation of ERV RNA. Z18 mutant human cell lines from the Cancer Cell Line Encyclopedia also expressed higher levels of ERV RNA. In engineered human melanoma cells, Z18trunc enhanced ERV RNA accumulation more than loss of one Z18 copy, indicating dominant negative activity. Z18trunc directly bound and stabilized ERV RNA. Notably, expression of ERV RNA was sufficient to expedite oncogenesis in a zebrafish model, which is the first evidence of which we are aware that ERV transcripts can play a functional role in cancer. Our work illuminates a mechanism for elevated ERV transcripts in cancer and supports that aberrant RNA accumulation is broadly oncogenic.

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.

DSIF factor Spt5 coordinates transcription, maturation and exoribonucleolysis of RNA polymerase II transcripts

Precursor messenger RNA (pre-mRNA) is processed into its functional form during RNA polymerase II (Pol II) transcription. Although functional coupling between transcription and pre-mRNA processing is established, the underlying mechanisms are not fully understood. We show that the key transcription termination factor, RNA exonuclease Xrn2 engages with Pol II forming a stable complex. Xrn2 activity is stimulated by Spt5 to ensure efficient degradation of nascent RNA leading to Pol II dislodgement from DNA. Our results support a model where Xrn2 first forms a stable complex with the elongating Pol II to achieve its full activity in degrading nascent RNA revising the current 'torpedo' model of termination, which posits that RNA degradation precedes Xrn2 engagement with Pol II. Spt5 is also a key factor that attenuates the expression of non-coding transcripts, coordinates pre-mRNA splicing and 3'-end processing. Our findings indicate that engagement with the transcribing Pol II is an essential regulatory step modulating the activity of RNA enzymes such as Xrn2, thus advancing our understanding of how RNA maturation is controlled during transcription.

Mapping Active RNA Polymerases in Proliferating and Quiescent Fission Yeast Cells Using Precision Run-On Sequencing

The development of next-generation sequencing (NGS) approaches to investigate the functioning of RNA polymerases has led to groundbreaking advances in the field of transcriptional regulation. One powerful method, Precision nuclear Run-On sequencing (PRO-seq), maps the locations of RNA polymerase active sites genome-wide at high resolution. PRO-seq provides a snapshot of strand-specific transcriptional activity and does not rely on immunoprecipitation of the polymerase of interest. Notably, this technique has been utilized to investigate the control of the RNA polymerase II transcription cycle in a variety of model systems. However, the initially published PRO-seq method required significant amounts of starting sample and was technically challenging, both of which were deterrents for its broader use. Recently, an improved and simplified version called qPRO-seq that reduced the length of the experiment and the quantity of necessary input sample was developed for human and Drosophila cell lines. Here we provide an updated, step-by-step protocol in which we have validated and optimized qPRO-seq for the fission yeast Schizosaccharomyces pombe. Importantly, we have implemented this method for assessing RNA polymerase activity in nutrient-limiting conditions, for both proliferating and nitrogen-depleted quiescent cells.

Mechanisms of RNA Polymerase II Termination at the 3'-End of Genes

RNA polymerase II (RNAPII) is responsible for the synthesis of a diverse set of RNA molecules, including protein-coding messenger RNAs (mRNAs) and many short non-coding RNAs (ncRNAs). For this purpose, RNAPII relies on a multitude of factors that regulate the transcription cycle, from initiation and promoter-proximal pausing, through elongation and finally termination. RNAPII transcription termination at the end of genes ensures the release of RNAPII from the DNA template and its efficient recycling for further rounds of transcription. Termination of RNAPII is tightly coupled to 3'-end mRNA processing, which constitutes an important trigger for the subsequent transcription termination event. In this review, we discuss the current understanding of RNAPII termination mechanisms, focusing on 'canonical' termination at the 3'-end of genes. We also integrate the allosteric and 'torpedo' models into a unified model of termination, and describe the different termination factors that have been identified to date, paying special attention to the human factors and their mechanism of action at the molecular level. Indeed, in recent years the development of novel approaches in structural biology, biochemistry and cell biology have together led to a more detailed comprehension of the different mechanisms of RNAPII termination, and a better understanding of their importance in regulating gene expression, especially under cellular stress and pathological situations.

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.

Fateful Decisions of Where to Cut the Line: Pathology Associated with Aberrant 3' End Processing and Transcription Termination

Aberrant gene expression lies at the heart of many pathologies. This review will point out how 3' end processing, the final mRNA-maturation step in the transcription cycle, is surprisingly prone to regulated as well as stochastic variations with a wide range of consequences. Whereas smaller variations contribute to the plasticity of gene expression, larger alternations to 3' end processing and coupled transcription termination can lead to pathological consequences. These can be caused by the local mutation of one gene or affect larger numbers of genes systematically, if aspects of the mechanisms of 3' end processing and transcription termination are altered.

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.

The phosphatase PP1 sustains global transcription by promoting RNA polymerase II pause release

RNA polymerase II progression from initiation to elongation is driven in part by a cascade of protein kinases acting on the core transcription machinery. Conversely, the corresponding phosphatases, notably PP2A and PP1-the most abundant serine-threonine phosphatases in cells-are thought to mainly impede polymerase progression, respectively restraining pause release at promoters and elongation at terminators. Here, we reveal an unexpected role of PP1, within the phosphatase 1 nuclear targeting subunit (PNUTS)-PP1 complex, in sustaining global transcriptional activation in human cells. Acute disruption of PNUTS-PP1 leads to severe defects in the release of paused polymerase and subsequent downregulation for the majority of transcribed genes. PNUTS-PP1 promotes pause release by dephosphorylating multiple substrates, including the 7SK small nuclear ribonucleoprotein particle (snRNP) subunit MEPCE, a known pausing regulator. PNUTS-PP1 exhibits antagonistic functions compared with Integrator-PP2A (INTAC) phosphatase, which generally inhibits pause release. Our research thus highlights opposing roles of PP1 and PP2A in modulating genome-wide transcriptional pausing and gene expression.

Distinct H3K9me3 heterochromatin maintenance dynamics govern different gene programmes and repeats in pluripotent cells

H3K9me3 heterochromatin, established by lysine methyltransferases (KMTs) and compacted by heterochromatin protein 1 (HP1) isoforms, represses alternative lineage genes and DNA repeats. Our understanding of H3K9me3 heterochromatin stability is presently limited to individual domains and DNA repeats. Here we engineered Suv39h2-knockout mouse embryonic stem cells to degrade remaining two H3K9me3 KMTs within 1 hour and found that both passive dilution and active removal contribute to H3K9me3 decay within 12-24 hours. We discovered four different H3K9me3 decay rates across the genome and chromatin features and transcription factor binding patterns that predict the stability classes. A 'binary switch' governs heterochromatin compaction, with HP1 rapidly dissociating from heterochromatin upon KMT depletion and a particular threshold level of HP1 limiting pioneer factor binding, chromatin opening and exit from pluripotency within 12 h. Unexpectedly, receding H3K9me3 domains unearth residual HP1β peaks enriched with heterochromatin-inducing proteins. Our findings reveal distinct H3K9me3 heterochromatin maintenance dynamics governing gene networks and repeats that together safeguard pluripotency.

PRC2-EZH1 contributes to circadian gene expression by orchestrating chromatin states and RNA polymerase II complex stability

Circadian rhythmicity of gene expression is a conserved feature of cell physiology. This involves fine-tuning between transcriptional and post-transcriptional mechanisms and strongly depends on the metabolic state of the cell. Together these processes guarantee an adaptive plasticity of tissue-specific genetic programs. However, it is unclear how the epigenome and RNA Pol II rhythmicity are integrated. Here we show that the PcG protein EZH1 has a gateway bridging function in postmitotic skeletal muscle cells. On the one hand, the circadian clock master regulator BMAL1 directly controls oscillatory behavior and periodic assembly of core components of the PRC2-EZH1 complex. On the other hand, EZH1 is essential for circadian gene expression at alternate Zeitgeber times, through stabilization of RNA Polymerase II preinitiation complexes, thereby controlling nascent transcription. Collectively, our data show that PRC2-EZH1 regulates circadian transcription both negatively and positively by modulating chromatin states and basal transcription complex stability.

RNA kinetics influence the response to transcriptional perturbation in leukaemia cell lines

Therapeutic targeting of dysregulated transcription has emerged as a promising strategy for the treatment of cancers, such as leukaemias. The therapeutic response to small molecule inhibitors of Bromodomain-Containing Proteins (BRD), such as BRD2 and BRD4, P300/cAMP-response element binding protein (CBP) and Cyclin Dependent Kinases (CDKs), is generally attributed to the selective disruption of oncogenic gene expression driven by enhancers, super-enhancers (SEs) and lineage-specific transcription factors (TFs), including the c-MYC oncogene. The selectivity of compounds targeting the transcriptional machinery may be further shaped by post-transcriptional processes. To quantitatively assess the contribution of post-transcriptional regulation in responses to transcription inhibition, we performed multi-omics analyses to accurately measure mRNA production and decay kinetics. We demonstrate that it is not only the selective disruption of mRNA production, but rather mRNA decay rates that largely influence the selectivity associated with transcriptional inhibition. Accordingly, genes down-regulated with transcriptional inhibitors are largely characterized by extremely rapid mRNA production and turnover. In line with this notion, stabilization of the c-MYC transcript through swapping of its 3' untranslated region (UTR) rendered c-MYC insensitive to transcriptional targeting. This failed to negate the impact on c-MYC downstream targets and did not abrogate therapeutic responses. Finally, we provide evidence that modulating post-transcriptional pathways, such as through ELAVL1 targeting, can sensitize long-lived mRNAs to transcriptional inhibition and be considered as a combination therapy approach in leukaemia. Taken together, these data demonstrate that mRNA kinetics influence the therapeutic response to transcriptional perturbation and can be modulated for novel therapeutic outcomes using transcriptional agents in leukaemia.

Uncovering the functions and mechanisms of regulatory elements-associated non-coding RNAs

Over the past decade, regulatory non-coding RNAs (ncRNAs) produced by RNA Pol II have been revealed as meaningful players in various essential cellular functions. In particular, thousands of ncRNAs are produced at transcriptional regulatory elements such as enhancers and promoters, where they may exert multiple functions to regulate proper development, cellular programming, transcription or genomic stability. Here, we review the mechanisms involving these regulatory element-associated ncRNAs, and particularly enhancer RNAs (eRNAs) and PROMoter uPstream Transcripts (PROMPTs). We contextualize the mechanisms described to the processing and degradation of these short lived RNAs. We summarize recent findings explaining how ncRNAs operate locally at promoters and enhancers, or further away, either shortly after their production by RNA Pol II, or through post-transcriptional stabilization. Such discoveries lead to a converging model accounting for how ncRNAs influence cellular fate, by acting on transcription and chromatin structure, which may further involve factors participating to 3D nuclear organization.

Transcriptome and Temporal Transcriptome Analyses in Single Cells

Transcriptome analysis in single cells, enabled by single-cell RNA sequencing, has become a prevalent approach in biomedical research, ranging from investigations of gene regulation to the characterization of tissue organization. Over the past decade, advances in single-cell RNA sequencing technology, including its underlying chemistry, have significantly enhanced its performance, marking notable improvements in methodology. A recent development in the field, which integrates RNA metabolic labeling with single-cell RNA sequencing, has enabled the profiling of temporal transcriptomes in individual cells, offering new insights into dynamic biological processes involving RNA kinetics and cell fate determination. In this review, we explore the chemical principles and design improvements that have enhanced single-molecule capture efficiency, improved RNA quantification accuracy, and increased cellular throughput in single-cell transcriptome analysis. We also illustrate the concept of RNA metabolic labeling for detecting newly synthesized transcripts and summarize recent advancements that enable single-cell temporal transcriptome analysis. Additionally, we examine data analysis strategies for the precise quantification of newly synthesized transcripts and highlight key applications of transcriptome and temporal transcriptome analyses in single cells.

Notch induces transcription by stimulating release of paused RNA polymerase II

Notch proteins undergo ligand-induced proteolysis to release a nuclear effector that influences a wide range of cellular processes by regulating transcription. Despite years of study, however, how Notch induces the transcription of its target genes remains unclear. Here, we comprehensively examine the response to human Notch1 across a time course of activation using high-resolution genomic assays of chromatin accessibility and nascent RNA production. Our data reveal that Notch induces target gene transcription primarily by releasing paused RNA polymerase II (RNAPII). Moreover, in contrast to prevailing models suggesting that Notch acts by promoting chromatin accessibility, we found that open chromatin was established at Notch-responsive regulatory elements prior to Notch signal induction through SWI/SNF-mediated remodeling. Together, these studies show that the nuclear response to Notch signaling is dictated by the pre-existing chromatin state and RNAPII distribution at the time of signal activation.

Human CCR4-NOT globally regulates gene expression and is a novel silencer of retrotransposon activation

CCR4-NOT regulates multiple steps in gene regulation and has been well studied in budding yeast, but much less is known about the human complex. Auxin-induced degradation was used to rapidly deplete the scaffold subunit CNOT1, and CNOT4, to characterize the functions of human CCR4-NOT in gene regulation. Depleting CNOT1 increased RNA levels and caused a widespread decrease in RNA decay. In contrast, CNOT4 depletion only modestly changed steady-state RNA levels and, surprisingly, led to a global acceleration in mRNA decay. Further, depleting either subunit resulted in a global increase in RNA synthesis. In contrast to most of the genome, the transcription of KRAB-Zinc-Finger-protein (KZNFs) genes, especially those on chromosome 19, was repressed. KZNFs are transcriptional repressors of retrotransposable elements (rTEs), and consistent with the decreased KZNFs expression, rTEs, mainly Long Interspersed Nuclear Elements (LINEs), were activated. These data establish CCR4-NOT as a global regulator of gene expression and a novel silencer of rTEs.

H3.3K122A results in a neomorphic phenotype in mouse embryonic stem cells

Canonical histone H3 and histone variant H3.3 are posttranslationally modified with the genomic distribution of these marks denoting different features and these modifications may influence transcription. While the majority of posttranslational modifications occur on histone tails, there are defined modifications within the globular domain, such as acetylation of H3K122/H3.3K122. To understand the function of the amino acid H3.3K122 in transcriptional regulation, we attempted to generate H3.3K122A mouse embryonic stem (mES) cells but were unsuccessful. Through multi-omic profiling of mutant cell lines harboring two or three of four H3.3 targeted alleles, we have uncovered that H3.3K122A is neomorphic and results in lethality. This is surprising as prior studies demonstrate H3.3-null mES cells are viable and pluripotent but exhibit a reduced differentiation capacity. Together, these studies have uncovered a novel dependence of a globular domain residue within H3.3 for viability and broadened our understanding of how histone variants contribute to transcription regulation and pluripotency in mES cells.

Expanding and improving analyses of nucleotide recoding RNA-seq experiments with the EZbakR suite

Nucleotide recoding RNA sequencing methods (NR-seq; TimeLapse-seq, SLAM-seq, TUC-seq, etc.) are powerful approaches for assaying transcript population dynamics. In addition, these methods have been extended to probe a host of regulated steps in the RNA life cycle. Current bioinformatic tools significantly constrain analyses of NR-seq data. To address this limitation, we developed EZbakR, an R package to facilitate a more comprehensive set of NR-seq analyses, and fastq2EZbakR, a Snakemake pipeline for flexible preprocessing of NR-seq datasets, collectively referred to as the EZbakR suite. Together, these tools generalize many aspects of the NR-seq analysis workflow. The fastq2EZbakR pipeline can assign reads to a diverse set of genomic features (e.g., genes, exons, splice junctions, etc.), and EZbakR can perform analyses on any combination of these features. EZbakR extends standard NR-seq mutational modeling to support multi-label analyses (e.g., s4U and s6G dual labeling), and implements an improved hierarchical model to better account for transcript-to-transcript variance in metabolic label incorporation. EZbakR also generalizes dynamical systems modeling of NR-seq data to support analyses of premature mRNA processing and flow between subcellular compartments. Finally, EZbakR implements flexible and well-powered comparative analyses of all estimated parameters via design matrix-specified generalized linear modeling. The EZbakR suite will thus allow researchers to make full, effective use of NR-seq data.

Transcriptome-Wide Profiling of Nascent RNA in Neurons with Enriched H3K27ac Signal Elevates eRNA Identification Efficiency

Growing evidence suggests that activity-dependent gene expression is crucial for neuronal plasticity and behavioral experience. Enhancer RNAs (eRNAs), a class of long noncoding RNAs, play a key role in these processes. However, eRNAs are highly dynamic and are often present at lower levels than their corresponding mRNAs, making them difficult to detect using total RNA-seq techniques. Nascent RNA sequencing, which separates nascent RNAs from the steady-state RNA population, has been shown to increase the ability to detect activity-induced eRNAs with a higher signal-to-noise ratio. However, there is a lack of bioinformatic tools or pipelines for detecting eRNAs utilizing nascent RNA-seq and other multiomics data sets. In this study, we addressed this gap by developing a novel bioinformatic framework, e-finder, for finding eRNAs and have made it available to the scientific community. Additionally, we reanalyzed our previous nascent RNA sequencing data and compared them with total RNA-seq data to identify activity-regulated RNAs in neuronal cell populations. Using H3K27 acetylome data, we characterized activity-dependent eRNAs that drive the transcriptional activity of the target genes. Our analysis identified a subset of eRNAs involved in mediating synapse organization, which showed increased activity-dependent transcription after the potassium chloride stimulation. Notably, our data suggest that nascent RNA-seq with an enriched H3K27ac signal exhibits high resolution to identify potential eRNAs in response to membrane depolarization. Our findings uncover the role of the eRNA-mediated gene activation network in neuronal systems, providing new insights into the molecular processes characterizing neurological diseases.

An RNA-centric view of transcription and genome organization

Foundational models of transcriptional regulation involve the assembly of protein complexes at DNA elements associated with specific genes. These assemblies, which can include transcription factors, cofactors, RNA polymerase, and various chromatin regulators, form dynamic spatial compartments that contribute to both gene regulation and local genome architecture. This DNA-protein-centric view has been modified with recent evidence that RNA molecules have important roles to play in gene regulation and genome structure. Here, we discuss evidence that gene regulation by RNA occurs at multiple levels that include assembly of transcriptional complexes and genome compartments, feedback regulation of active genes, silencing of genes, and control of protein kinases. We thus provide an RNA-centric view of transcriptional regulation that must reside alongside the more traditional DNA-protein-centric perspectives on gene regulation and genome architecture.

TRiC/CCT Chaperonin Governs RNA Polymerase II Activity in the Nucleus to Support RNA Homeostasis

The chaperonin TRiC/CCT is a large hetero-oligomeric ringed-structure that is essential in eukaryotes. While present in the nucleus, TRiC/CCT is typically considered to function in the cytosol where it mediates nascent polypeptide folding and the assembly/disassembly of protein complexes. Here, we investigated the nuclear role of TRiC/CCT. Inactivation of TRiC/CCT resulted in a significant increase in the production of nascent RNA leading to the accumulation of noncoding transcripts. The influence on transcription was not due to cytoplasmic TRiC/CCT-activities or other nuclear proteins as the effect was observed when TRiC/CCT was evicted from the nucleus and restricted to the cytoplasm. Rather, our data support a direct role of TRiC/CCT in regulating RNA polymerase II activity, as the chaperonin modulated nascent RNA production both in vivo and in vitro. Overall, our studies reveal a new avenue by which TRiC/CCT contributes to cell homeostasis by regulating the activity of nuclear RNA polymerase II.

A Compendium of G-Flipon Biological Functions That Have Experimental Validation

As with all new fields of discovery, work on the biological role of G-quadruplexes (GQs) has produced a number of results that at first glance are quite baffling, sometimes because they do not fit well together, but mostly because they are different from commonly held expectations. Like other classes of flipons, those that form G-quadruplexes have a repeat sequence motif that enables the fold. The canonical DNA motif (G3N1-7)3G3, where N is any nucleotide and G is guanine, is a feature that is under active selection in avian and mammalian genomes. The involvement of G-flipons in genome maintenance traces back to the invertebrate Caenorhabditis elegans and to ancient DNA repair pathways. The role of GQs in transcription is supported by the observation that yeast Rap1 protein binds both B-DNA, in a sequence-specific manner, and GQs, in a structure-specific manner, through the same helix. Other sequence-specific transcription factors (TFs) also engage both conformations to actuate cellular transactions. Noncoding RNAs can also modulate GQ formation in a sequence-specific manner and engage the same cellular machinery as localized by TFs, linking the ancient RNA world with the modern protein world. The coevolution of noncoding RNAs and sequence-specific proteins is supported by studies of early embryonic development, where the transient formation of G-quadruplexes coordinates the epigenetic specification of cell fate.

Multiple direct and indirect roles of the Paf1 complex in transcription elongation, splicing, and histone modifications

The polymerase-associated factor 1 (Paf1) complex (Paf1C) is a conserved protein complex with critical functions during eukaryotic transcription. Previous studies showed that Paf1C is multi-functional, controlling specific aspects of transcription ranging from RNA polymerase II (RNAPII) processivity to histone modifications. However, it is unclear how specific Paf1C subunits directly impact transcription and coupled processes. We have compared conditional depletion to steady-state deletion for each Paf1C subunit to determine the direct and indirect contributions to gene expression in Saccharomyces cerevisiae. Using nascent transcript sequencing, RNAPII profiling, and modeling of transcription elongation dynamics, we have demonstrated direct effects of Paf1C subunits on RNAPII processivity and elongation rate and indirect effects on transcript splicing and repression of antisense transcripts. Further, our results suggest that the direct transcriptional effects of Paf1C cannot be readily assigned to any particular histone modification. This work comprehensively analyzes both the immediate and the extended roles of each Paf1C subunit in transcription elongation and transcript regulation.

Distinct H3K9me3 heterochromatin maintenance dynamics govern different gene programs and repeats in pluripotent cells

H3K9me3-heterochromatin, established by lysine methyltransferases (KMTs) and compacted by HP1 isoforms, represses alternative lineage genes and DNA repeats. Our understanding of H3K9me3-heterochromatin stability is presently limited to individual domains and DNA repeats. We engineered Suv39h2 KO mouse embryonic stem cells to degrade remaining two H3K9me3-KMTs within one hour and found that both passive dilution and active removal contribute to H3K9me3 decay within 12-24 hours. We discovered four different H3K9me3 decay rates across the genome and chromatin features and transcription factor binding patterns that predict the stability classes. A "binary switch" governs heterochromatin compaction, with HP1 rapidly dissociating from heterochromatin upon KMTs' depletion and a particular threshold level of HP1 limiting pioneer factor binding, chromatin opening, and exit from pluripotency within 12 hr. Unexpectedly, receding H3K9me3 domains unearth residual HP1β peaks enriched with heterochromatin-inducing proteins. Our findings reveal distinct H3K9me3-heterochromatin maintenance dynamics governing gene networks and repeats that together safeguard pluripotency.

Study of the RNA splicing kinetics via in vivo 5-EU labeling

Splicing is an important step of gene expression in all eukaryotes. Splice sites might be used with different efficiency, giving rise to alternative splicing products. At the same time, splice sites might be used at a variable rate. We used 5-ethynyl uridine labeling to sequence a nascent transcriptome of HeLa cells and deduced the rate of splicing for each donor and acceptor splice site. The following correlation analysis showed a correspondence of primary transcript features with the rate of splicing. Some dependencies we revealed were anticipated, such as a splicing rate decrease with a decreased complementarity of the donor splice site to U1 and acceptor sites to U2 snRNAs. Other dependencies were more surprising, like a negative influence of a distance to the 5' end on the rate of the acceptor splicing site utilization, or the differences in splicing rate between long, short, and RBM17-dependent introns. We also observed a deceleration of last intron splicing with an increase of the distance to the poly(A) site, which might be explained by the cooperativity of the splicing and polyadenylation. Additional analysis of splicing kinetics of SF3B4 knockdown cells suggested the impairment of a U2 snRNA recognition step. As a result, we deconvoluted the effects of several examined features on the splicing rate into a single regression model. The data obtained here are useful for further studies in the field, as they provide general splicing rate dependencies as well as help to justify the existence of slowly removed splice sites.

The signaling cascade of induction and maintenance of ES cell diapause

Nutrient deficiency during pregnancy in numerous animal species can induce the state of embryonic diapause. Diapause is characterized by changes in protein and gene expression that minimize the organism's reliance on external energy sources and ensure survival. Remarkably, the systematic changes associated with diapause appear to spare the gene expression program that supports embryonic cells' maintenance in the pluripotent state. The phenomenon of the differentiation "freeze" during diapause can be reproduced in vitro. Mimicking nutrient deficiency by pharmacological inhibition of mTOR induces the diapause-like state in ES cells without affecting ES cell pluripotency. We discovered a connection between mTOR signaling and the chromatin-bound bromodomain and extra-terminal (BET) transcriptional regulator BRD4, showing a key role of BET-protein in the induction of diapause-like state in ES cells. mTOR inhibition rapidly and negatively impacts BRD4 binding to chromatin, which is associated with changes in gene expression that can contribute to diapause. Conversely, pharmacological inhibition of BET-protein circumvents the diapause dependence on mTOR inhibition and causes the diapause-like state. BET-repressed diapause-like ES cells retain the undifferentiated pluripotent state, which is associated with upregulation of a functionally linked group of genes encoding negative regulators of MAP kinase (MAPK) signaling and inactivation of MAP kinase. The transcriptional switch-off of MAP kinase following chronic BET inhibition imitates the transcriptional de-repression of MAP kinase negative regulators in response to mTOR inhibition. Mechanistically, suppression of mTOR or BET-protein leads to a profound decline in Capicua transcriptional repressor (CIC) at promoters of key negative regulators of MAP kinase. The discovered mTOR-BRD4 axis in the induction of diapause and the rapid transcriptional shut-off of differentiation program is likely to play a major role in the maintenance of embryonic diapause in vivo, as well as in controlling of the undifferentiated state of various types of stem cells during diapause-like metabolic dormancy.

Deep learning and direct sequencing of labeled RNA captures transcriptome dynamics

In eukaryotes, genes produce a variety of distinct RNA isoforms, each with potentially unique protein products, coding potential or regulatory signals such as poly(A) tail and nucleotide modifications. Assessing the kinetics of RNA isoform metabolism, such as transcription and decay rates, is essential for unraveling gene regulation. However, it is currently impeded by lack of methods that can differentiate between individual isoforms. Here, we introduce RNAkinet, a deep convolutional and recurrent neural network, to detect nascent RNA molecules following metabolic labeling with the nucleoside analog 5-ethynyl uridine and long-read, direct RNA sequencing with nanopores. RNAkinet processes electrical signals from nanopore sequencing directly and distinguishes nascent from pre-existing RNA molecules. Our results show that RNAkinet prediction performance generalizes in various cell types and organisms and can be used to quantify RNA isoform half-lives. RNAkinet is expected to enable the identification of the kinetic parameters of RNA isoforms and to facilitate studies of RNA metabolism and the regulatory elements that influence it.

YY1-controlled regulatory connectivity and transcription are influenced by the cell cycle

Few transcription factors have been examined for their direct roles in physically connecting enhancers and promoters. Here acute degradation of Yin Yang 1 (YY1) in erythroid cells revealed its requirement for the maintenance of numerous enhancer-promoter loops, but not compartments or domains. Despite its reported ability to interact with cohesin, the formation of YY1-dependent enhancer-promoter loops does not involve stalling of cohesin-mediated loop extrusion. Integrating mitosis-to-G1-phase dynamics, we observed partial retention of YY1 on mitotic chromatin, predominantly at gene promoters, followed by rapid rebinding during mitotic exit, coinciding with enhancer-promoter loop establishment. YY1 degradation during the mitosis-to-G1-phase interval revealed a set of enhancer-promoter loops that require YY1 for establishment during G1-phase entry but not for maintenance in interphase, suggesting that cell cycle stage influences YY1's architectural function. Thus, as revealed here for YY1, chromatin architectural functions of transcription factors can vary in their interplay with CTCF and cohesin as well as by cell cycle stage.

H3.3K122A results in a neomorphic phenotype in mouse embryonic stem cells

The histone variant H3.3 acts in coordination with histone posttranslational modifications and other chromatin features to facilitate appropriate transcription. Canonical histone H3 and histone variant H3.3 are post-translationally modified with the genomic distribution of these marks denoting different features and with more recent evidence suggesting that these modifications may influence transcription. While the majority of posttranslational modifications occur on histone tails, there are defined modifications within the globular domain, such as acetylation of H3K122/H3.3K122. To understand the function of the residue H3.3K122 in transcriptional regulation, we attempted to generate H3.3K122A mouse embryonic stem (mES) cells but were unsuccessful. Through multi-omic profiling of mutant cell lines harboring two or three of four H3.3 targeted alleles, we have uncovered that H3.3K122A is neomorphic and results in lethality. This is surprising as prior studies demonstrate H3.3-null mES cells are viable and pluripotent, albeit with reduced differentiation capacity. Together, these studies have uncovered a novel dependence of a globular domain residue of H3.3 for viability and broadened our understanding of how histone variants contribute to transcription regulation and pluripotency in mES cells.

LDB1 establishes multi-enhancer networks to regulate gene expression

How specific enhancer-promoter pairing is established is still mostly unclear. Besides the CTCF/cohesin machinery, only a few nuclear factors have been studied for a direct role in physically connecting regulatory elements. Here, we show via acute degradation experiments that LDB1 directly and broadly promotes enhancer-promoter loops. Most LDB1-mediated contacts, even those spanning hundreds of kb, can form in the absence of CTCF, cohesin, or YY1 as determined via the use of multiple degron systems. Moreover, an engineered LDB1-driven chromatin loop is cohesin independent. Cohesin-driven loop extrusion does not stall at LDB1 occupied sites but may aid the formation of a subset of LDB1 anchored loops. Leveraging the dynamic reorganization of nuclear architecture during the transition from mitosis to G1-phase, we establish a relationship between LDB1-dependent interactions in the context of TAD organization and gene activation. Lastly, Tri-C and Region Capture Micro-C reveal that LDB1 organizes multi-enhancer networks to activate transcription. This establishes LDB1 as a direct driver of regulatory network inter-connectivity.

A convenient single-cell newly synthesized transcriptome assay reveals gene expression dynamics during early-stage T-cell activation

Sequencing newly synthesized transcriptome alongside regular transcriptome in single cells enables the study of gene expression temporal dynamics during rapid chromatin and gene regulation processes. However, existing assays to profile single-cell newly synthesized transcriptome require in-house technical expertise to achieve high cellular throughput, limiting their widespread application. Here, we developed NOTE-seq, a method that simultaneously profiles regular and newly synthesized transcriptomes in single cells. NOTE-seq integrates 4-thiouridine labeling of newly synthesized RNA, thiol-alkylation-based chemical conversion, and a streamlined workflow on the 10X Genomics platform, offering high cellular throughput that is accessible and convenient for regular biology laboratories without specialized single-cell expertise. Using NOTE-seq, we characterized the temporal dynamics of gene expression during early-stage T-cell activation in Jurkat and naïve T cells, identified transcription factors and regulons, and discovered Fli-1 as a master transcription factor for gene regulation upon T-cell stimulation. Interestingly, chemotherapeutic topoisomerase inhibitor affects Fli-1 level in T cells, indicating potential complications for the immune system.

The Pseudo-Natural Product Tafbromin Selectively Targets the TAF1 Bromodomain 2

Phenotypic assays detect small-molecule bioactivity at functionally relevant cellular sites, and inherently cover a variety of targets and mechanisms of action. They can uncover new small molecule-target pairs and may give rise to novel biological insights. By means of an osteoblast differentiation assay which employs a Hedgehog (Hh) signaling agonist as stimulus and which monitors an endogenous marker for osteoblasts, we identified a pyrrolo[3,4-g]quinoline (PQ) pseudo-natural product (PNP) class of osteogenesis inhibitors. The most potent PQ, termed Tafbromin, impairs canonical Hh signaling and modulates osteoblast differentiation through binding to the bromodomain 2 of the TATA-box binding protein-associated factor 1 (TAF1). Tafbromin is the most selective TAF1 bromodomain 2 ligand and promises to be an invaluable tool for the study of biological processes mediated by TAF1(2) bromodomains.

Dissecting key regulators of transcriptome kinetics through scalable single-cell RNA profiling of pooled CRISPR screens

We present a combinatorial indexing method, PerturbSci-Kinetics, for capturing whole transcriptomes, nascent transcriptomes and single guide RNA (sgRNA) identities across hundreds of genetic perturbations at the single-cell level. Profiling a pooled CRISPR screen targeting various biological processes, we show the gene expression regulation during RNA synthesis, processing and degradation, miRNA biogenesis and mitochondrial mRNA processing, systematically decoding the genome-wide regulatory network that underlies RNA temporal dynamics at scale.

Transcription directionality is licensed by Integrator at active human promoters

A universal characteristic of eukaryotic transcription is that the promoter recruits RNA polymerase II (RNAPII) to produce both precursor mRNAs (pre-mRNAs) and short unstable promoter upstream transcripts (PROMPTs) toward the opposite direction. However, how the transcription machinery selects the correct direction to produce pre-mRNAs is largely unknown. Here, through multiple acute auxin-inducible degradation systems, we show that rapid depletion of an RNAPII-binding protein complex, Integrator, results in robust PROMPT accumulation throughout the genome. Interestingly, the accumulation of PROMPTs is compensated by the reduction of pre-mRNA transcripts in actively transcribed genes. Consistently, Integrator depletion alters the distribution of polymerase between the sense and antisense directions, which is marked by increased RNAPII-carboxy-terminal domain Tyr1 phosphorylation at PROMPT regions and a reduced Ser2 phosphorylation level at transcription start sites. Mechanistically, the endonuclease activity of Integrator is critical to suppress PROMPT production. Furthermore, our data indicate that the presence of U1 binding sites on nascent transcripts could counteract the cleavage activity of Integrator. In this process, the absence of robust U1 signal at most PROMPTs allows Integrator to suppress the antisense transcription and shift the transcriptional balance in favor of the sense direction.

Genome-wide quantification of RNA flow across subcellular compartments reveals determinants of the mammalian transcript life cycle

Dissecting the regulatory mechanisms controlling mammalian transcripts from production to degradation requires quantitative measurements of mRNA flow across the cell. We developed subcellular TimeLapse-seq to measure the rates at which RNAs are released from chromatin, exported from the nucleus, loaded onto polysomes, and degraded within the nucleus and cytoplasm in human and mouse cells. These rates varied substantially, yet transcripts from genes with related functions or targeted by the same transcription factors and RNA-binding proteins flowed across subcellular compartments with similar kinetics. Verifying these associations uncovered a link between DDX3X and nuclear export. For hundreds of RNA metabolism genes, most transcripts with retained introns were degraded by the nuclear exosome, while the remaining molecules were exported with stable cytoplasmic lifespans. Transcripts residing on chromatin for longer had extended poly(A) tails, whereas the reverse was observed for cytoplasmic mRNAs. Finally, machine learning identified molecular features that predicted the diverse life cycles of mRNAs.

KAP1 negatively regulates RNA polymerase II elongation kinetics to activate signal-induced transcription

Signal-induced transcriptional programs regulate critical biological processes through the precise spatiotemporal activation of Immediate Early Genes (IEGs); however, the mechanisms of transcription induction remain poorly understood. By combining an acute depletion system with several genomics approaches to interrogate synchronized, temporal transcription, we reveal that KAP1/TRIM28 is a first responder that fulfills the temporal and heightened transcriptional demand of IEGs. Acute KAP1 loss triggers an increase in RNA polymerase II elongation kinetics during early stimulation time points. This elongation defect derails the normal progression through the transcriptional cycle during late stimulation time points, ultimately leading to decreased recruitment of the transcription apparatus for re-initiation thereby dampening IEGs transcriptional output. Collectively, KAP1 plays a counterintuitive role by negatively regulating transcription elongation to support full activation across multiple transcription cycles of genes critical for cell physiology and organismal functions.

Notch induces transcription by stimulating release of paused RNA Polymerase II

Notch proteins undergo ligand-induced proteolysis to release a nuclear effector that influences a wide range of cellular processes by regulating transcription. Despite years of study, however, how Notch induces the transcription of its target genes remains unclear. Here, we comprehensively examined the response to human Notch1 across a time course of activation using high-resolution genomic assays of chromatin accessibility and nascent RNA production. Our data reveal that Notch induces target gene transcription primarily by releasing paused RNA polymerase II (RNAPII). Moreover, in contrast to prevailing models suggesting that Notch acts by promoting chromatin accessibility, we found that open chromatin was established at Notch-responsive regulatory elements prior to Notch signal induction, through SWI/SNF-mediated remodeling. Together, these studies show that the nuclear response to Notch signaling is dictated by the pre-existing chromatin state and RNAPII distribution at the time of signal activation.

Direct and selective pharmacological disruption of the YAP-TEAD interface by IAG933 inhibits Hippo-dependent and RAS-MAPK-altered cancers

The YAP-TEAD protein-protein interaction mediates YAP oncogenic functions downstream of the Hippo pathway. To date, available YAP-TEAD pharmacologic agents bind into the lipid pocket of TEAD, targeting the interaction indirectly via allosteric changes. However, the consequences of a direct pharmacological disruption of the interface between YAP and TEADs remain largely unexplored. Here, we present IAG933 and its analogs as potent first-in-class and selective disruptors of the YAP-TEAD protein-protein interaction with suitable properties to enter clinical trials. Pharmacologic abrogation of the interaction with all four TEAD paralogs resulted in YAP eviction from chromatin and reduced Hippo-mediated transcription and induction of cell death. In vivo, deep tumor regression was observed in Hippo-driven mesothelioma xenografts at tolerated doses in animal models as well as in Hippo-altered cancer models outside mesothelioma. Importantly this also extended to larger tumor indications, such as lung, pancreatic and colorectal cancer, in combination with RTK, KRAS-mutant selective and MAPK inhibitors, leading to more efficacious and durable responses. Clinical evaluation of IAG933 is underway.

Technology to the rescue: how to uncover the role of transposable elements in preimplantation development

Transposable elements (TEs) are highly expressed in preimplantation development. Preimplantation development is the phase when the cells of the early embryo undergo the first cell fate choice and change from being totipotent to pluripotent. A range of studies have advanced our understanding of TEs in preimplantation, as well as their epigenetic regulation and functional roles. However, many questions remain about the implications of TE expression during early development. Challenges originate first due to the abundance of TEs in the genome, and second because of the limited cell numbers in preimplantation. Here we review the most recent technological advancements promising to shed light onto the role of TEs in preimplantation development. We explore novel avenues to identify genomic TE insertions and improve our understanding of the regulatory mechanisms and roles of TEs and their RNA and protein products during early development.

CDK7 kinase activity promotes RNA polymerase II promoter escape by facilitating initiation factor release

Cyclin-dependent kinase 7 (CDK7), part of the general transcription factor TFIIH, promotes gene transcription by phosphorylating the C-terminal domain of RNA polymerase II (RNA Pol II). Here, we combine rapid CDK7 kinase inhibition with multi-omics analysis to unravel the direct functions of CDK7 in human cells. CDK7 inhibition causes RNA Pol II retention at promoters, leading to decreased RNA Pol II initiation and immediate global downregulation of transcript synthesis. Elongation, termination, and recruitment of co-transcriptional factors are not directly affected. Although RNA Pol II, initiation factors, and Mediator accumulate at promoters, RNA Pol II complexes can also proceed into gene bodies without promoter-proximal pausing while retaining initiation factors and Mediator. Further downstream, RNA Pol II phosphorylation increases and initiation factors and Mediator are released, allowing recruitment of elongation factors and an increase in RNA Pol II elongation velocity. Collectively, CDK7 kinase activity promotes the release of initiation factors and Mediator from RNA Pol II, facilitating RNA Pol II escape from the promoter.

De novo assembly of transcriptomes and differential gene expression analysis using short-read data from emerging model organisms - a brief guide

Many questions in biology benefit greatly from the use of a variety of model systems. High-throughput sequencing methods have been a triumph in the democratization of diverse model systems. They allow for the economical sequencing of an entire genome or transcriptome of interest, and with technical variations can even provide insight into genome organization and the expression and regulation of genes. The analysis and biological interpretation of such large datasets can present significant challenges that depend on the 'scientific status' of the model system. While high-quality genome and transcriptome references are readily available for well-established model systems, the establishment of such references for an emerging model system often requires extensive resources such as finances, expertise and computation capabilities. The de novo assembly of a transcriptome represents an excellent entry point for genetic and molecular studies in emerging model systems as it can efficiently assess gene content while also serving as a reference for differential gene expression studies. However, the process of de novo transcriptome assembly is non-trivial, and as a rule must be empirically optimized for every dataset. For the researcher working with an emerging model system, and with little to no experience with assembling and quantifying short-read data from the Illumina platform, these processes can be daunting. In this guide we outline the major challenges faced when establishing a reference transcriptome de novo and we provide advice on how to approach such an endeavor. We describe the major experimental and bioinformatic steps, provide some broad recommendations and cautions for the newcomer to de novo transcriptome assembly and differential gene expression analyses. Moreover, we provide an initial selection of tools that can assist in the journey from raw short-read data to assembled transcriptome and lists of differentially expressed genes.

FACT maintains chromatin architecture and thereby stimulates RNA polymerase II pausing during transcription in vivo

Facilitates chromatin transcription (FACT) is a histone chaperone that supports transcription through chromatin in vitro, but its functional roles in vivo remain unclear. Here, we analyze the in vivo functions of FACT with the use of multi-omics analysis after rapid FACT depletion from human cells. We show that FACT depletion destabilizes chromatin and leads to transcriptional defects, including defective promoter-proximal pausing and elongation, and increased premature termination of RNA polymerase II. Unexpectedly, our analysis revealed that promoter-proximal pausing depends not only on the negative elongation factor (NELF) but also on the +1 nucleosome, which is maintained by FACT.

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.