Identification of differential splicing genes in gliomas using exon expression profiling

Role of long noncoding RNAs in the regulation of alternative splicing in glioblastoma

Glioblastoma multiforme (GBM) is a highly severe primary brain tumor. Despite extensive research, effective treatments remain elusive. Long noncoding RNAs (lncRNAs) play a significant role in both cancer and normal biology. They influence alternative splicing (AS), which is crucial in cancer. Advances in lncRNA-specific microarrays and next-generation sequencing have enhanced understanding of AS. Abnormal AS contributes to cancer invasion, metastasis, apoptosis, therapeutic resistance, and tumor development, including glioma. lncRNA-mediated AS affects several cellular signaling pathways, promoting or suppressing cancer malignancy. This review discusses the lncRNAs regulating AS in glioblastoma and their mechanisms.

ATRT-SHH comprises three molecular subgroups with characteristic clinical and histopathological features and prognostic significance

Atypical teratoid/rhabdoid tumor (ATRT) is an aggressive central nervous system tumor characterized by loss of SMARCB1/INI1 protein expression and comprises three distinct molecular groups, ATRT–TYR, ATRT–MYC and ATRT–SHH. ATRT–SHH represents the largest molecular group and is heterogeneous with regard to age, tumor location and epigenetic profile. We, therefore, aimed to investigate if heterogeneity within ATRT–SHH might also have biological and clinical importance. Consensus clustering of DNA methylation profiles and confirmatory t-SNE analysis of 65 ATRT–SHH yielded three robust molecular subgroups, i.e., SHH-1A, SHH-1B and SHH-2. These subgroups differed by median age of onset (SHH-1A: 18 months, SHH-1B: 107 months, SHH-2: 13 months) and tumor location (SHH-1A: 88% supratentorial; SHH-1B: 85% supratentorial; SHH-2: 93% infratentorial, often extending to the pineal region). Subgroups showed comparable SMARCB1 mutational profiles, but pathogenic/likely pathogenic SMARCB1 germline variants were over-represented in SHH-2 (63%) as compared to SHH-1A (20%) and SHH-1B (0%). Protein expression of proneural marker ASCL1 (enriched in SHH-1B) and glial markers OLIG2 and GFAP (absent in SHH-2) as well as global mRNA expression patterns differed, but all subgroups were characterized by overexpression of SHH as well as Notch pathway members. In a Drosophila model, knockdown of Snr1 (the fly homologue of SMARCB1) in hedgehog activated cells not only altered hedgehog signaling, but also caused aberrant Notch signaling and formation of tumor-like structures. Finally, on survival analysis, molecular subgroup and age of onset (but not ASCL1 staining status) were independently associated with overall survival, older patients (> 3 years) harboring SHH-1B experiencing relatively favorable outcome. In conclusion, ATRT–SHH comprises three subgroups characterized by SHH and Notch pathway activation, but divergent molecular and clinical features. Our data suggest that molecular subgrouping of ATRT–SHH has prognostic relevance and might aid to stratify patients within future clinical trials.

OKN-007 Alters Protein Expression Profiles in High-Grade Gliomas: Mass Spectral Analysis of Blood Sera

Current therapies for high-grade gliomas, particularly glioblastomas (GBM), do not extend patient survival beyond 16–22 months. OKN-007 (OKlahoma Nitrone 007), which is currently in phase II (multi-institutional) clinical trials for GBM patients, and has demonstrated efficacy in several rodent and human xenograft glioma models, shows some promise as an anti-glioma therapeutic, as it affects most aspects of tumorigenesis (tumor cell proliferation, angiogenesis, migration, and apoptosis). Combined with the chemotherapeutic agent temozolomide (TMZ), OKN-007 is even more effective by affecting chemo-resistant tumor cells. In this study, mass spectrometry (MS) methodology ESI-MS, mass peak analysis (Leave One Out Cross Validation (LOOCV) and tandem MS peptide sequence analyses), and bioinformatics analyses (Ingenuity^® Pathway Analysis (IPA^®)), were used to identify up- or down-regulated proteins in the blood sera of F98 glioma-bearing rats, that were either untreated or treated with OKN-007. Proteins of interest identified by tandem MS-MS that were decreased in sera from tumor-bearing rats that were either OKN-007-treated or untreated included ABCA2, ATP5B, CNTN2, ITGA3, KMT2D, MYCBP2, NOTCH3, and VCAN. Conversely, proteins of interest in tumor-bearing rats that were elevated following OKN-007 treatment included ABCA6, ADAMTS18, VWA8, MACF1, and LAMA5. These findings, in general, support our previous gene analysis, indicating that OKN-007 may be effective against the ECM. These findings also surmise that OKN-007 may be more effective against oligodendrogliomas, other brain tumors such as medulloblastoma, and possibly other types of cancers.

