Characterization of alternative spliceoforms and the RNA splicing machinery in pancreatic cancer

Spliceosomic dysregulation in pancreatic cancer uncovers splicing factors PRPF8 and RBMX as novel candidate actionable targets

Pancreatic ductal adenocarcinoma (PDAC) is a highly lethal cancer, characterized by late diagnosis and poor treatment response. Surgery is the only curative approach, only available to early-diagnosed patients. Current therapies have limited effects, cause severe toxicities, and minimally improve overall survival. Understanding of splicing machinery alterations in PDAC remains incomplete. Here, we comprehensively examined 59 splicing machinery components, uncovering dysregulation in pre-mRNA processing factor 8 (PRPF8) and RNA-binding motif protein X-linked (RBMX). Their downregulated expression was linked to poor prognosis and malignancy features, including tumor stage, invasion and metastasis, and associated with poorer survival and the mutation of key PDAC genes. Experimental modulation of these splicing factors in pancreatic cancer cell lines reverted their expression to non-tumor levels and resulted in decreased key tumor-related features. These results provide evidence that the splicing machinery is altered in PDAC, wherein PRPF8 and RBMX emerge as candidate actionable therapeutic targets.

RSRC2 Expression Inhibits Malignant Progression of Triple-Negative Breast Cancer by Transcriptionally Regulating SCIN Expression

Triple-negative breast cancer (TNBC) has a shorter survival time and higher mortality rate than other molecular subtypes. RSRC2 is a newly discovered tumor suppressor gene. However, the potential functional mechanism of RSRC2 in TNBC remains unknown so far. Multiple bioinformatics databases were used. A Human Transcriptome Array 2.0 analysis, ChIP-seq analysis, ChIP-qPCR, RT-qPCR, Western blot, cell function assays in vitro and a metastatic mouse model in vivo were performed to demonstrate the role of RSRC2 in TNBC. Through the analysis of various databases, RSRC2 expression was the lowest in TNBC tissues compared to other molecular subtypes. The low expression of RSRC2 was associated with a worse prognosis for patients with breast cancer. The transcriptome array, ChIP-seq and bioinformatics analysis identified that GRHL2 and SCIN might have a close relationship with RSRC2. The functional bioinformatics enrichment analysis and functional cell experiments showed that RSRC2 was involved in cell adhesion, cell proliferation, cell migration and invasion. Furthermore, RSRC2 expression suppressed SCIN expression but not GRHL2 expression. SCIN re-expression in the RSRC2 overexpression cells or SCIN knockdown in the RSRC2 knockdown cells reversed the cellular function caused by RSRC2. Mechanistically, RSRC2 transcriptionally inhibited SCIN expression. In summary, our study reveals that RSRC2 acts as a tumor suppressor in TNBC development and progression through negatively regulating SCIN-mediated cell function, thus providing a potential target for TNBC treatment.

Splicing alterations in pancreatic ductal adenocarcinoma: a new molecular landscape with translational potential

Pancreatic ductal adenocarcinoma (PDAC) remains one of the most lethal cancers worldwide, mainly due to its late diagnosis and lack of effective therapies, translating into a low 5-year 12% survival rate, despite extensive clinical efforts to improve outcomes. International cooperative studies have provided informative multiomic landscapes of PDAC, but translation of these discoveries into clinical advances are lagging. Likewise, early diagnosis biomarkers and new therapeutic tools are sorely needed to tackle this cancer. The study of poorly explored molecular processes, such as splicing, can provide new tools in this regard. Alternative splicing of pre-RNA allows the generation of multiple RNA variants from a single gene and thereby contributes to fundamental biological processes by finely tuning gene expression. However, alterations in alternative splicing are linked to many diseases, and particularly to cancer, where it can contribute to tumor initiation, progression, metastasis and drug resistance. Splicing defects are increasingly being associated with PDAC, including both mutations or dysregulation of components of the splicing machinery and associated factors, and altered expression of specific relevant gene variants. Such disruptions can be a key element enhancing pancreatic tumor progression or metastasis, while they can also provide suitable tools to identify potential candidate biomarkers and discover new actionable targets. In this review, we aimed to summarize the current information about dysregulation of splicing-related elements and aberrant splicing isoforms in PDAC, and to describe their relationship with the development, progression and/or aggressiveness of this dismal cancer, as well as their potential as therapeutic tools and targets.

