Survival-associated alternative splicing signatures in non-small cell lung cancer

Novel prognostic alternative splicing events in colorectal Cancer: Impact on immune infiltration and therapy response

OBJECTIVE

This study aims to comprehensively analyze alternative splicing (AS) features in colorectal cancer (CRC) using integrative multi-omics and to elucidate their relationship with the CRC immune microenvironment.

METHODS

Transcriptomic data, clinical information, and Percent Spliced In (PSI) values of AS events for CRC patients were obtained from The Cancer Genome Atlas (TCGA) and TCGA SpliceSeq databases. Differentially expressed AS events were identified. Univariate Cox analysis was used to pinpoint prognosis-related AS events. A prognostic risk model was developed and validated using multivariate Cox analysis, patient survival analysis, and the area under the receiver operating characteristic (ROC) curve (AUC). Gene Set Enrichment Analysis (GSEA), immune infiltration, immunotherapy, chemotherapy sensitivity analyses, and regulatory relationships between AS events and splicing factors (SFs) were conducted. Single-cell sequencing was used to study the distribution of key factors. siRNA and overexpression vectors were utilized to silence/overexpress BCAS1 in CRC cells and evaluate their effects on cell growth, migration, and invasion. Furthermore, the interaction between BCAS1 and ANO7 pre-mRNA was investigated using RIP-PCR.

RESULTS

82 prognosis-related AS events were identified in CRC patients. A 15-AS prognostic model was constructed, which correlated with immune cell infiltration and showed differences in immunotherapy and chemotherapy sensitivity. BCAS1 was identified as a potential regulator of the ANO7|58341|AT splicing event in CRC. Single-cell sequencing analysis revealed the distribution of BCAS1 and ANO7 in cancer stem cells. In vitro experiments demonstrated that overexpression of BCAS1 and silencing of ANO7 inhibit the proliferation, migration, and invasion of CRC cells. Moreover, BCAS1 suppresses the progression of CRC by modulating ANO7 alternative splicing.

CONCLUSION

This study provides new insights into the role of alternative splicing in colorectal cancer, particularly the potential regulatory action of BCAS1 on the ANO7|58341|AT splicing event. It also identifies the impact of alternative splicing on the tumor microenvironment and potential implications for immunotherapy, highlighting its relevance for the in-depth study and treatment of CRC.

Alternative Splicing in Glioblastoma and its Clinical Implication in Outcome Prediction

BACKGROUND AND OBJECTIVE

Alternative splicing (AS) offers an important mechanism to form protein polymorphism. A growing body of evidence indicates the correlation between splicing abnormality and carcinoma. Nevertheless, an overall analysis of AS signatures in glioblastoma (GBM) is absent and urgently needed.

METHODS

TCGA SpliceSea data was used to evaluate the AS profiles and further classified into different AS events. The survival analysis was based on these AS events, and AS-related genes were identified and performed with enrichment analysis. At last, the splicing factor-AS regulatory network was established in Cytoscape.

RESULTS

Eight hundred forty-two splicing events were confirmed as prognostic molecular events in GBM. Furthermore, the final prognostic signature constructed by seven AS events gave good result with an area under the curve (AUC) of receiver operating characteristic (ROC) curve up to 0.935 for five years, showing high potency in predicting patients' outcome. We built the splicing regulatory network to show the internal relationship of splicing events in GBM. PC4 and SFRS1 interacting protein 1 (PSIP1) and histone H4 acetylation may play a significant part in the prognosis induced by splicing events.

CONCLUSION

In our study, a high-efficiency prognostic prediction model was built for GBM patients based on AS events, which could become potential prognostic biomarkers for GBM. Meanwhile, PSIP1 may be a critical target for pharmaceutical treatment.

Identification of Key Prognostic Alternative Splicing Events of Costimulatory Molecule-Related Genes in Colon Cancer

OBJECTIVE

This study aimed to explore the key alternative splicing events in costimulatory molecule-related genes in colon cancer and to determine their correlation with prognosis.

METHODS

Gene expression RNA-sequencing data, clinical data, and SpliceSeq data of colon cancer were obtained from The Cancer Genome Atlas. Differentially expressed alternative splicing events in genes were identified, Followed by correlation analysis of genes corresponding to differentially expressed alternative splicing events with costimulatory molecule-related genes. Survival analysis was conducted using differentially expressed alternative splicing events in these genes and a prognostic model was constructed. Functional enrichment, proteinprotein interaction network, and splicing factor analyses were performed.

