Transcriptome-wide identification and study of cancer-specific splicing events across multiple tumors

Transcriptomes of the tumor-adjacent normal tissues are more informative than tumors in predicting recurrence in colorectal cancer patients

BACKGROUND

Previous investigations of transcriptomic signatures of cancer patient survival and post-therapy relapse have focused on tumor tissue. In contrast, here we show that in colorectal cancer (CRC) transcriptomes derived from normal tissues adjacent to tumors (NATs) are better predictors of relapse.

RESULTS

Using the transcriptomes of paired tumor and NAT specimens from 80 Korean CRC patients retrospectively determined to be in recurrence or nonrecurrence states, we found that, when comparing recurrent with nonrecurrent samples, NATs exhibit a greater number of differentially expressed genes (DEGs) than tumors. Training two prognostic elastic net-based machine learning models-NAT-based and tumor-based in our Samsung Medical Center (SMC) cohort, we found that NAT-based model performed better in predicting the survival when the model was applied to the tumor-derived transcriptomes of an independent cohort of 450 COAD patients in TCGA. Furthermore, compositions of tumor-infiltrating immune cells in NATs were found to have better prognostic capability than in tumors. We also confirmed through Cox regression analysis that in both SMC-CRC as well as in TCGA-COAD cohorts, a greater proportion of genes exhibited significant hazard ratio when NAT-derived transcriptome was used compared to when tumor-derived transcriptome was used.

CONCLUSIONS

Taken together, our results strongly suggest that NAT-derived transcriptomes and immune cell composition of CRC are better predictors of patient survival and tumor recurrence than the primary tumor.

Broad misappropriation of developmental splicing profile by cancer in multiple organs

Oncogenesis mimics key aspects of embryonic development. However, the underlying mechanisms are incompletely understood. Here, we demonstrate that the splicing events specifically active during human organogenesis, are broadly reactivated in the organ-specific tumor. Such events are associated with key oncogenic processes and predict proliferation rates in cancer cell lines as well as patient survival. Such events preferentially target nitrosylation and transmembrane-region domains, whose coordinated splicing in multiple genes respectively affect intracellular transport and N-linked glycosylation. We infer critical splicing factors potentially regulating embryonic splicing events and show that such factors are potential oncogenic drivers and are upregulated specifically in malignant cells. Multiple complementary analyses point to MYC and FOXM1 as potential transcriptional regulators of critical splicing factors in brain and liver. Our study provides a comprehensive demonstration of a splicing-mediated link between development and cancer, and suggest anti-cancer targets including splicing events, and their upstream splicing and transcriptional regulators.

A widespread length-dependent splicing dysregulation in cancer

Dysregulation of alternative splicing is a key molecular hallmark of cancer. However, the common features and underlying mechanisms remain unclear. Here, we report an intriguing length-dependent splicing regulation in cancers. By systematically analyzing the transcriptome of thousands of cancer patients, we found that short exons are more likely to be mis-spliced and preferentially excluded in cancers. Compared to other exons, cancer-associated short exons (CASEs) are more conserved and likely to encode in-frame low-complexity peptides, with functional enrichment in GTPase regulators and cell adhesion. We developed a CASE-based panel as reliable cancer stratification markers and strong predictors for survival, which is clinically useful because the detection of short exon splicing is practical. Mechanistically, mis-splicing of CASEs is regulated by elevated transcription and alteration of certain RNA binding proteins in cancers. Our findings uncover a common feature of cancer-specific splicing dysregulation with important clinical implications in cancer diagnosis and therapies.

Genome-Wide Profiling of Alternative Splicing Signatures Associated with Prognosis and Immune Microenvironment of Hepatocellular Carcinoma

Background

The potential roles of alternative splicing (AS) in HCC remain unknown. This study aimed to identify AS signatures associated with the prognosis that influence the immune microenvironment of HCC.

Material/Methods

The SpliceSeq tool was employed for genome-wide profiling of 7 AS events in 361 HCC patients from The Cancer Genome Atlas (TCGA). A prognostic signature was built by integrating Cox regression and the least absolute shrinkage and selection operator (LASSO). The support vector machine (SVM) and receiver operating characteristic curve (ROC) were employed to analyze the AS events in the signatures to discriminate the immune microenvironment.

Results

