Genome-wide determination of RNA stability reveals hundreds of short-lived noncoding transcripts in mammals

Isoform and pathway-specific regulation of post-transcriptional RNA processing in human cells

Steady-state levels of RNA transcripts are controlled by their rates of synthesis and degradation. Here we used nascent RNA Bru-seq and BruChase-seq to profile RNA dynamics across 16 human cell lines as part of ENCODE4 Deeply Profiled Cell Lines collection. We show that RNA turnover dynamics differ widely between transcripts of different genes and between different classes of RNA. Gene set enrichment analysis (GSEA) revealed that transcripts encoding proteins belonging to the same pathway often show similar turnover dynamics. Furthermore, transcript isoforms show distinct dynamics suggesting that RNA turnover is important in regulating mRNA isoform choice. Finally, splicing across newly made transcripts appears to be cooperative with either all or none type splicing. These data sets generated as part of ENCODE4 illustrate the intricate and coordinated regulation of RNA dynamics in controlling gene expression to allow for the precise coordination of cellular functions.

Cellular control of protein levels: A systems biology perspective

How cells regulate protein levels is a central question of biology. Over the past decades, molecular biology research has provided profound insights into the mechanisms and the molecular machinery governing each step of the gene expression process, from transcription to protein degradation. Recent advances in transcriptomics and proteomics have complemented our understanding of these fundamental cellular processes with a quantitative, systems-level perspective. Multi-omic studies revealed significant quantitative, kinetic and functional differences between the genome, transcriptome and proteome. While protein levels often correlate with mRNA levels, quantitative investigations have demonstrated a substantial impact of translation and protein degradation on protein expression control. In addition, protein-level regulation appears to play a crucial role in buffering protein abundances against undesirable mRNA expression variation. These findings have practical implications for many fields, including gene function prediction and precision medicine.

Nascent transcription reveals regulatory changes in extremophile fishes inhabiting hydrogen sulfide-rich environments

Regulating transcription allows organisms to respond to their environment, both within a single generation (plasticity) and across generations (adaptation). We examined transcriptional differences in gill tissues of fishes in the Poecilia mexicana species complex (family Poeciliidae), which have colonized toxic springs rich in hydrogen sulfide (H2S) in southern Mexico. There are gene expression differences between sulfidic and non-sulfidic populations, yet regulatory mechanisms mediating this gene expression variation remain poorly studied. We combined capped-small RNA sequencing (csRNA-seq), which captures actively transcribed (i.e. nascent) transcripts, and messenger RNA sequencing (mRNA-seq) to examine how variation in transcription, enhancer activity, and associated transcription factor binding sites may facilitate adaptation to extreme environments. csRNA-seq revealed thousands of differentially initiated transcripts between sulfidic and non-sulfidic populations, many of which are involved in H2S detoxification and response. Analyses of transcription factor binding sites in promoter and putative enhancer csRNA-seq peaks identified a suite of transcription factors likely involved in regulating H2S-specific shifts in gene expression, including several key transcription factors known to respond to hypoxia. Our findings uncover a complex interplay of regulatory processes that reflect the divergence of extremophile populations of P. mexicana from their non-sulfidic ancestors and suggest shared responses among evolutionarily independent lineages.

Global, site-resolved analysis of ubiquitylation occupancy and turnover rate reveals systems properties

Ubiquitylation regulates most proteins and biological processes in a eukaryotic cell. However, the site-specific occupancy (stoichiometry) and turnover rate of ubiquitylation have not been quantified. Here we present an integrated picture of the global ubiquitylation site occupancy and half-life. Ubiquitylation site occupancy spans over four orders of magnitude, but the median ubiquitylation site occupancy is three orders of magnitude lower than that of phosphorylation. The occupancy, turnover rate, and regulation of sites by proteasome inhibitors are strongly interrelated, and these attributes distinguish sites involved in proteasomal degradation and cellular signaling. Sites in structured protein regions exhibit longer half-lives and stronger upregulation by proteasome inhibitors than sites in unstructured regions. Importantly, we discovered a surveillance mechanism that rapidly and site-indiscriminately deubiquitylates all ubiquitin-specific E1 and E2 enzymes, protecting them against accumulation of bystander ubiquitylation. The work provides a systems-scale, quantitative view of ubiquitylation properties and reveals general principles of ubiquitylation-dependent governance.

Cis-regulatory control of transcriptional timing and noise in response to estrogen

Cis-regulatory elements control transcription levels, temporal dynamics, and cell-cell variation or transcriptional noise. However, the combination of regulatory features that control these different attributes is not fully understood. Here, we used single-cell RNA-seq during an estrogen treatment time course and machine learning to identify predictors of expression timing and noise. We found that genes with multiple active enhancers exhibit faster temporal responses. We verified this finding by showing that manipulation of enhancer activity changes the temporal response of estrogen target genes. Analysis of transcriptional noise uncovered a relationship between promoter and enhancer activity, with active promoters associated with low noise and active enhancers linked to high noise. Finally, we observed that co-expression across single cells is an emergent property associated with chromatin looping, timing, and noise. Overall, our results indicate a fundamental tradeoff between a gene's ability to quickly respond to incoming signals and maintain low variation across cells.

Annotation of nuclear lncRNAs based on chromatin interactions

The human genome is pervasively transcribed and produces a wide variety of long non-coding RNAs (lncRNAs), constituting the majority of transcripts across human cell types. Some specific nuclear lncRNAs have been shown to be important regulatory components acting locally. As RNA-chromatin interaction and Hi-C chromatin conformation data showed that chromatin interactions of nuclear lncRNAs are determined by the local chromatin 3D conformation, we used Hi-C data to identify potential target genes of lncRNAs. RNA-protein interaction data suggested that nuclear lncRNAs act as scaffolds to recruit regulatory proteins to target promoters and enhancers. Nuclear lncRNAs may therefore play a role in directing regulatory factors to locations spatially close to the lncRNA gene. We provide the analysis results through an interactive visualization web portal at https://fantom.gsc.riken.jp/zenbu/reports/#F6_3D_lncRNA.

Time-resolved profiling of RNA binding proteins throughout the mRNA life cycle

mRNAs continually change their protein partners throughout their lifetimes, yet our understanding of mRNA-protein complex (mRNP) remodeling is limited by a lack of temporal data. Here, we present time-resolved mRNA interactome data by performing pulse metabolic labeling with photoactivatable ribonucleoside in human cells, UVA crosslinking, poly(A)+ RNA isolation, and mass spectrometry. This longitudinal approach allowed the quantification of over 700 RNA binding proteins (RBPs) across ten time points. Overall, the sequential order of mRNA binding aligns well with known functions, subcellular locations, and molecular interactions. However, we also observed RBPs with unexpected dynamics: the transcription-export (TREX) complex recruited posttranscriptionally after nuclear export factor 1 (NXF1) binding, challenging the current view of transcription-coupled mRNA export, and stress granule proteins prevalent in aged mRNPs, indicating roles in late stages of the mRNA life cycle. To systematically identify mRBPs with unknown functions, we employed machine learning to compare mRNA binding dynamics with Gene Ontology (GO) annotations. Our data can be explored at chronology.rna.snu.ac.kr.

N4-Allylcytidine: a new nucleoside analogue for RNA labelling and chemical sequencing

