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

Solid phase chemistry to covalently and reversibly capture thiolated RNA

Here, we describe an approach to enrich newly transcribed RNAs from primary mouse neurons using 4-thiouridine (s4U) metabolic labeling and solid phase chemistry. This one-step enrichment procedure captures s4U-RNA by using highly efficient methane thiosulfonate (MTS) chemistry in an immobilized format. Like solution-based methods, this solid-phase enrichment can distinguish mature RNAs (mRNA) with differential stability, and can be used to reveal transient RNAs such as enhancer RNAs (eRNAs) and primary microRNAs (pri-miRNAs) from short metabolic labeling. Most importantly, the efficiency of this solid-phase chemistry made possible the first large scale measurements of RNA polymerase II (RNAPII) elongation rates in mouse cortical neurons. Thus, our approach provides the means to study regulation of RNA metabolism in specific tissue contexts as a means to better understand gene expression in vivo.

Transcriptional Profiling of Quiescent Muscle Stem Cells In Vivo

Muscle stem cells (MuSCs) persist in a quiescent state and activate in response to specific stimuli. The quiescent state is both actively maintained and dynamically regulated. However, analyses of quiescence have come primarily from cells removed from their niche. Although these cells are still quiescent, biochemical changes certainly occur during the isolation process. Here, we analyze the transcriptome of MuSCs in vivo utilizing MuSC-specific labeling of RNA. Notably, labeling transcripts during the isolation procedure revealed very active transcription of specific subsets of genes. However, using the transcription inhibitor α-amanitin, we show that the ex vivo transcriptome remains largely reflective of the in vivo transcriptome. Together, these data provide perspective on the molecular regulation of the quiescent state at the transcriptional level, demonstrate the utility of these tools for probing transcriptional dynamics in vivo, and provide an invaluable resource for understanding stem cell state transitions.

Co-existing feedback loops generate tissue-specific circadian rhythms

The analysis of tissue-specific data-based models of the gene regulatory network of the mammalian circadian clock reveals organ-specific synergies of feedback loops.

Gene regulatory feedback loops generate autonomous circadian rhythms in mammalian tissues. The well-studied core clock network contains many negative and positive regulations. Multiple feedback loops have been discussed as primary rhythm generators but the design principles of the core clock and differences between tissues are still under debate. Here we use global optimization techniques to fit mathematical models to circadian gene expression profiles for different mammalian tissues. It turns out that for every investigated tissue multiple model parameter sets reproduce the experimental data. We extract for all model versions the most essential feedback loops and find auto-inhibitions of period and cryptochrome genes, Bmal1–Rev-erb-α loops, and repressilator motifs as possible rhythm generators. Interestingly, the essential feedback loops differ between tissues, pointing to specific design principles within the hierarchy of mammalian tissue clocks. Self-inhibitions of Per and Cry genes are characteristic for models of suprachiasmatic nucleus clocks, whereas in liver models many loops act in synergy and are connected by a repressilator motif. Tissue-specific use of a network of co-existing synergistic feedback loops could account for functional differences between organs.

Modeling circadian variability of core-clock and clock-controlled genes in four tissues of the rat

Circadian clocks, present in almost all cells of the body, are entrained to rhythmic changes in the environment (e.g. light/dark cycles). Genes responsible for this timekeeping are named core-clock genes, which through transcriptional feedback interactions mediated by transcription factor binding to Ebox/RRE/Dbox elements can generate oscillatory activity of their expression. By regulating the transcription of other clock-controlled genes (CCGs) circadian information is transmitted to tissue and organ levels. Recent studies have indicated that there is a considerable variability of clock-controlled gene expression between tissues both with respect to the circadian genes that are regulated and to their phase lags. In this work, a mathematical model was adapted to explore the dynamics of core-clock and clock-controlled genes measured in four tissues of the rat namely liver, muscle, adipose, and lung. The model efficiently described the synchronous rhythmicity of core-clock genes and further predicted that their phases are mainly regulated by Per2 and Cry1 transcriptional delays and Rev-Erba and Cry1 degradation rates. Similarly, after mining databases for potential Ebox/RRE/Dbox elements in the promoter region of clock-controlled genes, the phase variabilities of the same genes between different tissues were described. The analysis suggests that inter-tissue circadian variability of the same clock-controlled genes is an inherent component of homeostatic function and may arise due to different transcription factor activities on Ebox/RRE/Dbox elements.

Systematic pan-cancer analysis of somatic allele frequency

Imbalanced expression of somatic alleles in cancer can suggest functional and selective features, and can therefore indicate possible driving potential of the underlying genetic variants. To explore the correlation between allele frequency of somatic variants and total gene expression of their harboring gene, we used the unique data set of matched tumor and normal RNA and DNA sequencing data of 5523 distinct single nucleotide variants in 381 individuals across 10 cancer types obtained from The Cancer Genome Atlas (TCGA). We analyzed the allele frequency in the context of the variant and gene functional features and linked it with changes in the total gene expression. We documented higher allele frequency of somatic variants in cancer-implicated genes (Cancer Gene Census, CGC). Furthermore, somatic alleles bearing premature terminating variants (PTVs), when positioned in CGC genes, appeared to be less frequently degraded via nonsense-mediated mRNA decay, indicating possible favoring of truncated proteins by the tumor transcriptome. Among the genes with multiple PTVs with high allele frequency, ARID1, TP53 and NSD1 were known key cancer genes. All together, our analyses suggest that high allele frequency of tumor somatic variants can indicate driving functionality and can serve to identify potential cancer-implicated genes.

Acute Noise Exposure Is Associated With Intrinsic Apoptosis in Murine Central Auditory Pathway

Noise that is capable of inducing the hearing loss (NIHL) has a strong impact on the inner ear structures and causes early and most obvious pathophysiological changes in the auditory periphery. Several studies indicated that intrinsic apoptotic cell death mechanisms are the key factors inducing cellular degeneration immediately after noise exposure and are maintained for days or even weeks. In addition, studies demonstrated several changes in the central auditory system following noise exposure, consistent with early apoptosis-related pathologies. To clarify the underlying mechanisms, the present study focused on the noise-induced gene and protein expression of the pro-apoptotic protease activating factor-1 (APAF1) and the anti-apoptotic B-cell lymphoma 2 related protein a1a (BCL2A1A) in the cochlear nucleus (CN), inferior colliculus (IC) and auditory cortex (AC) of the murine central auditory pathway. The expression of Bcl2a1a mRNA was upregulated immediately after trauma in all tissues investigated, whereas the protein levels were significantly reduced at least in the auditory brainstem. Conversely, acute noise has decreased the expression of Apaf1 gene along the auditory pathway. The changes in APAF1 protein level were not statistically significant. It is tempting to speculate that the acoustic overstimulation leads to mitochondrial dysfunction and induction of apoptosis by regulation of proapoptotic and antiapoptotic proteins. The inverse expression pattern on the mRNA level of both genes might reflect a protective response to decrease cellular damage. Our results indicate the immediate presence of intrinsic apoptosis following noise trauma. This, in turn, may significantly contribute to the development of central structural deficits. Auditory pathway-specific inhibition of intrinsic apoptosis could be a therapeutic approach for the treatment of acute (noise-induced) hearing loss to prevent irreversible neuronal injury in auditory brain structures and to avoid profound deficits in complex auditory processing.

