Microarray analysis of newly synthesized RNA in cells and animals

Three-way junction structure-mediated reverse transcription-free exponential amplification reaction for pathogen RNA detection

Since RNA is an important biomarker of many infectious pathogens, RNA detection of pathogenic organisms is crucial for disease diagnosis and environmental and food safety. By simulating the base mismatch during DNA replication, this study presents a novel three-way junction structure-mediated reverse transcription-free exponential amplification reaction (3WJ-RTF-EXPAR) for the rapid and sensitive detection of pathogen RNA. The target RNA served as a switch to initiate the reaction by forming a three-way junction (3WJ) structure with the ex-trigger strand and the ex-primer strand. The generated trigger strand could be significantly amplified through EXPAR to open the stem-loop structure of the molecular beacon to emit fluorescence signal. The proofreading activity of Vent DNA polymerase, in combination with the unique structure of 2+1 bases at the 3'-end of the ex-primer strand, could enhance the role of target RNA as a reaction switch to reduce non-specific amplification and ensure excellent specificity to differentiate target pathogen from those causing similar symptoms. Furthermore, detection of target RNA showed a detection limit of 1.0×104 copies/mL, while the time consumption was only 20 min, outperforming qRT-LAMP and qRT-PCR, the most commonly used RNA detection methods in clinical practice. All those indicates the great application prospects of this method in clinical diagnostic.

Engineering of a Chimeric Template Triggers RNase H-Based Isothermal Amplification Approach for Sensitive Detection of Pathogen RNA

RNA-based detection of pathogenic organisms is an emerging field of research that is crucial for disease diagnosis and environmental and food safety. By rationally engineering an RNA-DNA tandem (RDT) structural template, we proposed a novel RNase H-based isothermal exponential amplification (RH-IEA) reaction to rapidly identify long-stranded RNA. In this strategy, the rigid and compact RDT template selectively recognized the target RNA and formed a stable hybrid with it. Upon site-specific cleavage of RNase H, the 3' overhang of the target RNA was cut off, and a free hydroxyl end at the hydrolysis site was generated to trigger an exponential amplification reaction (EXPAR). This method maintained the high efficiency and rapid amplification kinetics of EXPAR. As a result, the RH-IEA strategy was able to sensitively and specifically detect the characteristic sequence of Escherichia coli O157:H7 RNA, with a detection sensitivity of 1 fg/μL. Besides, the RDT template can be used as an RNA protector to prevent specific segments of the target RNA from being degraded by RNase enzymes, allowing the sample to be stored at room temperature for a long time. With this advantage, the practicality of RH-IEA will be more flexible than the reverse transcription polymerase chain reaction. It was successfully applied in the identification of E. coli O157:H7 in milk with a minimum detection concentration of 1.0 × 102 CFU/mL. Therefore, the RH-IEA method will serve as a powerful tool for detecting long-stranded RNA and will also shed light on the pathogen detection in food safety and molecular diagnosis.

Genome-wide Mapping of 5'-monophosphorylated Ends of Mammalian Nascent RNA Transcripts

In eukaryotic cells, RNA biogenesis generally requires processing of the nascent transcript as it is being synthesized by RNA polymerase. These processing events include endonucleolytic cleavage, exonucleolytic trimming, and splicing of the growing nascent transcript. Endonucleolytic cleavage events that generate an exposed 5'-monophosphorylated (5'-PO4) end on the growing nascent transcript occur in the maturation of rRNAs, tRNAs, and mRNAs. These 5'-PO4 ends can be a target of further processing or be subjected to 5'-3' exonucleolytic digestion that may result in termination of transcription. Here, we describe how to identify 5'-PO4 ends of intermediates in nascent RNA metabolism. We capture these species via metabolic labeling with bromouridine followed by immunoprecipitation and specific ligation of 5'-PO4 RNA ends with the 3'-hydroxyl group of a 5' adaptor (5'-PO4 Bru-Seq) using RNA ligase I. These ligation events are localized at single nucleotide resolution via highthroughput sequencing, which identifies the position of 5'-PO4 groups precisely. This protocol successfully detects the 5'monophosphorylated ends of RNA processing intermediates during production of mature ribosomal, transfer, and micro RNAs. When combined with inhibition of the nuclear 5'-3' exonuclease Xrn2, 5'-PO4 Bru-Seq maps the 5' splice sites of debranched introns and mRNA and tRNA 3' end processing sites cleaved by CPSF73 and RNaseZ, respectively. Key features • Metabolic labeling for brief periods with bromouridine focuses the analysis of 5'-PO4 RNA ends on the population of nascent transcripts that are actively transcribed. • Detects 5'-PO4 RNA ends on nascent transcripts produced by all RNA polymerases. • Detects 5'-PO4 RNA ends at single nucleotide resolution.

Improving the study of RNA dynamics through advances in RNA-seq with metabolic labeling and nucleotide-recoding chemistry

RNA metabolic labeling using 4-thiouridine (s4U) captures the dynamics of RNA synthesis and decay. The power of this approach is dependent on appropriate quantification of labeled and unlabeled sequencing reads, which can be compromised by the apparent loss of s4U-labeled reads in a process we refer to as dropout. Here we show that s4U-containing transcripts can be selectively lost when RNA samples are handled under sub-optimal conditions, but that this loss can be minimized using an optimized protocol. We demonstrate a second cause of dropout in nucleotide recoding and RNA sequencing (NR-seq) experiments that is computational and downstream of library preparation. NR-seq experiments involve chemically converting s4U from a uridine analog to a cytidine analog and using the apparent T-to-C mutations to identify the populations of newly synthesized RNA. We show that high levels of T-to-C mutations can prevent read alignment with some computational pipelines, but that this bias can be overcome using improved alignment pipelines. Importantly, kinetic parameter estimates are affected by dropout independent of the NR chemistry employed, and all chemistries are practically indistinguishable in bulk, short-read RNA-seq experiments. Dropout is an avoidable problem that can be identified by including unlabeled controls, and mitigated through improved sample handing and read alignment that together improve the robustness and reproducibility of NR-seq experiments.

Profiling transcription factor activity dynamics using intronic reads in time-series transcriptome data

Activities of transcription factors (TFs) are temporally modulated to regulate dynamic cellular processes, including development, homeostasis, and disease. Recent developments of bioinformatic tools have enabled the analysis of TF activities using transcriptome data. However, because these methods typically use exon-based target expression levels, the estimated TF activities have limited temporal accuracy. To address this, we proposed a TF activity measure based on intron-level information in time-series RNA-seq data, and implemented it to decode the temporal control of TF activities during dynamic processes. We showed that TF activities inferred from intronic reads can better recapitulate instantaneous TF activities compared to the exon-based measure. By analyzing public and our own time-series transcriptome data, we found that intron-based TF activities improve the characterization of temporal phasing of cycling TFs during circadian rhythm, and facilitate the discovery of two temporally opposing TF modules during T cell activation. Collectively, we anticipate that the proposed approach would be broadly applicable for decoding global transcriptional architecture during dynamic processes.

