Physical isolation of nascent RNA chains transcribed by RNA polymerase II: evidence for cotranscriptional splicing

The regulation and function of post-transcriptional RNA splicing

Eukaryotic RNA transcripts undergo extensive processing before becoming functional messenger RNAs, with splicing being a critical and highly regulated step that occurs both co-transcriptionally and post-transcriptionally. Recent analyses have revealed, with unprecedented spatial and temporal resolution, that up to 40% of mammalian introns are retained after transcription termination and are subsequently removed largely while transcripts remain chromatin-associated. Post-transcriptional splicing has emerged as a key layer of gene expression regulation during development, stress response and disease progression. The control of post-transcriptional splicing regulates protein production through delayed splicing and nuclear export, or nuclear retention and degradation of specific transcript isoforms. Here, we review current methodologies for detecting post-transcriptional splicing, discuss the mechanisms controlling the timing of splicing and examine how this temporal regulation affects gene expression programmes in healthy cells and in disease states.

RNA m5C oxidation by TET2 regulates chromatin state and leukaemogenesis

Mutation of tet methylcytosine dioxygenase 2 (encoded by TET2) drives myeloid malignancy initiation and progression1-3. TET2 deficiency is known to cause a globally opened chromatin state and activation of genes contributing to aberrant haematopoietic stem cell self-renewal4,5. However, the open chromatin observed in TET2-deficient mouse embryonic stem cells, leukaemic cells and haematopoietic stem and progenitor cells5 is inconsistent with the designated role of DNA 5-methylcytosine oxidation of TET2. Here we show that chromatin-associated retrotransposon RNA 5-methylcytosine (m5C) can be recognized by the methyl-CpG-binding-domain protein MBD6, which guides deubiquitination of nearby monoubiquitinated Lys119 of histone H2A (H2AK119ub) to promote an open chromatin state. TET2 oxidizes m5C and antagonizes this MBD6-dependent H2AK119ub deubiquitination. TET2 depletion thereby leads to globally decreased H2AK119ub, more open chromatin and increased transcription in stem cells. TET2-mutant human leukaemia becomes dependent on this gene activation pathway, with MBD6 depletion selectively blocking proliferation of TET2-mutant leukaemic cells and largely reversing the haematopoiesis defects caused by Tet2 loss in mouse models. Together, our findings reveal a chromatin regulation pathway by TET2 through retrotransposon RNA m5C oxidation and identify the downstream MBD6 protein as a feasible target for developing therapies specific against TET2 mutant malignancies.

Genome organization around nuclear speckles drives mRNA splicing efficiency

The nucleus is highly organized, such that factors involved in the transcription and processing of distinct classes of RNA are confined within specific nuclear bodies1,2. One example is the nuclear speckle, which is defined by high concentrations of protein and noncoding RNA regulators of pre-mRNA splicing3. What functional role, if any, speckles might play in the process of mRNA splicing is unclear4,5. Here we show that genes localized near nuclear speckles display higher spliceosome concentrations, increased spliceosome binding to their pre-mRNAs and higher co-transcriptional splicing levels than genes that are located farther from nuclear speckles. Gene organization around nuclear speckles is dynamic between cell types, and changes in speckle proximity lead to differences in splicing efficiency. Finally, directed recruitment of a pre-mRNA to nuclear speckles is sufficient to increase mRNA splicing levels. Together, our results integrate the long-standing observations of nuclear speckles with the biochemistry of mRNA splicing and demonstrate a crucial role for dynamic three-dimensional spatial organization of genomic DNA in driving spliceosome concentrations and controlling the efficiency of mRNA splicing.

The chromatin - triple helix connection

Mammalian genomes are extensively transcribed, producing a large number of coding and non-coding transcripts. A large fraction of the nuclear RNAs is physically associated with chromatin, functioning in gene activation and silencing, shaping higher-order genome organisation, such as involvement in long-range enhancer-promoter interactions, transcription hubs, heterochromatin, nuclear bodies and phase transitions. Different mechanisms allow the tethering of these chromatin-associated RNAs (caRNA) to chromosomes, including RNA binding proteins, the RNA polymerases and R-loops. In this review, we focus on the sequence-specific targeting of RNA to DNA by forming triple helical structures and describe its interplay with chromatin. It turns out that nucleosome positioning at triple helix target sites and the nucleosome itself are essential factors in determining the formation and stability of triple helices. The histone H3-tail plays a critical role in triple helix stabilisation, and the role of its epigenetic modifications in this process is discussed.

METTL14 is a chromatin regulator independent of its RNA N6-methyladenosine methyltransferase activity

METTL3 and METTL14 are two components that form the core heterodimer of the main RNA m6A methyltransferase complex (MTC) that installs m6A. Surprisingly, depletion of METTL3 or METTL14 displayed distinct effects on stemness maintenance of mouse embryonic stem cell (mESC). While comparable global hypo-methylation in RNA m6A was observed in Mettl3 or Mettl14 knockout mESCs, respectively. Mettl14 knockout led to a globally decreased nascent RNA synthesis, whereas Mettl3 depletion resulted in transcription upregulation, suggesting that METTL14 might possess an m6A-independent role in gene regulation. We found that METTL14 colocalizes with the repressive H3K27me3 modification. Mechanistically, METTL14, but not METTL3, binds H3K27me3 and recruits KDM6B to induce H3K27me3 demethylation independent of METTL3. Depletion of METTL14 thus led to a global increase in H3K27me3 level along with a global gene suppression. The effects of METTL14 on regulation of H3K27me3 is essential for the transition from self-renewal to differentiation of mESCs. This work reveals a regulatory mechanism on heterochromatin by METTL14 in a manner distinct from METTL3 and independently of m6A, and critically impacts transcriptional regulation, stemness maintenance, and differentiation of mESCs.

RBFOX2 recognizes N6-methyladenosine to suppress transcription and block myeloid leukaemia differentiation

N6-methyladenosine (m6A) methylation can be deposited on chromatin-associated RNAs (caRNAs) by the RNA methyltransferase complex (MTC) to regulate chromatin state and transcription. However, the mechanism by which MTC is recruited to distinct genomic loci remains elusive. Here we identify RBFOX2, a well-studied RNA-binding protein, as a chromatin factor that preferentially recognizes m6A on caRNAs. RBFOX2 can recruit RBM15, an MTC component, to facilitate methylation of promoter-associated RNAs. RBM15 also physically interacts with YTHDC1 and recruits polycomb repressive complex 2 (PRC2) to the RBFOX2-bound loci for chromatin silencing and transcription suppression. Furthermore, we found that this RBFOX2/m6A/RBM15/YTHDC1/PRC2 axis plays a critical role in myeloid leukaemia. Downregulation of RBFOX2 notably inhibits survival/proliferation of acute myeloid leukaemia cells and promotes their myeloid differentiation. RBFOX2 is also required for self-renewal of leukaemia stem/initiation cells and acute myeloid leukaemia maintenance. Our study presents a pathway of m6A MTC recruitment and m6A deposition on caRNAs, resulting in locus-selective chromatin regulation, which has potential therapeutic implications in leukaemia.