RNA labelling has become indispensable in studying RNA biology. Nucleoside analogues with a chemical sequencing power represent desirable RNA labelling molecules because precise labelling information at base resolution can be obtained. Here, we report a new nucleoside analogue, N4-allylcytidine (a4C), which is able to tag RNA through both in vitro and in vivo pathways and further specifically reacts with iodine to form 3, N4-cyclized cytidine (cyc-C) in a catalyst-free, fast and complete manner. Full spectroscopic characterization concluded that cyc-C consisted of paired diastereoisomers with opposite chiral carbon centers in the fused 3, N4-five-membered ring. During RNA reverse transcription into complementary DNA, cyc-C induces base misincorporation due to the disruption of canonical hydrogen bonding by the cyclized structure and thus can be accurately identified by sequencing at single base resolution. With the chemical sequencing rationale of a4C, successful applications have been performed including pinpointing N4-methylcytidine methyltransferases' substrate modification sites, metabolically labelling mammalian cellular RNAs, and mapping active cellular RNA polymerase locations with the chromatin run-on RNA sequencing technique. Collectively, our work demonstrates that a4C is a promising molecule for RNA labelling and chemical sequencing and expands the toolkit for studying sophisticated RNA biology.

Screening radiation-differentially expressed circular RNAs and establishing dose classification models in the human lymphoblastoid cell line AHH-1

PURPOSE

In the event of a large-scale radiological accident, rapid and high-throughput biodosimetry is the most vital basis in medical resource allocation for the prompt treatment of victims. However, the current biodosimeter is yet to be rapid and high-throughput. Studies have shown that ionizing radiation modulates expressions of circular RNAs (circRNAs) in healthy human cell lines and tumor tissue. circRNA expressions can be quantified rapidly and high-throughput. However, whether circRNAs are suitable for early radiation dose classification remains unclear.

METHODS

We employed transcriptome sequencing and bioinformatics analysis to screen for radiation-differentially expressed circRNAs in the human lymphoblastoid cell line AHH-1 at 4 h following exposure to 0, 2, and 5 Gy 60Co γ-rays. The dose-response relationships between differentially expressed circRNA expressions and absorbed doses were investigated using real-time polymerase chain reaction and linear regression analysis at 4 h, 24 h, and 48 h post-exposure to 0, 2, 4, 6, and 8 Gy. Six distinct dose classification models of circRNA panels were established and validated by receiver operating characteristic (ROC) curve analysis.

RESULTS

A total of 11 radiation-differentially expressed circRNAs were identified and validated. Based on dose-response effects, those circRNAs changed in a dose-responsive or dose-dependent manner were combined into panels A through F at 4 h, 24 h, and 48 h post-irradiation. ROC curve analysis showed that panels A through C had the potential to effectively classify exposed and non-exposed conditions, which area under the curve (AUC) of these three panels were all 1.000, and the associate p values were .009. Panels D through F excellently distinguished between different dose groups (AUC = 0.963-1.000, p < .05). The validation assay showed that panels A through F demonstrated consistent excellence in sensitivity and specificity in dose classification.

CONCLUSIONS

Ionizing radiation can indeed modulate the circRNA expression profile in the human lymphoblastoid cell line AHH-1. The differentially expressed circRNAs exhibit the potential for rapid and high-throughput dose classification.

Late consolidation of rRNA structure during co-transcriptional assembly in E. coli by time-resolved DMS footprinting

The production of new ribosomes requires proper folding of the rRNA and the addition of more than 50 ribosomal proteins. The structures of some assembly intermediates have been determined by cryo-electron microscopy, yet these structures do not provide information on the folding dynamics of the rRNA. To visualize the changes in rRNA structure during ribosome assembly in E. coli cells, transcripts were pulse-labeled with 4-thiouridine and the structure of newly made rRNA probed at various times by dimethyl sulfate modification and mutational profiling sequencing (4U-DMS-MaPseq). The in-cell DMS modification patterns revealed that many long-range rRNA tertiary interactions and protein binding sites through the 16S and 23S rRNA remain partially unfolded 1.5 min after transcription. By contrast, the active sites were continually shielded from DMS modification, suggesting that these critical regions are guarded by cellular factors throughout assembly. Later, bases near the peptidyl tRNA site exhibited specific rearrangements consistent with the binding and release of assembly factors. Time-dependent structure-probing in cells suggests that many tertiary interactions throughout the new ribosomal subunits remain mobile or unfolded until the late stages of subunit maturation.

Exploring IDI2-AS1, OIP5-AS1, and LITATS1: Changes in Long Non-coding RNAs Induced by the Poly I:C Stimulation

Long non-coding RNAs (lncRNAs) are sequences longer than 200 nucleotides, but they do not encode proteins. Nevertheless, they have significant roles in diverse biological functions. It remains unclear how viral infections trigger the expression of lncRNAs. In our previous research, we revealed a distinct type of lncRNAs with a lifespan under 4 h in human HeLa cells. These short-lived lncRNAs might be associated with numerous regulatory roles. Given their potential impact on human physiology, these short-lived lncRNAs could be key indicators to measure polyinosinic:polycytidylic acid (poly I:C) stimulation. In our recent work, we discovered three lncRNAs: IDI2-AS1, OIP5-AS1, and LITATS1. After exposure to poly I:C, imitating viral assault in human A549 cells, IDI2-AS1 levels dropped significantly while OIP5-AS1 and LITATS1 levels rose markedly. Our results indicate that short-lived lncRNAs respond to poly I:C stimulation, exhibiting substantial changes in expression. This indicates that the understanding the role of lncRNAs in the host response to viral infection and the potential for these molecules to serve as novel therapeutic targets.

Deep learning and direct sequencing of labeled RNA captures transcriptome dynamics

Quantification of the dynamics of RNA metabolism is essential for understanding gene regulation in health and disease. Existing methods rely on metabolic labeling of nascent RNAs and physical separation or inference of labeling through PCR-generated mutations, followed by short-read sequencing. However, these methods are limited in their ability to identify transient decay intermediates or co-analyze RNA decay with cis-regulatory elements of RNA stability such as poly(A) tail length and modification status, at single molecule resolution. Here we use 5-ethynyl uridine (5EU) to label nascent RNA followed by direct RNA sequencing with nanopores. We developed RNAkinet, a deep convolutional and recurrent neural network that processes the electrical signal produced by nanopore sequencing to identify 5EU-labeled nascent RNA molecules. RNAkinet demonstrates generalizability to distinct cell types and organisms and reproducibly quantifies RNA kinetic parameters allowing the combined interrogation of RNA metabolism and cis-acting RNA regulatory elements.

Profiling dynamic RNA-protein interactions using small-molecule-induced RNA editing

RNA-binding proteins (RBPs) play an important role in biology, and characterizing dynamic RNA-protein interactions is essential for understanding RBP function. In this study, we developed targets of RBPs identified by editing induced through dimerization (TRIBE-ID), a facile strategy for quantifying state-specific RNA-protein interactions upon rapamycin-mediated chemically induced dimerization and RNA editing. We performed TRIBE-ID with G3BP1 and YBX1 to study RNA-protein interactions during normal conditions and upon oxidative stress-induced biomolecular condensate formation. We quantified editing kinetics to infer interaction persistence and show that stress granule formation strengthens pre-existing RNA-protein interactions and induces new RNA-protein binding events. Furthermore, we demonstrate that G3BP1 stabilizes its targets under normal and oxidative stress conditions independent of stress granule formation. Finally, we apply our method to characterize small-molecule modulators of G3BP1-RNA binding. Taken together, our work provides a general approach to profile dynamic RNA-protein interactions in cellular contexts with temporal control.

