Skip to main content

Advertisement

Genome-wide analysis reveals the association between alternative splicing and DNA methylation across human solid tumors

Abstract

Background

Dysregulation of alternative splicing (AS) is a critical signature of cancer. However, the regulatory mechanisms of cancer-specific AS events, especially the impact of DNA methylation, are poorly understood.

Methods

By using The Cancer Genome Atlas (TCGA) SpliceSeq and TCGA data for ten solid tumor types, association analysis was performed to characterize the potential link between cancer-specific AS and DNA methylation. Functional and pathway enrichment analyses were performed, and the protein-protein interaction (PPI) network was constructed with the String website. The prognostic analysis was carried out with multivariate Cox regressions models.

Results

15,818 AS events in 3955 annotated genes were identified across ten solid tumor types. The different DNA methylation patterns between tumor and normal tissues at the corresponding alternative spliced exon boundaries were shown, and 51.3% of CpG sites (CpGs) revealed hypomethylated in tumors. Notably, 607 CpGs were found to be highly correlated with 369 cancer-specific AS events after permutation tests. Among them, the hypomethylated CpGs account for 52.7%, and the number of down-regulated exons was 173. Furthermore, we found 38 AS events in 35 genes could serve as new molecular biomarkers to predict patient survival.

Conclusions

Our study described the relationship between DNA methylation and AS events across ten human solid tumor types and provided new insights into intragenic DNA methylation and exon usage during the AS process.

Peer Review reports

Background

Alternative splicing (AS) is one of the conserved biological processes diversifying the transcriptome and proteome [1]. Pre-mRNA splicing is a key step of gene expression in which introns within nascent RNA are removed, and exons are ligated to form mature mRNA [2]. It is estimated that more than 90% of human genes with multiple exons undergo AS during pre-mRNA maturation, highlighting the importance of AS in determining gene function [3, 4]. Recently, AS has been implicated as an important signature to understand tumorigenesis, cancer progression, and resistance to therapy [5]. Several lines of evidence have shown that the disruption of AS is frequently associated with the inactivation of tumor suppressors and activation of oncogenes [6, 7]. For example, Bcl-x, an apoptosis regulator, is found to be switched from its pro-apoptotic into anti-apoptotic splicing isoforms in a number of cancer types [8, 9]. Another well-characterized example is CD44, whose different splicing isoforms have been associated with tumor evasion and metastasis [10, 11]. Moreover, Kahles and colleagues have recently provided a comprehensive landscape of AS events across different tumor types [12]. Other studies have also reported the prognostic value of AS events in multiple cancer types, such as ovarian cancer [13], breast cancer [14], glioblastma [15] and gastrointestinal adenocarcinomas [16], suggesting a predominant role of splicing dysregulation in cancers.

The splicing reaction is generally regulated by cis-elements within the pre-mRNA and trans-acting splicing factors that bind to these cis-elements [17]. Alterations of splicing factors, mutations in spliceosomal proteins, and regulation exerted by pre-mRNA cis-elements can all contribute to the splicing alterations in cancers [18,19,20]. Additionally, epigenetic modification has been recently proposed as another regulator of alternative pre-mRNA splicing patterns [11, 21, 22]. The genome-wide mapping analysis has unveiled that DNA methylation was a strong marker for exon boundaries and suggested the possible role of DNA methylation in exon definition and splicing regulation [23, 24]. Anastasiadou et al. have reported that the increased CpG methylation was frequent in alternatively spliced sites [25]. In support, by studying wild-type and methylation deficient embryonic stem cells, Yearim et al. have found that DNA methylation could affect the splicing of more than 20% of alternative exons [26].

The global change of splicing in tumor tissues in comparison to their normal tissue counterparts is increasingly appreciated. However, the regulatory mechanisms underlying cancer-specific AS events, especially the influence of DNA methylation on AS, remain poorly understood. To figure out the association between AS events and DNA methylation of alternatively spliced exons, we first compared the transcriptome-wide splicing between tumor and matched normal tissues across ten solid tumor types. Further, we compared the DNA methylation of CpG sites (CpGs) at the boundaries of alternatively spliced exons. Finally, we utilized correlation analysis and permutation test to assess whether there is an association between DNA methylation and cancer-specific AS events. By these means, we expect to explore the role of methylation in exon usage during transcriptional process and carcinogenesis.

Methods

Overall workflow

The overall workflow of the present study is shown in Fig. 1. Using TCGA SpliceSeq dataset for ten solid tumor types, AS events were analyzed in 580 paired tumor-normal tissues. Then, the methylation levels of CpGs at the boundaries of alternatively spliced exons were compared between 641 normal and tumor samples from TCGA. Subsequently, we correlated cancer-specific AS events with methylation features and observed a significant association between them. To further understand the function and mechanism of the genes containing methylation-associated cancer-specific AS events, functional enrichment, and network analyses were applied. Finally, we found that a set of cancer-specific AS events could serve as reliable prognostic biomarkers for cancers.