3D genome organization around nuclear speckles drives mRNA splicing efficiency

The nucleus is highly organized such that factors involved in transcription and processing of distinct classes of RNA are organized within specific nuclear bodies. One such nuclear body is the nuclear speckle, which is defined by high concentrations of protein and non-coding RNA regulators of pre-mRNA splicing. What functional role, if any, speckles might play in the process of mRNA splicing remains unknown. Here we show that genes localized near nuclear speckles display higher spliceosome concentrations, increased spliceosome binding to their pre-mRNAs, and higher co-transcriptional splicing levels relative to genes that are located farther from nuclear speckles. We show that directed recruitment of a pre-mRNA to nuclear speckles is sufficient to drive increased mRNA splicing levels. Finally, we show that gene organization around nuclear speckles is highly dynamic with differential localization between cell types corresponding to differences in Pol II occupancy. Together, our results integrate the longstanding observations of nuclear speckles with the biochemistry of mRNA splicing and demonstrate a critical role for dynamic 3D spatial organization of genomic DNA in driving spliceosome concentrations and controlling the efficiency of mRNA splicing.

A simple and robust method for isolating and analyzing chromatin-bound RNAs in Arabidopsis

BACKGROUND

Chromatin-bound RNAs are the primary product of transcription that undergo on-chromatin processing such as capping, splicing, and polyadenylation. These processing steps then determine the fate of RNAs. Albeit its vital importance, a simple and robust method for isolating different fractions of chromatin-bound RNAs is missing in plants.

RESULT

Here, we describe our updated method and the associated step-by-step protocol for chromatin-bound RNAs isolation in A. thaliana. The chromatin-bound RNAs isolation is based on the 1 M UREA wash that removes the majority of non-chromatin-associated proteins from the nucleus, as previously developed in mammalian cells. On-demand, the isolated chromatin-bound RNAs can be either used directly for gene-specific analysis or subject to further rRNA removal and also the optional polyadenylated RNA removal, followed by high-throughput sequencing. Detailed protocols for these procedures are also provided. Comparison of sequencing results of chromatin-bound RNAs with and without polyadenylated RNA removal revealed that a small fraction of CB-RNAs is polyadenylated but not yet fully spliced, representing RNA-processing intermediate on-chromatin.

CONCLUSION

This optimized chromatin-bound RNAs purification method is simple and robust and can be used to study transcription and its-coupled RNA processing in plants.

ZC3H4 restricts non-coding transcription in human cells

The human genome encodes thousands of non-coding RNAs. Many of these terminate early and are then rapidly degraded, but how their transcription is restricted is poorly understood. In a screen for protein-coding gene transcriptional termination factors, we identified ZC3H4. Its depletion causes upregulation and extension of hundreds of unstable transcripts, particularly antisense RNAs and those transcribed from so-called super-enhancers. These loci are occupied by ZC3H4, suggesting that it directly functions in their transcription. Consistently, engineered tethering of ZC3H4 to reporter RNA promotes its degradation by the exosome. ZC3H4 is predominantly metazoan –interesting when considering its impact on enhancer RNAs that are less prominent in single-celled organisms. Finally, ZC3H4 loss causes a substantial reduction in cell proliferation, highlighting its overall importance. In summary, we identify ZC3H4 as playing an important role in restricting non-coding transcription in multicellular organisms.

Tracking pre-mRNA maturation across subcellular compartments identifies developmental gene regulation through intron retention and nuclear anchoring

Steps of mRNA maturation are important gene regulatory events that occur in distinct cellular locations. However, transcriptomic analyses often lose information on the subcellular distribution of processed and unprocessed transcripts. We generated extensive RNA-seq data sets to track mRNA maturation across subcellular locations in mouse embryonic stem cells, neuronal progenitor cells, and postmitotic neurons. We find disparate patterns of RNA enrichment between the cytoplasmic, nucleoplasmic, and chromatin fractions, with some genes maintaining more polyadenylated RNA in chromatin than in the cytoplasm. We bioinformatically defined four regulatory groups for intron retention, including complete cotranscriptional splicing, complete intron retention in the cytoplasmic RNA, and two intron groups present in nuclear and chromatin transcripts but fully excised in cytoplasm. We found that introns switch their regulatory group between cell types, including neuronally excised introns repressed by polypyrimidine track binding protein 1 (PTBP1). Transcripts for the neuronal gamma-aminobutyric acid (GABA) B receptor, 1 (Gabbr1) are highly expressed in mESCs but are absent from the cytoplasm. Instead, incompletely spliced Gabbr1 RNA remains sequestered on chromatin, where it is bound by PTBP1, similar to certain long noncoding RNAs. Upon neuronal differentiation, Gabbr1 RNA becomes fully processed and exported for translation. Thus, splicing repression and chromatin anchoring of RNA combine to allow posttranscriptional regulation of Gabbr1 over development. For this and other genes, polyadenylated RNA abundance does not indicate functional gene expression. Our data sets provide a rich resource for analyzing many other aspects of mRNA maturation in subcellular locations and across development.

Intron exon boundary junctions in human genome have in-built unique structural and energetic signals

Precise identification of correct exon–intron boundaries is a prerequisite to analyze the location and structure of genes. The existing framework for genomic signals, delineating exon and introns in a genomic segment, seems insufficient, predominantly due to poor sequence consensus as well as limitations of training on available experimental data sets. We present here a novel concept for characterizing exon–intron boundaries in genomic segments on the basis of structural and energetic properties. We analyzed boundary junctions on both sides of all the exons (3 28 368) of protein coding genes from human genome (GENCODE database) using 28 structural and three energy parameters. Study of sequence conservation at these sites shows very poor consensus. It is observed that DNA adopts a unique structural and energy state at the boundary junctions. Also, signals are somewhat different for housekeeping and tissue specific genes. Clustering of 31 parameters into four derived vectors gives some additional insights into the physical mechanisms involved in this biological process. Sites of structural and energy signals correlate well to the positions playing important roles in pre-mRNA splicing.

Subcellular Localization of uc.8+ as a Prognostic Biomarker in Bladder Cancer Tissue

Simple Summary

DNA regions having high sequence similarity among human, rat and mouse genomes are defined as Ultraconserved Regions. Non-coding RNA transcripts originating by these regions may play relevant roles in the onset and progression of multiple cancer types. We recently found that ultra-conserved-transcript-8+ (uc.8+) levels correlate with the grading and staging of bladder cancer. The aim of this study is to systematically evaluate the expression of ultra-conserved-transcript-8+ (uc.8+) in biopsies and assess its intracellular localization. Furthermore, we aimed to correlate uc.8+ levels with clinical parameters and patient survival. Our analysis indicates that uc.8+ can localize both in the cytoplasm and nucleus of bladder cells at early stages of tumorigenesis, while in tumors at advanced stages, uc.8+ has a prevalent cytoplasmic localization. These data provide relevant information about uc.8+ localization as a hallmark of tumor stage. Finally, using advanced computer-based techniques, we predicted the binding of uc.8+ to RNA-binding proteins. Our study overall suggests that uc.8+ localization can be used as a prognostic biomarker for bladder cancer.

