Tumor-specific isoform switch of the fibroblast growth factor receptor 2 underlies the mesenchymal and malignant phenotypes of clear cell renal cell carcinomas

FGFR families: biological functions and therapeutic interventions in tumors

There are five fibroblast growth factor receptors (FGFRs), namely, FGFR1-FGFR5. When FGFR binds to its ligand, namely, fibroblast growth factor (FGF), it dimerizes and autophosphorylates, thereby activating several key downstream pathways that play an important role in normal physiology, such as the Ras/Raf/mitogen-activated protein kinase kinase/extracellular signal-regulated kinase, phosphoinositide 3-kinase (PI3K)/AKT, phospholipase C gamma/diacylglycerol/protein kinase c, and signal transducer and activator of transcription pathways. Furthermore, as an oncogene, FGFR genetic alterations were found in 7.1% of tumors, and these alterations include gene amplification, gene mutations, gene fusions or rearrangements. Therefore, FGFR amplification, mutations, rearrangements, or fusions are considered as potential biomarkers of FGFR therapeutic response for tyrosine kinase inhibitors (TKIs). However, it is worth noting that with increased use, resistance to TKIs inevitably develops, such as the well-known gatekeeper mutations. Thus, overcoming the development of drug resistance becomes a serious problem. This review mainly outlines the FGFR family functions, related pathways, and therapeutic agents in tumors with the aim of obtaining better outcomes for cancer patients with FGFR changes. The information provided in this review may provide additional therapeutic ideas for tumor patients with FGFR abnormalities.

Immunomodulatory Gene-Splicing Dysregulation in Tumorigenesis: Unmasking the Complexity

Cancer is a global health concern with rising incidence, morbidity, and mortality. The interaction between the tumor and immune cells within the tumor microenvironment is facilitated by signaling pathways driven by immunomodulatory proteins. Alternative splicing regulates the production of multiple immunomodulatory proteins with diverse functionality from a single mRNA transcript. Splicing factors are pivotal in modulating alternative splicing processes but are also subject to regulation. The dysregulation of alternative splicing may result from splicing factor (SF) abnormal expression levels and mutations in the cis and trans-acting elements and small nuclear RNA (snRNA) molecules. Aberrant splicing may generate abnormal mRNA transcripts encoding isoforms with altered functions that contribute to tumorigenesis or cancer progression. This review uncovers the complexity of immunomodulatory genes splicing dysregulation in oncogenesis. Identifying specific immunomodulatory splicing isoforms that contribute to cancer could be utilized to improve current immunotherapeutic drugs or develop novel therapeutic interventions for cancer.

CanIsoNet: a database to study the functional impact of isoform switching events in diseases

Motivation

Alternative splicing, as an essential regulatory mechanism in normal mammalian cells, is frequently disturbed in cancer and other diseases. Switches in the expression of most dominant alternative isoforms can alter protein interaction networks of associated genes giving rise to disease and disease progression. Here, we present CanIsoNet, a database to view, browse and search isoform switching events in diseases. CanIsoNet is the first webserver that incorporates isoform expression data with STRING interaction networks and ClinVar annotations to predict the pathogenic impact of isoform switching events in various diseases.

Results

Data in CanIsoNet can be browsed by disease or searched by genes or isoforms in annotation-rich data tables. Various annotations for 11 811 isoforms and 14 357 unique isoform switching events across 31 different disease types are available. The network density score for each disease-specific isoform, PFAM domain IDs of disrupted interactions, domain structure visualization of transcripts and expression data of switched isoforms for each sample is given. Additionally, the genes annotated in ClinVar are highlighted in interactive interaction networks.

Availability and implementation

CanIsoNet is freely available at https://www.caniso.net. The source codes can be found under a Creative Common License at https://github.com/kahramanlab/CanIsoNet_Web.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

[Expression and significance of fibroblast growth factor receptor 2 in clear cell renal cell carcinoma]

OBJECTIVE

To investigate the expression of fibroblast growth factor receptor 2 (FGFR2) in clear cell renal cell carcinoma (ccRCC; or kidney renal clear cell carcinoma, KIRC), to analyze the relationship between the expression of FGFR2 and the clinical pathological features and prognosis of ccRCC, to study the relationship between the expression of FGFR2 and other molecules, and to explore its role in the development of ccRCC.

METHODS

Gene expressional and clinical information of ccRCC patients were downloaded from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus(GEO) database. Next, the data were transformed and collated. In the study, 104 clinical ccRCC samples and corresponding paracancerous normal tissue samples were collected from Department of Urology, Peking University First Hospital. Immunohistochemistry (IHC) was performed and the staining results were scored, so as to compare the expression of FGFR2 in ccRCC and paracancerous normal tissues. Besides, quantify real-time polymerase chain reaction (qRT-PCR) was used to detect the mRNA expression level of FGFR2 in normal renal epithelial cell lines (293) and ccRCC cell lines (786-O, 769-P, OSRC-2, Caki-1, ACHN, and A498). In addition, the relationship between FGFR2 expression and clinical pathological characteristics (including TNM staging and pathological grading) and survival prognosis in ccRCC patients was further analyzed. Furthermore, the relationship between FGFR2 expression and B cells, T cells, natural killer (NK) cells and neutrophil infiltration in the ccRCC patients was analyzed, and the Biological General Repository for Interactionh Datasets (BioGRID) was used to builds protein-protein interaction (PPI) networks to study molecules that interacted with the FGFR2 protein.

RESULTS

In the TCGA database, the expression of FGFR2 was down-regulated in ccRCC tissue samples compared with normal tissue samples, and the expression in the GEO database also showed this differences. Furthermore, FGFR2 expression was downregulated in ccRCC clinical samples and ccRCC cell lines, compared with corresponding paracancerous normal tissue or normal renal epithelial cell lines. In addition, FGFR2 high expression was associated with earlier, lower-level ccRCC and was associated with a better prognosis in the patients with ccRCC. Moreover, FGFR2 expression was not significantly related to B cells, T cells, NK cells and neutrophil infiltration, and the PPI network showed that FGFR2 protein interacted with certain molecules.

CONCLUSION

Our work sheds light on the potential role of FGFR2 in the development of ccRCC, suggesting that FGFR2 may serve as a prognostic marker and potential therapeutic target for patients with ccRCC.

The Role of Alternative Splicing in Cancer: Regulatory Mechanism, Therapeutic Strategy, and Bioinformatics Application

[Formula: see text] Alternative splicing (AS) can generate distinct transcripts and subsequent isoforms that play differential functions from the same pre-mRNA. Recently, increasing numbers of studies have emerged, unmasking the association between AS and cancer. In this review, we arranged AS events that are closely related to cancer progression and presented promising treatments based on AS for cancer therapy. Obtaining proliferative capacity, acquiring invasive properties, gaining angiogenic features, shifting metabolic ability, and getting immune escape inclination are all splicing events involved in biological processes. Spliceosome-targeted and antisense oligonucleotide technologies are two novel strategies that are hopeful in tumor therapy. In addition, bioinformatics applications based on AS were summarized for better prediction and elucidation of regulatory routines mingled in. Together, we aimed to provide a better understanding of complicated AS events associated with cancer biology and reveal AS a promising target of cancer treatment in the future.

Spatial expression of the FGFR2b splice isoform and its prognostic significance in endometrioid endometrial carcinoma