There were 3546 AS events highly linked to the survival of patients with HCC. The AS signature could effectively stratify HCC patients. Clustering analysis revealed 3 different immune clusters characterized with significantly different prognoses and were significantly correlated with AS signatures. The AS events in the final prognostic signature classified the immune cluster with an average AUC of the ROC (0.88). Moreover, a potential regulatory network of splicing events in HCC is presented.

Conclusions

We established the prognostic signature based on AS, which can effectively stratify HCC patients and predict the immune subtypes. Moreover, novel RNA splicing patterns and splicing-regulatory networks involved in HCC were discovered.

SUVA: splicing site usage variation analysis from RNA-seq data reveals highly conserved complex splicing biomarkers in liver cancer

Most of the current alternative splicing (AS) analysis tools are powerless to analyse complex splicing. To address this, we developed SUVA (Splice sites Usage Variation Analysis) that decomposes complex splicing events into five types of splice junction pairs. By analysing real and simulated data, SUVA showed higher sensitivity and accuracy in detecting AS events than the compared methods. Notably, SUVA detected extensive complex AS events and screened out 69 highly conserved and dominant AS events associated with cancer. The cancer-associated complex AS events in FN1 and the co-regulated RNA-binding proteins were significantly correlated with patient survival.

Novel insights into clear cell renal cell carcinoma prognosis by comprehensive characterization of aberrant alternative splicing signature: a study based on large-scale sequencing data

Clear cell renal cell carcinoma (ccRCC) is the most common type with poor prognosis in kidney tumor. Growing evidence has indicated that aberrant alternative splicing (AS) events are efficacious signatures for tumor prognosis prediction and therapeutic targets. However, the detailed roles of AS events in ccRCC are largely unknown. In our study, level 3 RNA-seq data was acquired from The Cancer Genome Atlas dataset and corresponding AS profiles were detected with the assistance of SpliceSeq software. A total of 2100 aberrant survival-associated AS events were identified via differential expression and univariate cox regression analysis. The final prognostic panel formed by 17 specific events was developed by stepwise least absolute shrinkage and selection operator (LASSO) penalty, with the area under curve (AUC) values of receiver operator characteristic (ROC) curves keeping above 0.7 spanning 1 year to 5 years. And the results from functional enrichment analyses are unanimous that autophagy could be a potential mechanism of splicing regulation in ccRCC. Furthermore, splicing regulatory network was constructed via Spearman correlation between splicing factors and AS events. Finally, unsupervised clustering analysis revealed three clusters with distinct survival patterns, and associated with specific clinicopathological phenotypes. In overall, we developed a robust and individualized predictive model based on large-scale sequencing data. The identified AS events and splicing network may be valuable in deciphering the crucial posttranscriptional mechanisms on tumorigenesis of ccRCC.

Profiles of alternative splicing events in the diagnosis and prognosis of Gastric Cancer

Background: Gastric cancer (GC) is a heterogeneous disease, and alternative splicing (AS) is a powerful universal transcriptional regulatory mechanism that contributes to the occurrence and development of cancer. However, the systematic analysis of AS events in GC is lacking; therefore, further studies are needed. Methods: Genome-wide analysis of AS events was performed using RNA-Seq data to evaluate the difference between GC and adjacent tissues at the AS level. Prognostic signatures based on differentially expressed alternative splicing (DEAS) events and a correlation network between DEAS and genes were built. Results: We identified 48,141 AS events, of which 2325 showed differential expression patterns. The parental genes before DEAS events play an essential role in regulating GC-related processes such as ribosome (FDR < 0.0001) and thermogenesis (FDR = 0.0002). There were 76 survival-associated DEAS cases. Stratifying patients according to the percent spliced in index value of six types of splicing patterns formed significant Kaplan-Meier curves in the overall survival analysis. A prognostic feature based on DEAS performed well for stratification in patients with GC. Conclusion: The present study will enrich our understanding regarding the distinction of GC and provide a generous amount of biomarkers and potential targets for the treatment of GC.

Isoginkgetin derivative IP2 enhances the adaptive immune response against tumor antigens