Abstract

Non-coding RNA transcripts originating from Ultraconserved Regions (UCRs) have tissue-specific expression and play relevant roles in the pathophysiology of multiple cancer types. Among them, we recently identified and characterized the ultra-conserved-transcript-8+ (uc.8+), whose levels correlate with grading and staging of bladder cancer. Here, to validate uc.8+ as a potential biomarker in bladder cancer, we assessed its expression and subcellular localization by using tissue microarray on 73 human bladder cancer specimens. We quantified uc.8+ by in-situ hybridization and correlated its expression levels with clinical characteristics and patient survival. The analysis of subcellular localization indicated the simultaneous presence of uc.8+ in the cytoplasm and nucleus of cells from the Low-Grade group, whereas a prevalent cytoplasmic localization was observed in samples from the High-Grade group, supporting the hypothesis of uc.8+ nuclear-to-cytoplasmic translocation in most malignant tumor forms. Moreover, analysis of uc.8+ expression and subcellular localization in tumor-surrounding stroma revealed a marked down-regulation of uc.8+ levels compared to the paired (adjacent) tumor region. Finally, deep machine-learning approaches identified nucleotide sequences associated with uc.8+ localization in nucleus and/or cytoplasm, allowing to predict possible RNA binding proteins associated with uc.8+, recognizing also sequences involved in mRNA cytoplasm-translocation. Our model suggests uc.8+ subcellular localization as a potential prognostic biomarker for bladder cancer.

Chromatin regulatory dynamics of early human small intestinal development using a directed differentiation model

The establishment of the small intestinal (SI) lineage during human embryogenesis ensures functional integrity of the intestine after birth. The chromatin dynamics that drive SI lineage formation and regional patterning in humans are essentially unknown. To fill this knowledge void, we apply a cutting-edge genomic technology to a state-of-the-art human model of early SI development. Specifically, we leverage chromatin run-on sequencing (ChRO-seq) to define the landscape of active promoters, enhancers and gene bodies across distinct stages of directed differentiation of human pluripotent stem cells into SI spheroids with regional specification. Through comprehensive ChRO-seq analysis we identify candidate stage-specific chromatin activity states, novel markers and enhancer hotspots during the directed differentiation. Moreover, we propose a detailed transcriptional network associated with SI lineage formation or regional patterning. Our ChRO-seq analyses uncover a previously undescribed pattern of enhancer activity and transcription at HOX gene loci underlying SI regional patterning. We also validated this unique HOX dynamics by the analysis of single cell RNA-seq data from human fetal SI. Overall, the results lead to a new proposed working model for the regulatory underpinnings of human SI development, thereby adding a novel dimension to the literature that has relied almost exclusively on non-human models.

TGFβ promotes widespread enhancer chromatin opening and operates on genomic regulatory domains

The Transforming Growth Factor-β (TGFβ) signaling pathway controls transcription by regulating enhancer activity. How TGFβ-regulated enhancers are selected and what chromatin changes are associated with TGFβ-dependent enhancers regulation are still unclear. Here we report that TGFβ treatment triggers fast and widespread increase in chromatin accessibility in about 80% of the enhancers of normal mouse mammary epithelial-gland cells, irrespective of whether they are activated, repressed or not regulated by TGFβ. This enhancer opening depends on both the canonical and non-canonical TGFβ pathways. Most TGFβ-regulated genes are located around enhancers regulated in the same way, often creating domains of several co-regulated genes that we term TGFβ regulatory domains (TRD). CRISPR-mediated inactivation of enhancers within TRDs impairs TGFβ-dependent regulation of all co-regulated genes, demonstrating that enhancer targeting is more promiscuous than previously anticipated. The area of TRD influence is restricted by topologically associating domains (TADs) borders, causing a bias towards co-regulation within TADs.

Preparation of Mammalian Nascent RNA for Long Read Sequencing

Long read sequencing technologies now allow high-quality sequencing of RNAs (or their cDNAs) that are hundreds to thousands of nucleotides long. Long read sequences of nascent RNA provide single-nucleotide-resolution information about co-transcriptional RNA processing events-e.g., splicing, folding, and base modifications. Here, we describe how to isolate nascent RNA from mammalian cells through subcellular fractionation of chromatin-associated RNA, as well as how to deplete poly(A)+ RNA and rRNA, and, finally, how to generate a full-length cDNA library for use on long read sequencing platforms. This approach allows for an understanding of coordinated splicing status across multi-intron transcripts by revealing patterns of splicing or other RNA processing events that cannot be gained from traditional short read RNA sequencing. © 2020 Wiley Periodicals LLC. Basic Protocol 1: Subcellular fractionation Basic Protocol 2: Nascent RNA isolation and adapter ligation Basic Protocol 3: cDNA amplicon preparation.

METTL3 and N6-Methyladenosine Promote Homologous Recombination-Mediated Repair of DSBs by Modulating DNA-RNA Hybrid Accumulation

Double-strand breaks (DSBs) are the most deleterious DNA lesions, which, if left unrepaired, may lead to genome instability or cell death. Here, we report that, in response to DSBs, the RNA methyltransferase METTL3 is activated by ATM-mediated phosphorylation at S43. Phosphorylated METTL3 is then localized to DNA damage sites, where it methylates the N6 position of adenosine (m6A) in DNA damage-associated RNAs, which recruits the m6A reader protein YTHDC1 for protection. In this way, the METTL3-m6A-YTHDC1 axis modulates accumulation of DNA-RNA hybrids at DSBs sites, which then recruit RAD51 and BRCA1 for homologous recombination (HR)-mediated repair. METTL3-deficient cells display defective HR, accumulation of unrepaired DSBs, and genome instability. Accordingly, depletion of METTL3 significantly enhances the sensitivity of cancer cells and murine xenografts to DNA damage-based therapy. These findings uncover the function of METTL3 and YTHDC1 in HR-mediated DSB repair, which may have implications for cancer therapy.

Gene-Specific Control of tRNA Expression by RNA Polymerase II

Increasing evidence suggests that tRNA levels are dynamically and specifically regulated in response to internal and external cues to modulate the cellular translational program. However, the molecular players and the mechanisms regulating the gene-specific expression of tRNAs are still unknown. Using an inducible auxin-degron system to rapidly deplete RPB1 (the largest subunit of RNA Pol II) in living cells, we identified Pol II as a direct gene-specific regulator of tRNA transcription. Our data suggest that Pol II transcription robustly interferes with Pol III function at specific tRNA genes. This activity was further found to be essential for MAF1-mediated repression of a large set of tRNA genes during serum starvation, indicating that repression of tRNA genes by Pol II is dynamically regulated. Hence, Pol II plays a direct and central role in the gene-specific regulation of tRNA expression.

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.

