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

3' UTR-isoform choice has limited influence on the stability and translational efficiency of most mRNAs in mouse fibroblasts

Variation in protein output across the genome is controlled at several levels, but the relative contributions of different regulatory mechanisms remain poorly understood. Here, we obtained global measurements of decay and translation rates for mRNAs with alternative 3′ untranslated regions (3′ UTRs) in murine 3T3 cells. Distal tandem isoforms had slightly but significantly lower mRNA stability and greater translational efficiency than proximal isoforms on average. The diversity of alternative 3′ UTRs also enabled inference and evaluation of both positively and negatively acting cis-regulatory elements. The 3′ UTR elements with the greatest implied influence were microRNA complementary sites, which were associated with repression of 32% and 4% at the stability and translational levels, respectively. Nonetheless, both the decay and translation rates were highly correlated for proximal and distal 3′ UTR isoforms from the same genes, implying that in 3T3 cells, alternative 3′ UTR sequences play a surprisingly small regulatory role compared to other mRNA regions.

Inactivation of intergenic enhancers by EBNA3A initiates and maintains polycomb signatures across a chromatin domain encoding CXCL10 and CXCL9

Epstein-Barr virus (EBV) causes a persistent infection in human B cells by establishing specific transcription programs to control B cell activation and differentiation. Transcriptional reprogramming of EBV infected B cells is predominantly driven by the action of EBV nuclear antigens, among them the transcriptional repressor EBNA3A. By comparing gene expression profiles of wt and EBNA3A negative EBV infected B cells, we have previously identified a broad array of cellular genes controlled by EBNA3A. We now find that genes repressed by EBNA3A in these cells are significantly enriched for the repressive histone mark H3K27me3, which is installed by Polycomb group (PcG) proteins. This PcG-controlled subset of genes also carries H3K27me3 marks in a variety of other tissues, suggesting that the commitment to PcG silencing is an intrinsic feature of these gene loci that can be used by EBNA3A. In addition, EBNA3A targets frequently reside in co-regulated gene clusters. To study the mechanism of gene repression by EBNA3A and to evaluate the relative contribution of PcG proteins during this process, we have selected the genomic neighbors CXCL10 and CXCL9 as a model for co-repressed and PcG-controlled genes. We show that EBNA3A binds to CBF1 occupied intergenic enhancers located between CXCL10 and CXCL9 and displaces the transactivator EBNA2. This impairs enhancer activity, resulting in a rapid transcriptional shut-down of both genes in a CBF1-dependent manner and initiation of a delayed gain of H3K27me3 marks covering an extended chromatin domain. H3K27me3 marks increase gradually and are maintained by EBNA3A. Our study provides direct evidence that repression by EBNA3A requires CBF1 and that EBNA3A and EBNA2 compete for access to CBF1 at identical genomic sites. Most importantly, our results demonstrate that transcriptional silencing by EBNA3A precedes the appearance of repressive PcG marks and indicate that both events are triggered by loss of enhancer activity.

Epstein-Barr virus (EBV) is a γ-herpesvirus which establishes a latent infection in human B cells and is associated with the pathogenesis of several types of cancer. Here, we report that cellular genes repressed by the EBV nuclear antigen 3A (EBNA3A) in EBV infected B cells frequently form contiguous clusters in the human genome and are committed to epigenetic silencing by Polycomb group (PcG) proteins. The chemokine genes CXCL10 and CXCL9 and their receptors on NK and T cells are critical weapons of the infected host to control herpesvirus infections. CXCL10 and CXCL9 are close neighbors within an extended PcG-controlled domain. We show that EBNA3A binds to intergenic enhancers located between CXCL10 and CXCL9 and displaces the transactivator EBNA2. This process impairs enhancer activity, resulting in a rapid transcriptional shut-down of both genes followed by a delayed gain of PcG histone marks. These PcG marks increase within the following weeks and are maintained by EBNA3A. Our results show that rapid transcriptional shut-down of distal genes and domain-wide PcG silencing is triggered by loss of enhancer activity and suggest that EBNA3A can reprogram the cellular genome in order to escape the immune surveillance of the host.

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.

Minor introns are embedded molecular switches regulated by highly unstable U6atac snRNA

Eukaryotes have two types of spliceosomes, comprised of either major (U1, U2, U4, U5, U6) or minor (U11, U12, U4atac, U6atac; <1%) snRNPs. The high conservation of minor introns, typically one amidst many major introns in several hundred genes, despite their poor splicing, has been a long-standing enigma. Here, we discovered that the low abundance minor spliceosome's catalytic snRNP, U6atac, is strikingly unstable (t½<2 hr). We show that U6atac level depends on both RNA polymerases II and III and can be rapidly increased by cell stress-activated kinase p38MAPK, which stabilizes it, enhancing mRNA expression of hundreds of minor intron-containing genes that are otherwise suppressed by limiting U6atac. Furthermore, p38MAPK-dependent U6atac modulation can control minor intron-containing tumor suppressor PTEN expression and cytokine production. We propose that minor introns are embedded molecular switches regulated by U6atac abundance, providing a novel post-transcriptional gene expression mechanism and a rationale for the minor spliceosome's evolutionary conservation. DOI:http://dx.doi.org/10.7554/eLife.00780.001.

The viral chemokine MCK-2 of murine cytomegalovirus promotes infection as part of a gH/gL/MCK-2 complex

