psichomics: graphical application for alternative splicing quantification and analysis

betAS: intuitive analysis and visualization of differential alternative splicing using beta distributions

Next-generation RNA sequencing allows alternative splicing (AS) quantification with unprecedented resolution, with the relative inclusion of an alternative sequence in transcripts being commonly quantified by the proportion of reads supporting it as percent spliced-in (PSI). However, PSI values do not incorporate information about precision, proportional to the respective AS events' read coverage. Beta distributions are suitable to quantify inclusion levels of alternative sequences, using reads supporting their inclusion and exclusion as surrogates for the two distribution shape parameters. Each such beta distribution has the PSI as its mean value and is narrower when the read coverage is higher, facilitating the interpretability of its precision when plotted. We herein introduce a computational pipeline, based on beta distributions accurately modeling PSI values and their precision, to quantitatively and visually compare AS between groups of samples. Our methodology includes a differential splicing significance metric that compromises the magnitude of intergroup differences, the estimation uncertainty in individual samples, and the intragroup variability, being therefore suitable for multiple-group comparisons. To make our approach accessible and clear to both noncomputational and computational biologists, we developed betAS, an interactive web app and user-friendly R package for visual and intuitive differential splicing analysis from read count data.

An alternative splicing signature defines the basal-like phenotype and predicts worse clinical outcome in pancreatic cancer

Pancreatic ductal adenocarcinoma (PDAC) is characterized by extremely poor prognosis. PDAC presents with molecularly distinct subtypes, with the basal-like one being associated with enhanced chemoresistance. Splicing dysregulation contributes to PDAC; however, its involvement in subtype specification remains elusive. Herein, we uncover a subtype-specific splicing signature associated with prognosis in PDAC and the splicing factor Quaking (QKI) as a determinant of the basal-like signature. Single-cell sequencing analyses highlight QKI as a marker of the basal-like phenotype. QKI represses splicing events associated with the classical subtype while promoting basal-like events associated with shorter survival. QKI favors a plastic, quasi-mesenchymal phenotype that supports migration and chemoresistance in PDAC organoids and cell lines, and its expression is elevated in high-grade primary tumors and metastatic lesions. These studies identify a splicing signature that defines PDAC subtypes and indicate that QKI promotes an undifferentiated, plastic phenotype, which renders PDAC cells chemoresistant and adaptable to environmental changes.

Dysregulation of RNA splicing in early non-alcoholic fatty liver disease through hepatocellular carcinoma

While changes in RNA splicing have been extensively studied in hepatocellular carcinoma (HCC), no studies have systematically investigated changes in RNA splicing during earlier liver disease. Mouse studies have shown that disruption of RNA splicing can trigger liver disease and we have shown that the splicing factor SRSF3 is decreased in the diseased human liver, so we profiled RNA splicing in liver samples from twenty-nine individuals with no-history of liver disease or varying degrees of non-alcoholic fatty liver disease (NAFLD). We compared our results with three publicly available transcriptome datasets that we re-analyzed for splicing events (SEs). We found many changes in SEs occurred during early liver disease, with fewer events occurring with the onset of inflammation and fibrosis. Many of these early SEs were enriched for SRSF3-dependent events and were associated with SRSF3 binding sites. Mapping the early and late changes to gene ontologies and pathways showed that the genes harboring these early SEs were involved in normal liver metabolism, whereas those harboring late SEs were involved in inflammation, fibrosis and proliferation. We compared the SEs with HCC data from the TCGA and observed that many of these early disease SEs are found in HCC samples and, furthermore, are correlated with disease survival. Changes in splicing factor expression are also observed, which may be associated with distinct subsets of the SEs. The maintenance of these SEs through the multi-year oncogenic process suggests that they may be causative. Understanding the role of these splice variants in metabolic liver disease progression may shed light on the triggers of liver disease progression and the pathogenesis of HCC.

NLRP3 leucine-rich repeats control induced and spontaneous inflammasome activation in cryopyrin-associated periodic syndrome

BACKGROUND

