Alternative splicing and tumor progression

MAGOH promotes gastric cancer progression via hnRNPA1 expression inhibition-mediated RONΔ160/PI3K/AKT signaling pathway activation

BACKGROUND

Gastric cancer (GC) is associated with high mortality and heterogeneity and poses a great threat to humans. Gene therapies for the receptor tyrosine kinase RON and its spliceosomes are attracting increasing amounts of attention due to their unique characteristics. However, little is known about the mechanism involved in the formation of the RON mRNA alternative spliceosome RONΔ160.

METHODS

Fourteen human GC tissue samples and six normal gastric tissue samples were subjected to label-free relative quantitative proteomics analysis, and MAGOH was identified as a candidate protein for subsequent studies. The expression of MAGOH in clinical specimens was verified by quantitative real-time PCR and western blotting. We then determined the biological function of MAGOH in GC through in vitro and in vivo experiments. RNA pulldown, RNA sequencing and RNA immunoprecipitation (RIP) were subsequently conducted to uncover the underlying mechanism by which MAGOH regulated the formation of RONΔ160.

RESULTS

Proteomic analysis revealed that MAGOH, which is located at key nodes and participates in RNA processing and mRNA splicing, was upregulated in GC tissue and GC cell lines and was associated with poor prognosis. Functional analysis showed that MAGOH promoted the proliferation, migration and invasion of GC cells in vitro and in vivo. Mechanistically, MAGOH inhibited the expression of hnRNPA1 and reduced the binding of hnRNPA1 to RON mRNA, thereby promoting the formation of RONΔ160 to activate the PI3K/AKT signaling pathway and consequently facilitating GC progression.

CONCLUSIONS

Our study revealed that MAGOH could promote the formation of RONΔ160 and activate the PI3K/AKT signaling pathway through the inhibition of hnRNPA1 expression. We elucidate a novel mechanism and potential therapeutic targets for the growth and metastasis of GC based on the MAGOH-RONΔ160 axis, and these findings have important guiding significance and clinical value for the future development of effective therapeutic strategies for GC.

Plant homeodomain-finger protein 5A: A key player in cancer progression

PHF5A is a member of the zinc-finger proteins. To advance knowledge on their role in carcinogenesis, data from experimental studies, animal models and clinical studies in different tumorigenesis have been reviewed. Furthermore, PHF5A as an oncogenic function, is frequently high expressed in tumor cells and a potential prognostic marker for different cancers. PHF5A is implicated in the regulation of cancer cell proliferation, invasion, migration and metastasis. Knockdown of PHF5A prevented the invasion and metastasis of tumor cells. Here, the role of PHF5A in different cancers and their possible mechanism in relation to recent literature is reviewed and discussed. There is an open promising perspective to their therapeutic management for different cancer types.

Human histone H1 variants impact splicing outcome by controlling RNA polymerase II elongation

Histones shape chromatin structure and the epigenetic landscape. H1, the most diverse histone in the human genome, has 11 variants. Due to the high structural similarity between the H1s, their unique functions in transferring information from the chromatin to mRNA-processing machineries have remained elusive. Here, we generated human cell lines lacking up to five H1 subtypes, allowing us to characterize the genomic binding profiles of six H1 variants. Most H1s bind to specific sites, and binding depends on multiple factors, including GC content. The highly expressed H1.2 has a high affinity for exons, whereas H1.3 binds intronic sequences. H1s are major splicing regulators, especially of exon skipping and intron retention events, through their effects on the elongation of RNA polymerase II (RNAPII). Thus, H1 variants determine splicing fate by modulating RNAPII elongation.

Identification of Intron Retention in the Slc16a3 Gene Transcript Encoding the Transporter MCT4 in the Brain of Aged and Alzheimer-Disease Model (APPswePS1dE9) Mice

The monocarboxylate transporter 4 (MCT4; Slc16a3) is expressed in the central nervous system, notably by astrocytes. It is implicated in lactate release and the regulation of glycolytic flux. Whether its expression varies during normal and/or pathological aging is unclear. As the presence of its mature transcript in the brain of young and old mice was determined, an unexpectedly longer RT-PCR fragment was detected in the mouse frontal cortex and hippocampus at 12 vs. 3 months of age. Cultured astrocytes expressed the expected 516 base pair (bp) fragment but treatment with IL-1β to mimic inflammation as can occur during aging led to the additional expression of a 928 bp fragment like that seen in aged mice. In contrast, cultured pericytes (a component of the blood-brain barrier) only exhibited the 516 bp fragment. Intriguingly, cultured endothelial cells constitutively expressed both fragments. When RT-PCR was performed on brain subregions of an Alzheimer mouse model (APPswePS1dE9), no fragment was detected at 3 months, while only the 928 bp fragment was present at 12 months. Sequencing of MCT4 RT-PCR products revealed the presence of a remaining intron between exon 2 and 3, giving rise to the longer fragment detected by RT-PCR. These results unravel the existence of intron retention for the MCT4 gene in the central nervous system. Such alternative splicing appears to increase with age in the brain and might be prominent in neurodegenerative diseases such as Alzheimer's disease. Hence, further studies in vitro and in vivo of intron 2 retention in the Slc16a3 gene transcript are required for adequate characterization concerning the biological roles of Slc16a3 isoforms in the context of aging and Alzheimer's disease pathology.

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.

A Splicing Transcriptome-Wide Association Study Identifies Candidate Altered Splicing for Prostate Cancer Risk

Prostate cancer (PCa) represents a huge public health burden among men. Many susceptibility genetic factors for PCa still remain unknown. In this study, we performed a large splicing transcriptome-wide association study (spTWAS) using three modeling strategies to develop alternative splicing genetic prediction models for identifying novel susceptibility loci and splicing introns for PCa risk by assessing 79,194 cases and 61,112 controls of European ancestry in the PRACTICAL, CRUK, CAPS, BPC3, and PEGASUS consortia. We identified 120 splicing introns of 97 genes showing an association with PCa risk at false discovery rate (FDR)-corrected threshold (FDR <0.05). Of them, 33 genes were enriched in PCa-related diseases and function categories. Fine-mapping analysis suggested that 21 splicing introns of 19 genes were likely causally associated with PCa risk. Thirty-five splicing introns of 34 novel genes were identified to be related to PCa susceptibility for the first time, and 11 of the genes were enriched in a cancer-related network. Our study identified novel loci and splicing introns associated with PCa risk, which can improve our understanding of the etiology of this common malignancy.

Novel thrombospondin-1 transcript exhibits distinctive expression and activity in thyroid tumorigenesis

Thrombospondin 1 (TSP1) is known for its cell-specific functions in cancer progression, such as proliferation and migration. It contains 22 exons that may potentially produce several different transcripts. Here, we identified TSP1V as a novel TSP1-splicing variant produced by intron retention (IR) in human thyroid cancer cells and tissues. We observed that TSP1V functionally inhibited tumorigenesis contrary to TSP1 wild-type, as identified in vivo and in vitro. These activities of TSP1V are caused by inhibiting phospho-Smad and phospho-focal adhesion kinase. Reverse transcription polymerase chain reaction and minigene experiments revealed that some phytochemicals/non-steroidal anti-inflammatory drugs enhanced IR. We further found that RNA-binding motif protein 5 (RBM5) suppressed IR induced by sulindac sulfide treatment. Additionally, sulindac sulfide reduced phospho-RBM5 levels in a time-dependent manner. Furthermore, trans-chalcone demethylated TSP1V, thereby preventing methyl-CpG-binding protein 2 binding to TSP1V gene. In addition, TSP1V levels were significantly lower in patients with differentiated thyroid carcinoma than in those with benign thyroid nodule, indicating its potential application as a diagnostic biomarker in tumor progression.

