Unbalanced alternative splicing and its significance in cancer

MYCN and SNRPD3 cooperate to maintain a balance of alternative splicing events that drives neuroblastoma progression

Many of the pro-tumorigenic functions of the oncogene MYCN are attributed to its regulation of global gene expression programs. Alternative splicing is another important regulator of gene expression and has been implicated in neuroblastoma development, however, the molecular mechanisms remain unknown. We found that MYCN up-regulated the expression of the core spliceosomal protein, SNRPD3, in models of neuroblastoma initiation and progression. High mRNA expression of SNRPD3 in human neuroblastoma tissues was a strong, independent prognostic factor for poor patient outcome. Repression of SNRPD3 expression correlated with loss of colony formation in vitro and reduced tumorigenicity in vivo. The effect of SNRPD3 on cell viability was in part dependent on MYCN as an oncogenic co-factor. RNA-sequencing revealed a global increase in the number of genes being differentially spliced when MYCN was overexpressed. Surprisingly, depletion of SNRPD3 in the presence of overexpressed MYCN further increased differential splicing, particularly of cell cycle regulators, such as BIRC5 and CDK10. MYCN directly bound SNRPD3, and the protein arginine methyltransferase, PRMT5, consequently increasing SNRPD3 methylation. Indeed, the PRMT5 inhibitor, JNJ-64619178, reduced cell viability and SNRPD3 methylation in neuroblastoma cells with high SNRPD3 and MYCN expression. Our findings demonstrate a functional relationship between MYCN and SNRPD3, which maintains the fidelity of MYCN-driven alternative splicing in the narrow range required for neuroblastoma cell growth. SNRPD3 methylation and its protein-protein interface with MYCN represent novel therapeutic targets. Hypothetical model for SNRPD3 as a co-factor for MYCN oncogenesis. SNRPD3 and MYCN participate in a regulatory loop to balance splicing fidelity in neuroblastoma cells. First MYCN transactivates SNRPD3 to lead to high-level expression. Second, SNRPD3 and MYCN form a protein complex involving PRMT5. Third, this leads to balanced alterative splicing (AS) activitiy that is favorable to neuroblastoma. Together this forms as a therapeutic vulnerability where SNRPD3 perturbation or PRMT5 inhibitors are selectively toxic to neuroblastoma by conditionally disturbing splicing activity.

Introns and Their Therapeutic Applications in Biomedical Researches

Context

Although for a long time, it was thought that intervening sequences (introns) were junk DNA without any function, their critical roles and the underlying molecular mechanisms in genome regulation have only recently come to light. Introns not only carry information for splicing, but they also play many supportive roles in gene regulation at different levels. They are supposed to function as useful tools in various biological processes, particularly in the diagnosis and treatment of diseases. Introns can contribute to numerous biological processes, including gene silencing, gene imprinting, transcription, mRNA metabolism, mRNA nuclear export, mRNA localization, mRNA surveillance, RNA editing, NMD, translation, protein stability, ribosome biogenesis, cell growth, embryonic development, apoptosis, molecular evolution, genome expansion, and proteome diversity through various mechanisms.

Evidence Acquisition

In order to fulfill the objectives of this study, the following databases were searched: Medline, Scopus, Web of Science, EBSCO, Open Access Journals, and Google Scholar. Only articles published in English were included.

Results & Conclusions

The intervening sequences of eukaryotic genes have critical functions in genome regulation, as well as in molecular evolution. Here, we summarize recent advances in our understanding of how introns influence genome regulation, as well as their effects on molecular evolution. Moreover, therapeutic strategies based on intron sequences are discussed. According to the obtained results, a thorough understanding of intron functional mechanisms could lead to new opportunities in disease diagnosis and therapies, as well as in biotechnology applications.

NRF2 has a splicing regulatory function involving the survival of motor neuron (SMN) in non-small cell lung cancer

The nuclear factor erythroid 2-like 2 (NFE2L2; NRF2) signaling pathway is frequently deregulated in human cancers. The critical functions of NRF2, other than its transcriptional activation, in cancers remain largely unknown. Here, we uncovered a previously unrecognized role of NRF2 in the regulation of RNA splicing. Global splicing analysis revealed that NRF2 knockdown in non-small cell lung cancer (NSCLC) A549 cells altered 839 alternative splicing (AS) events in 485 genes. Mechanistic studies demonstrated that NRF2 transcriptionally regulated SMN mRNA expression by binding to two antioxidant response elements in the SMN1 promoter. Post-transcriptionally, NRF2 was physically associated with the SMN protein. The Neh2 domain of NRF2, as well as the YG box and the region encoded by exon 7 of SMN, were required for their interaction. NRF2 formed a complex with SMN and Gemin2 in nuclear gems and Cajal bodies. Furthermore, the NRF2-SMN interaction regulated RNA splicing by expressing SMN in NRF2-knockout HeLa cells, reverting some of the altered RNA splicing. Moreover, SMN overexpression was significantly associated with alterations in the NRF2 pathway in patients with lung squamous cell carcinoma from The Cancer Genome Atlas. Taken together, our findings suggest a novel therapeutic strategy for cancers involving an aberrant NRF2 pathway.

METTL3-Mediated m6A Modification Controls Splicing Factor Abundance and Contributes to Aggressive CLL

RNA splicing dysregulation underlies the onset and progression of cancers. In chronic lymphocytic leukemia (CLL), spliceosome mutations leading to aberrant splicing occur in ∼20% of patients. However, the mechanism for splicing defects in spliceosome-unmutated CLL cases remains elusive. Through an integrative transcriptomic and proteomic analysis, we discover that proteins involved in RNA splicing are posttranscriptionally upregulated in CLL cells, resulting in splicing dysregulation. The abundance of splicing complexes is an independent risk factor for poor prognosis. Moreover, increased splicing factor expression is highly correlated with the abundance of METTL3, an RNA methyltransferase that deposits N6-methyladenosine (m6A) on mRNA. METTL3 is essential for cell growth in vitro and in vivo and controls splicing factor protein expression in a methyltransferase-dependent manner through m6A modification-mediated ribosome recycling and decoding. Our results uncover METTL3-mediated m6A modification as a novel regulatory axis in driving splicing dysregulation and contributing to aggressive CLL.

SIGNIFICANCE

METTL3 controls widespread splicing factor abundance via translational control of m6A-modified mRNA, contributes to RNA splicing dysregulation and disease progression in CLL, and serves as a potential therapeutic target in aggressive CLL. See related commentary by Janin and Esteller, p. 176. This article is highlighted in the In This Issue feature, p. 171.

Differential effects of WRAP53 transcript variants on non-small cell lung cancer cell behaviors

BACKGROUND

The WD40-encoding RNA antisense to p53 (WRAP53) is an antisense gene of TP53 with three transcriptional start sites producing three transcript variants involved in the progression of non-small cell lung cancer. However, the mechanism by which these different transcript variants regulate non-small cell lung cancer cell behaviors is to be elucidated.

METHODS

Two non-small cell lung cancer cell lines, A549 cells with wild-type p53 and H1975 with mutated p53, were transfected with WRAP53-1α and WRAP53-1β siRNA. The biological effects were assessed via colony formation, cell viability, apoptosis, cell cycle, wound healing and cell invasion assays, as well as immunoblotting.

RESULTS

Knockdown of WRAP53-1α increased the mRNA and protein levels of p53; suppressed colony formation and proliferation of A549 cells but promoted them in H1975 cells; increased the proportion of cells in the G0/G1 phase in A549 cells but decreased that in H1975 cells; and suppressed migration and invasion in A549 cells but not in H1975 cells. Conversely, knockdown of WRAP53-1β had no effect on p53 expression; promoted the growth of A549 cells but not of H1975 cells; decreased the proportion of cells in the G0/G1 phase in A549 cells but not in H1975 cells; and promoted migration and invasion in A549 cells but not in H1975 cells. Knockdown of both WRAP53-1α and WRAP53-1β promoted apoptosis in A549 cells but not in H1975 cells.

CONCLUSIONS

WRAP53 transcript variants exerted different functions in non-small cell lung cancer cells and regulated non-small cell lung cancer cell behaviors depending on the p53 expression.