RESULTS

In total, 6504 differentially expressed alternative splicing events in 3949 genes were identified between tumor and normal tissues. Correlation analysis revealed 3499 differentially expressed alternative splicing events in 2168 costimulatory molecule-related genes. Moreover, 328 differentially expressed alternative splicing events in 288 costimulatory molecule-related genes were associated with overall survival. The prognostic models constructed using these showed considerable power in predicting survival. The ubiquitin A-52 residue ribosomal protein fusion product 1 and ribosomal protein S9 were the hub nodes in the protein-protein interaction network. Furthermore, one splicing factor, splicing factor proline and glutamine-rich, was significantly associated with patient prognosis. Four splicing factor-alternative splicing pairs were obtained from four alternative splicing events in three genes: TBC1 domain family member 8 B, complement factor H, and mitochondrial fission 1.

CONCLUSION

The identified differentially expressed alternative splicing events of costimulatory molecule-related genes may be used to predict patient prognosis and immunotherapy responses in colon cancer.

Novel, alternative splicing signature to detect lymph node metastasis in prostate adenocarcinoma with machine learning

Background

The presence of lymph node metastasis leads to a poor prognosis for prostate cancer (Pca). Recently, many studies have indicated that gene signatures may be able to predict the status of lymph nodes. The purpose of this study is to probe and validate a new tool to predict lymph node metastasis (LNM) based on alternative splicing (AS).

Methods

Gene expression profiles and clinical information of prostate adenocarcinoma cohort were retrieved from The Cancer Genome Atlas (TCGA) database, and the corresponding RNA-seq splicing events profiles were obtained from the TCGA SpliceSeq. Limma package was used to identify the differentially expressed alternative splicing (DEAS) events between LNM and non-LNM groups. Eight machine learning classifiers were built to train with stratified five-fold cross-validation. SHAP values was used to explain the model.

Results

333 differentially expressed alternative splicing (DEAS) events were identified. Using correlation filter and the least absolute shrinkage and selection operator (LASSO) method, a 96 AS signature was identified that had favorable discrimination in the training set and validated in the validation set. The linear discriminant analysis (LDA) was the best classifier after 100 iterations of training. The LDA classifier was able to distinguish between LNM and non-LNM with an area under the receiver operating curve of 0.962 ± 0.026 in the training set (D1 = 351) and 0.953 in the validation set (D2 = 62). The decision curve analysis plot proved the clinical application of the AS-based model.

Conclusion

Machine learning combined with AS data could robustly distinguish between LNM and non-LNM in Pca.

The Role of Alternative Splicing in Cancer: Regulatory Mechanism, Therapeutic Strategy, and Bioinformatics Application

[Formula: see text] Alternative splicing (AS) can generate distinct transcripts and subsequent isoforms that play differential functions from the same pre-mRNA. Recently, increasing numbers of studies have emerged, unmasking the association between AS and cancer. In this review, we arranged AS events that are closely related to cancer progression and presented promising treatments based on AS for cancer therapy. Obtaining proliferative capacity, acquiring invasive properties, gaining angiogenic features, shifting metabolic ability, and getting immune escape inclination are all splicing events involved in biological processes. Spliceosome-targeted and antisense oligonucleotide technologies are two novel strategies that are hopeful in tumor therapy. In addition, bioinformatics applications based on AS were summarized for better prediction and elucidation of regulatory routines mingled in. Together, we aimed to provide a better understanding of complicated AS events associated with cancer biology and reveal AS a promising target of cancer treatment in the future.

Identification of prognostic alternative splicing signature in gastric cancer

Background

Aberrant alternative splicing (AS) events could be viewed as prognostic indicators in a large number of malignancies. This study aims to identify prognostic AS events, illuminate the function of the splicing variants biomarkers and provide reliable evidence for formulating public health strategies for gastric cancer (GC) surveillance.

Methods

RNA-Seq data, clinical information and percent spliced in (PSI) values were available in The cancer genome atlas (TCGA) and TCGA SpliceSeq data portal. A three-step regression method was conducted to identify prognostic AS events and construct multi-AS-based signatures. The associations between prognostic AS events and splicing factors were also investigated.

Results