The cryopyrin-associated periodic syndromes (CAPS) comprise a group of rare autoinflammatory diseases caused by gain-of-function mutations in the NLRP3 gene. NLRP3 contains a leucine-rich repeats (LRR) domain with a highly conserved exonic organization that is subjected to extensive alternative splicing. Aberrant NLRP3 inflammasome assembly in CAPS causes chronic inflammation; however, the mechanisms regulating inflammasome function remain unclear.

OBJECTIVE

We aimed to elucidate the mechanisms regulating NLRP3-mediated autoinflammation in human disease, characterizing the role of LRR in inflammasome activation.

METHODS

We analyzed sequence read archive data to characterize the pattern of NLRP3 splicing in human monocytes and investigated the role of each LRR-coding exon in inflammasome regulation in genetically modified U937 cells representing CAPS and healthy conditions.

RESULTS

We detected a range of NLRP3 splice variants in human primary cells and monocytic cell lines, including 2 yet-undescribed splice variants. We observe that lipopolysaccharides affect the abundance of certain splice variants, suggesting that they may regulate NLRP3 activation by affecting alternative splicing. We showed that exons 4, 5, 7, and 9 are essential for inflammasome function, both in the context of wild-type NLRP3 activation by the agonist molecule nigericin and in a model of CAPS-mediated NLRP3 inflammasome assembly. Moreover, the SGT1-NLRP3 interaction is decreased in nonfunctional variants, suggesting that alternative splicing may regulate the recruitment of proteins that facilitate inflammasome assembly.

CONCLUSION

These findings demonstrate the contribution of the LRR domain in inflammasome function and suggest that navigating LRR exon usage within NLRP3 is sufficient to dampen inflammasome assembly in CAPS.

FOXA1 regulates alternative splicing in prostate cancer

Dysregulation of alternative splicing in prostate cancer is linked to transcriptional programs activated by AR, ERG, FOXA1, and MYC. Here, we show that FOXA1 functions as the primary orchestrator of alternative splicing dysregulation across 500 primary and metastatic prostate cancer transcriptomes. We demonstrate that FOXA1 binds to the regulatory regions of splicing-related genes, including HNRNPK and SRSF1. By controlling trans-acting factor expression, FOXA1 exploits an "exon definition" mechanism calibrating alternative splicing toward dominant isoform production. This regulation especially impacts splicing factors themselves and leads to a reduction of nonsense-mediated decay (NMD)-targeted isoforms. Inclusion of the NMD-determinant FLNA exon 30 by FOXA1-controlled oncogene SRSF1 promotes cell growth in vitro and predicts disease recurrence. Overall, we report a role for FOXA1 in rewiring the alternative splicing landscape in prostate cancer through a cascade of events from chromatin access, to splicing factor regulation, and, finally, to alternative splicing of exons influencing patient survival.

Hepatocyte Deletion of IGF2 Prevents DNA Damage and Tumor Formation in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is the fifth most common cancer worldwide. Serine-arginine rich splicing factor 3 (SRSF3) plays a critical role in hepatocyte function and its loss in mice promotes chronic liver damage and leads to HCC. Hepatocyte-specific SRSF3 knockout mice (SKO mice) also overexpress insulin-like growth factor 2 (IGF2). In the present study, double deletion of Igf2 and Srsf3 (DKO mice) prevents hepatic fibrosis and inflammation, and completely prevents tumor formation, and is associated with decreased proliferation, apoptosis and DNA damage, and restored DNA repair enzyme expression. This is confirmed in vitro, where IGF2 treatment of HepG2 hepatoma cells decreases DNA repair enzyme expression and causes DNA damage. Tumors from the SKO mice also show mutational signatures consistent with homologous recombination and mismatch repair defects. Analysis of frozen human samples shows that SRSF3 protein is decreased sixfold in HCC compared to normal liver tissue but SRSF3 mRNA is increased. Looking at public TCGA data, HCC patients having high SRSF3 mRNA expression show poor survival, as do patients with alterations in known SRSF3-dependent splicing events. The results indicate that IGF2 overexpression in conjunction with reduced SRSF3 splicing activity could be a major cause of DNA damage and driver of liver cancer.

DJExpress: An Integrated Application for Differential Splicing Analysis and Visualization

