Metabolic labeling of newly transcribed RNA for high resolution gene expression profiling of RNA synthesis, processing and decay in cell culture

Gene-specific RNA homeostasis revealed by perturbation of coactivator complexes

Transcript buffering entails the reciprocal modulation of mRNA synthesis and degradation rates to maintain stable RNA levels under varying cellular conditions. Current research supports a global, non-sequence-specific connection between mRNA synthesis and degradation, but the underlying mechanisms are still unclear. In this study, we investigated changes in RNA metabolism following acute depletion of TIP60/KAT5, the acetyltransferase subunit of the NuA4 transcriptional coactivator complex, in mouse embryonic stem cells. By combining RNA sequencing of nuclear, cytoplasmic, and newly synthesised transcript fractions with biophysical modelling, we demonstrate that TIP60 predominantly enhances transcription of numerous genes, while a smaller set of genes undergoes TIP60-dependent transcriptional repression. Surprisingly, transcription changes caused by TIP60 depletion were offset by corresponding changes in RNA nuclear export and cytoplasmic stability, indicating gene-specific buffering mechanisms. Similarly, disruption of the unrelated ATAC coactivator complex also resulted in gene-specific transcript buffering. These findings reveal that transcript buffering functions at a gene-specific level and suggest that cells dynamically adjust RNA splicing, export, and degradation in response to individual RNA synthesis alterations, thereby sustaining cellular homeostasis.

Development and Comparison of 4-Thiouridine to Cytidine Base Conversion Reaction

To study transcriptome dynamics without harming cells, it is essential to convert chemical bases. 4-Thiouridine (4sU) is a biocompatible uridine analogue that can be converted into a cytidine analogue. Although several reactions can convert 4sU into a cytidine analogue, few studies have compared the features of these reactions. In this study, we performed three reported base conversion reactions, including osmium tetroxide, iodoacetamide, and sodium periodate treatment, as well as a new reaction using 2,4-dinitrofluorobenzene. We compared the reaction time, conversion efficacy, and effects on reverse transcription. These reactions successfully converted 4sU into a cytidine analogue quantitatively using trinucleotides. However, the conversion efficacy and effect on reverse transcription vary depending on the reaction with the RNA transcript. OsO4 treatment followed by NH4Cl treatment showed the best base-conversion efficiency. Nevertheless, each reaction has its own advantages and disadvantages as a tool for studying the transcriptome. Therefore, it is crucial to select the appropriate reaction for the target of interest.

hGRAD: A versatile "one-fits-all" system to acutely deplete RNA binding proteins from condensates

Nuclear RNA binding proteins (RBPs) are difficult to study because they often belong to large protein families and form extensive networks of auto- and crossregulation. They are highly abundant and many localize to condensates with a slow turnover, requiring long depletion times or knockouts that cannot distinguish between direct and indirect or compensatory effects. Here, we developed a system that is optimized for the rapid degradation of nuclear RBPs, called hGRAD. It comes as a "one-fits-all" plasmid, and integration into any cell line with endogenously GFP-tagged proteins allows for an inducible, rapid, and complete knockdown. We show that the nuclear RBPs SRSF3, SRSF5, SRRM2, and NONO are completely cleared from nuclear speckles and paraspeckles within 2 h. hGRAD works in various cell types, is more efficient than previous methods, and does not require the expression of exogenous ubiquitin ligases. Combining SRSF5 hGRAD degradation with Nascent-seq uncovered transient transcript changes, compensatory mechanisms, and an effect of SRSF5 on transcript stability.

Dynamically regulated transcription factors are encoded by highly unstable mRNAs in the Drosophila larval brain

The level of each RNA species depends on the balance between its rates of production and decay. Although previous studies have measured RNA decay across the genome in tissue culture and single-celled organisms, few experiments have been performed in intact complex tissues and organs. It is therefore unclear whether the determinants of RNA decay found in cultured cells are preserved in an intact tissue, and whether they differ between neighboring cell types and are regulated during development. To address these questions, we measured RNA synthesis and decay rates genome wide via metabolic labeling of whole cultured Drosophila larval brains using 4-thiouridine. Our analysis revealed that decay rates span a range of more than 100-fold, and that RNA stability is linked to gene function, with mRNAs encoding transcription factors being much less stable than mRNAs involved in core metabolic functions. Surprisingly, among transcription factor mRNAs there was a clear demarcation between more widely used transcription factors and those that are expressed only transiently during development. mRNAs encoding transient transcription factors are among the least stable in the brain. These mRNAs are characterized by epigenetic silencing in most cell types, as shown by their enrichment with the histone modification H3K27me3. Our data suggest the presence of an mRNA destabilizing mechanism targeted to these transiently expressed transcription factors to allow their levels to be regulated rapidly with high precision. Our study also demonstrates a general method for measuring mRNA transcription and decay rates in intact organs or tissues, offering insights into the role of mRNA stability in the regulation of complex developmental programs.

Splicing activates transcription from weak promoters upstream of alternative exons

Transcription and splicing are intrinsically coupled. Alternative splicing of internal exons can fine-tune gene expression through a recently described phenomenon called exon-mediated activation of transcription starts (EMATS). However, the association of this phenomenon with human diseases remains unknown. Here, we develop a strategy to activate gene expression through EMATS and demonstrate its potential for treatment of genetic diseases caused by loss of expression of essential genes. We first identified a catalog of human EMATS genes and provide a list of their pathological variants. To test if EMATS can be used to activate gene expression, we constructed stable cell lines expressing a splicing reporter based on the alternative splicing of motor neuron 2 (SMN2) gene. Using small molecules and antisense oligonucleotides (ASOs) currently used for treatment of spinal muscular atrophy, we demonstrated that increase of inclusion of alternative exons can trigger an activation of gene expression up to 45-fold by enhancing transcription in EMATS-like genes. We observed the strongest effects in genes under the regulation of weak human promoters located proximal to highly included skipped exons.

Measuring the tail: Methods for poly(A) tail profiling

The 3'-end poly(A) tail is an important and potent feature of most mRNA molecules that affects mRNA fate and translation efficiency. Polyadenylation is a posttranscriptional process that occurs in the nucleus by canonical poly(A) polymerases (PAPs). In some specific instances, the poly(A) tail can also be extended in the cytoplasm by noncanonical poly(A) polymerases (ncPAPs). This epitranscriptomic regulation of mRNA recently became one of the most interesting aspects in the field. Advances in RNA sequencing technologies and software development have allowed the precise measurement of poly(A) tails, identification of new ncPAPs, expansion of the function of known enzymes, discovery and a better understanding of the physiological role of tail heterogeneity, and recognition of a correlation between tail length and RNA translatability. Here, we summarize the development of polyadenylation research methods, including classic low-throughput approaches, Illumina-based genome-wide analysis, and advanced state-of-art techniques that utilize long-read third-generation sequencing with Pacific Biosciences and Oxford Nanopore Technologies platforms. A boost in technical opportunities over recent decades has allowed a better understanding of the regulation of gene expression at the mRNA level. This article is categorized under: RNA Methods > RNA Analyses In Vitro and In Silico.

Parental mutations influence wild-type offspring via transcriptional adaptation

