Aberrant splicing in neuroblastoma generates RNA-fusion transcripts and provides vulnerability to spliceosome inhibitors

Spliceosomal vulnerability of MYCN-amplified neuroblastoma is contingent on PRMT5-mediated regulation of epitranscriptomic and metabolomic pathways

Approximately 50% of poor prognosis neuroblastomas arise due to MYCN over-expression. We previously demonstrated that MYCN and PRMT5 proteins interact and PRMT5 knockdown led to apoptosis of MYCN amplified (MNA) neuroblastoma. Here we evaluate the highly selective first-in-class PRMT5 inhibitor GSK3203591 and its in vivo analogue GSK3326593 as targeted therapeutics for MNA neuroblastoma. Cell-line analyses show MYCN-dependent growth inhibition and apoptosis, with approximately 200-fold greater sensitivity of MNA neuroblastoma lines. RNA sequencing of three MNA neuroblastoma lines treated with GSK3203591 reveal deregulated MYCN transcriptional programmes and altered mRNA splicing, converging on key regulatory pathways such as DNA damage response, epitranscriptomics and cellular metabolism. Stable isotope labelling experiments in the same cell lines demonstrate that glutamine metabolism is impeded following GSK3203591 treatment, linking with disruption of the MLX/Mondo nutrient sensors via intron retention of MLX mRNA. Interestingly, glutaminase (GLS) protein decreases after GSK3203591 treatment despite unchanged transcript levels. We demonstrate that the RNA methyltransferase METTL3 and cognate reader YTHDF3 proteins are lowered following their mRNAs undergoing GSK3203591-induced splicing alterations, indicating epitranscriptomic regulation of GLS; accordingly, we observe decreases of GLS mRNA m6A methylation following GSK3203591 treatment, and decreased GLS protein following YTHDF3 knockdown. In vivo efficacy of GSK3326593 is confirmed by increased survival of Th-MYCN mice, with drug treatment triggering splicing events and protein decreases consistent with in vitro data. Together our study demonstrates the PRMT5-dependent spliceosomal vulnerability of MNA neuroblastoma and identifies the epitranscriptome and glutamine metabolism as critical determinants of this sensitivity.

Identification of gene fusions associated with amyotrophic lateral sclerosis

INTRODUCTION/AIMS

Genetics is an important risk factor for amyotrophic lateral sclerosis (ALS), a neurodegenerative disease affecting motor neurons. Recent findings demonstrate that in addition to specific genetic mutations, structural variants caused by genetic instability can also play a causative role in ALS. Genomic instability can lead to deletions, duplications, insertions, inversions, and translocations in the genome, and these changes can sometimes lead to fusion of distinct genes into a single transcript. Gene fusion events have been studied extensively in cancer; however, they have not been thoroughly investigated in ALS. The aim of this study was to determine whether gene fusions are present in ALS.

METHODS

Gene fusions were identified using STAR Fusion v1.10.0 software in bulk RNA-Seq data from human postmortem samples from publicly available data sets from Target ALS and the New York Genome Center ALS Consortium.

RESULTS

We report the presence of gene fusion events in several brain regions as well as in spinal cord samples in ALS. Although most gene fusions were intra-chromosomal events between neighboring genes and present in both ALS and control samples, there was a significantly greater number of unique gene fusions in ALS compared to controls. Lastly, we identified specific gene fusions with a significant burden in ALS, that were absent from both control samples and known cancer gene fusion databases.

DISCUSSION

Collectively, our findings reveal an enrichment of gene fusions in ALS and suggest that these events may be an additional genetic cause linked to ALS pathogenesis.

RBM39 degrader invigorates natural killer cells to eradicate neuroblastoma despite cancer cell plasticity

The cellular plasticity of neuroblastoma is defined by a mixture of two major cell states, adrenergic (ADRN) and mesenchymal (MES), which may contribute to therapy resistance. However, how neuroblastoma cells switch cellular states during therapy remains largely unknown and how to eradicate neuroblastoma regardless of their cell states is a clinical challenge. To better understand the lineage switch of neuroblastoma in chemoresistance, we comprehensively defined the transcriptomic and epigenetic map of ADRN and MES types of neuroblastomas using human and murine models treated with indisulam, a selective RBM39 degrader. We showed that cancer cells not only undergo a bidirectional switch between ADRN and MES states, but also acquire additional cellular states, reminiscent of the developmental pliancy of neural crest cells. The lineage alterations are coupled with epigenetic reprogramming and dependency switch of lineage-specific transcription factors, epigenetic modifiers and targetable kinases. Through targeting RNA splicing, indisulam induces an inflammatory tumor microenvironment and enhances anticancer activity of natural killer cells. The combination of indisulam with anti-GD2 immunotherapy results in a durable, complete response in high-risk transgenic neuroblastoma models, providing an innovative, rational therapeutic approach to eradicate tumor cells regardless of their potential to switch cell states.

