Aberrant RNA Splicing in Cancer and Drug Resistance

Transcriptome complexity in intravascular NK/T-cell lymphoma

AIMS

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Differential roles of protease isoforms in the tumor microenvironment

Alternative splicing of precursor mRNA is a key mediator of gene expression regulation leading to greater diversity of the proteome in complex organisms. Systematic sequencing of the human genome and transcriptome has led to our understanding of how alternative splicing of critical genes leads to multiple pathological conditions such as cancer. For many years, proteases were known only for their roles as proteolytic enzymes, acting to regulate/process proteins associated with diverse cellular functions. However, the differential expression and altered function of various protease isoforms, such as (i) anti-apoptotic activities, (ii) mediating intercellular adhesion, and (iii) modifying the extracellular matrix, are evidence of their specific contribution towards shaping the tumor microenvironment. Revealing the alternative splicing of protease genes and characterization of their protein products/isoforms with distinct and opposing functions creates a platform to understand how protease isoforms contribute to specific cancer hallmarks. Here, in this review, we address cancer-specific isoforms produced by the alternative splicing of proteases and their distinctive roles in the tumor microenvironment.

In silico analysis of alternative splicing on drug-target gene interactions

Identifying and evaluating the right target are the most important factors in early drug discovery phase. Most studies focus on one protein ignoring the multiple splice-variant or protein-isoforms, which might contribute to unexpected therapeutic activity or adverse side effects. Here, we present computational analysis of cancer drug-target interactions affected by alternative splicing. By integrating information from publicly available databases, we curated 883 FDA approved or investigational stage small molecule cancer drugs that target 1,434 different genes, with an average of 5.22 protein isoforms per gene. Of these, 618 genes have ≥5 annotated protein-isoforms. By analyzing the interactions with binding pocket information, we found that 76% of drugs either miss a potential target isoform or target other isoforms with varied expression in multiple normal tissues. We present sequence and structure level alignments at isoform-level and make this information publicly available for all the curated drugs. Structure-level analysis showed ligand binding pocket architectures differences in size, shape and electrostatic parameters between isoforms. Our results emphasize how potentially important isoform-level interactions could be missed by solely focusing on the canonical isoform, and suggest that on- and off-target effects at isoform-level should be investigated to enhance the productivity of drug-discovery research.

Rockets, gauges, and pendulums: applying engineering principles to cell biology

From flight to radar to Velcro, biological form and function have inspired engineers for centuries. It is equally valuable to consider whether concepts in engineering might provide insights into core biological processes. To explore this idea, cell cycle checkpoints, biological clocks, and signaling pathways are viewed here from an engineering perspective. Engineering concepts covered include gauge error, the distinction between precision and accuracy, and the Taguchi method of robust design. Also discussed is the Pareto principle, which describes the observation that, in complex systems, a minority of the components (or inputs) are responsible for a majority of the outputs. These concepts enable engineers to manage complexity, both in system design and in operation. Thus, with new techniques and large data sets revealing ever-increasing levels of biological complexity, an engineering mindset may be particularly valuable for the study of living systems.

Discovery, characterization and potential roles of a novel NF-YAx splice variant in human neuroblastoma

BACKGROUND

Identification of novel cancer-associated splice variants is of potential diagnostic, prognostic and therapeutic importance. NF-Y transcription factor is comprised of NF-YA, NF-YB and NF-YC subunits, binds inverted CCAAT-boxes in ≈70% of gene promoters, regulates > 1000 cancer-associated genes and proteins involved in proliferation, staminality, differentiation, apoptosis, metabolism and is subject to component alternative splicing. RT-PCR evaluation of alternative NF-YA splicing in primary human neuroblastomas (NBs), led to discovery of a novel NF-YAx splice variant, also expressed during mouse embryo development and induced by doxorubicin in NB cells. Here, we report the discovery and characterisation of NF-YAx and discus its potential roles in NB.

METHODS