Transgenerational epigenetic inheritance (TEI) is mostly discussed in the context of physiological or environmental factors. Here, we show intergenerational and transgenerational inheritance of transcriptional adaptation (TA), a process whereby mutant messenger RNA (mRNA) degradation affects gene expression, in nematodes and zebrafish. Wild-type offspring of animals heterozygous for mRNA-destabilizing alleles display increased expression of adapting genes. Notably, offspring of animals heterozygous for nontranscribing alleles do not display this response. Germline-specific mutations are sufficient to induce TA in wild-type offspring, indicating that, at least for some genes, mutations in somatic tissues are not necessary for this process. Microinjecting total RNA from germ cells of TA-displaying heterozygous zebrafish can trigger TA in wild-type embryos and in their progeny, suggesting a model whereby mutant mRNAs in the germline trigger a TA response that can be epigenetically inherited. In sum, this previously unidentified mode of TEI reveals a means by which parental mutations can modulate the offspring's transcriptome.

SRSF1 acts as an IFN-I-regulated cellular dependency factor decisively affecting HIV-1 post-integration steps

Efficient HIV-1 replication depends on balanced levels of host cell components including cellular splicing factors as the family of serine/arginine-rich splicing factors (SRSF, 1-10). Type I interferons (IFN-I) play a crucial role in the innate immunity against HIV-1 by inducing the expression of IFN-stimulated genes (ISGs) including potent host restriction factors. The less well known IFN-repressed genes (IRepGs) might additionally affect viral replication by downregulating host dependency factors that are essential for the viral life cycle; however, so far, the knowledge about IRepGs involved in HIV-1 infection is very limited. In this work, we could demonstrate that HIV-1 infection and the associated ISG induction correlated with low SRSF1 levels in intestinal lamina propria mononuclear cells (LPMCs) and peripheral blood mononuclear cells (PBMCs) during acute and chronic HIV-1 infection. In HIV-1-susceptible cell lines as well as primary monocyte-derived macrophages (MDMs), expression levels of SRSF1 were transiently repressed upon treatment with specific IFNα subtypes in vitro. Mechanically, 4sU labeling of newly transcribed mRNAs revealed that IFN-mediated SRSF1 repression is regulated on early RNA level. SRSF1 knockdown led to an increase in total viral RNA levels, but the relative proportion of the HIV-1 viral infectivity factor (Vif) coding transcripts, which is essential to counteract APOBEC3G-mediated host restriction, was significantly reduced. In the presence of high APOBEC3G levels, however, increased LTR activity upon SRSF1 knockdown facilitated the overall replication, despite decreased vif mRNA levels. In contrast, SRSF1 overexpression significantly impaired HIV-1 post-integration steps including LTR transcription, alternative splice site usage, and virus particle production. Since balanced SRSF1 levels are crucial for efficient viral replication, our data highlight the so far undescribed role of SRSF1 acting as an IFN-modulated cellular dependency factor decisively regulating HIV-1 post-integration steps.

Metabolic labeling of cardiomyocyte-derived small extracellular-vesicle (sEV) miRNAs identifies miR-208a in cardiac regulation of lung gene expression

Toxoplasma gondii uracil phosphoribosyltransferase (UPRT) converts 4-thiouracil (4TUc) into 4-thiouridine (4TUd), which is incorporated into nascent RNAs and can be biotinylated, then labelled with streptavidin conjugates or isolated via streptavidin-affinity methods. Here, we generated mice that expressed T. gondii UPRT only in cardiomyocytes (CM UPRT mice) and tested our hypothesis that CM-derived miRNAs (CM miRs) are transferred into remote organs after myocardial infarction (MI) by small extracellular vesicles (sEV) that are released from the heart into the peripheral blood (PB sEV). We found that 4TUd was incorporated with high specificity and sensitivity into RNAs isolated from the hearts and PB sEV of CM UPRT mice 6 h after 4TUc injection. In PB sEV, 4TUd was incorporated into CM-specific/enriched miRs including miR-208a, but not into miRs with other organ or tissue-type specificities. 4TUd-labelled miR208a was also present in lung tissues, especially lung endothelial cells (ECs), and CM-derived miR-208a (CM miR-208a) levels peaked 12 h after experimentally induced MI in PB sEV and 24 h after MI in the lung. Notably, miR-208a is expressed from intron 29 of α myosin heavy chain (αMHC), but αMHC transcripts were nearly undetectable in the lung. When PB sEV from mice that underwent MI (MI-PB sEV) or sham surgery (Sham-PB sEV) were injected into intact mice, the expression of Tmbim6 and NLK, which are suppressed by miR-208a and cooperatively regulate inflammation via the NF-κB pathway, was lower in the lungs of MI-PB sEV-treated animals than the lungs of animals treated with Sham-PB sEV or saline. In MI mice, Tmbim6 and NLK were downregulated, whereas endothelial adhesion molecules and pro-inflammatory cells were upregulated in the lung; these changes were significantly attenuated when the mice were treated with miR-208a antagomirs prior to MI surgery. Thus, CM UPRT mice enables us to track PB sEV-mediated transport of CM miRs and identify an miR-208a-mediated mechanism by which myocardial injury alters the expression of genes and inflammatory response in the lung.

Slow RNAPII Transcription Elongation Rate, Low Levels of RNAPII Pausing, and Elevated Histone H1 Content at Promoters Associate with Higher m6A Deposition on Nascent mRNAs

N6-methyladenosine modification (m6A) fine-tunes RNA fate in a variety of ways, thus regulating multiple fundamental biological processes. m6A writers bind to chromatin and interact with RNA polymerase II (RNAPII) during transcription. To evaluate how the dynamics of the transcription process impact m6A deposition, we studied RNAPII elongation rates in mouse embryonic stem cells with altered chromatin configurations, due to reductions in linker histone H1 content. We found that genes transcribed at slow speed are preferentially methylated and display unique signatures at their promoter region, namely high levels of histone H1, together with marks of bivalent chromatin and low RNAPII pausing. They are also highly susceptible to m6A loss upon histone H1 reduction. These results indicate that RNAPII velocity links chromatin structure and the deposition of m6A, highlighting the intricate relationship between different regulatory layers on nascent mRNA molecules.

CDK11 regulates pre-mRNA splicing by phosphorylation of SF3B1

RNA splicing, the process of intron removal from pre-mRNA, is essential for the regulation of gene expression. It is controlled by the spliceosome, a megadalton RNA-protein complex that assembles de novo on each pre-mRNA intron through an ordered assembly of intermediate complexes^1,2. Spliceosome activation is a major control step that requires substantial protein and RNA rearrangements leading to a catalytically active complex^1-5. Splicing factor 3B subunit 1 (SF3B1) protein-a subunit of the U2 small nuclear ribonucleoprotein⁶-is phosphorylated during spliceosome activation^7-10, but the kinase that is responsible has not been identified. Here we show that cyclin-dependent kinase 11 (CDK11) associates with SF3B1 and phosphorylates threonine residues at its N terminus during spliceosome activation. The phosphorylation is important for the association between SF3B1 and U5 and U6 snRNAs in the activated spliceosome, termed the B^act complex, and the phosphorylation can be blocked by OTS964, a potent and selective inhibitor of CDK11. Inhibition of CDK11 prevents spliceosomal transition from the precatalytic complex B to the activated complex B^act and leads to widespread intron retention and accumulation of non-functional spliceosomes on pre-mRNAs and chromatin. We demonstrate a central role of CDK11 in spliceosome assembly and splicing regulation and characterize OTS964 as a highly selective CDK11 inhibitor that suppresses spliceosome activation and splicing.

Neuronal Cells Display Distinct Stability Controls of Alternative Polyadenylation mRNA Isoforms, Long Non-Coding RNAs, and Mitochondrial RNAs