Human cytomegalovirus (HCMV) forms two gH/gL glycoprotein complexes, gH/gL/gO and gH/gL/pUL(128,130,131A), which determine the tropism, the entry pathways and the mode of spread of the virus. For murine cytomegalovirus (MCMV), which serves as a model for HCMV, a gH/gL/gO complex functionally homologous to the HCMV gH/gL/gO complex has been described. Knock-out of MCMV gO does impair, but not abolish, virus spread indicating that also MCMV might form an alternative gH/gL complex. Here, we show that the MCMV CC chemokine MCK-2 forms a complex with the glycoprotein gH, a complex which is incorporated into the virion. We could additionally show that mutants lacking both, gO and MCK-2 are not able to produce infectious virus. Trans-complementation of these double mutants with either gO or MCK-2 showed that both proteins can promote infection of host cells, although through different entry pathways. MCK-2 has been extensively studied in vivo by others. It has been shown to be involved in attracting cells for virus dissemination and in regulating antiviral host responses. We now show that MCK-2, by forming a complex with gH, strongly promotes infection of macrophages in vitro and in vivo. Thus, MCK-2 may play a dual role in MCMV infection, as a chemokine regulating the host response and attracting specific target cells and as part of a glycoprotein complex promoting entry into cells crucial for virus dissemination.

Several human herpesviruses form alternative gH/gL complexes which determine the tropism for different cell types. For murine cytomegalovirus (MCMV), a gH/gL/gO complex has recently been characterized. Here, we present the identification and characterization of an alternative gH/gL/MCK-2 complex which promotes MCMV spread and is important for efficient infection of macrophages in vitro and in vivo. Association of the MCMV CC chemokine MCK-2 with a glycoprotein complex promoting virus entry is a novel function for the well-characterized MCK-2. Virus mutants lacking MCK-2 have been shown to exhibit a reduced capacity to attract leukocytes and a disregulated T cell control of the MCMV infection in vivo. These defects can be attributed to the chemokine function of MCK-2. Yet, the observation that MCK-2 knock-out mutants additionally are impaired in infecting leukocytes in vivo is consistent with our new finding that MCK-2 forms a glycoprotein complex promoting entry into monocytic cells. gH/gL complexes associating with multifunctional proteins add a new level of complexity to the interpretation of infection phenotypes of the respective knock-out herpesviruses.

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.

Exonuclease hDIS3L2 specifies an exosome-independent 3'-5' degradation pathway of human cytoplasmic mRNA

Turnover of mRNA in the cytoplasm of human cells is thought to be redundantly conducted by the monomeric 5'-3' exoribonuclease hXRN1 and the 3'-5' exoribonucleolytic RNA exosome complex. However, in addition to the exosome-associated 3'-5' exonucleases hDIS3 and hDIS3L, the human genome encodes another RNase II/R domain protein-hDIS3L2. Here, we show that hDIS3L2 is an exosome-independent cytoplasmic mRNA 3'-5' exonuclease, which exhibits processive activity on structured RNA substrates in vitro. hDIS3L2 associates with hXRN1 in an RNA-dependent manner and can, like hXRN1, be found on polysomes. The impact of hDIS3L2 on cytoplasmic RNA metabolism is revealed by an increase in levels of cytoplasmic RNA processing bodies (P-bodies) upon hDIS3L2 depletion, which also increases half-lives of investigated mRNAs. Consistently, RNA sequencing (RNA-seq) analyses demonstrate that depletion of hDIS3L2, like downregulation of hXRN1 and hDIS3L, causes changed levels of multiple mRNAs. We suggest that hDIS3L2 is a key exosome-independent effector of cytoplasmic mRNA metabolism.

Nontransformed, GM-CSF-dependent macrophage lines are a unique model to study tissue macrophage functions

Macrophages are diverse cell types in the first line of antimicrobial defense. Only a limited number of primary mouse models exist to study their function. Bone marrow-derived, macrophage-CSF-induced cells with a limited life span are the most common source. We report here a simple method yielding self-renewing, nontransformed, GM-CSF/signal transducer and activator of transcription 5-dependent macrophages (Max Planck Institute cells) from mouse fetal liver, which reflect the innate immune characteristics of alveolar macrophages. Max Planck Institute cells are exquisitely sensitive to selected microbial agents, including bacterial LPS, lipopeptide, Mycobacterium tuberculosis, cord factor, and adenovirus and mount highly proinflammatory but no anti-inflammatory IL-10 responses. They show a unique pattern of innate responses not yet observed in other mononuclear phagocytes. This includes differential LPS sensing and an unprecedented regulation of IL-1α production upon LPS exposure, which likely plays a key role in lung inflammation in vivo. In conclusion, Max Planck Institute cells offer an useful tool to study macrophage biology and for biomedical science.

Drosophila miR-277 controls branched-chain amino acid catabolism and affects lifespan

Development, growth and adult survival are coordinated with available metabolic resources, ascertaining that the organism responds appropriately to environmental conditions. MicroRNAs are short (21-23 nt) regulatory RNAs that confer specificity on the RNA-induced silencing complex (RISC) to inhibit a given set of mRNA targets. We profiled changes in miRNA expression during adult life in Drosophila melanogaster and determined that miR-277 is downregulated during adult life. Molecular analysis revealed that this miRNA controls branched-chain amino acid (BCAA) catabolism and as a result it can modulate the activity of the TOR kinase, a central growth regulator, in cultured cells. Metabolite analysis in cultured cells as well as flies suggests that the mechanistic basis may be an accumulation of branched-chain α-keto-acids (BCKA), rather than BCAAs, thus avoiding potentially detrimental consequences of increased branched chain amino acid levels on e.g., translational fidelity. Constitutive miR-277 expression shortens lifespan and is synthetically lethal with reduced insulin signaling, indicating that metabolic control underlies this phenotype. Transgenic inhibition with a miRNA sponge construct also shortens lifespan, in particular on protein-rich food. Thus, optimal metabolic adaptation appears to require tuning of cellular BCAA catabolism by miR-277.

Complex degradation processes lead to non-exponential decay patterns and age-dependent decay rates of messenger RNA