Endometrial carcinoma (EC) is the most common gynecological malignancy and fibroblast growth factor receptor 2 (FGFR2) is a frequently dysregulated receptor tyrosine kinase. FGFR2b and FGFR2c are the two main splice isoforms of FGFR2 and are normally localized in epithelial and mesenchymal cells, respectively. Previously, we demonstrated that FGFR2c mRNA expression was associated with aggressive tumor characteristics, shorter progression-free survival (PFS), and disease-specific survival (DSS) in endometrioid ECs (EECs). The objectives of this study were to investigate the spatial expression of FGFR2b in normal and hyperplasia with and without atypia of human endometrium and to assess the prognostic significance of FGFR2b expression in EC. FGFR2b and FGFR2c mRNA expression was evaluated in normal (proliferative [n = 10], secretory [n = 15], and atrophic [n = 10] endometrium), hyperplasia with and without atypia (n = 19) as well as two patient cohorts of EC samples (discovery [n = 78] and Vancouver [n = 460]) using isoform-specific BaseScope RNA in situ hybridization assays. Tumors were categorized based on FGFR2 isoform expression (one, both, or neither) and categories were correlated with clinicopathologic markers, molecular subtypes, and clinical outcomes. The FGFR2b splice isoform was exclusively expressed in the epithelial compartment of normal endometrium and hyperplasia without atypia. We observed FGFR2c expression at the basalis layer of glands in 33% (3/9) of hyperplasia with atypia. In patients with EEC, FGFR2b+/FGFR2c- expression was found in 48% of the discovery cohort and 35% of the validation Vancouver cohort. In univariate analyses, tumors with FGFR2b+/FGFR2c- expression had longer PFS (hazard ratio [HR] 0.265; 95% CI 0.145-0.423; log-rank p < 0.019) and DSS (HR 0.31; 95% CI 0.149-0.622; log-rank p < 0.001) compared to tumors with FGFR2b-/FGFR2c+ expression in the large EEC Vancouver cohort. In multivariable Cox regression analyses, tumors with FGFR2b+/FGFR2c- expression were significantly associated with longer DSS (HR 0.37; 95% CI 0.153-0.872; log-rank p < 0.023) compared to FGFR2b-/FGFR2c+ tumors. In conclusion, FGFR2b+/FGFR2c- expression is associated with favorable clinicopathologic markers and clinical outcomes suggesting that FGFR2b could play a role in tailoring the management of EEC patients in the clinic if these findings are confirmed in an independent cohort.

Abnormal GRHL2 Methylation Confers Malignant Progression to Acute Leukemia

Purpose

Abnormal methylation of Grainyhead-like 2 (GRHL2) is associated with a substantial role in the malignant phenotype of tumor patients. Our present research is aimed at studying the abnormal expression of GRHL2 and the association of methylation in patients with acute leukemia and its relationship with prognosis.

Materials and Methods

We used quantitative real-time polymerase chain reaction (qRT-PCR) for detecting the aberrant expression level of GRHL2 in 60 patients with acute leukemia and 60 normal controls. We analyzed the significant correlation between the expression level of GRHL2 with clinicopathological features and patients' prognosis in acute leukemia using the corresponding statistical methods. Secondly, we employed qRT-PCR and Western blotting to detect the mRNA and protein levels of GRHL2 in leukemia cell lines. Next, we used methylation-specific polymerase chain reaction (MSP) technology for detecting the methylation of GRHL2 in clinical samples with acute leukemia and cell lines. Then we investigated the demethylating effect of arsenic trioxide and 5-azacitidine on the mRNA and protein expression levels of GRHL2 in cell lines of acute leukemia. Finally, we studied the effects of arsenide trioxide and 5-azacitidine on the proliferation of leukemia cells and the TGF-β signaling pathway.

Results

We found a lower level of GRHL2 expression not only in acute leukemia patients but also in cell lines when compared with normal controls. At the same time, the expression level of GRHL2 in patients with acute leukemia was significantly correlated with leukocyte count, platelet count, and cytogenetic risk grouping. In addition, the lower GRHL2 expression group showed a significantly lower overall survival rate in acute leukemia patients than that of patients with a higher GRHL2 expression group. Univariate and multivariate analyses revealed that the expression of GRHL2 is an independent risk factor in acute leukemia patients. The methylation level of the GRHL2 promoter region in acute leukemia patients and cell lines was significantly higher than the normal control group, and we found the elevated mRNA and protein levels of GRHL2 in acute leukemia cell lines after the use of the demethylation drug arsenic trioxide and 5-azacitidine. At the same time, arsenide trioxide and 5-azacitidine are associated with the inhibition of cellular proliferation of acute leukemia cells and also promote the elevated expression of TGF-β signaling pathway-linked proteins, including TGF-β, Smad2, Smad3, and Smad4.

Conclusion

Increased expression and methylation level of GRHL2 are closely associated with the prognosis and malignant phenotype of acute leukemia patients and play an irreplaceable role in the occurrence and development of patients with acute leukemia.

A repositioning screen using an FGFR2 splicing reporter reveals compounds that regulate epithelial-mesenchymal transitions and inhibit growth of prostate cancer xenografts

Research in the area of hallmarks of cancer has opened the possibility of designing new therapies based on modulating these cancer properties. We present here a screen designed to find chemicals that modulate epithelial-mesenchymal transitions (EMTs) in prostate cancer. For screening, we used a repurposing library and, as a readout, an FGFR2-based splicing reporter, which has been shown previously to be a sensor for EMTs. Various properties of cancer cells were assessed, signaling pathways investigated, and in vivo experiments in nude mice xenografts performed. The screen yielded three hit compounds (a T-type Ca channel inhibitor, an L-type Ca channel inhibitor, and an opioid antagonist) that switch FGFR2 splicing and induce an epithelial phenotype in prostate cancer cells. The compounds affected differently various properties of cancer cells, but all of them decreased cell migration, which is in line with modulating EMTs. We further present mechanistic insights into one of the compounds, nemadipine-A. The administration of nemadipine-A intraperitoneally in a nude mouse xenograft model of prostate cancer slowed tumor growth. To conclude, we show that knowledge of the molecular mechanisms that connect alternative splicing and various cancer properties may be used as a platform for drug development.

METTL14-mediated Lnc-LSG1 m6A modification inhibits clear cell renal cell carcinoma metastasis via regulating ESRP2 ubiquitination

Clear cell renal cell carcinoma (ccRCC) is the most lethal urological cancer and is characterized by a high rate of metastasis and relapse. N6-Methyladenosine (m⁶A) is implicated in various stages of cancer development. However, a thorough understanding of m⁶A-modified lncRNAs in ccRCC is lacking. The results showed that METTL14 had decreased expression in ccRCC tissues. In addition, the expression of METTL14 was negatively correlated to the prognosis, stage, and ccRCC tumor grade. The silencing of METTL14 was shown to significantly increase metastasis in vitro and in vivo. High-throughput methylated RNA immunoprecipitation sequencing (MeRIP-seq) showed that the m⁶A levels of Lnc-LSG1 could be regulated by METTL14. Lnc-LSG1 can directly bind to ESRP2 protein and promote ESRP2 degradation via facilitating ESRP2 ubiquitination. However, m⁶A modification on Lnc-LSG1 can block the interaction between Lnc-LSG1 and ESRP2 via the m⁶A reader, YTHDC1. Taken together, our findings unraveled the novel mechanism of METTL14 inhibiting ccRCC progression, and explored the correlation between m⁶A and lncRNA in ccRCC for the first time.

Alternative RNA splicing in tumour heterogeneity, plasticity and therapy

Alternative splicing is a process by which a single gene is able to encode multiple different protein isoforms. It is regulated by the inclusion or exclusion of introns and exons that are joined in different patterns prior to protein translation, thus enabling transcriptomic and proteomic diversity. It is now widely accepted that alternative splicing is dysregulated across nearly all cancer types. This widespread dysregulation means that nearly all cellular processes are affected - these include processes synonymous with the hallmarks of cancer - evasion of apoptosis, tissue invasion and metastasis, altered cellular metabolism, genome instability and drug resistance. Emerging evidence indicates that the dysregulation of alternative splicing also promotes a permissive environment for increased tumour heterogeneity and cellular plasticity. These are fundamental regulators of a patient's response to therapy. In this Review, we introduce the mechanisms of alternative splicing and the role of aberrant splicing in cancer, with particular focus on newfound evidence of alternative splicing promoting tumour heterogeneity, cellular plasticity and altered metabolism. We discuss recent in vivo models generated to study alternative splicing and the importance of these for understanding complex tumourigenic processes. Finally, we review the effects of alternative splicing on immune evasion, cell death and genome instability, and how targeting these might enhance therapeutic efficacy.

Regulation of Alternative Splicing during Epithelial-Mesenchymal Transition

Alternative splicing is an essential mechanism of gene regulation, giving rise to remarkable protein diversity in higher eukaryotes. Epithelial-mesenchymal transition (EMT) is a developmental process that plays an essential role in metazoan embryogenesis. Recent studies have revealed that alternative splicing serves as a fundamental layer of regulation that governs cells to undergo EMT. In this review, we summarize recent findings on the functional impact of alternative splicing in EMT and EMT-associated activities. We then discuss the regulatory mechanisms that control alternative splicing changes during EMT.

Alternative RNA splicing in stem cells and cancer stem cells: Importance of transcript-based expression analysis

