A Detailed Protocol for Subcellular RNA Sequencing (subRNA-seq)

Control of endothelial cell function and arteriogenesis by MEG3:EZH2 epigenetic regulation of integrin expression

Epigenetic processes involving long non-coding RNAs regulate endothelial gene expression. However, the underlying regulatory mechanisms causing endothelial dysfunction remain to be elucidated. Enhancer of zeste homolog 2 (EZH2) is an important rheostat of histone H3K27 trimethylation (H3K27me3) that represses endothelial targets, but EZH2 RNA binding capacity and EZH2:RNA functional interactions have not been explored in post-ischemic angiogenesis. We used formaldehyde/UV-assisted crosslinking ligation and sequencing of hybrids and identified a new role for maternally expressed gene 3 (MEG3). MEG3 formed the predominant RNA:RNA hybrid structures in endothelial cells. Moreover, MEG3:EZH2 assists recruitment onto chromatin. By EZH2-chromatin immunoprecipitation, following MEG3 depletion, we demonstrated that MEG3 controls recruitment of EZH2/H3K27me3 onto integrin subunit alpha4 (ITGA4) promoter. Both MEG3 knockdown or EZH2 inhibition (A-395) promoted ITGA4 expression and improved endothelial cell migration and adhesion to fibronectin in vitro. The A-395 inhibitor re-directed MEG3-assisted chromatin remodeling, offering a direct therapeutic benefit by increasing endothelial function and resilience. This approach subsequently increased the expression of ITGA4 in arterioles following ischemic injury in mice, thus promoting arteriogenesis. Our findings show a context-specific role for MEG3 in guiding EZH2 to repress ITGA4. Novel therapeutic strategies could antagonize MEG3:EZH2 interaction for pre-clinical studies.

Post-transcriptional splicing can occur in a slow-moving zone around the gene

Splicing is the stepwise molecular process by which introns are removed from pre-mRNA and exons are joined together to form mature mRNA sequences. The ordering and spatial distribution of these steps remain controversial, with opposing models suggesting splicing occurs either during or after transcription. We used single-molecule RNA FISH, expansion microscopy, and live-cell imaging to reveal the spatiotemporal distribution of nascent transcripts in mammalian cells. At super-resolution levels, we found that pre-mRNA formed clouds around the transcription site. These clouds indicate the existence of a transcription-site-proximal zone through which RNA move more slowly than in the nucleoplasm. Full-length pre-mRNA undergo continuous splicing as they move through this zone following transcription, suggesting a model in which splicing can occur post-transcriptionally but still within the proximity of the transcription site, thus seeming co-transcriptional by most assays. These results may unify conflicting reports of co-transcriptional versus post-transcriptional splicing.

TOP1 and R-loops facilitate transcriptional DSBs at hypertranscribed cancer driver genes

DNA double-stranded breaks (DSBs) pose a significant threat to genomic integrity, and their generation during essential cellular processes like transcription remains poorly understood. In this study, we employ several techniques to map DSBs, R-loops, and topoisomerase 1 cleavage complex (TOP1cc) to comprehensively investigate the interplay between transcription, DSBs, topoisomerase 1 (TOP1), and R-loops. Our findings reveal the presence of DSBs at highly expressed genes enriched with TOP1 and R-loops. Remarkably, transcription-associated DSBs at these loci are significantly reduced upon depletion of R-loops and TOP1, uncovering the pivotal roles of TOP1 and R-loops in transcriptional DSB formation. By elucidating the intricate interplay between TOP1cc trapping, R-loops, and DSBs, our study provides insights into the mechanisms underlying transcription-associated genomic instability. Moreover, we establish a link between transcriptional DSBs and early molecular changes driving cancer development, highlighting the distinct etiology and molecular characteristics of driver mutations compared to passenger mutations.

The chromatin-associated lncREST ensures effective replication stress response by promoting the assembly of fork signaling factors

