Genome-wide profiling of RNA polymerase transcription at nucleotide resolution in human cells with native elongating transcript sequencing

Chemo-phosphoproteomic profiling with ATR inhibitors berzosertib and gartisertib uncovers new biomarkers and DNA damage response regulators

The ATR kinase protects cells against DNA damage and replication stress and represents a promising anti-cancer drug target. The ATR inhibitors (ATRi) berzosertib and gartisertib are both in clinical trials for the treatment of advanced solid tumours as monotherapy or in combination with genotoxic agents. We carried out quantitative phospho-proteomic screening for ATR biomarkers that are highly sensitive to berzosertib and gartisertib, using an optimized mass spectrometry pipeline. Screening identified a range of novel ATR-dependent phosphorylation events, which were grouped into three broad classes: i) targets whose phosphorylation is highly sensitive to ATRi and which could be the next generation of ATR biomarkers; ii) proteins with known genome maintenance roles not previously known to be regulated by ATR; iii) novel targets whose cellular roles are unclear. Class iii targets represent candidate DNA damage response proteins and, with this in mind, proteins in this class were subjected to secondary screening for recruitment to DNA damage sites. We show that one of the proteins recruited, SCAF1, interacts with RNAPII in a phospho-dependent manner and recruitment requires PARP activity and interaction with RNAPII. We also show that SCAF1 deficiency partly rescues RAD51 loading in cells lacking the BRCA1 tumour suppressor. Taken together these data reveal potential new ATR biomarkers and new genome maintenance factors.

Predicting RNA polymerase II transcriptional elongation pausing and associated histone code

RNA Polymerase II (Pol II) transcriptional elongation pausing is an integral part of the dynamic regulation of gene transcription in the genome of metazoans. It plays a pivotal role in many vital biological processes and disease progression. However, experimentally measuring genome-wide Pol II pausing is technically challenging and the precise governing mechanism underlying this process is not fully understood. Here, we develop RP3 (RNA Polymerase II Pausing Prediction), a network regularized logistic regression machine learning method, to predict Pol II pausing events by integrating genome sequence, histone modification, gene expression, chromatin accessibility, and protein-protein interaction data. RP3 can accurately predict Pol II pausing in diverse cellular contexts and unveil the transcription factors that are associated with the Pol II pausing machinery. Furthermore, we utilize a forward feature selection framework to systematically identify the combination of histone modification signals associated with Pol II pausing. RP3 is freely available at https://github.com/AMSSwanglab/RP3.

A kinetic dichotomy between mitochondrial and nuclear gene expression processes

Oxidative phosphorylation (OXPHOS) complexes, encoded by both mitochondrial and nuclear DNA, are essential producers of cellular ATP, but how nuclear and mitochondrial gene expression steps are coordinated to achieve balanced OXPHOS subunit biogenesis remains unresolved. Here, we present a parallel quantitative analysis of the human nuclear and mitochondrial messenger RNA (mt-mRNA) life cycles, including transcript production, processing, ribosome association, and degradation. The kinetic rates of nearly every stage of gene expression differed starkly across compartments. Compared with nuclear mRNAs, mt-mRNAs were produced 1,100-fold more, degraded 7-fold faster, and accumulated to 160-fold higher levels. Quantitative modeling and depletion of mitochondrial factors LRPPRC and FASTKD5 identified critical points of mitochondrial regulatory control, revealing that the mitonuclear expression disparities intrinsically arise from the highly polycistronic nature of human mitochondrial pre-mRNA. We propose that resolving these differences requires a 100-fold slower mitochondrial translation rate, illuminating the mitoribosome as a nexus of mitonuclear co-regulation.

Molecular robotic agents that survey molecular landscapes for information retrieval

DNA-based artificial motors have allowed the recapitulation of biological functions and the creation of new features. Here, we present a molecular robotic system that surveys molecular environments and reports spatial information in an autonomous and repeated manner. A group of molecular agents, termed 'crawlers', roam around and copy information from DNA-labeled targets, generating records that reflect their trajectories. Based on a mechanism that allows random crawling, we show that our system is capable of counting the number of subunits in example molecular complexes. Our system can also detect multivalent proximities by generating concatenated records from multiple local interactions. We demonstrate this capability by distinguishing colocalization patterns of three proteins inside fixed cells under different conditions. These mechanisms for examining molecular landscapes may serve as a basis towards creating large-scale detailed molecular interaction maps inside the cell with nanoscale resolution.

Splicing quality control mediated by DHX15 and its G-patch activator SUGP1

Pre-mRNA splicing is surveilled at different stages by quality control (QC) mechanisms. The leukemia-associated DExH-box family helicase hDHX15/scPrp43 is known to disassemble spliceosomes after splicing. Here, using rapid protein depletion and analysis of nascent and mature RNA to enrich for direct effects, we identify a widespread splicing QC function for DHX15 in human cells, consistent with recent in vitro studies. We find that suboptimal introns with weak splice sites, multiple branch points, and cryptic introns are repressed by DHX15, suggesting a general role in promoting splicing fidelity. We identify SUGP1 as a G-patch factor that activates DHX15's splicing QC function. This interaction is dependent on both DHX15's ATPase activity and on SUGP1's U2AF ligand motif (ULM) domain. Together, our results support a model in which DHX15 plays a major role in splicing QC when recruited and activated by SUGP1.

RNA polymerase II pausing temporally coordinates cell cycle progression and erythroid differentiation

Controlled release of promoter-proximal paused RNA polymerase II (RNA Pol II) is crucial for gene regulation. However, studying RNA Pol II pausing is challenging, as pause-release factors are almost all essential. In this study, we identified heterozygous loss-of-function mutations in SUPT5H, which encodes SPT5, in individuals with β-thalassemia. During erythropoiesis in healthy human cells, cell cycle genes were highly paused as cells transition from progenitors to precursors. When the pathogenic mutations were recapitulated by SUPT5H editing, RNA Pol II pause release was globally disrupted, and as cells began transitioning from progenitors to precursors, differentiation was delayed, accompanied by a transient lag in erythroid-specific gene expression and cell cycle kinetics. Despite this delay, cells terminally differentiate, and cell cycle phase distributions normalize. Therefore, hindering pause release perturbs proliferation and differentiation dynamics at a key transition during erythropoiesis, identifying a role for RNA Pol II pausing in temporally coordinating the cell cycle and erythroid differentiation.