Genome-wide profiling of alternative splicing in glioblastoma and their clinical value

Background

Alternative splicing (AS), one of the main post-transcriptional biological regulation mechanisms, plays a key role in the progression of glioblastoma (GBM). Systematic AS profiling in GBM is limited and urgently needed.

Methods

TCGA SpliceSeq data and the corresponding clinical data were downloaded from the TCGA data portal. Survival-related AS events were identified through Kaplan–Meier survival analysis and univariate Cox analysis. Then, splicing correlation network was constructed based on these AS events and associated splicing factors. LASSO regression followed by multivariate Cox analysis was performed to validate independent AS biomarkers and to construct a risk prediction model. Enrichment analysis was subsequently conducted to explore potential signaling pathways of these AS events.

Results

A total of 132 TCGA GBM samples and 45,610 AS events were included in our study, among which 416 survival-related AS events were identified. An AS correlation network, including 54 AS events and 94 splicing factors, was constructed, and further functional enrichment was performed. Moreover, the novel risk prediction model we constructed displayed moderate performance (the area under the curves were > 0.7) at both one, two and three years.

Conclusions

Survival-related AS events may be vital factors of both biological function and prognosis. Our findings in this study can deepen the understanding of the complicated mechanisms of AS in GBM and provide novel insights for further study. Moreover, our risk prediction model is ready for preliminary clinical applications. Further verification is required.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08681-z.

Alternative RNA Splicing-The Trojan Horse of Cancer Cells in Chemotherapy

Almost all transcribed human genes undergo alternative RNA splicing, which increases the diversity of the coding and non-coding cellular landscape. The resultant gene products might have distinctly different and, in some cases, even opposite functions. Therefore, the abnormal regulation of alternative splicing plays a crucial role in malignant transformation, development, and progression, a fact supported by the distinct splicing profiles identified in both healthy and tumor cells. Drug resistance, resulting in treatment failure, still remains a major challenge for current cancer therapy. Furthermore, tumor cells often take advantage of aberrant RNA splicing to overcome the toxicity of the administered chemotherapeutic agents. Thus, deciphering the alternative RNA splicing variants in tumor cells would provide opportunities for designing novel therapeutics combating cancer more efficiently. In the present review, we provide a comprehensive outline of the recent findings in alternative splicing in the most common neoplasms, including lung, breast, prostate, head and neck, glioma, colon, and blood malignancies. Molecular mechanisms developed by cancer cells to promote oncogenesis as well as to evade anticancer drug treatment and the subsequent chemotherapy failure are also discussed. Taken together, these findings offer novel opportunities for future studies and the development of targeted therapy for cancer-specific splicing variants.

Overlapping variants in the blood, tissues and cell lines for patients with intracranial meningiomas are predominant in stem cell-related genes

OBJECTIVE

Bulk tissue genomic analysis of meningiomas identified common somatic mutations, however, it often excluded blood-related variants. In contrast, genomic characterisation of primary cell lines that can provide critical information regarding growth and proliferation, have been rare. In our work, we identified the variants that are present in the blood, tissues and corresponding cell lines that are likely to be predictive, tumorigenic and progressive.

METHOD

Whole-exome sequencing was used to identify variants and distinguish related pathways that exist in 42 blood, tissues and corresponding cell lines (BTCs) samples for patients with intracranial meningiomas. Conventional sequencing was used for the confirmation of variants. Integrative analysis of the gene expression for the corresponding samples was utilised for further interpretations.

RESULTS

In total, 926 BTC variants were detected, implicating 845 genes. A pathway analysis of all BTC genes with damaging variants indicated the 'cell morphogenesis involved in differentiation' stem cell-related pathway to be the most frequently affected pathway. Concordantly, five stem cell-related genes, GPRIN2, ALDH3B2, ASPN, THSD7A and SIGLEC6, showed BTC variants in at least five of the patients. Variants that were heterozygous in the blood and homozygous in the tissues or the corresponding cell lines were rare (average: 1.3 ± 0.3%), and included variants in the RUNX2 and CCDC114 genes. An analysis comparing the variants detected only in tumours with aggressive features indicated a total of 240 BTC genes, implicating the 'homophilic cell adhesion via plasma membrane adhesion molecules' pathway, and identifying the stem cell-related transcription coactivator NCOA3/AIB1/SRC3 as the most frequent BTC gene. Further analysis of the possible impact of the poly-Q mutation present in the NCOA3 gene indicated associated deregulation of 15 genes, including the up-regulation of the stem cell related SEMA3D gene and the angiogenesis related VEGFA gene.