Besides the well-characterized protein network involved in the replication stress response, several regulatory RNAs have been shown to play a role in this critical process. However, it has remained elusive whether they act locally at the stressed forks. Here, by investigating the RNAs localizing on chromatin upon replication stress induced by hydroxyurea, we identified a set of lncRNAs upregulated in S-phase and controlled by stress transcription factors. Among them, we demonstrate that the previously uncharacterized lncRNA lncREST (long non-coding RNA REplication STress) is transcriptionally controlled by p53 and localizes at stressed replication forks. LncREST-depleted cells experience sustained replication fork progression and accumulate un-signaled DNA damage. Under replication stress, lncREST interacts with the protein NCL and assists in engaging its interaction with RPA. The loss of lncREST is associated with a reduced NCL-RPA interaction and decreased RPA on chromatin, leading to defective replication stress signaling and accumulation of mitotic defects, resulting in apoptosis and a reduction in tumorigenic potential of cancer cells. These findings uncover the function of a lncRNA in favoring the recruitment of replication proteins to sites of DNA replication.

ORC1 binds to cis-transcribed RNAs for efficient activation of replication origins

Cells must coordinate the activation of thousands of replication origins dispersed throughout their genome. Active transcription is known to favor the formation of mammalian origins, although the role that RNA plays in this process remains unclear. We show that the ORC1 subunit of the human Origin Recognition Complex interacts with RNAs transcribed from genes with origins in their transcription start sites (TSSs), displaying a positive correlation between RNA binding and origin activity. RNA depletion, or the use of ORC1 RNA-binding mutant, result in inefficient activation of proximal origins, linked to impaired ORC1 chromatin release. ORC1 RNA binding activity resides in its intrinsically disordered region, involved in intra- and inter-molecular interactions, regulation by phosphorylation, and phase-separation. We show that RNA binding favors ORC1 chromatin release, by regulating its phosphorylation and subsequent degradation. Our results unveil a non-coding function of RNA as a dynamic component of the chromatin, orchestrating the activation of replication origins.

SINEUP non-coding RNA activity depends on specific N6-methyladenosine nucleotides

SINEUPs are natural and synthetic antisense long non-coding RNAs (lncRNAs) selectively enhancing target mRNAs translation by increasing their association with polysomes. This activity requires two RNA domains: an embedded inverted SINEB2 element acting as effector domain, and an antisense region, the binding domain, conferring target selectivity. SINEUP technology presents several advantages to treat genetic (haploinsufficiencies) and complex diseases restoring the physiological activity of diseased genes and of compensatory pathways. To streamline these applications to the clinic, a better understanding of the mechanism of action is needed. Here we show that natural mouse SINEUP AS Uchl1 and synthetic human miniSINEUP-DJ-1 are N⁶-methyladenosine (m⁶A) modified by METTL3 enzyme. Then, we map m⁶A-modified sites along SINEUP sequence with Nanopore direct RNA sequencing and a reverse transcription assay. We report that m⁶A removal from SINEUP RNA causes the depletion of endogenous target mRNA from actively translating polysomes, without altering SINEUP enrichment in ribosomal subunit-associated fractions. These results prove that SINEUP activity requires an m⁶A-dependent step to enhance translation of target mRNAs, providing a new mechanism for m⁶A translation regulation and strengthening our knowledge of SINEUP-specific mode of action. Altogether these new findings pave the way to a more effective therapeutic application of this well-defined class of lncRNAs.

A genome-wide CRISPR activation screen identifies SCREEM a novel SNAI1 super-enhancer demarcated by eRNAs

The genome is pervasively transcribed to produce a vast array of non-coding RNAs (ncRNAs). Long non-coding RNAs (lncRNAs) are transcripts of >200 nucleotides and are best known for their ability to regulate gene expression. Enhancer RNAs (eRNAs) are subclass of lncRNAs that are synthesized from enhancer regions and have also been shown to coordinate gene expression. The biological function and significance of most lncRNAs and eRNAs remain to be determined. Epithelial to mesenchymal transition (EMT) is a ubiquitous cellular process that occurs during cellular migration, homeostasis, fibrosis, and cancer-cell metastasis. EMT-transcription factors, such as SNAI1 induce a complex transcriptional program that coordinates the morphological and molecular changes associated with EMT. Such complex transcriptional programs are often subject to coordination by networks of ncRNAs and thus can be leveraged to identify novel functional ncRNA loci. Here, using a genome-wide CRISPR activation (CRISPRa) screen targeting ∼10,000 lncRNA loci we identified ncRNA loci that could either promote or attenuate EMT. We discovered a novel locus that we named SCREEM (SNAI1 cis-regulatory eRNAs expressed in monocytes). The SCREEM locus contained a cluster of eRNAs that when activated using CRISPRa induced expression of the neighboring gene SNAI1, driving concomitant EMT. However, the SCREEM eRNA transcripts themselves appeared dispensable for the induction of SNAI1 expression. Interestingly, the SCREEM eRNAs and SNAI1 were co-expressed in activated monocytes, where the SCREEM locus demarcated a monocyte-specific super-enhancer. These findings suggest a potential role for SNAI1 in monocytes. Exploration of the SCREEM-SNAI axis could reveal novel aspects of monocyte biology.