Fig. 1
figure1

Flowchart of the present study

Data acquisition

Data on transcriptome alterations were downloaded from TCGA SpliceSeq portal (http://bioinformatics.mdanderson.org/TCGASpliceSeq/). TCGA SpliceSeq is a web-based bioinformatics resource containing the mRNA splicing patterns of 33 types of TCGA cancers and adjacent normal samples. SpliceSeq can be used to calculate the “percent spliced in” (PSI) value for each splicing event in each cancer sample [3], which provides a clear view of the splice junction and the proportion of exons included in different samples. PSI values were utilized to quantify seven types of AS events: exon skip (ES), alternate donor site (AD), alternate acceptor site (AA), retained intron (RI), mutually exclusive exons (ME), alternate terminator (AT) and alternate promoter (AP) [27]. In the present study, AS events (379,749 in total) that have PSI values in more than 75% of samples were included.

Data of somatic mutations, copy number variations, mRNA expression, and DNA methylations were downloaded from publicly released TCGA level 3 data (processed data from USCS Cancer Genome Browser). In this study, ten types of cancers were included: breast invasive carcinoma (BRCA), colon adenocarcinoma (COAD), head and neck squamous cell carcinoma (HNSC), kidney renal clear cell carcinoma (KIRC), kidney renal papillary carcinoma (KIRP), liver hepatocellular carcinoma (LIHC), lung adenocarcinoma (LUAD), lung squamous cell carcinoma (LUSC), prostate adenocarcinoma (PRAD) and thyroid carcinoma (THCA). Illumina Human Methylation 450 BeadChip (450 K array) was used for DNA methylation data. Patient clinical information was also downloaded from TCGA dataset. The sample size of each cancer is summarized in Additional file 1: Table S1.

Gene coordinates and RefSeq annotations were obtained from UCSC (Jul 2013 release, hg 19) [28]. All the features of exons and introns referred to RefSeq genes. The ‘exon boundary’ was determined as 200 nucleotides in size for intronic regions and 39 nucleotides for exonic regions [29].

Identification of cancer-specific AS events and methylated CpGs

Wilcoxon paired test was applied to compare the PSI values between tumor tissues and matched adjacent normal tissues. In each cancer type, we identified cancer-specific AS events that satisfied the following criteria: (1) the distributions of PSI values were significantly different between tumor and normal tissues (FDR < 0.05); (2) the change of PSI values (|PSI|) between tumor and normal tissues was greater than 0.1. Differentially methylated CpGs were detected by applying a paired t-test, setting the absolute average methylation change (|meth|) ≥ 0.1 and FDR < 0.05 as selection criteria. All statistical analyses were performed using R software (version 3.6.0).

Correlation analysis and permutation test

Patients with both splicing information and methylation data were used for correlation analysis. Covariance analysis model, after adjusting for copy number variants and somatic mutation in tumor, as well as Spearman correlation analysis in normal tissues, were performed to investigate the correlations between PSI and DNA methylation respectively. In order to infer their statistical significance, we further carried out the permutation test. The number of permutations was 1000, and we considered the permutation P < 0.05 was statistically significant. In the present study, the aovperm() function in R package “permuco” and spearman_test() function in package “coin” were used for permutation tests.

Functional enrichment analysis

Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were conducted on genes with methylation-associated cancer-specific AS events using “ClusterProfiler” package in R software [30]. GO and KEGG enrichment analyses were based on the threshold of FDR < 0.05.

Protein-protein interaction analysis

Genes with significantly methylation-associated AS events were analyzed through the protein-protein interaction (PPI) network. The PPI analysis was constructed with the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING, https://string-db.org/cgi/input.pl) [31]. The confidence score ≥ 0.4 and the maximum number of interactors = 0 were set as the selection criteria.

Survival analysis

To assess the prognostic value of cancer-specific AS events and methylated CpGs, we performed survival analysis by using the “survival” package in R. The Cox proportional hazard model was utilized to test the interactions between different variables and overall survival in multivariate analysis by adjusting for age, gender, TNM stage and adjuvant therapy (including chemotherapy and radiotherapy). P < 0.05 was considered as the threshold for significance.

Results

Identification of cancer-specific AS events in common solid tumors

To systematically characterize abnormal AS events in common solid tumors, we compared PSI values between tumors and matched adjacent normal tissues. We found that 15,818 AS events in 3955 annotated genes were significantly altered in cancer vs. normal tissues. Splicing alterations were most abundant in KIRC and LUAD, while least in COAD (Additional file 1: Table S2). For KIRC, a total of 2678 differential AS events in 1550 genes were identified, including 81 AAs in 74 genes, 91 ADs in 83 genes, 1190 APs in 643 genes, 535 ATs in 286 genes, 611 ESs in 474 genes, 12 MEs 12 genes, and 142 RIs in 137 genes. For LUAD, a total of 2212 differential AS events in 1369 genes were identified (Fig. 2a).

Fig. 2
figure2

Cancer-specific AS events across human solid tumors. a Bar plot showing the numbers of cancer-specific AS events in different tumors split according to types of event. b Bar plot showing the numbers of cancer-specific AS events in seven different types of AS. c The UpSet plot of interactions between the seven types of cancer-specific AS events. One gene may have several types of AS events

For the seven AS events types, we found that AP and ES took up 44.0 and 21.2% of the AS events, respectively, followed by AT (19.9%). In addition, we observed that the numbers of AS events with increased vs. decreased PSI values were similar, with the exception of RI and AD events. (Fig. 2b).

Notably, one gene could possess more than one cancer-specific AS events. The numbers of genes that possess differential AS events in varied cancer types are depicted in Fig. 2c.

Identification of differentially methylated CpGs at exon boundaries

To examine the differences of DNA methylation, we retrieved methylation data of CpGs at exon boundaries from TCGA data portal. After comparing between tumor and matched normal tissues, we identified a large amount of CpGs with differential methylation levels (Fig. 3a). Totally, 1180 CpGs were detected, among which KIRC displayed the highest number of CpGs, including 5 CpGs in AA, 3 CpGs in AD, 186 CpGs in AP, 23 CpGs in AT, 15 CpGs in ES and 18 CpGs in RI (Additional file 1: Table S3). We also observed that the proportions of hypomethylated CpGs in different solid tumors were similar, with the exception of COAD (32.4%) and PRAD (10.9%) (Fig. 3b). When comparing the differentially methylated CpGs in varied types of AS events, 833 CpGs were identified in AP type, of which 48.6% were hypomethylated. The proportions of hypomethylated CpGs in other types of AS were similar, ranging from 43.9% in ES to 66.7% in AD (Fig. 3c). Taken together, hypomethylation was a more frequently observed form of differential methylation in tumor tissues, compared to hypermethylation.

Fig. 3
figure3

Differentially methylated CpGs at alternatively spliced exon boundaries in multiple tumors. a Volcano plot showing a comparison of the methylation level of CpGs for tumor tissues versus adjacent normal tissues. This plot depicts the biological significance (|meth|) ≥ 10%) on the X axis and the statistical significance (−log10 FDR) on the Y axis. b The proportion of hypomethylated and hypermethylated CpGs at cancer-specific alternative spliced exon boundaries in different tumors. c The proportion of hypomethylated and hypermethylated CpGs at cancer-specific alternative spliced exon boundaries in different splicing types

