Systematic Profiling of Alternative mRNA Splicing Signature for Predicting Glioblastoma Prognosis

mRNA isoform balance in neuronal development and disease

Balanced mRNA isoform diversity and abundance are spatially and temporally regulated throughout cellular differentiation. The proportion of expressed isoforms contributes to cell type specification and determines key properties of the differentiated cells. Neurons are unique cell types with intricate developmental programs, characteristic cellular morphologies, and electrophysiological potential. Neuron-specific gene expression programs establish these distinctive cellular characteristics and drive diversity among neuronal subtypes. Genes with neuron-specific alternative processing are enriched in key neuronal functions, including synaptic proteins, adhesion molecules, and scaffold proteins. Despite the similarity of neuronal gene expression programs, each neuronal subclass can be distinguished by unique alternative mRNA processing events. Alternative processing of developmentally important transcripts alters coding and regulatory information, including interaction domains, transcript stability, subcellular localization, and targeting by RNA binding proteins. Fine-tuning of mRNA processing is essential for neuronal activity and maintenance. Thus, the focus of neuronal RNA biology research is to dissect the transcriptomic mechanisms that underlie neuronal homeostasis, and consequently, predispose neuronal subtypes to disease. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA in Disease and Development > RNA in Development.

The inhibitory effect of human umbilical cord mesenchymal stem cells expressing anti-HAAH scFv-sTRAIL fusion protein on glioma

Glioma is the most common malignant intracranial tumor with low 5-year survival rate. In this study, we constructed a plasmid expressing anti-HAAH single-chain antibody and sTRAIL fusion protein (scFv-sTRAIL), and explored the effects of the double gene modified human umbilical cord mesenchyreal stem cells (hucMSCs) on the growth of glioma in vitro and in vivo. The isolated hucMSCs were identified by detecting the adipogenic differentiation ability and the osteogenic differentiation ability. The phenotypes of hucMSCs were determined by the flow cytometry. The hucMSCs were infected with lentivirus expression scFv-sTRAIL fusion protein. The expression of sTRAIL in hucMSCs were detected by immunofluorescence staining, western blot and ELISA. The tropism of hucMSCs toward U87G cells was assessed by transwell assay. The inhibitory effect of hucMSCs on U87G cells were explored by CCK8 and apoptosis assay. The xenograft tumor was established by subcutaneously injection of U87G cells into the back of mice. The hucMSCs were injected via tail veins. The inhibitory effect of hucMSCs on glioma in vivo was assessed by TUNEL assay. The hucMSCs migrated into the xenograft tumor were revealed by detecting the green fluorescent. The results showed that the scFv-sTRAIL expression did not affect the phenotypes of hucMSCs. The scFv-sTRAIL expression promoted the tropism of hucMSCs toward U87G cells, enhanced the inhibitory effect and tumor killing effect of hucMSCs on U87G cells. The in vivo study showed that hucMSCs expressing scFv-sTRAIL demonstrated significantly higher inhibitory effect and tumor killing effect than hucMSCs expressing sTRAIL. The green fluorescence intensity in the mice injected with hucMSCs expressing scFv-sTRAIL was significantly higher than that injected with hucMSCs expressing sTRAIL. These data suggested that the scFv conferred the targeting effect of hucMSCs tropism towards the xenograft tumor. In conclusion, the hucMSCs expressing scFv-sTRAIL fusion protein gained the capability to target and kill gliomas cells in vitro and in vivo. These findings shed light on a potential therapy for glioma treatment.

Splicing is an alternate oncogenic pathway activation mechanism in glioma

High-grade diffuse glioma (HGG) is the leading cause of brain tumour death. While the genetic drivers of HGG have been well described, targeting these has thus far had little impact on survival suggesting other mechanisms are at play. Here we interrogate the alternative splicing landscape of pediatric and adult HGG through multi-omic analyses, uncovering an increased splicing burden compared with normal brain. The rate of recurrent alternative splicing in cancer drivers exceeds their mutation rate, a pattern that is recapitulated in pan-cancer analyses, and is associated with worse prognosis in HGG. We investigate potential oncogenicity by interrogating cancer pathways affected by alternative splicing in HGG; spliced cancer drivers include members of the RAS/MAPK pathway. RAS suppressor neurofibromin 1 is differentially spliced to a less active isoform in >80% of HGG downstream from REST upregulation, activating the RAS/MAPK pathway and reducing glioblastoma patient survival. Overall, our results identify non-mutagenic mechanisms by which cancers activate oncogenic pathways which need to accounted for in personalized medicine approaches.