Transcriptome analysis of clock disrupted cancer cells reveals differential alternative splicing of cancer hallmarks genes

Emerging evidence points towards a regulatory role of the circadian clock in alternative splicing (AS). Whether alterations in core-clock components may contribute to differential AS events is largely unknown. To address this, we carried out a computational analysis on recently generated time-series RNA-seq datasets from three core-clock knockout (KO) genes (ARNTL, NR1D1, PER2) and WT of a colorectal cancer (CRC) cell line, and time-series RNA-seq datasets for additional CRC and Hodgkin’s lymphoma (HL) cells, murine WT, Arntl KO, and Nr1d1/2 KO, and murine SCN WT tissue. The deletion of individual core-clock genes resulted in the loss of circadian expression in crucial spliceosome components such as SF3A1 (in ARNTL^KO), SNW1 (in NR1D1^KO), and HNRNPC (in PER2^KO), which led to a differential pattern of KO-specific AS events. All HCT116^KO cells showed a rhythmicity loss of a crucial spliceosome gene U2AF1, which was also not rhythmic in higher progression stage CRC and HL cancer cells. AS analysis revealed an increase in alternative first exon events specific to PER2 and NR1D1 KO in HCT116 cells, and a KO-specific change in expression and rhythmicity pattern of AS transcripts related to cancer hallmarks genes including FGFR2 in HCT116_ARNTL^KO, CD44 in HCT116_NR1D1^KO, and MET in HCT116_PER2^KO. KO-specific changes in rhythmic properties of known spliced variants of these genes (e.g. FGFR2 IIIb/FGFR2 IIIc) correlated with epithelial-mesenchymal-transition signalling. Altogether, our bioinformatic analysis highlights a role for the circadian clock in the regulation of AS, and reveals a potential impact of clock disruption in aberrant splicing in cancer hallmark genes.

How VEGF-A and its splice variants affect breast cancer development - clinical implications

BACKGROUND

Altered expression levels and structural variations in the vascular endothelial growth factor (VEGF) have been found to play important roles in cancer development and to be associated with the overall survival and therapy response of cancer patients. Particularly VEGF-A and its splice variants have been found to affect physiological and pathological angiogenic processes, including tumor angiogenesis, correlating with tumor progression, mostly caused by overexpression. This review focuses on the expression and impact of VEGF-A splice variants under physiologic conditions and in tumors and, in particular, the distribution and role of isoform VEGF165b in breast cancer.

CONCLUSIONS AND PERSPECTIVES

Many publications already highlighted the importance of VEGF-A and its splice variants in tumor therapy, especially in breast cancer, which are summarized in this review. Furthermore, we were able to demonstrate that cytoplasmatic VEGFA/165b expression is higher in invasive breast cancer tumor cells than in normal tissues or stroma. These examples show that the detection of VEGF splice variants can be performed also on the protein level in formalin fixed tissues. Although no quantitative conclusions can be drawn, these results may be the starting point for further studies at a quantitative level, which can be a major step towards the design of targeted antibody-based (breast) cancer therapies.

Nicotinamide N-Methyltransferase in Acquisition of Stem Cell Properties and Therapy Resistance in Cancer

The activity of nicotinamide N-methyltransferase (NNMT) is tightly linked to the maintenance of the nicotinamide adenine dinucleotide (NAD⁺) level. This enzyme catalyzes methylation of nicotinamide (NAM) into methyl nicotinamide (MNAM), which is either excreted or further metabolized to N1-methyl-2-pyridone-5-carboxamide (2-PY) and H₂O₂. Enzymatic activity of NNMT is important for the prevention of NAM-mediated inhibition of NAD⁺-consuming enzymes poly-adenosine -diphosphate (ADP), ribose polymerases (PARPs), and sirtuins (SIRTs). Inappropriately high expression and activity of NNMT, commonly present in various types of cancer, has the potential to disrupt NAD⁺ homeostasis and cellular methylation potential. Largely overlooked, in the context of cancer, is the inhibitory effect of 2-PY on PARP-1 activity, which abrogates NNMT's positive effect on cellular NAD⁺ flux by stalling liberation of NAM and reducing NAD⁺ synthesis in the salvage pathway. This review describes, and discusses, the mechanisms by which NNMT promotes NAD⁺ depletion and epigenetic reprogramming, leading to the development of metabolic plasticity, evasion of a major tumor suppressive process of cellular senescence, and acquisition of stem cell properties. All these phenomena are related to therapy resistance and worse clinical outcomes.

SNRPD1 confers diagnostic and therapeutic values on breast cancers through cell cycle regulation

Background

SNRPD1 is a spliceosome-associated protein and has previously been implicated with important roles in cancer development.

Methods

Through analyzing the differential expression patterns and clinical association of splicing associated genes among tumor and tumor adjacent samples across different tumors and among different breast cancer subtypes, we identify the tumor promotive role of SNRPD1 using multiple publicly available datasets. Through pathway, gene ontology enrichment analysis and network construction, we linked the onco-therapeutic role of SNRPD1 with cell cycle. Via a series of experimental studies including knockdown assay, qPCR, western blotting, cell cycle, drug response assay, we confirmed the higher expression of SNPRD1 at both gene and protein expression levels in triple negative breast cancer cells, as well as its roles in promoting cell cycle and chemotherapy response.

Results

Our study revealed that SNRPD1 over-expression was significantly associated with genes involved in cell cycle, cell mitosis and chromatin replication, and silencing SNRPD1 in breast cancer cells could lead to halted tumor cell growth and cell cycle arrest at the G₀/G₁ stage. We also found that triple negative breast cancer cells with reduced SNRPD1 expression lost certain sensitivity to doxorubicin whereas luminal cancer cells did not.

Conclusions

Our results suggested the prognostic value of SNRPD1 on breast cancer survival, its potential as the therapeutic target halting cell cycle progression for breast cancer control, and warranted special attention on the combined use of doxorubicin and drugs targeting SNRPD1.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-01932-w.

Systematic Identification of Survival-Associated Alternative Splicing Events in Kidney Renal Clear Cell Carcinoma

There is growing evidence that aberrant alternative splicing (AS) is highly correlated with driving tumorigenesis, but its function in kidney renal clear cell carcinoma (KIRC) remains to be discovered. In this study, we obtained the level-3 RNA sequencing and clinical data of KIRC from The Cancer Genome Atlas (TGCA). Combining with the splicing event detail information from TGCA SpliceSeq database, we established the independent prognosis signatures for KIRC with the univariate and multivariate Cox regression analyses. Then, we used the Kaplan-Meier analysis and receiver operating characteristic curves (ROCs) to assess the accuracy of prognosis signatures. We also constructed the regulatory network of splicing factors (SFs) and AS events. Our results showed that a total of 12029 survival-associated AS events of 5761 genes were found in 524 KIRC patients. All types of prognosis signatures displayed a satisfactory ability to reliably predict, especially in exon skip model which the area under curve of ROC was 0.802. Moreover, 18 splicing factors (SFs) highly correlated to AS events were identified. With the construction of the SF-AS interactive network, we found that SF powerfully promotes the occurrence of abnormal AS and may have a profound role in KIRC. Collectively, we screened survival-associated AS events and established prognosis signatures for KIRC, coupling with the SF-AS interactive network, which might provide a key perspective to clarify the potential mechanism of AS in KIRC.

microRNAs Promoting Growth of Gastric Cancer Xenografts and Correlation to Clinical Prognosis

The annual death toll for gastric cancer is in the range of 700,000 worldwide. Even in patients with early-stage gastric cancer recurrence within five years has been observed after surgical resection and following chemotherapy with therapy-resistant features. Therefore, the identification of new targets and treatment modalities for gastric cancer is of paramount importance. In this review we focus on the role of microRNAs with documented efficacy in preclinical xenograft models with respect to growth of human gastric cancer cells. We have identified 31 miRs (-10b, -19a, -19b, -20a, -23a/b, -25, -27a-3p, -92a, -93, -100, -106a, -130a, -135a, -135b-5p, -151-5p, -187, -199-3p, -215, -221-3p, -224, -340a, -382, -421, -425, -487a, -493, -532-3p, -575, -589, -664a-3p) covering 26 different targets which promote growth of gastric cancer cells in vitro and in vivo as xenografts. Five miRs (miRs -10b, 151-5p, -187, 532-3p and -589) additionally have an impact on metastasis. Thirteen of the identified miRs (-19b, -20a/b, -25, -92a, -106a, -135a, -187, -221-3p, -340a, -421, -493, -575 and -589) have clinical impact on worse prognosis in patients.