Correlation between cancer-specific AS events and DNA methylation

To explore possible regulatory mechanisms of cancer-specific AS, we further applied correlation analysis between the PSI values and methylation levels of CpGs at alternatively spliced exon boundaries. A significant association between PSI and DNA methylation was observed from the covariance analysis model after adjusting for covariates, including copy number variants and somatic mutation. After performing the permutation test, 891 CpGs were found to be highly correlated with 483 cancer-specific AS events (permuted P < 0.05) (Additional file 1: Table S4). To better determine whether the changes in AS resulted from changed methylation levels or vice versa, we also performed association analysis in the normal tissues. Results showed that the methylation levels of 607 CpGs were not significantly correlated to PSI in the normal tissues after permutation tests (permuted P < 0.05) (Additional file 1: Table S4). The number of PSI-associated CpGs is presented in Fig. 4a. Among them, 320 CpGs were hypomethylated, and 287 were hypermethylated in tumor tissues. Expect for COAD, LIHC, LUAD, LUSC, and PRAD, the hypomethylated CpGs account for more than half of AS-associated CpGs in cancers (ranging from 54.5 to 95.7%)(Additional file 2: Figure S1-A). Among the exons, 173 exons were down-regulated (ranging from 35.7% in THCA to57.6% in LUSC) in the tumor tissues (Additional file 2: Figure S1-B).

Fig. 4
figure4

A significant correlation between cancer-specific AS events and DNA methylation at exon boundaries. a Column plot of the number of alternative splicing-related CpGs in different cancer types split according to type pf event. b Circos plot showing the distribution of methylation-associated splicing events across the human genome in ten solid tumors. The track outside the chromosome track showed the distribution of gene density across the genome. Inside tracks denote the Spearman correlation coefficient (ranging from − 1 to 1) between PSI and methylation level across the genome in each human tumor (from outside to inside): BRAC, COAD, HNSC, KIRC, KIRP, LIHC, LUAD, LUSC, PRAD, THCA. c Common change of ZNF577 in five cancer types. Upper: Box plot of PSI values of exon 1.2 and 1.3 in ZNF577 in tumor and normal samples in five cancer types. Middle: Box plot of the methylation level of CpG site (cg11269599) in ZNF577 in tumor and normal samples in five cancer types. Lower: Scatter plot showing the significant relationship between PSI values and methylation levels in five cancer types