DRUID: a pipeline for transcriptome-wide measurements of mRNA stability

Control of messenger RNA (mRNA) stability is an important aspect of gene regulation. The gold standard for measuring mRNA stability transcriptome-wide uses metabolic labeling, biochemical isolation of labeled RNA populations, and high-throughput sequencing. However, difficult normalization procedures have inhibited widespread adoption of this approach. Here, we present DRUID (for determination of rates using intron dynamics), a new computational pipeline that is robust, easy to use, and freely available. Our pipeline uses endogenous introns to normalize time course data and yields reproducible half-lives, even with data sets that were otherwise unusable. DRUID can handle data sets from a variety of organisms, spanning yeast to humans, and we even applied it retroactively on published data sets. We anticipate that DRUID will allow broad application of metabolic labeling for studies of transcript stability.

A Multi-Parameter Analysis of Cellular Coordination of Major Transcriptome Regulation Mechanisms

To understand cellular coordination of multiple transcriptome regulation mechanisms, we simultaneously measured transcription rate (TR), mRNA abundance (RA) and translation activity (TA). This revealed multiple insights. First, the three parameters displayed systematic statistical differences. Sequentially more genes exhibited extreme (low or high) expression values from TR to RA, and then to TA; that is, cellular coordination of multiple transcriptome regulatory mechanisms leads to sequentially enhanced gene expression selectivity as the genetic information flow from the genome to the proteome. Second, contribution of the stabilization-by-translation regulatory mechanism to the cellular coordination process was assessed. The data enabled an estimation of mRNA stability, revealing a moderate but significant positive correlation between mRNA stability and translation activity. Third, the proportion of mRNA occupied by un-translated regions (UTR) exhibited a negative relationship with the level of this correlation, and was thus a major determinant of the mode of regulation of the mRNA. High-UTR-proportion mRNAs tend to defy the stabilization-by-translation regulatory mechanism, staying out of the polysome but remaining stable; mRNAs with little UTRs largely followed this regulation. In summary, we quantitatively delineated the relationship among multiple transcriptome regulation parameters, i.e., cellular coordination of corresponding regulatory mechanisms.

A mechanistic pan-cancer pathway model informed by multi-omics data interprets stochastic cell fate responses to drugs and mitogens

Most cancer cells harbor multiple drivers whose epistasis and interactions with expression context clouds drug and drug combination sensitivity prediction. We constructed a mechanistic computational model that is context-tailored by omics data to capture regulation of stochastic proliferation and death by pan-cancer driver pathways. Simulations and experiments explore how the coordinated dynamics of RAF/MEK/ERK and PI-3K/AKT kinase activities in response to synergistic mitogen or drug combinations control cell fate in a specific cellular context. In this MCF10A cell context, simulations suggest that synergistic ERK and AKT inhibitor-induced death is likely mediated by BIM rather than BAD, which is supported by prior experimental studies. AKT dynamics explain S-phase entry synergy between EGF and insulin, but simulations suggest that stochastic ERK, and not AKT, dynamics seem to drive cell-to-cell proliferation variability, which in simulations is predictable from pre-stimulus fluctuations in C-Raf/B-Raf levels. Simulations suggest MEK alteration negligibly influences transformation, consistent with clinical data. Tailoring the model to an alternate cell expression and mutation context, a glioma cell line, allows prediction of increased sensitivity of cell death to AKT inhibition. Our model mechanistically interprets context-specific landscapes between driver pathways and cell fates, providing a framework for designing more rational cancer combination therapy.

Cancer is a complex and diverse disease. Two people with the same cancer type often respond differently to the same treatment. These differences are primarily driven by the fact that two type-matched tumors can possess distinct sets of mutations and gene expression profiles, provoking differential sensitivity to drugs. Over the past few decades, we have seen a shift away from more broadly cytotoxic drugs to more targeted molecules therapies; but how to match a patient with a specific drug or drug cocktail remains a difficult problem. Here, we build a mechanistic ordinary differential equation model describing the interactions between commonly mutated pan-cancer signaling pathways—receptor tyrosine kinases, Ras/RAF/ERK, PI3K/AKT, mTOR, cell cycle, DNA damage, and apoptosis. We develop methods for how to tailor the model to multi-omics data from a specific biological context, devise a novel stochastic algorithm to induce non-genetic cell-to-cell fluctuations in mRNA and protein quantities over time, and train the model against a wealth of biochemical and cell fate data to gain insight into the systems-level, context-specific control of proliferation and death. One day, we hope models of this kind could be tailored to patient-derived tumor mRNA sequencing data and used to prioritize patient-specific drug regimens.

Effect of Postmortem Interval and Years in Storage on RNA Quality of Tissue at a Repository of the NIH NeuroBioBank

Brain tissue from 1068 donors was analyzed for RNA quality as a function of postmortem interval (PMI) and years in storage. Approximately 83% of the cortical and cerebellar samples had an RNA integrity number (RIN) of 6 or greater, indicating their likely suitability for real-time quantitative polymerase chain reaction research. The average RIN value was independent of the PMI, up to at least 36 hours. The RNA quality for specific donated brains could not be predicted based on the PMI. Individual samples with a low PMI could have a poor RIN value, while a sample with a PMI over 36 hours may have a high RIN value. The RIN values for control brain donors, all of whom died suddenly and unexpectedly, were marginally higher than for individuals with clinical brain disorders. Polymerase chain reaction (PCR) analysis of samples confirmed that RIN values were more critical than PMI for determining suitability of tissue for molecular biological studies and samples should be matched by their RIN values rather than PMI. Importantly, PCR analysis established that tissue stored up to 23 years at −80°C yielded high-quality RNA. These results confirm that postmortem human brain tissue collected by brain and tissue banks over decades can serve as high quality material for the study of human disorders.

TimeLapse-seq: adding a temporal dimension to RNA sequencing through nucleoside recoding

RNA sequencing (RNA-seq) offers a snapshot of cellular RNA populations, but not temporal information about the sequenced RNA. Here we report TimeLapse-seq, which uses oxidative-nucleophilic-aromatic substitution to convert 4-thiouridine into cytidine analogs, yielding apparent U-to-C mutations that mark new transcripts upon sequencing. TimeLapse-seq is a single-molecule approach that is adaptable to many applications and reveals RNA dynamics and induced differential expression concealed in traditional RNA-seq.

Long noncoding RNA growth arrest-specific 5 promotes proliferation and survival of female germline stem cells in vitro

Female germline stem cells (FGSCs) are proposed to be a key factor for ameliorating female infertility. Previously we have shown that neonatal and adult FGSCs could be isolated and purified from mouse ovarian tissues. The long noncoding (lnc) RNA growth arrest-specific 5 sequence (GAS5) transcribed from mammalian genomes plays important regulatory roles in various developmental processes. However, there is no study on the relationship between GAS5 and FGSC development in vitro. In this study, we showed that GAS5 was highly expressed in the neonatal mouse ovary and was located in both FGSCs and oocytes. GAS5 facilitated FGSC proliferation and promoted their survival in vitro. Moreover, GAS5 also inhibited apoptosis of cultured FGSCs. These findings indicate that GAS5 is a crucial regulator of FGSC development. This might serve as a foundation for a strategy of lncRNA-directed diagnosis or treatment of female infertility.