SERBP1 affects the apoptotic level by regulating the expression and alternative splicing of cellular and metabolic process genes in HeLa cells

Background

RNA-binding proteins (RBPs) have important roles in orchestrating posttranscriptional regulation and modulating many tumorigenesis events. SERBP1 has been recognized as an important regulator in multiple cancers, while it remains unclear whether SERBP1-regulated gene expression at the transcriptome-wide level is significantly correlated with tumorigenesis.

Methods

We overexpressed SERBP1 in HeLa cells and explored whether SERBP1 overexpression (SERBP1-OE) affects the proliferation and apoptosis of HeLa cells. We analyzed the transcriptome-wide gene expression changes and alternative splicing changes mediated by SERBP1-OE using the transcriptome sequencing method (RNA-seq). RT-qPCR was conducted to assay SERBP1-regulated alternative splicing.

Results

SERBP1-OE induced the apoptosis of HeLa cells. The downregulated genes were strongly enriched in the cell proliferation and apoptosis pathways according to the GO analysis, including FOS, FOSB, PAK6 and RAB26. The genes undergoing at least one SERBP1-regulated alternative splicing event were enriched in transcriptional regulation, suggesting a mechanism of the regulation of gene expression, and in pyruvate and fatty acid metabolic processes critical for tumorigenesis events. The SERBP1-regulated alternative splicing of ME3, LPIN3, CROT, PDP1, SLC27A1 and ALKBH7 was validated by RT-qPCR analysis.

Conclusions

We for the first time demonstrated the cellular function and molecular targets of SERBP1 in HeLa cells at transcriptional and post-transcriptional levels. The SERBP1-regulated gene expression and alternative splicing networks revealed by this study provide important information for exploring the functional roles and regulatory mechanisms of SERBP1 in cancer development and progression.

SRSF6 Regulates the Alternative Splicing of the Apoptotic Fas Gene by Targeting a Novel RNA Sequence

Simple Summary

Alternative splicing (AS) produces multiple mRNA isoforms from a gene to make a large number of proteins. Fas (Apo-1/CD95) pre-mRNA, a member of TNF receptor family that mediates apoptosis, can generate pro-apoptotic and anti-apoptotic proteins through AS. Here, we identified SRSF6 as an essential regulator protein in Fas AS. We further located a new functional target sequence of SRSF6 in Fas splicing. In addition, our large-scale RNA-seq analysis using GTEX and TCGA indicated that while SRSF6 expression was correlated with Fas expression in normal tissues, the correlation was disrupted in tumors. Our results suggest a novel regulatory mechanisms of Fas AS.

Abstract

Alternative splicing (AS) is a procedure during gene expression that allows the production of multiple mRNAs from a single gene, leading to a larger number of proteins with various functions. The alternative splicing (AS) of Fas (Apo-1/CD95) pre-mRNA can generate membrane-bound or soluble isoforms with pro-apoptotic and anti-apoptotic functions. SRSF6, a member of the Serine/Arginine-rich protein family, plays essential roles in both constitutive and alternative splicing. Here, we identified SRSF6 as an important regulatory protein in Fas AS. The cassette exon inclusion of Fas was decreased by SRSF6-targeting shRNA treatment, but increased by SRSF6 overexpression. The deletion and substitution mutagenesis of the Fas minigene demonstrated that the UGCCAA sequence in the cassette exon of the Fas gene causes the functional disruption of SRSF6, indicating that these sequences are essential for SRSF6 function in Fas splicing. In addition, biotin-labeled RNA-pulldown and immunoblotting analysis showed that SRSF6 interacted with these RNA sequences. Mutagenesis in the splice-site strength alteration demonstrated that the 5′ splice-site, but not the 3′ splice-site, was required for the SRSF6 regulation of Fas pre-mRNA. In addition, a large-scale RNA-seq analysis using GTEX and TCGA indicated that while SRSF6 expression was correlated with Fas expression in normal tissues, the correlation was disrupted in tumors. Furthermore, high SRSF6 expression was linked to the high expression of pro-apoptotic and immune activation genes. Therefore, we identified a novel RNA target with 5′ splice-site dependence of SRSF6 in Fas pre-mRNA splicing, and a correlation between SRSF6 and Fas expression.

The Role of Extracellular Matrix Proteins in Breast Cancer

The extracellular matrix is a structure composed of many molecules, including fibrillar (types I, II, III, V, XI, XXIV, XXVII) and non-fibrillar collagens (mainly basement membrane collagens: types IV, VIII, X), non-collagenous glycoproteins (elastin, laminin, fibronectin, thrombospondin, tenascin, osteopontin, osteonectin, entactin, periostin) embedded in a gel of negatively charged water-retaining glycosaminoglycans (GAGs) such as non-sulfated hyaluronic acid (HA) and sulfated GAGs which are linked to a core protein to form proteoglycans (PGs). This highly dynamic molecular network provides critical biochemical and biomechanical cues that mediate the cell–cell and cell–matrix interactions, influence cell growth, migration and differentiation and serve as a reservoir of cytokines and growth factors’ action. The breakdown of normal ECM and its replacement with tumor ECM modulate the tumor microenvironment (TME) composition and is an essential part of tumorigenesis and metastasis, acting as key driver for malignant progression. Abnormal ECM also deregulate behavior of stromal cells as well as facilitating tumor-associated angiogenesis and inflammation. Thus, the tumor matrix modulates each of the classically defined hallmarks of cancer promoting the growth, survival and invasion of the cancer. Moreover, various ECM-derived components modulate the immune response affecting T cells, tumor-associated macrophages (TAM), dendritic cells and cancer-associated fibroblasts (CAF). This review article considers the role that extracellular matrix play in breast cancer. Determining the detailed connections between the ECM and cellular processes has helped to identify novel disease markers and therapeutic targets.

PHF5A promotes colorectal cancerprogression by alternative splicing of TEAD2

Dysregulated alternative splicing (AS) plays critical roles in driving cancer progression, and the underlying mechanisms remain largely unknown. Here, we demonstrated that PHF5A, a component of U2 small nuclear ribonucleoproteins, was frequently upregulated in colorectal cancer (CRC) samples and associated with poor prognosis. PHF5A promoted proliferation and metastasis of CRC cells in vitro and in vivo. Transcriptomic analysis identified PHF5A-regulated AS targets and pathways. Particularly, PHF5A induced TEAD2 exon 2 inclusion to activate YAP signaling, and interference of TEAD2-L partially reversed the PHF5A-mediated tumor progression. Pharmacological inhibition of PHF5A using pladienolide B had potent antitumor activity. Collectively, these data revealed the oncogenic role of PHF5A in CRC through regulating AS and established PHF5A as potential therapeutic target.

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

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

Global Alternative Splicing Defects in Human Breast Cancer Cells

Simple Summary

Aberrant alternative splicing (AS) regulation plays a pivotal role in breast cancer development, progression, and resistance to therapeutical interventions. Indeed, cancer cells can adapt their own transcriptome by changing different AS programs, thus generating cancer-specific AS isoforms involved in every hallmark of cancer. Here, we investigated global AS errors occurring in human breast cancer cells by using RNA-mediated oligonucleotide annealing, selection, and ligation coupled with next-generation sequencing. Our results identified several dysregulated AS events potentially relevant for breast cancer-related biological processes and that provide a better comprehension of the molecular mechanisms that orchestrate the malignant transformation.

Abstract