A kinetic dichotomy between mitochondrial and nuclear gene expression drives OXPHOS biogenesis

Oxidative phosphorylation (OXPHOS) complexes, encoded by both mitochondrial and nuclear DNA, are essential producers of cellular ATP, but how nuclear and mitochondrial gene expression steps are coordinated to achieve balanced OXPHOS biogenesis remains unresolved. Here, we present a parallel quantitative analysis of the human nuclear and mitochondrial messenger RNA (mt-mRNA) life cycles, including transcript production, processing, ribosome association, and degradation. The kinetic rates of nearly every stage of gene expression differed starkly across compartments. Compared to nuclear mRNAs, mt-mRNAs were produced 700-fold higher, degraded 5-fold faster, and accumulated to 170-fold higher levels. Quantitative modeling and depletion of mitochondrial factors, LRPPRC and FASTKD5, identified critical points of mitochondrial regulatory control, revealing that the mitonuclear expression disparities intrinsically arise from the highly polycistronic nature of human mitochondrial pre-mRNA. We propose that resolving these differences requires a 100-fold slower mitochondrial translation rate, illuminating the mitoribosome as a nexus of mitonuclear co-regulation.

High-sensitive nascent transcript sequencing reveals BRD4-specific control of widespread enhancer and target gene transcription

Gene transcription by RNA polymerase II (Pol II) is under control of promoters and distal regulatory elements known as enhancers. Enhancers are themselves transcribed by Pol II correlating with their activity. How enhancer transcription is regulated and coordinated with transcription at target genes has remained unclear. Here, we developed a high-sensitive native elongating transcript sequencing approach, called HiS-NET-seq, to provide an extended high-resolution view on transcription, especially at lowly transcribed regions such as enhancers. HiS-NET-seq uncovers new transcribed enhancers in human cells. A multi-omics analysis shows that genome-wide enhancer transcription depends on the BET family protein BRD4. Specifically, BRD4 co-localizes to enhancer and promoter-proximal gene regions, and is required for elongation activation at enhancers and their genes. BRD4 keeps a set of enhancers and genes in proximity through long-range contacts. From these studies BRD4 emerges as a general regulator of enhancer transcription that may link transcription at enhancers and genes.

7SK methylation by METTL3 promotes transcriptional activity

A fundamental feature of cell signaling is the conversion of extracellular signals into adaptive transcriptional responses. The role of RNA modifications in this process is poorly understood. The small nuclear RNA 7SK prevents transcriptional elongation by sequestering the cyclin dependent kinase 9/cyclin T1 (CDK9/CCNT1) positive transcription elongation factor (P-TEFb) complex. We found that epidermal growth factor signaling induces phosphorylation of the enzyme methyltransferase 3 (METTL3), leading to METTL3-mediated methylation of 7SK. 7SK methylation enhanced its binding to heterogeneous nuclear ribonucleoproteins, causing the release of the HEXIM1 P-TEFb complex subunit1 (HEXIM1)/P-TEFb complex and inducing transcriptional elongation. Our findings establish the mechanism underlying 7SK activation and uncover a previously unknown function for the m6A modification in converting growth factor signaling events into a regulatory transcriptional response via an RNA methylation-dependent switch.

RNA Polymerase II pausing temporally coordinates cell cycle progression and erythroid differentiation

The controlled release of promoter-proximal paused RNA polymerase II (Pol II) into productive elongation is a major step in gene regulation. However, functional analysis of Pol II pausing is difficult because factors that regulate pause release are almost all essential. In this study, we identified heterozygous loss-of-function mutations in SUPT5H , which encodes SPT5, in individuals with β-thalassemia unlinked to HBB mutations. During erythropoiesis in healthy human cells, cell cycle genes were highly paused at the transition from progenitors to precursors. When the pathogenic mutations were recapitulated by SUPT5H editing, Pol II pause release was globally disrupted, and the transition from progenitors to precursors was delayed, marked by a transient lag in erythroid-specific gene expression and cell cycle kinetics. Despite this delay, cells terminally differentiate, and cell cycle phase distributions normalize. Therefore, hindering pause release perturbs proliferation and differentiation dynamics at a key transition during erythropoiesis, revealing a role for Pol II pausing in the temporal coordination between the cell cycle and differentiation.

Aging is associated with increased chromatin accessibility and reduced polymerase pausing in liver

Regulation of gene expression is linked to the organization of the genome. With age, chromatin alterations occur on all levels of genome organization, accompanied by changes in the gene expression profile. However, little is known about the changes in the level of transcriptional regulation. Here, we used a multi‐omics approach and integrated ATAC‐, RNA‐ and NET‐seq to identify age‐related changes in the chromatin landscape of murine liver and to investigate how these are linked to transcriptional regulation. We provide the first systematic inventory of the connection between aging, chromatin accessibility, and transcriptional regulation in a whole tissue. Aging in murine liver is characterized by an increase in chromatin accessibility at promoter regions, but not in an increase in transcriptional output. Instead, aging is accompanied by a decrease in promoter‐proximal pausing of RNA polymerase II (Pol II), while initiation of transcription is not decreased as assessed by RNA polymerase mapping using CUT&RUN. Based on the data reported, we propose that these age‐related changes in transcriptional regulation are due to a reduced stability of the pausing complex.

PAMP-triggered genetic reprogramming involves widespread alternative transcription initiation and an immediate transcription factor wave

Immune responses triggered by pathogen-associated molecular patterns (PAMPs) are key to pathogen defense, but drivers and stabilizers of the growth-to-defense genetic reprogramming remain incompletely understood in plants. Here, we report a time-course study of the establishment of PAMP-triggered immunity (PTI) using cap analysis of gene expression. We show that around 15% of all transcription start sites (TSSs) rapidly induced during PTI define alternative transcription initiation events. From these, we identify clear examples of regulatory TSS change via alternative inclusion of target peptides or domains in encoded proteins, or of upstream open reading frames in mRNA leader sequences. We also find that 60% of PAMP response genes respond earlier than previously thought. In particular, a cluster of rapidly and transiently PAMP-induced genes is enriched in transcription factors (TFs) whose functions, previously associated with biological processes as diverse as abiotic stress adaptation and stem cell activity, appear to converge on growth restriction. Furthermore, examples of known potentiators of PTI, in one case under direct mitogen-activated protein kinase control, support the notion that the rapidly induced TFs could constitute direct links to PTI signaling pathways and drive gene expression changes underlying establishment of the immune state.