RNA-seq analysis of alternative pre-mRNA splicing has facilitated an unprecedented understanding of transcriptome complexity in health and disease. However, despite the availability of countless bioinformatic pipelines for transcriptome-wide splicing analysis, the use of these tools is often limited to expert bioinformaticians. The need for high computational power, combined with computational outputs that are complicated to visualize and interpret present obstacles to the broader research community. Here we introduce DJExpress, an R package for differential expression analysis of transcriptomic features and expression-trait associations. To determine gene-level differential junction usage as well as associations between junction expression and molecular/clinical features, DJExpress uses raw splice junction counts as input data. Importantly, DJExpress runs on an average laptop computer and provides a set of interactive and intuitive visualization formats. In contrast to most existing pipelines, DJExpress can handle both annotated and de novo identified splice junctions, thereby allowing the quantification of novel splice events. Moreover, DJExpress offers a web-compatible graphical interface allowing the analysis of user-provided data as well as the visualization of splice events within our custom database of differential junction expression in cancer (DJEC DB). DJEC DB includes not only healthy and tumor tissue junction expression data from TCGA and GTEx repositories but also cancer cell line data from the DepMap project. The integration of DepMap functional genomics data sets allows association of junction expression with molecular features such as gene dependencies and drug response profiles. This facilitates identification of cancer cell models for specific splicing alterations that can then be used for functional characterization in the lab. Thus, DJExpress represents a powerful and user-friendly tool for exploration of alternative splicing alterations in RNA-seq data, including multi-level data integration of alternative splicing signatures in healthy tissue, tumors and cancer cell lines.

Differential Expression of CREM/ICER Isoforms Is Associated with the Spontaneous Control of HIV Infection

A rare subset of HIV-infected individuals, termed elite controllers (ECs), can maintain long-term control over HIV replication in the absence of antiretroviral therapy (ART). To elucidate the biological mechanism of resistance to HIV replication at the molecular and cellular levels, we performed RNA sequencing and identified alternative splicing variants from ECs, HIV-infected individuals undergoing ART, ART-naive HIV-infected individuals, and healthy controls. We identified differential gene expression patterns that are specific to ECs and may influence HIV resistance, including alternative RNA splicing and exon usage variants of the CREM/ICER gene (cyclic AMP [cAMP]-responsive element modulator/inducible cAMP early repressors). The knockout and knockdown of specific ICER exons that were found to be upregulated in ECs resulted in significantly increased HIV infection in a CD4⁺ T cell line and primary CD4⁺ T cells. Overexpression of ICER isoforms decreased HIV infection in primary CD4⁺ T cells. Furthermore, ICER regulated HIV long terminal repeat (LTR) promoter activity in a Tat-dependent manner. Together, these results suggest that ICER is an HIV host factor that may contribute to the HIV resistance of ECs. These findings will help elucidate the mechanisms of HIV control by ECs and may yield a new approach for treatment of HIV. IMPORTANCE A small group of HIV-infected individuals, termed elite controllers (ECs), display control of HIV replication in the absence of antiretroviral therapy (ART). However, the mechanism of ECs' resistance to HIV replication is not clear. In our work, we found an increased expression of specific, small isoforms of ICER in ECs. Further experiments proved that ICER is a robust host factor to regulate viral replication. Furthermore, we found that ICER regulates HIV LTR promoter activity in a Tat-dependent manner. These findings suggest that ICER is related to spontaneous control of HIV infection in ECs. This study may help elucidate a novel target for treatment of HIV.

Effects of alternative splicing events and transcriptome changes on kidney stone formation

During the development of urinary stone disease, the formation of tiny crystals that adhere to the renal tubular epithelium induces epithelial cell damage. This damage and repair of the epithelium is associated with the establishment of more crystal adhesion sites, which in turn stimulates further crystal adhesion and, eventually, stone formation. Deposited crystals typically cause changes in epithelial cell gene expression, such as transcriptome changes and alternative splicing events. Although considered important for regulating gene expression, alternative splicing has not been reported in studies related to kidney stones. To date, whether alternative splicing events are involved in the regulation of stone formation and whether crystallographic cell interactions are regulated by alternative splicing at the transcriptional level have remained unknown. Therefore, we conducted RNA sequencing and alternative splicing-related bioassays by modeling the in vitro stone environment. Many alternative splicing events were associated with crystallographic cell interactions. Moreover, these events regulated transcription and significantly affected the capacity of crystals to adhere to renal tubular epithelial cells and regulate apoptosis.