Circadian clock-dependent and -independent posttranscriptional regulation underlies temporal mRNA accumulation in mouse liver

Rhythms in gene expression propelled by the circadian clock and environmental signals are ubiquitous across cells and tissues. In particular, in mouse tissues, thousands of transcripts show oscillations with a period of 24 hours. Keys question are how such rhythms propagate and eventually exert functions, but also how these are generated. Here, we developed a mathematical model based on total RNA-seq to classify genes according to the respective contributions of transcriptional and posttranscriptional regulation toward mRNA expression profiles. We found that about one-third of rhythmically accumulating mRNA are under posttranscriptional regulation. Such regulation is only partially dependent on the circadian clock, showing that systemic pathways and feeding patterns contribute important posttranscriptional control of gene expression in liver.

The mammalian circadian clock coordinates physiology with environmental cycles through the regulation of daily oscillations of gene expression. Thousands of transcripts exhibit rhythmic accumulations across mouse tissues, as determined by the balance of their synthesis and degradation. While diurnally rhythmic transcription regulation is well studied and often thought to be the main factor generating rhythmic mRNA accumulation, the extent of rhythmic posttranscriptional regulation is debated, and the kinetic parameters (e.g., half-lives), as well as the underlying regulators (e.g., mRNA-binding proteins) are relatively unexplored. Here, we developed a quantitative model for cyclic accumulations of pre-mRNA and mRNA from total RNA-seq data, and applied it to mouse liver. This allowed us to identify that about 20% of mRNA rhythms were driven by rhythmic mRNA degradation, and another 15% of mRNAs regulated by both rhythmic transcription and mRNA degradation. The method could also estimate mRNA half-lives and processing times in intact mouse liver. We then showed that, depending on mRNA half-life, rhythmic mRNA degradation can either amplify or tune phases of mRNA rhythms. By comparing mRNA rhythms in wild-type and Bmal1^−/− animals, we found that the rhythmic degradation of many transcripts did not depend on a functional BMAL1. Interestingly clock-dependent and -independent degradation rhythms peaked at distinct times of day. We further predicted mRNA-binding proteins (mRBPs) that were implicated in the posttranscriptional regulation of mRNAs, either through stabilizing or destabilizing activities. Together, our results demonstrate how posttranscriptional regulation temporally shapes rhythmic mRNA accumulation in mouse liver.

Arabidopsis mRNA decay landscape arises from specialized RNA decay substrates, decapping-mediated feedback, and redundancy

The decay of mRNA plays a vital role in modulating mRNA abundance, which, in turn, influences cellular and organismal processes. In plants and metazoans, three distinct pathways carry out the decay of most cytoplasmic mRNAs: The mRNA decapping complex, which requires the scaffold protein VARICOSE (VCS), removes a protective 5' cap, allowing for 5' to 3' decay via EXORIBONUCLEASE4 (XRN4, XRN1 in metazoans and yeast), and both the exosome and SUPPRESSOR OF VCS (SOV)/DIS3L2 degrade RNAs in the 3' to 5' direction. However, the unique biological contributions of these three pathways, and whether they degrade specialized sets of transcripts, are unknown. In Arabidopsis, the participation of SOV in RNA homeostasis is also unclear, because Arabidopsis sov mutants have a normal phenotype. We carried out mRNA decay analyses in wild-type, sov, vcs, and vcs sov seedlings, and used a mathematical modeling approach to determine decay rates and quantify gene-specific contributions of VCS and SOV to decay. This analysis revealed that VCS (decapping) contributes to decay of 68% of the transcriptome, and, while it initiates degradation of mRNAs with a wide range of decay rates, it especially contributes to decay of short-lived RNAs. Only a few RNAs were clear SOV substrates in that they decayed more slowly in sov mutants. However, 4,506 RNAs showed VCS-dependent feedback in sov that modulated decay rates, and, by inference, transcription, to maintain RNA abundances, suggesting that these RNAs might also be SOV substrates. This feedback was shown to be independent of siRNA activity.

Determining Genome-wide Transcript Decay Rates in Proliferating and Quiescent Human Fibroblasts

Quiescence is a temporary, reversible state in which cells have ceased cell division, but retain the capacity to proliferate. Multiple studies, including ours, have demonstrated that quiescence is associated with widespread changes in gene expression. Some of these changes occur through changes in the level or activity of proliferation-associated transcription factors, such as E2F and MYC. We have demonstrated that mRNA decay can also contribute to changes in gene expression between proliferating and quiescent cells. In this protocol, we describe the procedure for establishing proliferating and quiescent cultures of human dermal foreskin fibroblasts. We then describe the procedures for inhibiting new transcription in proliferating and quiescent cells with Actinomycin D (ActD). ActD treatment represents a straightforward and reproducible approach to dissociating new transcription from transcript decay. A disadvantage of ActD treatment is that the time course must be limited to a short time frame because ActD affects cell viability. Transcript levels are monitored over time to determine transcript decay rates. This procedure allows for the identification of genes and isoforms that exhibit differential decay in proliferating versus quiescent fibroblasts.

The Dynamics of mRNA Turnover Revealed by Single-Molecule Imaging in Single Cells

RNA degradation plays a fundamental role in regulating gene expression. In order to characterize the spatiotemporal dynamics of RNA turnover in single cells, we developed a fluorescent biosensor based on dual-color, single-molecule RNA imaging that allows intact transcripts to be distinguished from stabilized degradation intermediates. Using this method, we measured mRNA decay in single cells and found that individual degradation events occur independently within the cytosol and are not enriched within processing bodies. We show that slicing of an mRNA targeted for endonucleolytic cleavage by the RNA-induced silencing complex can be observed in real time in living cells. This methodology provides a framework for investigating the entire life history of individual mRNAs from birth to death in single cells.

Stress and the nonsense-mediated RNA decay pathway

Cells respond to internal and external cellular stressors by activating stress-response pathways that re-establish homeostasis. If homeostasis is not achieved in a timely manner, stress pathways trigger programmed cell death (apoptosis) to preserve organism integrity. A highly conserved stress pathway is the unfolded protein response (UPR), which senses excessive amounts of unfolded proteins in the ER. While a physiologically beneficial pathway, the UPR requires tight regulation to provide a beneficial outcome and avoid deleterious consequences. Recent work has demonstrated that a conserved and highly selective RNA degradation pathway-nonsense-mediated RNA decay (NMD)-serves as a major regulator of the UPR pathway. NMD degrades mRNAs encoding UPR components to prevent UPR activation in response to innocuous ER stress. In response to strong ER stress, NMD is inhibited by the UPR to allow for a full-magnitude UPR response. Recent studies have indicated that NMD also has other stress-related functions, including promoting the timely termination of the UPR to avoid apoptosis; NMD also regulates responses to non-ER stressors, including hypoxia, amino-acid deprivation, and pathogen infection. NMD regulates stress responses in species across the phylogenetic scale, suggesting that it has conserved roles in shaping stress responses. Stress pathways are frequently constitutively activated or dysregulated in human disease, raising the possibility that "NMD therapy" may provide clinical benefit by downmodulating stress responses.

Alternative exon definition events control the choice between nuclear retention and cytoplasmic export of U11/U12-65K mRNA