EU-RNA-seq for in vivo labeling and high throughput sequencing of nascent transcripts

The protocol allows for labeling nascent RNA without isolating nuclei. The cell-permeable uridine analog, 5-ethynyluridine (EU), is added to media to allow in vivo labeling of nascent transcripts. Cells are lysed, total RNA is collected, and biotin is conjugated to EU-labeled RNAs. Custom biotin RNAs are added and biotinylated RNAs are isolated for generation of cDNA libraries. The sequencing data are normalized to controls for quantitative assessment of the nascent transcriptome. The protocol takes 4 days, not including sequencing and analysis.

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

•

Labeling nascent RNAs with EU does not require nuclear isolation

•

EU-RNA-Seq maps nascent transcripts across a large dynamic range

•

Control biotin RNAs are customizable

•

Custom biotin RNAs allow for global normalization following sequencing

Loss of Gene Information: Discrepancies between RNA Sequencing, cDNA Microarray, and qRT-PCR

Molecular analyses of normal and diseased cells give insight into changes in gene expression and help in understanding the background of pathophysiological processes. Years after cDNA microarrays were established in research, RNA sequencing (RNA-seq) became a key method of quantitatively measuring the transcriptome. In this study, we compared the detection of genes by each of the transcriptome analysis methods: cDNA array, quantitative RT-PCR, and RNA-seq. As expected, we found differences in the gene expression profiles of the aforementioned techniques. Here, we present selected genes that exemplarily demonstrate the observed differences and calculations to reveal that a strong RNA secondary structure, as well as sample preparation, can affect RNA-seq. In summary, this study addresses an important issue with a strong impact on gene expression analysis in general. Therefore, we suggest that these findings need to be considered when dealing with data from transcriptome analyses.

A BRD4-mediated elongation control point primes transcribing RNA polymerase II for 3'-processing and termination

Transcription elongation has emerged as a regulatory hub in gene expression of metazoans. A major control point occurs during early elongation before RNA polymerase II (Pol II) is released into productive elongation. Prior research has linked BRD4 with transcription elongation. Here, we use rapid BET protein and BRD4-selective degradation along with quantitative genome-wide approaches to investigate direct functions of BRD4 in Pol II transcription regulation. Notably, as an immediate consequence of acute BRD4 loss, promoter-proximal pause release is impaired, and transcriptionally engaged Pol II past this checkpoint undergoes readthrough transcription. An integrated proteome-wide analysis uncovers elongation and 3'-RNA processing factors as core BRD4 interactors. BRD4 ablation disrupts the recruitment of general 3'-RNA processing factors at the 5'-control region, which correlates with RNA cleavage and termination defects. These studies, performed in human cells, reveal a BRD4-mediated checkpoint and begin to establish a molecular link between 5'-elongation control and 3'-RNA processing.

RNAnue: efficient data analysis for RNA-RNA interactomics

RNA–RNA inter- and intramolecular interactions are fundamental for numerous biological processes. While there are reasonable approaches to map RNA secondary structures genome-wide, understanding how different RNAs interact to carry out their regulatory functions requires mapping of intermolecular base pairs. Recently, different strategies to detect RNA–RNA duplexes in living cells, so called direct duplex detection (DDD) methods, have been developed. Common to all is the Psoralen-mediated in vivo RNA crosslinking followed by RNA Proximity Ligation to join the two interacting RNA strands. Sequencing of the RNA via classical RNA-seq and subsequent specialised bioinformatic analyses the result in the prediction of inter- and intramolecular RNA–RNA interactions. Existing approaches adapt standard RNA-seq analysis pipelines, but often neglect inherent features of RNA–RNA interactions that are useful for filtering and statistical assessment. Here we present RNAnue, a general pipeline for the inference of RNA–RNA interactions from DDD experiments that takes into account hybridisation potential and statistical significance to improve prediction accuracy. We applied RNAnue to data from different DDD studies and compared our results to those of the original methods. This showed that RNAnue performs better in terms of quantity and quality of predictions.

A heat shock-responsive lncRNA Heat acts as a HSF1-directed transcriptional brake via m6A modification

Long noncoding RNAs (lncRNAs) are key regulators of gene expression in diverse cellular contexts and biological processes. Given the surprising range of shapes and sizes, how distinct lncRNAs achieve functional specificity remains incompletely understood. Here, we identified a heat shock-inducible lncRNA, Heat, in mouse cells that acts as a transcriptional brake to restrain stress gene expression. Functional characterization reveals that Heat directly binds to heat shock transcription factor 1 (HSF1), thereby targeting stress genes in a trans-acting manner. Intriguingly, Heat is heavily methylated in the form of m⁶A. Although dispensable for HSF1 binding, Heat methylation is required for silencing stress genes to attenuate heat shock response. Consistently, m⁶A depletion results in prolonged activation of stress genes. Furthermore, Heat mediates these effects via the nuclear m⁶A reader YTHDC1, forming a transcriptional silencing complex for stress genes. Our study reveals a crucial role of nuclear epitranscriptome in the transcriptional regulation of heat shock response.

Conserved DNA sequence features underlie pervasive RNA polymerase pausing

Pausing of transcribing RNA polymerase is regulated and creates opportunities to control gene expression. Research in metazoans has so far mainly focused on RNA polymerase II (Pol II) promoter-proximal pausing leaving the pervasive nature of pausing and its regulatory potential in mammalian cells unclear. Here, we developed a pause detecting algorithm (PDA) for nucleotide-resolution occupancy data and a new native elongating transcript sequencing approach, termed nested NET-seq, that strongly reduces artifactual peaks commonly misinterpreted as pausing sites. Leveraging PDA and nested NET-seq reveal widespread genome-wide Pol II pausing at single-nucleotide resolution in human cells. Notably, the majority of Pol II pauses occur outside of promoter-proximal gene regions primarily along the gene-body of transcribed genes. Sequence analysis combined with machine learning modeling reveals DNA sequence properties underlying widespread transcriptional pausing including a new pause motif. Interestingly, key sequence determinants of RNA polymerase pausing are conserved between human cells and bacteria. These studies indicate pervasive sequence-induced transcriptional pausing in human cells and the knowledge of exact pause locations implies potential functional roles in gene expression.

