Identification of a prognostic alternative splicing signature in oral squamous cell carcinoma

Alternative splicing (AS) is critically associated with tumorigenesis and patient's prognosis. Here, we systematically analyzed survival-associated AS signatures in oral squamous cell carcinoma (OSCC) and evaluated their prognostic predictive values. Survival-related AS events were identified by univariate and multivariate Cox regression analyses using OSCC data from the TCGA head neck squamous cell carcinoma data set. The Percent Spliced In calculated by SpliceSeq from 0 to 1 was used to quantify seven types of AS events. A predictive model based on AS events was constructed by least absolute shrinkage and selection operator Cox regression assay and further validated using a training-testing cohort design. Patient survival was estimated using the Kaplan–Meier method and compared with Log-rank test. The receiver operating characteristics curve area under the curves was used to evaluate the predictive abilities of these predictive models. Furthermore, gene–gene interaction networks and the splicing factors (SFs)-AS regulatory network was generated by Cytoscape. A total of 825 survival-related AS events within 719 genes were identified in OSCC samples. The integrative predictive model was better at predicting outcomes of patients as compared to those models built with the individual AS event. The predictive model based on three AS-related genes also effectively predicted patients’ survival. Moreover, seven survival-related SFs were detected in OSCC including RBM4, HNRNPD, and HNRNPC, which have been linked to tumorigenesis. The SF-AS network revealed a significant correlation between survival-related AS genes and these SFs. Our findings revealed a systemic portrait of survival-associated AS events and the splicing network in OSCC, suggesting that AS events might serve as novel prognostic biomarkers and therapeutic targets for OSCC.     

Systematic analysis of alternative splicing signature unveils prognostic predictor for kidney renal clear cell carcinoma

There is growing evidence that alternative splicing (AS) plays an important role in cancer development. However, a comprehensive analysis of AS signatures in kidney renal clear cell carcinoma (KIRC) is lacking and urgently needed. It remains unclear whether AS acts as diagnostic biomarkers in predicting the prognosis of KIRC patients. In the work, gene expression and clinical data of KIRC were obtained from The Cancer Genome Atlas (TCGA), and profiles of AS events were downloaded from the SpliceSeq database. The RNA sequence/AS data and clinical information were integrated, and we conducted the Cox regression analysis to screen survival-related AS events and messenger RNAs (mRNAs). Correlation between prognostic AS events and gene expression were analyzed using the Pearson correlation coefficient. Protein-protein interaction analysis was conducted for the prognostic AS-related genes, and a potential regulatory network was built using Cytoscape (version 3.6.1). Meanwhile, functional enrichment analysis was conducted. A prognostic risk score model is then established based on seven hub genes (KRT222, LENG8, APOB, SLC3A1, SCD5, AQP1, and ADRA1A) that have high performance in the risk classification of KIRC patients. A total 46,415 AS events including 10,601 genes in 537 patients with KIRC were identified. In univariate Cox regression analysis, 13,362 survival associated AS events and 8,694 survival-specific mRNAs were detected. Common 3,105 genes were screen by overlapping 13,362 survival associated AS events and 8,694 survival-specific mRNAs. The Pearson correlation analysis suggested that 13 genes were significantly correlated with AS events (Pearson correlation coefficient >0.8 or <−0.8). Then, We conducted multivariate Cox regression analyses to select the potential prognostic AS genes. Seven genes were identified to be significantly related to OS. A prognostic model based on seven genes was constructed. The area under the ROC curve was 0.767. In the current study, a robust prognostic prediction model was constructed for KIRC patients, and the findings revealed that the AS events could act as potential prognostic biomarkers for KIRC.     

Downstream neighbor of SON (DONSON) is associated with unfavorable survival across diverse cancers with oncogenic properties in clear cell renal cell carcinoma