We identified a total of 1,318 survival-related AS events in GC, parent genes of which were implicated in numerous oncogenic pathways. The final prognostic signatures stratified by seven types of AS events or not stratified performed well in risk prediction for GC patients. Moreover, five signatures based on AA, AD, AT, ES and RI events function as independent prognostic indicators after multivariate adjustment of other clinical variables. Splicing network also showed marked correlation between the expression of splicing factors and PSI value of AS events in GC patients.

Conclusion

Our findings provide a landscape of AS events and regulatory network in GC, indicating that AS events might serve as prognostic biomarkers and therapeutic targets for GC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13690-022-00894-3.

Ferroptosis‑related long non‑coding RNAs and the roles of LASTR in stomach adenocarcinoma

Ferroptosis is a form of programmed cell death that participates in diverse physiological processes. Increasing evidence suggests that long noncoding RNAs (lncRNAs) regulate ferroptosis in tumors, including stomach adenocarcinoma (STAD). In the present study, RNA-sequencing data from The Cancer Genome Atlas database and ferroptosis-related markers from the FerrDb data resource were analyzed to select differentially expressed lncRNAs. Univariate and multivariate Cox regression analyses were performed on these differentially expressed lncRNAs to screen 12 lncRNAs linked with overall survival (OS) and 13 associated with progression-free survival (PFS). Subsequently, two signatures for predicting OS and PFS were established based on these lncRNAs. Kaplan-Meier analyses indicated that the high-risk group of patients with STAD had relatively poor prognosis. The areas under the receiver operating characteristic curves of the two signatures indicated their excellent efficacy in predicting STAD prognosis. In addition, the effect of the lncRNA LASTR on proliferation and migration in gastric cancer was confirmed and the relationship between LASTR and ferroptosis was initially explored through experiments. These results provide potential novel targets for tumor treatment and promote personalized medicine.

A Novel Nomogram Combining Alternative Splicing Events and Clinical Factors for Prognosis Prediction in Head and Neck Squamous Cell Carcinoma

Due to limitations of sensitive biomarkers, the clinical prognosis of patients with head and neck squamous cell carcinoma (HNSCC) remains poor. Alternative splicing (AS) is the basis of both transcriptome and proteome richness, so more and more evidence indicates an important relationship between AS and tumor progression. The aim of this study was to offer a comprehensive analysis on AS events and then investigate its potentials as a new biomarker for patients with squamous cell carcinoma of the head and neck. In this study, univariate assays were conducted to examine the prognosis-associated AS events, and we screened 4068 survival-related AS events in 2573 genes. Then, the AS events related to survival were further determined and analyzed using LASSO regression and multivariate assays, and an eleven-AS signature was developed. Kaplan–Meier assays indicated patients with high-risk scores exhibited a shorter OS than those with low-risk scores. Multivariate assays further demonstrated that the signature’s risk score was independent of HNSCC survivals. Meanwhile, we analyzed the clinical association of AS-based prognostic signature in HNSCC patients and observed that tumor specimens with advanced stages and grades exhibited a high risk score. In addition, the results of survival nomogram revealed that predicted outcomes and actual outcomes were highly consistent. Overall, our group showed an eleven-AS signature of HNSCC, which could be regarded as a separate prognostic factor.

Switched alternative splicing events as attractive features in lung squamous cell carcinoma

Background

Alternative splicing (AS) plays important roles in transcriptome and proteome diversity. Its dysregulation has a close affiliation with oncogenic processes. This study aimed to evaluate AS-based biomarkers by machine learning algorithms for lung squamous cell carcinoma (LUSC) patients.

Method

The Cancer Genome Atlas (TCGA) database and TCGA SpliceSeq database were utilized. After data composition balancing, Boruta feature selection and Spearman correlation analysis were used for differentially expressed AS events. Random forests and a nested fivefold cross-validation were applied for lymph node metastasis (LNM) classifier building. Random survival forest combined with Cox regression model was performed for a prognostic model, based on which a nomogram was developed. Functional enrichment analysis and Spearman correlation analysis were also conducted to explore underlying mechanisms. The expression of some switch-involved AS events along with parent genes was verified by qRT-PCR with 20 pairs of normal and LUSC tissues.

Results

We found 16 pairs of splicing events from same parent genes which were strongly related to the splicing switch (intrapair correlation coefficient = − 1). Next, we built a reliable LNM classifier based on 13 AS events as well as a nice prognostic model, in which switched AS events behaved prominently. The qRT-PCR presented consistent results with previous bioinformatics analysis, and some AS events like ITIH5-10715-AT and QKI-78404-AT showed remarkable detection efficiency for LUSC.