Cellular homeostasis of the minor spliceosome is regulated by a negative feed-back loop that targets U11-48K and U11/U12-65K mRNAs encoding essential components of the U12-type intron-specific U11/U12 di-snRNP. This involves interaction of the U11 snRNP with an evolutionarily conserved splicing enhancer giving rise to unproductive mRNA isoforms. In the case of U11/U12-65K, this mechanism controls the length of the 3′ untranslated region (3′UTR). We show that this process is dynamically regulated in developing neurons and some other cell types, and involves a binary switch between translation-competent mRNAs with a short 3′UTR to non-productive isoforms with a long 3′UTR that are retained in the nucleus or/and spliced to the downstream amylase locus. Importantly, the choice between these alternatives is determined by alternative terminal exon definition events regulated by conserved U12- and U2-type 5′ splice sites as well as sequence signals used for pre-mRNA cleavage and polyadenylation. We additionally show that U11 snRNP binding to the U11/U12-65K mRNA species with a long 3′UTR is required for their nuclear retention. Together, our studies uncover an intricate molecular circuitry regulating the abundance of a key spliceosomal protein and shed new light on the mechanisms limiting the export of non-productively spliced mRNAs from the nucleus to the cytoplasm.

The cellular homeostasis of many components of the eukaryotic RNA processing machinery is regulated via negative feed-back pathways that result in the formation of both productive and non-productive mRNA species. Typically, the formation of non-productive mRNAs species results from changes in alternative splicing that disrupt the reading frame of the protein coding region and leads to destabilization of the mRNA. Here, we have investigated the homeostasis regulation of the U11/U12-65K mRNA that encodes an essential protein component of the minor (U12-dependent) spliceosome intron recognition complex. We show that homeostasis is regulated at the level of nuclear mRNA export and mRNA 3′-end formation, and that it can be further regulated during neuronal differentiation. We describe a multilayered regulatory system utilizing alternative exon definition interactions that use the input from both spliceosomes and the polyadenylation machinery to decide between productive and non-productive mRNA formation. Because the 65K protein is an essential component of the minor spliceosome, this regulatory pathway can potentially affect the expression of ~700 genes containing U12-type introns.

The lncRNA VELUCT strongly regulates viability of lung cancer cells despite its extremely low abundance

Little is known about the function of most non-coding RNAs (ncRNAs). The majority of long ncRNAs (lncRNAs) is expressed at very low levels and it is a matter of intense debate whether these can be of functional relevance. Here, we identified lncRNAs regulating the viability of lung cancer cells in a high-throughput RNA interference screen. Based on our previous expression profiling, we designed an siRNA library targeting 638 lncRNAs upregulated in human cancer. In a functional siRNA screen analyzing the viability of lung cancer cells, the most prominent hit was a novel lncRNA which we called Viability Enhancing LUng Cancer Transcript (VELUCT). In silico analyses confirmed the non-coding properties of the transcript. Surprisingly, VELUCT was below the detection limit in total RNA from NCI-H460 cells by RT-qPCR as well as RNA-Seq, but was robustly detected in the chromatin-associated RNA fraction. It is an extremely low abundant lncRNA with an RNA copy number of less than one copy per cell. Blocking transcription with actinomycin D revealed that VELUCT RNA was highly unstable which may partially explain its low steady-state concentration. Despite its extremely low abundance, loss-of-function of VELUCT with three independent experimental approaches in three different lung cancer cell lines led to a significant reduction of cell viability: Next to four individual siRNAs, also two complex siPOOLs as well as two antisense oligonucleotides confirmed the strong and specific phenotype. In summary, the extremely low abundant lncRNA VELUCT is essential for regulation of cell viability in several lung cancer cell lines. Hence, VELUCT is the first example for a lncRNA that is expressed at a very low level, but has a strong loss-of-function phenotype. Thus, our study proves that at least individual low-abundant lncRNAs can play an important functional role.

Rapid Genome-wide Recruitment of RNA Polymerase II Drives Transcription, Splicing, and Translation Events during T Cell Responses

Activation of immune cells results in rapid functional changes, but how such fast changes are accomplished remains enigmatic. By combining time courses of 4sU-seq, RNA-seq, ribosome profiling (RP), and RNA polymerase II (RNA Pol II) ChIP-seq during T cell activation, we illustrate genome-wide temporal dynamics for ∼10,000 genes. This approach reveals not only immediate-early and posttranscriptionally regulated genes but also coupled changes in transcription and translation for >90% of genes. Recruitment, rather than release of paused RNA Pol II, primarily mediates transcriptional changes. This coincides with a genome-wide temporary slowdown in cotranscriptional splicing, even for polyadenylated mRNAs that are localized at the chromatin. Subsequent splicing optimization correlates with increasing Ser-2 phosphorylation of the RNA Pol II carboxy-terminal domain (CTD) and activation of the positive transcription elongation factor (pTEFb). Thus, rapid de novo recruitment of RNA Pol II dictates the course of events during T cell activation, particularly transcription, splicing, and consequently translation.

An immediate-late gene expression module decodes ERK signal duration

The RAF‐MEK‐ERK signalling pathway controls fundamental, often opposing cellular processes such as proliferation and apoptosis. Signal duration has been identified to play a decisive role in these cell fate decisions. However, it remains unclear how the different early and late responding gene expression modules can discriminate short and long signals. We obtained both protein phosphorylation and gene expression time course data from HEK293 cells carrying an inducible construct of the proto‐oncogene RAF. By mathematical modelling, we identified a new gene expression module of immediate–late genes (ILGs) distinct in gene expression dynamics and function. We find that mRNA longevity enables these ILGs to respond late and thus translate ERK signal duration into response amplitude. Despite their late response, their GC‐rich promoter structure suggested and metabolic labelling with 4SU confirmed that transcription of ILGs is induced immediately. A comparative analysis shows that the principle of duration decoding is conserved in PC12 cells and MCF7 cells, two paradigm cell systems for ERK signal duration. Altogether, our findings suggest that ILGs function as a gene expression module to decode ERK signal duration.

DDA1, a novel oncogene, promotes lung cancer progression through regulation of cell cycle

Lung cancer is globally widespread and associated with high morbidity and mortality. DDA1 (DET1 and DDB1 associated 1) was first discovered and registered in the GenBank database by our colleagues. DDA1, an evolutionarily conserved gene, might have significant functions. Recent reports have demonstrated that DDA1 is linked to the ubiquitin–proteasome pathway and facilitates the degradation of target proteins. However, the function of DDA1 in lung cancer was previously unknown. This study aimed to investigate whether DDA1 contributes to tumorigenesis and progression of lung cancer. We found that the expression of DDA1 in normal lung cells and tissue was significantly lower than that in lung cancer and was associated with poor prognosis. DDA1 overexpression promoted proliferation of lung tumour cells and facilitated cell cycle progression in vitro and subcutaneous xenograft tumour progression in vivo. Mechanistically, this was associated with the regulation of S phase and cyclins including cyclin D1/D3/E1. These results indicate that DDA1 promotes lung cancer progression, potentially through promoting cyclins and cell cycle progression. Therefore, DDA1 may be a potential novel target for lung cancer treatment, and a biomarker for tumour prognosis.

Widespread changes in mRNA stability contribute to quiescence-specific gene expression patterns in a fibroblast model of quiescence

Background