Experimental studies on mRNA stability have established several, qualitatively distinct decay patterns for the amount of mRNA within the living cell. Furthermore, a variety of different and complex biochemical pathways for mRNA degradation have been identified. The central aim of this paper is to bring together both the experimental evidence about the decay patterns and the biochemical knowledge about the multi-step nature of mRNA degradation in a coherent mathematical theory. We first introduce a mathematical relationship between the mRNA decay pattern and the lifetime distribution of individual mRNA molecules. This relationship reveals that the mRNA decay patterns at steady state expression level must obey a general convexity condition, which applies to any degradation mechanism. Next, we develop a theory, formulated as a Markov chain model, that recapitulates some aspects of the multi-step nature of mRNA degradation. We apply our theory to experimental data for yeast and explicitly derive the lifetime distribution of the corresponding mRNAs. Thereby, we show how to extract single-molecule properties of an mRNA, such as the age-dependent decay rate and the residual lifetime. Finally, we analyze the decay patterns of the whole translatome of yeast cells and show that yeast mRNAs can be grouped into three broad classes that exhibit three distinct decay patterns. This paper provides both a method to accurately analyze non-exponential mRNA decay patterns and a tool to validate different models of degradation using decay data.

Mechanisms coordinating ELAV/Hu mRNA regulons

The 5' and 3' untranslated regions (UTRs) of messenger RNAs (mRNAs) function as platforms that can determine the fate of each mRNA individually and in aggregate. Multiple mRNAs that encode proteins that are functionally related often interact with RNA-binding proteins (RBPs) and noncoding RNAs (ncRNAs) that coordinate their expression in time and space as RNA regulons within the ribonucleoprotein (RNP) infrastructure we term the ribonome. Recent ribonomic methods have emerged that can determine which mRNAs are bound and regulated by RBPs and ncRNAs, some of which act in combination to determine global outcomes. ELAV/Hu proteins bind to AU-rich elements (ARE) in mRNAs and regulate their stability from splicing to translation, and the ubiquitous HuR protein has been implicated in cancerous cell growth. Recent work is focused on mechanistic models of how ELAV/Hu proteins increase mRNA stability and translation by repressing microRNAs (miRs) and the RNA induced silencing complex (RISC) via ARE-based ribonucleosomes that may affect global functions of mRNA regulons.

CDK8 kinase phosphorylates transcription factor STAT1 to selectively regulate the interferon response

Gene regulation by cytokine-activated transcription factors of the signal transducer and activator of transcription (STAT) family requires serine phosphorylation within the transactivation domain (TAD). STAT1 and STAT3 TAD phosphorylation occurs upon promoter binding by an unknown kinase. Here, we show that the cyclin-dependent kinase 8 (CDK8) module of the Mediator complex phosphorylated regulatory sites within the TADs of STAT1, STAT3, and STAT5, including S727 within the STAT1 TAD in the interferon (IFN) signaling pathway. We also observed a CDK8 requirement for IFN-γ-inducible antiviral responses. Microarray analyses revealed that CDK8-mediated STAT1 phosphorylation positively or negatively regulated over 40% of IFN-γ-responsive genes, and RNA polymerase II occupancy correlated with gene expression changes. This divergent regulation occurred despite similar CDK8 occupancy at both S727 phosphorylation-dependent and -independent genes. These data identify CDK8 as a key regulator of STAT1 and antiviral responses and suggest a general role for CDK8 in STAT-mediated transcription. As such, CDK8 represents a promising target for therapeutic manipulation of cytokine responses.

Coordinated regulation of synthesis and stability of RNA during the acute TNF-induced proinflammatory response

Steady-state gene expression is a coordination of synthesis and decay of RNA through epigenetic regulation, transcription factors, micro RNAs (miRNAs), and RNA-binding proteins. Here, we present bromouride labeling and sequencing (Bru-Seq) and bromouridine pulse-chase and sequencing (BruChase-Seq) to assess genome-wide changes to RNA synthesis and stability in human fibroblasts at homeostasis and after exposure to the proinflammatory tumor necrosis factor (TNF). The inflammatory response in human cells involves rapid and dramatic changes in gene expression, and the Bru-Seq and BruChase-Seq techniques revealed a coordinated and complex regulation of gene expression both at the transcriptional and posttranscriptional levels. The combinatory analysis of both RNA synthesis and stability using Bru-Seq and BruChase-Seq allows for a much deeper understanding of mechanisms of gene regulation than afforded by the analysis of steady-state total RNA and should be useful in many biological settings.

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.

Human mitochondrial RNA decay mediated by PNPase-hSuv3 complex takes place in distinct foci

RNA decay is usually mediated by protein complexes and can occur in specific foci such as P-bodies in the cytoplasm of eukaryotes. In human mitochondria nothing is known about the spatial organization of the RNA decay machinery, and the ribonuclease responsible for RNA degradation has not been identified. We demonstrate that silencing of human polynucleotide phosphorylase (PNPase) causes accumulation of RNA decay intermediates and increases the half-life of mitochondrial transcripts. A combination of fluorescence lifetime imaging microscopy with Förster resonance energy transfer and bimolecular fluorescence complementation (BiFC) experiments prove that PNPase and hSuv3 helicase (Suv3, hSuv3p and SUPV3L1) form the RNA-degrading complex in vivo in human mitochondria. This complex, referred to as the degradosome, is formed only in specific foci (named D-foci), which co-localize with mitochondrial RNA and nucleoids. Notably, interaction between PNPase and hSuv3 is essential for efficient mitochondrial RNA degradation. This provides indirect evidence that degradosome-dependent mitochondrial RNA decay takes place in foci.

Inhibition of polyadenylation reduces inflammatory gene induction