RBFOX2 modulates a metastatic signature of alternative splicing in pancreatic cancer

Pancreatic ductal adenocarcinoma (PDA) is characterized by aggressive local invasion and metastatic spread, leading to high lethality. Although driver gene mutations during PDA progression are conserved, no specific mutation is correlated with the dissemination of metastases1-3. Here we analysed RNA splicing data of a large cohort of primary and metastatic PDA tumours to identify differentially spliced events that correlate with PDA progression. De novo motif analysis of these events detected enrichment of motifs with high similarity to the RBFOX2 motif. Overexpression of RBFOX2 in a patient-derived xenograft (PDX) metastatic PDA cell line drastically reduced the metastatic potential of these cells in vitro and in vivo, whereas depletion of RBFOX2 in primary pancreatic tumour cell lines increased the metastatic potential of these cells. These findings support the role of RBFOX2 as a potent metastatic suppressor in PDA. RNA-sequencing and splicing analysis of RBFOX2 target genes revealed enrichment of genes in the RHO GTPase pathways, suggesting a role of RBFOX2 splicing activity in cytoskeletal organization and focal adhesion formation. Modulation of RBFOX2-regulated splicing events, such as via myosin phosphatase RHO-interacting protein (MPRIP), is associated with PDA metastases, altered cytoskeletal organization and the induction of focal adhesion formation. Our results implicate the splicing-regulatory function of RBFOX2 as a tumour suppressor in PDA and suggest a therapeutic approach for metastatic PDA.

Dysregulated splicing factor SF3B1 unveils a dual therapeutic vulnerability to target pancreatic cancer cells and cancer stem cells with an anti-splicing drug

Background

Pancreatic ductal adenocarcinoma (PDAC) is a highly lethal cancer, requiring novel treatments to target both cancer cells and cancer stem cells (CSCs). Altered splicing is emerging as both a novel cancer hallmark and an attractive therapeutic target. The core splicing factor SF3B1 is heavily altered in cancer and can be inhibited by Pladienolide-B, but its actionability in PDAC is unknown. We explored the presence and role of SF3B1 in PDAC and interrogated its potential as an actionable target.

Methods

SF3B1 was analyzed in PDAC tissues, an RNA-seq dataset, and publicly available databases, examining associations with splicing alterations and key features/genes. Functional assays in PDAC cell lines and PDX-derived CSCs served to test Pladienolide-B treatment effects in vitro, and in vivo in zebrafish and mice.

Results

SF3B1 was overexpressed in human PDAC and associated with tumor grade and lymph-node involvement. SF3B1 levels closely associated with distinct splicing event profiles and expression of key PDAC players (KRAS, TP53). In PDAC cells, Pladienolide-B increased apoptosis and decreased multiple tumor-related features, including cell proliferation, migration, and colony/sphere formation, altering AKT and JNK signaling, and favoring proapoptotic splicing variants (BCL-XS/BCL-XL, KRASa/KRAS, Δ133TP53/TP53). Importantly, Pladienolide-B similarly impaired CSCs, reducing their stemness capacity and increasing their sensitivity to chemotherapy. Pladienolide-B also reduced PDAC/CSCs xenograft tumor growth in vivo in zebrafish and in mice.

Conclusion

SF3B1 overexpression represents a therapeutic vulnerability in PDAC, as altered splicing can be targeted with Pladienolide-B both in cancer cells and CSCs, paving the way for novel therapies for this lethal cancer.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13046-021-02153-9.

Alternative Splicing Events in Tumor Immune Infiltration in Colorectal Cancer