A Roadmap Toward the Definition of Actionable Tumor-Specific Antigens

The search for tumor-specific antigens (TSAs) has considerably accelerated during the past decade due to the improvement of proteogenomic detection methods. This provides new opportunities for the development of novel antitumoral immunotherapies to mount an efficient T cell response against one or multiple types of tumors. While the identification of mutated antigens originating from coding exons has provided relatively few TSA candidates, the possibility of enlarging the repertoire of targetable TSAs by looking at antigens arising from non-canonical open reading frames opens up interesting avenues for cancer immunotherapy. In this review, we outline the potential sources of TSAs and the mechanisms responsible for their expression strictly in cancer cells. In line with the heterogeneity of cancer, we propose that discrete families of TSAs may be enriched in specific cancer types.

Genome-Wide Analysis of Prognostic Alternative Splicing Signature and Splicing Factors in Lung Adenocarcinoma

Analysis of The Cancer Genome Atlas data revealed that alternative splicing (AS) events could serve as prognostic biomarkers in various cancer types. This study examined lung adenocarcinoma (LUAD) tissues for AS and assessed AS events as potential indicators of prognosis in our cohort. RNA sequencing and bioinformatics analysis were performed. We used SUPPA2 to analyze the AS profiles. Using univariate Cox regression analysis, overall survival (OS)-related AS events were identified. Genes relating to the OS-related AS events were imported into Cytoscape, and the CytoHubba application was run. OS-related splicing factors (SFs) were explored using the log-rank test. The relationship between the percent spliced-in value of the OS-related AS events and SF expression was identified by Spearman correlation analysis. We found 1957 OS-related AS events in 1151 genes, and most were protective factors. Alternative first exon splicing was the most frequent type of splicing event. The hub genes in the gene network of the OS-related AS events were FBXW11, FBXL5, KCTD7, UBB and CDC27. The area under the curve of the MIX prediction model was 0.847 for 5-year survival based on seven OS-related AS events. Overexpression of SFs CELF2 and SRSF5 was associated with better OS. We constructed a correlation network between SFs and OS-related AS events. In conclusion, we identified prognostic predictors using AS events that stratified LUAD patients into high- and low-risk groups. The discovery of the splicing networks in this study provides an insight into the underlying mechanisms.

Transcriptomic Analysis Revealed an Emerging Role of Alternative Splicing in Embryonal Tumor with Multilayered Rosettes

Embryonal tumor with multilayered rosettes (ETMR) is an aggressive and rare pediatric embryonal brain tumor. Amplification of C19MC microRNA cluster and expression of LIN28 are distinctive features of ETMR. Despite the increasing efforts to decipher ETMR, the biology remains poorly understood. To date, the role of aberrant alternative splicing in ETMR has not been thoroughly investigated. In the current study, a comprehensive analysis was performed on published unprocessed RNA-seq reads of tissue-matched ETMR and fetal controls datasets. Gene expression was quantified in samples using Kallisto/sleuth pipeline. For the alternative splicing analysis, STAR, SplAdder and rMATS were used. Functional enrichment analysis was subsequently performed using Metascape. The expression analysis identified a total of 3622 differentially expressed genes (DEGs) between ETMR and fetal controls while 1627 genes showed differential alternative splicing patterns. Interestingly, genes with significant alternative splicing events in ETMR were identified to be involved in signaling pathways such as ErbB, mTOR and MAPK pathways as well as ubiquitin-mediated proteolysis, cell cycle and autophagy. Moreover, up-regulated DEGs with alternative splicing events were involved in important biological processes including nuclear transport, regulation of cell cycle and regulation of Wnt signaling pathway. These findings highlight the role of aberrant alternative splicing in shaping the ETMR tumor landscape, and the identified pathways constitute potential therapeutic targets.

NONO Inhibits Lymphatic Metastasis of Bladder Cancer via Alternative Splicing of SETMAR

Bladder cancer patients with lymph node (LN) metastasis have an extremely poor prognosis and no effective treatment. The alternative splicing of precursor (pre-)mRNA participates in the progression of various tumors. However, the precise mechanisms of splicing factors and cancer-related variants in LN metastasis of bladder cancer remain largely unknown. The present study identified a splicing factor, non-POU domain-containing octamer-binding protein (NONO), that was significantly downregulated in bladder cancer tissues and correlated with LN metastasis status, tumor stage, and prognosis. Functionally, NONO markedly inhibited bladder cancer cell migration and invasion in vitro and LN metastasis in vivo. Mechanistically, NONO regulated the exon skipping of SETMAR by binding to its motif, mainly through the RRM2 domain. NONO directly interacted with splicing factor proline/glutamine rich (SFPQ) to regulate the splicing of SETMAR, and it induced metastasis suppression of bladder cancer cells. SETMAR-L overexpression significantly reversed the metastasis of NONO-knockdown bladder cancer cells, both in vitro and in vivo. The further analysis revealed that NONO-mediated SETMAR-L can induce H3K27me3 at the promotor of metastatic oncogenes and inhibit their transcription, ultimately resulting in metastasis suppression. Therefore, the present findings uncover the molecular mechanism of lymphatic metastasis in bladder cancer, which may provide novel clinical markers and therapeutic strategies for LN-metastatic bladder cancer.

VEGF-A in Cardiomyocytes and Heart Diseases

The vascular endothelial growth factor (VEGF), a homodimeric vasoactive glycoprotein, is the key mediator of angiogenesis. Angiogenesis, the formation of new blood vessels, is responsible for a wide variety of physio/pathological processes, including cardiovascular diseases (CVD). Cardiomyocytes (CM), the main cell type present in the heart, are the source and target of VEGF-A and express its receptors, VEGFR1 and VEGFR2, on their cell surface. The relationship between VEGF-A and the heart is double-sided. On the one hand, VEGF-A activates CM, inducing morphogenesis, contractility and wound healing. On the other hand, VEGF-A is produced by CM during inflammation, mechanical stress and cytokine stimulation. Moreover, high concentrations of VEGF-A have been found in patients affected by different CVD, and are often correlated with an unfavorable prognosis and disease severity. In this review, we summarized the current knowledge about the expression and effects of VEGF-A on CM and the role of VEGF-A in CVD, which are the most important cause of disability and premature death worldwide. Based on clinical studies on angiogenesis therapy conducted to date, it is possible to think that the control of angiogenesis and VEGF-A can lead to better quality and span of life of patients with heart disease.

Identification and expression analysis of novel splice variants of the human carcinoembryonic antigen-related cell adhesion molecule 19 (CEACAM19) gene using a high-throughput sequencing approach

Alternative splicing is commonly involved in carcinogenesis, being highly implicated in differential expression of cancer-related genes. Recent studies have shown that the human CEACAM19 gene is overexpressed in malignant breast and ovarian tumors, possessing significant biomarker attributes. In the present study, 3' rapid amplification of cDNA ends (3' RACE) and next-generation sequencing (NGS) were used for the detection and identification of novel CEACAM19 transcripts. Bioinformatical analysis of our NGS data revealed novel splice junctions between previously annotated exons and ultimately new exons. Next, fifteen novel CEACAM19 transcripts were identified with Sanger sequencing. Additionally, their expression profile was investigated in a wide panel of human cell lines, using nested PCR with variant-specific primers. The broad expression pattern of the CEACAM19 gene, along with the fact that its overexpression has previously been associated with ovarian and breast cancer progression, indicate the potential of novel CEACAM19 transcripts as putative diagnostic and/or prognostic biomarkers.

Anti-VEGF Drugs in the Treatment of Multiple Myeloma Patients

