Dataset Information

Alternative splicing associated with cancer stemness in kidney renal clear cell carcinoma.

ABSTRACT:

Backgroud

Cancer stemness is associated with metastases in kidney renal clear cell carcinoma (KIRC) and negatively correlates with immune infiltrates. Recent stemness evaluation methods based on the absolute expression have been proposed to reveal the relationship between stemness and cancer. However, we found that existing methods do not perform well in assessing the stemness of KIRC patients, and they overlooked the impact of alternative splicing. Alternative splicing not only progresses during the differentiation of stem cells, but also changes during the acquisition of the stemness features of cancer stem cells. There is an urgent need for a new method to predict KIRC-specific stemness more accurately, so as to provide help in selecting treatment options.

Methods

The corresponding RNA-Seq data were obtained from the The Cancer Genome Atlas (TCGA) data portal. We also downloaded stem cell RNA sequence data from the Progenitor Cell Biology Consortium (PCBC) Synapse Portal. Independent validation sets with large sample size and common clinic pathological characteristics were obtained from the Gene Expression Omnibus (GEO) database. we constructed a KIRC-specific stemness prediction model using an algorithm called one-class logistic regression based on the expression and alternative splicing data to predict stemness indices of KIRC patients, and the model was externally validated. We identify stemness-associated alternative splicing events (SASEs) by analyzing different alternative splicing event between high- and low- stemness groups. Univariate Cox and multivariable logistic regression analysisw as carried out to detect the prognosis-related SASEs respectively. The area under curve (AUC) of receiver operating characteristic (ROC) was performed to evaluate the predictive values of our model.

Results

Here, we constructed a KIRC-specific stemness prediction model with an AUC of 0.968,and to provide a user-friendly interface of our model for KIRC stemness analysis, we have developed KIRC Stemness Calculator and Visualization (KSCV), hosted on the Shiny server, can most easily be accessed via web browser and the url https://jiang-lab.shinyapps.io/kscv/ . When applied to 605 KIRC patients, our stemness indices had a higher correlation with the gender, smoking history and metastasis of the patients than the previous stemness indices, and revealed intratumor heterogeneity at the stemness level. We identified 77 novel SASEs by dividing patients into high- and low- stemness groups with significantly different outcome and they had significant correlations with expression of 17 experimentally validated splicing factors. Both univariate and multivariate survival analysis demonstrated that SASEs closely correlated with the overall survival of patients.

Conclusions

Basing on the stemness indices, we found that not only immune infiltration but also alternative splicing events showed significant different at the stemness level. More importantly, we highlight the critical role of these differential alternative splicing events in poor prognosis, and we believe in the potential for their further translation into targets for immunotherapy.

SUBMITTER: Xiao L

PROVIDER: S-EPMC8204412 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Background: Kidney renal clear cell carcinoma (KIRC) is the malignancy originated from the renal epithelium, with a high rate of distant metastasis. Aberrant alternative splicing (AS) of pre-mRNA are widely reported to be involved in the tumorigenesis and metastasis of multiple cancers. The aim of this study is to explore the mechanism of alternative splicing events (ASEs) underlying tumorigenesis and metastasis of KIRC. Methods: RNA-seq of 537 KIRC samples downloaded from the TCGA database and ASEs data from the TCGASpliceSeq database were used to identify ASEs in patients with KIRC. The univariate and Lasso regression analysis were used to screen the most significant overall survival-related ASEs (OS-SEs). Based on those, the OS-SEs model was proposed. The interaction network of OS-SEs and splicing factors (SFs) with absolute value of correlation coefficient value >0.750 was constructed by Pearson correlation analysis. The OS-SEs significantly related to distant metastasis and clinical stage were identified by non-parametric test, and those were also integrated into co-expression analysis with prognosis-related Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways identified by Gene Set Variation Analysis (GSVA). ASEs with significance were selected for multiple online database validation. Results: A total of prognostic 6,081 overall survival-related ASEs (OS-SEs) were identified by univariate Cox regression analysis and a prediction model was constructed based on 5 OS-SEs screened by Lasso regression with the Area Under Curve of 0.788. Its risk score was also illustrated to be an independent predictor, which the good reliability of the model. Among 390 identified candidate SFs, DExD-Box Helicase 39B (DDX39B) was significantly correlated with OS and metastasis. After external database validation, Retained Intron of Ras Homolog Family Member T2 (RHOT2) and T-Cell Immune Regulator 1 (TCIRG1) were identified. In the co-expression analysis, overlapped co-expression signal pathways for RHOT2 and TCIRG1 were sphingolipid metabolism and N-glycan biosynthesis. Conclusions: Based on the results of comprehensive bioinformatic analysis, we proposed that aberrant DDX39B regulated RHOT2-32938-RI and TCIRG1-17288-RI might be associated with the tumorigenesis, metastasis, and poor prognosis of KIRC via sphingolipid metabolism or N-glycan biosynthesis pathway.