Targeting genetic drivers of high grade diffuse glioma (HGG) has not improved patient survival, suggesting the involvement of other mechanisms. Here, across cancer types, the authors identify increased alternative splicing burden in cancer drivers compared to mutation rate as an alternative mechanism for activation of oncogenic pathways such as RAS/MAPK.

Systematic analysis of the relationship between ovarian cancer prognosis and alternative splicing

Background

Ovarian cancer(OC) is the gynecological tumor with the highest mortality rate, effective biomarkers are of great significance in improving its prognosis. In recent years, there have been many studies on alternative splicing (AS) events, and the role of AS events in tumor has become a focus of attention.

Methods

Data were downloaded from the TCGA database and Univariate Cox regression analysis was performed to determine AS events associated with OC prognosis.Eight prognostic models of OC were constructed in R package, and the accuracy of the models were evaluated by the time-dependent receiver operating characteristic (ROC) curves.Eight types of survival curves were drawn to evaluate the differences between the high and low risk groups.Independent prognostic factors of OC were analyzed by single factor independent analysis and multi-factor independent prognostic analysis.Again, Univariate Cox regression analysis was used to analyze the relationship between splicing factors(SF) and AS events, and Gene Ontology(GO) and Kyoto Encyclopedia of Genes and Genomes(KEGG) enrichment analysis were performed on OS-related SFs to understand the pathways.

Results

Univariate Cox regression analysis showed that among the 15,278 genes, there were 31,286 overall survival (OS) related AS events, among which 1524 AS events were significantly correlated with OS. The area under the time-dependent receiver operating characteristic curve (AUC) of AT and ME were the largest and the RI was the smallest,which were 0.757 and 0.68 respectively. The constructed models have good value for the prognosis assessment of OC patients. Among the eight survival curves, AP was the most significant difference between the high and low risk groups, with a P value of 1.61e − 1.The results of single factor independent analysis and multi-factor independent prognostic analysis showed that risk score calculated by the model and age could be used as independent risk factors.According to univariate COX regression analysis,109 SFs were correlated with AS events and adjusted in two ways: positive and negative.

Conclusions

SFs and AS events can directly or indirectly affect the prognosis of OC patients. It is very important to find effective prognostic markers to improve the survival rate of OC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13048-021-00866-1.

Elucidation of the BMI1 interactome identifies novel regulatory roles in glioblastoma

Glioblastoma (GBM) is the most common and aggressive intrinsic brain tumour in adults. Epigenetic mechanisms controlling normal brain development are often dysregulated in GBM. Among these, BMI1, a structural component of the Polycomb Repressive Complex 1 (PRC1), which promotes the H2AK119ub catalytic activity of Ring1B, is upregulated in GBM and its tumorigenic role has been shown in vitro and in vivo. Here, we have used protein and chromatin immunoprecipitation followed by mass spectrometry (MS) analysis to elucidate the protein composition of PRC1 in GBM and transcriptional silencing of defining interactors in primary patient-derived GIC lines to assess their functional impact on GBM biology. We identify novel regulatory functions in mRNA splicing and cholesterol transport which could represent novel targetable mechanisms in GBM.

A cell-to-patient machine learning transfer approach uncovers novel basal-like breast cancer prognostic markers amongst alternative splice variants

BACKGROUND

Breast cancer is amongst the 10 first causes of death in women worldwide. Around 20% of patients are misdiagnosed leading to early metastasis, resistance to treatment and relapse. Many clinical and gene expression profiles have been successfully used to classify breast tumours into 5 major types with different prognosis and sensitivity to specific treatments. Unfortunately, these profiles have failed to subclassify breast tumours into more subtypes to improve diagnostics and survival rate. Alternative splicing is emerging as a new source of highly specific biomarkers to classify tumours in different grades. Taking advantage of extensive public transcriptomics datasets in breast cancer cell lines (CCLE) and breast cancer tumours (TCGA), we have addressed the capacity of alternative splice variants to subclassify highly aggressive breast cancers.