Alternative ribonucleic acid (RNA) splicing can lead to the assembly of different protein isoforms with distinctive functions. The outcome of alternative splicing (AS) can result in a complete loss of function or the acquisition of new functions. There is a gap in knowledge of abnormal RNA splice variants promoting cancer stem cells (CSCs), and their prospective contribution in cancer progression. AS directly regulates the self-renewal features of stem cells (SCs) and stem-like cancer cells. Notably, octamer-binding transcription factor 4A spliced variant of octamer-binding transcription factor 4 contributes to maintaining stemness properties in both SCs and CSCs. The epithelial to mesenchymal transition pathway regulates the AS events in CSCs to maintain stemness. The alternative spliced variants of CSCs markers, including cluster of differentiation 44, aldehyde dehydrogenase, and doublecortin-like kinase, α6β1 integrin, have pivotal roles in increasing self-renewal properties and maintaining the pluripotency of CSCs. Various splicing analysis tools are considered in this study. LeafCutter software can be considered as the best tool for differential splicing analysis and identification of the type of splicing events. Additionally, LeafCutter can be used for efficient mapping splicing quantitative trait loci. Altogether, the accumulating evidence re-enforces the fact that gene and protein expression need to be investigated in parallel with alternative splice variants.

Emerging Multi-cancer Regulatory Role of ESRP1: Orchestration of Alternative Splicing to Control EMT

RNA binding proteins (RBPs) associate with nascent and mature RNAs to perform biological functions such as alternative splicing and RNA stability. Having unique RNA recognition binding motifs, RBPs form complexes with RNA in a sequence- and structure-based manner. Aberrant expressions of several RBPs have been identified in tumorigenesis and cancer progression. These uncontrolled RBPs affect several mechanisms, including cell proliferation, tumor growth, invasion, metastasis and chemoresistance. Epithelial splicing regulatory protein 1 (ESRP1) is a member of the hnRNP family of proteins that play a crucial role in regulating numerous cellular processes, including alternative splicing and translation of multiple genes during organogenesis. Abnormal expression of ESRP1 alters the cell morphology, and leads to cell proliferation and tumor growth during cancer progression. ESRP1 mediated alternative splicing of target genes, including CD44, FGFR, PTBP1, LYN, ENAH, SPAG1 and ZMYND8, results in cancer progression. In addition, ESRP1 also regulates circularization and biogenesis of circular RNAs such as circUHRF1, circNOL10 and circANKS1B, whose expressions have been identified as key factors in various cancers. This multi-functional protein is also involved in imposing stability of target mRNAs such as cyclin A2, and thereby cell cycle regulation. The scope of this review is to examine recent scientific data, outcomes of the up- and down-regulated proteins, and the role of ESRP1 in various cancers. We conclude by summarizing ESRP1 dysregulation and its consequences on target genes in various human cancers. Collectively, the consequences of ESRP1 mediated splicing in cancer cells suggest the role of ESRP1 in cell proliferation and chemoresistance via apoptosis and autophagy modulation, which could, therefore, be potential targets for cancer therapeutics.

2b or Not 2b: How Opposing FGF Receptor Splice Variants Are Blocking Progress in Precision Oncology

More than ten thousand peer-reviewed studies have assessed the role of fibroblast growth factors (FGFs) and their receptors (FGFRs) in cancer, but few patients have yet benefited from drugs targeting this molecular family. Strategizing how best to use FGFR-targeted drugs is complicated by multiple variables, including RNA splicing events that alter the affinity of ligands for FGFRs and hence change the outcomes of stromal-epithelial interactions. The effects of splicing are most relevant to FGFR2; expression of the FGFR2b splice isoform can restore apoptotic sensitivity to cancer cells, whereas switching to FGFR2c may drive tumor progression by triggering epithelial-mesenchymal transition. The differentiating and regulatory actions of wild-type FGFR2b contrast with the proliferative actions of FGFR1 and FGFR3, and may be converted to mitogenicity either by splice switching or by silencing of tumor suppressor genes such as CDH1 or PTEN. Exclusive use of small-molecule pan-FGFR inhibitors may thus cause nonselective blockade of FGFR2 isoforms with opposing actions, undermining the rationale of FGFR2 drug targeting. This splice-dependent ability of FGFR2 to switch between tumor-suppressing and -driving functions highlights an unmet oncologic need for isoform-specific drug targeting, e.g., by antibody inhibition of ligand-FGFR2c binding, as well as for more nuanced molecular pathology prediction of FGFR2 actions in different stromal-tumor contexts.

Epithelial and Mesenchymal Markers in Adrenocortical Tissues: How Mesenchymal Are Adrenocortical Tissues?

A clinically relevant proportion of adrenocortical carcinoma (ACC) cases shows a tendency to metastatic spread. The objective was to determine whether the epithelial to mesenchymal transition (EMT), a mechanism associated with metastasizing in several epithelial cancers, might play a crucial role in ACC. 138 ACC, 29 adrenocortical adenomas (ACA), three normal adrenal glands (NAG), and control tissue samples were assessed for the expression of epithelial (E-cadherin and EpCAM) and mesenchymal (N-cadherin, SLUG and SNAIL) markers by immunohistochemistry. Using real-time RT-PCR we quantified the alternative isoform splicing of FGFR 2 and 3, another known indicator of EMT. We also assessed the impact of these markers on clinical outcome. Results show that both normal and neoplastic adrenocortical tissues lacked expression of epithelial markers but strongly expressed mesenchymal markers N-cadherin and SLUG. FGFR isoform splicing confirmed higher similarity of adrenocortical tissues to mesenchymal compared to epithelial tissues. In ACC, higher SLUG expression was associated with clinical markers indicating aggressiveness, while N-cadherin expression inversely associated with these markers. In conclusion, we could not find any indication of EMT as all adrenocortical tissues lacked expression of epithelial markers and exhibited closer similarity to mesenchymal tissues. However, while N-cadherin might play a positive role in tissue structure upkeep, SLUG seems to be associated with a more aggressive phenotype.

Alternative splicing modulates cancer aggressiveness: role in EMT/metastasis and chemoresistance

Enhanced metastasis and disease recurrence accounts for the high mortality rates associated with cancer. The process of Epithelial-Mesenchymal Transition (EMT) contributes towards the augmentation of cancer invasiveness along with the gain of stem-like and the subsequent drug-resistant behavior. Apart from the well-established transcriptional regulation, EMT is also controlled post-transcriptionally by virtue of alternative splicing (AS). Numerous genes including Fibroblast Growth Factor receptor (FGFR) as well as CD44 are differentially spliced during this trans-differentiation process which, in turn, governs cancer progression. These splicing alterations are controlled by various splicing factors including ESRP, RBFOX2 as well as hnRNPs. Here, we have depicted the mechanisms governing the splice isoform switching of FGFR and CD44. Moreover, the role of the splice variants generated by AS of these gene transcripts in modulating the metastatic potential and stem-like/chemoresistant behavior of cancer cells has also been highlighted. Additionally, the involvement of splicing factors in regulating EMT/invasiveness along with drug-resistance as well as the metabolic properties of the cells has been emphasized. Tumorigenesis is accompanied by a remodeling of the cellular splicing profile generating diverse protein isoforms which, in turn, control the cancer-associated hallmarks. Therefore, we have also briefly discussed about a wide variety of genes which are differentially spliced in the tumor cells and promote cancer progression. We have also outlined different strategies for targeting the tumor-associated splicing events which have shown promising results and therefore this approach might be useful in developing therapies to reduce cancer aggressiveness in a more specific manner.

Deletion of Yy1 in mouse lung epithelium unveils molecular mechanisms governing pleuropulmonary blastoma pathogenesis