a6A-seq: N 6-allyladenosine-based cellular messenger RNA metabolic labelling and sequencing

The integration of RNA metabolic labelling by nucleoside analogues with high-throughput RNA sequencing has been harnessed to study RNA dynamics. The immunoprecipitation purification or chemical pulldown technique is generally required to enrich the analogue-labelled RNAs. Here we developed an a6A-seq method, which takes advantage of N6-allyladenosine (a6A) metabolic labelling on cellular mRNAs and profiles them in an immunoprecipitation-free and mutation-based manner. a6A plays a role as a chemical sequencing tag in that the iodination of a6A in mRNAs results in 1,N 6-cyclized adenosine (cyc-A), which induces base misincorporation during RNA reverse transcription, thus making a6A-labelled mRNAs detectable by sequencing. A nucleic acid melting assay was utilized to investigate why cyc-A prefers to be paired with guanine. a6A-seq was utilized to study cellular gene expression changes under a methionine-free stress condition. Compared with regular RNA-seq, a6A-seq could more sensitively detect the change of mRNA production over a time scale. The experiment of a6A-containing mRNA immunoprecipitation followed by qPCR successfully validated the high-throughput a6A-seq data. Together, our results show a6A-seq is an effective tool to study RNA dynamics.

Simultaneous measurement of nascent transcriptome and translatome using 4-thiouridine metabolic RNA labeling and translating ribosome affinity purification

Regulation of gene expression in response to various biological processes, including extracellular stimulation and environmental adaptation requires nascent RNA synthesis and translation. Analysis of the coordinated regulation of dynamic RNA synthesis and translation is required to determine functional protein production. However, reliable methods for the simultaneous measurement of nascent RNA synthesis and translation at the gene level are limited. Here, we developed a novel method for the simultaneous assessment of nascent RNA synthesis and translation by combining 4-thiouridine (4sU) metabolic RNA labeling and translating ribosome affinity purification (TRAP) using a monoclonal antibody against evolutionarily conserved ribosomal P-stalk proteins. The P-stalk-mediated TRAP (P-TRAP) technique recovered endogenous translating ribosomes, allowing easy translatome analysis of various eukaryotes. We validated this method in mammalian cells by demonstrating that acute unfolded protein response (UPR) in the endoplasmic reticulum (ER) induces dynamic reprogramming of nascent RNA synthesis and translation. Our nascent P-TRAP (nP-TRAP) method may serve as a simple and powerful tool for analyzing the coordinated regulation of transcription and translation of individual genes in various eukaryotes.

Nuclear SRP9/SRP14 heterodimer transcriptionally regulates 7SL and BC200 RNA expression

The SRP9/SRP14 heterodimer is a central component of signal recognition particle (SRP) RNA (7SL) processing and Alu retrotransposition. In this study, we sought to establish the role of nuclear SRP9/SRP14 in the transcriptional regulation of 7SL and BC200 RNA. 7SL and BC200 RNA steady-state levels, rate of decay, and transcriptional activity were evaluated under SRP9/SRP14 knockdown conditions. Immunofluorescent imaging, and subcellular fractionation of MCF-7 cells, revealed a distinct nuclear localization for SRP9/SRP14. The relationship between this localization and transcriptional activity at 7SL and BC200 genes was also examined. These findings demonstrate a novel nuclear function of SRP9/SRP14 establishing that this heterodimer transcriptionally regulates 7SL and BC200 RNA expression. We describe a model in which SRP9/SRP14 cotranscriptionally regulate 7SL and BC200 RNA expression. Our model is also a plausible pathway for regulating Alu RNA transcription and is consistent with the hypothesized roles of SRP9/SRP14 transporting 7SL RNA into the nucleolus for posttranscriptional processing, and trafficking of Alu RNA for retrotransposition.

Translation reinitiation after uORFs does not fully protect mRNAs from nonsense-mediated decay

It is estimated that nearly 50% of mammalian transcripts contain at least one upstream open reading frame (uORF), which are typically one to two orders of magnitude smaller than the downstream main ORF. Most uORFs are thought to be inhibitory as they sequester the scanning ribosome, but in some cases allow for translation reinitiation. However, termination in the 5' UTR at the end of uORFs resembles premature termination that is normally sensed by the nonsense-mediated mRNA decay (NMD) pathway. Translation reinitiation has been proposed as a method for mRNAs to prevent NMD. Here, we test how uORF length influences translation reinitiation and mRNA stability in HeLa cells. Using custom 5' UTRs and uORF sequences, we show that reinitiation can occur on heterologous mRNA sequences, favors small uORFs, and is supported when initiation occurs with more initiation factors. After determining reporter mRNA half-lives in HeLa cells and mining available mRNA half-life data sets for cumulative predicted uORF length, we conclude that translation reinitiation after uORFs is not a robust method for mRNAs to prevent NMD. Together, these data suggest that the decision of whether NMD ensues after translating uORFs occurs before reinitiation in mammalian cells.

Regulation of mature mRNA levels by RNA processing efficiency

Transcription and co-transcriptional processes, including pre-mRNA splicing and mRNA cleavage and polyadenylation, regulate the production of mature mRNAs. The carboxyl terminal domain (CTD) of RNA polymerase (pol) II, which comprises 52 repeats of the Tyr1Ser2Pro3Thr4Ser5Pro6Ser7 peptide, is involved in the coordination of transcription with co-transcriptional processes. The pol II CTD is dynamically modified by protein phosphorylation, which regulates recruitment of transcription and co-transcriptional factors. We have investigated whether mature mRNA levels from intron-containing protein-coding genes are related to pol II CTD phosphorylation, RNA stability, and pre-mRNA splicing and mRNA cleavage and polyadenylation efficiency. We find that genes that produce a low level of mature mRNAs are associated with relatively high phosphorylation of the pol II CTD Thr4 residue, poor RNA processing, increased chromatin association of transcripts, and shorter RNA half-life. While these poorly-processed transcripts are degraded by the nuclear RNA exosome, our results indicate that in addition to RNA half-life, chromatin association due to a low RNA processing efficiency also plays an important role in the regulation of mature mRNA levels.

RNA degradome analysis reveals DNE1 endoribonuclease is required for the turnover of diverse mRNA substrates in Arabidopsis

In plants, cytoplasmic mRNA decay is critical for posttranscriptionally controlling gene expression and for maintaining cellular RNA homeostasis. Arabidopsis DCP1-ASSOCIATED NYN ENDORIBONUCLEASE 1 (DNE1) is a cytoplasmic mRNA decay factor that interacts with proteins involved in mRNA decapping and nonsense-mediated mRNA decay (NMD). There is limited information on the functional role of DNE1 in RNA turnover, and the identities of its endogenous targets are unknown. In this study, we utilized RNA degradome approaches to globally investigate DNE1 substrates. Monophosphorylated 5' ends, produced by DNE1, should accumulate in mutants lacking the cytoplasmic exoribonuclease XRN4, but be absent from DNE1 and XRN4 double mutants. In seedlings, we identified over 200 such transcripts, most of which reflect cleavage within coding regions. While most DNE1 targets were NMD-insensitive, some were upstream ORF (uORF)-containing and NMD-sensitive transcripts, indicating that this endoribonuclease is required for turnover of a diverse set of mRNAs. Transgenic plants expressing DNE1 cDNA with an active-site mutation in the endoribonuclease domain abolished the in planta cleavage of transcripts, demonstrating that DNE1 endoribonuclease activity is required for cleavage. Our work provides key insights into the identity of DNE1 substrates and enhances our understanding of DNE1-mediated mRNA decay.