NF-YAx cDNA was RT-PCR-cloned from a stage 3 NB (provided by the Italian Association of Haematology and Paediatric Oncology, Genova, IT), sequenced and expressed as a protein using standard methods and compared to known fully-spliced NF-YAl and exon B-skipped NF-YAs isoforms in: EMSAs for capacity to form NF-Y complexes; by co-transfection, co-immunoprecipitation and Western blotting for capacity to bind Sp1; by IF for localisation; in AO/EtBr cell-death and colony formation assays for relative cytotoxicity, and by siRNA knockdown, use of inhibitors and Western blotting for potential mechanisms of action. Stable SH-SY5Y transfectants of all three NF-YA isoforms were also propagated and compared by RT-PCR and Western blotting for differences in cell-death and stem cell (SC)-associated gene expression, in cell-death assays for sensitivity to doxorubicin and in in vitro proliferation, substrate-independent growth and in vivo tumour xenograft assays for differences in growth and tumourigenic capacity.

RESULTS

NF-YAx was characterized as a novel variant with NF-YA exons B, D and partial F skipping, detected in 20% of NF-YA positive NBs, was the exclusive isoform in a stage 3 NB, expressed in mouse stage E11.5-14 embryos and induced by doxorubicin in SH-SY5Y NB cells. The NF-YAx protein exhibited nuclear localisation, competed with other isoforms in CCAAT box-binding NF-Y complexes but, in contrast to other isoforms, did not bind Sp1. NF-YAx expression in neural-related progenitor and NB cells repressed Bmi1 expression, induced KIF1Bβ expression and promoted KIF1Bβ-dependent necroptosis but in NB cells also selected tumourigenic, doxorubicin-resistant, CSC-like sub-populations, resistant to NF-YAx cytotoxicity.

CONCLUSIONS

The discovery of NF-YAx in NBs, its expression in mouse embryos and induction by doxorubicin in NB cells, unveils a novel NF-YA splice mechanism and variant, regulated by and involved in development, genotoxic-stress and NB. NF-YAx substitution of other isoforms in NF-Y complexes and loss of capacity to bind Sp1, characterises this novel isoform as a functional modifier of NF-Y and its promotion of KIF1Bβ-dependent neural-lineage progenitor and NB cell necroptosis, association with doxorubicin-induced necroptosis and expression in mouse embryos coinciding with KIF1Bβ-dependent sympathetic neuroblast-culling, confirm a cytotoxic function and potential role in suppressing NB initiation. On the other hand, the in vitro selection of CSC-like NB subpopulations resistant to NF-YAx cytotoxicity not only helps to explain high-level exclusive NF-YAx expression in a stage 3 NB but also supports a role for NF-YAx in disease progression and identifies a potential doxorubicin-inducible mechanism for post-therapeutic relapse.

Role of alternative splicing signatures in the prognosis of glioblastoma

BACKGROUND

Increasing evidence has validated the crucial role of alternative splicing (AS) in tumors. However, comprehensive investigations on the entirety of AS and their clinical value in glioblastoma (GBM) are lacking.

METHODS

The AS profiles and clinical survival data related to GBM were obtained from The Cancer Genome Atlas database. Univariate and multivariate Cox regression analyses were performed to identify survival-associated AS events. A risk score was calculated, and prognostic signatures were constructed using seven different types of independent prognostic AS events, respectively. The Kaplan-Meier estimator was used to display the survival of GBM patients. The receiver operating characteristic curve was applied to compare the predictive efficacy of each prognostic signature. Enrichment analysis and protein interactive networks were conducted using the gene symbols of the AS events to investigate important processes in GBM. A splicing network between splicing factors and AS events was constructed to display the potential regulatory mechanism in GBM.

RESULTS

A total of 2355 survival-associated AS events were identified. The splicing prognostic model revealed that patients in the high-risk group have worse survival rates than those in the low-risk group. The predictive efficacy of each prognostic model showed satisfactory performance; among these, the Alternate Terminator (AT) model showed the best performance at an area under the curve (AUC) of 0.906. Enrichment analysis uncovered that autophagy was the most enriched process of prognostic AS gene symbols in GBM. The protein network revealed that UBC, VHL, KCTD7, FBXL19, RNF7, and UBE2N were the core genes in GBM. The splicing network showed complex regulatory correlations, among which ELAVL2 and SYNE1_AT_78181 were the most correlated (r = -.506).