The success of cancer immunotherapy relies on the induction of an immunoprotective response targeting tumor antigens (TAs) presented on MHC-I molecules. We demonstrated that the splicing inhibitor isoginkgetin and its water-soluble and non-toxic derivative IP2 act at the production stage of the pioneer translation products (PTPs). We showed that IP2 increases PTP-derived antigen presentation in cancer cells in vitro and impairs tumor growth in vivo. IP2 action is long-lasting and dependent on the CD8+ T cell response against TAs. We observed that the antigen repertoire displayed on MHC-I molecules at the surface of MCA205 fibrosarcoma is modified upon treatment with IP2. In particular, IP2 enhances the presentation of an exon-derived epitope from the tumor suppressor nischarin. The combination of IP2 with a peptide vaccine targeting the nischarin-derived epitope showed a synergistic antitumor effect in vivo. These findings identify the spliceosome as a druggable target for the development of epitope-based immunotherapies.

Roles and mechanisms of alternative splicing in cancer - implications for care

Removal of introns from messenger RNA precursors (pre-mRNA splicing) is an essential step for the expression of most eukaryotic genes. Alternative splicing enables the regulated generation of multiple mRNA and protein products from a single gene. Cancer cells have general as well as cancer type-specific and subtype-specific alterations in the splicing process that can have prognostic value and contribute to every hallmark of cancer progression, including cancer immune responses. These splicing alterations are often linked to the occurrence of cancer driver mutations in genes encoding either core components or regulators of the splicing machinery. Of therapeutic relevance, the transcriptomic landscape of cancer cells makes them particularly vulnerable to pharmacological inhibition of splicing. Small-molecule splicing modulators are currently in clinical trials and, in addition to splice site-switching antisense oligonucleotides, offer the promise of novel and personalized approaches to cancer treatment.

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

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.

Classification of glioma based on prognostic alternative splicing

BACKGROUND

Previously developed classifications of glioma have provided enormous advantages for the diagnosis and treatment of glioma. Although the role of alternative splicing (AS) in cancer, especially in glioma, has been validated, a comprehensive analysis of AS in glioma has not yet been conducted. In this study, we aimed at classifying glioma based on prognostic AS.

METHODS

Using the TCGA glioblastoma (GBM) and low-grade glioma (LGG) datasets, we analyzed prognostic splicing events. Consensus clustering analysis was conducted to classified glioma samples and correlation analysis was conducted to characterize regulatory network of splicing factors and splicing events.

RESULTS

We analyzed prognostic splicing events and proposed novel splicing classifications across pan-glioma samples (labeled pST1-7) and across GBM samples (labeled ST1-3). Distinct splicing profiles between GBM and LGG were observed, and the primary discriminator for the pan-glioma splicing classification was tumor grade. Subtype-specific splicing events were identified; one example is AS of zinc finger proteins, which is involved in glioma prognosis. Furthermore, correlation analysis of splicing factors and splicing events identified SNRPB and CELF2 as hub splicing factors that upregulated and downregulated oncogenic AS, respectively.

CONCLUSION

A comprehensive analysis of AS in glioma was conducted in this study, shedding new light on glioma heterogeneity and providing new insights into glioma diagnosis and treatment.

psichomics: graphical application for alternative splicing quantification and analysis

Alternative pre-mRNA splicing generates functionally distinct transcripts from the same gene and is involved in the control of multiple cellular processes, with its dysregulation being associated with a variety of pathologies. The advent of next-generation sequencing has enabled global studies of alternative splicing in different physiological and disease contexts. However, current bioinformatics tools for alternative splicing analysis from RNA-seq data are not user-friendly, disregard available exon-exon junction quantification or have limited downstream analysis features. To overcome such limitations, we have developed psichomics, an R package with an intuitive graphical interface for alternative splicing quantification and downstream dimensionality reduction, differential splicing and gene expression and survival analyses based on The Cancer Genome Atlas, the Genotype-Tissue Expression project, the Sequence Read Archive project and user-provided data. These integrative analyses can also incorporate clinical and molecular sample-associated features. We successfully used psichomics in a laptop to reveal alternative splicing signatures specific to stage I breast cancer and associated novel putative prognostic factors.

Genetic basis of functional variability in adhesion G protein-coupled receptors

