Pharmacologic modulation of RNA splicing enhances anti-tumor immunity.

Although mutations in DNA are the best-studied source of neoantigens that determine response to immune checkpoint blockade, alterations in RNA splicing within cancer cells could similarly result in neoepitope production. However, the endogenous antigenicity and clinical potential of such splicing-derived epitopes have not been tested. Here, we demonstrate that pharmacologic modulation of splicing via specific drug classes generates bona fide neoantigens and elicits anti-tumor immunity, augmenting checkpoint immunotherapy. Splicing modulation inhibited tumor growth and enhanced checkpoint blockade in a manner dependent on host T cells and peptides presented on tumor MHC class I. Splicing modulation induced stereotyped splicing changes across tumor types, altering the MHC I-bound immunopeptidome to yield splicing-derived neoepitopes that trigger an anti-tumor T cell response in vivo. These data definitively identify splicing modulation as an untapped source of immunogenic peptides and provide a means to enhance response to checkpoint blockade that is readily translatable to the clinic.

ERα is an RNA-binding protein sustaining tumor cell survival and drug resistance.

Estrogen receptor α (ERα) is a hormone receptor and key driver for over 70% of breast cancers that has been studied for decades as a transcription factor. Unexpectedly, we discover that ERα is a potent non-canonical RNA-binding protein. We show that ERα RNA binding function is uncoupled from its activity to bind DNA and critical for breast cancer progression. Employing genome-wide cross-linking immunoprecipitation (CLIP) sequencing and a functional CRISPRi screen, we find that ERα-associated mRNAs sustain cancer cell fitness and elicit cellular responses to stress. Mechanistically, ERα controls different steps of RNA metabolism. In particular, we demonstrate that ERα RNA binding mediates alternative splicing of XBP1 and translation of the eIF4G2 and MCL1 mRNAs, which facilitates survival upon stress conditions and sustains tamoxifen resistance of cancer cells. ERα is therefore a multifaceted RNA-binding protein, and this activity transforms our knowledge of post-transcriptional regulation underlying cancer development and drug response.

Spliceosome-targeted therapies trigger an antiviral immune response in triple-negative breast cancer.

Many oncogenic insults deregulate RNA splicing, often leading to hypersensitivity of tumors to spliceosome-targeted therapies (STTs). However, the mechanisms by which STTs selectively kill cancers remain largely unknown. Herein, we discover that mis-spliced RNA itself is a molecular trigger for tumor killing through viral mimicry. In MYC-driven triple-negative breast cancer, STTs cause widespread cytoplasmic accumulation of mis-spliced mRNAs, many of which form double-stranded structures. Double-stranded RNA (dsRNA)-binding proteins recognize these endogenous dsRNAs, triggering antiviral signaling and extrinsic apoptosis. In immune-competent models of breast cancer, STTs cause tumor cell-intrinsic antiviral signaling, downstream adaptive immune signaling, and tumor cell death. Furthermore, RNA mis-splicing in human breast cancers correlates with innate and adaptive immune signatures, especially in MYC-amplified tumors that are typically immune cold. These findings indicate that dsRNA-sensing pathways respond to global aberrations of RNA splicing in cancer and provoke the hypothesis that STTs may provide unexplored strategies to activate anti-tumor immune pathways.

Structures of the Catalytically Activated Yeast Spliceosome Reveal the Mechanism of Branching.

Pre-mRNA splicing is executed by the spliceosome. Structural characterization of the catalytically activated complex (B∗) is pivotal for understanding the branching reaction. In this study, we assembled the B∗ complexes on two different pre-mRNAs from Saccharomyces cerevisiae and determined the cryo-EM structures of four distinct B∗ complexes at overall resolutions of 2.9-3.8 Å. The duplex between U2 small nuclear RNA (snRNA) and the branch point sequence (BPS) is discretely away from the 5'-splice site (5'SS) in the three B∗ complexes that are devoid of the step I splicing factors Yju2 and Cwc25. Recruitment of Yju2 into the active site brings the U2/BPS duplex into the vicinity of 5'SS, with the BPS nucleophile positioned 4 Å away from the catalytic metal M2. This analysis reveals the functional mechanism of Yju2 and Cwc25 in branching. These structures on different pre-mRNAs reveal substrate-specific conformations of the spliceosome in a major functional state.

