Spliceosomal snRNA Epitranscriptomics

m6A modification of U6 snRNA modulates usage of two major classes of pre-mRNA 5' splice site

Alternative splicing of messenger RNAs is associated with the evolution of developmentally complex eukaryotes. Splicing is mediated by the spliceosome, and docking of the pre-mRNA 5' splice site into the spliceosome active site depends upon pairing with the conserved ACAGA sequence of U6 snRNA. In some species, including humans, the central adenosine of the ACAGA box is modified by N6 methylation, but the role of this m6A modification is poorly understood. Here, we show that m6A modified U6 snRNA determines the accuracy and efficiency of splicing. We reveal that the conserved methyltransferase, FIONA1, is required for Arabidopsis U6 snRNA m6A modification. Arabidopsis fio1 mutants show disrupted patterns of splicing that can be explained by the sequence composition of 5' splice sites and cooperative roles for U5 and U6 snRNA in splice site selection. U6 snRNA m6A influences 3' splice site usage. We generalise these findings to reveal two major classes of 5' splice site in diverse eukaryotes, which display anti-correlated interaction potential with U5 snRNA loop 1 and the U6 snRNA ACAGA box. We conclude that U6 snRNA m6A modification contributes to the selection of degenerate 5' splice sites crucial to alternative splicing.

Implications of m6A-associated snRNAs in the prognosis and immunotherapeutic responses of hepatocellular carcinoma

Background

Hepatocellular carcinoma (HCC) is the most prevalent pathological type of liver cancer worldwide with high mortality and poor prognosis. N6-methyladenosine (m6A) can modify RNAs such as mRNA, lncRNA, miRNA, and tRNA, thereby playing a critical role in the pathogenesis of HCC. However, the role of m6A-associated small nuclear RNA (snRNA) in the prognostic value and immunotherapeutic response in HCC remains unclear.

Materials and methods

In this study, snRNA expression data, gene mutation data, and clinical data of HCC patients were acquired from The Cancer Genome Atlas (TCGA) database. We used the least absolute shrinkage and selection operator (LASSO) Cox regression analysis to identify significant prognostic m6A-associated snRNAs, and then developed a multivariate Cox model based on the selected snRNAs. HCC patients were split into low- and high-risk groups based on the median risk score. We subsequently performed Kaplan-Meier curve analysis to estimate overall survival (OS) by clinicopathological characteristics and tumor mutational burden (TMB) status in low- and high-risk HCC patients. Finally, we compared the immunotherapeutic response as represented by tumor immune dysfunction and exclusion (TIDE) scores between the two risk groups.

Results

Eight m6A-associated snRNAs were selected as independent predictors to develop the risk model. Our results revealed that the OS of HCC patients in the high-risk group was significantly worse than that in the low-risk group on clinicopathologic characteristics, including age (≤65 years and >65 years), gender (male), grade (G I-II and G III-IV) and TNM staging (Stage I-II and Stage III-IV). In addition, the OS of low-TMB and low-risk group was longer than that of high-TMB and high-risk group. The TIDE score indicated that HCC patients in the high-risk group were more susceptible to immunotherapy.

Conclusion

Our study suggests that m6A-associated snRNAs may be useful biomarkers for the prognosis of HCC and that m6A-associated snRNA models can predict the effect of immunotherapy in HCC patients.

The coilin N-terminus mediates multivalent interactions between coilin and Nopp140 to form and maintain Cajal bodies

Cajal bodies (CBs) are ubiquitous nuclear membraneless organelles (MLOs) that concentrate and promote efficient biogenesis of snRNA-protein complexes involved in splicing (snRNPs). Depletion of the CB scaffolding protein coilin disperses snRNPs, making CBs a model system for studying the structure and function of MLOs. Although it is assumed that CBs form through condensation, the biomolecular interactions responsible remain elusive. Here, we discover the unexpected capacity of coilin’s N-terminal domain (NTD) to form extensive fibrils in the cytoplasm and discrete nuclear puncta in vivo. Single amino acid mutational analysis reveals distinct molecular interactions between coilin NTD proteins to form fibrils and additional NTD interactions with the nuclear Nopp140 protein to form puncta. We provide evidence that Nopp140 has condensation capacity and is required for CB assembly. From these observations, we propose a model in which coilin NTD–NTD mediated assemblies make multivalent contacts with Nopp140 to achieve biomolecular condensation in the nucleus.