CONCLUSIONS

Applying the prognostic signatures constructed by independent AS events shows promise for predicting the survival of GBM patients. A splicing regulatory network might be the potential splicing mechanism in GBM.

Resistance Mechanisms to CAR T-Cell Therapy and Overcoming Strategy in B-Cell Hematologic Malignancies

Chimeric antigen receptor (CAR) T-cell therapy has shown promising clinical impact against hematologic malignancies. CD19 is a marker on the surface of normal B cells as well as most B-cell malignancies, and thus has a role as an effective target for CAR T-cell therapy. In numerous clinical data, successes with cell therapy have provided anticancer therapy as a potential therapeutic option for patients who are resistant to standard chemotherapies. However, recent growing evidence showed the limitations of the treatment such as antigen-positive relapse due to poor CAR T-cell persistence and antigen-negative relapses associated with CAR-driven mutations, alternative splicing, epitope masking, low antigen density, and lineage switching. The understanding of the resistance mechanisms to the cell therapy has developed novel potential treatment strategies, including dual-targeting therapy (dual and tandem CAR), and armored and universal CAR T-cell therapies. In this review, we provide an overview of resistance mechanisms to CD19 CAR T-cell therapy in B-cell malignancies and also review therapeutic strategies to overcome these resistances.

FJ, Liu Z, Tang W. Systematic profiling of alternative splicing events and splicing factors in left- and right-sided colon cancer

Left- and right-sided colon cancer (LC and RC) differ substantially in their molecular characteristics and prognoses, and are thus treated using different strategies. We systematically analyzed alternative splicing (AS) events and splicing factors in LC and RC. RNA-seq data were used for genome-wide profiling of AS events that could distinguish LC from RC. The Exon Skip splicing pattern was more common in RC, while the Retained Intron pattern was more common in LC. The AS events that were upregulated in RC were enriched for genes in the axon guidance pathway, while those that were upregulated in LC were enriched for genes in immune-related pathways. Prognostic models based on differentially expressed AS events were built, and a prognostic signature based on these AS events performed well for risk stratification in colon cancer patients. A correlation network of differentially expressed AS events and differentially expressed splicing factors was constructed, and RBM25 was identified as the hub gene in the network. In conclusion, large differences in AS events may contribute to the phenotypic differences between LC and RC. The differentially expressed AS events reported herein could be used as biomarkers and treatment targets for colon cancer.

Systematic profiling of a novel prognostic alternative splicing signature in hepatocellular carcinoma

Alternative splicing (AS) is a pervasive and vital mechanism involved in the progression of cancer by expanding genomic encoding capacity and increasing protein complexity. However, the systematic analysis of AS in hepatocellular carcinoma (HCC) is lacking and urgently required. In the present study, genome‑wide AS events with corresponding clinical information were profiled in 290 patients with HCC from the Cancer Genome Atlas and SpliceSeq software. Functional enrichment analyses revealed the pivotal biological process of AS regulation. Univariate Cox regression analyses were performed, followed by stepwise forward multivariate analysis to develop the prognostic signatures. Spearman's correlation analyses were also used to construct potential regulatory network between the AS events and aberrant splicing factors. A total of 34,163 AS events were detected, among which 1,805 AS events from 1,314 parent genes were significantly associated with the overall survival (OS) of patients with HCC, and their parent genes serve crucial roles in HCC‑related oncogenic processes, including the p53 signaling pathway, AMPK signaling pathway and HIF‑1 signaling pathway. A prognostic AS signature was established that was found to be an independent prognostic factor for OS in stratified cohorts, harboring a noteworthy ability to distinguish between the distinct prognoses of patients with HCC (high‑risk vs. low‑risk, 827 vs. 3,125 days, P<2e‑16). Time‑dependent receiver‑-operator characteristic curves confirmed its robustness and clinical efficacy, with the area under the curves maintained >0.9 for short‑term and long‑term prognosis prediction. The splicing correlation network suggested a trend in the interactions between splicing factors and prognostic AS events, further revealing the underlying mechanism of AS in the oncogenesis of HCC. In conclusion, the present study provides a comprehensive portrait of global splicing alterations involved in the progression and HCC in addition to valuable prognostic factors for patients, which may represent as underappreciated hallmark and provide novel clues of therapeutic targets in HCC.