RNA stability plays an important role in gene expression. Here, using 3' end sequencing of newly made and pre-existing poly(A)+ RNAs, we compare transcript stability in multiple human cell lines, including HEK293T, HepG2, and SH-SY5Y. We show that while mRNA stability is generally conserved across the cell lines, specific transcripts having a high GC content and possibly more stable secondary RNA structures are relatively more stable in SH-SY5Y cells compared to the other 2 cell lines. These features also differentiate stability levels of alternative polyadenylation (APA) 3'UTR isoforms in a cell type-specific manner. Using differentiation of a neural stem cell line as a model, we show that mRNA stability difference could contribute to gene expression changes in neurogenesis and confirm the neuronal identity of SH-SY5Y cells at both gene expression and APA levels. In addition, compared to transcripts using 3'-most exon cleavage/polyadenylation sites (PASs), those using intronic PASs are generally less stable, especially when the PAS-containing intron is large and has a strong 5' splice site, suggesting that intronic polyadenylation mostly plays a negative role in gene expression. Interestingly, the differential mRNA stability among APA isoforms appears to buffer PAS choice in these cell lines. Moreover, we found that several other poly(A)+ RNA species, including promoter-associated long noncoding RNAs and transcripts encoded by the mitochondrial genome, are more stable in SH-SY5Y cells than the other 2 cell lines, further highlighting distinct RNA metabolism in neuronal cells. Together, our results indicate that distinct RNA stability control in neuronal cells may contribute to the gene expression and APA programs that define their cell identity.

Enhancer RNA Expression in Response to Glucocorticoid Treatment in Murine Macrophages

Glucocorticoids are potent anti-inflammatory drugs; however, their molecular mode of action remains complex and elusive. They bind to the glucocorticoid receptor (GR), a nuclear receptor that controls gene expression in almost all tissues in a cell type-specific manner. While GR's transcriptional targets mediate beneficial reactions in immune cells, they also harbor the potential of adverse metabolic effects in other cell types such as hepatocytes. Here, we have profiled nascent transcription upon glucocorticoid stimulation in LPS-activated primary murine macrophages using 4sU-seq. We compared our results to publicly available nascent transcriptomics data from murine liver and bioinformatically identified non-coding RNAs transcribed from intergenic GR binding sites in a tissue-specific fashion. These tissue-specific enhancer RNAs (eRNAs) correlate with target gene expression, reflecting cell type-specific glucocorticoid responses. We further associate GR-mediated eRNA expression with changes in H3K27 acetylation and BRD4 recruitment in inflammatory macrophages upon glucocorticoid treatment. In summary, we propose a common mechanism by which GR-bound enhancers regulate target gene expression by changes in histone acetylation, BRD4 recruitment and eRNA expression. We argue that local eRNAs are potential therapeutic targets downstream of GR signaling which may modulate glucocorticoid response in a cell type-specific way.

Genome-wide identification and characterization of circular RNA m6A modification in pancreatic cancer

BACKGROUND

N⁶-methyladenosine (m⁶A) is the most abundant modification of RNA in eukaryotic cells and play critical roles in cancer. While most related studies focus on m⁶A modifications in linear RNA, transcriptome-wide profiling and exploration of m⁶A modification in circular RNAs in cancer is still lacking.

METHODS

For the detection of m⁶A modification in circRNAs, we developed a new bioinformatics tools called Circm6A and applied it to the m⁶A-seq data of 77 tissue samples from 58 individuals with pancreatic ductal adenocarcinoma (PDAC).

RESULTS

Circm6A performs better than the existing circRNA identification tools, which achieved highest F1 score among these tools in the detection of circRNAs with m⁶A modifications. By using Circm6A, we identified a total of 8807 m⁶A-circRNAs from our m⁶A-seq data. The m⁶A-circRNAs tend to be hypermethylated in PDAC tumor tissues compared with normal tissues. The hypermethylated m⁶A-circRNAs were associated with a significant gain of circRNA-mRNA coexpression network, leading to the dysregulation of many important cancer-related pathways. Moreover, we found the cues that hypermethylated m⁶A-circRNAs may promote the circularization and translation of circRNAs.

CONCLUSIONS

These comprehensive findings further bridged the knowledge gaps between m⁶A modification and circRNAs fields by depicting the m⁶A-circRNAs genomic landscape of PDAC patients and revealed the emerging roles played by m⁶A-circRNAs in pancreatic cancer. Circm6A is available at https://github.com/canceromics/circm6a .

RN7SK small nuclear RNA controls bidirectional transcription of highly expressed gene pairs in skin

Pausing of RNA polymerase II (Pol II) close to promoters is a common regulatory step in RNA synthesis, and is coordinated by a ribonucleoprotein complex scaffolded by the noncoding RNA RN7SK. The function of RN7SK-regulated gene transcription in adult tissue homoeostasis is currently unknown. Here, we deplete RN7SK during mouse and human epidermal stem cell differentiation. Unexpectedly, loss of this small nuclear RNA specifically reduces transcription of numerous cell cycle regulators leading to cell cycle exit and differentiation. Mechanistically, we show that RN7SK is required for efficient transcription of highly expressed gene pairs with bidirectional promoters, which in the epidermis co-regulated cell cycle and chromosome organization. The reduction in transcription involves impaired splicing and RNA decay, but occurs in the absence of chromatin remodelling at promoters and putative enhancers. Thus, RN7SK is directly required for efficient Pol II transcription of highly transcribed bidirectional gene pairs, and thereby exerts tissue-specific functions, such as maintaining a cycling cell population in the epidermis.

The noncoding RNA RN7SK regulates RNA polymerase II pausing and splicing. Here the authors deplete RN7SK in mouse and human during epidermal stem cell differentiation and reveal a novel role in orchestrating bidirectional transcription of highly expressed gene pairs.

The related coactivator complexes SAGA and ATAC control embryonic stem cell self-renewal through acetyltransferase-independent mechanisms

SAGA (Spt-Ada-Gcn5 acetyltransferase) and ATAC (Ada-two-A-containing) are two related coactivator complexes, sharing the same histone acetyltransferase (HAT) subunit. The HAT activities of SAGA and ATAC are required for metazoan development, but the role of these complexes in RNA polymerase II transcription is less understood. To determine whether SAGA and ATAC have redundant or specific functions, we compare the effects of HAT inactivation in each complex with that of inactivation of either SAGA or ATAC core subunits in mouse embryonic stem cells (ESCs). We show that core subunits of SAGA or ATAC are required for complex assembly and mouse ESC growth and self-renewal. Surprisingly, depletion of HAT module subunits causes a global decrease in histone H3K9 acetylation, but does not result in significant phenotypic or transcriptional defects. Thus, our results indicate that SAGA and ATAC are differentially required for self-renewal of mouse ESCs by regulating transcription through different pathways in a HAT-independent manner.

Clip Chemistry: Diverse (Bio)(macro)molecular and Material Function through Breaking Covalent Bonds

In the two decades since the introduction of the "click chemistry" concept, the toolbox of "click reactions" has continually expanded, enabling chemists, materials scientists, and biologists to rapidly and selectively build complexity for their applications of interest. Similarly, selective and efficient covalent bond breaking reactions have provided and will continue to provide transformative advances. Here, we review key examples and applications of efficient, selective covalent bond cleavage reactions, which we refer to herein as "clip reactions." The strategic application of clip reactions offers opportunities to tailor the compositions and structures of complex (bio)(macro)molecular systems with exquisite control. Working in concert, click chemistry and clip chemistry offer scientists and engineers powerful methods to address next-generation challenges across the chemical sciences.

Dynamic imaging of nascent RNA reveals general principles of transcription dynamics and stochastic splice site selection