Lactate modulates RNA splicing to promote CTLA-4 expression in tumor-infiltrating regulatory T cells.

RNA splicing is involved in cancer initiation and progression, but how it influences host antitumor immunity in the metabolically abnormal tumor microenvironment (TME) remains unclear. Here, we demonstrate that lactate modulates Foxp3-dependent RNA splicing to maintain the phenotypic and functional status of tumor-infiltrating regulatory T (Treg) cells via CTLA-4. RNA splicing in Treg cells was correlated with the Treg cell signatures in the TME. Ubiquitin-specific peptidase 39 (USP39), a component of the RNA splicing machinery, maintained RNA-splicing-mediated CTLA-4 expression to control Treg cell function. Mechanistically, lactate promoted USP39-mediated RNA splicing to facilitate CTLA-4 expression in a Foxp3-dependent manner. Moreover, the efficiency of CTLA-4 RNA splicing was increased in tumor-infiltrating Treg cells from patients with colorectal cancer. These findings highlight the immunological relevance of RNA splicing in Treg cells and provide important insights into the environmental mechanism governing CTLA-4 expression in Treg cells.

Post-transcriptional Regulation of De Novo Lipogenesis by mTORC1-S6K1-SRPK2 Signaling.

mTORC1 is a signal integrator and master regulator of cellular anabolic processes linked to cell growth and survival. Here, we demonstrate that mTORC1 promotes lipid biogenesis via SRPK2, a key regulator of RNA-binding SR proteins. mTORC1-activated S6K1 phosphorylates SRPK2 at Ser494, which primes Ser497 phosphorylation by CK1. These phosphorylation events promote SRPK2 nuclear translocation and phosphorylation of SR proteins. Genome-wide transcriptome analysis reveals that lipid biosynthetic enzymes are among the downstream targets of mTORC1-SRPK2 signaling. Mechanistically, SRPK2 promotes SR protein binding to U1-70K to induce splicing of lipogenic pre-mRNAs. Inhibition of this signaling pathway leads to intron retention of lipogenic genes, which triggers nonsense-mediated mRNA decay. Genetic or pharmacological inhibition of SRPK2 blunts de novo lipid synthesis, thereby suppressing cell growth. These results thus reveal a novel role of mTORC1-SRPK2 signaling in post-transcriptional regulation of lipid metabolism and demonstrate that SRPK2 is a potential therapeutic target for mTORC1-driven metabolic disorders.

CAG repeat expansions create splicing acceptor sites and produce aberrant repeat-containing RNAs.

Expansions of CAG trinucleotide repeats cause several rare neurodegenerative diseases. The disease-causing repeats are translated in multiple reading frames and without an identifiable initiation codon. The molecular mechanism of this repeat-associated non-AUG (RAN) translation is not known. We find that expanded CAG repeats create new splice acceptor sites. Splicing of proximal donors to the repeats produces unexpected repeat-containing transcripts. Upon splicing, depending on the sequences surrounding the donor, CAG repeats may become embedded in AUG-initiated open reading frames. Canonical AUG-initiated translation of these aberrant RNAs may account for proteins that have been attributed to RAN translation. Disruption of the relevant splice donors or the in-frame AUG initiation codons is sufficient to abrogate RAN translation. Our findings provide a molecular explanation for the abnormal translation products observed in CAG trinucleotide repeat expansion disorders and add to the repertoire of mechanisms by which repeat expansion mutations disrupt cellular functions.