Project description:BackgroundRenal cell carcinoma (RCC) accounts for 90% of renal cancers, of which clear cell carcinoma (ccRCC) is the most usual histological type. The process of alternative splicing (AS) contributes to protein diversity, and the dysregulation of protein diversity may have a great influence on tumorigenesis. We developed a prognostic signature and comprehensively analyzed the role of tumor immune microenvironment (TIME) and immune checkpoint blocking (ICB) treatment in ccRCC.MethodsTo identify prognosis-related AS events, univariate Cox regression was used and functional annotation was performed using gene set enrichment analysis (GSEA). In this study, prognostic signatures were developed based on multivariate Cox, univariate Cox, and LASSO regression models. Moreover, to assess the prognostic value, the proportional hazards model, Kruskal-Wallis analysis, and ROC curves were used. To obtain a better understanding of TIME in ccRCC, the ESTIMATE R package, single sample gene set enrichment analysis (ssGSEA) algorithm, CIBERSORT method, and the tumor immune estimation resource (TIMER) were applied. The database was searched to verify the expression of C4OF19 in tumor and normal samples. Regulatory networks for AS-splicing factors (SFs) were visualized using Cytoscape 3.9.1.ResultsThere were 9,347 AS cases associated with the survival of ccRCC patients screened. A total of eight AS prognostic signatures were developed with stable prognostic predictive accuracy based on splicing subtypes. In addition, a qualitative prognostic nomogram was developed, and the prognostic prediction showed high effectiveness. In addition, we found that the combined signature was significantly associated with the diversity of TIME and ICB treatment-related genes. C4ORF19 might become an important prognostic factor for ccRCC. Finally, the AS-SF regulatory network was established to clearly reveal the potential function of SFs.ConclusionWe found novel and robust indicators (i.e., risk signature, prognostic nomogram, etc.) for the prognostic prediction of ccRCC. A new and reliable prognostic nomogram was established to quantitatively predict the clinical outcome. The AS-SF networks could provide a new way for the study of potential regulatory mechanisms, and the important roles of AS events in the context of TIME and immunotherapy efficiency were exhibited. C4ORF19 was found to be a vital gene in TIME and ICB treatment.

Project description:BackgroundAberrant alternative splicing events play critical roles in carcinogenesis and progression of many cancers, while sparse studies regarding to alternative splicing are available for clear cell renal cell carcinoma (ccRCC). We identified that alternative splicing of coiled-coil domain containing 50 (CCDC50) was dysregulated in ccRCC, whereas the clinical significance of this splicing event and its splicing regulation mechanisms were still elusive.MethodsBioinformatic algorithm was utilized to identify significant exon skipping events in ccRCC via exon sequencing data from The Cancer Genome Atlas. Semi-quantitative real-time polymerase chain reaction and western blot were used to validate the aberrant expression of different transcripts in renal cancer tissues, cell lines and corresponding noncancerous controls. Short hairpin RNA targeting CCDC50 and overexpressing plasmids for each transcript were introduced into ccRCC cell lines, followed by a series of in vitro and in vivo functional experiments. Moreover, a panel of splicing factors were identified and their roles on splicing regulation of CCDC50 precursor mRNA (pre-mRNA) were studied. Furthermore, RNAseq data were analyzed to elucidate downstream molecules of CCDC50. Two-way analysis of variance and unpaired Student t test were used in statistical analysis.ResultsPre-mRNA of CCDC50 generated two transcripts, full-length transcript (CCDC50-FL) and truncated transcript (CCDC50-S) with exon 6 skipped. CCDC50-S was overexpressed in ccRCC tissues and cell lines compared to noncancerous counterparts, but CCDC50-FL was only detected in noncancerous tissues and normal renal epithelial cells. Higher percent spliced-in index was associated with better survival in ccRCC patients. In vitro and in vivo functional experiments indicated that CCDC50-S transcript promoted the proliferation, migration, invasion and tumorigenesis of ccRCC, while CCDC50-FL exerted opposite tumor suppressive functions. Besides, we identified that heterogeneous nuclear ribonucleoprotein A1 (HnRNP A1) could promote the skipping of exon 6, which resulted in higher portion of CCDC50-S and oncogenic transformation. Moreover, zinc finger protein 395 (ZNF395) was identified as a downstream protein of CCDC50-S, and the interaction initiated oncogenic pathways which were involved in ccRCC progression.ConclusionsAberrant alternative splicing of CCDC50 is regulated by HnRNP A1 in ccRCC. This splicing event contributes to cancer progression through the downstream pathway involving ZNF395.