The interaction between the bone marrow microenvironment and plasma cells plays an essential role in multiple myeloma progression and drug resistance. The vascular endothelial growth factor (VEGF)/VEGF receptor (VEGFR) pathway in vascular endothelial cells activates and promotes angiogenesis. Moreover, VEGF activates and promotes vasculogenesis and vasculogenic mimicry when it interacts with VEGF receptors expressed in precursor cells and inflammatory cells, respectively. In myeloma bone marrow, VEGF and VEGF receptor expression are upregulated and hyperactive in the stromal and tumor cells. It has been demonstrated that several antiangiogenic agents can effectively target VEGF-related pathways in the preclinical phase. However, they are not successful in treating multiple myeloma, probably due to the vicarious action of other cytokines and signaling pathways. Thus, the simultaneous blocking of multiple cytokine pathways, including the VEGF/VEGFR pathway, may represent a valid strategy to treat multiple myeloma. This review aims to summarize recent advances in understanding the role of the VEGF/VEGFR pathway in multiple myeloma, and mainly focuses on the transcription pathway and on strategies that target this pathway.

In Silico Analysis of a Highly Mutated Gene in Cancer Provides Insight into Abnormal mRNA Splicing: Splicing Factor 3B Subunit 1K700E Mutant

Cancer is the second leading cause of death worldwide. The etiology of the disease has remained elusive, but mutations causing aberrant RNA splicing have been considered one of the significant factors in various cancer types. The association of aberrant RNA splicing with drug/therapy resistance further increases the importance of these mutations. In this work, the impact of the splicing factor 3B subunit 1 (SF3B1) K700E mutation, a highly prevalent mutation in various cancer types, is investigated through molecular dynamics simulations. Based on our results, K700E mutation increases flexibility of the mutant SF3B1. Consequently, this mutation leads to i) disruption of interaction of pre-mRNA with SF3B1 and p14, thus preventing proper alignment of mRNA and causing usage of abnormal 3' splice site, and ii) disruption of communication in critical regions participating in interactions with other proteins in pre-mRNA splicing machinery. We anticipate that this study enhances our understanding of the mechanism of functional abnormalities associated with splicing machinery, thereby, increasing possibility for designing effective therapies to combat cancer at an earlier stage.

Comparative proteogenomics profiling of non-small and small lung carcinoma cell lines using mass spectrometry

Background

Evidences indicated that non-small-cell lung cancer (NSCLC) and small-cell lung cancer (SCLC) might originate from the same cell type, which however ended up to be two different subtypes of lung carcinoma, requiring different therapeutic regimens. We aimed to identify the differences between these two subtypes of lung cancer by using integrated proteome and genome approaches.

Methods and Materials

Two representative cell lines for each lung cancer subtype were comparatively analysed by quantitative proteomics, and their corresponding transcriptomics data were obtained from the Gene Expression Omnibus database. The integrated analyses of proteogenomic data were performed to determine key differentially expressed proteins that were positively correlated between proteomic and transcriptomic data.

Result

The proteomics analysis revealed 147 differentially expressed proteins between SCLC and NSCLC from a total of 3,970 identified proteins. Combined with available transcriptomics data, we further confirmed 14 differentially expressed proteins including six known and eight new lung cancer related proteins that were positively correlated with their transcriptomics data. These proteins are mainly involved in cell migration, proliferation, and invasion.

Conclusion

The proteogenomic data on both NSCLC and SCLC cell lines presented in this manuscript is complementary to existing genomic and proteomic data related to lung cancers and will be crucial for a systems biology-level understanding of the molecular mechanism of lung cancers. The raw mass spectrometry data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD015270.

A Novel Prognostic Index Based on Alternative Splicing in Papillary Renal Cell Carcinoma

Background

Papillary renal cell carcinoma (pRCC) is a heterogeneous multifocal or isolated tumor with an invasive phenotype. Previous studies presented that alternative splicing, as a crucial posttranscriptional regulator in gene expression, is associated with tumorigenesis. However, the association between alternative splicing and pRCC has not been clarified

Methods

The RNA sequencing data and clinical information were downloaded from The Cancer Genome Atlas database and mRNA splicing profiles from TCGASpliceSeq. The percent spliced in data of alternative splicing merged with survival information was firstly calculated by univariate Cox regression analysis to screen for survival‐associated alternative splicing events, and survival‐associated alternative splicing events were then analyzed by Gene Ontology categories using Kyoto Encyclopedia of Genes and Genomes. Meanwhile, the least absolute shrinkage and selection operator Cox analysis and multivariate Cox analysis were performed to calculate the prognostic index for each alternative splicing type. In addition, clinical factors were introduced to assess the performance of prognostic index.

Results

A total of 4,084 candidate survival-associated alternative splicing events in 2,558 genes were screened out. Patients were divided into the low-risk group and the high-risk group based on the median prognostic index value. The Kaplan-Meier survival analysis (p < 0.05) and receiver operating characteristics curves (AUC>0.9) indicated that prognostic index was effective and stable for predicting the prognosis of pRCC patients. Furthermore, a regulatory network was constructed incorporating alternative splicing events and survival-associated splicing factors.

Conclusion

Our study provides new insights into the mechanism of alternative splicing events in tumorigenesis and their clinical potential for pRCC.

Survival-Associated Alternative Splicing Events in Pan-Renal Cell Carcinoma

Alternative splicing is an important modification process for the genome to generate mature mRNA by transcription, which has been found associated with survival in some tumors. However, systematic analysis of AS events in pan-renal cell carcinoma at the genome-wide level has been seldom conducted yet. In the current study, Upset plot and Venn plot were utilized to present the distribution characteristics of AS events. Those SREs were screened out with multivariate COX regression analyses, and functional enrichment analysis was performed to figure out potential pathways. ROC model was conducted to compare the efficiency of those potential SREs. A total of 2,169, 1,671, and 1,414 SREs were found in renal clear cell carcinoma (KIRC), renal chromophobe cell carcinoma (KICH), and renal papillary cell carcinoma (KIRP), respectively. Functional enrichment analysis results suggested possible mechanism such as changes in the branched-chain amino acid catabolic process due to SREs might play a key role in KIRC. The binary logistic regression equation based on the SREs had a good performance in each model compared to the single factor. The 5 year survival model presented that the AUC of the predicted probabilities in KIRC, KICH, and KIRP were 0.754, 1 and 0.841, and in the diagnostic model were 0.988, 0.970, and 0.999, respectively. Some AS types that were significantly different in pan-RCC and paracancerous tissues have also been discovered to play a role in carcinoma screening. To sum up, alternative splicing events significantly interfere with the prognosis of patients with pan-RCC and are capable as biomarkers for prognosis.

PTBP1-mediated regulation of AXL mRNA stability plays a role in lung tumorigenesis

AXL is expressed in many types of cancer and promotes cancer cell survival, metastasis and drug resistance. Here, we focus on identifying modulators that regulate AXL at the mRNA level. We have previously observed that the AXL promoter activity is inversely correlated with the AXL expression levels, suggesting that post-transcriptional mechanisms exist that down-regulate the expression of AXL mRNA. Here we show that the RNA binding protein PTBP1 (polypyrimidine tract-binding protein) directly targets the 5′-UTR of AXL mRNA in vitro and in vivo. Moreover, we also demonstrate that PTBP1, but not PTBP2, inhibits the expression of AXL mRNA and the RNA recognition motif 1 (RRM1) of PTBP1 is crucial for this interaction. To clarify how PTBP1 regulates AXL expression at the mRNA level, we found that, while the transcription rate of AXL was not significantly different, PTBP1 decreased the stability of AXL mRNA. In addition, over-expression of AXL may counteract the PTBP1-mediated apoptosis. Knock-down of PTBP1 expression could enhance tumor growth in animal models. Finally, PTBP1 was found to be negatively correlated with AXL expression in lung tumor tissues in Oncomine datasets and in tissue micro-array (TMA) analysis. In conclusion, we have identified a molecular mechanism of AXL expression regulation by PTBP1 through controlling the AXL mRNA stability. These findings may represent new thoughts alternative to current approaches that directly inhibit AXL signaling and may eventually help to develop novel therapeutics to avoid cancer metastasis and drug resistance.