RESULTS

Transcriptomics analysis of alternative splicing events between luminal, basal A and basal B breast cancer cell lines identified a unique splicing signature for a subtype of tumours, the basal B, whose classification is not in use in the clinic yet. Basal B cell lines, in contrast with luminal and basal A, are highly metastatic and express epithelial-to-mesenchymal (EMT) markers, which are hallmarks of cell invasion and resistance to drugs. By developing a semi-supervised machine learning approach, we transferred the molecular knowledge gained from these cell lines into patients to subclassify basal-like triple negative tumours into basal A- and basal B-like categories. Changes in splicing of 25 alternative exons, intimately related to EMT and cell invasion such as ENAH, CD44 and CTNND1, were sufficient to identify the basal-like patients with the worst prognosis. Moreover, patients expressing this basal B-specific splicing signature also expressed newly identified biomarkers of metastasis-initiating cells, like CD36, supporting a more invasive phenotype for this basal B-like breast cancer subtype.

CONCLUSIONS

Using a novel machine learning approach, we have identified an EMT-related splicing signature capable of subclassifying the most aggressive type of breast cancer, which are basal-like triple negative tumours. This proof-of-concept demonstrates that the biological knowledge acquired from cell lines can be transferred to patients data for further clinical investigation. More studies, particularly in 3D culture and organoids, will increase the accuracy of this transfer of knowledge, which will open new perspectives into the development of novel therapeutic strategies and the further identification of specific biomarkers for drug resistance and cancer relapse.

Transcriptome analysis and prognostic model construction based on splicing profiling in glioblastoma

Glioblastoma (GBM) is the most aggressive malignant brain tumour, with high morbidity and mortality rates. Currently, there is a lack of systematic and comprehensive analysis on the prognostic significance of alternative splicing (AS) profiling for GBM. The GBM data, including RNA-sequencing, corresponding clinical information and the expression levels of splicing factor genes, were downloaded from The Cancer Genome Atlas and the SpliceAid2 database. The prognostic models were assessed by the least absolute shrinkage and selection operator Cox regression analysis. The correlation network between survival-associated AS events and splicing factors was plotted. Prognostic models were built for every AS event type and performed well for risk stratification in patients with GBM. The final prognostic signature served as an independent prognostic factor [hazard ratio (HR), 4.61; 95% confidence interval (CI), 2.97-7.16; P=9.66×10-12] for several clinical parameters, including age, sex, isocitrate dehydrogenase mutation, O6-methylguanine-DNA methyltransferase promoter methylation and risk score. The HR for risk score with GBM was 1.0063 (95% CI, 1.0024-1.0103). The splicing regulatory network indicated that heat shock protein b-1, protein arginine N-methyltransferase 5, protein FAM50B and endoplasmic reticulum chaperone BiP genes were independent prognostic factors for GBM. The results of the present study support the ongoing effort in developing novel genomic models and providing potentially more effective treatment options for patients with GBM.

Feasibility of phosphoproteomics to uncover oncogenic signalling in secreted extracellular vesicles using glioblastoma-EGFRVIII cells as a model