Despite extensive research, the exact mechanisms involved in colorectal cancer (CRC) etiology and pathogenesis remain unclear. This study aimed to examine the correlation between tumor-associated alternative splicing (AS) events and tumor immune infiltration (TII) in CRC. We analyzed transcriptome profiling and clinical CRC data from The Cancer Genome Atlas (TCGA) database and lists of AS-related and immune-related signatures from the SpliceSeq and Innate databases, respectively to develop and validate a risk model of differential AS events and subsequently a TII risk model. We then conducted a two-factor survival analysis to study the association between TII and AS risk and evaluated the associations between immune signatures and six types of immune cells based on the TIMER database. Subsequently, we studied the distribution of six types of TII cells in high- and low-risk groups for seven AS events and in total. We obtained the profiles of AS events/genes for 484 patients, which included 473 CRC tumor samples and 41 corresponding normal samples, and detected 22581 AS events in 8122 genes. Exon Skip (ES) (8446) and Mutually Exclusive Exons (ME) (74) exhibited the most and fewest AS events, respectively. We then classified the 433 patients with CRC into low-risk (n = 217) and high-risk (n = 216) groups based on the median risk score in different AS events. Compared with patients with low-risk scores (mortality = 11.8%), patients with high-risk scores were associated with poor overall survival (mortality = 27.6%). The risk score, cancer stage, and pathological stage (T, M, and N) were closely correlated with prognosis in patients with CRC (P < 0.001). We identified 6479 differentially expressed genes from the transcriptome profiles of CRC and intersected 468 differential immune-related signatures. High-AS-risk and high-TII-risk predicted a poor prognosis in CRC. Different AS types were associated with different TII risk characteristics. Alternate Acceptor site (AA) and Alternate Promoter (AP) events directly affected the concentration of CD4T cells, and the level of CD8T cells was closely correlated with Alternate Terminator (AT) and Exon Skip (ES) events. Thus, the concentration of CD4T and CD8T cells in the CRC immune microenvironment was not specifically modulated by AS. However, B cell, dendritic cell, macrophage, and neutrophilic cell levels were strongly correlated with AS events. These results indicate adverse associations between AS event risk levels and immune cell infiltration density. Taken together, our findings show a clear association between tumor-associated alternative splicing and immune cell infiltration events and patient outcome and could form a basis for the identification of novel markers and therapeutic targets for CRC and other cancers in the future.

Survival-Associated Alternative Messenger RNA Splicing Signatures in Pancreatic Ductal Adenocarcinoma: A Study Based on RNA-Sequencing Data

Multiple studies have shown that cancer-specific alternative splicing (AS) alterations are associated with clinical outcome. In this study, we aimed to profile prognostic AS signatures for pancreatic ductal adenocarcinoma (PDAC). We integrated the percent-spliced-in (PSI) data of AS in 140 PDAC patients based on the Cancer Genome Atlas (TCGA) dataset. We identified overall survival (OS)-associated AS events using univariate Cox regression analysis. Then, prognostic AS signatures were constructed for OS and chemoresistance prediction using the least absolute shrinkage and selection operator (LASSO) method. We also analyzed splicing factors (SFs) regulatory networks by Pearson's correlation. We detected 677 OS-related AS events in 485 genes by profiling 10,354 AS events obtained from 140 PDAC patients. Gene functional enrichment analysis demonstrated the pathways enriched by survival-associated AS. The AS signatures constructed with significant survival-associated AS events revealed high performance in predicting PDAC survival and gemcitabine chemoresistance. The area under the receiver operator characteristic curve was 0.937 in training cohort and 0.748 in validation cohort at 2000 days of OS. Furthermore, we identified prognostic SFs (e.g., ESRP1 and HNRNPC) to build the AS regulatory network. We constructed AS signatures for OS and gemcitabine chemoresistance in PDAC patients, which may provide clues for further experiment-based mechanism study.

Splice variants as novel targets in pancreatic ductal adenocarcinoma