Systematic Profiling of Alternative mRNA Splicing Signature for Predicting Glioblastoma Prognosis

Emerging evidence suggests that alternative splicing (AS) is modified in cancer and is associated with cancer progression. Systematic analysis of AS signature in glioblastoma (GBM) is lacking and is greatly needed. We profiled genome-wide AS events in 498 GBM patients in TCGA using RNA-seq data, and splicing network and prognostic predictor were built by integrated bioinformatics analysis. Among 45,610 AS events in 10,434 genes, we detected 1,829 AS events in 1,311 genes, and 1,667 AS events in 1,146 genes that were significantly associated with overall survival and disease-free survival of GBM patients, respectively. Five potential feature genes, S100A4, ECE2, CAST, ASPH, and LY6K, were discovered after network mining as well as correlation analysis between AS and gene expression, most of which were related to carcinogenesis and development. Multivariate survival model analysis indicated that these five feature genes could classify the prognosis at AS event and gene expression level. This report opens up a new avenue for exploration of the pathogenesis of GBM through AS, thus more precisely guiding clinical treatment and prognosis judgment.

Structural and functional characterisation of human RNA helicase DHX8 provides insights into the mechanism of RNA-stimulated ADP release

DHX8 is a crucial DEAH-box RNA helicase involved in splicing and required for the release of mature mRNA from the spliceosome. Here, we report the biochemical characterisation of full-length human DHX8 and the catalytically active helicase core DHX8Δ547, alongside crystal structures of DHX8Δ547 bound to ADP and a structure of DHX8Δ547 bound to poly(A)6 single-strand RNA. Our results reveal that DHX8 has an in vitro binding preference for adenine-rich RNA and that RNA binding triggers the release of ADP through significant conformational flexibility in the conserved DEAH-, P-loop and hook-turn motifs. We demonstrate the importance of R620 and both the hook-turn and hook-loop regions for DHX8 helicase activity and propose that the hook-turn acts as a gatekeeper to regulate the directional movement of the 3' end of RNA through the RNA-binding channel. This study provides an in-depth understanding of the activity of DHX8 and contributes insights into the RNA-unwinding mechanisms of the DEAH-box helicase family.

Downregulation of SRSF3 by antisense oligonucleotides sensitizes oral squamous cell carcinoma and breast cancer cells to paclitaxel treatment

PURPOSE

Paclitaxel (PTX) is widely used in the chemotherapy of many cancers, including breast cancer and oral squamous cell carcinoma (OSCC). However, many patients respond poorly to PTX treatment. The SRSF3 oncogene and several splicing factors play important roles in OSCC tumorigenesis. This study aimed to understand the function of splicing factors in PTX treatment and improve the therapeutic effects of PTX treatment.

METHODS

Splicing factors regulated by PTX treatment were screened in CAL 27 cell by reverse transcription polymerase chain reaction. The function of SRSF3 in PTX treatment was analyzed by gain-of-function or loss-of-function assay in OSCC cell lines CAL 27 and SCC-9 and breast cancer cell line MCF-7. Alternative splicing of SRSF3 exon 4 in cancer tissues or cells was analyzed by RT-PCR and online program TSVdb. SRSF3-specific antisense oligonucleotide (ASO) SR-3 was used to downregulate SRSF3 expression and enhance the effect of PTX treatment.

RESULTS

PTX treatment decreased SRSF3 expression, and SRSF3 overexpression rescued the growth inhibition caused by PTX in both OSCC and breast cancer cells. Moreover, we found that PTX treatment could repress SRSF3 exon 4 (containing an in-frame stop codon) exclusion and then decrease the SRSF3 protein expression. Increased exclusion of SRSF3 exon 4 is correlated with poor survival in OSCC and breast cancer patients. SR-3 downregulated SRSF3 protein expression and significantly increased the sensitivity of cancer cells to PTX treatment.