Pleuropulmonary blastoma (PPB) is a very rare pediatric lung disease. It can progress from abnormal epithelial cysts to an aggressive sarcoma with poor survival. PPB is difficult to diagnose as it can be confounded with other cystic lung disorders, such as congenital pulmonary airway malformation (CPAM). PPB is associated with mutations in DICER1 that perturb the microRNA (miRNA) profile in lung. How DICER1 and miRNAs act during PPB pathogenesis remains unsolved. Lung epithelial deletion of the Yin Yang1 (Yy1) gene in mice causes a phenotype mimicking the cystic form of PPB and affects the expression of key regulators of lung development. Similar changes in expression were observed in PPB but not in CPAM lung biopsies, revealing a distinctive PPB molecular signature. Deregulation of molecules promoting epithelial-mesenchymal transition (EMT) was detected in PPB specimens, suggesting that EMT might participate in tumor progression. Changes in miRNA expression also occurred in PPB lung biopsies. miR-125a-3p, a candidate to regulate YY1 expression and lung branching, was abnormally highly expressed in PPB samples. Together, these findings support the concept that reduced expression of YY1, due to the abnormal miRNA profile resulting from DICER1 mutations, contributes to PPB development via its impact on the expression of key lung developmental genes.This article has an associated First Person interview with the joint first authors of the paper.

Role of PKCε in the epithelial-mesenchymal transition induced by FGFR2 isoform switch

Background

The epithelial isoform of the fibroblast growth factor receptor 2 (FGFR2b) controls the entire program of keratinocyte differentiation via the sequential involvement of protein kinase C (PKC) δ and PKCα. In contrast, the FGFR2 isoform switch and the aberrant expression of the mesenchymal FGFR2c isoform leads to impairment of differentiation, epithelial-mesenchymal transition (EMT) and tumorigenic features. Aim of our present study was to contribute in clarifying the complex network of signaling pathways involved in the FGFR2c-mediated oncogenic outcomes focusing on PKCε, which appears to be involved in the induction of EMT and tumorigenesis in several epithelial contexts.

Methods

Biochemical and molecular analysis, as well as in vitro invasion assays, combined with the use of specific small interfering RNA (siRNA), were performed in human keratinocytes stably expressing FGFR2c or FGFR2b isoforms.

Results

Our results showed that aberrant expression and signaling of FGFR2c, but not those of FGFR2b, in human keratinocytes induced a strong phosphorylation/activation of PKCε. The use of siRNA approach showed that PKCε is the hub signaling downstream FGFR2c responsible for the modulation of EMT markers and for the induction of the EMT-related transcription factors STAT3, Snail1 and FRA1, as well as for the acquisition of the invasive behavior. Moreover, experiments of depletion of ESRP1, responsible for FGFR2 splicing in epithelial cells, indicated that the activation of PKCε is the key molecular event triggered by FGFR2 isoform switch and underlying EMT induction.

Conclusions

Overall, our results point to the identification of the downstream PKC isoform responsible for the FGFR signaling deregulation occurring in epithelial tissues from the physiological oncosoppressive to the pathological oncogenic profile.

FGFR2c Mesenchymal Isoform Expression Is Associated with Poor Prognosis and Further Refines Risk Stratification within Endometrial Cancer Molecular Subtypes

PURPOSE

The two most common molecular subtypes of endometrial cancers, mismatch repair deficient (MMRd) and p53 wild-type (p53wt) comprise the majority of endometrial cancers and have intermediate prognoses where additional risk stratification biomarkers are needed. Isoform switching of FGFR2 from FGFR2b to FGFR2c (normally expressed in mesenchymal cells), has been reported in other solid carcinomas. The objective of this study was to investigate the role of FGFR2c in risk stratification of endometrial cancer.

EXPERIMENTAL DESIGN

We have developed and optimized a BaseScope RNA ISH assay to detect FGFR2c. FGFR2c expression was determined in a preliminary screening cohort of 78 endometrial cancers and a clinically and molecularly annotated Vancouver cohort (n = 465). Cox regression model analyses were performed to assess the prognostic value of FGFR2c.

RESULTS

Univariate and multivariate analyses revealed FGFR2c expression was significantly associated with shorter disease-specific survival (DSS) and progression-free survival (PFS) in endometrioid endometrial cancer (EEC, n = 302). Notably, FGFR2c expression was significantly associated with shorter PFS and DSS in patients with grade 3 EECs (P < 0.003 and P < 0.002) and the European Society Medical Oncology (ESMO) high-risk group (P < 0.0001 and P < 0.002), respectively. Moreover, within the MMRd subtype, FGFR2c expression was significantly associated with shorter PFS (P < 0.048) and DSS (P < 0.001).

CONCLUSIONS

FGFR2c expression appears an independent prognostic biomarker in patients with EEC and further discerns the outcomes within grade 3 tumors, ESMO high-risk groups, as well as within the MMRd and p53wt subtypes. FGFR2c inclusion into future molecular subtyping can further refine risk stratification of EEC.

Transcriptome complexity in intravascular NK/T-cell lymphoma

AIMS

Intravascular NK/T-cell lymphoma (IVNKTCL) is a rare disease, which is characterised by exclusive growth of large cells within the lumen of small vessels, Epstein-Barr virus infection and somatic mutations in epigenetic regulator genes. Here, we elucidate the transcriptomic complexity of IVNKTCL.

METHODS

IVNKTCL cases were retrieved from a single-centre cohort of 25 intravascular lymphomas. RNA-seq and whole exome sequencing (WES) were performed to analyse transcriptomic abnormalities and mutations in splicing factors.

RESULTS

Approximately 88% of the total reads from the RNA-seq were considered exonic, while the remaining reads (12%) were mapped to intronic or intergenic regions. We detected 28,941 alternative splicing events, some of which would produce abnormal proteins rarely found in normal cells. The detected events also included tumour-specific splicing alterations in oncogenes and tumour suppressors (e.g., HRAS, MDM2 and VEGFA). WES identified premature termination mutations or copy number losses in a total of 15 splicing regulator genes, including SF3B5, SRSF12 and TNPO3.

CONCLUSIONS

This study raises the possibility that IVNKTCL may be driven by multiple complex regulatory loops, including non-exonic expression and aberrant splicing, in addition to defects in epigenetic regulation.

The RNA-binding protein AKAP8 suppresses tumor metastasis by antagonizing EMT-associated alternative splicing

Alternative splicing has been shown to causally contribute to the epithelial–mesenchymal transition (EMT) and tumor metastasis. However, the scope of splicing factors that govern alternative splicing in these processes remains largely unexplored. Here we report the identification of A-Kinase Anchor Protein (AKAP8) as a splicing regulatory factor that impedes EMT and breast cancer metastasis. AKAP8 not only is capable of inhibiting splicing activity of the EMT-promoting splicing regulator hnRNPM through protein–protein interaction, it also directly binds to RNA and alters splicing outcomes. Genome-wide analysis shows that AKAP8 promotes an epithelial cell state splicing program. Experimental manipulation of an AKAP8 splicing target CLSTN1 revealed that splice isoform switching of CLSTN1 is crucial for EMT. Moreover, AKAP8 expression and the alternative splicing of CLSTN1 predict breast cancer patient survival. Together, our work demonstrates the essentiality of RNA metabolism that impinges on metastatic breast cancer.

Splice isoform switching regulated by the heterogeneous nuclear ribonucleoprotein M (hnRNPM) induces EMT and metastasis. Here, the authors report that AKAP8 is a metastasis suppressor that inhibits the splicing activity of hnRNPM and antagonizes genome-wide EMT-associated alternative splicing to maintain epithelial cell state.

Genome-wide Analysis of the Alternative Splicing Profiles Revealed Novel Prognostic Index for Kidney Renal Cell Clear Cell Carcinoma

Alternative splicing (AS) is a major mechanism that greatly enhanced the diversity of proteome. Mounting evidence demonstrated that aberration of AS are important steps for the initiation and progression of human cancers. Here, we comprehensively investigated the association between whole landscape of AS profiles and the survival outcome of renal cell carcinoma (RCC) patients using RNA-seq data from TCGA SpliceSeq. Because of the limited number size of deaths in kidney chromophobe renal cell carcinoma (KICH) and papillary renal cell carcinoma (KIRP) TCGA cohorts, we only conducted survival analysis in kidney clear renal cell carcinoma (KIRC). We further constructed prognostic index (PI) based on prognosis-related AS events and built correlation network for splicing factors and prognosis-related AS events. According to the results, a total of 5351 AS events in 3522 genes were significantly correlated with the overall survival (OS) of kidney clear cell renal cell carcinoma (KIRC) patients. Seven of the PI models exhibited preferable prognosis-predicting capacity for KIRC with PI-ALL reaching the highest area under curve value of 0.875. The splicing regulatory network between splicing factors and prognosis-related AS events depicted a tangled web of relationships between them. One of the splicing factors: KHDRBS3 was validated by immunohistochemistry to be down-regulated in KIRC tissues. In conclusion, the powerful efficiency of risk stratification of PI models indicated the potential of AS signature as promising prognostic markers for KIRC and the splicing regulation network provided possible genetic mechanism of KIRC.