Preanalytic and Analytic Quality System Considerations in Noncoding RNA Biomarker Development for Clinical Diagnostics

A frequent topic of biomedical research is the potential clinical use of non-coding (nc) RNAs as quantitative biomarkers for a broad spectrum of health and disease. However, ncRNA analyses have not been pressed into widespread diagnostic use. Strong preclinical evidence suggests obstacles in the translation and reproducibility of this type of biomarker which may result from preanalytical and analytical variation in the non-standardized processes used to collect, process, and store samples, as well as the substantive differences between small and long ncRNA. We performed a narrative review of selected literature, through the lens of key laboratory-developed test (LDT) regulations under the Clinical Laboratory Improvement Amendments (CLIA) in the United States, to study critical gaps in ncRNA validation studies. This review describes the leading candidate ncRNA subclasses, their biogenesis and cellular function, and identifies specific pre-analytical variables with disproportionate impact on testing performance. We summarize these findings with strategic recommendations to clinicians and biomedical scientists involved in the design, conduct, and translation of ncRNA biomarker development.

Cell-to-cell variability in Myc dynamics drives transcriptional heterogeneity in cancer cells

Cell-to-cell heterogeneity is vital for tumor evolution and survival. How cancer cells achieve and exploit this heterogeneity remains an active area of research. Here, we identify c-Myc as a highly heterogeneously expressed transcription factor and an orchestrator of transcriptional and phenotypic diversity in cancer cells. By monitoring endogenous c-Myc protein in individual living cells, we report the surprising pulsatile nature of c-Myc expression and the extensive cell-to-cell variability in its dynamics. We further show that heterogeneity in c-Myc dynamics leads to variable target gene transcription and that timing of c-Myc expression predicts cell-cycle progression rates and drug sensitivities. Together, our data advocate for a model in which cancer cells increase the heterogeneity of functionally diverse transcription factors such as c-Myc to rapidly survey transcriptional landscapes and survive stress.

Differential nuclear import regulates nuclear RNA inheritance following mitosis

Mitosis results in a dramatic reorganization of chromatin structure to promote chromosome compaction and segregation to daughter cells. Consequently, mitotic entry is accompanied by transcriptional silencing and removal of most chromatin-bound RNA from chromosomes. As cells exit mitosis, chromatin rapidly decondenses and transcription restarts as waves of differential gene expression. However, little is known about the fate of chromatin-bound RNAs following cell division. Here we explored whether nuclear RNA from the previous cell cycle is present in G1 nuclei following mitosis. We found that half of all nuclear RNA is inherited in a transcription-independent manner following mitosis. Interestingly, the snRNA U2 is efficiently inherited by G1 nuclei, while the lncRNAs NEAT1 and MALAT1 show no inheritance following mitosis. We found that the nuclear protein SAF-A, which is hypothesized to tether RNA to DNA, did not play a prominent role in nuclear RNA inheritance, indicating that the mechanism for RNA inheritance may not involve RNA chaperones that have chromatin-binding activity. Instead, we observe that the timing of RNA inheritance indicates that a select group of nuclear RNAs are reimported into the nucleus after the nuclear envelope has reassembled. Our work demonstrates that there is a fraction of nuclear RNA from the previous cell cycle that is reimported following mitosis and suggests that mitosis may serve as a time to reset the interaction of lncRNAs with chromatin.

Higher-order modular regulation of the human proteome

Operons are transcriptional modules that allow bacteria to adapt to environmental changes by coordinately expressing the relevant set of genes. In humans, biological pathways and their regulation are more complex. If and how human cells coordinate the expression of entire biological processes is unclear. Here, we capture 31 higher-order co-regulation modules, which we term progulons, by help of supervised machine-learning on proteomics data. Progulons consist of dozens to hundreds of proteins that together mediate core cellular functions. They are not restricted to physical interactions or co-localisation. Progulon abundance changes are primarily controlled at the level of protein synthesis and degradation. Implemented as a web app at www.proteomehd.net/progulonFinder, our approach enables the targeted search for progulons of specific cellular processes. We use it to identify a DNA replication progulon and reveal multiple new replication factors, validated by extensive phenotyping of siRNA-induced knockdowns. Progulons provide a new entry point into the molecular understanding of biological processes.

Adolescent binge ethanol impacts H3K36me3 regulation of synaptic genes

Adolescence is marked in part by the ongoing development of the prefrontal cortex (PFC). Binge ethanol use during this critical stage in neurodevelopment induces significant structural changes to the PFC, as well as cognitive and behavioral deficits that can last into adulthood. Previous studies showed that adolescent binge ethanol causes lasting deficits in working memory, decreases in the expression of chromatin remodeling genes responsible for the methylation of histone 3 lysine 36 (H3K36), and global decreases in H3K36 in the PFC. H3K36me3 is present within the coding region of actively-transcribed genes, and safeguards against aberrant, cryptic transcription by RNA Polymerase II. We hypothesize that altered methylation of H3K36 could play a role in adolescent binge ethanol-induced memory deficits. To investigate this at the molecular level, ethanol (4 g/kg, i.g.) or water was administered intermittently to adolescent mice. RNA-and ChIP-sequencing were then performed within the same tissue to determine gene expression changes and identify genes and loci where H3K36me3 was disrupted by ethanol. We further assessed ethanol-induced changes at the transcription level with differential exon-use and cryptic transcription analysis - a hallmark of decreased H3K36me3. Here, we found ethanol-induced changes to the gene expression and H3K36me3-regulation of synaptic-related genes in all our analyses. Notably, H3K36me3 was differentially trimethylated between ethanol and control conditions at synaptic-related genes, and Snap25 and Cplx1 showed evidence of cryptic transcription in males and females treated with ethanol during adolescence. Our results provide preliminary evidence that ethanol-induced changes to H3K36me3 during adolescent neurodevelopment may be linked to synaptic dysregulation at the transcriptional level, which may explain the reported ethanol-induced changes to PFC synaptic function.

RNA out of the mist

RNA has long been regarded primarily as the intermediate between genes and proteins. It was a surprise then to discover that eukaryotic genes are mosaics of mRNA sequences interrupted by large tracts of transcribed but untranslated sequences, and that multicellular organisms also express many long 'intergenic' and antisense noncoding RNAs (lncRNAs). The identification of small RNAs that regulate mRNA translation and half-life did not disturb the prevailing view that animals and plant genomes are full of evolutionary debris and that their development is mainly supervised by transcription factors. Gathering evidence to the contrary involved addressing the low conservation, expression, and genetic visibility of lncRNAs, demonstrating their cell-specific roles in cell and developmental biology, and their association with chromatin-modifying complexes and phase-separated domains. The emerging picture is that most lncRNAs are the products of genetic loci termed 'enhancers', which marshal generic effector proteins to their sites of action to control cell fate decisions during development.

Using protein-per-mRNA differences among human tissues in codon optimization

BACKGROUND