The enormous sizes of adhesion G protein-coupled receptors (aGPCRs) go along with complex genomic exon-intron architectures giving rise to multiple mRNA variants. There is a need for a comprehensive catalog of aGPCR variants for proper evaluation of the complex functions of aGPCRs found in structural, in vitro and animal model studies. We used an established bioinformatics pipeline to extract, quantify and visualize mRNA variants of aGPCRs from deeply sequenced transcriptomes. Data analysis showed that aGPCRs have multiple transcription start sites even within introns and that tissue-specific splicing is frequent. On average, 19 significantly expressed transcript variants are derived from a given aGPCR gene. The domain architecture of the N terminus encoded by transcript variants often differs and N termini without or with an incomplete seven-helix transmembrane anchor as well as separate seven-helix transmembrane domains are frequently derived from aGPCR genes. Experimental analyses of selected aGPCR transcript variants revealed marked functional differences. Our analysis has an impact on a rational design of aGPCR constructs for structural analyses and gene-deficient mouse lines and provides new support for independent functions of both, the large N terminus and the transmembrane domain of aGPCRs.

Survival-associated alternative splicing signatures in esophageal carcinoma

Alternative splicing (AS), a major mechanism for the enhancement of transcriptome and proteome diversity, has been widely demonstrated to be involved in the full spectrum of oncogenic processes. High-throughput sequencing technology and the rapid accumulation of clinical data sets have provided an opportunity to systemically analyze the association between messenger RNA AS variants and patient clinical outcomes. Here, we compared differentially spliced AS transcripts between esophageal carcinoma (ESCA) and non-tumor tissues, profiled genome-wide survival-associated AS events in 87 patients with esophageal adenocarcinoma (EAC) and 79 patients with esophageal squamous cell carcinoma (ESCC) using The Cancer Genome Atlas (TCGA) RNA-seq data set, and constructed predictive models as well as splicing regulation networks by integrated bioinformatic analysis. A total of 2326 AS events in 1738 genes and 1812 AS events in 1360 genes were determined to be significantly associated with overall survival (OS) of patients in the EAC and ESCC cohorts, respectively, including some essential participants in the oncogenic process. The predictive model of each splice type performed reasonably well in distinguishing good and poor outcomes of patients with esophageal cancer, and values for the area under curve reached 0.942 and 0.815 in the EAC exon skip predictive model and the ESCC alternate acceptor site predictive model, respectively. The splicing regulation networks revealed an interesting correlation between survival-associated splicing factors and prognostic AS genes. In summary, we created prognostic models for patients with esophageal cancer based on AS signatures and constructed novel splicing correlation networks.

Alternative RNA Splicing as a Potential Major Source of Untapped Molecular Targets in Precision Oncology and Cancer Disparities

Studies of alternative RNA splicing (ARS) have the potential to provide an abundance of novel targets for development of new biomarkers and therapeutics in oncology, which will be necessary to improve outcomes for patients with cancer and mitigate cancer disparities. ARS, a key step in gene expression enabling individual genes to encode multiple proteins, is emerging as a major driver of abnormal phenotypic heterogeneity. Recent studies have begun to identify RNA splicing-related genetic and genomic variation in tumors, oncogenes dysregulated by ARS, RNA splice variants driving race-related cancer aggressiveness and drug response, spliceosome-dependent transformation, and RNA splicing-related immunogenic epitopes in cancer. In addition, recent studies have begun to identify and test, preclinically and clinically, approaches to modulate and exploit ARS for therapeutic application, including splice-switching oligonucleotides, small molecules targeting RNA splicing or RNA splice variants, and combination regimens with immunotherapies. Although ARS data hold such promise for precision oncology, inclusion of studies of ARS in translational and clinical cancer research remains limited. Technologic developments in sequencing and bioinformatics are being routinely incorporated into clinical oncology that permit investigation of clinically relevant ARS events, yet ARS remains largely overlooked either because of a lack of awareness within the clinical oncology community or perceived barriers to the technical complexity of analyzing ARS. This perspective aims to increase such awareness, propose immediate opportunities to improve identification and analysis of ARS, and call for bioinformaticians and cancer researchers to work together to address the urgent need to incorporate ARS into cancer biology and precision oncology.

Transcriptome-wide analysis of alternative mRNA splicing signature in the diagnosis and prognosis of stomach adenocarcinoma