Cancer cells secrete extracellular vesicles (EVs) that contain molecular information, including proteins and RNA. Oncogenic signalling can be transferred via the cargo of EVs to recipient cells and may influence the behaviour of neighbouring cells or cells at a distance. This cargo may contain cancer drivers, such as EGFR, and also phosphorylated (activated) components of oncogenic signalling cascades. Till date, the cancer EV phosphoproteome has not been studied in great detail. In the present study, we used U87 and U87EGFRvIII cells as a model to explore EV oncogenic signalling components in comparison to the cellular profile. EVs were isolated using the VN96 ME-kit and subjected to LC-MS/MS based phosphoproteomics and dedicated bioinformatics. Expression of (phosphorylated)-EGFR was highly increased in EGFRvIII overexpressing cells and their secreted EVs. The increased phosphorylated proteins in both cells and EVs were associated with activated components of the EGFR-signalling cascade and included EGFR, AKT2, MAPK8, SMG1, MAP3K7, DYRK1A, RPS6KA3 and PAK4 kinases. In conclusion, EVs harbour oncogenic signalling networks including multiple activated kinases including EGFR, AKT and mTOR. SIGNIFICANCE: Extracellular vesicles (EVs) are biomarker treasure troves and are widely studied for their biomarker content in cancer. However, little research has been done on the phosphorylated protein profile within cancer EVs. In the current study, we demonstrate that EVs that are secreted by U87-EGFRvIII mutant glioblastoma cells contain high levels of oncogenic signalling networks. These networks contain multiple activated (phosphorylated) kinases, including EGFR, MAPK, AKT and mTOR.

SRRM4 Expands the Repertoire of Circular RNAs by Regulating Microexon Inclusion

High-throughput RNA sequencing (RNA-seq) and dedicated bioinformatics pipelines have synergized to identify an expansive repertoire of unique circular RNAs (circRNAs), exceeding 100,000 variants. While the vast majority of these circRNAs comprise canonical exonic and intronic sequences, microexons (MEs)-which occur in 30% of functional mRNA transcripts-have been entirely overlooked. CircRNAs which contain these known MEs (ME-circRNAs) could be identified with commonly utilized circRNA prediction pipelines, CIRCexplorer2 and CIRI2, but were not previously recognized as ME-circRNAs. In addition, when employing a bespoke bioinformatics pipeline for identifying RNA chimeras, called Hyb, we could also identify over 2000 ME-circRNAs which contain novel MEs at their backsplice junctions, that are uncalled by either CIRCexplorer2 or CIRI2. Analysis of circRNA-seq datasets from gliomas of varying clinical grades compared with matched control tissue has shown circRNAs have potential as prognostic markers for stratifying tumor from healthy tissue. Furthermore, the abundance of microexon-containing circRNAs (ME-circRNAs) between tumor and normal tissues is correlated with the expression of a splicing associated factor, Serine/arginine repetitive matrix 4 (SRRM4). Overexpressing SRRM4, known for regulating ME inclusion in mRNAs critical for neural differentiation, in human HEK293 cells resulted in the biogenesis of over 2000 novel ME-circRNAs, including ME-circEIF4G3, and changes in the abundance of many canonical circRNAs, including circSETDB2 and circLBRA. This shows SRRM4, in which its expression is correlated with poor prognosis in gliomas, acts as a bona fide circRNA biogenesis factor. Given the known roles of MEs and circRNAs in oncogenesis, the identification of these previously unrecognized ME-circRNAs further increases the complexity and functional purview of this non-coding RNA family.

Clusterization in acute myeloid leukemia based on prognostic alternative splicing signature to reveal the clinical characteristics in the bone marrow microenvironment

Background

Alternative splicing (AS), a crucial post-transcriptional regulatory mechanism in expanding the coding capacities of genomes and increasing the diversity of proteins, still faces various challenges in the splicing regulation mechanism of acute myeloid leukemia (AML) and microenvironmental changes.

Results

A total of 27,833 AS events were detected in 8337 genes in 178 AML patients, with exon skip being the predominant type. Approximately 11% of the AS events were significantly related to prognosis, and the prediction models based on various events demonstrated high classification efficiencies. Splicing factors correlation networks further altered the diversity of AS events through epigenetic regulation and clarified the potential mechanism of the splicing pathway. Unsupervised cluster analysis revealed significant correlations between AS and immune features, molecular mutations, immune checkpoints and clinical outcome. The results suggested that AS clusters could be used to identify patient subgroups with different survival outcomes in AML, among which C1 was both associated with good outcome in overall survival. Interestingly, C1 was associated with lower immune scores compared with C2 and C3, and favorable-risk cytogenetics was rarely distributed in C2, but much more common in C1.

Conclusions