Cajal bodies are membraneless organelles scaffolded by coilin protein. Here, coilin–coilin and coilin–Nopp140 interaction sites are identified and perturbed, revealing coilin’s capacity to form long fibrils and be remodeled into spherical structures.

The PAICE suite reveals circadian posttranscriptional timing of noncoding RNAs and spliceosome components in Mus musculus macrophages

Circadian rhythms broadly regulate physiological functions by tuning oscillations in the levels of mRNAs and proteins to the 24-h day/night cycle. Globally assessing which mRNAs and proteins are timed by the clock necessitates accurate recognition of oscillations in RNA and protein data, particularly in large omics data sets. Tools that employ fixed-amplitude models have previously been used to positive effect. However, the recognition of amplitude change in circadian oscillations required a new generation of analytical software to enhance the identification of these oscillations. To address this gap, we created the Pipeline for Amplitude Integration of Circadian Exploration suite. Here, we demonstrate the Pipeline for Amplitude Integration of Circadian Exploration suite's increased utility to detect circadian trends through the joint modeling of the Mus musculus macrophage transcriptome and proteome. Our enhanced detection confirmed extensive circadian posttranscriptional regulation in macrophages but highlighted that some of the reported discrepancy between mRNA and protein oscillations was due to noise in data. We further applied the Pipeline for Amplitude Integration of Circadian Exploration suite to investigate the circadian timing of noncoding RNAs, documenting extensive circadian timing of long noncoding RNAs and small nuclear RNAs, which control the recognition of mRNA in the spliceosome complex. By tracking oscillating spliceosome complex proteins using the PAICE suite, we noted that the clock broadly regulates the spliceosome, particularly the major spliceosome complex. As most of the above-noted rhythms had damped amplitude changes in their oscillations, this work highlights the importance of the PAICE suite in the thorough enumeration of oscillations in omics-scale datasets.

Current Status of Regulatory Non-Coding RNAs Research in the Tritryp

Trypanosomatids are protozoan parasites that cause devastating vector-borne human diseases. Gene expression regulation of these organisms depends on post-transcriptional control in responding to diverse environments while going through multiple developmental stages of their complex life cycles. In this scenario, non-coding RNAs (ncRNAs) are excellent candidates for a very efficient, quick, and economic strategy to regulate gene expression. The advent of high throughput RNA sequencing technologies show the presence and deregulation of small RNA fragments derived from canonical ncRNAs. This review seeks to depict the ncRNA landscape in trypanosomatids, focusing on the small RNA fragments derived from functional RNA molecules observed in RNA sequencing studies. Small RNA fragments derived from canonical ncRNAs (tsRNAs, snsRNAs, sdRNAs, and sdrRNAs) were identified in trypanosomatids. Some of these RNAs display changes in their levels associated with different environments and developmental stages, demanding further studies to determine their functional characterization and potential roles. Nevertheless, a comprehensive and detailed ncRNA annotation for most trypanosomatid genomes is still needed, allowing better and more extensive comparative and functional studies.

PAMP-triggered genetic reprogramming involves widespread alternative transcription initiation and an immediate transcription factor wave

Immune responses triggered by pathogen-associated molecular patterns (PAMPs) are key to pathogen defense, but drivers and stabilizers of the growth-to-defense genetic reprogramming remain incompletely understood in plants. Here, we report a time-course study of the establishment of PAMP-triggered immunity (PTI) using cap analysis of gene expression. We show that around 15% of all transcription start sites (TSSs) rapidly induced during PTI define alternative transcription initiation events. From these, we identify clear examples of regulatory TSS change via alternative inclusion of target peptides or domains in encoded proteins, or of upstream open reading frames in mRNA leader sequences. We also find that 60% of PAMP response genes respond earlier than previously thought. In particular, a cluster of rapidly and transiently PAMP-induced genes is enriched in transcription factors (TFs) whose functions, previously associated with biological processes as diverse as abiotic stress adaptation and stem cell activity, appear to converge on growth restriction. Furthermore, examples of known potentiators of PTI, in one case under direct mitogen-activated protein kinase control, support the notion that the rapidly induced TFs could constitute direct links to PTI signaling pathways and drive gene expression changes underlying establishment of the immune state.