Conclusion

AS events, especially switched ones from the same parent genes, could provide new insights into the molecular diagnosis and therapeutic drug design of LUSC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-02429-2.

A Prediction Model Using Alternative Splicing Events and the Immune Microenvironment Signature in Lung Adenocarcinoma

Background

Alternative splicing (AS) is a gene regulatory mechanism that drives protein diversity. Dysregulation of AS is thought to play an essential role in cancer initiation and development. This study aimed to construct a prognostic signature based on AS and explore the role in the tumor immune microenvironment (TIME) in lung adenocarcinoma.

Methods

We analyzed transcriptome profiling and clinical lung adenocarcinoma data from The Cancer Genome Atlas (TCGA) database and lists of AS-related and immune-related signatures from the SpliceSeq. Prognosis-related AS events were analyzed by univariate Cox regression analysis. Gene set enrichment analyses (GSEA) were performed for functional annotation. Prognostic signatures were identified and validated using univariate and multivariate Cox regression, LASSO regression, Kaplan–Meier survival analyses, and proportional hazards model. The context of TIME in lung adenocarcinoma was also analyzed. Gene and protein expression data of Cyclin-Dependent Kinase Inhibitor 2A (CDKN2A) were obtained from ONCOMINE and Human Protein Atlas. Splicing factor (SF) regulatory networks were visualized.

Results

A total of 19,054 survival-related AS events in lung adenocarcinoma were screened in 1,323 genes. Exon skip (ES) and mutually exclusive exons (ME) exhibited the most and fewest AS events, respectively. Based on AS subtypes, eight AS prognostic signatures were constructed. Patients with high-risk scores were associated with poor overall survival. A nomogram with good validity in prognostic prediction was generated. AUCs of risk scores at 1, 2, and 3 years were 0.775, 0.736, and 0.759, respectively. Furthermore, the prognostic signatures were significantly correlated with TIME diversity and immune checkpoint inhibitor (ICI)-related genes. Low-risk patients had a higher StromalScore, ImmuneScore, and ESTIMATEScore. AS-based risk score signature was positively associated with CD8+ T cells. CDKN2A was also found to be a prognostic factor in lung adenocarcinoma. Finally, potential functions of SFs were determined by regulatory networks.

Conclusion

Taken together, our findings show a clear association between AS and immune cell infiltration events and patient outcome, which could provide a basis for the identification of novel markers and therapeutic targets for lung adenocarcinoma. SF networks provide information of regulatory mechanisms.

The pattern of alternative splicing in lung adenocarcinoma shows novel events correlated with tumorigenesis and immune microenvironment

Lung adenocarcinoma (LUAD) is the leading cause of cancer deaths worldwide due to the lack of early diagnostic markers and specific drugs. Previous studies have shown the association of LUAD growth with aberrant alternative splicing (AS). Herein, clinical data of 535 tumor tissues and 59 normal tissues were extracted from The Cancer Genome Atlas (TCGA) database. Each sample was analyzed using the ESTIMATE algorithm; a comparison between higher and lower score groups (stromal or immune) was made to determine the overall- and progression-free survival-related differentially expressed AS (DEAS) events. We then performed unsupervised clustering of these DEASs, followed by determining their relationship with survival rate, immune cells, and the tumor microenvironment (TME). Next, two prognostic signatures were developed using bioinformatics tools to explore the prognosis of cases with LUAD. Five OS- and six PFS-associated DEAS events were implemented to establish a prognostic risk score model. When compared to the high-risk group (HRG), the PFS and OS of the low-risk group (LRG) were found to be considerable. Additionally, a better prognosis was found considerably associated with the ESTIMATE score of the patients as well as immune cells infiltration. Our analysis of AS events in LUAD not only helps to clarify the tumorigenesis mechanism of AS but also provides ideas for revealing potential prognostic biomarkers and therapeutic targets.

Prognostic Signatures of Alternative Splicing Events in Esophageal Carcinoma Based on TCGA Splice-Seq Data