Coding RNA Sequencing of Equine Endometrium during Maternal Recognition of Pregnancy

Equine maternal recognition of pregnancy (MRP) is a process whose signal remains unknown. During MRP the conceptus and endometrium communicate to attenuate prostaglandin F_2α (PGF) secretion, sparing the corpus luteum and maintaining progesterone production. Recognition of a mobile conceptus by the endometrium is critical by days 14–16 post-ovulation (PO), when endometrium produces PGF, initiating luteolysis. The objective of this study was to evaluate endometrial gene expression changes based upon pregnancy status via RNA sequencing. This experiment utilized a cross-over design with each mare serving as both a pregnant and non-mated control on days nine, 11, and 13 PO (n = 3/status/day). Mares were randomly assigned to collection day and pregnancy confirmed by terminal uterine lavage at the time of endometrial biopsy. Total RNA was isolated and libraries prepared using Illumina TruSeq RNA sample preparation kit. Reads were mapped and annotated using HISAT2 and Stringtie. Expression values were evaluated with DESEQ2 (P ≤ 0.05 indicated significance). On day nine, 11, and 13 there were 1435, 1435 and 916 significant transcripts, respectively. Multiple genes with splice variants had different expression patterns within the same day. These are the first data to evaluate the endometrial transcriptome during MRP on days nine, 11, and 13.

Identification of FUBP1 as a Long Tail Cancer Driver and Widespread Regulator of Tumor Suppressor and Oncogene Alternative Splicing

Comprehensive sequencing approaches have allowed for the identification of the most frequent contributors to cancer, known as drivers. They have also revealed a class of mutations in understudied, infrequently altered genes, referred to as "long tail" (LT) drivers. A key challenge has been to find clinically relevant LT drivers and to understand how they cooperate to drive disease. Here, we identified far upstream binding protein 1 (FUBP1) as an LT driver using an in vivo CRISPR screen. FUBP1 cooperates with other tumor suppressor genes to transform mammary epithelial cells by disrupting cellular differentiation and tissue architecture. Mechanistically, FUBP1 participates in regulating N6-methyladenosine (m6A) RNA methylation, and its loss leads to global changes in RNA splicing and widespread expression of aberrant driver isoforms. These findings suggest that somatic alteration of a single gene involved in RNA splicing and m6A methylation can produce the necessary panoply of contributors for neoplastic transformation.

Identification of novel alternative splice variants of the human L-DOPA decarboxylase (DDC) gene in human cancer cells, using high-throughput sequencing approaches

The human L-DOPA decarboxylase (DDC) is a gene that has been in the center of research attention in many laboratories the last decades, due to its major implication in various disorders, including many types of cancer. In the current work, we used in-house developed RACE and high-throughput sequencing approaches, in order to detect and identify novel DDC transcripts. Bioinformatic analysis revealed new alternative splicing events that support the existence of novel DDC transcripts. As a result, a total of 14 DDC splice variants were identified and their expression profile was investigated in a wide panel of human cancer cell lines. From all 14 novel DDC transcripts that were identified, 9 transcripts are predicted to encode new protein isoforms, while the remaining 5 are nonsense-mediated mRNA decay (NMD) candidates. Our results demonstrate that the human DDC gene undergoes complex processing leading to the figuration of multiple mRNA isoforms in cancer cells.

Novel splice variants of the human kallikrein-related peptidases 11 (KLK11) and 12 (KLK12), unraveled by next-generation sequencing technology

Tissue kallikrein, kallikrein-related peptidases (KLKs), and plasma kallikrein form the largest group of serine proteases in the human genome, sharing many structural and functional characteristics. In this study, we describe the molecular cloning of four novel splice variants of the human KLK11 and KLK12 genes, discovered by combining 3' rapid amplification of cDNA ends (3' RACE), next-generation sequencing (NGS) technology, advanced bioinformatic analysis and Sanger sequencing. Expression analysis of these new transcripts in cell lines originating from 17 cancerous and two normal tissues revealed the expression pattern of each transcript. These novel KLK11 and KLK12 splice variants represent new potential cancer biomarkers.

Hypoxia alters splicing of the cancer associated Fas gene

The removal of introns from mRNA precursors (pre-mRNAs) is an essential step in eukaryotic gene expression. The splicing machinery heavily contributes to biological complexity and especially to the ability of cells to adapt to altered cellular conditions. Hypoxia also plays a key role in the pathophysiology of many disease states. Recent studies have revealed that tumorigenesis and hypoxia involve large-scale alterations in alternative pre-mRNA splicing. Cancer associated Fas protein plays a central role in the physiological regulation of programmed cell death and has been implicated in the pathogenesis of various malignancies and diseases of the immune system. Fas pre-mRNA is alternatively spliced by excluding exon 6 to produce soluble Fas (sFas) protein that lacks a transmembrane domain and acts by inhibiting Fas mediated apoptosis. Another cancer related protein Rac1 plays an important regulatory role specifically in cells' motility, growth and survival. Rac pre-mRNA is alternatively spliced to produce Rac1b protein, which is upregulated in metastatic diseases. In the present study we, for the first time, show that anti-apoptotic Fas mRNA isoform formation is regulated by cellular microenvironment - hypoxia. Hypoxic microenvironment, however, does not influence Rac1 pre-mRNAs alternative splicing. Also our presented results indicate that splicing factors hnRNP A1 and SPF45, previously shown to regulate Fas alternative splicing in normoxic cells, are not involved in hypoxia dependent alternative Fas pre-mRNA splicing regulation in an amount dependent manner. Our observations on hypoxia dependent alternative Fas pre-mRNA splicing regulation indicate a probable involvement of other, yet unidentified splicing factors. Presented data also shows the contribution of pre-mRNA splicing to cell survival under unfavorable conditions.

Differential expression and functions of Ehm2 transcript variants in lung adenocarcinoma

Ehm2 [also known as erythrocyte membrane protein band 4.1‑like protein 4B (EPB41L4B)] is a member of the NF2/ERM/4.1 superfamily. The overexpression of Ehm2 has been observed in metastatic cancer cells. Through alternative splicing, the Ehm2 gene produces two transcript variants that encode the two different isoforms, Ehm2/1 and Ehm2/2. The biological functions of these different Ehm2 transcript variants remain unclear. The present study aimed to determine the expression of the Ehm2 variants in lung adenocarcinoma and their involvement in the disease progression of the patients. The expression of Ehm2 transcript variants in human lung adenocarcinoma tissues was analyzed using immunohistochemistry and western blot analysis. Ehm2 variants were overexpressed or knocked down in A549 human lung adenocarcinoma cells. The consequent effects of the genetic modifications on the cellular functions of lung cancer cells were then examined using in vitro cell viability, invasion and migration assays. The expression of epithelial‑mesenchymal transition (EMT)‑related markers was evaluated by western blot analysis in the cell models. The association of Ehm2 variant expression with patient survival was analyzed using Kaplan‑Meier survival analysis. The expression of Ehm2/1 was significantly decreased in lung cancers compared with the paired normal lung tissues (P<0.05), while the Ehm2/2 protein levels were higher in the tumors than in the paired normal lung tissues, although this was not statistically significant. The overexpression of Ehm2/1 exerted inhibitory effects, while the knockdown of Ehm2/1 promoted the growth, invasion and migration of A549 cells in vitro. Ehm2/2 was expressed at low levels in the A549 cells and the enforced expression of Ehm2/2 significantly increased the invasiveness and migration of the A549 cells. Immunofluorescence staining revealed that Ehm2/1 was confined to the plasma membrane, while Ehm2/2 was observed at both the plasma membrane and cytoplasm. The overexpression of Ehm2/1 resulted in the upregulation of the epithelial marker, E‑cadherin, and in the decreased expression of the mesenchymal markers, N‑cadherin and Snail1, while the knockdown of Ehm2/1 and the enforced expression of Ehm2/2 had the opposite effects on the protein levels of EMT‑related markers. Kaplan‑Meier survival analysis revealed that higher Ehm2/1 transcript levels were associated with the longer survival of patients with lung adenocarcinoma, while the lower expression of Ehm2/2 exhibited a similar association with patient survival. Taken together, the two Ehm2 variants appear to be differentially expressed in lung adenocarcinoma. Ehm2/1 may function as a putative tumor suppressor in the disease progression of lung adenocarcinoma, while Ehm2/2 may have an opposite function.