PTPN21-CDSlong isoform inhibits the response of acute lymphoblastic leukemia cells to NK-mediated lysis via the KIR/HLA-I axis

Protein tyrosine phosphatase non-receptor type 21 (PTPN21) is a member of the non-receptor tyrosine phosphatase family. We have found that PTPN21 is mutated in relapsed Philadelphia chromosome-negative acute lymphoblastic leukemia (ALL) after allogeneic hematopoietic stem cell transplantation. PTPN21 consists of three types of isoforms according to the length of the protein encoded. However, the roles of different isoforms in leukemic cells have not been elucidated. In the study, PTPN21 isoform constitution in five ALL cell lines were identified by transcriptome polymerase chain reaction combined with Sanger sequencing, and the relationship between PTPN21 isoforms and sensitivity to natural killer (NK) cells mediated killing in ALL cell lines were further assessed by knock-out of different isoforms of PTPN21 using CRISPR-Cas9 technique. Subsequently, we explored the functional mechanisms through RNA sequencing and confirmatory testing. The results showed that there was no significant change when all PTPN21 isoforms were knocked out in ALL cells, but the sensitivity of NALM6 cells with PTPN21-CDS_long knock-out (NALM6-PTPN21_lk ) to NK-mediated killing was significantly increased. Whole transcriptome sequencing and further validation testing showed that human leukocyte antigen class I (HLA-I) molecules were significantly decreased, accompanied by a significantly downregulated expression of antigen presenting-related chaperones in NALM6-PTPN21_lk cells. Our results uncovered a previously unknown mechanism that PTPN21-CDS_long and CDS_short isoforms may play opposite roles in NK-mediated killing in ALL cells, and showed that the endogenous PTPN21-CDS_long isoform inhibited ALL cells to NK cell-mediated lysis by regulating the KIR-HLA-I axis.

Integrated Analysis of Structural Variation and RNA Expression of FGFR2 and Its Splicing Modulator ESRP1 Highlight the ESRP1amp-FGFR2norm-FGFR2-IIIchigh Axis in Diffuse Gastric Cancer

Gastric Cancer (GC) is one of the most common and deadliest types of cancer in the world. To improve GC prognosis, increasing efforts are being made to develop new targeted therapies. Although FGFR2 genetic amplification and protein overexpression in GC have been targeted in clinical trials, so far no improvement in patient overall survival has been found. To address this issue, we studied genetic and epigenetic events affecting FGFR2 and its splicing regulator ESRP1 in GC that could be used as new therapeutic targets or predictive biomarkers. We performed copy number variation (CNV), DNA methylation, and RNA expression analyses of FGFR2/ESRP1 across several cohorts. We discovered that both genes were frequently amplified and demethylated in GC, resulting in increased ESRP1 expression and of a specific FGFR2 isoform: FGFR2-IIIb. We also showed that ESRP1 amplification in GC correlated with a significant decreased expression of FGFR2-IIIc, an alternative FGFR2 splicing isoform. Furthermore, when we performed a survival analysis, we observed that patients harboring diffuse-type tumors with low FGFR2-IIIc expression revealed a better overall survival than patients with FGFR2-IIIc high-expressing diffuse tumors. Our results encourage further studies on the role of ESRP1 in GC and support FGFR2-IIIc as a relevant biomarker in GC.

Emerging roles of histone modifications and HDACs in RNA splicing

Histone modifications and RNA splicing, two seemingly unrelated gene regulatory processes, greatly increase proteome diversity and profoundly influence normal as well as pathological eukaryotic cellular functions. Like many histone modifying enzymes, histone deacetylases (HDACs) play critical roles in governing cellular behaviors and are indispensable in numerous biological processes. While the association between RNA splicing and histone modifications is beginning to be recognized, a lack of knowledge exists regarding the role of HDACs in splicing. Recent studies however, reveal that HDACs interact with spliceosomal and ribonucleoprotein complexes, actively control the acetylation states of splicing-associated histone marks and splicing factors, and thereby unexpectedly could modulate splicing. Here, we review the role of histone/protein modifications and HDACs in RNA splicing and discuss the convergence of two parallel fields, which supports the argument that HDACs, and perhaps most histone modifying enzymes, are much more versatile and far more complicated than their initially proposed functions. Analogously, an HDAC-RNA splicing connection suggests that splicing is regulated by additional upstream factors and pathways yet to be defined or not fully characterized. Some human diseases share common underlying causes of aberrant HDACs and dysregulated RNA splicing and, thus, further support the potential link between HDACs and RNA splicing.

Cancer the'RBP'eutics-RNA-binding proteins as therapeutic targets for cancer

RNA-binding proteins (RBPs) play a critical role in the regulation of various RNA processes, including splicing, cleavage and polyadenylation, transport, translation and degradation of coding RNAs, non-coding RNAs and microRNAs. Recent studies indicate that RBPs not only play an instrumental role in normal cellular processes but have also emerged as major players in the development and spread of cancer. Herein, we review the current knowledge about RNA binding proteins and their role in tumorigenesis as well as the potential to target RBPs for cancer therapeutics.

The Aberrant Expression of the Mesenchymal Variant of FGFR2 in the Epithelial Context Inhibits Autophagy

Signaling of the epithelial splice variant of fibroblast growth factor receptor 2 (FGFR2b) triggers both differentiation and autophagy, while the aberrant expression of the mesenchymal FGFR2c isoform in epithelial cells induces impaired differentiation, epithelial mesenchymal transition (EMT) and tumorigenic features. Here we analyzed in the human keratinocyte cell line, as well as in primary cultured cells, the possible impact of FGFR2c forced expression on the autophagic process. Biochemical and quantitative immunofluorescence analysis, coupled to the use of autophagic flux sensors, specific substrate inhibitors or silencing approaches, showed that ectopic expression and the activation of FGFR2c inhibit the autophagosome formation and that AKT/MTOR is the downstream signaling mainly involved. Interestingly, the selective inhibition of AKT or MTOR substrates caused a reversion of the effects of FGFR2c on autophagy, which could also arise from the imbalance of the interplay between AKT/MTOR pathway and JNK1 signaling in favor of JNK1 activation, BCL-2 phosphorylation and possibly phagophore nucleation. Finally, silencing experiments of depletion of ESRP1, responsible for FGFR2 splicing and consequent FGFR2b expression, indicated that the switching from FGFR2b to FGFR2c isoform could represent the key event underlying the inhibition of the autophagic process in the epithelial context. Our results provide the first evidence of a negative impact of the out-of-context expression of FGFR2c on autophagy, suggesting a possible role of this receptor in the modulation of the recently proposed negative loop between autophagy and EMT during carcinogenesis.

Cause-and-Effect relationship between FGFR1 expression and epithelial-mesenchymal transition in EGFR-mutated non-small cell lung cancer cells

OBJECTIVES

Increased FGFR1 expression is associated with resistance to tyrosine kinase inhibitors (TKIs) in EGFR-mutated NSCLC cells and often concomitant with epithelial to mesenchymal transition (EMT). However, the cause-and-effect relationship between increased FGFR1 expression and EMT in the genetic background of EGFR-mutated non-small cell lung cancer (NSCLC) cells is not clear. Previous studies have specifically addressed the relationship between EMT and increased FGFR1 expression in the context of simultaneous TKI-mediated blocking of EGFR-signaling. Here, in the context of EGFR-mutated NSCLC cells with active EGFR-signaling, we have examined whether increased FGFR1 expression drives EMT or is an EMT passenger event.

MATERIALS AND METHODS

For cause-and-effect analyses between EMT and FGFR1 expression, including expression of alternative spliced FGFR1 isoforms, we used CRISPR-dCAS9-SAM-mediated induction of the endogenous FGFR1 and ZEB1 genes, as well as biochemical EMT-induction, in PC9 and HCC827 NSCLC cell lines harboring activating EGFR-mutations.

RESULTS

We find that FGFR1 expression correlates with a ZEB1-associated EMT gene expression profile in NSCLC cells. In experiments using NSCLC cell lines harboring activating EGFR-mutations we show that CRISPR-dCAS9-SAM-mediated induction of FGFR1 expression is neither driving an increase in ZEB1 expression nor EMT characteristics. However, CRISPR-dCAS9-SAM-mediated induction of ZEB1 expression drives EMT characteristics and an increase in FGFR1 expression. Biochemical induction of EMT also drives an increase in FGFR1 expression.

