Comprehensive Analysis of Prognostic Alternative Splicing Signatures in Oral Squamous Cell Carcinoma

Splicing Machinery Is Impaired in Oral Squamous Cell Carcinomas and Linked to Key Pathophysiological Features

Alternative splicing dysregulation is an emerging cancer hallmark, potentially serving as a source of novel diagnostic, prognostic, or therapeutic tools. Inhibitors of the activity of the splicing machinery can exert antitumoral effects in cancer cells. We aimed to characterize the splicing machinery (SM) components in oral squamous cell carcinoma (OSCC) and to evaluate the direct impact of the inhibition of SM-activity on OSCC-cells. The expression of 59 SM-components was assessed using a prospective case-control study of tumor and healthy samples from 37 OSCC patients, and the relationship with clinical and histopathological features was assessed. The direct effect of pladienolide-B (SM-inhibitor) on the proliferation rate of primary OSCC cell cultures was evaluated. A significant dysregulation in several SM components was found in OSCC vs. adjacent-healthy tissues [i.e., 12 out of 59 (20%)], and their expression was associated with clinical and histopathological features of less aggressiveness and overall survival. Pladienolide-B treatment significantly decreased OSCC-cell proliferation. Our data reveal a significantly altered expression of several SM-components and link it to pathophysiological features, reinforcing a potential clinical and pathophysiological relevance of the SM dysregulation in OSCC. The inhibition of SM-activity might be a therapeutic avenue in OSCC, offering a clinically relevant opportunity to be explored.

An anti-PD-1 antisense oligonucleotide promotes the expression of soluble PD-1 by blocking the interaction between SRSF3 and an exonic splicing enhancer of PD-1 exon 3

PD-1 is a key immune checkpoint molecule. Anti-PD-1 immunotherapy is encouraging in cancer treatment. However, it still needs to be improved. PD-1 has at least five isoforms generated by alternative splicing. An isoform without exon 3 encoding soluble PD-1 (sPD-1) can activate anti-tumor immunity by inhibiting the interaction between cellular surface full-length PD-1 (flPD-1) and PD-L1. However, the regulatory mechanism of exon 3 splicing remains largely unknown. Here, we screened the exon 3 sequence by mutation and searched corresponding splicing factors by SpliceAid database and pulldown assay. The alternative splicing of PD-1 exon 3 was analyzed by RT-PCR. The expression levels of flPD-1 and sPD-1 were analyzed by Western blot, flow cytometry, and ELISA. We discovered that an exonic splicing enhancer (ESE) of exon 3 is essential for its inclusion. Moreover, SRSF3 can bind to this ESE and enhance exon 3 inclusion and flPD-1 expression. We designed and screened out an antisense oligonucleotide (ASO) targeting PD-1 to block the interaction between SRSF3 and ESE, and significantly increase exon 3 skipping and sPD-1 expression, which was verified in various tumor cells in addition to oral cancer cells. Altogether, our results uncovered the regulatory mechanism of human PD-1 exon 3 splicing and sPD-1 expression and further designed a novel anti-PD-1 ASO, which are useful for developing a new method of anti-cancer immunotherapy.

Upregulation of miRNA-196a and miRNA-196b correlates with Bryne's prognostic score in oral squamous cell carcinoma

BACKGROUND

microRNA(miRNA)-196a and miRNA-196b expression has been found to be dysregulated and involved in tumorigenesis and tumor progression in array of cancers through different targets. The role of these miRNAs together in clinical application is not always consistent and, its prognostic value in oral squamous cell carcinoma (OSCC) is still elusive. This study was performed to investigate the correlation of these miRNAs expression with histological grades of OSCC according to Bryne's histological grading system, to predict prognosis and to evaluate their relationship with clinico-pathological data.

METHODS

Real-time quantitative Reverse Transcription Polymerase Chain Reaction (qRT-PCR) was done to evaluate the expression levels of miRNA-196a and miRNA-196b in 75 pairs of OSCC tissue matched with adjacent normal mucosa, used as a control.

RESULTS