Codon usage and nucleotide composition of coding sequences have profound effects on protein expression. However, while it is recognized that different tissues have distinct tRNA profiles and codon usages in their transcriptomes, the effect of tissue-specific codon optimality on protein synthesis remains elusive.

RESULTS

We leverage existing state-of-the-art transcriptomics and proteomics datasets from the GTEx project and the Human Protein Atlas to compute the protein-to-mRNA ratios of 36 human tissues. Using this as a proxy of translational efficiency, we build a machine learning model that identifies codons enriched or depleted in specific tissues. We detect two clusters of tissues with an opposite pattern of codon preferences. We then use these identified patterns for the development of CUSTOM, a codon optimizer algorithm which suggests a synonymous codon design in order to optimize protein production in a tissue-specific manner. In human cell-line models, we provide evidence that codon optimization should take into account particularities of the translational machinery of the tissues in which the target proteins are expressed and that our approach can design genes with tissue-optimized expression profiles.

CONCLUSIONS

We provide proof-of-concept evidence that codon preferences exist in tissue-specific protein synthesis and demonstrate its application to synthetic gene design. We show that CUSTOM can be of benefit in biological and biotechnological applications, such as in the design of tissue-targeted therapies and vaccines.

Distinct dosage compensations of ploidy-sensitive and -insensitive X chromosome genes during development and in diseases

The active X chromosome in mammals is upregulated to balance its dosage to autosomes during evolution. However, it is elusive why the known dosage compensation machinery showed uneven and small influence on X genes. Here, based on >20,000 transcriptomes, we identified two X gene groups (ploidy-sensitive [PSX] and ploidy-insensitive [PIX]), showing distinct but evolutionarily conserved dosage compensations (termed XAR). We demonstrated that XAR-PIX was downregulated whereas XAR-PSX upregulated at both RNA and protein levels across cancer types, in contrast with their trends during stem cell differentiation. XAR-PIX, but not XAR-PSX, was lower and correlated with autoantibodies and inflammation in patients of lupus, suggesting that insufficient dosage of PIX genes contribute to lupus pathogenesis. We further identified and experimentally validated two XAR regulators, TP53 and ATRX. Collectively, we provided insights into X dosage compensation in mammals and demonstrated different regulation of PSX and PIX and their pathophysiological roles in human diseases.

Rapid nuclear deadenylation of mammalian messenger RNA

Poly(A) tails protect RNAs from degradation and their deadenylation rates determine RNA stability. Although poly(A) tails are generated in the nucleus, deadenylation of tails has mostly been investigated within the cytoplasm. Here, we combined long-read sequencing with metabolic labeling, splicing inhibition and cell fractionation experiments to quantify, separately, the genesis and trimming of nuclear and cytoplasmic tails in vitro and in vivo. We present evidence for genome-wide, nuclear synthesis of tails longer than 200 nt, which are rapidly shortened after transcription. Our data suggests that rapid deadenylation is a nuclear process, and that different classes of transcripts and even transcript isoforms have distinct nuclear tail lengths. For example, many long-noncoding RNAs retain long poly(A) tails. Modeling deadenylation dynamics predicts nuclear deadenylation about 10 times faster than cytoplasmic deadenylation. In summary, our data suggests that nuclear deadenylation might be a key mechanism for regulating mRNA stability, abundance, and subcellular localization.

The Emerging Role of uORF-Encoded uPeptides and HLA uLigands in Cellular and Tumor Biology

Recent technological advances have facilitated the detection of numerous non-canonical human peptides derived from regulatory regions of mRNAs, long non-coding RNAs, and other cryptic transcripts. In this review, we first give an overview of the classification of these novel peptides and summarize recent improvements in their annotation and detection by ribosome profiling, mass spectrometry, and individual experimental analysis. A large fraction of the novel peptides originates from translation at upstream open reading frames (uORFs) that are located within the transcript leader sequence of regular mRNA. In humans, uORF-encoded peptides (uPeptides) have been detected in both healthy and malignantly transformed cells and emerge as important regulators in cellular and immunological pathways. In the second part of the review, we focus on various functional implications of uPeptides. As uPeptides frequently act at the transition of translational regulation and individual peptide function, we describe the mechanistic modes of translational regulation through ribosome stalling, the involvement in cellular programs through protein interaction and complex formation, and their role within the human leukocyte antigen (HLA)-associated immunopeptidome as HLA uLigands. We delineate how malignant transformation may lead to the formation of novel uORFs, uPeptides, or HLA uLigands and explain their potential implication in tumor biology. Ultimately, we speculate on a potential use of uPeptides as peptide drugs and discuss how uPeptides and HLA uLigands may facilitate translational inhibition of oncogenic protein messages and immunotherapeutic approaches in cancer therapy.

The genetic and biochemical determinants of mRNA degradation rates in mammals

BACKGROUND

Degradation rate is a fundamental aspect of mRNA metabolism, and the factors governing it remain poorly characterized. Understanding the genetic and biochemical determinants of mRNA half-life would enable more precise identification of variants that perturb gene expression through post-transcriptional gene regulatory mechanisms.

RESULTS

We establish a compendium of 39 human and 27 mouse transcriptome-wide mRNA decay rate datasets. A meta-analysis of these data identified a prevalence of technical noise and measurement bias, induced partially by the underlying experimental strategy. Correcting for these biases allowed us to derive more precise, consensus measurements of half-life which exhibit enhanced consistency between species. We trained substantially improved statistical models based upon genetic and biochemical features to better predict half-life and characterize the factors molding it. Our state-of-the-art model, Saluki, is a hybrid convolutional and recurrent deep neural network which relies only upon an mRNA sequence annotated with coding frame and splice sites to predict half-life (r=0.77). The key novel principle learned by Saluki is that the spatial positioning of splice sites, codons, and RNA-binding motifs within an mRNA is strongly associated with mRNA half-life. Saluki predicts the impact of RNA sequences and genetic mutations therein on mRNA stability, in agreement with functional measurements derived from massively parallel reporter assays.

CONCLUSIONS

Our work produces a more robust ground truth for transcriptome-wide mRNA half-lives in mammalian cells. Using these revised measurements, we trained Saluki, a model that is over 50% more accurate in predicting half-life from sequence than existing models. Saluki succinctly captures many of the known determinants of mRNA half-life and can be rapidly deployed to predict the functional consequences of arbitrary mutations in the transcriptome.

Application of nucleoside or nucleotide analogues in RNA dynamics and RNA-binding protein analysis

Cellular RNAs undergo dynamic changes during RNA biological processes, which are tightly orchestrated by RNA-binding proteins (RBPs). Yet, the investigation of RNA dynamics is hurdled by highly abundant steady-state RNAs, which make the signals of dynamic RNAs less detectable. Notably, the exert of nucleoside or nucleotide analogue-based RNA technologies has provided a remarkable platform for RNA dynamics research, revealing diverse unnoticed features in RNA metabolism. In this review, we focus on the application of two types of analogue-based RNA sequencing, antigen-/antibody- and click chemistry-based methodologies, and summarize the RNA dynamics features revealed. Moreover, we discuss emerging single-cell newly transcribed RNA sequencing methodologies based on nucleoside analogue labeling, which provides novel insights into RNA dynamics regulation at single-cell resolution. On the other hand, we also emphasize the identification of RBPs that interact with polyA, non-polyA RNAs, or newly transcribed RNAs and also their associated RNA-binding domains at genomewide level through ultraviolet crosslinking and mass spectrometry in different contexts. We anticipated that further modification and development of these analogue-based RNA and RBP capture technologies will aid in obtaining an unprecedented understanding of RNA biology. This article is categorized under: RNA Interactions with Proteins and Other Molecules &gt; Protein-RNA Recognition RNA Structure and Dynamics &gt; RNA Structure, Dynamics and Chemistry RNA Methods &gt; RNA Analyses in Cells.

