Landscape of Alternative Splicing Events Related to Prognosis and Immune Infiltration in Glioma: A Data Analysis and Basic Verification

Su, Hong-xin; Yang, Gang; Su, Fei; Hu, Chen-xiao; Zhang, Tao; Ran, Jun-tao; Guan, Quan-lin

doi:https://doi.org/10.1155/2022/2671891

Journal of Immunology Research

On this page

Abstract Introduction Methods Results Discussion Data Availability Consent Conflicts of Interest Authors’ Contributions References Copyright Related Articles

Special Issue

Interaction between Innate Immunity and Adaptive Immunity in Human Diseases

View this Special Issue

Research Article | Open Access

Volume 2022 | Article ID 2671891 | https://doi.org/10.1155/2022/2671891

Landscape of Alternative Splicing Events Related to Prognosis and Immune Infiltration in Glioma: A Data Analysis and Basic Verification

Hong-xin Su,¹Gang Yang,²Fei Su,³Chen-xiao Hu,⁴Tao Zhang,³Jun-tao Ran,¹and Quan-lin Guan⁵

Academic Editor: Dawei Cui

Received16 Feb 2022

Revised04 May 2022

Accepted12 May 2022

Published04 Jul 2022

Abstract

Background. Glioma is a prevalent primary brain cancer with high invasiveness and typical local diffuse infiltration. Alternative splicing (AS), as a pervasive transcriptional regulatory mechanism, amplifies the coding capacity of the genome and promotes the progression of malignancies. This study was aimed at identifying AS events and novel biomarkers associated with survival for glioma. Methods. RNA splicing patterns were collected from The Cancer Genome Atlas SpliceSeq database, followed by calculating the percentage of splicing index. Expression profiles and related clinical information of glioma were integrated based on the UCSC Xena database. The AS events in glioma were further analyzed, and glioma prognosis-related splicing factors were identified with the use of bioinformatics analysis and laboratory techniques. Further immune infiltration analysis was performed. Results. Altogether, 9028 AS events were discovered. Upon univariate Cox analysis, 425 AS events were found to be related to the survival of patients with glioma, and 42 AS events were further screened to construct the final prognostic model (). Additionally, decreased expression of the splicing factors including Neuro-Oncological Ventral Antigen 1 (NOVA1), heterogeneous nuclear ribonucleoprotein C (HNRNPC), heterogeneous nuclear ribonucleoprotein L-like protein (HNRNPLL), and RNA-Binding Motif Protein 4 (RBM4) contributed to the poor survival in glioma. The immune infiltration analysis demonstrated that AS events were related to the proportion of immune cells infiltrating in glioma. Conclusions. It is of great value for comprehensive consideration of AS events, splicing networks, and related molecular subtype clusters in revealing the underlying mechanism and immune microenvironment remodeling for glioma, which provides clues for the further verification of related therapeutic targets.

1. Introduction

Genomic analysis has become an important tool for developing new therapeutic options, which has facilitated the opening of a new era of tumor genomics research. The mRNA precursors are spliced, and the remaining exons can be reconnected in different ways at the posttranslational level. The resultant diverse mature mRNAs are then translated into different protein variants. The process is defined as alternative splicing (AS), and over 90% of human multiexon genes undergo AS [1]. Protein factors involved in the splicing process of RNA precursors are considered splicing factors (SFs). Consequently, a gene can encode multiple proteins, which contributes to proteome diversity and cell complexity [2].

AS plays a critical role in numerous vital biological processes, including tissue differentiation; its dysfunction may induce multiple diseases, such as neurodegenerative diseases [3], malignant tumors [4], bladder cancer [5], lung cancer [6], ovarian cancer [7], gastric cancer [8], colorectal carcinoma [9], breast cancer [10], oral squamous cell carcinoma [11], and lymphoma [12]. Previous studies [2] have confirmed that a number of genes have tumor-specific splice variants, which share different functions from the typical variants expressed in normal tissues. These gene variants participate in multiple key pathways related to tumorigenesis, such as apoptosis regulation, DNA repair, proliferation, and migration. Splicing defects can also lead to genomic instability, which is a common feature of tumors, as well as being the driver of tumor progression. Recently, accumulating evidence demonstrated that AS contributed to the tumor microenvironment remodeling by affecting the proportion of immune infiltrating cells and regulating the activity of immune cells [13, 14]. Thus, cancer-specific AS might serve as a potential biomarker.