Depletion of Tip60/Kat5 affects the hepatic antioxidant system in mice

Tat-interactive protein 60 kDa (TIP60, also known as lysine acetyltransferase 5 [KAT5]) is a member of the MYST protein family with histone acetyltransferase activity. Recent studies have reported that TIP60 has multiple functions in many signal transduction mechanisms, especially p53-mediated apoptosis. Although the activation of apoptosis signaling pathways requires the presence of cellular reactive oxygen species (ROS) at a certain level, an imbalance between the production and consumption of ROS in cells results in oxidative stress (OS). In this study, we investigated for the first time how the absence of the Tip60 gene in the liver affects gene expression, enzyme activity, and protein expression of the hepatic antioxidant members localized in the cytoplasm, including superoxide dismutase (SOD), catalase (CAT), glutathione reductase (GR), glutathione peroxidase (GPx), and glutathione S-transferase (GST). First, we successfully generated liver-specific Tip60 knockout mice (mutants) using Cre/LoxP recombination. The reduced glutathione level and nicotinamide adenine dinucleotide phosphate oxidase 4 (Nox4) expression, a marker of OS, increased significantly in the Tip60 mutant liver. Gene expression, activity, and protein expression of the enzymatic antioxidant system, including SOD, CAT, GR, GPx, and GST were investigated in mutants and control groups. Despite a significant correlation between the gene, enzyme activity, and protein content for CAT and GR, this was not true for SOD and GPx. The overall results suggest that TIP60 acts on the hepatic antioxidant system both at the gene and protein levels, but the actual effect of the deletion of Tip60 is observed at the protein level, especially for SOD and GPx.

Non-coding RNAs associated with Prader-Willi syndrome regulate transcription of neurodevelopmental genes in human induced pluripotent stem cells

Mutations and aberrant gene expression during cellular differentiation lead to neurodevelopmental disorders, such as Prader-Willi syndrome (PWS), which results from the deletion of an imprinted locus on paternally inherited chromosome 15. We analyzed chromatin-associated RNA in human induced pluripotent cells (iPSCs) upon depletion of hybrid small nucleolar long non-coding RNAs (sno-lncRNAs) and 5' snoRNA capped and polyadenylated long non-coding RNAs (SPA-lncRNAs) transcribed from the locus deleted in PWS. We found that rapid ablation of these lncRNAs affects transcription of specific gene classes. Downregulated genes contribute to neurodevelopment and neuronal maintenance, while upregulated genes are predominantly involved in the negative regulation of cellular metabolism and apoptotic processes. Our data reveal the importance of SPA-lncRNAs and sno-lncRNAs in controlling gene expression in iPSCs and provide a platform for synthetic experimental approaches in PWS studies. We conclude that ncRNAs transcribed from the PWS locus are critical regulators of a transcriptional signature, which is important for neuronal differentiation and development.

SUBCELLULAR TRANSCRIPTOMICS & PROTEOMICS: A COMPARATIVE METHODS REVIEW