Though the regions of splicing-related CpGs showed differently across varied chromosomes and regions (Fig. 4b), some common features in different tumor types were sill discovered. The common AS events having similar relationships with CpGs, which changed in at least three cancer types are shown in Table 1. This result suggested that there were indeed a significant number of AS-associated CpGs shared by multiple cancer types. Taking ZNF577 and cg11269599 as examples, we found that the exon 1.2 and exon 1.3 of ZNF577 showed increased PSI values in five cancer types when comparing tumors to the cognate normal tissues. The CpG site, cg11269599, also showed significantly hypermethylated in HNSC, KIRC, KIRP, LUSC, and PRAD. The correlations between both aforementioned were also observed consistently in the five cancer types (Fig. 4c). The preference of a specific splice variant by several cancer types could help to confirm the link between AS and DNA methylation [32].

Table 1 Splicing events and DNA methylation altered in the same direction in at least three cancer types

Gene functional enrichment analysis and network construction

We further conducted GO analysis on genes containing methylation-related cancer-specific AS events for different cancer types. Results showed that these genes were enriched in different GO terms across different cancers. In COAD, high enrichment of genes was associated with cellular components, including sarcolemma and basal part. In THCA, genes were mainly associated with molecular functions, including hormone binding and acyl-CoA ligase activity (Additional file 1: Table S5). KEGG pathway analysis was also performed. Only the pathway “Glycerophospholipid metabolism” was found significantly enriched by genes in LIHC (Additional file 1: Table S6). These results suggested that the activity of these genes might be in an independent manner in different cancer types.

Moreover, we constructed a PPI network using the STRING database. In the PPI network, NCK1 and DOK1 were the top hub genes in the HNSC network. PIK3R1 PRKACA and LCK were the top three hub genes in the LUSC network. However, in other cancers, no PPI networks were constructed with statistical significance (Additional file 3: Figure S2).

Survival analysis of AS events and CpGs

We then explored whether the associated AS events and CpGs could be effectively used as a prognostic signature for each cancer. Cox regression analysis was performed adjusting for age, gender, TNM stage and adjuvant therapy, and results showed that 10.3% (38/369) AS events could significantly distinguish patients with longer versus shorter survival prognoses (Table 2). Because splicing of introns could be regulated under a co-transcriptional mechanism [33], an alteration in gene expression may affect AS of corresponding genes. Thus, we examined mRNA expression of genes containing cancer-specific AS events and found that the majority of them (85.71%, 30/35) were not significantly associated with survival, suggesting that these AS events are partially influenced by their gene expression (Table 2). Additionally, we evaluated DNA methylation of CpGs at the survival-related AS exons boundaries. 14 (20.6%) CpGs in 8 (21.1%) AS exons were found significantly associated with the overall survival of cancer patients. (Table 2).

Table 2 Multivariate Cox proportional hazard model analyses of associations between alternative exons, methylation of corresponding CpG and cancer patients’ survival

Discussion

Collectively, we performed a comprehensive analysis of the relationship between cancer-specific AS events and DNA methylation. We found that approximately half of the AS events were AP and ES types in human solid tumors. The boundaries of alternatively spliced exons were more likely to be hypomethylated in tumor tissues compared with adjacent normal tissues. Association analyses and permutation tests revealed that cancer-specific AS events were significantly correlated with DNA methylation. Our research provided a novel perspective of the regulatory mechanism of cancer-specific AS.

It is well accepted that AS plays an important role in multiple cellular processes and development programs, as well as contributing to tumorigenesis [34, 35]. Nonetheless, to our knowledge, few studies have investigated the role of DNA methylation as a regulator of the cancer-specific AS. Recently, Laurent et al. have demonstrated that methylation differences were more prominent at the exon-intron boundaries [36]. Other studies also found that DNA methylation links with the inclusion rate of alternative exons [26, 37]. These findings supported the hypothesis that AS could be functional by cis-regulation of DNA methylation. With the advantage of high-throughput data, TCGA data portal provides opportunities for the integration analyses of multi-omics data. TCGA Splice Seq, web-based bioinformatics, providing a clear view of the mRNA splicing patterns of 33 tumor types, across a dataset of more than 10,000 TCGA samples. Ryan et al. identified and calculated each potential splicing event across 33 types of cancer to establish TCGA SpliceSeq database, while did not evaluate the potential mechanism and clinical usage of AS events [27]. In the present study, we integrated AS events from SpliceSeq and TCGA data together to comprehensively explore potential regulatory mechanisms of DNA methylation for cancer-specific AS. We performed an association analysis between cancer-specific AS events and CpGs at their exon boundaries. The parameter selected in this study was based on previous reports [38,39,40], especially Castle’s study, in which they examined the exon neighborhoods in the size of 200 nucleotides for intronic regions and 39 nucleotides for exonic regions to identify splicing cis-regulatory elements in sequences [29].