Cancer-Specifically Re-Spliced TSG101 mRNA Promotes Invasion and Metastasis of Nasopharyngeal Carcinoma

TSG101 (Tumor susceptibility 101) gene and its aberrantly spliced isoform, termed TSG101∆154-1054, are tightly linked to tumorigenesis in various cancers. The aberrant TSG101∆154-1054 mRNA is generated from cancer-specific re-splicing of mature TSG101 mRNA. The TSG101∆154-1054 protein protects the full-length TSG101 protein from ubiquitin-mediated degradation, implicating TSG101∆154-1054 protein in the progression of cancer. Here, we confirmed that the presence of TSG101∆154-1054 mRNA indeed caused an accumulation of the TSG101 protein in biopsies of human nasopharyngeal carcinoma (NPC), which was recapitulated by the overexpression of TSG101∆154-1054 in the NPC cell line TW01. We demonstrate the potential function of the TSG101∆154-1054 protein in the malignancy of human NPC with scratch-wound healing and transwell invasion assays. By increasing the stability of the TSG101 protein, TSG101∆154-1054 specifically enhanced TSG101-mediated TW01 cell migration and invasion, suggesting the involvement in NPC metastasis in vivo. This finding sheds light on the functional significance of TSG101∆154-1054 generation via re-splicing of TSG101 mRNA in NPC metastasis and hints at its potential importance as a therapeutic target.

Polypyrimidine tract binding protein 1 promotes lymphatic metastasis and proliferation of bladder cancer via alternative splicing of MEIS2 and PKM

Lymph node (LN) metastasis is the leading cause of bladder cancer-related mortality. Splicing factors facilitate cancer progression by modulating oncogenic variants, but it is unclear whether and how splicing factors regulate bladder cancer LN metastasis. In this study, Polypyrimidine tract binding protein 1 (PTBP1) expression was found to relate to bladder cancer LN metastasis, and was positively correlated with LN metastasis status, tumor stage, histological grade, and poor patient prognosis. Functional assays demonstrated that PTBP1 promoted bladder cancer cell migration, invasion, and proliferation in vitro, as well as LN metastasis and tumor growth in vivo. Mechanistic investigations revealed that PTBP1 upregulated MEIS2-L variant to promote metastasis and increased expression of PKM2 variant to enhance proliferation by modulating alternative mRNA splicing. Moreover, overexpression of MEIS2-L or PKM2 could rescue the oncogenic abilities of bladder cancer cells and the expression of MMP9 or CCND1 respectively after PTBP1 knockdown. In conclusion, our data demonstrate that PTBP1 induces bladder cancer LN metastasis and proliferation through an alternative splicing mechanism. PTBP1 may serve as a novel prognostic marker and therapeutic target for LN-metastatic bladder cancer.

CD44 splice isoform switching determines breast cancer stem cell state

Zhang et al. show that manipulating the splicing regulator ESRP1 to shift alternative splicing from splice isoform CD44v to CD44s leads to an induction of cancer stem cell properties.

Although changes in alternative splicing have been observed in cancer, their functional contributions still remain largely unclear. Here we report that splice isoforms of the cancer stem cell (CSC) marker CD44 exhibit strikingly opposite functions in breast cancer. Bioinformatic annotation in patient breast cancer in The Cancer Genome Atlas (TCGA) database reveals that the CD44 standard splice isoform (CD44s) positively associates with the CSC gene signatures, whereas the CD44 variant splice isoforms (CD44v) exhibit an inverse association. We show that CD44s is the predominant isoform expressed in breast CSCs. Elimination of the CD44s isoform impairs CSC traits. Conversely, manipulating the splicing regulator ESRP1 to shift alternative splicing from CD44v to CD44s leads to an induction of CSC properties. We further demonstrate that CD44s activates the PDGFRβ/Stat3 cascade to promote CSC traits. These results reveal CD44 isoform specificity in CSC and non-CSC states and suggest that alternative splicing provides functional gene versatility that is essential for distinct cancer cell states and thus cancer phenotypes.

Networks of mRNA Processing and Alternative Splicing Regulation in Health and Disease

mRNA processing events introduce an intricate layer of complexity into gene expression processes, supporting a tremendous level of diversification of the genome's coding and regulatory potential, particularly in vertebrate species. The recent development of massive parallel sequencing methods and their adaptation to the identification and quantification of different RNA species and the dynamics of mRNA metabolism and processing has generated an unprecedented view over the regulatory networks that are established at this level, which contribute to sustain developmental, tissue specific or disease specific gene expression programs. In this chapter, we provide an overview of the recent evolution of transcriptome profiling methods and the surprising insights that have emerged in recent years regarding distinct mRNA processing events - from the 5' end to the 3' end of the molecule.

Extensive Alternative Splicing of KIR Transcripts

The killer-cell Ig-like receptors (KIR) form a multigene entity involved in modulating immune responses through interactions with MHC class I molecules. The complexity of the KIR cluster is reflected by, for instance, abundant levels of allelic polymorphism, gene copy number variation, and stochastic expression profiles. The current transcriptome study involving human and macaque families demonstrates that KIR family members are also subjected to differential levels of alternative splicing, and this seems to be gene dependent. Alternative splicing may result in the partial or complete skipping of exons, or the partial inclusion of introns, as documented at the transcription level. This post-transcriptional process can generate multiple isoforms from a single KIR gene, which diversifies the characteristics of the encoded proteins. For example, alternative splicing could modify ligand interactions, cellular localization, signaling properties, and the number of extracellular domains of the receptor. In humans, we observed abundant splicing for KIR2DL4, and to a lesser extent in the lineage III KIR genes. All experimentally documented splice events are substantiated by in silico splicing strength predictions. To a similar extent, alternative splicing is observed in rhesus macaques, a species that shares a close evolutionary relationship with humans. Splicing profiles of Mamu-KIR1D and Mamu-KIR2DL04 displayed a great diversity, whereas Mamu-KIR3DL20 (lineage V) is consistently spliced to generate a homolog of human KIR2DL5 (lineage I). The latter case represents an example of convergent evolution. Although just a single KIR splice event is shared between humans and macaques, the splicing mechanisms are similar, and the predicted consequences are comparable. In conclusion, alternative splicing adds an additional layer of complexity to the KIR gene system in primates, and results in a wide structural and functional variety of KIR receptors and its isoforms, which may play a role in health and disease.

Aberrant RNA Splicing in Cancer and Drug Resistance

More than 95% of the 20,000 to 25,000 transcribed human genes undergo alternative RNA splicing, which increases the diversity of the proteome. Isoforms derived from the same gene can have distinct and, in some cases, opposing functions. Accumulating evidence suggests that aberrant RNA splicing is a common and driving event in cancer development and progression. Moreover, aberrant splicing events conferring drug/therapy resistance in cancer is far more common than previously envisioned. In this review, aberrant splicing events in cancer-associated genes, namely BCL2L1, FAS, HRAS, CD44, Cyclin D1, CASP2, TMPRSS2-ERG, FGFR2, VEGF, AR and KLF6, will be discussed. Also highlighted are the functional consequences of aberrant splice variants (BCR-Abl35INS, BIM-γ, IK6, p61 BRAF V600E, CD19-∆2, AR-V7 and PIK3CD-S) in promoting resistance to cancer targeted therapy or immunotherapy. To overcome drug resistance, we discuss opportunities for developing novel strategies to specifically target the aberrant splice variants or splicing machinery that generates the splice variants. Therapeutic approaches include the development of splice variant-specific siRNAs, splice switching antisense oligonucleotides, and small molecule inhibitors targeting splicing factors, splicing factor kinases or the aberrant oncogenic protein isoforms.

Insulin Receptor Isoforms in Cancer

