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.     

Systematic analysis of survival-associated alternative splicing signatures in clear cell renal cell carcinoma

Alternative splicing (AS) constitutes a major reason for messenger RNA (mRNA) and protein diversity. Increasing studies have shown a link to splicing dysfunction associated with malignant neoplasia. Systematic analysis of AS events in kidney cancer remains poorly reported. Therefore, we generated AS profiles in 533 kidney renal clear cell carcinoma (KIRC) patients in The Cancer Genome Atlas (TCGA) database using RNA-seq data. Then, prognostic models were developed in a primary cohort (N = 351) and validated in a validation cohort (N = 182). In addition, splicing networks were built by integrating bioinformatics analyses. A total of 11 268 and 8083 AS variants were significantly associated with patient overall survival time in the primary and validation KIRC cohorts, respectively, including STAT1, DAZAP1, IDS, NUDT7, and KLHDC4. The AS events in the primary KIRC cohorts served as candidate AS events to screen the independent risk factors associated with survival in the primary cohort and to develop prognostic models. The area under the curve of the receiver-operator characteristic curve for prognostic prediction in the primary and validation KIRC cohorts was 0.84 and 0.82 at 2500 days of overall survival, respectively. In addition, splicing correlation networks revealed key splicing factors (SFs) in KIRC, such as HNRNPH1, HNRNPU, KHDBS1, KHDBS3, SRSF9, RBMX, SFQ, SRP54, HNRNPA0, and SRSF6. In this study, we analyzed the AS landscape in the TCGA KIRC cohort and detected predictors (prognostic) based on AS variants with high performance for risk stratification of the KIRC cohort and revealed key SFs in splicing networks, which could act as underlying mechanisms.     

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.     

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.     

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.     

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.     

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.

     

Global genome splicing analysis reveals an increased number of alternatively spliced genes with aging

Alternative splicing (AS) is a key regulatory mechanism for the development of different tissues; however, not much is known about changes to alternative splicing during aging. Splicing events may become more frequent and widespread genome‐wide as tissues age and the splicing machinery stringency decreases. Using skin, skeletal muscle, bone, thymus, and white adipose tissue from wild‐type C57BL6/J male mice (4 and 18 months old), we examined the effect of age on splicing by AS analysis of the differential exon usage of the genome. The results identified a considerable number of AS genes in skeletal muscle, thymus, bone, and white adipose tissue between the different age groups (ranging from 27 to 246 AS genes corresponding to 0.3–3.2% of the total number of genes analyzed). For skin, skeletal muscle, and bone, we included a later age group (28 months old) that showed that the number of alternatively spliced genes increased with age in all three tissues (P < 0.01). Analysis of alternatively spliced genes across all tissues by gene ontology and pathway analysis identified 158 genes involved in RNA processing. Additional analysis of AS in a mouse model for the premature aging disease Hutchinson–Gilford progeria syndrome was performed. The results show that expression of the mutant protein, progerin, is associated with an impaired developmental splicing. As progerin accumulates, the number of genes with AS increases compared to in wild‐type skin. Our results indicate the existence of a mechanism for increased AS during aging in several tissues, emphasizing that AS has a more important role in the aging process than previously known.     

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.     

Identification of prognostic alternative splicing signatures in hepatitis B or/and C viruses related hepatocellular carcinoma

BackgroundAlternative splicing (AS) takes a crucial part in tumor process. We aim to analyze AS in Hepatitis B virus (HBV) or/and hepatitis C virus (HCV) related hepatocellular carcinoma (HCC).MethodsCox regression analysis was conducted to screen survival-associated AS events. The receiver operating characteristic curve used to evaluate the predictive accuracy. Splicing network was built to investigate the relationship between splicing factors and AS events.ResultsNinety-six survival-associated AS events were obtained by univariate Cox regression. Final prognostic model could significantly distinguish the prognosis. We identified RBFOX2 as the hub gene in splicing network based on differentially expressed splicing factors, and obtained MAP3K13_AT as the key AS event in survival-related splicing network.ConclusionOur results highlight the AS signatures in HCC patients with HBV or/and HCV infection. Meanwhile, AS events and splicing factors in different virus-infected HCC subgroups can provide novel perspectives as biomarkers and individualized therapeutic targets.     