A precise stratification of our patients is essential and can support clinicians to determine the right therapy. The aim of this study was to identify clinically relevant genes using The Cancer Genome Atlas (TCGA) datasets.A comprehensive pan-cancer analysis of 30 distinct tumor entities (N = 9022) identified the largely unknown gene Downstream neighbor of SON (DONSON) to be particularly associated with unfavorable overall survival in clear cell renal cell carcinoma (KIRC). This prognostic potential of DONSON was validated in an independent KIRC cohort via quantitative real-time PCR (n = 152). Further, DONSON protein expression was evaluated via immunohistochemical staining followed by quantitative image analysis using the image analysis software QuPath on a renal cancer tissue microarray (n = 270).Interestingly, DONSON overexpression was preferentially associated with poor survival in 9 of the 30 entities, suggesting tumor-independent oncogenic properties of this largely unknown gene. A particularly strong association of DONSON to an aggressive phenotype was evident in KIRC and proved to be a strong independent predictor of unfavorable overall survival in two additional cohorts on the mRNA and protein level. In our KIRC cell culture model, we observed a substantial attenuation of proliferative activity and migration capacity of the KIRC cells Caki1 and 769p.In conclusion, we identified DONSON as a robust biomarker for risk stratification in KIRC in three independent cohorts and provide evidence that DONSON is linked to a malignant phenotype in the KIRC cell culture model.     

Identification of novel prognostic alternative splicing signature in papillary renal cell carcinoma

Papillary renal cell carcinoma (pRCC) is a heterogeneous disease containing multifocal or solitary tumors with an aggressive phenotype. Increasing evidence has indicated the involvement of aberrant splicing variants in renal cell cancer, while systematic profiling of aberrant alternative splicing (AS) in pRCC was lacking and largely unknown. In the current study, comprehensive profiling of AS events were performed based on the integration of pRCC cohort from the Cancer Genome Atlas database and SpliceSeq software. With rigorous screening and univariate Cox analysis, a total of 2077 prognoses AS events from 1642 parent genes were identified. Then, stepwise least absolute shrinkage and selection operator method-penalized Cox regression analyses with 10-fold cross-validation followed by multivariate Cox regression were used to construct the prognostic AS signatures within each AS type. And a final 21 AS event-based signature was proposed which showed potent prognostic capability in stratifying patients into low- and high-risk subgroups (P < .0001). Furthermore, time-dependent receiver operating characteristics curves confirmed that the final AS signature was effective and robust in predicting overall survival for pRCC patients with the area under the curve above 0.9 from 1 to 5 years. In addition, splicing correlation network was built to uncover the potential regulatory pattern among prognostic splicing factors and candidate AS events. Besides, gene set enrichment analysis revealed the involvement of these candidates AS events in tumor-related pathways including extracellular matrix organization, oxidative phosphorylation, and P53 signaling pathways. Taken together, our results could contribute to elucidating the underlying mechanism of AS in the oncogenesis process and broaden the novel field of prognostic and clinical application of molecule-targeted approaches in pRCC.     

Unique protein expression signatures of survival time in kidney renal clear cell carcinoma through a pan-cancer screening

Background

In 2016, it is estimated that there will be 62,700 new cases of kidney cancer in the United States, and 14,240 patients will die from the disease. Because the incidence of kidney renal clear cell carcinoma (KIRC), the most common type of kidney cancer, is expected to continue to increase in the US, there is an urgent need to find effective diagnostic biomarkers for KIRC that could help earlier detection of and customized treatment strategies for the disease. Accordingly, in this study we systematically investigated KIRC’s prognostic biomarkers for survival using the reverse phase protein array (RPPA) data and the high throughput sequencing data from The Cancer Genome Atlas (TCGA).

Results

With comprehensive data available in TCGA, we systematically screened protein expression based survival biomarkers in 10 major cancer types, among which KIRC presented many protein prognostic biomarkers of survival time. This is in agreement with a previous report that expression level changes (mRNAs, microRNA and protein) may have a better performance for prognosis of KIRC. In this study, we also identified 52 prognostic genes for KIRC, many of which are involved in cell-cycle and cancer signaling, as well as 15 tumor-stage-specific prognostic biomarkers. Notably, we found fewer prognostic biomarkers for early-stage than for late-stage KIRC. Four biomarkers (the RPPA protein IDs: FASN, ACC1, Cyclin_B1 and Rad51) were found to be prognostic for survival based on both protein and mRNA expression data.