The insulin receptor (IR) mediates both metabolic and mitogenic effects especially when overexpressed or in clinical conditions with compensatory hyperinsulinemia, due to the metabolic pathway resistance, as obesity diabetes. In many cancers, IR is overexpressed preferentially as IR-A isoform, derived by alternative splicing of exon 11. The IR-A overexpression, and the increased IR-A:IR-B ratio, are mechanisms that promote the mitogenic response of cancer cells to insulin and IGF-2, which is produced locally by both epithelial and stromal cancer cells. In cancer IR-A, isoform predominance may occur for dysregulation at both mRNA transcription and post-transcription levels, including splicing factors, non-coding RNAs and protein degradation. The mechanisms that regulate IR isoform expression are complex and not fully understood. The IR isoform overexpression may play a role in cancer cell stemness, in tumor progression and in resistance to target therapies. From a clinical point of view, the IR-A overexpression in cancer may be a determinant factor for the resistance to IGF-1R target therapies for this issue. IR isoform expression in cancers may have the meaning of a predictive biomarker and co-targeting IGF-1R and IR-A may represent a new more efficacious treatment strategy.

Roles of Smads Family and Alternative Splicing Variants of Smad4 in Different Cancers

Transforming Growth Factor β (TGF-β) is one of the most common secretory proteins which are recognized by membrane receptors joined to transcription regulatory factor. TGF-β signals are transduced by the Smads family that regulate differentiation, proliferation, early growth, apoptosis, homeostasis, and tumor development. Functional study of TGF-β signaling pathway and Smads role is vital for certain diseases such as cancer. Alternative splicing produces a diverse range of protein isoforms with unique function and the ability to react differently with various pharmaceutical products. This review organizes to describe the general study of Smads family, the process of alternative splicing, the general aspect of alternative splicing of Smad4 in cancer and the possible use of spliceoforms for the diagnosis and therapeutic purpose. The main aim and objective of this article are to highlight some particular mechanisms involving in alternatives splicing of cancer and also to demonstrate new evidence about alternative splicing in different steps given cancer initiation and progression.

Mutually exclusive acetylation and ubiquitylation of the splicing factor SRSF5 control tumor growth

Most tumor cells take up more glucose than normal cells. Splicing dysregulation is one of the molecular hallmarks of cancer. However, the role of splicing factor in glucose metabolism and tumor development remains poorly defined. Here, we show that upon glucose intake, the splicing factor SRSF5 is specifically induced through Tip60-mediated acetylation on K125, which antagonizes Smurf1-mediated ubiquitylation. SRSF5 promotes the alternative splicing of CCAR1 to produce CCAR1S proteins, which promote tumor growth by enhancing glucose consumption and acetyl-CoA production. Conversely, upon glucose starvation, SRSF5 is deacetylated by HDAC1, and ubiquitylated by Smurf1 on the same lysine, resulting in proteasomal degradation of SRSF5. The CCAR1L proteins accumulate to promote apoptosis. Importantly, SRSF5 is hyperacetylated and upregulated in human lung cancers, which correlates with increased CCAR1S expression and tumor progression. Thus, SRSF5 responds to high glucose to promote cancer development, and SRSF5-CCAR1 axis may be valuable targets for cancer therapeutics.

PTBP1 promotes the growth of breast cancer cells through the PTEN/Akt pathway and autophagy

Invasion and migration is the hallmark of malignant tumors as well as the major cause for breast cancer death. The polypyrimidine tract binding, PTB, protein serves as an important model for understanding how RNA binding proteins affect proliferation and invasion and how changes in the expression of these proteins can control complex programs of tumorigenesis. We have investigated some roles of polypyrimidine tract binding protein 1 (PTBP1) in human breast cancer. We found that PTBP1 was upregulated in breast cancer tissues compared with normal tissues and the same result was confirmed in breast cancer cell lines. Knockdown of PTBP1 substantially inhibited tumor cell growth, migration, and invasion. These results suggest that PTBP1 is associated with breast tumorigenesis and appears to be required for tumor cell growth and maintenance of metastasis. We further analyzed the relationship between PTBP1 and clinicopathological parameters and found that PTBP1 was correlated with her‐2 expression, lymph node metastasis, and pathological stage. This will be a novel target for her‐2(⁺) breast cancer. PTBP1 exerts these effects, in part, by regulating the phosphatase and tensin homolog‐phosphatidylinositol‐4,5‐bisphosphate 3‐kinase/protein kinase B (PTEN‐PI3K/Akt) pathway and autophagy, and consequently alters cell growth and contributes to the invasion and metastasis.

Discovery and expression analysis of novel transcripts of the human SR-related CTD-associated factor 1 (SCAF1) gene in human cancer cells using Next-Generation Sequencing

The human SR-related CTD associated factor 1 (SCAF1) gene is a new member of the human SR (Ser/Arg-rich) superfamily of pre-mRNA splicing factors, which has been discovered and cloned by members of our lab. SCAF1 interacts with the CTD domain of the RNA polymerase II polypeptide A and is firmly involved in pre-mRNA splicing. Although it was found to be expressed widely in multiple human tissues, its mRNA levels vary a lot. The significant relation of SCAF1 with cancer has been confirmed by many studies, since SCAF1 mRNA transcript was found to be overexpressed in breast and ovarian tumors, confirming its significant prognostic value as a cancer biomarker in both these human malignancies. In this study, we describe the discovery and cloning of fifteen novel transcripts of the human SCAF1 gene (SCAF1 v.2 - v.16), using nested PCR and NGS technology. In detail, extensive bioinformatic analysis revealed that these novel SCAF1 splice variants comprise a total of nine novel alternative splicing events between the annotated exons of the gene, thus producing seven novel SCAF1 transcripts with open-reading frames, which are predicted to encode novel SCAF1 isoforms and eight novel SCAF1 transcripts with premature termination codons that are likely long non-coding RNAs. Additionally, a novel 3' UTR was discovered and cloned using nested 3' RACE and was validated with Sanger sequencing. In order to validate the NGS findings as well as to investigate the expression profile of each novel transcript, RT-PCR experiments were carried out with the use of variant-specific primers. Since SCAF1 is implicated in many human malignancies, qualifying as a potential biomarker, the quantification of the presented novel transcripts in human samples may have clinical applications in different types of cancer.

Modulation of alternative splicing induced by paclitaxel in human lung cancer

Paclitaxel is utilized as the first-line chemotherapeutic regimen for the majority of advanced non-small-cell lung carcinoma. However, whether paclitaxel could suppress cancer progression through modulating RNA alternative splicing remains largely unknown. Here, we demonstrated the effects of paclitaxel on cell proliferation inhibition, cell cycle arrest, and apoptosis. Mechanistically, paclitaxel leads to transcriptional alteration of networks involved in DNA replication and repair, chromosome segregation, chromatin silencing at rDNA, and mitosis at the transcriptional level. Moreover, paclitaxel regulates a number of cancer-associated RNA alternative splicing events, including genes involved in cellular response to DNA damage stimulus, preassembly of GPI anchor in ER membrane, transcription, and DNA repair. In particular, paclitaxel modulates the splicing of ECT2, a key factor involved in the regulation of cytokinesis. Briefly, paclitaxel favors the production of ECT2-S, the short splicing isoforms of ECT2, thereby inhibiting cancer cell proliferation. Our study provides mechanistic insights of paclitaxel on RNA alternative splicing regulation, thus to offer a potential novel route for paclitaxel to inhibit cancer progression.

Molecular Pharmacology of VEGF-A Isoforms: Binding and Signalling at VEGFR2

Vascular endothelial growth factor-A (VEGF-A) is a key mediator of angiogenesis, signalling via the class IV tyrosine kinase receptor family of VEGF Receptors (VEGFRs). Although VEGF-A ligands bind to both VEGFR1 and VEGFR2, they primarily signal via VEGFR2 leading to endothelial cell proliferation, survival, migration and vascular permeability. Distinct VEGF-A isoforms result from alternative splicing of the Vegfa gene at exon 8, resulting in VEGFxxxa or VEGFxxxb isoforms. Alternative splicing events at exons 5⁻7, in addition to recently identified posttranslational read-through events, produce VEGF-A isoforms that differ in their bioavailability and interaction with the co-receptor Neuropilin-1. This review explores the molecular pharmacology of VEGF-A isoforms at VEGFR2 in respect to ligand binding and downstream signalling. To understand how VEGF-A isoforms have distinct signalling despite similar affinities for VEGFR2, this review re-evaluates the typical classification of these isoforms relative to the prototypical, “pro-angiogenic” VEGF165a. We also examine the molecular mechanisms underpinning the regulation of VEGF-A isoform signalling and the importance of interactions with other membrane and extracellular matrix proteins. As approved therapeutics targeting the VEGF-A/VEGFR signalling axis largely lack long-term efficacy, understanding these isoform-specific mechanisms could aid future drug discovery efforts targeting VEGF receptor pharmacology.