Expression of the miR-302/367 microRNA cluster is regulated by a conserved long non-coding host-gene

MicroRNAs are important regulators of cellular functions. MiR-302/367 is a polycistronic miRNA cluster that can induce and maintain pluripotency. Here we investigate the transcriptional control and the processing of the miR-302 host-gene in mice. Our results indicate that the mmu-miR-302 host-gene is alternatively spliced, polyadenylated and exported from the nucleus. The regulatory sequences extend at least 2 kb upstream of the transcription start site and contain several conserved binding sites for both transcriptional activators and repressors. The gene structure and regulatory elements are highly conserved between mouse and human. So far, regulating miR-302 expression is the only known function of the miR-302 host-gene. Even though we here only provide one example, regulation of microRNA transcription might be a so far little recognized function of long non-coding RNA genes.

A dual role for H2A.Z.1 in modulating the dynamics of RNA polymerase II initiation and elongation

RNA polymerase II (RNAPII) pausing immediately downstream of the transcription start site is a critical rate-limiting step for the expression of most metazoan genes. During pause release, RNAPII encounters a highly conserved +1 H2A.Z nucleosome, yet how this histone variant contributes to transcription is poorly understood. Here, using an inducible protein degron system combined with genomic approaches and live cell super-resolution microscopy, we show that H2A.Z.1 modulates RNAPII dynamics across most genes in murine embryonic stem cells. Our quantitative analysis shows that H2A.Z.1 slows the rate of RNAPII pause release and consequently impacts negative elongation factor dynamics as well as nascent transcription. Consequently, H2A.Z.1 also impacts re-loading of the pre-initiation complex components TFIIB and TBP. Altogether, this work provides a critical mechanistic link between H2A.Z.1 and the proper induction of mammalian gene expression programs through the regulation of RNAPII dynamics and pause release.

Variability in mRNA translation: a random matrix theory approach

The rate of mRNA translation depends on the initiation, elongation, and termination rates of ribosomes along the mRNA. These rates depend on many "local" factors like the abundance of free ribosomes and tRNA molecules in the vicinity of the mRNA molecule. All these factors are stochastic and their experimental measurements are also noisy. An important question is how protein production in the cell is affected by this considerable variability. We develop a new theoretical framework for addressing this question by modeling the rates as identically and independently distributed random variables and using tools from random matrix theory to analyze the steady-state production rate. The analysis reveals a principle of universality: the average protein production rate depends only on the of the set of possible values that the random variable may attain. This explains how total protein production can be stabilized despite the overwhelming stochasticticity underlying cellular processes.

Revealing nascent RNA processing dynamics with nano-COP

During maturation, eukaryotic precursor RNAs undergo processing events including intron splicing, 3'-end cleavage, and polyadenylation. Here we describe nanopore analysis of co-transcriptional processing (nano-COP), a method for probing the timing and patterns of RNA processing. An extension of native elongating transcript sequencing, which quantifies transcription genome-wide through short-read sequencing of nascent RNA 3' ends, nano-COP uses long-read nascent RNA sequencing to observe global patterns of RNA processing. First, nascent RNA is stringently purified through a combination of 4-thiouridine metabolic labeling and cellular fractionation. In contrast to cDNA or short-read-based approaches relying on reverse transcription or amplification, the sample is sequenced directly through nanopores to reveal the native context of nascent RNA. nano-COP identifies both active transcription sites and splice isoforms of single RNA molecules during synthesis, providing insight into patterns of intron removal and the physical coupling between transcription and splicing. The nano-COP protocol yields data within 3 d.

Molecular motor traffic with a slow binding site

We discuss how the presence of a slow binding site in molecular motor traffic gives rise to defect-induced "traffic jams" that have properties different from those of the well-studied boundary-induced jams that originate from an imbalance between initiation and termination. To this end we analyze in detail the stationary distribution of a lattice gas model for traffic of molecular motors with a defect. In particular, we obtain analytically the exact spatial distribution of motors, the probability distribution of the random position of the molecular traffic jam and we report unexpected spatial anticorrelations between local molecular motor densities near the defect.

New tools for 'ZEBRA-FISHING'

Zebrafish has been established as a classical model for developmental studies, yet in the past years, with the explosion of novel technological methods, the use of zebrafish as a model has expanded. One of the prominent fields that took advantage of zebrafish as a model organism early on is hematopoiesis, the process of blood cell generation from hematopoietic stem and progenitor cells (HSPCs). In zebrafish, HSPCs are born early during development in the aorta-gonad-mesonephros region and then translocate to the caudal hematopoietic tissue, where they expand and finally take residence in the kidney marrow. This journey is tightly regulated at multiple levels from extracellular signals to chromatin. In order to delineate the mechanistic underpinnings of this process, next-generation sequencing techniques could be an important ally. Here, we describe genome-wide approaches that have been undertaken to delineate zebrafish hematopoiesis.

Divergence in alternative polyadenylation contributes to gene regulatory differences between humans and chimpanzees

While comparative functional genomic studies have shown that inter-species differences in gene expression can be explained by corresponding inter-species differences in genetic and epigenetic regulatory mechanisms, co-transcriptional mechanisms, such as alternative polyadenylation (APA), have received little attention. We characterized APA in lymphoblastoid cell lines from six humans and six chimpanzees by identifying and estimating the usage for 44,432 polyadenylation sites (PAS) in 9518 genes. Although APA is largely conserved, 1705 genes showed significantly different PAS usage (FDR 0.05) between species. Genes with divergent APA also tend to be differentially expressed, are enriched among genes showing differences in protein translation, and can explain a subset of observed inter-species protein expression differences that do not differ at the transcript level. Finally, we found that genes with a dominant PAS, which is used more often than other PAS, are particularly enriched for differentially expressed genes.

Persistence of RNA transcription during DNA replication delays duplication of transcription start sites until G2/M