miRNA-196a and 196b expression in OSCC was significantly higher than that in corresponding adjacent normal tissues (p > 0.001). Also, a significant differential correlation was found in between the expression levels of these two miRNAs (Pearson correlation test r = 0.676, p-value<0.0001). The increased expression of these miRNAs was more frequently observed in OSCC tissues with advanced clinical and pathological TNM stages (IVa and IVb, pIVb respectively, p-value<0.0001). Significant correlation was found between miRNA-196a upregulation and moderate prognostic score (p-value<0.0001) in comparison with good and poor prognostic score of histological grades of OSCC. Sensitivity and specificity for miRNA-196a were 95 % and 85 %, respectively (AUC = 1, 95 % CI = 0.617-0.850; p 0.001), while for miRNA-196b were 94 % and 86 %, respectively (AUC = 0.808, 95 % CI = 0.701-0.916; p0.0001).

CONCLUSIONS

These findings suggest that the increased expression of miRNA-196a and 196b may play an important role in tumor progression in OSCC. miRNA-196a might be a useful marker for predicting the clinical outcome of OSCC, especially for advanced stages. In conclusion, our data demonstrate for the first time that these miRNAs may serve as a potent prognostic marker for tumor progression. We further highlight miRNA-196a and miRNA-196b as a promising predictor of prognostic assessment in OSCC.

CD44 occurring alternative splicing promotes cisplatin resistance and evokes tumor immune response in oral squamous cell carcinoma cells

BACKGROUND

Oral squamous cell carcinoma (OSCC) is the most prevalent malignant tumor in head and neck region. Platinum drug resistance limits the clinical application of chemotherapy regardless of medical development. The aim of our study is to identify cisplatin-resistant genes which can be used as new therapeutic targets and investigate the functional mechanism of OSCC chemoresistance.

METHODS

The OSCC Cal27 and HSC4 cisplatin-resistant cell lines were constructed to screen the differential genes/transcripts expression. GO, KEGG and GSEA were performed to reveal the relevant signaling pathways. Alternative splicing (AS) software rMATs was applied to explore AS events in chemoresistance. R package and TIMER tools were used to evaluate the linear correlation between CD44 and immune cell subpopulations. The co-culture model of dendritic cells (DCs) and OSCC cells was applied to explore the effect of CD44 on immune microenvironment and cisplatin resistance.

RESULTS

Our results showed that CD44 was differentially expressed in cisplatin-resistant OSCC cells. Through bioinformatics prediction and experimental verification, we confirmed that CD44 occurring AS was involved in tumor progression and cisplatin resistance. Moreover, CD44 could further enhance the cisplatin resistance of OSCC by activating DCs, making CD44 to be a potential intervention target. We also identified DC as a new target for platinum drugs to stimulate the growth of OSCC.

CONCLUSION

Our findings not only make it possible to explore new therapeutic methods, such as CD44 inhibitors or antisense oligonucleotides, but also provide insights into the new mechanisms of cisplatin resistance to chemotherapy.

COL3A1: Potential prognostic predictor for head and neck cancer based on immune-microenvironment alternative splicing

We aimed to identify a novel prognostic biomarker for head and neck squamous cell carcinoma (HNSCC) based on tumor immunology-related alternative splicing (AS). Data for 502 HNSCC and 44 normal samples were obtained from the TCGA database and used to establish an AS-related risk model through univariate, least absolute shrinkage, and selection operator Cox regression analyses. Fresh HNSCC and normal oral tissues were surgically obtained from 44 HNSCC patients. Western blotting and quantitative reverse transcription-PCR were used to assess gene expression levels. Kaplan-Meier was performed to evaluate patients' overall survival (OS) rate. The CIBERSORT algorithm, single-sample gene set enrichment analysis, and immune checkpoint analyses were performed to compare immune activities between subgroups. The risk model was established using 10 pivotal AS events first. Collagen Type III Alpha 1 Chain (COL3A1) were screened based on |log2FC| ≥ 1 and FDR < 0.05 criteria. COL3A1 expression levels in HNSCC tissues were elevated relative to normal tissues (p < 0.001). Moreover, COL3A1 was a reliable biomarker for HNSCC patients' prognostic prediction in both cohorts (p < 0.001, p = 0.0085, respectively). COL3A1 protein (p = 0.0054) and mRNA (p < 0.0001) levels were correlated with HNSCC differentiation. Furthermore, the T stage was correlated with COL3A1 expression (p = 0.043), and COL3A1 expression was an independent prognostic predictor for HNSCC patients (p = 0.006). Compared with the risk model, COL3A1 was better at evaluating immune cell infiltrations, immune activities, and immune checkpoint gene expressions of HNSCC lesions.