Discovery of novel transcripts of the human tissue kallikrein (KLK1) and kallikrein-related peptidase 2 (KLK2) in human cancer cells, exploiting Next-Generation Sequencing technology

Tissue kallikrein, kallikrein-related peptidases (KLKs), and plasma kallikrein form the largest group of serine proteases in the human genome, sharing many structural and functional properties. Several KLK transcripts have been found aberrantly expressed in numerous human malignancies, confirming their prognostic or/and diagnostic values. However, the process of alternative splicing can now be studied in-depth due to the development of Next-Generation Sequencing (NGS). In the present study, we used NGS to discover novel transcripts of the KLK1 and KLK2 genes, after nested touchdown PCR. Bioinformatics analysis and PCR experiments revealed a total of eleven novel KLK transcripts (two KLK1 and nine KLK2 transcripts). In addition, the expression profiles of each novel transcript were investigated with nested PCR experiments using variant-specific primers. Since KLKs are implicated in human malignancies, qualifying as potential biomarkers, the quantification of the presented novel transcripts in human samples may have clinical applications in different types of cancer.

High expression of PTBP1 promote invasion of colorectal cancer by alternative splicing of cortactin

Polypyrimidine tract-binding protein 1 (PTBP1) involving in almost all steps of mRNA regulation including alternative splicing metabolism during tumorigenesis due to its RNA-binding activity. Initially, we found that high expressed PTBP1 and poor prognosis was interrelated in colorectal cancer (CRC) patients with stages II and III CRC, which widely different in prognosis and treatment, by immunohistochemistry. PTBP1 was also upregulated in colon cancer cell lines. In our study, knockdown of PTBP1 by siRNA transfection decreased cell proliferation and invasion in vitro. Denovirus shRNA knockdown of PTBP1 inhibited colorectal cancer growth in vivo. Furthermore, PTBP1 regulates alternative splicing of many target genes involving in tumorgenesis in colon cancer cells. We confirmed that the splicing of cortactin exon 11 which was only contained in cortactin isoform-a, as a PTBP1 target. Knockdown of PTBP1 decreased the expression of cortactin isoform-a by exclusion of exon 11. Also the mRNA levels of PTBP1 and cortactin isoform-a were cooperatively expressed in colorectal cancer tissues. Knocking down cortactin isoform-a significantly decreased cell migration and invasion in colorectal cancer cells. Overexpression of cortactin isoform-a could rescue PTBP1-knockdown effect of cell motility. In summary the study revealed that PTBP1 facilitates colorectal cancer migration and invasion activities by inclusion of cortactin exon 11.

Signaling Pathways Driving Aberrant Splicing in Cancer Cells

Aberrant profiles of pre-mRNA splicing are frequently observed in cancer. At the molecular level, an altered profile results from a complex interplay between chromatin modifications, the transcriptional elongation rate of RNA polymerase, and effective binding of the spliceosome to the generated transcripts. Key players in this interplay are regulatory splicing factors (SFs) that bind to gene-specific splice-regulatory sequence elements. Although mutations in genes of some SFs were described, a major driver of aberrant splicing profiles is oncogenic signal transduction pathways. Signaling can affect either the transcriptional expression levels of SFs or the post-translational modification of SF proteins, and both modulate the ratio of nuclear versus cytoplasmic SFs in a given cell. Here, we will review currently known mechanisms by which cancer cell signaling, including the mitogen-activated protein kinases (MAPK), phosphatidylinositol 3 (PI3)-kinase pathway (PI3K) and wingless (Wnt) pathways but also signals from the tumor microenvironment, modulate the activity or subcellular localization of the Ser/Arg rich (SR) proteins and heterogeneous nuclear ribonucleoproteins (hnRNPs) families of SFs.

Molecular cloning of novel transcripts of human kallikrein-related peptidases 5, 6, 7, 8 and 9 (KLK5 - KLK9), using Next-generation sequencing

Alternative splicing of cancer-related genes is a common cellular mechanism accounting for cancer cell transcriptome complexity and affecting cell cycle control, proliferation, apoptosis, angiogenesis, invasion, and metastasis. In this study, we describe the discovery and molecular cloning of thirty novel transcripts of the human KLK5, KLK6, KLK7, KLK8 and KLK9 genes, using 3′ rapid amplification of cDNA ends (3′ RACE) and NGS technology, as well as their expression analysis in many established cell lines, originating from several distinct cancerous and normal tissues. Extensive bioinformatic analysis revealed novel splice variants of these five members of the KLK family, comprising entirely new exons, previously unknown boundaries of the already annotated exons (extensions and truncations) as well as alternative splicing events between these exons. Nested RT-PCR in a panel of human cell lines originating from seventeen cancerous and two normal tissues with the use of variant-specific pairs of primers was carried out for expression analysis of these novel splice variants, and Sanger sequencing of the respective amplicons confirmed our NGS results. Given that some splice variants of KLK family members possess clinical value, novel alternatively spliced transcripts appear as new candidate biomarkers for diagnostic and/or prognostic purposes and as targets for therapeutic strategies.

Insulin Receptor Isoforms in Physiology and Disease: An Updated View

The insulin receptor (IR) gene undergoes differential splicing that generates two IR isoforms, IR-A and IR-B. The physiological roles of IR isoforms are incompletely understood and appear to be determined by their different binding affinities for insulin-like growth factors (IGFs), particularly for IGF-2. Predominant roles of IR-A in prenatal growth and development and of IR-B in metabolic regulation are well established. However, emerging evidence indicates that the differential expression of IR isoforms may also help explain the diversification of insulin and IGF signaling and actions in various organs and tissues by involving not only different ligand-binding affinities but also different membrane partitioning and trafficking and possibly different abilities to interact with a variety of molecular partners. Of note, dysregulation of the IR-A/IR-B ratio is associated with insulin resistance, aging, and increased proliferative activity of normal and neoplastic tissues and appears to sustain detrimental effects. This review discusses novel information that has generated remarkable progress in our understanding of the physiology of IR isoforms and their role in disease. We also focus on novel IR ligands and modulators that should now be considered as an important strategy for better and safer treatment of diabetes and cancer and possibly other IR-related diseases.

Alternative splicing discriminates molecular subtypes and has prognostic impact in diffuse large B-cell lymphoma

Effect of alternative splicing (AS) on diffuse large B-cell lymphoma (DLBCL) pathogenesis and survival has not been systematically addressed. Here, we compared differentially expressed genes and exons in association with survival after chemoimmunotherapy, and between germinal center B-cell like (GCB) and activated B-cell like (ABC) DLBCLs. Genome-wide exon array-based screen was performed from samples of 38 clinically high-risk patients who were treated in a Nordic phase II study with dose-dense chemoimmunotherapy and central nervous system prophylaxis. The exon expression profile separated the patients according to molecular subgroups and survival better than the gene expression profile. Pathway analyses revealed enrichment of AS genes in inflammation and adhesion-related processes, and in signal transduction, such as phosphatidylinositol signaling system and adenosine triphosphate binding cassette transporters. Altogether, 49% of AS-related exons were protein coding, and domain prediction showed 28% of such exons to include a functional domain, such as transmembrane helix domain or phosphorylation sites. Validation in an independent cohort of 92 DLBCL samples subjected to RNA-sequencing confirmed differential exon usage of selected genes and association of AS with molecular subtypes and survival. The results indicate that AS events are able to discriminate GCB and ABC DLBCLs and have prognostic impact in DLBCL.