Conclusions

Through pan-cancer screening, we found that many protein biomarkers were prognostic for patients’ survival in KIRC. Stage-specific survival biomarkers in KIRC were also identified. Our study indicated that these protein biomarkers might have potential clinical value in terms of predicting survival in KIRC patients and developing individualized treatment strategies. Importantly, we found many biomarkers in KIRC at both the mRNA expression level and the protein expression level. These biomarkers shared a significant overlap, indicating that they were technically replicable.

     

Identification of survival-related alternative splicing signatures in acute myeloid leukemia

Aberrant RNA alternative splicing (AS) variants play critical roles in tumorigenesis and prognosis in human cancers. Here, we conducted a comprehensive profiling of aberrant AS events in acute myeloid leukemia (AML). RNA AS profile, including seven AS types, and the percent spliced in (PSI) value for each patient were generated by SpliceSeq using RNA-seq data from TCGA. Univariate followed by multivariate Cox regression analysis were used to identify survival-related AS events and develop the AS signatures. A nomogram was developed, and its predictive efficacy was assessed. About 27,892 AS events and 3,178 events were associated with overall survival (OS) after strict filtering. Parent genes of survival-associated AS events were mainly enriched in leukemia-associated processes including chromatin modification, autophagy, and T-cell receptor signaling pathway. The 10 AS signature based on seven types of AS events showed better efficacy in predicting OS of patients than those built on a single AS event type. The area under curve (AUC) value of the 10 AS signature for 3-year OS was 0.91. Gene set enrichment analysis (GSEA) confirmed that these survival-related AS events contribute to AML progression. Moreover, the nomogram showed good predictive performance for patient''s prognosis. Finally, the correlation network of AS variants with splicing factor genes found potential important regulatory genes in AML. The present study presented a systematic analysis of survival-related AS events and developed AS signatures for predicting the patient’s survival. Further studies are needed to validate the signatures in independent AML cohorts and might provide a promising perspective for developing therapeutic targets.     

Comprehensively analysis of splicing factors to construct prognosis prediction classifier in prostate cancer

Splicing factors (SFs) are proteins that control the alternative splicing (AS) of RNAs, which have been recognized as new cancer hallmarks. Their dysregulation has been found to be involved in many biological processes of cancer, such as carcinogenesis, proliferation, metastasis and senescence. Dysregulation of SFs has been demonstrated to contribute to the progression of prostate cancer (PCa). However, a comprehensive analysis of the prognosis value of SFs in PCa is limited. In this work, we systematically analysed 393 SFs to deeply characterize the expression patterns, clinical relevance and biological functions of SFs in PCa. We identified 53 survival-related SFs that can stratify PCa into two de nove molecular subtypes with distinct mRNA expression and AS-event expression patterns and displayed significant differences in pathway activity and clinical outcomes. An SF-based classifier was established using LASSO-COX regression with six key SFs (BCAS1, LSM3, DHX16, NOVA2, RBM47 and SNRPN), which showed promising prognosis-prediction performance with a receiver operating characteristic (ROC) >0.700 in both the training and testing datasets, as well as in three external PCa cohorts (DKFZ, GSE70769 and GSE21035). CRISPR/CAS9 screening data and cell-level functional analysis suggested that LSM3 and DHX16 are essential factors for the proliferation and cell cycle progression in PCa cells. This study proposes that SFs and AS events are potential multidimensional biomarkers for the diagnosis, prognosis and treatment of PCa.     

Construction and validation of a seven-gene signature for predicting overall survival in patients with kidney renal clear cell carcinoma via an integrated bioinformatics analysis

ABSTRACT