Seq'ing identity and function in a repeat-derived noncoding RNA world

Innovations in high-throughout sequencing approaches are being marshaled to both reveal the composition of the abundant and heterogeneous noncoding RNAs that populate cell nuclei and lend insight to the mechanisms by which noncoding RNAs influence chromosome biology and gene expression. This review focuses on some of the recent technological developments that have enabled the isolation of nascent transcripts and chromatin-associated and DNA-interacting RNAs. Coupled with emerging genome assembly and analytical approaches, the field is poised to achieve a comprehensive catalog of nuclear noncoding RNAs, including those derived from repetitive regions within eukaryotic genomes. Herein, particular attention is paid to the challenges and advances in the sequence analyses of repeat and transposable element-derived noncoding RNAs and in ascribing specific function(s) to such RNAs.

The Features and Regulation of Co-transcriptional Splicing in Arabidopsis

Precursor mRNA (pre-mRNA) splicing is essential for gene expression in most eukaryotic organisms. Previous studies from mammals, Drosophila, and yeast show that the majority of splicing events occurs co-transcriptionally. In plants, however, the features of co-transcriptional splicing (CTS) and its regulation still remain largely unknown. Here, we used chromatin-bound RNA sequencing to study CTS in Arabidopsis thaliana. We found that CTS is widespread in Arabidopsis seedlings, with a large proportion of alternative splicing events determined co-transcriptionally. CTS efficiency correlated with gene expression level, the chromatin landscape and, most surprisingly, the number of introns and exons of individual genes, but is independent of gene length. In combination with enhanced crosslinking and immunoprecipitation sequencing analysis, we further showed that the hnRNP-like proteins RZ-1B and RZ-1C promote efficient CTS globally through direct binding, frequently to exonic sequences. Notably, this general effect of RZ-1B/1C on splicing promotion is mainly observed at the chromatin level, not at the mRNA level. RZ-1C promotes CTS of multiple-exon genes in association with its binding to regions both proximal and distal to the regulated introns. We propose that RZ-1C promotes efficient CTS of genes with multiple exons through cooperative interactions with many exons, introns, and splicing factors. Our work thus reveals important features of CTS in plants and provides methodologies for the investigation of CTS and RNA-binding proteins in plants.

N 6-methyladenosine of chromosome-associated regulatory RNA regulates chromatin state and transcription

N ⁶-methyladenosine (m⁶A) regulates stability and translation of messenger RNA (mRNA) in various biological processes. In this work, we show that knockout of the m⁶A writer Mettl3 or the nuclear reader Ythdc1 in mouse embryonic stem cells increases chromatin accessibility and activates transcription in an m⁶A-dependent manner. We found that METTL3 deposits m⁶A modifications on chromosome-associated regulatory RNAs (carRNAs), including promoter-associated RNAs, enhancer RNAs, and repeat RNAs. YTHDC1 facilitates the decay of a subset of these m⁶A-modified RNAs, especially elements of the long interspersed element-1 family, through the nuclear exosome targeting-mediated nuclear degradation. Reducing m⁶A methylation by METTL3 depletion or site-specific m⁶A demethylation of selected carRNAs elevates the levels of carRNAs and promotes open chromatin state and downstream transcription. Collectively, our results reveal that m⁶A on carRNAs can globally tune chromatin state and transcription.

Directed RNase H Cleavage of Nascent Transcripts Causes Transcription Termination

An attractive approach to reduce gene expression is via the use of antisense oligonucleotides (ASOs) that harness the RNase H1 mechanism. Here we show that RNase H ASOs targeted to introns or exons robustly reduce the level of spliced RNA associated with chromatin. Surprisingly, intron-targeted ASOs reduce the level of pre-mRNA associated with chromatin to a greater extent than exon-targeted ASOs. This indicates that exon-targeted ASOs achieve full activity after the pre-mRNA has undergone splicing, but before the mRNA is released from chromatin. Even though RNase H ASOs can reduce the level of RNA associated with chromatin, the effect of ASO-directed RNA degradation on transcription has never been documented. Here we show that intron-targeted ASOs and, to a lesser extent, exon-targeted ASOs cause RNA polymerase II (Pol II) transcription termination in cultured cells and mice. Furthermore, ASO-directed transcription termination is mediated by the nuclear exonuclease XRN2.

Histone lysine demethylase KDM4B regulates the alternative splicing of the androgen receptor in response to androgen deprivation

Alternative splicing is emerging as an oncogenic mechanism. In prostate cancer, generation of constitutively active forms of androgen receptor (AR) variants including AR-V7 plays an important role in progression of castration-resistant prostate cancer (CRPC). AR-V7 is generated by alternative splicing that results in inclusion of cryptic exon CE3 and translation of truncated AR protein that lacks the ligand binding domain. Whether AR-V7 can be a driver for CRPC remains controversial as the oncogenic mechanism of AR-V7 activation remains elusive. Here, we found that KDM4B promotes AR-V7 and identified a novel regulatory mechanism. KDM4B is phosphorylated by protein kinase A under conditions that promote castration-resistance, eliciting its binding to the splicing factor SF3B3. KDM4B binds RNA specifically near the 5'-CE3, upregulates the chromatin accessibility, and couples the spliceosome to the chromatin. Our data suggest that KDM4B can function as a signal responsive trans-acting splicing factor and scaffold that recruits and stabilizes the spliceosome near the alternative exon, thus promoting its inclusion. Genome-wide profiling of KDM4B-regulated genes also identified additional alternative splicing events implicated in tumorigenesis. Our study defines KDM4B-regulated alternative splicing as a pivotal mechanism for generating AR-V7 and a contributing factor for CRPC, providing insight for mechanistic targeting of CRPC.

Epigenetic regulation of TIMP1 expression by 8-oxoguanine DNA glycosylase-1 binding to DNA:RNA hybrid

8-Oxoguanine DNA glycosylase-1 (OGG1)-initiated base excision repair pathway is primarily responsible for 7, 8-dihydro-8-oxoguanine (8-oxoG) removal from DNA. Recent studies, however, have shown that 8-oxoG in gene regulatory elements may serve as an epigenetic mark, and OGG1 has distinct functions in modulating gene expression. Genome-wide mapping of oxidative stress-induced OGG1 enrichment within introns was documented, but its significance has not yet been fully characterized. Here, we explored whether OGG1 recruited to intron 1 of tissue inhibitor of metalloproteinase-1 (TIMP1) gene and modulated its expression. Using chromatin and DNA:RNA hybrid immunoprecipitation assays, we report recruitment of OGG1 to the DNA:RNA hybrid in intron 1, where it increases nascent RNA but lowers mRNA levels in O₃-exposed human airway epithelial cells and mouse lungs. Decrease in TIMP1 expression is alleviated by antioxidant administration, small interfering RNA depletion, or inhibition of OGG1 binding to its genomic substrate. In vitro studies revealed direct interaction between OGG1 and 8-oxoG containing DNA:RNA hybrid, without excision of its substrate. Inhibition of OGG1 binding to DNA:RNA hybrid translated into an increase in TIMP1 expression and a decrease in oxidant-induced lung inflammatory responses as well as airway remodeling. Data documented here reveal a novel molecular link between OGG1 at damaged sites and transcription dynamics that may contribute to oxidative stress-induced cellular and tissue responses.-Pan, L., Wang, H., Luo, J., Zeng, J., Pi, J., Liu, H., Liu, C., Ba, X., Qu, X., Xiang, Y., Boldogh, I., Qin, X. Epigenetic regulation of TIMP1 expression by 8-oxoguanine DNA glycosylase-1 binding to DNA:RNA hybrid.