An alternative splicing (AS) event is a highly complex process that plays an essential role in post-transcriptional gene expression. Several studies have suggested that abnormal AS events were the primary element in the pathological process of cancer. However, few works are dedicated to the study of AS events in esophageal carcinoma (EC). In the present study, clinical information and RNA-seq data of EC patients were downloaded from The Cancer Genome Atlas (TCGA) database. The percent spliced in (PSI) values of AS events were acquired from the TCGA Splice-seq. A total of 183 EC patients were enrolled in this study, and 2,212 AS events were found significantly associated with the overall survival of these patients by univariate Cox regression analysis. The prognostic signatures based on AS events were built by multivariate Cox analysis. Receiver operating characteristic (ROC) curves displayed that the area under the curve (AUC) of the following prognostic signatures, including exon skip (ES), alternate terminator (AT), alternate acceptor site (AA), alternate promoter (AP), alternate donor site (AD), retained intron (RI), and total events, was greater than 0.8, suggesting that these seven signatures had valuable prognosis prediction capacity. Finally, the risk score of prognostic signatures was indicated as an independent risk factor of survival. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to explore the function of splicing factors (SFs) that were associated with AS events. Also, the interactive network between AS events and SFs identified several hub genes and AS events which need further study. This was a comprehensive study that explored prognosis-related AS events and established valuable prognosis signatures in EC patients. The network of interactions between AS events and SFs might offer novel insights into the fundamental mechanisms of tumorigenesis and progression of EC.

Controversial roles of cold‑inducible RNA‑binding protein in human cancer (Review)

Cold‑inducible RNA‑binding protein (CIRBP) is a cold‑shock protein comprised of an RNA‑binding motif that is induced by several stressors, such as cold shock, UV radiation, nutrient deprivation, reactive oxygen species and hypoxia. CIRBP can modulate post‑transcriptional regulation of target mRNA, which is required to control DNA repair, circadian rhythms, cell growth, telomere integrity and cardiac physiology. In addition, the crucial function of CIRBP in various human diseases, including cancers and inflammatory disease, has been reported. Although CIRBP is primarily considered to be an oncogene, it may also serve a role in tumor suppression. In the present study, the controversial roles of CIRBP in various human cancers is summarized, with a focus on the interconnectivity between CIRBP and its target mRNAs involved in tumorigenesis. CIRBP may represent an important prognostic marker and therapeutic target for cancer therapy.

Systemic characterization of alternative splicing related to prognosis and immune infiltration in malignant mesothelioma

BACKGROUND

Malignant mesothelioma (MM) is a relatively rare and highly lethal tumor with few treatment options. Thus, it is important to identify prognostic markers that can help clinicians diagnose mesothelioma earlier and assess disease activity more accurately. Alternative splicing (AS) events have been recognized as critical signatures for tumor diagnosis and treatment in multiple cancers, including MM.

METHODS

We systematically examined the AS events and clinical information of 83 MM samples from TCGA database. Univariate Cox regression analysis was used to identify AS events associated with overall survival. LASSO analyses followed by multivariate Cox regression analyses were conducted to construct the prognostic signatures and assess the accuracy of these prognostic signatures by receiver operating characteristic (ROC) curve and Kaplan-Meier survival analyses. The ImmuCellAI and ssGSEA algorithms were used to assess the degrees of immune cell infiltration in MM samples. The survival-related splicing regulatory network was established based on the correlation between survival-related AS events and splicing factors (SFs).

RESULTS

A total of 3976 AS events associated with overall survival were identified by univariate Cox regression analysis, and ES events accounted for the greatest proportion. We constructed prognostic signatures based on survival-related AS events. The prognostic signatures proved to be an efficient predictor with an area under the curve (AUC) greater than 0.9. Additionally, the risk score based on 6 key AS events proved to be an independent prognostic factor, and a nomogram composed of 6 key AS events was established. We found that the risk score was significantly decreased in patients with the epithelioid subtype. In addition, unsupervised clustering clearly showed that the risk score was associated with immune cell infiltration. The abundances of cytotoxic T (Tc) cells, natural killer (NK) cells and T-helper 17 (Th17) cells were higher in the high-risk group, whereas the abundances of induced regulatory T (iTreg) cells were lower in the high-risk group. Finally, we identified 3 SFs (HSPB1, INTS1 and LUC7L2) that were significantly associated with MM patient survival and then constructed a regulatory network between the 3 SFs and survival-related AS to reveal potential regulatory mechanisms in MM.

CONCLUSION

Our study provided a prognostic signature based on 6 key events, representing a better effective tumor-specific diagnostic and prognostic marker than the TNM staging system. AS events that are correlated with the immune system may be potential therapeutic targets for MM.

Identification of Survival-Associated Gene Signature in Lung Cancer Coexisting With COPD