Interrogating the transcriptome with metabolically incorporated ribonucleosides

RNA is a central player in biological processes, but there remain major gaps in our understanding of transcriptomic processes and the underlying biochemical mechanisms regulating RNA in cells. A powerful strategy to facilitate molecular analysis of cellular RNA is the metabolic incorporation of chemical probes. In this review, we discuss current approaches for RNA metabolic labeling with modified ribonucleosides and their integration with Next-Generation Sequencing, mass spectrometry-based proteomics, and fluorescence microscopy in order to interrogate RNA behavior in its native context.

A comparison of metabolic labeling and statistical methods to infer genome-wide dynamics of RNA turnover

Metabolic labeling of newly transcribed RNAs coupled with RNA-seq is being increasingly used for genome-wide analysis of RNA dynamics. Methods including standard biochemical enrichment and recent nucleotide conversion protocols each require special experimental and computational treatment. Despite their immediate relevance, these technologies have not yet been assessed and benchmarked, and no data are currently available to advance reproducible research and the development of better inference tools. Here, we present a systematic evaluation and comparison of four RNA labeling protocols: 4sU-tagging biochemical enrichment, including spike-in RNA controls, SLAM-seq, TimeLapse-seq and TUC-seq. All protocols are evaluated based on practical considerations, conversion efficiency and wet lab requirements to handle hazardous substances. We also compute decay rate estimates and confidence intervals for each protocol using two alternative statistical frameworks, pulseR and GRAND-SLAM, for over 11 600 human genes and evaluate the underlying computational workflows for their robustness and ease of use. Overall, we demonstrate a high inter-method reliability across eight use case scenarios. Our results and data will facilitate reproducible research and serve as a resource contributing to a fuller understanding of RNA biology.

Characterizing RNA stability genome-wide through combined analysis of PRO-seq and RNA-seq data

BACKGROUND

The concentrations of distinct types of RNA in cells result from a dynamic equilibrium between RNA synthesis and decay. Despite the critical importance of RNA decay rates, current approaches for measuring them are generally labor-intensive, limited in sensitivity, and/or disruptive to normal cellular processes. Here, we introduce a simple method for estimating relative RNA half-lives that is based on two standard and widely available high-throughput assays: Precision Run-On sequencing (PRO-seq) and RNA sequencing (RNA-seq).

RESULTS

Our method treats PRO-seq as a measure of transcription rate and RNA-seq as a measure of RNA concentration, and estimates the rate of RNA decay required for a steady-state equilibrium. We show that this approach can be used to assay relative RNA half-lives genome-wide, with good accuracy and sensitivity for both coding and noncoding transcription units. Using a structural equation model (SEM), we test several features of transcription units, nearby DNA sequences, and nearby epigenomic marks for associations with RNA stability after controlling for their effects on transcription. We find that RNA splicing-related features are positively correlated with RNA stability, whereas features related to miRNA binding and DNA methylation are negatively correlated with RNA stability. Furthermore, we find that a measure based on U1 binding and polyadenylation sites distinguishes between unstable noncoding and stable coding transcripts but is not predictive of relative stability within the mRNA or lincRNA classes. We also identify several histone modifications that are associated with RNA stability.

CONCLUSION

We introduce an approach for estimating the relative half-lives of individual RNAs. Together, our estimation method and systematic analysis shed light on the pervasive impacts of RNA stability on cellular RNA concentrations.

Nascent RNA analyses: tracking transcription and its regulation

The programmes that direct an organism's development and maintenance are encoded in its genome. Decoding of this information begins with regulated transcription of genomic DNA into RNA. Although transcription and its control can be tracked indirectly by measuring stable RNAs, it is only by directly measuring nascent RNAs that the immediate regulatory changes in response to developmental, environmental, disease and metabolic signals are revealed. Multiple complementary methods have been developed to quantitatively track nascent transcription genome-wide at nucleotide resolution, all of which have contributed novel insights into the mechanisms of gene regulation and transcription-coupled RNA processing. Here we critically evaluate the array of strategies used for investigating nascent transcription and discuss the recent conceptual advances they have provided.

Chemo-enzymatic treatment of RNA to facilitate analyses

Labeling RNA is a recurring problem to make RNA compatible with state-of-the-art methodology and comes in many flavors. Considering only cellular applications, the spectrum still ranges from site-specific labeling of individual transcripts, for example, for live-cell imaging of mRNA trafficking, to metabolic labeling in combination with next generation sequencing to capture dynamic aspects of RNA metabolism on a transcriptome-wide scale. Combining the specificity of RNA-modifying enzymes with non-natural substrates has emerged as a valuable strategy to modify RNA site- or sequence-specifically with functional groups suitable for subsequent bioorthogonal reactions and thus label RNA with reporter moieties such as affinity or fluorescent tags. In this review article, we will cover chemo-enzymatic approaches (a) for in vitro labeling of RNA for application in cells, (b) for treatment of total RNA, and (c) for metabolic labeling of RNA. This article is categorized under: RNA Processing < RNA Editing and Modification RNA Methods < RNA Analyses in vitro and In Silico RNA Methods < RNA Analyses in Cells.

Gaining insight into transcriptome-wide RNA population dynamics through the chemistry of 4-thiouridine

Cellular RNA levels are the result of a juggling act between RNA transcription, processing, and degradation. By tuning one or more of these parameters, cells can rapidly alter the available pool of transcripts in response to stimuli. While RNA sequencing (RNA-seq) is a vital method to quantify RNA levels genome-wide, it is unable to capture the dynamics of different RNA populations at steady-state or distinguish between different mechanisms that induce changes to the steady-state (i.e., altered rate of transcription vs. degradation). The dynamics of different RNA populations can be studied by targeted incorporation of noncanonical nucleosides. 4-Thiouridine (s⁴ U) is a commonly used and versatile RNA metabolic label that allows the study of many properties of RNA metabolism from synthesis to degradation. Numerous experimental strategies have been developed that leverage the power of s⁴ U to label newly transcribed RNA in whole cells, followed by enrichment with activated disulfides or chemistry to induce C mutations at sites of s⁴ U during sequencing. This review presents existing methods to study RNA population dynamics genome-wide using s⁴ U metabolic labeling, as well as a discussion of considerations and challenges when designing s⁴ U metabolic labeling experiments. This article is categorized under: RNA Methods > RNA Analyses in Cells RNA Turnover and Surveillance > Regulation of RNA Stability.

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.

Contributions of regulated transcription and mRNA decay to the dynamics of gene expression