The oncogenic kinase NEK2 regulates an RBFOX2-dependent pro-mesenchymal splicing program in triple-negative breast cancer cells

BACKGROUND

Triple-negative breast cancer (TNBC) is the most heterogeneous and malignant subtype of breast cancer (BC). TNBC is defined by the absence of expression of estrogen, progesterone and HER2 receptors and lacks efficacious targeted therapies. NEK2 is an oncogenic kinase that is significantly upregulated in TNBC, thereby representing a promising therapeutic target. NEK2 localizes in the nucleus and promotes oncogenic splice variants in different cancer cells. Notably, alternative splicing (AS) dysregulation has recently emerged as a featuring trait of TNBC that contributes to its aggressive phenotype.

METHODS

To investigate whether NEK2 modulates TNBC transcriptome we performed RNA-sequencing analyses in a representative TNBC cell line (MDA-MB-231) and results were validated in multiple TNBC cell lines. Bioinformatics and functional analyses were carried out to elucidate the mechanism of splicing regulation by NEK2. Data from The Cancer Genome Atlas were mined to evaluate the potential of NEK2-sensitive exons as markers to identify the TNBC subtype and to assess their prognostic value.

RESULTS

Transcriptome analysis revealed a widespread impact of NEK2 on the transcriptome of TNBC cells, with 1830 AS events that are susceptible to its expression. NEK2 regulates the inclusion of cassette exons in splice variants that discriminate TNBC from other BC and that correlate with poor prognosis, suggesting that this kinase contributes to the TNBC-specific splicing program. NEK2 elicits its effects by modulating the expression of the splicing factor RBFOX2, a well-known regulator of epithelial to mesenchymal transition (EMT). Accordingly, NEK2 splicing-regulated genes are enriched in functional terms related to cell adhesion and contractile cytoskeleton and NEK2 depletion in mesenchymal TNBC cells induces phenotypic and molecular traits typical of epithelial cells. Remarkably, depletion of select NEK2-sensitive splice-variants that are prognostic in TNBC patients is sufficient to interfere with TNBC cell morphology and motility, suggesting that NEK2 orchestrates a pro-mesenchymal splicing program that modulates migratory and invasive properties of TNBC cells.

CONCLUSIONS

Our study uncovers an extensive splicing program modulated by NEK2 involving splice variants that confer an invasive phenotype to TNBCs and that might represent, together with NEK2 itself, valuable therapeutic targets for this disease.

Transcript-targeted analysis reveals isoform alterations and double-hop fusions in breast cancer

Although transcriptome alteration is an essential driver of carcinogenesis, the effects of chromosomal structural alterations on the cancer transcriptome are not yet fully understood. Short-read transcript sequencing has prevented researchers from directly exploring full-length transcripts, forcing them to focus on individual splice sites. Here, we develop a pipeline for Multi-Sample long-read Transcriptome Assembly (MuSTA), which enables construction of a transcriptome from long-read sequence data. Using the constructed transcriptome as a reference, we analyze RNA extracted from 22 clinical breast cancer specimens. We identify a comprehensive set of subtype-specific and differentially used isoforms, which extended our knowledge of isoform regulation to unannotated isoforms including a short form TNS3. We also find that the exon-intron structure of fusion transcripts depends on their genomic context, and we identify double-hop fusion transcripts that are transcribed from complex structural rearrangements. For example, a double-hop fusion results in aberrant expression of an endogenous retroviral gene, ERVFRD-1, which is normally expressed exclusively in placenta and is thought to protect fetus from maternal rejection; expression is elevated in several TCGA samples with ERVFRD-1 fusions. Our analyses provide direct evidence that full-length transcript sequencing of clinical samples can add to our understanding of cancer biology and genomics in general.

MDM4 Isoform Expression in Melanoma Supports an Oncogenic Role for MDM4-A