CONCLUSIONS

SRSF3 downregulation by ASO sensitizes cancer cells to PTX treatment.

Whole Genome Analysis and Prognostic Model Construction Based on Alternative Splicing Events in Endometrial Cancer

Objectives

A growing body of evidence has shown that aberrant alternative splicing (AS) is closely related to the occurrence and development of cancer. However, prior studies mainly have concentrated on a few genes that exhibit aberrant AS. This study aimed to determine AS events through whole genome analysis and construct a prognostic model of endometrial cancer (EC).

Methods

We downloaded gene expression RNAseq data from UCSC Xena, and seven types of AS events from TCGA SpliceSeq. Univariate Cox regression was employed to analyze the prognostic-related alternative splicing events (PASEs) and splicing factors; multivariate Cox regression was conducted to analyze the effect of risk score (All) and clinicopathological parameters on EC prognosis. An underlying interaction network of PASEs of EC was constructed by Cytoscape Reactome FI, GO, and KEGG pathway enrichment was performed by DAVID. ROC curves and Kaplan-Meir analysis were used to assess the diagnostic value of prognostic model. The correlation between PASEs and splicing factors was analyzed by GraphPad Prism; then a network was constructed using Cytoscape.

Results

In total, 28,281 AS events in EC were identified, which consisted of 1166 PASEs. RNPS1, NEK2, and CTNNB1 were the hub genes in the network of the top 600 PASEs. The area under the curve (AUC) of risk score (All) reached 0.819. Risk score (All) together with FIGO stage, cancer status, and primary therapy outcome success was risk factors of the prognosis of EC patients. Splicing factors YBX1, HNRNPDL, and HNRNPA1 were significantly related to the overall survival (OS). The splicing network indicated that the expression of splicing factors was significantly correlated with percent-splice-in (PSI) value of PASEs.

Conclusion

We constructed a model for predicting the prognosis of EC patients based on PASEs using whole genome analysis of AS events and thereby provided a reliable theoretical basis for EC clinical prognosis evaluation.

A Novel FGFR3 Splice Variant Preferentially Expressed in African American Prostate Cancer Drives Aggressive Phenotypes and Docetaxel Resistance

Alternative splicing (AS) has been shown to participate in prostate cancer development and progression; however, a link between AS and prostate cancer health disparities has been largely unexplored. Here we report on the cloning of a novel splice variant of FGFR3 that is preferentially expressed in African American (AA) prostate cancer. This novel variant (FGFR3-S) omits exon 14, comprising 123 nucleotides that encode the activation loop in the intracellular split kinase domain. Ectopic overexpression of FGFR3-S in European American (EA) prostate cancer cell lines (PC-3 and LNCaP) led to enhanced receptor autophosphorylation and increased activation of the downstream signaling effectors AKT, STAT3, and ribosomal S6 compared with FGFR3-L (retains exon 14). The increased oncogenic signaling imparted by FGFR3-S was associated with a substantial gain in proliferative and antiapoptotic activities, as well as a modest but significant gain in cell motility. Moreover, the FGFR3-S-conferred proliferative and motility gains were highly resistant to the pan-FGFR small-molecule inhibitor dovitinib and the antiapoptotic gain was insensitive to the cytotoxic drug docetaxel, which stands in marked contrast with dovitinib- and docetaxel-sensitive FGFR3-L. In an in vivo xenograft model, mice injected with PC-3 cells overexpressing FGFR3-S exhibited significantly increased tumor growth and resistance to dovitinib treatment compared with cells overexpressing FGFR3-L. In agreement with our in vitro and in vivo findings, a high FGFR3-S/FGFR3-L expression ratio in prostate cancer specimens was associated with poor patient prognosis. IMPLICATIONS: This work identifies a novel FGFR3 splice variant and supports the hypothesis that differential AS participates in prostate cancer health disparities.

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

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

microRNA-1908-5p inhibits proliferation and promotes apoptosis by targeting PP5 in NSCLC