Organisms have acquired sophisticated regulatory networks that control gene expression in response to cellular perturbations. Understanding of the mechanisms underlying the coordinated changes in gene expression in response to external and internal stimuli is a fundamental issue in biology. Recent advances in high-throughput technologies have enabled the measurement of diverse biological information, including gene expression levels, kinetics of gene expression, and interactions among gene expression regulatory molecules. By coupling these technologies with quantitative modeling, we can now uncover the biological roles and mechanisms of gene regulation at the system level. This review consists of two parts. First, we focus on the methods using uridine analogs that measure synthesis and decay rates of RNAs, which demonstrate how cells dynamically change the regulation of gene expression in response to both internal and external cues. Second, we discuss the underlying mechanisms of these changes in kinetics, including the functions of transcription factors and RNA-binding proteins. Overall, this review will help to clarify a system-level view of gene expression programs in cells. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Turnover and Surveillance > Regulation of RNA Stability RNA Methods > RNA Analyses in vitro and In Silico.

Genome-wide real-time in vivo transcriptional dynamics during Plasmodium falciparum blood-stage development

Genome-wide analysis of transcription in the malaria parasite Plasmodium falciparum has revealed robust variation in steady-state mRNA abundance throughout the 48-h intraerythrocytic developmental cycle (IDC), suggesting that this process is highly dynamic and tightly regulated. Here, we utilize rapid 4-thiouracil (4-TU) incorporation via pyrimidine salvage to specifically label, capture, and quantify newly-synthesized RNA transcripts at every hour throughout the IDC. This high-resolution global analysis of the transcriptome captures the timing and rate of transcription for each newly synthesized mRNA in vivo, revealing active transcription throughout all IDC stages. Using a statistical model to predict the mRNA dynamics contributing to the total mRNA abundance at each timepoint, we find varying degrees of transcription and stabilization for each mRNA corresponding to developmental transitions. Finally, our results provide new insight into co-regulation of mRNAs throughout the IDC through regulatory DNA sequence motifs, thereby expanding our understanding of P. falciparum mRNA dynamics.

Transcriptomic analysis often doesn’t differentiate between newly synthesized and stabilized mRNAs. Using rapid 4-thiouracil incorporation, Painter et al. here define genome-wide active transcription throughout Plasmodium blood-stage developmental stages and identify associated regulatory DNA sequence motifs.

EC-tagging allows cell type-specific RNA analysis

Purification of cell type-specific RNAs remains a significant challenge. One solution involves biosynthetic tagging of target RNAs. RNA tagging via incorporation of 4-thiouracil (TU) in cells expressing transgenic uracil phosphoribosyltransferase (UPRT), a method known as TU-tagging, has been used in multiple systems but can have limited specificity due to endogenous pathways of TU incorporation. Here, we describe an alternative method that requires the activity of two enzymes: cytosine deaminase (CD) and UPRT. We found that the sequential activity of these enzymes converts 5-ethynylcytosine (EC) to 5-ethynyluridine monophosphate that is subsequently incorporated into nascent RNAs. The ethynyl group allows efficient detection and purification of tagged RNAs. We show that ‘EC-tagging’ occurs in tissue culture cells and Drosophila engineered to express CD and UPRT. Additional control can be achieved through a split-CD approach in which functional CD is reconstituted from independently expressed fragments. We demonstrate the sensitivity and specificity of EC-tagging by obtaining cell type-specific gene expression data from intact Drosophila larvae, including transcriptome measurements from a small population of central brain neurons. EC-tagging provides several advantages over existing techniques and should be broadly useful for investigating the role of differential RNA expression in cell identity, physiology and pathology.

Capturing in vivo RNA transcriptional dynamics from the malaria parasite Plasmodium falciparum

To capture the transcriptional dynamics within proliferating cells, methods to differentiate nascent transcription from preexisting mRNAs are desired. One approach is to label newly synthesized mRNA transcripts in vivo through the incorporation of modified pyrimidines. However, the human malaria parasite, Plasmodium falciparum, is incapable of pyrimidine salvage for mRNA biogenesis. To capture cellular mRNA dynamics during Plasmodium development, we engineered parasites that can salvage pyrimidines through the expression of a single bifunctional yeast fusion gene, cytosine deaminase/uracil phosphoribosyltransferase (FCU). We show that expression of FCU allows for the direct incorporation of thiol-modified pyrimidines into nascent mRNAs. Using developmental stage-specific promoters to express FCU-GFP enables the biosynthetic capture and in-depth analysis of mRNA dynamics from subpopulations of cells undergoing differentiation. We demonstrate the utility of this method by examining the transcriptional dynamics of the sexual gametocyte stage transition, a process that is essential to malaria transmission between hosts. Using the pfs16 gametocyte-specific promoter to express FCU-GFP in 3D7 parasites, we found that sexual stage commitment is governed by transcriptional reprogramming and stabilization of a subset of essential gametocyte transcripts. We also measured mRNA dynamics in F12 gametocyte-deficient parasites and demonstrate that the transcriptional program required for sexual commitment and maturation is initiated but likely aborted due to the absence of the PfAP2-G transcriptional regulator and a lack of gametocyte-specific mRNA stabilization. Biosynthetic labeling of Plasmodium mRNAs is incredibly versatile, can be used to measure transcriptional dynamics at any stage of parasite development, and will allow for future applications to comprehensively measure RNA-protein interactions in the malaria parasite.

The regulation of transcriptional repression in hypoxia

A sufficient supply molecular oxygen is essential for the maintenance of physiologic metabolism and bioenergetic homeostasis for most metazoans. For this reason, mechanisms have evolved for eukaryotic cells to adapt to conditions where oxygen demand exceeds supply (hypoxia). These mechanisms rely on the modification of pre-existing proteins, translational arrest and transcriptional changes. The hypoxia inducible factor (HIF; a master regulator of gene induction in response to hypoxia) is responsible for the majority of induced gene expression in hypoxia. However, much less is known about the mechanism(s) responsible for gene repression, an essential part of the adaptive transcriptional response. Hypoxia-induced gene repression leads to a reduction in energy demanding processes and the redirection of limited energetic resources to essential housekeeping functions. Recent developments have underscored the importance of transcriptional repressors in cellular adaptation to hypoxia. To date, at least ten distinct transcriptional repressors have been reported to demonstrate sensitivity to hypoxia. Central among these is the Repressor Element-1 Silencing Transcription factor (REST), which regulates over 200 genes. In this review, written to honor the memory and outstanding scientific legacy of Lorenz Poellinger, we provide an overview of our existing knowledge with respect to transcriptional repressors and their target genes in hypoxia.

High-Resolution Gene Expression Profiling of RNA Synthesis, Processing, and Decay by Metabolic Labeling of Newly Transcribed RNA Using 4-Thiouridine