CONCLUSION

Stem cell-related pathways and genes showed high prevalence in the BTC variants, and novel variants in stem cell-related genes were identified for meningioma. These variants can potentially be used as predictive, tumorigenic and progressive biomarkers for meningioma.

Prognostic alternative splicing signature reveals the landscape of immune infiltration in Pancreatic Cancer

Background: Pancreatic cancer (PC) is an aggressive cancer with worse survival in the world. Emerging evidence suggested that the imbalance of alternative splicing (AS) is a hallmark of cancer and indicated poor prognosis of patients. Genes-derived splicing events can produce neoepitopes for immunotherapy. However, the profound study of splicing profiling in PC is still elusive. We aimed to identification of novel prognostic signature across a comprehensive splicing landscape and reveal their relationship with tumor-infiltrating immune cells in pancreatic cancer microenvironment.

Methods: Based on integrated analysis of splicing profiling and clinical data, differentially splicing events were filtered out. Then, stepwise Cox regression analysis was applied to identify survival-related splicing events and construct prognostic signature. Functional enrichment analysis was performed to explore biology function. Kaplan-Meier curves and receiver operating characteristic (ROC) curves were performed to validate the predictive effect of predictive signature. We also verified the clinical value of prognostic signature under the influence of different clinical parameters. For deeper analysis, we evaluated the correlation between prognostic signature and infiltrating immune cells by CIBERSORT.

Results: According to systematic analyzing, a final six splicing events were identified and validated the good prognostic capability in entire TCGA dataset, validation set 1 and validation set 2 by Kaplan-Meier curves (P < 0.0001). The area under the curve (AUC) of ROC curves were also confirmed the high predictive efficiency of the prognostic signature in these three cohorts (AUC = 0.857, 0.895 and 0.788). In order to validate whether prognostic signature highlights a correlation between AS and immune contexture, CIBERSORT was performed to analyze the proportion of tumor-infiltrating immune cells in PC. Based on prognostic signature, we identified survival-related immune cells including CD8 T cells (P = 0.0111), activated CD4 memory T cells (P = 0.0329) and resting mast cells (P = 0.0352).

Conclusion: In conclusion, our study contribute to provide a promising prognostic signature based on six splicing events and revealed prognosis-related immune cells which indeed represented novel tumor drivers and provide potential targets for personalized therapeutic.

hsa_circ_0001947 suppresses acute myeloid leukemia progression via targeting hsa-miR-329-5p/CREBRF axis

Aim: Accumulating evidence has indicated that circular RNAs (circRNAs) are involved in cancer biology. However, their roles in acute myeloid leukemia (AML) remain unclear. Therefore, we aimed to define novel circRNAs involved the development and progression of AML. Materials & methods: We used circRNAs microarray to determine the differential expression profile. Quantitative reverse transcription PCR analyzed the expression of hsa_circ_0001947. The siRNA assesses the function of hsa_circ_0001947 in vitro and in vivo. A dual-luciferase and mimics/inhibitor were to determine the target gene relationship. Results: hsa_circ_0001947 functions as a tumor inhibitor to suppress AML cell proliferation through hsa-miR-329-5p/ CREBRF axis. Conclusion: hsa_circ_0001947 may be as a novel potential biomarker for the treatment of AML.

RNA processing genes characterize RNA splicing and further stratify lower-grade glioma

BACKGROUND

Aberrant expression of RNA processing genes may drive the alterative RNA profile in lower-grade gliomas (LGGs). Thus, we aimed to further stratify LGGs based on the expression of RNA processing genes.

METHODS

This study included 446 LGGs from The Cancer Genome Atlas (TCGA, training set) and 171 LGGs from the Chinese Glioma Genome Atlas (CGGA, validation set). The least absolute shrinkage and selection operator (LASSO) Cox regression algorithm was conducted to develop a risk-signature. The receiver operating characteristic (ROC) curves and Kaplan-Meier curves were used to study the prognosis value of the risk-signature.

RESULTS

Among the tested 784 RNA processing genes, 276 were significantly correlated with the OS of LGGs. Further LASSO Cox regression identified a 19-gene risk-signature, whose risk score was also an independently prognosis factor (P<0.0001, multiplex Cox regression) in the validation dataset. The signature had better prognostic value than the traditional factors "age", "grade" and "WHO 2016 classification" for 3- and 5-year survival both two datasets (AUCs > 85%). Importantly, the risk-signature could further stratify the survival of LGGs in specific subgroups of WHO 2016 classification. Furthermore, alternative splicing events for genes such as EGFR and FGFR were found to be associated with the risk score. mRNA expression levels for genes, which participated in cell proliferation and other processes, were significantly correlated to the risk score.