Alternative mRNA splicing (AS) contributes greatly to expanding the diversity and function of the proteome. Increasing evidence has suggested that dysregulation of mRNA splicing may be associated with various types of cancer. In the present study, RNA sequencing data were used to investigate alterations to the global mRNA splicing landscape of cellular genes from 452 stomach adenocarcinoma (STAD) tissues available in The Cancer Genome Atlas. Seven types of AS events, including the profiles of exon skipping events, were analyzed using SpliceSeq software. A total of 60,754 AS events in 10,611 genes were detected, more than half of which were exon skipping events. The AS events were compared between 415 STAD tissues and 37 normal tissues, and 3,895 differentially spliced cancer-specific events were identified. In addition, the association of the AS events with the overall survival of 373 STAD patients was analyzed. Multivariate Cox regression analysis revealed that prognosis prediction models based on the AS events with clinical parameters had an excellent performance in predicting the survival of STAD patients. This study provides a comprehensive portrait of global changes in mRNA splicing signatures that occur in gastric cancer. These results allowed the identification of a core set of AS in gastric cancer and indicated that AS events may serve as prognostic indicators.

Therapy-induced stress response is associated with downregulation of pre-mRNA splicing in cancer cells

BACKGROUND

Abnormal pre-mRNA splicing regulation is common in cancer, but the effects of chemotherapy on this process remain unclear.

METHODS

To evaluate the effect of chemotherapy on slicing regulation, we performed meta-analyses of previously published transcriptomic, proteomic, phosphoproteomic, and secretome datasets. Our findings were verified by LC-MS/MS, western blotting, immunofluorescence, and FACS analyses of multiple cancer cell lines treated with cisplatin and pladienolide B.

RESULTS

Our results revealed that different types of chemotherapy lead to similar changes in alternative splicing by inducing intron retention in multiple genes. To determine the mechanism underlying this effect, we analyzed gene expression in 101 cell lines affected by ɣ-irradiation, hypoxia, and 10 various chemotherapeutic drugs. Strikingly, оnly genes involved in the cell cycle and pre-mRNA splicing regulation were changed in a similar manner in all 335 tested samples regardless of stress stimuli. We revealed significant downregulation of gene expression levels in these two pathways, which could be explained by the observed decrease in splicing efficiency and global intron retention. We showed that the levels of active spliceosomal proteins might be further post-translationally decreased by phosphorylation and export into the extracellular space. To further explore these bioinformatics findings, we performed proteomic analysis of cisplatin-treated ovarian cancer cells. Finally, we demonstrated that the splicing inhibitor pladienolide B impairs the cellular response to DNA damage and significantly increases the sensitivity of cancer cells to chemotherapy.

CONCLUSIONS

Decreased splicing efficiency and global intron retention is a novel stress response mechanism that may promote survival of malignant cells following therapy. We found that this mechanism can be inhibited by pladienolide B, which significantly increases the sensitivity of cancer cells to cisplatin which makes it a good candidate drug for improving the efficiency of cancer therapy.

SRSF10-mediated IL1RAP alternative splicing regulates cervical cancer oncogenesis via mIL1RAP-NF-κB-CD47 axis

High-risk human papillomavirus oncoproteins E6 and E7 are the major etiological factors of cervical cancer but are insufficient for malignant transformation of cervical cancer. Dysregulated alternative splicing, mainly ascribed to aberrant splicing factor levels and activities, contributes to most cancer hallmarks. However, do E6 and E7 regulate the expression of splicing factors? Does alternative splicing acts as an “accomplice” of E6E7 to promote cervical cancer progression? Here, we identified that the splicing factor SRSF10, which promotes tumorigenesis of cervix, was upregulated by E6E7 via E2F1 transcriptional activation. SRSF10 modulates the alternate terminator of interleukin-1 receptor accessory protein exon 13 to increase production of the membrane form of interleukin-1 receptor accessory protein. SRSF10-mediated mIL1RAP upregulates the expression of the “don’t eat me” signal CD47 to inhibit macrophage phagocytosis by promoting nuclear factor-κB activation, which is pivotal in inflammatory, immune, and tumorigenesis processes. Altogether, these data reveal a close relationship among HPV infection, alternative splicing and tumor immune evasion, and also suggests that the SRSF10-mIL1RAP-CD47 axis could be an attractive therapeutic target for the treatment of cervical cancer.

Computer Analysis of Glioma Transcriptome Profiling: Alternative Splicing Events