Cellular RNA levels are orchestrated by highly regulated processes involving RNA synthesis (transcription), processing (e.g., splicing, polyadenylation, transport), and degradation. Profiling these changes provides valuable information on the regulation of gene expression. Total cellular RNA is a poor template for revealing short-term changes in gene expression, alterations in RNA decay rates, and the kinetics of RNA processing as well as the differentiation thereof. Here, we describe the metabolic labeling and purification of newly transcribed RNA with 4-thiouridine, by which these limitations are overcome.

A Nonradioactive Assay to Measure Production and Processing of Ribosomal RNA by 4sU-Tagging

In vivo metabolic pulse labeling is a classical approach to assess production and processing of ribosomal RNA (rRNA). However, conventional labeling techniques can be indirect and require work with radioactivity. Here, we describe in detail a protocol for in vivo metabolic labeling, purification, and readout of nascent rRNA by 4-thiouridine (4sU). We propose 4sU labeling as standard nonradioactive technique for the analysis of rRNA metabolism during ribosome biogenesis.

Temporally and spatially restricted gene expression profiling

Identifying gene function in specific cells is critical for understanding the processes that make cells unique. Several different methods are available to isolate actively transcribed RNA or actively translated RNA in specific cells at a chosen time point. Cell-specific mRNA isolation can be accomplished by the expression of transgenes in cells of interest, either directly from a specific promoter or using a modular system such as Gal4/UAS or Cre/lox. All of the methods described in this review, namely thiol-labeling of RNA (TU-tagging or RABT), TRAP (translating ribosome affinity purification) and INTACT (isolation of nuclei tagged in specific cell types), allow next generation sequencing, permitting the identification of enriched gene transcripts within the specific cell-type. We describe here the general concept of each method, include examples, evaluate possible problems related to each technique, and suggest the types of questions for which each method is best suited.

Basonuclin-1 modulates epithelial plasticity and TGF-β1-induced loss of epithelial cell integrity

Transforming growth factor-β1 (TGF-β1) is a multifunctional cytokine and critically involved in the progression of a variety of cancers. TGF-β1 signaling can impair tumor development by its anti-proliferative and pro-apoptotic features. In contrast, it may actively promote tumor progression and cancer cell dissemination by inducing a gradual switch from epithelial towards mesenchymal-like cell features (EMT-like), including decreased intercellular adhesion. Here, we show that expression of the transcription factor Basonuclin-1 (Bnc1) modulates TGF-β1-induced epithelial dedifferentiation of mammary epithelial cells. RNAi-mediated repression of Bnc1 resulted in enhanced intercellular adhesion and strongly impaired TGF-β1-dependent sheet disintegration and cell scattering. In contrast, forced expression of Bnc1 modifies plasma membrane/cytoskeletal dynamics and seemingly interferes with the initiation of sustainable cell-cell contacts. Follow-up analyses revealed that Bnc1 affects the expression of numerous TGF-β1-responsive genes including distinct EMT-related transcription factors, some of which modulate the expression of Bnc1 themselves. These results suggest that Bnc1 is part of a transcription factor network related to epithelial plasticity with reciprocal feedback-loop connections on which Smad-factors integrate TGF-β1 signaling. Our study demonstrates that Bnc1 regulates epithelial plasticity of mammary epithelial cells and influences outcome of TGF-β1 signaling.

SnapShot-Seq: a method for extracting genome-wide, in vivo mRNA dynamics from a single total RNA sample

mRNA synthesis, processing, and destruction involve a complex series of molecular steps that are incompletely understood. Because the RNA intermediates in each of these steps have finite lifetimes, extensive mechanistic and dynamical information is encoded in total cellular RNA. Here we report the development of SnapShot-Seq, a set of computational methods that allow the determination of in vivo rates of pre-mRNA synthesis, splicing, intron degradation, and mRNA decay from a single RNA-Seq snapshot of total cellular RNA. SnapShot-Seq can detect in vivo changes in the rates of specific steps of splicing, and it provides genome-wide estimates of pre-mRNA synthesis rates comparable to those obtained via labeling of newly synthesized RNA. We used SnapShot-Seq to investigate the origins of the intrinsic bimodality of metazoan gene expression levels, and our results suggest that this bimodality is partly due to spillover of transcriptional activation from highly expressed genes to their poorly expressed neighbors. SnapShot-Seq dramatically expands the information obtainable from a standard RNA-Seq experiment.

Applying thiouracil tagging to mouse transcriptome analysis

Transcriptional profiling is a powerful approach for studying mouse development, physiology and disease models. Here we describe a protocol for mouse thiouracil tagging (TU tagging), a transcriptome analysis technology that includes in vivo covalent labeling, purification and analysis of cell type-specific RNA. TU tagging enables the isolation of RNA from a given cell population of a complex tissue, avoiding transcriptional changes induced by cell isolation trauma, as well as the identification of actively transcribed RNAs and not preexisting transcripts. Therefore, in contrast to other cell-specific transcriptional profiling methods based on the purification of tagged ribosomes or nuclei, TU tagging provides a direct examination of transcriptional regulation. We describe how to (i) deliver 4-thiouracil to transgenic mice to thio-label cell lineage-specific transcripts, (ii) purify TU-tagged RNA and prepare libraries for Illumina sequencing and (iii) follow a straightforward bioinformatics workflow to identify cell type-enriched or differentially expressed genes. Tissue containing TU-tagged RNA can be obtained in 1 d, RNA-seq libraries can be generated within 2 d and, after sequencing, an initial bioinformatics analysis can be completed in 1 additional day.

4-thiouridine inhibits rRNA synthesis and causes a nucleolar stress response

High concentrations (> 100 µM) of the ribonucleoside analog 4-thiouridine (4sU) is widely used in methods for RNA analysis like photoactivatable-ribonucleoside-enhanced crosslinking and immunoprecipitation (PAR-CLIP) and nascent messenger (m)RNA labeling (4sU-tagging). Here, we show that 4sU-tagging at low concentrations ≤ 10 µM can be used to measure production and processing of ribosomal (r)RNA. However, elevated concentrations of 4sU (> 50 µM), which are usually used for mRNA labeling experiments, inhibit production and processing of 47S rRNA. The inhibition of rRNA synthesis is accompanied by nucleoplasmic translocation of nucleolar nucleophosmin (NPM1), induction of the tumor suppressor p53, and inhibition of proliferation. We conclude that metabolic labeling of RNA by 4sU triggers a nucleolar stress response, which might influence the interpretation of results. Therefore, functional ribosome biogenesis, nucleolar integrity, and cell cycle should be addressed in 4sU labeling experiments.

Deciphering the modulation of gene expression by type I and II interferons combining 4sU-tagging, translational arrest and in silico promoter analysis