FAM120A couples SREBP-dependent transcription and splicing of lipogenesis enzymes downstream of mTORC1.

The mechanistic target of rapamycin complex 1 (mTORC1) is a master regulator of cell growth that stimulates macromolecule synthesis through transcription, RNA processing, and post-translational modification of metabolic enzymes. However, the mechanisms of how mTORC1 orchestrates multiple steps of gene expression programs remain unclear. Here, we identify family with sequence similarity 120A (FAM120A) as a transcription co-activator that couples transcription and splicing of de novo lipid synthesis enzymes downstream of mTORC1-serine/arginine-rich protein kinase 2 (SRPK2) signaling. The mTORC1-activated SRPK2 phosphorylates splicing factor serine/arginine-rich splicing factor 1 (SRSF1), enhancing its binding to FAM120A. FAM120A directly interacts with a lipogenic transcription factor SREBP1 at active promoters, thereby bridging the newly transcribed lipogenic genes from RNA polymerase II to the SRSF1 and U1-70K-containing RNA-splicing machinery. This mTORC1-regulated, multi-protein complex promotes efficient splicing and stability of lipogenic transcripts, resulting in fatty acid synthesis and cancer cell proliferation. These results elucidate FAM120A as a critical transcription co-factor that connects mTORC1-dependent gene regulation programs for anabolic cell growth.

SRSF2 plays an unexpected role as reader of m5C on mRNA, linking epitranscriptomics to cancer.

A common mRNA modification is 5-methylcytosine (m5C), whose role in gene-transcript processing and cancer remains unclear. Here, we identify serine/arginine-rich splicing factor 2 (SRSF2) as a reader of m5C and impaired SRSF2 m5C binding as a potential contributor to leukemogenesis. Structurally, we identify residues involved in m5C recognition and the impact of the prevalent leukemia-associated mutation SRSF2P95H. We show that SRSF2 binding and m5C colocalize within transcripts. Furthermore, knocking down the m5C writer NSUN2 decreases mRNA m5C, reduces SRSF2 binding, and alters RNA splicing. We also show that the SRSF2P95H mutation impairs the ability of the protein to read m5C-marked mRNA, notably reducing its binding to key leukemia-related transcripts in leukemic cells. In leukemia patients, low NSUN2 expression leads to mRNA m5C hypomethylation and, combined with SRSF2P95H, predicts poor outcomes. Altogether, we highlight an unrecognized mechanistic link between epitranscriptomics and a key oncogenesis driver.

Mechanism of exon ligation by human spliceosome.

Pre-mRNA splicing involves two sequential reactions: branching and exon ligation. The C complex after branching undergoes remodeling to become the C∗ complex, which executes exon ligation. Here, we report cryo-EM structures of two intermediate human spliceosomal complexes, pre-C∗-I and pre-C∗-II, both at 3.6 Å. In both structures, the 3' splice site is already docked into the active site, the ensuing 3' exon sequences are anchored on PRP8, and the step II factor FAM192A contacts the duplex between U2 snRNA and the branch site. In the transition of pre-C∗-I to pre-C∗-II, the step II factors Cactin, FAM32A, PRKRIP1, and SLU7 are recruited. Notably, the RNA helicase PRP22 is positioned quite differently in the pre-C∗-I, pre-C∗-II, and C∗ complexes, suggesting a role in 3' exon binding and proofreading. Together with information on human C and C∗ complexes, our studies recapitulate a molecular choreography of the C-to-C∗ transition, revealing mechanistic insights into exon ligation.

Translating recent advances in the pathogenesis of acute myeloid leukemia to the clinic.

Despite FDA approval of nine new drugs for patients with acute myeloid leukemia (AML) in the United States over the last 4 years, AML remains a major area of unmet medical need among hematologic malignancies. In this review, we discuss the development of promising new molecular targeted approaches for AML, including menin inhibition, novel IDH1/2 inhibitors, and preclinical means to target TET2, ASXL1, and RNA splicing factor mutations. In addition, we review progress in immune targeting of AML through anti-CD47, anti-SIRPα, and anti-TIM-3 antibodies; bispecific and trispecific antibodies; and new cellular therapies in development for AML.