The internal environment of cells is molecularly crowded, which requires spatial organization via subcellular compartmentalization. These compartments harbor specific conditions for molecules to perform their biological functions, such as coordination of the cell cycle, cell survival, and growth. This compartmentalization is also not static, with molecules trafficking between these subcellular neighborhoods to carry out their functions. For example, some biomolecules are multifunctional, requiring an environment with differing conditions or interacting partners, and others traffic to export such molecules. Aberrant localization of proteins or RNA species has been linked to many pathological conditions, such as neurological, cancer, and pulmonary diseases. Differential expression studies in transcriptomics and proteomics are relatively common, but the majority have overlooked the importance of subcellular information. In addition, subcellular transcriptomics and proteomics data do not always colocate because of the biochemical processes that occur during and after translation, highlighting the complementary nature of these fields. In this review, we discuss and directly compare the current methods in spatial proteomics and transcriptomics, which include sequencing- and imaging-based strategies, to give the reader an overview of the current tools available. We also discuss current limitations of these strategies as well as future developments in the field of spatial -omics.

•

Subcellular information of protein and RNA give insights into molecular function.

•

This review discusses strategies available to measure subcellular information.

•

Hybridization of methods shows promise for exploring the composition of organelles.

•

Advances are aiding understanding of the organisation and dynamics of cells.

Optimization of library preparation based on SMART for ultralow RNA-seq in mice brain tissues

BACKGROUND

Single-cell RNA sequencing (scRNA-seq) provides new insights to address biological and medical questions, and it will benefit more from the ultralow input RNA or subcellular sequencing.

RESULTS

Here, we present a highly sensitive library construction protocol for ultralow input RNA sequencing (ulRNA-seq). We systematically evaluate experimental conditions of this protocol, such as reverse transcriptase, template-switching oligos (TSO), and template RNA structure. It was found that Maxima H Minus reverse transcriptase and rN modified TSO, as well as all RNA templates capped with m7G improved the sequencing sensitivity and low abundance gene detection ability. RNA-seq libraries were successfully prepared from total RNA samples as low as 0.5 pg, and more than 2000 genes have been identified.

CONCLUSIONS

The ability of low abundance gene detection and sensitivity were largely enhanced with this optimized protocol. It was also confirmed in single-cell sequencing, that more genes and cell markers were identified compared to conventional sequencing method. We expect that ulRNA-seq will sequence and transcriptome characterization for the subcellular of disease tissue, to find the corresponding treatment plan.

RNALocate v2.0: an updated resource for RNA subcellular localization with increased coverage and annotation

Resolving the spatial distribution of the transcriptome at a subcellular level can increase our understanding of biology and diseases. To facilitate studies of biological functions and molecular mechanisms in the transcriptome, we updated RNALocate, a resource for RNA subcellular localization analysis that is freely accessible at http://www.rnalocate.org/ or http://www.rna-society.org/rnalocate/. Compared to RNALocate v1.0, the new features in version 2.0 include (i) expansion of the data sources and the coverage of species; (ii) incorporation and integration of RNA-seq datasets containing information about subcellular localization; (iii) addition and reorganization of RNA information (RNA subcellular localization conditions and descriptive figures for method, RNA homology information, RNA interaction and ncRNA disease information) and (iv) three additional prediction tools: DM3Loc, iLoc-lncRNA and iLoc-mRNA. Overall, RNALocate v2.0 provides a comprehensive RNA subcellular localization resource for researchers to deconvolute the highly complex architecture of the cell.

RNA-bound PGC-1α controls gene expression in liquid-like nuclear condensates

Plasticity of cells, tissues, and organs is controlled by the coordinated transcription of biological programs. However, the mechanisms orchestrating such context-specific transcriptional networks mediated by the dynamic interplay of transcription factors and coregulators are poorly understood. The peroxisome proliferator-activated receptor γ coactivator 1α (PGC-1α) is a prototypical master regulator of adaptive transcription in various cell types. We now uncovered a central function of the C-terminal domain of PGC-1α to bind RNAs and assemble multiprotein complexes including proteins that control gene transcription and RNA processing. These interactions are important for PGC-1α recruitment to chromatin in transcriptionally active liquid-like nuclear condensates. Notably, such a compartmentalization of active transcription mediated by liquid-liquid phase separation was observed in mouse and human skeletal muscle, revealing a mechanism by which PGC-1α regulates complex transcriptional networks. These findings provide a broad conceptual framework for context-dependent transcriptional control of phenotypic adaptations in metabolically active tissues.

Nascent RNA scaffolds contribute to chromosome territory architecture and counter chromatin compaction