The plant epitranscriptome: revisiting pseudouridine and 2'-O-methyl RNA modifications

There is growing evidence that post-transcriptional RNA modifications are highly dynamic and can be used to improve crop production. Although more than 172 unique types of RNA modifications have been identified throughout the kingdom of life, we are yet to leverage upon the understanding to optimize RNA modifications in crops to improve productivity. The contributions of internal mRNA modifications such as N6-methyladenosine (m⁶ A) and 5-methylcytosine (m⁵ C) methylations to embryonic development, root development, leaf morphogenesis, flowering, fruit ripening and stress response are sufficiently known, but the roles of the two most abundant RNA modifications, pseudouridine (Ψ) and 2'-O-methylation (Nm), in the cell remain unclear due to insufficient advances in high-throughput technologies in plant development. Therefore, in this review, we discuss the latest methods and insights gained in mapping internal Ψ and Nm and their unique properties in plants and other organisms. In addition, we discuss the limitations that remain in high-throughput technologies for qualitative and quantitative mapping of these RNA modifications and highlight future challenges in regulating the plant epitranscriptome.

Structure Characterization of Escherichia coli Pseudouridine Kinase PsuK

Pseudouridine (Ψ) is one of the most abundant RNA modifications in cellular RNAs that post-transcriptionally impact many aspects of RNA. However, the metabolic fate of modified RNA nucleotides has long been a question. A pseudouridine kinase (PsuK) and a pseudouridine monophosphate glycosylase (PsuG) in Escherichia coli were first characterized as involved in pseudouridine degradation by catalyzing the phosphorylation of pseudouridine to pseudouridine 5′-phosphate (ΨMP) and further hydrolyzing 5′-ΨMP to produce uracil and ribose 5′-phosphate. Recently, their homolog proteins in eukaryotes were also identified, which were named PUKI and PUMY in Arabidopsis. Here, we solved the crystal structures of apo-EcPsuK and its binary complex with Ψ or N¹-methyl-pseudouridine (m1Ψ). The structure of EcPsuK showed a homodimer conformation assembled by its β-thumb region. EcPsuK has an appropriate binding site with a series of hydrophilic and hydrophobic interactions for Ψ. Moreover, our complex structure of EcPsuK-m1Ψ suggested the binding pocket has an appropriate capacity for m1Ψ. We also identified the monovalent ion-binding site and potential ATP-binding site. Our studies improved the understanding of the mechanism of Ψ turnover.

Identification and functional implications of pseudouridine RNA modification on small noncoding RNAs in the mammalian pathogen Trypanosoma brucei

Trypanosoma brucei, the parasite that causes sleeping sickness, cycles between an insect and a mammalian host. However, the effect of RNA modifications such as pseudouridinylation on its ability to survive in these two different host environments is unclear. Here, two genome-wide approaches were applied for mapping pseudouridinylation sites (Ψs) on small nucleolar RNA (snoRNA), 7SL RNA, vault RNA, and tRNAs from T. brucei. We show using HydraPsiSeq and RiboMeth-seq that the Ψ on C/D snoRNA guiding 2'-O-methylation increased the efficiency of the guided modification on its target, rRNA. We found differential levels of Ψs on these noncoding RNAs in the two life stages (insect host and mammalian host) of the parasite. Furthermore, tRNA isoform abundance and Ψ modifications were characterized in these two life stages demonstrating stage-specific regulation. We conclude that the differential Ψ modifications identified here may contribute to modulating the function of noncoding RNAs involved in rRNA processing, rRNA modification, protein synthesis, and protein translocation during cycling of the parasite between its two hosts.

The Progress and Promise of RNA Medicine─An Arsenal of Targeted Treatments