Poison Exon Splicing Regulates a Coordinated Network of SR Protein Expression during Differentiation and Tumorigenesis.

The RNA isoform repertoire is regulated by splicing factor (SF) expression, and alterations in SF levels are associated with disease. SFs contain ultraconserved poison exon (PE) sequences that exhibit greater identity across species than nearby coding exons, but their physiological role and molecular regulation is incompletely understood. We show that PEs in serine-arginine-rich (SR) proteins, a family of 14 essential SFs, are differentially spliced during induced pluripotent stem cell (iPSC) differentiation and in tumors versus normal tissues. We uncover an extensive cross-regulatory network of SR proteins controlling their expression via alternative splicing coupled to nonsense-mediated decay. We define sequences that regulate PE inclusion and protein expression of the oncogenic SF TRA2ß using an RNA-targeting CRISPR screen. We demonstrate location dependency of RS domain activity on regulation of TRA2ß-PE using CRISPR artificial SFs. Finally, we develop splice-switching antisense oligonucleotides to reverse the increased skipping of TRA2ß-PE detected in breast tumors, altering breast cancer cell viability, proliferation, and migration.

Splicing-specific transcriptome-wide association uncovers genetic mechanisms for schizophrenia.

Recent studies have highlighted the essential role of RNA splicing, a key mechanism of alternative RNA processing, in establishing connections between genetic variations and disease. Genetic loci influencing RNA splicing variations show considerable influence on complex traits, possibly surpassing those affecting total gene expression. Dysregulated RNA splicing has emerged as a major potential contributor to neurological and psychiatric disorders, likely due to the exceptionally high prevalence of alternatively spliced genes in the human brain. Nevertheless, establishing direct associations between genetically altered splicing and complex traits has remained an enduring challenge. We introduce Spliced-Transcriptome-Wide Associations (SpliTWAS) to integrate alternative splicing information with genome-wide association studies to pinpoint genes linked to traits through exon splicing events. We applied SpliTWAS to two schizophrenia (SCZ) RNA-sequencing datasets, BrainGVEX and CommonMind, revealing 137 and 88 trait-associated exons (in 84 and 67 genes), respectively. Enriched biological functions in the associated gene sets converged on neuronal function and development, immune cell activation, and cellular transport, which are highly relevant to SCZ. SpliTWAS variants impacted RNA-binding protein binding sites, revealing potential disruption of RNA-protein interactions affecting splicing. We extended the probabilistic fine-mapping method FOCUS to the exon level, identifying 36 genes and 48 exons as putatively causal for SCZ. We highlight VPS45 and APOPT1, where splicing of specific exons was associated with disease risk, eluding detection by conventional gene expression analysis. Collectively, this study supports the substantial role of alternative splicing in shaping the genetic basis of SCZ, providing a valuable approach for future investigations in this area.

Spliceosome component PHD finger 5A is essential for early B lymphopoiesis.

The spliceosome, a multi-megadalton ribonucleoprotein complex, is essential for pre-mRNA splicing in the nucleus and ensuring genomic stability. Its precise and dynamic assembly is pivotal for its function. Spliceosome malfunctions can lead to developmental abnormalities and potentially contribute to tumorigenesis. The specific role of the spliceosome in B cell development is poorly understood. Here, we reveal that the spliceosomal U2 snRNP component PHD finger protein 5A (Phf5a) is vital for early B cell development. Loss of Phf5a results in pronounced defects in B cell development, causing an arrest at the transition from pre-pro-B to early pro-B cell stage in the bone marrow of mutant mice. Phf5a-deficient B cells exhibit impaired immunoglobulin heavy (IgH) chain expression due to defective V-to-DJ gene rearrangement. Mechanistically, our findings suggest that Phf5a facilitates IgH gene rearrangement by regulating the activity of recombination-activating gene endonuclease and influencing chromatin interactions at the Igh locus.