Quiescence, reversible exit from the cell division cycle, is characterized by large-scale changes in steady-state gene expression, yet mechanisms controlling these changes are in need of further elucidation. In order to characterize the effects of post-transcriptional control on the quiescent transcriptome in human fibroblasts, we determined mRNA decay rates for over 10,000 genes using a transcription shut-off time-course.

Results

We found that ~500 of the genes monitored exhibited significant changes in decay rate upon quiescence induction. Genes involved in RNA processing and ribosome biogenesis were destabilized with quiescence, while genes involved in the developmental process were stabilized with quiescence. Moreover, extracellular matrix genes demonstrated an upregulation of gene expression that corresponded with a stabilization of these transcripts. Additionally, targets of a quiescence-associated microRNA (miR-29) were significantly enriched in the fraction of transcripts that were stabilized during quiescence.

Conclusion

Coordinated stability changes in clusters of genes with important functions in fibroblast quiescence maintenance are highly correlated with quiescence gene expression patterns. Analysis of miR-29 target decay rates suggests that microRNA-induced changes in RNA stability are important contributors to the quiescence gene expression program in fibroblasts. The identification of multiple stability-related gene clusters suggests that other posttranscriptional regulators of transcript stability may contribute to the coordination of quiescence gene expression. Such regulators may ultimately prove to be valuable targets for therapeutics that target proliferative cells, for instance, in cancer or fibrosis.

Electronic supplementary material

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

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

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

Synthetic recording and in situ readout of lineage information in single cells

Reconstructing the lineage relationships and dynamic event histories of individual cells within their native spatial context is a long-standing challenge in biology. Many biological processes of interest occur in optically opaque or physically inaccessible contexts, necessitating approaches other than direct imaging. Here we describe a synthetic system that enables cells to record lineage information and event histories in the genome in a format that can be subsequently read out of single cells in situ. This system, termed memory by engineered mutagenesis with optical in situ readout (MEMOIR), is based on a set of barcoded recording elements termed scratchpads. The state of a given scratchpad can be irreversibly altered by CRISPR/Cas9-based targeted mutagenesis, and later read out in single cells through multiplexed single-molecule RNA fluorescence hybridization (smFISH). Using MEMOIR as a proof of principle, we engineered mouse embryonic stem cells to contain multiple scratchpads and other recording components. In these cells, scratchpads were altered in a progressive and stochastic fashion as the cells proliferated. Analysis of the final states of scratchpads in single cells in situ enabled reconstruction of lineage information from cell colonies. Combining analysis of endogenous gene expression with lineage reconstruction in the same cells further allowed inference of the dynamic rates at which embryonic stem cells switch between two gene expression states. Finally, using simulations, we show how parallel MEMOIR systems operating in the same cell could enable recording and readout of dynamic cellular event histories. MEMOIR thus provides a versatile platform for information recording and in situ, single-cell readout across diverse biological systems.

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.

Rapid monitoring of RNA degradation activity in vivo for mammalian cells

We have developed a rapid fluorescence assay based on fluorescence resonance energy transfer (FRET) for the monitoring of RNA degradation activity in mammalian cells. In this technique, double-stranded RNA (dsRNA) fluorescent probes are used. The dsRNA fluorescent probes consist of a 5' fluorophore-labeled strand hybridized to a 3' quencher-labeled strand, and the fluorescent dye is quenched by a quencher dye. When the dsRNA is degraded by nascent RNases in cells, the fluorescence emission of the fluorophore is induced following the degradation of the double strands. The degradation rates of the dsRNA are decelerated in response to chemical or environmental toxicity; therefore, in the case of cellular toxicity, the dsRNA is not degraded and remains intact, thus quenching the fluorescence. Unlike in conventional cell-counting assays, this new assay eliminates time-consuming steps, and can be used to simply evaluate the cellular toxicity via a single reaction. Our results demonstrate that this assay can rapidly quantify the RNA degradation rates in vivo within 4 h for three model chemicals. We propose that this assay will be useful for monitoring cellular toxicity in high-throughput applications.

Evaluating the Stability of mRNAs and Noncoding RNAs

Changes in RNA stability have an important impact in the gene expression regulation. Different methods based on the transcription blockage with RNA polymerase inhibitors or metabolic labeling of newly synthesized RNAs have been developed to evaluate RNA decay rates in cultured cell. Combined with techniques to measure transcript abundance genome-wide, these methods have been used to reveal novel features of the eukaryotic transcriptome. The stability of protein-coding mRNAs is in general closely associated to the physiological function of their encoded proteins, with short-lived mRNAs being significantly enriched among regulatory genes whereas genes associated with housekeeping functions are predominantly stable. Likewise, the stability of noncoding RNAs (ncRNAs) seems to reflect their functional role in the cell. Thus, investigating RNA stability can provide insights regarding the function of yet uncharacterized regulatory ncRNAs. In this chapter, we discuss the methodologies currently used to estimate RNA decay and outline an experimental protocol for genome-wide estimation of RNA stability of protein-coding and lncRNAs. This protocol details the transcriptional blockage of cultured cells with actinomycin D, followed by RNA isolation at different time points, the determination of transcript abundance by qPCR/DNA oligoarray hybridization, and the calculation of individual transcript half-lives.

Feedback Loops of the Mammalian Circadian Clock Constitute Repressilator

Mammals evolved an endogenous timing system to coordinate their physiology and behaviour to the 24h period of the solar day. While it is well accepted that circadian rhythms are generated by intracellular transcriptional feedback loops, it is still debated which network motifs are necessary and sufficient for generating self-sustained oscillations. Here, we systematically explore a data-based circadian oscillator model with multiple negative and positive feedback loops and identify a series of three subsequent inhibitions known as "repressilator" as a core element of the mammalian circadian oscillator. The central role of the repressilator motif is consistent with time-resolved ChIP-seq experiments of circadian clock transcription factors and loss of rhythmicity in core clock gene knockouts.

Control of Gene Expression by RNA Binding Protein Action on Alternative Translation Initiation Sites

Transcript levels do not faithfully predict protein levels, due to post-transcriptional regulation of gene expression mediated by RNA binding proteins (RBPs) and non-coding RNAs. We developed a multivariate linear regression model integrating RBP levels and predicted RBP-mRNA regulatory interactions from matched transcript and protein datasets. RBPs significantly improved the accuracy in predicting protein abundance of a portion of the total modeled mRNAs in three panels of tissues and cells and for different methods employed in the detection of mRNA and protein. The presence of upstream translation initiation sites (uTISs) at the mRNA 5’ untranslated regions was strongly associated with improvement in predictive accuracy. On the basis of these observations, we propose that the recently discovered widespread uTISs in the human genome can be a previously unappreciated substrate of translational control mediated by RBPs.

Gene expression is a dynamic program by which the information stored in the genome is rendered functional by production and degradation of two types of macromolecules, RNAs and proteins. mRNAs are templates for proteins; therefore we expect correspondence between quantities of mRNAs and proteins. Genome-wide studies instead indicate a marked discrepancy between them, when considering their steady-state levels or their variations across different conditions. We employed linear regression approaches with paired mRNA/protein datasets in order to develop a model predicting the protein level of a gene from both the mRNA level and the protein levels of RBPs inferred to bind the mRNA untranslated regions. The results of our analyses restricted the utility of RBPs to improve accuracy of predicted protein abundance to a small fraction of the total modelled genes, and identified a novel association of the improvement induced by RBPs with the presence of upstream translation sites. This finding suggests a new avenue of experimental studies aimed at exploring the hypothesis that RBPs could influence protein abundance by changing the preference for certain translation initiation sites.