CONCLUSION

From our findings concerning the cause-and-effect relationship in the genetic background of EGFR-mutated NSCLC cells, we conclude that an increase in ZEB1 expression is a driver of EMT resulting in concomitant increased FGFR1 expression, whereas an increase in FGFR1 expression is insufficient to drive concomitant EMT.

Overexpression of grainyhead-like transcription factor 2 is associated with poor prognosis in human pancreatic carcinoma

Recent studies have demonstrated that the abnormal expression of the grainyhead-like transcription factor 2 (GRHL2) gene contributes to the progression and poor prognosis of cancer through multiple mechanisms, but little is known about its expression status and prognostic value in pancreatic carcinoma (PC). The aim of the present study was to investigate the expression of GRHL2 in PC and to evaluate its clinicopathological and prognostic significance. Immunohistochemistry and western blotting were used to detect the expression of GRHL2 in PC tissues and cell lines, respectively. The expression of GRHL2 was investigated in 92 PC tissue samples by immunohistochemistry. High expression of GRHL2 was significantly associated with histological differentiation (P=0.018) and lymphatic metastasis (P=0.024). A Kaplan-Meier analysis revealed that high expression of GRHL2 was associated with worsened overall survival time (P<0.001). Multivariate analysis indicated that GRHL2 may be an independent prognostic factor for poor overall survival time (P=0.001). Additionally, western blot analysis demonstrated that the GRHL2 protein was highly expressed in PC cell lines. GRHL2 may serve an important role in the tumourigenesis of PC and serve as a potential therapeutic target for the prevention of PC progression.

The Cancer Spliceome: Reprograming of Alternative Splicing in Cancer

Alternative splicing allows for the expression of multiple RNA and protein isoforms from one gene, making it a major contributor to transcriptome and proteome diversification in eukaryotes. Advances in next generation sequencing technologies and genome-wide analyses have recently underscored the fact that the vast majority of multi-exon genes under normal physiology engage in alternative splicing in tissue-specific and developmental-specific manner. On the other hand, cancer cells exhibit remarkable transcriptome alterations partly by adopting cancer-specific splicing isoforms. These isoforms and their encoded proteins are not insignificant byproducts of the abnormal physiology of cancer cells, but either drivers of cancer progression or small but significant contributors to specific cancer hallmarks. Thus, it is paramount that the pathways that regulate alternative splicing in cancer, including the splicing factors that bind to pre-mRNAs and modulate spliceosome recruitment. In this review, we present a few distinct cases of alternative splicing in cancer, with an emphasis on their regulation as well as their contribution to cancer cell phenotype. Several categories of splicing aberrations are highlighted, including alterations in cancer-related genes that directly affect their pre-mRNA splicing, mutations in genes encoding splicing factors or core spliceosomal subunits, and the seemingly mutation-free disruptions in the balance of the expression of RNA-binding proteins, including components of both the major (U2-dependent) and minor (U12-dependent) spliceosomes. Given that the latter two classes cause global alterations in splicing that affect a wide range of genes, it remains a challenge to identify the ones that contribute to cancer progression. These challenges necessitate a systematic approach to decipher these aberrations and their impact on cancer. Ultimately, a sufficient understanding of splicing deregulation in cancer is predicted to pave the way for novel and innovative RNA-based therapies.

Transcriptomic signature associated with carcinogenesis and aggressiveness of papillary thyroid carcinoma

Papillary thyroid carcinoma (PTC) is the fastest-growing disease caused by numerous molecular alterations in addition to previously reported DNA mutations. There is a compelling need to identify novel transcriptomic alterations that are associated with the pathogenesis of PTC with potential diagnostic and prognostic implications. Methods: We gathered and compared 242 expression profiles between paired PTC and adjacent normal tissues and identified and validated the coding and long non-coding RNAs (lncRNAs) associated with the extrathyroidal extension (ETE) of 655 PTC patients in two independent cohorts, followed by predicting their interactions with drugs. Co-expression, RNA interaction, Kaplan-Meier survival and multivariate Cox proportional regression analyses were performed to identify dysregulated lncRNAs and genes that correlated with clinical outcomes of PTC. Alternative splicing (AS), RNA circularization, and editing were also compared between transcriptomes to expand the repertoire of molecular alterations in PTC. Results: Numerous genes related to cellular microenvironment and steroid hormone response were associated with the ETE of PTC. Drug susceptibility predictions of the expression signature revealed two highly ranked compounds, 6-bromoindirubin-3'-oxime and lovastatin. Co-expression and RNA interaction analysis revealed the essential role of lncRNAs in PTC pathogenesis by modulating extracellular matrix and cell adhesion. Eight genes and two novel lncRNAs were identified that correlated with the aggressive nature and disease-free survival of PTC. Furthermore, this study provided the transcriptome-wide landscape of circRNAs in PTC and uncovered dissimilar expression profiles among circRNAs originating from the same host gene, suggesting the functional complexity of circRNAs in PTC carcinogenesis. The newly identified AS events in the SERPINA1 and FN1 genes may improve the sensitivity and specificity of these diagnostic biomarkers. Conclusions: Our study uncovered a comprehensive transcriptomic signature associated with the carcinogenesis and aggressive behavior of PTC, as well as presents a catalog of 10 potential biomarkers, which would facilitate PTC prognosis and development of new therapeutic strategies for this cancer.

Alternative-splicing defects in cancer: Splicing regulators and their downstream targets, guiding the way to novel cancer therapeutics

Defects in alternative splicing are frequently found in human tumors and result either from mutations in splicing-regulatory elements of specific cancer genes or from changes in the regulatory splicing machinery. RNA splicing regulators have emerged as a new class of oncoproteins and tumor suppressors, and contribute to disease progression by modulating RNA isoforms involved in the hallmark cancer pathways. Thus, dysregulation of alternative RNA splicing is fundamental to cancer and provides a potentially rich source of novel therapeutic targets. Here, we review the alterations in splicing regulatory factors detected in human tumors, as well as the resulting alternatively spliced isoforms that impact cancer hallmarks, and discuss how they contribute to disease pathogenesis. RNA splicing is a highly regulated process and, as such, the regulators are themselves tightly regulated. Differential transcriptional and posttranscriptional regulation of splicing factors modulates their levels and activities in tumor cells. Furthermore, the composition of the tumor microenvironment can also influence which isoforms are expressed in a given cell type and impact drug responses. Finally, we summarize current efforts in targeting alternative splicing, including global splicing inhibition using small molecules blocking the spliceosome or splicing-factor-modifying enzymes, as well as splice-switching RNA-based therapeutics to modulate cancer-specific splicing isoforms. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Processing > Splicing Regulation/Alternative Splicing.

Global isoform-specific transcript alterations and deregulated networks in clear cell renal cell carcinoma

Extensive genome-wide analyses of deregulated gene expression have now been performed for many types of cancer. However, most studies have focused on deregulation at the gene-level, which may overlook the alterations of specific transcripts for a given gene. Clear cell renal cell carcinoma (ccRCC) is one of the best-characterized and most pervasive renal cancers, and ccRCCs are well-documented to have aberrant RNA processing. In the present study, we examine the extent of aberrant isoform-specific RNA expression by reporting a comprehensive transcript-level analysis, using the new kallisto-sleuth-RATs pipeline, investigating coding and non-coding differential transcript expression in ccRCC. We analyzed 50 ccRCC tumors and their matched normal samples from The Cancer Genome Altas datasets. We identified 7,339 differentially expressed transcripts and 94 genes exhibiting differential transcript isoform usage in ccRCC. Additionally, transcript-level coexpression network analyses identified vasculature development and the tricarboxylic acid cycle as the most significantly deregulated networks correlating with ccRCC progression. These analyses uncovered several uncharacterized transcripts, including lncRNAs FGD5-AS1 and AL035661.1, as potential regulators of the tricarboxylic acid cycle associated with ccRCC progression. As ccRCC still presents treatment challenges, our results provide a new resource of potential therapeutics targets and highlight the importance of exploring alternative methodologies in transcriptome-wide studies.

Role of FGFR2b expression and signaling in keratinocyte differentiation: sequential involvement of PKCδ and PKCα