Breast cancer is the most frequently occurred cancer type and the second cause of death in women worldwide. Alternative splicing (AS) is the process that generates more than one mRNA isoform from a single gene, and it plays a major role in expanding the human protein diversity. Aberrant AS contributes to breast cancer metastasis and resistance to chemotherapeutic interventions. Therefore, identifying cancer-specific isoforms is the prerequisite for therapeutic interventions intended to correct aberrantly expressed AS events. Here, we performed RNA-mediated oligonucleotide annealing, selection, and ligation coupled with next-generation sequencing (RASL-seq) in breast cancer cells, to identify global breast cancer-specific AS defects. By RT-PCR validation, we demonstrate the high accuracy of RASL-seq results. In addition, we analyzed identified AS events using the Cancer Genome Atlas (TCGA) database in a large number of non-pathological and breast tumor specimens and validated them in normal and breast cancer samples. Interestingly, aberrantly regulated AS cassette exons in cancer tissues do not encode for known functional domains but instead encode for amino acids constituting regions of intrinsically disordered protein portions characterized by high flexibility and prone to be subjected to post-translational modifications. Collectively, our results reveal novel AS errors occurring in human breast cancer, potentially affecting breast cancer-related biological processes.

Microarray analysis reveals ONC201 mediated differential mechanisms of CHOP gene regulation in metastatic and nonmetastatic colorectal cancer cells