Interferons (IFN) play a pivotal role in innate immunity, orchestrating a cell-intrinsic anti-pathogenic state and stimulating adaptive immune responses. The complex interplay between the primary response to IFNs and its modulation by positive and negative feedback loops is incompletely understood. Here, we implement the combination of high-resolution gene-expression profiling of nascent RNA with translational inhibition of secondary feedback by cycloheximide. Unexpectedly, this approach revealed a prominent role of negative feedback mechanisms during the immediate (≤60 min) IFNα response. In contrast, a more complex picture involving both negative and positive feedback loops was observed on IFNγ treatment. IFNγ-induced repression of genes associated with regulation of gene expression, cellular development, apoptosis and cell growth resulted from cycloheximide-resistant primary IFNγ signalling. In silico promoter analysis revealed significant overrepresentation of SP1/SP3-binding sites and/or GC-rich stretches. Although signal transducer and activator of transcription 1 (STAT1)-binding sites were not overrepresented, repression was lost in absence of STAT1. Interestingly, basal expression of the majority of these IFNγ-repressed genes was dependent on STAT1 in IFN-naïve fibroblasts. Finally, IFNγ-mediated repression was also found to be evident in primary murine macrophages. IFN-repressed genes include negative regulators of innate and stress response, and their decrease may thus aid the establishment of a signalling perceptive milieu.

Genome-wide survey of interindividual differences of RNA stability in human lymphoblastoid cell lines

The extent to which RNA stability differs between individuals and its contribution to the interindividual expression variation remain unknown. We conducted a genome-wide analysis of RNA stability in seven human HapMap lymphoblastoid cell lines (LCLs) and analyzed the effect of DNA sequence variation on RNA half-life differences. Twenty-six percent of the expressed genes exhibited RNA half-life differences between LCLs at a false discovery rate (FDR) < 0.05, which accounted for ~ 37% of the gene expression differences between individuals. Nonsense polymorphisms were associated with reduced RNA half-lives. In genes presenting interindividual RNA half-life differences, higher coding GC3 contents (G and C percentages at the third-codon positions) were correlated with increased RNA half-life. Consistently, G and C alleles of single nucleotide polymorphisms (SNPs) in protein coding sequences were associated with enhanced RNA stability. These results suggest widespread interindividual differences in RNA stability related to DNA sequence and composition variation.

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

The development of whole-transcriptome microarrays and next-generation sequencing has revolutionized our understanding of the complexity of cellular gene expression. Along with a better understanding of the involved molecular mechanisms, precise measurements of the underlying kinetics have become increasingly important. Here, these powerful methodologies face major limitations due to intrinsic properties of the template samples they study, i.e. total cellular RNA. In many cases changes in total cellular RNA occur either too slowly or too quickly to represent the underlying molecular events and their kinetics with sufficient resolution. In addition, the contribution of alterations in RNA synthesis, processing, and decay are not readily differentiated. We recently developed high-resolution gene expression profiling to overcome these limitations. Our approach is based on metabolic labeling of newly transcribed RNA with 4-thiouridine (thus also referred to as 4sU-tagging) followed by rigorous purification of newly transcribed RNA using thiol-specific biotinylation and streptavidin-coated magnetic beads. It is applicable to a broad range of organisms including vertebrates, Drosophila, and yeast. We successfully applied 4sU-tagging to study real-time kinetics of transcription factor activities, provide precise measurements of RNA half-lives, and obtain novel insights into the kinetics of RNA processing. Finally, computational modeling can be employed to generate an integrated, comprehensive analysis of the underlying molecular mechanisms.

Global analysis of mRNA half-lives and de novo transcription in a dinoflagellate, Karenia brevis

Dinoflagellates possess many physiological processes that appear to be under post-transcriptional control. However, the extent to which their genes are regulated post-transcriptionally remains unresolved. To gain insight into the roles of differential mRNA stability and de novo transcription in dinoflagellates, we biosynthetically labeled RNA with 4-thiouracil to isolate newly transcribed and pre-existing RNA pools in Karenia brevis. These isolated fractions were then used for analysis of global mRNA stability and de novo transcription by hybridization to a K. brevis microarray. Global K. brevis mRNA half-lives were calculated from the ratio of newly transcribed to pre-existing RNA for 7086 array features using the online software HALO (Half-life Organizer). Overall, mRNA half-lives were substantially longer than reported in other organisms studied at the global level, ranging from 42 minutes to greater than 144 h, with a median of 33 hours. Consistent with well-documented trends observed in other organisms, housekeeping processes, including energy metabolism and transport, were significantly enriched in the most highly stable messages. Shorter-lived transcripts included a higher proportion of transcriptional regulation, stress response, and other response/regulatory processes. One such family of proteins involved in post-transcriptional regulation in chloroplasts and mitochondria, the pentatricopeptide repeat (PPR) proteins, had dramatically shorter half-lives when compared to the arrayed transcriptome. As transcript abundances for PPR proteins were previously observed to rapidly increase in response to nutrient addition, we queried the newly synthesized RNA pools at 1 and 4 h following nitrate addition to N-depleted cultures. Transcriptome-wide there was little evidence of increases in the rate of de novo transcription during the first 4 h, relative to that in N-depleted cells, and no evidence for increased PPR protein transcription. These results lend support to the growing consensus of post-transcriptional control of gene expression in dinoflagellates.

Eukaryotic mRNA decay: methodologies, pathways, and links to other stages of gene expression

mRNA concentration depends on the balance between transcription and degradation rates. On both sides of the equilibrium, synthesis and degradation show, however, interesting differences that have conditioned the evolution of gene regulatory mechanisms. Here, we discuss recent genome-wide methods for determining mRNA half-lives in eukaryotes. We also review pre- and posttranscriptional regulons that coordinate the fate of functionally related mRNAs by using protein- or RNA-based trans factors. Some of these factors can regulate both transcription and decay rates, thereby maintaining proper mRNA homeostasis during eukaryotic cell life.

High resolution gene expression profiling of RNA synthesis, processing, and decay by metabolic labeling of newly transcribed RNA using 4-thiouridine

Gene expression profiling of changes in total RNA levels has provided invaluable knowledge on the regulation of gene expression. Studies on the kinetics of this regulation, however, have been limited by the fact that total cellular RNA is a poor template for revealing short-term changes in gene expression, alterations in RNA decay rates and the kinetics of RNA processing as well as the differentiation thereof. Here, we describe the metabolic labeling and purification of newly transcribed RNA with 4-thiouridine to study the molecular mechanisms governing RNA synthesis, processing, and decay.

Genome-wide technology for determining RNA stability in mammalian cells: historical perspective and recent advantages based on modified nucleotide labeling

