- Research article
- Open Access
Characterization of alternative splicing events and prognostic signatures in breast cancer
BMC Cancer volume 21, Article number: 587 (2021)
Breast cancer (BRCA) is one of the most common cancers worldwide. Abnormal alternative splicing (AS) frequently observed in cancers. This study aims to demonstrate AS events and signatures that might serve as prognostic indicators for BRCA.
Original data for all seven types of splice events were obtained from TCGA SpliceSeq database. RNA-seq and clinical data of BRCA cohorts were downloaded from TCGA database. Survival-associated AS events in BRCA were analyzed by univariate COX proportional hazards regression model. Prognostic signatures were constructed for prognosis prediction in patients with BRCA based on survival-associated AS events. Pearson correlation analysis was performed to measure the correlation between the expression of splicing factors (SFs) and the percent spliced in (PSI) values of AS events. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were conducted to demonstrate pathways in which survival-associated AS event is enriched.
A total of 45,421 AS events in 21,232 genes were identified. Among them, 1121 AS events in 931 genes significantly correlated with survival for BRCA. The established AS prognostic signatures of seven types could accurately predict BRCA prognosis. The comprehensive AS signature could serve as independent prognostic factor for BRCA. A SF-AS regulatory network was therefore established based on the correlation between the expression levels of SFs and PSI values of AS events.
This study revealed survival-associated AS events and signatures that may help predict the survival outcomes of patients with BRCA. Additionally, the constructed SF-AS networks in BRCA can reveal the underlying regulatory mechanisms in BRCA.
BRCA ranks among the top most common female malignancies in China and worldwide [1, 2]. BRCA is treated with a multidisciplinary approach including surgical resection, chemotherapy, hormonotherapy, molecular targeting treatments, and radiotherapy . In the past three decades, breast cancer death rates decreased by 39%, which translates to more than 300,000 averted breast cancer deaths in the United States . However, the heterogeneity and complexity of BRCA still led to poor prognosis in patients with certain types of BRCA. Thus, it is urgent to identify the potential molecular mechanisms underlying BRCA and to improve the prognosis of BRCA patients.
Alternative splicing (AS) is an important post-transcriptional process through which multiple transcripts are generated from a single gene . Dysregulation of AS is known to be implicated in multiple human malignancies [6, 7]. There are seven types of AS events including exon skip (ES), alternate donor site (AD), alternate acceptor site (AA), mutually exclusive exon (ME), alternate terminator (AT), alternate promoter (AP), and retained intron (RI). Increasing evidence indicates that AS is related to cancer development and metastasis [8,9,10,11]. In addition, AS events were reported to serve as prognostic predictors for multiple types of cancer .
Advances in the next-generation sequencing technologies have pushed forward cancer epigenetic researches. SpliceSeq  database provides AS profiles across 33 tumors and enables researchers to study the global profiling of AS events in most predominant human malignant tumors.
Our study aimed to identify AS events in BRCA and its prognostic significance in BRCA patients using data downloaded from TCGA database. AS events that significantly associated with prognosis of BRCA were identified, and AS prognostic signatures based on AS events were then generated. Moreover, a SF-AS regulatory network was also established to reveal the underlying correlations among SFs and AS in BRCA.
AS data download and process
The percent spliced in (PSI) value was calculated and processed by SpliceSeq to quantify AS events. The PSI value indicates the inclusion of a transcript element divided by the total number of reads for that AS event. PSI values range from 0 to 100%, which indicates a shift percentage in AS events. AS events with PSI value of larger than 75% in were downloaded from SpliceSeq database. An UpSet plot was created to show the intersection and distribution among all types of AS. Then clinical data of BRCA patients were obtained from TCGA-BRCA database. The primary endpoint in our study is overall survival (OS).
Patients with an OS of less than 90 days were excluded. The follow-up periods of the BRCA cohort ranged from 91 to 8605 days. Univariate proportional hazards regression model was performed to evaluate the correlation between the PSI value of each AS event and prognosis of BRCA patients.
Construction of prognostic signatures
A least absolute shrinkage and selection operator (LASSO) analysis to develop prognostic signatures based on the top 20 most significant AS events selected from the univariate Cox analysis. The coefficients and partial likelihood deviance were also calculated in LASSO analysis. The prognostic AS signatures were generated by multiplying the PSI values of AS events with the coefficients assigned by LASSO analysis. Then AS prognostic signatures along with multiple clinical parameters were included into multivariable Cox regression analysis to identify independent predictors. A time-dependent receiver-operator characteristic (ROC) curve was used to evaluate the prognostic prediction efficacy of the AS signatures. The risk score of each AS event was calculated to assess the performance of the prognostic signatures. Kaplan-Meier survival analysis with log-rank test was performed to evaluate the prognostic difference between high- and low-risk groups.
SF-AS regulatory network
A list of four hundred and four SFs was obtained from the study of Seiler et al. (Table S3) . The expression data of SFs were downloaded from TCGA-BRCA dataset. Correlations among the expression of SFs and the PSI values of prognostic AS events were assessed by Pearson’s correlation analysis. SFs-AS interactomes with a correlation coefficient more than 0.5 and P value less than 0.05 were selected to create the SF-AS network by Cytoscape software.
KEGG and GO analysis were conducted to investigate the functional categories of genes with prognostic AS events by the package “clusterProfiler” in R software.
Landscape of AS events in BRCA
TCGA splice-seq data and clinical information in TCGA-BRCA project were downloaded, and 1098 patients were enrolled in this study. In total, 45,421 AS events in 21,232 genes were identified in BRCA, suggesting that AS is a common biological process in BRCA. To be specific, 3731 AA events in 2628 genes, 3246 AD events in 2278 genes, 9112 AP events in 3654 genes, 8595 AT events in 3755 genes, 17,702 ES events in 6812 genes, 233 ME events in 227 genes, and 2802 RI events in 1878 genes were observed in preliminary analysis (Fig. 1a). Figure 1b showed the distribution and intersection among seven types of AS. ES was the predominant type, and one gene may have multiple types of AS.
Prognostic AS events
The differential AS events between normal and BRCA tissues were identified using “limma” package in R software and demonstrated with a volcano plot. However, no significant difference of AS event was detected (Table S1). In addition, the differential expressed genes between normal and BRCA tissues were analyzed. In comparison with normal tissue, there are 57 genes upregulated and 215 genes downregulated in BRCA tissues (Table S6). A univariate Cox analysis was performed to evaluate the prognostic significance of AS events in BRCA patients. Results suggested that 1121 AS events in 931 genes were significantly correlated with the prognosis of patients with BRCA (Table S2). Specifically, 88 AA events in 85 genes, 83 AD events in 82 genes, 281 AP events in 190 genes, 62 AT events in 42 genes, 526 ES events in 453 genes, 8 ME events in 7 genes, and 74 RI events in 71 genes were identified as prognostic AS events (Fig. 1c). Additionally, one gene could present multiple AS events that were significantly associated with the OS of patients with BRCA. Figure 1d showed that ES was also the predominant AS type. To verify if a particular splicing-pattern is more enriched of prognostic isoforms, we calculate the enrichment ratio for each type of AS by divide the number of prognostic AS events by the number of overall AS events of each AS type. The ratios of AA, AD, AP, AT, ES, ME, and RI are 2.36, 2.56, 3.08, 0.72, 2.97, 3.43, and 2.64% (Table S7). With the exception of AT, which is lowest (0.72%), there was no significant difference among the ratios of the rest of 6 types of AS (2.36-3.43%). ES does not have a significantly higher ratio than other types of AS.
Prognostic AS signatures
Figure 2a-g showed the 20 most significant prognostic AS events of each of the seven types. Because ME type accounts for the smallest percentage of AS, the only seven ME events were significant associated with survival. Seven types of prognostic signatures based on prognostic AA, AD, AP, AT, ES, ME, and RI events were developed using LASSO Cox analysis (Fig. 3a-g). Moreover, an integrated analysis of all the seven types of AS events was performed to create a comprehensive prognostic signature (abbreviated as “ALL”), which consist of PARPBP-24031-ES, NCOR1–39424-ES, COPZ1–22159-RI, ANK3–11845-AP, ITGB5–100223-ES, PHTF1–4284-RI, HSPBP1–52052-AP, TCF12–30783-AP, RPS6KA1–1282-AP, CNST-10497-ES, TMEM25–19023-AA, TMEM25–19017-AA and BTN3A2–75,630-ES (Fig. 3h, Table 1). Kaplan-Meier survival analysis showed that the eight signatures could effectively separate the survival curves of low-risk groups from those of the high-risk groups (Fig. 4a-h). Figure 5a-h showed the risk scores of eight signatures which ranked from low to high (upper panel). The median was used as a cut-off to divide high- and low-risk groups. Patients with a low-risk score had longer survival time (lower panel). Next, the efficacy of these eight prognostic signatures in prognosis prediction was evaluated by ROC curves. The area under the cure (AUC) of eight signatures was larger than 0.6. It is worth noting that the AUC of comprehensive signatures reached 0.801 (Fig. 6a). Univariate Cox regression analysis showed that the eight signatures were significantly associated with survival of patients with BRCA (Fig. 6b). Additionally, all eight signatures were identified as independent prognostic predictors for BRCA in multivariate COX analyses (Fig. 6c-j). Moreover, the correlations between ER, PR, HER2 status and prognostic risk defined by AS signatures were analyzed (Table S4). Positive ER status was associated with low risk defined by AA and ALL signatures, and high risk defined by AT signature. Positive PR status was associated with low risk defined by AA, AP, ME and ALL signatures, and high risk defined by AT signature. Positive HER2 status was associated with high risk defined by AP, AT and ALL signatures.
Prognostic SF-AS network
The Pearson correlation analysis showed that there were 30 SFs were negatively correlated with 20 AS events, whereas 34 SFs positively correlated with 35 AS events (Table S8). A regulatory network was generated based on the correlation between AS and SFs, which consist of 29 protective AS events (associated with good prognosis), 9 risk AS events (associated with poor prognosis) and 38 SFs, (Fig. 7a). Among the SFs, CCDC12, CLASRP and LUC7L were significantly correlated with more than 10 AS events, therefore they were considered as a core SF. AS events INPP5F-13,276-RI and NRBP2–85507-RI were both regulated by more 13 SFs, they might play important roles in AS of BRCA. Moreover, the prognostic significance of the correlating SFs we analyzed. There were 8 correlating SFs (HSPA8, U2AF1L4, SNRNP70, SRSF5, CLASRP, CCDC12, SART1, and WDR83) significantly associated with survival of BRCA patients (Table S9).
Functional enrichment analysis
The results of GO analysis suggested that AS genes were implicated in carcinogenesis associated biological processes such as “negative regulation of mitotic cell cycle”, “cell cycle G2/M phase transition”, “cell−matrix adhesion” (Fig. 7b, Table S5). In the KEGG analysis, AS genes were enriched in pathways associated with cancer, such as “MAPK signaling pathway”, “Small cell lung cancer”, and “Cell cycle” (Fig. 7c, Table S5).
Alternative splicing is a vital process involved in the RNA transcription and modification of mRNA isoforms [5, 15]. Increasing evidence has demonstrated that abnormal AS is associated with the carcinogenesis of multiple cancers [16,17,18,19]. Hence, exploration of AS mechanisms deepens our understanding of posttranscriptional regulatory patterns.
With the rapid development of the next-generation sequencing technology, progress has been made in the field of bioinformatics. TCGA and SpliceSeq database provides researchers with a great amount of high quality RNA sequencing data, which enabled the studies of AS patterns in various cancer types [20, 21]. To our knowledge, several studies reported AS profiles and established prognostic prediction model for several cancers, including kidney renal clear cell carcinoma , hepatocellular carcinoma , esophageal carcinoma , prostate adenocarcinoma , colorectal cancer , and soft tissue sarcoma .
Our study demonstrated that 45,421 AS events in 21,232 mRNAs were found in BRCA, and 1121 AS events in 931 genes are significantly correlated with the survival of BRCA patients. Zhang et al. identified 3071 AS events in breast cancer patients as significant for prognoses . In both Zhang’s and ours studies, overall survival associated AS events were analyzed by univariate Cox proportional hazard regression analysis. In Zhang’s study, BRCA cohort were first divided into two groups (categorical variable) by a median cut of PSI value before survival analysis. In our study, we conduct Cox survival analysis with the original PSI value (continuous variable), instead of deliberately separate the cohort into two groups. Therefore, the number of identified AS events is less in our study. By integrating all seven types of AS, a comprehensive prognostic signature was generated, which included PARPBP , NCOR1 , COPZ1 , ANK3, ITGB5, PHTF1 , HSPBP1, TCF12 , RPS6KA1, CNST, TMEM25  and BTN3A2 , which are play essential roles in carcinogenesis. Claudia et al. reported that PARPBP inhibits activation of the NF-κB pathway, which can initiate p21-mediated differentiation and proliferation arrest . Wang et al. suggested that NCoR1 may act as tumor suppressors in GIST cells through the Smad signaling pathway . Maria et al. indicated that COPZ1 represents an example of non-oncogene addiction in thyroid tumor cells, COPZ1 depletion affects thyroid tumor cell viability in vivo and in vitro . Huang et al. found that PHTF1 may be a tumor-suppressor like gene and a therapeutic target for triggering the PHTF1-FEM1b-Apaf-1 apoptosis pathway . Yang et al. revealed that TCF12 promotes the tumorigenesis and metastasis of hepatocellular carcinoma via upregulation of CXCR4 expression . Moreover, a recent study suggested that BTN3A2 serves as a prognostic marker and favors immune infiltration in triple-negative breast cancer .
The comprehensive signature can serve as a useful tool to predict the survival outcomes of patients with BRCA with an AUC value of 0.801. Besides, we found there were correlation between breast cancer subtype (HER2, ER/PR) and prognostic risk defined by AS signatures. Accordingly, we speculated that HER2, ER and PR status might affect alternative splicing of certain genes that associated with cancer progression and survival outcome. More in-depth researches are warranted to provide novel insights into the molecular mechanism of BRCA. Additionally, an SFs-AS network was created and we found that CCDC12, CLASRP and LUC7L might serve as core SFs on account of their significant correlation with multiple AS events.
The AS prognostic signatures accurately predict survival outcomes of BRCA patients, suggesting that AS signatures might act as ideal prognostic indicators. The SFs-AS regulatory network demonstrated the molecular mechanisms of AS in BRCA. Our study may provide potential therapeutic targets for future BRCA management.
Availability of data and materials
All data are included in this article and supplementary files. The original data are available upon reasonable request to the corresponding author. The data of alternative splicing were downloaded from TCGA SpliceSeq database (https://bioinformatics.mdanderson.org/TCGASpliceSeq/). The expression and clinical data were downloaded from TCGA Genomic Data Commons Data Portal (https://portal.gdc.cancer.gov/). The public access to these databases is open.
The Cancer Genome Atlas
Percent Spliced In
Alternate acceptor site
Alternate donor site
Mutually exclusive exons
Least absolute shrinkage and selection operator
Receiver operating characteristic
Area under the curve
Kyoto Encyclopedia of Genes and Genomes
Chen W, Zheng R, Baade PD, Zhang S, Zeng H, Bray F, et al. Cancer statistics in China, 2015. CA Cancer J Clin. 2016;66(2):115–32. https://doi.org/10.3322/caac.21338.
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424. https://doi.org/10.3322/caac.21492.
Brentnall AR, Cuzick J, Buist DSM, Bowles EJA. Long-term accuracy of breast Cancer risk assessment combining classic risk factors and breast density. JAMA Oncol. 2018;4(9):e180174. https://doi.org/10.1001/jamaoncol.2018.0174.
DeSantis CE, Ma J, Goding Sauer A, Newman LA, Jemal A. Breast cancer statistics, 2017, racial disparity in mortality by state. CA Cancer J Clin. 2017;67(6):439–48. https://doi.org/10.3322/caac.21412.
Climente-González H, Porta-Pardo E, Godzik A, Eyras E. The functional impact of alternative splicing in Cancer. Cell Rep. 2017;20(9):2215–26. https://doi.org/10.1016/j.celrep.2017.08.012.
Nilsen TW, Graveley BR. Expansion of the eukaryotic proteome by alternative splicing. Nature. 2010;463(7280):457–63. https://doi.org/10.1038/nature08909.
El Marabti E, Younis I. The Cancer Spliceome: reprograming of alternative splicing in Cancer. Front Mol Biosci. 2018;5:80. https://doi.org/10.3389/fmolb.2018.00080.
Lee SC, Abdel-Wahab O. Therapeutic targeting of splicing in cancer. Nat Med. 2016;22(9):976–86. https://doi.org/10.1038/nm.4165.
Kozlovski I, Siegfried Z, Amar-Schwartz A, Karni R. The role of RNA alternative splicing in regulating cancer metabolism. Hum Genet. 2017;136(9):1113–27. https://doi.org/10.1007/s00439-017-1803-x.
Sciarrillo R, Wojtuszkiewicz A, Assaraf YG, Jansen G, Kaspers GJL, Giovannetti E, et al. The role of alternative splicing in cancer: from oncogenesis to drug resistance. Drug Resist Updat. 2020;53:100728. https://doi.org/10.1016/j.drup.2020.100728.
Bonnal SC, López-Oreja I, Valcárcel J. Roles and mechanisms of alternative splicing in cancer - implications for care. Nat Rev Clin Oncol. 2020;17(8):457–74. https://doi.org/10.1038/s41571-020-0350-x.
Singh B, Eyras E. The role of alternative splicing in cancer. Transcription. 2017;8(2):91–8. https://doi.org/10.1080/21541264.2016.1268245.
Ryan M, Wong WC, Brown R, Akbani R, Su X, Broom B, et al. TCGASpliceSeq a compendium of alternative mRNA splicing in cancer. Nucleic Acids Res. 2016;44(D1):D1018–22. https://doi.org/10.1093/nar/gkv1288.
Seiler M, Peng S, Agrawal AA, Palacino J, Teng T, Zhu P, et al. Somatic Mutational Landscape of Splicing Factor Genes and Their Functional Consequences across 33 Cancer Types. Cell Rep. 2018;23(1):282–96 e284.
Yang X, Coulombe-Huntington J, Kang S, Sheynkman GM, Hao T, Richardson A, et al. Widespread expansion of protein interaction capabilities by alternative splicing. Cell. 2016;164(4):805–17. https://doi.org/10.1016/j.cell.2016.01.029.
Wong ACH, Rasko JEJ, Wong JJL. We skip to work: alternative splicing in normal and malignant myelopoiesis. Leukemia. 2018;32(5):1081–93. https://doi.org/10.1038/s41375-018-0021-4.
Jin Y, Dong H, Shi Y, Bian L. Mutually exclusive alternative splicing of pre-mRNAs. Wiley Interdiscip Rev RNA. 2018;9(3):e1468. https://doi.org/10.1002/wrna.1468.
Ule J, Blencowe BJ. Alternative splicing regulatory networks: functions, mechanisms, and evolution. Mol Cell. 2019;76(2):329–45. https://doi.org/10.1016/j.molcel.2019.09.017.
Zhang J, Manley JL. Misregulation of pre-mRNA alternative splicing in cancer. Cancer Discov. 2013;3(11):1228–37. https://doi.org/10.1158/2159-8290.CD-13-0253.
Lin P, He R-Q, Ma F-C, Liang L, He Y, Yang H, et al. Systematic analysis of survival-associated alternative splicing signatures in gastrointestinal pan-adenocarcinomas. EBioMedicine. 2018;34:46–60. https://doi.org/10.1016/j.ebiom.2018.07.040.
Zhang Y, Yan L, Zeng J, Zhou H, Liu H, Yu G, et al. Pan-cancer analysis of clinical relevance of alternative splicing events in 31 human cancers. Oncogene. 2019;38(40):6678–95. https://doi.org/10.1038/s41388-019-0910-7.
Song J, Liu YD, Su J, Yuan D, Sun F, Zhu J. Systematic analysis of alternative splicing signature unveils prognostic predictor for kidney renal clear cell carcinoma. J Cell Physiol. 2019;234(12):22753–64. https://doi.org/10.1002/jcp.28840.
Chen Q-F, Li W, Wu P, Shen L, Huang Z-L. Alternative splicing events are prognostic in hepatocellular carcinoma. Aging (Albany NY). 2019;11(13):4720–35. https://doi.org/10.18632/aging.102085.
Mao S, Li Y, Lu Z, Che Y, Sun S, Huang J, et al. Survival-associated alternative splicing signatures in esophageal carcinoma. Carcinogenesis. 2019;40(1):121–30. https://doi.org/10.1093/carcin/bgy123.
Huang Z-G, He R-Q, Mo Z-N. Prognostic value and potential function of splicing events in prostate adenocarcinoma. Int J Oncol. 2018;53(6):2473–87. https://doi.org/10.3892/ijo.2018.4563.
Xiong Y, Deng Y, Wang K, Zhou H, Zheng X, Si L, et al. Profiles of alternative splicing in colorectal cancer and their clinical significance: a study based on large-scale sequencing data. EBioMedicine. 2018;36:183–95. https://doi.org/10.1016/j.ebiom.2018.09.021.
Yang X, Huang W-T, He R-Q, Ma J, Lin P, Xie Z-C, et al. Determining the prognostic significance of alternative splicing events in soft tissue sarcoma using data from the Cancer genome atlas. J Transl Med. 2019;17(1):283. https://doi.org/10.1186/s12967-019-2029-6.
Zhang D, Duan Y, Cun J, Yang Q. Identification of prognostic alternative splicing signature in breast carcinoma. Front Genet. 2019;10:278. https://doi.org/10.3389/fgene.2019.00278.
Nicolae CM, O'Connor MJ, Schleicher EM, Song C, Gowda R, Robertson G, et al. PARI (PARPBP) suppresses replication stress-induced myeloid differentiation in leukemia cells. Oncogene. 2019;38(27):5530–40. https://doi.org/10.1038/s41388-019-0810-x.
Wang W, Song XW, Bu XM, Zhang N, Zhao CH. PDCD2 and NCoR1 as putative tumor suppressors in gastric gastrointestinal stromal tumors. Cell Oncol. 2016;39(2):129–37. https://doi.org/10.1007/s13402-015-0258-0.
Anania MC, Cetti E, Lecis D, Todoerti K, Gulino A, Mauro G, et al. Targeting COPZ1 non-oncogene addiction counteracts the viability of thyroid tumor cells. Cancer Lett. 2017;410:201–11. https://doi.org/10.1016/j.canlet.2017.09.024.
Huang X, Geng S, Weng J, Lu Z, Zeng L, Li M, et al. Analysis of the expression of PHTF1 and related genes in acute lymphoblastic leukemia. Cancer Cell Int. 2015;15(1):93. https://doi.org/10.1186/s12935-015-0242-9.
Wang F, Wei XL, Wang FH, Xu N, Shen L, Dai GH, et al. Safety, efficacy and tumor mutational burden as a biomarker of overall survival benefit in chemo-refractory gastric cancer treated with toripalimab, a PD-1 antibody in phase Ib/II clinical trial NCT02915432. Ann Oncol. 2019;30(9):1479–86. https://doi.org/10.1093/annonc/mdz197.
Hrašovec S, Hauptman N, Glavač D, Jelenc F, Ravnik-Glavač M. TMEM25 is a candidate biomarker methylated and down-regulated in colorectal cancer. Dis Markers. 2013;34(2):93–104. https://doi.org/10.1155/2013/427890.
Cai P, Lu Z, Wu J, Qin X, Wang Z, Zhang Z, et al. BTN3A2 serves as a prognostic marker and favors immune infiltration in triple-negative breast cancer. J Cell Biochem. 2020;121(3):2643–54. https://doi.org/10.1002/jcb.29485.
Ethics approval and consent to participate
This study was conducted following the TCGA publication guidelines, and approval from a local Ethics Committee were unnecessary.
Consent for publication
The authors declare that they have no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Differential AS eventsR3.
AS events significantly associated with OSR3.
List of 404 splice factorsR3.
ER, PR, and HER2 correlationR3.
Detailed results of GO and KEGG analysesR3.
Differential expressed genesR3.
AS events Ratio for 7 typesR3.
Correlation between AS events and SFsR3.
About this article
Cite this article
Han, P., Zhu, J., Feng, G. et al. Characterization of alternative splicing events and prognostic signatures in breast cancer. BMC Cancer 21, 587 (2021). https://doi.org/10.1186/s12885-021-08305-6
- Alternative splicing
- Breast cancer
- The Cancer genome atlas