In the present study, to ensure the AS difference resulted from DNA methylation, we excluded CpGs showing evidence of associations with AS in normal tissues. By using such analysis, we identified a number of PSI-related CpGs, whose locations appeared differently across varied chromosomes and regions in different cancers. Moreover, the higher proportion of decreased methylation levels of AS-associated CpGs were observed in the majority of solid cancer types, except for COAD, LIHC, LUAD, LUSC, and PRAD. We speculated that this was due to the heterogeneity among tumor samples. It is well known that the splicing process is predominantly tissue- and cancer-specific [12, 41]. Meanwhile, DNA methylation is also demonstrated to be varied across tumor types [42]. It is worth to note that several cases showing the same pattern of changes were identified across different cancer types. These signatures were potential candidates for actively selected variants, which might drive tumorigenesis. For example, we found the consistent changes of AS event and cg11269599 in gene ZNF577 (zinc finger protein 577) across several tumor types. Zinc finger proteins are commonly involved in transcriptional regulation of genes, but the mechanism of how DNA methylation affects AS of ZNF577 has not been clarified and needs to be determined in further studies [43].

Recently, several studies have given insight into the regulation of DNA methylation in AS. It has been suggested that DNA methylation at exonic regions could affect the binding of methyl-sensitive DNA-binding proteins, such as MeCp2, CTCF, and HP1 [37, 44,45,46]. MeCp2 was found to be enriched in spliced exons, which could facilitate the recruitment of histone deacetylase, thus causing Pol II pause and exon inclusion [46]. In 2011, Shukla et al. found the DNA methylation-mediated splicing of CD45 pre-mRNA, which inhibited CTCF binding [37]. HP1 could bring several splicing factors to methylated exons [26]. In addition, DNA methylation had also been reported to determine nucleosome positioning and affected the transcriptional elongation rate [47].

We also investigated the prognostic values of the cancer-specific AS events by conducting Cox regression. We found 38 methylation-related cancer-specific AS events could significantly distinguish patients with longer versus shorter survival prognoses. It should be noted that these identified signatures were overrepresented in renal carcinoma. Renal carcinoma has various histological and molecular subtypes and is characterized by poor prognosis and high recurrence. Until now, several molecular biomarkers have been investigated for renal carcinoma, while none of them have been used in clinical practice [48, 49]. Therefore, the identification of novel and effective prognostic biomarkers is important for patients suffering from renal carcinoma. The signatures identified in the present study in renal cancers were not well reported previously and could provide the basis for further study into the pathomechanisms and serve as potential therapeutic targets. Additionally, 14 CpGs at the exon boundaries of 8 survival-related AS also showed significant association with overall survival, even after adjusting for age, gender, TNM stage, and adjuvant therapy. These results led us to speculate that the identified CpGs might be the specific cancer drivers, and serve as prognostic biomarkers for cancers.

By the integrated analyses of multi-dimensional data, our study revealed a deeper understanding of cancer-specific AS events and their relationships with DNA methylation across different solid tumor types. Another strength of our study was the pan-cancer analysis, which could reveal similarities and specificities across different cancer types, as well as enabled the elimination of false-positive and false-negative calls made in several single-tumor-type projects [50]. Finally, we provided novel mechanisms and therapeutic perspectives for either one or several cancer types. The present study has several limitations. First, our study investigated the methylation-splicing features across multiple cancer types instead of the difference within cancer. In addition, we only studied the differentially methylated CpGs between tumor and normal tissues, which could miss the CpGs with hyper- or hypo-methylation in tumors that were associated with either down- or up-regulation of PSI. Finally, our research was based on the public data which lacked several important clinical features, such as a history of drug use, and has not been studied in survival analysis. Therefore, further research is needed to clarify and demonstrate the nature of the regulation of methylation in AS.

Conclusion

Overall, our analyses identified a significant association between DNA methylation and cancer-specific AS events. The present study contributes to the understanding of the role of methylation in exon usage during transcriptional process and carcinogenesis. Finally, this relatively small set of AS events and CpGs will facilitate the discovery of critical regulators, which are responsible for splicing dysregulation in cancers and thus can be used as new therapeutic biomarkers.

Availability of data and materials