Metabolic reprogramming of cancer cells by JMJD6-mediated pre-mRNA splicing associated with therapeutic response to splicing inhibitor

Dysregulated pre-mRNA splicing and metabolism are two hallmarks of MYC-driven cancers. Pharmacological inhibition of both processes has been extensively investigated as potential therapeutic avenues in preclinical and clinical studies. However, how pre-mRNA splicing and metabolism are orchestrated in response to oncogenic stress and therapies is poorly understood. Here, we demonstrate that jumonji domain containing 6, arginine demethylase, and lysine hydroxylase, JMJD6, acts as a hub connecting splicing and metabolism in MYC-driven human neuroblastoma. JMJD6 cooperates with MYC in cellular transformation of murine neural crest cells by physically interacting with RNA binding proteins involved in pre-mRNA splicing and protein homeostasis. Notably, JMJD6 controls the alternative splicing of two isoforms of glutaminase (GLS), namely kidney-type glutaminase (KGA) and glutaminase C (GAC), which are rate-limiting enzymes of glutaminolysis in the central carbon metabolism in neuroblastoma. Further, we show that JMJD6 is correlated with the anti-cancer activity of indisulam, a 'molecular glue' that degrades splicing factor RBM39, which complexes with JMJD6. The indisulam-mediated cancer cell killing is at least partly dependent on the glutamine-related metabolic pathway mediated by JMJD6. Our findings reveal a cancer-promoting metabolic program is associated with alternative pre-mRNA splicing through JMJD6, providing a rationale to target JMJD6 as a therapeutic avenue for treating MYC-driven cancers.

Metabolic reprogramming of cancer cells by JMJD6-mediated pre-mRNA splicing is associated with therapeutic response to splicing inhibitor

Dysregulated pre-mRNA splicing and metabolism are two hallmarks of MYC-driven cancers. Pharmacological inhibition of both processes has been extensively investigated as potential therapeutic avenues in preclinical and clinical studies. However, how pre-mRNA splicing and metabolism are orchestrated in response to oncogenic stress and therapies is poorly understood. Here, we demonstrate that Jumonji Domain Containing 6, Arginine Demethylase and Lysine Hydroxylase, JMJD6, acts as a hub connecting splicing and metabolism in MYC-driven neuroblastoma. JMJD6 cooperates with MYC in cellular transformation by physically interacting with RNA binding proteins involved in pre-mRNA splicing and protein homeostasis. Notably, JMJD6 controls the alternative splicing of two isoforms of glutaminase (GLS), namely kidney-type glutaminase (KGA) and glutaminase C (GAC), which are rate-limiting enzymes of glutaminolysis in the central carbon metabolism in neuroblastoma. Further, we show that JMJD6 is correlated with the anti-cancer activity of indisulam, a "molecular glue" that degrades splicing factor RBM39, which complexes with JMJD6. The indisulam-mediated cancer cell killing is at least partly dependent on the glutamine-related metabolic pathway mediated by JMJD6. Our findings reveal a cancer-promoting metabolic program is associated with alternative pre-mRNA splicing through JMJD6, providing a rationale to target JMJD6 as a therapeutic avenue for treating MYC-driven cancers.

Aurora kinase A regulates cancer-associated RNA aberrant splicing in breast cancer

The contribution of oncogenes to tumor-associated RNA splicing and the relevant molecular mechanisms therein require further elaboration. Here, we show that oncogenic Aurora kinase A (AURKA) promotes breast cancer-related RNA aberrant splicing in a context-dependent manner. AURKA regulated pan-breast cancer-associated RNA splicing events including GOLGA4, RBM4 and UBQLN1. Aberrant splicing of GOLGA4 and RBM4 was closely related to breast cancer development. Mechanistically, AURKA interacted with the splicing factor YBX1 and promoted AURKA-YBX1 complex-mediated GOLGA4 exon inclusion. AURKA binding to the splicing factor hnRNPK promoted AURKA-hnRNPK complex-mediated RBM4 exon skipping. Analysis of clinical data identified an association between the AURKA-YBX1/hnRNPK complex and poor prognosis in breast cancer. Blocking AURKA nuclear translocation with small molecule drugs partially reversed the oncogenic splicing of RBM4 and GOLGA4 in breast cancer cells. In summary, oncogenic AURKA executes its function on modulating breast cancer-related RNA splicing, and nuclear AURKA is distinguished as a hopeful target in the case of treating breast cancer.