IQGAP1 and RNA Splicing in the Context of Head and Neck via Phosphoproteomics

IQGAP1 (IQ motif-containing GTPase-activating protein 1) scaffolds several signaling pathways in mammalian cells that are implicated in carcinogenesis, including the RAS and PI3K pathways that involve multiple protein kinases. IQGAP1 has been shown to promote head and neck squamous cell carcinoma (HNSCC); however, the underlying mechanism(s) remains unclear. Here, we report a mass spectrometry-based analysis identifying differences in phosphorylation of cellular proteins in vivo and in vitro in the presence or absence of IQGAP1. By comparing the esophageal phosphoproteome profiles between Iqgap1+/+ and Iqgap1-/- mice, we identified RNA splicing as one of the most altered cellular processes. Serine/arginine-rich splicing factor 6 (SRSF6) was the protein with the most downregulated levels of phosphorylation in Iqgap1-/- tissue. We confirmed that the absence of IQGAP1 reduced SRSF6 phosphorylation both in vivo and in vitro. We then expanded our analysis to human normal oral keratinocytes. Again, we found factors involved in RNA splicing to be highly altered in the phosphoproteome profile upon genetic disruption of IQGAP1. Both the Clinical Proteomic Tumor Analysis Consortium (CPTAC) and the Cancer Genome Atlas (TCGA) data sets indicate that phosphorylation of splicing-related proteins is important in HNSCC prognosis. The Biological General Repository for Interaction Datasets (BioGRID) repository also suggested multiple interactions between IQGAP1 and splicing-related proteins. Based on these collective observations, we propose that IQGAP1 regulates the phosphorylation of splicing proteins, which potentially affects their splicing activities and, therefore, contributes to HNSCC. Raw data are available from the MassIVE database with identifier MSV000087770.

High expression of serine and arginine-rich splicing factor 9 (SRSF9) is associated with hepatocellular carcinoma progression and a poor prognosis

BACKGROUND

Serine and arginine-rich splicing factor 9 (SRSF9) has been linked to the occurrence and progression of various cancers; however, its effects and mechanism of action hepatocellular carcinoma (HCC) have not been reported. In this study, we used a bioinformatics approach and in vitro assays to evaluate the expression of SRSF9 in HCC, its prognostic value, and its underlying regulatory mechanisms, including analyses of related pathways and the role of methylation.

METHODS

Transcriptomic and DNA methylation data for 357 HCC cases and 50 paratumor tissues in The Cancer Genome Atlas database were obtained. Additionally, protein expression data for cell lines and tissue samples were obtained from the Human Protein Atlas. The CMap databased was used to predict candidate drugs targeting SRSF9. Various cell lines were used for in vitro validation.

RESULTS

SRSF9 expression was significantly elevated in HCC and was negatively regulated by its methylation site cg06116271. The low expression of SRSF9 and hypermethylation of cg06116271 were both associated with a longer overall survival time. A correlation analysis revealed ten genes that were co-expressed with SRSF9; levels of CDK4, RAN, DENR, RNF34, and ANAPC5 were positively correlated and levels of RBP4, APOC1, MASP2, HP, and HPX were negatively correlated with SRSF9 expression. The knockdown of SRSF9 in vitro inhibited the proliferation and migration of HCC cells and significantly reduced the expression of proteins in the Wnt signaling pathway (DVL2 and β-catenin) and cell cycle pathway (Cyclin D and Cyclin E). A CMap analysis identified two drugs, camptothecin and apigenin, able to target and inhibit the expression of SRSF9.

CONCLUSIONS

This study expands our understanding of the molecular biological functions of SRSF9 and cg06116271 and provides candidate diagnostic and therapeutic targets for HCC.

CELF Family Proteins in Cancer: Highlights on the RNA-Binding Protein/Noncoding RNA Regulatory Axis