The p53 tumor suppressor integrates upstream signals such as DNA damage and active oncogenes to initiate cell cycle arrest or apoptosis. This response is critical to halting inappropriate growth signals. As such, p53 activity is lost in cancer. In melanoma, however, the p53 gene is intact in a reported 94% of human cases. Rather than direct mutation, p53 is held inactive through interaction with inhibitory proteins. Here, we examine the expression of the two primary inhibitors of p53, MDM2 and MDM4, in genomic databases and biopsy specimens. We find that MDM4 is frequently overexpressed. Moreover, changes in splicing of MDM4 occur frequently and early in melanomagenesis. These changes in splicing must be considered in the design of therapeutic inhibitors of the MDM2/4 proteins for melanoma.

The scaffold protein IQGAP1 links heat-induced stress signals to alternative splicing regulation in gastric cancer cells

In response to oncogenic signals, Alternative Splicing (AS) regulators such as SR and hnRNP proteins show altered expression levels, subnuclear distribution and/or post-translational modification status, but the link between signals and these changes remains unknown. Here, we report that a cytosolic scaffold protein, IQGAP1, performs this task in response to heat-induced signals. We show that in gastric cancer cells, a nuclear pool of IQGAP1 acts as a tethering module for a group of spliceosome components, including hnRNPM, a splicing factor critical for the response of the spliceosome to heat-shock. IQGAP1 controls hnRNPM's sumoylation, subnuclear localisation and the relevant response of the AS machinery to heat-induced stress. Genome-wide analyses reveal that IQGAP1 and hnRNPM co-regulate the AS of a cell cycle-related RNA regulon in gastric cancer cells, thus favouring the accelerated proliferation phenotype of gastric cancer cells. Overall, we reveal a missing link between stress signals and AS regulation.

HIV-1 Infection Transcriptomics: Meta-Analysis of CD4+ T Cells Gene Expression Profiles

HIV-1 infection elicits a complex dynamic of the expression various host genes. High throughput sequencing added an expressive amount of information regarding HIV-1 infections and pathogenesis. RNA sequencing (RNA-Seq) is currently the tool of choice to investigate gene expression in a several range of experimental setting. This study aims at performing a meta-analysis of RNA-Seq expression profiles in samples of HIV-1 infected CD4+ T cells compared to uninfected cells to assess consistently differentially expressed genes in the context of HIV-1 infection. We selected two studies (22 samples: 15 experimentally infected and 7 mock-infected). We found 208 differentially expressed genes in infected cells when compared to uninfected/mock-infected cells. This result had moderate overlap when compared to previous studies of HIV-1 infection transcriptomics, but we identified 64 genes already known to interact with HIV-1 according to the HIV-1 Human Interaction Database. A gene ontology (GO) analysis revealed enrichment of several pathways involved in immune response, cell adhesion, cell migration, inflammation, apoptosis, Wnt, Notch and ERK/MAPK signaling.

The RNA-binding protein MEX3A is a prognostic factor and regulator of resistance to gemcitabine in pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is a highly aggressive cancer. Most patients present with advanced disease at diagnosis, which only permits palliative chemotherapeutic treatments. RNA dysregulation is a hallmark of most human cancers, including PDAC. To test the impact of RNA processing dysregulation on PDAC pathology, we performed a bioinformatics analysis to identify RNA-binding proteins (RBPs) associated with prognosis. Among the 12 RBPs associated with progression-free survival, we focused on MEX3A because it was recently shown to mark an intestinal stem cell population that is refractory to chemotherapeutic treatments, a typical feature of PDAC. Increased expression of MEX3A was correlated with higher disease stage in PDAC patients and with tumor development in a mouse model of PDAC. Depletion of MEX3A in PDAC cells enhanced sensitivity to chemotherapeutic treatment with gemcitabine, whereas its expression was increased in PDAC cells selected upon chronic exposure to the drug. RNA-sequencing analyses highlighted hundreds of genes whose expression is sensitive to MEX3A expression, with significant enrichment in cell cycle genes. MEX3A binds to its target mRNAs, like cyclin-dependent kinase 6 (CDK6), and promotes their stability. Accordingly, knockdown of MEX3A caused a significant reduction in PDAC cell proliferation and in progression to the S phase of the cell cycle. These findings uncover a novel role for MEX3A in the acquisition and maintenance of chemoresistance by PDAC cells, suggesting that it may represent a novel therapeutic target for PDAC.