Post-transcriptional regulation of cytokine and growth factor signaling in cancer

Cytokines and growth factors regulate cell proliferation, differentiation, migration and apoptosis, and play important roles in coordinating growth signal responses during development. The expression of cytokine genes and the signals transmitted through cytokine receptors are tightly regulated at several levels, including transcriptional and post-transcriptional levels. A majority of cytokine mRNAs, including growth factor transcripts, contain AU-rich elements (AREs) in their 3' untranslated regions that control gene expression by regulating mRNA degradation and changing translational rates. In addition, numerous proteins involved in transmitting signals downstream of cytokine receptors are regulated at the level of mRNA degradation by GU-rich elements (GREs) found in their 3' untranslated regions. Abnormal stabilization and overexpression of ARE or GRE-containing transcripts had been observed in many malignancies, which is a consequence of the malfunction of RNA-binding proteins. In this review, we briefly summarize the role of AREs and GREs in regulating mRNA turnover to coordinate cytokine and growth factor expression, and we describe how dysregulation of mRNA degradation mechanisms contributes to the development and progression of cancer.

Systematic Analysis of Drug Targets Confirms Expression in Disease-Relevant Tissues

It is commonly assumed that drug targets are expressed in tissues relevant to their indicated diseases, even under normal conditions. While multiple anecdotal cases support this hypothesis, a comprehensive study has not been performed to verify it. We conducted a systematic analysis to assess gene and protein expression for all targets of marketed and phase III drugs across a diverse collection of normal human tissues. For 87% of gene-disease pairs, the target is expressed in a disease-affected tissue under healthy conditions. This result validates the importance of confirming expression of a novel drug target in an appropriate tissue for each disease indication and strengthens previous findings showing that targets of efficacious drugs should be expressed in relevant tissues under normal conditions. Further characterization of the remaining 13% of gene-disease pairs revealed that most genes are expressed in a different tissue linked to another disease. Our analysis demonstrates the value of extensive tissue specific expression resources.both in terms of tissue and cell diversity as well as techniques used to measure gene expression.

Cytoplasmic RNA decay pathways - Enzymes and mechanisms

RNA decay plays a crucial role in post-transcriptional regulation of gene expression. Work conducted over the last decades has defined the major mRNA decay pathways, as well as enzymes and their cofactors responsible for these processes. In contrast, our knowledge of the mechanisms of degradation of non-protein coding RNA species is more fragmentary. This review is focused on the cytoplasmic pathways of mRNA and ncRNA degradation in eukaryotes. The major 3' to 5' and 5' to 3' mRNA decay pathways are described with emphasis on the mechanisms of their activation by the deprotection of RNA ends. More recently discovered 3'-end modifications such as uridylation, and their relevance to cytoplasmic mRNA decay in various model organisms, are also discussed. Finally, we provide up-to-date findings concerning various pathways of non-coding RNA decay in the cytoplasm.

The rules and impact of nonsense-mediated mRNA decay in human cancers

Premature termination codons (PTCs) cause a large proportion of inherited human genetic diseases. PTC-containing transcripts can be degraded by an mRNA surveillance pathway termed nonsense-mediated mRNA decay (NMD). However, the efficiency of NMD varies; it is inefficient when a PTC is located downstream of the last exon junction complex (EJC). We used matched exome and transcriptome data from 9,769 human tumors to systematically elucidate the rules of NMD targeting in human cells. An integrated model incorporating multiple rules beyond the canonical EJC model explains approximately three-quarters of the non-random variance in NMD efficiency across thousands of PTCs. We also show that dosage compensation may mask the effects of NMD. Applying the NMD model identifies signatures of both positive and negative selection on NMD-triggering mutations in human tumors and provides a classification of tumor suppressor genes.

Linking Core Promoter Classes to Circadian Transcription

Circadian rhythms in transcription are generated by rhythmic abundances and DNA binding activities of transcription factors. Propagation of rhythms to transcriptional initiation involves the core promoter, its chromatin state, and the basal transcription machinery. Here, I characterize core promoters and chromatin states of genes transcribed in a circadian manner in mouse liver and in Drosophila. It is shown that the core promoter is a critical determinant of circadian mRNA expression in both species. A distinct core promoter class, strong circadian promoters (SCPs), is identified in mouse liver but not Drosophila. SCPs are defined by specific core promoter features, and are shown to drive circadian transcriptional activities with both high averages and high amplitudes. Data analysis and mathematical modeling further provided evidence for rhythmic regulation of both polymerase II recruitment and pause release at SCPs. The analysis provides a comprehensive and systematic view of core promoters and their link to circadian mRNA expression in mouse and Drosophila, and thus reveals a crucial role for the core promoter in regulated, dynamic transcription.

Functions and mechanisms of long noncoding RNAs in lung cancer

Lung cancer is a heterogeneous disease, and there is a lack of adequate biomarkers for diagnosis. Long noncoding RNAs (lncRNAs) are emerging as an important set of molecules because of their roles in various key pathophysiological pathways, including cell growth, apoptosis, and metastasis. We review the current knowledge of the lncRNAs in lung cancer. In-depth analyses of lncRNAs in lung cancer have increased the number of potential effective biomarkers, thus providing options to increase the therapeutic benefit. In this review, we summarize the functions, mechanisms, and regulatory networks of lncRNAs in lung cancer, providing a basis for further research in this field.

The SLE variant Ala71Thr of BLK severely decreases protein abundance and binding to BANK1 through impairment of the SH3 domain function

The B-lymphocyte kinase (BLK) gene is associated genetically with several human autoimmune diseases including systemic lupus erythematosus. We recently described that the genetic risk is given by two haplotypes: one covering several strongly linked single-nucleotide polymorphisms within the promoter of the gene that correlated with low transcript levels, and a second haplotype that includes a rare nonsynonymous variant (Ala71Thr). Here we show that this variant, located within the BLK SH3 domain, is a major determinant of protein levels. In vitro analyses show that the 71Thr isoform is hyperphosphorylated and promotes kinase activation. As a consequence, BLK is ubiquitinated, its proteasomal degradation enhanced and the average life of the protein is reduced by half. Altogether, these findings suggest that an intrinsic autoregulatory mechanism previously unappreciated in BLK is disrupted by the 71Thr substitution. Because the SH3 domain is also involved in protein interactions, we sought for differences between the two isoforms in trafficking and binding to protein partners. We found that binding of the 71Thr variant to the adaptor protein BANK1 is severely reduced. Our study provides new insights on the intrinsic regulation of BLK activation and highlights the dominant role of its SH3 domain in BANK1 binding.

Degradation dynamics of microRNAs revealed by a novel pulse-chase approach

The regulation of miRNAs is critical to the definition of cell identity and behavior in normal physiology and disease. To date, the dynamics of miRNA degradation and the mechanisms involved in remain largely obscure, in particular, in higher organisms. Here, we developed a pulse-chase approach based on metabolic RNA labeling to calculate miRNA decay rates at genome-wide scale in mammalian cells. Our analysis revealed heterogeneous miRNA half-lives, with many species behaving as stable molecules (T_1/2 > 24 h), while others, including passenger miRNAs and a number (25/129) of guide miRNAs, are quickly turned over (T_1/2 = 4–14 h). Decay rates were coupled with other features, including genomic organization, transcription rates, structural heterogeneity (isomiRs), and target abundance, measured through quantitative experimental approaches. This comprehensive analysis highlighted functional mechanisms that mediate miRNA degradation, as well as the importance of decay dynamics in the regulation of the miRNA pool under both steady-state conditions and during cell transitions.