Here we present the analysis of alternative splicing events on an example of glioblastoma cell culture samples using a set of computer tools in combination with database integration. The gene expression profiles of glioblastoma were obtained from cell culture samples of primary glioblastoma which were isolated and processed for RNA extraction. Transcriptome profiling of normal brain samples and glioblastoma were done by Illumina sequencing. The significant differentially expressed exon-level probes and their corresponding genes were identified using a combination of the splicing index method. Previous studies indicated that tumor-specific alternative splicing is important in the regulation of gene expression and corresponding protein functions during cancer development. Multiple alternative splicing transcripts have been identified as progression markers, including generalized splicing abnormalities and tumor- and stage-specific events. We used a set of computer tools which were recently applied to analysis of gene expression in laboratory animals to study differential splicing events. We found 69 transcripts that are differentially alternatively spliced. Three cancer-associated genes were considered in detail, in particular: APP (amyloid beta precursor protein), CASC4 (cancer susceptibility candidate 4) and TP53. Such alternative splicing opens new perspectives for cancer research.

Alternative splicing in cancers: From aberrant regulation to new therapeutics

Alternative splicing is one of the most common mechanisms for gene regulation in humans, and plays a vital role to increase the complexity of functional proteins. In this article, we seek to provide a general review on the relationships between alternative splicing and tumorigenesis. We briefly introduce the basic rules for regulation of alternative splicing, and discuss recent advances on dynamic regulation of alternative splicing in cancers by highlighting the roles of a variety of RNA splicing factors in tumorigenesis. We further discuss several important questions regarding the splicing of long noncoding RNAs and back-splicing of circular RNAs in cancers. Finally, we discuss the current technologies that can be used to manipulate alternative splicing and serve as potential cancer treatment.

RNA processing as an alternative route to attack glioblastoma

Genomic analyses have become an important tool to identify new avenues for therapy. This is especially true for cancer types with extremely poor outcomes, since our lack of effective therapies offers no tangible clinical starting point to build upon. The highly malignant brain tumor glioblastoma (GBM) exemplifies such a refractory cancer, with only 15 month average patient survival. Analyses of several hundred GBM samples compiled by the TCGA (The Cancer Genome Atlas) have produced an extensive transcriptomic map, identified prevalent chromosomal alterations, and defined important driver mutations. Unfortunately, clinical trials based on these results have not yet delivered an improvement on outcome. It is, therefore, necessary to characterize other regulatory routes known for playing a role in tumor relapse and response to treatment. Alternative splicing affects more than 90% of the human coding genes and it is an important source for transcript variation and gene regulation. Mutations and alterations in splicing factors are highly prevalent in multiple cancers, demonstrating the potential for splicing to act as a tumor driver. As a result, numerous genes are expressed as cancer-specific splicing isoforms that are functionally distinct from the canonical isoforms found in normal tissue. These include genes that regulate cancer-critical pathways such as apoptosis, DNA repair, cell proliferation, and migration. Splicing defects can even induce genomic instability, a common characteristic of cancer, and a driver of tumor evolution. Importantly, components of the splicing machinery are targetable; multiple drugs can inhibit splicing factors or promote changes in splicing which could be exploited to begin improving clinical outcomes. Here, we review the current literature and present a case for exploring RNA processing as therapeutic route for the treatment of GBM.

Allelic imbalance of somatic mutations in cancer genomes and transcriptomes

Somatic mutations in cancer genomes often show allelic imbalance (AI) of mutation abundance between the genome and transcriptome, but there is not yet a systematic understanding of AI. In this study, we performed large-scale DNA and RNA AI analyses of >100,000 somatic mutations in >2,000 cancer specimens across five tumor types using the exome and transcriptome sequencing data of the Cancer Genome Atlas consortium. First, AI analysis of nonsense mutations and frameshift indels revealed that nonsense-mediated decay is typical in cancer genomes, and we identified the relationship between the extent of AI and the location of mutations in addition to the well-recognized 50-nt rules. Second, the AI with splice site mutations may reflect the extent of intron retention and is frequently observed in known tumor suppressor genes. For missense mutations, we observed that mutations frequently subject to AI are enriched to genes related to cancer, especially those of apoptosis and the extracellular matrix, and C:G > A:T transversions. Our results suggest that mutations in known cancer-related genes and their transcripts are subjected to different levels of transcriptional or posttranscriptional regulation compared to wildtype alleles and may add an additional regulatory layer to the functions of cancer-relevant genes.

Alternative splicing: the pledge, the turn, and the prestige : The key role of alternative splicing in human biological systems