Mutations primarily alter the inclusion of alternatively spliced exons

Genetic analyses and systematic mutagenesis have revealed that synonymous, non-synonymous and intronic mutations frequently alter the inclusion levels of alternatively spliced exons, consistent with the concept that altered splicing might be a common mechanism by which mutations cause disease. However, most exons expressed in any cell are highly-included in mature mRNAs. Here, by performing deep mutagenesis of highly-included exons and by analysing the association between genome sequence variation and exon inclusion across the transcriptome, we report that mutations only very rarely alter the inclusion of highly-included exons. This is true for both exonic and intronic mutations as well as for perturbations in trans. Therefore, mutations that affect splicing are not evenly distributed across primary transcripts but are focussed in and around alternatively spliced exons with intermediate inclusion levels. These results provide a resource for prioritising synonymous and other variants as disease-causing mutations.

The changing paradigm of intron retention: regulation, ramifications and recipes

Intron retention (IR) is a form of alternative splicing that has long been neglected in mammalian systems although it has been studied for decades in non-mammalian species such as plants, fungi, insects and viruses. It was generally assumed that mis-splicing, leading to the retention of introns, would have no physiological consequence other than reducing gene expression by nonsense-mediated decay. Relatively recent landmark discoveries have highlighted the pivotal role that IR serves in normal and disease-related human biology. Significant technical hurdles have been overcome, thereby enabling the robust detection and quantification of IR. Still, relatively little is known about the cis- and trans-acting modulators controlling this phenomenon. The fate of an intron to be, or not to be, retained in the mature transcript is the direct result of the influence exerted by numerous intrinsic and extrinsic factors at multiple levels of regulation. These factors have altered current biological paradigms and provided unexpected insights into the transcriptional landscape. In this review, we discuss the regulators of IR and methods to identify them. Our focus is primarily on mammals, however, we broaden the scope to non-mammalian organisms in which IR has been shown to be biologically relevant.

Right and left-sided colon cancers - specificity of molecular mechanisms in tumorigenesis and progression

Background

Given the differences in embryonic origin, vascular and nervous supplies, microbiotic burden, and main physiological functions of left and right colons, tumor location is increasingly suggested to dictate tumor behavior affecting pathology, progression and prognosis. Right-sided colon cancers arise in the cecum, ascending colon, hepatic flexure and/or transverse colon, while left-sided colon cancers arise in the splenic flexure, descending, and/or sigmoid colon. In contrast to prior reports, we attempt to delineate programs of tumorigenesis independently for each side.

Methods

Four hundred and eleven samples were extracted from The Cancer Genome Atlas-COAD cohort, based on a conservative sample inclusion criterion. Each side was independently analyzed with respect to their respective normal tissue, at the level of transcription, post-transcription, miRNA control and methylation in both a stage specific and stage-agnostic manner.

Results

Our results indicate a suppression of enzymes involved in various stages of carcinogen breakdown including CYP2C8, CYP4F12, GSTA1, and UGT1A within right colon tumors. This implies its reduced capacity to detoxify carcinogens, contributing to a genotoxic tumor environment, and subsequently a more aggressive phenotype. Additionally, we highlight a crucial nexus between calcium homeostasis (sensing, mobilization and absorption) and immune/GPCR signaling within left-sided tumors, possibly contributing to its reduced proliferative and metastatic potential. Interestingly, two genes SLC6A4 and HOXB13 show opposing regulatory trends within right and left tumors. Post-transcriptional regulation mediated by both RNA-binding proteins (e.g. NKRF (in left) and MSI2 (in right)) and miRNAs (e.g. miR-29a (in left); miR-155, miR181-d, miR-576 and miR23a (in right)) appear to exhibit side-specificity in control of their target transcripts and is pronounced in right colon tumors. Additionally, methylation results depict location-specific differences, with increased hypomethylation in open seas within left tumors, and increased hypermethylation of CpG islands within right tumors.

Conclusions

Differences in molecular mechanisms captured here highlight distinctions in tumorigenesis and progression between left and right colon tumors, which will serve as the basis for future studies, influencing the efficacies of existing and future diagnostic, prognostic and therapeutic interventions.

ASCOT identifies key regulators of neuronal subtype-specific splicing