Promoter-proximal pausing of RNA polymerase II: a nexus of gene regulation

In this review, Core et al. discuss the recent advances in our understanding of the early steps in Pol II transcription, highlighting the events and factors involved in the establishment and release of paused Pol II. They also discuss a number of unanswered questions about the regulation and function of Pol II pausing.

Precise spatio–temporal control of gene activity is essential for organismal development, growth, and survival in a changing environment. Decisive steps in gene regulation involve the pausing of RNA polymerase II (Pol II) in early elongation, and the controlled release of paused polymerase into productive RNA synthesis. Here we describe the factors that enable pausing and the events that trigger Pol II release into the gene. We also discuss open questions in the field concerning the stability of paused Pol II, nucleosomes as obstacles to elongation, and potential roles of pausing in defining the precision and dynamics of gene expression.

RNA m6 A modification enzymes shape innate responses to DNA by regulating interferon β

In this study, Rubio et al. researched how the dynamic genome-wide landscape of m⁶A-modified mRNAs impacts virus infection and host immune responses. They show that type I interferon (IFN) production triggered by dsDNA or human cytomegalovirus (HCMV) is controlled by the cellular m⁶A methyltrasferase subunit METTL14 and ALKBH5 demethylase, and their results demonstrate that responses to nonmicrobial dsDNA in uninfected cells are regulated by enzymes controlling m⁶A epitranscriptomic changes.

Modification of mRNA by N⁶-adenosine methylation (m⁶A) on internal bases influences gene expression in eukaryotes. How the dynamic genome-wide landscape of m⁶A-modified mRNAs impacts virus infection and host immune responses remains poorly understood. Here, we show that type I interferon (IFN) production triggered by dsDNA or human cytomegalovirus (HCMV) is controlled by the cellular m⁶A methyltrasferase subunit METTL14 and ALKBH5 demethylase. While METTL14 depletion reduced virus reproduction and stimulated dsDNA- or HCMV-induced IFNB1 mRNA accumulation, ALKBH5 depletion had the opposite effect. Depleting METTL14 increased both nascent IFNB1 mRNA production and stability in response to dsDNA. In contrast, ALKBH5 depletion reduced nascent IFNB1 mRNA production without detectably influencing IFN1B mRNA decay. Genome-wide transcriptome profiling following ALKBH5 depletion identified differentially expressed genes regulating antiviral immune responses, while METTL14 depletion altered pathways impacting metabolic reprogramming, stress responses, and aging. Finally, we determined that IFNB1 mRNA was m⁶A-modified within both the coding sequence and the 3′ untranslated region (UTR). This establishes that the host m⁶A modification machinery controls IFNβ production triggered by HCMV or dsDNA. Moreover, it demonstrates that responses to nonmicrobial dsDNA in uninfected cells, which shape host immunity and contribute to autoimmune disease, are regulated by enzymes controlling m⁶A epitranscriptomic changes.

Chromatin run-on and sequencing maps the transcriptional regulatory landscape of glioblastoma multiforme

The human genome encodes a variety of poorly understood RNA species that remain challenging to identify using existing genomic tools. We developed chromatin run-on and sequencing (ChRO-seq) to map the location of RNA polymerase for almost any input sample, including samples with degraded RNA that are intractable to RNA sequencing. We used ChRO-seq to map nascent transcription in primary human glioblastoma (GBM) brain tumors. Enhancers identified in primary GBMs resemble open chromatin in the normal human brain. Rare enhancers that are activated in malignant tissue drive regulatory programs similar to the developing nervous system. We identified enhancers that regulate groups of genes that are characteristic of each known GBM subtype and transcription factors that drive them. Finally we discovered a core group of transcription factors that control the expression of genes associated with clinical outcomes. This study characterizes the transcriptional landscape of GBM and introduces ChRO-seq as a method to map regulatory programs that contribute to complex diseases.

Pre-mRNA modifications and their role in nuclear processing

Background

Cellular non-coding RNAs are extensively modified post-transcriptionally, with more than 100 chemically distinct nucleotides identified to date. In the past five years, new sequencing based methods have revealed widespread decoration of eukaryotic messenger RNA with diverse RNA modifications whose functions in mRNA metabolism are only beginning to be known.

Results

Since most of the identified mRNA modifying enzymes are present in the nucleus, these modifications have the potential to function in nuclear pre-mRNA processing including alternative splicing. Here we review recent progress towards illuminating the role of pre-mRNA modifications in splicing and highlight key areas for future investigation in this rapidly growing field.

Conclusions

Future studies to identify which modifications are added to nascent pre-mRNA and to interrogate the direct effects of individual modifications are likely to reveal new mechanisms by which nuclear pre-mRNA processing is regulated.

Microprocessor-dependent processing of splice site overlapping microRNA exons does not result in changes in alternative splicing

MicroRNAs are found throughout the genome and are processed by the microprocessor complex (MPC) from longer precursors. Some precursor miRNAs overlap intron:exon junctions. These splice site overlapping microRNAs (SO-miRNAs) are mostly located in coding genes. It has been intimated, in the rarer examples of SO-miRNAs in noncoding RNAs, that the competition between the spliceosome and the MPC modulates alternative splicing. However, the effect of this overlap on coding transcripts is unknown. Unexpectedly, we show that neither Drosha silencing nor SF3b1 silencing changed the inclusion ratio of SO-miRNA exons. Two SO-miRNAs, located in genes that code for basal membrane proteins, are known to inhibit proliferation in primary keratinocytes. These SO-miRNAs were up-regulated during differentiation and the host mRNAs were down-regulated, but again there was no change in inclusion ratio of the SO-miRNA exons. Interestingly, Drosha silencing increased nascent RNA density, on chromatin, downstream from SO-miRNA exons. Overall our data suggest a novel mechanism for regulating gene expression in which MPC-dependent cleavage of SO-miRNA exons could cause premature transcriptional termination of coding genes rather than affecting alternative splicing.

Nuclear PTEN safeguards pre-mRNA splicing to link Golgi apparatus for its tumor suppressive role