Long non-coding RNA and RNA-binding protein interactions in cancer: Experimental and machine learning approaches

Understanding the complex and specific roles played by non-coding RNAs (ncRNAs), which comprise the bulk of the genome, is important for understanding virtually every hallmark of cancer. This large group of molecules plays pivotal roles in key regulatory mechanisms in various cellular processes. Regulatory mechanisms, mediated by long non-coding RNA (lncRNA) and RNA-binding protein (RBP) interactions, are well documented in several types of cancer. Their effects are enabled through networks affecting lncRNA and RBP stability, RNA metabolism including N⁶-methyladenosine (m⁶A) and alternative splicing, subcellular localization, and numerous other mechanisms involved in cancer. In this review, we discuss the reciprocal interplay between lncRNAs and RBPs and their involvement in epigenetic regulation via histone modifications, as well as their key role in resistance to cancer therapy. Other aspects of RBPs including their structural domains, provide a deeper knowledge on how lncRNAs and RBPs interact and exert their biological functions. In addition, current state-of-the-art knowledge, facilitated by machine and deep learning approaches, unravels such interactions in better details to further enhance our understanding of the field, and the potential to harness RNA-based therapeutics as an alternative treatment modality for cancer are discussed.

Sequence features around cleavage sites are highly conserved among different species and a critical determinant for RNA cleavage position across eukaryotes

RNA degradation is one of the critical steps for control of gene expression, and endonucleolytic cleavage-dependent RNA degradation is conserved among eukaryotes. Some cleavage sites are secondarily capped in the cytoplasm and identified using the Cap analysis of gene expression (CAGE) method. Although uncapped cleavage sites are widespread in eukaryotes, comparatively little information has been obtained about these sites using CAGE-based degradome analysis. Previously, we developed the truncated RNA-end sequencing (TREseq) method in plant species and used it to acquire comprehensive information about uncapped cleavage sites; we observed G-rich sequences near cleavage sites. However, it remains unclear whether this finding is general to other eukaryotes. In this study, we conducted TREseq analyses in fruit flies (Drosophila melanogaster) and budding yeast (Saccharomyces cerevisiae). The results revealed specific sequence features related to RNA cleavage in D. melanogaster and S. cerevisiae that were similar to sequence patterns in Arabidopsis thaliana. Although previous studies suggest that ribosome movements are important for determining cleavage position, feature selection using a random forest classifier showed that sequences around cleavage sites were major determinant for cleaved or uncleaved sites. Together, our results suggest that sequence features around cleavage sites are critical for determining cleavage position, and that sequence-specific endonucleolytic cleavage-dependent RNA degradation is highly conserved across eukaryotes.

Control of immediate early gene expression by CPEB4-repressor complex-mediated mRNA degradation

Background

Cytoplasmic polyadenylation element-binding protein 4 (CPEB4) is known to associate with cytoplasmic polyadenylation elements (CPEs) located in the 3′ untranslated region (UTR) of specific mRNAs and assemble an activator complex promoting the translation of target mRNAs through cytoplasmic polyadenylation.

Results

Here, we find that CPEB4 is part of an alternative repressor complex that mediates mRNA degradation by associating with the evolutionarily conserved CCR4-NOT deadenylase complex. We identify human CPEB4 as an RNA-binding protein (RBP) with enhanced association to poly(A) RNA upon inhibition of class I histone deacetylases (HDACs), a condition known to cause widespread degradation of poly(A)-containing mRNA. Photoactivatable ribonucleoside-enhanced crosslinking and immunoprecipitation (PAR-CLIP) analysis using endogenously tagged CPEB4 in HeLa cells reveals that CPEB4 preferentially binds to the 3′UTR of immediate early gene mRNAs, at G-containing variants of the canonical U- and A-rich CPE located in close proximity to poly(A) sites. By transcriptome-wide mRNA decay measurements, we find that the strength of CPEB4 binding correlates with short mRNA half-lives and that loss of CPEB4 expression leads to the stabilization of immediate early gene mRNAs. Akin to CPEB4, we demonstrate that CPEB1 and CPEB2 also confer mRNA instability by recruitment of the CCR4-NOT complex.

Conclusions

While CPEB4 was previously known for its ability to stimulate cytoplasmic polyadenylation, our findings establish an additional function for CPEB4 as the RNA adaptor of a repressor complex that enhances the degradation of short-lived immediate early gene mRNAs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02760-5.

Live-Cell RNA Imaging with Metabolically Incorporated Fluorescent Nucleosides

Fluorescence imaging is a powerful method for probing macromolecular dynamics in biological systems; however, approaches for cellular RNA imaging are limited to the investigation of individual RNA constructs or bulk RNA labeling methods compatible primarily with fixed samples. Here, we develop a platform for fluorescence imaging of bulk RNA dynamics in living cells. We show that fluorescent bicyclic and tricyclic cytidine analogues can be metabolically incorporated into cellular RNA by overexpression of uridine-cytidine kinase 2. In particular, metabolic feeding with the tricyclic cytidine-derived nucleoside tC combined with confocal imaging enables the investigation of RNA synthesis, degradation, and trafficking at single-cell resolution. We apply our imaging modality to study RNA metabolism and localization during the oxidative stress response and find that bulk RNA turnover is greatly accelerated upon NaAsO2 treatment. Furthermore, we identify cytoplasmic RNA granules containing RNA transcripts generated during oxidative stress that are distinct from canonical stress granules and P-bodies and co-localize with the RNA helicase DDX6. Taken together, our work provides a powerful approach for live-cell RNA imaging and reveals how cells reshape RNA transcriptome dynamics in response to oxidative stress.

Comprehensive analysis of the circadian nuclear and cytoplasmic transcriptome in mouse liver

In eukaryotes, RNA is synthesised in the nucleus, spliced, and exported to the cytoplasm where it is translated and finally degraded. Any of these steps could be subject to temporal regulation during the circadian cycle, resulting in daily fluctuations of RNA accumulation and affecting the distribution of transcripts in different subcellular compartments. Our study analysed the nuclear and cytoplasmic, poly(A) and total transcriptomes of mouse livers collected over the course of a day. These data provide a genome-wide temporal inventory of enrichment in subcellular RNA, and revealed specific signatures of splicing, nuclear export and cytoplasmic mRNA stability related to transcript and gene lengths. Combined with a mathematical model describing rhythmic RNA profiles, we could test the rhythmicity of export rates and cytoplasmic degradation rates of approximately 1400 genes. With nuclear export times usually much shorter than cytoplasmic half-lives, we found that nuclear export contributes to the modulation and generation of rhythmic profiles of 10% of the cycling nuclear mRNAs. This study contributes to a better understanding of the dynamic regulation of the transcriptome during the day-night cycle.

Arabidopsis RNA processing body components LSM1 and DCP5 aid in the evasion of translational repression during Cauliflower mosaic virus infection