Cell-type-specific dysregulation of RNA alternative splicing in short tandem repeat mouse knockin models of myotonic dystrophy.

Short tandem repeats (STRs) are prone to expansion mutations that cause multiple hereditary neurological and neuromuscular diseases. To study pathomechanisms using mouse models that recapitulate the tissue specificity and developmental timing of an STR expansion gene, we used rolling circle amplification and CRISPR/Cas9-mediated genome editing to generate Dmpk CTG expansion (CTGexp) knockin models of myotonic dystrophy type 1 (DM1). We demonstrate that skeletal muscle myoblasts and brain choroid plexus epithelial cells are particularly susceptible to Dmpk CTGexp mutations and RNA missplicing. Our results implicate dysregulation of muscle regeneration and cerebrospinal fluid homeostasis as early pathogenic events in DM1.

Exon-independent recruitment of SRSF1 is mediated by U1 snRNP stem-loop 3.

SRSF1 protein and U1 snRNPs are closely connected splicing factors. They both stimulate exon inclusion, SRSF1 by binding to exonic splicing enhancer sequences (ESEs) and U1 snRNPs by binding to the downstream 5' splice site (SS), and both factors affect 5' SS selection. The binding of U1 snRNPs initiates spliceosome assembly, but SR proteins such as SRSF1 can in some cases substitute for it. The mechanistic basis of this relationship is poorly understood. We show here by single-molecule methods that a single molecule of SRSF1 can be recruited by a U1 snRNP. This reaction is independent of exon sequences and separate from the U1-independent process of binding to an ESE. Structural analysis and cross-linking data show that SRSF1 contacts U1 snRNA stem-loop 3, which is required for splicing. We suggest that the recruitment of SRSF1 to a U1 snRNP at a 5'SS is the basis for exon definition by U1 snRNP and might be one of the principal functions of U1 snRNPs in the core reactions of splicing in mammals.

Using the ACMG/AMP framework to capture evidence related to predicted and observed impact on splicing: Recommendations from the ClinGen SVI Splicing Subgroup.

The American College of Medical Genetics and Genomics (ACMG)/Association for Molecular Pathology (AMP) framework for classifying variants uses six evidence categories related to the splicing potential of variants: PVS1, PS3, PP3, BS3, BP4, and BP7. However, the lack of guidance on how to apply such codes has contributed to variation in the specifications developed by different Clinical Genome Resource (ClinGen) Variant Curation Expert Panels. The ClinGen Sequence Variant Interpretation Splicing Subgroup was established to refine recommendations for applying ACMG/AMP codes relating to splicing data and computational predictions. We utilized empirically derived splicing evidence to (1) determine the evidence weighting of splicing-related data and appropriate criteria code selection for general use, (2) outline a process for integrating splicing-related considerations when developing a gene-specific PVS1 decision tree, and (3) exemplify methodology to calibrate splice prediction tools. We propose repurposing the PVS1_Strength code to capture splicing assay data that provide experimental evidence for variants resulting in RNA transcript(s) with loss of function. Conversely, BP7 may be used to capture RNA results demonstrating no splicing impact for intronic and synonymous variants. We propose that the PS3/BS3 codes are applied only for well-established assays that measure functional impact not directly captured by RNA-splicing assays. We recommend the application of PS1 based on similarity of predicted RNA-splicing effects for a variant under assessment in comparison with a known pathogenic variant. The recommendations and approaches for consideration and evaluation of RNA-assay evidence described aim to help standardize variant pathogenicity classification processes when interpreting splicing-based evidence.

MERISTEM-DEFECTIVE regulates the balance between stemness and differentiation in the root meristem through RNA splicing control.