Investigating SH-SY5Y Neuroblastoma Cell Surfaceome as a Model for Neuronal-Targeted Novel Therapeutic Modalities

The SH-SY5Y neuroblastoma cells are a widely used in vitro model approximating neurons for testing the target engagement of therapeutics designed for neurodegenerative diseases and pain disorders. However, their potential as a model for receptor-mediated delivery and uptake of novel modalities, such as antibody-drug conjugates, remains understudied. Investigation of the SH-SY5Y cell surfaceome will aid in greater in vitro to in vivo correlation of delivery and uptake, thereby accelerating drug discovery. So far, the majority of studies have focused on total cell proteomics from undifferentiated and differentiated SH-SY5Y cells. While some studies have investigated the expression of specific proteins in neuroblastoma tissue, a global approach for comparison of neuroblastoma cell surfaceome to the brain and dorsal root ganglion (DRG) neurons remains uninvestigated. Furthermore, an isoform-specific evaluation of cell surface proteins expressed on neuroblastoma cells remains unexplored. In this study, we define a bioinformatic workflow for the identification of high-confidence surface proteins expressed on brain and DRG neurons using tissue proteomic and transcriptomic data. We then delineate the SH-SY5Y cell surfaceome by surface proteomics and show that it significantly overlaps with the human brain and DRG neuronal surface proteome. We find that, for 32% of common surface proteins, SH-SY5Y-specific major isoforms are alternatively spliced, maintaining their protein-coding ability, and are predicted to localize to the cell surface. Validation of these isoforms using surface proteomics confirms a SH-SY5Y-specific alternative NRCAM (neuron-glia related cell adhesion molecule) isoform, which is absent in typical brain neurons, but present in neuroblastomas, making it a receptor of interest for neuroblastoma-specific therapeutics.

A signature based on five immune-related genes to predict the survival and immune characteristics of neuroblastoma

BACKGROUND

MYCN amplification (MNA) has been proved to be related to poor prognosis in neuroblastoma (NBL), but the MYCN-related immune signatures and genes remain unclear.

METHODS

Enrichment analysis was used to identify the significant enrichment pathways of differentially expressed immune-related genes (DEIRGs). Weight gene coexpression network analysis (WGCNA) was applied to reveal the correlation between these DEIRGs and MYCN status. Univariate and multivariate Cox analyses were used to construct risk model. The relevant fractions of immune cells were evaluated by CIBERSORT and single-sample gene set enrichment analysis (ssGSEA).

RESULTS

Five genes, including CHGA, PTGER1, SHC3, PLXNC1, and TRIM55 were enrolled into the risk model. Kaplan-Meier survival analysis and receiver operating characteristic (ROC) curve showed that our model performed well in predicting the outcomes of NBL (3-years AUC = 0.720, 5-year AUC = 0.775, 10-years AUC = 0.782), which has been validated in the GSE49711 dataset and the E-MTAB-8248 dataset. By comparing with the tumor immune dysfunction and exclusion (TIDE) and tumor inflammation signature (TIS), we further proved that our model is reliable. Univariate and multivariate Cox regression analyses indicated that the risk score, age, and MYCN can serve as independent prognostic factors in the E-MATB-8248. Functional enrichment analysis showed the DEIRGs were enriched in leukocyte adhesion-related signaling pathways. Gene set enrichment analysis (GSEA) revealed the significantly enriched pathways of the five MYCN-related DEIRGs. The risk score was negatively correlated with the immune checkpoint CD274 (PD-L1) but no significant difference with the TMB. We also confirmed the prognostic value of our model in predicting immunotherapeutics.

CONCLUSION

We constructed and verified a signature based on DEIRG that related to MNA and predicted the survival of NBL based on relevant immune signatures. These findings could provide help for predicting prognosis and developing immunotherapy in NBL.

Non-Canonical Splicing and Its Implications in Brain Physiology and Cancer