In the past decade, there has been a shift in research, clinical development, and commercial activity to exploit the many physiological roles of RNA for use in medicine. With the rapid success in the development of lipid-RNA nanoparticles for mRNA vaccines against COVID-19 and with several approved RNA-based drugs, RNA has catapulted to the forefront of drug research. With diverse functions beyond the role of mRNA in producing antigens or therapeutic proteins, many classes of RNA serve regulatory roles in cells and tissues. These RNAs have potential as new therapeutics, with RNA itself serving as either a drug or a target. Here, based on the CAS Content Collection, we provide a landscape view of the current state and outline trends in RNA research in medicine across time, geography, therapeutic pipelines, chemical modifications, and delivery mechanisms.

Changes in the Small Noncoding RNAome During M1 and M2 Macrophage Polarization

Macrophages belong to a special phagocytic subgroup of human leukocytes and are one of the important cells of the human immune system. Small noncoding RNAs are a group of small RNA molecules that can be transcribed without the ability to encode proteins but could play a specific function in cells. SncRNAs mainly include microRNAs (miRNAs) and piwi-interacting RNAs (piRNAs), small nucleolar RNAs (snoRNAs), small nuclear RNAs (snRNAs) and repeat RNAs. We used high-throughput sequencing analysis and qPCR to detect the expression changes of the small noncoding RNAome during macrophage polarization. Our results showed that 84 miRNAs and 47 miRNAs with were downregulated during M1 macrophage polarization and that 11 miRNAs were upregulated and 19 miRNAs were downregulated during M2 macrophage polarization. MiR-novel-3-nature and miR-27b-5p could promote expression of TNF-α which was marker gene of M1 macrophages. The piRNA analysis results showed that 69 piRNAs were upregulated and 61 piRNAs were downregulated during M1 macrophage polarization and that 3 piRNAs were upregulated and 10 piRNAs were downregulated during M2 macrophage polarization. DQ551351 and DQ551308 could promote the mRNA expression of TNF-α and DQ551351overexpression promoted the antitumor activity of M1 macrophages. SnoRNA results showed that 62 snoRNAs were upregulated and 59 snoRNAs were downregulated during M1 macrophage polarization, whereas 6 snoRNAs were upregulated and 10 snoRNAs were downregulated during M2 macrophage polarization. Overexpression of snoRNA ENSMUST00000158683.2 could inhibit expression of TNF-α. For snRNA, we found that 12 snRNAs were upregulated and 15 snRNAs were downregulated during M1 macrophage polarization and that 2 snRNAs were upregulated during M2 macrophage polarization. ENSMUSG00000096786 could promote expression of IL-1 and iNOS and ENSMUSG00000096786 overexpression promoted the antitumor activity of M1 macrophages. Analysis of repeat RNAs showed that 7 repeat RNAs were upregulated and 9 repeat RNAs were downregulated during M1 macrophage polarization and that 2 repeat RNAs were downregulated during M2 macrophage polarization. We first reported the expression changes of piRNA, snoRNA, snRNA and repeat RNA during macrophage polarization, and preliminarily confirmed that piRNA, snoRNA and snRNA can regulate the function of macrophages.

Guide RNA acrobatics: the one-for-two shuffle

RNA modifications are crucial for the proper function of the RNAs. The sites of pseudouridines are often specified by dual hairpin guide RNAs, with one or both hairpins identifying a target uridine. In this issue of Genes & Development, Jády and colleagues (pp. 70-83) identify a novel mechanism by which a single guide RNA hairpin can specify two uridines adjacent to each other or separated by 1 nt; i.e., one for two or guide RNA acrobatics.

Epitranscriptomic regulation of cognitive development and decline

Functional genomics and systems biology have opened new doors to previously inaccessible genomic information and holistic approaches to study complex networks of genes and proteins in the central nervous system. The advances are revolutionizing our understanding of the genetic underpinning of cognitive development and decline by facilitating identifications of novel molecular regulators and physiological pathways underlying brain function, and by associating polymorphism and mutations to cognitive dysfunction and neurological diseases. However, our current understanding of these complex gene regulatory mechanisms has yet lacked sufficient mechanistic resolution for further translational breakthroughs. Here we review recent findings from the burgeoning field of epitranscriptomics in association of cognitive functions with a special focus on the epitranscritomic regulation in subcellular locations such as chromosome, synapse, and mitochondria. Although there are important gaps in knowledge, current evidence is suggesting that this layer of RNA regulation may be of particular interest for the spatiotemporally coordinated regulation of gene networks in developing and maintaining brain function that underlie cognitive changes.