Background: Chronic obstructive pulmonary disease (COPD) and lung cancer often coexist, which is associated with a worse prognosis. Thousands of biomarkers related to the survival of lung cancer have been investigated. However, those which can predict the survival of lung cancer coexisting with COPD are currently lacking. The present study aimed to identify novel gene signatures to predict the survival of patients with lung cancer coexisting COPD.

Method: RNA-sequence data of lung cancer and control accompanying with matched clinical information were retrieved from the Cancer Genome Atlas (TCGA). Differently expressed genes (DEGs) associated with lung cancer coexisting COPD were screened. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were performed. Univariate and multivariate Cox regression analyses were applied to identify survival-associated DEGs and to construct survival-associated gene signature. Kaplan-Meier survival analysis and calibration plots of the nomogram were performed to test the predictive accuracy of the gene signature. qPCR was performed to validate the genes in the prognostic signature.

Results: Sequence data from 70 patients with lung cancer coexisting COPD, 127 with lung cancer alone and 108 control tissues were included for analysis. A total of 2424 DEGs were identified when comparing lung cancer coexisting COPD with controls. The biological process was primarily associated with DNA-binding transcription activator activity, peptidase inhibitor activity, endopeptidase inhibitor activity, et al. KEGG pathways were mainly enriched in neuroactive ligand-receptor interaction, cell cycle, and Staphylococcus aureus infection. A survival-associated gene signature consisting of CEACAM5, RASAL1, CSTL1, CNGB1, and SLC4A3 was identified and represented as risk score. The high-risk score group had significantly worse survival than the low-risk score group (P < 0.001). Areas under receiver operating characteristic curves were 0.943, 0.773, 0.888 for predicting overall survival at 1-, 3-, and 5-year, respectively. The risk score was an independent predictor of survival, independent of clinical factors. High conformity of the actual survival and the nomogram–predicted probability of survival by applying the risk score. Upregulation of the five genes in patients with lung cancer coexisting COPD were confirmed by qPCR in an independent cohort.

Conclusion: Our study constructed and validated a novel prognostic gene signature for predicting survival of patient with lung cancer coexisting COPD, which may contribute to the clinical treatment decisions.

N6-Methyladenosine RNA Methylation Regulator-Related Alternative Splicing (AS) Gene Signature Predicts Non-Small Cell Lung Cancer Prognosis

Aberrant N6-methyladenosine (m6A) RNA methylation regulatory genes and related gene alternative splicing (AS) could be used to predict the prognosis of non-small cell lung carcinoma. This study focused on 13 m6A regulatory genes (METTL3, METTL14, WTAP, KIAA1429, RBM15, ZC3H13, YTHDC1, YTHDC2, YTHDF1, YTHDF2, HNRNPC, FTO, and ALKBH5) and expression profiles in TCGA-LUAD (n = 504) and TCGA-LUSC (n = 479) datasets from the Cancer Genome Atlas database. The data were downloaded and bioinformatically and statistically analyzed, including the gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses. There were 43,948 mRNA splicing events in lung adenocarcinoma (LUAD) and 46,020 in lung squamous cell carcinoma (LUSC), and the data suggested that m6A regulators could regulate mRNA splicing. Differential HNRNPC and RBM15 expression was associated with overall survival (OS) of LUAD and HNRNPC and METTL3 expression with the OS of LUSC patients. Furthermore, the non-small cell lung cancer prognosis-related AS events signature was constructed and divided patients into high- vs. low-risk groups using seven and 14 AS genes in LUAD and LUSC, respectively. The LUAD risk signature was associated with gender and T, N, and TNM stages, but the LUSC risk signature was not associated with any clinical features. In addition, the risk signature and TNM stage were independent prognostic predictors in LUAD and the risk signature and T stage were independent prognostic predictors in LUSC after the multivariate Cox regression and receiver operating characteristic analyses. In conclusion, this study revealed the AS prognostic signature in the prediction of LUAD and LUSC prognosis.

Exploration of predictive and prognostic alternative splicing signatures in lung adenocarcinoma using machine learning methods

Background

Alternative splicing (AS) plays critical roles in generating protein diversity and complexity. Dysregulation of AS underlies the initiation and progression of tumors. Machine learning approaches have emerged as efficient tools to identify promising biomarkers. It is meaningful to explore pivotal AS events (ASEs) to deepen understanding and improve prognostic assessments of lung adenocarcinoma (LUAD) via machine learning algorithms.

Method