Viral infections impose extraordinary RNA stress, triggering cellular RNA surveillance pathways such as RNA decapping, nonsense-mediated decay, and RNA silencing. Viruses need to maneuver among these pathways to establish infection and succeed in producing high amounts of viral proteins. Processing bodies (PBs) are integral to RNA triage in eukaryotic cells, with several distinct RNA quality control pathways converging for selective RNA regulation. In this study, we investigated the role of Arabidopsis thaliana PBs during Cauliflower mosaic virus (CaMV) infection. We found that several PB components are co-opted into viral factories that support virus multiplication. This pro-viral role was not associated with RNA decay pathways but instead, we established that PB components are helpers in viral RNA translation. While CaMV is normally resilient to RNA silencing, dysfunctions in PB components expose the virus to this pathway, which is similar to previous observations for transgenes. Transgenes, however, undergo RNA quality control-dependent RNA degradation and transcriptional silencing, whereas CaMV RNA remains stable but becomes translationally repressed through decreased ribosome association, revealing a unique dependence among PBs, RNA silencing, and translational repression. Together, our study shows that PB components are co-opted by the virus to maintain efficient translation, a mechanism not associated with canonical PB functions.

Long-term imaging of individual mRNA molecules in living cells

Single-cell imaging of individual mRNAs has revealed core mechanisms of the central dogma. However, most approaches require cell fixation or have limited sensitivity for live-cell applications. Here, we describe SunRISER (SunTag-based reporter for imaging signal-enriched mRNA), a computationally and experimentally optimized approach for unambiguous detection of single mRNA molecules in living cells. When viewed by epifluorescence microscopy, SunRISER-labeled mRNAs show strong signal to background and resistance to photobleaching, which together enable long-term mRNA imaging studies. SunRISER variants, using 8× and 10× stem-loop arrays, demonstrate effective mRNA detection while significantly reducing alterations to target mRNA sequences. We characterize SunRISER to observe mRNA inheritance during mitosis and find that stressors enhance diversity among post-mitotic sister cells. Taken together, SunRISER enables a glimpse into living cells to observe aspects of the central dogma and the role of mRNAs in rare and dynamical trafficking events.

Global SLAM-seq for accurate mRNA decay determination and identification of NMD targets

Gene expression analysis requires accurate measurements of global RNA degradation rates, earlier problematic with methods disruptive to cell physiology. Recently, metabolic RNA labeling emerged as an efficient and minimally invasive technique applied in mammalian cells. Here, we have adapted SH-linked alkylation for the metabolic sequencing of RNA (SLAM-seq) for a global mRNA stability study in yeast using 4-thiouracil pulse-chase labeling. We assign high-confidence half-life estimates for 67.5% of expressed ORFs, and measure a median half-life of 9.4 min. For mRNAs where half-life estimates exist in the literature, their ranking order was in good agreement with previous data, indicating that SLAM-seq efficiently classifies stable and unstable transcripts. We then leveraged our yeast protocol to identify targets of the nonsense-mediated decay (NMD) pathway by measuring the change in RNA half-lives, instead of steady-state RNA level changes. With SLAM-seq, we assign 580 transcripts as putative NMD targets, based on their measured half-lives in wild-type and upf3Δ mutants. We find 225 novel targets, and observe a strong agreement with previous reports of NMD targets, 61.2% of our candidates being identified in previous studies. This indicates that SLAM-seq is a simpler and more economic method for global quantification of mRNA half-lives. Our adaptation for yeast yielded global quantitative measures of the NMD effect on transcript half-lives, high correlation with RNA half-lives measured previously with more technically challenging protocols, and identification of novel NMD regulated transcripts that escaped prior detection.

Mechanisms of cellular mRNA transcript homeostasis

For most genes, mRNA transcript abundance scales with cell size to ensure a constant concentration. Scaling of mRNA synthesis rates with cell size plays an important role, with regulation of the activity and abundance of RNA polymerase II (Pol II) now emerging as a key point of control. However, there is also considerable evidence for feedback mechanisms that kinetically couple the rates of mRNA synthesis, nuclear export, and degradation to allow cells to compensate for changes in one by adjusting the others. Researchers are beginning to integrate results from these different fields to reveal the mechanisms underlying transcript homeostasis. This will be crucial for moving beyond our current understanding of relative gene expression towards an appreciation of how absolute transcript levels are linked to other aspects of the cellular phenotype.

Feedback from nuclear RNA on transcription promotes robust RNA concentration homeostasis in human cells

RNA concentration homeostasis involves coordinating RNA abundance and synthesis rates with cell size. Here, we study this in human cells by combining genome-wide perturbations with quantitative single-cell measurements. Despite relative ease in perturbing RNA synthesis, we find that RNA concentrations generally remain highly constant. Perturbations that would be expected to increase nuclear mRNA levels, including those targeting nuclear mRNA degradation or export, result in downregulation of RNA synthesis. This is associated with reduced abundance of transcription-associated proteins and protein states that are normally coordinated with RNA production in single cells, including RNA polymerase II (RNA Pol II) itself. Acute perturbations, elevation of nuclear mRNA levels, and mathematical modeling indicate that mammalian cells achieve robust mRNA concentration homeostasis by the mRNA-based negative feedback on transcriptional activity in the nucleus. This ultimately acts to coordinate RNA Pol II abundance with nuclear mRNA degradation and export rates and may underpin the scaling of mRNA abundance with cell size.

Modulation of RNA stability regulates gene expression in two opposite ways: through buffering of RNA levels upon global perturbations and by supporting adapted differential expression

The steady state levels of RNAs, often referred to as expression levels, result from a well-balanced combination of RNA transcription and decay. Alterations in RNA levels will therefore result from tight regulation of transcription rates, decay rates or both. Here, we explore the role of RNA stability in achieving balanced gene expression and present genome-wide RNA stabilities in Drosophila melanogaster male and female cells as well as male cells depleted of proteins essential for dosage compensation. We identify two distinct RNA-stability mediated responses involved in regulation of gene expression. The first of these responds to acute and global changes in transcription and thus counteracts potentially harmful gene mis-expression by shifting the RNA stability in the direction opposite to the transcriptional change. The second response enhances inter-individual differential gene expression by adjusting the RNA stability in the same direction as a transcriptional change. Both mechanisms are global, act on housekeeping as well as non-housekeeping genes and were observed in both flies and mammals. Additionally, we show that, in contrast to mammals, modulation of RNA stability does not detectably contribute to dosage compensation of the sex-chromosomes in D. melanogaster.

Sorting and packaging of RNA into extracellular vesicles shape intracellular transcript levels

BACKGROUND

Extracellular vesicles (EVs) are released by nearly every cell type and have attracted much attention for their ability to transfer protein and diverse RNA species from donor to recipient cells. Much attention has been given so far to the features of EV short RNAs such as miRNAs. However, while the presence of mRNA and long noncoding RNA (lncRNA) transcripts in EVs has also been reported by multiple different groups, the properties and function of these longer transcripts have been less thoroughly explored than EV miRNA. Additionally, the impact of EV export on the transcriptome of exporting cells has remained almost completely unexamined. Here, we globally investigate mRNA and lncRNA transcripts in endothelial EVs in multiple different conditions.

RESULTS