Combined exposure to polychlorinated biphenyls and high-fat diet modifies the global epitranscriptomic landscape in mouse liver

Exposure to a single dose of polychlorinated biphenyls (PCBs) and a 12-week high-fat diet (HFD) results in nonalcoholic steatohepatitis (NASH) in mice by altering intracellular signaling and inhibiting epidermal growth factor receptor signaling. Post-transcriptional chemical modification (PTM) of RNA regulates biological processes, but the contribution of epitranscriptomics to PCB-induced steatosis remains unknown. This study tested the hypothesis that PCB and HFD exposure alters the global RNA epitranscriptome in male mouse liver. C57BL/6J male mice were fed a HFD for 12 weeks and exposed to a single dose of Aroclor 1260 (20 mg/kg), PCB 126 (20 µg/kg), both Aroclor 1260 and PCB 126 or vehicle control after 2 weeks on HFD. Chemical RNA modifications were identified at the nucleoside level by liquid chromatography-mass spectrometry. From 22 PTM global RNA modifications, we identified 10 significant changes in RNA modifications in liver with HFD and PCB 126 exposure. Only two modifications were significantly different from HFD control liver in all three PCB exposure groups: 2'-O-methyladenosine (Am) and N(6)-methyladenosine (m6A). Exposure to HFD + PCB 126 + Aroclor 1260 increased the abundance of N(6), O(2)-dimethyladenosine (m6Am), which is associated with the largest number of transcript changes. Increased m6Am and pseudouridine were associated with increased protein expression of the writers of these modifications: Phosphorylated CTD Interacting Factor 1 (PCIF1) and Pseudouridine Synthase 10 (PUS10), respectively, in HFD + PCB 126- + Aroclor 1260-exposed mouse liver. Increased N1-methyladenosine (m1A) and m6A were associated with increased transcript levels of the readers of these modifications: YTH N6-Methyladenosine RNA Binding Protein 2 (YTHDF2), YTH Domain Containing 2 (YTHDC2), and reader FMRP Translational Regulator 1 (FMR1) transcript and protein abundance. The results demonstrate that PCB exposure alters the global epitranscriptome in a mouse model of NASH; however, the mechanism for these changes requires further investigation.

Pseudouridines in RNAs: switching atoms means shifting paradigms

The structure, stability, and function of various coding and noncoding RNAs are influenced by chemical modifications. Pseudouridine (Ψ) is one of the most abundant post‐transcriptional RNA base modifications and has been detected at individual positions in tRNAs, rRNAs, mRNAs, and snRNAs, which are referred to as Ψ‐sites. By allowing formation of additional bonds with neighboring atoms, Ψ strengthens RNA–RNA and RNA–protein interactions. Although many aspects of the underlying modification reactions remain unclear, the advent of new transcriptome‐wide methods to quantitatively detect Ψ‐sites has recently changed our perception of the functional roles and importance of Ψ. For instance, it is now clear that the occurrence of Ψs appears to be directly linked to the lifetime and the translation efficiency of a given mRNA molecule. Furthermore, the administration of Ψ‐containing RNAs reduces innate immune responses, which appears strikingly advantageous for the development of generations of mRNA‐based vaccines. In this review, we aim to comprehensively summarize recent discoveries that highlight the impact of Ψ on various types of RNAs and outline possible novel biomedical applications of Ψ.

Nopp140-chaperoned 2'-O-methylation of small nuclear RNAs in Cajal bodies ensures splicing fidelity

In this study, Bizarro et al. sought to understand the function and subcellular site of snRNA modification, and found that Cajal body (CB) localization of the protein Nopp140 is essential for concentration of small Cajal body-specific ribonucleoproteins (scaRNPs) in nuclear condensate and that phosphorylation by casein kinase 2 (CK2) at ∼80 serines targets Nopp140 to CBs. Nopp140 knockdown-mediated release of scaRNPs from CBs severely compromises 2′-O-methylation of spliceosomal snRNAs, identifying CBs as the site of scaRNP catalysis.