Kidney renal clear cell carcinoma (KIRC) remains a significant challenge worldwide because of its poor prognosis and high mortality rate, and accurate prognostic gene signatures are urgently required for individual therapy. This study aimed to construct and validate a seven-gene signature for predicting overall survival (OS) in patients with KIRC. The mRNA expression profile and clinical data of patients with KIRC were obtained from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC). Prognosis-associated genes were identified, and a prognostic gene signature was constructed. Then, the prognostic efficiency of the gene signature was assessed. The results obtained using data from the TCGA were validated using those from the ICGC and other online databases. Gene set enrichment analyses (GSEA) were performed to explore potential molecular mechanisms. A seven-gene signature (PODXL, SLC16A12, ZIC2, ATP2B3, KRT75, C20orf141, and CHGA) was constructed, and it was found to be effective in classifying KIRC patients into high- and low-risk groups, with significantly different survival based on the TCGA and ICGC validation data set. Cox regression analysis revealed that the seven-gene signature had an independent prognostic value. Then, we established a nomogram, including the seven-gene signature, which had a significant clinical net benefit. Interestingly, the seven-gene signature had a good performance in distinguishing KIRC from normal tissues. GSEA revealed that several oncological signatures and GO terms were enriched. This study developed a novel seven-gene signature and nomogram for predicting the OS of patients with KIRC, which may be helpful for clinicians in establishing individualized treatments.     

Characterization and prognostic significance of alternative splicing events in lower‐grade diffuse gliomas

Alternative splicing (AS) is assumed to play important roles in the progression and prognosis of cancer. Currently, the comprehensive analysis and clinical relevance of AS in lower‐grade diffuse gliomas have not been systematically addressed. Here, we gathered alternative splicing data of lower‐grade diffuse gliomas from SpliceSeq. Based on the Percent Spliced In (PSI) values of 515 lower‐grade diffuse glioma patients from the Cancer Genome Atlas (TCGA), we performed subtype‐differential AS analysis and consensus clustering to determine robust clusters of patients. A total of 48 050 AS events in 10 787 genes in lower‐grade diffuse gliomas were profiled. Subtype‐differential splicing analysis and functional annotation revealed that spliced genes were significantly enriched in numerous cancer‐related biological phenotypes and signalling pathways. Consensus clustering using AS events identified three robust clusters of patients with distinguished pathological and prognostic features. Moreover, each cluster was also associated with distinct genomic alterations. Finally, we developed and validated an AS‐related signature with Cox proportional hazards model. The signature, significantly associated with clinical and molecular features, could serve as an independent prognostic factor for lower‐grade diffuse gliomas. Thus, our results indicated that AS events could discriminate molecular subtypes and have prognostic impact in lower‐grade diffuse gliomas.     

Long noncoding RNA DGCR5 involves in tumorigenesis of esophageal squamous cell carcinoma via SRSF1-mediated alternative splicing of Mcl-1

Long noncoding RNAs (lncRNAs) emerge as essential roles in the regulation of alternative splicing (AS) in various malignancies. Serine- and arginine-rich splicing factor 1 (SRSF1)-mediated AS events are the most important molecular hallmarks in cancer. Nevertheless, the biological mechanism underlying tumorigenesis of lncRNAs correlated with SRSF1 in esophageal squamous cell carcinoma (ESCC) remains elusive. In this study, we found that lncRNA DiGeorge syndrome critical region gene 5 (DGCR5) was upregulated in ESCC clinical samples, which associated with poor prognosis. Through RNA interference and overexpression approaches, we confirmed that DGCR5 contributed to promote ESCC cell proliferation, migration, and invasion while inhibited apoptosis in vitro. Mechanistically, DGCR5 could directly bind with SRSF1 to increase its stability and thus stimulate alternative splicing events. Furthermore, we clarified that SRSF1 regulated the aberrant splicing of myeloid cell leukemia-1 (Mcl-1) and initiated a significant Mcl-1L (antiapoptotic) isoform switch, which contributed to the expression of the full length of Mcl-1. Moreover, the cell-derived xenograft (CDX) model was validated that DGCR5 could facilitate the tumorigenesis of ESCC in vivo. Collectively, our findings identified that the key biological role of lncRNA DGCR5 in alternative splicing regulation and emphasized DGCR5 as a potential biomarker and therapeutic target for ESCC.Subject terms: Tumour biomarkers, Apoptosis     