This study revealed a comprehensive landscape of AS events, and provides new insight into molecular targeted therapy and immunotherapy strategy for AML.

Diverse molecular functions of aspartate β‑hydroxylase in cancer (Review)

Aspartate/asparagine β-hydroxylase (AspH) is a type II transmembrane protein that catalyzes the post-translational hydroxylation of definite aspartyl and asparaginyl residues in epidermal growth factor-like domains of substrates. In the last few decades, accumulating evidence has indicated that AspH expression is upregulated in numerous types of human malignant cancer and is associated with poor survival and prognosis. The AspH protein aggregates on the surface of tumor cells, which contributes to inducing tumor cell migration, infiltration and metastasis. However, small-molecule inhibitors targeting hydroxylase activity can markedly block these processes, both in vitro and in vivo. Immunization of tumor-bearing mice with a phage vaccine fused with the AspH protein can substantially delay tumor growth and progression. Additionally, AspH antigen-specific CD4⁺ and CD8⁺ T cells were identified in the spleen of tumor-bearing mice. Therefore, these agents may be used as novel strategies for cancer treatment. The present review summarizes the current progress on the underlying mechanisms of AspH expression in cancer development.

Transcriptome-wide analysis and modelling of prognostic alternative splicing signatures in invasive breast cancer: a prospective clinical study

Aberrant alternative splicing (AS) has been highly involved in the tumorigenesis and progression of most cancers. The potential role of AS in invasive breast cancer (IBC) remains largely unknown. In this study, RNA sequencing of IBC samples from The Cancer Genome Atlas was acquired. AS events were screened by conducting univariate and multivariate Cox analysis and least absolute shrinkage and selection operator regression. In total, 2146 survival-related AS events were identified from 1551 parental genes, of which 93 were related to prognosis, and a prognostic marker model containing 14 AS events was constructed. We also constructed the regulatory network of splicing factors (SFs) and AS events, and identified DDX39B as the node SF gene, and verified the accuracy of the network through experiments. Next, we performed quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) in triple negative breast cancer patients with different responses to neoadjuvant chemotherapy, and found that the exon-specific expression of EPHX2, C6orf141, and HERC4 was associated with the different status of patients that received neoadjuvant chemotherapy. In conclusion, this study found that DDX39B, EPHX2 (exo7), and HERC4 (exo23) can be used as potential targets for the treatment of breast cancer, which provides a new idea for the treatment of breast cancer.

LY6K promotes glioblastoma tumorigenicity via CAV-1-mediated ERK1/2 signaling enhancement

BACKGROUND

Lymphocyte antigen 6 complex, locus K (LY6K) is a putative oncogene in various cancers. Elevated expression of LY6K is correlated with poor patient prognosis in glioblastoma (GBM). The aim of this study is to advance our understanding of the mechanism by which LY6K contributes to GBM tumor biology.

METHODS

Bioinformatic data mining was used to investigate LY6K expression in relation to GBM clinical outcome. To understand the role of LY6K in GBM, we utilized patient-derived glioma stemlike cells (GSCs) and U87 cells and employed immunoblotting, immunofluorescent staining, radiation treatment, and orthotopic GBM xenograft models.

RESULTS

Our results show that increased expression of LY6K inversely correlates with GBM patient survival. LY6K promotes tumorigenicity in GBM cells both in vitro and in vivo. The mechanism underlying this tumorigenic behavior is enhancement of extracellular signal-regulated kinase 1 and 2 (ERK1/2) signaling. Interestingly, we observed that tumor-promoting LY6K-ERK1/2 signaling is mediated by the interaction of LY6K with caveolin-1, rather than through oncogenic receptor tyrosine kinase-mediated signaling. Moreover, association of LY6K with the cell membrane is crucial for its tumorigenic functions. Finally, DNA methylation maintains LY6K silencing, and hypomethylation of the LY6K promoter increases its expression. In GSCs, ionizing radiation leads to demethylation of the LY6K promoter, thereby increasing LY6K expression and GSC resistance to radiation.

CONCLUSIONS

Our study highlights the importance of the contribution of LY6K to GBM tumor biology and suggests LY6K as a potential membrane target for treating GBM.