MicroRNA-1908 is involved in the occurrence and development of various tumors. However, the mechanism of microRNA-1908-5p in the pathogenesis of non-small cell lung cancer (NSCLC) is not thoroughly studied. Protein phosphatase 5 catalytic subunit (PP5), a member of the protein phosphatase catalytic subunit family, may be a target of the microRNA-1908-5p. In order to further explore the mechanism of microRNA-1908-5p, real-time PCR was used to detect the expression of microRNA-1908-5p in non-small cell lung cancer tissues, and analyze the relationship between the expression of microRNA-1908-5p and clinical characteristics of lung cancer patients. The target of microRNA-1908-5p was predicted by bioinformatics and verified by Dual-luciferase assay. The effects of microRNA-1908-5p on the proliferation and apoptosis of lung cancer cells were examined at the cellular level. Nude mice tumorigenesis experiment was used to study the effect of microRNA-1908-5p on cancer cells. Western blot was used to detect the expression of related proteins. The results showed that the expression of microRNA-1908-5p in lung cancer tissues was significantly lower than that in adjacent tissues. The expression of microRNA-1908-5p in the non-metastatic lung cancer tissues was significantly higher than that in the metastatic lung cancer tissues, and the expression of microRNA-1908-5p was closely related to the survival rate of patients. Bioinformatics analysis combined with double luciferase assay showed that PP5 was a significant target of microRNA-1908-5p. Our results suggest that microRNA-1908-5p can regulate the pathogenesis of NSCLC by inhibiting PP5.

Signature microRNAs and long noncoding RNAs in laryngeal cancer recurrence identified using a competing endogenous RNA network

The aim of the present study was to identify novel microRNA (miRNA) or long noncoding RNA (lncRNA) signatures of laryngeal cancer recurrence and to investigate the regulatory mechanisms associated with this malignancy. Datasets of recurrent and nonrecurrent laryngeal cancer samples were downloaded from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus database (GSE27020 and GSE25727) to examine differentially expressed miRNAs (DE-miRs), lncRNAs (DE-lncRs) and mRNAs (DEGs). miRNA-mRNA and lncRNA-miRNA networks were constructed by investigating the associations among these RNAs in various databases. Subsequently, the interactions identified were combined into a competing endogenous RNA (ceRNA) regulatory network. Feature genes in the miRNA-mRNA network were identified via topological analysis and a recursive feature elimination algorithm. A support vector machine (SVM) classifier was established using the betweenness centrality values in the miRNA-mRNA network, consisting of 32 optimal feature-coding genes. The classification effect was tested using two validation datasets. Furthermore, coding genes in the ceRNA network were examined via pathway enrichment analyses. In total, 21 DE-lncRs, 507 DEGs and 55 DE-miRs were selected. The SVM classifier exhibited an accuracy of 94.05% (79/84) for sample classification prediction in the TCGA dataset, and 92.66 and 91.07% in the two validation datasets. The ceRNA regulatory network comprised 203 nodes, corresponding to mRNAs, miRNAs and lncRNAs, and 346 lines, corresponding to the interactions among RNAs. In particular, the interactions with the highest scores were HLA complex group 4 (HCG4)-miR-33b, HOX transcript antisense RNA (HOTAIR)-miR-1-MAGE family member A2 (MAGEA2), EMX2 opposite strand/antisense RNA (EMX2OS)-miR-124-calcitonin related polypeptide α (CALCA) and EMX2OS-miR-124-γ-aminobutyric acid type A receptor γ2 subunit (GABRG2). Gene enrichment analysis of the genes in the ceRNA network identified that 11 pathway terms and 16 molecular function terms were significantly enriched. The SVM classifier based on 32 feature coding genes exhibited high accuracy in the classification of laryngeal cancer samples. miR-1, miR-33b, miR-124, HOTAIR, HCG4 and EMX2OS may be novel biomarkers of recurrent laryngeal cancer, and HCG4-miR-33b, HOTAIR-miR-1-MAGEA2 and EMX2OS-miR-124-CALCA/GABRG2 may be associated with the molecular mechanisms regulating recurrent laryngeal cancer.