The activities of RNA polymerase and the spliceosome are responsible for the heterogeneity in the abundance and isoform composition of mRNA in human cells. However, the dynamics of these megadalton enzymatic complexes working in concert on endogenous genes have not been described. Here, we establish a quasi-genome-scale platform for observing synthesis and processing kinetics of single nascent RNA molecules in real time. We find that all observed genes show transcriptional bursting. We also observe large kinetic variation in intron removal for single introns in single cells, which is inconsistent with deterministic splice site selection. Transcriptome-wide footprinting of the U2AF complex, nascent RNA profiling, long-read sequencing, and lariat sequencing further reveal widespread stochastic recursive splicing within introns. We propose and validate a unified theoretical model to explain the general features of transcription and pervasive stochastic splice site selection.

ELAV/Hu RNA binding proteins determine multiple programs of neural alternative splicing

ELAV/Hu factors are conserved RNA binding proteins (RBPs) that play diverse roles in mRNA processing and regulation. The founding member, Drosophila Elav, was recognized as a vital neural factor 35 years ago. Nevertheless, little was known about its impacts on the transcriptome, and potential functional overlap with its paralogs. Building on our recent findings that neural-specific lengthened 3' UTR isoforms are co-determined by ELAV/Hu factors, we address their impacts on splicing. While only a few splicing targets of Drosophila are known, ectopic expression of each of the three family members (Elav, Fne and Rbp9) alters hundreds of cassette exon and alternative last exon (ALE) splicing choices. Reciprocally, double mutants of elav/fne, but not elav alone, exhibit opposite effects on both classes of regulated mRNA processing events in larval CNS. While manipulation of Drosophila ELAV/Hu RBPs induces both exon skipping and inclusion, characteristic ELAV/Hu motifs are enriched only within introns flanking exons that are suppressed by ELAV/Hu factors. Moreover, the roles of ELAV/Hu factors in global promotion of distal ALE splicing are mechanistically linked to terminal 3' UTR extensions in neurons, since both processes involve bypass of proximal polyadenylation signals linked to ELAV/Hu motifs downstream of cleavage sites. We corroborate the direct action of Elav in diverse modes of mRNA processing using RRM-dependent Elav-CLIP data from S2 cells. Finally, we provide evidence for conservation in mammalian neurons, which undergo broad programs of distal ALE and APA lengthening, linked to ELAV/Hu motifs downstream of regulated polyadenylation sites. Overall, ELAV/Hu RBPs orchestrate multiple broad programs of neuronal mRNA processing and isoform diversification in Drosophila and mammalian neurons.

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.

Overlapping Activities of ELAV/Hu Family RNA Binding Proteins Specify the Extended Neuronal 3' UTR Landscape in Drosophila

The tissue-specific deployment of highly extended neural 3' UTR isoforms, generated by alternative polyadenylation (APA), is a broad and conserved feature of metazoan genomes. However, the factors and mechanisms that control neural APA isoforms are not well understood. Here, we show that three ELAV/Hu RNA binding proteins (Elav, Rbp9, and Fne) have similar capacities to induce a lengthened 3' UTR landscape in an ectopic setting. These factors promote accumulation of chromatin-associated, 3' UTR-extended, nascent transcripts, through inhibition of proximal polyadenylation site (PAS) usage. Notably, Elav represses an unannotated splice isoform of fne, switching the normally cytoplasmic Fne toward the nucleus in elav mutants. We use genomic profiling to reveal strong and broad loss of neural APA in elav/fne double mutant CNS, the first genetic background to largely abrogate this distinct APA signature. Overall, we demonstrate how regulatory interplay and functionally overlapping activities of neural ELAV/Hu RBPs drives the neural APA landscape.

Statistical inference of the rate of RNA polymerase II elongation by total RNA sequencing

Motivation

Sequencing total RNA without poly-A selection enables us to obtain a transcriptomic profile of nascent RNAs undergoing transcription with co-transcriptional splicing. In general, the RNA-seq reads exhibit a sawtooth pattern in a gene, which is characterized by a monotonically decreasing gradient across introns in the 5’–3’ direction, and by substantially higher levels of RNA-seq reads present in exonic regions. Such patterns result from the process of underlying transcription elongation by RNA polymerase II, which traverses the DNA strand in a 5’–3’ direction as it performs a complex series of mRNA synthesis and processing. Therefore, data of sequenced total RNAs could be utilized to infer the rate of transcription elongation by solving the inverse problem.

Results

Though solving the inverse problem in total RNA-seq has the great potential, statistical methods have not yet been fully developed. We demonstrate what extent the newly developed method can be useful. The objective is to reconstruct the spatial distribution of transcription elongation rates in a gene from a given noisy, sawtooth-like profile. It is necessary to recover the signal source of the elongation rates separately from several types of nuisance factors, such as unobserved modes of co-transcriptionally occurring mRNA splicing, which exert significant influences on the sawtooth shape. The present method was tested using published total RNA-seq data derived from mouse embryonic stem cells. We investigated the spatial characteristics of the estimated elongation rates, focusing especially on the relation to promoter-proximal pausing of RNA polymerase II, nucleosome occupancy and histone modification patterns.

Availability and implementation

A C implementation of PolSter and sample data are available at https://github.com/yoshida-lab/PolSter.

Supplementary information

Supplementary data are available at Bioinformatics online.

Differential Stability of miR-9-5p and miR-9-3p in the Brain Is Determined by Their Unique Cis- and Trans-Acting Elements

microRNAs (miRs) are fundamental regulators of protein coding genes. In the CNS, miR-9 is highly enriched and critical for neuronal development and function. Mature miRs are derived from a duplex precursor, and the -5p strand ("guide") is preferentially incorporated into an RNA-induced silencing complex (RISC) to exert its regulatory functions, while the complementary -3p strand ("passenger") is thought to be rapidly degraded. By contrast, both strands of the miR-9 duplex have unique functions critical for neuronal physiology, yet their respective degradation rates and mechanisms governing degradation are not well understood. Therefore, we determined the degradation kinetics of miR-9-5p and miR-9-3p and investigated the cis and trans elements that affected their stability in the brain. Using a combination of homogeneous neuronal/astrocyte cell models and heterogeneous brain tissue lysate, we demonstrate the novel finding that miR-9-3p was more stable than the miR-9-5p guide strand in all models tested. Moreover, the degradation kinetics of both miR-9-5p and miR-9-3p were brain-region specific, suggesting that each brain region was differentially enriched for specific degradation factors. We also determined that the 3' nucleotides harbor important cis elements required to not only maintain stability, but also to recruit potential protein degradation factors. We used mass spectrometry to assess the miR-9 interacting proteins and found that the -5p and -3p strands were associated with functionally distinct proteins. Overall, these studies revealed unique miR-9-5p and miR-9-3p degradation kinetics in the brain and proposed critical nucleotide sequences and protein partners that could contribute to this differential stability.

YTHDF2 destabilizes m6A-modified neural-specific RNAs to restrain differentiation in induced pluripotent stem cells

N⁶-methyladenosine (m⁶A) is an abundant post-transcriptional modification that can impact RNA fate via interactions with m⁶A-specific RNA binding proteins. Despite accumulating evidence that m⁶A plays an important role in modulating pluripotency, the influence of m⁶A reader proteins in pluripotency is less clear. Here, we report that YTHDF2, an m⁶A reader associated with mRNA degradation, is highly expressed in induced pluripotent stem cells (iPSCs) and down-regulated during neural differentiation. Through RNA sequencing, we identified a group of m⁶A-modified transcripts associated with neural development that are directly regulated by YTDHF2. Depletion of YTHDF2 in iPSCs leads to stabilization of these transcripts, loss of pluripotency, and induction of neural-specific gene expression. Collectively, our results suggest YTHDF2 functions to restrain expression of neural-specific mRNAs in iPSCs and facilitate their rapid and coordinated up-regulation during neural induction. These effects are both achieved by destabilization of the targeted transcripts.