CONCLUSIONS

Our results highlight the role of RNA processing genes for further stratifying the survival of patients with LGGs and provide insight into the alternative splicing events underlying this role.

Identification of a Novel Splice Variant of Neural Cell Adhesion Molecule in Glioblastoma Through Proteogenomics Analysis

Splice variants are known to be important in the pathophysiology of tumors, including the brain cancers. We applied a proteogenomics pipeline to identify splice variants in glioblastoma (GBM, grade IV glioma), a highly malignant brain tumor, using in-house generated mass spectrometric proteomic data and public domain RNASeq dataset. Our analysis led to the identification of a novel exon that maps to the long isoform of Neural cell adhesion molecule 1 (NCAM1), expressed on the surface of glial cells and neurons, important for cell adhesion and cell signaling. The presence of the novel exon is supported with the identification of five peptides spanning it. Additional peptides were also detected in sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) gel separated proteins from GBM patient tissue, underscoring the presence of the novel peptides in the intact brain protein. The novel exon was detected in the RNASeq dataset in 18 of 25 GBM samples and separately validated in additional 10 GBM tumor tissues using quantitative real-time-polymerase chain reaction (qRT-PCR). Both transcriptomic and proteomic data indicate downregulation of NCAM1, including the novel variant, in GBM. Domain analysis of the novel NCAM1 sequence indicates that the insertion of the novel exon contributes extra low-complexity region in the protein that may be important for protein-protein interactions and hence for cell signaling associated with tumor development. Taken together, the novel NCAM1 variant reported in this study exemplifies the importance of future multiomics research and systems biology applications in GBM.

Current Challenges and Opportunities in Treating Glioblastoma

Glioblastoma multiforme (GBM), the most common and aggressive primary brain tumor, has a high mortality rate despite extensive efforts to develop new treatments. GBM exhibits both intra- and intertumor heterogeneity, lending to resistance and eventual tumor recurrence. Large-scale genomic and proteomic analysis of GBM tumors has uncovered potential drug targets. Effective and “druggable” targets must be validated to embark on a robust medicinal chemistry campaign culminating in the discovery of clinical candidates. Here, we review recent developments in GBM drug discovery and delivery. To identify GBM drug targets, we performed extensive bioinformatics analysis using data from The Cancer Genome Atlas project. We discovered 20 genes, BOC, CLEC4GP1, ELOVL6, EREG, ESR2, FDCSP, FURIN, FUT8-AS1, GZMB, IRX3, LITAF, NDEL1, NKX3-1, PODNL1, PTPRN, QSOX1, SEMA4F, TH, VEGFC, and C20orf166AS1 that are overexpressed in a subpopulation of GBM patients and correlate with poor survival outcomes. Importantly, nine of these genes exhibit higher expression in GBM versus low-grade glioma and may be involved in disease progression. In this review, we discuss these proteins in the context of GBM disease progression. We also conducted computational multi-parameter optimization to assess the blood-brain barrier (BBB) permeability of small molecules in clinical trials for GBM treatment. Drug delivery in the context of GBM is particularly challenging because the BBB hinders small molecule transport. Therefore, we discuss novel drug delivery methods, including nanoparticles and prodrugs. Given the aggressive nature of GBM and the complexity of targeting the central nervous system, effective treatment options are a major unmet medical need. Identification and validation of biomarkers and drug targets associated with GBM disease progression present an exciting opportunity to improve treatment of this devastating disease.

Computer Analysis of Glioma Transcriptome Profiling: Alternative Splicing Events

Here we present the analysis of alternative splicing events on an example of glioblastoma cell culture samples using a set of computer tools in combination with database integration. The gene expression profiles of glioblastoma were obtained from cell culture samples of primary glioblastoma which were isolated and processed for RNA extraction. Transcriptome profiling of normal brain samples and glioblastoma were done by Illumina sequencing. The significant differentially expressed exon-level probes and their corresponding genes were identified using a combination of the splicing index method. Previous studies indicated that tumor-specific alternative splicing is important in the regulation of gene expression and corresponding protein functions during cancer development. Multiple alternative splicing transcripts have been identified as progression markers, including generalized splicing abnormalities and tumor- and stage-specific events. We used a set of computer tools which were recently applied to analysis of gene expression in laboratory animals to study differential splicing events. We found 69 transcripts that are differentially alternatively spliced. Three cancer-associated genes were considered in detail, in particular: APP (amyloid beta precursor protein), CASC4 (cancer susceptibility candidate 4) and TP53. Such alternative splicing opens new perspectives for cancer research.