Public archives of next-generation sequencing data are growing exponentially, but the difficulty of marshaling this data has led to its underutilization by scientists. Here, we present ASCOT, a resource that uses annotation-free methods to rapidly analyze and visualize splice variants across tens of thousands of bulk and single-cell data sets in the public archive. To demonstrate the utility of ASCOT, we identify novel cell type-specific alternative exons across the nervous system and leverage ENCODE and GTEx data sets to study the unique splicing of photoreceptors. We find that PTBP1 knockdown and MSI1 and PCBP2 overexpression are sufficient to activate many photoreceptor-specific exons in HepG2 liver cancer cells. This work demonstrates how large-scale analysis of public RNA-Seq data sets can yield key insights into cell type-specific control of RNA splicing and underscores the importance of considering both annotated and unannotated splicing events.

A stroma-corrected ZEB1 transcriptional signature is inversely associated with antitumor immune activity in breast cancer

The epithelial-to-mesenchymal transition (EMT) is an essential developmental process which can be hijacked by cancer cells, leading to enhanced metastasis and chemoresistance in experimental models. Recent studies have linked gene expression of EMT-associated gene signatures to increased inflammatory immune response in multiple cancer types. However, these studies did not account for the potential confounding effects of gene expression by tumor-infiltrating mesenchymal stromal cells. In this study, we comprehensively dissect the associations between multiple EMT transcription factors and EMT markers with stromal and immune tumor infiltration. We find that EMT-related genes are highly correlated with intratumoral stromal cell abundance and identify a specific relationship between stroma-corrected ZEB1 expression and decreased immune activity in multiple cancer types. We derive a stroma-corrected ZEB1-activated transcriptional signature and demonstrate that this signature includes several known inhibitors of inflammation, including BMPR2. Finally, multivariate survival analysis reveals that ZEB1 and its expression signature are significantly associated with reduced overall survival in breast cancer patients. In conclusion, this study identifies a novel association between stroma-adjusted ZEB1 expression and tumor immune activity and addresses the critical issue of confounding between EMT-associated genes and tumor stromal content.

Alternative splicing in lung cancer

Lung cancer has the highest mortality rate of all cancers worldwide. Lung cancer is a very heterogeneous disease that is often diagnosed at later stages which have a poor prognosis. Aberrant alternative splicing patterns found in lung cancer contribute to important cell functions. These include changes in splicing for the BCL2L1, MDM2, MDM4, NUMB and MET genes during lung tumourigenesis, to affect pathways involved in apoptosis, cell proliferation and cellular cohesion. Global analyses of RNASeq datasets suggest there may be many more potentially influential aberrant splicing events that need to be investigated in lung cancer. Changes in expression of the splicing factors that regulate alternative splicing events have also been identified in lung cancer. Of these, changes in expression of QKI, RBM4, RBM5, RBM6, RBM10 and SRSF1 proteins regulate many of the most frequently referenced aberrant splicing events in lung cancer. The expanding list of genes known to be aberrantly spliced in lung cancer along with the altered expression of splicing factors that regulate them are providing new clues as to how lung cancer develops, and how these events can be exploited for better treatment. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

Androgen-regulated transcription of ESRP2 drives alternative splicing patterns in prostate cancer

Prostate is the most frequent cancer in men. Prostate cancer progression is driven by androgen steroid hormones, and delayed by androgen deprivation therapy (ADT). Androgens control transcription by stimulating androgen receptor (AR) activity, yet also control pre-mRNA splicing through less clear mechanisms. Here we find androgens regulate splicing through AR-mediated transcriptional control of the epithelial-specific splicing regulator ESRP2. Both ESRP2 and its close paralog ESRP1 are highly expressed in primary prostate cancer. Androgen stimulation induces splicing switches in many endogenous ESRP2-controlled mRNA isoforms, including splicing switches correlating with disease progression. ESRP2 expression in clinical prostate cancer is repressed by ADT, which may thus inadvertently dampen epithelial splice programmes. Supporting this, treatment with the AR antagonist bicalutamide (Casodex) induced mesenchymal splicing patterns of genes including FLNB and CTNND1. Our data reveals a new mechanism of splicing control in prostate cancer with important implications for disease progression.