As transcription and replication use DNA as substrate, conflicts between transcription and replication can occur, leading to genome instability with direct consequences for human health. To determine how the two processes are coordinated throughout S phase, we characterize both processes together at high resolution. We find that transcription occurs during DNA replication, with transcription start sites (TSSs) not fully replicated along with surrounding regions and remaining under-replicated until late in the cell cycle. TSSs undergo completion of DNA replication specifically when cells enter mitosis, when RNA polymerase II is removed. Intriguingly, G2/M DNA synthesis occurs at high frequency in unperturbed cell culture, but it is not associated with increased DNA damage and is fundamentally separated from mitotic DNA synthesis. TSSs duplicated in G2/M are characterized by a series of specific features, including high levels of antisense transcription, making them difficult to duplicate during S phase.

A machine learning-based framework for modeling transcription elongation

RNA polymerase II (Pol II) generally pauses at certain positions along gene bodies, thereby interrupting the transcription elongation process, which is often coupled with various important biological functions, such as precursor mRNA splicing and gene expression regulation. Characterizing the transcriptional elongation dynamics can thus help us understand many essential biological processes in eukaryotic cells. However, experimentally measuring Pol II elongation rates is generally time and resource consuming. We developed PEPMAN (polymerase II elongation pausing modeling through attention-based deep neural network), a deep learning-based model that accurately predicts Pol II pausing sites based on the native elongating transcript sequencing (NET-seq) data. Through fully taking advantage of the attention mechanism, PEPMAN is able to decipher important sequence features underlying Pol II pausing. More importantly, we demonstrated that the analyses of the PEPMAN-predicted results around various types of alternative splicing sites can provide useful clues into understanding the cotranscriptional splicing events. In addition, associating the PEPMAN prediction results with different epigenetic features can help reveal important factors related to the transcription elongation process. All these results demonstrated that PEPMAN can provide a useful and effective tool for modeling transcription elongation and understanding the related biological factors from available high-throughput sequencing data.

Rapid and Scalable Profiling of Nascent RNA with fastGRO

Genome-wide profiling of nascent RNA has become a fundamental tool to study transcription regulation. Unlike steady-state RNA-sequencing (RNA-seq), nascent RNA profiling mirrors real-time activity of RNA polymerases and provides an accurate readout of transcriptome-wide variations. Some species of nuclear RNAs (i.e., large intergenic noncoding RNAs [lincRNAs] and eRNAs) have a short half-life and can only be accurately gauged by nascent RNA techniques. Furthermore, nascent RNA-seq detects post-cleavage RNA at termination sites and promoter-associated antisense RNAs, providing insights into RNA polymerase II (RNAPII) dynamics and processivity. Here, we present a run-on assay with 4-thio ribonucleotide (4-S-UTP) labeling, followed by reversible biotinylation and affinity purification via streptavidin. Our protocol allows streamlined sample preparation within less than 3 days. We named the technique fastGRO (fast Global Run-On). We show that fastGRO is highly reproducible and yields a more complete and extensive coverage of nascent RNA than comparable techniques can. Importantly, we demonstrate that fastGRO is scalable and can be performed with as few as 0.5 × 10⁶ cells.

Barbieri et al. developed fastGRO, a nascent RNA-sequencing technique based on nuclear run-on. Using a streamlined, under-3-days protocol, fastGRO tracks the activity of RNA polymerase for differential gene expression analysis, polymerase kinetic studies, and profiling of lowly expressed and unstable RNA species. A low-input fastGRO protocol profiles nascent RNA in as little as 0.5 × 10⁶ cells.

Alternative polyadenylation mediates genetic regulation of gene expression

Little is known about co-transcriptional or post-transcriptional regulatory mechanisms linking noncoding variation to variation in organismal traits. To begin addressing this gap, we used 3’ Seq to study the impact of genetic variation on alternative polyadenylation (APA) in the nuclear and total mRNA fractions of 52 HapMap Yoruba human lymphoblastoid cell lines. We mapped 602 APA quantitative trait loci (apaQTLs) at 10% FDR, of which 152 were nuclear specific. Effect sizes at intronic apaQTLs are negatively correlated with eQTL effect sizes. These observations suggest genetic variants can decrease mRNA expression levels by increasing usage of intronic PAS. We also identified 24 apaQTLs associated with protein levels, but not mRNA expression. Finally, we found that 19% of apaQTLs can be associated with disease. Thus, our work demonstrates that APA links genetic variation to variation in gene expression, protein expression, and disease risk, and reveals uncharted modes of genetic regulation.

Metabolic Labeling of RNAs Uncovers Hidden Features and Dynamics of the Arabidopsis Transcriptome

Transcriptome analysis by RNA sequencing (RNA-seq) has become an indispensable research tool in modern plant biology. Virtually all RNA-seq studies provide a snapshot of the steady state transcriptome, which contains valuable information about RNA populations at a given time but lacks information about the dynamics of RNA synthesis and degradation. Only a few specialized sequencing techniques, such as global run-on sequencing, have been used to provide information about RNA synthesis rates in plants. Here, we demonstrate that RNA labeling with the modified, nontoxic uridine analog 5-ethynyl uridine (5-EU) in Arabidopsis (Arabidopsis thaliana) seedlings provides insight into plant transcriptome dynamics. Pulse labeling with 5-EU revealed nascent and unstable RNAs, RNA processing intermediates generated by splicing, and chloroplast RNAs. Pulse-chase experiments with 5-EU allowed us to determine RNA stabilities without the need for chemical transcription inhibitors such as actinomycin and cordycepin. Inhibitor-free, genome-wide analysis of polyadenylated RNA stability via 5-EU pulse-chase experiments revealed RNAs with shorter half-lives than those reported after chemical inhibition of transcription. In summary, our results indicate that the Arabidopsis nascent transcriptome contains unstable RNAs and RNA processing intermediates and suggest that polyadenylated RNAs have low stability in plants. Our technique lays the foundation for easy, affordable, nascent transcriptome analysis and inhibitor-free analysis of RNA stability in plants.

Splicing Kinetics and Coordination Revealed by Direct Nascent RNA Sequencing through Nanopores