Spliceosomal small nuclear RNAs (snRNAs) are modified by small Cajal body (CB)-specific ribonucleoproteins (scaRNPs) to ensure snRNP biogenesis and pre-mRNA splicing. However, the function and subcellular site of snRNA modification are largely unknown. We show that CB localization of the protein Nopp140 is essential for concentration of scaRNPs in that nuclear condensate; and that phosphorylation by casein kinase 2 (CK2) at ∼80 serines targets Nopp140 to CBs. Transiting through CBs, snRNAs are apparently modified by scaRNPs. Indeed, Nopp140 knockdown-mediated release of scaRNPs from CBs severely compromises 2′-O-methylation of spliceosomal snRNAs, identifying CBs as the site of scaRNP catalysis. Additionally, alternative splicing patterns change indicating that these modifications in U1, U2, U5, and U12 snRNAs safeguard splicing fidelity. Given the importance of CK2 in this pathway, compromised splicing could underlie the mode of action of small molecule CK2 inhibitors currently considered for therapy in cholangiocarcinoma, hematological malignancies, and COVID-19.

From Antisense RNA to RNA Modification: Therapeutic Potential of RNA-Based Technologies

Therapeutic oligonucleotides interact with a target RNA via Watson-Crick complementarity, affecting RNA-processing reactions such as mRNA degradation, pre-mRNA splicing, or mRNA translation. Since they were proposed decades ago, several have been approved for clinical use to correct genetic mutations. Three types of mechanisms of action (MoA) have emerged: RNase H-dependent degradation of mRNA directed by short chimeric antisense oligonucleotides (gapmers), correction of splicing defects via splice-modulation oligonucleotides, and interference of gene expression via short interfering RNAs (siRNAs). These antisense-based mechanisms can tackle several genetic disorders in a gene-specific manner, primarily by gene downregulation (gapmers and siRNAs) or splicing defects correction (exon-skipping oligos). Still, the challenge remains for the repair at the single-nucleotide level. The emerging field of epitranscriptomics and RNA modifications shows the enormous possibilities for recoding the transcriptome and repairing genetic mutations with high specificity while harnessing endogenously expressed RNA processing machinery. Some of these techniques have been proposed as alternatives to CRISPR-based technologies, where the exogenous gene-editing machinery needs to be delivered and expressed in the human cells to generate permanent (DNA) changes with unknown consequences. Here, we review the current FDA-approved antisense MoA (emphasizing some enabling technologies that contributed to their success) and three novel modalities based on post-transcriptional RNA modifications with therapeutic potential, including ADAR (Adenosine deaminases acting on RNA)-mediated RNA editing, targeted pseudouridylation, and 2′-O-methylation.

Nopp140-chaperoned 2'-O-methylation of small nuclear RNAs in Cajal bodies ensures splicing fidelity

Spliceosomal small nuclear RNAs (snRNAs) are modified by small Cajal body (CB) specific ribonucleoproteins (scaRNPs) to ensure snRNP biogenesis and pre-mRNA splicing. However, the function and subcellular site of snRNA modification are largely unknown. We show that CB localization of the protein Nopp140 is essential for concentration of scaRNPs in that nuclear condensate; and that phosphorylation by casein kinase 2 (CK2) at some 80 serines targets Nopp140 to CBs. Transiting through CBs, snRNAs are apparently modified by scaRNPs. Indeed, Nopp140 knockdown-mediated release of scaRNPs from CBs severely compromises 2'-O-methylation of spliceosomal snRNAs, identifying CBs as the site of scaRNP catalysis. Additionally, alternative splicing patterns change indicating that these modifications in U1, U2, U5, and U12 snRNAs safeguard splicing fidelity. Given the importance of CK2 in this pathway, compromised splicing could underlie the mode of action of small molecule CK2 inhibitors currently considered for therapy in cholangiocarcinoma, hematological malignancies, and COVID-19.