Glioma is a tumor that originates from neuroglial cells in the brain. According to the standard formulated by the World Health Organization (WHO), glioma is classified as grade I, II, III, or IV [15]. Typically, low-grade gliomas (grades I and II) have low invasiveness and relatively favorable prognosis, but they usually cause damage to the functional areas of the brain, whereas high-grade gliomas (grades III and IV) are known to be aggressive, with poor survival. The average survival of patients with glioma is about 15 months, and the 5-year overall survival (OS) rate is <5% (WHO 2016). Currently, there are only a few reliable biomarkers for predicting the course or prognosis of patients with glioma [16]. Due to the limited treatment methods and poor curative effect of glioma, especially high-grade glioma, there is an increased need to find more efficient targets.

Glioma exhibits significant heterogeneity at microscopic and molecular levels. The Cancer Genome Atlas (TCGA) can provide an integrated somatic cell atlas for glioma based on molecular and clinical data, which is a valuable resource for investigating genomic disorder in glioma. RNA processing is a major contributor to transcript variation and gene expression regulation and plays a pivotal role in many diseases and cancers. Therefore, it is extremely necessary to systemically analyze AS and its regulations for tumor immune microenvironment in glioma.

2. Methods

2.1. AS Events in TCGA RNA-seq Database

The mRNA splicing profiles were extracted from TCGA SpliceSeq (https://bioinformatics.mdanderson.org/TCGASpliceSeq/) database. The percentage splicing index (PSI) (range, 0–100%) is used to quantify AS events and calculate seven AS types: Mutually Exclusive Exons (ME), Exon Skip (ES), Retained Intron (RI), Alternate Terminator (AT), Alternate Promoter (AP), Alternate Acceptor site (AA), and Alternate Donor site (AD). In this study, and were established as the thresholds to screen the splicing patterns of related protein-coding genes in patients with glioma. In addition, UCSC Xena (http://xena.ucsc.edu/) was utilized to obtain the glioma expression profile data and related clinical information, which were later integrated with AS data for further investigation.

2.2. Survival Analysis

Patients were classified into two groups based on the median threshold of every parameter. AS events were subjected to univariate and multivariate Cox regression analyses to identify their association with OS, and was used to indicate significant associations. A prognostic risk score was calculated by multiplying the AS PSI linear combination by the related regression coefficient () that indicated related weight. In this study, the regression coefficient was determined using the multivariate Cox proportional hazard regression model. The risk score formula is shown below:

Finally, the AS events that were shown to independently predict prognosis were integrated to construct a prognostic prediction model. We also plotted the survival curve and receiver operating characteristic (ROC) curve for the establishment of prognostic models. Kaplan-Meier analysis was performed for survival analysis, which was compared using the log-rank test.

2.3. Construction of a Splicing-Related Network

Altogether, 67 human SF lists were generated from manual literature and database screening. The gene expression profile data of SFs were extracted from UCSC Xena for subsequent analysis. Survival analysis () was applied for the identification of SFs related to prognosis. Spearman correlation analysis was performed to analyze the correlation between survival-associated SF gene expression and survival-associated AS PSI value, and and were set to indicate the significant difference. Moreover, Cytoscape (3.6.0 National Resource for Network Biology) was used to construct an AS event-SF interaction network.

2.4. Prognostic- and Subtype Analysis-Related Clusters

AS events varied considerably among individual samples. Thus, to obtain a reliable classification, this study utilized the unsupervised consensus method with R package Consensus Cluster Plus to identify the glioma molecular subtypes. We assessed the prognostic value of different clusters to determine relationships of subtypes with survival. To investigate the potential regulation in biological processes, cellular components, and molecular functions, the ClusterProfiler package was employed for Gene Ontology (GO) functional enrichment analysis with genes involved in survival-associated AS events. Additionally, the ReactomePA package was utilized for Reactome pathway enrichment analysis to identify the significant pathways.

2.5. Tissue Samples

To validate the expression of key SFs in glioma, 10 tumors were collected from patients diagnosed with glioma who underwent surgery at the First Hospital of Lanzhou University. Of the 10 patients with glioma, there were 6 men and 4 women. They ranged in age from 35 to 74 years old, with an average age of years. Normal brain tissues were collected from 10 patients whose brain tissues were partially removed due to head injury. The samples were obtained with informed consent from all patients. This study was approved by the Ethics Committee of the First Hospital of Lanzhou University.

2.6. Western Blotting

Total protein was extracted from frozen tissues, and protein concentration was measured. The protein sample was isolated by sodium dodecyl sulfate polyacrylamide gel electrophoresis and then transferred onto a polyvinylidene fluoride membrane. After the membranes were blocked in BSA (10%) for 2 h, the target antibody was applied for the incubation with membranes overnight at 4°C. Subsequently, TBST was used to wash the membrane three times for 10 min, followed by incubation with a goat anti-IgG secondary antibody at ambient temperature. Thereafter, the membranes treated with ECL photoluminescence solution were analyzed using ImageJ software. More details about the antibodies were shown below: the anti-Nova1 antibody (1 : 1000), anti-RBM4 antibody (1 : 800), and anti-HNRNPLL antibody (1 : 1200) were purchased from Abcam (Cambridge, MA, USA).

2.7. Immune Cell Infiltration Analysis

The TIMER database, as a platform, could provide in-depth analysis and visualization of tumor 5infiltrating immune cells for various cancers [17]. The correlation between key SF gene expression and infiltration levels of B cells, CD4+ T cells, CD8+ T cells, macrophages, neutrophils, and dendritic cells in glioma was calculated with the TIMER database. To further assess the relationship between the key SF genes and immune cell infiltration, the CBERSORT algorithm was applied, which could calculate the proportion of 22 immune cells.

3. Results

3.1. Overview of AS Events

In total, 9028 AS events were observed among 4117 genes in the 146 glioma samples based on TCGA database (Figure 1(a)). The findings indicated that a single gene was associated with multiple mRNA splicing event types, which underwent as many as five types of AS events. ES was the predominant AS event: over 1/3 AS events were ESs (Figure 1(b)).

(a)

(b)

3.2. Survival-Related AS Events in Glioma

To explore the significance of AS events in predicting the prognosis of glioma, univariate Cox regression analysis was performed. A total of 425 AS events were identified as significantly associated with survival (), including 24 AAs, 26 ADs, 188 APs, 57 ATs, 158 ESs, 1 ME, and 41 RIs (Figure 2(a)). As shown in Figure 2(b), the Circos graph displays the survival-associated AS events and the involved genes.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 2

AS events and prognosis in glioma. (a) The number of survival-related AS events and related genes in each type of differential AS. (b) Circos diagram of survival-related AS events and related genes: the Circos panel from outside to inside is expressed as follows: chromosome number, genome axis, and survival-related AS event-related genes; the number of related genes in the overall AS event (showing 1-10 different heights, more than 10 calculated as 10); correlation; the number of AS types of genes in the overall event; the values of related genes in the single-factor COX regression analysis (expressed by the conversion value of ; the higher the height, the more significant the value); the HR value of correlation genes in the univariate COX regression analysis (in which red represents and black represents ); and the correlation between genes. (c–h) Forest diagram of top 15 AS events related to survival in each AS type.

3.3. Construction of the Glioma Prognostic Models

The forest plot showed the top 15 prognostic AS events of the univariate analysis in each AS type (Figures 2(c)–2(h)) The top 15 AS events in each AS type were selected for multivariate Cox regression analysis and obtained 42 independent prognostic AS events, including AA events (), AD events (), AP events (), AT events (), ES () events, and RI () events. In total, the independent variables included 7 AA events, 8 AD events, 6 AP events, 5 AT events, 10 ES events, and 6 RI events. The prognostic AS events among the different AS types were incorporated to construct a prognosis model. Based on the analysis for every AS pattern, patients with a high risk score had lower survival rate compared to patients with a low risk score, as presented in Figures 3(a)–3(f), of which the model with AD events had the best accuracy in predicting 5-year survival () (Figure 3(h)). Moreover, when all possible AS events were included for independent prognosis prediction from the different types (total 42 AS events), the model displayed favorable ability in predicting the 5-year survival (), as presented in Figures 3(g) and 3(h).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

3.4. Survival-Related AS-SF Interaction Network

To determine the prognostic value of SFs, survival analysis was applied for SFs according to gene expression. As shown in Figure 4(a), four SFs exhibited a significant correlation with OS, and decreased levels led to poor prognosis. We also clarified the relationships between SFs and survival-associated AS events and established a prognostic model. Most SF gene levels were negatively correlated with the PSI (Figure 4(b)).

(a)

(b)

(c)

Figure 4

SFs in glioma. (a) Survival curve of SFs. (b) Correlation analysis of the splicing factor and AS prognostic predictor. The upper right figure shows the significance of the correlation between the correlation coefficient and SF expression and the PSI value of prognostic-related AS events, where , , and ; the lower left figure shows a pie chart of the correlation between SF expression and the PSI value of survival-related AS events, where the size of the pie chart represents the correlation, the green line represents the negative correlation, and the red line represents the positive correlation. (c) The AS network, where the green node represents the gene involved in the differential AS event and the red node represents the AS factor.

Additionally, the correlations between PSI values of survival-associated AS event and survival-related SF expression were investigated by Spearman correlation analysis (Figure 4(a)). Four survival-related SFs (red dots) showed a significant correlation between 58 genes (green dots) and survival-associated AS events. There were 59 associations of the SF NOVA1 and 35 associations of the SF RBM4 with AS events. In addition, as shown in Figure 4(c), CCDC121, MAD2L2, TCF12, and ZNF138 interacted with all four SFs.

3.5. Prognosis-Related Molecular Subtype Clusters

To further identify different AS patterns, unsupervised analysis was conducted for all samples according to prognosis-related AS events. According to the results of the Consensus Cluster Plus analysis and the distribution of consensus values from 0 to 1, cluster 1 (, 47.3%) and cluster 2 (, 52.7%) were finally identified by the elbow method, as displayed in Figures 5(a) and 5(b). Subsequently, survival analysis was performed to evaluate the effect of these clusters on prognosis. The results showed that cluster 2 was related to poor survival, whereas cluster 1 was associated with favorable survival (Figure 5(c)). The related genes were significantly enriched in multiple GO terms and Reactome pathways, as presented in Figures 5(c)–5(g).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

Figure 5

AS clusters associated with prognosis and molecular subtypes: (a) the elbow method identifies the optimal number of clusters; (b) the consensus matrix heat map defines two sample clusters, the consensus value of which ranges from 0 (white, samples never clustered together) to 1 (dark blue, samples always cluster together); (c) survival analysis in the two identified sample clusters; (d–f) top 10 significantly enriched biological processes, cell components, and molecular functions in GO analysis; (g) top 10 significantly enriched Reactome pathways.

3.6. The Expression of Key SFs in Glioma

To confirm the expression of key SFs in glioma, we conducted western blotting based on glioma and nontumor tissues. The results indicated that the expression levels of these four SFs, including NOVA1, RBM4, HNRNPC, and HNRPLL, in tumor tissues were significantly decreased, compared to those in normal tissues (Figure 6).

(a)

(b)

(c)

(d)

3.7. The Correlation between Key SFs and Immune Cell Infiltration

To calculate the correlation between the SF expression and immune infiltrates, the TIMER database was applied with the key SF expression in TCGA dataset. As presented in Figure 7, the expressions of NOVA1 and RBM4 were positively correlated with the infiltration of CD8+ T cells. Besides, NOVA1 expression was significantly related to the infiltration level of CD4+ T cells, macrophages, neutrophils, and dendritic cells. High expression levels of HNRNPC contributed to the infiltration of neutrophils and dendritic cells. Although there was no significant difference between HNRPLL expression and the infiltration of immune cells, HNRPLL was closely related to immune purity. Besides, the results of CBERSORT analysis were consistent with the results of the TIMER database (Figure 8). The above findings suggested that there might be an underlying correlation between the AS events and immune infiltration.

(a)

(b)

(c)

(d)

4. Discussion

In recent years, accumulated evidence suggested that AS occurred in multiple tumors and played a critical role in the progression of malignancies, including glioma. Our study found that genes involved in AS events affected glioma survival including NOVA1, HNRNPC, HNRNPLL, RBM4. Other scholars have discovered more influencing factors. Li et al. [18] verified that CX-5461 suppressed telomerase activity and induces cell apoptosis via the regulation of hTERT in glioma cells. In addition, it has been reported that the expression level of ELK1 exhibited the positive correlation with tumor malignant progression, which was regulated by AS [19]. Shao et al. [20] found that two variants of ITSN1 via AS exercised distinctly different functions in glioma. Specifically, the short variant of ITSN1 promoted the development of glioma, whereas the long variant exerted a tumor-suppressive effect.

Previous studies examining AS have usually focused on a single gene or SF. Additionally, the prognostic value and immune infiltration analysis of AS in glioma have not yet been studied. To our knowledge, this study is the first to systematically identify and analyze AS events related to survival and immune cell infiltration in glioma. Based on our findings, 425 survival-related AS events were revealed, and SF genes related to survival-associated AS events were identified, including NOVA1, HNRNPC, HNRNPLL, and RBM4. These key SF genes were significantly related to the infiltration of immune cells in the tumor microenvironment. These factors exert vital roles in the genesis and development of multiple tumors [11, 21–32].

Among these four SFs, HNRNPC, a small nuclear ribonucleoprotein particle protein factor, is the most extensively investigated. Fischl et al. [23] reported that the HNRNPC-dependent alternative cleavage and polyadenylation (APA) profile changes in colonic tumors compared to noncarcinoma colonic epithelial cells. They also verified that HNRNPC is a key regulator of physiology-related APA events and postulated that HNRNPC might have facilitated gene expression related to proliferation and metastasis and promote carcinogenesis. The findings of Zhang et al. [24] demonstrated that HNRNPC interacted with the lncRNA LBX2-AS1 and transcription factors ZEB1 and ZEB2 related to the epithelial-mesenchymal transition (EMT), which further facilitated the progression of EMT. Wen et al. [10] discovered that HNRNPC could serve as a predictor for tumor burden and prognosis. Additionally, high expression level of HNRNPC contributed to poor survival of patients with gastric cancer receiving chemotherapy and resulted in chemoresistance [33].

HNRNPLL, a tissue-specific heterogeneous nuclear ribonucleoprotein, has been shown to be related to progression of colorectal cancer (CRC). Sakuma et al. [29] found that HNRNPLL suppressed metastasis via EMT. Interestingly, another one of their studies [27] suggested that HNRNPLL stabilized the transcripts of DNA replication regulators (PCNA, RFC3, and FEN1), thus accelerating the proliferation of CRC cells. It has been reported that HNRNPLL could reduce T cell accumulation in lymphoid tissues [34], which was consistent with the findings of our study. This indicates that an SF may impact tumor development and prognosis in multiple mechanisms. We demonstrated that NOVA1, HNRNPC, HNRNPLL, and RBM4 were independent prognostic factors for glioma, and their expression was verified by western blotting analysis. In contrast, NOVA1 is considered to be an oncogene and a prognostic factor in multiple malignant tumors, including CRC [35], thyroid cancer [36], breast cancer [37], melanoma [30], osteosarcoma [38], and gastric cancer [31, 39]. This phenomenon has also been observed for other SFs. We speculate that, unlike the effect of HNRNPLL in CRC, the overexpression of these four SFs may facilitate glioma metastasis but suppress glioma cell proliferation. However, considering the biological characteristics of glioma, it is difficult to study carcinogenesis and construct models to better predict patient prognosis. Notably, SFs suppressing cell proliferation also exist in other types of tumors [40–42]. However, the innate significant heterogeneity of glioma samples may lead to different results. Consequently, more studies are warranted to further verify this result and uncover the mechanism behind this phenomenon.

In conclusion, this study screened and analyzed AS events in glioma using TCGA database. We created prognostic signatures associated with 42 AS events, which displayed favorable performance in predicting prognosis in glioma patients. We also discovered that the SFs NOVA1, HNRNPC, HNRNPLL, and RBM4 are independent risk factors that affect the prognosis and immune cell infiltration of patients with glioma. Certain limitations should also be noted. First, our research was based on TCGA and UCSC Xena databases and was not verified using other independent databases. Second, the four prognosis-related SFs identified in our results may differ in terms of their expression in other tumor types, and the real mechanism behind their association with survival and immune infiltration should be further verified. In summary, our research illustrates the value of AS events and AS-related genes in glioma and presents a survival-related molecular subtype-SF interaction network, which provides a bioinformatics basis for subsequent study of related mechanisms.

Data Availability

The data used in this study are available from the corresponding authors upon request.

All the patients signed the informed consent.

Conflicts of Interest

The authors have declared no conflict of interest.

Authors’ Contributions

GQL, RJT, and ZT were responsible for the resources, project administration, funding acquisition, and supervision. SHX and SF were responsible for the writing and validation. YG and HCX were responsible for the methodology and visualization. Hong-xin Su, Gang Yang, and Fei Su contributed equally to this work.

References

H. Feng, Z. Qin, and X. Zhang, “Opportunities and methods for studying alternative splicing in cancer with RNA-Seq,” Cancer Letters, vol. 340, no. 2, pp. 179–191, 2013.
View at: Publisher Site | Google Scholar
H. Climente-Gonzalez, E. Porta-Pardo, A. Godzik, and E. Eyras, “The functional impact of alternative splicing in cancer,” Cell Reports, vol. 20, no. 9, pp. 2215–2226, 2017.
View at: Publisher Site | Google Scholar
M. M. Scotti and M. S. Swanson, “RNA mis-splicing in disease,” Nature Reviews Genetics, vol. 17, no. 1, pp. 19–32, 2016.
View at: Google Scholar
L. Shkreta, B. Bell, T. Revil et al., “Cancer-associated perturbations in alternative pre-messenger RNA splicing,” Cancer Treatment and Research, vol. 158, pp. 41–94, 2013.
View at: Publisher Site | Google Scholar
R. Q. He, X. G. Zhou, Q. Y. Yi et al., “Prognostic signature of alternative splicing events in bladder urothelial carcinoma based on Spliceseq data from 317 cases,” Cellular Physiology and Biochemistry, vol. 48, no. 3, pp. 1355–1368, 2018.
View at: Publisher Site | Google Scholar
Y. Li, N. Sun, Z. Lu et al., “Prognostic alternative mRNA splicing signature in non-small cell lung cancer,” Cancer Letters, vol. 393, pp. 40–51, 2017.
View at: Publisher Site | Google Scholar
J. Zhu, Z. Chen, and L. Yong, “Systematic profiling of alternative splicing signature reveals prognostic predictor for ovarian cancer,” Gynecologic Oncology, vol. 148, no. 2, pp. 368–374, 2018.
View at: Publisher Site | Google Scholar
P. Lin, R. Q. He, F. C. Ma et al., “Systematic analysis of survival-associated alternative splicing signatures in gastrointestinal pan-adenocarcinomas,” eBioMedicine, vol. 34, pp. 46–60, 2018.
View at: Publisher Site | Google Scholar
A. Bisognin, S. Pizzini, L. Perilli et al., “An integrative framework identifies alternative splicing events in colorectal cancer development,” Molecular Oncology, vol. 8, no. 1, pp. 129–141, 2014.
View at: Publisher Site | Google Scholar
J. Wen, K. H. Toomer, Z. Chen, and X. Cai, “Genome-wide analysis of alternative transcripts in human breast cancer,” Breast Cancer Research and Treatment, vol. 151, no. 2, pp. 295–307, 2015.
View at: Publisher Site | Google Scholar
S. Zhang, X. Wu, P. Diao et al., “Identification of a prognostic alternative splicing signature in oral squamous cell carcinoma,” Journal of cellular physiology, vol. 235, no. 5, pp. 4804–4813, 2020.
View at: Google Scholar
R. Zhang, P. Lin, X. Yang et al., “Survival associated alternative splicing events in diffuse large B-cell lymphoma,” American Journal of Translational Research, vol. 10, no. 8, pp. 2636–2647, 2018.
View at: Google Scholar
Z. X. Li, Z. Q. Zheng, Z. H. Wei et al., “Comprehensive characterization of the alternative splicing landscape in head and neck squamous cell carcinoma reveals novel events associated with tumorigenesis and the immune microenvironment,” Theranostics, vol. 9, no. 25, pp. 7648–7665, 2019.
View at: Publisher Site | Google Scholar
T. Tang, X. Cheng, B. Truong, L. Sun, X. Yang, and H. Wang, “Molecular basis and therapeutic implications of CD40/CD40L immune checkpoint,” Pharmacology & Therapeutics, vol. 219, article 107709, 2021.
View at: Publisher Site | Google Scholar
Q. T. Ostrom, G. Cioffi, H. Gittleman et al., “CBTRUS statistical report: primary brain and other central nervous system tumors diagnosed in the United States in 2012-2016,” Neuro-Oncology, vol. 21, Supplement_5, pp. v1–v100, 2019.
View at: Publisher Site | Google Scholar
P. Diamandis and K. D. Aldape, “Insights from molecular profiling of adult glioma,” Journal of Clinical Oncology, vol. 35, no. 21, pp. 2386–2393, 2017.
View at: Publisher Site | Google Scholar
T. Li, J. Fan, B. Wang et al., “TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells,” Cancer Research, vol. 77, no. 21, pp. e108–e110, 2017.
View at: Publisher Site | Google Scholar
G. Li, J. Shen, J. Cao et al., “Alternative splicing of human telomerase reverse transcriptase in gliomas and its modulation mediated by CX-5461,” Journal of Experimental & Clinical Cancer Research, vol. 37, no. 1, pp. 1–13, 2018.
View at: Google Scholar
J. Li, Y. Wang, X. Meng, and H. Liang, “Modulation of transcriptional activity in brain lower grade glioma by alternative splicing,” PeerJ, vol. 6, article e4686, 2018.
View at: Google Scholar
Y. Shao, W. Chong, X. Liu et al., “Alternative splicing-derived intersectin1-L and intersectin1-S exert opposite function in glioma progression,” Cell Death & Disease, vol. 10, no. 6, p. 431, 2019.
View at: Publisher Site | Google Scholar
Y. J. Zhou, G. Q. Zhu, Q. W. Zhang et al., “Survival-associated alternative messenger RNA splicing signatures in pancreatic ductal adenocarcinoma: a study based on RNA-sequencing data,” DNA and Cell Biology, vol. 38, no. 11, pp. 1207–1222, 2019.
View at: Publisher Site | Google Scholar
X. Liu, L. Liu, Z. Dong et al., “Expression patterns and prognostic value of m₆A-related genes in colorectal cancer,” American Journal of Translational Research, vol. 11, no. 7, pp. 3972–3991, 2019.
View at: Google Scholar
H. Fischl, J. Neve, Z. Wang et al., “hnRNPC regulates cancer-specific alternative cleavage and polyadenylation profiles,” Nucleic Acids Research, vol. 47, no. 14, pp. 7580–7591, 2019.
View at: Publisher Site | Google Scholar
Y. Zhang, W. Chen, T. Pan, H. Wang, Y. Zhang, and C. Li, “LBX2-AS1 is activated by ZEB1 and promotes the development of esophageal squamous cell carcinoma by interacting with HNRNPC to enhance the stability of ZEB1 and ZEB2 mRNAs,” Biochemical and Biophysical Research Communications, vol. 511, no. 3, pp. 566–572, 2019.
View at: Publisher Site | Google Scholar
M. P. Tremblay, V. E. Armero, A. Allaire et al., “Global profiling of alternative RNA splicing events provides insights into molecular differences between various types of hepatocellular carcinoma,” BMC Genomics, vol. 17, no. 1, p. 683, 2016.
View at: Publisher Site | Google Scholar
Y. A. Zhang, J. M. Zhu, J. Yin et al., “High expression of neuro-oncological ventral antigen 1 correlates with poor prognosis in hepatocellular carcinoma,” PLoS One, vol. 9, no. 3, article e90955, 2014.
View at: Publisher Site | Google Scholar
K. Sakuma, E. Sasaki, K. Kimura et al., “HNRNPLL stabilizes mRNA for DNA replication proteins and promotes cell cycle progression in colorectal cancer cells,” Cancer Science, vol. 109, no. 8, pp. 2458–2468, 2018.
View at: Publisher Site | Google Scholar
S. Baranwal, “Commentary: HNRNPLL, a newly identified colorectal cancer metastasis suppressor, modulates alternative splicing of CD44 during epithelial-mesenchymal transition,” Frontiers in Cell and Development Biology, vol. 5, p. 91, 2017.
View at: Publisher Site | Google Scholar
K. Sakuma, E. Sasaki, K. Kimura et al., “HNRNPLL, a newly identified colorectal cancer metastasis suppressor, modulates alternative splicing of CD44 during epithelial-mesenchymal transition,” Gut, vol. 67, no. 6, pp. 1103–1111, 2018.
View at: Publisher Site | Google Scholar
X. Yu, H. Zheng, M. T. V. Chan, and W. K. K. Wu, “NOVA1 acts as an oncogene in melanoma via regulating FOXO3a expression,” Journal of Cellular and Molecular Medicine, vol. 22, no. 5, pp. 2622–2630, 2018.
View at: Publisher Site | Google Scholar
E. K. Kim, S. O. Yoon, W. Y. Jung et al., “Implications of NOVA1 suppression within the microenvironment of gastric cancer: association with immune cell dysregulation,” Gastric Cancer, vol. 20, no. 3, pp. 438–447, 2017.
View at: Publisher Site | Google Scholar
M. A. Markus, Y. H. Yang, and B. J. Morris, “Transcriptome-wide targets of alternative splicing by RBM4 and possible role in cancer,” Genomics, vol. 107, no. 4, pp. 138–144, 2016.
View at: Publisher Site | Google Scholar
H. Huang, Y. Han, C. Zhang et al., “HNRNPC as a candidate biomarker for chemoresistance in gastric cancer,” Tumour Biology, vol. 37, no. 3, pp. 3527–3534, 2016.
View at: Publisher Site | Google Scholar
M. van den Buuse, P. Halley, and G. F. Hoyne, “Behavioural phenotyping of thunder mice with a hypomorphic mutation of heterogeneous nuclear ribonuclear protein L-like (hnRNPLL) and reduced T cell function,” Neuroscience Letters, vol. 740, article 135469, 2021.
View at: Publisher Site | Google Scholar
Y. G. Hong, G. S. Xu, G. Y. Yu et al., “The RNA binding protein neuro-oncological ventral antigen 1 (NOVA1) regulates IL-6 mRNA stability to enhance JAK2-STAT3 signaling in CRC,” Surgical Oncology, vol. 31, pp. 67–74, 2019.
View at: Google Scholar
Y. Luo, T. Hao, J. Zhang, M. Zhang, P. Sun, and L. Wu, “MicroRNA-592 suppresses the malignant phenotypes of thyroid cancer by regulating lncRNA NEAT1 and downregulating NOVA1,” International Journal of Molecular Medicine, vol. 44, no. 3, pp. 1172–1182, 2019.
View at: Publisher Site | Google Scholar
X. Y. Yu, J. R. Tian, D. Yang, and H. R. Tan, “TATDN1 promotes the development and progression of breast cancer by targeting microRNA-140-3p,” European Review for Medical and Pharmacological Sciences, vol. 23, no. 12, pp. 5293–5300, 2019.
View at: Publisher Site | Google Scholar
C. Li, Y. He, H. Ma, and S. Han, “NOVA1 acts as an oncogene in osteosarcoma,” American Journal of Translational Research, vol. 9, no. 10, pp. 4450–4457, 2017.
View at: Google Scholar
S. O. Yoon, E. K. Kim, M. Lee et al., “NOVA1 inhibition by miR-146b-5p in the remnant tissue microenvironment defines occult residual disease after gastric cancer removal,” Oncotarget, vol. 7, no. 3, pp. 2475–2495, 2016.
View at: Publisher Site | Google Scholar
H. Yong, W. Zhao, X. Zhou et al., “RNA-binding motif 4 (RBM4) suppresses tumor growth and metastasis in human gastric cancer,” Medical Science Monitor, vol. 25, pp. 4025–4034, 2019.
View at: Publisher Site | Google Scholar
Y. Wang, D. Chen, H. Qian et al., “The splicing factor RBM4 controls apoptosis, proliferation, and migration to suppress tumor progression,” Cancer Cell, vol. 26, no. 3, pp. 374–389, 2014.
View at: Publisher Site | Google Scholar
“RBM4-regulated alternative splicing suppresses tumorigenesis,” Cancer Discovery, vol. 4, no. 11, pp. 1253.1–121253, 2014.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Hong-xin Su et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

268

Downloads

539

Citations