Evaluation of serum exosomal LncRNA‐based biomarker panel for diagnosis and recurrence prediction of bladder cancer

Exosomes are small membrane vesicles released by many cells. These vesicles can mediate cellular communications by transmitting active molecules including long non‐coding RNAs (lncRNAs). In this study, our aim was to identify a panel of lncRNAs in serum exosomes for the diagnosis and recurrence prediction of bladder cancer (BC). The expressions of 11 candidate lncRNAs in exosome were investigated in training set (n = 200) and an independent validation set (n = 320) via quantitative real‐time PCR. A three‐lncRNA panel (PCAT‐1, UBC1 and SNHG16) was finally identified by multivariate logistic regression model to provide high diagnostic accuracy for BC with an area under the receiver‐operating characteristic curve (AUC) of 0.857 and 0.826 in training set and validation set, respectively, which was significantly higher than that of urine cytology. The corresponding AUCs of this panel for patients with Ta, T1 and T2‐T4 were 0.760, 0.827 and 0.878, respectively. In addition, Kaplan‐Meier analysis showed that non‐muscle‐invasive BC (NMIBC) patients with high UBC1 expression had significantly lower recurrence‐free survival (P = 0.01). Multivariate Cox analysis demonstrated that UBC1 was independently associated with tumour recurrence of NMIBC (P = 0.018). Our study suggested that lncRNAs in serum exosomes may serve as considerable diagnostic and prognostic biomarkers of BC.     

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.     

Exploring Splicing Variants and Novel Genes in Sacred Lotus Based on RNA-seq Data

Sacred lotus (Nelumbo nucifera) is a typical aquatic plant, belonging to basal eudicot plant, which is ideal for genome and genetic evolutionary study. Understanding lotus gene diversity is important for the study of molecular genetics and breeding. In this research, public RNA-seq data and the annotated reference genome were used to identify the genes in lotus. A total of 26,819 consensus and 1,081 novel genes were identified. Meanwhile, a comprehensive analysis of gene alternative splicing events was conducted, and a total of 19,983 “internal” alternative splicing (AS) events and 14,070 “complete” AS events were detected in 5,878 and 5,881 multi-exon expression genes, respectively. Observations made from the AS events show the predominance of intron retention (IR) subtype of AS events representing 33%. IR is followed by alternative acceptor (AltA), alternative donor (AltD) and exon skipping (ES), highlighting the universality of the intron definition model in plants. In addition, functional annotations of the gene with AS indicated its relationship to a number of biological processes such as cellular process and metabolic process, showing the key role for alternative splicing in influencing the growth and development of lotus. The results contribute to a better understanding of the current gene diversity in lotus, and provide an abundant resource for future functional genome analysis in lotus.     

Proteogenomic analysis reveals alternative splicing and translation as part of the abscisic acid response in Arabidopsis seedlings

In eukaryotes, mechanisms such as alternative splicing (AS) and alternative translation initiation (ATI) contribute to organismal protein diversity. Specifically, splicing factors play crucial roles in responses to environment and development cues; however, the underlying mechanisms are not well investigated in plants. Here, we report the parallel employment of short‐read RNA sequencing, single molecule long‐read sequencing and proteomic identification to unravel AS isoforms and previously unannotated proteins in response to abscisic acid (ABA) treatment. Combining the data from the two sequencing methods, approximately 83.4% of intron‐containing genes were alternatively spliced. Two AS types, which are referred to as alternative first exon (AFE) and alternative last exon (ALE), were more abundant than intron retention (IR); however, by contrast to AS events detected under normal conditions, differentially expressed AS isoforms were more likely to be translated. ABA extensively affects the AS pattern, indicated by the increasing number of non‐conventional splicing sites. This work also identified thousands of unannotated peptides and proteins by ATI based on mass spectrometry and a virtual peptide library deduced from both strands of coding regions within the Arabidopsis genome. The results enhance our understanding of AS and alternative translation mechanisms under normal conditions, and in response to ABA treatment.