TCGA SpliceSeq portal (http://projects.insilico.us.com/TCGASpliceSeq/PSIdownload.jsp) applies alternative splicing data. Information on somatic mutations, copy number variations, and DNA methylations of the human solid tumors in this study are from the UCSC Cancer Genomics Browser (https://xenabrowser.net/datapages/).

Abbreviations

AA:

Alternate acceptor site

AD:

Alternate donor site

AP:

Alternate promoter

AS:

Alternative splicing

AT:

Alternate terminator

BRCA:

Breast invasive carcinoma

COAD:

Colon adenocarcinoma

CpGs:

CpG site

ES:

Exon skip

FDR:

False discovery rate

GO:

Gene Ontology

HNSC:

Head and neck squamous cell carcinoma

KEGG:

Kyoto Encyclopedia of Genes and Genomes

KIRC:

Kidney renal clear cell carcinoma

KIRP:

Kidney renal papillary carcinoma

LIHC:

Liver hepatocellular carcinoma

LUAD:

Lung adenocarcinoma

LUSC:

Lung squamous cell carcinoma

ME:

Mutually exclusive exons

PRAD:

Prostate adenocarcinoma

PSI:

Percent spliced in

RI:

Retained intron

TCGA:

The Cancer Genome Atlas

THCA:

Thyroid carcinoma

References

  1. 1.

    Nilsen TW, Graveley BR. Expansion of the eukaryotic proteome by alternative splicing. Nature. 2010;463(7280):457–63.

  2. 2.

    Black DL. Mechanisms of alternative pre-messenger RNA splicing. Annu Rev Biochem. 2003;72:291–336.

  3. 3.

    Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, Mayr C, Kingsmore SF, Schroth GP, Burge CB. Alternative isoform regulation in human tissue transcriptomes. Nature. 2008;456(7221):470–6.

  4. 4.

    Pan Q, Shai O, Lee LJ, Frey BJ, Blencowe BJ. Deep surveying of alternative splicing complexity in the human transcriptome by high-throughput sequencing. Nat Genet. 2008;40(12):1413–5.

  5. 5.

    Grosso AR, Carmo-Fonseca M. The potential of targeting splicing for cancer therapy. In: Nuclear signaling pathways and targeting transcription in cancer. New York: Humana Press; 2014. p. 313–36.

  6. 6.

    Gallego-Paez LM, Bordone MC, Leote AC, Saraiva-Agostinho N, Ascensao-Ferreira M, Barbosa-Morais NL. Alternative splicing: the pledge, the turn, and the prestige: the key role of alternative splicing in human biological systems. Hum Genet. 2017;136(9):1015–42.

  7. 7.

    Sveen A, Kilpinen S, Ruusulehto A, Lothe RA, Skotheim RI. Aberrant RNA splicing in cancer; expression changes and driver mutations of splicing factor genes. Oncogene. 2016;35(19):2413–27.

  8. 8.

    Takehara T, Liu X, Fujimoto J, Friedman SL, Takahashi H. Expression and role of Bcl-xL in human hepatocellular carcinomas. Hepatology. 2001;34(1):55–61.

  9. 9.

    Olopade OI, Adeyanju MO, Safa AR, Hagos F, Mick R, Thompson CB, Recant WM. Overexpression of BCL-x protein in primary breast cancer is associated with high tumor grade and nodal metastases. Cancer J Sci Am. 1997;3(4):230–7.

  10. 10.

    Yae T, Tsuchihashi K, Ishimoto T, Motohara T, Yoshikawa M, Yoshida GJ, Wada T, Masuko T, Mogushi K, Tanaka H, et al. Alternative splicing of CD44 mRNA by ESRP1 enhances lung colonization of metastatic cancer cell. Nat Commun. 2012;3:883.

  11. 11.

    Brown RL, Reinke LM, Damerow MS, Perez D, Chodosh LA, Yang J, Cheng C. CD44 splice isoform switching in human and mouse epithelium is essential for epithelial-mesenchymal transition and breast cancer progression. J Clin Invest. 2011;121(3):1064–74.

  12. 12.

    Kahles A, Lehmann KV, Toussaint NC, Huser M, Stark SG, Sachsenberg T, Stegle O, Kohlbacher O, Sander C, Cancer Genome Atlas Research N, et al. Comprehensive Analysis of Alternative Splicing Across Tumors from 8,705 Patients. Cancer Cell. 2018;34(2):211–24 e216.

  13. 13.

    Zhu J, Chen Z, Yong L. Systematic profiling of alternative splicing signature reveals prognostic predictor for ovarian cancer. Gynecol Oncol. 2018;148(2):368–74.

  14. 14.

    Bjorklund SS, Panda A, Kumar S, Seiler M, Robinson D, Gheeya J, Yao M, Alnaes GIG, Toppmeyer D, Riis M, et al. Widespread alternative exon usage in clinically distinct subtypes of invasive ductal carcinoma. Sci Rep. 2017;7(1):5568.

  15. 15.

    Meliso FM, Hubert CG, Galante PAF, Penalva LO. RNA processing as an alternative route to attack glioblastoma. Hum Genet. 2017;136(9):1129–41.

  16. 16.

    Lin P, He RQ, Ma FC, Liang L, He Y, Yang H, Dang YW, Chen G. Systematic analysis of survival-associated alternative splicing signatures in gastrointestinal Pan-adenocarcinomas. Ebiomedicine. 2018;34:46–60.

  17. 17.

    Narayanan SP, Singh S, Shukla S. A saga of cancer epigenetics: linking epigenetics to alternative splicing. Biochem J. 2017;474(6):885–96.

  18. 18.

    Misquitta-Ali CM, Cheng E, O'Hanlon D, Liu N, McGlade CJ, Tsao MS, Blencowe BJ. Global profiling and molecular characterization of alternative splicing events misregulated in lung cancer. Mol Cell Biol. 2011;31(1):138–50.

  19. 19.

    Fu XD, Ares M Jr. Context-dependent control of alternative splicing by RNA-binding proteins. Nat Rev Genet. 2014;15(10):689–701.

  20. 20.

    Coltri PP, Dos Santos MGP, da Silva GHG. Splicing and cancer: challenges and opportunities. Wiley Interdiscip Rev RNA. 2019;10(3):e1527.

  21. 21.

    Zhou HL, Luo G, Wise JA, Lou H. Regulation of alternative splicing by local histone modifications: potential roles for RNA-guided mechanisms. Nucleic Acids Res. 2014;42(2):701–13.

  22. 22.

    Kornblihtt AR. Chromatin, transcript elongation and alternative splicing. Nat Struct Mol Biol. 2006;13(1):5–7.

  23. 23.

    Curado J, Iannone C, Tilgner H, Valcarcel J, Guigo R. Promoter-like epigenetic signatures in exons displaying cell type-specific splicing. Genome Biol. 2015;16:236.

  24. 24.

    Chodavarapu RK FS, Bernatavichute YV, Chen PY, Stroud H, Yu Y, Hetzel JA, Kuo F, Kim J, Cokus SJ, Casero D, Bernal M, Huijser P, Clark AT, Krämer U, Merchant SS, Zhang X, Jacobsen SE, Pellegrini M. Relationship between nucleosome positioning and DNA methylation. Nature. 2010;466(7304):388-92.

  25. 25.

    Anastasiadou C, Malousi A, Maglaveras N, Kouidou S. Human epigenome data reveal increased CpG methylation in alternatively spliced sites and putative exonic splicing enhancers. DNA Cell Biol. 2011;30(5):267–75.

  26. 26.

    Ahuvi Yearim SG, Shayevitch R, Melcer S, Glaich O, Mallm J-P, Nissim-Rafinia M, Cohen A-HS, Rippe K, Meshorer E, Ast G. HP1 Is Involved in Regulating the Global Impact of DNA Methylation on Alternative Splicing. Cell Rep. 2015;10(7):1122-34.

  27. 27.

    Ryan M, Wong WC, Brown R, Akbani R, Su X, Broom B, Melott J, Weinstein J. TCGASpliceSeq a compendium of alternative mRNA splicing in cancer. Nucleic Acids Res. 2016;44(D1):D1018–22.

  28. 28.

    Li S, Zhang J, Huang S, He X. Genome-wide analysis reveals that exon methylation facilitates its selective usage in the human transcriptome. Brief Bioinform. 2017;19(5):754-764.

  29. 29.

    Castle JC, Zhang C, Shah JK, Kulkarni AV, Kalsotra A, Cooper TA, Johnson JM. Expression of 24,426 human alternative splicing events and predicted cis regulation in 48 tissues and cell lines. Nat Genet. 2008;40(12):1416–25.

  30. 30.

    Yu GC, Wang LG, Han YY, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics-a J Integr Biol. 2012;16(5):284–7.

  31. 31.

    Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, Simonovic M, Doncheva NT, Morris JH, Bork P, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–13.

  32. 32.

    David CJ, Manley JL. Alternative pre-mRNA splicing regulation in cancer: pathways and programs unhinged. Genes Dev. 2010;24(21):2343–64.

  33. 33.

    Kornblihtt AR, Schor IE, Allo M, Dujardin G, Petrillo E, Munoz MJ. Alternative splicing: a pivotal step between eukaryotic transcription and translation. Nat Rev Mol Cell Biol. 2013;14(3):153–65.

  34. 34.

    Chen M, Manley JL. Mechanisms of alternative splicing regulation: insights from molecular and genomics approaches. Nat Rev Mol Cell Biol. 2009;10(11):741–54.

  35. 35.

    Singh RK, Cooper TA. Pre-mRNA splicing in disease and therapeutics. Trends Mol Med. 2012;18(8):472–82.

  36. 36.

    Laurent L, Wong E, Li G, Huynh T, Tsirigos A, Ong CT, Low HM, Kin Sung KW, Rigoutsos I, Loring J, et al. Dynamic changes in the human methylome during differentiation. Genome Res. 2010;20(3):320–31.

  37. 37.

    Shukla S, Kavak E, Gregory M, Imashimizu M, Shutinoski B, Kashlev M, Oberdoerffer P, Sandberg R, Oberdoerffer S. CTCF-promoted RNA polymerase II pausing links DNA methylation to splicing. Nature. 2011;479(7371):74–9.

  38. 38.

    Sugnet CW, Srinivasan K, Clark TA, O'Brien G, Cline MS, Wang H, Williams A, Kulp D, Blume JE, Haussler D, et al. Unusual intron conservation near tissue-regulated exons found by splicing microarrays. PLoS Comput Biol. 2006;2(1):22–35.

  39. 39.

    Ule J, Stefani G, Mele A, Ruggiu M, Wang XN, Taneri B, Gaasterland T, Blencowe BJ, Darnell RB. An RNA map predicting Nova-dependent splicing regulation. Nature. 2006;444(7119):580–6.

  40. 40.

    Fairbrother WG, Holste D, Burge CB, Sharp PA. Single nucleotide polymorphism-based validation of exonic splicing enhancers. PLoS Biol. 2004;2(9):1388–95.

  41. 41.

    Tsai YS, Dominguez D, Gomez SM, Wang Z. Transcriptome-wide identification and study of cancer-specific splicing events across multiple tumors. Oncotarget. 2015;6(9):6825–39.

  42. 42.

    Saghafinia S, Mina M, Riggi N, Hanahan D, Ciriello G. Pan-Cancer landscape of aberrant DNA methylation across human tumors. Cell Rep. 2018;25(4):1066–80 e1068.

  43. 43.

    Sehrawat B, Sridharan M, Ghosh S, Robson P, Cass CE, Mackey JR, Greiner R, Damaraju S. Potential novel candidate polymorphisms identified in genome-wide association study for breast cancer susceptibility. Hum Genet. 2011;130(4):529–37.

  44. 44.

    Osakabe A, Adachi F, Arimura Y, Maehara K, Ohkawa Y, Kurumizaka H. Influence of DNA methylation on positioning and DNA flexibility of nucleosomes with pericentric satellite DNA. Open Biol. 2015;5(10):150128.

  45. 45.

    Jimenez-Useche I, Ke J, Tian Y, Shim D, Howell SC, Qiu X, Yuan C. DNA methylation regulated nucleosome dynamics. Sci Rep. 2013;3:2121.

  46. 46.

    Maunakea AK, Chepelev I, Cui K, Zhao K. Intragenic DNA methylation modulates alternative splicing by recruiting MeCP2 to promote exon recognition. Cell Res. 2013;23(11):1256–69.

  47. 47.

    Gelfman S, Cohen N, Yearim A, Ast G. DNA-methylation effect on cotranscriptional splicing is dependent on GC architecture of the exon–intron structure. Genome Res. 2013;23(5):789–99.

  48. 48.

    Czarnecka AM, Kukwa W, Kornakiewicz A, Lian F, Szczylik C. Clinical and molecular prognostic and predictive biomarkers in clear cell renal cell cancer. Future Oncol. 2014;10(15):2493–508.

  49. 49.

    Maroto P, Rini B. Molecular biomarkers in advanced renal cell carcinoma. Clin Cancer Res. 2014;20(8):2060–71.

  50. 50.

    Lawrence MS, Stojanov P, Polak P, Kryukov GV, Cibulskis K, Sivachenko A, Carter SL, Stewart C, Mermel CH, Roberts SA, et al. Mutational heterogeneity in cancer and the search for new cancer-associated genes. Nature. 2013;499(7457):214–8.

Download references

Acknowledgments

The authors thank the TCGA Research Network for their public high-throughput data of diverse cancer types. Also, the authors thank TCGASplice Seq for their public data.

Funding

There is no funding in the present study.

Author information

YZ and XS designed this study. XS and YT analyzed and wrote the manuscript with assistance from JW and ZS. All authors read and approved the final manuscript.

Correspondence to Yimin Zhu.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1. The number of samples from TCGA SpliceSeq and TCGA we used in the analyses. Table S2. The number and distribution of cancer-specific AS events across different cancers. Table S3. The number and distribution of AS-associated CpGs across different cancers. Table S4. Correlation of cancer-specific AS and DNA methylation in tumor tissues and normal tissues. Table S5. GO biological processes that were enriched by genes harboring differentially regulated exon in ten human cancers. Table S6. KEGG pathway enrichment of genes with methylation-related cancer-specific AS events across ten human cancer types.

Additional file 2: Figure S1. (A) The proportion of hypomethylated and hypermethylated CpGs at cancer-specific AS exon boundaries which showed association with DNA methylation in different tumors. (B) The proportion of up-regulated and down-regulated cancer-specific AS exons in different tumors.

Additional file 3: Figure S2. Protein-protein interaction networks of genes with methylation-associated cancer-specific AS events in different cancers without statistical significance.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Sun, X., Tian, Y., Wang, J. et al. Genome-wide analysis reveals the association between alternative splicing and DNA methylation across human solid tumors. BMC Med Genomics 13, 4 (2020). https://doi.org/10.1186/s12920-019-0654-9

Download citation

Keywords

  • Alternative splicing
  • DNA methylation
  • Exon boundaries
  • Cancers
  • Association