DNA Methylation Regulates Alternative Polyadenylation via CTCF and the Cohesin Complex

Dysregulation of DNA methylation and mRNA alternative cleavage and polyadenylation (APA) are both prevalent in cancer and have been studied as independent processes. We discovered a DNA methylation-regulated APA mechanism when we compared genome-wide DNA methylation and polyadenylation site usage between DNA methylation-competent and DNA methylation-deficient cells. Here, we show that removal of DNA methylation enables CTCF binding and recruitment of the cohesin complex, which, in turn, form chromatin loops that promote proximal polyadenylation site usage. In this DNA demethylated context, either deletion of the CTCF binding site or depletion of RAD21 cohesin complex protein can recover distal polyadenylation site usage. Using data from The Cancer Genome Atlas, we authenticated the relationship between DNA methylation and mRNA polyadenylation isoform expression in vivo. This DNA methylation-regulated APA mechanism demonstrates how aberrant DNA methylation impacts transcriptome diversity and highlights the potential sequelae of global DNA methylation inhibition as a cancer treatment.

A Protocol for Transcriptome-Wide Inference of RNA Metabolic Rates in Mouse Embryonic Stem Cells

The relative ease of mouse Embryonic Stem Cells (mESCs) culture and the potential of these cells to differentiate into any of the three primary germ layers: ectoderm, endoderm and mesoderm (pluripotency), makes them an ideal and frequently used ex vivo system to dissect how gene expression changes impact cell state and differentiation. These efforts are further supported by the large number of constitutive and inducible mESC mutants established with the aim of assessing the contributions of different pathways and genes to cell homeostasis and gene regulation. Gene product abundance is controlled by the modulation of the rates of RNA synthesis, processing, and degradation. The ability to determine the relative contribution of these different RNA metabolic rates to gene expression control using standard RNA-sequencing approaches, which only capture steady state abundance of transcripts, is limited. In contrast, metabolic labeling of RNA with 4-thiouridine (4sU) coupled with RNA-sequencing, allows simultaneous and reproducible inference of transcriptome wide synthesis, processing, and degradation rates. Here we describe, a detailed protocol for 4sU metabolic labeling in mESCs that requires short 4sU labeling times at low concentration and minimally impacts cellular homeostasis. This approach presents a versatile method for in-depth characterization of the gene regulatory strategies governing gene steady state abundance in mESC.

A-Tail of Telomerase RNA Maturation

In this issue of Molecular Cell, Roake et al. (2019) define a feedforward kinetic pathway consisting of a cycle of oligoadenylation and deadenylation that regulates the production of mature human telomerase RNA.

Eukaryotic translation initiation factor 4AII contributes to microRNA-122 regulation of hepatitis C virus replication

Hepatitis C virus (HCV) is a positive sense RNA virus that persistently infects human liver, leading to cirrhosis and hepatocellular carcinoma. HCV replication requires the liver-specific microRNA-122 (miR-122). In contrast to canonical miRNA-mediated repression via 3′UTR sites, miR-122 positively regulates HCV replication by a direct interaction with the 5′ untranslated region (UTR) of the viral RNA. The protein factor requirements for this unusual miRNA regulation remain poorly understood. Here, we identify eIF4AII, previously implicated in miRNA-mediated repression via 3′UTR sites, as a host factor that is important for HCV replication. We demonstrate that eIF4AII interacts with HCV RNA and that this interaction is miR-122-dependent. We show that effective miR-122 binding to, and regulation of, HCV RNA are reduced following eIF4AII depletion. We find that the previously identified HCV co-factor CNOT1, which has also been implicated in miRNA-mediated repression via 3′UTR sites, contributes to regulation of HCV by eIF4AII. Finally, we show that eIF4AI knockdown alleviates the inhibition of HCV replication mediated by depletion of either eIF4AII or CNOT1. Our results suggest a competition effect between the eIF4A proteins to influence HCV replication by modulation of miR-122 function.

P-TEFb Activation by RBM7 Shapes a Pro-survival Transcriptional Response to Genotoxic Stress

DNA damage response (DDR) involves dramatic transcriptional alterations, the mechanisms of which remain ill defined. Here, we show that following genotoxic stress, the RNA-binding motif protein 7 (RBM7) stimulates RNA polymerase II (Pol II) transcription and promotes cell viability by activating the positive transcription elongation factor b (P-TEFb) via its release from the inhibitory 7SK small nuclear ribonucleoprotein (7SK snRNP). This is mediated by activation of p38MAPK, which triggers enhanced binding of RBM7 with core subunits of 7SK snRNP. In turn, P-TEFb relocates to chromatin to induce transcription of short units, including key DDR genes and multiple classes of non-coding RNAs. Critically, interfering with the axis of RBM7 and P-TEFb provokes cellular hypersensitivity to DNA-damage-inducing agents due to activation of apoptosis. Our work uncovers the importance of stress-dependent stimulation of Pol II pause release, which enables a pro-survival transcriptional response that is crucial for cell fate upon genotoxic insult.

A slow transcription rate causes embryonic lethality and perturbs kinetic coupling of neuronal genes

The rate of RNA polymerase II (RNAPII) elongation has an important role in the control of alternative splicing (AS); however, the in vivo consequences of an altered elongation rate are unknown. Here, we generated mouse embryonic stem cells (ESCs) knocked in for a slow elongating form of RNAPII We show that a reduced transcriptional elongation rate results in early embryonic lethality in mice. Focusing on neuronal differentiation as a model, we observed that slow elongation impairs development of the neural lineage from ESCs, which is accompanied by changes in AS and in gene expression along this pathway. In particular, we found a crucial role for RNAPII elongation rate in transcription and splicing of long neuronal genes involved in synapse signaling. The impact of the kinetic coupling of RNAPII elongation rate with AS is greater in ESC-differentiated neurons than in pluripotent cells. Our results demonstrate the requirement for an appropriate transcriptional elongation rate to ensure proper gene expression and to regulate AS during development.

Disruption of Telomerase RNA Maturation Kinetics Precipitates Disease

Mutations in RNA-processing enzymes are increasingly linked to human disease. Telomerase RNA and related noncoding RNAs require 3' end-processing steps, including oligoadenylation. Germline mutations in poly(A)ribonuclease (PARN) cause accumulation of extended human telomerase RNA (hTR) species and precipitate dyskeratosis congenita and pulmonary fibrosis. Here, we develop nascent RNAend-seq to measure processing rates of RNA precursors. We find that mature hTR derives from extended precursors but that in PARN-mutant cells hTR maturation kinetically stalls and unprocessed precursors are degraded. Loss of poly(A)polymerase PAPD5 in PARN-mutant cells accelerates hTR maturation and restores hTR processing, indicating that oligoadenylation and deadenylation set rates of hTR maturation. The H/ACA domain mediates hTR maturation by precisely defining the 3' end, recruiting poly(A)polymerase activity, and conferring sensitivity to PARN regulation. These data reveal a feedforward circuit in which post-transcriptional oligoadenylation controls RNA maturation kinetics. Similar alterations in RNA processing rates may contribute to mechanisms of RNA-based human disease.

Regulation of Intronic Polyadenylation by PCF11 Impacts mRNA Expression of Long Genes