The advance of experimental and computational techniques has allowed us to highlight the existence of numerous different mechanisms of RNA maturation, which have been so far unknown. Besides canonical splicing, consisting of the removal of introns from pre-mRNA molecules, non-canonical splicing events may occur to further increase the regulatory and coding potential of the human genome. Among these, splicing of microexons, recursive splicing and biogenesis of circular and chimeric RNAs through back-splicing and trans-splicing processes, respectively, all contribute to expanding the repertoire of RNA transcripts with newly acquired regulatory functions. Interestingly, these non-canonical splicing events seem to occur more frequently in the central nervous system, affecting neuronal development and differentiation programs with important implications on brain physiology. Coherently, dysregulation of non-canonical RNA processing events is associated with brain disorders, including brain tumours. Herein, we summarize the current knowledge on molecular and regulatory mechanisms underlying canonical and non-canonical splicing events with particular emphasis on cis-acting elements and trans-acting factors that all together orchestrate splicing catalysis reactions and decisions. Lastly, we review the impact of non-canonical splicing on brain physiology and pathology and how unconventional splicing mechanisms may be targeted or exploited for novel therapeutic strategies in cancer.

Targeting the spliceosome through RBM39 degradation results in exceptional responses in high-risk neuroblastoma models

Aberrant alternative pre-mRNA splicing plays a critical role in MYC-driven cancers and therefore may represent a therapeutic vulnerability. Here, we show that neuroblastoma, a MYC-driven cancer characterized by splicing dysregulation and spliceosomal dependency, requires the splicing factor RBM39 for survival. Indisulam, a “molecular glue” that selectively recruits RBM39 to the CRL4-DCAF15 E3 ubiquitin ligase for proteasomal degradation, is highly efficacious against neuroblastoma, leading to significant responses in multiple high-risk disease models, without overt toxicity. Genetic depletion or indisulam-mediated degradation of RBM39 induces significant genome-wide splicing anomalies and cell death. Mechanistically, the dependency on RBM39 and high-level expression of DCAF15 determine the exquisite sensitivity of neuroblastoma to indisulam. Our data indicate that targeting the dysregulated spliceosome by precisely inhibiting RBM39, a vulnerability in neuroblastoma, is a valid therapeutic strategy.

Ral GEF with the PH Domain and SH3 Binding Motif 1 Regulated by Splicing Factor Junction Plakoglobin and Pyrimidine Metabolism Are Prognostic in Uterine Carcinosarcoma

Uterine carcinosarcoma (UCS) is a highly invasive malignant tumor that originated from the uterine epithelium. Many studies suggested that the abnormal changes of alternative splicing (AS) of pre-mRNA are related to the occurrence and metastasis of the tumor. This study investigates the mechanism of alternative splicing events (ASEs) in the tumorigenesis and metastasis of UCS. RNA-seq of UCS samples and alternative splicing event (ASE) data of UCS samples were downloaded from The Cancer Genome Atlas (TCGA) and TCGASpliceSeq databases, several times. Firstly, we performed the Cox regression analysis to identify the overall survival-related alternative splicing events (OSRASEs). Secondly, a multivariate model was applied to approach the prognostic values of the risk score. Afterwards, a coexpressed network between splicing factors (SFs) and OSRASEs was constructed. In order to explore the relationship between the potential prognostic signaling pathways and OSRASEs, we fabricated a network between these pathways and OSRASEs. Finally, validations from multidimension platforms were used to explain the results unambiguously. 1,040 OSRASEs were identified by Cox regression. Then, 6 OSRASEs were incorporated in a multivariable model by Lasso regression. The area under the curve (AUC) of the receiver operator characteristic (ROC) curve was 0.957. The risk score rendered from the multivariate model was corroborated to be an independent prognostic factor (P < 0.001). In the network of SFs and ASEs, junction plakoglobin (JUP) noteworthily regulated RALGPS1-87608-AT (P < 0.001, R = 0.455). Additionally, RALGPS1-87608-AT (P = 0.006) showed a prominent relationship with distant metastasis. KEGG pathways related to prognosis of UCS were selected by gene set variation analysis (GSVA). The pyrimidine metabolism (P < 0.001, R = -0.470) was the key pathway coexpressed with RALGPS1. We considered that aberrant JUP significantly regulated RALGPS1-87608-AT and the pyrimidine metabolism pathway might play a significant part in the metastasis and prognosis of UCS.