Alternative pre-mRNA splicing is a tightly controlled process conducted by the spliceosome, with the assistance of several regulators, resulting in the expression of different transcript isoforms from the same gene and increasing both transcriptome and proteome complexity. The differences between alternative isoforms may be subtle but enough to change the function or localization of the translated proteins. A fine control of the isoform balance is, therefore, needed throughout developmental stages and adult tissues or physiological conditions and it does not come as a surprise that several diseases are caused by its deregulation. In this review, we aim to bring the splicing machinery on stage and raise the curtain on its mechanisms and regulation throughout several systems and tissues of the human body, from neurodevelopment to the interactions with the human microbiome. We discuss, on one hand, the essential role of alternative splicing in assuring tissue function, diversity, and swiftness of response in these systems or tissues, and on the other hand, what goes wrong when its regulatory mechanisms fail. We also focus on the possibilities that splicing modulation therapies open for the future of personalized medicine, along with the leading techniques in this field. The final act of the spliceosome, however, is yet to be fully revealed, as more knowledge is needed regarding the complex regulatory network that coordinates alternative splicing and how its dysfunction leads to disease.

Prognostic alternative mRNA splicing signature in non-small cell lung cancer

Alternative splicing provides a major mechanism to generate protein diversity. Increasing evidence suggests a link of dysregulation of splicing associated with cancer. Genome-wide alternative splicing profiling in lung cancer remains largely unstudied. We generated alternative splicing profiles in 491 lung adenocarcinoma (LUAD) and 471 lung squamous cell carcinoma (LUSC) patients in TCGA using RNA-seq data, prognostic models and splicing networks were built by integrated bioinformatics analysis. A total of 3691 and 2403 alternative splicing events were significantly associated with patient survival in LUAD and LUSC, respectively, including EGFR, CD44, PIK3C3, RRAS2, MAPKAP1 and FGFR2. The area under the curve of the receiver-operator characteristic curve for prognostic predictor in NSCLC was 0.817 at 2000 days of overall survival which were also over 0.8 in LUAD and LUSC, separately. Interestingly, splicing correlation networks uncovered opposite roles of splicing factors in LUAD and LUSC. We created prognostic predictors based on alternative splicing events with high performances for risk stratification in NSCLC patients and uncovered interesting splicing networks in LUAD and LUSC which could be underlying mechanisms.

A splicing isoform of TEAD4 attenuates the Hippo-YAP signalling to inhibit tumour proliferation

Aberrant splicing is frequently found in cancer, yet the biological consequences of such alterations are mostly undefined. Here we report that the Hippo–YAP signalling, a key pathway that regulates cell proliferation and organ size, is under control of a splicing switch. We show that TEAD4, the transcription factor that mediates Hippo–YAP signalling, undergoes alternative splicing facilitated by the tumour suppressor RBM4, producing a truncated isoform, TEAD4-S, which lacks an N-terminal DNA-binding domain, but maintains YAP interaction domain. TEAD4-S is located in both the nucleus and cytoplasm, acting as a dominant negative isoform to YAP activity. Consistently, TEAD4-S is reduced in cancer cells, and its re-expression suppresses cancer cell proliferation and migration, inhibiting tumour growth in xenograft mouse models. Furthermore, TEAD4-S is reduced in human cancers, and patients with elevated TEAD4-S levels have improved survival. Altogether, these data reveal a splicing switch that serves to fine tune the Hippo–YAP pathway.

The Hippo/Yap signalling pathway is found deregulated in several cancers. Here, the authors uncover an additional mechanism of YAP regulation that occurs via alternately spliced isoform of TEAD4, which acts as a dominant negative regulator of YAP-TEAD signalling.

Functional genomics analyses of RNA-binding proteins reveal the splicing regulator SNRPB as an oncogenic candidate in glioblastoma

BACKGROUND

Glioblastoma (GBM) is the most common and aggressive type of brain tumor. Currently, GBM has an extremely poor outcome and there is no effective treatment. In this context, genomic and transcriptomic analyses have become important tools to identify new avenues for therapies. RNA-binding proteins (RBPs) are master regulators of co- and post-transcriptional events; however, their role in GBM remains poorly understood. To further our knowledge of novel regulatory pathways that could contribute to gliomagenesis, we have conducted a systematic study of RBPs in GBM.

RESULTS