Post-transcriptional modifications to coding and non-coding RNAs are unquestionably a pivotal way in which human mRNA and protein diversity can influence the different phases of a transcript's life cycle. CELF (CUGBP Elav-like family) proteins are RBPs (RNA-binding proteins) with pleiotropic capabilities in RNA processing. Their responsibilities extend from alternative splicing and transcript editing in the nucleus to mRNA stability, and translation into the cytoplasm. In this way, CELF family members have been connected to global alterations in cancer proliferation and invasion, leading to their identification as potential tumor suppressors or even oncogenes. Notably, genetic variants, alternative splicing, phosphorylation, acetylation, subcellular distribution, competition with other RBPs, and ultimately lncRNAs, miRNAs, and circRNAs all impact CELF regulation. Discoveries have emerged about the control of CELF functions, particularly via noncoding RNAs, and CELF proteins have been identified as competing, antagonizing, and regulating agents of noncoding RNA biogenesis. On the other hand, CELFs are an intriguing example through which to broaden our understanding of the RBP/noncoding RNA regulatory axis. Balancing these complex pathways in cancer is undeniably pivotal and deserves further research. This review outlines some mechanisms of CELF protein regulation and their functional consequences in cancer physiology.

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.

Identification of Prognostic alternative splicing signatures and their clinical significance in uveal melanoma

As a posttranscriptional regulatory mechanism, alternative splicing (AS) has the potential to generate a large amount of protein diversity from limited genes. The purpose of our study was to assess the usefulness of prognostic splicing events as novel diagnostic and therapeutic signatures for uveal melanoma (UM). The datasets, clinical traits and AS data of UM were obtained from The Cancer Genome Atlas (TCGA) database and TCGA SpliceSeq database. Using bioinformatics analysis, we identified 1047 AS events as candidate AS events closely related to prognosis from 920 parent genes. The gene enrichment analysis indicated that these genes were mainly enriched in cellular components (CC) including cytosol, nucleoplasm, cytoplasm and ribosome, and in molecular functions (MF), including protein binding and poly(A) RNA binding. Furthermore, we selected all survival-associated splicing events to generate prognostic signatures, which included 4 exon skip (ES) events (DNASE1L1-90581-ES, NUDT1-78611-ES, BIN1-55198-ES, SEPN1-1195-ES) and 1 alternate promoter (AP) event (DPYSL2-83132-AP). The AS prognostic model was confirmed as independent overall survival (OS)-related factors (p = 0.014). A total of 17 splicing factors (SFs) involved in the regulation of AS were identified as related to the OS of UM patients. Our pooled data highlighted the usefulness and importance of AS biomarkers, which provided a potential strategy for the diagnosis and treatment of UM.

Genome-Wide Analyses of Prognostic and Therapeutic Alternative Splicing Signatures in Bladder Urothelial Carcinoma

Background

Alternative splicing (AS) is an indispensable post-transcriptional modification applied during the maturation of mRNA, and AS defects have been associated with many cancers. This study was designed to thoroughly analyze AS events in bladder urothelial carcinoma (BLCA) at the genome-wide level.

Methods

We adopted a gap analysis to screen for significant differential AS events (DASEs) associated with BLCA. DASEs with prognostic value for OS and the disease-free interval (DFI) were identified by Cox analysis. In addition, a differential AS network and AS clusters were identified using unsupervised cluster analysis. We examined differences in the sensitivity to chemotherapy and immunotherapy between BLCA patients with high and low overall survival (OS) risk.

Results

An extensive number of DASEs (296) were found to be clinically relevant in BLCA. A prognosis model was established based prognostic value of OS and DFI. CUGBP elav-like family member 2 (CELF2) was identified as a hub splicing factor for AS networks. We also identified AS clusters associated with OS using unsupervised cluster analysis, and we predicted that the effects of cisplatin and gemcitabine chemotherapy would be different between high- and low-risk groups based on OS prognosis.

Conclusion

We completed a comprehensive analysis of AS events in BLCA at the genome-wide level. The present findings revealed that DASEs and splicing factors tended to impact BLCA patient survival and sensitivity to chemotherapy drugs, which may provide novel prospects for BLCA therapies.