The imipramine ONC201 has antiproliferative effects in several cancer cell types and activates integrated stress response pathway associated with the induction of Damage Inducible Transcript 3 (DDIT3, also known as C/EBP homologous protein or CHOP). We investigated the signaling pathways through which ONC201/CHOP crosstalk is regulated in ONC201-treated nonmetastatic and metastatic cancer cell lines (Dukes' type B colorectal adenocarcinoma nonmetastatic SW480 and metastatic LS-174T cells, respectively). Cell proliferation and apoptosis were evaluated by MTT assays and flow cytometry, gene expression was assessed by Affymetrix microarray, signaling pathway perturbations were assessed in silico, and key regulatory proteins were validated by Western blotting. Unlike LS-174T cells, SW480 cells were resistant to ONC201 treatment; Gene Ontology analysis of differentially expressed genes showed that cellular responsiveness to ONC201 treatment also differed substantially. In both ONC201-treated cell lines, CHOP expression was upregulated; however, its upstream regulatory mechanisms were perturbed. Although, PERK, ATF6 and IRE1 ER-stress pathways upregulated CHOP in both cell types, the Bak/Bax pathway regulated CHOP only LS-174T cells. Additionally, CHOP RNA splicing profiles varied between cell lines; these were further modified by ONC201 treatment. In conclusion, we delineated the signaling mechanisms by which CHOP expression is regulated in ONC201-treated non-metastatic and metastatic colorectal cell lines. The observed differences could be related to cellular plasticity and metabolic reprogramming, nevertheless, detailed mechanistic studies are required for further validations.

Prediction of disease-associated functional variants in noncoding regions through a comprehensive analysis by integrating datasets and features

One of the greatest challenges in human genetics is deciphering the link between functional variants in noncoding sequences and the pathophysiology of complex diseases. To address this issue, many methods have been developed to sort functional single-nucleotide variants (SNVs) for neutral SNVs in noncoding regions. In this study, we integrated well-established features and commonly used datasets and merged them into large-scale datasets based on a random forest model, which yielded promising performance and outperformed some cutting-edge approaches. Our analyses of feature importance and data coverage also provide certain clues for future research in enhancing the prediction of functional noncoding SNVs.

Widespread Alternative Splicing Changes in Metastatic Breast Cancer Cells

Aberrant alternative splicing (AS) is a hallmark of cancer and a potential target for novel anti-cancer therapeutics. Breast cancer-associated AS events are known to be linked to disease progression, metastasis, and survival of breast cancer patients. To identify altered AS programs occurring in metastatic breast cancer, we perform a global analysis of AS events by using RNA-mediated oligonucleotide annealing, selection, and ligation coupled with next-generation sequencing (RASL-seq). We demonstrate that, relative to low-metastatic, high-metastatic breast cancer cells show different AS choices in genes related to cancer progression. Supporting a global reshape of cancer-related splicing profiles in metastatic breast cancer we found an enrichment of RNA-binding motifs recognized by several splicing regulators, which have aberrant expression levels or activity during breast cancer progression, including SRSF1. Among SRSF1-regulated targets we found DCUN1D5, a gene for which skipping of exon 4 in its pre-mRNA introduces a premature termination codon (PTC), thus generating an unstable transcript degraded by nonsense-mediated mRNA decay (NMD). Significantly, distinct breast cancer subtypes show different DCUN1D5 isoform ratios with metastatic breast cancer expressing the highest level of the NMD-insensitive DCUN1D5 mRNA, thus showing high DCUN1D5 expression levels, which are ultimately associated with poor overall and relapse-free survival in breast cancer patients. Collectively, our results reveal global AS features of metastatic breast tumors, which open new possibilities for the treatment of these aggressive tumor types.

lncRNA TRMP-S directs dual mechanisms to regulate p27-mediated cellular senescence

Long noncoding RNAs (lncRNAs) undergo extensive alternative splicing, but little is known about isoform functions. A prior investigation of lncRNA RP11-369C8.1 reported that its splice variant TRMP suppressed p27 translation through PTBP1. Here we characterize a second major splice variant, TRMP-S (short variant), whose enforced loss promotes cancer cell-cycle arrest and p27-dependent entry into cellular senescence. Remarkably, despite sharing a single common exon with TRMP, TRMP-S restrains p27 expression through distinct mechanisms. First, TRMP-S stabilizes UHRF1 protein levels, an epigenetic inhibitor of p27, by promoting interactions between UHRF1 and its deubiquitinating enzyme USP7. Alternatively, binding interactions between TRMP-S and FUBP3 prevent p53 mRNA interactions with RPL26 ribosomal protein, the latter essential for promoting p53 translation with ensuing suppression of p53 translation limiting p27 expression. Significantly, as TRMP-S is itself transactivated by p53, this identifies negative feedback regulation between p53 and TRMP-S. Different splicing variants of the RP11-369C8.1 gene thereby exert distinct roles that converge on the homeostatic control of p27 expression, providing an important precedent for understanding the actions of alternatively spliced lncRNAs.

Alternative splicing in Alzheimer's disease

Alzheimer's disease (AD) is the most frequent neurodegenerative disorder in the elderly, occurring in approximately 20% of people older than 80. The molecular causes of AD are still poorly understood. However, recent studies have shown that Alternative Splicing (AS) is involved in the gene expression reprogramming associated with the functional changes observed in AD patients. In particular, mutations in cis-acting regulatory sequences as well as alterations in the activity and sub-cellular localization of trans-acting splicing factors and components of the spliceosome machinery are associated with splicing abnormalities in AD tissues, which may influence the onset and progression of the disease. In this review, we discuss the current molecular understanding of how alterations in the AS process contribute to AD pathogenesis. Finally, recent therapeutic approaches targeting aberrant AS regulation in AD are also reviewed.

SRSF9 Regulates Cassette Exon Splicing of Caspase-2 by Interacting with Its Downstream Exon

Alternative splicing (AS) is an important posttranscriptional regulatory process. Damaged or unnecessary cells need to be removed though apoptosis to maintain physiological processes. Caspase-2 pre-mRNA produces pro-apoptotic long mRNA and anti-apoptotic short mRNA isoforms through AS. How AS of Caspase-2 is regulated remains unclear. In the present study, we identified a novel regulatory protein SRSF9 for AS of Caspase-2 cassette exon 9. Knock-down (KD) of SRSF9 increased inclusion of cassette exon and on the other hand, overexpression of SRSF9 decreased inclusion of this exon. Deletion mutagenesis demonstrated that exon 9, parts of intron 9, exon 8 and exon 10 were not required for the role of SRSF9 in Caspase-2 AS. However, deletion and substitution mutation analysis revealed that AGGAG sequence located at exon 10 provided functional target for SRSF9. In addition, RNA-pulldown mediated immunoblotting analysis showed that SRSF9 interacted with this sequence. Gene ontology analysis of RNA-seq from SRSF9 KD cells demonstrates that SRSF9 could regulate AS of a subset of apoptosis related genes. Collectively, our results reveal a basis for regulation of Caspase-2 AS.

Hypoxia-induced alternative splicing in human diseases: the pledge, the turn, and the prestige

Maintenance of oxygen homeostasis is an indispensable criterion for the existence of multicellular life-forms. Disruption of this homeostasis due to inadequate oxygenation of the respiring tissues leads to pathological hypoxia, which acts as a significant stressor in several pathophysiological conditions including cancer, cardiovascular defects, bacterial infections, and neurological disorders. Consequently, the hypoxic tissues develop necessary adaptations both at the tissue and cellular level. The cellular adaptations involve a dramatic alteration in gene expression, post-transcriptional and post-translational modification of gene products, bioenergetics, and metabolism. Among the key responses to oxygen-deprivation is the skewing of cellular alternative splicing program. Herein, we discuss the current concepts of oxygen tension-dependent alternative splicing relevant to various pathophysiological conditions. Following a brief description of cellular response to hypoxia and the pre-mRNA splicing mechanism, we outline the impressive number of hypoxia-elicited alternative splicing events associated with maladies like cancer, cardiovascular diseases, and neurological disorders. Furthermore, we discuss how manipulation of hypoxia-induced alternative splicing may pose promising strategies for novel translational diagnosis and therapeutic interventions.

Splicing factors control triple-negative breast cancer cell mitosis through SUN2 interaction and sororin intron retention

BACKGROUND

Triple negative breast cancer (TNBC) is an aggressive subtype of breast cancer with limited therapeutic opportunities. Recently, splicing factors have gained attention as potential targets for cancer treatment. Here we systematically evaluated the role of RNA splicing factors in TNBC cell proliferation.

METHODS

In this study, we performed an RNAi screen targeting 244 individual splicing factors to systematically evaluate their role in TNBC cell proliferation. For top candidates, mechanistic insight was gained using amongst others western blot, PCR, FACS, molecular imaging and cloning. Pulldown followed by mass spectrometry were used to determine protein-protein interactions and patient-derived RNA sequencing data was used relate splicing factor expression levels to proliferation markers.

RESULTS

We identified nine splicing factors, including SNRPD2, SNRPD3 and NHP2L1, of which depletion inhibited proliferation in two TNBC cell lines by deregulation of sister chromatid cohesion (SCC) via increased sororin intron 1 retention and down-regulation of SMC1, MAU2 and ESPL1. Protein-protein interaction analysis of SNRPD2, SNRPD3 and NHP2L1 identified that seven out of the nine identified splicing factors belong to the same spliceosome complex including novel component SUN2 that was also critical for efficient sororin splicing. Finally, sororin transcript levels are highly correlated to various proliferation markers in BC patients.

CONCLUSION

We systematically determined splicing factors that control proliferation of breast cancer cells through a mechanism that involves effective sororin splicing and thereby appropriate sister chromatid cohesion. Moreover, we identified SUN2 as an important new spliceosome complex interacting protein that is critical in this process. We anticipate that deregulating sororin levels through targeting of the relevant splicing factors might be a potential strategy to treat TNBC.

U2AF65-Dependent SF3B1 Function in SMN Alternative Splicing

Splicing factor 3b subunit 1 (SF3B1) is an essential protein in spliceosomes and mutated frequently in many cancers. While roles of SF3B1 in single intron splicing and roles of its cancer-linked mutant in aberrant splicing have been identified to some extent, regulatory functions of wild-type SF3B1 in alternative splicing (AS) are not well-understood yet. Here, we applied RNA sequencing (RNA-seq) to analyze genome-wide AS in SF3B1 knockdown (KD) cells and to identify a large number of skipped exons (SEs), with a considerable number of alternative 5′ splice-site selection, alternative 3′ splice-site selection, mutually exclusive exons (MXE), and retention of introns (RI). Among altered SEs by SF3B1 KD, survival motor neuron 2 (SMN2) pre-mRNA exon 7 splicing was a regulatory target of SF3B1. RT-PCR analysis of SMN exon 7 splicing in SF3B1 KD or overexpressed HCT116, SH-SY5Y, HEK293T, and spinal muscular atrophy (SMA) patient cells validated the results. A deletion mutation demonstrated that the U2 snRNP auxiliary factor 65 kDa (U2AF65) interaction domain of SF3B1 was required for its function in SMN exon 7 splicing. In addition, mutations to lower the score of the polypyrimidine tract (PPT) of exon 7, resulting in lower affinity for U2AF65, were not able to support SF3B1 function, suggesting the importance of U2AF65 in SF3B1 function. Furthermore, the PPT of exon 7 with higher affinity to U2AF65 than exon 8 showed significantly stronger interactions with SF3B1. Collectively, our results revealed SF3B1 function in SMN alternative splicing.

Alternative Splicing: Expanding the Landscape of Cancer Biomarkers and Therapeutics

Alternative splicing (AS) is a critical post-transcriptional regulatory mechanism used by more than 95% of transcribed human genes and responsible for structural transcript variation and proteome diversity. In the past decade, genome-wide transcriptome sequencing has revealed that AS is tightly regulated in a tissue- and developmental stage-specific manner, and also frequently dysregulated in multiple human cancer types. It is currently recognized that splicing defects, including genetic alterations in the spliced gene, altered expression of both core components or regulators of the precursor messenger RNA (pre-mRNA) splicing machinery, or both, are major drivers of tumorigenesis. Hence, in this review we provide an overview of our current understanding of splicing alterations in cancer, and emphasize the need to further explore the cancer-specific splicing programs in order to obtain new insights in oncology. Furthermore, we also discuss the recent advances in the identification of dysregulated splicing signatures on a genome-wide scale and their potential use as biomarkers. Finally, we highlight the therapeutic opportunities arising from dysregulated splicing and summarize the current approaches to therapeutically target AS in cancer.

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.

Alternative splicing of leptin receptor overlapping transcript in osteosarcoma

IMPACT STATEMENT

Osteosarcoma (OS, also known as osteogenic sarcoma) is the most common primary malignancy of bone in children and adolescents. The molecular mechanisms of OS are extremely complicated and its molecular mediators remain to be elucidated. We sequenced total RNA from 18 OS bone samples (paired normal-tumor biopsies). We found statistically significant (FDR <0.05) 26 differentially expressed transcript variants of LEPROT gene with different expressions in normal and tumor samples. These findings contribute to the understanding of molecular mechanisms of OS development and provide encouragement to pursue further research.

Overexpression of Extradomain-B Fibronectin is Associated with Invasion of Breast Cancer Cells

Breast tumor heterogeneity is a major impediment to oncotherapy. Cancer cells undergo rapid clonal evolution, thereby acquiring significant growth and invasive advantages. The absence of specific markers of these high-risk populations precludes efficient therapeutic and diagnostic management of the disease. Given the critical function of tumor microenvironment in the oncogenic circuitry, we sought to determine the expression profile of the extracellular matrix oncoprotein, extradomain-B fibronectin (EDB-FN) in invasive breast cancer. Analyses of TCGA/GTEx databases and immunostaining of clinical samples found a significant overexpression of EDB-FN in breast tumors, which correlated with poor overall survival. Significant upregulation of EDB-FN was observed in invasive cell populations generated from relatively less invasive MCF7 and MDA-MB-468 cells by long-term TGF-β treatment and acquired chemoresistance. Treatment of the invasive cell populations with an AKT inhibitor (MK2206-HCl) reduced their invasive potential, with a concomitant decrease in their EDB-FN expression, partly through the phosphoAKT-SRp55 pathway. EDB-FN downregulation, with direct RNAi of EDB-FN or indirectly through RNAi of SRp55, also resulted in reduced motility of the invasive cell populations, validating the correlation between EDB-FN expression and invasion of breast cancer cells. These data establish EDB-FN as a promising molecular marker for non-invasive therapeutic surveillance of aggressive breast cancer.

HnRNP A1 - mediated alternative splicing of CCDC50 contributes to cancer progression of clear cell renal cell carcinoma via ZNF395

BACKGROUND

Aberrant alternative splicing events play critical roles in carcinogenesis and progression of many cancers, while sparse studies regarding to alternative splicing are available for clear cell renal cell carcinoma (ccRCC). We identified that alternative splicing of coiled-coil domain containing 50 (CCDC50) was dysregulated in ccRCC, whereas the clinical significance of this splicing event and its splicing regulation mechanisms were still elusive.

METHODS

Bioinformatic algorithm was utilized to identify significant exon skipping events in ccRCC via exon sequencing data from The Cancer Genome Atlas. Semi-quantitative real-time polymerase chain reaction and western blot were used to validate the aberrant expression of different transcripts in renal cancer tissues, cell lines and corresponding noncancerous controls. Short hairpin RNA targeting CCDC50 and overexpressing plasmids for each transcript were introduced into ccRCC cell lines, followed by a series of in vitro and in vivo functional experiments. Moreover, a panel of splicing factors were identified and their roles on splicing regulation of CCDC50 precursor mRNA (pre-mRNA) were studied. Furthermore, RNAseq data were analyzed to elucidate downstream molecules of CCDC50. Two-way analysis of variance and unpaired Student t test were used in statistical analysis.

RESULTS

Pre-mRNA of CCDC50 generated two transcripts, full-length transcript (CCDC50-FL) and truncated transcript (CCDC50-S) with exon 6 skipped. CCDC50-S was overexpressed in ccRCC tissues and cell lines compared to noncancerous counterparts, but CCDC50-FL was only detected in noncancerous tissues and normal renal epithelial cells. Higher percent spliced-in index was associated with better survival in ccRCC patients. In vitro and in vivo functional experiments indicated that CCDC50-S transcript promoted the proliferation, migration, invasion and tumorigenesis of ccRCC, while CCDC50-FL exerted opposite tumor suppressive functions. Besides, we identified that heterogeneous nuclear ribonucleoprotein A1 (HnRNP A1) could promote the skipping of exon 6, which resulted in higher portion of CCDC50-S and oncogenic transformation. Moreover, zinc finger protein 395 (ZNF395) was identified as a downstream protein of CCDC50-S, and the interaction initiated oncogenic pathways which were involved in ccRCC progression.

CONCLUSIONS

Aberrant alternative splicing of CCDC50 is regulated by HnRNP A1 in ccRCC. This splicing event contributes to cancer progression through the downstream pathway involving ZNF395.

Insight into the Development of PET Radiopharmaceuticals for Oncology

While the development of positron emission tomography (PET) radiopharmaceuticals closely follows that of traditional drug development, there are several key considerations in the chemical and radiochemical synthesis, preclinical assessment, and clinical translation of PET radiotracers. As such, we outline the fundamentals of radiotracer design, with respect to the selection of an appropriate pharmacophore. These concepts will be reinforced by exemplary cases of PET radiotracer development, both with respect to their preclinical and clinical evaluation. We also provide a guideline for the proper selection of a radionuclide and the appropriate labeling strategy to access a tracer with optimal imaging qualities. Finally, we summarize the methodology of their evaluation in in vitro and animal models and the road to clinical translation. This review is intended to be a primer for newcomers to the field and give insight into the workflow of developing radiopharmaceuticals.

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.

RBM4 regulates M1 macrophages polarization through targeting STAT1-mediated glycolysis

M1/M2 macrophages polarization play important roles in regulating tissue homeostasis. Recently, RNA-binding motif 4 (RBM4) has been reported to modulate the proliferation and expression of inflammatory factors in HeLa cells. However, whether RBM4 is involved in regulating macrophage polarization and inflammatory factor expression are still unknown. In this study, RAW264.7, a mouse macrophage cell line, were stimulated with interferon γ (IFN-γ) or interleukin-4 (IL-4) to induce M1/M2 macrophages polarization. We found that IFN-γ, but not IL-4, stimulation decreased RBM4 expression in macrophages, and RBM4 overexpression inhibits IFN-γ-induced M1 macrophage polarization. Furthermore, RNA-Sequencing, protein immunoprecipitation accompanied with mass spectrometry, and extracellular acidification rate analysis showed that RBM4 suppresses IFN-γ-induced M1 macrophage polarization though inhibiting glycolysis. Moreover, RBM4 knockdown promoted IFN-γ-induced signal transducer and activator of transcription 1 (STAT1) activation via increasing STAT1 mRNA stability, leading to the increase of glycolysis-related gene transcripts regulated by STAT1. Finally, we find that RBM4 interacts with YTH N6-methyladenosine RNA binding protein 2 (YTHDF2) to degrade m6A modified STAT1 mRNA, thereby regulating glycolysis and M1 macrophage polarization. Collectively, the current study firstly reports that RBM4 regulates M1 macrophages polarization through targeting STAT1-mediated glycolysis and shows that RBM4 is a possible candidate for regulating macrophage M1 polarization and inflammatory responses.

Differential splicing analysis based on isoforms expression with NBSplice

Alternative splicing alterations have been widely related to several human diseases revealing the importance of their study for the success of translational medicine. Differential splicing (DS) occurrence has been mainly analyzed through exon-based approaches over RNA-seq data. Although these strategies allow identifying differentially spliced genes, they ignore the identity of the affected gene isoforms which is crucial to understand the underlying pathological processes behind alternative splicing changes. Moreover, despite several isoform quantification tools for RNA-seq data have been recently developed, DS tools have not taken advantage of them. Here, the NBSplice R package for differential splicing analysis by means of isoform expression data is presented. It estimates differences on relative expressions of gene transcripts between experimental conditions to infer changes in gene alternative splicing patterns. The developed tool was evaluated using a synthetic RNA-seq dataset with controlled differential splicing. NBSplice accurately predicted DS occurrence, outperforming current methods in terms of accuracy, sensitivity, F-score, and false discovery rate control. The usefulness of our development was demonstrated by the analysis of a real cancer dataset, revealing new differentially spliced genes that could be studied pursuing new colorectal cancer biomarkers discovery.

Gene Expression Profiles Controlled by the Alternative Splicing Factor Nova2 in Endothelial Cells

Alternative splicing (AS) plays an important role in expanding the complexity of the human genome through the production of specialized proteins regulating organ development and physiological functions, as well as contributing to several pathological conditions. How AS programs impact on the signaling pathways controlling endothelial cell (EC) functions and vascular development is largely unknown. Here we identified, through RNA-seq, changes in mRNA steady-state levels in ECs caused by the neuro-oncological ventral antigen 2 (Nova2), a key AS regulator of the vascular morphogenesis. Bioinformatics analyses identified significant enrichment for genes regulated by peroxisome proliferator-activated receptor-gamma (Ppar-γ) and E2F1 transcription factors. We also showed that Nova2 in ECs controlled the AS profiles of Ppar-γ and E2F dimerization partner 2 (Tfdp2), thus generating different protein isoforms with distinct function (Ppar-γ) or subcellular localization (Tfdp2). Collectively, our results supported a mechanism whereby Nova2 integrated splicing decisions in order to regulate Ppar-γ and E2F1 activities. Our data added a layer to the sequential series of events controlled by Nova2 in ECs to orchestrate vascular biology.

AIDE: annotation-assisted isoform discovery with high precision

Genome-wide accurate identification and quantification of full-length mRNA isoforms is crucial for investigating transcriptional and posttranscriptional regulatory mechanisms of biological phenomena. Despite continuing efforts in developing effective computational tools to identify or assemble full-length mRNA isoforms from second-generation RNA-seq data, it remains a challenge to accurately identify mRNA isoforms from short sequence reads owing to the substantial information loss in RNA-seq experiments. Here, we introduce a novel statistical method, annotation-assisted isoform discovery (AIDE), the first approach that directly controls false isoform discoveries by implementing the testing-based model selection principle. Solving the isoform discovery problem in a stepwise and conservative manner, AIDE prioritizes the annotated isoforms and precisely identifies novel isoforms whose addition significantly improves the explanation of observed RNA-seq reads. We evaluate the performance of AIDE based on multiple simulated and real RNA-seq data sets followed by PCR-Sanger sequencing validation. Our results show that AIDE effectively leverages the annotation information to compensate the information loss owing to short read lengths. AIDE achieves the highest precision in isoform discovery and the lowest error rates in isoform abundance estimation, compared with three state-of-the-art methods Cufflinks, SLIDE, and StringTie. As a robust bioinformatics tool for transcriptome analysis, AIDE enables researchers to discover novel transcripts with high confidence.

Splicing regulatory factors in breast cancer hallmarks and disease progression

By regulating transcript isoform expression levels, alternative splicing provides an additional layer of protein control. Recent studies show evidence that cancer cells use different splicing events to fulfill their requirements in order to develop, progress and metastasize. However, there has been less attention for the role of the complex catalyzing the complicated multistep splicing reaction: the spliceosome. The spliceosome consists of multiple sub-complexes in total comprising 244 proteins or splice factors and 5 associated RNA molecules. Here we discuss the role of splice factors in the oncogenic processes tumors cells need to fulfill their oncogenic properties (the so-called the hallmarks of cancer). Despite the fact that splice factors have been investigated only recently, they seem to play a prominent role in already five hallmarks of cancer: angiogenesis, resisting cell death, sustaining proliferation, deregulating cellular energetics and invasion and metastasis formation by affecting major signaling pathways such as epithelial-to-mesenchymal transition, the Warburg effect, DNA damage response and hormone receptor dependent proliferation. Moreover, we could relate expression of representative genes of four other hallmarks (enabling replicative mortality, genomic instability, avoiding immune destruction and evading growth suppression) to splice factor levels in human breast cancer tumors, suggesting that also these hallmarks could be regulated by splice factors. Since many splice factors are involved in multiple hallmarks of cancer, inhibiting splice factors might provide a new layer of oncogenic control and a powerful method to combat breast cancer progression.

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

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

A bioinformatic analysis identifies circadian expression of splicing factors and time-dependent alternative splicing events in the HD-MY-Z cell line

The circadian clock regulates key cellular processes and its dysregulation is associated to several pathologies including cancer. Although the transcriptional regulation of gene expression by the clock machinery is well described, the role of the clock in the regulation of post-transcriptional processes, including splicing, remains poorly understood. In the present work, we investigated the putative interplay between the circadian clock and splicing in a cancer context. For this, we applied a computational pipeline to identify oscillating genes and alternatively spliced transcripts in time-course high-throughput data sets from normal cells and tissues, and cancer cell lines. We investigated the temporal phenotype of clock-controlled genes and splicing factors, and evaluated their impact in alternative splice patterns in the Hodgkin Lymphoma cell line HD-MY-Z. Our data points to a connection between clock-controlled genes and splicing factors, which correlates with temporal alternative splicing in several genes in the HD-MY-Z cell line. These include the genes DPYD, SS18, VIPR1 and IRF4, involved in metabolism, cell cycle, apoptosis and proliferation. Our results highlight a role for the clock as a temporal regulator of alternative splicing, which may impact malignancy in this cellular model.

Events of alternative splicing in head and neck cancer via RNA sequencing - an update

Background

Alternative splicing (AS) is a regulatory mechanism used to create many forms of mature messengers RNAs (mRNAs) from the same gene. Sequencing of RNA (RNA-Seq) is an advanced technology, which has been utilized by different studies to find AS mechanisms in head and neck cancer (HNC). Hitherto, there is no available review that could inform us of the major findings from these studies. Hence, we aim to perform a systematic literature search following PRISMA guidelines to study AS events in HNC identified through RNA-Seq studies.

Results

A total of five records were identified that utilized RNA-Seq data for identifying AS events in HNC. Five software was used in these records to identify AS events. Two genes influenced by AS i.e. MLL3 and RPS9 were found to be common in 4 out of 5 records. Likewise, 38 genes were identified to be similar in at least 3 records.

Conclusions

Alternative splicing in HNC is a multifaceted regulatory mechanism of gene expression. It can be studied via RNA-Seq using different bioinformatics tools. Genes MLL3, as well as RPS9, were repeatedly found to be associated with HNC, however needs further functional validation.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5794-y) contains supplementary material, which is available to authorized users.