RNA processing as an alternative route to attack glioblastoma

Genomic analyses have become an important tool to identify new avenues for therapy. This is especially true for cancer types with extremely poor outcomes, since our lack of effective therapies offers no tangible clinical starting point to build upon. The highly malignant brain tumor glioblastoma (GBM) exemplifies such a refractory cancer, with only 15 month average patient survival. Analyses of several hundred GBM samples compiled by the TCGA (The Cancer Genome Atlas) have produced an extensive transcriptomic map, identified prevalent chromosomal alterations, and defined important driver mutations. Unfortunately, clinical trials based on these results have not yet delivered an improvement on outcome. It is, therefore, necessary to characterize other regulatory routes known for playing a role in tumor relapse and response to treatment. Alternative splicing affects more than 90% of the human coding genes and it is an important source for transcript variation and gene regulation. Mutations and alterations in splicing factors are highly prevalent in multiple cancers, demonstrating the potential for splicing to act as a tumor driver. As a result, numerous genes are expressed as cancer-specific splicing isoforms that are functionally distinct from the canonical isoforms found in normal tissue. These include genes that regulate cancer-critical pathways such as apoptosis, DNA repair, cell proliferation, and migration. Splicing defects can even induce genomic instability, a common characteristic of cancer, and a driver of tumor evolution. Importantly, components of the splicing machinery are targetable; multiple drugs can inhibit splicing factors or promote changes in splicing which could be exploited to begin improving clinical outcomes. Here, we review the current literature and present a case for exploring RNA processing as therapeutic route for the treatment of GBM.

Contactin-2/TAG-1, active on the front line for three decades

Contactin-2/transiently expressed axonal surface glycoprotein-1 (TAG-1) is a cell adhesion molecule belonging to the immunoglobulin superfamily (IgSF). It has six immunoglobulin-like extracellular domains and four fibronectin III-like ones, with anchoring to the cell membrane through glycosylphosphatidyl inositol. Contactin-2/TAG-1 is expressed in specific neurons transiently on the axonal surface during the fetal period. In postnatal stages, Contactin-2/TAG-1 is expressed in cerebellar granule cells, hippocampal pyramidal cells, and the juxtaparanodal regions of myelinated nerve fibers. In the embryonic nervous system, Contactin-2/TAG-1 plays important roles in axonal elongation, axonal guidance, and cellular migration. In the postnatal nervous system, it also plays an essential role in the formation of myelinated nerve fibers. Moreover, Contactin-2/TAG-1 has been linked to autoimmune diseases of the human nervous system. Taken together, Contactin-2/TAG-1 plays a central role in a variety of functions from development to disease.

RACK1 affects glioma cell growth and differentiation through the CNTN2-mediated RTK/Ras/MAPK pathway

Receptor for activated C kinase 1 (RACK1) and contactin-2 (CNTN2) are known to be abnormally expressed in gliomas; however, the association between RACK1 and CNTN2, and the effects of RACK1 and CNTN2 on glioma cell differentiation and the related molecular mechanisms remain largely unknown. The present study aimed to investigate the interaction between RACK1 and CNTN2, and to examine whether RACK1/CNTN2/receptor tyrosine kinase (RTK)/Ras/mitogen-activated protein kinase (MAPK) axis plays a role in glioma growth and differentiation. The results from western blot analysis revealed that the protein expression levels of RACK1 and CNTN2 were higher in high‑grade glioma tissues and cells, and lower in low-grade glioma tissues and cells. A co-immunoprecipitation assay demonstrated that RACK1 interacts with CNTN2, and RACK1 upregulated the expression of CNTN2. Gain-of‑function and loss-of‑function experiments indicated that both RACK1 and CNTN2 promoted glioma cell proliferation, inhibited glioma cell differentiation and activated the RTK/Ras/MAPK pathway. However, the effects of RACK1 on glioma cell proliferation, differentiation and the activation of the RTK/Ras/MAPK signaling pathway were abolished by the knockdown of CNTN2 using siRNA. In Therefore, the findings of this study firstly demonstrate that RACK1 interacts with CNTN2, and that the effects of RACK1 on glioma cell growth and differentiation are mediated by CNTN2. The RACK1/CNTN2/RTK/Ras/MAPK axis exists in glioma cells, and it may be a potential therapeutic target in gliomas.