The tumor suppressor epithelial isoform of the fibroblast growth factor receptor 2 (FGFR2b) induces human keratinocyte early differentiation. Moreover, protein kinases C (PKCs) are known to regulate the differentiation program in several cellular contexts, including keratinocytes. Therefore, in this paper we propose to clarify if FGFR2b could play a role also in the late steps of keratinocyte differentiation and to assess if this receptor-induced process would sequentially involve PKCδ and PKCα isoforms. Immunofluorescence, biochemical, and molecular approaches, performed on 2D cultures or 3D organotypic rafts of human keratinocytes overexpressing FGFR2b by stable transduction, showed that receptor signaling induced the precocious onset and an accelerated progression of keratinocyte differentiation, indicating that FGFR2b is a crucial regulator of the entire program of keratinocyte differentiation. In addition, the use of specific inhibitors and gene silencing approaches through specific siRNA demonstrated that PKCδ controls the onset of FGFR2b-triggered differentiation, while PKCα plays a role restricted to the terminal stages of the process. Molecular analysis revealed that the two PKC isoforms sequentially act via induction of KLF4 and DLX3, two transcription factors linked by negative loops to p63, suggesting that p63 would represent the hub molecule at the crossroad of an intricate signaling network downstream FGFR2b, involving multiple PKC-induced transcription factors.

Epithelial splicing regulatory protein 1 and 2 paralogues correlate with splice signatures and favorable outcome in human colorectal cancer

ESRPs are master splice regulators implicated in alternative mRNA splicing programs important for epithelial-mesenchymal transition (EMT) and tumor progression. ESRP1 was identified in some tumors as good or worse predictor of outcome, but in colorectal cancer (CRC) the prognostic value of ESRPs and relation with mesenchymal splice variants is not clear. Here, we studied 68 CRC cases, compared tissue expression of ESRPs with clinical data and with EMT gene splice patterns of conditional CRC cells with deficient ESRP1 expression.Around 72% of patients showed global decreased transcript expression of both ESRPs in tumor as compared to matched non-neoplastic colorectal epithelium. Reduction of ESRP1 in tumor cells was evaluated by immunohistochemistry, associated with microsatellite stability and switch to mesenchymal splice signatures of FGFRs, CD44, ENAH and CTNND1(p120-catenin). Expression of ESRPs was significantly associated with favorable overall survival (log-rank test, P=0.0186 and 0.0408), better than prognostic stratification by tumor staging; and for ESRP1 confirmed with second TCGA cohort (log-rank test, P=0.0435). Prognostic value is independent of the pathological stage and microsatellite instability (ESRP1: HR=0.36, 95%CI 0.15-0.91, P=0.032; ESRP2: HR=0.23, 95%CI 0.08-0.65, P=0.006).Our study supports the role of ESRP1 as tumor suppressor and strongly suggests that ESRPs are candidate markers for early detection, diagnosis, and prognosis of CRC.

Human PDE4D isoform composition is deregulated in primary prostate cancer and indicative for disease progression and development of distant metastases

Phosphodiesterase 4D7 was recently shown to be specifically over-expressed in localized prostate cancer, raising the question as to which regulatory mechanisms are involved and whether other isoforms of this gene family (PDE4D) are affected under the same conditions.We investigated PDE4D isoform composition in prostatic tissues using a total of seven independent expression datasets and also included data on DNA methylation, copy number and AR and ERG binding in PDE4D promoters to gain insight into their effect on PDE4D transcription.We show that expression of PDE4D isoforms is consistently altered in primary human prostate cancer compared to benign tissue, with PDE4D7 being up-regulated while PDE4D5 and PDE4D9 are down-regulated. Disease progression is marked by an overall down-regulation of long PDE4D isoforms, while short isoforms (PDE4D1/2) appear to be relatively unaffected. While these alterations seem to be independent of copy number alterations in the PDE4D locus and driven by AR and ERG binding, we also observed increased DNA methylation in the promoter region of PDE4D5, indicating a long lasting alteration of the isoform composition in prostate cancer tissues.We propose two independent metrics that may serve as diagnostic and prognostic markers for prostate disease: (PDE4D7 - PDE4D5) provides an effective means for distinguishing PCa from normal adjacent prostate, whereas PDE4D1/2 - (PDE4D5 + PDE4D7 + PDE4D9) offers strong prognostic potential to detect aggressive forms of PCa and is associated with metastasis free survival. Overall, our findings highlight the relevance of PDE4D as prostate cancer biomarker and potential drug target.

Integrative analysis of protein-coding and non-coding RNAs identifies clinically relevant subtypes of clear cell renal cell carcinoma

Protein-coding genes and non-coding RNAs cooperate mutually in cells. Integrative analysis of protein-coding and non-coding RNAs may facilitate characterizing tumor heterogeneity. We introduced integrated consensus clustering (ICC) method to integrate mRNA, miRNA and lncRNA expression profiles of 431 primary clear cell renal cell carcinomas (ccRCCs). We identified one RCC subgroup easily misdiagnosed as ccRCC in clinic and four robust ccRCC subtypes associated with distinct clinicopathologic and molecular features. In subtype R1, AMPK signaling pathway is significantly upregulated, which may improve the oncologic-metabolic shift and partially account for its best prognosis. Subtype R2 has more chromosomal abnormities, higher expression of cell cycle genes and less expression of genes in various metabolism pathways, which may explain its more aggressive characteristic and the worst prognosis. Moreover, much more miRNAs and lncRNAs are significantly upregulated in R2 and R4 respectively, suggesting more important roles of miRNAs in R2 and lncRNAs in R4. Triple-color co-expression network analysis identified 28 differentially expressed modules, indicating the importance of cooperative regulation of mRNAs, miRNAs and lncRNAs in ccRCC. This study establishes an integrated transcriptomic classification which may contribute to understanding the heterogeneity and implicating the treatment of ccRCC.

Expression of the FGFR2c mesenchymal splicing variant in human keratinocytes inhibits differentiation and promotes invasion

The altered isoform switching of the fibroblast growth factor receptor 2 (FGFR2) and aberrant expression of the mesenchymal FGFR2c isoform in epithelial cells is involved in cancer progression. We have recently described that the ectopic expression of FGFR2c in normal human keratinocytes induces epithelial‐mesenchymal transition and leads to invasiveness and anchorage‐independent growth. Here, we extended our analysis to the effects of this FGFR2c forced expression on human keratinocyte differentiation and stratification. Our findings demonstrated that, differently from cells overexpressing the epithelial splicing variant FGFR2b, keratinocytes ectopically expressing FGFR2c are not able to form a monolayer and display decreased expression of early differentiation markers. This impaired ability to enter the differentiation program is related to the up‐modulation of the transcription factor ΔNp63. In addition, FGFR2c‐expressing keratinocytes undergo defective stratification and invasion of the collagen matrix in 3D organotypic cultures, further suggesting their tumorigenic potential. Taken together, our results support the hypothesis that the receptor switching and the consequent appearance of the mesenchymal FGFR2c variant in the epithelial context would drive early steps of carcinogenesis, unbalancing the p63/FGFR interplay, and altering the paracrine response to the microenvironment.

Regulation of splicing and circularisation of RNA in epithelial mesenchymal plasticity

Interconversions between epithelial and mesenchymal states, often referred to as epithelial mesenchymal transition (EMT) and its reverse MET, play important roles in embryonic development and are recapitulated in various adult pathologies including cancer progression. These conversions are regulated by complex transcriptional and post-transcriptional mechanisms including programs of alternative splicing which are orchestrated by specific splicing factors. This review will focus on the latest developments in our understanding of the splicing factors regulating epithelial mesenchymal plasticity associated with cancer progression and the induction of pluripotency, including potential roles for circular RNAs (circRNAs) which have been recently implicated in these processes.

RNA-Binding Proteins in Cancer: Old Players and New Actors

RNA-binding proteins (RBPs) are key players in post-transcriptional events. The combination of versatility of their RNA-binding domains with structural flexibility enables RBPs to control the metabolism of a large array of transcripts. Perturbations in RBP-RNA networks activity have been causally associated with cancer development, but the rational framework describing these contributions remains fragmented. We review here the evidence that RBPs modulate multiple cancer traits, emphasize their functional diversity, and assess future trends in the study of RBPs in cancer.