In basal conditions, long RNA transcripts enriched in EVs have longer than average half-lives and distinctive stability-related sequence and structure characteristics including shorter transcript length, higher exon density, and fewer 3' UTR A/U-rich elements. EV-enriched long RNA transcripts are also enriched in HNRNPA2B1 binding motifs and are impacted by HNRNPA2B1 depletion, implicating this RNA-binding protein in the sorting of long RNA to EVs. After signaling-dependent modification of the cellular transcriptome, we observed that, unexpectedly, the rate of EV enrichment relative to cells was altered for many mRNA and lncRNA transcripts. This change in EV enrichment was negatively correlated with intracellular abundance, with transcripts whose export to EVs increased showing decreased abundance in cells and vice versa. Correspondingly, after treatment with inhibitors of EV secretion, levels of mRNA and lncRNA transcripts that are normally highly exported to EVs increased in cells, indicating a measurable impact of EV export on the long RNA transcriptome of the exporting cells. Compounds with different mechanisms of inhibition of EV secretion affected the cellular transcriptome differently, suggesting the existence of multiple EV subtypes with different long RNA profiles.

CONCLUSIONS

We present evidence for an impact of EV physiology on the characteristics of EV-producing cell transcriptomes. Our work suggests a new paradigm in which the sorting and packaging of transcripts into EVs participate, together with transcription and RNA decay, in controlling RNA homeostasis and shape the cellular long RNA abundance profile.

Genetic variants associated mRNA stability in lung

Background

Expression quantitative trait loci (eQTLs) analyses have been widely used to identify genetic variants associated with gene expression levels to understand what molecular mechanisms underlie genetic traits. The resultant eQTLs might affect the expression of associated genes through transcriptional or post-transcriptional regulation. In this study, we attempt to distinguish these two types of regulation by identifying genetic variants associated with mRNA stability of genes (stQTLs).

Results

Here, we presented a computational framework that takes advantage of recently developed methods to infer the mRNA stability of genes based on RNA-seq data and performed association analysis to identify stQTLs. Using the Genotype-Tissue Expression (GTEx) lung RNA-Seq data, we identified a total of 142,801 stQTLs for 3942 genes and 186,132 eQTLs for 4751 genes from 15,122,700 genetic variants for 13,476 genes on the autosomes, respectively. Interestingly, our results indicated that stQTLs were enriched in the CDS and 3’UTR regions, while eQTLs are enriched in the CDS, 3’UTR, 5’UTR, and upstream regions. We also found that stQTLs are more likely than eQTLs to overlap with RNA binding protein (RBP) and microRNA (miRNA) binding sites. Our analyses demonstrate that simultaneous identification of stQTLs and eQTLs can provide more mechanistic insight on the association between genetic variants and gene expression levels.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08405-y.

Data-driven identification of inherent features of eukaryotic stress-responsive genes

Living organisms are continuously challenged by changes in their environment that can propagate to stresses at the cellular level, such as rapid changes in osmolarity or oxygen tension. To survive these sudden changes, cells have developed stress-responsive mechanisms that tune cellular processes. The response of Saccharomyces cerevisiae to osmostress includes a massive reprogramming of gene expression. Identifying the inherent features of stress-responsive genes is of significant interest for understanding the basic principles underlying the rewiring of gene expression upon stress. Here, we generated a comprehensive catalog of osmostress-responsive genes from 5 independent RNA-seq experiments. We explored 30 features of yeast genes and found that 25 (83%) were distinct in osmostress-responsive genes. We then identified 13 non-redundant minimal osmostress gene traits and used statistical modeling to rank the most stress-predictive features. Intriguingly, the most relevant features of osmostress-responsive genes are the number of transcription factors targeting them and gene conservation. Using data on HeLa samples, we showed that the same features that define yeast osmostress-responsive genes can predict osmostress-responsive genes in humans, but with changes in the rank-ordering of feature-importance. Our study provides a holistic understanding of the basic principles of the regulation of stress-responsive gene expression across eukaryotes.

Molecular and cellular functions of long non-coding RNAs in prostate and breast cancer

Long noncoding RNAs (lncRNAs) are defined as noncoding RNA transcripts with a length greater than 200 nucleotides. Research over the last decade has made great strides in our understanding of lncRNAs, especially in the biology of their role in cancer. In this article, we will briefly discuss the biogenesis and characteristics of lncRNAs, then review their molecular and cellular functions in cancer by using prostate and breast cancer as examples. LncRNAs are abundant, diverse, and evolutionarily, less conserved than protein-coding genes. They are often expressed in a tumor and cell-specific manner. As a key epigenetic factor, lncRNAs can use a wide variety of molecular mechanisms to regulate gene expression at each step of the genetic information flow pathway. LncRNAs display widespread effects on cell behavior, tumor growth, and metastasis. They act intracellularly and extracellularly in an autocrine, paracrine and endocrine fashion. Increased understanding of lncRNA's role in cancer has facilitated the development of novel biomarkers for cancer diagnosis, led to greater understanding of cancer prognosis, enabled better prediction of therapeutic responses, and promoted identification of potential targets for cancer therapy.

A method for simultaneous detection of small and long RNA biotypes by ribodepleted RNA-Seq

RNA sequencing offers unprecedented access to the transcriptome. Key to this is the identification and quantification of many different species of RNA from the same sample at the same time. In this study we describe a novel protocol for simultaneous detection of coding and non-coding transcripts using modifications to the Ion Total RNA-Seq kit v2 protocol, with integration of QIASeq FastSelect rRNA removal kit. We report highly consistent sequencing libraries can be produced from both frozen high integrity mouse hippocampal tissue and the more challenging post-mortem human tissue. Removal of rRNA using FastSelect was extremely efficient, resulting in less than 1.5% rRNA content in the final library. We identified > 30,000 unique transcripts from all samples, including protein-coding genes and many species of non-coding RNA, in biologically-relevant proportions. Furthermore, the normalized sequencing read count for select genes significantly negatively correlated with Ct values from qRT-PCR analysis from the same samples. These results indicate that this protocol accurately and consistently identifies and quantifies a wide variety of transcripts simultaneously. The highly efficient rRNA depletion, coupled with minimized sample handling and without complicated and high-loss size selection protocols, makes this protocol useful to researchers wishing to investigate whole transcriptomes.

Massively parallel analysis of human 3' UTRs reveals that AU-rich element length and registration predict mRNA destabilization

AU-rich elements (AREs) are 3' UTR cis-regulatory elements that regulate the stability of mRNAs. Consensus ARE motifs have been determined, but little is known about how differences in 3' UTR sequences that conform to these motifs affect their function. Here, we use functional annotation of sequences from 3' UTRs (fast-UTR), a massively parallel reporter assay (MPRA), to investigate the effects of 41,288 3' UTR sequence fragments from 4653 transcripts on gene expression and mRNA stability in Jurkat and Beas2B cells. Our analyses demonstrate that the length of an ARE and its registration (the first and last nucleotides of the repeating ARE motif) have significant effects on gene expression and stability. Based on this finding, we propose improved ARE classification and concomitant methods to categorize and predict the effect of AREs on gene expression and stability. Finally, to investigate the advantages of our general experimental design we examine other motifs including constitutive decay elements (CDEs), where we show that the length of the CDE stem-loop has a significant impact on steady-state expression and mRNA stability. We conclude that fast-UTR, in conjunction with our analytical approach, can produce improved yet simple sequence-based rules for predicting the activity of human 3' UTRs.