Understanding how splicing events are coordinated across numerous introns in metazoan RNA transcripts requires quantitative analyses of transient RNA processing events in living cells. We developed nanopore analysis of co-transcriptional processing (nano-COP), in which nascent RNAs are directly sequenced through nanopores, exposing the dynamics and patterns of RNA splicing without biases introduced by amplification. Long nano-COP reads reveal that, in human and Drosophila cells, splicing occurs after RNA polymerase II transcribes several kilobases of pre-mRNA, suggesting that metazoan splicing transpires distally from the transcription machinery. Inhibition of the branch-site recognition complex SF3B rapidly diminished global co-transcriptional splicing. We found that splicing order does not strictly follow the order of transcription and is associated with cis-acting elements, alternative splicing, and RNA-binding factors. Further, neighboring introns in human cells tend to be spliced concurrently, implying that splicing of these introns occurs cooperatively. Thus, nano-COP unveils the organizational complexity of RNA processing.

Antisense-Mediated Transcript Knockdown Triggers Premature Transcription Termination

Antisense oligonucleotides (ASOs) that trigger RNase-H-mediated cleavage are commonly used to knock down transcripts for experimental or therapeutic purposes. In particular, ASOs are frequently used to functionally interrogate long noncoding RNAs (lncRNAs) and discriminate lncRNA loci that produce functional RNAs from those whose activity is attributable to the act of transcription. Transcription termination is triggered by cleavage of nascent transcripts, generally during polyadenylation, resulting in degradation of the residual RNA polymerase II (Pol II)-associated RNA by XRN2 and dissociation of elongating Pol II. Here, we show that ASOs act upon nascent transcripts and, consequently, induce premature transcription termination downstream of the cleavage site in an XRN2-dependent manner. Targeting the transcript 3' end with ASOs, however, allows transcript knockdown while preserving Pol II association with the gene body. These results demonstrate that the effects of ASOs on transcription must be considered for appropriate experimental and therapeutic use of these reagents.

Transcription factor binding at Ig enhancers is linked to somatic hypermutation targeting

Secondary diversification of the Ig repertoire occurs through somatic hypermutation (SHM), gene conversion (GCV), and class switch recombination (CSR)-three processes that are initiated by activation-induced cytidine deaminase (AID). AID targets Ig genes at orders of magnitude higher than the rest of the genome, but the basis for this specificity is poorly understood. We have previously demonstrated that enhancers and enhancer-like sequences from Ig genes are capable of stimulating SHM of neighboring genes in a capacity distinct from their roles in increasing transcription. Here, we use an in vitro proteomics approach to identify E-box, MEF2, Ets, and Ikaros transcription factor family members as potential binders of these enhancers. ChIP assays in the hypermutating Ramos B cell line confirmed that many of these factors bound the endogenous Igλ enhancer and/or the IgH intronic enhancer (Eμ) in vivo. Further investigation using SHM reporter assays identified binding sites for E2A and MEF2B in Eμ and demonstrated an association between loss of factor binding and decreases in the SHM stimulating activity of Eμ mutants. Our results provide novel insights into trans-acting factors that dictate SHM targeting and link their activity to specific DNA binding sites within Ig enhancers.

NET-CAGE characterizes the dynamics and topology of human transcribed cis-regulatory elements

Promoters and enhancers are key cis-regulatory elements, but how they operate to generate cell type-specific transcriptomes is not fully understood. We developed a simple and robust method, native elongating transcript-cap analysis of gene expression (NET-CAGE), to sensitively detect 5' ends of nascent RNAs in diverse cells and tissues, including unstable transcripts such as enhancer-derived RNAs. We studied RNA synthesis and degradation at the transcription start site level, characterizing the impact of differential promoter usage on transcript stability. We quantified transcription from cis-regulatory elements without the influence of RNA turnover, and show that enhancer-promoter pairs are generally activated simultaneously on stimulation. By integrating NET-CAGE data with chromatin interaction maps, we show that cis-regulatory elements are topologically connected according to their cell type specificity. We identified new enhancers with high sensitivity, and delineated primary locations of transcription within super-enhancers. Our NET-CAGE dataset derived from human and mouse cells expands the FANTOM5 atlas of transcribed enhancers, with broad applicability to biomedical research.

Nascent RNA and the Coordination of Splicing with Transcription

At each active protein-encoding gene, nascent RNA is tethered to the DNA axis by elongating RNA polymerase II (Pol II) and is continuously altered by splicing and other processing events during its synthesis. This review discusses the development of three major methods that enable us to track the conversion of precursor messenger RNA (pre-mRNA) to messenger RNA (mRNA) products in vivo: live-cell imaging, metabolic labeling of RNA, and RNA-seq of purified nascent RNA. These approaches are complementary, addressing distinct issues of transcription rates and intron lifetimes alongside spatial information regarding the gene position of Pol II at which spliceosomes act. The findings will be placed in the context of active transcription units, each of which-because of the presence of nascent RNA, Pol II, and features of the chromatin environment-will recruit a potentially gene-specific constellation of RNA binding proteins and processing machineries.

Methods for the analysis of transcriptome dynamics

The transcriptome is the complete set of transcripts in a cell or tissue and includes ribosomal RNA (rRNA), messenger RNA (mRNA), transfer RNA (tRNA), and regulatory noncoding RNA. At steady-state, the transcriptome results from a compensatory variation of the transcription and decay rate to maintain the RNA concentration constant. RNA transcription constitutes the first stage in gene expression, and thus is a major and primary mode of gene expression control. Nevertheless, regulation of RNA decay is also a key factor in gene expression control, involving either selective RNA stabilization or enhanced degradation. Transcriptome analysis allows the identification of gene expression alterations, providing new insights regarding the pathways and mechanisms involved in physiological and pathological processes. Upon perturbation of cell homeostasis, rapid changes in gene expression are required to adapt to new conditions. Thus, to better understand the regulatory mechanisms associated with gene expression alterations, it is vital to acknowledge the relative contribution of RNA synthesis and decay to the transcriptome. To the toxicology field, the study of gene expression regulation mechanisms can help identify the early and mechanistic relevant cellular events associated with a particular response. This review aims to provide a critical comparison of the available methods used to analyze the contribution of RNA transcription and decay to gene expression dynamics. Notwithstanding, an integration of the data obtained is necessary to understand the entire repercussions of gene transcription changes at a system-level. Thus, a brief overview of the methods available for the integration and analysis of the data obtained from transcriptome analysis will also be provided.

NET-prism enables RNA polymerase-dedicated transcriptional interrogation at nucleotide resolution