Dysregulation of pre-mRNA alternative splicing (AS) is closely associated with cancers. However, the relationships between the AS and classic oncogenes/tumor suppressors are largely unknown. Here we show that the deletion of tumor suppressor PTEN alters pre-mRNA splicing in a phosphatase-independent manner, and identify 262 PTEN-regulated AS events in 293T cells by RNA sequencing, which are associated with significant worse outcome of cancer patients. Based on these findings, we report that nuclear PTEN interacts with the splicing machinery, spliceosome, to regulate its assembly and pre-mRNA splicing. We also identify a new exon 2b in GOLGA2 transcript and the exon exclusion contributes to PTEN knockdown-induced tumorigenesis by promoting dramatic Golgi extension and secretion, and PTEN depletion significantly sensitizes cancer cells to secretion inhibitors brefeldin A and golgicide A. Our results suggest that Golgi secretion inhibitors alone or in combination with PI3K/Akt kinase inhibitors may be therapeutically useful for PTEN-deficient cancers.

Mechanism of N6-methyladenosine modification and its emerging role in cancer

N⁶-methyladenosine (m⁶A), the most prevalent internal methylation in messenger RNA (mRNA) that is deposited by m⁶A methyltransferases, removed by m⁶A demethylases and recognized by different RNA-binding proteins, distinguishes the transcripts through multilayer interactions with mRNA processing, export, degradation and translation machineries. m⁶A plays an important role in regulation of gene expression for fundamental cellular processes and diverse physiological functions. Aberrant m⁶A decorations lead to cancer but also have the potential to yield new therapies. This review outlines the evolution of the m⁶A field, formation of key concepts, important open questions and also discusses the molecular basis of mRNA m⁶A modification and its effect in cancer, highlighting the potential of demethylase as a therapeutic target for cancer treatment.

Human TFIIH Kinase CDK7 Regulates Transcription-Associated Chromatin Modifications

CDK7 phosphorylates the RNA polymerase II (pol II) C-terminal domain CTD and activates the P-TEFb-associated kinase CDK9, but its regulatory roles remain obscure. Here, using human CDK7 analog-sensitive (CDK7as) cells, we observed reduced capping enzyme recruitment, increased pol II promoter-proximal pausing, and defective termination at gene 3' ends upon CDK7 inhibition. We also noted that CDK7 regulates chromatin modifications downstream of transcription start sites. H3K4me3 spreading was restricted at gene 5' ends and H3K36me3 was displaced toward gene 3' ends in CDK7as cells. Mass spectrometry identified factors that bound TFIIH-phosphorylated versus P-TEFb-phosphorylated CTD (versus unmodified); capping enzymes and H3K4 methyltransferase complexes, SETD1A/B, selectively bound phosphorylated CTD, and the H3K36 methyltransferase SETD2 specifically bound P-TEFb-phosphorylated CTD. Moreover, TFIIH-phosphorylated CTD stimulated SETD1A/B activity toward nucleosomes, revealing a mechanistic basis for CDK7 regulation of H3K4me3 spreading. Collectively, these results implicate a CDK7-dependent "CTD code" that regulates chromatin marks in addition to RNA processing and pol II pausing.

Pre-mRNA processing includes N6 methylation of adenosine residues that are retained in mRNA exons and the fallacy of "RNA epigenetics"

By using a cell fraction technique that separates chromatin-associated nascent RNA, newly completed nucleoplasmic mRNA and cytoplasmic mRNA, we have shown in a previous study that residues in exons are methylated (m⁶A) in nascent pre-mRNA and remain methylated in the same exonic residues in nucleoplasmic and cytoplasmic mRNA. Thus, there is no evidence of a substantial degree of demethylation in mRNA exons that would correspond to so-called "epigenetic" demethylation. The turnover rate of mRNA molecules is faster, depending on m⁶A content in HeLa cell mRNA, suggesting that specification of mRNA stability may be the major role of m⁶A exon modification. In mouse embryonic stem cells (mESCs) lacking Mettl3, the major mRNA methylase, the cells continue to grow, making the same mRNAs with unchanged splicing profiles in the absence (>90%) of m⁶A in mRNA, suggesting no common obligatory role of m⁶A in splicing. All these data argue strongly against a commonly used "reversible dynamic methylation/demethylation" of mRNA, calling into question the concept of "RNA epigenetics" that parallels the well-established role of dynamic DNA epigenetics.

Measuring Nascent Transcripts by Nascent-seq

A complete understanding of transcription and co-transcriptional RNA processing events by polymerase requires precise and robust approaches to visualize polymerase progress and quantify nascent transcripts on a genome-wide scale. Here, we present a transcriptome-wide method to measure the level of nascent transcribing RNA in a fast and unbiased manner.

m6A mRNA modifications are deposited in nascent pre-mRNA and are not required for splicing but do specify cytoplasmic turnover

Understanding the biologic role of N⁶-methyladenosine (m⁶A) RNA modifications in mRNA requires an understanding of when and where in the life of a pre-mRNA transcript the modifications are made. We found that HeLa cell chromatin-associated nascent pre-mRNA (CA-RNA) contains many unspliced introns and m⁶A in exons but very rarely in introns. The m⁶A methylation is essentially completed upon the release of mRNA into the nucleoplasm. Furthermore, the content and location of each m⁶A modification in steady-state cytoplasmic mRNA are largely indistinguishable from those in the newly synthesized CA-RNA or nucleoplasmic mRNA. This result suggests that quantitatively little methylation or demethylation occurs in cytoplasmic mRNA. In addition, only ∼10% of m⁶As in CA-RNA are within 50 nucleotides of 5' or 3' splice sites, and the vast majority of exons harboring m⁶A in wild-type mouse stem cells is spliced the same in cells lacking the major m⁶A methyltransferase Mettl3. Both HeLa and mouse embryonic stem cell mRNAs harboring m⁶As have shorter half-lives, and thousands of these mRNAs have increased half-lives (twofold or more) in Mettl3 knockout cells compared with wild type. In summary, m⁶A is added to exons before or soon after exon definition in nascent pre-mRNA, and while m⁶A is not required for most splicing, its addition in the nascent transcript is a determinant of cytoplasmic mRNA stability.

Promiscuous DNA-binding of a mutant zinc finger protein corrupts the transcriptome and diminishes cell viability

The rules of engagement between zinc finger transcription factors and DNA have been partly defined by in vitro DNA-binding and structural studies, but less is known about how these rules apply in vivo. Here, we demonstrate how a missense mutation in the second zinc finger of Krüppel-like factor-1 (KLF1) leads to degenerate DNA-binding specificity in vivo, resulting in ectopic transcription and anemia in the Nan mouse model. We employed ChIP-seq and 4sU-RNA-seq to identify aberrant DNA-binding events genome wide and ectopic transcriptional consequences of this binding. We confirmed novel sequence specificity of the mutant recombinant zinc finger domain by performing biophysical measurements of in vitro DNA-binding affinity. Together, these results shed new light on the mechanisms by which missense mutations in DNA-binding domains of transcription factors can lead to autosomal dominant diseases.

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.