SRSF5 functions as a novel oncogenic splicing factor and is upregulated by oncogene SRSF3 in oral squamous cell carcinoma

Alternative splicing of precursor messenger RNA has been increasingly associated with tumorigenesis. The serine/arginine-rich protein (SR) family plays key roles in the regulation of pre-mRNA alternative splicing. Increasing evidence has demonstrated that the SR protein family is involved in tumorigenesis. However, the functions and mechanisms of SR proteins in tumourigenesis remain largely unknown. In the present study, we discovered that serine/arginine-rich splicing factor 5 (SRSF5) is a novel oncogenic splicing factor that is overexpressed in oral squamous cell carcinoma (OSCC) tissues and cells, being crucial for OSCC cell proliferation and tumor formation. Overexpression of SRSF5 transformed immortal rodent fibroblasts to form tumors in nude mice, while downregulation of SRSF5 in oral squamous cell lines retarded cell growth, cell cycle progression, and tumor growth. The expression of SRSF5 is controlled by an autoregulation mechanism. Serine/arginine-rich splicing factor 3 (SRSF3) has been identified as an oncogene. We found that SRSF5 is a novel target of SRSF3. SRSF3 impairs the autoregulation of SRSF5 and promotes SRSF5 overexpression in cancer cells. Altogether, the present study demonstrated that SRSF5 is a novel oncogene that is upregulated by SRSF3 in OSCC cells.     

Expression Profile of Three Splicing Factors in Pleural Cells Based on the Underlying Etiology and Its Clinical Values in Patients with Pleural Effusion

Splicing factors (SFs) are involved in oncogenesis or immune modulation, the common underlying processes giving rise to pleural effusion (PE). The expression profiles of three SFs (HNRNPA1, SRSF1, and SRSF3) and their clinical values have never been assessed in PE. The three SFs (in pellets of PE) and conventional tumor markers were analyzed using PE samples in patients with PE (N = 336). The sum of higher–molecular weight (Mw) forms of HNRNPA1 (Sum-HMws-HNRNPA1) and SRSF1 (Sum-HMws-SRSF1) and SRSF3 levels were upregulated in malignant PE (MPE) compared to benign PE (BPE); they were highest in cytology-positive MPE, followed by tuberculous PE and parapneumonic PE. Meanwhile, the lowest-Mw HNRNPA1 (LMw-HNRNPA1) and SRSF1 (LMw-SRSF1) levels were not upregulated in MPE. Sum-HMws-HNRNPA1, Sum-HMws-SRSF1, and SRSF3, but neither LMw-HNRNPA1 nor LMw-SRSF1, showed positive correlations with cancer cell percentages in MPE. The detection accuracy for MPE was high in the order of carcinoembryonic antigen (CEA, 85%), Sum-HMws-HNRNPA1 (76%), Sum-HMws-SRSF1 (68%), SRSF3, cytokeratin-19 fragments (CYFRA 21-1), LMw-HNRNPA1, and LMw-SRSF1. Sum-HMws-HNRNPA1 detected more than half of the MPE cases that were undetected by cytology and CEA. Sum-HMws-HNRNPA1, but not other SFs or conventional tumor markers, showed an association with longer overall survival among patients with MPE receiving chemotherapy. Our results demonstrated different levels of the three SFs with their Mw-specific profiles depending on the etiology of PE. We suggest that Sum-HMws-HNRNPA1 is a supplementary diagnostic marker for MPE and a favorable prognostic indicator for patients with MPE receiving chemotherapy.     

Immune-pyroptosis-related genes predict the prognosis of kidney renal clear cell carcinoma