RNA-Binding Motif 4 (RBM4) Suppresses Tumor Growth and Metastasis in Human Gastric Cancer

BACKGROUND Dysregulation of the splicing activator, RNA-binding motif 4 (RBM4), has recently been reported to be involved in the progression of several cancers. However, the mechanisms that underpin the activity of RBM4 in gastric cancer (GC) remain unknown. The purpose of our study was to explore how RBM4 affects the biological behavior of GC through in vivo and in vitro experiments. MATERIAL AND METHODS Western blot and flow cytometry analyses were used to investigate the RBM4 protein levels in normal gastric epithelial cells and 5 types of GC cells. Cell Counting Kit-8 assay, flow cytometry analysis, wound-healing, and migration and invasion assays were evaluated in vitro in BGC823 and MGC803 GC cells. A xenograft tumor model was used to assess whether RBM4 inhibits GC growth in vivo. Mitogen-activated protein kinase (MAPK) protein levels were determined using western blot analyses. RESULTS Our study revealed that RBM4 protein was downregulated in GC cells. Re-expression of RBM4 inhibited the proliferation, migration, and invasion of GC cells, while promoting apoptosis. Thus, the overexpression of RBM4 can inhibit tumor growth in GC mouse models. We also report that RBM4 was involved in the activation of MAPK-dependent signaling pathways in human GC. CONCLUSIONS It is hoped that these findings will improve our understanding of GC pathogenesis while also helping us to explore the feasibility of RBM4-targeted therapy for GC treatment.