RNA sequencing data and AS data were extracted from The Cancer Genome Atlas (TCGA) database and TCGA SpliceSeq database. Using several machine learning methods, we identified 24 pairs of LUAD-related ASEs implicated in splicing switches and a random forest-based classifiers for identifying lymph node metastasis (LNM) consisting of 12 ASEs. Furthermore, we identified key prognosis-related ASEs and established a 16-ASE-based prognostic model to predict overall survival for LUAD patients using Cox regression model, random survival forest analysis, and forward selection model. Bioinformatics analyses were also applied to identify underlying mechanisms and associated upstream splicing factors (SFs).

Results

Each pair of ASEs was spliced from the same parent gene, and exhibited perfect inverse intrapair correlation (correlation coefficient = − 1). The 12-ASE-based classifier showed robust ability to evaluate LNM status of LUAD patients with the area under the receiver operating characteristic (ROC) curve (AUC) more than 0.7 in fivefold cross-validation. The prognostic model performed well at 1, 3, 5, and 10 years in both the training cohort and internal test cohort. Univariate and multivariate Cox regression indicated the prognostic model could be used as an independent prognostic factor for patients with LUAD. Further analysis revealed correlations between the prognostic model and American Joint Committee on Cancer stage, T stage, N stage, and living status. The splicing network constructed of survival-related SFs and ASEs depicts regulatory relationships between them.

Conclusion

In summary, our study provides insight into LUAD researches and managements based on these AS biomarkers.

A SNP involved in alternative splicing of ABCB1 is associated with clopidogrel resistance in coronary heart disease in Chinese population

Although many scientists are studying the association between genetic polymorphism of ABCB1 and CR in patients, the molecular mechanism has not been further studied in patients with CHD. This study investigated the relationship between SNP of the ABCB1 gene in patients with CHD and CR, and whether the polymorphism of the ABCB1 gene affects the AS of the gene. 741 patients were enrolled in the study, 316 CR cases and 425 NCR cases. The correlation between CR risk and clinical-pathological characteristics were studied. Additionally, the five SNPs were analysed by PCR and Mass Array genotyping methods. Furthermore, silicon analysis was used to predict whether the polymorphism affects the process of AS. Results showed that there was a significant correlation between rs1045642 polymorphism and CR in genotyping and allele analysis. The rs1045642 polymorphism of the ABCB1 gene of CHD patients carrying the A allele are more likely to develop CR. Silicon analysis showed that rs1045642 generated a new ESE sequence which might affect AS of ABCB1 gene. We hypothesize that the mechanism of CR might be caused by a change in the AS caused by the polymorphism of the gene. Thus, this work provides guidance for the clinical use of clopidogrel.

QKI-5 regulates the alternative splicing of cytoskeletal gene ADD3 in lung cancer

Accumulating evidence indicates that the alternative splicing program undergoes extensive changes during cancer development and progression. The RNA-binding protein QKI-5 is frequently downregulated and exhibits anti-tumor activity in lung cancer. Howeve-r, little is known about the functional targets and regulatory mechanism of QKI-5. Here, we report that upregulation of exon 14 inclusion of cytoskeletal gene Adducin 3 (ADD3) significantly correlates with a poor prognosis in lung cancer. QKI-5 inhibits cell proliferation and migration in part through suppressing the splicing of ADD3 exon 14. Through genome-wide mapping of QKI-5 binding sites in vivo at nucleotide resolution by iCLIP-seq analysis, we found that QKI-5 regulates alternative splicing of its target mRNAs in a binding position-dependent manner. By binding to multiple sites in an upstream intron region, QKI-5 represses the splicing of ADD3 exon 14. We also identified several QKI mutations in tumors, which cause dysregulation of the splicing of QKI targets ADD3 and NUMB. Taken together, our results reveal that QKI-mediated alternative splicing of ADD3 is a key lung cancer-associated splicing event, which underlies in part the tumor suppressor function of QKI.

Genome-Wide Analysis of Prognostic Alternative Splicing Signature and Splicing Factors in Lung Adenocarcinoma