Role of CD56 in Normal Kidney Development and Wilms Tumorigenesis

The cell-surface glycoprotein CD56 has three major isoforms that play important roles in cell adhesion and signaling, which may promote cell proliferation, differentiation, survival, or migration. It is an important molecule in normal kidney development and acts as a key marker in Wilms tumor stem and progenitor cells. Here, we review the structural and genetic features of the CD56 glycoprotein, and summarize its roles in the normal versus diseased metanephric blastema. We discuss areas of CD56-related research that may complement or improve existing Wilms tumor treatment strategies, including the antibody-drug conjugate lorvotuzumab mertansine that binds to CD56.

EMT and stemness: flexible processes tuned by alternative splicing in development and cancer progression

Epithelial-to-mesenchymal transition (EMT) is associated with metastasis formation as well as with generation and maintenance of cancer stem cells. In this way, EMT contributes to tumor invasion, heterogeneity and chemoresistance. Morphological and functional changes involved in these processes require robust reprogramming of gene expression, which is only partially accomplished at the transcriptional level. Alternative splicing is another essential layer of gene expression regulation that expands the cell proteome. This step in post-transcriptional regulation of gene expression tightly controls cell identity between epithelial and mesenchymal states and during stem cell differentiation. Importantly, dysregulation of splicing factor function and cancer-specific splicing isoform expression frequently occurs in human tumors, suggesting the importance of alternative splicing regulation for cancer biology.

In this review, we briefly discuss the role of EMT programs in development, stem cell differentiation and cancer progression. Next, we focus on selected examples of key factors involved in EMT and stem cell differentiation that are regulated post-transcriptionally through alternative splicing mechanisms. Lastly, we describe relevant oncogenic splice-variants that directly orchestrate cancer stem cell biology and tumor EMT, which may be envisioned as novel targets for therapeutic intervention.

Expression of the FGFR2 mesenchymal splicing variant in epithelial cells drives epithelial-mesenchymal transition

The FGFRs are receptor tyrosine kinases expressed by tissue-specific alternative splicing in epithelial IIIb or mesenchymal IIIc isoforms. Deregulation of FGF/FGFR signaling unbalances the epithelial-stromal homeostasis and may lead to cancer development. In the epithelial-context, while FGFR2b/KGFR acts as tumor suppressor, FGFR2c appears to play an oncogenic role. Based on our recent observation that the switching of FGFR2b versus FGFR2c induces EMT, here we investigated the biological outcome of the ectopic expression of FGFR2c in normal human keratinocytes. Morphological analysis showed that, differently from FGFR2b overexpression, the forced expression and activation of FGFR2c drive the epithelial cells to acquire a mesenchymal-like shape and actin reorganization. Moreover, the appearance of invasiveness and anchorage-independent growth ability in FGFR2c transfected keratinocytes was consistent with the potential tumorigenic role proposed for this receptor variant. Biochemical and molecular approaches revealed that the observed phenotypic changes were accompanied by modulation of EMT biomarkers and indicated the involvement of EMT transcription factors and miRs. Finally, the analysis of the expression pattern of discriminating markers strongly suggested that activation of FGFR2c triggers a process corresponding to the initiation of the pathological type III EMT, but not to the more physiological type II EMT occurring during FGFR2b-mediated wound healing.

ISOexpresso: a web-based platform for isoform-level expression analysis in human cancer

Background

Alternative splicing events that result in the production of multiple gene isoforms reveals important molecular mechanisms. Gene isoforms are often differentially expressed across organs and tissues, developmental stages, and disease conditions. Specifically, recent studies show that aberrant regulation of alternative splicing frequently occurs in cancer to affect tumor cell transformation and growth. While analysis of isoform expression is important for discovering tumor-specific isoform signatures and interpreting relevant genomic mutations, there is currently no web-based, easy-to-use, and publicly available platform for this purpose.

Description

We developed ISOexpresso to provide information regarding isoform existence and expression, which can be grouped by cancer vs. normal conditions, cancer types, and tissue types. ISOexpresso implements two main functions: First, the Isoform Expression View function creates visualizations for condition-specific RNA/isoform expression patterns upon query of a gene of interest. With this function, users can easily determine the major isoform (the most expressed isoform in a sample) of a gene with respect to the condition and check whether it matches the known canonical isoform. ISOexpresso outputs expression levels of all known transcripts to check alterations of expression landscape and to find potential tumor-specific isoforms. Second, the User Data Annotation function supports annotation of genomic variants to determine the most plausible consequence of a variation (e.g., an amino acid change) among many possible interpretations. As most coding sequence mutations are effective through the subsequent transcription and translation, ISOexpresso automatically prioritizes transcripts that act as backbones for mutation effect prediction by their relative expression. By employing ISOexpresso, we could investigate the consistency between the most expressed and known canonical/principal isoforms, as well as infer candidate tumor-specific isoforms based on their expression levels. In addition, we confirmed that ISOexpresso could easily reproduce previously known isoform expression patterns: recurrent observation of a major isoform across tissues, differential isoform expression patterns in a given tissue, and switching of major isoform during tumorigenesis.

Conclusions

ISOexpresso serves as a web-based, easy-to-use platform for isoform expression and alteration analysis based on large-scale cancer database. We anticipate that ISOexpresso will expedite formulation and confirmation of novel hypotheses by providing isoform-level perspectives on cancer research. The ISOexpresso database is available online at http://wiki.tgilab.org/ISOexpresso/.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2852-6) contains supplementary material, which is available to authorized users.

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.

Altered Expression and Splicing of ESRP1 in Malignant Melanoma Correlates with Epithelial-Mesenchymal Status and Tumor-Associated Immune Cytolytic Activity

Melanoma is one of the major cancer types for which new immune-based cancer treatments have achieved promising results. However, anti-PD-1 and anti-CTLA-4 therapies are effective only in some patients. Hence, predictive molecular markers for the development of clinical strategies targeting immune checkpoints are needed. Using The Cancer Genome Atlas (TCGA) RNAseq data, we found that expression of ESRP1, encoding a master splicing regulator in the epithelial-mesenchymal transition (EMT), was inversely correlated with tumor-associated immune cytolytic activity. That association holds up across multiple TCGA tumor types, suggesting a link between tumor EMT status and infiltrating lymphocyte activity. In melanoma, ESRP1 mainly exists in a melanocyte-specific truncated form transcribed from exon 13. This was validated by analyzing CCLE cell line data, public CAGE data, and RT-PCR in primary cultured melanoma cell lines. Based on ESRP1 expression, we divided TCGA melanoma cases into ESRP1-low, -truncated, and -full-length groups. ESRP1-truncated tumors comprise approximately two thirds of melanoma samples and reside in an apparent transitional state between epithelial and mesenchymal phenotypes. ESRP1 full-length tumors express epithelial markers and constitute about 5% of melanoma samples. In contrast, ESRP1-low tumors express mesenchymal markers and are high in immune cytolytic activity as well as PD-L2 and CTLA-4 expression. Those tumors are associated with better patient survival. Results from our study suggest a path toward the use of ESRP1 and other EMT markers as informative biomarkers for immunotherapy. Cancer Immunol Res; 4(6); 552-61. ©2016 AACR.

jSplice: a high-performance method for accurate prediction of alternative splicing events and its application to large-scale renal cancer transcriptome data

MOTIVATION

Alternative splicing represents a prime mechanism of post-transcriptional gene regulation whose misregulation is associated with a broad range of human diseases. Despite the vast availability of transcriptome data from different cell types and diseases, bioinformatics-based surveys of alternative splicing patterns remain a major challenge due to limited availability of analytical tools that combine high accuracy and rapidity.

RESULTS

We describe here a novel junction-centric method, jSplice, that enables de novo extraction of alternative splicing events from RNA-sequencing data with high accuracy, reliability and speed. Application to clear cell renal carcinoma (ccRCC) cell lines and 65 ccRCC patients revealed experimentally validatable alternative splicing changes and signatures able to prognosticate ccRCC outcome. In the aggregate, our results propose jSplice as a key analytic tool for the derivation of cell context-dependent alternative splicing patterns from large-scale RNA-sequencing datasets.

AVAILABILITY AND IMPLEMENTATION

jSplice is a standalone Python application freely available at http://www.mhs.biol.ethz.ch/research/krek/jsplice

CONTACT

wilhelm.krek@biol.ethz.ch

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.