Changing the abundance of transcripts by regulated RNA degradation is a critical step in the control of various biological pathways. Recently, genome-wide inhibitor-free technologies for determining RNA stabilities in mammalian cells have been developed. In these methods, endogenous RNAs are pulse labeled by uridine analogs [e.g., 4-thiouridine (4sU), 5-etyniluridine (EU) and 5'-bromo-uridine (BrU)], followed by purification of labeled de novo RNAs. These technologies have revealed that the specific half-life of each mRNA is closely related to its physiological function. Genes with short-lived mRNAs are significantly enriched among regulatory genes, while genes with long-lived mRNAs are enriched among housekeeping genes. This review describes the recent progress of experimental procedures for measuring RNA stability.

Real-time transcriptional profiling of cellular and viral gene expression during lytic cytomegalovirus infection

During viral infections cellular gene expression is subject to rapid alterations induced by both viral and antiviral mechanisms. In this study, we applied metabolic labeling of newly transcribed RNA with 4-thiouridine (4sU-tagging) to dissect the real-time kinetics of cellular and viral transcriptional activity during lytic murine cytomegalovirus (MCMV) infection. Microarray profiling on newly transcribed RNA obtained at different times during the first six hours of MCMV infection revealed discrete functional clusters of cellular genes regulated with distinct kinetics at surprising temporal resolution. Immediately upon virus entry, a cluster of NF-κB- and interferon-regulated genes was induced. Rapid viral counter-regulation of this coincided with a very transient DNA-damage response, followed by a delayed ER-stress response. Rapid counter-regulation of all three clusters indicated the involvement of novel viral regulators targeting these pathways. In addition, down-regulation of two clusters involved in cell-differentiation (rapid repression) and cell-cycle (delayed repression) was observed. Promoter analysis revealed all five clusters to be associated with distinct transcription factors, of which NF-κB and c-Myc were validated to precisely match the respective transcriptional changes observed in newly transcribed RNA. 4sU-tagging also allowed us to study the real-time kinetics of viral gene expression in the absence of any interfering virion-associated-RNA. Both qRT-PCR and next-generation sequencing demonstrated a sharp peak of viral gene expression during the first two hours of infection including transcription of immediate-early, early and even well characterized late genes. Interestingly, this was subject to rapid gene silencing by 5–6 hours post infection. Despite the rapid increase in viral DNA load during viral DNA replication, transcriptional activity of some viral genes remained remarkably constant until late-stage infection, or was subject to further continuous decline. In summary, this study pioneers real-time transcriptional analysis during a lytic herpesvirus infection and highlights numerous novel regulatory aspects of virus-host-cell interaction.

Cytomegaloviruses are large DNA viruses, which establish life-long latent infections, leaving the infected individual at risk of reactivation and disease. Here, we applied 4-thiouridine-(4sU)-tagging of newly transcribed RNA to monitor the real-time kinetics of transcriptional activity of both cellular and viral genes during lytic murine CMV (MCMV) infection. We observed a cascade of MCMV-induced signaling events including a rapid inflammatory/interferon-response, a transient DNA-damage-response and a delayed ER-stress-response. All of these were heavily counter-regulated by viral gene expression. Besides dramatically increasing temporal resolution, our approach provides the unique opportunity to study viral transcriptional activity in absence of any interfering virion-associated-RNA. Virion-associated-RNA consists of transcripts that are unspecifically incorporated into the virus particles thereby resembling the cellular RNA profile of late stage infection. A clear picture of which viral genes are expressed, particularly at very early times of infection, could thus not be obtained. By overcoming this problem, we provide intriguing insights into the regulation of viral gene expression, namely 1) a peak of viral gene expression during the first two hours of infection including the expression of well-characterized late genes and 2) remarkably constant or even continuously declining expression of some viral genes despite the onset of rapid viral DNA replication.

Detection and correction of probe-level artefacts on microarrays

Background

A recent large-scale analysis of Gene Expression Omnibus (GEO) data found frequent evidence for spatial defects in a substantial fraction of Affymetrix microarrays in the GEO. Nevertheless, in contrast to quality assessment, artefact detection is not widely used in standard gene expression analysis pipelines. Furthermore, although approaches have been proposed to detect diverse types of spatial noise on arrays, the correction of these artefacts is mostly left to either summarization methods or the corresponding arrays are completely discarded.

Results

We show that state-of-the-art robust summarization procedures are vulnerable to artefacts on arrays and cannot appropriately correct for these. To address this problem, we present a simple approach to detect artefacts with high recall and precision, which we further improve by taking into account the spatial layout of arrays. Finally, we propose two correction methods for these artefacts that either substitute values of defective probes using probeset information or filter corrupted probes. We show that our approach can identify and correct defective probe measurements appropriately and outperforms existing tools.

Conclusions

While summarization is insufficient to correct for defective probes, this problem can be addressed in a straightforward way by the methods we present for identification and correction of defective probes. As these methods output CEL files with corrected probe values that serve as input to standard normalization and summarization procedures, they can be easily integrated into existing microarray analysis pipelines as an additional pre-processing step. An R package is freely available from http://www.bio.ifi.lmu.de/artefact-correction.

Ultrashort and progressive 4sU-tagging reveals key characteristics of RNA processing at nucleotide resolution

RNA synthesis and decay rates determine the steady-state levels of cellular RNAs. Metabolic tagging of newly transcribed RNA by 4-thiouridine (4sU) can reveal the relative contributions of RNA synthesis and decay rates. The kinetics of RNA processing, however, had so far remained unresolved. Here, we show that ultrashort 4sU-tagging not only provides snapshot pictures of eukaryotic gene expression but, when combined with progressive 4sU-tagging and RNA-seq, reveals global RNA processing kinetics at nucleotide resolution. Using this method, we identified classes of rapidly and slowly spliced/degraded introns. Interestingly, each class of splicing kinetics was characterized by a distinct association with intron length, gene length, and splice site strength. For a large group of introns, we also observed long lasting retention in the primary transcript, but efficient secondary splicing or degradation at later time points. Finally, we show that processing of most, but not all small nucleolar (sno)RNA-containing introns is remarkably inefficient with the majority of introns being spliced and degraded rather than processed into mature snoRNAs. In summary, our study yields unparalleled insights into the kinetics of RNA processing and provides the tools to study molecular mechanisms of RNA processing and their contribution to the regulation of gene expression.

Genome-wide studies of mRNA synthesis and degradation in eukaryotes

In recent years, the use of genome-wide technologies has revolutionized the study of eukaryotic transcription producing results for thousands of genes at every step of mRNA life. The statistical analyses of the results for a single condition, different conditions, different transcription stages, or even between different techniques, is outlining a totally new landscape of the eukaryotic transcription process. Although most studies have been conducted in the yeast Saccharomyces cerevisiae as a model cell, others have also focused on higher eukaryotes, which can also be comparatively analyzed. The picture which emerges is that transcription is a more variable process than initially suspected, with large differences between genes at each stage of the process, from initiation to mRNA degradation, but with striking similarities for functionally related genes, indicating that all steps are coordinately regulated. This article is part of a Special Issue entitled: Nuclear Transport and RNA Processing.