Nuclear chromosomes transcribe far more RNA than required to encode protein. Here we investigate whether non-coding RNA broadly contributes to cytological-scale chromosome territory architecture. We develop a procedure that depletes soluble proteins, chromatin, and most nuclear RNA from the nucleus but does not delocalize XIST, a known architectural RNA, from an insoluble chromosome "scaffold." RNA-seq analysis reveals that most RNA in the nuclear scaffold is repeat-rich, non-coding, and derived predominantly from introns of nascent transcripts. Insoluble, repeat-rich (C₀T-1) RNA co-distributes with known scaffold proteins including scaffold attachment factor A (SAF-A), and distribution of these components inversely correlates with chromatin compaction in normal and experimentally manipulated nuclei. We further show that RNA is required for SAF-A to interact with chromatin and for enrichment of structurally embedded "scaffold attachment regions" prevalent in euchromatin. Collectively, the results indicate that long nascent transcripts contribute a dynamic structural role that promotes the open architecture of active chromosome territories.

The splicing factor XAB2 interacts with ERCC1-XPF and XPG for R-loop processing

RNA splicing, transcription and the DNA damage response are intriguingly linked in mammals but the underlying mechanisms remain poorly understood. Using an in vivo biotinylation tagging approach in mice, we show that the splicing factor XAB2 interacts with the core spliceosome and that it binds to spliceosomal U4 and U6 snRNAs and pre-mRNAs in developing livers. XAB2 depletion leads to aberrant intron retention, R-loop formation and DNA damage in cells. Studies in illudin S-treated cells and Csb^m/m developing livers reveal that transcription-blocking DNA lesions trigger the release of XAB2 from all RNA targets tested. Immunoprecipitation studies reveal that XAB2 interacts with ERCC1-XPF and XPG endonucleases outside nucleotide excision repair and that the trimeric protein complex binds RNA:DNA hybrids under conditions that favor the formation of R-loops. Thus, XAB2 functionally links the spliceosomal response to DNA damage with R-loop processing with important ramifications for transcription-coupled DNA repair disorders.

XPA-binding protein (XAB)-2 is the human homologue of the yeast pre-mRNA splicing factor Syf1. Here the authors use an in vivo biotinylation tagging approach to show XAB2’s role in DNA repair, RNA splicing and transcription during mammalian development.

Co-transcriptional splicing regulates 3' end cleavage during mammalian erythropoiesis

Pre-mRNA processing steps are tightly coordinated with transcription in many organisms. To determine how co-transcriptional splicing is integrated with transcription elongation and 3' end formation in mammalian cells, we performed long-read sequencing of individual nascent RNAs and precision run-on sequencing (PRO-seq) during mouse erythropoiesis. Splicing was not accompanied by transcriptional pausing and was detected when RNA polymerase II (Pol II) was within 75-300 nucleotides of 3' splice sites (3'SSs), often during transcription of the downstream exon. Interestingly, several hundred introns displayed abundant splicing intermediates, suggesting that splicing delays can take place between the two catalytic steps. Overall, splicing efficiencies were correlated among introns within the same transcript, and intron retention was associated with inefficient 3' end cleavage. Remarkably, a thalassemia patient-derived mutation introducing a cryptic 3'SS improved both splicing and 3' end cleavage of individual β-globin transcripts, demonstrating functional coupling between the two co-transcriptional processes as a determinant of productive gene output.

Global Positioning System: Understanding Long Noncoding RNAs through Subcellular Localization

The localization of long noncoding RNAs (lncRNAs) within the cell is the primary determinant of their molecular functions. LncRNAs are often thought of as chromatin-restricted regulators of gene transcription and chromatin structure. However, a rich population of cytoplasmic lncRNAs has come to light, with diverse roles including translational regulation, signaling, and respiration. RNA maps of increasing resolution and scope are revealing a subcellular world of highly specific localization patterns and hint at sequence-based address codes specifying lncRNA fates. We propose a new framework for analyzing sequencing-based data, which suggests that numbers of cytoplasmic lncRNA molecules rival those in the nucleus. New techniques promise to create high-resolution, transcriptome-wide maps associated with all organelles of the mammalian cell. Given its intimate link to molecular roles, subcellular localization provides a means of unlocking the mystery of lncRNA functions.