Analysis of The Cancer Genome Atlas data revealed that alternative splicing (AS) events could serve as prognostic biomarkers in various cancer types. This study examined lung adenocarcinoma (LUAD) tissues for AS and assessed AS events as potential indicators of prognosis in our cohort. RNA sequencing and bioinformatics analysis were performed. We used SUPPA2 to analyze the AS profiles. Using univariate Cox regression analysis, overall survival (OS)-related AS events were identified. Genes relating to the OS-related AS events were imported into Cytoscape, and the CytoHubba application was run. OS-related splicing factors (SFs) were explored using the log-rank test. The relationship between the percent spliced-in value of the OS-related AS events and SF expression was identified by Spearman correlation analysis. We found 1957 OS-related AS events in 1151 genes, and most were protective factors. Alternative first exon splicing was the most frequent type of splicing event. The hub genes in the gene network of the OS-related AS events were FBXW11, FBXL5, KCTD7, UBB and CDC27. The area under the curve of the MIX prediction model was 0.847 for 5-year survival based on seven OS-related AS events. Overexpression of SFs CELF2 and SRSF5 was associated with better OS. We constructed a correlation network between SFs and OS-related AS events. In conclusion, we identified prognostic predictors using AS events that stratified LUAD patients into high- and low-risk groups. The discovery of the splicing networks in this study provides an insight into the underlying mechanisms.

Identification of Survival-Associated Alternative Splicing Signatures in Lung Squamous Cell Carcinoma

Purpose: Alternative splicing (AS) is a post-transcriptional process that plays a significant role in enhancing the diversity of transcription and protein. Accumulating evidences have demonstrated that dysregulation of AS is associated with oncogenic processes. However, AS signature specifically in lung squamous cell carcinoma (LUSC) remains unknown. This study aimed to evaluate the prognostic values of AS events in LUSC patients.

Methods: The RNA-seq data, AS events data and corresponding clinical information were obtained from The Cancer Genome Atlas (TCGA) database. Univariate Cox regression analysis was performed to identify survival-related AS events and survival-related parent genes were subjected to Gene Ontology enrichment analysis and gene network analysis. The least absolute shrinkage and selection operator (LASSO) method and multivariate Cox regression analysis were used to construct prognostic prediction models, and their predictive values were assessed by Kaplan-Meier analysis and receiver operating characteristic (ROC) curves. Then a nomogram was established to predict the survival of LUSC patients. And the interaction network of splicing factors (SFs) and survival-related AS events was constructed by Spearman correlation analysis and visualized by Cytoscape.

Results: Totally, 467 LUSC patients were included in this study and 1,991 survival-related AS events within 1,433 genes were identified. SMAD4, FOS, POLR2L, and RNPS1 were the hub genes in the gene interaction network. Eight prognostic prediction models (seven types of AS and all AS) were constructed and all exhibited high efficiency in distinguishing good or poor survival of LUSC patients. The final integrated prediction model including all types of AS events exhibited the best prognostic power with the maximum AUC values of 0.778, 0.816, 0.814 in 1, 3, 5 years ROC curves, respectively. Meanwhile, the nomogram performed well in predicting the 1-, 3-, and 5-year survival of LUSC patients. In addition, the SF-AS regulatory network uncovered a significant correlation between SFs and survival-related AS events.

Conclusion: This is the first comprehensive study to analyze the role of AS events in LUSC specifically, which improves our understanding of the prognostic value of survival-related AS events for LUSC. And these survival-related AS events might serve as novel prognostic biomarkers and drug therapeutic targets for LUSC.

The Construction of Bone Metastasis-Specific Prognostic Model and Co-expressed Network of Alternative Splicing in Breast Cancer

Background

Breast cancer (BRCA) ranks among the top most common female malignancies and was regarded as incurable when combined with bone and distant metastasis. Alternative splicing events (ASEs) together with splicing factors (SFs) were considered responsible for the development and progression of tumors.

Methods

Datasets including RNA sequencing and ASEs of BRCA samples were achieved from TCGA and TCGASpliceSeq databases. Then, a survival model was built including 15 overall-survival-associated splicing events (OS-SEs) by Cox regression and Lasso regression. The co-expressed SFs of each bone-and-distant-metastasis-related OS-SE were discovered by Pearson correlation analysis. Additionally, Gene Set Variation Analysis (GSVA) was performed to identify the downstream mechanisms of the key OS-SEs. Finally, the results were validated in different online platforms.

Results

A reliable survival model was established (the area under ROC = 0.856), and CIRBP was found co-expressed with FAM110B (R = 0.320, P < 0.001) associated with the fatty acid metabolism pathway.

Conclusion

Aberrant SF, CIRBP, regulated a specific ASE, exon skip (ES) of FAM110B, during which the fatty acid metabolism pathway played an essential part in tumorigenesis and prognosis of BRCA.