Genetic oscillators in development

In development, morphogenetic processes are strictly coordinated in time. Cells in a developing tissue would need mechanisms for time-keeping. One such time-keeping mechanism is to use oscillations of gene expression. Oscillatory gene expression can be generated by transcriptional/translational feedback loops, usually referred to as a genetic oscillator. In this review article, we discuss genetic oscillators in the presence of developmental processes such as cell division, cell movement and cell differentiation. We first introduce the gene regulatory network for generating a rhythm of gene expression. We then discuss how developmental processes influence genetic oscillators. Examples include vertebrate somitogenesis and neural progenitor cell differentiation, as well as the circadian clock for comparison. To understand the behaviors of genetic oscillators in development, it is necessary to consider both gene expression dynamics and cellular behaviors simultaneously. Theoretical modeling combined with live imaging at single-cell resolution will be a powerful tool to analyze genetic oscillators in development.

Circular RNA enrichment in platelets is a signature of transcriptome degradation

In platelets, splicing and translation occur in the absence of a nucleus. However, the integrity and stability of mRNAs derived from megakaryocyte progenitor cells remain poorly quantified on a transcriptome-wide level. As circular RNAs (circRNAs) are resistant to degradation by exonucleases, their abundance relative to linear RNAs can be used as a surrogate marker for mRNA stability in the absence of transcription. Here we show that circRNAs are enriched in human platelets 17- to 188-fold relative to nucleated tissues and 14- to 26-fold relative to samples digested with RNAse R to selectively remove linear RNA. We compare RNAseq read depths inside and outside circRNAs to provide in silico evidence of transcript circularity, show that exons within circRNAs are enriched on average 12.7 times in platelets relative to nucleated tissues and identify 3162 genes significantly enriched for circRNAs, including some where all RNAseq reads appear to be derived from circular molecules. We also confirm that this is a feature of other anucleate cells through transcriptome sequencing of mature erythrocytes, demonstrate that circRNAs are not enriched in cultured megakaryocytes, and demonstrate that linear RNAs decay more rapidly than circRNAs in platelet preparations. Collectively, these results suggest that circulating platelets have lost >90% of their progenitor mRNAs and that translation in platelets occurs against the backdrop of a highly degraded transcriptome. Finally, we find that transcripts previously classified as products of reverse transcriptase template switching are both enriched in platelets and resistant to decay, countering the recent suggestion that up to 50% of rearranged RNAs are artifacts.

Circular RNAs are long-lived and display only minimal early alterations in response to a growth factor

Circular RNAs (circRNAs) are widespread circles of non-coding RNAs with largely unknown function. Because stimulation of mammary cells with the epidermal growth factor (EGF) leads to dynamic changes in the abundance of coding and non-coding RNA molecules, and culminates in the acquisition of a robust migratory phenotype, this cellular model might disclose functions of circRNAs. Here we show that circRNAs of EGF-stimulated mammary cells are stably expressed, while mRNAs and microRNAs change within minutes. In general, the circRNAs we detected are relatively long-lived and weakly expressed. Interestingly, they are almost ubiquitously co-expressed with the corresponding linear transcripts, and the respective, shared promoter regions are more active compared to genes producing linear isoforms with no detectable circRNAs. These findings imply that altered abundance of circRNAs, unlike changes in the levels of other RNAs, might not play critical roles in signaling cascades and downstream transcriptional networks that rapidly commit cells to specific outcomes.

The whole-genome expression analysis of peripheral blood mononuclear cells from aspirin sensitive asthmatics versus aspirin tolerant patients and healthy donors after in vitro aspirin challenge

Background

Up to 30 % of adults with severe asthma are hypersensitive to aspirin and no unambiguous theory exists which provides a satisfactory explanation for the occurrence of aspirin-induced asthma (AIA) in some asthmatic patients. Therefore, the aim of this study was to compare the AIA expression profile against aspirin tolerant asthma (ATA) and healthy volunteers (HV) profile in peripheral blood mononuclear cells (PBMCs) after in vitro aspirin challenge in Caucasian population.

Methods

PBMCs were separated from blood of three groups of subjects - 11 AIA, 7 ATA and 15 HV and then stimulated by either 2 μM lysine aspirin or 20 μM lysine as a control. Subsequently, RNA was isolated, transcribed into cDNA and subjected to microarray and qPCR studies. Simultaneously, protein was extracted from PBMCs and used in further immunoblotting analysis.

Results

The validation of results at mRNA level has shown only three genes, whose expression was significantly altered between comprising groups. mRNA expression of CNPY3 in PBMCs in AIA was significantly lower (-0.41 ± 2.67) than in HV (1.04 ± 2.69), (p = 0.02); mRNA expression of FOSL1 in PBMCs in AIA was also significantly decreased (-0.66 ± 2.97) as opposed to HV (0.31 ± 4.83), (p = 0.02). While mRNA expression of ERAS in PBMCs was increased (1.15 ± 0.23) in AIA in comparison to HV (-1.32 ± 0.41), (p = 0.03). At protein level the changed expression of one protein was confirmed. Protein expression of FOSL1 in PBMCs in AIA was both significantly lower (-0.86 ± 0.08) than in ATA (0.39 ± 0.42), (p = 0.046) and in HV (0.9 ± 0.27), (p = 0.007).

Conclusions

This pilot study implies a positive association between CNPY3, ERAS, FOSL1 and aspirin-intolerant asthma, suggesting that these findings would be useful for further investigations of NSAIDs mechanism.

Electronic supplementary material

The online version of this article (doi:10.1186/s12931-015-0305-4) contains supplementary material, which is available to authorized users.

Functional coordination and HuR-mediated regulation of mRNA stability during T cell activation

Global mRNA abundance depends on the balance of synthesis and decay of a population of mRNAs. To account for this balance during activation of T cells, we used metabolic labeling to quantify the contributions of RNA transcription and decay over a 4 h time course during activation of leukemia-derived Jurkat T cells. While prior studies suggested more than half of the changes in mRNA abundance were due to RNA stability, we found a smaller but more interesting population of mRNAs changed stability. These mRNAs clustered into functionally related subpopulations that included replicative histones, ribosomal biogenesis and cell motility functions. We then applied a novel analysis based on integrating global protein-RNA binding with concurrent changes in RNA stability at specific time points following activation. This analysis demonstrated robust stabilization of mRNAs by the HuR RNA-binding protein 4 h after activation. Our unexpected findings demonstrate that the temporal regulation of mRNA stability coordinates vital cellular pathways and is in part controlled by the HuR RNA binding protein in Jurkat T cells following activation.

Next-generation analysis of gene expression regulation--comparing the roles of synthesis and degradation