CRKL regulates alternative splicing of cancer-related genes in cervical cancer samples and HeLa cell

BACKGROUND

Aberrant spliced isoforms are specifically associated with cancer progression and metastasis. The cytoplasmic adaptor CRKL (v-crk avian sarcoma virus CT10 oncogene homolog-like) is a CRK like proto-oncogene, which encodes a SH2 and SH3 (src homology) domain-containing adaptor protein. CRKL is tightly linked to leukemia via its binding partners BCR-ABL and TEL-ABL, upregulated in multiple types of human cancers, and induce cancer cell proliferation and invasion. However, it remains unclear whether signaling adaptors such as CRKL could regulate alternative splicing.

METHODS

We analyzed the expression level of CRKL in 305 cervical cancer tissue samples available in TCGA database, and then selected two groups of cancer samples with CRKL differentially expressed to analyzed potential CRKL-regulated alternative splicing events (ASEs). CRKL was knocked down by shRNA to further study CRKL-regulated alternative splicing and the activity of SR protein kinases in HeLa cells using RNA-Seq and Western blot techniques. We validated 43 CRKL-regulated ASEs detected by RNA-seq in HeLa cells, using RT-qPCR analysis of HeLa cell samples and using RNA-seq data of the two group of clinical cervical samples.

RESULTS