By measuring expression levels of 1542 human RBPs in GBM samples and glioma stem cell samples, we identified 58 consistently upregulated RBPs. Survival analysis revealed that increased expression of 21 RBPs was also associated with a poor prognosis. To assess the functional impact of those RBPs, we modulated their expression in GBM cell lines and performed viability, proliferation, and apoptosis assays. Combined results revealed a prominent oncogenic candidate, SNRPB, which encodes core spliceosome machinery components. To reveal the impact of SNRPB on splicing and gene expression, we performed its knockdown in a GBM cell line followed by RNA sequencing. We found that the affected genes were involved in RNA processing, DNA repair, and chromatin remodeling. Additionally, genes and pathways already associated with gliomagenesis, as well as a set of general cancer genes, also presented with splicing and expression alterations.

CONCLUSIONS

Our study provides new insights into how RBPs, and specifically SNRPB, regulate gene expression and directly impact GBM development.

SURVIV for survival analysis of mRNA isoform variation

The rapid accumulation of clinical RNA-seq data sets has provided the opportunity to associate mRNA isoform variations to clinical outcomes. Here we report a statistical method SURVIV (Survival analysis of mRNA Isoform Variation), designed for identifying mRNA isoform variation associated with patient survival time. A unique feature and major strength of SURVIV is that it models the measurement uncertainty of mRNA isoform ratio in RNA-seq data. Simulation studies suggest that SURVIV outperforms the conventional Cox regression survival analysis, especially for data sets with modest sequencing depth. We applied SURVIV to TCGA RNA-seq data of invasive ductal carcinoma as well as five additional cancer types. Alternative splicing-based survival predictors consistently outperform gene expression-based survival predictors, and the integration of clinical, gene expression and alternative splicing profiles leads to the best survival prediction. We anticipate that SURVIV will have broad utilities for analysing diverse types of mRNA isoform variation in large-scale clinical RNA-seq projects.

An extensive program of periodic alternative splicing linked to cell cycle progression

Progression through the mitotic cell cycle requires periodic regulation of gene function at the levels of transcription, translation, protein-protein interactions, post-translational modification and degradation. However, the role of alternative splicing (AS) in the temporal control of cell cycle is not well understood. By sequencing the human transcriptome through two continuous cell cycles, we identify ~1300 genes with cell cycle-dependent AS changes. These genes are significantly enriched in functions linked to cell cycle control, yet they do not significantly overlap genes subject to periodic changes in steady-state transcript levels. Many of the periodically spliced genes are controlled by the SR protein kinase CLK1, whose level undergoes cell cycle-dependent fluctuations via an auto-inhibitory circuit. Disruption of CLK1 causes pleiotropic cell cycle defects and loss of proliferation, whereas CLK1 over-expression is associated with various cancers. These results thus reveal a large program of CLK1-regulated periodic AS intimately associated with cell cycle control.

Functional divergence and convergence between the transcript network and gene network in lung adenocarcinoma

INTRODUCTION

Alternative RNA splicing is a critical regulatory mechanism during tumorigenesis. However, previous oncological studies mainly focused on the splicing of individual genes. Whether and how transcript isoforms are coordinated to affect cellular functions remain underexplored. Also of great interest is how the splicing regulome cooperates with the transcription regulome to facilitate tumorigenesis. The answers to these questions are of fundamental importance to cancer biology.

RESULTS

Here, we report a comparative study between the transcript-based network (TN) and the gene-based network (GN) derived from the transcriptomes of paired tumor-normal tissues from 77 lung adenocarcinoma patients. We demonstrate that the two networks differ significantly from each other in terms of patient clustering and the number and functions of network modules. Interestingly, the majority (89.5%) of multi-transcript genes have their transcript isoforms distributed in at least two TN modules, suggesting regulatory and functional divergences between transcript isoforms. Furthermore, TN and GN modules share onlŷ50%-60% of their biological functions. TN thus appears to constitute a regulatory layer separate from GN. Nevertheless, our results indicate that functional convergence and divergence both occur between TN and GN, implying complex interactions between the two regulatory layers. Finally, we report that the expression profiles of module members in both TN and GN shift dramatically yet concordantly during tumorigenesis. The mechanisms underlying this coordinated shifting remain unclear yet are worth further explorations.

CONCLUSION

We show that in lung adenocarcinoma, transcript isoforms per se are coordinately regulated to conduct biological functions not conveyed by the network of genes. However, the two networks may interact closely with each other by sharing the same or related biological functions. Unraveling the effects and mechanisms of such interactions will significantly advance our understanding of this deadly disease.