Plants respond to environmental stresses through controlled stem cell maintenance and meristem activity. One level of gene regulation is RNA alternative splicing. However, the mechanistic link between stress, meristem function and RNA splicing is poorly understood. The MERISTEM-DEFECTIVE (MDF) Arabidopsis gene encodes an SR-related family protein, required for meristem function and leaf vascularization, and is the likely orthologue of the human SART1 and yeast Snu66 splicing factors. MDF is required for the correct splicing and expression of key transcripts associated with root meristem function. We identified RSZ33 and ACC1, both known to regulate cell patterning, as splicing targets required for MDF function in the meristem. MDF expression is modulated by osmotic and cold stress, associated with differential splicing and specific isoform accumulation and shuttling between nucleus and cytosol, and acts in part via a splicing target SR34. We propose a model in which MDF controls splicing in the root meristem to promote stemness and to repress stress response, cell differentiation and cell death pathways.

Self-supervised learning on millions of primary RNA sequences from 72 vertebrates improves sequence-based RNA splicing prediction.

Language models pretrained by self-supervised learning (SSL) have been widely utilized to study protein sequences, while few models were developed for genomic sequences and were limited to single species. Due to the lack of genomes from different species, these models cannot effectively leverage evolutionary information. In this study, we have developed SpliceBERT, a language model pretrained on primary ribonucleic acids (RNA) sequences from 72 vertebrates by masked language modeling, and applied it to sequence-based modeling of RNA splicing. Pretraining SpliceBERT on diverse species enables effective identification of evolutionarily conserved elements. Meanwhile, the learned hidden states and attention weights can characterize the biological properties of splice sites. As a result, SpliceBERT was shown effective on several downstream tasks: zero-shot prediction of variant effects on splicing, prediction of branchpoints in humans, and cross-species prediction of splice sites. Our study highlighted the importance of pretraining genomic language models on a diverse range of species and suggested that SSL is a promising approach to enhance our understanding of the regulatory logic underlying genomic sequences.

IRIS: Discovery of cancer immunotherapy targets arising from pre-mRNA alternative splicing.

Alternative splicing (AS) is prevalent in cancer, generating an extensive but largely unexplored repertoire of novel immunotherapy targets. We describe Isoform peptides from RNA splicing for Immunotherapy target Screening (IRIS), a computational platform capable of discovering AS-derived tumor antigens (TAs) for T cell receptor (TCR) and chimeric antigen receptor T cell (CAR-T) therapies. IRIS leverages large-scale tumor and normal transcriptome data and incorporates multiple screening approaches to discover AS-derived TAs with tumor-associated or tumor-specific expression. In a proof-of-concept analysis integrating transcriptomics and immunopeptidomics data, we showed that hundreds of IRIS-predicted TCR targets are presented by human leukocyte antigen (HLA) molecules. We applied IRIS to RNA-seq data of neuroendocrine prostate cancer (NEPC). From 2,939 NEPC-associated AS events, IRIS predicted 1,651 epitopes from 808 events as potential TCR targets for two common HLA types (A*02:01 and A*03:01). A more stringent screening test prioritized 48 epitopes from 20 events with "neoantigen-like" NEPC-specific expression. Predicted epitopes are often encoded by microexons of ≤30 nucleotides. To validate the immunogenicity and T cell recognition of IRIS-predicted TCR epitopes, we performed in vitro T cell priming in combination with single-cell TCR sequencing. Seven TCRs transduced into human peripheral blood mononuclear cells (PBMCs) showed high activity against individual IRIS-predicted epitopes, providing strong evidence of isolated TCRs reactive to AS-derived peptides. One selected TCR showed efficient cytotoxicity against target cells expressing the target peptide. Our study illustrates the contribution of AS to the TA repertoire of cancer cells and demonstrates the utility of IRIS for discovering AS-derived TAs and expanding cancer immunotherapies.