The expression of CRKL was mostly up-regulated in stage I cervical cancer samples. Knock-down of CRKL led to a reduced cell proliferation. CRKL-regulated alternative splicing of a large number of genes were enriched in cancer-related functional pathways, among which DNA repair and G2/M mitotic cell cycle, GnRH signaling were shared among the top 10 enriched GO terms and KEGG pathways by results from clinical samples and HeLa cell model. We showed that CRKL-regulated ASEs revealed by computational analysis using ABLas software in HeLa cell were highly validated by RT-qPCR, and also validated by cervical cancer clinical samples.

CONCLUSIONS

This is the first report of CRKL-regulation of the alternative splicing of a number of genes critical in tumorigenesis and cancer progression, which is consistent with CRKL reported role as a signaling adaptor and a kinase. Our results underline that the signaling adaptor CRKL might integrate the external and intrinsic cellular signals and coordinate the dynamic activation of cellular signaling pathways including alternative splicing regulation.

Pan-Cancer Analysis of TCGA Data Revealed Promising Reference Genes for qPCR Normalization

Quantitative PCR (qPCR) remains the most widely used technique for gene expression evaluation. Obtaining reliable data using this method requires reference genes (RGs) with stable mRNA level under experimental conditions. This issue is especially crucial in cancer studies because each tumor has a unique molecular portrait. The Cancer Genome Atlas (TCGA) project provides RNA-Seq data for thousands of samples corresponding to dozens of cancers and presents the basis for assessment of the suitability of genes as reference ones for qPCR data normalization. Using TCGA RNA-Seq data and previously developed CrossHub tool, we evaluated mRNA level of 32 traditionally used RGs in 12 cancer types, including those of lung, breast, prostate, kidney, and colon. We developed an 11-component scoring system for the assessment of gene expression stability. Among the 32 genes, PUM1 was one of the most stably expressed in the majority of examined cancers, whereas GAPDH, which is widely used as a RG, showed significant mRNA level alterations in more than a half of cases. For each of 12 cancer types, we suggested a pair of genes that are the most suitable for use as reference ones. These genes are characterized by high expression stability and absence of correlation between their mRNA levels. Next, the scoring system was expanded with several features of a gene: mutation rate, number of transcript isoforms and pseudogenes, participation in cancer-related processes on the basis of Gene Ontology, and mentions in PubMed-indexed articles. All the genes covered by RNA-Seq data in TCGA were analyzed using the expanded scoring system that allowed us to reveal novel promising RGs for each examined cancer type and identify several "universal" pan-cancer RG candidates, including SF3A1, CIAO1, and SFRS4. The choice of RGs is the basis for precise gene expression evaluation by qPCR. Here, we suggested optimal pairs of traditionally used RGs for 12 cancer types and identified novel promising RGs that demonstrate high expression stability and other features of reliable and convenient RGs (high expression level, low mutation rate, non-involvement in cancer-related processes, single transcript isoform, and absence of pseudogenes).

Interplay between post-translational cyclooxygenase-2 modifications and the metabolic and proteomic profile in a colorectal cancer cohort

BACKGROUND

Colorectal cancer (CRC) is the second most common cause of cancer death worldwide. It is broadly described that cyclooxygenase-2 (COX-2) is mainly overexpressed in CRC but less is known regarding post-translational modifications of this enzyme that may regulate its activity, intracellular localization and stability. Since metabolic and proteomic profile analysis is essential for cancer prognosis and diagnosis, our hypothesis is that the analysis of correlations between these specific parameters and COX-2 state in tumors of a high number of CRC patients could be useful for the understanding of the basis of this cancer in humans.

AIM

To analyze COX-2 regulation in colorectal cancer and to perform a detailed analysis of their metabolic and proteomic profile.

METHODS

Biopsies from both healthy and pathological colorectal tissues were taken under informed consent from patients during standard colonoscopy procedure in the University Hospital of Bellvitge (Barcelona, Spain) and Germans Trias i Pujol University Hospital (Campus Can Ruti) (Barcelona, Spain). Western blot analysis was used to determine COX-2 levels. Deglycosylation assays were performed in both cells and tumor samples incubating each sample with peptide N-glycosidase F (PNGase F). Prostaglandin E₂ (PGE₂) levels were determined using a specific ELISA. ¹H high resolution magic angle spinning (HRMAS) analysis was performed using a Bruker AVIII 500 MHz spectrometer and proteomic analysis was performed in a nano-liquid chromatography-tandem mass spectrometer (nano LC-MS/MS) using a QExactive HF orbitrap MS.

RESULTS

Our data show that COX-2 has a differential expression profile in tumor tissue of CRC patients vs the adjacent non-tumor area, which correspond to a glycosylated and less active state of the protein. This fact was associated to a lesser PGE₂ production in tumors. These results were corroborated in vitro performing deglycosylation assays in HT29 cell line where COX-2 protein profile was modified after PNGase F incubation, showing higher PGE₂ levels. Moreover, HRMAS analysis indicated that tumor tissue has altered metabolic features vs non-tumor counterparts, presenting increased levels of certain metabolites such as taurine and phosphocholine and lower levels of lactate. In proteomic experiments, we detected an enlarged number of proteins in tumors that are mainly implicated in basic biological functions like mitochondrial activity, DNA/RNA processing, vesicular trafficking, metabolism, cytoskeleton and splicing.

CONCLUSION

In our colorectal cancer cohort, tumor tissue presents a differential COX-2 expression pattern with lower enzymatic activity that can be related to an altered metabolic and proteomic profile.

Aberrant mRNA splicing generates oncogenic RNA isoforms and contributes to the development and progression of cholangiocarcinoma

Cholangiocarcinoma is a lethal biliary cancer, with an unclear molecular pathogenesis. Alternative splicing is a post-transcriptional modification that generates mature mRNAs, which are subsequently translated into proteins. Aberrant alternative splicing has been reported to serve a role in tumor initiation, maintenance and metastasis in several types of human cancer, including cholangiocarcinoma. In this review, the aberrant splicing of genes and the functional contributions of the spliced genes, in the carcinogenesis, progression and aggressiveness of cholangiocarcinoma are summarized. In addition, factors that influence this aberrant splicing that may be relevant as therapeutic targets or prognosis markers for cholangiocarcinoma are discussed.

FGFR1β is a driver isoform of FGFR1 alternative splicing in breast cancer cells

Abnormal FGFR1 alternative splicing is correlated with tumorigenicity and poor prognosis in several tumor types. We sought to determine the roles of FGFR1α and FGFR1β variants in breast cancer. TCGA samples and cell lines were analyzed for FGFR1α/FGFR1β expression. MCF-10A cells were used to overexpress these variants. Cell growth and transformation were assessed by SRB, colony formation, 3D-Matrigel, soft agar, cell motility assays. In TCGA, compared to FGFR1 non-amplified samples, FGFR1-amplified samples had significantly higher FGFR1α but not FGFR1β levels. FGFR1β expression levels and FGFR1β/FGFR1α ratio were higher in basal subtype samples than in ER-positive/luminal samples in both TCGA and breast cancer cell lines. Both FGFR1α and FGFR1β induced transformation of MCF-10A cells. However, only FGFR1β-expressing cells, not FGFR1α, enhanced cell growth and cell motility. Cells with higher FGFR1β levels and FGFR1β/FGFR1α ratio were more sensitive to FGFR inhibitor BGJ-398. Interestingly, in ER-negative cells, FGFR inhibitors decreased FGFR1β levels, likely by increasing expression of splicing repressor PTBP1. In ER-positive cells, estrogen treatment increased FGFR1β levels by decreasing PTBP1 expression, which was blocked by 4-OHT. Lastly, combination treatment with BGJ-398 and 4-OHT synergistically inhibited cell survival. These findings suggest that FGFR1 alternative FGFR1α/FGFR1β splicing plays an important role in breast cancer.

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.