Regulation of cleavage and polyadenylation (CPA) affects gene expression and polyadenylation site (PAS) choice. Here, we report that the CPA and termination factor PCF11 modulates gene expression on the basis of gene size. Although downregulation of PCF11 leads to inhibition of short gene expression, long genes are upregulated because of suppressed intronic polyadenylation (IPA) enriched in large introns. We show that this regulatory scheme, named PCF11-mediated expression regulation through IPA (PEIPA), takes place in cell differentiation, during which downregulation of PCF11 is coupled with upregulation of long genes with functions in cell morphology, adhesion, and migration. PEIPA targets distinct gene sets in different cell contexts with similar rules. Furthermore, PCF11 is autoregulated through a conserved IPA site, the removal of which leads to global activation of PASs close to gene promotors. Therefore, PCF11 uses distinct mechanisms to regulate genes of different sizes, and its autoregulation maintains homeostasis of PAS usage in the cell.

Saccharomyces cerevisiae Metabolic Labeling with 4-thiouracil and the Quantification of Newly Synthesized mRNA As a Proxy for RNA Polymerase II Activity

Global defects in RNA polymerase II transcription might be overlooked by transcriptomic studies analyzing steady-state RNA. Indeed, the global decrease in mRNA synthesis has been shown to be compensated by a simultaneous decrease in mRNA degradation to restore normal steady-state levels. Hence, the genome-wide quantification of mRNA synthesis, independently from mRNA decay, is the best direct reflection of RNA polymerase II transcriptional activity. Here, we discuss a method using non-perturbing metabolic labeling of nascent RNAs in Saccharomyces cerevisiae (S. cerevisiae). Specifically, the cells are cultured for 6 min with a uracil analog, 4-thiouracil, and the labeled newly transcribed RNAs are purified and quantified to determine the synthesis rates of all individual mRNA. Moreover, using labeled Schizosaccharomyces pombe cells as internal standard allows comparing mRNA synthesis in different S. cerevisiae strains. Using this protocol and fitting the data with a dynamic kinetic model, the corresponding mRNA decay rates can be determined.

Suppression of Wnt/β-catenin signaling by EGF receptor is required for hair follicle development

Mice that lack the epidermal growth factor receptor (EGFR) fail to develop a hair coat, but the mechanism responsible for this deficit is not completely understood. Here, we show that EGFR plays a critical role to attenuate wingless-type MMTV integration site family member (Wnt)/β-catenin signaling during postnatal hair follicle development. Genetic ablation of EGFR in mice resulted in increased mitotic activity in matrix cells, apoptosis in hair follicles, and impaired differentiation of epithelial lineages that form hair. EGFR is activated in wild-type hair follicle stem cells marked with SOX9 or NFATc1 and is essential to restrain proliferation and support stem cell numbers and their quiescence. We observed elevated levels of Wnt4, 6, 7b, 10a, 10b, and 16 transcripts and hyperactivation of the β-catenin pathway in EGFR knockout follicles. Using primary keratinocytes, we linked ligand-induced EGFR activation to suppression of nascent mRNA synthesis of Wnt genes. Overexpression of the Wnt antagonist sFRP1 in mice lacking EGFR demonstrated that elevated Wnts are a major cause for the hair follicle defects. Colocalization of transforming growth factor α and Wnts regulated by EGFR in stem cells and progeny indicates that EGFR autocrine loops control Wnts. Our findings define a novel mechanism that integrates EGFR and Wnt/β-catenin pathways to coordinate the delicate balance between proliferation and differentiation during development.

Small RNAs Are Trafficked from the Epididymis to Developing Mammalian Sperm

The biogenesis of the RNA payload of mature sperm is of great interest, because RNAs delivered to the zygote at fertilization can affect early development. Here, we tested the hypothesis that small RNAs are trafficked to mammalian sperm during the process of post-testicular maturation in the epididymis. By characterizing small RNA dynamics during germ cell maturation in mice, we confirm and extend prior observations that sperm undergo a dramatic switch in the RNA payload from piRNAs to tRNA fragments (tRFs) upon exiting the testis and entering the epididymis. Small RNA delivery to sperm could be recapitulated in vitro by incubating testicular spermatozoa with caput epididymosomes. Finally, tissue-specific metabolic labeling of RNAs in intact mice definitively shows that mature sperm carry RNAs that were originally synthesized in the epididymal epithelium. These data demonstrate that soma-germline RNA transfer occurs in male mammals, most likely via vesicular transport from the epididymis to maturing sperm.

Decoupling the impact of microRNAs on translational repression versus RNA degradation in embryonic stem cells

Translation and mRNA degradation are intimately connected, yet the mechanisms that link them are not fully understood. Here, we studied these mechanisms in embryonic stem cells (ESCs). Transcripts showed a wide range of stabilities, which correlated with their relative translation levels and that did not change during early ESC differentiation. The protein DHH1 links translation to mRNA stability in yeast; however, loss of the mammalian homolog, DDX6, in ESCs did not disrupt the correlation across transcripts. Instead, the loss of DDX6 led to upregulated translation of microRNA targets, without concurrent changes in mRNA stability. The Ddx6 knockout cells were phenotypically and molecularly similar to cells lacking all microRNAs (Dgcr8 knockout ESCs). These data show that the loss of DDX6 can separate the two canonical functions of microRNAs: translational repression and transcript destabilization. Furthermore, these data uncover a central role for translational repression independent of transcript destabilization in defining the downstream consequences of microRNA loss.

The Nucleosome Remodeling and Deacetylation Complex Modulates Chromatin Structure at Sites of Active Transcription to Fine-Tune Gene Expression

Chromatin remodeling complexes play essential roles in metazoan development through widespread control of gene expression, but the precise molecular mechanisms by which they do this in vivo remain ill defined. Using an inducible system with fine temporal resolution, we show that the nucleosome remodeling and deacetylation (NuRD) complex controls chromatin architecture and the protein binding repertoire at regulatory regions during cell state transitions. This is primarily exerted through its nucleosome remodeling activity while deacetylation at H3K27 follows changes in gene expression. Additionally, NuRD activity influences association of RNA polymerase II at transcription start sites and subsequent nascent transcript production, thereby guiding the establishment of lineage-appropriate transcriptional programs. These findings provide a detailed molecular picture of genome-wide modulation of lineage-specific transcription by an essential chromatin remodeling complex as well as insight into the orchestration of molecular events involved in transcriptional transitions in vivo.

Video Abstract

•

NuRD increases nucleosome density, expelling TFs or inhibiting recruitment

•

NuRD displaces RNA Pol II from TSSs, reducing nascent transcription

•

Local gains in TF and Mediator occupancy can be indirect effects of NuRD activity

•

Resetting protein binding at regulatory elements can promote or suppress transcription

Systematic Detection of Poly(A)+ RNA-Interacting Proteins and Their Differential Binding

RNA-binding proteins are dynamic posttranscriptional regulators of gene expression. Identification of mRNA-binding proteins in a given experimental setting is thus of great importance. We describe a procedure to enrich for direct poly(A)⁺ RNA protein binders by 4-thiouridine-enhanced UV cross-linking and oligo(dT) purification. Subsequent nuclease-mediated release of RNA-binding proteins (RBPs) from mRNA allows for detection of eluted proteins by mass spectrometry. In addition, we provide a comparative approach to detect differences in RBP binding activity upon a biological stimulus.

Cellular stress alters 3'UTR landscape through alternative polyadenylation and isoform-specific degradation