Genomic insights into the different layers of gene regulation in yeast

The model organism Saccharomyces cerevisiae has allowed the development of new functional genomics techniques devoted to the study of transcription in all its stages. With these techniques, it has been possible to find interesting new mechanisms to control gene expression that act at different levels and for different gene sets apart from the known cis-trans regulation in the transcription initiation step. Here we discuss a method developed in our laboratory, Genomic Run-On, and other new methods that have recently appeared with distinct technical features. A comparison between the datasets generated by them provides interesting genomic insights into the different layers of gene regulation in eukaryotes.

Regulatory mechanisms and networks couple the different phases of gene expression

Gene expression comprises multiple stages, from transcription to protein degradation. Although much is known about the regulation of each stage separately, an understanding of the regulatory coupling between the different stages is only beginning to emerge. For example, there is a clear crosstalk between translation and transcription, and the localization and stability of an mRNA in the cytoplasm could already be determined during transcription in the nucleus. We review a diversity of mechanisms discovered in recent years that couple the different stages of gene expression. We then speculate on the functional and evolutionary significance of this coupling and suggest certain systems-level functionalities that might be optimized via the various coupling modes. In particular, we hypothesize that coupling is often an economic strategy that allows biological systems to respond robustly and precisely to genetic and environmental perturbations.

Metabolic labeling of RNA uncovers principles of RNA production and degradation dynamics in mammalian cells

Cellular RNA levels are determined by the interplay of RNA production, processing and degradation. However, because most studies of RNA regulation do not distinguish the separate contributions of these processes, little is known about how they are temporally integrated. Here we combine metabolic labeling of RNA at high temporal resolution with advanced RNA quantification and computational modeling to estimate RNA transcription and degradation rates during the response of mouse dendritic cells to lipopolysaccharide. We find that changes in transcription rates determine the majority of temporal changes in RNA levels, but that changes in degradation rates are important for shaping sharp 'peaked' responses. We used sequencing of the newly transcribed RNA population to estimate temporally constant RNA processing and degradation rates genome wide. Degradation rates vary significantly between genes and contribute to the observed differences in the dynamic response. Certain transcripts, including those encoding cytokines and transcription factors, mature faster. Our study provides a quantitative approach to study the integrative process of RNA regulation.

Synthesis and decay of varicella zoster virus transcripts

Varicella zoster virus (VZV) is highly cell-associated. At least 68 VZV open reading frames (ORFs) are transcribed in varying amounts that increase as infection progresses. Using reverse transcriptase PCR, quantification of total and newly synthesized mRNA showed that ongoing VZV DNA replication is required for continued accumulation of VZV ORF 63, 9, and 40 transcripts. Analysis of stability of 4-thiouridine-labeled transcripts of nine VZV ORFs revealed a similar half-life for all VZV ORFs tested. Thus, difference in mRNA synthesis, and not mRNA decay, is the major factor contributing to the difference in the relative abundance of VZV transcripts in infected cells.

Dynamic transcriptome analysis measures rates of mRNA synthesis and decay in yeast

To obtain rates of mRNA synthesis and decay in yeast, we established dynamic transcriptome analysis (DTA). DTA combines non-perturbing metabolic RNA labeling with dynamic kinetic modeling. DTA reveals that most mRNA synthesis rates are around several transcripts per cell and cell cycle, and most mRNA half-lives range around a median of 11 min. DTA can monitor the cellular response to osmotic stress with higher sensitivity and temporal resolution than standard transcriptomics. In contrast to monotonically increasing total mRNA levels, DTA reveals three phases of the stress response. During the initial shock phase, mRNA synthesis and decay rates decrease globally, resulting in mRNA storage. During the subsequent induction phase, both rates increase for a subset of genes, resulting in production and rapid removal of stress-responsive mRNAs. During the recovery phase, decay rates are largely restored, whereas synthesis rates remain altered, apparently enabling growth at high salt concentration. Stress-induced changes in mRNA synthesis rates are predicted from gene occupancy with RNA polymerase II. DTA-derived mRNA synthesis rates identified 16 stress-specific pairs/triples of cooperative transcription factors, of which seven were known. Thus, DTA realistically monitors the dynamics in mRNA metabolism that underlie gene regulatory systems.

HALO--a Java framework for precise transcript half-life determination

Recent improvements in experimental technologies now allow measurements of de novo transcription and/or RNA decay at whole transcriptome level and determination of precise transcript half-lives. Such transcript half-lives provide important insights into the regulation of biological processes and the relative contributions of RNA decay and de novo transcription to differential gene expression. In this article, we present HALO (Half-life Organizer), the first software for the precise determination of transcript half-lives from measurements of RNA de novo transcription or decay determined with microarrays or RNA-seq. In addition, methods for quality control, filtering and normalization are supplied. HALO provides a graphical user interface, command-line tools and a well-documented Java application programming interface (API). Thus, it can be used both by biologists to determine transcript half-lives fast and reliably with the provided user interfaces as well as software developers integrating transcript half-life analysis into other gene expression profiling pipelines.

AVAILABILITY

Source code, executables and documentation are available at http://www.bio.ifi.lmu.de/software/halo.

Time-dependent c-Myc transactomes mapped by Array-based nuclear run-on reveal transcriptional modules in human B cells

Background

The definition of transcriptional networks through measurements of changes in gene expression profiles and mapping of transcription factor binding sites is limited by the moderate overlap between binding and gene expression changes and the inability to directly measure global nuclear transcription (coined “transactome”).

Methodology/Principal Findings

We developed a method to measure nascent nuclear gene transcription with an Array-based Nuclear Run-On (ANRO) assay using commercial microarray platforms. This strategy provides the missing component, the transactome, to fully map transcriptional networks. ANRO measurements in an inducible c-Myc expressing human P493-6 B cell model reveals time-dependent waves of transcription, with a transactome early after c-Myc induction that does not persist at a late, steady-state phase, when genes that are regulated by c-Myc and E2F predominate. Gene set matrix analysis further uncovers functionally related groups of genes putatively regulated by waves of transcription factor motifs following Myc induction, starting with AP1 and CREB that are followed by EGR1, NFkB and STAT, and ending with E2F, Myc and ARNT/HIF motifs.

Conclusions/Significance

By coupling ANRO with previous global mapping of c-Myc binding sites by chromatin immunoprecipitation (ChIP) in P493-6 cells, we define a set of transcriptionally regulated direct c-Myc target genes and pave the way for the use of ANRO to comprehensively map any transcriptional network.

Metabolic tagging and purification of nascent RNA: implications for transcriptomics