Characterization of kinase gene expression and splicing profile in prostate cancer with RNA-Seq data

Background

Alternative splicing is a ubiquitous post-transcriptional regulation mechanism in most eukaryotic genes. Aberrant splicing isoforms and abnormal isoform ratios can contribute to cancer development. Kinase genes are key regulators of multiple cellular processes. Many kinases are found to be oncogenic and have been intensively investigated in the study of cancer and drugs. RNA-Seq provides a powerful technology for genome-wide study of alternative splicing in cancer besides the conventional gene expression profiling. But this potential has not been fully demonstrated yet.

Methods

We characterized the transcriptome profile of prostate cancer using RNA-Seq data from viewpoints of both differential expression and differential splicing, with an emphasis on kinase genes and their splicing variations. We built a pipeline to conduct differential expression and differential splicing analysis, followed by functional enrichment analysis. We performed kinase domain analysis to identify the functionally important candidate kinase gene in prostate cancer, and calculated the expression levels of isoforms to explore the function of isoform switching of kinase genes in prostate cancer.

Results

We identified distinct gene groups from differential expression and splicing analyses, which suggested that alternative splicing adds another level to gene expression regulation. Enriched GO terms of differentially expressed and spliced kinase genes were found to play different roles in regulation of cellular metabolism. Function analysis on differentially spliced kinase genes showed that differentially spliced exons of these genes are significantly enriched in protein kinase domains. Among them, we found that gene CDK5 has isoform switching between prostate cancer and benign tissues, which may affect cancer development by changing androgen receptor (AR) phosphorylation. The observation was validated in another RNA-Seq dataset of prostate cancer cell lines.

Conclusions

Our work characterized the expression and splicing profiles of kinase genes in prostate cancer and proposed a hypothetical model on isoform switching of CDK5 and AR phosphorylation in prostate cancer. These findings bring new understanding to the role of alternatively spliced kinases in prostate cancer and also demonstrate the use of RNA-Seq data in studying alternative splicing in cancer.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4925-1) contains supplementary material, which is available to authorized users.

SRSF2 directly inhibits intron splicing to suppresses cassette exon inclusion

SRSF2, a Serine-Arginine rich (SR) protein, is a splicing activator that mediates exon inclusion and exclusion events equally well. Here we show SRSF2 directly suppresses intron splicing to suppress cassette exon inclusion in SMN pre-mRNA. Through a serial mutagenesis, we demonstrate that a 10 nt RNA sequence surrounding the branch-point (BP), is important for SRSF2-mediated inhibition of cassette exon inclusion through directly interacting with SRSF2. We conclude that SRSF2 inhibits intron splicing to promote exon exclusion.

In1-ghrelin splicing variant is associated with reduced disease-free survival of breast cancer patients and increases malignancy of breast cancer cells lines

Ghrelin gene generates several variants that regulate multiple pathophysiological functions, including tumor-related processes. In1-ghrelin is a splicing variant that was previously shown to be overexpressed in breast cancer (BCa), where it correlated with proliferation markers; however, its possible association with clinical outcome of BCa patients and underlying mechanisms are still unknown. To address this issue, expression levels and clinical associations of In1-ghrelin were analyzed in a cohort of 117 BCa samples. Additionally, a battery of cellular and molecular assays was implemented using two BCa cell lines (MCF-7 and MDA-MB-231), wherein the role of In1-ghrelin on proliferation, migration, dedifferentiation and signaling pathways was explored. The results generated revealed that high expression of In1-ghrelin in BCa samples was associated with lymph node metastasis and reduced disease-free survival. Indeed, In1-ghrelin overexpression stimulated proliferation and migration in MCF-7 and MDA-MB-231 cells. Similar results were found by treating MDA-MB-231 and MCF-7 with In1-ghrelin-derived peptides. Conversely, In1-ghrelin silencing decreased proliferation and migration capacities of MDA-MB-231. Furthermore, In1-ghrelin (but not ghrelin) overexpression increased the capacity to form mammospheres in both cell lines. These effects could be associated with activation of MAPK-ERK, Jag1/Notch, Wnt/β-catenin and/or TGF-β1 pathways. Altogether, our data indicate that In1-ghrelin could play relevant functional roles in the regulation of BCa development and progression and may provide insights to identify novel biomarkers and new therapeutic approaches for this pathology.

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.

Collateral Damage Intended-Cancer-Associated Fibroblasts and Vasculature Are Potential Targets in Cancer Therapy

After oncogenic transformation, tumor cells rewire their metabolism to obtain sufficient energy and biochemical building blocks for cell proliferation, even under hypoxic conditions. Glucose and glutamine become their major limiting nutritional demands. Instead of being autonomous, tumor cells change their immediate environment not only by their metabolites but also by mediators, such as juxtacrine cell contacts, chemokines and other cytokines. Thus, the tumor cells shape their microenvironment as well as induce resident cells, such as fibroblasts and endothelial cells (ECs), to support them. Fibroblasts differentiate into cancer-associated fibroblasts (CAFs), which produce a qualitatively and quantitatively different extracellular matrix (ECM). By their contractile power, they exert tensile forces onto this ECM, leading to increased intratumoral pressure. Moreover, along with enhanced cross-linkage of the ECM components, CAFs thus stiffen the ECM. Attracted by tumor cell- and CAF-secreted vascular endothelial growth factor (VEGF), ECs sprout from pre-existing blood vessels during tumor-induced angiogenesis. Tumor vessels are distinct from EC-lined vessels, because tumor cells integrate into the endothelium or even mimic and replace it in vasculogenic mimicry (VM) vessels. Not only the VM vessels but also the characteristically malformed EC-lined tumor vessels are typical for tumor tissue and may represent promising targets in cancer therapy.

MCM2-regulated functional networks in lung cancer by multi-dimensional proteomic approach

DNA replication control is vital for maintaining genome stability and the cell cycle, perhaps most notably during cell division. Malignancies often exhibit defective minichromosome maintenance protein 2 (MCM2), a cancer proliferation biomarker that serves as a licensing factor in the initiation of DNA replication. MCM2 is also known to be one of the ATPase active sites that facilitates conformational changes and drives DNA unwinding at the origin of DNA replication. However, the biological networks of MCM2 in lung cancer cells via protein phosphorylation remain unmapped. The RNA-seq datasets from The Cancer Genome Atlas (TCGA) revealed that MCM2 overexpression is correlated with poor survival rate in lung cancer patients. To uncover MCM2-regulated functional networks in lung cancer, we performed multi-dimensional proteomic approach by integrating analysis of the phosphoproteome and proteome, and identified a total of 2361 phosphorylation sites on 753 phosphoproteins, and 4672 proteins. We found that the deregulation of MCM2 is involved in lung cancer cell proliferation, the cell cycle, and migration. Furthermore, HMGA1^S99 phosphorylation was found to be differentially expressed under MCM2 perturbation in opposite directions, and plays an important role in regulating lung cancer cell proliferation. This study therefore enhances our capacity to therapeutically target cancer-specific phosphoproteins.