Targeting Chemoresistant Tumors: Could TRIM Proteins-p53 Axis Be a Possible Answer?

Chemosensitivity is a crucial feature for all tumours so that they can be successfully treated, but the huge heterogeneity of these diseases, to be intended both inter- and intra-tumour, makes it a hard-to-win battle. Indeed, this genotypic and phenotypic variety, together with the adaptability of tumours, results in a plethora of chemoresistance acquisition mechanisms strongly affecting the effectiveness of treatments at different levels. Tripartite motif (TRIM) proteins are shown to be involved in some of these mechanisms thanks to their E3-ubiquitin ligase activity, but also to other activities they can exert in several cellular pathways. Undoubtedly, the ability to regulate the stability and activity of the p53 tumour suppressor protein, shared by many of the TRIMs, represents the preeminent link between this protein family and chemoresistance. Indeed, they can modulate p53 degradation, localization and subset of transactivated target genes, shifting the cellular response towards a cytoprotective or cytotoxic reaction to whatever damage induced by therapy, sometimes in a cellular-dependent way. The involvement in other chemoresistance acquisition mechanisms, independent by p53, is known, affecting pivotal processes like PI3K/Akt/NF-κB signalling transduction or Wnt/beta catenin pathway, to name a few. Hence, the inhibition or the enhancement of TRIM proteins functionality could be worth investigating to better understand chemoresistance and as a strategy to increase effectiveness of anticancer therapies.

Alternative splicing, RNA-seq and drug discovery

Alternative splicing, hereafter referred to as AS, is an essential component of gene expression regulation that contributes to the diversity of proteomes. Recent developments in RNA sequencing (RNA-seq) technologies, combined with the advent of computational tools, have enabled transcriptome-wide studies of AS at an unprecedented scale and resolution. RNA mis-splicing can cause human disease, and to target alternative splicing has led to the development of novel therapeutics. Splice variants diversify the repertoire of biomarkers and functionally contribute to drug resistance. Our expanding knowledge of AS variation in human populations holds great promise for improving disease diagnoses and ultimately patient care in the era of sequencing and precision medicine.

CircSMARCA5 Regulates VEGFA mRNA Splicing and Angiogenesis in Glioblastoma Multiforme Through the Binding of SRSF1

Circular RNAs are a large group of RNAs whose cellular functions are still being investigated. We recently proposed that circSMARCA5 acts as sponge for the splicing factor Serine and Arginine Rich Splicing Factor 1 (SRSF1) in glioblastoma multiforme (GBM). After demonstrating by RNA immunoprecipitation a physical interaction between SRFS1 and circSMARCA5, we assayed by real-time PCR in a cohort of 31 GBM biopsies and 20 unaffected brain parenchyma controls (UC) the expression of total, pro-angiogenic (Iso8a) and anti-angiogenic (Iso8b) mRNA isoforms of Vascular Endothelial Growth Factor A (VEGFA), a known splicing target of SRSF1. The Iso8a to Iso8b ratio: (i) increased in GBM biopsies with respect to UC (p-value < 0.00001); (ii) negatively correlated with the expression of circSMARCA5 (r-value = −0.46, p-value = 0.006); (iii) decreased in U87-MG overexpressing circSMARCA5 with respect to negative control (p-value = 0.0055). Blood vascular microvessel density, estimated within the same biopsies, negatively correlated with the expression of circSMARCA5 (r-value = −0.59, p-value = 0.00001), while positively correlated with that of SRSF1 (r-value = 0.38, p-value = 0.00663) and the Iso8a to Iso8b ratio (r-value = 0.41, p-value = 0.0259). Kaplan-Meier survival analysis showed that GBM patients with low circSMARCA5 expression had lower overall and progression free survival rates than those with higher circSMARCA5 expression (p-values = 0.033, 0.012, respectively). Our data convincingly suggest that circSMARCA5 is an upstream regulator of pro- to anti-angiogenic VEGFA isoforms ratio within GBM cells and a highly promising GBM prognostic and prospective anti-angiogenic molecule.