Cordycepin (3' deoxyadenosine) has long been used in the study of in vitro assembled polyadenylation complexes, because it terminates the poly(A) tail and arrests the cleavage complex. It is derived from caterpillar fungi, which are highly prized in Chinese traditional medicine. Here we show that cordycepin specifically inhibits the induction of inflammatory mRNAs by cytokines in human airway smooth muscle cells without affecting the expression of control mRNAs. Cordycepin treatment results in shorter poly(A) tails, and a reduction in the efficiency of mRNA cleavage and transcription termination is observed, indicating that the effects of cordycepin on 3' processing in cells are similar to those described in in vitro reactions. For the CCL2 and CXCL1 mRNAs, the effects of cordycepin are post-transcriptional, with the mRNA disappearing during or immediately after nuclear export. In contrast, although the recruitment of RNA polymerase II to the IL8 promoter is also unaffected, the levels of nascent transcript are reduced, indicating a defect in transcription elongation. We show that a reporter construct with 3' sequences from a histone gene is unaffected by cordycepin, while CXCL1 sequences confer cordycepin sensitivity to the reporter, demonstrating that polyadenylation is indeed required for the effect of cordycepin on gene expression. In addition, treatment with another polyadenyation inhibitor and knockdown of poly(A) polymerase α also specifically reduced the induction of inflammatory mRNAs. These data demonstrate that there are differences in the 3' processing of inflammatory and housekeeping genes and identify polyadenylation as a novel target for anti-inflammatory drugs.

Genome-wide RNA polymerase II profiles and RNA accumulation reveal kinetics of transcription and associated epigenetic changes during diurnal cycles

Genome-wide rhythms in RNA polymerase II loading and dynamic chromatin remodeling underlie periodic gene expression during diurnal cycles in the mouse liver.

Interactions of cell-autonomous circadian oscillators with diurnal cycles govern the temporal compartmentalization of cell physiology in mammals. To understand the transcriptional and epigenetic basis of diurnal rhythms in mouse liver genome-wide, we generated temporal DNA occupancy profiles by RNA polymerase II (Pol II) as well as profiles of the histone modifications H3K4me3 and H3K36me3. We used these data to quantify the relationships of phases and amplitudes between different marks. We found that rhythmic Pol II recruitment at promoters rather than rhythmic transition from paused to productive elongation underlies diurnal gene transcription, a conclusion further supported by modeling. Moreover, Pol II occupancy preceded mRNA accumulation by 3 hours, consistent with mRNA half-lives. Both methylation marks showed that the epigenetic landscape is highly dynamic and globally remodeled during the 24-hour cycle. While promoters of transcribed genes had tri-methylated H3K4 even at their trough activity times, tri-methylation levels reached their peak, on average, 1 hour after Pol II. Meanwhile, rhythms in tri-methylation of H3K36 lagged transcription by 3 hours. Finally, modeling profiles of Pol II occupancy and mRNA accumulation identified three classes of genes: one showing rhythmicity both in transcriptional and mRNA accumulation, a second class with rhythmic transcription but flat mRNA levels, and a third with constant transcription but rhythmic mRNAs. The latter class emphasizes widespread temporally gated posttranscriptional regulation in the mouse liver.

In mammalian organs such as the liver, many metabolic and physiological processes occur preferentially at specific times during the 24-hour daily cycle. The timing of these rhythmic functions depends on a complex interplay between the endogenous circadian clock and environmental timing cues relayed through the master circadian clock in the suprachiasmatic nucleus, or via feeding rhythms. These rhythms can be implemented on several regulatory levels, and here we aimed at a better understanding of the transcriptional and epigenetic changes that regulate diurnal rhythms. We performed genome-wide analysis of the locations of RNA polymerase II (Pol II) and the epigenetic histone modifications H3K4me3 and H3K36me3 at specific times of day, relating these data to mRNA expression levels. Our analyses show that Pol II transcriptional rhythms are biphasic in mouse liver, having predominant peak activities in the morning and evening. Moreover, dynamic changes in histone marks lag transcription rhythms genome-wide, indicating that the epigenetic landscape can be remodeled during the 24-hour cycle. Finally, a quantitative analysis of temporal Pol II and mRNA accumulation profiles indicates that posttranscriptional regulation significantly contributes to the amplitude and phase of mRNA accumulation profiles. While many studies have analyzed how transcription and chromatin states are modified during irreversible cell differentiation processes, our work highlights how these states can evolve reversibly in a system exhibiting periodicity in time.

Analysis of interferon-beta mRNA stability control after poly(I:C) stimulation using RNA metabolic labeling by ethynyluridine

Interferon-beta (IFN-β) is a critical antiviral cytokine and is essential for innate and acquired immune responses to pathogens. Treatment with polyinosinic:polycytidylic acid (poly(I:C)) induces transient accumulation of IFN-β mRNA, which involves an increase and a decrease of IFN-β mRNA. This phenomenon has been extensively analyzed as a model for understanding the mechanisms of transient gene induction in response to external stimuli. Using a new RNA metabolic labeling method with ethynyluridine to directly measure de novo RNA synthesis and RNA stability, we reassessed both de novo synthesis and degradation of IFN-β mRNA. We found that transcriptional activity is maintained after the maximum accumulation of IFN-β mRNA following poly(I:C) treatment on immortalized human bronchial epithelial cells. We also observed an unexpected change in the stability of IFN-β mRNA before and after the maximum accumulation. The results indicate that this method of RNA metabolic labeling provides a general approach for the simultaneous analysis of transcriptional activity and mRNA stability coupled with transcriptional timing.

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.

The mRNA-bound proteome and its global occupancy profile on protein-coding transcripts

Protein-RNA interactions are fundamental to core biological processes, such as mRNA splicing, localization, degradation, and translation. We developed a photoreactive nucleotide-enhanced UV crosslinking and oligo(dT) purification approach to identify the mRNA-bound proteome using quantitative proteomics and to display the protein occupancy on mRNA transcripts by next-generation sequencing. Application to a human embryonic kidney cell line identified close to 800 proteins. To our knowledge, nearly one-third were not previously annotated as RNA binding, and about 15% were not predictable by computational methods to interact with RNA. Protein occupancy profiling provides a transcriptome-wide catalog of potential cis-regulatory regions on mammalian mRNAs and showed that large stretches in 3' UTRs can be contacted by the mRNA-bound proteome, with numerous putative binding sites in regions harboring disease-associated nucleotide polymorphisms. Our observations indicate the presence of a large number of mRNA binders with diverse molecular functions participating in combinatorial posttranscriptional gene-expression networks.

Minireview: applications of next-generation sequencing on studies of nuclear receptor regulation and function

Next-generation sequencing technologies have expanded the experimental possibilities for studying the genome-wide regulation of transcription by nuclear receptors, their collaborating transcription factors, and coregulators. These technologies allow investigators to obtain abundance and DNA sequence information in a single experiment. In this review, we highlight proven and potential uses of next-generation sequencing in the study of gene regulation by nuclear receptors. We also provide suggestions on how to effectively leverage this technology in a collaborative environment.

C16orf57, a gene mutated in poikiloderma with neutropenia, encodes a putative phosphodiesterase responsible for the U6 snRNA 3' end modification

C16orf57 encodes a human protein of unknown function, and mutations in the gene occur in poikiloderma with neutropenia (PN), which is a rare, autosomal recessive disease. Interestingly, mutations in C16orf57 were also observed among patients diagnosed with Rothmund-Thomson syndrome (RTS) and dyskeratosis congenita (DC), which are caused by mutations in genes involved in DNA repair and telomere maintenance. A genetic screen in Saccharomyces cerevisiae revealed that the yeast ortholog of C16orf57, USB1 (YLR132C), is essential for U6 small nuclear RNA (snRNA) biogenesis and cell viability. Usb1 depletion destabilized U6 snRNA, leading to splicing defects and cell growth defects, which was suppressed by the presence of multiple copies of the U6 snRNA gene SNR6. Moreover, Usb1 is essential for the generation of a unique feature of U6 snRNA; namely, the 3'-terminal phosphate. RNAi experiments in human cells followed by biochemical and functional analyses confirmed that, similar to yeast, C16orf57 encodes a protein involved in the 2',3'-cyclic phosphate formation at the 3' end of U6 snRNA. Advanced bioinformatics predicted that C16orf57 encodes a phosphodiesterase whose putative catalytic activity is essential for its function in vivo. Our results predict an unexpected molecular basis for PN, DC, and RTS and provide insight into U6 snRNA 3' end formation.

Chromatin measurements reveal contributions of synthesis and decay to steady-state mRNA levels

Histone modification, polymerase binding, mRNA half-life, and miRNA abundance measurements in mouse cells are used to dissect the relative contribution of each to mRNA levels, revealing control primarily at the level of transcription, with minor contributions from post-transcriptional processes.

A linear model of three histone modifications and RNAP II occupancy can predict >80% of the variance in mRNA levels.

mRNA half-life explains an additional 1.4% variance in mRNA levels.

miRNA-mediated silencing does not explain any variance on a genome-wide scale.

H3K36me3 has different predictive power in dividing and non-dividingcells.

Messenger RNA levels in eukaryotes are controlled by multiple consecutive regulatory processes, which can be classified into two layers: primary transcriptional regulation at the chromosomal level and secondary, co- and post-transcriptional regulation of the mRNA. To identify the individual contribution of these layers to steady-state RNA levels requires separate quantification. Using mouse as a model organism, we show that chromatin features are sufficient to model RNA levels but with different sensitivities in dividing versus postmitotic cells. In both cases, chromatin-derived transcription rates explain over 80% of the observed variance in measured RNA levels. Further inclusion of measurements of mRNA half-life and microRNA expression data enabled the identification of a low quantitative contribution of RNA decay by either microRNA or general differential turnover to final mRNA levels. Together, this establishes a chromatin-based quantitative model for the contribution of transcriptional and post-transcriptional processes to steady-state levels of messenger RNA.

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.

Transient phenomena in gene expression after induction of transcription

When transcription of a gene is induced by a stimulus, the number of its mRNA molecules changes with time. Here we discuss how this time evolution depends on the shape of the mRNA lifetime distribution. Analysis of the statistical properties of this change reveals transient effects on polysomes, ribosomal profiles, and rate of protein synthesis. Our studies reveal that transient phenomena in gene expression strongly depend on the specific form of the mRNA lifetime distribution.

Computational inference of mRNA stability from histone modification and transcriptome profiles

Histone modifications play important roles in regulating eukaryotic gene expression and have been used to model expression levels. Here, we present a regression model to systematically infer mRNA stability by comparing transcriptome profiles with ChIP-seq of H3K4me3, H3K27me3 and H3K36me3. The results from multiple human and mouse cell lines show that the inferred unstable mRNAs have significantly longer 3′Untranslated Regions (UTRs) and more microRNA binding sites within 3′UTR than the inferred stable mRNAs. Regression residuals derived from RNA-seq, but not from GRO-seq, are highly correlated with the half-lives measured by pulse-labeling experiments, supporting the rationale of our inference. Whereas, the functions enriched in the inferred stable and unstable mRNAs are consistent with those from pulse-labeling experiments, we found the unstable mRNAs have higher cell-type specificity under functional constraint. We conclude that the systematical use of histone modifications can differentiate non-expressed mRNAs from unstable mRNAs, and distinguish stable mRNAs from highly expressed ones. In summary, we represent the first computational model of mRNA stability inference that compares transcriptome and epigenome profiles, and provides an alternative strategy for directing experimental measurements.

Systematic discovery of structural elements governing stability of mammalian messenger RNAs

Decoding post-transcriptional regulatory programs in RNA is a critical step in the larger goal to develop predictive dynamical models of cellular behavior. Despite recent efforts^1–3, the vast landscape of RNA regulatory elements remain largely uncharacterized. A longstanding obstacle is the contribution of local RNA secondary structure in defining interaction partners in a variety of regulatory contexts, including but not limited to transcript stability³, alternative splicing⁴ and localization³. There are many documented instances where the presence of a structural regulatory element dictates alternative splicing patterns (e.g. human cardiac troponin T) or affects other aspects of RNA biology⁵. Thus, a full characterization of post-transcriptional regulatory programs requires capturing information provided by both local secondary structures and the underlying sequence^3,6. We have developed a computational framework based on context-free grammars^3,7 and mutual information² that systematically explores the immense space of small structural elements and reveals motifs that are significantly informative of genome-wide measurements of RNA behavior. The application of this framework to genome-wide mammalian mRNA stability data revealed eight highly significant elements with substantial structural information, for the strongest of which we showed a major role in global mRNA regulation. Through biochemistry, mass-spectrometry, and in vivo binding studies, we identified HNRPA2B1 as the key regulator that binds this element and stabilizes a large number of its target genes. Ultimately, we created a global post-transcriptional regulatory map based on the identity of the discovered linear and structural cis-regulatory elements, their regulatory interactions and their target pathways. This approach can also be employed to reveal the structural elements that modulate other aspects of RNA behavior.

The 3' untranslated region of the rabies virus glycoprotein mRNA specifically interacts with cellular PCBP2 protein and promotes transcript stability

Viral polymerase entry and pausing at intergenic junctions is predicted to lead to a defined polarity in the levels of rhabdovirus gene expression. Interestingly, we observed that the rabies virus glycoprotein mRNA is differentially over-expressed based on this model relative to other transcripts during infection of 293T cells. During infection, the rabies virus glycoprotein mRNA also selectively interacts with the cellular poly(rC)-binding protein 2 (PCBP2), a factor known to influence mRNA stability. Reporter assays performed both in electroporated cells and in a cell-free RNA decay system indicate that the conserved portion of the 3' UTR of the rabies virus glycoprotein mRNA contains an RNA stability element. PCBP2 specifically interacts with reporter transcripts containing this 72 base 3' UTR sequence. Furthermore, the PCBP2 interaction is directly associated with the stability of reporter transcripts. Therefore, we conclude that PCBP2 specifically and selectively interacts with the rabies virus glycoprotein mRNA and that this interaction may contribute to the post-transcriptional regulation of glycoprotein expression.

Genome-wide analysis of long noncoding RNA stability

Transcriptomic analyses have identified tens of thousands of intergenic, intronic, and cis-antisense long noncoding RNAs (lncRNAs) that are expressed from mammalian genomes. Despite progress in functional characterization, little is known about the post-transcriptional regulation of lncRNAs and their half-lives. Although many are easily detectable by a variety of techniques, it has been assumed that lncRNAs are generally unstable, but this has not been examined genome-wide. Utilizing a custom noncoding RNA array, we determined the half-lives of ∼800 lncRNAs and ∼12,000 mRNAs in the mouse Neuro-2a cell line. We find only a minority of lncRNAs are unstable. LncRNA half-lives vary over a wide range, comparable to, although on average less than, that of mRNAs, suggestive of complex metabolism and widespread functionality. Combining half-lives with comprehensive lncRNA annotations identified hundreds of unstable (half-life < 2 h) intergenic, cis-antisense, and intronic lncRNAs, as well as lncRNAs showing extreme stability (half-life > 16 h). Analysis of lncRNA features revealed that intergenic and cis-antisense RNAs are more stable than those derived from introns, as are spliced lncRNAs compared to unspliced (single exon) transcripts. Subcellular localization of lncRNAs indicated widespread trafficking to different cellular locations, with nuclear-localized lncRNAs more likely to be unstable. Surprisingly, one of the least stable lncRNAs is the well-characterized paraspeckle RNA Neat1, suggesting Neat1 instability contributes to the dynamic nature of this subnuclear domain. We have created an online interactive resource (http://stability.matticklab.com) that allows easy navigation of lncRNA and mRNA stability profiles and provides a comprehensive annotation of ∼7200 mouse lncRNAs.

In vivo and transcriptome-wide identification of RNA binding protein target sites

Animal mRNAs are regulated by hundreds of RNA binding proteins (RBPs). The identification of RBP targets is crucial for understanding their function. A recent method, PAR-CLIP, uses photoreactive nucleosides to crosslink RBPs to target RNAs in cells prior to immunoprecipitation. Here, we establish iPAR-CLIP (in vivo PAR-CLIP) to determine, at nucleotide resolution, transcriptome-wide binding sites of GLD-1, a conserved, germline-specific translational repressor in C. elegans. We identified 439 reproducible target mRNAs and demonstrate an excellent dynamic range of target detection by iPAR-CLIP. Upon GLD-1 knockdown, protein but not mRNA expression of the 439 targets was specifically upregulated, demonstrating functionality. Finally, we discovered strongly conserved GLD-1 binding sites near the start codon of target genes. These sites are functional in vitro and likely confer strong repression in vivo. We propose that GLD-1 interacts with the translation machinery near the start codon, a so-far-unknown mode of gene regulation in eukaryotes.

Transcriptional activation of the adenoviral genome is mediated by capsid protein VI

Gene expression of DNA viruses requires nuclear import of the viral genome. Human Adenoviruses (Ads), like most DNA viruses, encode factors within early transcription units promoting their own gene expression and counteracting cellular antiviral defense mechanisms. The cellular transcriptional repressor Daxx prevents viral gene expression through the assembly of repressive chromatin remodeling complexes targeting incoming viral genomes. However, it has remained unclear how initial transcriptional activation of the adenoviral genome is achieved. Here we show that Daxx mediated repression of the immediate early Ad E1A promoter is efficiently counteracted by the capsid protein VI. This requires a conserved PPxY motif in protein VI. Capsid proteins from other DNA viruses were also shown to activate the Ad E1A promoter independent of Ad gene expression and support virus replication. Our results show how Ad entry is connected to transcriptional activation of their genome in the nucleus. Our data further suggest a common principle for genome activation of DNA viruses by counteracting Daxx related repressive mechanisms through virion proteins.

To initiate infection, DNA viruses deliver their genome to the nucleus and express viral genes required for genome replication. Efficient transport is achieved by packing the viral genome as a condensed, transcriptionally inactive nucleo-protein complex. However, for most DNA viruses, including Adenoviruses (Ads), it remains unclear how the viral genome is decondensed and how transcription is initiated inside the nucleus. Cells control unwanted gene expression by chromatin modification mediated through transcriptionally repressive complexes. A key factor in repressive complex assemblies is the transcriptional repressor Daxx. The Ad structural capsid protein VI is required for endosomal escape and nuclear transport. Here we show that protein VI also activates the Ad E1A promoter to initiate Ad gene expression. This is achieved through the removal of Daxx repression from the E1A promoter, which requires a conserved ubiquitin ligase interacting motif (PPxY-motif) in protein VI. We further show that capsid proteins from other unrelated DNA viruses also activate the Ad E1A promoter and support Ad replication by counteracting Daxx repression, functionally replacing protein VI. Our data suggest that reversal of Daxx repression by virion proteins is a widespread mechanism among DNA viruses that is not restricted to a single virus family.

Synergistic and multidimensional regulation of plasminogen activator inhibitor type 1 expression by transforming growth factor type β and epidermal growth factor

The major physiological inhibitor of plasminogen activator, type I plasminogen activator inhibitor (PAI-1), controls blood clotting and tissue remodeling events that involve cell migration. Transforming growth factor type β (TGFβ) and epidermal growth factor (EGF) interact synergistically to increase PAI-1 mRNA and protein levels in human HepG2 and mink Mv1Lu cells. Other growth factors that activate tyrosine kinase receptors can substitute for EGF. EGF and TGFβ regulate PAI-1 by synergistically activating transcription, which is further amplified by a decrease in the rate of mRNA degradation, the latter being regulated only by EGF. The combined effect of transcriptional activation and mRNA stabilization results in a rapid 2-order of magnitude increase in the level of PAI-1. TGFβ also increases the sensitivity of the cells to EGF, thereby recruiting the cooperation of EGF at lower than normally effective concentrations. The contribution of EGF to the regulation of PAI-1 involves the MAPK pathway, and the synergistic interface with the TGFβ pathway is downstream of MEK1/2 and involves phosphorylation of neither ERK1/2 nor Smad2/3. Synergism requires the presence of both Smad and AP-1 recognition sites in the promoter. This work demonstrates the existence of a multidimensional cellular mechanism by which EGF and TGFβ are able to promote large and rapid changes in PAI-1 expression.

Regulation of lymphocyte development and function by RNA-binding proteins

Lymphocyte development requires cells to progress through a series of stages, each associated with changes in gene expression. Intense effort has been invested into characterising the dynamic networks of transcription factors underlying these regulated changes. Whilst transcription factors determine the tempo at which mRNA is produced, recent results highlight the importance of the selective regulation of mRNA decay and translation in regulating gene expression. These processes are regulated by sequence-specific RNA-binding proteins (RBP) as well as noncoding RNA such as microRNAs. RNA-binding proteins are emerging as important regulators of cell fate and function in both developing and mature lymphocytes. At the molecular level the function of RNA-binding proteins is integrated with signal transduction pathways that also govern gene transcription.

Measurement of genome-wide RNA synthesis and decay rates with Dynamic Transcriptome Analysis (DTA)

Standard transcriptomics measures total cellular RNA levels. Our understanding of gene regulation would be greatly improved if we could measure RNA synthesis and decay rates on a genome-wide level. To that end, the Dynamic Transcriptome Analysis (DTA) method has been developed. DTA combines metabolic RNA labeling with standard transcriptomics to measure RNA synthesis and decay rates in a precise and non-perturbing manner. Here, we present the open source R/Bioconductor software package DTA. It implements all required bioinformatics steps that allow the accurate absolute quantification and comparison of RNA turnover.

Genome-wide assessment of AU-rich elements by the AREScore algorithm

In mammalian cells, AU-rich elements (AREs) are well known regulatory sequences located in the 3′ untranslated region (UTR) of many short-lived mRNAs. AREs cause mRNAs to be degraded rapidly and thereby suppress gene expression at the posttranscriptional level. Based on the number of AUUUA pentamers, their proximity, and surrounding AU-rich regions, we generated an algorithm termed AREScore that identifies AREs and provides a numerical assessment of their strength. By analyzing the AREScore distribution in the transcriptomes of 14 metazoan species, we provide evidence that AREs were selected for in several vertebrates and Drosophila melanogaster. We then measured mRNA expression levels genome-wide to address the importance of AREs in SL2 cells derived from D. melanogaster hemocytes. Tis11, a zinc finger RNA–binding protein homologous to mammalian tristetraprolin, was found to target ARE–containing reporter mRNAs for rapid degradation in SL2 cells. Drosophila mRNAs whose expression is elevated upon knock down of Tis11 were found to have higher AREScores. Moreover high AREScores correlate with reduced mRNA expression levels on a genome-wide scale. The precise measurement of degradation rates for 26 Drosophila mRNAs revealed that the AREScore is a very good predictor of short-lived mRNAs. Taken together, this study introduces AREScore as a simple tool to identify ARE–containing mRNAs and provides compelling evidence that AREs are widespread regulatory elements in Drosophila.

Many genes are regulated at the posttranscriptional level by factors that influence the stability of the messenger RNA. In mammals, AU-rich elements are known to cause rapid degradation of messenger RNAs and thereby suppress gene expression. In order to identify such elements on a genome-wide scale, we developed a bioinformatic tool with which we can score messenger RNAs for the presence of AU-rich elements. Using the AREScore algorithm, we observe that AU-rich elements correlate with reduced messenger RNA stability and expression levels. We then used the AREScore to compare the transcriptomes of 14 metazoan species and found that messenger RNAs with high AREScores are enriched in several vertebrates and the fruit fly Drosophila melanogaster. We identified messenger RNAs whose levels are regulated by the Drosophila Tis11 protein, which binds to AU-rich elements. Our study introduces the AREScore as a means to globally assess AU-rich elements and predict short-lived messenger RNAs. Furthermore, it demonstrates the regulatory role of AU-rich elements in suppressing gene expression by accelerating messenger RNA degradation in D. melanogaster cells.

An emerging role of RNA-binding proteins as multifunctional regulators of lymphocyte development and function

Sequence-specific RNA-binding proteins (RBP) and the regulation of RNA decay have long been recognized as important regulators of the inflammatory response. RBP influence gene expression throughout the lifespan of the mRNA by regulating splicing, polyadenylation, cellular localization, translation, and decay. Increasing evidence now indicates that these proteins, together with the RNA decay machinery that they recruit, also regulate the development and activation of lymphocytes. The activity of RBP is regulated by the same signal transduction pathways that govern lymphocyte development and differentiation in response to antigen and cytokine receptor engagement. Roles for these proteins in regulating the diverse functions of lymphocytes are becoming increasingly apparent.

Identification of RNA-protein interaction networks using PAR-CLIP

All mRNA molecules are subject to some degree of post-transcriptional gene regulation (PTGR) involving sequence-dependent modulation of splicing, cleavage and polyadenylation, editing, transport, stability, and translation. The recent introduction of deep-sequencing technologies enabled the development of new methods for broadly mapping interaction sites between RNA-binding proteins (RBPs) and their RNA target sites. In this article, we review crosslinking and immunoprecipitation (CLIP) methods adapted for large-scale identification of target RNA-binding sites and the respective RNA recognition elements. CLIP methods have the potential to detect hundreds of thousands of binding sites in single experiments although the separation of signal from noise can be challenging. As a consequence, each CLIP method has developed different strategies to distinguish true targets from background. We focus on photoactivatable ribonucleoside-enhanced CLIP, which relies on the intracellular incorporation of photoactivatable ribonucleoside analogs into nascent transcripts, and yields characteristic sequence changes upon crosslinking that facilitate the separation of signal from noise. The precise knowledge of the position and distribution of binding sites across mature and primary mRNA transcripts allows critical insights into cellular localization and regulatory function of the examined RBP. When coupled with other systems-wide approaches measuring transcript and protein abundance, the generation of high-resolution RBP-binding site maps across the transcriptome will broaden our understanding of PTGR and thereby lead to new strategies for therapeutic treatment of genetic diseases perturbing these processes.

Strategies to discover regulatory circuits of the mammalian immune system

Recent advances in technologies for genome- and proteome-scale measurements and perturbations promise to accelerate discovery in every aspect of biology and medicine. Although such rapid technological progress provides a tremendous opportunity, it also demands that we learn how to use these tools effectively. One application with great potential to enhance our understanding of biological systems is the unbiased reconstruction of genetic and molecular networks. Cells of the immune system provide a particularly useful model for developing and applying such approaches. Here, we review approaches for the reconstruction of signalling and transcriptional networks, with a focus on applications in the mammalian innate immune system.

3'UTR-mediated gene silencing of the Mixed Lineage Leukemia (MLL) gene

Translocations involving the Mixed Lineage Leukemia (MLL) gene generate in-frame fusions of MLL with more than 50 different partner genes (PGs). Common to all MLL translocations is the exchange not only of coding regions, but also of MLL and PG 3′-untranslated regions (3′UTRs). As a result, the MLL-PG fusion is normally highly expressed and considered the main driver of leukemia development, whereas the function of the PG-MLL fusions in leukemic disease is unclear. As 3′UTRs have been recognized as determinant regions for regulation of gene expression, we hypothesized that loss of the MLL 3′UTR could have a role in generating high MLL-PG levels and leukemia development. Here, we first tested the MLL-PG and PG-MLL mRNA levels in different leukemic cells and tumours and uncovered differential expression that indicates strong repression by the MLL-3′UTR. Reporter assays confirmed that the 3′UTR of MLL, but not of its main PGs, harbours a region that imposes a strong gene silencing effect. Gene suppression by the MLL 3′UTR was largely microRNA independent and did not affect mRNA stability, but inhibited transcription. This effect can at least partially be attributed to a tighter interaction of the MLL 3′UTR with RNA polymerase II than PG 3′UTRs, affecting its phosphorylation state. Altogether, our findings indicate that MLL translocations relieve oncogenic MLL-PG fusions from the repressive MLL 3′UTR, contributing to higher activity of these genes and leukaemia development.

Reversible inhibition of murine cytomegalovirus replication by gamma interferon (IFN-γ) in primary macrophages involves a primed type I IFN-signaling subnetwork for full establishment of an immediate-early antiviral state

Activated macrophages play a central role in controlling inflammatory responses to infection and are tightly regulated to rapidly mount responses to infectious challenge. Type I interferon (alpha/beta interferon [IFN-α/β]) and type II interferon (IFN-γ) play a crucial role in activating macrophages and subsequently restricting viral infections. Both types of IFNs signal through related but distinct signaling pathways, inducing a vast number of interferon-stimulated genes that are overlapping but distinguishable. The exact mechanism by which IFNs, particularly IFN-γ, inhibit DNA viruses such as cytomegalovirus (CMV) is still not fully understood. Here, we investigate the antiviral state developed in macrophages upon reversible inhibition of murine CMV by IFN-γ. On the basis of molecular profiling of the reversible inhibition, we identify a significant contribution of a restricted type I IFN subnetwork linked with IFN-γ activation. Genetic knockout of the type I-signaling pathway, in the context of IFN-γ stimulation, revealed an essential requirement for a primed type I-signaling process in developing a full refractory state in macrophages. A minimal transient induction of IFN-β upon macrophage activation with IFN-γ is also detectable. In dose and kinetic viral replication inhibition experiments with IFN-γ, the establishment of an antiviral effect is demonstrated to occur within the first hours of infection. We show that the inhibitory mechanisms at these very early times involve a blockade of the viral major immediate-early promoter activity. Altogether our results show that a primed type I IFN subnetwork contributes to an immediate-early antiviral state induced by type II IFN activation of macrophages, with a potential further amplification loop contributed by transient induction of IFN-β.

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.