Technological advances now enable routine measurement of mRNA and protein abundances, and estimates of their rates of synthesis and degradation that inform on their values and the degree of change in response to stimuli. Importantly, more and more data on time-series experiments are emerging, e.g. of cells responding to stress, enabling first insights into a new dimension of gene expression regulation - its dynamics and how it allows for very different response signals across genes. This review discusses recently published methods and datasets, their impact on what we now know about the relationships between concentrations and synthesis rates of mRNAs and proteins in yeast and mammalian cells, their evolution, and new hypotheses on translation regulatory mechanisms generated by approaches that involve ribosome footprinting.

mRIN for direct assessment of genome-wide and gene-specific mRNA integrity from large-scale RNA-sequencing data

The volume of RNA-Seq data sets in public repositories has been expanding exponentially, providing unprecedented opportunities to study gene expression regulation. Because degraded RNA samples, such as those collected from post-mortem tissues, can result in distinct expression profiles with potential biases, a particularly important step in mining these data is quality control. Here we develop a method named mRIN to directly assess mRNA integrity from RNA-Seq data at the sample and individual gene level. We systematically analyse large-scale RNA-Seq data sets of the human brain transcriptome generated by different consortia. Our analysis demonstrates that 3′ bias resulting from partial RNA fragmentation in post-mortem tissues has a marked impact on global expression profiles, and that mRIN effectively identifies samples with different levels of mRNA degradation. Unexpectedly, this process has a reproducible and gene-specific component, and transcripts with different stabilities are associated with distinct functions and structural features reminiscent of mRNA decay in living cells.

With the rapid increase in the volume of publically available RNA-seq data, quality control is an increasingly important consideration. Here Feng et al. develop mRIN, a method to directly assess mRNA integrity, and show that RNA degradation in post-mortem samples has a strong impact on global expression profiles.

Rapid proteasomal elimination of 3-hydroxy-3-methylglutaryl-CoA reductase by interferon-γ in primary macrophages requires endogenous 25-hydroxycholesterol synthesis

Interferons (IFNs) play a central role in immunity and emerging evidence suggests that IFN-signalling coordinately regulates sterol biosynthesis in macrophages, via Sterol Regulatory Element-Binding Protein (SREBP) dependent and independent pathways. However, the precise mechanisms and kinetic steps by which IFN controls sterol biosynthesis are as yet not fully understood. Here, we elucidate the molecular circuitry governing how IFN controls the first regulated step in the mevalonate-sterol pathway, 3-hydroxy-3-methylglutaryl-CoA reductase (HMGCR), through the synthesis of 25-Hydroxycholesterol (25-HC) from cholesterol by the IFN-inducible Cholesterol-25-Hydroxylase (CH25H). We show for the first 30-min of IFN stimulation of macrophages the rate of de novo synthesis of the Ch25h transcript is markedly increased but by 120-min becomes transcriptionally curtailed, coincident with induction of the Activating Transcription Factor 3 (ATF3) repressor. We demonstrate ATF3 induction by Toll-like receptors is strictly dependent on IFN-signalling. While the SREBP-pathway dependent rates of de novo transcription of Hmgcr are relatively unchanged in the first 90-min of IFN treatment, we find HMGCR enzyme levels undergo a rapid proteasomal-mediated degradation, defining a previously unappreciated SREBP-independent mechanism for IFN-action. These events precede a sustained marked reduction in Hmgcr RNA levels involving SREBP-dependent mechanisms. We demonstrate that HMGCR proteasomal-degradation by IFN strictly requires the synthesis of endogenous 25-HC and functionally couples HMGCR to CH25H to coordinately suppress sterol biosynthesis. In conclusion, we quantitatively delineate proteomic and transcriptional levels of IFN-mediated control of HMGCR, the primary enzymatic step of the mevalonate-sterol biosynthesis pathway, providing a foundational framework for mathematically modelling the therapeutic outcome of immune-metabolic pathways.

Dynamic regulation of mRNA decay during neural development

Background

Gene expression patterns are determined by rates of mRNA transcription and decay. While transcription is known to regulate many developmental processes, the role of mRNA decay is less extensively defined. A critical step toward defining the role of mRNA decay in neural development is to measure genome-wide mRNA decay rates in neural tissue. Such information should reveal the degree to which mRNA decay contributes to differential gene expression and provide a foundation for identifying regulatory mechanisms that affect neural mRNA decay.

Results

We developed a technique that allows genome-wide mRNA decay measurements in intact Drosophila embryos, across all tissues and specifically in the nervous system. Our approach revealed neural-specific decay kinetics, including stabilization of transcripts encoding regulators of axonogenesis and destabilization of transcripts encoding ribosomal proteins and histones. We also identified correlations between mRNA stability and physiologic properties of mRNAs; mRNAs that are predicted to be translated within axon growth cones or dendrites have long half-lives while mRNAs encoding transcription factors that regulate neurogenesis have short half-lives. A search for candidate cis-regulatory elements identified enrichment of the Pumilio recognition element (PRE) in mRNAs encoding regulators of neurogenesis. We found that decreased expression of the RNA-binding protein Pumilio stabilized predicted neural mRNA targets and that a PRE is necessary to trigger reporter-transcript decay in the nervous system.

Conclusions

We found that differential mRNA decay contributes to the relative abundance of transcripts involved in cell-fate decisions, axonogenesis, and other critical events during Drosophila neural development. Neural-specific decay kinetics and the functional specificity of mRNA decay suggest the existence of a dynamic neurodevelopmental mRNA decay network. We found that Pumilio is one component of this network, revealing a novel function for this RNA-binding protein.

Electronic supplementary material

The online version of this article (doi:10.1186/s13064-015-0038-6) contains supplementary material, which is available to authorized users.

Post-transcriptional regulation of cytokine signaling by AU-rich and GU-rich elements

Cytokines are necessary for cell communication to enable responses to external stimuli that are imperative for the survival and maintenance of homeostasis. Dysfunction of the cytokine network has detrimental effects on intra- and extracellular environments. Thus, it is critical that the expression of cytokines and the signals transmitted by cytokines to target cells are tightly regulated at numerous levels, including transcriptional and post-transcriptional levels. Here, we briefly summarize the role of AU-rich elements (AREs) in the regulation of cytokine gene expression at the post-transcriptional level and describe a role for GU-rich elements (GREs) in coordinating the regulation of cytokine signaling. GREs function as post-transcriptional regulators of proteins that control cellular activation, growth, and apoptosis. GREs and AREs work in concert to coordinate cytokine signal transduction pathways. The precise regulation of cytokine signaling is particularly important, because its dysregulation can lead to human diseases.

Increased expression of X-linked genes in mammals is associated with a higher stability of transcripts and an increased ribosome density

Mammalian sex chromosomes evolved from the degeneration of one homolog of a pair of ancestral autosomes, the proto-Y. This resulted in a gene dose imbalance that is believed to be restored (partially or fully) through upregulation of gene expression from the single active X-chromosome in both sexes by a dosage compensatory mechanism. We analyzed multiple genome-wide RNA stability data sets and found significantly longer average half-lives for X-chromosome transcripts than for autosomal transcripts in various human cell lines, both male and female, and in mice. Analysis of ribosome profiling data shows that ribosome density is higher on X-chromosome transcripts than on autosomal transcripts in both humans and mice, suggesting that the higher stability is causally linked to a higher translation rate. Our results and observations are in accordance with a dosage compensatory upregulation of expressed X-linked genes. We therefore propose that differential mRNA stability and translation rates of the autosomes and sex chromosomes contribute to an evolutionarily conserved dosage compensation mechanism in mammals.