The advent of quantitative approaches that enable interrogation of transcription at single nucleotide resolution has allowed a novel understanding of transcriptional regulation previously undefined. However, little is known, at such high resolution, how transcription factors directly influence RNA Pol II pausing and directionality. To map the impact of transcription/elongation factors on transcription dynamics genome-wide at base pair resolution, we developed an adapted NET-seq protocol called NET-prism (Native Elongating Transcription by Polymerase-Regulated Immunoprecipitants in the Mammalian genome). Application of NET-prism on elongation factors (Spt6, Ssrp1), splicing factors (Sf1), and components of the pre-initiation complex (PIC) (TFIID, and Mediator) reveals their inherent command on transcription dynamics, with regards to directionality and pausing over promoters, splice sites, and enhancers/super-enhancers. NET-prism will be broadly applicable as it exposes transcription factor/Pol II dependent topographic specificity and thus, a new degree of regulatory complexity during gene expression.

An Overview of Methodologies in Studying lncRNAs in the High-Throughput Era: When Acronyms ATTACK!

The discovery of pervasive transcription in eukaryotic genomes provided one of many surprising (and perhaps most surprising) findings of the genomic era and led to the uncovering of a large number of previously unstudied transcriptional events. This pervasive transcription leads to the production of large numbers of noncoding RNAs (ncRNAs) and thus opened the window to study these diverse, abundant transcripts of unclear relevance and unknown function. Since that discovery, recent advances in high-throughput sequencing technologies have identified a large collection of ncRNAs, from microRNAs to long noncoding RNAs (lncRNAs). Subsequent discoveries have shown that many lncRNAs play important roles in various eukaryotic processes; these discoveries have profoundly altered our understanding of the regulation of eukaryotic gene expression. Although the identification of ncRNAs has become a standard experimental approach, the functional characterization of these diverse ncRNAs remains a major challenge. In this chapter, we highlight recent progress in the methods to identify lncRNAs and the techniques to study the molecular function of these lncRNAs and the application of these techniques to the study of plant lncRNAs.

Mettl3 Deficiency Sustains Long-Chain Fatty Acid Absorption through Suppressing Traf6-Dependent Inflammation Response

A better understanding of the molecular mechanism of intestinal fatty acid absorption could lead to novel approaches to treatment and prevention of fatty acid-related metabolic diseases. Although it is confirmed that absorption of long-chain fatty acids (LCFAs) decreases during the pathological processes, the genetic basis and molecular mechanisms remain largely unknown. N ⁶-methyladenosine (m⁶A) is the most prevalent internal modification on eukaryotic mRNA. Recently, m⁶A has been found to play important roles in inflammation and antiviral responses. In this study, we show that deficiency of Mettl3, the core methyltransferase of m⁶A, exerts antimalabsorption of LCFA activity in vitro through inhibiting the inflammation response mediated by LPS. To substantiate this finding further, we found the levels of triglycerides were also sustained in cells with depleted Mettl3, which were cultured in Transwell to polarize with villus formation to simulate the situation in vivo. Mechanistically, depletion of Mettl3 decreases the m⁶A level of Traf6 mRNA, thereby its transcripts are entrapped in the nucleus, followed by the decreased expression of Traf6, leading to the suppression of NF-κB and MAPK signaling pathway. Thus, the inflammation response was suppressed, resulting in the sustained absorption of LCFA. Moreover, we found that ectopic expression of Traf6 largely abolishes the sustained absorption LCFA in Mettl3 depletion cells. Collectively, silencing Mettl3 could sustain LCFA absorption through blocking the TRAF6-dependent inflammation response. Our work uncovers a critical function of m⁶A methylation and provides insight into critical roles of Mettl3 in LCFA absorption and inflammatory disease.

Transcriptional repression by FACT is linked to regulation of chromatin accessibility at the promoter of ES cells

Depletion of FACT in murine embryonic stem cells show mild changes on the nucleosomal landscape but widespread changes in the transcriptome, faster proliferation, and neuronal differentiation rates.

The conserved and essential histone chaperone, facilitates chromatin transcription (FACT), reorganizes nucleosomes during DNA transcription, replication, and repair and ensures both efficient elongation of RNA Pol II and nucleosome integrity. In mammalian cells, FACT is a heterodimer, consisting of SSRP1 and SUPT16. Here, we show that in contrast to yeast, FACT accumulates at the transcription start site of genes reminiscent of RNA polymerase II profile. Depletion of FACT in mouse embryonic stem cells leads to deregulation of developmental and pro-proliferative genes concomitant with hyper-proliferation of mES cells. Using MNase-seq, Assay for Transposase-Accessible Chromatin sequencing, and nascent elongating transcript sequencing, we show that up-regulation of genes coincides with loss of nucleosomes upstream of the transcription start site and concomitant increase in antisense transcription, indicating that FACT impacts the promoter architecture to regulate the expression of these genes. Finally, we demonstrate a role for FACT in cell fate determination and show that FACT depletion primes embryonic stem cells for the neuronal lineage.

Considerations on Experimental Design and Data Analysis of Chromatin Immunoprecipitation Experiments

Arguably one of the most valuable techniques to study chromatin organization, ChIP is the method of choice to map the contacts established between proteins and genomic DNA. Ever since its inception, more than 30 years ago, ChIP has been constantly evolving, improving, and expanding its capabilities and reach. Despite its widespread use by many laboratories across a wide variety of disciplines, ChIP assays can be sometimes challenging to design, and are often sensitive to variations in practical implementation.In this chapter, we provide a general overview of the ChIP method and its most common variations, with a special focus on ChIP-seq. We try to address some of the most important aspects that need to be taken into account in order to design and perform experiments that generate the most reproducible, high-quality data. Some of the main topics covered include the use of properly characterized antibodies, alternatives to chromatin preparation, the need for proper controls, and some recommendations about ChIP-seq data analysis.

Pho dynamically interacts with Spt5 to facilitate transcriptional switches at the hsp70 locus

Background

Numerous target genes of the Polycomb group (PcG) are transiently activated by a stimulus and subsequently repressed. However, mechanisms by which PcG proteins regulate such target genes remain elusive.

Results