BackgroundKidney renal clear cell carcinoma (KIRC) is a common cancer of the adult urological system. Recent developments in tumor immunology and pyroptosis biology have provided new directions for kidney cancer treatment. Therefore, there is an urgent need to identify potential targets and prognostic biomarkers for the combination of immunotherapy and pyroptosis-targeted therapy.MethodsThe expression of immune-pyroptosis-related differentially expressed genes (IPR-DEGs) between KIRC and healthy tissues was examined using the Gene Expression Omnibus datasets. The GSE168845 dataset was selected for subsequent analyses. Data of 1793 human immune-related genes were downloaded from the ImmPort database (https://www.immport.org./home), while those of 33 pyroptosis-related genes were extracted from previous reviews. The independent prognostic value of IPR-DEGs was determined using differential expression, prognostic, and univariate and multivariate Cox regression analyses. The GSE53757 dataset was used to further verify the GSDMB and PYCARD levels. In our cohorts, the association among DEGs and clinicopathological features and overall survival was analyzed. The least absolute shrinkage and selection operator Cox regression model was established to evaluate the correlation of IPR-DEGs with the immune score, immune checkpoint gene expression, and one-class logistic regression (OCLR) score. KIRC cells and clinical tissue samples were subjected to quantitative real-time polymerase chain reaction to examine the GSDMB and PYCARD mRNA levels. The GSDMB and PYCARD levels in a healthy kidney cell line (HK-2 cells) and two KIRC cell lines (786-O and Caki-1 cells) were verified. The tissue levels of GSDMB and PYCARD were evaluated using immunohistochemical analysis. GSDMB and PYCARD were knocked down in 786-O cells using short-interfering RNA. Cell proliferation was examined using the cell counting kit-8 assay. Cell migration was measured by transwell migration assaysResultsGSDMB and PYCARD were determined to be IPR-DEGs with independent prognostic values. A risk prognostic model based on GSDMB and PYCARD was successfully established. In the GSE53757 dataset, the GSDMB and PYCARD levels in KIRC tissues were significantly higher than those in healthy tissues. The GSDMB and PYCARD expression was related to T stage and OS in our cohort. The GSDMB and PYCARD levels were significantly correlated with the immune score, immune checkpoint gene expression, and OCLR score. The results of experimental studies were consistent with those of bioinformatics analysis. The GSDMB and PYCARD levels in KIRC cells were significantly upregulated when compared with those in healthy kidney cells. Consistently, GSDMB and PYCARD in KIRC tissues were significantly upregulated when compared with those in adjacent healthy kidney tissues. GSDMB and PYCARD knockdown significantly decreased 786-O cell proliferation (p < 0.05). Transwell migration result reflects that silencing GSDMB and PYCARD inhibited 786-O cell migration (p < 0.05) .ConclusionsGSDMB and PYCARD are potential targets and effective prognostic biomarkers for the combination of immunotherapy and pyroptosis-targeted therapy in KIRC.     

Identification of a universal 6-lncRNA prognostic signature for three pathologic subtypes of renal cell carcinoma

Renal cell carcinoma (RCC) is the most common adult renal epithelial cancer susceptible to metastasis and patients with irresectable RCC always have a poor prognosis. Long noncoding RNAs (lncRNAs) have recently been documented as having critical roles in the etiology of RCC. Nevertheless, the prognostic significance of lncRNA-based signature for outcome prediction in patients with RCC has not been well investigated. Therefore, it is essential to identify a lncRNA-based signature for predicting RCC prognosis. In the current study, we comprehensively analyzed the RNA sequencing data of the three main pathological subtypes of RCC (kidney renal clear cell carcinoma [KIRC], kidney renal papillary cell carcinoma [KIRP], and kidney chromophobe carcinoma [KICH]) from The Cancer Genome Atlas (TCGA) database, and identified a 6-lncRNA prognostic signature with the help of a step-wise multivariate Cox regression model. The 6-lncRNA signature stratified the patients into low- and high-risk groups with significantly different prognosis. Multivariate Cox regression analysis showed that predictive value of the 6-lncRNA signature was independent of other clinical or pathological factors in the entire cohort and in each cohort of RCC subtypes. In addition, the three independent prognostic clinical factors (including age, pathologic stage III, and stage IV) was also stratified into low- and high-risk groups basis on the risk score, and the stratification analyses demonstrated that the high-risk score was a poor prognostic factor. In conclusion, these findings indicate that the 6-lncRNA signature is a novel prognostic biomarker for all three subtypes of RCC, and can increase the accuracy of predicting overall survival.