Despite a wealth of genomic information, a comprehensive alternative splicing (AS) analysis of pancreatic ductal adenocarcinoma (PDAC) has not been performed yet. In the present study, we assessed whole exome-based transcriptome and AS profiles of 43 pancreas tissues using Affymetrix exon array. The AS analysis of PDAC indicated on average two AS probe-sets (ranging from 1-28) in 1,354 significantly identified protein-coding genes, with skipped exon and alternative first exon being the most frequently utilised. In addition to overrepresented extracellular matrix (ECM)-receptor interaction and focal adhesion that were also seen in transcriptome differential expression (DE) analysis, Fc gamma receptor-mediated phagocytosis and axon guidance AS genes were also highly represented. Of note, the highest numbers of AS probe-sets were found in collagen genes, which encode the characteristically abundant stroma seen in PDAC. We also describe a set of 37 'hypersensitive' genes which were frequently targeted by somatic mutations, copy number alterations, DE and AS, indicating their propensity for multidimensional regulation. We provide the most comprehensive overview of the AS landscape in PDAC with underlying changes in the spliceosomal machinery. We also collate a set of AS and DE genes encoding cell surface proteins, which present promising diagnostic and therapeutic targets in PDAC.

SF3B4 is decreased in pancreatic cancer and inhibits the growth and migration of cancer cells

Splicing factor 3b subunit 4, a critical component of pre-message RNA splicing complex, has been reported to play an important part in the tumorigenesis. However, the expression pattern and biological role of splicing factor 3b subunit 4 in pancreatic cancer have never been investigated. In this study, we found that both the messenger RNA ( p < 0.001) and protein level of splicing factor 3b subunit 4 were decreased significantly in pancreatic cancer specimens compared with their adjacent normal tissues. Overexpression of splicing factor 3b subunit 4 in pancreatic cancer cells inhibited cell growth and motility in vitro, while suppressing splicing factor 3b subunit 4 expression promoted the proliferation and migration of pancreatic cancer cells. In addition, splicing factor 3b subunit 4 was found to inhibit the activity of signal transducer and activator of transcription 3 signaling via downregulating the phosphorylation of signal transducer and activator of transcription 3 on a tyrosine residue at position 705. Taken together, these findings demonstrated that splicing factor 3b subunit 4 acted as a suppressive role in pancreatic cancer and indicated that restoring the function of splicing factor 3b subunit 4 might be a strategy for cancer therapy.

Genome-wide Profiling of RNA splicing in prostate tumor from RNA-seq data using virtual microarrays

Background

Second generation RNA sequencing technology (RNA-seq) offers the potential to interrogate genome-wide differential RNA splicing in cancer. However, since short RNA reads spanning spliced junctions cannot be mapped contiguously onto to the chromosomes, there is a need for methods to profile splicing from RNA-seq data. Before the invent of RNA-seq technologies, microarrays containing probe sequences representing exon-exon junctions of known genes have been used to hybridize cellular RNAs for measuring context-specific differential splicing. Here, we extend this approach to detect tumor-specific splicing in prostate cancer from a RNA-seq dataset.

Method

A database, SPEventH, representing probe sequences of under a million non-redundant splice events in human is created with exon-exon junctions of optimized length for use as virtual microarray. SPEventH is used to map tens of millions of reads from matched tumor-normal samples from ten individuals with prostate cancer. Differential counts of reads mapped to each event from tumor and matched normal is used to identify statistically significant tumor-specific splice events in prostate.

Results

We find sixty-one (61) splice events that are differentially expressed with a p-value of less than 0.0001 and a fold change of greater than 1.5 in prostate tumor compared to the respective matched normal samples. Interestingly, the only evidence, EST (BF372485), in the public database for one of the tumor-specific splice event joining one of the intron in KLK3 gene to an intron in KLK2, is also derived from prostate tumor-tissue. Also, the 765 events with a p-value of less than 0.001 is shown to cluster all twenty samples in a context-specific fashion with few exceptions stemming from low coverage of samples.

Conclusions

We demonstrate that virtual microarray experiments using a non-redundant database of splice events in human is both efficient and sensitive way to profile genome-wide splicing in biological samples and to detect tumor-specific splicing signatures in datasets from RNA-seq technologies. The signature from the large number of splice events that could cluster tumor and matched-normal samples into two tight separate clusters, suggests that differential splicing is yet another RNA phenotype, alongside gene expression and SNPs, that can be exploited for tumor stratification.