Mot1, Ino80C, and NC2 Function Coordinately to Regulate Pervasive Transcription in Yeast and Mammals

Pervasive transcription initiates from cryptic promoters and is observed in eukaryotes ranging from yeast to mammals. The Set2-Rpd3 regulatory system prevents cryptic promoter function within expressed genes. However, conserved systems that control pervasive transcription within intergenic regions have not been well established. Here we show that Mot1, Ino80 chromatin remodeling complex (Ino80C), and NC2 co-localize on chromatin and coordinately suppress pervasive transcription in S. cerevisiae and murine embryonic stem cells (mESCs). In yeast, all three proteins bind subtelomeric heterochromatin through a Sir3-stimulated mechanism and to euchromatin via a TBP-stimulated mechanism. In mESCs, the proteins bind to active and poised TBP-bound promoters along with promoters of polycomb-silenced genes apparently lacking TBP. Depletion of Mot1, Ino80C, or NC2 by anchor away in yeast or RNAi in mESCs leads to near-identical transcriptome phenotypes, with new subtelomeric transcription in yeast, and greatly increased pervasive transcription in both yeast and mESCs.

Post-transcriptional regulation of mouse neurogenesis by Pumilio proteins

Despite extensive studies on mammalian neurogenesis, its post-transcriptional regulation remains under-explored. Here we report that neural-specific inactivation of two murine post-transcriptional regulators, Pumilio 1 (Pum1) and Pum2, severely reduced the number of neural stem cells (NSCs) in the postnatal dentate gyrus (DG), drastically increased perinatal apoptosis, altered DG cell composition, and impaired learning and memory. Consistently, the mutant DG neurospheres generated fewer NSCs with defects in proliferation, survival, and differentiation, supporting a major role of Pum1 and Pum2 in hippocampal neurogenesis and function. Cross-linking immunoprecipitation revealed that Pum1 and Pum2 bind to thousands of mRNAs, with at least 694 common targets in multiple neurogenic pathways. Depleting Pum1 and/or Pum2 did not change the abundance of most target mRNAs but up-regulated their proteins, indicating that Pum1 and Pum2 regulate the translation of their target mRNAs. Moreover, Pum1 and Pum2 display RNA-dependent interaction with fragile X mental retardation protein (FMRP) and bind to one another's mRNA. This indicates that Pum proteins might form collaborative networks with FMRP and possibly other post-transcriptional regulators to regulate neurogenesis.

Chromatin-enriched lncRNAs can act as cell-type specific activators of proximal gene transcription

We recently described a new class of long noncoding RNA defined by especially tight chromatin association, whose presence is strongly correlated with expression of nearby genes in HEK293 cells. Here we critically examine the generality and cis-enhancer mechanism of this class of chromatin enriched RNA (cheRNA). CheRNA are largely cell-type specific, and remain the most effective chromatin signature for predicting cis-gene transcription in all cell types examined. Targeted depletion of three cheRNAs decreases gene expression of their neighbors, indicating potential co-activator function. Single-molecule FISH of one cheRNA-distal target gene pair suggests spatial overlap consistent with a role in chromosome looping. In another example, the cheRNA HIDALGO stimulates the fetal hemoglobin HBG1 gene during erythroid differentiation by promoting contacts to a downstream enhancer. Our results suggest that many cheRNAs activate proximal, lineage-specific gene transcription.

m6A Demethylase ALKBH5 Maintains Tumorigenicity of Glioblastoma Stem-like Cells by Sustaining FOXM1 Expression and Cell Proliferation Program

The dynamic and reversible N⁶-methyladenosine (m⁶A) RNA modification installed and erased by N⁶-methyltransferases and demethylases regulates gene expression and cell fate. We show that the m⁶A demethylase ALKBH5 is highly expressed in glioblastoma stem-like cells (GSCs). Silencing ALKBH5 suppresses the proliferation of patient-derived GSCs. Integrated transcriptome and m⁶A-seq analyses revealed altered expression of certain ALKBH5 target genes, including the transcription factor FOXM1. ALKBH5 demethylates FOXM1 nascent transcripts, leading to enhanced FOXM1 expression. Furthermore, a long non-coding RNA antisense to FOXM1 (FOXM1-AS) promotes the interaction of ALKBH5 with FOXM1 nascent transcripts. Depleting ALKBH5 and FOXM1-AS disrupted GSC tumorigenesis through the FOXM1 axis. Our work uncovers a critical function for ALKBH5 and provides insight into critical roles of m⁶A methylation in glioblastoma.

Pharmacological inhibition of the spliceosome subunit SF3b triggers exon junction complex-independent nonsense-mediated decay

Spliceostatin A, meayamycin, and pladienolide B are small molecules that target the SF3b subunit of the spliceosomal U2 small nuclear ribonucleoprotein (snRNP). These compounds are attracting much attention as tools to manipulate splicing and for use as potential anti-cancer drugs. We investigated the effects of these inhibitors on mRNA transport and stability in human cells. Upon splicing inhibition, unspliced pre-mRNAs accumulated in the nucleus, particularly within enlarged nuclear speckles. However, a small fraction of the pre-mRNA molecules were exported to the cytoplasm. We identified the export adaptor ALYREF as being associated with intron-containing transcripts and show its requirement for the nucleo-cytoplasmic transport of unspliced pre-mRNA. In contrast, the exon junction complex (EJC) core protein eIF4AIII failed to form a stable complex with intron-containing transcripts. Despite the absence of EJC, unspliced transcripts in the cytoplasm were degraded by nonsense-mediated decay (NMD), suggesting that unspliced transcripts are degraded by an EJC-independent NMD pathway. Collectively, our results indicate that although blocking the function of SF3b elicits a massive accumulation of unspliced pre-mRNAs in the nucleus, intron-containing transcripts can still bind the ALYREF export factor and be transported to the cytoplasm, where they trigger an alternative NMD pathway.

Cellular Fractionation and Isolation of Chromatin-Associated RNA

In eukaryotic cells, the synthesis, processing, and functions of RNA molecules are confined to distinct subcellular compartments. Biochemical fractionation of cells prior to RNA isolation thus enables the analysis of distinct steps in the lifetime of individual RNA molecules that would be masked in bulk RNA preparations from whole cells. Here, we describe a simple two-step differential centrifugation protocol for the isolation of cytoplasmic, nucleoplasmic, and chromatin-associated RNA that can be used in downstream applications such as qPCR or deep sequencing. We discuss various aspects of this fractionation protocol, which can be readily applied to many mammalian cell types. For the study of long noncoding RNAs and enhancer RNAs in regulation of transcription especially the preparation of chromatin-associated RNA can contribute significantly to further developments.

Rbfox Proteins Regulate Splicing as Part of a Large Multiprotein Complex LASR