Most eukaryotic genes express alternative polyadenylation (APA) isoforms with different 3'UTR lengths, production of which is influenced by cellular conditions. Here, we show that arsenic stress elicits global shortening of 3'UTRs through preferential usage of proximal polyadenylation sites during stress and enhanced degradation of long 3'UTR isoforms during recovery. We demonstrate that RNA-binding protein TIA1 preferentially interacts with alternative 3'UTR sequences through U-rich motifs, correlating with stress granule association and mRNA decay of long 3'UTR isoforms. By contrast, genes with shortened 3'UTRs due to stress-induced APA can evade mRNA clearance and maintain transcript abundance post stress. Furthermore, we show that stress causes distinct 3'UTR size changes in proliferating and differentiated cells, highlighting its context-specific impacts on the 3'UTR landscape. Together, our data reveal a global, 3'UTR-based mRNA stability control in stressed cells and indicate that APA can function as an adaptive mechanism to preserve mRNAs in response to stress.

Real-time Analysis of Transcription Factor Binding, Transcription, Translation, and Turnover to Display Global Events During Cellular Activation

Upon activation, cells rapidly change their functional programs and, thereby, their gene expression profile. Massive changes in gene expression occur, for example, during cellular differentiation, morphogenesis, and functional stimulation (such as activation of immune cells), or after exposure to drugs and other factors from the local environment. Depending on the stimulus and cell type, these changes occur rapidly and at any possible level of gene regulation. Displaying all molecular processes of a responding cell to a certain type of stimulus/drug is one of the hardest tasks in molecular biology. Here, we describe a protocol that enables the simultaneous analysis of multiple layers of gene regulation. We compare, in particular, transcription factor binding (Chromatin-immunoprecipitation-sequencing (ChIP-seq)), de novo transcription (4-thiouridine-sequencing (4sU-seq)), mRNA processing, and turnover as well as translation (ribosome profiling). By combining these methods, it is possible to display a detailed and genome-wide course of action.

Sequencing newly transcribed RNA is especially recommended when analyzing rapidly adapting or changing systems, since this depicts the transcriptional activity of all genes during the time of 4sU exposure (irrespective of whether they are up- or downregulated). The combinatorial use of total RNA-seq and ribosome profiling additionally allows the calculation of RNA turnover and translation rates. Bioinformatic analysis of high-throughput sequencing results allows for many means for analysis and interpretation of the data. The generated data also enables tracking co-transcriptional and alternative splicing, just to mention a few possible outcomes.

The combined approach described here can be applied for different model organisms or cell types, including primary cells. Furthermore, we provide detailed protocols for each method used, including quality controls, and discuss potential problems and pitfalls.

Drosha drives the formation of DNA:RNA hybrids around DNA break sites to facilitate DNA repair

The error-free and efficient repair of DNA double-stranded breaks (DSBs) is extremely important for cell survival. RNA has been implicated in the resolution of DNA damage but the mechanism remains poorly understood. Here, we show that miRNA biogenesis enzymes, Drosha and Dicer, control the recruitment of repair factors from multiple pathways to sites of damage. Depletion of Drosha significantly reduces DNA repair by both homologous recombination (HR) and non-homologous end joining (NHEJ). Drosha is required within minutes of break induction, suggesting a central and early role for RNA processing in DNA repair. Sequencing of DNA:RNA hybrids reveals RNA invasion around DNA break sites in a Drosha-dependent manner. Removal of the RNA component of these structures results in impaired repair. These results show how RNA can be a direct and critical mediator of DNA damage repair in human cells.

The kinetics of pre-mRNA splicing in the Drosophila genome and the influence of gene architecture

Production of most eukaryotic mRNAs requires splicing of introns from pre-mRNA. The splicing reaction requires definition of splice sites, which are initially recognized in either intron-spanning (‘intron definition’) or exon-spanning (‘exon definition’) pairs. To understand how exon and intron length and splice site recognition mode impact splicing, we measured splicing rates genome-wide in Drosophila, using metabolic labeling/RNA sequencing and new mathematical models to estimate rates. We found that the modal intron length range of 60–70 nt represents a local maximum of splicing rates, but that much longer exon-defined introns are spliced even faster and more accurately. We observed unexpectedly low variation in splicing rates across introns in the same gene, suggesting the presence of gene-level influences, and we identified multiple gene level variables associated with splicing rate. Together our data suggest that developmental and stress response genes may have preferentially evolved exon definition in order to enhance the rate or accuracy of splicing.

mRNA pseudouridylation affects RNA metabolism in the parasite Toxoplasma gondii

RNA contains over 100 modified nucleotides that are created post-transcriptionally, among which pseudouridine (Ψ) is one of the most abundant. Although it was one of the first modifications discovered, the biological role of this modification is still not fully understood. Recently, we reported that a pseudouridine synthase (TgPUS1) is necessary for differentiation of the single-celled eukaryotic parasite Toxoplasma gondii from active to chronic infection. To better understand the biological role of pseudouridylation, we report here gel-based and deep-sequencing methods to identify TgPUS1-dependent Ψ's in Toxoplasma RNA, and the use of TgPUS1 mutants to examine the effect of this modification on mRNAs. In addition to identifying conserved sites of pseudouridylation in Toxoplasma rRNA, tRNA, and snRNA, we also report extensive pseudouridylation of Toxoplasma mRNAs, with the Ψ's being relatively depleted in the 3'-UTR but enriched at position 1 of codons. We show that many Ψ's in tRNA and mRNA are dependent on the action of TgPUS1 and that TgPUS1-dependent mRNA Ψ's are enriched in developmentally regulated transcripts. RNA-seq data obtained from wild-type and TgPUS1-mutant parasites shows that genes containing a TgPUS1-dependent Ψ are relatively more abundant in mutant parasites, while pulse/chase labeling of RNA with 4-thiouracil shows that mRNAs containing TgPUS1-dependent Ψ have a modest but statistically significant increase in half-life in the mutant parasites. These data are some of the first evidence suggesting that mRNA Ψ's play an important biological role.

BET-Bromodomain Inhibitors Engage the Host Immune System and Regulate Expression of the Immune Checkpoint Ligand PD-L1

BET inhibitors (BETi) target bromodomain-containing proteins and are currently being evaluated as anti-cancer agents. We find that maximal therapeutic effects of BETi in a Myc-driven B cell lymphoma model required an intact host immune system. Genome-wide analysis of the BETi-induced transcriptional response identified the immune checkpoint ligand Cd274 (Pd-l1) as a Myc-independent, BETi target-gene. BETi directly repressed constitutively expressed and interferon-gamma (IFN-γ) induced CD274 expression across different human and mouse tumor cell lines and primary patient samples. Mechanistically, BETi decreased Brd4 occupancy at the Cd274 locus without any change in Myc occupancy, resulting in transcriptional pausing and rapid loss of Cd274 mRNA production. Finally, targeted inhibition of the PD-1/PD-L1 axis by combining anti-PD-1 antibodies and the BETi JQ1 caused synergistic responses in mice bearing Myc-driven lymphomas. Our data uncover an interaction between BETi and the PD-1/PD-L1 immune-checkpoint and provide mechanistic insight into the transcriptional regulation of CD274.

•

BETi require an intact host immune system to promote robust anti-tumor responses

•

BRD4 inhibition inhibits PD-L1 transcription independently from MYC expression

•

BRD4 and IRF1 co-regulate interferon-induced PD-L1 transcription

•

Combinations of BET inhibitor and immune modulating therapy are efficacious in vivo

Identification of RNA-binding domains of RNA-binding proteins in cultured cells on a system-wide scale with RBDmap