Gene expression profiling to analyze cellular responses against different stimuli or conditions is generally performed at the total cellular RNA level. This results in poor resolution of the temporal kinetics of the cellular response and a bias towards detecting up-regulation of short-lived transcripts. Furthermore, changes in transcription rate and RNA stability cannot be distinguished. These problems can be addressed by analyzing nascent RNA instead of total cellular RNA. Throughout the last few years methods have been developed for metabolic tagging and purification of nascent RNA. In this article, we review these experimental procedures and discuss their implications for large-scale gene expression profiling.

Endo-siRNAs depend on a new isoform of loquacious and target artificially introduced, high-copy sequences

Colonization of genomes by a new selfish genetic element is detrimental to the host species and must lead to an efficient, repressive response. In vertebrates as well as in Drosophila, piRNAs repress transposons in the germ line, whereas endogenous siRNAs take on this role in somatic cells. We show that their biogenesis depends on a new isoform of the Drosophila TRBP homologue loquacious, which arises by alternative polyadenylation and is distinct from the one that functions during the biogenesis of miRNAs. For endo-siRNAs and piRNAs, it is unclear how an efficient response can be initiated de novo. Our experiments establish that the endo-siRNA pathway will target artificially introduced sequences without the need for a pre-existing template in the genome. This response is also triggered in transiently transfected cells, thus genomic integration is not essential. Deep sequencing showed that corresponding endo-siRNAs are generated throughout the sequence, but preferentially from transcribed regions. One strand of the dsRNA precursor can come from spliced mRNA, whereas the opposite strand derives from independent transcripts in antisense orientation.

Conserved principles of mammalian transcriptional regulation revealed by RNA half-life

RNA levels in a cell are regulated by the relative rates of RNA synthesis and decay. We recently developed a new approach for measuring both RNA synthesis and decay in a single experimental setting by biosynthetic labeling of newly transcribed RNA. Here, we show that this provides measurements of RNA half-lives from microarray data with a so far unreached accuracy. Based on such measurements of RNA half-lives for human B-cells and mouse fibroblasts, we identified conserved regulatory principles for a large number of biological processes. We show that different regulatory patterns between functionally similar proteins are characterized by differences in the half-life of the corresponding transcripts and can be identified by measuring RNA half-life. We identify more than 100 protein families which show such differential regulatory patterns in both species. Additionally, we provide strong evidence that the activity of protein complexes consisting of subunits with overall long transcript half-lives can be regulated by transcriptional regulation of individual key subunits with short-lived transcripts. Based on this observation, we predict more than 100 key regulatory subunits for human complexes of which 28% could be confirmed in mice (P < 10⁻⁹). Therefore, this atlas of transcript half-lives provides new fundamental insights into many cellular processes.

An upstream open reading frame regulates translation of GADD34 during cellular stresses that induce eIF2alpha phosphorylation

Cellular stress such as endoplasmic reticulum stress, hypoxia, and viral infection activates an integrated stress response, which includes the phosphorylation of the eukaryotic initiation factor 2alpha (eIF2alpha) to inhibit overall protein synthesis. Paradoxically, this leads to translation of a subset of mRNAs, like transcription factor ATF4, which in turn induces transcription of downstream stress-induced genes such as growth arrest DNA-inducible gene 34 (GADD34). GADD34 interacts with protein phosphatase 1 to dephosphorylate eIF2alpha, resulting in a negative feedback loop to recover protein synthesis and allow translation of stress-induced transcripts. Here, we show that GADD34 is not only transcriptionally induced but also translationally regulated to ensure maximal expression during eIF2alpha phosphorylation. GADD34 mRNAs are preferentially associated with polysomes during eIF2alpha phosphorylation, which is mediated by its 5'-untranslated region (5'UTR). The human GADD34 5'UTR contains two non-overlapping upstream open reading frames (uORFs), whereas the mouse version contains two overlapping and out of frame uORFs. Using 5'UTR GADD34 reporter constructs, we show that the downstream uORF mediates repression of basal translation and directs translation during eIF2alpha phosphorylation. Furthermore, we show that the upstream uORF is poorly translated and that a proportion of scanning ribosomes bypasses the upstream uORF to recognize the downstream uORF. These findings suggest that GADD34 translation is regulated by a unique 5'UTR uORF mechanism to ensure proper GADD34 expression during eIF2alpha phosphorylation. This mechanism may serve as a model for understanding how other 5'UTR uORF-containing mRNAs are regulated during cellular stress.

High-resolution gene expression profiling for simultaneous kinetic parameter analysis of RNA synthesis and decay

RNA levels in a cell are determined by the relative rates of RNA synthesis and decay. State-of-the-art transcriptional analyses only employ total cellular RNA. Therefore, changes in RNA levels cannot be attributed to RNA synthesis or decay, and temporal resolution is poor. Recently, it was reported that newly transcribed RNA can be biosynthetically labeled for 1-2 h using thiolated nucleosides, purified from total cellular RNA and subjected to microarray analysis. However, in order to study signaling events at molecular level, analysis of changes occurring within minutes is required. We developed an improved approach to separate total cellular RNA into newly transcribed and preexisting RNA following 10-15 min of metabolic labeling. Employing new computational tools for array normalization and half-life determination we simultaneously study short-term RNA synthesis and decay as well as their impact on cellular transcript levels. As an example we studied the response of fibroblasts to type I and II interferons (IFN). Analysis of RNA transcribed within 15-30 min at different times during the first three hours of interferon-receptor activation resulted in a >10-fold increase in microarray sensitivity and provided a comprehensive profile of the kinetics of IFN-mediated changes in gene expression. We identify a previously undisclosed highly connected network of short-lived transcripts selectively down-regulated by IFNgamma in between 30 and 60 min after IFN treatment showing strong associations with cell cycle and apoptosis, indicating novel mechanisms by which IFNgamma affects these pathways.

Novel DNA microarray system for analysis of nascent mRNAs

Transcriptional activation and repression are a key step in the regulation of all cellular activities. The development of comprehensive analysis methods such as DNA microarray has advanced our understanding of the correlation between the regulation of transcription and that of cellular mechanisms. However, DNA microarray analysis based on steady-state mRNA (total mRNA) does not always correspond to transcriptional activation or repression. To comprehend these transcriptional regulations, the detection of nascent RNAs is more informative. Although the nuclear run-on assay can detect nascent RNAs, it has not been fully applied to DNA microarray analysis. In this study, we have developed a highly efficient method for isolating bromouridine-labeled nascent RNAs that can be successfully applied to DNA microarray analysis. This method can linearly amplify small amounts of mRNAs with little bias. Furthermore, we have applied this method to DNA microarray analysis from mouse G₂-arrested cells and have identified several genes that exhibit novel expression profiles. This method will provide important information in the field of transcriptome analysis of various cellular processes.