Rbfox proteins control alternative splicing and posttranscriptional regulation in mammalian brain and are implicated in neurological disease. These proteins recognize the RNA sequence (U)GCAUG, but their structures and diverse roles imply a variety of protein-protein interactions. We find that nuclear Rbfox proteins are bound within a large assembly of splicing regulators (LASR), a multimeric complex containing the proteins hnRNP M, hnRNP H, hnRNP C, Matrin3, NF110/NFAR-2, NF45, and DDX5, all approximately equimolar to Rbfox. We show that splicing repression mediated by hnRNP M is stimulated by Rbfox. Virtually all the intron-bound Rbfox is associated with LASR, and hnRNP M motifs are enriched adjacent to Rbfox crosslinking sites in vivo. These findings demonstrate that Rbfox proteins bind RNA with a defined set of cofactors and affect a broader set of exons than previously recognized. The function of this multimeric LASR complex has implications for deciphering the regulatory codes controlling splicing networks.

TRIBE: Hijacking an RNA-Editing Enzyme to Identify Cell-Specific Targets of RNA-Binding Proteins

RNA transcripts are bound and regulated by RNA-binding proteins (RBPs). Current methods for identifying in vivo targets of an RBP are imperfect and not amenable to examining small numbers of cells. To address these issues, we developed TRIBE (targets of RNA-binding proteins identified by editing), a technique that couples an RBP to the catalytic domain of the Drosophila RNA-editing enzyme ADAR and expresses the fusion protein in vivo. RBP targets are marked with novel RNA editing events and identified by sequencing RNA. We have used TRIBE to identify the targets of three RBPs (Hrp48, dFMR1, and NonA). TRIBE compares favorably to other methods, including CLIP, and we have identified RBP targets from as little as 150 specific fly neurons. TRIBE can be performed without an antibody and in small numbers of specific cells.

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

Many features of how gene transcription occurs in human cells remain unclear, mainly because of a lack of quantitative approaches to follow genome transcription with nucleotide precision in vivo. Here we present a robust genome-wide approach for studying RNA polymerase II (Pol II)-mediated transcription in human cells at single-nucleotide resolution by native elongating transcript sequencing (NET-seq). Elongating RNA polymerase and the associated nascent RNA are prepared by cell fractionation, avoiding immunoprecipitation or RNA labeling. The 3' ends of nascent RNAs are captured through barcode linker ligation and converted into a DNA sequencing library. The identity and abundance of the 3' ends are determined by high-throughput sequencing, which reveals the exact genomic locations of Pol II. Human NET-seq can be applied to the study of the full spectrum of Pol II transcriptional activities, including the production of unstable RNAs and transcriptional pausing. By using the protocol described here, a NET-seq library can be obtained from human cells in 5 d.

Mitotic Transcriptional Activation: Clearance of Actively Engaged Pol II via Transcriptional Elongation Control in Mitosis

Although it is established that some general transcription factors are inactivated at mitosis, many details of mitotic transcription inhibition (MTI) and its underlying mechanisms are largely unknown. We have identified mitotic transcriptional activation (MTA) as a key regulatory step to control transcription in mitosis for genes with transcriptionally engaged RNA polymerase II (Pol II) to activate and transcribe until the end of the gene to clear Pol II from mitotic chromatin, followed by global impairment of transcription reinitiation through MTI. Global nascent RNA sequencing and RNA fluorescence in situ hybridization demonstrate the existence of transcriptionally engaged Pol II in early mitosis. Both genetic and chemical inhibition of P-TEFb in mitosis lead to delays in the progression of cell division. Together, our study reveals a mechanism for MTA and MTI whereby transcriptionally engaged Pol II can progress into productive elongation and finish transcription to allow proper cellular division.

RNA Binding Proteins RZ-1B and RZ-1C Play Critical Roles in Regulating Pre-mRNA Splicing and Gene Expression during Development in Arabidopsis

Nuclear-localized RNA binding proteins are involved in various aspects of RNA metabolism, which in turn modulates gene expression. However, the functions of nuclear-localized RNA binding proteins in plants are poorly understood. Here, we report the functions of two proteins containing RNA recognition motifs, RZ-1B and RZ-1C, in Arabidopsis thaliana. RZ-1B and RZ-1C were localized to nuclear speckles and interacted with a spectrum of serine/arginine-rich (SR) proteins through their C termini. RZ-1C preferentially bound to purine-rich RNA sequences in vitro through its N-terminal RNA recognition motif. Disrupting the RNA binding activity of RZ-1C with SR proteins through overexpression of the C terminus of RZ-1C conferred defective phenotypes similar to those observed in rz-1b rz-1c double mutants, including delayed seed germination, reduced stature, and serrated leaves. Loss of function of RZ-1B and RZ-1C was accompanied by defective splicing of many genes and global perturbation of gene expression. In addition, we found that RZ-1C directly targeted FLOWERING LOCUS C (FLC), promoting efficient splicing of FLC introns and likely also repressing FLC transcription. Our findings highlight the critical role of RZ-1B/1C in regulating RNA splicing, gene expression, and many key aspects of plant development via interaction with proteins including SR proteins.

Nuclear Fractionation Reveals Thousands of Chromatin-Tethered Noncoding RNAs Adjacent to Active Genes

A number of long noncoding RNAs (lncRNAs) have been reported to regulate transcription via recruitment of chromatin modifiers or bridging distal enhancer elements to gene promoters. However, the generality of these modes of regulation and the mechanisms of chromatin attachment for thousands of unstudied human lncRNAs remain unclear. To address these questions, we performed stringent nuclear fractionation coupled to RNA sequencing. We provide genome-wide identification of human chromatin-associated lncRNAs and demonstrate tethering of RNA to chromatin by RNAPII is a pervasive mechanism of attachment. We also uncovered thousands of chromatin-enriched RNAs (cheRNAs) that share molecular properties with known lncRNAs. Although distinct from eRNAs derived from active prototypical enhancers, the production of cheRNAs is strongly correlated with the expression of neighboring protein-coding genes. This work provides an updated framework for nuclear RNA organization that includes a large chromatin-associated transcript population correlated with active genes and may prove useful in de novo enhancer annotation.

Native elongating transcript sequencing reveals human transcriptional activity at nucleotide resolution

Major features of transcription by human RNA polymerase II (Pol II) remain poorly defined due to a lack of quantitative approaches for visualizing Pol II progress at nucleotide resolution. We developed a simple and powerful approach for performing native elongating transcript sequencing (NET-seq) in human cells that globally maps strand-specific Pol II density at nucleotide resolution. NET-seq exposes a mode of antisense transcription that originates downstream and converges on transcription from the canonical promoter. Convergent transcription is associated with a distinctive chromatin configuration and is characteristic of lower-expressed genes. Integration of NET-seq with genomic footprinting data reveals stereotypic Pol II pausing coincident with transcription factor occupancy. Finally, exons retained in mature transcripts display Pol II pausing signatures that differ markedly from skipped exons, indicating an intrinsic capacity for Pol II to recognize exons with different processing fates. Together, human NET-seq exposes the topography and regulatory complexity of human gene expression.