We employed the heat shock-responsive hsp70 locus in Drosophila to study the chromatin dynamics of PRC1 and its interplay with known regulators of the locus before, during and after heat shock. We detected mutually exclusive binding patterns for HSF and PRC1 at the hsp70 locus. We found that Pleiohomeotic (Pho), a DNA-binding PcG member, dynamically interacts with Spt5, an elongation factor. The dynamic interaction switch between Pho and Spt5 is triggered by the recruitment of HSF to chromatin. Mutation in the protein–protein interaction domain (REPO domain) of Pho interferes with the dynamics of its interaction with Spt5. The transcriptional kinetics of the heat shock response is negatively affected by a mutation in the REPO domain of Pho.

Conclusions

We propose that a dynamic interaction switch between PcG proteins and an elongation factor enables stress-inducible genes to efficiently switch between ON/OFF states in the presence/absence of the activating stimulus.

Electronic supplementary material

The online version of this article (10.1186/s13072-017-0166-9) contains supplementary material, which is available to authorized users.

A Detailed Protocol for Subcellular RNA Sequencing (subRNA-seq)

In eukaryotic cells, RNAs at various maturation and processing levels are distributed across cellular compartments. The standard approach to determine transcript abundance and identity in vivo is RNA sequencing (RNA-seq). RNA-seq relies on RNA isolation from whole-cell lysates and thus mainly captures fully processed, stable, and more abundant cytoplasmic RNAs over nascent, unstable, and nuclear RNAs. Here, we provide a step-by-step protocol for subcellular RNA-seq (subRNA-seq). subRNA-seq allows the quantitative measurement of RNA polymerase II-generated RNAs from the chromatin, nucleoplasm, and cytoplasm of mammalian cells. This approach relies on cell fractionation prior to RNA isolation and sequencing library preparation. High-throughput sequencing of the subcellular RNAs can then be used to reveal the identity, abundance, and subcellular distribution of transcripts, thus providing insights into RNA processing and maturation. Deep sequencing of the chromatin-associated RNAs further offers the opportunity to study nascent RNAs. Subcellular RNA-seq libraries are obtained within 5 days. © 2017 by John Wiley & Sons, Inc.

Context-specific effects of genetic variants associated with autoimmune disease

Autoimmune diseases such as rheumatoid arthritis and coeliac disease are typical examples of complex genetic diseases caused by a combination of genetic and non-genetic risk factors. Insight into the genetic risk factors (single nucleotide polymorphisms (SNPs)) has increased since genome-wide association studies (GWAS) became possible in 2007 and, for individual diseases, SNPs can now explain some 15-50% of genetic risk. GWAS have also shown that some 50% of the genetic risk factors for individual autoimmune diseases overlap between different diseases. Thus, shared risk factors may converge to pathways that, when perturbed by genetic variation, predispose to autoimmunity in general. This raises the question of what determines disease specificity, and suggests that identical risk factors may have different effects in various autoimmune diseases. Addressing this question requires translation of genetic risk factors to causal genes and then to molecular and cellular pathways. Since >90% of the genetic risk factors are found in the non-coding part of the genome (i.e. outside the exons of protein-coding genes) and can have an impact on gene regulation, there is an urgent need to better understand the non-coding part of the genome. Here, we will outline the methods being used to unravel the gene regulatory networks perturbed in autoimmune diseases and the importance of doing this in the relevant cell types. We will highlight findings in coeliac disease, which manifests in the small intestine, to demonstrate how cell type and disease context can impact on the consequences of genetic risk factors.

The Ground State and Evolution of Promoter Region Directionality

Eukaryotic promoter regions are frequently divergently transcribed in vivo, but it is unknown whether the resultant antisense RNAs are a mechanistic by-product of RNA polymerase II (Pol II) transcription or biologically meaningful. Here, we use a functional evolutionary approach that involves nascent transcript mapping in S. cerevisiae strains containing foreign yeast DNA. Promoter regions in foreign environments lose the directionality they have in their native species. Strikingly, fortuitous promoter regions arising in foreign DNA produce equal transcription in both directions, indicating that divergent transcription is a mechanistic feature that does not imply a function for these transcripts. Fortuitous promoter regions arising during evolution promote bidirectional transcription and over time are purged through mutation or retained to enable new functionality. Similarly, human transcription is more bidirectional at newly evolved enhancers and promoter regions. Thus, promoter regions are intrinsically bidirectional and are shaped by evolution to bias transcription toward coding versus non-coding RNAs.

Pause & go: from the discovery of RNA polymerase pausing to its functional implications

The synthesis of nascent RNA is a discontinuous process in which phases of productive elongation by RNA polymerase are interrupted by frequent pauses. Transcriptional pausing was first observed decades ago, but was long considered to be a special feature of transcription at certain genes. This view was challenged when studies using genome-wide approaches revealed that RNA polymerase II pauses at promoter-proximal regions in large sets of genes in Drosophila and mammalian cells. High-resolution genomic methods uncovered that pausing is not restricted to promoters, but occurs globally throughout gene-body regions, implying the existence of key-rate limiting steps in nascent RNA synthesis downstream of transcription initiation. Here, we outline the experimental breakthroughs that led to the discovery of pervasive transcriptional pausing, discuss its emerging roles and regulation, and highlight the importance of pausing in human development and disease.

Total RNA-seq to identify pharmacological effects on specific stages of mRNA synthesis

Pharmacological perturbation is a powerful tool for understanding mRNA synthesis, but identification of the specific steps of this multi-step process that are targeted by small molecules remains challenging. Here we applied strand-specific total RNA sequencing (RNA-seq) to identify and distinguish specific pharmacological effects on transcription and pre-mRNA processing in human cells. We found unexpectedly that the natural product isoginkgetin, previously described as a splicing inhibitor, inhibits transcription elongation. Compared to well-characterized elongation inhibitors that target CDK9, isoginkgetin caused RNA polymerase accumulation within a broader promoter-proximal band, indicating that elongation inhibition by isoginkgetin occurs after release from promoter-proximal pause. RNA-seq distinguished isoginkgetin and CDK9 inhibitors from topoisomerase I inhibition, which alters elongation across gene bodies. We were able to detect these and other specific defects in mRNA synthesis at low sequencing depth using simple metagene-based metrics. These metrics now enable total-RNA-seq-based screening for high-throughput identification of pharmacological effects on individual stages of mRNA synthesis.