This protocol is an extension to: Nat. Protoc. 8, 491-500 (2013); doi:10.1038/nprot.2013.020; published online 14 February 2013RBDmap is a method for identifying, in a proteome-wide manner, the regions of RNA-binding proteins (RBPs) engaged in native interactions with RNA. In brief, cells are irradiated with UV light to induce protein-RNA cross-links. Following stringent denaturing washes, the resulting covalently linked protein-RNA complexes are purified with oligo(dT) magnetic beads. After elution, RBPs are subjected to partial proteolysis, in which the protein regions still bound to the RNA and those released to the supernatant are separated by a second oligo(dT) selection. After sample preparation and mass-spectrometric analysis, peptide intensity ratios between the RNA-bound and released fractions are used to determine the RNA-binding regions. As a Protocol Extension, this article describes an adaptation of an existing Protocol and offers additional applications. The earlier protocol (for the RNA interactome capture method) describes how to identify the active RBPs in cultured cells, whereas this Protocol Extension also enables the identification of the RNA-binding domains of RBPs. The experimental workflow takes 1 week plus 2 additional weeks for proteomics and data analysis. Notably, RBDmap presents numerous advantages over classic methods for determining RNA-binding domains: it produces proteome-wide, high-resolution maps of the protein regions contacting the RNA in a physiological context and can be adapted to different biological systems and conditions. Because RBDmap relies on the isolation of polyadenylated RNA via oligo(dT), it will not provide RNA-binding information on proteins interacting exclusively with nonpolyadenylated transcripts. Applied to HeLa cells, RBDmap uncovered 1,174 RNA-binding sites in 529 proteins, many of which were previously unknown.

Widespread activation of antisense transcription of the host genome during herpes simplex virus 1 infection

Background

Herpesviruses can infect a wide range of animal species. Herpes simplex virus 1 (HSV-1) is one of the eight herpesviruses that can infect humans and is prevalent worldwide. Herpesviruses have evolved multiple ways to adapt the infected cells to their needs, but knowledge about these transcriptional and post-transcriptional modifications is sparse.

Results

Here, we show that HSV-1 induces the expression of about 1000 antisense transcripts from the human host cell genome. A subset of these is also activated by the closely related varicella zoster virus. Antisense transcripts originate either at gene promoters or within the gene body, and they show different susceptibility to the inhibition of early and immediate early viral gene expression. Overexpression of the major viral transcription factor ICP4 is sufficient to turn on a subset of antisense transcripts. Histone marks around transcription start sites of HSV-1-induced and constitutively transcribed antisense transcripts are highly similar, indicating that the genetic loci are already poised to transcribe these novel RNAs. Furthermore, an antisense transcript overlapping with the BBC3 gene (also known as PUMA) transcriptionally silences this potent inducer of apoptosis in cis.

Conclusions

We show for the first time that a virus induces widespread antisense transcription of the host cell genome. We provide evidence that HSV-1 uses this to downregulate a strong inducer of apoptosis. Our findings open new perspectives on global and specific alterations of host cell transcription by viruses.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1329-5) contains supplementary material, which is available to authorized users.

DDX54 regulates transcriptome dynamics during DNA damage response

The cellular response to genotoxic stress is mediated by a well-characterized network of DNA surveillance pathways. The contribution of post-transcriptional gene regulatory networks to the DNA damage response (DDR) has not been extensively studied. Here, we systematically identified RNA-binding proteins differentially interacting with polyadenylated transcripts upon exposure of human breast carcinoma cells to ionizing radiation (IR). Interestingly, more than 260 proteins, including many nucleolar proteins, showed increased binding to poly(A)+ RNA in IR-exposed cells. The functional analysis of DDX54, a candidate genotoxic stress responsive RNA helicase, revealed that this protein is an immediate-to-early DDR regulator required for the splicing efficacy of its target IR-induced pre-mRNAs. Upon IR exposure, DDX54 acts by increased interaction with a well-defined class of pre-mRNAs that harbor introns with weak acceptor splice sites, as well as by protein-protein contacts within components of U2 snRNP and spliceosomal B complex, resulting in lower intron retention and higher processing rates of its target transcripts. Because DDX54 promotes survival after exposure to IR, its expression and/or mutation rate may impact DDR-related pathologies. Our work indicates the relevance of many uncharacterized RBPs potentially involved in the DDR.

Widespread Influence of 3'-End Structures on Mammalian mRNA Processing and Stability

The physiological relevance of structures within mammalian mRNAs has been elusive, as these mRNAs are less folded in cells than in vitro and have predicted secondary structures no more stable than those of random sequences. Here, we investigate the possibility that mRNA structures facilitate the 3'-end processing of thousands of human mRNAs by juxtaposing poly(A) signals (PASs) and cleavage sites that are otherwise too far apart. We find that RNA structures are predicted to be more prevalent within these extended 3'-end regions than within PAS-upstream regions and indeed are substantially more folded within cells, as determined by intracellular probing. Analyses of thousands of ectopically expressed variants demonstrate that this folding both enhances processing and increases mRNA metabolic stability. Even folds with predicted stabilities resembling those of random sequences can enhance processing. Structure-controlled processing can also regulate neighboring gene expression. Thus, RNA structure has widespread roles in mammalian mRNA biogenesis and metabolism.

Metabolic labeling and recovery of nascent RNA to accurately quantify mRNA stability

Changes in the rate of mRNA decay are closely coordinated with transcriptional changes and together these events have profound effects on gene expression during development and disease. Traditional approaches to assess mRNA decay have relied on inhibition of transcription, which can alter mRNA decay rates and confound interpretation. More recently, metabolic labeling combined with chemical modification and fractionation of labeled RNAs has allowed the isolation of nascent transcripts and the subsequent calculation of mRNA decay rates. This approach has been widely adopted for measuring mRNA half-lives on a global scale, but has proven challenging to use for analysis of single genes. We present a series of normalization and quality assurance steps to be used in combination with 4-thiouridine pulse labeling of cultured eukaryotic cells. Importantly, we demonstrate how the relative amount of 4sU-labeled nascent RNA influences accurate quantification. The approach described facilitates reproducible measurement of individual mRNA half-lives using 4-thiouridine and could be adapted for use with other nucleoside analogs.

The CELF1 RNA-Binding Protein Regulates Decay of Signal Recognition Particle mRNAs and Limits Secretion in Mouse Myoblasts

We previously identified several mRNAs encoding components of the secretory pathway, including signal recognition particle (SRP) subunit mRNAs, among transcripts associated with the RNA-binding protein CELF1. Through immunoprecipitation of RNAs crosslinked to CELF1 in myoblasts and in vitro binding assays using recombinant CELF1, we now provide evidence that CELF1 directly binds the mRNAs encoding each of the subunits of the SRP. Furthermore, we determined the half-lives of the Srp transcripts in control and CELF1 knockdown myoblasts. Our results indicate CELF1 is a destabilizer of at least five of the six Srp transcripts and that the relative abundance of the SRP proteins is out of balance when CELF1 is depleted. CELF1 knockdown myoblasts exhibit altered secretion of a luciferase reporter protein and are impaired in their ability to migrate and close a wound, consistent with a defect in the secreted extracellular matrix. Importantly, similar defects in wound healing are observed when SRP subunit imbalance is induced by over-expression of SRP68. Our studies support the existence of an RNA regulon containing Srp mRNAs that is controlled by CELF1. One implication is that altered function of CELF1 in myotonic dystrophy may contribute to changes in the extracellular matrix of affected muscle through defects in secretion.