Aspartate β-hydroxylase as a target for cancer therapy

As metastasis is a major cause of death in cancer patients, new anti-metastatic strategies are needed to improve cancer therapy outcomes. Numerous pathways have been shown to contribute to migration and invasion of malignant tumors. Aspartate β-hydroxylase (ASPH) is a key player in the malignant transformation of solid tumors by enhancing cell proliferation, migration, and invasion. ASPH also promotes tumor growth by stimulation of angiogenesis and immunosuppression. These effects are mainly achieved via the activation of Notch and SRC signaling pathways. ASPH expression is upregulated by growth factors and hypoxia in different human tumors and its inactivation may have broad clinical impact. Therefore, small molecule inhibitors of ASPH enzymatic activity have been developed and their anti-metastatic effect confirmed in preclinical mouse models. ASPH can also be targeted by monoclonal antibodies and has also been used as a tumor-associated antigen to induce both cluster of differentiation (CD) 8+ and CD4+ T cells in mice. The PAN-301-1 vaccine against ASPH has already been tested in a phase 1 clinical trial in patients with prostate cancer. In summary, ASPH is a promising target for anti-tumor and anti-metastatic therapy based on inactivation of catalytic activity and/or immunotherapy.

Systematic Profiling of Alternative Splicing for Sarcoma Patients Reveals Novel Prognostic Biomarkers Associated with Tumor Microenvironment and Immune Cells

Background

Alternative splicing (AS) events is a novel biomarker of tumor prognosis, but the role of AS events in sarcoma patients remains unclear.

Material/Methods

RNA-seq and clinicopathologic data of the sarcoma cohort were extracted from the TCGA database and data on AS events were downloaded from the TCGASpliceSeq database. Univariate Cox analysis, LASSO regression analysis, and multivariate Cox analysis were performed to determine the overall survival (OS)- and disease-free survival (DFS)-related AS events. Two nomograms were developed based on the independent variables, and subgroup analysis was performed. The area under the curve (AUC), calibration curve, and decision curve analysis (DCA) were used to evaluate the nomograms. Then, we used the CIBERSORT and ESTIMATE package to determine the immune cell proportion and tumor microenvironment (TME) score, respectively. The associations between AS events-based clusters and TME and immune cells were studied.

Results

We identified 1945 and 1831 AS events as OS- and DFS-related AS events, respectively. Two nomograms based on the AS events and clinical data were established and the AUCs of nomograms ranged from 0.807 to 0.894. The calibration curve and DCA showed excellent performance of nomograms. In addition, the results indicated the distinct relationships between AS events-based clusters and OS, DFS, immune score, stromal score, and 10 immune cells.

Conclusions

Our study indicated that AS events are novel prognostic biomarkers for sarcoma patients that may be associated with the TME and immune cells.

Malignant Evaluation and Clinical Prognostic Values of m6A RNA Methylation Regulators in Glioblastoma

N6-methyladenosine (m6A) RNA methylation, the most common form of mRNA modification and regulated by the m6A RNA methylation regulators ("writers," "erasers," and "readers"), has been reported to be associated with the progression of the malignant tumor. However, its role in glioblastoma (GBM) has been poorly known. This study aimed to identify the expression, potential functions, and prognostic values of m6A RNA methylation regulators in GBM. Here, we revealed that the 13 central m6A RNA methylation regulators were firmly related to the clinical and molecular phenotype of GBM. Taking advantage of consensus cluster analysis, we obtained two categories of GBM samples and found malignancy-related processes of m6A methylation regulators and compounds that specifically targeted the malignant processes. Besides, we also obtained a list of genes with poor prognosis in GBM. Finally, we derived a risk-gene signature with three selected m6A RNA methylation regulators, which allowed us to extend the in-depth study and dichotomized the OS of patients with GBM into high- and low-risk subgroups. Notably, this risk-gene signature could be used as independent prognostic markers and accurate clinicopathological parameter predictors. In conclusion, m6A RNA methylation regulators are a type of vital participant in the malignant progression of GBM, with a critical potential in the prognostic stratification and treatment strategies of GBM.